Transmission, Infection, Diversification and Evolution Group, Max Planck Institute for the Science of Human History, Kahlaische Straße 10, Germany.
Organismal and Evolutionary Biology Research Programme, Faculty of Biological Sciences, University of Helsinki, P.O. Box 56, FI-00014 Helsinki, Finland.
Genome Biol Evol. 2022 Jul 2;14(7). doi: 10.1093/gbe/evac097.
Resolving the absolute timescale of phylogenetic trees stipulates reliable estimates for the rate of DNA sequence evolution. For this end, various calibration methods have been developed and studied intensively. Intraspecific rate variation among distinct genetic lineages, however, has gained less attention. Here, we have assessed lineage-specific molecular rates of human mitochondrial DNA (mtDNA) by performing tip-calibrated Bayesian phylogenetic analyses. Tip-calibration, as opposed to traditional nodal time stamps from dated fossil evidence or geological events, is based on sample ages and becoming ever more feasible as ancient DNA data from radiocarbon-dated samples accumulate. We focus on subhaplogroups U2, U4, U5a, and U5b, the data including ancient mtDNA genomes from 14C-dated samples (n = 234), contemporary genomes (n = 301), and two outgroup sequences from haplogroup R. The obtained molecular rates depended on the data sets (with or without contemporary sequences), suggesting time-dependency. More notable was the rate variation between haplogroups: U4 and U5a stand out having a substantially higher rate than U5b. This is also reflected in the divergence times obtained (U5a: 17,700 years and U5b: 29,700 years), a disparity not reported previously. After ruling out various alternative causes (e.g., selection, sampling, and sequence quality), we propose that the substitution rates have been influenced by demographic histories, widely different among populations where U4/U5a or U5b are frequent. As with the Y-chromosomal subhaplogroup R1b, the mitochondrial U4 and U5a have been associated with remarkable range extensions of the Yamnaya culture in the Bronze Age.
解决系统发育树的绝对时间尺度规定了 DNA 序列进化率的可靠估计。为此,已经开发并深入研究了各种校准方法。然而,不同遗传谱系之间的种内速率变化受到的关注较少。在这里,我们通过执行尖峰校准贝叶斯系统发育分析来评估人类线粒体 DNA(mtDNA)的谱系特异性分子速率。与基于化石证据或地质事件的传统节点时间戳不同,尖峰校准基于样本年龄,随着放射性碳定年样本的古老 DNA 数据的积累,变得越来越可行。我们专注于亚单倍群 U2、U4、U5a 和 U5b,数据包括来自 14C 定年样本的古老 mtDNA 基因组(n=234)、当代基因组(n=301)和来自单倍群 R 的两个外群序列。获得的分子速率取决于数据集(是否包含当代序列),表明存在时间依赖性。更值得注意的是单倍群之间的速率变化:U4 和 U5a 的速率明显高于 U5b。这也反映在获得的分歧时间上(U5a:17700 年和 U5b:29700 年),这是以前没有报道过的差异。在排除了各种替代原因(例如选择、采样和序列质量)之后,我们提出替代率受到人口之间广泛不同的人口统计学历史的影响,U4/U5a 或 U5b 频繁的人群。与 Y 染色体亚单倍群 R1b 一样,线粒体 U4 和 U5a 与青铜时代 Yamnaya 文化的显著范围扩展有关。