整合生理化学、转录组学与代谢组学揭示类黄酮（FLVs）生物合成途径在紫大麦干旱响应中的潜在作用

《Plant Stress》：Integrated physio-chemical, transcriptomics and metabolomics analyses reveal the potential roles of flavonoids biosynthesis pathway in drought response of barley

干旱胁迫显著降低大麦产量和品质。类黄酮（FLVs）是与作物干旱适应相关的重要抗氧化代谢物。然而，与大麦FLVs介导的干旱响应相关的分子机制仍未完全理解，需进一步功能研究。本类研究导向综述整合生理学、转录组学、代谢组学和计算证据，将紫大麦中FLVs相关干旱响应的当前知识转化为未来研究方向。实验证据表明，干旱胁迫增加FLVs积累、抗氧化酶活性和渗透调节物质含量，同时优化与有效保水相关的气孔参数。转录组分析鉴定出826个上调基因，富集于苯丙烷类和FLVs生物合成途径。代谢组分析检测到3,335个差异积累代谢物，FLVs合成途径显著富集。整合多组学分析提示，干旱胁迫下参与FLVs生物合成的酶（CHS、CHI、F3H、FLS、DFR）表达与代表性FLVs代谢物发生协同变化。二十个高度互作蛋白由FLVs生物合成途径的十个酶和十个转录因子（ERF、MYB、bZIP、WRKY、NAC）组成，在干旱胁迫期间显著上调，其启动子区域含有显著的干旱响应顺式元件。与关键FLVs代谢物的分子对接也预测出良好结合亲和力，并鉴定出候选蛋白质–配体相互作用以供未来功能验证。总体而言，这些发现为FLVs介导的大麦干旱响应分子机制提供了有价值见解，可指导未来开发耐旱大麦品种的研究。

大麦（Hordeum vulgare L.）是全球第四大谷物，广泛用于食物、饲料、制麦和酿造。紫大麦积累较高水平的花青素和类黄酮（flavonoids，FLVs），并表现出对包括干旱在内的非生物胁迫较强耐受性。干旱严重限制全球作物生产力，严重条件下可造成作物减产达50%。尽管大麦相对许多谷类作物更耐旱，但FLVs介导的干旱适应分子基础仍不完全清楚；已有研究识别了生理、生化和遗传响应，但如何将生理响应、转录组重编程和代谢物积累整合起来仍有限。Morex大麦基因组已测序，使调控大麦干旱信号通路的分子基础研究可行。研究人员开展的这项研究导向综述整合了当前关于FLVs介导大麦耐旱性的知识，并辅以生理、转录组、代谢组和计算分析，以支持生物学解释并提供干旱适应分子机制的整合视角。大麦中FLVs主要包括黄烷醇、花青素和原花青素，黄烷醇和花青素多定位于籽粒果皮；蓝粒和紫粒大麦FLVs含量高于其他类型，且粒色与FLVs含量成正比。苯丙烷途径中关键酶包括查尔酮合酶（CHS）、查尔酮异构酶（CHI）、黄烷酮-3-羟化酶（F3H）、黄酮醇合酶（FLS）和二氢黄酮醇4-还原酶（DFR）等。干旱胁迫可增加大麦FLVs含量，并通过清除活性氧（reactive oxygen species，ROS）参与耐旱。大麦在灌浆期尤其易受干旱影响。研究人员此前发现108个FLVs生物合成途径酶，其外显子-内含子结构、基序、顺式调控元件、染色体分布和基因复制存在显著变异。因此，本研究进一步探索紫大麦在干旱胁迫下的生理、转录组和代谢组变化，以指导未来耐旱大麦育种。

研究人员以紫大麦品种Ziguangmangluoerleng为材料，对21 d幼苗采用15%聚乙二醇（polyethylene glycol，PEG）处理24 h模拟干旱，并开展生理、转录组、代谢组和计算分析。研究结论显示：干旱胁迫增加FLVs、可溶性糖、脯氨酸和抗氧化酶活性，并优化气孔参数；转录组鉴定出826个上调差异表达基因（differentially expressed genes，DEGs），富集于苯丙烷类和FLVs生物合成途径；代谢组检测到3,335个差异积累代谢物（differentially accumulated metabolites，DAMs）；20个高度互作蛋白由10个FLVs生物合成酶和10个转录因子（transcription factors，TFs）组成，在干旱胁迫下显著上调，其启动子含有关键干旱响应顺式元件；分子对接预测了候选蛋白质–配体结合。该研究为FLVs介导的大麦干旱响应分子机制提供了整合性见解，并为开发耐旱大麦品种提供了候选靶点。

