光温耦合调控雪茄烟草外包皮叶次生脉木质化的代谢与转录重编程

《Environmental and Experimental Botany》：Combined light and temperature regimes coordinate metabolic and transcriptional reprogramming underlying secondary vein lignification in cigar tobacco wrapper leaves

　　

光温组合协调雪茄烟草外包皮叶次生脉木质化的代谢与转录重编程。光温耦合调控次生脉木质化；低光高温（LLHT）抑制苯丙烷代谢通量及单木质醇积累；在LLHT下，木质化降低与发育抑制共同限制次生脉直径；LLHT下的激素变化与维管分化及木质素沉积受抑相关；拮抗基因模块编码木质化的环境调控。研究人员通过整合次生脉解剖结构、细胞壁组成、苯丙烷代谢物、激素谱、时间分辨转录组学及加权基因共表达网络分析（WGCNA），揭示了光温条件调控雪茄烟草（Nicotiana tabacum L.）外包皮叶次生脉木质化的分子与代谢机制。研究发现，低光高温（LLHT）条件下，次生脉直径、维管束面积及木质素含量显著降低（分别较HLHT减少11.5%、20.9%和22.3%），苯丙烷代谢通量受抑，单木质醇前体（如p-香豆酸、松柏醇、p-香豆醇）积累减少。激素分析显示，LLHT下细胞分裂素（如tZ、tZR、IP）水平降低，乙烯前体ACC积累增加，ABA在中间阶段短暂升高，茉莉酸信号减弱。转录组分析表明，LLHT诱导了阶段特异性转录重编程，差异表达基因（DEGs）在21天时大幅增加（2343个上调，1197个下调）。WGCNA鉴定出两个与木质素含量高度相关的模块：Honeydew模块（正相关，r=0.95）在LLHT下下调，包含蔗糖合酶（NtaSUS）、咖啡酸O-甲基转移酶（NtaCOMT）及细胞分裂素响应调节因子（NtaARR-A）等枢纽基因；Greenyellow模块（负相关，r=-0.92）在LLHT下上调，包含乙烯响应因子（NtaEREBP）、ABA受体（NtaPYL）、谷氧还蛋白（NtaGrxC）等应激信号基因。研究结论：光温联合调控叶片维管木质化，低光高温驱动一种低木质素发育状态，该状态通过苯丙烷代谢抑制、激素重编程及拮抗基因网络调控实现，低温可缓解低光的抑制效应，为优化雪茄外包皮叶质量提供了理论框架。

**论文解读：光温耦合调控雪茄烟草外包皮叶次生脉木质化的机制研究**

**一、研究背景与问题**

雪茄烟草（*Nicotiana tabacum* L.）的外包皮叶是决定雪茄品质的关键部位，其质量由叶片大小、表面均匀性、弹性、燃烧性及叶脉特征等复杂性状综合决定。与卷烟用烟叶不同，雪茄外包皮叶要求叶片大而薄、柔韧且结构均匀，这对叶片发育和环境管理提出了独特要求。在传统和现代雪茄生产中，适度遮荫（降低光照强度）被广泛用于促进叶片扩展、减少厚度并改善表面光滑度，从而提升包裹性能。然而，光照降低同时限制碳同化并改变发育信号，可能影响叶片内部结构，尤其是维管系统的发育。次生脉（secondary vein）的木质化程度直接影响叶片的机械性能和视觉外观：过粗的次生脉会降低弹性、破坏表面均匀性并影响卷制和发酵过程。木质素（lignin）作为次生细胞壁的主要结构聚合物，赋予刚性和机械强度，但过度沉积会损害叶片柔韧性。木质素生物合成依赖于苯丙烷代谢（phenylpropanoid metabolism），该途径对环境因子（如光、温）敏感，并与发育和激素信号紧密协调。尽管对茎秆和木质组织中木质素调控的研究已较深入，但叶片维管木质化的环境调控机制，尤其在雪茄烟草中，仍缺乏系统研究。此外，叶片维管发育受生长素（auxin）、细胞分裂素（cytokinin）、乙烯（ethylene）和茉莉酸（jasmonate）等激素网络调控，但光温条件如何整合这些信号以重塑次生脉木质化，尚不清楚。因此，本研究旨在阐明光温组合调控雪茄烟草次生脉木质化的分子与代谢机制，为优化叶片结构性状和质量提供理论依据。

