This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Snps In Linkage Disequilibrium For Ncbi Build 37 (Hg19)

Hi all

I'd like to know if anyone has or might know where to get a list of SNPs in high LD for the last human genome build (GRCh37). I'm currently working with the Brazilian population. I'm trying to do PCA for population structure analysis and, therefore, need to remove high LD regions from my dataset to get a low biased PC axes.

Thanks a lot in advance for the attention!

Best regards, Pedro. PS: Just --indep-pairwise option in PLINK won't solve my problem

ld snp gwas plink
• 7,142 views • link
updated by Biostar Biostar • written by cruzpedro cruzpedro

2 answers

The LD data for hapmap / hg18 are available here: http://hapmap.ncbi.nlm.nih.gov/downloads/ld_data/?N=D

you could run liftOver to map those positions to hg19.

• 47 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thanks a lot for you answer, really shed some light to my problem!

• 1 views • link
written by cruzpedro cruzpedro

Do you know if there is the same data from the newest release of the 1000G project?

• 0 views • link
written by Jimbou Jimbou

Thank you very much. But is there also a ftp or whole genome download site?

• 0 views • link
written by Jimbou Jimbou

Had a chat with my colleague Laura Clarke from the 1000 Genomes here at EMBL-EBI and this is what she said

There are no bulk downloads for LD information. Doing pairwise comparisons of 80m sites between 2500 individuals is impossible. The thing to recommend is to convert the files to plink format using our VCF to PED tool or vcftools

http://www.1000genomes.org/faq/can-i-convert-vcf-files-plinkped-format

http://vcftools.sourceforge.net/documentation.html#plink

Then look at the region of interest in haploview or equivalent. There is no feasible way to get this in bulk for any large quantity of 1000 genomes sites or individuals

• 0 views • link
updated by Ram Ram • written by Denise CS Denise CS

Well I guess what is true for other populations might not be same for Brazilian population, therefore you may need to calculate the pairwise LD in you sample set and construct LD blocks ( and scan the consecutive window to locate the high LD region) Please refer to this paper for more information: A Price et a;., Long-Range LD Can Confound Genome Scansin Admixed Populations, AJHG 2008 (https://www.sciencedirect.com/science/article/pii/S0002929708003534)

• 49 views • link
written by Tky Tky

Hey! Thanks a lot for your answer! I was planning to do this as well later, since I was first geting a "test" set for running EIGENSTRAT. Well, I think I'll use for my sample the "--ld-window-kb" option in PLINK: http://pngu.mgh.harvard.edu/~purcell/plink/ld.shtml#ld2. Do you think this might suffice for LD controlling in a 60 individuals dataset? I have some data on other 50 individuals I can include to infer patterns of LD in my population.

Best regards!

• 1 views • link
written by cruzpedro cruzpedro

Yeah, you should check the eigenstrat plot first (against other HapMap populations) And more samples give more accurate LD estimations.

• 1 views • link
written by Tky Tky

Log in to answer this question.

More posts like this