The matrix coalescent and an application to human single-nucleotide polymorphisms.

Area of Science:

• Population genetics
• Computational biology
• Evolutionary dynamics

Background:

• The coalescent process is fundamental to population genetics for inferring evolutionary history.
• Traditional coalescent models can be complex to implement, especially with changing population sizes.
• Understanding past population dynamics is crucial for interpreting genetic variation.

Purpose of the Study:

• To introduce a matrix coalescent model as a reformulation of the standard coalescent process.
• To demonstrate the utility of this model in calculating the site frequency spectrum.
• To apply the model to genetic data for testing hypotheses about human population size history.

Main Methods:

• Reformulated the coalescent process as a Markov process using a matrix formulation.
• Focused on the decay of ancestral lineages backward in time, ignoring gene tree topology.
• Calculated the expectation of the site frequency spectrum under different population history scenarios.

Main Results:

• The matrix formulation simplifies calculations for populations with changing sizes by isolating effects in eigenvalues.
• The model successfully calculates the expected site frequency spectrum.
• Applied to single nucleotide polymorphism (SNP) data to infer human population size changes.

Conclusions:

• The matrix coalescent offers a computationally convenient approach for population genetics inference.
• This method provides a robust framework for analyzing genetic data to reconstruct population histories.
• The study demonstrates the practical application of the matrix coalescent in human evolutionary studies.