整合转录组学与代谢组学揭示不同施氮水平下头花蓼黄酮生物合成的调控机制

《Metabolites》：Integrated Transcriptomics and Metabolomics to Reveal the Regulatory Mechanisms of Flavonoid Biosynthesis in Persicaria capitata Under Different Nitrogen Fertilization Levels

　　

摘要
背景：头花蓼（*Persicaria capitata*）为多年生草本植物，是贵州省十大苗药之一，黄酮类化合物（flavonoids）为其主要活性成分。农业生产实践表明，施用氮肥可降低头花蓼中黄酮的积累，但其潜在的分子调控机制尚不清楚。方法：本研究采用盆栽试验，探究不同氮处理对头花蓼生长表现和产量的影响，并系统分析了相关转录组特征和代谢谱的变化。结果：施氮显著促进了茎 elongation、分枝数和产量形成。在0.3 g N/kg处理组中，最大茎长、分枝数和单株产量分别达到80.74 cm、37.30 和123.64 g/株。转录组筛选共获得10,880个差异表达基因（differentially expressed genes，DEGs），且所有比较组中下调DEGs的数量始终高于上调DEGs。代谢组分析共鉴定出3356个差异积累代谢物（differentially accumulated metabolites），其中包括168个差异黄酮代谢物。随着施氮量的增加，大多数黄酮类化合物（如绿原酸、5-羟基阿魏酸和4-羟基肉桂酸）表现出显著下调的积累模式。转录组与代谢组联合分析显示，差异黄酮代谢物和功能基因主要富集于苯丙烷生物合成（phenylpropanoid biosynthesis）、黄酮生物合成（flavonoid biosynthesis）以及黄酮和黄酮醇代谢（flavone and flavonol metabolism）通路。关键结构基因（包括 *PAL*、*4CL*、*CHS*、*FLS*、*F3H* 和 *F3′H*）表达下调，抑制了核心黄酮组分（如绿原酸、柚皮素和木犀草素）的生物合成与积累，从而导致头花蓼药用品质下降。结论：施氮提高了头花蓼的产量，但抑制了黄酮的生物合成。黄酮代谢途径中关键基因的下调表达是导致黄酮含量降低的原因。本研究为优化氮素管理、协调头花蓼种植实践中产量形成与品质提升提供了理论支持。

论文解读文章

头花蓼（*Persicaria capitata*）是贵州省特色苗药品种之一，其干燥全草或地上部分入药，具有清热、利尿、解毒、镇痛等功效。黄酮类化合物是头花蓼的核心药理活性成分。然而，在实际农业生产中，为追求高产，种植户往往过量施用氮肥，导致药材中黄酮含量显著下降，药材质量不达标，造成经济损失。已有研究表明，氮素水平能够显著影响药用植物次生代谢产物的积累，例如过量氮输入会降低苹果果实中总黄酮含量，高氮供应也会减少积雪草和淫羊藿中的黄酮产量。但针对头花蓼，氮肥如何调控其黄酮生物合成的分子机制尚不明确。因此，阐明施氮条件下头花蓼黄酮生物合成的调控规律，对于优化其栽培管理、协调产量与品质具有重要意义。

研究人员设置了四个氮素浓度（0、0.10、0.20、0.30 g N/kg干土）进行盆栽对比试验，测定了不同氮处理下头花蓼的生长形态指标和产量，并选取对照组（CK）、低氮组（T1）和高氮组（T3）三个关键梯度，结合转录组测序和广泛靶向代谢组学技术，系统分析了氮素调控黄酮生物合成的分子网络。

研究获得的主要结果如下。

**不同氮肥处理对头花蓼生长性状和产量的影响**：施氮显著促进了头花蓼的茎伸长、分枝增加和产量形成。随着施氮量增加，株高、分枝数和叶片氮含量均呈持续上升趋势。T3处理（0.30 g N/kg）下茎长最大（80.74 cm）、分枝数最多（37.30）、单株产量最高（123.64 g/株），较CK（14.2 g/株）增产约8.7倍。这表明氮肥是提高头花蓼生物量的关键营养因子。

