高时间分辨率基因表达分析揭示佛罗里达群岛星芒滨珊瑚热耐受性的对比动态

《Molecular Ecology》：High Temporal Resolution Gene Expression Analysis Reveals Contrasting Dynamics of Thermal Tolerance in the Coral Porites astreoides From the Florida Keys

摘要 海洋温度上升威胁着珊瑚礁，然而珊瑚物种和群体的热耐受性各不相同。利用珊瑚白化自动应激系统（CBASS），研究人员检测了来自两个热环境截然不同的栖息地（一个耐受性更强的近岸群体与一个更敏感的离岸群体）的野生星芒滨珊瑚群体的热阈值和转录组动态。在为期19小时的急性热胁迫实验中的六个时间点进行RNA测序，并结合光合效率测量（Fv/Fm），揭示了近岸珊瑚（37.77°C）相较于离岸珊瑚（36.94°C）具有更高的标准化热耐受阈值，即有效剂量（ED50），这支持了先前在长期热胁迫实验中观察到的差异。重要的是，在施加的最低热胁迫温度34°C下，离岸群体始终表现出更多差异表达基因（DEGs），且其变化幅度高于近岸群体。例如，在胁迫高峰期，离岸珊瑚中有1136个DEGs以更高幅度表达，而近岸珊瑚仅有50个。主成分分析（PCA）显示温度是转录变异的主要驱动因素（解释了PC1上30%的方差），而群体起源则占较小但显著的比例（PC2上21%）。近岸珊瑚表现出平缓的转录组响应和快速恢复，而离岸珊瑚则表现出更大、延迟的响应。研究人员将这些响应差异解读为转录组弹性和驯化，这是与珊瑚热耐受性增强相关的两种新兴机制模式。未检测到基因表达前置的证据，表明表观遗传或转录后调控可能是所观察到的弹性差异的基础。研究结果证明了弹性机制中长期存在的、精细尺度的群体差异，这表明稳定的热耐受性策略可为未来的保护工作提供参考。

论文解读

研究背景

海洋温度持续升高对全球珊瑚礁生态系统构成严重威胁，导致珊瑚白化事件频发，进而引发大规模珊瑚死亡。尽管面临相同的气候变化压力，不同珊瑚物种乃至同一物种的不同地理群体却展现出迥异的热耐受性。这种差异背后的分子调控机制，特别是宿主珊瑚的基因表达响应动力学，是理解珊瑚能否适应未来变暖海洋的关键。以往的研究多局限于少数时间点的采样，难以捕捉热胁迫发生初期及恢复期关键的转录动态，如转录组弹性、基因表达前置和转录组驯化等过程。因此，迫切需要采用高时间分辨率的实验设计，深入探究不同热历史背景的珊瑚群体在应对急性热胁迫时的转录组响应差异及其潜在机制。

研究目的与意义

本研究旨在利用珊瑚白化自动应激系统（CBASS），对来自佛罗里达群岛热环境对比鲜明的两个星芒滨珊瑚群体（近岸耐受群体与离岸敏感群体）进行急性热胁迫实验，并通过高时间分辨率的转录组测序，系统比较两者在热胁迫及恢复过程中宿主基因表达的时序和幅度差异。研究成果不仅验证了短时CBASS实验能够复现长期实验观察到的群体间热耐受性差异，还揭示了近岸珊瑚可能通过“转录组弹性”和“转录组驯化”等机制实现更高效的胁迫响应与恢复，而非依赖于共生藻群落组成改变或基因表达前置。这些发现强调了宿主调控在珊瑚热耐受性中的核心作用，并为识别天然耐热群体作为未来珊瑚礁保护与修复的种质资源提供了重要理论依据。

关键技术方法

