Inferring population history from molecular phylogenies

Sean Nee, Eddie C Holmes, Andrew Rambaut, Paul H Harvey

1995-07-29

Abstract:Variable molecular sequences sampled from a population can be used to infer its dynamic history. Graphical methods are developed and applied to real data, illustrating ways of navigating through hypothesis space with two landmarks for reference: constant population size and exponentially growing population size.

What problem does this paper attempt to address?

GENIE: estimating demographic history from molecular phylogenies

Oliver G Pybus, Andrew Rambaut

2002-10-01

Abstract:Summary: GENIE implements a statistical framework for inferring the demographic history of a population from phylogenies that have been reconstructed from sampled DNA sequences. The methods are based on population genetic models known collectively as coalescent theory. Availability: GENIE is available from http://evolve.zoo.ox.ac.uk. All popular operating systems are supported. Contact: oliver.pybus@zoo.ox.ac.uk
An integrated framework for the inference of viral population history from reconstructed genealogies

Oliver G Pybus, Andrew Rambaut, Paul H Harvey

2000-07-01

Abstract:We describe a unified set of methods for the inference of demographic history using genealogies reconstructed from gene sequence data. We introduce the skyline plot, a graphical, nonparametric estimate of demographic history. We discuss both maximum-likelihood parameter estimation and demographic hypothesis testing. Simulations are carried out to investigate the statistical properties of maximum-likelihood estimates of demographic parameters. The simulations reveal that (i) the performance of exponential growth model estimates is determined by a simple function of the true parameter values and (ii) under some conditions, estimates from reconstructed trees perform as well as estimates from perfect trees. We apply our methods to HIV-1 sequence data and find strong evidence that subtypes A and B have different demographic histories. We also provide the first (albeit tentative) genetic evidence for a …
Bayesian coalescent inference of past population dynamics from molecular sequences

Alexei J Drummond, Andrew Rambaut, Beth Shapiro, Oliver G Pybus

2005-02-09

Abstract:We introduce the Bayesian skyline plot, a new method for estimating past population dynamics through time from a sample of molecular sequences without dependence on a prespecified parametric model of demographic history. We describe a Markov chain Monte Carlo sampling procedure that efficiently samples a variant of the generalized skyline plot, given sequence data, and combines these plots to generate a posterior distribution of effective population size through time. We apply the Bayesian skyline plot to simulated data sets and show that it correctly reconstructs demographic history under canonical scenarios. Finally, we compare the Bayesian skyline plot model to previous coalescent approaches by analyzing two real data sets (hepatitis C virus in Egypt and mitochondrial DNA of Beringian bison) that have been previously investigated using alternative coalescent methods. In the bison analysis, we …
Jointly inferring the dynamics of population size and sampling intensity from molecular sequences

Kris V Parag, Louis du Plessis, Oliver G Pybus

2020-08-01

Abstract:Estimating past population dynamics from molecular sequences that have been sampled longitudinally through time is an important problem in infectious disease epidemiology, molecular ecology, and macroevolution. Popular solutions, such as the skyline and skygrid methods, infer past effective population sizes from the coalescent event times of phylogenies reconstructed from sampled sequences but assume that sequence sampling times are uninformative about population size changes. Recent work has started to question this assumption by exploring how sampling time information can aid coalescent inference. Here, we develop, investigate, and implement a new skyline method, termed the epoch sampling skyline plot (ESP), to jointly estimate the dynamics of population size and sampling rate through time. The ESP is inspired by real-world data collection practices and comprises a flexible model in …
Inferring the historical patterns of biological evolution

Mark Pagel

DOI: https://doi.org/10.1038/44766

IF: 64.8

1999-10-01

Nature

Abstract:Phylogenetic trees describe the pattern of descent amongst a group of species. With the rapid accumulation of DNA sequence data, more and more phylogenies are being constructed based upon sequence comparisons. The combination of these phylogenies with powerful new statistical approaches for the analysis of biological evolution is challenging widely held beliefs about the history and evolution of life on Earth.

multidisciplinary sciences
Extinction rates can be estimated from molecular phylogenies

Sean Nee, Edward C Holmes, Robert Mccredie May, Paul H Harvey

1994-04-29

