《Metabolomics》：Comparative metabolomic profiling reveals salinity tolerance mechanisms in a rice introgression line
编辑推荐：
摘要：水稻（Oryza sativa）对盐度极敏感，但盐胁迫下耐盐的代谢机制仍不完全清楚。在本研究中，研究人员对耐盐渐渗系JN100（JN）、供体亲本Nona Bokra（NB）和轮回亲本Jupiter（JU）进行叶片组织特异性非靶向代谢组学（untarget
摘要：水稻（Oryza sativa）对盐度极敏感，但盐胁迫下耐盐的代谢机制仍不完全清楚。在本研究中，研究人员对耐盐渐渗系JN100（JN）、供体亲本Nona Bokra（NB）和轮回亲本Jupiter（JU）进行叶片组织特异性非靶向代谢组学（untargeted metabolomics）分析，以表征盐胁迫下的代谢响应。比较分析鉴定出跨越氨基酸、糖与碳水化合物、脂质、有机酸、辅因子、电子载体和核苷酸等多种化学类的差异积累代谢物（differentially accumulated metabolites，DAMs）。在盐胁迫（salt stress，SS）下，JN相对于JU检测到201个DAMs（89个上调、112个下调）。显著差积代谢物包括尿囊素、大豆黄素（glycitin）、烟酰胺核糖苷（nicotinamide ribotide）、D-阿拉伯糖酸-1,4-内酯（D-arabinono-1,4-lactone）、堇黄素（violanthin）、L-甲硫氨酸亚砜（L-methionine S-oxide）、核糖醇（ribitol）、赖氨酸、芦丁（rutin）、谷氨酰胺、泛酸（pantothenic acid）和奎尼酸（quinic acid）。通路富集分析显示精氨酸生物合成、嘌呤代谢、丙氨酸-天冬氨酸-谷氨酸代谢显著富集，表明盐胁迫下氮与能量相关代谢通路广泛重编程。将SS实验的转录组与代谢组数据集整合，进一步在JN vs. JU比较中鉴定出10个与代谢物网络相关的差异表达基因（differentially expressed genes，DEGs）。其中OsDHQDT/SDH、OsFd-GOGAT、苯丙氨酰-tRNA合成酶、OsP5CS1、OsP5CS2和吡哆醛磷酸依赖转移酶与莽草酸、氨基酸和脯氨酸代谢相关代谢物关联。
研究背景方面，水稻作为甘薯型植物对盐胁迫高度敏感，耐盐阈值仅约3 dS/m，而土壤盐渍化阈值高于4 dS/m。已有研究提示脯氨酸、海藻糖、半乳糖醇、棉子糖、黄酮类、谷胱甘肽等参与渗透调节与活性氧（reactive oxygen species，ROS）清除，但学界仍缺乏在渐渗背景下把耐盐相关代谢特征与遗传背景差异区分开的完整代谢快照。沿海地方品种虽耐盐却常伴低产，Nona Bokra耐盐但具落粒、光敏等驯化劣势，Jupiter高产但盐敏感。为弥补背景基因型变异掩盖耐盐信号这一空白，研究人员以JN100渐渗系、NB供体亲本、JU轮回亲本为对象，在叶组织水平开展非靶向代谢组学并整合已有转录组，以界定耐盐代谢签名及系统层面机制。该研究发表于《Metabolomics》。
主要技术方法上，样本来自路易斯安那州立大学农业中心温室培养的JN100、Jupiter、Nona Bokra苗期叶片，盐处理以6 dS/m两步升至12 dS/m电导率（electrical conductivity，EC），每处理三生物学重复。研究人员采用液相色谱-四极杆飞行时间质谱（liquid chromatography-mass spectrometry，LC-MS；Q-TOF）的HILIC酰胺柱广靶检测极性代谢物，用Progenesis QI做峰拾取与注释，METLIN等库做二级谱注释；随后在MetaboAnalyst 6.0做对数变换、自动标度、单变量与正交偏最小二乘判别分析（orthogonal partial least squares discriminant analysis，OPLS-DA）、KEGG通路富集与代谢物集富集分析（metabolite set enrichment analysis，MSEA）；并以MetaboAnalyst联合通路模块和STITCH数据库整合既往JN100转录组差异表达基因与本研究差异积累代谢物，置信度≥0.400保留基因-代谢物互作。
研究结果如下。3.1 全局代谢谱揭示盐胁迫的基因型特异响应：通过OPLS-DA与UpSet分析，JN-JU盐迫下201个DAMs（89上/112下），JN、JU、NB各自盐迫DAMs分别为134、140、156个；氨基酸类最多（25个），说明叶代谢按基因型重编程，DAMs涵盖氨基酸、碳水、脂肪酸、嘌呤嘧啶等。3.2 多变量分析鉴定盐下关键判别代谢物：OPLS-DA将对照与盐胁迫清晰分簇，JN-JU预测分量在盐迫下解释91.2%方差；VIP>1的TOP代谢物中D-阿拉伯糖酸-1,4-内酯、堇黄素、L-甲硫氨酸亚砜、核糖醇、赖氨酸、芦丁、谷氨酰胺、泛酸、奎尼酸、乳糖醇在JN-JU盐迫特有富集，Berberine、水杨酸、L-精氨酸等在JN特有，JU特有L-天冬酰胺等，NB特有N5-甲基谷氨酸等。3.3 JN100与Jupiter盐迫下差异积累代谢物：通过火山图与倍数变化分析，大豆黄素在JN-JU对照20倍、盐迫25倍累积且JU单独盐迫下未检出；D-阿拉伯糖酸-1,4-内酯（抗坏血酸前体）盐迫23倍、核糖醇14—19倍上调；尿囊素、堇黄素、L-甲硫氨酸亚砜、赖氨酸、芦丁下调，反映ROS清除、氧化还原与渗透调整协同改变。3.4 适应代谢性可塑性赋予水稻耐盐性：KEGG拓扑分析显示JN-JU盐迫下精氨酸生物合成（pathway impact 0.494）、丙氨酸-天冬氨酸-谷氨酸代谢（0.787）影响最大，还富集乙醛酸与二羧酸代谢、磷酸戊糖途径（pentose phosphate pathway，PPP）、谷胱甘肽代谢、嘌呤代谢；JU侧重精氨酸与三氨基酸代谢重排，NB侧重精氨酸、半乳糖、谷胱甘肽、精氨酸-脯氨酸代谢，说明渐渗区段带来氮同化与抗氧化重构。3.5 联合通路分析：把转录组DEGs与代谢DAMs映射到Oryza sativa KEGG，JN-JU盐迫下精氨酸生物合成、甘氨酸-丝氨酸-苏氨酸代谢、谷胱甘肽代谢、D-氨基酸代谢、嘌呤代谢共富集；赖氨酸生物合成、亚油酸代谢、果糖-甘露糖代谢仅盐迫激活，硫代谢、泛酸与CoA生物合成、β-丙氨酸代谢也盐迫特异，JU则富乙醛酸与二羧酸代谢。3.6 整合通路分析揭示转录-代谢协同响应：STITCH网络里OsGLT1、OsGLT2、OsFd-GOGAT三个谷氨酸合酶基因居核心，将谷氨酰胺+2-氧戊二酸转为谷氨酸；OsP5CS1与OsP5CS2把谷氨酸转向脯氨酸，OsNAGS2催化精氨酸起始步，OsDHQDT/SDH连莽草酸/奎尼酸/赤藓糖-4-磷酸，苯丙氨酰-tRNA合成酶连苯丙氨酸，OsGMPS1连嘌呤核苷酸，吡哆醛磷酸依赖转移酶跨氨基酸与辅因子；所有保留互作合并分值0.988—0.999，说明氨基酸、渗透调节、辅因子代谢在耐盐中收敛。
讨论与结论部分总结：研究人员指出，渐渗系JN100相对JU的耐盐性不是单一代谢物决定，而是氨基酸代谢（精氨酸、脯氨酸、谷氨酸节点）、氮同化（GOGAT簇）、氧化还原防御（谷胱甘肽、黄酮、抗坏血酸前体）、碳-氮平衡（PPP、乙醛酸循环、嘌呤）与莽草酸分支共同重编程的结果。Nona Bokra供体背景提供耐盐代谢骨架，Jupiter背景代表敏感代谢表型，JN100通过渐渗区段把供体耐盐签名固定到高产背景中。该工作把叶组织非靶向代谢组与前期转录组衔接，提供了可区分背景变异的盐迫代谢地图，并为水稻耐盐marker与候选基因（OsFd-GOGAT、OsP5CS1/2、OsDHQDT/SDH、OsNAGS2）挖掘提供依据。
结论翻译：比较代谢组学分析揭示，水稻耐盐渐渗系JN100在盐胁迫下通过精氨酸生物合成、丙氨酸-天冬氨酸-谷氨酸代谢、嘌呤代谢、谷胱甘肽与磷酸戊糖途径、莽草酸及脯氨酸代谢的协同重编程实现适应；转录组-代谢组整合定位到OsFd-GOGAT、OsP5CS1、OsP5CS2、OsDHQDT/SDH、OsNAGS2、苯丙氨酰-tRNA合成酶、OsGMPS1与吡哆醛磷酸依赖转移酶等节点，证明氮同化、渗透调节、ROS清除和辅因子代谢在渐渗背景中收敛为耐盐表型，为高产耐盐水稻设计的代谢-基因枢纽提供证据。
