Abstract:Motivation: Reconstructing evolutionary histories of biological entities, such as genes, cells, organisms, populations, and species, from phenotypic and molecular sequencing data is central to many biological, palaeontological, and biomedical disciplines. Typically, due to uncertainties and incompleteness in data, the true evolutionary history (phylogeny) is challenging to estimate. Statistical modelling approaches address this problem by introducing and studying probability distributions over all possible evolutionary histories, but can also introduce uncertainties due to misspecification. In practice, computational methods are deployed to learn those distributions typically by sampling them. This approach, however, is fundamentally challenging as it requires designing and implementing various statistical methods over a space of phylogenetic trees (or treespace).Although the problem of developing statistics over a treespace has received substantial attention in the literature and numerous breakthroughs have been made, it remains largely unsolved. The challenge of solving this problem is two-fold: A treespace has non-trivial often counter-intuitive geometry implying that much of classical Euclidean statistics does not immediately apply; many parametrisations of treespace with promising statistical properties are computationally hard, so they cannot be used in data analyses. As a result, there is no single conventional method for estimating even the most fundamental statistics over any treespace, such as mean and variance, and various heuristics are used in practice. Despite the existence of numerous tree summary methods to approximate means of probability distributions over a treespace based on its geometry, and the theoretical promise of this idea, none of the attempts resulted in a practical method for summarising tree samples. Results: In this paper we present a tree summary method along with useful properties of our chosen treespace while focusing on its impact on phylogenetic analyses of real datasets. We perform an extensive benchmark study and demonstrate that our method outperforms currently most popular methods with respect to a number of important "quality" statistics. Further, we apply our method to three empirical datasets ranging from cancer evolution to linguistics and find novel insights into corresponding evolutionary problems in all of them. We hence conclude that this treespace is a promising candidate to serve as a foundation for developing statistics over phylogenetic trees analytically, as well as new computational tools for evolutionary data analyses. Availability and implementation: An implementation is available at https://github.com/bioDS/Centroid-Code. Supplementary information: Supplementary data are available at Bioinformatics online.

Phylogenetic least squares estimation without genetic distances

A consistent least-squares criterion for calibrating edge lengths in phylogenetic networks

Maximum likelihood estimates of pairwise rearrangement distances

A Distance-Based Method for Inferring Phylogenetic Networks in the Presence of Incomplete Lineage Sorting

Combinatorics of least squares trees

Performance comparison between k-tuple distance and four model-based distances in phylogenetic tree reconstruction.

Bayesian inference of phylogenetic distances: revisiting the eigenvalue approach

Distangsd: Fast and Accurate Inference of Genetic Distances for Next-Generation Sequencing Data

Distance-based phylogenetic methods around a polytomy

Scalable distance-based phylogeny inference using divide-and-conquer

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

Bayesian Least-Squares Supertrees (BLeSS): flexible inference of large time-calibrated phylogenies

Genetic Programming For Maximum-Likelihood Phylogeny Inference

Distance to the stochastic part of phylogenetic varieties

Toward a semi-supervised learning approach to phylogenetic estimation

Genetic Distance for a General Non-Stationary Markov Substitution Process

Measuring Fit of Sequence Data to Phylogenetic Model: Gain of Power Using Marginal Tests

A mixed integer linear programming model to reconstruct phylogenies from single nucleotide polymorphism haplotypes under the maximum parsimony criterion

Estimating the mean in the space of ranked phylogenetic trees

Fast phylogeny reconstruction through learning of ancestral sequences

IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies