DOI: 10.1016/j.tpb.2019.11.005
论文题名: Fast likelihood calculation for multivariate Gaussian phylogenetic models with shifts
作者: Mitov V. ; Bartoszek K. ; Asimomitis G. ; Stadler T.
刊名: Theoretical Population Biology
ISSN: 405809
出版年: 2020
卷: 131 语种: 英语
英文关键词: Drift
; Jumps
; Lévy process
; Pruning
; Punctuated equilibrium
; Selection
Scopus关键词: algorithm
; Brownian motion
; evolution
; Gaussian method
; heterogeneity
; maximum likelihood analysis
; natural selection
; phylogenetics
; Animalia
; algorithm
; article
; calculation
; controlled study
; evolutionary rate
; expectation
; heart
; maximum likelihood method
; memory
; microorganism
; motion
; nonhuman
; phylogenetic tree
; probability
; quantitative analysis
; quantitative trait
英文摘要: Phylogenetic comparative methods (PCMs) have been used to study the evolution of quantitative traits in various groups of organisms, ranging from micro-organisms to animal and plant species. A common approach has been to assume a Gaussian phylogenetic model for the trait evolution along the tree, such as a branching Brownian motion (BM) or an Ornstein–Uhlenbeck (OU) process. Then, the parameters of the process have been inferred based on a given tree and trait data for the sampled species. At the heart of this inference lie multiple calculations of the model likelihood, that is, the probability density of the observed trait data, conditional on the model parameters and the tree. With the increasing availability of big phylogenetic trees, spanning hundreds to several thousand sampled species, this approach is facing a two-fold challenge. First, the assumption of a single Gaussian process governing the entire tree is not adequate in the presence of heterogeneous evolutionary forces acting in different parts of the tree. Second, big trees present a computational challenge, due to the time and memory complexity of the model likelihood calculation. Here, we explore a sub-family, denoted GLInv, of the Gaussian phylogenetic models, with the transition density exhibiting the properties that the expectation depends Linearly on the ancestral trait value and the variance is Invariant with respect to the ancestral value. We show that GLInv contains the vast majority of Gaussian models currently used in PCMs, while supporting an efficient (linear in the number of nodes) algorithm for the likelihood calculation. The algorithm supports scenarios with missing data, as well as different types of trees, including trees with polytomies and non-ultrametric trees. To account for the heterogeneity in the evolutionary forces, the algorithm supports models with “shifts” occurring at specific points in the tree. Such shifts can include changes in some or all parameters, as well as the type of the model, provided that the model remains within the GLInv family. This contrasts with most of the current implementations where, due to slow likelihood calculation, the shifts are restricted to specific parameters in a single type of model, such as the long-term selection optima of an OU process, assuming that all of its other parameters, such as evolutionary rate and selection strength, are global for the entire tree. We provide an implementation of this likelihood calculation algorithm in an accompanying R-package called PCMBase. The package has been designed as a generic library that can be integrated with existing or novel maximum likelihood or Bayesian inference tools. © 2019 The Author(s) Citation statistics:
资源类型: 期刊论文
标识符: http://119.78.100.158/handle/2HF3EXSE/159973
Appears in Collections: 气候变化与战略
There are no files associated with this item.
作者单位: Department of Biosystems Science and Engineering, Eidgenössische Technische Hochschule Zürich, Basel, Switzerland; Swiss Institute of Bioinformatics, Lausanne, Switzerland; Department of Computer and Information Science, Linköping University, Linköping, Sweden
Recommended Citation:
Mitov V.,Bartoszek K.,Asimomitis G.,et al. Fast likelihood calculation for multivariate Gaussian phylogenetic models with shifts[J]. Theoretical Population Biology,2020-01-01,131