Asymptotic Distributions of Coalescence Times and Ancestral Lineage Numbers for Populations with Temporally Varying Size
Open Access
- 1 July 2013
- journal article
- Published by Oxford University Press (OUP) in Genetics
- Vol. 194 (3), 721-736
- https://doi.org/10.1534/genetics.113.151522
Abstract
The distributions of coalescence times and ancestral lineage numbers play an essential role in coalescent modeling and ancestral inference. Both exact distributions of coalescence times and ancestral lineage numbers are expressed as the sum of alternating series, and the terms in the series become numerically intractable for large samples. More computationally attractive are their asymptotic distributions, which were derived in Griffiths (1984) for populations with constant size. In this article, we derive the asymptotic distributions of coalescence times and ancestral lineage numbers for populations with temporally varying size. For a sample of size n, denote by Tm the mth coalescent time, when m + 1 lineages coalesce into m lineages, and An(t) the number of ancestral lineages at time t back from the current generation. Similar to the results in Griffiths (1984), the number of ancestral lineages, An(t), and the coalescence times, Tm, are asymptotically normal, with the mean and variance of these distributions depending on the population size function, N(t). At the very early stage of the coalescent, when t → 0, the number of coalesced lineages n − An(t) follows a Poisson distribution, and as m → n, n(n−1)Tm/2N(0) follows a gamma distribution. We demonstrate the accuracy of the asymptotic approximations by comparing to both exact distributions and coalescent simulations. Several applications of the theoretical results are also shown: deriving statistics related to the properties of gene genealogies, such as the time to the most recent common ancestor (TMRCA) and the total branch length (TBL) of the genealogy, and deriving the allele frequency spectrum for large genealogies. With the advent of genomic-level sequencing data for large samples, the asymptotic distributions are expected to have wide applications in theoretical and methodological development for population genetic inference.Keywords
This publication has 38 references indexed in Scilit:
- Intercoalescence Time Distribution of Incomplete Gene Genealogies in Temporally Varying Populations, and Applications in Population Genetic InferenceAnnals of Human Genetics, 2013
- A Simple Method for Finding Explicit Analytic Transition Densities of Diffusion Processes with General Diploid SelectionGenetics, 2012
- Non-equilibrium allele frequency spectra via spectral methodsTheoretical Population Biology, 2011
- Deep resequencing reveals excess rare recent variants consistent with explosive population growthNature Communications, 2010
- A map of human genome variation from population-scale sequencingNature, 2010
- Inferring the Joint Demographic History of Multiple Populations from Multidimensional SNP Frequency DataPLoS Genetics, 2009
- The Joint Allele-Frequency Spectrum in Closely Related SpeciesGenetics, 2007
- ESTIMATION OF THE NUMBER OF INDIVIDUALS FOUNDING COLONIZED POPULATIONSEvolution, 2007
- The age of a mutation in a general coalescent treeCommunications in Statistics. Stochastic Models, 1998
- Sampling theory for neutral alleles in a varying environmentPhilosophical Transactions B, 1994