Approximating the coalescent with recombination - PubMed (original) (raw)
Comparative Study
Approximating the coalescent with recombination
Gilean A T McVean et al. Philos Trans R Soc Lond B Biol Sci. 2005.
Abstract
The coalescent with recombination describes the distribution of genealogical histories and resulting patterns of genetic variation in samples of DNA sequences from natural populations. However, using the model as the basis for inference is currently severely restricted by the computational challenge of estimating the likelihood. We discuss why the coalescent with recombination is so challenging to work with and explore whether simpler models, under which inference is more tractable, may prove useful for genealogy-based inference. We introduce a simplification of the coalescent process in which coalescence between lineages with no overlapping ancestral material is banned. The resulting process has a simple Markovian structure when generating genealogies sequentially along a sequence, yet has very similar properties to the full model, both in terms of describing patterns of genetic variation and as the basis for statistical inference.
Figures
Figure 1
The ratio of the average number of recombination events in the ARG for the standard coalescent to the average number of recombination events in the SMC model for _n_=2. The average number of recombination events in the SMC is equal to ρ.
Figure 2
The sequentially Markov coalescent with recombination. The point of the recombination event (indicated by a crossmark) is placed uniformly on the tree. The branch above it is removed and the lineage coalesces back to the remaining tree at a rate proportional to the number of lineages present.
Figure 3
The decay of LD as a function of genetic distance (ρ) as approximated by _σ_d2 under the standard coalescent process (black) and the sequentially Markov version (grey).
Figure 4
Log likelihood surface for θ and ρ under the standard coalescent and sequentially Markov processes. Although the maximum likelihood estimates of both parameters are very similar under the two models, it should also be noted that the likelihood surfaces are very flat.
Figure 5
Expected values of the number of recombination events and TMRCA at each position under the standard coalescent (black) and sequentially Markov processes (grey). Maximum likelihood estimates of θ and ρ were used for each. Open circles represent the position of the mutations in the sequences.
Similar articles
- Topological linkage disequilibrium calculated from coalescent genealogies.
Wirtz J, Rauscher M, Wiehe T. Wirtz J, et al. Theor Popul Biol. 2018 Dec;124:41-50. doi: 10.1016/j.tpb.2018.09.001. Epub 2018 Sep 19. Theor Popul Biol. 2018. PMID: 30243857 - The Promise of Inferring the Past Using the Ancestral Recombination Graph.
Brandt DYC, Huber CD, Chiang CWK, Ortega-Del Vecchyo D. Brandt DYC, et al. Genome Biol Evol. 2024 Feb 1;16(2):evae005. doi: 10.1093/gbe/evae005. Genome Biol Evol. 2024. PMID: 38242694 Free PMC article. - Markovian approximation to the finite loci coalescent with recombination along multiple sequences.
Hobolth A, Jensen JL. Hobolth A, et al. Theor Popul Biol. 2014 Dec;98:48-58. doi: 10.1016/j.tpb.2014.01.002. Epub 2014 Jan 28. Theor Popul Biol. 2014. PMID: 24486389 - Linkage disequilibrium: what history has to tell us.
Nordborg M, Tavaré S. Nordborg M, et al. Trends Genet. 2002 Feb;18(2):83-90. doi: 10.1016/s0168-9525(02)02557-x. Trends Genet. 2002. PMID: 11818140 Review. - Ancestral Population Genomics.
Dutheil JY, Hobolth A. Dutheil JY, et al. Methods Mol Biol. 2019;1910:555-589. doi: 10.1007/978-1-4939-9074-0_18. Methods Mol Biol. 2019. PMID: 31278677 Review.
Cited by
- Exact Decoding of a Sequentially Markov Coalescent Model in Genetics.
Ki C, Terhorst J. Ki C, et al. J Am Stat Assoc. 2024;119(547):2242-2255. doi: 10.1080/01621459.2023.2252570. Epub 2023 Oct 3. J Am Stat Assoc. 2024. PMID: 39323740 - Improved inference of population histories by integrating genomic and epigenomic data.
Sellinger T, Johannes F, Tellier A. Sellinger T, et al. Elife. 2024 Sep 12;12:RP89470. doi: 10.7554/eLife.89470. Elife. 2024. PMID: 39264367 Free PMC article. - A general and efficient representation of ancestral recombination graphs.
Wong Y, Ignatieva A, Koskela J, Gorjanc G, Wohns AW, Kelleher J. Wong Y, et al. Genetics. 2024 Sep 4;228(1):iyae100. doi: 10.1093/genetics/iyae100. Genetics. 2024. PMID: 39013109 Free PMC article. - Biases in ARG-Based Inference of Historical Population Size in Populations Experiencing Selection.
Marsh JI, Johri P. Marsh JI, et al. Mol Biol Evol. 2024 Jul 3;41(7):msae118. doi: 10.1093/molbev/msae118. Mol Biol Evol. 2024. PMID: 38874402 Free PMC article. - TRAILS: Tree reconstruction of ancestry using incomplete lineage sorting.
Rivas-González I, Schierup MH, Wakeley J, Hobolth A. Rivas-González I, et al. PLoS Genet. 2024 Feb 8;20(2):e1010836. doi: 10.1371/journal.pgen.1010836. eCollection 2024 Feb. PLoS Genet. 2024. PMID: 38330138 Free PMC article.
References
- Crawford D.C, Bhangale T, Li N, Hellenthal G, Rieder M.J, Nickerson D.A, Stephens M. Evidence for substantial fine-scale variation in recombination rates across the human genome. Nat. Genet. 2004;36:700–706. - PubMed
- Fearnhead P, Donnelly P. Approximate likelihood methods for estimating local recombination rates. J. R. Stat. Soc. B. 2002;64:1–24.
Publication types
MeSH terms
LinkOut - more resources
Full Text Sources
Other Literature Sources