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.

Distances between Extension Spaces of Phylogenetic Trees

A combinatorial method for connecting BHV spaces representing different numbers of taxa

Genomic data analysis in tree spaces

Branch Length Transforms using Optimal Tree Metric Matching

A Metric on the Space of Rooted Phylogenetic Trees

Distances Between Phylogenetic Trees: A Survey

Comparison of phylogenetic trees defined on different but mutually overlapping sets of taxa: A review

Normalizing Kernels in the Billera-Holmes-Vogtmann Treespace

Distance-based phylogenetic methods around a polytomy

Estimating the mean in the space of ranked phylogenetic trees

Species, Clusters and the 'Tree of Life': A graph-theoretic perspective

Variational Bayesian Phylogenies through Matrix Representation of Tree Space

Novel metric for hyperbolic phylogenetic tree embeddings

Tropical Density Estimation of Phylogenetic Trees

Ranked Subtree Prune and Regraft

A New Distance Metric and Its Application in Phylogenetic Tree Construction

A distance-based model for convergent evolution

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

Uncertainty in phylogenetic tree estimates

Do branch lengths help to locate a tree in a phylogenetic network?