研究人员于2021年6月采集了来自佛罗里达群岛近岸和离岸两个位点的星芒滨珊瑚群体各5个基因型。利用珊瑚白化自动应激系统（CBASS）进行19小时急性热胁迫实验，设置30°C基线、34°C、37°C和39°C四个温度处理，并在实验开始前（TF）、基线稳定后（T0）、升温结束（T1）、高温维持结束（T2）、降温至基线后（T3）、夜间恢复（T4）和次日早晨（T5）共6个时间点取样进行mRNA测序（mRNA-Seq）。同时，利用脉冲振幅调制（PAM）荧光仪测量光系统II最大量子产量（F_v/F_m），用于拟合剂量-响应曲线计算半数有效剂量（ED50）。通过DESeq2进行差异表达基因（DEGs）分析和交互作用对比，利用主成分分析（PCA）评估全局基因表达模式，并计算转录组偏移量衡量响应与恢复动态。此外，通过ITS2扩增子测序鉴定共生藻（Symbiodiniaceae）群落组成。

研究结果

3.1 近岸群体比离岸群体具有更高的热阈值

通过将F_v/F_m测量值拟合三参数对数-逻辑斯蒂剂量-响应模型，研究发现近岸群体的ED50为37.77°C，显著高于离岸群体的36.94°C，差值达0.83°C。这表明近岸珊瑚的共生藻光化学系统在更高温度下才出现同等程度的损伤，验证了其更强的热耐受性。

3.2 RNA测序概况

经过滤和去重后，共获得216个可用样本（每个群体108个），平均每个样本获得约390万条去重复读数，用于后续基因表达定量分析。

3.3 近岸与离岸珊瑚基因表达响应的差异

交互作用分析显示，离岸群体在多数时间点和温度处理下拥有更多且变化幅度更大的DEGs。尤其在34°C处理的T3时间点（恢复初期），离岸群体有1136个DEGs变化幅度大于近岸群体，而近岸仅有50个，差异近23倍。随着恢复进行（T4、T5），两群体的DEGs数量均下降，但离岸群体始终更多，表明其恢复较慢。共享DEGs占绝大多数（>90.9%），说明两群体激活的核心胁迫响应通路相似，差异主要体现在响应幅度上。初步的基因本体（GO）富集分析在34°C处理的T2和T4时间点发现了少量显著富集的条目，涉及钙离子结合、细胞外区域和Wnt信号通路等。

3.4 近岸珊瑚表现出更快的胁迫响应恢复速度

主成分分析（PCA）显示，PC1（30%方差）主要由温度驱动，PC2（21%方差）主要由群体来源分离。PERMANOVA分析证实，在多个时间点和温度下，两群体的基因表达谱存在显著分离。转录组偏移量分析表明，在34°C处理下，近岸群体的偏移量在不同时间点间无显著差异，暗示其表达谱在整个实验期间相对稳定；而离岸群体的偏移量在T2和T3与T5相比显著更高，表明其经历了显著的表达波动后才逐渐恢复。在39°C极端胁迫下，两群体在恢复后期（T5）均表现出极高的偏移量，表明转录组弹性丧失。

3.5 群体间具有不同的ITS2类型图谱但群落组成相似

ITS2测序结果显示，两群体均以共生藻属（Symbiodinium）的A4/A4a类型为主导，但具体ITS2类型组合略有不同。离岸群体中的一个样本检测到少量Gerakladium（G）类型。总体而言，两群体的共生藻群落结构高度相似，不支持共生藻差异是导致热耐受性不同的主要原因。

总结与讨论

本研究证实，来自热波动较大环境的近岸星芒滨珊瑚群体，其热耐受性（以ED50衡量）显著高于来自热稳定环境的离岸群体，且这一差异在过去近十年中持续存在，表明其驱动因素可能是稳定的。转录组分析揭示，近岸群体在面对亚致死热胁迫时，表现出一种“平缓”的转录组响应，即DEGs数量少、变化幅度小，且能快速恢复到基线水平，符合“转录组弹性”和“转录组驯化”的特征。相反，离岸群体则表现出强烈且延迟的转录组响应，恢复也较慢。由于未检测到基因表达前置现象，且两群体共生藻组成相似，研究人员推测，近岸珊瑚更高效的胁迫管理能力可能源于表观遗传修饰（如DNA甲基化）或转录后调控（如选择性剪接）等机制，而非简单的基因表达水平改变。这些发现强调了宿主基因调控在塑造珊瑚热耐受性中的核心角色，并为利用CBASS结合高时间分辨率转录组学快速评估和筛选天然耐热珊瑚群体提供了有效的实验范式。