How to obtain the chromosomal location for a list of 8M SNPs?
0
0
Entering edit mode
3.8 years ago
aroso491 • 0

Hello, I have a list of 8M SNPs from a meta-GWAS which does not include the chromosome ID and chromosome location in bp. I thought I could do this using the grch37.ensembl.org mart of bioMart (I have never done this before and I am still learning), but I have read in other posts that there is a 500 line limit.

My list of rsIDs is of about 8.5M and I have read in other posts that there are other ways to do this like Parse the VCF, Use the APIs or Use the VEP.

The thing is that I do not know how to do any of these and I have not managed to find any sort of tutorial or anything to learn to do that. Do any of you have any sort of indication about where can I learn to retrieve the chromosomal location and ID for my list, please?

Thanks in advance.

VCF SNP R ensembl bioMart • 1.2k views
ADD COMMENT
2
Entering edit mode

An alternate way would be to use galaxy https://usegalaxy.org/

  • In the left tools search bar searching for "UCSC"
  • "UCSC Main table browser" click redirects to ucsc table browser
  • Here group "variation" and dbSNP table can be selected
  • Then in "get output" -> "send query to galaxy"
  • The job might take a while but the the bed of dbSNP with rsids will be uploaded to galaxy
  • Then one can upload the file with rsids and use "join two datasets side by side" tool to get the join of the dnsnp location of the uploaded SNPs
ADD REPLY
0
Entering edit mode

Thank you so much! This sounds exactly like what I needed. I am giving it a try right now, will let you know if it worked. :)

Edit: It worked just fine, although it's only been able to retrieve 3.5 million out of my 8 million file. I assume this is normal as some of the SNPs might just not be in the database or identified yet, but I was wondering if there'd be a way to double check this.

ADD REPLY
1
Entering edit mode
ADD REPLY

Login before adding your answer.

Traffic: 2713 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6