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 Phylogeny and Its Impact on Evolutionary Biology

Bayesian inference of phylogenetic distances: revisiting the eigenvalue approach

Efficient Bayesian Inference of General Gaussian Models on Large Phylogenetic Trees

Sequential Bayesian Phylogenetic Inference

Inferring the historical patterns of biological evolution

Bayesian Inference of Evolutionary Histories under Time-Dependent Substitution Rates

Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10

Emerging Concepts of Data Integration in Pathogen Phylodynamics

Bayesian Inference of Species Trees from Multilocus Data

MRBAYES: Bayesian inference of phylogenetic trees

Bayesian inference in ecology

A Bayesian approach to inferring the phylogenetic structure of communities from metagenomic data

Online Bayesian phylodynamic inference in BEAST with application to epidemic reconstruction

Evaluating the Performance of Probabilistic Algorithms for Phylogenetic Analysis of Big Morphological Datasets: A Simulation Study

A Variational Approach to Bayesian Phylogenetic Inference

On the importance of assessing topological convergence in Bayesian phylogenetic inference

Bayesian Inference of Pathogen Phylogeography using the Structured Coalescent Model

Approximate Bayesian computation for Markovian binary trees in phylogenetics

Accelerating Bayesian inference of dependency between complex biological traits

Bayesian inference of mixed Gaussian phylogenetic models

MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice Across a Large Model Space