Question: extract SNPs of certain gene
1
gravatar for marwa.cs91
3.4 years ago by
marwa.cs9120
Egypt
marwa.cs9120 wrote:

Hello all,I am new in plink analysis

I have Map file and I know that location (start /end) of CD2AP gene on chromosome 6 has base pairs 47,477,745 to 47,627,262

How can I use physical position column in MAP file to extract SNPs of this gene 

 

Any help will be appreciated

map location plink snps gene • 1.9k views
ADD COMMENTlink modified 3.4 years ago by chrchang5234.9k • written 3.4 years ago by marwa.cs9120
1
gravatar for Ibrahim Tanyalcin
3.4 years ago by
Belgium
Ibrahim Tanyalcin930 wrote:

Dear,

First you need to make sure which transcript you are interested. Keep in mind that the same SNP will have different cDNA coordinates between transcripts. Also you need to watch out for the genome version. Because g coordinates will also vary between different genome versions. Moving from your example, you can go to www.ensembl.org and write in the search CD2AP. You will see that the golden transcript is ENST00000359314 with 629 aminoacids. Now go to biomart at www.ensembl.org/biomart . Here select ensembl variation database, than select SNPs and Indels. Next from the left click on filters and go to gene associated filters. Here paste your gene ID which is ENSG00000198087. Next go to the attributes tab and select the coordinates you want (genomic, mRNA, coding sequence..etc) Download your results and you will have your SNPs. 

I hope this helps,

ADD COMMENTlink written 3.4 years ago by Ibrahim Tanyalcin930

Thanks alot for your help I really appreciate

ADD REPLYlink written 3.4 years ago by marwa.cs9120
1
gravatar for chrchang523
3.4 years ago by
chrchang5234.9k
United States
chrchang5234.9k wrote:

plink cannot process isolated .map files; do you have an accompanying .ped?  (If you don't, it might be simpler to write a short script for this job.)

If you do, "plink --file [.ped/.map filename prefix] --chr 6 --from-bp 47477745 --to-bp 47627262 --write-snplist --out [output filename prefix]" should give you the list you're looking for.

ADD COMMENTlink modified 3.4 years ago • written 3.4 years ago by chrchang5234.9k

Thanks alot for your help I really appreciate

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by marwa.cs9120
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1691 users visited in the last hour