This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Get Coding Snps For 1000 Genomes Data

Edit. How can I get the coding SNPs (coding synonymous and non-synonymous polymorphisms) for a particular gene (say, BRCA1) from 1000 genomes data?

BioMart provides a nice interface to 1000 genomes data, but it just takes forever to output. Is there an alternative where I can simply look for the coding SNPs by providing the gene id or the corresponding genomic region?

Thanks

genome gene snp retrieval
• 5,320 views • link
updated by Biostar Biostar • written by Korban • Korban

Why would you want to retrieve this from 1000 Genomes data rather than Ensembl? Do you want to retrieve reads mapping to a gene or the sequence itself?

• 0 views • link
written by User 59 User 59

Hi Daniel, I want to check for the polymorphisms in a gene of interest.

• 0 views • link
written by Korban • Korban

2 answers

To find polymorphisms in a gene of interest, I would get the gene's general region as chromosome and start/stop location (plus a kilobase or ten on each end for regulatory regions). Then go to Ensembl or Biomart or dbSNP and ask for known variants in that region.

• 0 views • link
written by Karl • Karl

You can also search by gene name / symbol in Ensembl (e.g. 'BRCA2'). To get a list of the variants in a gene go to the 'Variant Table' page (e.g. http://www.ensembl.org/Homo_sapiens/Gene/Variation_Gene/Table?g=ENSG00000139618;r=13:32889611-32973805). You can configure the page using [Configure this page] in the side menu. Note that by default we show only variants that are located in the exons and up till 100 bp away from the exons. To see all intronic variants as well as the variants up to 5 kb up- and downstream of the gene, set 'Intron Context' on the configuration page to 'Full Introns'.

• 0 views • link
updated by Ram Ram • written by Bert Overduin Bert Overduin

You could use one of the SIFT tools to filter out non-coding variants. The name of the tool is "Restrict to coding variants" here is the link http://sift.bii.a-star.edu.sg/www/SIFT_intersect_coding_submit.html

Remember to choose the genome build and the input format. It can take in multiple formats like VCF4, pileup, MAQ or SIFT's own simple comma delimited formats.

• 2 views • link
updated by Ram Ram • written by Prateek Prateek

Log in to answer this question.

More posts like this