Performing inference on contemporary samples of DNA sequence data is an important and challenging task. Computationally intensive methods such as importance sampling (IS) are attractive because they make full use of the available data, but in the presence of recombination the large state space of genealogies can be prohibitive. In this article, we make progress by developing an efficient IS proposal distribution for a two-locus model of sequence data. We show that the proposal developed here leads to much greater efficiency, outperforming existing IS methods that could be adapted to this model. Among several possible applications, the algorithm can be used to find maximum likelihood estimates for mutation and crossover rates, and to perform ancestral inference. We illustrate the method on previously reported sequence data covering two loci either side of the well-studied TAP2 recombination hotspot. The two loci are themselves largely non-recombining, so we obtain a gene tree at each locus and are able to infer in detail the effect of the hotspot on their joint ancestry. We summarize this joint ancestry by introducing the gene graph, a summary of the well-known ancestral recombination graph.
Get full access to this article
View all access options for this article.
References
1.
BahloM., GriffithsR.2000. Inference from gene trees in a subdivided population. Theor. Popul. Biol., 57:79–95.
2.
BeaumontM.A., ZhangW., BaldingD.J.2002. Approximate Bayesian computation in population genetics. Genetics, 162:2025–2035.
3.
ChenY., XieJ., LiuJ.S.2005. Stopping-time resampling for sequential Monte Carlo methods. J. R. Stat. Soc. Ser. B Stat. Methodol., 67:199–217.
4.
CrawfordD.C., BhangaleT., LiN.et al.2004. Evidence for substantial fine-scale variation in recombination rates across the human genome. Nat. Genet., 36:700–706.
5.
DalyM.J., RiouxJ.D., SchaffnerS.F.et al.2001. High-resolution haplotype structure in the human genome. Nat. Genet., 29:229–232.
6.
De IorioM., GriffithsR.C.2004a. Importance sampling on coalescent histories I. Adv. Appl. Prob., 36:417–433.
7.
De IorioM., GriffithsR.C.2004b. Importance sampling on coalescent histories II. Adv. Appl. Prob., 36:434–454.
8.
EthierS.N., GriffithsR.C.1990. On the two-locus sampling distribution. J. Math. Biol., 29:131–159.
9.
FearnheadP.2003. Consistency of estimators of the population-scaled recombination rate. Theor. Popul. Biol., 64:67–79.
10.
FearnheadP., DonnellyP.2001. Estimating recombination rates from population genetic data. Genetics, 159:1299–1318.
11.
FrisseL., HudsonR.R., BartoszewiczA.et al.2001. Gene conversion and different population histories may explain the contrast between polymorphism and linkage disequilibrium levels. Am. J. Hum. Genet., 69:831–843.
12.
GabrielS.B., SchaffnerS.F., NguyenH.et al.2002. The structure of haplotype blocks in the human genome. Science, 296:2225–2229.
GriffithsR.C.1989. Genealogical-tree probabilities in the infinitely-many-site model. J. Math. Biol., 27:667–680.
15.
GriffithsR.C., MarjoramP.1996. Ancestral inference from samples of DNA sequences with recombination. J. Comput. Biol., 3:479–502.
16.
GriffithsR.C., MarjoramP.1997. An ancestral recombination graph, 257–270. DonnellyP., TavaréS.Progress in Population Genetics and Human Evolution, 87. Springer-Verlag: Berlin.
17.
GriffithsR.C., TavaréS.1994. Simulating probability distributions in the coalescent. Theor. Popul. Biol., 46:131–159.
18.
GriffithsR.C., JenkinsP.A., SongY.S.2008. Importance sampling and the two-locus model with subdivided population structure. Adv. Appl. Prob., 40:473–500.
19.
HardingR.M., FullertonS.M., GriffithsR.C.et al.1997. Archaic African and Asian lineages in the genetic ancestry of modern humans. Am. J. Hum. Genet., 60:772–789.
20.
HellenthalG., AutonA., FalushD.2008. Inferring human colonization using a copying model. PLoS Genet., 4:e1000078.
21.
HindsD. A., StuveL. L., NilsenG. B.et al.2005. Whole-genome patterns of common DNA variation in three human populations. Science, 307:1072–1079.
22.
HudsonR.R.1983. Properties of a neutral allele model with intragenic recombination. Theor. Popul. Biol., 23:183–201.
23.
HudsonR.R.2001. Two-locus sampling distributions and their application. Genetics, 159:1805–1817.
24.
HudsonR.R.2002. Generating samples under a Wright-Fisher neutral model of genetic variation. Bioinformatics, 18:337–338.
25.
HudsonR.R., KaplanN.L.1985. Statistical properties of the number of recombination events in the history of a sample of DNA sequences. Genetics, 111:147–164.
26.
InnanH., ZhangK., MarjoramP.et al.2005. Statistical tests of the coalescent model based on the haplotype frequency distribution and the number of segregating sites. Genetics, 169:1763–1777.
27.
JeffreysA.J., RitchieA., NeumannR.2000. High-resolution analysis of haplotype diversity and meiotic crossover in the human TAP2 recombination hotspot. Hum. Mol. Genet., 9:725–733.
28.
JenkinsP.A.2008. Importance sampling on the coalescent with recombination. [Ph.D. dissertation]University of Oxford: Oxford, UK.
29.
KingmanJ.F.C.1982. The coalescent. Stoch. Proc. Appl., 13:235–248.
30.
KuhnerM.K., YamatoJ., FelsensteinJ.2000. Maximum likelihood estimation of recombination rates from population data. Genetics, 156:1393–1401.
31.
LiN., StephensM.2003. Modeling linkage disequilibrium and identifying recombination hotspots using single-nucleotide polymorphism data. Genetics, 165:2213–2233.
32.
LiuJ.S.2001. Monte Carlo Strategies in Scientific Computing. Springer-Verlag: New York.
33.
MarjoramP., MolitorJ., PlagnolV.et al.2003. Markov chain Monte Carlo without likelihoods. Proc. Natl. Acad. Sci. USA, 100:15324–15328.
34.
McVeanG., AwadallaP., FearnheadP.2002. A coalescent-based method for detecting and estimating recombination from gene sequences. Genetics, 160:1231–1241.
35.
McVeanG., SpencerC.C.A., ChaixR.2005. Perspectives on human genetic variation from the HapMap Project. PLoS Genet., 1:e54.
36.
McVeanG.A.T., MyersS.R., HuntS.et al.2004. The fine-scale structure of recombination rate variation in the human genome. Science, 304:581–584.
37.
MyersS., BottoloL., FreemanC.et al.2005. A fine-scale map of recombination rates and hotspots across the human genome. Science, 310:321–324.
38.
NachmanM.W., CrowellS.L.2000. Estimate of the mutation rate per nucleotide in humans. Genetics, 156:297–304.
39.
NielsenR.2000. Estimation of population parameters and recombination rates from single nucleotide polymorphisms. Genetics, 154:931–942.
40.
PtakS.E., VoelpelK., PrzeworskiM.2004. Insight into recombination from patterns of linkage disequilibrium in humans. Genetics, 167:387–397.
41.
ReichD.E., CargillM., BolkS.et al.2001. Linkage disequilibrium in the human genome. Nature, 411:199–204.
42.
ReichD.E., SchaffnerS.F., DalyM.J.et al.2002. Human genome sequence variation and the influence of gene history, mutation and recombination. Nat. Genet., 32:135–142.
43.
ScheetP., StephensM.2006. A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am. J. Hum. Genet., 78:629–644.
44.
SimonsenK.L., ChurchillG.A.1997. A Markov chain model of coalescence with recombination. Theor. Popul. Biol., 52:43–59.
45.
SousaV.C., FritzM., BeaumontM.A.et al.2009. Approximate Bayesian Computation (ABC) without summary statistics: the case of admixture. Genetics, 181:1507–1519.
46.
StephensM.2001. Inference under the coalescent. BaldingD., BishopM., CanningsC.Handbook of Statistical Genetics. Wiley: Chichester, UK.
47.
StephensM., DonnellyP.2000. Inference in molecular population genetics. J. R. Stat. Soc. Ser. B Stat. Methodol., 62:605–655.
48.
StephensM., ScheetP.2005. Accounting for decay of linkage disequilibrium in haplotype inference and missing-data imputation. Am. J. Hum. Genet., 76:449–462.
49.
TajimaF.1989. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics, 123:585–595.
50.
The International HapMap Consortium. 2007. A second generation human haplotype map of over 3.1 million SNPs. Nature, 449:851–861.
WangY., RannalaB.2008. Bayesian inference of fine-scale recombination rates using population genomic data. Philos. Trans. R. Soc. Lond. B Biol. Sci., 363:3921–3930.
53.
WiufC., HeinJ.1999. The ancestry of a sample of sequences subject to recombination. Genetics, 151:1217–1228.
