整合多组学分析揭示苯丙烷介导的木质素沉积在增强烟草耐冷性中的作用

《Plant Stress》：Integrated Multi-omics Analysis Reveals the Role of Phenylpropanoid-Mediated Lignin Deposition in Enhancing Cold Tolerance in Tobacco

　　

冷胁迫显著制约作物产量和空间分布。尽管已报道了通用的冷响应途径，但烟草中的特定调控结构仍不清楚。研究人员对耐冷品种（YM2020）和冷敏感品种（YY99）的生理和多组学特征进行了系统比较。结果表明，YM2020表现出优越的耐冷性，表现为叶片萎蔫减轻、抗氧化活性增强和有效的渗透调节。转录组分析揭示了YM2020中冷响应转录因子（TF）家族的选择性上调。加权基因共表达网络分析（WGCNA）鉴定出一个与耐冷性相关的关键模块，突出显示WRKY58、EIN3和REVEILLE7为候选枢纽基因。整合分析表明，苯丙烷代谢相关通路是烟草响应冷胁迫的关键通路，其中关键基因（如PAL- Nta09g12930）和代谢物（如松柏醇）与上述枢纽转录因子显著相关。此外，YM2020表现出显著升高的木质素含量和苯丙氨酸解氨酶（PAL）活性，并通过增加的木质素沉积得到证实。本研究结果表明，烟草耐冷性受协调的转录因子（TF）介导的木质素生物合成途径增强所调控，为分子育种提供了靶点。

**论文解读：整合多组学分析揭示苯丙烷介导的木质素沉积在增强烟草耐冷性中的作用**

**研究背景**
冷胁迫是限制作物地理分布、生长与产量的主要非生物胁迫之一。烟草（*Nicotiana tabacum*）作为重要经济作物，在幼苗期对低温极为敏感，冷害导致生长抑制、品质下降及经济损失。不同烟草品种间耐冷性存在显著差异，但其分子调控机制尚不明确。尽管已有研究报道了冷胁迫的一般响应途径，但烟草中特定调控网络，尤其是苯丙烷代谢与木质素沉积在耐冷性中的角色，缺乏系统解析。因此，本研究旨在通过整合生理、转录组与代谢组方法，揭示烟草耐冷性的分子基础，为分子育种提供靶点。论文发表在《Plant Stress》。

**主要技术方法**
研究人员以耐冷品种YM2020和冷敏感品种YY99（来自河南农业大学烟草学院育种实验室）为材料，在6–8叶期进行4°C冷胁迫处理，分别于0小时和2小时采样。关键技术包括：非靶向代谢组学（UHPLC-Q Exactive HF-X系统，LC-MS/MS）分析差异积累代谢物（DAMs）；转录组测序（RNA-Seq，Illumina平台）鉴定差异表达基因（DEGs），并利用Plant Transcription Factor Database预测转录因子（TF）家族；加权基因共表达网络分析（WGCNA，R包）构建共表达模块并识别枢纽基因；整合Spearman相关性分析及KEGG通路富集；苯丙氨酸解氨酶（PAL）活性测定、木质素含量（重铬酸钾氧化滴定法）及组织染色（Wiesner反应）验证。

**研究结果**
**3.1 冷胁迫下烟草的表型与生理响应**
通过形态观察和生理指标测定，发现冷胁迫2小时后，YM2020叶片萎蔫程度显著轻于YY99。组织化学染色（DAB和NBT）显示YY99中积累更多超氧阴离子（O₂^?）和过氧化氢（H₂O₂）。定量分析表明，YM2020具有更低的相对电导率、丙二醛（MDA）和H₂O₂含量，更高的超氧化物歧化酶（SOD）、过氧化物酶（POD）和过氧化氢酶（CAT）活性，以及更多的可溶性糖和可溶性蛋白积累，证实其更强的抗氧化能力和渗透调节能力。

**3.2 代谢组与差异代谢物分析**
主成分分析（PCA）和偏最小二乘判别分析（PLS-DA）表明冷胁迫显著改变了两品种的代谢谱。YM2020中，能量代谢（丙酮酸代谢、氧化磷酸化、三羧酸循环）相关DAMs以上调为主，而YY99中以下调为主。在次生代谢方面，YM2020中黄酮和苯丙烷生物合成通路的DAMs显著上调，而YY99中这些通路普遍下调。抗氧化相关通路（如抗坏血酸代谢）在YM2020中上调，在YY99中下调。

