How To Get A List Of All Correlated Snps (In Ld) With Another List Of Snps?
6
1
Entering edit mode
9.6 years ago

Hello. I have a list of 54 tagging SNPs. And I would like to get a list of all SNPs in LD with each one of my 54 SNPs. I think I need some kind of software that would allow me to get all correlated SNPs using r2>=0.6 or 0.8.. I cannot find a way to do this sistematically for 54 SNPs.. Can anyone help? Thanks

snp linkage haplotype • 12k views
ADD COMMENT
3
Entering edit mode
6.6 years ago
cedlund ▴ 30

rAggr can find SNPs and indels in LD with your list of SNPs. See: http://raggr.usc.edu/

It now has data from 1000 Genomes Phase 3 and can search up to 5Mb in each direction from the queried SNP.

ADD COMMENT
0
Entering edit mode

I am getting "A server error occurred. Try submitting again. If you still receive an error, contact the site administrator" while trying to query LD for marker SNP's .

ADD REPLY
0
Entering edit mode

Were you able to resolve this ?

ADD REPLY
0
Entering edit mode

Hi, I used this https://analysistools.nci.nih.gov/LDlink/ for calculating LD. Thanks

ADD REPLY
2
Entering edit mode
8.7 years ago

SNAP Proxy search is the one I have been using, it is easy to use, they should update to hg19!

ADD COMMENT
2
Entering edit mode
8.7 years ago
Erik Garrison ★ 2.3k

If you can pull the SNPs out of the 1000 Genomes VCF, you can then run vcf-tools --hap-r2 or something suitable. If the SNPs are close together, you can use the fully-phased VCF and --hap-r2. Otherwise you may need to use --geno-r2, which is probably what you want anyway. Good luck.

ADD COMMENT
1
Entering edit mode
9.6 years ago

PLINK can do this:

"to obtain all LD values from a group of SNPs with other SNPs, use the command --ld-snp-list mysnps.txt where mysnps.txt is a list of SNPs."

http://pngu.mgh.harvard.edu/~purcell/plink/ld.shtml

ADD COMMENT
0
Entering edit mode
9.6 years ago
Joey ▴ 430

May be SNAP (http://www.broadinstitute.org/mpg/snap/ldsearch.php) is the tool you need.

ADD COMMENT
0
Entering edit mode
9.6 years ago

The problem using PLINK is that you need to have the full list of genotypes for all possible SNPs. From the question (correct me if I am wrong) I assume you have several SNP names (probably with rs as they are tag SNPs so they should be known and previously genotyped), and you want a list of possible correlated SNPs in certain population similar to yours, and you have not genotyped them yet.

For this kind of task you have several option, the ones I have used before when I was at Sanger are:

GLIDERS (unfortunately it is not up to date with 1000 genomes, only HAPMAP projects)
[http://www.sanger.ac.uk/resources/software/gliders/][1]

EnsemblAPI (I think that you can do that now in the web, probably you could ask the help-desk team, they are very helpful)
[http://www.ensembl.org/info/docs/api/variation/variation_tutorial.html#ld][2]
ADD COMMENT

Login before adding your answer.

Traffic: 2277 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