Abstract:Phylogenetic analyses which include fossils or molecular sequences that are sampled through time require models that allow one sample to be a direct ancestor of another sample. As previously available phylogenetic inference tools assume that all samples are tips, they do not allow for this possibility. We have developed and implemented a Bayesian Markov Chain Monte Carlo (MCMC) algorithm to infer what we call sampled ancestor trees, that is, trees in which sampled individuals can be direct ancestors of other sampled individuals. We use a family of birth-death models where individuals may remain in the tree process after sampling, in particular we extend the birth-death skyline model [Stadler et al., 2013] to sampled ancestor trees. This method allows the detection of sampled ancestors as well as estimation of the probability that an individual will be removed from the process when it is sampled. We show that even if sampled ancestors are not of specific interest in an analysis, failing to account for them leads to significant bias in parameter estimates. We also show that sampled ancestor birth-death models where every sample comes from a different time point are non-identifiable and thus require one parameter to be known in order to infer other parameters. We apply our phylogenetic inference accounting for sampled ancestors to epidemiological data, where the possibility of sampled ancestors enables us to identify individuals that infected other individuals after being sampled and to infer fundamental epidemiological parameters. We also apply the method to infer divergence times and diversification rates when fossils are included along with extant species samples, so that fossilisation events are modelled as a part of the tree branching process. Such modelling has many advantages as argued in the literature. The sampler is available as an open-source BEAST2 package (https://github.com/CompEvol/sampled-ancestors).

Bayesian coalescent inference of past population dynamics from molecular sequences

Comparison of Bayesian Coalescent Skyline Plot Models for Inferring Demographic Histories

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

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

Accelerated Bayesian inference of population size history from recombining sequence data

An Efficient Coalescent Epoch Model for Bayesian Phylogenetic Inference

The effects of sampling strategy on the quality of reconstruction of viral population dynamics using Bayesian skyline family coalescent methods: A simulation study

Bayesian Estimation of Past Population Dynamics in BEAST 1.10 Using the Skygrid Coalescent Model

Understanding Past Population Dynamics: Bayesian Coalescent-Based Modeling with Covariates

Bayesian Inference of Dependent Population Dynamics in Coalescent Models

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

Inferring the demographic history from DNA sequences: An importance sampling approach based on non-homogeneous processes

An Efficient Bayesian Inference Framework for Coalescent-Based Nonparametric Phylodynamics

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

Bayesian Inference of Sampled Ancestor Trees for Epidemiology and Fossil Calibration

Bayesian Inference of Pathogen Phylogeography using the Structured Coalescent Model

Bayesian Inference of Species Trees from Multilocus Data

Inferring population history from molecular phylogenies

Multiple merger coalescent inference of effective population size

Bayesian Inference of Natural Selection from Allele Frequency Time Series

A Bayesian phylogenetic method to estimate unknown sequence ages