**二、研究内容与结论**

研究人员在青岛即墨实验基地（中国农业科学院烟草研究所）开展实验，以雪茄烟草品种HY103为材料，设置四种光温处理：高光高温（HLHT，600±30 μmol·m?2·s?1，30±2°C）、低光高温（LLHT，300±30 μmol·m?2·s?1，30±2°C）、高光低温（HLLT，600±30 μmol·m?2·s?1，20±2°C）、低光低温（LLLT，300±30 μmol·m?2·s?1，20±2°C）。所有测量均针对第10叶进行。研究得出以下结论：（1）低光高温（LLHT）显著抑制叶片生长，表现为叶干物质含量（LDMC）和收获指数（HI）最低，次生脉直径较HLHT减少11.5%，维管束面积减少20.9%，木质素含量降低22.3%。（2）LLHT下苯丙烷代谢通量受抑，单木质醇前体（p-香豆酸、p-香豆醇、松柏醇）积累减少，而低温（LLLT）可部分缓解低光的抑制效应。（3）激素分析显示，LLHT下细胞分裂素（tZ、tZR、IP）水平降低，乙烯前体ACC积累增加，ABA在14天短暂升高，JA-Ile水平下降，表明LLHT创建了不利于次生壁沉积的激素信号环境。（4）转录组分析揭示LLHT诱导了阶段特异性转录重编程，差异表达基因（DEGs）在21天时大幅增加（2343个上调，1197个下调），且苯丙烷生物合成途径基因在多个节点（PAL、C4H、4CL、HCT、CCR、CAD、COMT等）受到广泛抑制。（5）加权基因共表达网络分析（WGCNA）鉴定了两个与木质素含量高度相关的模块：Honeydew模块（正相关，r=0.95）包含蔗糖合酶（NtaSUS）、咖啡酸O-甲基转移酶（NtaCOMT）及细胞分裂素响应调节因子（NtaARR-A）等枢纽基因，在LLHT下表达下调；Greenyellow模块（负相关，r=-0.92）包含乙烯响应因子（NtaEREBP）、ABA受体（NtaPYL）、谷氧还蛋白（NtaGrxC）等应激信号基因，在LLHT下表达上调，表明木质化可塑性由两个拮抗转录网络编码。该研究为理解光温信号如何通过代谢与调控网络整合以塑造叶片维管发育和木质化提供了重要机制见解，对优化雪茄外包皮叶的受控栽培具有实践指导意义。论文发表在《Environmental and Experimental Botany》。

**三、主要技术方法**

本研究采用以下关键技术方法：（1）可控环境生长室光温处理：使用红（660 nm）、蓝（450 nm）、绿（525 nm）LED光源，保持红蓝绿比（R:B:G）为7:1:1，设置四种光温组合，约40株/处理。（2）次生脉解剖学分析：21天时取样，石蜡包埋、甲苯胺蓝O染色，光学显微镜观察并测量维管束面积和次生脉直径（ImageJ）。（3）细胞壁组分定量：依据中国农业行业标准NY/T 3494-2019，通过两步硫酸水解法测定纤维素、半纤维素和木质素含量，酸溶木质素用分光光度法（320 nm）测定，酸不溶木质素用重量法（575°C灰化）测定。（4）靶向代谢物与激素分析：采用超高效液相色谱-串联质谱（UPLC-MS/MS）多反应监测（MRM）模式，检测苯丙烷途径中间体（p-香豆酸、咖啡酸、阿魏酸、p-香豆醇、松柏醇）及激素（生长素、细胞分裂素、茉莉酸、水杨酸、脱落酸、乙烯前体ACC）。（5）RNA提取与转录组测序：7、14、21天时提取次生脉总RNA，Illumina测序，比对至雪茄烟草参考基因组，DESeq2鉴定差异表达基因（DEGs，FDR<0.05，|log?FC|≥1），GO和KEGG富集分析。（6）加权基因共表达网络分析（WGCNA）：构建共表达网络，鉴定与木质素及代谢物性状显著相关的模块，提取枢纽基因。（7）qRT-PCR验证：对WGCNA枢纽基因进行定量PCR验证，以Actin为内参，2?ΔΔCt法计算相对表达量。

