近未来气候变化下变温动物胚胎对预测热胁迫的调控响应

《Molecular Ecology》：Regulatory Responses of Ectothermic Embryos to Predicted Heat Stresses Under Near-Future Climate Change

　　

热胁迫会破坏支撑正常发育的协调调控过程，导致胚胎畸形与不可存活。因此，胚胎对热挑战的缓冲作用在个体生存与种群韧性应对热胁迫中可能起关键作用。然而，支撑胚胎抵御具有生态意义热胁迫的分子机制仍知之甚少。研究人员以三种在近未来气候变暖下将会增加的热胁迫模式，探究了褐安乐蜥（Anolis sagrei）胚胎的调控响应。在慢性、持续性与急性热处理 across 中，胚胎表现出快速的转录组广泛扰动，随后在24 h内恢复至正常发育轨迹。尽管存在这种广泛扰动，一个 governing 前脑与颅面发育的核心调控模块在各热胁迫中仍保持强烈保守。热胁迫诱导了广泛的转录抑制，同时转录后调控（post-transcriptional regulation）依赖性增加，以及RNA通量（RNA flux）呈剂量依赖性减慢，这与蛋白质毒性应激（proteotoxic stress）的缓解一致。与此同时，胚胎通过招募模块外应激响应基因与加强发育通路内共表达，重构调控网络。这些结果共同揭示了一种层级调控策略，胚胎在其中缓冲必需发育程序，同时灵活重构外周网络以管理细胞应激。这种系统水平的韧性强调了早期生命适应性调控可塑性在塑造生物体响应气候变暖中的重要性。

研究背景方面，环境温度是普遍存在的压力因子，约束着从生长、表现到适合度各生物学层级的极限。全球变暖、城市热岛效应及热浪频率与强度增加，正对野生动物施加日益增加的生理胁迫。所有生物依赖动态体温调节机制缓冲热应激，维持正常生长、生存与繁殖所需的稳态，而这种适应性可塑性的有效性取决于热剂量、进化史与生物生活史阶段。胚胎发育期间热是一种尤其 potent 的致畸原，热胁迫可诱导细胞死亡并破坏支撑正常发育的关键信号通路，导致后期生活中众多不适应的生物学后果。仅超出正常形态发生所需范围2°C–4°C的温度，就可在包括哺乳类、鸟类、鱼类与爬行类的脊椎动物中导致颅面畸形、神经缺陷、脑容量减少、行为异常与存活率下降。现有物种对气候敏感性模型偏重于理解成体对慢性热挑战的响应，而早期发育期间的体温调节损伤可能显著影响面临近未来气候变化的野生动物的健康、存活力与进化。卵生物种产卵后很少或无亲代照料，胚胎热孵化环境由母体巢址选择决定，发育中胚胎需应对日常温度波动、可变平均温度、长期热胁迫与急性热激暴露以在现代生态系统生存，其行为缓冲各类热胁迫能力有限，适应性调控可塑性可能在卵绑胚胎生命维持热损伤期间正常颅面发育模式中起关键作用。尽管特定有序转录程序是正确发育图式形成所需，但这些程序外基因表达更具变异性并可响应环境线索差异调制，细胞热休克响应（heat-shock response）是进化上高度保守的保护过程以减轻热挑战下蛋白质错误折叠与聚集诱导的细胞应激，为启动这种关键补偿响应，细胞须显著重组其调控程序，迅速上调分子伴侣 suite 协助蛋白质折叠并帮助清除受损大分子，同时下调数千个调控细胞周期、转录与代谢的基因，若该热休克响应不足以保护细胞免受降解蛋白、受损DNA与染色体畸变积累，细胞进入程序性细胞死亡。胚胎细胞须在维持正确发育轨迹与防御热诱导稳态威胁间权衡，尽管特定分子伴侣如热休克蛋白（heat shock proteins, HSPs）在许多物种中被充分研究，但对胚胎热应激响应的更广义系统生物学、发育细胞如何校准对不同幅度与暴露时间热胁迫的投入，以及适应性调控可塑性在近自然环境经历的热挑战下保护维持正确形态发生所需敏感细胞信号通路的程度仍知之甚少。理解热胁迫下面维持正常发育模式的胚胎机制，将为热韧性与气候改变下不易在后期成体阶段观察到的热诱导选择潜在靶标提供关键信息。为此，研究人员评估了古巴褐安乐蜥（Anolis sagrei）前脑与面部胚胎发育期间适应性调控可塑性的模式，其为广泛的温热广义者与新兴的脊椎动物热应激发育模型系统。

