Département de Biochimie, Centre Robert Cedergren, Université de Montréal, Montréal, Québec, Canada.
Mol Biol Evol. 2011 Jan;28(1):729-44. doi: 10.1093/molbev/msq244. Epub 2010 Oct 6.
The comparative approach is routinely used to test for possible correlations between phenotypic or life-history traits. To correct for phylogenetic inertia, the method of independent contrasts assumes that continuous characters evolve along the phylogeny according to a multivariate Brownian process. Brownian diffusion processes have also been used to describe time variations of the parameters of the substitution process, such as the rate of substitution or the ratio of synonymous to nonsynonymous substitutions. Here, we develop a probabilistic framework for testing the coupling between continuous characters and parameters of the molecular substitution process. Rates of substitution and continuous characters are jointly modeled as a multivariate Brownian diffusion process of unknown covariance matrix. The covariance matrix, the divergence times and the phylogenetic variations of substitution rates and continuous characters are all jointly estimated in a Bayesian Monte Carlo framework, imposing on the covariance matrix a prior conjugate to the Brownian process so as to achieve a greater computational efficiency. The coupling between rates and phenotypes is assessed by measuring the posterior probability of positive or negative covariances, whereas divergence dates and phenotypic variations are marginally reconstructed in the context of the joint analysis. As an illustration, we apply the model to a set of 410 mammalian cytochrome b sequences. We observe a negative correlation between the rate of substitution and mass and longevity, which was previously observed. We also find a positive correlation between ω = dN/dS and mass and longevity, which we interpret as an indirect effect of variations of effective population size, thus in partial agreement with the nearly neutral theory. The method can easily be extended to any parameter of the substitution process and to any continuous phenotypic or environmental character.
比较方法通常用于测试表型或生活史特征之间可能存在的相关性。为了校正系统发育惯性,独立对比方法假设连续性状沿着系统发育按照多元布朗过程进化。布朗扩散过程也被用于描述替代过程参数的时间变化,例如替代率或同义替代与非同义替代的比值。在这里,我们开发了一个用于检验连续性状与分子替代过程参数之间耦合关系的概率框架。替代率和连续性状共同建模为未知协方差矩阵的多元布朗扩散过程。协方差矩阵、分歧时间以及替代率和连续性状的系统发育变异在贝叶斯蒙特卡罗框架中共同估计,协方差矩阵的先验共轭于布朗过程,以实现更高的计算效率。通过测量正或负协方差的后验概率来评估速率与表型之间的耦合关系,而分歧时间和表型变异则在联合分析的背景下进行边际重建。作为说明,我们将模型应用于 410 个哺乳动物细胞色素 b 序列。我们观察到替代率与质量和寿命之间存在负相关,这与之前的观察结果一致。我们还发现ω=dN/dS 与质量和寿命之间存在正相关,我们将其解释为有效种群大小变化的间接影响,因此在一定程度上与近中性理论一致。该方法可以很容易地扩展到替代过程的任何参数和任何连续的表型或环境特征。