关键技术方法概括（不超过250字）：研究人员以紫大麦品种Ziguangmangluoerleng的21 d幼苗为样本，用15% PEG模拟干旱24 h；测定气孔、FLVs、可溶性糖、脯氨酸及超氧化物歧化酶（SOD）、过氧化物酶（POD）、过氧化氢酶（CAT）活性；进行RNA-Seq转录组分析及差异表达、GO和京都基因与基因组百科全书（KEGG）富集、可变剪接、新基因注释和主成分分析（PCA）；进行代谢组检测，以PCA和正交偏最小二乘判别分析（OPLS-DA）鉴定DAMs；开展转录组-代谢组关联和加权基因共表达网络分析（WGCNA）；对20个候选基因进行染色体定位、基因结构、2,000 bp启动子顺式元件、三维蛋白结构、TFs和miRNA互作预测及分子对接；用相对mRNA表达验证候选基因。

研究结果如下。

3.1 干旱胁迫改变大麦多条信号通路的生理。通过气孔和生理生化测定，干旱使紫大麦气孔长度、宽度和密度降低，同时增加总FLVs、脯氨酸、可溶性糖含量及SOD、POD、CAT活性，表明多条干旱信号通路协调参与缓解胁迫。

3.2 干旱胁迫下转录组特征改变，包括FLVs生物合成途径主要基因表达。RNA-Seq数据总体映射率为87.82%–90.62%，唯一映射率为79.99%–83.18%。DEGs参与多种生物过程、细胞组分和分子功能，KEGG富集显示其参与苯丙烷生物合成、激素信号和胁迫信号通路。干旱还改变可变剪接（alternative splicing，AS）模式，TSS/TTS型AS较丰富；鉴定出新基因和功能注释基因。826个上调DEGs中，苯丙烷途径有80个DEGs，其中42个上调；FLVs生物合成途径有20个DEGs。启动子顺式元件在调控中发挥重要作用。

3.3 干旱胁迫下FLVs及其他信号通路的代谢物积累改变。代谢组分析鉴定出3,335个DAMs，其中1,827个上调、1,508个下调；层次聚类显示对照与干旱处理组清晰分离。KEGG富集显示氨基酸tRNA生物合成、氰基氨基酸代谢、甘油磷脂代谢、组氨酸代谢等显著相关，ABC转运蛋白相关代谢物最多，flavone、flavonol和anthocyanin生物合成途径也富集。鉴定出luteolin、quercitrin、isovitexin、isorhamnetin、epigallocatechin、procyanidin B2、kaempferol 3-O-gentiobioside和peonidin 3-O-glucoside等FLVs相关代谢物，并可映射到CHS、CHI、F3H、FLS、DFR。该结果支持FLVs相关代谢参与大麦干旱响应，但作者指出仍需靶向定量、代谢通量或遗传方法进行功能验证。

3.4 差异基因表达引导相应代谢物积累。整合转录组和代谢组后，对照与干旱处理样本清晰分离，基因-代谢物关联显示苯丙烷生物合成和flavone/flavonol生物合成途径的协调响应。WGCNA鉴定出八个共表达模块，模块-性状关系支持FLVs生物合成属于更广泛的干旱响应调控网络。十个干旱响应基因与高积累FLVs代谢物显著相关，提示FLVs生物合成是干旱响应网络的重要功能模块。

