Dataset Information

Phylogenetic Curved Optimal Regression for Adaptive Trait Evolution.

ABSTRACT: Regression analysis using line equations has been broadly applied in studying the evolutionary relationship between the response trait and its covariates. However, the characteristics among closely related species in nature present abundant diversities where the nonlinear relationship between traits have been frequently observed. By treating the evolution of quantitative traits along a phylogenetic tree as a set of continuous stochastic variables, statistical models for describing the dynamics of the optimum of the response trait and its covariates are built herein. Analytical representations for the response trait variables, as well as their optima among a group of related species, are derived. Due to the models' lack of tractable likelihood, a procedure that implements the Approximate Bayesian Computation (ABC) technique is applied for statistical inference. Simulation results show that the new models perform well where the posterior means of the parameters are close to the true parameters. Empirical analysis supports the new models when analyzing the trait relationship among kangaroo species.

SUBMITTER: Jhwueng DC

PROVIDER: S-EPMC7916804 | biostudies-literature | 2021 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Phylogenetic Curved Optimal Regression for Adaptive Trait Evolution.

Jhwueng Dwueng-Chwuan DC Wang Chih-Ping CP

Entropy (Basel, Switzerland) 20210210 2

Regression analysis using line equations has been broadly applied in studying the evolutionary relationship between the response trait and its covariates. However, the characteristics among closely related species in nature present abundant diversities where the nonlinear relationship between traits have been frequently observed. By treating the evolution of quantitative traits along a phylogenetic tree as a set of continuous stochastic variables, statistical models for describing the dynamics o ...[more]

PMID: 33579023

Similar Datasets

Project description:A number of methods have been developed for modeling the evolution of a quantitative trait on a phylogeny. These methods have received renewed interest in the context of genome-wide studies of gene expression, in which the expression levels of many genes can be modeled as quantitative traits. We here develop a new method for joint analyses of quantitative traits within- and between species, the Expression Variance and Evolution (EVE) model. The model parameterizes the ratio of population to evolutionary expression variance, facilitating a wide variety of analyses, including a test for lineage-specific shifts in expression level, and a phylogenetic ANOVA that can detect genes with increased or decreased ratios of expression divergence to diversity, analogous to the famous Hudson Kreitman Aguadé (HKA) test used to detect selection at the DNA level. We use simulations to explore the properties of these tests under a variety of circumstances and show that the phylogenetic ANOVA is more accurate than the standard ANOVA (no accounting for phylogeny) sometimes used in transcriptomics. We then apply the EVE model to a mammalian phylogeny of 15 species typed for expression levels in liver tissue. We identify genes with high expression divergence between species as candidates for expression level adaptation, and genes with high expression diversity within species as candidates for expression level conservation and/or plasticity. Using the test for lineage-specific expression shifts, we identify several candidate genes for expression level adaptation on the catarrhine and human lineages, including genes putatively related to dietary changes in humans. We compare these results to those reported previously using a model which ignores expression variance within species, uncovering important differences in performance. We demonstrate the necessity for a phylogenetic model in comparative expression studies and show the utility of the EVE model to detect expression divergence, diversity, and branch-specific shifts.

Project description:Understanding the processes that give rise to quantitative measurements associated with molecular sequence data remains an important issue in statistical phylogenetics. Examples of such measurements include geographic coordinates in the context of phylogeography and phenotypic traits in the context of comparative studies. A popular approach is to model the evolution of continuously varying traits as a Brownian diffusion process acting on a phylogenetic tree. However, standard Brownian diffusion is quite restrictive and may not accurately characterize certain trait evolutionary processes. Here, we relax one of the major restrictions of standard Brownian diffusion by incorporating a nontrivial estimable mean into the process. We introduce a relaxed directional random walk (RDRW) model for the evolution of multivariate continuously varying traits along a phylogenetic tree. Notably, the RDRW model accommodates branch-specific variation of directional trends while preserving model identifiability. Furthermore, our development of a computationally efficient dynamic programming approach to compute the data likelihood enables scaling of our method to large data sets frequently encountered in phylogenetic comparative studies and viral evolution. We implement the RDRW model in a Bayesian inference framework to simultaneously reconstruct the evolutionary histories of molecular sequence data and associated multivariate continuous trait data, and provide tools to visualize evolutionary reconstructions. We demonstrate the performance of our model on synthetic data, and we illustrate its utility in two viral examples. First, we examine the spatiotemporal spread of HIV-1 in central Africa and show that the RDRW model uncovers a clearer, more detailed picture of the dynamics of viral dispersal than standard Brownian diffusion. Second, we study antigenic evolution in the context of HIV-1 resistance to three broadly neutralizing antibodies. Our analysis reveals evidence of a continuous drift at the HIV-1 population level towards enhanced resistance to neutralization by the VRC01 monoclonal antibody over the course of the epidemic. [Brownian Motion; Diffusion Processes; Phylodynamics; Phylogenetics; Phylogeography; Trait Evolution.].

Project description:Orchidaceae is the 3rd largest family of angiosperms, an evolved young branch of monocotyledons. This family contains a number of economically-important horticulture and flowering plants. However, the limited availability of genomic information largely hindered the study of molecular evolution and phylogeny of Orchidaceae. In this study, we determined the evolutionary characteristics of whole chloroplast (cp) genomes and the phylogenetic relationships of the family Orchidaceae. We firstly characterized the cp genomes of four orchid species: Cremastra appendiculata, Calanthe davidii, Epipactis mairei, and Platanthera japonica. The size of the chloroplast genome ranged from 153,629 bp (C. davidi) to 160,427 bp (E. mairei). The gene order, GC content, and gene compositions are similar to those of other previously-reported angiosperms. We identified that the genes of ndhC, ndhI, and ndhK were lost in C. appendiculata, in that the ndh I gene was lost in P. japonica and E. mairei. In addition, the four types of repeats (forward, palindromic, reverse, and complement repeats) were examined in orchid species. E. mairei had the highest number of repeats (81), while C. davidii had the lowest number (57). The total number of Simple Sequence Repeats is at least 50 in C. davidii, and, at most, 78 in P. japonica. Interestingly, we identified 16 genes with positive selection sites (the psbH, petD, petL, rpl22, rpl32, rpoC1, rpoC2, rps12, rps15, rps16, accD, ccsA, rbcL, ycf1, ycf2, and ycf4 genes), which might play an important role in the orchid species' adaptation to diverse environments. Additionally, 11 mutational hotspot regions were determined, including five non-coding regions (ndhB intron, ccsA-ndhD, rpl33-rps18, ndhE-ndhG, and ndhF-rpl32) and six coding regions (rps16, ndhC, rpl32, ndhI, ndhK, and ndhF). The phylogenetic analysis based on whole cp genomes showed that C. appendiculata was closely related to C. striata var. vreelandii, while C. davidii and C. triplicate formed a small monophyletic evolutionary clade with a high bootstrap support. In addition, five subfamilies of Orchidaceae, Apostasioideae, Cypripedioideae, Epidendroideae, Orchidoideae, and Vanilloideae, formed a nested evolutionary relationship in the phylogenetic tree. These results provide important insights into the adaptive evolution and phylogeny of Orchidaceae.

Dataset Information

Phylogenetic Curved Optimal Regression for Adaptive Trait Evolution.

Publications

Phylogenetic Curved Optimal Regression for Adaptive Trait Evolution.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets