环境驱动因素对大星珊瑚Montastraea cavernosa隐存谱系内与谱系间代谢组特征的影响

《Metabolomics》：Environmental drivers of metabolomic profiles within and between cryptic lineages of Montastraea cavernosa, the great star coral

　　

环境驱动因素对大星珊瑚Montastraea cavernosa隐存谱系内和谱系间代谢组特征的影响

摘要

引言：珊瑚礁恢复实践者旨在通过保护珊瑚礁和在珊瑚苗圃中培育多种基因型来保护珊瑚遗传多样性。然而，大多数珊瑚中的隐存遗传谱系使恢复策略复杂化，因为谱系间遗传分化在适应中的作用仍不清楚。在Montastraea cavernosa中，研究人员已鉴定出隐存谱系，其中一些按深度强烈分离。

目标：研究人员旨在量化对环境谱系特异性反应，以解释两个谱系之间的遗传和生态分化。

方法：研究人员进行了一项为期十周的双向移植实验，使用两个局限于浅水（小于10米深度）的隐存谱系，其中一个谱系在近岸珊瑚礁更为常见，另一个在离岸珊瑚礁。移植后，研究人员分析了珊瑚的生长和对移植的代谢组响应。

结果：令人惊讶的是，最强的响应并非谱系特异性。相反，两个谱系均表现出显著且相似的生长和代谢组特征变化，这取决于它们被移植到的最终栖息地。

结论：这些结果表明，尽管隐存谱系具有遗传差异，它们在代谢物水平上采用相似的适应和驯化机制来应对环境挑战。

随着人为胁迫加剧，珊瑚种群持续衰退，遗传多样性降低，削弱了珊瑚适应未来环境的能力。许多名义珊瑚物种内存在形态相似但基因交流受限的隐存遗传谱系（cryptic genetic lineages），这给珊瑚礁恢复中的遗传多样性管理带来复杂挑战。局部适应可能是驱动隐存遗传结构的环境因素，但维持这些谱系分布的环境参数尚不清楚。基因型×环境（G×E）互作能揭示遗传背景与环境条件如何共同影响珊瑚表型，包括生长、胁迫耐受性和代谢组特征。双向移植实验（reciprocal transplant experiment）是解析G×E互作的重要手段。大星珊瑚Montastraea cavernosa是分布于大西洋、加勒比海和墨西哥湾的常见造礁珊瑚，在佛罗里达群岛存在多个隐存谱系，其中两个浅水特化谱系MC4和MC7分别主导近岸和离岸环境。研究人员利用二者常同域分布的特点，设计双向移植实验，旨在量化这两个谱系对环境的特异性响应，从而解释其遗传与生态分化。

研究人员从佛罗里达群岛的Carly’s Patch和Haslun’s Heads两处珊瑚礁采集样本，经基因分型鉴定谱系后，选取MC4和MC7各15个群体，每群体分裂为两份，分别移植至两个地点，进行为期十周的相互移植。通过浮重测量评估生长速率，并利用非靶向代谢组学分析移植后的代谢组变化。结果发现，最强的响应并非谱系特异性，而是取决于最终移植地点：两个谱系在生长和代谢组特征上均发生显著且相似的变化。该发现表明，尽管隐存谱系存在遗传差异，它们在代谢水平上可能采用相似的适应与驯化机制，这为恢复实践中珊瑚基因型和移植地点的匹配提供了新认识：短期内在热应激条件下，移植地点的环境条件可能比谱系身份更关键。

关键技术方法如下：研究样本来自佛罗里达群岛Carly’s Patch和Haslun’s Heads两处珊瑚礁，最初采样160份组织，经2bRAD测序和参考基因组比对鉴定隐存谱系，剔除克隆和杂交个体后，每个谱系每地选取15个群体用于移植实验。移植装置使用PVC支架和防食网罩，HOBO温度记录仪连续监测水温。代谢组分析采用改良的100%乙醇保存液提取方案，以超高效液相色谱-加热电喷雾电离-串联质谱（UHPLC-HESI-MS/MS）进行非靶向检测；数据处理使用MZMine、GNPS分子网络、SIRIUS及NPClassifier等工具。统计推断基于线性混合模型、置换多元方差分析（PERMANOVA）和主坐标分析（PCoA）等方法。

环境监测：HOBO记录显示，Carly’s Patch的周平均温度略高于Haslun’s Heads，日平均温度差异显著（p<2×10^-16），表明两移植地点存在明显的热环境差异。

浮重（生长）：在控制初始大小影响后，三组克隆比较显示移植到Haslun’s Heads的片段平均增重显著大于移植到Carly’s Patch的片段，第四组也呈相同趋势但未显著。线性混合模型表明最终地点显著影响生长（F=19.12，p<0.001），而谱系、初始地点及其交互作用均不显著。估计边际均值进一步确认Haslun’s Heads的生长优势。

化学组成：置换多元方差分析显示，最终地点解释代谢组变异的24%（R²=0.24，p≤1×10^-6），谱系解释9%（R²=0.09，p≤1×10^-6），初始地点解释3%，处理效应不显著。火山图显示谱系比较有203个差异丰度代谢特征，最终地点比较有195个，初始地点比较72个，处理比较仅2个。热图显示样本主要按最终地点聚类。主坐标分析中PCo2主要由最终地点分离，PCo1主要由谱系分离，方差分解表明最终地点对PCo2的解释量远大于谱系，且无谱系×最终地点交互作用。因此，最终栖息地是代谢组特征最强的统计驱动因素。

讨论部分指出，生长结果与部分物种的局部驯化现象不同，而与Acropora millepora的移植结果类似，即移植到更有利栖息地的个体生长更好。Haslun’s Heads邻近佛罗里达海流，温度较低，可能作为热应激避难所，但实验期间正值热浪和区域白化事件，故需跨季节的长期实验验证。代谢组分析虽显示谱系和最终地点均有影响，但无法区分宿主、共生体或微生物组的贡献，建议未来结合微生物组基因表达研究。最终地点的代谢组差异与温度差异一致，但由于实验与急性热事件重合，尚不能判断是稳定的栖息地特征还是热暴露差异所致，应将代谢组响应解释为对短时热应激环境的反应。

与最初假设相反，无论初始来源地如何，移植到离岸环境的个体比移植到近岸环境的个体生长更快。此外，虽然谱系对珊瑚代谢组有显著影响，但最终栖息地也有显著影响。这些结果表明，该珊瑚物种当前所处的位置影响其化学表型和适应性。然而，由于实验与急性热应激事件重合且地点间存在温度差异，这些地点效应应被解释为热浪期间环境条件的响应，而非稳定栖息地特化的证据。从恢复角度而言，这可能意味着珊瑚移植的栖息地可能比其隐存谱系身份更重要，至少在热胁迫条件下的短期内如此。长期实验或可揭示谱系与栖息地互作对适应性的影响。