草铵膦残留对油菜的慢性植物毒性：转录组-代谢组解偶联机制研究

《Plant Stress》：Transcriptome–Metabolome Uncoupling Underlies Chronic Phytotoxicity of Glufosinate Residues in Rapeseed

摘要专业翻译 草铵膦（Glufosinate-ammonium, GLA）在全球范围内的使用日益广泛，但其土壤残留对后茬作物的慢性植物毒性仍知之甚少。为填补这一知识空白，研究人员整合了生理学、超微结构、转录组学及代谢组学分析，研究了油菜（Brassica napus L.）在三种土壤GLA残留浓度及五个移栽间隔条件下的响应，所有测定均在移栽后30天（days after transplanting, DAT）进行。叶片和根系中的GLA残留水平随处理浓度升高而增加，并随移栽间隔延长而降低。GLA残留与可见的生长抑制及光合作用下降相关，且呈浓度和时间依赖性，表现为叶绿素含量和净光合速率（net photosynthetic rate, Pn）的降低。叶绿体超微结构损伤随浓度升高而逐渐加重。谷氨酰胺合成酶（glutamine synthetase, GS）活性受到抑制，导致铵离子积累、可溶性蛋白含量下降，值得注意的是，在30 DAT时观察到谷氨酸积累，这一模式与经典的急性毒性范式——谷氨酸耗竭——形成鲜明对比。相比之下，核酮糖-1,5-二磷酸羧化酶（Ribulose-1,5-bisphosphate carboxylase, Rubisco）活性在所有GLA处理中均被一致抑制至相似的低水平。至第30天，抗氧化系统表现出功能性损伤，可能与其能量限制有关：尽管抗氧化相关基因在转录水平上仍保持响应性，但H2O2水平无显著变化，而丙二醛（malondialdehyde, MDA）积累，抗坏血酸-谷胱甘肽（ascorbate-glutathione, ASA-GSH）循环组分下降。转录组分析显示，三个处理组间共有1,967个共同上调的差异表达基因（differentially expressed genes, DEGs），但代谢组学仅鉴定出56个共同上调的代谢物，且高度集中于以谷氨酸为中心的氮代谢——这表明存在深度的转录组-代谢组解偶联。值得注意的是，M3T8组的DEG数量超过了急性胁迫组M2T1，表明在本盆栽试验的高剂量条件下，8天移栽间隔可能不足以确保安全。综上所述，这些结果表明，在30 DAT时，GLA暴露与一种独特的“慢性衰退”模式相关——叶绿体解体、谷氨酸积累、抗氧化系统损伤以及转录组-代谢组解偶联——这一模式与急性氨驱动的致死机制截然不同。

论文解读文章

一、研究背景与科学问题

油菜（Brassica napus L.）是全球最重要的油料作物之一，为食用油和饲料蛋白生产做出重要贡献。有效的杂草管理是保障油菜稳产的前提，而随着百草枯在多个国家被逐步禁用，草铵膦（Glufosinate-ammonium, GLA）凭借其速效性和在抗除草剂种植体系中的适用性，迅速成为主流的广谱非选择性除草剂。由此，GLA的施用量及在农田土壤中的残留负荷大幅增加，其持久性和对后茬作物的残留效应引发广泛关注。

在分子层面，GLA通过抑制谷氨酰胺合成酶（glutamine synthetase, GS）——铵同化途径的关键酶——阻断氮素循环，进而引发一系列植物毒性事件，包括细胞内铵离子（NH₄⁺）积累、光呼吸受损、叶绿素降解和光合作用抑制。急性叶面施用GLA已被充分证实可诱导膜脂过氧化（MDA升高）、可溶性蛋白丧失以及活性氧（reactive oxygen species, ROS）激增并压倒抗氧化防御系统。这些急性响应已被较为系统地表征。

