Barba-Montoya Jose, Dos Reis Mario, Yang Ziheng
Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E 6BT, UK.
Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E 6BT, UK; School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, UK.
Mol Phylogenet Evol. 2017 Sep;114:386-400. doi: 10.1016/j.ympev.2017.07.005. Epub 2017 Jul 11.
Fossil calibrations are the utmost source of information for resolving the distances between molecular sequences into estimates of absolute times and absolute rates in molecular clock dating analysis. The quality of calibrations is thus expected to have a major impact on divergence time estimates even if a huge amount of molecular data is available. In Bayesian molecular clock dating, fossil calibration information is incorporated in the analysis through the prior on divergence times (the time prior). Here, we evaluate three strategies for converting fossil calibrations (in the form of minimum- and maximum-age bounds) into the prior on times, which differ according to whether they borrow information from the maximum age of ancestral nodes and minimum age of descendent nodes to form constraints for any given node on the phylogeny. We study a simple example that is analytically tractable, and analyze two real datasets (one of 10 primate species and another of 48 seed plant species) using three Bayesian dating programs: MCMCTree, MrBayes and BEAST2. We examine how different calibration strategies, the birth-death process, and automatic truncation (to enforce the constraint that ancestral nodes are older than descendent nodes) interact to determine the time prior. In general, truncation has a great impact on calibrations so that the effective priors on the calibration node ages after the truncation can be very different from the user-specified calibration densities. The different strategies for generating the effective prior also had considerable impact, leading to very different marginal effective priors. Arbitrary parameters used to implement minimum-bound calibrations were found to have a strong impact upon the prior and posterior of the divergence times. Our results highlight the importance of inspecting the joint time prior used by the dating program before any Bayesian dating analysis.
在分子钟定年分析中,化石校准是将分子序列间的距离解析为绝对时间和绝对速率估计值的最重要信息来源。因此,即使有大量分子数据,校准质量预计仍会对分歧时间估计产生重大影响。在贝叶斯分子钟定年中,化石校准信息通过分歧时间的先验(时间先验)纳入分析。在此,我们评估了三种将化石校准(以最小和最大年龄界限的形式)转换为时间先验的策略,这些策略根据是否从祖先节点的最大年龄和后代节点的最小年龄借用信息来为系统发育树上的任何给定节点形成约束而有所不同。我们研究了一个易于分析处理的简单示例,并使用三个贝叶斯定年程序(MCMCTree、MrBayes和BEAST2)分析了两个真实数据集(一个是10种灵长类物种,另一个是48种种子植物物种)。我们研究了不同的校准策略、出生 - 死亡过程和自动截断(以强制祖先节点比后代节点更古老的约束)如何相互作用来确定时间先验。一般来说,截断对校准有很大影响,以至于截断后校准节点年龄的有效先验可能与用户指定的校准密度非常不同。生成有效先验的不同策略也有相当大的影响,导致边际有效先验非常不同。发现用于实施最小界限校准的任意参数对分歧时间的先验和后验有强烈影响。我们的结果强调了在任何贝叶斯定年分析之前检查定年程序使用的联合时间先验的重要性。