**3.3 转录组与差异表达基因分析**
共有767,871,202个有效读段，比对到烟草参考基因组，Q20≥98.57%，唯一比对率91.33%–93.35%。差异表达分析共鉴定263个DEGs，以上调为主。GO富集显示DEGs显著富集于非生物胁迫应答、防御反应、质膜及激酶活性等条目。KEGG分析表明，苯丙烷生物合成、植物激素信号转导和黄酮生物合成通路在两品种中均显著富集，但YM2020特有硫代谢、玉米素生物合成等通路，YY99特有苯乙烯生物合成等通路。

**3.4 转录因子分析**
对35个转录因子家族的系统分析显示，在YM2020中，冷胁迫后ERF、WRKY、MYB和NAC家族成员上调数量显著多于YY99。例如，WRKY家族在YM2020中有49个成员上调，仅3个下调；而在YY99中仅3个上调，26个下调。MYB家族在YM2020中39个上调，YY99中仅6个上调。这些结果表明转录因子差异表达是耐冷性差异的重要分子基础。

**3.5 WGCNA与枢纽基因鉴定**
WGCNA将7,520个DEGs划分为12个共表达模块。MEblue模块与POD、CAT、可溶性糖、可溶性蛋白等生理指标显著正相关（|r|≥0.8），MEturquoise模块与相对电导率、H₂O₂显著负相关。在MEturquoise模块中，以WRKY58（Nta03g32330）和EIN3（Nta15g09050、Nta14g17750）为核心的调控网络被构建；在MEblue模块中，以REVEILLE7（RVE7）家族基因（Nta03g30760、Nta04g20210、Nta03g06750）为核心的网络被构建。qRT-PCR验证了这些枢纽基因的表达趋势与转录组数据一致。

**3.6 冷胁迫下转录组与代谢组的协同响应**
KEGG联合分析表明，苯丙氨酸代谢、苯丙氨酸/酪氨酸/色氨酸生物合成和苯丙烷生物合成是两品种共有的显著富集通路，且也是YM2020与YY99在冷胁迫下的差异通路。聚焦于这些通路，发现在YM2020中，木质素生物合成关键基因（PAL、C4H、4CL、C3H、COMT、CCR、CAD、POD）均显著上调，相关代谢物（L-苯丙氨酸、肉桂酸、对香豆酸、松柏醇）也呈上升趋势。Mantel检验显示PAL与松柏醇显著正相关，COMT与咖啡酰辅酶A正相关等。相关性网络表明，RVE7与木质素合成基因和代谢物高度正相关，而EIN3和WRKY58仅与L-苯丙氨酸负相关。木质素含量测定和PAL活性测定显示，冷胁迫后YM2020的木质素含量和PAL活性显著高于YY99，组织染色（间苯三酚-HCl）证实YM2020茎基部维管束周围木质素沉积更广泛、颜色更深。

**讨论与结论**
讨论部分总结了：冷胁迫下烟草的生理生化响应、转录因子及枢纽基因的作用、以及苯丙烷代谢介导的木质素沉积的关键角色。研究人员提出，RVE7可能通过直接或间接调控*PAL*基因表达驱动木质素积累，而EIN3和WRKY58可能参与其他调控分支。该模型强调耐冷性涉及多个防御层的协同（转录调控、苯丙烷代谢激活、抗氧化系统增强、渗透调节物质积累）。研究局限性包括酶活性测定不全面、代谢物定量不足、以及枢纽基因与木质素基因的直接调控关系尚需分子实验验证（如酵母单杂交、双荧光素酶报告、染色质免疫共沉淀）。

**结论翻译**：本研究通过整合生理、转录组和代谢组数据集，阐明了烟草耐冷性的差异机制。研究人员的研究表明，烟草品种YM2020的优越耐冷性源于木质素生物合成途径的协同增强。低温胁迫下，生理评估显示YM2020比YY99维持更好的细胞稳态。转录组分析表明耐冷品种中多个冷响应转录因子（TF）家族选择性上调。通过WGCNA，研究人员鉴定出与耐冷性相关的特定模块，并确定WRKY58、EIN3和REVEILLE7为候选枢纽调控因子。整合组学分析将苯丙烷代谢确定为主要代谢响应，其中关键结构节点（特别是PAL（Nta09g12930）和松柏醇）与这些枢纽转录因子表现出显著协调。后续的木质素特异性染色和酶学测定证实，YM2020在早期胁迫阶段保持较高的PAL活性和木质素积累。总之，这些发现表明烟草耐冷性由枢纽转录因子（TF）调控的木质素生物合成途径增强所介导。本研究为烟草耐冷性机制及木质素在此过程中的功能作用提供了新的理论见解。