arrow
Return

CherryML: scalable maximum likelihood estimation of phylogenetic models

delete2023-06-29
delete2
delete
OA
AI
S
Sebastian Prillo
Y
Yun Deng
P
Pierre Boyeau
X
Xingyu Li
P
Po‐Yen Chen
Y
Yun S. Song *
DOI:10.1038/s41592-023-01917-9delete
deleteOriginal
deleteOriginal request for help
deleteShare
deleteSave
Abstract

Abstract

En 中文
Phylogenetic models of molecular evolution are central to numerous biological applications spanning diverse timescales, from hundreds of millions of years involving orthologous proteins to just tens of days relating to single cells within an organism. A fundamental problem in these applications is estimating model parameters, for which maximum likelihood estimation is typically employed. Unfortunately, maximum likelihood estimation is a computationally expensive task, in some cases prohibitively so. To address this challenge, we here introduce CherryML, a broadly applicable method that achieves several orders of magnitude speedup by using a quantized composite likelihood over cherries in the trees. The massive speedup offered by our method should enable researchers to consider more complex and biologically realistic models than previously possible. Here we demonstrate CherryML's utility by applying it to estimate a general 400 x 400 rate matrix for residue-residue coevolution at contact sites in three-dimensional protein structures; we estimate that using current state-of-the-art methods such as the expectation-maximization algorithm for the same task would take >100,000 times longer. CherryML is a method to scale up maximum likelihood estimation for general phylogenetic models of molecular evolution, providing several orders of magnitude speedup over traditional methods.
Keywords:
PROTEIN EVOLUTION
RATES

Journal

Nature Methods cover
Nature Methods
IF:
32.1
Papers:
7.2K
Citations:
12.7W

Organization

U
University of California Berkeley
Scholars:
3.5W
Papers: 2.8W
Citations: 11.3W
University of California System cover
University of California System
Scholars:
37.5W
Papers: 33.7W
Citations: 6.6K