多组学解析氮介导盐胁迫下苦豆子器官特异性黄酮类与生物碱代谢通路

《Plant Stress》：Multi-omics elucidates organ-specific flavonoid and alkaloid metabolic pathways in Sophora alopecuroides under nitrogen-mediated salt stress

　　

摘要
土壤盐渍化严重威胁农业生产，而苦豆子（Sophora alopecuroides）等盐生药用植物通过胁迫诱导的生物活性化合物积累，为边际土地的可持续利用提供了解决方案。补充氮（N）可能增强耐盐性和代谢产物产量，但组织特异性机制仍不清楚。本研究确定了盐胁迫下苦豆子的最佳氮 regime，并利用整合表型分析、代谢组学和转录组学解析了器官分化的代谢适应。在200 mmol/L NaCl胁迫下的表型筛选表明，10 mmol/L混合氮（NH??:NO??比例为1:0）最大化生长和耐盐性，显著降低地下部盐敏感指数（SSI）81.1%（相对于N0，p < 0.05），同时使耐盐指数（STI）增加56.6%。观察到器官特异性响应：在此处理下地上部SSI降低22.5%，在所有氮 regime中持续低于地下部SSI。代谢组学分析检测到796种代谢物，其中黄酮类（103种化合物）占主导；根部显示含氧有机物、喹啉脂质和关键黄酮类（橙皮素、紫云英苷）积累增加，而叶片富集生物碱（血根碱）。转录组学表明根部响应更强（7477个差异表达基因 vs. 叶片4027个），其中3936个根部基因上调。KEGG分析突出了不同的初生代谢（叶片中氰氨基酸/酪氨酸；根中谷胱甘肽/色氨酸）但共同汇聚于苯丙烷生物合成。关键的是，多组学整合阐明了两个核心适应通路：根部类黄酮代谢物生物合成与SaCHS1～3表达正相关，叶片中异喹啉生物碱生物合成与SaTAT1/SaAAT1～2相关。这些发现为改良苦豆子耐盐性和药用化合物积累鉴定了候选基因位点。

基于多组学整合的苦豆子氮介导盐胁迫适应性代谢机制研究解读

**研究背景与科学问题**
土壤盐渍化是全球农业面临的严峻挑战，超过14亿公顷土地受其影响，预计至2050年全球50%的灌溉土地可能盐渍化，严重威胁粮食安全。盐生药用植物因兼具耐盐特性和药用价值，成为边际盐碱地开发利用的理想对象。苦豆子（Sophora alopecuroides）作为分布于中亚和西亚干旱半干旱区的豆科多年生草本植物，其生物碱和黄酮类化合物具有重要药用价值，同时可作为生物农药和生物肥料，还是潜在的优质饲草资源。然而，苦豆子对盐分敏感，当土壤总盐含量超过0.3%时生长受抑，达0.4%时无法正常生长，幼苗耐盐阈值约为200 mmol/L NaCl。氮素作为植物生长的必需大量元素，不仅调控初生代谢，还通过影响关键酶活性和相关基因表达间接调控黄酮、萜类等次生代谢产物的合成，从而增强植物抗逆性。然而，氮素形态和浓度如何组织特异性调控苦豆子在盐胁迫下次生代谢物积累的机制尚不清楚。基于此，本研究提出核心科学问题：合理施用氮肥能否通过调控苦豆子的代谢过程，在盐胁迫下促进植株生长和药用活性成分积累？

**研究设计与技术方法**
研究人员以新疆伊犁地区采集的苦豆子种子（千粒重23.70 g，发芽率>98%）为材料，在石英砂盆栽体系中设置200 mmol/L NaCl胁迫，比较3种NH??-N:NO??-N摩尔比（1:0、0:1、1:1）和4个氮浓度水平（0、2、10、50 mmol/L）共10个处理，处理60天后测定生长和生理指标。基于表型筛选结果，选取最适处理组（S组：无氮对照；SN组：10 mmol/L氮，NH??:NO??=1:1）的叶片和根部样本进行代谢组学和转录组学分析。代谢组学采用超高效液相色谱-串联质谱（UPLC-MS/MS）技术，转录组学采用Illumina平台测序并从头组装，通过qRT-PCR验证关键差异表达基因的表达模式。多组学数据整合分析采用双向聚类、主成分分析（PCA）和典范对应分析（CCA）等方法。