然而，移栽油菜所面临的GLA暴露农业生态情境与单次急性剂量存在本质差异。土壤施用后，GLA以可测水平持续残留，短期内移栽的幼苗可被其吸收。现行推荐的移栽前安全间隔期——通常基于急性毒性阈值或单时间点测定——假定延迟移栽即可恢复作物安全性，但这一假设在现实的浓度-时间组合情境下很少被检验。土壤残留GLA的慢性、亚致死及分子水平效应仍缺乏系统研究。此外，以往工作多依赖孤立的生理终点指标，无法捕捉慢性胁迫下从调控到执行的动态过程。整合转录组与代谢组分析能够揭示转录重编程如何转化——或未能转化——为代谢调整，为解析慢性胁迫响应机制提供强有力手段。

研究人员据此提出假设：在慢性土壤残留GLA暴露下，转录与代谢层面的协调可能发生解偶联——植物启动强烈的转录“ distress ”响应，但代谢执行因能量耗竭和靶蛋白失活而崩溃。

二、研究设计与技术方法

研究人员建立了全因子浓度×时间试验设计：在GLA预处理的土壤中，以三个残留水平（M1为225 g a.i. ha^?1、M2为450 g a.i. ha^?1、M3为900 g a.i. ha^?1）和五个移栽间隔（喷药后0、1、3、5、8天，分别记为T0、T1、T3、T5、T8）移栽60日龄油菜幼苗，以清水处理为对照（CK），共20个处理组合、5次重复，总计100盆。盆栽试验于2025年12月18日在浙江省农业科学院作物与核技术利用研究所的天然光照玻璃温室内进行，土壤采自休耕田0–20 cm土层。所有样品均在移栽后30天（DAT）统一采集。

主要关键技术方法包括：（1）采用液相色谱-串联质谱法（LC-ESI-MS/MS）按国家标准GB23200.122-2026测定叶片和根系中GLA残留含量；（2）使用透射电子显微镜（TEM）观察叶绿体超微结构；（3）利用商业试剂盒测定叶绿素含量、GS和Rubisco酶活性、铵态氮/谷氨酸/可溶性蛋白含量以及抗氧化酶活性和氧化损伤指标；（4）基于Illumina NovaSeq 6000平台进行转录组测序（RNA-seq），以|log₂FC| ≥ 1且FDR < 0.01为标准鉴定差异表达基因（DEGs）；（5）采用超高效液相色谱-高分辨Q-Exactive Orbitrap质谱联用（UPLC-MS/MS）进行非靶向代谢组学分析，以|log₂FC| ≥ 1、VIP ≥ 1且p < 0.05为标准鉴定差异代谢物（DEMs）。转录组数据已提交至NCBI SRA数据库（BioProject编号PRJNA1514845）。

三、主要研究结果

3.1 油菜幼苗对GLA胁迫的生理响应

在T0、T1、T3和T5时，较高GLA浓度与更严重的可见生长抑制相关；至T8时，各处理组与对照间无显著生长差异。叶绿素含量在T0–T3呈显著浓度依赖性下降，T5部分恢复（M1和M2与对照相当，M3仍降低），T8完全恢复。Rubisco活性在所有GLA处理组中均显著低于对照，且三个浓度组间无显著差异，但延长移栽间隔从T1至T8使Rubisco活性逐步恢复。净光合速率（Pn）和气孔导度（Gs）随GLA浓度升高和移栽间隔缩短而下降，T8部分恢复。胞间CO₂浓度（Ci）从T0至T5逐渐升高，至T8各处理组恢复至对照水平。Ci升高伴随Pn和Gs降低排除了气孔限制——CO₂固定已受损，表明叶肉同化能力受损是光合抑制的主导因素。

3.2 GLA对油菜幼苗叶绿体的损伤

透射电镜观察显示，对照叶绿体具有发育良好、紧密堆叠的基粒和规则排列的基质片层。M1处理下叶绿体类囊体膜轻度松弛但基本完整；M2处理下类囊体片层明显肿胀分散，基粒堆叠显著紊乱；M3T8处理下类囊体膜广泛崩解为不规则碎片或囊泡结构，大量质体小球聚集成簇，表明严重膜损伤和叶绿体衰老。叶绿体进行性解体直接解释了叶绿素含量的显著下降。

3.3 GLA对油菜叶片氮代谢的影响

