Genetics, Vol. 165, 2213-2233, December 2003, Copyright © 2003

Modeling Linkage Disequilibrium and Identifying Recombination Hotspots Using Single-Nucleotide Polymorphism Data

Na Lia and Matthew Stephensb
a Department of Biostatistics, University of Washington, Seattle, Washington 98195
b Department of Statistics, University of Washington, Seattle, Washington 98195

Corresponding author: Matthew Stephens, University of Washington, Padelford B-313, Stevens Way, Seattle, WA 98195., stephens{at}stat.washington.edu (E-mail)

Communicating editor: S. TAVARE

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.





This article has been cited by other articles:


Home page
BioinformaticsHome page
B. Pasaniuc, S. Sankararaman, G. Kimmel, and E. Halperin
Inference of locus-specific ancestry in closely related populations
Bioinformatics, June 15, 2009; 25(12): i213 - i221.
[Abstract] [Full Text] [PDF]


Home page
BioinformaticsHome page
S. Kim, K.-A. Sohn, and E. P. Xing
A multivariate regression approach to association analysis of a quantitative trait network
Bioinformatics, June 15, 2009; 25(12): i204 - i212.
[Abstract] [Full Text] [PDF]


Home page
BioinformaticsHome page
J. Yin, M. I. Jordan, and Y. S. Song
Joint estimation of gene conversion rates and mean conversion tract lengths from population SNP data
Bioinformatics, June 15, 2009; 25(12): i231 - i239.
[Abstract] [Full Text] [PDF]


Home page
Eur Heart JHome page
U. H. Frey, M. Adamzik, E. Kottenberg-Assenmacher, H. Jakob, I. Manthey, M. Broecker-Preuss, L. Bergmann, G. Heusch, W. Siffert, J. Peters, et al.
A novel functional haplotype in the human GNAS gene alters G{alpha}s expression, responsiveness to {beta}-adrenoceptor stimulation, and peri-operative cardiac performance
Eur. Heart J., June 1, 2009; 30(11): 1402 - 1410.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
K. E. Lohmueller, C. D. Bustamante, and A. G. Clark
Methods for Human Demographic Inference Using Haplotype Patterns From Genomewide Single-Nucleotide Polymorphism Data
Genetics, May 1, 2009; 182(1): 217 - 231.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
Y. Wang and B. Rannala
Population genomic inference of recombination rates and hotspots
PNAS, April 14, 2009; 106(15): 6215 - 6219.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
V. C. Sousa, M. Fritz, M. A. Beaumont, and L. Chikhi
Approximate Bayesian Computation Without Summary Statistics: The Case of Admixture
Genetics, April 1, 2009; 181(4): 1507 - 1519.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
S. Boitard, C. Schlotterer, and A. Futschik
Detecting Selective Sweeps: A New Approach Based on Hidden Markov Models
Genetics, April 1, 2009; 181(4): 1567 - 1578.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
G. K. Chen, P. Marjoram, and J. D. Wall
Fast and flexible simulation of DNA sequence data
Genome Res., January 1, 2009; 19(1): 136 - 142.
[Abstract] [Full Text] [PDF]


Home page
Phil Trans R Soc BHome page
Y. Wang and B. Rannala
Bayesian inference of fine-scale recombination rates using population genomic data
Phil Trans R Soc B, December 27, 2008; 363(1512): 3921 - 3930.
[Abstract] [Full Text] [PDF]


Home page
Mol Biol EvolHome page
B. C. Verrelli, C. M. Lewis Jr, A. C. Stone, and G. H. Perry
Different Selective Pressures Shape the Molecular Evolution of Color Vision in Chimpanzee and Human Populations
Mol. Biol. Evol., December 1, 2008; 25(12): 2735 - 2743.
[Abstract] [Full Text] [PDF]


Home page
Schizophr BullHome page
J. van Os, B. P. Rutten, and R. Poulton
Gene-Environment Interactions in Schizophrenia: Review of Epidemiological Findings and Future Directions
Schizophr Bull, November 1, 2008; 34(6): 1066 - 1082.
[Abstract] [Full Text] [PDF]


