随着全球气候变暖加剧，预测生物如何通过进化来适应不断变化的环境，成为生态与进化生物学领域的核心挑战。在这个过程中，表型可塑性——即同一基因型在不同环境下产生不同表型的能力——被认为是生物快速响应环境变化的关键策略。然而，我们对表型可塑性背后的遗传和分子机制仍知之甚少，这使得预测其在气候变化下的进化模式困难重重。海洋无脊椎动物，尤其是生活在潮间带、经受剧烈温度波动的牡蛎，是研究热适应与可塑性的理想模型。在近缘的牡蛎物种中，适应中国北方冷水的长牡蛎（Crassostrea gigas）和适应南方暖水的葡萄牙牡蛎（Crassostrea angulata），为探索温度适应如何塑造基因表达模式的分歧提供了绝佳的对比体系。先前研究表明，热休克蛋白90（Hsp90）作为关键的分子伴侣，在热胁迫响应中扮演核心角色，并且在C. angulata的上游非编码区显示出强烈的选择性清除信号，提示其可能是温度适应与可塑性进化的关键靶点。为了深入解析表型可塑性分歧的遗传基础，一项发表于《Marine Life Science》的研究，以Hsp90为研究对象，开展了一项整合功能验证与遗传机制探索的综合性工作。
为开展此项研究，研究人员运用了多种关键技术方法。在功能验证层面，通过RNA干扰（RNAi）在活体牡蛎中敲低Hsp90表达，并结合热激实验评估其对耐热性的影响；同时，在人类细胞系（HEK293T和HeLa细胞）中进行Hsp90的过表达实验，通过流式细胞术和共聚焦显微镜分析其在热胁迫下的抗凋亡作用与亚细胞定位。在遗传机制探索方面，研究结合了正向与反向遗传学策略：首先通过对两个物种野生个体启动子区进行测序，筛选出显著差异的SNP和Indel位点；然后利用双荧光素酶报告基因检测系统，在细胞水平系统评估了这些位点突变对启动子活性的影响，特别是在热胁迫（42°C）下的反应，以鉴定具有基因型（G）和基因型-环境互作（G×E）效应的位点。此外，研究利用了一个已有的C. angulata (♀) × C. gigas (♂) 杂交F2代群体，在活体水平验证了这些位点与Hsp90表达模式的关联。为了阐明变异位点的作用机制，研究人员进行了DNA pull-down结合质谱分析，以寻找差异结合的反式作用因子，并通过电泳迁移率变动分析（EMSA）和报告基因共转染实验，验证了特定转录因子与变异位点的特异性结合及功能调控。样本主要来源于两个物种的自然栖息地（中国北方和南方海域）以及实验室构建的杂交F2群体。
研究结果部分通过一系列实验逐步揭示了Hsp90的功能及其表达可塑性分歧的机制。
1. 功能性表征证实牡蛎Hsp90是耐热性的理想代理指标
研究人员首先在细胞和个体水平验证了Hsp90在热适应中的核心作用。在HeLa细胞中，热胁迫导致带有mCherry标签的CgHsp90蛋白的红色荧光信号显著增强。在HEK293T细胞中，过表达Hsp90能显著降低热胁迫（42°C，1小时）诱导的细胞凋亡比例。相反，在活体牡蛎中，通过RNAi敲低CgHsp90表达后，个体在42°C热激下的存活率显著下降。这些结果直接证明了Hsp90对于增强牡蛎耐热性至关重要。进一步的免疫共沉淀结合质谱分析鉴定出491个与Hsp90相互作用的蛋白质，其中295个为特异性结合。基因本体（GO）和京都基因与基因组百科全书（KEGG）富集分析显示，这些互作蛋白主要参与能量代谢（如氧化磷酸化、TCA循环）、蛋白质翻译和细胞骨架调控等通路，提示Hsp90可能通过协调能量供应和蛋白质稳态来增强耐热性。
2. 筛选具有潜在致病作用的遗传变异
比较两个物种的Hsp90表达模式发现，C. angulata的Hsp90基础转录水平显著高于C. gigas（6.3倍）。在热胁迫下（37°C），C. gigas的Hsp90表达在12小时达到峰值（16.2倍），然后迅速下降；而C. angulata的表达上调缓慢，在24小时才达到峰值（2.6倍）。这种模式分别与它们栖息地温度的波动性（北方）和均值（南方）相适应。对Hsp90启动子区的测序分析，在两个物种间鉴定出六个SNP和三个Indel存在显著差异。
3. 鉴定所检测遗传变异上的G×E效应
通过双荧光素酶报告基因实验，研究人员发现四个变异位点（三个SNP: g.-2291 T>G, g.-1736A>T, g.-821C>T 和一个Indel: g.-1511ATT>A）在非胁迫条件下能显著降低转录活性，其中g.-2291 T>G和g.-1736A>T在热胁迫下还显著增加了转录幅度。双因素方差分析（ANOVA）确认了这四个位点具有G和G×E效应。在F2杂交群体中的验证进一步支持了这些位点与Hsp90表达模式的关联，其中C. angulata纯合等位基因在24小时热胁迫下维持高表达，而C. gigas纯合等位基因则在6或12小时达到表达峰值，表现出更敏感的模式。
4. 顺式变异中G×E效应的功能分析
为阐明变异位点的作用机制，研究针对g.-2291位点进行了深入探究。DNA pull-down结合质谱分析发现，嘌呤丰富元件结合蛋白B（Purine-rich element binding protein B, PURB）能特异性结合C. angulata的等位基因（g.-2291G）探针，而不结合C. gigas的等位基因（g.-2291 T）探针。EMSA实验证实了这种结合的特异性。报告基因实验显示，过表达PURB能显著提升含有g.-2291G等位基因的启动子活性，但对含有g.-2291 T等位基因的启动子无影响。此外，两个物种的PURB基因在热胁迫下的表达趋势相反：在C. angulata中，PURB表达先显著下降后上升；而在C. gigas中则呈现相反模式。这解释了C. angulata中Hsp90表达可塑性较低的现象。
在结论与讨论部分，研究系统地总结了发现并阐述了其重要意义。本研究发现，Hsp90的基础（组成型）表达和热诱导可塑性分别与适应栖息地温度的均值和波动性相关联。适应温暖环境的C. angulata表现出更高的基础表达和更低的可塑性，这是一种优化能量配置以适应稳定高温环境的策略；而适应寒冷但温度多变的北方环境的C. gigas则表现出更高的表达可塑性，以应对剧烈的温度波动。这种表达模式的分歧受到Hsp90启动子区紧密聚集的四个顺式变异位点的调控，这些位点被证实具有G+G×E效应。其中，位于g.-2291的关键SNP（T>G）通过创建（G等位基因）或删除（T等位基因）转录因子PURB的结合位点，直接调控了Hsp90的基础表达水平。同时，PURB本身对热胁迫的响应（反式变异）则决定了Hsp90可塑性表达的大小。因此，可塑性表达模式的分歧，是由热胁迫诱导的顺式变异与反式变异之间的相互作用共同调控的。
这项研究的意义深远。首先，它在非模式生物中建立了一个研究范式，成功地将正向遗传学（群体变异筛选与关联分析）与反向遗传学（细胞水平的功能验证与机制解析）相结合，将G×E效应的遗传基础解析精度提升至单核苷酸水平。这为在其他生物中开展类似研究提供了可借鉴的技术路线。其次，研究揭示了表型可塑性进化的精细遗传调控机制，表明环境选择的热点（顺式变异）可以通过与反式作用因子的互作来共同塑造生物对环境变化的反应（G×E效应），深化了我们对表型可塑性在适应性响应中重要作用的理解。最后，这项研究具有明确的应用潜力。所鉴定的关键遗传变异（如g.-2291G）为通过分子标记辅助选择或基因编辑技术，培育具有增强耐热性的牡蛎新品种提供了直接靶点，有助于应对全球变暖导致的夏季大规模死亡现象，保障牡蛎养殖产业的可持续发展。总之，该工作通过一个生动的案例，阐明了温度变化（均值和波动）如何塑造可塑性反应（G×E）的改变，为预测海洋生物在持续气候变化下的适应潜力提供了新的见解。
