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

TL;DR

This paper extends principal component analysis to the non-Euclidean space of phylogenetic trees by introducing a geometric object based on the Fréchet mean, enabling better visualization and analysis of complex biological data.

Contribution

It proposes a novel geometric approach for higher-order principal components in tree-space using the locus of the weighted Fréchet mean, with an efficient projection algorithm.

Findings

Algorithms perform well in simulations

Applied to gene trees and genome data

Reveals structure within complex phylogenetic data

Abstract

Most biological data are multidimensional, posing a major challenge to human comprehension and computational analysis. Principal component analysis is the most popular approach to rendering two- or three-dimensional representations of the major trends in such multidimensional data. The problem of multidimensionality is acute in the rapidly growing area of phylogenomics. Evolutionary relationships are represented by phylogenetic trees, and very typically a phylogenomic analysis results in a collection of such trees, one for each gene in the analysis. Principal component analysis offers a means of quantifying variation and summarizing a collection of phylogenies by dimensional reduction. However, the space of all possible phylogenies on a fixed set of species does not form a Euclidean vector space, so principal component analysis must be reformulated in the geometry of tree-space, which is a CAT(0) geodesic metric space. Previous work has focused on construction of the first principal component, or principal geodesic. Here we propose a geometric object which represents a $k$-th order principal component: the locus of the weighted Fr\'echet mean of $k+1$ points in tree-space, where the weights vary over the standard $k$-dimensional simplex. We establish basic properties of these objects, in particular that locally they generically have dimension $k$, and we propose an efficient algorithm for projection onto these surfaces. Combined with a stochastic optimization algorithm, this projection algorithm gives a procedure for constructing a principal component of arbitrary order in tree-space. Simulation studies confirm these algorithms perform well, and they are applied to data sets of Apicomplexa gene trees and the African coelacanth genome. The results enable visualizations of slices of tree-space, revealing structure within these complex data sets.