Home page
Proc R Soc BHome page
J. A Andres, L. S Maroja, and R. G Harrison
Searching for candidate speciation genes using a proteomic approach: seminal proteins in field crickets
Proc R Soc B, September 7, 2008; 275(1646): 1975 - 1983.
[Abstract] [Full Text] [PDF]


Home page
Alcohol AlcoholHome page
C. A. Hodgkinson, Q. Yuan, K. Xu, P.-H. Shen, E. Heinz, E. A. Lobos, E. B. Binder, J. Cubells, C. L. Ehlers, J. Gelernter, et al.
Addictions Biology: Haplotype-Based Analysis for 130 Candidate Genes on a Single Array
Alcohol Alcohol., September 1, 2008; 43(5): 505 - 515.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
P. F. O'Reilly, E. Birney, and D. J. Balding
Confounding between recombination and selection, and the Ped/Pop method for detecting selection
Genome Res., August 1, 2008; 18(8): 1304 - 1313.
[Abstract] [Full Text] [PDF]


Home page
Clin. Chem.Home page
J. A. G. Agundez, K. Golka, C. Martinez, S. Selinski, M. Blaszkewicz, and E. Garcia-Martin
Unraveling Ambiguous NAT2 Genotyping Data
Clin. Chem., August 1, 2008; 54(8): 1390 - 1394.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
J.-F. Lefebvre and D. Labuda
Fraction of Informative Recombinations: A Heuristic Approach to Analyze Recombination Rates
Genetics, April 1, 2008; 178(4): 2069 - 2079.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
I. J. Tsai, D. Bensasson, A. Burt, and V. Koufopanou
Population genomics of the wild yeast Saccharomyces paradoxus: Quantifying the life cycle
PNAS, March 25, 2008; 105(12): 4957 - 4962.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
A. Ojeda, L.-S. Huang, J. Ren, A. Angiolillo, I.-C. Cho, H. Soto, C. Lemus-Flores, S. M. Makuza, J. M. Folch, and M. Perez-Enciso
Selection in the Making: A Worldwide Survey of Haplotypic Diversity Around a Causative Mutation in Porcine IGF2
Genetics, March 1, 2008; 178(3): 1639 - 1652.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
J.-B. Veyrieras, L. Camus-Kulandaivelu, B. Gouesnard, D. Manicacci, and A. Charcosset
Bridging Genomics and Genetic Diversity: Linkage Disequilibrium Structure and Association Mapping in Maize and Other Cereals
Crop Sci., December 18, 2007; 47(Supplement_3): S-60 - S-71.
[Abstract] [Full Text] [PDF]


Home page
Mol Biol EvolHome page
F. Rousset and R. Leblois
Likelihood and Approximate Likelihood Analyses of Genetic Structure in a Linear Habitat: Performance and Robustness to Model Mis-Specification
Mol. Biol. Evol., December 1, 2007; 24(12): 2730 - 2745.
[Abstract] [Full Text] [PDF]


Home page
J HeredHome page
C. M. Seabury, C. A. Gill, J. W. Templeton, J. B. Dyar, J. N. Derr, D. L. Adelson, E. Owens, D. S. Davis, D. C. Kraemer, and J. E. Womack
Molecular Characterization of the Rocky Mountain Elk (Cervus elaphus nelsoni) PRNP Putative Promoter
J. Hered., November 21, 2007; (2007) esm091v1.
[Abstract] [Full Text] [PDF]


Home page
BioinformaticsHome page
B. Kirkpatrick, C. S. Armendariz, R. M. Karp, and E. Halperin
HAPLOPOOL: improving haplotype frequency estimation through DNA pools and phylogenetic modeling
Bioinformatics, November 15, 2007; 23(22): 3048 - 3055.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
J. Gay, S. Myers, and G. McVean
Estimating Meiotic Gene Conversion Rates From Population Genetic Data
Genetics, October 1, 2007; 177(2): 881 - 894.
[Abstract] [Full Text] [PDF]


