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.

Gradients Do Grow on Trees: A Linear-Time O(N)-Dimensional Gradient for Statistical Phylogenetics

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

Many-core algorithms for high-dimensional gradients on phylogenetic trees

Automatic differentiation is no panacea for phylogenetic gradient computation

Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion

Fast phylogeny reconstruction through learning of ancestral sequences

Many-core algorithms for statistical phylogenetics

Harnessing machine learning to guide phylogenetic-tree search algorithms

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

Random-effects substitution models for phylogenetics via scalable gradient approximations

Genetic Programming For Maximum-Likelihood Phylogeny Inference

FastTree: Computing Large Minimum Evolution Trees with Profiles instead of a Distance Matrix

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

Genomic data analysis in tree spaces

A simulation approach for change-points on phylogenetic trees

PhyloGen: Language Model-Enhanced Phylogenetic Inference via Graph Structure Generation

Recovering the tree-like trend of evolution despite extensive lateral genetic transfer: A probabilistic analysis

GeoPhy: Differentiable Phylogenetic Inference via Geometric Gradients of Tree Topologies

Estimating the mean in the space of ranked phylogenetic trees

Principal component analysis and the locus of the Frechet mean in the space of phylogenetic trees