为开展研究，研究人员主要运用了以下关键技术方法：首先基于公开历史气候数据与未来降尺度气候投影，在中等变暖情景下量化A. sagrei分布范围内超过安全筑巢温度的最大日均环境温度比例，并用比例相等性检验比较近未来与历史时期差异；其次使用生物物理模型NicheMapper整合环境影响辐射、对流与传导的变量，以2.5°空间与月度时间分辨率估算未来百年巢址微气候条件，并用广义线性混合效应模型检验未来巢址超过临界热敏感窗口的小时数变化；动物方面收集佛罗里达州珊瑚阁市入侵种群妊娠雌蜥所产卵，随机分配至对照、慢性轻度胁迫、持续中度胁迫与急性热激四种孵化处理，在热挑战后立即与24小时后采集分期匹配的胚胎前脑与面部组织；提取总RNA并制备mRNA文库进行测序，原始表达数据修剪后与卡洛琳安乐蜥（Anolis carolinensis）参考基因组比对，过滤低表达基因得到终末转录组；随后依次运用主成分分析（principal component analysis, PCA）评估整体基因表达模式，用edgeR包估计对照与热处理间对数倍差异表达并做功能富集；通过外显子-内含子分裂分析（exon-intron split analysis, EISA）评估转录后调控贡献，量化成熟与未成熟RNA读数变化；采用加权基因共表达网络分析（weighted gene correlation network analysis, WGCNA）识别端脑转录组发育进展相关的从头共表达网络并检验模块保存性；运用个体水平乘积相关（correlation by individual level product, CILP）评估热应激对调控连贯性（regulatory coherence）的影响，解析应激响应、潜在恢复与整体发育轨迹的调控影响并分类基因-基因连接。

研究结果部分，3.1 近未来气候变化情景下的胚胎应激（Embryonic Stress Under Near-Future Climate Change Scenarios）：研究人员首先评估A. sagrei巢址在未来百年可能变化的热胁迫情景，发现到2100年A. sagrei胚胎在发育期间可能面临显著更具挑战性的热条件，其原生范围日均最高环境温度与模拟巢址温度较历史条件显著增加，过去十年仅18%原生范围年均日最高温超过安全孵化温度33°C，而在SSP370中等变暖情景下到本世纪末超85%原生范围将经历危及胚胎发育的环境温度，平均巢址每日超过33°C的小时数较近年历史条件增加超3小时，典型巢址在2081–2100年平均每日有10小时超33°C、7小时超36°C、3小时超40°C，既往实验显示此类热挑战使胚胎存活率降低超10%、发育畸形增加四倍，因此世纪末气候可能大幅降低A. sagrei在其原生范围大部份的繁殖成功与幼体适合度。

3.2 不同热胁迫情景下发散性热诱导单基因调控模式（Patterns of Divergent Heat-Induced Single Gene Regulation Under Diverse Heat Stress Scenarios）：研究人员发现热胁迫暴露导致形态发生早期胚胎端脑组织全转录组调控扰动及剂量依赖性向正常表达的恢复，热挑战后立即发育中前脑与面部组织表达显著偏离对照孵化条件模式，持久中度热胁迫与急性热激暴露胚胎总体表达在对照孵化温度恢复24小时后大体恢复，慢性轻度胁迫条件下孵化第二日热应激胚胎基因表达模式继续显著偏离对照水平。单基因水平直接比较对照与热处理表达模式发现，热挑战后立即（t₀）有76个基因跨处理普遍上调、25个普遍下调，24小时后（t₁）较少基因差异表达，31个普遍上调、5个普遍下调，这些共享差异表达基因关联细胞在热胁迫下补偿损伤的核心机制，t₀上调关联微管运动、纤毛组装与蛋白结合以维持对端脑分化关键的Hedgehog信号通路，t₀转录相关基因下调提示热应激下快速转录抑制以缓解蛋白质毒性应激，t₁热应激胚胎脑与面部细胞显示转录后调控相关基因上调与细胞应激管理相关基因与通路下调，提示生理恢复与回归稳态模式。