Home page
BioinformaticsHome page
F. A. Wright, H. Huang, X. Guan, K. Gamiel, C. Jeffries, W. T. Barry, F. Pardo-Manuel de Villena, P. F. Sullivan, K. C. Wilhelmsen, and F. Zou
Simulating association studies: a data-based resampling method for candidate regions or whole genome scans
Bioinformatics, October 1, 2007; 23(19): 2581 - 2588.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
A. Auton and G. McVean
Recombination rate estimation in the presence of hotspots
Genome Res., August 1, 2007; 17(8): 1219 - 1227.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
T. H. E. Meuwissen and M. E. Goddard
Multipoint Identity-by-Descent Prediction Using Dense Markers to Map Quantitative Trait Loci and Estimate Effective Population Size
Genetics, August 1, 2007; 176(4): 2551 - 2560.
[Abstract] [Full Text] [PDF]


Home page
BiostatisticsHome page
T. Johnson
Bayesian method for gene detection and mapping, using a case and control design and DNA pooling
Biostat., July 1, 2007; 8(3): 546 - 565.
[Abstract] [Full Text] [PDF]


Home page
BioinformaticsHome page
K.-A. Sohn and E. P. Xing
Spectrum: joint bayesian inference of population structure and recombination events
Bioinformatics, July 1, 2007; 23(13): i479 - i489.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
A. RoyChoudhury and M. Stephens
Fast and Accurate Estimation of the Population-Scaled Mutation Rate, {theta}, From Microsatellite Genotype Data
Genetics, June 1, 2007; 176(2): 1363 - 1366.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
P. Calabrese
A population genetics model with recombination hotspots that are heterogeneous across the population
PNAS, March 13, 2007; 104(11): 4748 - 4752.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
V. J. Clark, S. E. Ptak, I. Tiemann, Y. Qian, G. Coop, A. C. Stone, M. Przeworski, N. Arnheim, and A. D. Rienzo
Combining Sperm Typing and Linkage Disequilibrium Analyses Reveals Differences in Selective Pressures or Recombination Rates Across Human Populations
Genetics, February 1, 2007; 175(2): 795 - 804.
[Abstract] [Full Text] [PDF]


Home page
J. Pharmacol. Exp. Ther.Home page
S. B. Koukouritaki, M. T. Poch, M. C. Henderson, L. K. Siddens, S. K. Krueger, J. E. VanDyke, D. E. Williams, N. M. Pajewski, T. Wang, and R. N. Hines
Identification and Functional Analysis of Common Human Flavin-Containing Monooxygenase 3 Genetic Variants
J. Pharmacol. Exp. Ther., January 1, 2007; 320(1): 266 - 273.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
S. Eyheramendy, J. Marchini, G. McVean, S. Myers, and P. Donnelly
A model-based approach to capture genetic variation for future association studies
Genome Res., January 1, 2007; 17(1): 88 - 95.
[Abstract] [Full Text] [PDF]


Home page
BioinformaticsHome page
P. Fearnhead
SequenceLDhot: detecting recombination hotspots
Bioinformatics, December 15, 2006; 22(24): 3061 - 3066.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
A. Ojeda, J. Rozas, J. M. Folch, and M. Perez-Enciso
Unexpected High Polymorphism at the FABP4 Gene Unveils a Complex History for Pig Populations
Genetics, December 1, 2006; 174(4): 2119 - 2127.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
B. Padhukasahasram, J. D. Wall, P. Marjoram, and M. Nordborg
Estimating Recombination Rates From Single-Nucleotide Polymorphisms Using Summary Statistics
Genetics, November 1, 2006; 174(3): 1517 - 1528.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
Y.-K. Yoo, X. Ke, S. Hong, H.-Y. Jang, K. Park, S. Kim, T. Ahn, Y.-D. Lee, O. Song, N.-Y. Rho, et al.
Fine-Scale Map of Encyclopedia of DNA Elements Regions in the Korean Population
Genetics, September 1, 2006; 174(1): 491 - 497.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
R. Jiang, P. Marjoram, J. O. Borevitz, and S. Tavare
Inferring Population Parameters From Single-Feature Polymorphism Data
Genetics, August 1, 2006; 173(4): 2257 - 2267.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
P. M. Visscher and W. G. Hill
Estimation of Recombination Rate and Detection of Recombination Hotspots From Dense Single-Nucleotide Polymorphism Trio Data
Genetics, August 1, 2006; 173(4): 2415 - 2417.
[Full Text] [PDF]


