《Life》：Transcriptomic and Metabolomic Analyses Reveal the Phenylpropanoid Metabolism Regulatory Network in Arabidopsis Response to Colletotrichum higginsianum Infection
编辑推荐：
炭疽菌（Colletotrichum higginsianum, Ch）是一种典型的半活体营养型子囊菌真菌，其引起的病害常导致全球十字花科作物生产遭受重大经济损失。然而，目前的知识仍不足以使研究人员更深入地理解宿主植物在Ch侵染过程中转录水平上的响应机制。本研
炭疽菌（Colletotrichum higginsianum, Ch）是一种典型的半活体营养型子囊菌真菌，其引起的病害常导致全球十字花科作物生产遭受重大经济损失。然而，目前的知识仍不足以使研究人员更深入地理解宿主植物在Ch侵染过程中转录水平上的响应机制。本研究中，研究人员对Mock组和Ch侵染组样本进行了转录组学和代谢组学分析。结果表明，Ch侵染显著抑制了宿主植株的地上部鲜重和初生根长度。此外，与苯丙烷代谢相关的基因本体（gene ontology, GO）术语和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）通路高度富集。研究鉴定到防御激素水杨酸（salicylic acid, SA）及三种属于苯丙烷代谢通路的代谢物。同时，研究人员筛选出5个相关的酶编码基因和78个转录因子（transcription factors, TFs）。其中，5个WRKY、3个MYB和2个NAC转录因子表现出显著的表达变化，并与酶基因高度相关。该研究结果丰富了十字花科植物病原响应的调控网络，为后续解析Ch抗性机制提供了潜在的分子候选靶标。
论文解读：苯丙烷代谢在拟南芥应答炭疽菌侵染中的调控网络
研究背景与问题：炭疽菌属（Colletotrichum）是一类具有高度多样性和强致病性的丝状真菌植物病原菌，其分类体系随研究技术发展不断更新，目前包含280个种、16个主要种复合群和15个单种系。其中，Colletotrichum higginsianum（Ch）是典型的半活体营养型子囊菌，可侵染包括重要芸薹属作物和模式植物拟南芥在内的多种十字花科植物，造成炭疽病，导致全球十字花科作物生产蒙受重大经济损失。Ch的侵染过程包括分生孢子附着、萌发、形成附着胞、穿透寄主角质层和细胞壁、进入活体营养阶段并最终转为死体营养阶段，引起宿主组织死亡和病斑形成。植物通过模式触发免疫（PTI）和效应子触发免疫（ETI）抵御病原，并由水杨酸（SA）、茉莉酸（JA）和乙烯（ET）等激素信号协调。苯丙烷代谢途径是植物次生代谢的核心枢纽，可将苯丙氨酸转化为800余种代谢物，其产物如木质素、类黄酮和酚酸在植物防御中发挥关键作用。然而，目前对于十字花科植物在Ch侵染过程中苯丙烷代谢途径在转录和代谢水平的系统 reprogramming 仍缺乏全面理解。以往研究多集中于单一层面，如特定致病基因、单一防御组分或小RNA动态，而整合转录组和代谢组数据专门解析苯丙烷代谢调控网络的研究仍不充分。因此，开展本研究具有必要性。
研究内容与结论：研究人员以拟南芥Col-0为材料，采用根接种Ch的方法，在侵染后14天（dpi）进行转录组（RNA-seq）和广泛靶向代谢组学分析。结果显示，Ch侵染显著抑制拟南芥地上部鲜重和初生根长度。转录组分析共筛选到2202个差异表达基因（DEGs），其中上调1023个、下调1179个。GO富集分析显示，生物过程类别中富集了“对真菌的响应”、“对几丁质的响应”、“防御反应”、“对SA的响应”等条目；细胞组分中富集于“质外体”、“细胞壁”、“胞外区”等；分子功能中富集于“蛋白激酶活性”、“氧化还原酶活性”、“过氧化物酶活性”等。KEGG富集分析显示“植物-病原互作”、“苯丙烷生物合成”、“MAPK信号通路-植物”、“谷胱甘肽代谢”等通路显著富集。代谢组分析共鉴定到834个差异积累代谢物（DAMs），其中上调564个、下调270个，主要包括香豆素、酮、醛、酯、萜类、糖类、生物碱、氨基酸、类黄酮和多酚等。整合分析发现苯丙烷生物合成通路中的酪氨酸、松柏醛和芥子碱三种代谢物在Ch侵染后发生显著变化，SA含量显著升高（log2FC=16.97）。进一步分析发现，苯丙烷途径中参与芥子碱合成的关键结构基因如AT4G36220（阿魏酸5-羟化酶1，F5H1）、AT1G77520和AT1G77530（O-甲基转移酶基因）显著上调，而最终产物芥子碱含量却下降，同时松柏醛大量积累，提示代谢流可能被重新导向木质素前体合成。转录因子分析共鉴定到78个差异表达转录因子（DE-TFs），主要属于WRKY、MYB和NAC家族，其中5个WRKY、3个MYB和2个NAC与酶基因高度相关。相关性网络分析显示这些转录因子可能通过调控苯丙烷生物合成基因的表达，改变碳分配，促进松柏醛积累，从而参与植物免疫防御。
主要关键技术方法：研究采用了以下主要技术方法：（1）转录组测序（RNA-seq）：利用Illumina NovaSeq 6000平台进行双端150 bp测序，经fastp质控、HISAT2比对到拟南芥参考基因组TAIR10，DESeq筛选差异表达基因（DEGs）；（2）广泛靶向代谢组学分析：采用UPLC-ESI-MS/MS系统（Waters Acquity I-Class PLUS，Applied Biosystems QTRAP 6500+）检测3000种预设代谢物，结合OPLS-DA模型和VIP值筛选差异积累代谢物（DAMs）；（3）基因表达验证：采用qRT-PCR对15个候选基因进行表达验证；（4）相关性网络分析：基于Pearson相关系数（FDR<0.05，|r|>0.8）构建基因表达与代谢物丰度之间的相关网络，并识别枢纽节点；（5）真菌鉴定：通过四基因座（ITS、TUB2、ACT、GAPDH）系统发育分析结合形态学特征将菌株鉴定为Ch。样本来源为韶关大学田间病萝卜上分离的菌株，拟南芥种子为Col-0野生型。
研究结果分述：
3.1 Ch侵染根部抑制拟南芥生长：通过含菌丝琼脂块接种拟南芥根部，14 dpi观察表型发现Ch侵染幼苗叶片卷曲、根变短；定量数据显示Mock组地上部鲜重和初生根长度均显著高于Ch组，表明Ch侵染导致拟南芥生长受抑。
3.2 转录组数据概览：RNA-seq数据质量高（Q30>96.87%，比对率>98.28%），PCA显示PC1和PC2共解释64.8%变异，组内相关性高（>0.992）。以fold change≥2且FDR<0.01筛选到2202个DEGs（上调1023个，下调1179个），表达模式聚类为Mock高表达/Ch低表达和Ch高表达/Mock低表达两类。
3.3 DEGs的GO和KEGG富集分析：GO分析显示DEGs富集于防御反应、SA响应、细胞壁、质外体、过氧化物酶活性等条目；KEGG分析显示上调通路包括“植物-病原互作”、“苯丙烷生物合成”、“MAPK信号通路-植物”等，下调通路包括“植物激素信号转导”、“类黄酮生物合成”、“花青素生物合成”、“角质、木栓质和蜡生物合成”等。
3.4 GO和KEGG富集分离的DEGs分析：去除交集后获得94个DEGs，再以|log2FC|>5筛选得到22个DEGs，其中除AT2G33790（AGP30）下调外，其余均被Ch诱导上调，包括RmlC样cupin基因、植物防御素基因、CRK受体基因、VANGUARD样蛋白基因、果胶修饰抑制蛋白基因、LCR67、过氧化物酶超家族蛋白基因、谷胱甘肽转移酶tau 3、几丁质酶家族蛋白基因和蛋白激酶家族蛋白基因等。
3.5 DEGs表达验证：通过qRT-PCR检测15个候选基因的相对表达量，结果显示其表达趋势与转录组FPKM数据一致，验证了转录组数据的可靠性。
3.6 代谢组比较分析：PCA显示PC1和PC2累计解释86.11%变异，Ch组与Mock组明显分离。共鉴定到834个DAMs（上调564个，下调270个），主要类别包括香豆素、酮、醛、酯、萜类、糖类、生物碱、氨基酸、类黄酮和多酚。KEGG富集显示DAMs显著富集于ABC转运蛋白、多种植物次生代谢物生物合成、苯丙烷生物合成、苯丙氨酸/酪氨酸/色氨酸生物合成、淀粉和蔗糖代谢、半乳糖代谢、氨基酸代谢和植物激素信号转导等通路。
3.7 转录组与代谢组整合分析：结合DEGs和DAMs的KEGG富集结果，聚焦苯丙烷生物合成通路。发现酪氨酸、松柏醛和芥子碱三种DAMs在Ch侵染后分别呈现上调和下调变化，SA含量显著升高。转录组显示F5H1（AT4G36220）和两个O-甲基转移酶基因（AT1G77520、AT1G77530）显著上调，两个NAD(P)结合Rossmann折叠超家族蛋白基因（AT1G09480、AT1G76470）表达增加。尽管转录水平上芥子碱合成基因被激活，但终产物芥子碱显著减少，同时松柏醛大量积累，表明代谢流可能被重新导向松柏醛/木质素合成方向。
3.8 转录因子分析：共鉴定到78个DE-TFs，主要属于WRKY、MYB、NAC家族，其中WRKY家族以上调为主，MYB家族多数下调。TF-代谢物相关网络显示许多TF与酪氨酸和松柏醛显著正相关，与芥子碱负相关。进一步选取|log2FC|大于3的10个DE-TFs与5个酶编码基因进行相关性分析，鉴定出与芥子碱和松柏醛生物合成关键结构基因（AT4G36220、AT1G77520、AT1G77530）显著相关的枢纽TF。结果表明TF可能通过调控苯丙烷生物合成基因的转录来重塑次生代谢中的碳分配。
讨论与结论：讨论部分指出，Ch侵染显著重编程宿主转录组和代谢组，激活防御反应相关基因，同时抑制生长相关进程，体现了生长-防御权衡。苯丙烷代谢在Ch侵染中发挥核心作用，SA含量升高与苯丙烷途径激活有关。松柏醛作为木质素前体的积累和芥子碱的减少提示代谢流可能被重定向至木质素合成，但该假说需进一步实验验证。转录因子（WRKY、MYB、NAC、bZIP、C2H2）可能作为上游开关调控苯丙烷代谢网络。结论总结：本研究通过整合转录组和代谢组分析，证实了苯丙烷代谢在Ch侵染中的关键作用；酪氨酸、松柏醛和芥子碱的差异积累可能与Ch侵染和宿主生长抑制有关；筛选到5个酶编码基因作为潜在抗性改良靶点；78个TFs（属于WRKY、MYB、NAC、C2H2和bZIP家族）中多个TF与苯丙烷途径酶基因呈显著正/负相关，可能参与调控苯丙烷代谢。该研究为植物-病原互作提供了有价值的遗传资源和分子靶标。研究局限性包括：仅采用单一时间点采样，无法捕捉动态变化；生长-防御权衡假说未经验证；同一培养皿中多株幼苗共享培养基和接种源，可能存在平板水平混杂效应。