3.3 胚胎热胁迫的转录后响应（Posttranscriptional Response to Embryonic Heat Stress）：研究人员通过外显子-内含子分裂分析比较成熟mRNA与premRNA产物生产的时间依赖性变化，发现热应激下RNA通量普遍减慢，热暴露持续时间越长减慢越明显，慢性轻度胁迫相对对照孵化温度RNA通量降低最大，急性热激最小，若累积热暴露对RNA代谢约束强于短时极端温度，慢性变暖可能对发育稳态构成比急性热损伤更大挑战。个体基因水平发现持久中度胁迫与热激处理较对照与慢性孵化条件更易出现RNA稳态的转录后贡献，对照条件下4个关联早期脑发育与调控RNA通量的基因受转录后调控，其中mid1的转录后调控在持久中度胁迫与热激下保留，这些热处理关联转录后调控基因数增加至20与17个，关联脑发育、RNA通量调控与细胞应激响应，热胁迫下稳态RNA水平的富集转录后变化与跨处理剂量依赖性RNA通量减慢共同表明转录后机制在缓冲中短时程热挑战对发育必需基因表达影响中起重要作用。

3.4 跨转录组热诱导调控架构改变（Heat-Induced Alterations of Regulatory Architecture Across the Transcriptome）：研究人员检验此发育窗口基因调控网络如何在热应激下扰动，发现对照孵化条件48小时实验期内单一由1087个基因组成的调控模块显著关联神经发育进展（发育模块），富集胶质细胞分化生物学过程与前脑及颅面形态发生调控的先验候选基因且含sox9、dlx2、bmp7等模块枢纽基因，该发育模块在各热处理中强烈保存（Z_summary>10，medianRank<3），而其余转录组广泛调控重构，30个整体架构模块中16–20个在热处理下无保存证据、7–11个模块分散度大于偶然预期，热处理中出现富集热休克蛋白、蛋白质折叠与细胞代谢的新共表达模块。发育模块内基因在热应激下加强模块内共调控动态并形成新的显著模块外共调控互作，热应激后立即所有处理中发育模块内模块内连贯性（intramodular coherence）较对照条件更显著，次日通常弱于对照胚胎，对照条件下此窗口显著改变连贯性的基因对中大部份在热处理中维持或恢复连贯性模式，维持关系共同富集形态发生、神经发育与有丝分裂相关位点，而延迟或破坏的模块内基因对关系跨热剂量无共性。模块内更强共表达与显著基因-基因关系维持由招募模块外应激响应基因支持，热应激立即使发育模块与模块外基因的成对关系平均强度较理想孵化条件增加超两倍，两日后慢性轻度与热激处理仍显著更强，热处理两日增加发育模块与转录组其余基因的新关系，跨处理t₀常见获得关系由关联细胞过程负调控、细胞应激响应与蛋白结合的模块外基因组成，t₁常见获得关系关联转录后调控、发育与生长，表明热应激下通过补偿性招募额外基因维持协调发育程序以强力防御脑与面部形态发生所需转录程序，模块内协调维持由与负调控转录及转录后调控相关更强模块外关系支持，后者是热挑战的快速响应。

总结讨论部分，研究人员结论为胚胎对热胁迫的韧性源于多生物学组织层面的协调调控，热应激触发层级响应而非均匀破坏发育：胚胎保留形态发生必需的核心调控模块，同时动态抑制转录、调节RNA通量并招募辅助应激响应通路，这种稳健性与灵活性组合使胚胎组织在否则会损害存活力的条件下维持发育通路。这些发现推进了理解环境应激下发育耐受的系统水平框架，随着气候变化加剧热变异与极端，脆弱生活史阶段的整合调控机制可能在物种存续中起核心作用，确定这些调控策略跨分类群、环境与进化史的变异，是预测全球变化下自然种群适应能力与轨迹的关键下一步。未来工作应确定这些模式是否在物种范围内具普遍性及是否反映维持形态发生的功能机制，必要下一步包括：在A. sagrei入侵与原生范围及别的卵生变温动物物种复制实验以理解观察模式普适性及测试胚胎热应激的调控响应与适合度成本在范围扩张与定殖中是否分化；对关键发育调控因子与新招募应激响应基因做蛋白水平验证以确立此处观察的转录协调是否转化为蛋白丰度与活性的协调变化；通过CRISPR介导基因编辑或药理抑制蛋白稳态与RNA加工通路靶向候选枢纽基因或应激响应通路的扰动实验，检验招募这些基因是否为热挑战下保留前脑与颅面发育所必需，此类实验将在网络重构、发育稳健性与变暖环境中适应性调控可塑性间提供关键机制联系。论文发表在《Molecular Ecology》。