Home page
GeneticsHome page
P. L. Morrell, D. M. Toleno, K. E. Lundy, and M. T. Clegg
Estimating the Contribution of Mutation, Recombination and Gene Conversion in the Generation of Haplotypic Diversity
Genetics, July 1, 2006; 173(3): 1705 - 1723.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
V. Plagnol, B. Padhukasahasram, J. D. Wall, P. Marjoram, and M. Nordborg
Relative Influences of Crossing Over and Gene Conversion on the Pattern of Linkage Disequilibrium in Arabidopsis thaliana
Genetics, April 1, 2006; 172(4): 2441 - 2448.
[Abstract] [Full Text] [PDF]


Home page
BioinformaticsHome page
M. Perez-Enciso
Multiple association analysis via simulated annealing (MASSA)
Bioinformatics, March 1, 2006; 22(5): 573 - 580.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
D. J. Wilson and G. McVean
Estimating Diversifying Selection and Functional Constraint in the Presence of Recombination
Genetics, March 1, 2006; 172(3): 1411 - 1425.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
F. A. Reed and S. A. Tishkoff
Positive Selection Can Create False Hotspots of Recombination
Genetics, March 1, 2006; 172(3): 2011 - 2014.
[Abstract] [Full Text] [PDF]


Home page
HypertensionHome page
E. L. Riddle, B. K. Rana, K. K. Murthy, F. Rao, E. Eskin, D. T. O'Connor, and P. A. Insel
Polymorphisms and Haplotypes of the Regulator of G Protein Signaling-2 Gene in Normotensives and Hypertensives
Hypertension, March 1, 2006; 47(3): 415 - 420.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
D. M. Greenawalt, X. Cui, Y. Wu, Y. Lin, H.-Y. Wang, M. Luo, I. V. Tereshchenko, G. Hu, J. Y. Li, Y. Chu, et al.
Strong correlation between meiotic crossovers and haplotype structure in a 2.5-Mb region on the long arm of chromosome 21
Genome Res., February 1, 2006; 16(2): 208 - 214.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
N. G. C. Smith and P. Fearnhead
A Comparison of Three Estimators of the Population-Scaled Recombination Rate: Accuracy and Robustness
Genetics, December 1, 2005; 171(4): 2051 - 2062.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
G. M. Clarke and L. R. Cardon
Disentangling Linkage Disequilibrium and Linkage From Dense Single-Nucleotide Polymorphism Trio Data
Genetics, December 1, 2005; 171(4): 2085 - 2095.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
A. V. Smith, D. J. Thomas, H. M. Munro, and G. R. Abecasis
Sequence features in regions of weak and strong linkage disequilibrium
Genome Res., November 1, 2005; 15(11): 1519 - 1534.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
P. Hallast, L. Nagirnaja, T. Margus, and M. Laan
Segmental duplications and gene conversion: Human luteinizing hormone/chorionic gonadotropin {beta} gene cluster
Genome Res., November 1, 2005; 15(11): 1535 - 1546.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
M. T. Hamblin, M. G. Salas Fernandez, A. M. Casa, S. E. Mitchell, A. H. Paterson, and S. Kresovich
Equilibrium Processes Cannot Explain High Levels of Short- and Medium-Range Linkage Disequilibrium in the Domesticated Grass Sorghum bicolor
Genetics, November 1, 2005; 171(3): 1247 - 1256.
[Abstract] [Full Text] [PDF]


Home page
J HeredHome page
E. Perez-Enciso, P. G. Garcia-Bernal, and M. Perez-Enciso
Atlas: A Java-Based Tool for Managing Genotypes
J. Hered., September 1, 2005; 96(5): 623 - 625.
[Abstract] [Full Text] [PDF]


