《Plants》：Integrated Transcriptomics and Metabolomics with Machine Learning Identify Flavonoids as Key Effectors in Wheat Root Thermotolerance
Wenyuan Shen,
Qingming Ren,
Yiyang Dai,
Yu Zhang and
Fei Xiong
编辑推荐：
为解决全球变暖背景下小麦根系耐热性分子机制不明的关键问题，研究人员结合表型组、转录组、代谢组学与机器学习算法，对比了耐热（YM158）与热敏感（YM15）小麦品种。研究发现，耐热品种通过激活内质网蛋白加工、类黄酮生物合成途径（特别是柚皮素等）及抗氧化酶系统，维持蛋白稳态与氧化还原平衡，其中类黄酮介导的抗氧化系统是关键的耐热效应机制。该研究为培育气候适应性小麦品种提供了重要理论基础与靶点。
随着全球变暖加剧，极端高温事件日益频繁，对全球粮食安全构成重大威胁。小麦作为全球约20%人口的主粮，其生产尤其易受高温胁迫的影响。高温会抑制光合作用、阻碍籽粒灌浆、引发氧化损伤并导致代谢紊乱，造成严重减产和品质下降。尽管对小麦地上部分的耐热性已有较多研究，但作为感知土壤温度变化、吸收水分养分“第一线”的根系，其耐热机制却长期被忽视。事实上，由于土壤的温度缓冲能力有限，根系往往比地上部分承受更长时间的热胁迫。因此，深入揭示小麦根系响应高温的分子调控网络，对于培育具有气候韧性的作物品种具有重要的战略意义。然而，调控小麦根系耐热性的核心分子机制与关键效应通路仍不明确。
为了回答上述问题，沈文源、任清明、代一阳、张宇和熊飞等研究人员在《Plants》期刊上发表了一项研究。他们选取了在田间表现出显著耐热性差异的两个小麦品种——热敏感型YM15和耐热型YM158，通过模拟不同梯度的温度胁迫，系统分析了其根系的表型、生理、转录组和代谢组响应。该研究整合了多组学数据与机器学习算法，旨在精确解析耐热品种维持根系结构与氧化还原稳态的特有优势，并鉴定决定耐热性的关键代谢通路与调控网络。
研究人员运用了多项关键技术方法。在实验设计上，他们对YM15和YM158小麦品种进行了水培，并设置了25°C、30°C和35°C的温度梯度处理。在表型与生理分析中，使用了根系扫描仪和植物根系分析系统评估根长、侧根数等形态指标，并采用ELISA试剂盒和分光光度法测定了丙二醛(MDA)、过氧化氢(H2O2)、根系活力以及超氧化物歧化酶(SOD)、过氧化物酶(POD)、过氧化氢酶(CAT)等抗氧化酶活性。在组学分析层面，利用非靶向代谢组学(LC-MS)和转录组测序(RNA-seq)技术分别获取了根系在高温下的代谢物与基因表达谱。数据分析则融合了先进的生物信息学工具：采用随机森林(Random Forest)算法从高维代谢数据中筛选关键生物标志物；运用加权基因共表达网络分析(WGCNA)构建基因-代谢物关联网络以识别核心调控模块；并通过基因本体(GO)和京都基因与基因组百科全书(KEGG)富集分析揭示差异基因的生物学功能与通路。最后，通过实时定量PCR(qRT-PCR)对关键枢纽基因的表达趋势进行了验证。
2.1. 小麦根系在高温胁迫下的表型和生理响应
研究发现，高温（35°C）严重抑制了小麦生长。与25°C对照相比，热敏感品种YM15的总根长下降了85.7%，而耐热品种YM158仅下降68.3%，表现出更强的生长维持能力。在生理层面，极端高温下YM15的H2O2和MDA含量急剧上升，表明其细胞膜系统遭受严重氧化损伤。相反，YM158的增幅显著更低，且能维持较高的根系活力。YM158的SOD、POD、CAT等抗氧化酶活性在高温下显著升高，而YM15的酶活性则被抑制甚至下降。这表明YM158通过维持高效的抗氧化酶系统来缓解氧化损伤，是其耐热性的重要生理基础。
2.2. 代谢组学分析
代谢组分析显示，高温诱导了小麦根系显著的代谢重编程。随着胁迫加剧，两个品种的差异积累代谢物(DAMs)数量均增加。在35°C极端高温下，YM158的上调代谢物数量(1208个)显著多于YM15(997个)，表明其代谢动员能力更强。为了从海量数据中精准定位决定耐热性差异的关键代谢物，研究引入了随机森林特征选择算法。通过分析特征重要性，最终筛选出9个关键生物标志物，主要涉及激素调节、次生代谢防御和信号转导三类。其中包括细胞分裂素前体核糖基玉米素磷酸、生长素结合物IAA-Glu、类黄酮关键中间体柚皮素、没食子单宁ginnalin B、木脂素balanophonin以及氧化脂质信号分子9-羟基-10-氧代-12Z-十八碳烯酸等。
2.3. 转录组分析
转录组测序表明，温度是引起转录组变异的主要因素。在30°C中度胁迫下，YM158的差异表达基因(DEGs)数量(4022个)少于YM15(5568个)，显示出更强的转录稳态维持能力。而在35°C下，两者均发生全基因组范围的重编程。功能富集分析揭示了两品种响应策略的差异。在30°C时，YM158特异性地富集于乙烯激活信号通路、肉桂酸生物合成过程、MAPK信号通路和植物激素信号转导等，表明其能早期启动信号传导与次生代谢防御。而YM15则主要富集于谷胱甘肽代谢等简单的抗氧化反应。在35°C时，YM158特异性地显著激活了内质网蛋白质加工通路，而YM15则没有，提示蛋白稳态维持对极端高温下的细胞存活至关重要。
2.4. 加权基因共表达网络分析(WGCNA)
为构建系统性的调控图谱，研究整合转录组与9个关键代谢标志物进行了WGCNA分析，将差异表达基因聚类为14个共表达模块。其中，深绿色模块和浅黄色模块与YM158中特异积累的关键代谢物关联最强。深绿色模块功能富集于类黄酮生物合成、苯丙烷生物合成、过氧化物酶活性等，其蛋白互作网络(PPI)的核心枢纽基因包括查尔酮合成酶2(CHS)、柚皮素2-氧戊二酸3-双加氧酶以及阳离子过氧化物酶SPC4等，该模块与柚皮素积累呈显著负相关。浅黄色模块则富集于MAPK信号通路、植物激素信号转导以及多种激素响应过程，其PPI网络中的枢纽基因包括热激转录因子HSF B-1、HSF C-1b、受体样蛋白激酶FER等，该模块与激素相关代谢物（如核糖基玉米素磷酸）显著相关。
2.5. 根系响应高温胁迫的重要通路
基于转录组KEGG富集、代谢组机器学习和WGCNA三方面的证据均汇聚于同一通路，研究确认类黄酮生物合成是YM158抵御高温胁迫的核心防御机制。对通路中关键酶基因表达模式的分析显示，在YM158中，查尔酮合成酶(CHS)、查尔酮异构酶(CHI)等基因在高温下（特别是30°C）被迅速且显著诱导表达，这与代谢组中观测到的柚皮素等类黄酮物质的大量积累在时相和数量上高度吻合。
2.6. 实时定量PCR验证
通过qRT-PCR对筛选出的6个关键候选基因（包括CHS2、CHI、F3H、ANS、HSF B-1、HSF C-1b）的表达趋势进行验证，结果与RNA-seq数据高度一致，证实了转录组数据的可靠性。
研究结论与讨论
本研究通过多组学与机器学习整合分析，揭示了一个由酶促与非酶促抗氧化防御、蛋白稳态维持以及激素-次生代谢信号互作构成的多维耐热调控网络。研究人员据此提出了一个工作模型来总结这些复杂机制。
核心结论包括：
- 1.
类黄酮的核心作用：类黄酮生物合成途径，特别是其中间产物柚皮素，是耐热品种YM158应对高温的关键效应通路。类黄酮不仅作为强效的非酶促抗氧化剂，还能与抗氧化酶系统(SOD, POD, CAT)协同作用，高效清除过量活性氧(ROS)，恢复氧化还原稳态。
- 2.
信号感知与转导：耐热品种可能通过细胞壁完整性传感器（如受体激酶FER）感知高温引起的细胞壁流变学变化，进而激活下游复杂的激素信号网络（涉及ABA、乙烯、生长素等）。这些激素信号与类黄酮合成之间存在精细的反馈调节环路，例如ABA可促进类黄酮合成，而类黄酮又能调节ABA信号和生长素的运输与稳态，从而在防御与生长间取得平衡。
- 3.
蛋白稳态维持：在极端高温(35°C)下，YM158特异性地激活了内质网蛋白质加工通路，通过上调热激转录因子(HSFs)和分子伴侣来修复或清除错误折叠的蛋白，这一蛋白稳态维持机制对细胞在致命高温下的存活至关重要，而热敏感品种YM15缺乏有效的蛋白修复响应。
- 4.
方法学创新：研究引入随机森林等机器学习算法，成功从数千个代谢特征中筛选出具有明确生物学意义的9个关键生物标志物，克服了传统差异分析仅关注差异倍数的局限，为解析复杂农艺性状的机制提供了有力的数据驱动策略。
这项研究的意义在于，它不仅首次通过整合多组学与机器学习，系统描绘了小麦根系适应高温的分子图谱，还将类黄酮介导的抗氧化系统定位为培育气候适应性作物的潜在靶点。所鉴定的关键代谢标志物（如特定类黄酮）和核心枢纽基因，为耐热小麦种质资源的快速筛选与分子设计育种提供了宝贵的理论依据和遗传资源。尽管仍需通过基因编辑等功能实验进行最终因果验证，但该研究无疑深化了人们对作物根系耐热机制的理解，并为应对全球变暖下的粮食安全挑战提供了新的思路与工具。