Abstract:Molecular phylogenies can be used to reject null models of the way we think evolution occurred, including patterns of lineage extinction. They can also be used to provide maximum likelihood estimates of parameters associated with lineage birth and death rates. We illustrate: (i) how molecular phylogenies provide information about the extent to which particular clades are likely to be under threat from extinction; (ii) how cursory analyses of molecular phylogenies can lead to incorrect conclusions about the evolutionary processes that have been at work; and (iii) how different evolutionary processes leave distinctive marks on the structure of reconstructed phylogenies.
Exact Limits of Inference in Coalescent Models

James E. Johndrow,Julia A. Palacios

DOI: https://doi.org/10.48550/arXiv.1711.05724

2018-01-16

Abstract:Recovery of population size history from molecular sequence data is an important problem in population genetics. Inference commonly relies on a coalescent model linking the population size history to genealogies. The high computational cost of estimating parameters from these models usually compels researchers to select a subset of the available data or to rely on non-sufficient summary statistics for statistical inference. We consider the problem of recovering the true population size history from two possible alternatives on the basis of coalescent time data. We give exact expressions for the probability of selecting the correct alternative in a variety of biologically interesting cases as a function of the separation between the alternative size histories, the number of individuals, loci, and the sampling times. The results are applied to human population history. This work has significant implications for optimal design when the inferential goal is to test scientific hypotheses about population size trajectories in coalescent models with and without recombination.

Statistics Theory
Phylogeny Estimation by Integration over Isolation with Migration Models

Jody Hey,Yujin Chung,Arun Sethuraman,Joseph Lachance,Sarah Tishkoff,Vitor C Sousa,Yong Wang

DOI: https://doi.org/10.1093/molbev/msy162

IF: 10.7

2018-08-20

Molecular Biology and Evolution

Abstract:Phylogeny estimation is difficult for closely related populations and species, especially if they have been exchanging genes. We present a hierarchical Bayesian, Markov-chain Monte Carlo method with a state space that includes all possible phylogenies in a full Isolation-with-Migration model framework. The method is based on a new type of genealogy augmentation called a "hidden genealogy" that enables efficient updating of the phylogeny. This is the first likelihood-based method to fully incorporate directional gene flow and genetic drift for estimation of a species or population phylogeny. Application to human hunter-gatherer populations from Africa revealed a clear phylogenetic history, with strong support for gene exchange with an unsampled ghost population, and relatively ancient divergence between a ghost population and modern human populations, consistent with human/archaic divergence. In contrast, a study of five chimpanzee populations reveals a clear phylogeny with several pairs of populations having exchanged DNA, but does not support a history with an unsampled ghost population.

genetics & heredity,biochemistry & molecular biology,evolutionary biology
End-Epi: an application for inferring phylogenetic and population dynamical processes from molecular sequences

Andrew Rambaut, Paul H Harvey, Sean Nee

1997-06-01

Abstract:Motivation. Phylogenetic trees constructed from molecular sequences contain information about the evolutionary or population dynamical processes that created them. Here we describe a computer package (End-Epi) that uses graphical methods to allow researchers to make inferences about these processes from their data. Statistical analyses can be performed to test the consistency of the data with various competing hypotheses. Availability. End-Epi can be obtained by WWW from http://evolve.zoo.ox.ac.uk/ and by anonymous FTP from ftp://evolve.zoo.ox.ac.uk/packages/End-Epi10.hqx. This file contains the compiled application, the manual and a test tree.
Inferring heterogeneous evolutionary processes through time: from sequence substitution to phylogeography

Filip Bielejec, Philippe Lemey, Guy Baele, Andrew Rambaut, Marc A Suchard

2014-07-01

Abstract:Molecular phylogenetic and phylogeographic reconstructions generally assume time-homogeneous substitution processes. Motivated by computational convenience, this assumption sacrifices biological realism and offers little opportunity to uncover the temporal dynamics in evolutionary histories. Here, we propose an evolutionary approach that explicitly relaxes the time-homogeneity assumption by allowing the specification of different infinitesimal substitution rate matrices across different time intervals, called epochs, along the evolutionary history. We focus on an epoch model implementation in a Bayesian inference framework that offers great modeling flexibility in drawing inference about any discrete data type characterized as a continuous-time Markov chain, including phylogeographic traits. To alleviate the computational burden that the additional temporal heterogeneity imposes, we adopt a massively …
Exploring the demographic history of DNA sequences using the generalized skyline plot

Korbinian Strimmer, Oliver G Pybus

2001-12-01

