《Plant Physiology and Biochemistry》：Integrated Transcriptomic and Metabolomic Analyses Reveal Hormone-Mediated Crosstalk During Potato Virus Y and Potato Spindle Tuber Viroid Co-Infection
编辑推荐：
本文聚焦植物在面对马铃薯Y病毒（PVY）与马铃薯纺锤块茎类病毒（PSTVd）共侵染时的复杂分子响应机制。研究人员通过整合田间调查、转录组学和代谢组学分析，揭示了PVY在共侵染中的主导作用及其对宿主免疫和代谢通路的深刻重编程。研究发现，PVY驱动广泛的转录和代谢变化，尤其是细胞分裂素和茉莉酸途径的显著扰动，而PSTVd的影响则相对局限。该研究为理解病毒与类病毒之间的相互作用、宿主防御网络的整合以及制定针对复合病原体侵染的病害管理策略提供了新的分子见解。
在农业生产中，作物常常面临多种病原体的“组团”攻击，这使得病害防控变得异常复杂。想象一下，一株马铃薯可能同时遭受病毒和类病毒的侵扰，它们如何“狼狈为奸”或“相互竞争”，最终导致作物减产甚至绝收？这背后的分子“暗战”一直是科学家们想要破解的谜题。特别是马铃薯Y病毒（Potato virus Y， PVY）和马铃薯纺锤块茎类病毒（Potato spindle tuber viroid， PSTVd），两者都是极具破坏性的病原体，PVY是编码蛋白质的RNA病毒，而PSTVd则是更简单的、不编码蛋白质的环状RNA分子。当它们同时入侵植物时，植物细胞内部会发生怎样的“化学战”和“基因调控战”？宿主又是如何整合来自不同病原体的信号，做出最终的防御决策？为了回答这些问题，一项发表在《Plant Physiology and Biochemistry》上的研究，采用了多组学整合分析的方法，深入揭示了PVY与PSTVd共侵染下，植物宿主内部激素介导的复杂分子对话。
为了系统揭示PVY与PSTVd共侵染的分子机制，研究人员综合运用了多种关键技术。首先，他们在重庆马铃薯主产区开展了田间调查，采集了256份具有症状的马铃薯叶片样本，通过逆转录聚合酶链式反应（Reverse Transcription PCR， RT-PCR）系统筛查了包括PVY和PSTVd在内的多种病毒和类病毒的流行情况。其次，在可控环境下，利用本氏烟（Nicotiana benthamiana）和马铃薯“克新18号”品种构建了模拟单感染（PVY或PSTVd）与共感染（PVY-PSTVd）的实验体系，进行了系统的表型观察和病原积累量定量分析。研究的核心是组学分析：通过RNA测序（RNA-seq）技术对感染21天后的本氏烟叶片进行转录组测序，利用加权基因共表达网络分析（Weighted Gene Co-expression Network Analysis， WGCNA）挖掘与不同感染状态相关的基因模块；同时，采用超高效液相色谱-质谱联用（UPLC-MS）技术对相同样本进行广泛靶向代谢组学分析，重点关注激素类代谢物的变化。最后，通过整合转录组与代谢组数据，构建基因-代谢物关联网络，并利用实时荧光定量PCR（RT-qPCR）对关键差异表达基因进行验证。
研究结果层层递进，揭示了从田间现象到分子机制的完整故事：
3.1. 田间调查揭示PVY高流行且常与PSTVd共感染
对重庆地区256份马铃薯样本的检测发现，PVY是田间最流行的病原体，而PSTVd极少单独发生，却频繁与PVY共同出现，尤其是在多重感染组合中。这提示在自然条件下，PSTVd可能依赖于与PVY的共感染来维持其存在或传播，两者之间存在潜在的生物学关联。
3.2. 表型和分子响应显示PVY在共感染中占主导地位
在本氏烟和马铃薯中的实验表明，PVY单感染引起严重的植株矮化、叶片黄化和卷曲，而PSTVd单感染症状轻微。共感染植株的表型与PVY单感染几乎无法区分，PSTVd的症状被完全掩盖。分子水平上，共感染时PSTVd的积累量反而比单感染时略有下降，而PVY的积累则维持在高水平。这证实了PVY在共感染中无论是表型还是分子水平上都占据主导地位，甚至可能对PSTVd的复制产生一定的竞争或抑制。
3.3. 转录组分析揭示PVY与PSTVd引发截然不同的基因表达模式
RNA-seq分析显示，PVY感染（无论单感染还是共感染）引发了宿主基因组的剧烈重编程，差异表达基因（Differentially Expressed Genes， DEGs）数量远多于PSTVd单感染。主成分分析（Principal Component Analysis， PCA）清晰地将不同处理组分开，共感染的转录谱更接近于PVY单感染，但仍具有独特性。
3.4. 通路富集分析揭示主要的转录组变化由PVY和PSTVd单感染及共感染驱动
基因本体（Gene Ontology， GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes， KEGG）通路富集分析显示，PSTVd单感染主要富集于蛋白磷酸化、激酶活性以及与植物-病原互作相关的通路。而PVY单感染和共感染则强烈富集于光合作用、碳代谢、次级代谢物生物合成等全局性代谢通路，以及植物激素信号转导通路。这表明PSTVd的影响更侧重于特定的信号调控，而PVY则对宿主的整体代谢和能量平衡造成了广泛冲击。
3.5. 参与RNA沉默、免疫、细胞壁重塑和光合作用宿主基因的差异表达
具体到功能基因，PVY强烈上调了RNA沉默通路的核心组件（如AGO1， AGO2， RDR1， RDR6）和免疫相关基因（如PR1， EDS1， PAD4）。PSTVd对这些通路的激活作用较弱。两者都抑制了光合作用相关基因，但PVY的抑制作用更强。在细胞壁重塑和次级代谢方面，PVY和共感染引起了更广泛的基因上调。
3.6. 转录因子介导的调控响应和PSTVd与PVY感染下的通路富集
转录因子（Transcription Factor， TF）分析发现，PVY激活了包括WRKY、NAC、MYB、bHLH在内的多个TF家族的大量成员，而PSTVd激活的TF数量较少且更特异。共感染条件下，PVY的主导作用压制了部分PSTVd诱导的TF，但同时协同激活了与茉莉酸（Jasmonic Acid， JA）和赤霉素（Gibberellin， GA）信号相关的TF（如TIFY10A， SCL3），展现了复杂的转录调控交叉对话。
3.7. 共表达网络分析揭示响应PSTVd、PVY和共感染的感染特异性基因模块
WGCNA分析鉴定出与不同感染状态特异性相关的基因共表达模块。例如，棕色模块与PSTVd单感染高度相关，黄色模块与PVY单感染相关，蓝色模块则与共感染高度相关。蓝色模块的基因在PVY单感染中已有表达，在共感染中表达进一步增强，但PSTVd单感染无法诱导它们，这再次印证了PVY是驱动该共感染特异性转录程序的主要因素。
3.8. 共表达网络分析揭示响应PSTVd、PVY和共感染的独特且协同的调控模块
对上述特异性模块的深入分析发现，PSTVd相关模块富集于代谢和氧化还原过程；PVY相关模块富集于免疫信号和硫代谢；而共感染相关模块则显示出最复杂的富集模式，包括氨基酸生物合成、碳代谢和多激素信号通路，表明在共感染压力下宿主的转录活性和代谢协调性显著增强。
3.9. 代谢组分析揭示PSTVd、PVY及共感染植株的独特代谢特征
代谢组学分析检测了88种激素相关代谢物。PVY单感染引起了最广泛的激素代谢物变化（30个上调，9个下调），而共感染引起了最大数量的变化（35个上调，9个下调）。细胞分裂素（Cytokinin， CK）和生长素（Auxin）是受影响最显著的激素类别。特别值得注意的是，在共感染中，一些在PSTVd单感染中积累的CK代谢物（如2CltZ， tZOG）水平大幅下降，显示了PVY对PSTVd诱导的激素变化的拮抗作用。
3.10. 差异代谢物的KEGG富集分析
KEGG富集分析显示，玉米素（Zeatin）生物合成通路在PSTVd、PVY单感染及共感染中均被显著富集，成为共有的关键代谢扰动节点。PSTVd还特异性地富集了多种生物碱生物合成通路，而PVY则额外富集了二萜类生物合成通路。
3.11. 激素谱分析揭示单感染和共感染植株的不同调控模式
对具体激素的定量分析进一步证实了复杂的互作关系：PSTVd单感染特异性地积累了某些CK；而PVY和共感染则显著降低了这些CK的水平，但提升了赤霉素前体GA20和色氨酸（Tryptophan， TRP）的水平。茉莉酸途径的前体OPDA在单感染中升高，但在共感染中被抑制。这些结果表明，PVY和PSTVd在调控宿主激素网络上既存在拮抗，也存在协同效应。
3.12. 整合的基因-代谢物网络揭示对PSTVd、PVY及其共感染的独特分子响应
整合分析发现了基因表达与代谢物积累之间的显著相关性网络。在共感染中，一些关键基因（如WRKY26， PAP1）与激素代谢物（如JA-Ile， SA， ABA）呈现强相关，勾勒出一个整合了防御与生长信号的多层次调控网络。
3.13. 整合KEGG分析揭示PVY-PSTVd共感染下本氏烟中激素信号和基因调控的改变
深入的整合通路分析将焦点锁定在“植物激素信号转导”通路上。该通路在共感染中显著富集，涉及173个DEGs和3个差异积累代谢物（DAMs）：二氢玉米素核苷（Dihydrozeatin riboside， DHZR）、脱落酸（Abscisic Acid， ABA）和茉莉酸-异亮氨酸（Jasmonoyl-isoleucine， JA-Ile）。研究发现，共感染植株中DHZR和ABA水平下降，而JA-Ile水平上升。与此对应，细胞分裂素信号通路中的AHP基因、ABA受体PYL基因被下调，而茉莉酸受体COI1基因被上调。这些基因表达与代谢物变化的耦合，清晰地描绘了共感染如何重编程核心激素信号网络：抑制与生长维持相关的CK和ABA信号，同时增强与防御相关的JA信号。
研究结论与讨论部分强调，这项研究通过多组学整合，系统阐明了PVY与PSTVd共侵染的分子互作机制。核心结论是：PVY在共侵染中占据绝对主导地位，其强烈的转录和代谢重编程能力掩盖了PSTVd的影响，甚至可能通过抑制宿主防御为PSTVd的存续创造条件。具体表现为：1）PVY驱动了大规模的宿主基因表达重编程，广泛影响光合作用、初级和次级代谢；2）PSTVd则引发更特异、更温和的响应，主要靶向蛋白修饰和特定信号通路；3）共感染并非简单的叠加，而是产生了独特的、常表现为协同效应的转录和代谢特征，尤其是在激素信号网络（如JA、CK、ABA）的整合与重构上；4）玉米素生物合成是两者共同影响的关键代谢节点。这些发现不仅深化了对病毒-类病毒这种特殊互作模式的理解，揭示了植物整合多重生物胁迫信号的复杂策略，而且为未来通过干扰关键激素通路或利用分子标记辅助育种来管理马铃薯等作物的复合病害提供了重要的理论依据和潜在的靶点。研究结果凸显了在田间复杂环境下，针对优势病原体（如PVY）的防控对于遏制其协同病原体（如PSTVd）可能具有事半功倍的效果。