**不同施氮水平下的转录组分析**：对CK、T1、T3三组叶片进行转录组测序，共获得103.79 Gb clean data，Q30比例均高于97%，组装完整性达96.50%。主成分分析显示不同氮处理间转录谱明显分离。共鉴定出10,880个差异表达基因（DEGs），其中T1_vs_CK有5356个（上调2404，下调2952），T3_vs_CK有6766个（上调2688，下调4098），T3_vs_T1有3832个（上调1508，下调2315）。所有比较组中下调DEGs数量均多于上调DEGs，且随施氮量增加，DEGs总数增加。GO富集分析显示DEGs主要富集于细胞过程、代谢过程、结合和催化活性等功能条目。KEGG通路富集分析发现，苯丙烷生物合成（ko00940）、黄酮生物合成（ko00941）以及黄酮和黄酮醇生物合成（ko00944）三条与黄酮代谢密切相关的通路显著富集，且随氮素升高，大多数相关DEGs呈现下调趋势。

**不同处理下的代谢组分析**：采用LC-MS/MS方法对叶片代谢物进行非靶向检测，共鉴定出3356种代谢物，其中黄酮类及其衍生物168种。差异代谢物（DAMs）分析显示，T1_vs_CK有924个DAMs，T3_vs_CK有1528个，T3_vs_T1有1096个。随施氮量增加，下调DAMs数量逐渐增多。KEGG富集分析同样指向黄酮生物合成相关通路。进一步聚焦差异黄酮代谢物（DAFs），在168种黄酮中，三组间共有105种显著差异黄酮。T3_vs_CK中63种黄酮下调，仅28种上调；T3_vs_T1中62种下调，仅13种上调。关键黄酮如绿原酸、5-羟基阿魏酸、4-羟基肉桂酸、柚皮素、木犀草素等均随施氮增加而显著下降。值得注意的是，头花蓼指标性成分槲皮苷（quercitrin）在高氮下也显著下调（log2FC = ?1.59）。

**转录组与代谢组联合分析**：将差异表达基因与差异代谢物映射到黄酮生物合成通路中，构建了氮素调控网络图。研究发现，苯丙烷途径中编码PAL（苯丙氨酸解氨酶）和4CL（4-香豆酸-CoA连接酶）的基因随施氮增加而表达下调，与下游产物4-羟基肉桂酸、绿原酸和5-羟基阿魏酸积累减少的趋势一致。在黄酮生物合成途径中，多数CHS（查尔酮合酶）基因表达下调，F3H（黄烷酮3-双加氧酶）和FLS（黄酮醇合酶）基因几乎全部下调，F3′H（类黄酮3′-单加氧酶）基因也随氮素增加而表达下降。这些关键酶基因的下调导致柚皮素、儿茶素、表儿茶素、花旗松素、木犀草素等核心黄酮组分含量降低。部分UGT（黄酮醇3-O-葡萄糖基转移酶）基因上调，可能促进槲皮素向杨梅素和芦丁的转化，进一步降低槲皮素含量。

**讨论与结论**：讨论部分指出，施氮增加头花蓼生物量但降低黄酮含量，其机制符合碳素营养平衡假说。高氮条件降低了植物C/N比，使光合碳骨架优先用于合成氨基酸、蛋白质等含氮化合物，抑制了莽草酸途径向苯丙烷代谢流的分配，导致黄酮等碳基次生代谢产物合成减少。本研究通过多组学联合分析，阐明了氮素调控头花蓼黄酮生物合成的分子机制，鉴定出PAL、4CL、CHS、FLS、F3H、F3′H和UGTs等关键基因及其下游代谢物的变化规律，为合理施氮、协调头花蓼高产与优质提供了理论依据。但关键基因的生物学功能仍需进一步通过转基因或基因编辑等手段验证，未来应结合田间试验确定最优施氮量，以促进头花蓼产业的可持续发展。