各移栽时间下，叶片和根系中GLA残留量随处理浓度显著增加，同一浓度下延长移栽间隔使残留量逐渐下降。GS活性在所有处理组中均受显著抑制，呈浓度和时间依赖性。GS抑制导致其底物铵态氮和谷氨酸均显著积累，且浓度越高、间隔越短积累越多。可溶性蛋白含量呈浓度和时间依赖性下降，与GS抑制程度一致，证实GLA胁迫下蛋白质合成受损。

3.4 GLA对油菜氧化胁迫及抗氧化系统的影响

H₂O₂含量在所有GLA处理下均无显著变化。MDA含量在移栽间隔2天内的处理中显著升高，2天后的处理无显著变化。SOD、POD和GST活性在各处理组中均与对照无显著差异。APX活性显著降低，且浓度越高、间隔越短降低越明显。GPX和GR活性在2天移栽间隔内显著降低，其后无显著差异。GSH含量与对照相当，而GSSG、ASA和DHA含量均显著低于对照，且降低幅度随GLA浓度增加而增大。抗氧化系统表现为功能性损伤而非应激性增强。

3.5 油菜对GLA残留胁迫的转录组响应

对M2T1、M2T8和M3T8三个代表性处理组进行RNA-seq分析。与对照相比，M2T1组有8,135个DEGs（4,913上调、3,222下调）；M2T8组DEGs降至5,278个（4,075上调、1,203下调），较M2T1减少约35.1%；M3T8组DEGs最多，达10,148个（6,872上调、3,276下调）。三组间共有1,967个上调基因和319个下调基因。共同上调基因显著富集于氮代谢、丙氨酸/天冬氨酸/谷氨酸代谢、谷胱甘肽代谢和苯丙氨酸代谢等与GLA毒性直接相关的通路；共同下调基因主要富集于光合作用-天线蛋白、光合作用和脂肪酸代谢等通路。

3.6 油菜在草铵膦残留胁迫下的代谢组学特征

与对照相比，M2T1组有837个上调、637个下调DEMs；M2T8组有203个上调、141个下调DEMs；M3T8组有520个上调、362个下调DEMs。与转录组中1,967个共同上调基因形成鲜明对比，代谢组中仅鉴定出56个共同上调代谢物，富集于丙氨酸/天冬氨酸/谷氨酸代谢和D-氨基酸代谢两条通路，前者与谷氨酸积累的生理观察一致。55个共同下调代谢物富集于叶酸生物合成和α-亚麻酸代谢通路，分别与光呼吸受损和膜脂过氧化/叶绿体解体相关联。氨基糖和核苷酸糖代谢及硫代葡萄糖苷生物合成在所有三个处理组中均被富集，是油菜对GLA残留胁迫最保守的代谢响应。间隔效应比较显示，M2T8较M2T1缺失了丙酮酸代谢、TCA循环、糖酵解/糖异生、氧化磷酸化等能量代谢通路，表明延长间隔至8天有效缓解了急性能量代谢紊乱。浓度效应比较显示，M3T8重新富集了TCA循环、丙酮酸代谢、糖酵解/糖异生和氧化磷酸化等能量代谢通路，并新增赖氨酸降解、硫代谢、嘌呤代谢和乙醛酸/二羧酸代谢等通路，表明高浓度处理即使在8天降解后仍触发与急性胁迫相当的代谢紊乱，且抑制了恢复期代谢重编程。

四、讨论与结论

讨论部分要点：

关于安全移栽间隔的再评估：GLA残留对移栽油菜的植物毒性呈明确的浓度和时间依赖性模式。M3T8的转录组扰动程度（10,148个DEGs）超过急性胁迫组M2T1（8,135个DEGs），表明高初始浓度延长了“危险窗口”，8天移栽间隔对标准剂量可显著缓解转录组损伤（DEGs减少35.1%），但对高剂量仍不充分。浓度与时间的交互并非简单加和，慢性亚致死暴露持续驱动应激信号，导致比急性胁迫更广泛的转录响应和更深的代谢崩溃。

