《Food and Energy Security》：Time-Resolved Transcriptomics Reveal Core and Genotype-Specific Pathways Underlying Thermotolerance in Wheat
编辑推荐：
这篇研究通过评估30个小麦品种的幼苗耐热性，并选择耐热基因型Huaimai211和热敏感基因型Chuke316进行五个时间点的转录组测序，系统解析了小麦应对高温胁迫的分子机制。研究揭示了涉及蛋白质稳态维持（如HSP90结合、未折叠蛋白反应）和光合作用持续抑制的核心响应通路，并发现了与基因型特异耐热性相关的代谢可塑性（如氨基酸生物合成）。关键枢纽基因TaRbcS-2B.3的鉴定为理解光合调节与胁迫适应的联系提供了新见解。这些发现为通过分子标记辅助选择和基因编辑培育耐热小麦品种提供了重要靶点。
引言
全球气候变化及极端高温事件频发对作物生产构成严重威胁。小麦作为全球最重要的主粮作物之一，在开花和灌浆期对高温尤为敏感。因此，鉴定耐热基因型并理解其分子机制，对于培育能够稳定生产的高温适应型品种至关重要。以往研究虽然鉴定出许多与热胁迫相关的基因，但大多局限于单个时间点或基因型，缺乏动态、整合的网络分析。本研究旨在通过时间分辨的转录组学，系统阐明小麦幼苗期耐热性的核心与基因型特异调控框架。
材料与方法
2.1 植物材料、生长条件与热胁迫处理
本研究使用了30个主要源自中国湖北麦区的小麦品种。种子在滤纸上萌发后，移栽至培养盘中，在可控生长室内培养8天后，置于40°C下进行48小时热胁迫处理。胁迫结束后恢复48小时，记录存活率、株高和地上部鲜重。基于表型差异，选择耐热基因型Huaimai211和热敏感基因型Chuke316用于转录组分析，分别采集0、6、12、24和48小时热胁迫后的叶片样品。
2.2 RNA提取与文库构建
使用商业试剂盒提取总RNA，并通过NanoDrop和Bioanalyzer评估其纯度与完整性。使用Illumina建库试剂盒构建文库，经质量检测合格后进行测序。
2.3 转录组测序与读数比对
在Illumina平台上进行双端150bp测序。获得的高质量clean reads使用HISAT2比对至中国春参考基因组，并使用StringTie进行转录本组装。基因表达水平通过FPKM进行标准化，用于后续的主成分分析和层次聚类。
2.4 差异表达基因分析
使用DESeq2软件包识别差异表达基因。将|log2FC| ≥ 1且错误发现率FDR < 0.01的基因定义为DEGs。分析包括同一基因型内不同时间点间的比较，以及不同基因型在同一时间点间的比较。
2.5 GO与KEGG分析
使用GOseq R包进行基因本体富集分析，使用KOBAS软件进行KEGG通路富集分析。校正后的p值<0.05被认为显著富集。
2.6 加权基因共表达网络分析
使用WGCNA R包对所有DEGs进行共表达网络分析，以识别与热响应性相关的基因模块。通过模块特征基因与表型性状（耐热性）的相关性鉴定关键模块。
2.7 亚细胞定位
将枢纽基因TaRbcS-2B.3的编码序列克隆至带有GFP标签的表达载体，通过聚乙二醇转化法瞬时转化小麦叶肉原生质体，利用共聚焦显微镜观察GFP荧光信号，确定其亚细胞定位。
2.8 选定DEGs的表达验证
通过qRT-PCR对6个选定的DEGs进行表达验证，以内参基因Actin进行标准化，使用2?ΔΔCT法计算相对表达量。
结果
3.1 30个小麦种质资源间耐热性变异
对30个小麦品种的幼苗期耐热性评估显示，其存活率、株高和鲜重降低存在显著的自然变异。其中，Huaimai211表现出最高的存活率和最小的生长抑制，而Chuke316最为敏感。相关性分析表明，存活率与株高、鲜重降低呈显著负相关。
3.2 两个对比基因型在热胁迫下的转录组分析
对HR和HS基因型在五个时间点进行转录组测序，获得了高质量数据。主成分分析显示，样本按热胁迫时间和基因型清晰聚类，PC1代表时间进程，PC2区分基因型。样本间相关性热图证实了重复样本间的高度一致性。
3.3 共有与基因型特异DEGs的鉴定
在HR基因型中，DEGs数量在48小时达到峰值。相比之下，HS基因型在所有时间点的DEGs数量都更少。共鉴定出8517个在HR中持续热响应的基因，以及4155个在HS中持续响应的基因，其中3028个为两基因型共享的核心DEGs。此外，HR在多个时间点拥有236个自身特异的持续响应基因。
3.4 核心、时序及基因型特异热响应DEGs的富集分析
对3028个核心DEGs的GO分析显示，上调基因显著富集于“Hsp90蛋白结合”、“未折叠蛋白结合”等蛋白折叠与胁迫响应通路，而下调基因富集于“核酮糖-1,5-二磷酸羧化酶活性”、“叶绿素结合”等光合作用相关通路。KEGG分析进一步证实，上调基因涉及“内质网蛋白质加工”、“谷胱甘肽代谢”，下调基因涉及“光合作用-天线蛋白”等通路。时序分析揭示了从早期胁迫感应（如离子结合）到后期适应机制（如GTP酶活性）的转变。对HR特异DEGs的分析发现，其显著富集于氨基酸代谢、能量调控和蛋白质折叠相关通路。
3.5 共表达簇的鉴定
对核心DEGs的共表达分析将其分为四个上调簇和四个下调簇。上调簇（如C1、C5、C7）主要参与蛋白折叠和胁迫响应，下调簇（如C2、C4、C8）则与光合作用等代谢过程相关，表现出随着胁迫时间延长而持续下调的模式。
3.6 共表达网络中热响应模块的鉴定
WGCNA分析鉴定出15个共表达模块，其中firebrick3模块与热响应性的相关性最强。该模块中的基因模块内连接度与基因显著性高度相关。网络可视化揭示TaRbcS-2B.3为一个关键的枢纽基因，与众多共表达的DEGs紧密相连，这些基因主要富集于光合作用和碳固定相关功能。亚细胞定位实验证实TaRbcS-2B.3定位于叶绿体。
3.7 以光合作用和蛋白质稳态转换为中心的热响应模型
基于一致的时序表达模式，筛选出89个核心DEGs，其表达热图清晰区分了对照与热胁迫样本。qRT-PCR验证了包括TaHSP90、TaGST在内的6个基因的表达模式与RNA-seq结果高度一致。枢纽基因TaRbcS-2B.3属于RbcS家族，是卡尔文循环的关键酶，其表达量在耐热品种中显著高于热敏感品种。整合以上结果，研究提出了一个小麦幼苗热响应的工作模型：热胁迫触发光合作用组分（如CABs、Rubisco/RbcS、天线蛋白）的持续抑制，同时激活HSF-HSP分子伴侣网络以维持蛋白质稳态。这种协调的转换将资源从碳同化重新分配到蛋白质保护，可能构成了HR和HS基因型间耐热性差异的基础。
讨论
本研究整合了多品种表型评价和双基因型时序转录组分析，全面揭示了小麦幼苗在长期热胁迫下的细胞功能重组。共享的3028个核心热响应基因及其富集的通路，与其他物种中报道的热适应策略一致，构成了维持蛋白质稳态和限制氧化损伤的早期保护层。然而，本研究发现光合作用相关基因在所有胁迫时长中均表现出持续下调，这不同于温和或短期处理下的短暂下降，暗示了一种从生长导向碳同化向保护功能进行资源战略性重分配的适应策略。
基因型特异分析阐明了HR基因型具有更优越热性能的原因。HR中更多的DEGs反映了其更强的转录可塑性与调控能力，而非胁迫损伤，这由其能量代谢、蛋白质稳态和抗氧化通路的协调富集所支持。时序分析揭示了从早期分子感应到后期调控过程的阶段性适应模式。共表达和网络分析支持对热响应的系统水平理解，将核心DEGs解析为功能一致的表达簇，表明小麦幼苗采用了一种有组织的调控策略。firebrick3模块与热响应性的强关联，及其在光合和代谢通路中的一致性富集，强调了模块化基因调控在塑造基因型差异中的重要性。该模块中的枢纽基因TaRbcS-2B.3，可能有助于协调光合抑制与胁迫适应过程。
尽管存在研究局限，但本整合性方法揭示了小麦热响应的关键特征，包括光合作用的持续抑制、蛋白质稳态的强化以及基因型特异的代谢可塑性。这些发现，特别是firebrick3模块和RbcS家族内的基因，为分子标记辅助选择、基因编辑和基因组预测提供了有潜力的靶标。增强早期耐热性，是稳定小麦在日益频繁的极端高温事件下生产的一条重要途径。
