长寿造礁珊瑚Acropora palmata体细胞遗传变异镶嵌累积与遗传年龄估算研究

《Molecular Ecology》：Mosaic Accumulation of Somatic Genetic Variation and Estimates of Age in the Long-Lived Reef-Building Coral Acropora palmata

　　

体细胞遗传变异（SoGV）在生物体生命周期中持续累积，传统上因其与癌症和衰老的关联而被视为有害而非适应性特征。然而，在珊瑚等模块化生物中，有害突变可在细胞或多肽水平被清除，而适应性突变可能随着多肽形成遗传上独特的模块而提高频率。量化珊瑚中的体细胞遗传景观对于理解这些突变在珊瑚及克隆动物发育和进化中的作用至关重要。本研究对来自库拉索的8个Acropora palmata群体进行了体细胞遗传变异编目。通过全基因组测序（70–100×深度）记录了突变变异等位基因频率（VAF）随遗传体年龄的变化。在6至10岁的群体中观察到大量SoGV，并利用推断的突变率将一株年龄不确定的遗传体年龄估算至近一个世纪。突变未在多肽或分支水平固定，即其频率始终低于0.5，符合突变杂合位点的预期。然而，其等位基因频率遵循幂律分布，类似于衰老的人类组织。未发现正选择信号；相反，SoGV处于纯化选择之下。在一个群体中，利用SNP微阵列分析了沿分支28个样本中的突变。与预期相反，遗传距离与物理距离无关。结合观察到的缺乏固定现象，这一模式可能由大型干细胞池、干细胞去分化或休眠、强纯化选择的作用或上述因素的组合所解释。本研究结果为测试珊瑚模块水平选择的遗传变异提供了中性框架，探索了群体内物理距离与遗传距离的关系，并将体细胞遗传时钟应用于Acropora palmata群体。这项工作为造礁珊瑚体细胞突变景观提供了基础性见解，强调了研究体细胞遗传变异的重要性，因其可能有助于克隆动物的遗传多样性和适应性。

论文解读：造礁珊瑚Acropora palmata体细胞遗传变异的镶嵌累积与年龄估算

一、研究背景与科学问题

体细胞遗传变异（SoGV）在生物体生命周期中持续累积，传统观点认为其与癌症和衰老相关，因而被视为有害而非适应性因素。然而，在模块化生物中，情况可能有所不同：有害突变可在细胞或多肽水平被清除，而适应性突变可能随着多肽形成遗传上独特的模块而提高频率。珊瑚作为典型的模块化生物，其群体由大量多肽组成，多肽之间通过组织层共享资源与信号，但每个多肽可通过分裂产生新的多肽，从而在群体内形成遗传镶嵌。量化珊瑚中的体细胞遗传景观，对于理解这些突变在珊瑚及克隆动物发育和进化中的作用至关重要。此前，基于微卫星标记的逐步突变模型曾被用于估算珊瑚遗传体年龄，但该方法存在较大不确定性，且无法提供群体内体细胞突变动态的信息。此外，在植物（如海草）中，固定体细胞突变随时间线性增加，可用于构建体细胞分子钟；但珊瑚中是否存在固定突变尚不清楚，干细胞池大小对突变固定的影响也未知。因此，本研究旨在利用已知年龄的Acropora palmata群体，检验通用体细胞遗传时钟在珊瑚中的适用性，量化珊瑚体细胞突变的等位基因频率谱，并识别选择或中性的信号。

二、研究内容与主要结论

研究人员对来自库拉索两个礁点的7个已知年龄（6、8、10岁）的Acropora palmata群体及1个年龄不确定的“老”群体进行了全基因组测序（70–100×深度），每个群体采样5–6个模块（每个模块含3–4个多肽）。通过Mutect2鉴定体细胞突变，计算每个模块的突变负荷和突变率，并利用已知年龄群体校准体细胞分子钟，估算未知年龄群体的年龄。同时，研究人员分析了等位基因频率谱的分布特征，使用MOBSTER进行亚克隆去卷积，利用neutralitytestR检验中性偏离，并通过dndscv注释突变的功能影响。此外，还利用SNP微阵列分析了沿一个大型群体分支的28个样本，检验遗传距离与物理距离的相关性。

主要结论包括：（1）所有已知年龄群体的平均替换率为2.798e^?7 ± 1.977e^?7突变 bp^?1 module^?1 year^?1，据此估算未知年龄群体的年龄为72.09–92.30岁，中位数约80.77岁，远低于此前基于微卫星的估算（228–1924岁）。（2）在所有群体中均未检测到固定突变（VAF ≥ 0.5），突变频率均低于0.5，但等位基因频率谱随年龄发生变化：高频率突变密度随年龄增加，老群体的等位基因频率谱显著偏离中性1/f分布，更接近1/f²分布。（3）大多数SoGV为私有等位基因，但部分突变在多个模块间共享，共享突变频率随共享模块数增加而升高。（4）SNP微阵列分析显示，遗传距离与物理距离之间无显著线性关系（R² = 0.0714），邻接树拓扑也不一致。（5）编码区突变显著减少（比例最高0.0267，远低于基因组编码区比例0.1300），表明存在强纯化选择。

