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.

Non-bifurcating Phylogenetic Tree Inference Via the Adaptive LASSO

A LASSO-based approach to sample sites for phylogenetic tree search

Reliable estimation of tree branch lengths using deep neural networks

Harnessing machine learning to guide phylogenetic-tree search algorithms

`Lassoing' a phylogenetic tree I: Basic properties, shellings, and covers

Fast phylogeny reconstruction through learning of ancestral sequences

Variational Bayesian Phylogenetic Inference with Semi-implicit Branch Length Distributions

Defining Loci in Restriction-Based Reduced Representation Genomic Data from Nonmodel Species: Sources of Bias and Diagnostics for Optimal Clustering

The length of haplotype blocks and signals of structural variation in reconstructed genealogies

Estimating phylogenetic trees from genome-scale data

Leaping through tree space: continuous phylogenetic inference for rooted and unrooted trees

Sequential Bayesian Phylogenetic Inference

Sapling: Inferring and Summarizing Tumor Phylogenies from Bulk Data using Backbone Trees

Data Requirement for Phylogenetic Inference from Multiple Loci: A New Distance Method

TreeShrink: fast and accurate detection of outlier long branches in collections of phylogenetic trees

A Fast and Scalable Method for Inferring Phylogenetic Networks from Trees by Aligning Lineage Taxon Strings

Inferring clonal evolution of tumors from single nucleotide somatic mutations

Accelerating Bayesian inference of dependency between complex biological traits

Bayesian Inference of Species Trees from Multilocus Data

Gradients do grow on trees: a linear-time ${\cal O}\hspace{-0.2em}\left( N \right)$-dimensional gradient for statistical phylogenetics

Predicting locus phylogenetic utility using machine learning