Abstract:The ideal approach to Bayesian phylogenetic inference is to estimate all parameters of interest jointly in a single hierarchical model. However, this is often not feasible in practice due to the high computational cost. Instead, phylogenetic pipelines generally consist of sequential analyses, whereby a single point estimate from a given analysis is used as input for the next analysis (e.g., a single multiple sequence alignment is used to estimate a gene tree). In this framework, uncertainty is not propagated from step to step, which can lead to inaccurate or spuriously confident results. Here, we formally develop and test a sequential inference approach for Bayesian phylogenetic inference, which uses importance sampling to generate observations for the next step of an analysis pipeline from the posterior distribution produced in the previous step. Our sequential inference approach presented here not only accounts for uncertainty between analysis steps, but also allows for greater flexibility in software choice (and hence model availability) and can be computationally more efficient than the traditional joint inference approach when multiple models are being tested. We show that our sequential inference approach is identical in practice to the joint inference approach only if sufficient information in the data is present (a narrow posterior distribution) and/or sufficiently many importance samples are used. Conversely, we show that the common practice of using a single point estimate can be biased, e.g., a single phylogeny estimate to transform an unrooted phylogeny into a time-calibrated phylogeny. We demonstrate the theory of sequential Bayesian inference using both a toy example and an empirical case study of divergence-time estimation in insects using a relaxed clock model from transcriptome data. In the empirical example, we estimate three posterior distributions of branch lengths from the same data (DNA character matrix with a GTR+Γ+I substitution model, an amino acid data matrix with empirical substitution models, and an amino acid data matrix with the PhyloBayes CAT-GTR model). Finally, we apply three different node-calibration strategies and show that divergence-time estimates are affected by both the data source and underlying substitution process to estimate branch lengths as well as the node-calibration strategies. Thus, our new sequential Bayesian phylogenetic inference provides the opportunity to efficiently test different approaches for divergence time estimation, including branch-length estimation from other software.

Bayesian inference of evolutionary histories under time-dependent substitution rates

Bayesian Inference of Evolutionary Histories under Time-Dependent Substitution Rates

Estimating evolutionary rates using time-structured data: a general comparison of phylogenetic methods

Model-based inference of punctuated molecular evolution

Scalable Bayesian divergence time estimation with ratio transformations

Examining the molecular clock hypothesis for the contemporary evolution of the rabies virus

Bayesian inference of phylogenetic distances: revisiting the eigenvalue approach

Dating the time of viral subtype divergence

Modeling Substitution Rate Evolution across Lineages and Relaxing the Molecular Clock

Online Bayesian phylodynamic inference in BEAST with application to epidemic reconstruction

Bayesian Selection of Relaxed-clock Models: Distinguishing Between Independent and Autocorrelated Rates

Bayesian random local clocks, or one rate to rule them all

Sequential Bayesian Phylogenetic Inference

Bayesian Inference of Pathogen Phylogeography using the Structured Coalescent Model

Bayesian mixture models for phylogenetic source attribution from consensus sequences and time since infection estimates

Random-effects substitution models for phylogenetics via scalable gradient approximations

Rate variation and recurrent sequence errors in pandemic-scale phylogenetics

Accelerating Bayesian inference of dependency between complex biological traits

Bayesian Inference of Phylogeny and Its Impact on Evolutionary Biology

Efficient Bayesian Inference of General Gaussian Models on Large Phylogenetic Trees

Integrating Transmission Dynamics and Pathogen Evolution Through a Bayesian Approach