**四、研究结果**

**3.1. 光温处理对叶片形态和生物量的影响**
通过21天生长观察，LLHT植株表现出更高的株高和更薄的叶片，而HLLT植株矮小、叶片更厚更绿。叶长在LLHT和LLLT下最长；鲜重处理间无显著差异，但干重在LLHT和HLHT下显著低于HLLT和LLLT。总叶生物量和茎生物量在低光处理（LLHT、LLLT）下显著高于高光处理（HLHT、HLLT）。叶干物质含量（LDMC）在LLHT下最低，HLLT下最高；收获指数（HI）在LLHT下最低，LLLT下最高，表明低光高温下生物量分配效率受损。

**3.2. 光温对次生脉解剖结构、细胞壁组成及木质素相关代谢物的影响**
次生脉横切面观察显示，LLHT下维管束面积显著小于其他处理，次生脉直径较HLHT、HLLT、LLLT分别减少11.5%、4.7%和2.9%。细胞壁组分定量表明，LLHT下纤维素含量相对较高，半纤维素含量下降，木质素含量降低最为显著（较HLHT减少22.3%，较HLLT减少26.4%，较LLLT减少23.4%），木质素含量与次生脉直径及维管束面积呈平行变化。

**3.3. 木质素含量与次生脉直径紧密相关**
木质素含量、次生脉直径和维管束面积三者之间的对应关系表明，LLHT下次生脉增厚受限由木质化受损和维管组织发育抑制共同驱动，而非细胞壁多糖的整体减少。

**3.4. LLHT下单木质醇积累受抑反映苯丙烷输出减少**
21天时，苯丙烷途径底物L-苯丙氨酸（Phe）在HLLT下最低，LLHT下与HLHT相当；但下游单木质醇前体（p-香豆酸、p-香豆醇、松柏醇）在LLHT下均显著降低，表明途径通量受限。

**3.5. 光温处理重塑激素谱与调控代谢物**
差异代谢物分析显示，LLHT vs HLHT比较中，独特差异代谢物在14天最多（11个），21天消失；LLLT vs HLLT比较中，独特差异代谢物在三个时间点均有分布。激素定量分析表明，LLHT下细胞分裂素（tZ、tZR、IP）水平降低，乙烯前体ACC在14天显著升高，JA-Ile在后期降低，ABA在14天升高。生长素响应基因表达复杂，但细胞分裂素代谢物持续降低。低温处理（LLLT）下，上述激素变化被部分缓解。

**3.6. 光温组合处理下次生脉的转录组重编程**
主成分分析（PCA）显示样品按处理和时间点分离，LLHT在14和21天与HLHT明显分离。差异表达分析表明，LLHT vs HLHT的DEGs从7天（112个上调，238个下调）急剧增加至21天（2343个上调，1197个下调），且不同时间点DEGs重叠有限。LLLT vs HLLT的DEGs在所有时间点均较少。KEGG富集分析中，苯丙烷生物合成途径虽未达到多重检验校正后的显著性，但富集因子较高（~1.98），提示该途径参与。

**3.7. 整合转录组-代谢组分析揭示LLHT下苯丙烷通量的协同抑制**
将转录组和代谢组数据联合映射到苯丙烷生物合成途径，发现LLHT下从途径入口（PAL）到下游单木质醇合成酶（CCR、CAD、COMT）等基因均广泛下调，且p-香豆酸、p-香豆醇、松柏醇在21天时积累减少，表明途径通量在多个节点受到协同抑制。

