高羊茅响应低氮胁迫的整合生理与转录组分析：碳氮代谢重编程及bHLH候选调控因子

《Plants》：Integrated Physiological and Transcriptomic Profiling of Tall Fescue in Response to Low-Nitrogen Stress

　　

低氮（low nitrogen, LN）可利用性是限制植物营养生长和生殖发育的主要约束。为解析高羊茅（Festuca arundinacea）低氮胁迫的生理与分子机制，研究人员开展了表型、生理和转录响应的综合分析。结果表明，低氮胁迫显著降低株高、分蘖数、叶长和地上部生物量。此外，根系构型明显改变，总根长增加但根数和根体积降低。生理测定显示，总叶绿素（total chlorophyll, Chl）含量、光化学效率（photochemical efficiency, F_v/F_m）和净光合速率（net photosynthetic rate, P_n）显著降低。同时，硝酸还原酶（nitrate reductase, NR）活性、NO₃^?浓度和氮含量显著降低，而可溶性糖含量显著升高。转录组分析共鉴定1185个差异表达基因（differential expression genes, DEGs）（585个上调、600个下调），涉及氮转运（N transport）、氮同化（N assimilation）、糖酵解（glycolysis）和三羧酸循环（tricarboxylic acid cycle, TCA）的基因发生显著变化。此外，7个不同家族的27个转录因子（transcription factor, TF）基因显著差异表达。基于低氮胁迫下的表达模式，鉴定出3个bHLH转录因子FaIRO2、FabHLH69和FabHLH062为高羊茅低氮响应的候选调控因子。综上，该研究提供了高羊茅响应低氮胁迫的整合生理和转录见解，并突出可用于遗传改良具有高氮利用效率的多年生牧草和草坪草的潜在靶基因。

氮（N）是植物必需的大量元素，也是氨基酸、核酸、叶绿素（chlorophyll, Chl）和多种次生代谢物的重要组分。土壤中NO₃^?易因淋溶和反硝化而损失，导致低氮（low nitrogen, LN）胁迫成为限制植物营养生长和生殖发育的重要因素。低氮会降低光合色素合成、光合电子传递链效率和碳固定关键酶活性，造成碳氮代谢失衡、叶色变淡、分蘖减少和生物量下降。以往低氮胁迫机制研究多集中于水稻、玉米、小麦等主粮作物，而牧草和草坪草研究相对有限；同时从地上部和根系两个角度，并结合碳氮代谢途径解析草类低氮响应的研究仍较少。高羊茅（Festuca arundinacea）是多年生冷季型草坪草和牧草，前期研究初步阐明其低氮生理响应，但转录调控网络仍不清楚，关键上游调控基因和相关分子途径有待探索。基于此，研究人员整合生理、生化和转录组分析，研究高羊茅地上部和根系对低氮胁迫的响应机制，重点关注碳氮代谢途径，以期为耐低氮牧草和草坪草育种提供理论依据和候选基因。

在关键技术方法方面，研究人员采用水培体系，以高羊茅‘K31’种子萌发并培养21 d的幼苗为材料，设置低氮（0.25 mM NO₃^?）与中等氮（moderate nitrogen, MN；4 mM NO₃^?）处理；处理35 d测定表型、根系构型和生理指标，处理14 d取叶片进行转录组测序（每处理3个生物学重复），并用qRT-PCR验证差异基因及0—48 h转录因子时序表达；根系经扫描和WinRHIZO分析，光合参数用LCproT和IMAGING-PAM测定。