Abstract:We present an intuitive visual framework, the generalized skyline plot, to explore the demographic history of sampled DNA sequences. This approach is based on a genealogy inferred from the sequences and provides a nonparametric estimate of effective population size through time. In contrast to previous related procedures, the generalized skyline plot is more applicable to cases where the underlying tree is not fully resolved and the data is not highly variable. This is achieved by the grouping of adjacent coalescent intervals. We employ a small-sample Akaike information criterion to objectively choose the optimal grouping strategy. We investigate the performance of our approach using simulation and subsequently apply it to HIV-1 sequences from central Africa and mtDNA sequences from red pandas.
Inferring demographic history from a spectrum of shared haplotype lengths

Kelley Harris,Rasmus Nielsen

DOI: https://doi.org/10.1371/journal.pgen.1003521

Abstract:There has been much recent excitement about the use of genetics to elucidate ancestral history and demography. Whole genome data from humans and other species are revealing complex stories of divergence and admixture that were left undiscovered by previous smaller data sets. A central challenge is to estimate the timing of past admixture and divergence events, for example the time at which Neanderthals exchanged genetic material with humans and the time at which modern humans left Africa. Here, we present a method for using sequence data to jointly estimate the timing and magnitude of past admixture events, along with population divergence times and changes in effective population size. We infer demography from a collection of pairwise sequence alignments by summarizing their length distribution of tracts of identity by state (IBS) and maximizing an analytic composite likelihood derived from a Markovian coalescent approximation. Recent gene flow between populations leaves behind long tracts of identity by descent (IBD), and these tracts give our method power by influencing the distribution of shared IBS tracts. In simulated data, we accurately infer the timing and strength of admixture events, population size changes, and divergence times over a variety of ancient and recent time scales. Using the same technique, we analyze deeply sequenced trio parents from the 1000 Genomes project. The data show evidence of extensive gene flow between Africa and Europe after the time of divergence as well as substructure and gene flow among ancestral hominids. In particular, we infer that recent African-European gene flow and ancient ghost admixture into Europe are both necessary to explain the spectrum of IBS sharing in the trios, rejecting simpler models that contain less population structure.
Molecular Phylogenetic Analyses and Real Life Data

Kerstin Hoef-Emden

DOI: https://doi.org/10.48550/arXiv.q-bio/0501030

2005-01-30

Abstract:In molecular phylogeny, relationships among organisms are reconstructed using DNA or protein sequences and are displayed as trees. A linear increase in the number of sequences results in an exponential increase of possible trees. Thus, inferring trees from molecular data was shown to be NP-hard. This causes problems, if large data sets are used. This review gives an introduction to molecular phylogenetic methods and to the problems biologists are facing in molecular phylogenetic analyses.

Genomics
Inference of viral evolutionary rates from molecular sequences

Alexei Drummond, Oliver G Pybus, Andrew Rambaut

2003-01-01

Abstract:1. Introduction 1.1. Ancestral Diversity and Evolutionary Non-independence. 2. General Linear Regression and Other Distance-Based Methods. 2.1. Root-to-tip Linear Regression
How Many Subpopulations is Too Many? Exponential Lower Bounds for Inferring Population Histories

Younhun Kim,Frederic Koehler,Ankur Moitra,Elchanan Mossel,Govind Ramnarayan

DOI: https://doi.org/10.48550/arXiv.1811.03177

2019-05-08

Abstract:Reconstruction of population histories is a central problem in population genetics. Existing coalescent-based methods, like the seminal work of Li and Durbin (Nature, 2011), attempt to solve this problem using sequence data but have no rigorous guarantees. Determining the amount of data needed to correctly reconstruct population histories is a major challenge. Using a variety of tools from information theory, the theory of extremal polynomials, and approximation theory, we prove new sharp information-theoretic lower bounds on the problem of reconstructing population structure -- the history of multiple subpopulations that merge, split and change sizes over time. Our lower bounds are exponential in the number of subpopulations, even when reconstructing recent histories. We demonstrate the sharpness of our lower bounds by providing algorithms for distinguishing and learning population histories with matching dependence on the number of subpopulations. Along the way and of independent interest, we essentially determine the optimal number of samples needed to learn an exponential mixture distribution information-theoretically, proving the upper bound by analyzing natural (and efficient) algorithms for this problem.

Populations and Evolution,Statistics Theory,Quantitative Methods
Improving Bayesian population dynamics inference: a coalescent-based model for multiple loci

Mandev S Gill, Philippe Lemey, Nuno R Faria, Andrew Rambaut, Beth Shapiro, Marc A Suchard