**3.8. WGCNA鉴定木质素相关基因模块与枢纽网络**
WGCNA分析将基因聚类为不同模块，其中Honeydew模块与木质素含量呈强正相关（r=0.95），包含蔗糖合酶（NtaSUS）、咖啡酸O-甲基转移酶（NtaCOMT）、细胞分裂素响应调节因子（NtaARR-A）等枢纽基因，在LLHT下表达下调；Greenyellow模块与木质素含量呈强负相关（r=-0.92），包含乙烯响应因子（NtaEREBP）、ABA受体（NtaPYL）、谷氧还蛋白（NtaGrxC）等信号相关基因，在LLHT下表达上调，表明木质化可塑性由两个拮抗模块控制。

**3.9. qRT-PCR验证选定的WGCNA枢纽基因**
qRT-PCR验证了Honeydew模块基因（NtaCOMT、NtaSUS、NtaARR-A、NtaGAE）在LLHT下表达低于其他处理，而Greenyellow模块基因（NtaMYB、NtaGALS、NtaPYL）在LLHT下表达较高，但NtaEREBP和NtaGrxC在LLLT下最高，NtaTMK在HLLT下最高，整体上Honeydew模块基因的抑制模式在LLHT下得到一致验证。

**五、讨论与结论**

**讨论部分**：
（1）光温组合是次生脉发育的关键限制因素：光照强度是次生脉木质化的主导因子，温度调节光限制的幅度。LLHT下叶干物质含量、次生脉直径和木质素含量均最低，反映次生脉发育需要充足的碳供应和维管分化代谢能力。LLLT虽与LLHT光照相同，但低温下呼吸碳损失减少，部分维持了结构组织的碳分配。
（2）LLHT下次生脉直径受限与木质化降低和维管发育抑制相关：木质素含量与次生脉直径紧密对应，但维管束面积的减少也表明细胞组织发育差异贡献了脉径减小。半纤维素在LLHT下降低，进一步反映次生壁形成的整体受限。
（3）LLHT下苯丙烷通量抑制为木质化降低提供生化基础：LLHT降低了苯丙烷中间体和单木质醇丰度，低温可缓解低光对途径输出的抑制。
（4）LLHT下激素重编程主动限制维管分化：细胞分裂素降低、乙烯信号增强、ABA短暂升高、茉莉酸活动减弱，共同构成不利于次生壁沉积的信号环境。高温放大了激素转录重编程（24个激素靶基因差异表达，低温仅6个）。
（5）苯丙烷代谢作为生物学响应途径：尽管KEGG富集未达统计显著性，但该途径富集因子高，且DEGs包含经典木质素生物合成酶，从入口到分支点均受抑制，结合代谢物数据支持其在LLHT下的功能重要性。
（6）整合转录组-代谢组分析揭示LLHT下木质化的协同抑制：转录抑制从苯丙氨酸入口延伸至木质素分支酶，伴随单木质醇前体减少，表明途径水平抑制是累积性木质素沉积降低的基础。
（7）对比基因网络编码木质化可塑性的环境调控：Honeydew模块（正相关）促进碳分配和次生壁形成，Greenyellow模块（负相关）主导应激信号，两个模块的差异激活解释了LLHT下的低木质素状态。
（8）对雪茄外包皮叶质量和叶片功能生物学的影响：次生脉结构决定叶片机械性能和加工性能，LLHT作为高度限制条件对受控栽培有指导意义。
（9）局限性与未来方向：需功能验证优先调控因子，并结合酶活性测定和碳状态测量；LLHT下木质素减少可能源于发育延迟而非永久抑制，需后期时间点采样区分。

**结论部分翻译**：
光与温联合调控叶片维管木质化，低光高温驱动一种低木质素发育状态。研究人员表明，21天时次生脉直径减小与该发育阶段木质素积累减少相关，但这是否反映永久抑制或发育延迟的木质化，需通过后期时间点采样进一步研究。这些变化由苯丙烷代谢的协同抑制、激素重编程及拮抗基因网络调控支撑。重要的是，低温缓解了低光的抑制效应，突显温度是木质化可塑性的关键调节因子。通过整合解剖学、代谢组学、转录组学及网络分析，这项工作建立了一个连接环境信号与维管结构的多层次调控框架，为叶片发育可塑性提供机制见解，并为受控栽培下优化外包皮叶质量提供实践指导。