研究结果如下。2.1 低氮胁迫对高羊茅地上部生长的影响：通过21日龄幼苗低氮与中等氮处理35 d后的表型测定发现，低氮显著抑制地上部生长，株高、分蘖数、叶片长度、叶鞘长度和地上部生物量分别较中等氮降低34.9%、56.6%、45.8%、61.6%和77.9%。2.2 低氮胁迫对高羊茅根系生长的影响：根系扫描分析表明，低氮促进根伸长但抑制根生物量，最长根长增加24.7%，根生物量降低26.9%；冠根数、表面积和体积分别降低56.1%、25.6%和50.9%，冠根长无显著变化；侧根数降低42.4%，侧根长和体积分别增加47.1%和28.3%，侧根表面积无显著变化；总根长增加68.7%，总根数和总根体积分别降低42.6%和40.0%。2.3 低氮胁迫对高羊茅光合性能的影响：总叶绿素含量、P_n和F_v/F_m分别显著降低60.4%、37.2%和31.6%。2.4 低氮胁迫对高羊茅氮代谢的影响：地上部NO₃^?浓度、氮含量和NR活性分别降低65.1%、50.7%和99.4%；根系分别降低80.5%、49.1%和98.4%。2.5 低氮胁迫对高羊茅可溶性糖含量的影响：地上部葡萄糖、蔗糖和果糖含量平均增加462.8%，根系平均增加约92.1%。2.6 差异表达基因分析：转录组测序获得超过4200万条clean reads，Q30比例超过97.18%；以|log₂ fold change|≥1且P_adj≤0.05筛选到1185个DEGs，其中585个上调、600个下调；主成分分析和Pearson相关分析显示处理内生物学重复高度相似、处理间明显分离。2.7 DEGs的GO和KEGG富集分析：GO富集涉及蛋白质折叠、碳固定和α-氨基酸生物合成等生物学过程，分子功能富集于磷酸烯醇式丙酮酸羧激酶活性、磷酸烯醇式丙酮酸羧化酶活性和蛋白质折叠伴侣等；KEGG富集前五位包括Calvin循环碳固定、丙氨酸/天冬氨酸/谷氨酸代谢、辅因子生物合成、类胡萝卜素生物合成和氮代谢。2.8 碳氮代谢相关基因对不同氮处理的表达特征：RNA-Seq分析显示，糖酵解和TCA循环相关基因中GAPDH1和PEPC1表达显著升高，ME、MDH、OGDH2、ALT、GAPDH2和PEPC2表达显著降低；氮转运和同化相关基因中NRT1、NRT2.4和GS表达上调，NR和NIR表达显著下调。2.9 低氮胁迫下转录因子表达：以|log₂ fold change|≥2且P_adj≤0.05筛选到7个家族的27个差异表达TF，其中bHLH家族最多（3个上调、6个下调），其次为HSF家族（4个上调、4个下调）。2.10 qRT-PCR验证低氮胁迫下TF时序表达模式：对9个bHLH和5个HSF基因在0、1、3、6、12和48 h的表达分析显示，bHLH基因可分为持续上调型（FaIRO2-G1、FaIRO2-P、FabHLH69-P）、不受低氮诱导型（FaHEC2-G1）、持续下调型（FabHLH062-G2、FabHLH062-P）、早期上调型（FaILR3-G1、FaILR3-P）和波动表达型（FabHLH69-G2）；HSF基因分为早期上调型（FaHSFA2c-P、FaHSFA2a-P、FaHSFA2a-G2、FaHSFB2c-G1）和不受低氮诱导型（FaHSFA2a-G1）。2.11 qRT-PCR验证DEGs：选取6个碳氮代谢基因和3个TF进行验证，FaGS、FaNRT2.4、FaNFYA7和FaERF1表达显著上调，FaNR、FaME、FaMDH、FaOGDH和FaLBD37表达下调；9个DEGs的qRT-PCR与RNA-Seq log₂ fold change高度相关（R²=0.8387），表明RNA-Seq数据可靠。

讨论部分总结认为，低氮胁迫严重抑制高羊茅地上部生长，根系则表现出根长增加等觅食反应；地上部氮含量和生物量与叶绿素含量、P_n和F_v/F_m呈强正相关，说明低氮导致的氮状态下降与光合能力降低密切耦合，进而限制生物量积累。转录组结果表明低氮诱导碳氮代谢相关基因的转录重编程：NRT1和NRT2.4上调提示低亲和与高亲和硝酸盐吸收系统被激活，NR和NIR下调、GS上调提示硝酸还原受抑而内源氮再循环增强；NR活性降低预期减少硝酸同化对碳骨架和还原力的需求，可能促进可溶性糖积累和碳氮比升高；糖酵解和TCA循环关键基因表达改变，反映碳代谢重排。差异表达TF中bHLH家族最丰富，FaIRO2、FabHLH69-P和FabHLH062因其表达模式被鉴定为候选低氮响应基因。作者提出低氮响应机制模型：根系感知低氮后，NO₃^?转运和同化基因上调；地上部光合相关基因下调并与叶绿素、P_n和F_v/F_m下降一致；NR活性降低减少碳骨架和还原力消耗，可能引起可溶性糖积累和碳氮比升高；碳盈余可能分配至根系，支持根系构型可塑性；NO₃^?转运和GS基因上调提示根系氮获取和内源氮再循环潜力增强，可能部分缓解氮饥饿并形成适应性反馈环。研究结论翻译如下：综上，低氮（LN）胁迫严重抑制高羊茅地上部生长、改变根系构型、降低光合性能和氮代谢能力，同时提高可溶性糖含量。转录组分析揭示低氮胁迫触发碳代谢和氮代谢途径的协同改变，具体影响氮转运、氮同化、糖酵解和三羧酸循环（TCA）相关基因。最后，研究人员聚焦低氮胁迫下差异表达转录因子（TF）的表达模式，并鉴定出3个可能参与高羊茅低氮胁迫响应的候选基因。综上，该研究提供了高羊茅响应低氮胁迫的整合生理和转录见解，并突出可用于遗传改良具有高氮利用效率的多年生牧草和草坪草的潜在靶基因。