Insights into recombination from population genetic variation

Current Opinion in Genetics and Development

Volumn 16, Issue 6, 2006, Pages 565-572

  Hellenthal, Garrett ^a   Stephens, Matthew ^a  

^a UNIVERSITY OF WASHINGTON   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

CORRELATION ANALYSIS; DNA SEQUENCE; GENE CONVERSION; GENE SEQUENCE; GENETIC IDENTIFICATION; GENETIC RECOMBINATION; GENETIC VARIABILITY; HUMAN; INFORMATION PROCESSING; MOLECULAR EVOLUTION; NONHUMAN; POPULATION GENETICS; PRIORITY JOURNAL; REVIEW; SINGLE NUCLEOTIDE POLYMORPHISM; STATISTICAL ANALYSIS;

ANIMALS; CROSSING OVER, GENETIC; DATABASES, GENETIC; EVOLUTION, MOLECULAR; GENETICS, POPULATION; HUMANS; POPULATION; RECOMBINATION, GENETIC; VARIATION (GENETICS);

EID: 33751026170     PISSN: 0959437X     EISSN: None     Source Type: Journal    
DOI: 10.1016/j.gde.2006.10.001     Document Type: Review

References (50)

1. Pritchard J.K., and Przeworski M. Linkage disequilibrium in humans: models and data. Am J Hum Genet 69 (2001) 1-14
2. Chakravarti A., Buetow K.H., Antonarakis S.E., Waber P.G., Boehm C.D., and Kazazian H.H. Nonuniform recombination within the human β-globin gene cluster. Am J Hum Genet 36 (1984) 1239-1258
3. Hudson R.R. Estimating the recombination parameter of a finite population model without selection. Genet Res 50 (1987) 245-250
4. The International HapMap Consortium. A haplotype map of the human genome. Nature 437 (2005) 1299-1320
5. Hinds D.A., Stuve L.L., Nilsen G.B., Halperin E., Eskin E., Ballinger D.G., Frazer K.A., and Cox D.R. Whole-genome patterns of common DNA variation in three human populations. Science 307 (2005) 1072-1079
6. Hudson R.R. Two-locus sampling distributions and their application. Genetics 159 (2001) 1805-1817
7. Fearnhead P., and Donnelly P. Estimating recombination rates from population genetic data. Genetics 159 (2001) 1299-1318
8. Li N., and Stephens M. Modeling linkage disequilibrium and identifying recombination hotspots using single-nucleotide polymorphism data. Genetics 165 (2003) 2213-2233
9. Broman K.W., Murray J.C., Sheffield V.C., White R.L., and Weber J.L. Comprehensive human genetic map: individual and sex-specific variation in recombination. Am J Hum Genet 63 (1998) 861-869
10. Kong A., Gudbjartsson D.F., Sainz J., Jonsdottir G.M., Gudjonsson S.A., Richardsson B., Sigurdardottir S., Barnard J., Hallbeck B., Masson G., et al. A high-resolution recombination map of the human genome. Nat Genet 31 (2002) 241-247
11. Kingman J.F.C. The coalescent. Stochastic Process Appl 13 (1982) 235-248
12. Hudson R.R. Properties of a neutral allele model with intragenic recombination. Theor Popul Biol 23 (1983) 183-201
13. Griffiths R.C., and Marjoram P. Ancestral inference from samples of DNA sequences with recombination. J Comput Biol 3 (1996) 479-502
14. Kuhner M.K., Yamato J., and Felsenstein J. Maximum likelihood estimation of recombination rates from population data. Genetics 156 (2000) 1393-1401
15. McVean G., Awadalla P., and Fearnhead P. A coalescent-based method for detecting and estimating recombination from gene sequences. Genetics 160 (2002) 1231-1241
16. Frisse L., Hudson R.R., Bartoszewicz A., Wall J.D., Donfack J., and Di Rienzo A. Gene conversion and different population histories may explain the contrast between polymorphism and linkage disequilibrium levels. Am J Hum Genet 69 (2001) 831-843
17. Ptak S.E., Voelpel K., and Przeworski M. Insights into recombination from patterns of linkage disequilibrium in humans. Genetics 167 (2004) 387-397. The authors adapt Hudson's [6] composite-likelihood approach to jointly estimate crossover and gene conversion in the SeattleSNPs dataset, assuming crossover rates to be constant within genes, and f to be constant across the genome. Using a profile-likelihood approach that simultaneously maximizes ρ per gene, and f across genes, they estimate genome-wide average f to be ∼1.0 in the African Americans (using 84 genes) and ∼0.3 in the Centre d'Etude du Polymorphisme Humain (CEPH) population (77 genes), assuming a tract length of 500 bp.
18. McVean G.A., Myers S.R., Hunt S., Deloukas P., Bentley D.R., and Donnelly P. The fine-scale structure of recombination rate variation in the human genome. Science 304 (2004) 581-584. The authors update Hudson's [6] composite-likelihood method to estimate variable crossover rates along a sequence - in this case of kilobase scale. They examine rates in several different datasets, including sampled chromosomes from European and African American individuals (average marker spacing 2.3 kb), and sampled chromosomes of CEPH and Yoruban descent (average marker spacing 5.8 kb). Also, finding hotspots to be narrow, they estimate that about 50% of all crossover events take place in <10% of the sequence. They find crossover rates to be diminished in genic regions and note that although G + C content and kilobase-scale crossover rates appear to be correlated, this correlation disappears when looking at finer scales such as hotspots.
19. Fearnhead P., Harding R.M., Schneider J.A., Myers S., and Donnelly P. Application of coalescent methods to reveal fine-scale rate variation and recombination hotspots. Genetics 167 (2004) 2067-2081
20. Fearnhead P., and Smith N.G. A novel method with improved power to detect recombination hotspots from polymorphism data reveals multiple hotspots in human genes. Am J Hum Genet 77 (2005) 781-794. The authors extend the importance-sampling, simulation approach of Fearnhead et al. [19] to locate probable crossover hotspots using LD data, considering overlapping sub-regions (assuming them independent) of approximately five SNPs each. Analyzing 89 genes from SeattleSNPs in African Americans and Europeans (both jointly and separately), they find hotspot frequencies of about 1 per 30-40 kb, as well as some evidence that several hotspots (perhaps as many as 30%) do not overlap between the two populations. They also find evidence for increased SNP density and G + C content in hotspots relative to outside them.
21. Smith N.G., and Fearnhead P. A comparison of three estimators of the population-scaled recombination rate: accuracy and robustness. Genetics 171 (2005) 2051-2062. The authors use coalescent-based simulations with 50 haplotypes, θ = 1.0/kb, and sequence lengths of 2, 10, 25 and 100 kb, to test the crossover estimation methods of Hudson [6], Fearnhead and Donnelly [7], and Li and Stephens [8]. They find that all methods behave similarly, and that crossover estimates in all three are biased upwards in the presence of gene conversion (simulated in rates of 1, 5 or 10/kb). Some extreme SNP ascertainment schemes, when left unaccounted for, can also be a problem.
22. Crawford D.C., Bhangale T., Li N., Hellenthal G., Rieder M.J., Nickerson D.A., and Stephens M. Evidence for substantial fine-scale variation in recombination rates across the human genome. Nat Genet 36 (2004) 700-706. The authors find evidence for crossover hotspots in 35 of 74 genes from the SeattleSNPs dataset. Crossover rates on a kilobase scale appear to be correlated with sequence features such as G + C content and nucleotide diversity, although this correlation does not appear to extend to hotspots.
23. Jeffreys A.J., Ritchie A., and Neumann R. High resolution analysis of haplotype diversity and meiotic crossover in the human TAP2 recombination hotspot. Hum Mol Genet 9 (2000) 725-733
24. Stumpf M.P.H., and McVean G.A.T. Estimating recombination rates from population genetic data. Nat Rev Genet 4 (2003) 959-968
25. Gabriel S.B., Schaffner S.F., Nguyen H., Moore J.M., Roy J., Blumenstiel B., Higgins J., DeFelice M., Lochner A., Faggart M., et al. The structure of haplotype blocks in the human genome. Science 296 (2002) 2225-2229
26. Jeffreys A.J., Kauppi L., and Neumann R. Intensely punctate meiotic recombination in the class II region of the major histocompatibility complex. Nat Genet 29 (2001) 217-222
27. ••], whereas the authors' population-based crossover rates (on Mb scales) and pedigree-based estimates [10] appear to be strongly correlated. This leads the authors to speculate that crossover rates are constrained over large (i.e. Mb) regions but are perhaps rapidly evolving at finer scales.
28. Jeffreys A.J., and Neumann R. Reciprocal crossover asymmetry and meiotic drive in a human recombination hot spot. Nat Genet 31 (2002) 267-271
29. Jeffreys A.J., and Neumann R. Factors influencing recombination frequency and distribution in a human meiotic crossover hotspot. Hum Mol Genet 14 (2005) 2277-2287. The authors examine crossover rates in the hotspot NID1 of chromosome 1 (near the minisatellite MS32 hotspot) in seven men (with a mean peak of crossover rate = 70 cM/Mb). They find evidence of crossover asymmetry (the 'crossover-suppressing' allele perhaps reduces crossover by ≥2.8-fold) and transmission distortion; at an allele near the center of the hotspot, the crossover-suppressing allele was over-transmitted by about 3:1. They also looked at rates of gene conversion for one man; rates were only a quarter of the observed crossovers for this man, and transmission distortion rates were similar to those of crossovers.
30. Myers S., Spencer C.C., Auton A., Bottolo L., Freeman C., Donnelly P., and McVean G. The distribution and causes of meiotic recombination in the human genome. Biochem Soc Trans 34 (2006) 526-530
31. Rattray A.J., Shafer B.K., McGill C.B., and Strathern J.N. The roles of rev3 and rad57 in double-strand-break-repair-induced mutagenesis of Saccharomyces cerevisiae. Genetics 162 (2002) 1063-1077
32. Vander Molen J., Frisse L.M., Fullerton S.M., Qian Y., del Bosque-Plata L., Hudson R.R., and Di Rienzo A. Population genetics of capn10 and gpr35: implications for the evolution of type 2 diabetes variants. Am J Hum Genet 76 (2005) 548-560
33. Reed F.A., and Tishkoff S.A. Positive selection can create false hotspots of recombination. Genetics 172 (2006) 2011-2014
34. Wall J.D., Frisse L.A., Hudson R.R., and Di Rienzo A. Comparative linkage-disequilibrium analysis of the β-globin hotspot in primates. Am J Hum Genet 73 (2003) 1330-1340
35. Ptak S.E., Roeder A.D., Stephens M., Gilad Y., Paabo S., and Przeworski M. Absence of the TAP2 human recombination hotspot in chimpanzees. PLoS Biol 2 (2004) 849-855. The authors compare LD patterns between 30 UK humans and 24 resequenced western chimpanzees in the TAP2 region. In particular, they use the methods of Li and Stephens [8] to estimate ρ across the fine-scale region. They find substantial evidence for a hotspot in the human data, with crossover activity at least 10 times greater than rates outside the hotspot, whereas there is evidence against a hotspot in the chimpanzee data. They also highlight an example of positive association between rates of crossover and mutation, noting that diversity levels are significantly higher than average for the human data in the hotspot region, although this does not appear to be the case in chimpanzees.
36. Ptak S.E., Hinds D.A., Koehler K., Nickel B., Patil N., Ballinger D.G., Przeworski M., Frazer K.A., and Paabo S. Fine-scale recombination patterns differ between chimpanzees and humans. Nat Genet 37 (2005) 429-434. The authors compare LD patterns between the DNA of 23 humans and 8 central chimpanzees, in 14 Mb of sequence from two disjoint regions of chromosome 21 (and its chimp homologue), to determine what proportion of hotspots appear to be shared between the two species. Using the methods of Li and Stephens [8], they find that only 3 of the 39 hotspots inferred in the chimpanzees appeared in humans; lack of power to find chimpanzee hotpots makes it more difficult to infer how many characterized human hotspots appear in the chimps. Through simulation, they find that this observation is consistent with the hypothesis that hotspots are independently distributed along the genome of the two species. This suggests relatively rapid hotspot evolution. The authors also find that background crossover rates - rates outside of hotspots - for 50-kb windows are moderately correlated between the species, suggesting that perhaps recombination rates are conserved at large scales but not at fine scales.
37. Winckler W., Myers S.R., Richter D.J., Onofrio R.C., McDonald G.J., Bontrop R.E., McVean G.A., Gabriel S.B., Reich D., Donnelly P., and Altshuler D. Comparison of fine-scale recombination rates in humans and chimpanzees. Science 308 (2005) 107-111. The authors inferred rates of crossover from resequencing polymorphism data at orthologous loci in western chimpanzees and in CEPH and Yoruban HapMap ENCODE regions. They found statistical support for 18 hotspots in humans and 3 in chimps (ρ < 0.01), but with a striking lack of overlap: a crossover hotspot in one species was typically found to have crossover rate 10 to 60 times larger than that of the same region of the other species, which was typically not 'hot' at all relative to its background crossover rates.
38. e = 10000), which they compare against that found by sperm analysis. In three of the hotspots, the estimated values from LD analysis were significantly different (e.g. r for NID2a was estimated as 55 cM/Mb from LD analysis, compared with 10 cM/Mb from sperm analysis).
39. Tiemann-Boege I., Calabrese P., Cochran D.M., Sokol R., and Arnheim N. High-resolution recombination patterns in a region of human chromosome 21 measured by sperm typing. PLoS Genet 2 (2006) e70. The authors examine crossover activity across a 103 kb segment of chromosome 21, using the sperm of 240 individuals. They find two hotspots in this region, each of width of 1-2 kb, and estimate that 71% of the total crossover activity occurs in ∼12% of the sequence. They find that the patterns of estimated crossover activity from sperm analysis correspond well to those of various LD methods. Furthermore, studying one of the hotspots in more fine-scale detail in three men showed that the majority of the crossover activity for two of the men mapped to one ∼2 kb interval, whereas in the third man the majority of activity mapped to an apparently distinct region ∼2 kb away. The authors speculate that when one hotspot dies away, a new hotspot arises relatively nearby.
40. Jeffreys A.J., and May C.A. Intense and highly localized gene conversion activity in human meiotic crossover hot spots. Nat Genet 36 (2004) 151-156. The authors examine the DNA3 crossover hotspot of the human MHC by sperm analysis of two men, finding DNA3 to be a hotspot for gene conversion as well; they speculate that DNA3 is a hotspot for double-strand-breaks overall. They find that gene conversions occur 4-15 times more often than crossovers in the two men studied. They estimate tract length to be, on average, 55-290 bp per gene conversion event, and note that tracts for conversion accompanying crossover appear to be somewhat larger (∼460 bp).
41. Holloway K., Lawson V.E., and Jeffreys A.J. Allelic recombination and de novo deletions in sperm in the human β-globin gene region. Hum Mol Genet 17 (2006) 1099-1111. The authors use sperm analysis to characterize the crossover activity of 13.5 kb of sequence in the β-globin gene region. Using two donors, they find an intense, narrow crossover hotspot, ∼1.2 kb wide, with crossover activity of ∼150-260 cM/Mb. Examining the hotspot of one of these men for gene conversion alone, however, detected no gene conversions. This finding does not appear to be attributable to lack of power or marker density. The authors conclude that the relative rate of gene conversion to crossover in this region must be <1/12. Furthermore, they do not find any evidence of association between deletion breakpoints and high crossover activity.
42. Szostak J.W., Orr-Weaver T.L., Rothstein R.J., and Stahl F.W. The double-strand-break repair model for recombination. Cell 33 (1983) 25-35
43. Ardlie K., Liu-Cordero S.N., Eberle M.A., Daly M., Barrett J., Winchester E., Lander E.S., and Kruglyak L. Lower-than-expected linkage disequilibrium between tightly linked markers in humans suggests a role for gene conversion. Am J Hum Genet 69 (2001) 582-589
44. Padhukasahasram B., Marjoram P., and Nordborg M. Estimating the rate of gene conversion on human chromosome 21. Am J Hum Genet 75 (2004) 386-397. Using simple multilocus summary statistics, the authors estimate rates of gene conversion and crossover in humans. The rates are based on patterns of LD, using the DNA from a worldwide sample of 20 chromosomes representing a 28 Mb region of chromosome 21. They estimate f to be ∼1.6, assuming a tract length of 500 bp. Estimating rates over overlapping 2 Mb windows, they find rates of crossover and gene conversion to be uncorrelated.
45. a) regions of the Drosophila melanogaster X chromosome. Genetics 156 (2000) 1837-1852
46. Andolfatto P., and Wall J.D. Linkage disequilibrium patterns across a recombination gradient in African Drosophila melanogaster. Genetics 165 (2003) 1289-1305
47. Plagnol V., Padhukasahasram B., Wall J.D., Marjoram P., and Nordborg M. Relative influences of crossing over and gene conversion on the pattern of linkage disequilibrium in Arabidopsis thaliana. Genetics 172 (2006) 2441-2448. The authors examine rates of gene conversion in Arabidopsis thaliana, using a genome-wide survey consisting of 1347 fragments (500-600 bp each) sequenced in 96 accessions. To do so, they perform coalescent simulations with growth parameters that match the observed allele frequency spectrum, and a crossover rate estimated using the method of Hudson [6]. They then find gene conversion parameters that match the observed LD patterns in their data. They find, on average, that f = 1, or that the crossover rate and gene conversion rate are both ∼0.3/kb, assuming a tract length of 100 bp. They also estimate rates of crossover and gene conversion across 14 genome-wide 2-Mb windows - each window having, on average, 18 fragments - and find that variable rates of f, crossover, and gene conversion across these windows fit the data significantly better than either fixing f across regions or taking an identical crossover and gene conversion rate in all regions. Furthermore, they find estimated rates of crossover and gene conversion across 43 such 2-Mb windows to be uncorrelated.
48. Wall J.D. Estimating recombination rates using three-site likelihoods. Genetics 167 (2004) 1461-1473
49. Hudson R.R. Linkage disequilibrium and recombination. In: Balding D.J., Bishop M., and Cannings C. (Eds). Handbook of Statistical Genetics (2001), Wiley 309-323
50. Stephens M., Smith N.J., and Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet 68 (2001) 978-989