《Industrial Crops and Products》：Multi-omics analysis of leaf and root-nodule responses to combined heat and drought in R4-stage soybean
编辑推荐：
本研究对R4期（full-pod stage）大豆叶片与根瘤响应高温干旱复合胁迫（combined heat and drought, HD）进行了多组学整合分析。研究人员比较了耐性品种Heinong 65（HN65）与敏感品种Qihuang 34（QH34）
本研究对R4期（full-pod stage）大豆叶片与根瘤响应高温干旱复合胁迫（combined heat and drought, HD）进行了多组学整合分析。研究人员比较了耐性品种Heinong 65（HN65）与敏感品种Qihuang 34（QH34）在HD胁迫3天（HD 3 d）及恢复7天（re 7 d）后的表型与分子响应。结果显示，HN65表现出更强的光合恢复能力，其净光合速率、气孔导度、蒸腾速率及胞间CO2浓度均恢复至对照以上水平，而QH34上述指标仍显著受损；HN65的产量与种子品质性状相对稳定。在氧化还原层面，HN65在HD 3 d时活性氧（ROS）产生速率增幅较小（2.1倍 vs. QH34的6.3倍），超氧化物歧化酶（SOD）与过氧化物酶（POD）活性维持较高水平，丙二醛（MDA）含量在恢复期降至对照以下。miRNA分析鉴定出7个候选差异表达模块，其中保守型模块（如miR156-SPL6、miR168-AGO1）主要在HD 3 d发挥作用，而大豆或豆科特异的模块（如miR1510a/b-3p-RUN1/ROQ1、miR5559-SEOB）主要在re 7 d被招募。转录组与蛋白质组整合分析在HN65叶片HD 3 d时筛选出14个效应基因，包括3个上调基因（HSP90-5、CPN10-1及此前未表征的LHDI1）和11个下调基因，这些基因在加权基因共表达网络分析（WGCNA）中均为枢纽（hub）或候选枢纽基因，涉及蛋白质质量控制、叶绿体稳态及氧化还原调控。功能验证表明，在大豆毛状根复合植株中过表达LHDI1可提高HD胁迫下的存活率与相对含水量，降低叶片H2O2与O2-积累。综上，本研究揭示了R4期大豆叶片与根瘤应对HD胁迫的互补性组织响应程序，并鉴定了与耐性相关的关键分子节点。
研究背景与科学问题
大豆是全球重要的植物油与植物蛋白来源，栽培面积超过1.25亿公顷，生产正日益向半干旱与亚热带地区扩展，而这些区域干旱与高温事件常同时发生。气候变化预计将进一步增加此类复合胁迫事件的频率与持续时间。大豆在生殖发育阶段对高温干旱复合胁迫（combined heat and drought, HD）尤为敏感，因为开花、结荚、籽粒灌浆及种子成分依赖于持续的光合作用与源-库转运。R4期（full-pod stage）是连接结荚与籽粒灌浆的关键时期，此时生殖库对同化物供应的依赖性增强，然而针对该阶段HD胁迫及恢复期的分子研究仍然有限。已有研究表明，单一热胁迫或干旱胁迫下的表现无法可靠预测HD胁迫下的基因型表现，且HD胁迫会重组激素、免疫、氧化还原及代谢网络，因此亟需从多组学层面系统解析大豆R4期响应HD胁迫的分子机制，为培育气候适应性品种提供理论依据。
研究内容与主要发现
研究人员选取相对耐性品种Heinong 65（HN65）与敏感品种Qihuang 34（QH34），在R4期模拟HD胁迫处理3天（HD 3 d）后恢复正常条件恢复7天（re 7 d），系统比较了两品种在生理、产量、种子品质及多组学层面的差异。主要结论包括：（1）HN65在恢复期光合参数全面恢复至对照以上水平，而QH34仍显著受损；HN65的产量与种子成分保持稳定，表明耐性更体现在恢复能力而非单一胁迫指标的幅度。（2）HN65在HD 3 d时活性氧（ROS）积累幅度远小于QH34（2.1倍 vs. 6.3倍），抗氧化酶活性维持较高水平，且在恢复期丙二醛（MDA）含量降至对照以下，显示出更优的氧化还原稳态。（3）miRNA分析揭示保守型模块（如miR156-SPL6、miR168-AGO1）主要在HD 3 d参与防御与发育调控，而大豆或豆科特异的模块（如miR1510a/b-3p-RUN1/ROQ1、miR5559-SEOB）主要在re 7 d被招募，双荧光素酶实验证实了miR1510b-3p对RUN1与ROQ1的直接调控。（4）转录组与蛋白质组整合分析在HN65叶片HD 3 d时鉴定出14个效应基因，其中CPN10-1、HSP90-5及此前未表征的LHDI1上调，其余11个基因下调；所有效应基因在WGCNA中均为枢纽（hub）或候选枢纽基因，涉及蛋白质质量控制、叶绿体稳态、氧化还原调控及次生代谢等通路。（5）大豆毛状根复合植株实验表明，过表达LHDI1可显著提高HD胁迫下的存活率与相对含水量，降低叶片H2O2与O2-积累。
关键技术方法
研究采用的主要技术方法包括：（1）生理生化测定：利用便携式光合仪测量气体交换参数，使用荧光探针法与比色法测定ROS产生速率、MDA含量及SOD、CAT、POD活性与脯氨酸含量；（2）多组学测序：对叶片与根瘤样本进行miRNA-seq与mRNA-seq（NovaSeq 6000平台），对叶片样本进行基于Orbitrap Astral高分辨质谱的超深度定量蛋白质组学分析（数据非依赖采集模式，DIA），以Glycine max Wm82.a6.v1为参考基因组；（3）生物信息学分析：运用DESeq2进行品种×处理互作分析，采用加权基因共表达网络分析（WGCNA）构建共表达网络并识别枢纽基因，结合STRING数据库构建蛋白-蛋白互作（PPI）网络，利用PlantCARE分析启动子顺式元件；（4）分子验证：通过qRT-PCR验证表达趋势，采用拟南芥原生质体双荧光素酶报告系统验证miRNA-靶基因调控关系，在烟草叶片中瞬时表达进行亚细胞定位，利用发根农杆菌K599介导的大豆毛状根复合植株体系进行LHDI1过表达的功能评估。
研究结果分述
3.1 HN65维持光合恢复与农艺表现：HD 3 d时QH34出现严重萎蔫与茎弯曲，而HN65仅中度叶片下垂。re 7 d时HN65的净光合速率、气孔导度、蒸腾速率及胞间CO2浓度均恢复至对照以上水平，而QH34仍显著受损。成熟期QH34单株种子重量显著降低，百粒重、种子体积、粗蛋白、粗脂肪及水分含量均下降，而HN65这些性状相对稳定。
3.2 HN65限制急性ROS积累：HD 3 d时QH34的ROS产生速率升高6.3倍，HN65仅升高2.1倍。HN65基础SOD与POD活性更高，且SOD在HD 3 d与re 7 d均可诱导，而QH34的SOD活性下降。HN65的MDA在re 7 d降至对照以下，QH34则仍略高于对照。脯氨酸积累在QH34中更为显著（2.8倍 vs. 1.36倍）。
3.3 保守与谱系限制性miRNA模块的阶段依赖性模式：鉴定出7个候选HD响应的miRNA-靶基因模块。保守型模块（miR156-SPL6、miR168-AGO1、miR156-MYB33）主要在HD 3 d活跃，而谱系限制性模块（miR5559-SEOB、miR1514a-NAC014、miR1510a/b-RUN1/ROQ1）主要在re 7 d被招募。双荧光素酶实验证实miR1510b-3p可直接调控RUN1与ROQ1。
3.4 叶片mRNA-蛋白质组整合与组织转录组揭示互补响应：HN65叶片在HD 3 d富集热响应、蛋白折叠修复、抗氧化及光保护通路，同时抑制基础代谢与生长；根瘤则大部分代谢通路受抑，但选择性诱导CTP生物合成、SAM生物合成及MAPK信号通路。re 7 d时叶片转向大分子降解，根瘤重新激活代谢与蛋白加工通路。
3.5 品种拮抗响应基因的功能特征：DESeq2互作分析鉴定出81个在HN65与QH34中呈相反响应趋势的基因，其中HD 3 d时HN65叶片优先诱导胁迫响应与抗氧化基因，re 7 d时所有47个拮抗基因均在HN65中受抑而在QH34中诱导，主要涉及胁迫响应、免疫调节与激素信号。PPI网络显示HD 3 d时HSC-2蛋白为核心连接点，re 7 d时网络围绕CPK28、SNAP33及CAMTA3等免疫相关蛋白展开。
3.6 叶片跨组学整合鉴定HD效应层：从367个叶片严格差异丰度蛋白（SDEP）中，经全长序列比对保留36个蛋白-基因匹配，进一步整合转录组获得14个mRNA与蛋白变化一致的基因。所有14个基因在mRNA与蛋白WGCNA中均为枢纽或候选枢纽，跨组学模块映射显示多数基因保留了至少一个共模块伙伴。共表达网络将诱导臂与抑制臂分开，诱导效应子关联蛋白质质量控制与叶绿体稳态，抑制效应子关联氨基酸生物合成与光合作用。
3.7 14个效应子的分子特征与LHDI1过表达的表型评估：CPN10-1与HSP90-5定位于叶绿体，LHDI1呈弥散的核质分布。所有14个效应子启动子富含ABA响应元件及MYB/MYC结合位点。大豆毛状根复合植株中过表达LHDI1可减轻叶片萎蔫，提高存活率与相对含水量，减弱DAB与NBT染色信号。
讨论与结论总结
讨论部分指出，生殖期HD耐性更应通过恢复能力而非单一胁迫指标来定义。HN65在恢复期完全恢复气体交换并维持产量与种子成分，而QH34光合恢复不完全且产量损失更大。氧化还原行为同样支持这一观点：QH34在HD 3 d时大量积累脯氨酸并未阻止re 7 d时的MDA积累，表明高渗透调节物质或酶活性本身并非耐性证据。HN65的特征在于HD 3 d时限制ROS产生并适度积累脯氨酸，配合快速诱导伴侣蛋白与抗氧化转录响应，并在re 7 d时衰减免疫相关转录。
叶片与根瘤呈现互补的功能程序：叶片采取主动保护策略（抑制生长、激活抗氧化与光保护），根瘤则广泛降低代谢但保留核苷酸与信号通路；恢复期叶片转向损伤清除，根瘤重新激活代谢。miRNA数据为此分工增加了阶段依赖性调控层次，保守型模块在胁迫高峰发挥作用，谱系限制性模块在恢复期被招募。
14个跨组学效应基因构成一个与蛋白质质量控制、叶绿体稳态及氧化还原调控紧密关联的效应网络。LHDI1作为此前未表征的基因，其过表达可改善HD胁迫下的生存率与水分状态，减少氧化损伤，提示其可能通过泛素介导的蛋白降解途径参与胁迫适应。
研究结论：HN65的相对耐性与快速防御相关转录、恢复期免疫信号的及时衰减、受限的ROS积累以及更有效的光合恢复密切相关。叶片与根瘤的时空特异性程序以及保守与谱系限制性miRNA模块的阶段依赖性招募共同构成了大豆R4期应对HD胁迫的协调调控网络。
