Approximating the coalescent with recombination

Gil McVean, Niall J. Cardin

Philosophical Transactions of the Royal Society B Biological Sciences · 2005 · 429 citations · 54 references

DOIFull text

Open access

Concepts

TL;DR

The coalescent with recombination models genealogical histories and genetic variation, yet its likelihood is computationally prohibitive for inference. We examine why this model is difficult and investigate whether simpler, more tractable models can support genealogy‑based inference. We propose a simplified coalescent that forbids coalescence between lineages lacking overlapping ancestral material. This simplified process retains a Markovian structure and reproduces the full model’s patterns of genetic variation and inference performance.

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.

References

54