3.5 20个基因和转录因子显著互作FLVs生物合成途径的分子机制。对20个高互作蛋白（10个FLVs途径酶和10个TFs）进行结构表征，发现氨基酸长度从119（WRKY50）至581（HORVU4Hr1G062830），分子量从13.01 kDa至63.96 kDa，理论等电点（pI）从4.85至10.25；总平均疏水性（GRAVY）有正有负，亚细胞定位涉及细胞核和胞外。20个基因分布于8条染色体。外显子-内含子结构存在显著变异。2,000 bp启动子区顺式元件中，MeJA响应元件最多（132），其次为ABA响应元件（83）、干旱诱导元件（18）和防御/胁迫响应元件（7）。三维结构显示α螺旋、β链、跨膜螺旋和卷曲结构。TFs调控网络预测1,561个TFs来自32个家族，其中ERF 409、MYB 90、bZIP 85、WRKY 83、NAC 73、bHLH 61。miRNA–基因网络预测10个miRNA与12个基因互作。分子对接显示L-phenylalanine平均结合自由能最低（ΔG = ?5.22 kcal/mol）；HORVU0Hr1G021680与4-hydroxy-3-methoxycinnamaldehyde和L-phenylalanine预测结合最强（ΔG = ?6.3 kcal/mol）；WRKY50与quercetin（ΔG = ?6.7 kcal/mol）、kaempferol（?6.2 kcal/mol）和myricetin（?6.1 kcal/mol）预测结合较强。20个基因在干旱下均显著上调，但预测互作仍需实验验证。

3.6 干旱胁迫上调不同干旱信号通路标记基因。研究人员检测了HvCPK2a、HvCaM1、HvCAMTA4、HvABI5、HvSNAC1、HvPP2C59、HvRBOHF2、HvCAT2、HvFC1和HvFC2的相对mRNA表达，发现其在干旱下均显著上调。这些基因参与钙信号、ABA信号、ROS稳态、渗透调节和叶绿体逆行信号等过程，支持大麦干旱信号通路的协调响应。

3.7 协调的干旱信号通路整合于大麦耐旱性。研究人员提出两个整合模型。模型一描述膜感知、Ca²⁺/ROS/激素信号、叶绿体逆行信号、转录和转录后调控及下游生理输出，包括气孔关闭、抗氧化酶诱导、渗透调节和FLVs积累。模型二强调Ca²⁺/NO信号、SUMO化、DRP释放、ER未折叠蛋白反应以及线粒体和叶绿体ROS/基因组稳定性互作。模型中H₂O₂、O₂^?、CO₂、Na⁺和H⁺相关节点与本研究基因表达和生理数据相整合，说明钙、ROS、激素、TFs和亚细胞动态协同缓解干旱胁迫。

讨论与结论总结：未来展望部分提出，尽管研究鉴定了候选基因、代谢物和调控途径，但其生物学功能仍需实验验证。合成生物学应作为未来研究方向，包括CRISPR基因组工程、合成启动子、可编程基因电路和诱导调控系统；实际应用需体外、体内和田间验证。研究人员提出合成遗传电路包含传感器、整合器和执行器模块，并展示两输入NOR门电路。局限性部分指出，当前证据主要为相关性，预测的调控关系、转录组-代谢组整合和分子对接仅鉴定候选基因、TFs和代谢物，未建立直接因果机制；未来需基因敲除或过表达、靶向代谢物定量、生化检测、蛋白质–配体结合实验和生理验证。

结论翻译：大麦耐旱性对于缓解气候变化对农业生产力的不利影响日益重要。常规育种和单/多基因转基因已取得显著进展；然而进一步提高耐旱性仍具挑战。本类研究导向综述将生理、转录组、代谢组、调控网络和计算证据整合为统一框架，突出FLVs生物合成作为大麦干旱响应的核心组分。关键生理标志物包括FLVs、可溶性糖、脯氨酸、抗氧化活性和气孔关闭增加，均与干旱适应一致相关。FLVs生物合成及其他途径中的差异表达基因和代谢物提示协调的干旱响应调控网络。FLVs途径中20个显著上调基因、其相应代谢物及分子对接分析鉴定出可能促进干旱适应的候选调控组分。总体而言，整合分析一致支持FLVs生物合成与大麦干旱响应分子网络之间的关联。根据当前研究导向综述的发现，建议未来采用基因敲除、过表达、酶活性测定、代谢物饲喂实验和转基因验证研究，以建立FLVs介导的大麦耐旱性直接因果机制。此类经实验验证的知识最终可能通过精准育种、基因组工程和合成生物学方法促进耐旱大麦品种开发。