Modeling linkage disequilibrium and identifying recombination hotspots using single-nucleotide polymorphism data.
We introduce a new statistical model for patterns of linkage disequilibrium (LD) among multiple SNPs in a population sample. The model overcomes limitations of existing approaches to understanding, summarizing, and interpreting LD by (i) relating patterns of LD directly to the underlying recombination process; (ii) considering all loci simultaneously, rather than pairwise; (iii) avoiding the assumption that LD necessarily has a "block-like" structure; and (iv) being computationally tractable for huge genomic regions (up to complete chromosomes). We examine in detail one natural application of the model: estimation of underlying recombination rates from population data. Using simulation, we show that in the case where recombination is assumed constant across the region of interest, recombination rate estimates based on our model are competitive with the very best of current available methods. More importantly, we demonstrate, on real and simulated data, the potential of the model to help identify and quantify fine-scale variation in recombination rate from population data. We also outline how the model could be useful in other contexts, such as in the development of more efficient haplotype-based methods for LD mapping.
Li
Na
N
Department of Biostatistics, University of Washington, Seattle, Washington 98195, USA.
Stephens
Matthew
M
D000071450:ATP Binding Cassette Transporter, Subfamily B, Member 3; D018528:ATP-Binding Cassette Transporters / Q000235:genetics; D000465:Algorithms; D002874:Chromosome Mapping; D003198:Computer Simulation; D018740:Genetic Heterogeneity; D005828:Genetics, Population; D006239:Haplotypes; D006801:Humans; D015810:Linkage Disequilibrium; D008957:Models, Genetic; D020641:Polymorphism, Single Nucleotide; D011156:Population Density; D011995:Recombination, Genetic
D000071450:ATP Binding Cassette Transporter, Subfamily B, Member 3; D018528:ATP-Binding Cassette Transporters; C071692:TAP2 protein, human