English

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

Computation 2019-05-30 v1 Populations and Evolution Methodology

Abstract

Calculation of the log-likelihood stands as the computational bottleneck for many statistical phylogenetic algorithms. Even worse is its gradient evaluation, often used to target regions of high probability. Order O(N){\cal O}\hspace{-0.2em}\left( N \right)-dimensional gradient calculations based on the standard pruning algorithm require O(N2){\cal O}\hspace{-0.2em}\left( N^2 \right) operations where N is the number of sampled molecular sequences. With the advent of high-throughput sequencing, recent phylogenetic studies have analyzed hundreds to thousands of sequences, with an apparent trend towards even larger data sets as a result of advancing technology. Such large-scale analyses challenge phylogenetic reconstruction by requiring inference on larger sets of process parameters to model the increasing data heterogeneity. To make this tractable, we present a linear-time algorithm for O(N){\cal O}\hspace{-0.2em}\left( N \right)-dimensional gradient evaluation and apply it to general continuous-time Markov processes of sequence substitution on a phylogenetic tree without a need to assume either stationarity or reversibility. We apply this approach to learn the branch-specific evolutionary rates of three pathogenic viruses: West Nile virus, Dengue virus and Lassa virus. Our proposed algorithm significantly improves inference efficiency with a 126- to 234-fold increase in maximum-likelihood optimization and a 16- to 33-fold computational performance increase in a Bayesian framework.

Keywords

Cite

@article{arxiv.1905.12146,
  title  = {Gradients do grow on trees: a linear-time ${\cal O}\hspace{-0.2em}\left( N \right)$-dimensional gradient for statistical phylogenetics},
  author = {Xiang Ji and Zhenyu Zhang and Andrew Holbrook and Akihiko Nishimura and Guy Baele and Andrew Rambaut and Philippe Lemey and Marc A. Suchard},
  journal= {arXiv preprint arXiv:1905.12146},
  year   = {2019}
}