关于叶绿体解体驱动光合崩溃：叶绿体是慢性GLA暴露的首要靶细胞器，其超微结构进行性恶化与叶绿素含量和Rubisco活性下降密切平行。Rubisco在所有GLA处理组中均被一致抑制且无浓度差异，表明其催化功能对GLA胁迫极为敏感。Pn和Gs下降伴随Ci升高是典型非气孔限制的诊断特征，证明叶肉羧化能力受损——叶绿体功能障碍和Rubisco失活——是光合抑制的主导因素。

关于谷氨酸积累的双相转变：与经典急性毒性模型中GS抑制导致谷氨酸耗竭不同，30天慢性暴露实验中谷氨酸显著积累。可能机制包括：胁迫诱导的蛋白质降解释放游离谷氨酸；GS不可逆抑制阻断谷氨酸→谷氨酰胺转化，且能量限制或酶失活损害转氨作用、GABA支路等替代消耗途径；叶绿体解体导致区室化丧失，基质谷氨酸释放至胞质；谷氨酸脱氢酶（GDH）在铵积累下可能被激活直接合成谷氨酸。由于仅单时间点采样，谷氨酸先降后升的动态过程或合成/降解/消耗间的动态平衡均无法排除。

关于慢性能量限制下抗氧化系统的功能性损伤：MDA显著升高表明持续脂质过氧化，而H₂O₂无显著变化可能归因于叶绿体作为主要ROS来源的丧失。SOD、POD和GST活性不变反映能量耗竭下无法诱导新酶合成。APX活性下降及ASA/DHA/GSSG含量降低表明ASA-GSH循环再生能力受损，与ATP和NADPH供应不足一致。抗氧化基因在转录水平仍具响应性但酶活性未平行变化，支持能量依赖的功能性损伤模型。

关于转录组-代谢组解偶联：三组间1,967个基因共同上调而仅56个代谢物共同上调，存在35倍差异。M3T8 DEGs最多但差异代谢物未相应增加，表明代谢组存在响应幅度上限。解偶联表现为两种形式：转录激活而无代谢后果（MAPK信号、激素信号、脂质重塑等通路转录上调但代谢物无对应积累）；代谢下降独立于转录调控（叶酸生物合成基因无显著转录变化但叶酸含量显著下降）。其生物学基础可能涉及GS蛋白水平不可逆抑制和叶绿体——细胞能量工厂——的解体导致ATP和NADPH短缺，但ATP/NADPH未直接测定，该机制仍属假说。作者强调该现象仅在一种除草剂-作物-环境组合中观察到，不宜作为普适性标志，需多除草剂作用模式、多作物物种、多发育阶段和多环境条件下的验证。

结论部分翻译：

本研究通过整合生理学、超微结构、转录组学和代谢组学分析，系统研究了土壤残留GLA对油菜的慢性植物毒性。GLA残留与浓度和时间依赖性损伤相关，高剂量效应持续超过8天移栽前间隔。在机制上，GLA破坏叶绿体超微结构——从类囊体松弛到完全膜崩解——最终导致光合崩溃。GS抑制引起铵积累和谷氨酸显著增加，与急性毒性中的耗竭形成对比。Rubisco在各浓度下均被一致抑制。30天后，ASA-GSH组分下降，抗氧化酶活性基本无响应，表明抗氧化系统在慢性能量限制下发生功能性损伤。关键的是，研究人员发现了严重的转录组-代谢组解偶联：1,967个基因上调，却仅有56个代谢物积累。转录激活的胁迫通路（苯丙烷代谢、谷胱甘肽代谢、氮代谢）与代谢改变的通路（赖氨酸降解、硫代谢、TCA循环）完全分离，揭示胁迫信号传导与代谢执行之间的系统性错位。总之，慢性GLA植物毒性通过一种独特的“慢性衰退”模式运作——叶绿体解体、谷氨酸积累、抗氧化损伤和转录组-代谢组解偶联——与急性氨驱动的致死机制根本不同。研究人员提出转录组-代谢组解偶联作为油菜慢性GLA胁迫相关的一种现象，可作为该情境下系统性代谢功能障碍的候选指标。在多种除草剂、作物物种和环境条件下验证后方可推广至更广泛的普适性结论。