**研究结果**
**1. 氮素对盐胁迫下苦豆子幼苗生长特性的影响**
形态学指标测定显示，1:1铵硝比下10 mmol/L氮处理效果最佳，株高、茎粗、叶面积分别较无氮对照显著增加71.0%、18.6%和99.3%（p < 0.05）。生物量积累呈现相同趋势，1:1比例N10处理使地上部、地下部和总生物量分别显著增加60.9%、54.7%和56.7%。盐胁迫指数分析表明，1:1比例N10处理使地下部SSI显著降低81.1%，同时STI增加56.6%，表明平衡的铵硝混合供氮最有效增强根系耐盐性。双因素方差分析证实铵硝比和氮浓度间存在显著交互效应。

**2. 氮素对叶片和根部代谢组的影响**
共鉴定796种代谢物，黄酮类（103种，22.01%）为最丰富类别，其次为有机氧化合物（56种）、异戊烯醇脂质（53种）等。叶片中鉴定241个差异积累代谢物（DAMs），根部鉴定304个。根部上调DAMs以黄酮类（24种）、醌类（17种）和有机氧化合物（12种）为主，而叶片特异性上调DAMs主要是有机氧化合物（10种）、醌类（9种）和黄酮类（6种）。KEGG富集分析显示，叶片DAMs显著富集于苯丙烷生物合成、异喹啉生物碱生物合成、丙氨酸/天冬氨酸/谷氨酸代谢途径；根部DAMs则显著富集于黄酮生物合成、黄酮和黄酮醇生物合成及咖啡因代谢途径。根部DAMs数量和化学多样性均高于叶片，表明根部经历更广泛的代谢重编程。

**3. 氮素对叶片和根部转录组的影响**
转录组测序共获得8775个注释unigenes。表达分析显示，PC1（98%方差）完全分离叶片和根部样本，PC2（2%方差）区分氮处理差异。叶片中鉴定4027个DEGs（1529上调，2498下调），根部鉴定7477个DEGs（3936上调，3541下调），根部DEGs数量为叶片的1.86倍。GO富集分析显示叶片上调DEGs主要富集于代谢过程、细胞过程和刺激响应，而根部上调DEGs主要富集于酶活性和代谢过程相关功能。KEGG分析表明，叶片DEGs显著富集于氰氨基酸代谢、酪氨酸代谢以及异喹啉生物碱生物合成等途径，根部DEGs显著富集于谷胱甘肽代谢、色氨酸代谢等途径，两个组织共同富集于苯丙烷生物合成途径。

**4. 关键代谢通路分析**
通过DAMs-DEGs整合分析，研究人员鉴定了两个核心适应性通路。黄酮类代谢通路中，查尔酮合成酶（CHS）基因家族呈现组织特异性表达模式：SaCHS1与叶片中橙皮素积累正相关，SaCHS2与紫云英苷负相关，SaCHS3与根部飞燕草素（FC=5.81）和山柰素积累强正相关。异喹啉生物碱通路中，血根碱在叶片中显著上调（FC=4.83）而在根部下调（FC=0.17），其与SaAAT1在叶片中强相关，而与SaAAT2在根部相关；酪氨酸转氨酶基因SaTAT1则与叶片中血根碱积累呈负相关。典范对应分析（CCA）进一步验证了这些基因-代谢物关联的特异性。

**5. qRT-PCR验证**
选取13个关键DEGs进行qRT-PCR验证，结果显示qRT-PCR与RNA-seq数据高度正相关（叶片R2=0.9535，根部R2=0.9661，p < 0.001），证实了转录组数据的可靠性。

**讨论与结论**
本研究表明，氮素形态和浓度对盐胁迫下苦豆子的生长和代谢具有显著交互作用。平衡的铵硝混合供氮（1:1比例，10 mmol/L）通过协同效应最有效增强苦豆子幼苗的耐盐性和生物量积累，效应优于单一氮形态。代谢组分析揭示黄酮类化合物（103种）是苦豆子代谢组中最丰富的类别，氮素添加通过促进含氧有机物和喹啉脂质积累增强组织抗氧化能力。转录组分析表明根部对氮素添加的转录响应更强烈（7477个DEGs），与根部直接暴露于盐分胁迫的生理角色一致。关键的是，多组学整合阐明了组织特异性的代谢适应策略：根部通过SaCHS1～3激活黄酮生物合成通路实现盐胁迫适应，叶片通过SaTAT1/SaAAT1～2激活异喹啉生物碱通路完成响应。该研究鉴定的SaCHS1～3、SaTAT1和SaAAT1～2等关键基因为今后通过基因编辑或转基因技术改良苦豆子的耐盐性和药用成分定向合成提供了重要靶点，也为盐碱地药用植物资源的开发利用提供了理论依据。本研究发表于《Plant Stress》。