三、主要技术方法

研究人员主要采用了以下关键技术方法：（1）全基因组测序（WGS）：对每个群体采样5–6个模块，每个模块含3–4个多肽，进行配对末端测序（150 bp插入片段，目标深度80×）。（2）体细胞突变鉴定：使用Mutect2（GATK）以肿瘤-正常配对方式对所有样本进行两两比较，筛选SNP，过滤条件包括深度≥60、作图质量≥20、碱基质量≥25。（3）突变率计算与体细胞分子钟校准：基于每个模块的独特突变负荷和可调用碱基数计算突变率，利用已知年龄群体校准，估算未知年龄群体年龄。（4）等位基因频率谱分析：使用MOBSTER进行亚克隆去卷积，使用neutralitytestR检验中性偏离，使用dndscv注释突变功能影响。（5）SNP微阵列分析：使用Applied Biosystems Axiom Coral Genotyping Array—550962对沿分支的28个样本进行基因分型，计算Prevosti遗传距离并与物理距离进行线性回归。样本来源包括库拉索的Sea Aquarium、Marriott和Blue Bay三个礁点。

四、研究结果

4.1 校准替换负荷以估算Acropora palmata遗传体年龄

通过计算每个模块每碱基对的突变发生率，得到所有已知年龄群体的平均替换率为2.798e^?7 ± 1.977e^?7突变 bp^?1 module^?1 year^?1。利用该平均突变率，估算未知年龄群体各模块的年龄范围为72.09–92.30岁，中位数约80.77岁，显著低于此前微卫星估算的228–1924岁。

4.2 等位基因频率分布

在给定过滤条件下，未检测到固定SoGV（VAF ≥ 0.5且参考等位基因频率 ≤ 0.01）。MOBSTER分析检测到每个群体均存在两个亚克隆簇：低频率突变簇和高频率突变簇。老群体的高频率簇均值最高（0.25），但未检测到第三个簇，表明无正选择信号。已知年龄群体（6、8、10岁）的等位基因频率谱未显著偏离中性1/f分布，而老群体的等位基因频率谱显著偏离（AUC = 0.1423，K-S距离 = 0.2671），在高频率端更接近1/f²分布。此外，镶嵌区域（VAF 0.1–0.5）的突变密度随年龄显著增加（K-S检验，p < 0.05）。

4.3 群体内遗传多样性同质性

大多数SoGV为私有等位基因，但除S15群体外，所有群体均存在多模块共享突变。共享突变频率随共享模块数增加而升高。共现网络分析显示，M145群体中不同模块的连通性存在差异，M145_1模块连通性最高，而浅蓝色模块连通性最低。

4.4 等位基因频率的空间分布

私有突变主要为低频率变异；两个模块共享的突变包含镶嵌突变和低频率突变；三个或四个模块共享的突变主要为低频率和镶嵌突变。SNP微阵列分析显示，遗传距离与物理距离之间无显著线性关系（R² = 0.0714），邻接树拓扑不一致，路径评分差异为19.570。

4.5 体细胞突变累积的功能后果

所有检测到编码区突变的样本均表现出编码区突变的显著减少（χ²检验，p < 0.05），编码区突变比例最高为0.0267，远低于基因组编码区比例（0.1300）。CpG岛中突变比例平均为0.25%，且不随年龄显著变化（p > 0.05，ANOVA）。

五、讨论与结论

本研究表明，Acropora palmata的体细胞遗传变异随遗传体年龄增加而累积。基于平均SoGV累积速率2.798e^?7 ± 1.977e^?7突变 bp^?1 module^?1 year^?1，研究人员应用体细胞分子钟将未知年龄遗传体估算为约80.77岁，低于此前基于微卫星的估算（228–1194岁）。在6至10岁的群体中，80×以上深度的重复测序揭示了所有遗传体基因组中高水平的体细胞遗传变异。尽管突变未在多肽或分支水平固定，其等位基因频率遵循幂律分布，类似于衰老人类组织的模式，且未检测到正选择。相反，编码区突变的减少表明存在强纯化选择。遗传距离与物理距离无关的观察结果支持珊瑚可能从大型干细胞池中产生新模块，从而长期延迟突变固定。这些发现对珊瑚分子生物学之外具有更广泛的意义：理解体细胞突变如何在珊瑚群体内累积和维持，可改进受威胁礁系统中种群年龄结构、补充历史和适应潜力的估算，尤其是在长期种群数据有限的情况下。