Home page
Phil Trans R Soc BHome page
R. W Lawrence, D. M Evans, and L. R Cardon
Prospects and pitfalls in whole genome association studies
Phil Trans R Soc B, August 29, 2005; 360(1460): 1589 - 1595.
[Abstract] [Full Text] [PDF]


Home page
Phil Trans R Soc BHome page
M. De Iorio, E. de Silva, and M. P.H Stumpf
Recombination hotspots as a point process
Phil Trans R Soc B, August 29, 2005; 360(1460): 1597 - 1603.
[Abstract] [Full Text] [PDF]


Home page
Phil Trans R Soc BHome page
G. A.T McVean and N. J Cardin
Approximating the coalescent with recombination
Phil Trans R Soc B, July 29, 2005; 360(1459): 1387 - 1393.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
L. Zhu and C. D. Bustamante
A Composite-Likelihood Approach for Detecting Directional Selection From DNA Sequence Data
Genetics, July 1, 2005; 170(3): 1411 - 1421.
[Abstract] [Full Text] [PDF]


Home page
Mol Biol EvolHome page
A. D. Tsaousis, D. P. Martin, E. D. Ladoukakis, D. Posada, and E. Zouros
Widespread Recombination in Published Animal mtDNA Sequences
Mol. Biol. Evol., April 1, 2005; 22(4): 925 - 933.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
S. Zollner and J. K. Pritchard
Coalescent-Based Association Mapping and Fine Mapping of Complex Trait Loci
Genetics, February 1, 2005; 169(2): 1071 - 1092.
[Abstract] [Full Text] [PDF]


Home page
Mol Biol EvolHome page
K. Ding, K. Zhou, J. Zhang, J. Knight, X. Zhang, and Y. Shen
The Effect of Haplotype-Block Definitions on Inference of Haplotype-Block Structure and htSNPs Selection
Mol. Biol. Evol., January 1, 2005; 22(1): 148 - 159.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
N. Nikoh, A. Duty, and G. Gibson
Effects of Population Structure and Sex on Association Between Serotonin Receptors and Drosophila Heart Rate
Genetics, December 1, 2004; 168(4): 1963 - 1974.
[Abstract] [Full Text] [PDF]


Home page
Hum Mol GenetHome page
J. Zhang, F. Li, J. Li, M. Q. Zhang, and X. Zhang
Evidence and characteristics of putative human {alpha} recombination hotspots
Hum. Mol. Genet., November 15, 2004; 13(22): 2823 - 2828.
[Abstract] [Full Text] [PDF]


Home page
DiabetesHome page
F. Gibson and P. Froguel
Genetics of the APM1 Locus and Its Contribution to Type 2 Diabetes Susceptibility in French Caucasians
Diabetes, November 1, 2004; 53(11): 2977 - 2983.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
P. Fearnhead, R. M. Harding, J. A. Schneider, S. Myers, and P. Donnelly
Application of Coalescent Methods to Reveal Fine-Scale Rate Variation and Recombination Hotspots
Genetics, August 1, 2004; 167(4): 2067 - 2081.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
S. Lin, A. Chakravarti, and D. J. Cutler
Haplotype and Missing Data Inference in Nuclear Families
Genome Res., August 1, 2004; 14(8): 1624 - 1632.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
J. D. Wall
Estimating Recombination Rates Using Three-Site Likelihoods
Genetics, July 1, 2004; 167(3): 1461 - 1473.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
S. E. Ptak, K. Voelpel, and M. Przeworski
Insights Into Recombination From Patterns of Linkage Disequilibrium in Humans
Genetics, May 1, 2004; 167(1): 387 - 397.
[Abstract] [Full Text] [PDF]


Home page
ScienceHome page
G. A. T. McVean, S. R. Myers, S. Hunt, P. Deloukas, D. R. Bentley, and P. Donnelly
The Fine-Scale Structure of Recombination Rate Variation in the Human Genome
Science, April 23, 2004; 304(5670): 581 - 584.
[Abstract] [Full Text] [PDF]