《Journal of Molecular Evolution》：Co-mutation Based Genetic Networks to Infer Temporal Mutation Dynamics in Ancient Human Mitochondrial Genomes
编辑推荐：
智人（Homo sapiens）的进化历史以环境、文化与遗传因素之间的复杂相互作用为特征。为探究这些过程的分子印记，研究人员利用次要等位基因共现（定义为共突变，co-mutation）原理，通过可变位点的时空共突变网络，分析了跨时间与地理背景的古代线粒体DNA
智人（Homo sapiens）的进化历史以环境、文化与遗传因素之间的复杂相互作用为特征。为探究这些过程的分子印记，研究人员利用次要等位基因共现（定义为共突变，co-mutation）原理，通过可变位点的时空共突变网络，分析了跨时间与地理背景的古代线粒体DNA（mtDNA）。基于单倍群（haplogroup）的可变位点评估揭示了铜器–青铜时代（Copper–Bronze Age）从采集狩猎向农耕生活方式的主要转变。遗传网络分析表明，COX和CYB基因座表现出不同的时间动态，其相互作用以依赖地质年代的方式受到NADH脱氢酶基因的调控。为补充网络方法，研究人员构建了基于系统发育的基因相互作用网络，并评估了与黑猩猩相比的多态性–分歧度比值。基于树的网络所呈现的拓扑结构与共突变分析一致，但基因–基因连通性降低。多态性/分歧度分析进一步表明，CYB基因长期处于净化选择（purifying selection）之下，而ATP6、COX及NADH脱氢酶基因则经历了与不同历史阶段相对应的间歇性净化选择。总体而言，研究结果表明，基于网络的古代mtDNA分析为早期人类生活方式转变及单倍群多样化提供了见解，有助于理解现代人群的进化基础。
论文解读：基于共突变网络的古代人类线粒体基因组时间动态研究
一、研究背景与问题
线粒体基因组（mtDNA）作为研究真核生物进化现象的紧凑而有力的模型，仅含37个基因，其中13个参与氧化磷酸化过程。由于其高拷贝数及相对易于从保存较差的化石标本中提取，古代mtDNA在重建人类历史——包括早期人类迁徙与混血事件——中发挥着关键作用。然而，从化石中提取和纯化DNA样品历来面临化学降解、微量含量及有机污染等困难。近年来技术进展使得从考古遗址获取更高质量的DNA序列成为可能，其中大部分为线粒体DNA。
人类线粒体基因组如同核基因组一样，获得了对不同环境的遗传适应，这为揭示时间和生物地理进化过程提供了线索。线粒体基因在代谢中的核心作用使得研究者能够利用古代mtDNA变异来推断营养可利用性如何影响人类迁徙及后续进化。人类亚群可按多种线粒体单倍群分类，这些单倍群代表人类线粒体系统发育树从中石器时代（Mesolithic）经新石器时代（Neolithic）、铜器时代（Copper）、青铜时代（Bronze）、铁器时代（Iron）至中世纪（Middle）的主要分支点。在该时期（公元前9755年至公元1264年），人类经历了由气候变化和农业发展加速所介导的文化与饮食转变。
一个关键的科学问题是：不同时代之间单倍群频率的转变是否反映了由向农耕生活方式转变所驱动的饮食选择效应，还是仅仅反映了有效群体大小波动导致的遗传漂变。此前尚不清楚从欧洲单倍群U的主导地位向当前主导的单倍群H的转变是否由农业介导的选择驱动。为探究这些问题，研究人员利用古代人类线粒体基因组数据库（AmtDB）中的完整mtDNA组装数据，结合新颖的共突变网络分析、常规系统发育重建以及人类多态性与人–黑猩猩分歧度比值分析，调查了过去约一万年人类历史中mtDNA变异的模式。
二、主要关键技术方法
研究人员从AmtDB数据库下载了812条具有C14测年数据的完整mtDNA FASTA序列，按地质年代分组。使用Clustal Omega进行序列比对，去除含插入/缺失的序列，以人类线粒体参考序列（rCRS）为基准确定可变位点。构建共突变网络：节点代表可变位点，边代表共突变频率，通过10,000次随机置换评估统计显著性（Pij<0.05），并映射至基因水平构建基因–基因相互作用（GGI）网络。计算网络度量指标包括平均度、聚类系数、模块性（Louvain算法）及共突变位点比例。利用MEGA11以最大似然法构建13个蛋白编码基因的基因树，通过加权Robinson-Foulds距离构建系统发育基因网络，并与共突变网络进行NetSimile、谱相似性和图结构相似性比较。采用McDonald-Kreitman（MK）检验计算中性指数（NI），以黑猩猩panTro6参考线粒体为外群评估选择压力。此外，使用gargammel模拟工具生成模拟古代mtDNA数据作为对照分析。
三、研究结果
可变位点确定：共获得812条古代mtDNA序列，涵盖六个时代。中石器时代可变位点中祖先SNP比例最高（219个位点中占6.8%），中世纪最低（580个位点中占1.0%）。
共突变网络：各时代共突变网络初始均为非连通状态，分析聚焦于最大连通分量（LCC）。参与共突变的可变位点比例从中石器时代到中世纪逐渐增加。ND4L基因的参与度变化显著——从中石器时代近50%的非共突变位点降至中世纪的不足10%。除铜器时代和中世纪外，各时代均观察到节点度与聚类系数之间的负相关。所有时代均呈现高模块性，表明网络组织为紧密连接的簇。
单倍群时间动态：中石器时代超过80%的样本属于单倍群U。新石器时代单倍群U仍为最频繁（27%），伴随10个其他单倍群的出现，反映mtDNA多样性增加。中世纪检测到14个主要单倍群。单倍群U在中石器、新石器及青铜时代占主导，而单倍群H在铜器、铁器及中世纪成为优势单倍群。
枢纽/全局节点及基因互作的时间变化：共鉴定出8个枢纽（高连接度）节点。中石器时代的G9477A变异（CO3基因）为非同义改变，其余时代的枢纽节点均为复合体III和IV中的同义改变。鉴定出25个全局节点（存在于所有时代网络中的节点）。CO1、CO3、CYB、ND2和ND4五个基因的度在不同时代呈现显著变异。新石器时代度分布广泛，铜器时代出现显著降低并保持相对稳定，但CYB基因例外——其度在铜器时代后急剧增加。
遗传互作网络与系统发育基因网络：各时期共突变减少型遗传网络中的互作数量均多于共突变增强型。共突变减少型网络与系统发育树基网络表现出显著更高的相似性。ND基因（复合体I）在各时代均作为中心枢纽出现，常与新石器及铜器时代的CO基因（复合体IV）组合。中世纪未观察到任何基因组合的强共突变增强或减少关联。基于Jaccard距离的聚类树状图显示，中石器与新石器时代聚为一组，其余四个时代聚为另一组；铁器和中世纪距中石器和新石器时代最远，铜器与青铜时代彼此最为接近。
时间动态选择与中性指数：中世纪是唯一检测到显著正选择（positive selection）的时期，涉及ATP8、CYB、ND4L和ND6基因，且该时期任何蛋白编码基因均未检测到显著净化选择。中世纪之前，ATP8、ND4L和ND6为中性进化，CYB大多处于净化选择之下。ATP6、CO1、ND2和ND5在多数时代处于净化选择之下，CO3和ND1多为中性伴部分净化选择，CO2、ND3和ND4始终为中性。
四、讨论与结论
讨论部分指出，中石器时代单倍群U的主导地位与其与狩猎采集群体的紧密关联一致。新石器及铜器时代单倍群H和K的出现标志着向农耕的转变。中世纪祖先SNP比例最低反映了更强的遗传漂变和群体扩张。共突变位点比例的逐步增加表明农业社会崛起伴随更大的遗传同质化。共突变减少型网络与系统发育树基网络的相似性表明其反映共享祖先信息，而共突变增强型网络可能捕获系统发育上独立的、与时间和生活方式相关的功能性耦合。CYB基因从中石器时代的树基和共突变减少型高连接节点，到铁器时代转变为共突变增强型高连接节点并持续至中世纪，这一转变与正选择的检测相关联。ATP6基因虽不含枢纽节点变异，但在多个时代均显示显著且强烈的净化选择。
研究结论总结如下：共突变网络有效捕获了人类进化各时期线粒体基因组变异的时间动态。CO1和CO3持续作为高连接度节点出现，凸显细胞色素氧化酶复合体在塑造饮食和生活方式相关进化轨迹中的核心作用。从中石器到中世纪共突变可变位点的逐步增加提示向农耕生活方式的转变伴随遗传同质化趋势。减少型共突变网络与谱系特异性进化约束一致，而增强型网络可能捕获功能性相关耦合。模拟对照分析表明这些网络模式不能仅由古代DNA损伤所重现。这些发现确立了共突变网络作为一种强大框架，与系统发育方法整合后，可为理解塑造线粒体基因组动态的功能、上位性和进化力量之间的相互作用提供细致见解。该研究发表于《Journal of Molecular Evolution》。
