在農業生產中，作物常常面臨多種病原體的“組團”攻擊，這使得病害防控變得異常復雜。想象一下，一株馬鈴薯可能同時遭受病毒和類病毒的侵擾，它們如何“狼狽為奸”或“相互競爭”，最終導致作物減產甚至絕收？這背后的分子“暗戰”一直是科學家們想要破解的謎題。特別是馬鈴薯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）對關鍵差異表達基因進行驗證。

研究結果層層遞進，揭示了從田間現象到分子機制的完整故事：

對重慶地區256份馬鈴薯樣本的檢測發現，PVY是田間最流行的病原體，而PSTVd極少單獨發生，卻頻繁與PVY共同出現，尤其是在多重感染組合中。這提示在自然條件下，PSTVd可能依賴于與PVY的共感染來維持其存在或傳播，兩者之間存在潛在的生物學關聯。

在本氏煙和馬鈴薯中的實驗表明，PVY單感染引起嚴重的植株矮化、葉片黃化和卷曲，而PSTVd單感染癥狀輕微。共感染植株的表型與PVY單感染幾乎無法區分，PSTVd的癥狀被完全掩蓋。分子水平上，共感染時PSTVd的積累量反而比單感染時略有下降，而PVY的積累則維持在高水平。這證實了PVY在共感染中無論是表型還是分子水平上都占據主導地位，甚至可能對PSTVd的復制產生一定的競爭或抑制。

RNA-seq分析顯示，PVY感染（無論單感染還是共感染）引發了宿主基因組的劇烈重編程，差異表達基因（Differentially Expressed Genes， DEGs）數量遠多于PSTVd單感染。主成分分析（Principal Component Analysis， PCA）清晰地將不同處理組分開，共感染的轉錄譜更接近于PVY單感染，但仍具有獨特性。

基因本體（Gene Ontology， GO）和京都基因與基因組百科全書（Kyoto Encyclopedia of Genes and Genomes， KEGG）通路富集分析顯示，PSTVd單感染主要富集于蛋白磷酸化、激酶活性以及與植物-病原互作相關的通路。而PVY單感染和共感染則強烈富集于光合作用、碳代謝、次級代謝物生物合成等全局性代謝通路，以及植物激素信號轉導通路。這表明PSTVd的影響更側重于特定的信號調控，而PVY則對宿主的整體代謝和能量平衡造成了廣泛沖擊。

具體到功能基因，PVY強烈上調了RNA沉默通路的核心組件（如AGO1， AGO2， RDR1， RDR6）和免疫相關基因（如PR1， EDS1， PAD4）。PSTVd對這些通路的激活作用較弱。兩者都抑制了光合作用相關基因，但PVY的抑制作用更強。在細胞壁重塑和次級代謝方面，PVY和共感染引起了更廣泛的基因上調。

轉錄因子（Transcription Factor， TF）分析發現，PVY激活了包括WRKY、NAC、MYB、bHLH在內的多個TF家族的大量成員，而PSTVd激活的TF數量較少且更特異。共感染條件下，PVY的主導作用壓制了部分PSTVd誘導的TF，但同時協同激活了與茉莉酸（Jasmonic Acid， JA）和赤霉素（Gibberellin， GA）信號相關的TF（如TIFY10A， SCL3），展現了復雜的轉錄調控交叉對話。

WGCNA分析鑒定出與不同感染狀態特異性相關的基因共表達模塊。例如，棕色模塊與PSTVd單感染高度相關，黃色模塊與PVY單感染相關，藍色模塊則與共感染高度相關。藍色模塊的基因在PVY單感染中已有表達，在共感染中表達進一步增強，但PSTVd單感染無法誘導它們，這再次印證了PVY是驅動該共感染特異性轉錄程序的主要因素。

對上述特異性模塊的深入分析發現，PSTVd相關模塊富集于代謝和氧化還原過程；PVY相關模塊富集于免疫信號和硫代謝；而共感染相關模塊則顯示出最復雜的富集模式，包括氨基酸生物合成、碳代謝和多激素信號通路，表明在共感染壓力下宿主的轉錄活性和代謝協調性顯著增強。

代謝組學分析檢測了88種激素相關代謝物。PVY單感染引起了最廣泛的激素代謝物變化（30個上調，9個下調），而共感染引起了最大數量的變化（35個上調，9個下調）。細胞分裂素（Cytokinin， CK）和生長素（Auxin）是受影響最顯著的激素類別。特別值得注意的是，在共感染中，一些在PSTVd單感染中積累的CK代謝物（如2CltZ， tZOG）水平大幅下降，顯示了PVY對PSTVd誘導的激素變化的拮抗作用。

KEGG富集分析顯示，玉米素（Zeatin）生物合成通路在PSTVd、PVY單感染及共感染中均被顯著富集，成為共有的關鍵代謝擾動節點。PSTVd還特異性地富集了多種生物堿生物合成通路，而PVY則額外富集了二萜類生物合成通路。

對具體激素的定量分析進一步證實了復雜的互作關系：PSTVd單感染特異性地積累了某些CK；而PVY和共感染則顯著降低了這些CK的水平，但提升了赤霉素前體GA20和色氨酸（Tryptophan， TRP）的水平。茉莉酸途徑的前體OPDA在單感染中升高，但在共感染中被抑制。這些結果表明，PVY和PSTVd在調控宿主激素網絡上既存在拮抗，也存在協同效應。

整合分析發現了基因表達與代謝物積累之間的顯著相關性網絡。在共感染中，一些關鍵基因（如WRKY26， PAP1）與激素代謝物（如JA-Ile， SA， ABA）呈現強相關，勾勒出一個整合了防御與生長信號的多層次調控網絡。

深入的整合通路分析將焦點鎖定在“植物激素信號轉導”通路上。該通路在共感染中顯著富集，涉及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）可能具有事半功倍的效果。