《Biology》：Temporal Transcriptomic Responses of Turbinaria peltata to Repeated Marine Heatwaves and Recovery Cycles
编辑推荐：
海洋热浪（MHWs）正对珊瑚礁生态系统构成日益严重的威胁，然而珊瑚宿主在重复热胁迫下产生适应的转录组机制仍知之甚少。研究人员利用从头RNA测序技术，追踪了石珊瑚盾形陀螺珊瑚（Turbinaria peltata）在两个连续的MHW-恢复周期中的时间转录组响应。
海洋热浪（MHWs）正对珊瑚礁生态系统构成日益严重的威胁，然而珊瑚宿主在重复热胁迫下产生适应的转录组机制仍知之甚少。研究人员利用从头RNA测序技术，追踪了石珊瑚盾形陀螺珊瑚（Turbinaria peltata）在两个连续的MHW-恢复周期中的时间转录组响应。第二次热浪引发的差异表达基因（DEGs）数量少于第一次（279个对379个），而第二个恢复阶段则表现出最大的DEG集合（399个）。在功能上，第一次热浪触发了急性应激防御和信号抑制，而第一次恢复阶段涉及细胞外基质重塑和钙介导的组织修复，表明恢复不完全。第二个恢复阶段出现最大的DEG集合表明，反复的热胁迫施加了累积的生理成本，需要更广泛的分子重组——涵盖代谢重编程、mTOR和Wnt信号通路以及持续的氧化应激管理。结合先前对该物种的生理学和微生物组证据，研究结果支持一个模型，即重复热暴露逐渐将珊瑚宿主从急性休克防御反应转向代偿性代谢和结构修复。因此，在重复MHW模式下，恢复成为一种活跃但代价日益高昂的分子调整过程。这些发现确定了用于评估在不断加剧的MHW条件下珊瑚适应能力的候选分子标记物。
研究背景、现存问题与研究目的
全球气候变化导致海洋热浪（MHWs）的频率和强度不断增加，对珊瑚礁生态系统构成了严重威胁。尽管已有大量研究关注珊瑚在单次热胁迫下的分子响应，但关于珊瑚在重复热浪事件之间的恢复期如何动态调节其转录组，以及恢复是完全的还是逐渐受损的，目前认知十分有限。特别是对于位于边缘和高纬度、历史上被视为热避难所但如今正经历更频繁热浪的珊瑚礁系统而言，阐明恢复期的动态至关重要。本研究旨在填补这一空白，以南海北部徐闻珊瑚礁国家级自然保护区的优势造礁珊瑚——盾形陀螺珊瑚（Turbinaria peltata）为模型，探究其在两次连续热浪与恢复周期中的宿主转录组时间动态变化，揭示重复热胁迫下珊瑚适应性的分子机制。
研究内容与结论
研究人员通过模拟两个连续的MHW-恢复周期，对T. peltata进行RNA测序分析。研究发现，第二次热浪引发的差异表达基因（DEGs）数量（279个）少于第一次（379个），表明转录响应发生了重组而非单纯减弱。然而，第二次恢复阶段出现了最多的DEGs（399个），表明重复热胁迫产生了累积的生理成本，需要更广泛的分子重组，包括代谢重编程、mTOR和Wnt信号通路的参与以及持续的氧化应激管理。这揭示了恢复并非被动回归基线，而是一个活跃且代价递增的过程。结合前期生理学与微生物组证据，该研究支持一个模型：重复热暴露使珊瑚宿主从急性休克防御反应逐步转向代偿性代谢与结构修复，且这种适应性伴随着累积的恢复成本。
关键技术与方法
- 1.
实验设计与样品处理：将采自徐闻珊瑚礁国家级自然保护区的T. peltata样本，在实验室驯化后分为对照组（C）和热浪组（H）。实验包含两个连续的MHW-恢复周期：第1-5天为第一次热浪（升温至32℃），第6-10天为第一次恢复（降温至26℃），第11-15天为第二次热浪（再次升温至32℃），第16-20天为第二次恢复（再次降温至26℃）。分别在第5、10、15、20天（记为A、B、C、D阶段）采集样本。
- 2.
转录组测序与分析：提取珊瑚组织总RNA，采用Illumina NovaSeq 6000平台进行de novo RNA测序。使用Trinity进行转录本组装，DIAMOND软件比对NR、Swiss-Prot、COG、KOG、eggNOG和KEGG数据库进行功能注释。利用DESeq2鉴定差异表达基因（DEGs，筛选标准为|log2FC| ≥ 1且p < 0.01），并使用TBtools进行GO和KEGG富集分析。
- 3.
共表达趋势分析：利用BMKCloud平台的共表达模块，基于欧氏距离对四个实验阶段的共有基因进行K-means聚类，以识别不同的表达模式。
- 4.
RT-qPCR验证：选取8个代表性DEGs，以β-actin为内参基因，采用2?ΔΔCt法进行定量验证，以确保转录组数据的可靠性。
研究结果
3.1 转录组测序质量与功能注释：测序获得153.32 Gb高质量数据，经从头组装得到99,972条unigenes，N50长度为1029 bp。通过与多个公共数据库比对，成功注释了大量unigenes。
3.2 重复热胁迫与恢复过程中共有基因的表达趋势：通过K-means聚类将共有基因分为四种趋势。趋势1（169个基因）在恢复期高表达、热浪期低表达，富集于色氨酸代谢和视黄醇代谢通路。趋势2（196个基因）仅在第二次热浪期达到峰值，富集于Wnt和mTOR信号通路。趋势3（218个基因）表达稳定，富集于ECM-受体相互作用和TGF-beta信号通路。趋势4（315个基因）从第一次恢复期开始上调并持续至第二次恢复期，富集于ECM-受体相互作用。
3.3 差异表达基因（DEGs）统计：第一次热浪后共鉴定出379个DEGs（229个下调，150个上调）；第一次恢复后减少至240个（105个下调，135个上调）；第二次热浪后为279个（116个下调，163个上调）；第二次恢复后增至最多，达399个（197个下调，202个上调）。
3.4 DEGs的GO富集分析：第一次热浪显著富集15个GO条目，包括谷氨酰胺生物合成过程和生长因子活性等，多数与下调基因相关。第一次恢复期虽DEGs最少，但富集到最多GO条目（171个），涉及细胞粘附、突触囊泡回收和细胞-细胞信号传导等。第二次热浪仅富集4个分子功能GO条目，如清道夫受体活性和血红素结合。第二次恢复期富集70个GO条目，包括生长因子活性和钙离子结合等。
3.5 DEGs的KEGG通路富集分析：第一次热浪显著富集6条通路，包括ECM-受体相互作用和神经活性配体-受体相互作用，多数基因下调。第一次恢复期富集9条通路，包括ECM-受体相互作用和亚油酸代谢，多数基因上调。第二次热浪仅显著富集唾液分泌通路。第二次恢复期富集16条通路，包括视黄醇代谢、Wnt和mTOR信号通路。此外，55个在所有阶段共有的DEGs富集于唾液分泌通路，其中大部分基因在第二次热浪和第二次恢复期表现出更强的响应幅度。
3.6 代表性DEGs分析：对32个与免疫、能量代谢、氧化应激和钙结合相关的代表基因进行分析。热休克蛋白Hsp20在所有阶段持续差异表达。钙离子结合、免疫和营养调节相关基因主要在第一次恢复期上调。含EGF样结构域、EF-hand基序和血栓反应蛋白1型结构域的基因，以及编码交替氧化酶、细胞色素P450等的基因主要在第二次热浪期间上调。含氧化还原酶NAD结合结构域的基因在两个热浪事件中均显著表达。
3.7 RT-qPCR验证：对8个随机选择的DEGs进行RT-qPCR验证，结果显示其表达趋势与转录组测序数据基本一致，证实了测序结果的可靠性。
讨论总结与结论翻译
讨论总结：讨论部分指出，第二次热浪引发的DEGs少于第一次，表明转录响应发生了重组，这与“生态记忆”概念相符，区别于经典的启动模型。恢复期在此过程中扮演着主动角色，而非被动回归基线。第一次热浪引发了剧烈的急性应激反应，伴随信号通路抑制和生理损伤。第二次热浪虽然生理损伤减轻，但转录谱呈现出重新编程的调控状态，转向代偿性代谢调整。第一次恢复期显示部分转录恢复但不完全，涉及细胞外基质重塑和解毒等耗能过程。值得注意的是，第二次恢复期DEGs数量超过所有热浪阶段，表明重复热胁迫的累积成本超过了首次恢复建立的能力，需要更大规模的转录重组。研究还识别了Hsp20、钙信号、细胞外基质重塑和氧化还原调节等相关基因作为热适应性的潜在分子组分。与红海耐热珊瑚相比，T. peltata的响应不符合经典启动模型，而是表现为一种重新组织的转录状态，强调了珊瑚热响应的物种特异性和区域热历史依赖性。
结论翻译：总之，本研究表明，重复海洋热浪会诱导T. peltata产生阶段特异性的转录重编程。初次响应涉及广泛应激防御激活和不完全的生理恢复，而再次暴露则触发了一种重新组织的转录响应，伴随生理损伤减轻。值得注意的是，最大数量的差异表达基因出现在第二次恢复阶段，这强调了恢复是一个活跃的、消耗资源的过程，涉及细胞修复、代谢重组和全生物体调整，而非被动的基线回归。结合先前的生理学和微生物组证据，包括虫黄藻密度、叶绿素a含量和微生物群落的不完全恢复，这些发现表明，T. peltata的热适应性是通过连续热浪-恢复周期中的协调调控实现的。这项工作加深了对珊瑚应对不断加剧的热干扰的理解，并强调了代谢和离子稳态通路可作为评估未来气候情景下珊瑚适应能力的潜在指标。
