Cookies on this website

We use cookies to ensure that we give you the best experience on our website. If you click 'Accept all cookies' we'll assume that you are happy to receive all cookies and you won't see this message again. If you click 'Reject all non-essential cookies' only necessary cookies providing core functionality such as security, network management, and accessibility will be enabled. Click 'Find out more' for information on how to change your cookie settings.

Effective population size is fundamental in population genetics and characterizes genetic diversity. To infer past population dynamics from molecular sequence data, coalescent-based models have been developed for Bayesian nonparametric estimation of effective population size over time. Among the most successful is a Gaussian Markov random field (GMRF) model for a single gene locus. Here, we present a generalization of the GMRF model that allows for the analysis of multilocus sequence data. Using simulated data, we demonstrate the improved performance of our method to recover true population trajectories and the time to the most recent common ancestor (TMRCA). We analyze a multilocus alignment of HIV-1 CRF02_AG gene sequences sampled from Cameroon. Our results are consistent with HIV prevalence data and uncover some aspects of the population history that go undetected in Bayesian parametric estimation. Finally, we recover an older and more reconcilable TMRCA for a classic ancient DNA data set.

Original publication

DOI

10.1093/molbev/mss265

Type

Journal article

Journal

Mol Biol Evol

Publication Date

03/2013

Volume

30

Pages

713 - 724

Keywords

Algorithms, Bayes Theorem, Computer Simulation, Evolution, Molecular, Genes, Viral, Genetic Loci, Genetic Speciation, HIV-1, Humans, Markov Chains, Models, Genetic, Monte Carlo Method, Mutation, Population Density, Population Dynamics, Statistics, Nonparametric