《Frontiers in Plant Science》：Integrated transcriptome and metabolome analysis reveals the molecular mechanism underlying differences in Psa resistance between Actinidia valvata and Actinidia chinensis
編輯推薦：
本研究通過整合轉錄組與代謝組學分析，揭示了抗病物種綿毛獼猴桃(Actinidia valvata)與感病品種中華獼猴桃‘紅陽’(A. chinensis ‘Hongyang’)在感染丁香假單胞菌獼猴桃致病變種(Psa)后響應的分子機制。研究發現，A. valvata通過強烈激活植物-病原互作、MAPK和植物激素信號轉導等防御通路，并特異性上調木質素合成關鍵基因(CCR、CAD、POD)，促進成熟木質素聚合物積累以加固細胞壁，從而展現出卓越的Psa抗性。而感病品種‘紅陽’則表現出微弱的轉錄響應和代謝失調，其木質素合成途徑受阻于中間產物階段。該研究為解析獼猴桃潰瘍病抗性機制提供了多組學層面的見解，并為抗病育種提供了關鍵候選基因和代謝靶點。
獼猴桃是一種具有高經濟價值的水果作物，但極易受到由丁香假單胞菌獼猴桃致病變種(Pseudomonas syringae pv. actinidiae, Psa)引起的潰瘍病的侵害。為剖析其抗性機制，本研究對抗病物種綿毛獼猴桃(Actinidia valvata)和感病品種中華獼猴桃‘紅陽’(A. chinensis ‘Hongyang’， HY)在Psa侵染后進行了整合的轉錄組和代謝組分析。
表型觀察揭示抗感性差異
研究首先觀察了兩種材料接種Psa后的表型變化。高抗的A. valvata枝條水漬狀病斑平均長度為12.89毫米，而高感的HY枝條病斑平均長度達45.31毫米，兩者存在顯著差異，明確了A. valvata的高抗性和HY的高感性表型。
轉錄組分析揭示分子層面防御差異
通過對枝條進行高通量RNA-Seq分析，發現A. valvata在接種Psa后展現出強烈的防御相關轉錄響應，其差異表達基因(DEGs)數量約為HY的5倍。具體而言，在A. valvata接種與對照的比較中，鑒定出7931個DEGs；而在HY的比較中僅鑒定出1622個DEGs。對A. valvata特有的765個DEGs以及A. valvata與HY共有的1781個DEGs進行KEGG富集分析，發現它們顯著富集在植物-病原互作、MAPK信號通路和植物激素信號轉導等核心防御通路上。相比之下，HY特有的106個DEGs僅富集于淀粉和蔗糖代謝、苯丙烷生物合成等基礎代謝途徑，未顯示出與防御信號相關的顯著富集。此外，在A. valvata中上調的DEGs中包含3個病程相關(PR)蛋白基因和29個抗性受體激酶(RK)基因，這些都是已知的植物免疫關鍵組分。
代謝組分析揭示代謝物積累模式
利用UPLC-MS/MS技術分析枝條韌皮部代謝物，共鑒定出1229種代謝物，其中黃酮類、酚酸類占主導地位。主成分分析(PCA)顯示，Psa感染誘導了顯著的代謝變化，尤其是在A. valvata中。K-均值聚類分析將差異累積代謝物(DAMs)分為9個簇，其中簇1和簇6的代謝物（主要包括黃酮類、萜類、生物堿、酚酸等）在A. valvata接種Psa后特異性誘導積累，而在HY中無明顯變化。這些物質在植物生物脅迫響應中扮演重要角色。值得注意的是，在A. valvata接種后，代謝物總體呈現凈上調趨勢，而HY的代謝物則呈現凈下調，反映出兩者在應對病原侵染時截然不同的代謝重編程策略。
整合分析構建系統性調控網絡
通過聯合分析轉錄組和代謝組數據，發現A. valvata在感染后顯著富集的通路包括淀粉和蔗糖代謝、碳代謝、氨基糖和核苷酸糖代謝以及苯丙烷生物合成等。這些通路的激活為防御反應提供了能量和碳骨架前體。相比之下，HY僅微弱富集了苯丙烷生物合成和淀粉/蔗糖代謝兩條通路。加權基因共表達網絡分析(WGCNA)進一步將基因聚類為17個模塊，其中MEbrown模塊與苯丙烷生物合成代謝物L-苯丙氨酸(pme0021)的積累顯著相關。該模塊包含3281個基因，其中2763個被定義為核心基因，與1781個共有DEGs有766個重疊。
木質素合成通路是抗性的關鍵
深入分析發現，苯丙烷生物合成通路，特別是木質素合成分支，是A. valvata抗Psa的關鍵。與HY相比，A. valvata在接種后特異性積累了關鍵的木質素相關前體代謝物L-苯丙氨酸。同時，木質素合成途徑中的關鍵基因，包括肉桂酰-CoA還原酶(CCR)基因(AVb02g00135)、肉桂醇脫氫酶(CAD)基因(AVa06g01023, AVb06g00946)和過氧化物酶(POD)基因(AVa01g00801, AVb01g00780)在A. valvata中顯著上調。其中，兩個POD基因的表達量在接種5天后相較于HY分別上調了15倍和30倍。POD在木質素單體的氧化聚合形成木質素大分子的過程中起著至關重要的作用。對關鍵基因啟動子的順式作用元件預測發現，它們含有WRKY和MYB轉錄因子特異性結合的W-box和MYB結合位點，暗示這些基因可能受WRKY/MYB轉錄因子的直接調控，從而協同驅動木質素的高效合成。
表型與分子證據的相互印證
木質素含量的測定進一步證實了上述分子發現。抗病品種A. valvata在接種Psa后，其枝條中的木質素含量顯著增加。綜合來看，A. valvata通過激活以WRKY/MYB等轉錄因子為核心的調控網絡，上調CCR、CAD、POD等關鍵酶基因的表達，從而高效地將L-苯丙氨酸等前體轉化為成熟的木質素聚合物，加固細胞壁，形成抵御Psa入侵的物理屏障。相反，感病品種HY的木質素生物合成途徑受阻于中間產物階段，未能有效合成成熟的木質素，導致細胞壁脆弱，病原菌得以在韌皮部定殖和擴展。
討論與展望
研究揭示了A. valvata與HY在應對Psa侵染時存在根本性的分子與代謝策略差異。A. valvata具備一個預先配置或可快速誘導的免疫調控網絡，能夠協調激活廣泛的防御通路和代謝重編程，特別是通過強化木質素合成來加固細胞壁。而HY則缺乏有效的防御轉錄響應，其代謝狀態受到抑制，木質素合成通路不完整。本研究鑒定出的核心抗性相關基因（如CCR、CAD、POD）及共表達模塊，為獼猴桃抗潰瘍病分子育種提供了寶貴的候選靶點。未來的研究應通過CRISPR/Cas9基因編輯或過表達等遺傳手段對這些核心基因進行功能驗證，并通過田間試驗評估其抗性的穩定性，同時深入解析WRKY/MYB等轉錄因子在調控該網絡中的具體作用，為培育抗Psa的商業化獼猴桃品種奠定理論基礎。