2013-03-01

Abstract: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 …
A Bayesian phylogenetic method to estimate unknown sequence ages

Beth Shapiro, Simon YW Ho, Alexei J Drummond, Marc A Suchard, Oliver G Pybus, Andrew Rambaut

2011-02-01

Abstract:Heterochronous data sets comprise molecular sequences sampled at different points in time. If the temporal range of the sampled sequences is large relative to the rate of mutation, the sampling times can directly calibrate evolutionary rates to calendar time. Here, we extend this calibration process to provide a full probabilistic method that utilizes temporal information in heterochronous data sets to estimate sampling times (leaf-ages) for sequenced for which this information unavailable. Our method is similar to relaxing the constraints of the molecular clock on specific lineages within a phylogenetic tree. Using a combination of synthetic and empirical data sets, we demonstrate that the method estimates leaf-ages reliably and accurately. Potential applications of our approach include incorporating samples of uncertain or radiocarbon-infinite age into ancient DNA analyses, evaluating the temporal signal in a …
Using Genetic Data to Build Intuition about Population History

Alan R. Rogers

DOI: https://doi.org/10.48550/arXiv.2201.02668

2022-01-08

Abstract:Genetic data are now routinely used to study the history of population size, subdivision, and gene flow. A variety of formal statistical methods is available for testing hypotheses and fitting models to data. Yet it is often unclear which hypotheses are worth testing, which models worth fitting. There is a need for less formal methods that can be used in exploratory analysis of genetic data. One approach to this problem uses *nucleotide site patterns*, which provide a simple summary of the pattern in genetic data. This article shows how to use them in exploratory data analysis.

Populations and Evolution
Ancestral inference from haplotypes and mutations

Robert C. Griffiths,Simon Tavaré

DOI: https://doi.org/10.48550/arXiv.1705.09485

2018-03-01

Abstract:We consider inference about the history of a sample of DNA sequences, conditional upon the haplotype counts and the number of segregating sites observed at the present time. After deriving some theoretical results in the coalescent setting, we implement rejection sampling and importance sampling schemes to perform the inference. The importance sampling scheme addresses an extension of the Ewens Sampling Formula for a configuration of haplotypes and the number of segregating sites in the sample. The implementations include both constant and variable population size models. The methods are illustrated by two human Y chromosome data sets.

Statistics Theory,Populations and Evolution
New inferences from tree shape: numbers of missing taxa and population growth rates

OG Pybus, A Rambaut, EC Holmes, PH Harvey

2002-12-01

Abstract:The relative positions of branching events in a phylogeny contain information about evolutionary and population dynamic processes. We provide new summary statistics of branching event times and describe how these statistics can be used to infer rates of species diversification from interspecies trees or rates of population growth from intraspecies trees. We also introduce a phylogenetic method for estimating the level of taxon sampling in a clade. Different evolutionary models and different sampling regimes can produce similar patterns of branching events, so it is important to consider explicitly the model assumptions involved when making evolutionary inferences. Results of an analysis of the phylogeny of the mosquito-borne flaviviruses suggest that there could be several thousand currently unidentified viruses in this clade.

Inferring population history from molecular phylogenies

GENIE: estimating demographic history from molecular phylogenies

An integrated framework for the inference of viral population history from reconstructed genealogies

Bayesian coalescent inference of past population dynamics from molecular sequences

Jointly inferring the dynamics of population size and sampling intensity from molecular sequences

Inferring the historical patterns of biological evolution

Extinction rates can be estimated from molecular phylogenies

Exact Limits of Inference in Coalescent Models

Phylogeny Estimation by Integration over Isolation with Migration Models

End-Epi: an application for inferring phylogenetic and population dynamical processes from molecular sequences

Inferring heterogeneous evolutionary processes through time: from sequence substitution to phylogeography

Exploring the demographic history of DNA sequences using the generalized skyline plot

Inferring demographic history from a spectrum of shared haplotype lengths

Molecular Phylogenetic Analyses and Real Life Data

Inference of viral evolutionary rates from molecular sequences

How Many Subpopulations is Too Many? Exponential Lower Bounds for Inferring Population Histories

Improving Bayesian population dynamics inference: a coalescent-based model for multiple loci

A Bayesian phylogenetic method to estimate unknown sequence ages

Using Genetic Data to Build Intuition about Population History

Ancestral inference from haplotypes and mutations

New inferences from tree shape: numbers of missing taxa and population growth rates