馬蜂柑響應黃龍病菌不同侵染時期的生物學和轉錄組學分析

滕彩玲，鐘晰，吳昊娣，胡燕，周常勇，王雪峰

馬蜂柑響應黃龍病菌不同侵染時期的生物學和轉錄組學分析

滕彩玲1，鐘晰1，吳昊娣1，胡燕2，周常勇1，王雪峰1

（1中國農業科學院柑桔研究所國家柑桔工程技術研究中心，重慶 400712；2贛州市果業局，江西贛州 341000）

【】柑橘黃龍病（Huanglongbing，HLB）是制約柑橘產業發展的重大病害，我國發生的黃龍病由韌皮部桿菌亞洲種（‘Liberibacter asiaticus’，CLas）引起。本研究通過分析潛在耐黃龍病柑橘品種馬蜂柑（）在感染黃龍病后不同時期的生物學癥狀、顯微結構和轉錄組學差異，揭示馬蜂柑不同時期對黃龍病菌具有耐性的分子機理，為進一步深入挖掘抗性基因及開展抗病育種打下基礎。以嫁接在兩年生卡里佐枳橙砧木上的馬蜂柑作為供試材料，使用只感染CLas、其他常見柑橘病毒類病原均呈陰性的毒源為接毒材料對馬蜂柑進行接種，并在接種毒源后每隔15 d進行定期的熒光定量PCR檢測。將馬蜂柑接種毒源4個月（最早檢測到CLas陽性）作為感病前期處理，接種毒源14個月作為感病后期處理，以健康植株為對照，進行生物學癥狀和顯微結構觀察，分析其感病后不同時期的結構變化。結合比較轉錄組學分析與熒光定量PCR驗證，解析其耐病相關機制。癥狀觀察發現，馬蜂柑在感病前期和感病后期，植株葉片幾乎不顯癥；顯微結構觀察表明，馬蜂柑在感病前期中脈組織結構清晰，細胞形態正常，無淀粉粒積累的現象，僅在后期出現韌皮部極少數的篩管被堵塞；對轉錄組測序結果進行篩選鑒定，馬蜂柑感病前期共篩選鑒定到181個差異表達基因，感病后期共篩選鑒定到1 384個差異表達基因；比較轉錄組分析表明，馬蜂柑感病前期和后期的差異表達基因主要涉及細胞壁代謝、防御反應、淀粉與蔗糖代謝、胼胝質合成以及信號轉導等方面，其中在感病前期，淀粉與蔗糖代謝、細胞壁代謝相關的調控基因下調表達，在感病后期，水楊酸代謝及其信號轉導途徑、病程相關蛋白和谷胱甘肽轉移酶相關的調控基因上調表達。馬蜂柑響應CLas早期侵染主要表現為物理結構穩定、淀粉合成和光合作用等途徑不受干擾；而水楊酸介導的抗性信號轉導、效應蛋白觸發免疫反應（effector-triggered immunity，ETI）激活的病程相關蛋白和谷胱甘肽轉移酶參與的解毒作用是感病后期的主要耐病機制。

柑橘黃龍病；韌皮部桿菌亞洲種；馬蜂柑；耐病性；轉錄組

0 引言

【研究意義】柑橘黃龍病（Huanglongbing，HLB）是柑橘生產上的一種毀滅性病害，能侵染幾乎所有的柑橘栽培品種[1]，造成植株經濟壽命縮短，產量銳減，果品品質變劣，并可在短期內導致植株死亡，從而造成巨大的經濟損失[2]。自1919年Reinking首次在我國潮汕地區發現黃龍病以來，該病在中國的發生已有100多年，目前世界范圍內仍未研發出治療性的藥劑[3]。我國發生的柑橘黃龍病由韌皮部桿菌亞洲種（‘Liberibacter asiaticus’，CLas）引起，由于CLas純培養難題尚未攻克，病原致病機理尚不明確，極大程度地限制了病原生物學、病原-植物寄主互作及抗病基因工程育種的研究。以耐病品種馬蜂柑（）為研究材料，從生物學和轉錄水平上研究馬蜂柑響應不同時期CLas侵染的應答機制，可為抗/耐黃龍病常規和轉基因育種提供新思路。【前人研究進展】解剖學和生物學研究表明感染CLas引起的癥狀與淀粉累積和韌皮部組織壞死相關[4-7]。對不同柑橘品種及其近緣屬對黃龍病的敏感性評價發現澳指檬（）、來檬（）、檸檬（）、枳（）等品種具有一定的抗/耐病性[8-10]。近年來，利用轉錄組學、蛋白質組學和代謝組學方法鑒定柑橘中黃龍病應答基因已取得了較大進展，柑橘感染黃龍病后，不同部位中的基因表達譜均已被解析[11-16]。CLas侵染后柑橘嫩葉、成熟葉和果實中的糖代謝和淀粉代謝相關基因的表達均受到了影響，其中葉片中蔗糖代謝和淀粉合成相關基因均被誘導上調表達，顯癥果實中蔗糖和棉子糖代謝相關基因被誘導上調表達，而淀粉合成相關基因被抑制表達[17-18]。光合作用相關基因在帶病根、莖、葉中均下調表達[19-20]，而在帶病果實中卻上調表達[10,21]。運用轉錄組學或蛋白組學技術比較了不同敏感度柑橘品種的不同部位響應CLas侵染后的差異，發現病菌侵染可導致柑橘蔗糖代謝[22]、淀粉代謝[17-18]、光合作用[19-21]、激素代謝[23-24]的平衡失調。CLas對寄主植物體內糖類的偏好性分解利用（如果糖）會導致寄主糖類代謝平衡失調，糖類代謝的失衡可能導致染病植株呈現生理及病理方面的癥狀[18]。柑橘黃龍病帶病葉片的黃化癥狀與韌皮部堵塞和光合作用受抑制相關，而CLas侵染后寄主體內淀粉累積和胼胝質沉積是韌皮部堵塞的重要原因。淀粉累積的主要原因是淀粉分解過程受抑制[19]。感染CLas的柑橘中，茉莉酸（JA）和水楊酸（SA）的含量明顯高于未感染CLas的柑橘[25]。筆者實驗室前期對耐病品種馬蜂柑和感病品種甜橙（）早期感染黃龍病進行了比較轉錄組學研究，發現CLas侵染后甜橙的淀粉代謝和光合作用受到一定影響，但對馬蜂柑無顯著影響，馬蜂柑中細胞壁代謝的激活和部分抗病蛋白相關基因的高水平表達與其防御能力相關[26]。【本研究切入點】近年來，在不同類型柑橘響應黃龍病菌侵染的轉錄組和蛋白組方面開展了大量研究，但尚無針對同一柑橘品種在感染黃龍病不同時期的比較轉錄組學研究。因此，有必要將感病前期和后期的馬蜂柑進行生物學和轉錄組學分析，進一步了解感病不同時期馬蜂柑顯微結構及基因表達的變化。【擬解決的關鍵問題】以對柑橘黃龍病具有一定耐性的馬蜂柑為材料，分別設置感病前期和感病后期兩個時間節點。通過生物學和轉錄組學對兩個時間節點進行比較研究，觀察二者在相關代謝途徑和防御反應上的差異，揭示馬蜂柑在感病前期和后期的不同耐病調節機制。

1 材料與方法

試驗于2017—2019年在中國農業科學院柑桔研究所國家柑桔工程技術研究中心完成。

1.1 供試毒源及苗木

黃龍病毒源：筆者實驗室保存毒源采自江西贛州八鏡臺的琯溪蜜柚（cv. Guanximiyou）。該樣品中柑橘衰退病毒（，CTV）、柑橘碎葉病毒（Citrus tatter leaf virus，CTLV）、溫州蜜柑萎縮病毒（，SDV）、柑橘裂皮病類病毒（，CEVd）等常見柑橘病毒類病害檢測結果均呈陰性。試驗苗木：馬蜂柑接穗嫁接于兩年生卡里佐枳橙（×）砧木，苗木放置于國家柑橘苗木脫毒中心溫室，溫度為（25±3）℃。試驗處理：通過芽接的方式將毒源嫁接在試驗苗木上，兩個月后每隔15 d進行定期qPCR檢測。第一批次接種毒源傳毒并放置到14個月的馬蜂柑苗木作為感病后期處理；第二批次接種毒源并首次檢測到（約4個月）帶毒的馬蜂柑作為感病前期處理；同時期接種無毒琯溪蜜柚的植株作為健康對照。兩個時期分別選3株植株作為生物學重復，健康對照取2株作為生物學重復。

1.2 Spurr樹脂樣品制備及切片觀察

取感病前期和后期的馬蜂柑葉片及相應的健康葉片的中脈，切取2 mm×1 mm大小置于1.5 mL離心管中；加入固定液，固定過夜；用鉀鹽緩沖液漂洗樣品；用不同梯度的丙酮溶液浸漬樣品去除樣品水分；加入Spurr樹脂與純丙酮溶液（體積比1﹕3）滲透，每4 h追加一次的Spurr樹脂，3次后純樹脂過夜；在包埋模中放上滲透好的樣品，加入純樹脂后，放入70℃恒溫箱中聚合過夜；得到樹脂樣品，切片染色觀察。

1.3 RNA提取及轉錄組建庫、測序

取感病前期、后期和健康植株的完全展開嫩葉，使用Trizol法提取總RNA。用Agilent RNA 6000 Nano試劑盒處理樣品后，利用NanoDrop One檢測總RNA樣品的純度，用Agilent 2100檢測總RNA樣品的濃度、28S/18S 以及RIN值。質檢合格后的樣品用于文庫構建，再利用BGISEQ-500平臺進行雙末端測序，并進行3次技術重復（深圳華大基因股份有限公司）。

1.4 序列分析

測序的原始數據去除低質量、接頭污染和未知堿基N含量過高的reads后，得到clean reads，使用HISAT將clean reads比對到甜橙參考基因組序列（http://www. ncbi.nlm.nih.gov/genome/10702）。使用Bowtie 2將clean reads比對到參考序列，再使用RSEM計算基因和轉錄本的表達水平。使用PossionDis算法進行差異基因檢測，參數為fold change（FC）≥2.00，-value ≤0.005，FDR≤0.001。差異表達基因（differentially expressed gene，DEG）功能分析，將差異表達基因通過Gene Ontology（GO）數據庫（http://www.geneontology. org/）進行功能分類，同時基于NCBI mercator.results數據，利用MapMan軟件對差異基因進行代謝途徑分析。

1.5 差異表達基因的RT-qPCR驗證

選擇不同代謝途徑差異表達基因，同時選取在測序樣品中表達量一致的（甘油醛-3-磷酸脫氫酶）作為內參基因。以不同時期馬蜂柑葉片總RNA為模板，反轉錄合成cDNA第一鏈。利用GoTaq? qPCR Master Mix（Promega）熒光定量試劑盒進行實時定量PCR。使用Primer Premier 5.0 應用軟件設計RT-qPCR引物，由深圳華大基因股份有限公司合成（表1）。利用Bio-Rad iQ5進行2-ΔΔCT分析，并換算成log2FC，計算基因相對表達水平。

2 結果

2.1 馬蜂柑感染黃龍病前期、后期的生物學分析

觀察感染黃龍病前期（接種CLas 4個月，病原菌檢測陽性）和后期（接種CLas 14個月，病原菌檢測陽性）的馬蜂柑植株，發現馬蜂柑接種CLas后，未表現出明顯的黃龍病癥狀，生長也未受到顯著影響（圖1-A、1-B、1-C）。

表1 RT-qPCR基因及引物

通過光學顯微鏡觀察馬蜂柑健康和感病樣品的中脈組織發現，健康和感病樣品的大部分組織結構都比較清晰，皮層薄壁細胞排列整齊，木質部細胞大而疏松，韌皮部細胞排列緊密（圖1-D、1-E、1-F）。感病樣品的大多結構無明顯異常，進一步放大觀察可以發現感病前期樣品的韌皮部有極少數的篩管被堵塞，感病后期樣品篩管堵塞量比前期略多（圖1-H、1-I），未觀察到淀粉粒積累的現象。

2.2 差異表達基因及功能注釋

馬蜂柑感病前期，與健康對照相比，有149個基因上調表達，32個基因下調表達，共計出現181個差異表達基因；馬蜂柑感病后期，與健康對照相比，有386個基因上調表達，998個基因下調表達，共計出現1 384個差異表達基因。馬蜂柑感病前期與后期比較，有542個基因上調表達，1 115個基因下調表達，共計1 657個差異表達基因。

將差異表達基因進行GO功能分類，主要分為生物過程（biological process，BP）、分子功能（molecular function，MF）和細胞組分（cell component，CC）三大類。對三大功能類單獨進行進一步的分類和富集，得到不同時期、不同類別差異表達基因富集最多的前5個GO分類。馬蜂柑感病前期，生物過程共富集差異基因63個，富集最多的GO分類是光合作用和信號過程等；分子功能共富集差異基因68個，富集最多GO分類是轉移酶活性和電子載體活性等；細胞組分共富集差異基因37個，富集最多的是膜和膜部分相關等；馬蜂柑感病后期，生物過程共富集差異基因536個，富集最多GO分類是多糖代謝過程和外部封裝結構組織等；分子功能共富集差異基因623個，富集最多的是催化活性和轉移酶活性等；細胞組分共富集差異基因298個，富集最多是膜和膜部分等。

2.3 差異表達基因比較分析

利用MapMan軟件進行Pageman分析，將馬蜂柑感染CLas后兩個時間節點的差異表達基因代謝途徑富集可視化（圖2）。在馬蜂柑感染黃龍病前期，生物脅迫反應和核苷酸代謝下調表達，蛋白質降解（半胱氨酸蛋白酶）和LRR XII受體激酶上調表達；在馬蜂柑感染黃龍病后期，細胞壁相關、脂類代謝、過氧化物酶代謝、蛋白質降解（天冬氨酸蛋白酶）和氨基酸代謝下調表達，生物脅迫、激素代謝（水楊酸和茉莉酸）、細胞色素P450、RNA轉錄調控（MADS-box轉錄因子家族和WRKY轉錄因子家族）、蛋白質降解（泛素相關）和受體激酶等代謝途徑上調表達。

為進一步明確與耐病機制相關代謝途徑中差異基因表達的具體情況，利用Mapman軟件將馬蜂柑兩個時間節點感染黃龍病菌后的差異表達基因按功能進行分類，差異基因主要集中在糖代謝、光合作用、細胞壁代謝、激素代謝、PR-蛋白、氧化還原反應等代謝途徑。

2.3.1 淀粉和蔗糖代謝與光合作用相關差異基因的表達分析 CLas侵染寄主植物后，葉片出現黃化現象與CLas造成寄主植物葉片中淀粉過度累積進而導致葉綠體分解有關。在馬蜂柑感病前期，與蔗糖分解有關的-呋喃果糖苷酶轉化酶A基因（）和己糖激酶基因（）均下調表達；淀粉合成途徑中，在馬蜂柑感病前期并未出現相關基因的差異表達。馬蜂柑感病后期，與蔗糖分解有關的蔗糖合酶基因（）和己糖激酶基因（）呈下調表達；淀粉合成途徑中，葡萄糖-1-磷酸腺苷酰轉移酶基因（）下調表達（圖3-A）。

CLas的侵染在一定程度上抑制了馬蜂柑寄主的光合作用。馬蜂柑感病前期，與光系統Ⅱ相關的PsbD 和PsbC蛋白相關基因下調表達；感病后期，與光系統Ⅱ相關的PsbW和PsbP蛋白相關基因和多個捕光復合體II葉綠素a/b結合蛋白基因（）均下調表達（圖3-B）。

A：健康植株Mock plants；B：感病前期CLas-infection at the early stage；C：感病后期CLas-infection at the late stage；D：健康對照葉中脈leaf midrib of healthy control；E：感病前期葉中脈CLas-infected leaf midrib at the early stage；F：感病后期的葉中脈CLas-infected leaf midrib at the late stage；G：D局部放大照片magnification of D；H：E局部放大照magnification of E；I：F局部放大照片magnification of F；Co：皮層細胞cortex；Fi：纖維細胞fiber；Ph：韌皮部phloem；Pi：髓心細胞pith；X：木質部xylem；Sc：分泌腔secretory cavity；圖中箭頭所指的藍色斑點代表堵塞的篩管The blue spot indicated by arrow represents phloem plugging

CH-M-VS-CH-HLB-E：馬蜂柑染病前期與健康對照相比DEGs in CLas-infected C. hystrix at the early stagecompared withhealthy control； CH-M-VS-CH-HLB-L：馬蜂柑染病后期與健康對照相比DEGs in CLas-infected C. hystrix at the late stage compared withhealthy control。下同The same as below

彩色方塊代表上調或下調表達的差異基因Colored squares indicate up- or down-regulated genes。圖4同The same as Fig. 4

2.3.2 細胞壁代謝相關差異基因的表達分析 在馬蜂柑感病前期，編碼與纖維素合成相關的類纖維素合酶C基因家族基因（）、類纖維素合成酶D基因家族基因（）和果膠酯酶合成相關的果膠酯酶抑制因子基因（）均下調表達。在馬蜂柑感病后期，細胞壁代謝相關基因均呈下調表達，且下調表達倍數較大，包括纖維素合成相關的類纖維素合酶A基因家族基因（纖維素合酶A基因（）和COBRA-like蛋白基因（）；半纖維素合成相關的半乳糖醛酸轉移酶1樣蛋白基因（）、木聚糖-葡糖醛酸糖基轉移酶基因（）、-1,4-木糖基轉移酶基因（）；細胞壁蛋白相關的簇狀阿拉伯膠原蛋白基因（）；細胞壁降解相關的木葡聚糖內糖基酶基因（）、內切葡聚糖酶基因（）、鼠李糖半乳糖酸裂解酶B基因）；細胞壁修飾相關的擴張蛋白基因（、）、果膠酯酶基因（）和果膠酯酶抑制因子基因（）以及胼胝質合成相關的胞間連絲與胼胝質結合蛋白3基因（）和胼胝質合成酶3基因（）（圖4）。

2.3.3 激素相關差異基因的表達分析 在馬蜂柑感病前期，茉莉酸途徑中，多個具有ZIM結構域的轉錄抑制因子基因（）下調表達。在馬蜂柑感病后期，茉莉酸途徑中，與茉莉酸合成相關的脂氧合酶2基因（）上調表達、與信號轉導相關的轉錄因子基因（）下調表達；水楊酸途徑中，水楊酸合成相關的異分質酸合成酶基因（）、水楊酸甲酯化相關的苯甲酸羧甲基轉移酶基因（）均下調表達，水楊酸甲酯化相關的水楊酸羧甲基轉移酶基因（）與水楊酸糖基化相關的病原菌誘導的UDP葡萄糖基轉移酶基因（）以及水楊酸下游信號轉導相關的病程相關蛋白1基因（）均上調表達（圖4）。

2.3.4 免疫防御反應相關差異基因的表達分析 在馬蜂柑感病前期，病原相關分子模式觸發的免疫反應（PAMP-triggered immunity，PTI）系統中，與Ca2+-CaM信號系統相關的鈣依賴性蛋白激酶基因（）、鈣調蛋白基因（）、呼吸爆發氧化酶基因（）均表現為下調表達，與識別相關的多個受體蛋白基因（）上調表達但下游相關轉錄因子未出現差異表達；效應蛋白觸發的免疫反應（effector-triggered immunity，ETI）系統中，無相關基因的差異表達。而在馬蜂柑感病后期，PTI系統中，與Ca2+-CaM信號系統相關的環核苷酸門控通道蛋白基因（）的多個相關基因上調表達，鈣依賴性蛋白激酶基因（）、鈣調蛋白基因（）和呼吸爆發氧化酶基因（）輕微上調表達，與識別相關的多個受體蛋白基因（）上調表達且個數增多但下游相關轉錄因子依舊未出現差異表達；ETI系統中，相關的抗病蛋白基因（1.8倍）（1.7倍）（2.1倍）出現了上調表達（圖4）。

2.3.5 病程相關蛋白與抗病相關基因的表達分析 在馬蜂柑感病前期，多個抗病蛋白基因被CLas誘導下調表達，主要是TMV抗性蛋白。而在感病后期，前期下調的TMV抗性蛋白基因表現出回升趨勢，同時也新增多個抗病蛋白基因的差異表達，上調表達的基因包括抗病性蛋白（CC-NBS-LRR類）家族基因的等，含NB-ARC結構域抗病性蛋白基因的等；以及同時受SA和PTI調控的蛋白基因；屬于Kunitz家族胰蛋白酶和蛋白酶抑制劑蛋白的神秘果素樣蛋白基因（）和21 kD種子樣蛋白基因（）在馬蜂柑感病后期下調表達。在馬蜂柑感病后期，5個谷胱甘肽S轉移酶基因（）以及2個熱休克蛋白基因（）均上調表達（圖4）。

圖4 馬蜂柑感染黃龍病菌不同時期生物脅迫相關差異表達基因

2.4 RT-qPCR驗證

從不同的基因功能分類中挑選了12個基因進行RT-qPCR驗證，10個基因RT-qPCR結果與RNA-seq測序結果的相對表達量變化趨勢基本保持一致，證明RNA-seq測序結果具有較高的準確性。其中細胞壁代謝相關的在前期二者趨勢一致；在后期RNA-seq結果中下調1.6倍，而RT-qPCR結果則輕微上調，后期與前期相比上調幅度較大。胼胝質沉積相關的和抗病相關的在前期RT-qPCR結果與RNA-seq測序結果有細微偏差，但后期變化趨勢保持一致。值得關注的是：蔗糖和淀粉代謝相關的，激素代謝相關的以及抗病相關的的RT-qPCR結果均高于RNA-seq測序結果（圖5）。

3 討論

3.1 馬蜂柑響應CLas早期侵染的耐病機制

在感病柑橘的葉片中脈、莖和根中觀察到不同程度的淀粉積累和韌皮部堵塞現象[5,27]。Folimonova等[8]研究了30 種不同基因型的柑橘對黃龍病的響應，根據癥狀發展和持續生長的能力，將不同基因型分為敏感、中度耐受和耐受3類。對耐受型粗檸檬與敏感型甜橙的葉、莖和根組織進行解剖學比較分析，發現在粗檸檬中觀察到的破壞性解剖結構變化（韌皮部細胞塌陷、淀粉沉積、韌皮部纖維缺乏）比甜橙少。同樣的現象在Fan等[28]的研究中也得到了印證。Deng等[29]比較觀察敏感型‘Valencia’甜橙和耐病型‘LB8-9 Sugar Belle’雜柑、‘Bearss’檸檬的葉片和葉中脈顯微結構，發現與敏感型‘Valencia’甜橙的韌皮部相比，耐病型‘LB8-9 Sugar Belle’雜柑、‘Bearss’檸檬的韌皮部有大量的再生；韌皮部分裂程度越低，韌皮部再生能力越強，是影響不同柑橘品種HLB耐受性的兩個關鍵因素。馬蜂柑感病前期、后期均未明顯的黃龍病癥狀；比較觀察馬蜂柑健康和感病樣品的中脈組織切片，二者結構無明顯異常，且到感病后期所受影響依舊較小，表明馬蜂柑韌皮部結構穩定，使其對CLas具有一定的耐病性。

A: CH-M-VS-CH-HLB-E; B: CH-M-VS-CH-HLB-L; C: CH-HLB-E-VS-CH-HLB-L

CLas的侵染可導致韌皮部損傷和篩孔堵塞，這可能限制礦物質營養素和各種有機化合物從根際到植物或從老葉到新芽的吸收和運輸[13]。纖維素酶、擴張蛋白和果膠酯酶均與細胞壁的破壞有關，這些基因的高表達可能是導致敏感寄主表現癥狀的原因之一[30]。甜橙感染黃龍病后，多個編碼CESA、CSLA和COBL蛋白的基因顯著下調表達，使得寄主的適應能力減弱，從而有利于病原菌的侵染[26]。擴張蛋白、PP2和胼胝質與CLas感染引起的韌皮部損傷有關[31-32]。相關研究表明擴張蛋白可能導致CLas影響的葉片種類變硬，韌皮部細胞壁扭曲[32]。Wang等[30]研究發現編碼該蛋白的基因在感病品種‘Marsh’葡萄柚中表達水平高于耐病品種‘Jackson’葡萄柚，推測其可能與品種對黃龍病的敏感性相關。本研究中，相關基因的下調表達比例、數量隨病原菌侵染時間的推移而不斷增加，說明馬蜂柑細胞壁對CLas侵染的限制能力隨時間的推移在不斷減弱。而在馬蜂柑感病后期，擴張蛋白和木葡聚糖內糖基酶的表達下調，推測可能與抵抗病原菌侵染相關，但其具體調節模式有待進一步研究。

韌皮部壞死導致運輸受阻是蔗糖、淀粉等碳水化合物在葉片中大量積累的主要原因[33]。早期研究發現，甜橙感染黃龍病后，與淀粉合成有關的酶基因上調表達[27]。本研究中，蔗糖分解在感病前期受到抑制且隨病原菌侵染時間的推移而不斷加強；而淀粉合成途徑在感病前期不受影響、隨病原菌侵染時間的推移開始抑制；光合作用中光的吸收在馬蜂柑侵染過程中會受到持續的干擾，但差異表達相關基因較少。由此推測在蔗糖分解和淀粉合成通路上相關基因的下調表達和光合作用所受影響較小，是馬蜂柑對CLas侵染作出的耐病性調節反應。

3.2 馬蜂柑響應CLas后期侵染的主要耐病機制

植物通過激活病原物相關分子模式（pathogen- associated molecular pattern，PAMP）誘發免疫反應（PTI）或效應蛋白觸發免疫反應（ETI），保護自身免受病原體的攻擊[34-35]。植物PAMP的信號在細胞膜表面，而CLas存在于細胞內，由于這種物理隔離的存在，推測寄主對CLas進行防御反應時可能無PTI的參與[36-37]。Shi等[38]通過對耐病型‘Sun Chu Sha’雜柑和敏感型‘Duncan’葡萄柚的研究表明，兩種基因型間PTI水平的差異與HLB耐受性水平相關。高度或中度耐病的柑橘基因型中沒有發現PTI反應的跡象，由此推測其耐病性可能與抗性相關基因的表達、病程相關蛋白（PRs）的分泌和防御型植物激素（如水楊酸、茉莉酸）的激活相關[39]。本研究中，隨病原菌侵入時間的增加PTI系統中相關的上游調控基因種類和數目不斷增加，差異表達也由原來的下調表達轉變為上調表達模式，但未發現其下游調控基因的變化。因此，推測馬蜂柑在響應CLas的侵染過程中PTI免疫系統未能被激活。隨著侵染時間延長，在后期3種抗性蛋白的相關基因均出現了上調表達，且上調幅度較大，說明ETI系統參與了馬蜂柑后期防御CLas的過程。

植物激素在植物與病原菌互作的過程中扮演重要角色，參與植物抗病機制。以水楊酸、茉莉酸和乙烯為代表的免疫反應得到了廣泛的討論。植物受到病原菌侵染時，水楊酸介導產生的過敏性反應（HR）和系統獲得性抗性（SAR）使植物獲得一定的抗/耐病能力[40]。水楊酸羧甲基轉移酶（SAMT）是植物在遠端葉片累積水楊酸或水楊酸葡糖糖苷（SAG），轉化水楊酸的關鍵基因[41]。MeSA與柑橘黃龍病的耐受性有關，CsSAMT1和CsSABP2-1參與了MeSA與水楊酸的相互轉化，與耐病酸柚中MeSA的累積有關，耐病酸柚具有較強的SAR反應[42]。水楊酸依賴的通路，通常導致SAR標記致病相關（PR）基因（和編碼基因等）在局部和系統植物組織中的誘導表達[43]。根部中的上調表達使得施文格枳柚（×）與Cleopatra桔（）相比，其能更有效地抵抗CLas[44]。馬蜂柑感病后期，水楊酸羧甲基轉移酶基因（）的上調表達與甲基化水楊酸，遠距離傳輸防御信號，觸發相關防御反應；UDP葡萄糖基轉移酶基因（）的上調表達將水楊酸糖基化轉化為非生理活性的水楊酸酰葡萄糖酯（SGE）和SAG兩種水楊酸苷儲存在液泡中，為后續持續激活水楊酸觸發的反應做儲備；水楊酸途徑的標志基因的上調表達與抵抗CLas、提高馬蜂柑抗性相關。

不同類別的病程相關蛋白在CLas侵染后被誘導差異表達，包括受體蛋白激酶、含TIR-NBS-LRR、CC-NBS-LRR和NB-ARC等結構域的抗性蛋白、Kunitz家族胰蛋白酶和蛋白酶抑制劑蛋白等[16,18,45-46]。在本研究中，從前期抗病蛋白的下調表達，再到后期抗病蛋白表達的回升和新增抗病蛋白的上調，表明馬蜂柑后期防衛反應更強烈。研究表明，在CLas侵染檸檬后，某些防御蛋白反而下調表達，推測該類蛋白不是CLas誘導的特異蛋白，而是參與非寄主抗性的蛋白，這些蛋白的降低表達有助于CLas侵染下的寄主保留能量[19]。因此，在馬蜂柑中下調表達的病程相關蛋白，是否與馬蜂柑耐病性相關有待進一步研究。

Martinelli等[47]通過比較感病品種臍橙和耐病品種‘Volkameriana’檸檬感染黃龍病后的蛋白組學差異，發現與植物解毒途徑有關的4個谷胱甘肽-S-轉移酶（GSTs）僅在耐病品種‘Volkameriana’檸檬上調表達，推測解毒途徑可能在提高黃龍病抗性中發揮重要作用，并提出的高表達是HLB耐受基因型的潛在候選標記。本研究中，在馬蜂柑感病后期，5個均上調表達，可能與提高馬蜂柑對CLas的耐受性相關。

4 結論

通過比較耐病品種馬蜂柑在感染黃龍病不同時期的生物學和轉錄組學差異，發現馬蜂柑響應CLas早期侵染主要表現為韌皮部結構穩定，細胞壁代謝、淀粉合成和光合作用等途徑不受干擾；而水楊酸介導的抗性信號、效應蛋白觸發免疫反應激活病程相關蛋白、谷胱甘肽轉移酶參與的解毒作用是感病后期的主要耐病機制。

[1] BOVé J M. Huanglongbing: a destructive, newly-emerging, century-old disease of citrus., 2006, 88(1): 7-37.

[2] 胡燕,王雪峰,周常勇.柑橘黃龍病菌亞洲種、蟲媒及植物寄主互作研究進展.園藝學報,2016,43(9):1688-1698.

HU Y, WANG X F, ZHOU C Y. Recent advances in interactions among ‘Liberibacter asiaticus’, insect vector and plant host.,2016,43(9):1688-1698. (in Chinese)

[3] Reinking O A. Diseases of economic plants in Southern China.,1919, 8(4): 109-134.

[4] Schneider H. Anatomy of greening-diseased sweet orange shoots., 1968, 58(1): 1155-1160.

[5] Etxeberria E, Gonzalez P, Achor D, ALBRIGO G. Anatomical distribution of abnormally high levels of starch in HLB-affected Valencia orange trees., 2009, 74(1):76-83.

[6] Albrecht U, Bowman K D. Transcriptional response of susceptible and tolerant citrus to infection with‘Liberibacter asiaticus’., 2012, 185/186:118-130.

[7] Boava L P, Cristofani-yaly M,Machado M A. Physiologic, anatomic, and gene expression changes in,and their hybrids after ‘Liberibacter asiaticus’ infection., 2017, 107(5):590-599.

[8] Folimonova S Y, Robertson C J, Garnsey S M, GOWDA S, DAWSON W O. Examination of the responses of different genotypes of citrus to huanglongbing (citrus greening) under different conditions., 2009, 99(12):1346-1354.

[9] SHOKROLLAH H,ABDULLAH T L, SIJAM K, ABDULLAH S NA, ABDULLAH NAP. Differential reaction of citrus species in Malaysia to Huanglongbing (HLB) disease using grafting method., 2009, 4(1): 338.

[10] Martinelli F, URATSU S L, Albrecht U, REAGAN R L, PHU M L, Britton M, Buffalo V, Fass J, Leicht E, Zhao W, Lin D, D’Souza R, Davis C E, Bowman K D, DandekarA M. Transcriptome profiling of citrus fruit response to huanglongbing disease., 2012, 7(5): e38039.

[11] ALBRECHT U, BOWMAN K D. Gene expression in(L.) Osbeck following infection with the bacterial pathogen ‘Liberibacter asiaticus’ causing Huanglongbing in Florida., 2008, 175(3): 291-306.

[12] KIM J S, SAGARAM U S, BURNS J K, LI J L, WANG N. Response of sweet orange () to ‘Liberibacter asiaticus’ infection: microscopy and microarray analyses., 2009, 99(1): 50-57.

[13] FAN J, CHEN C, BRLANSKY R H, GMITTER F G, LI Z G. Changes in carbohydrate metabolism ininfected with ‘Liberibacter asiaticus’., 2010, 59(6): 1037-1043.

[14] Fan J, Chen C, Yu Q, BRLANSKY R H, LI Z G, GMITTER F G. Comparative iTRAQ proteome and transcriptome analyses of sweet orange infected by ‘Liberibacter asiaticus’.,2011, 143(3): 235-245.

[15] Albrecht U, Fiehn O, BOWMAN K D. Metabolic variations in different citrus root stock cultivars associated with different responses to Huanglongbing., 2016, 107: 33-44.

[16] ARITUA V, ACHOR D, Gmitter F G, ALBRIGO G, WANG N. Transcriptional and microscopic analyses of citrus stem and root responses to ‘Liberibacter asiaticus’ infection., 2013, 8(9):e73742.

[17] Martinelli F, Reagan R L, Uratsu S L, PHU M L, ALBRECHT U, ZHAO W, DAVIS C E, BOWMAN K D, DANDEKAR A M. Gene regulatory networks elucidating huanglongbing disease mechanisms., 2013, 8(9):e74256.

[18] Zhong Y, Cheng C Z, Jiang B, JIANG N H, ZHANG Y Y, HU M L, ZHONG G Y. Digital gene expression analysis of Ponkan mandarin (Blanco) in response to Asia citrus psyllid-vectored Huanglongbing infection., 2016, 17(7):1063.

[19] Fan J, Chen C, Yu Q, Khalaf A, Achor D S, Brlansky R H, Moore G A, Li Z G, Gmitter F G. Comparative transcriptional and anatomical analyses of tolerant rough lemon and susceptible sweet orange in response to ‘Liberibacter asiaticus’ infection., 2012, 25(11):1396-1407.

[20] Zhong Y, Cheng C Z, Jiang N H, JIANG B, ZHANG Y Y, WU B, HU M L, ZENG J W, YAN H X, YI G J, ZHONG G Y. Comparative transcriptome and iTRAQ proteome analyses of citrus root responses to ‘Liberibacter asiaticus’ infection., 2015,10(6):e0126973.

[21] Liao H L, Burns J K. Gene expression infruit tissues harvested from huanglongbing-infected trees: comparison with girdled fruit., 2012, 63(8):3307-3319.

[22] Duan Y P, Zhou L J, Hall D G, LI W B,Doddapaneni H, Lin H, Liu L, Vahling C M, Gabriel D W, Williams K P, Dickerman A, Sun Y J, COTTWALD T. Complete genome sequence of citrus huanglongbing bacterium, ‘Liberibacter asiaticus’ obtained through metagenomics., 2009, 22(8):1011-1020.

[23] Lu H, Zhang C, Albrecht U, Shimizu R, Wang G, BOWMAN K D. Overexpression of a citrus NDR1 ortholog increases disease resistance in., 2013, 4:157.

[24] Rosales R, Burns J K. Phytohormone changes and carbohydrate status in sweet orange fruit from Huanglongbing-infected trees., 2011, 30(3):312-321.

[25] Nehela Y, Hijaz F, Elzaawely A A, EL-ZAHABY H M, KILLINY N. Citrus phytohormonal response to ‘Liberibacter asiaticus’ and its vector., 2018, 102: 24-35.

[26] Hu Y, Zhong X, Liu X L, LOU B H, ZHOU C Y, WANG X F. Comparative transcriptome analysis unveils the tolerance mechanisms ofin response to ‘Liberibacter asiaticus’ infection., 2017, 12(12):e0189229.

[27] Koh E J, Zhou L, Williams D S, PARK J, DING N, DUAN Y P, KANG B H. Callose deposition in the phloem plasmodesmata and inhibition of phloem transport in citrus leaves infected with ‘Liberibacter asiaticus’., 2012, 249(3):687-697.

[28] Fan J, Chen C, Achor D S, BRLANSKY R H, LI Z G, GMITTER F G. Differential anatomical responses of tolerant and susceptible citrus species to the infection of ‘Liberibacter asiaticus’., 2013, 83:69-74.

[29] Deng H, Achor D S, Exteberria E, Yu Q, Du D, Stanton D, Liang G, Gmitter FG. Phloem regeneration is a mechanism for Huanglongbing-tolerance of ‘Bearss’ lemon and ‘LB8-9’ sugar Belle? mandarin., 2019, 10: 277.

[30] Wang Y S, Zhou L J, Yu X Y, STOVER E, LUO F, DUAN Y P. Transcriptome profiling of Huanglongbing (HLB) tolerant and susceptible citrus plants reveals the role of basal resistance in HLB tolerance., 2016, 7: 933.

[31] Lionetti V. PECTOPLATE: the simultaneous phenotyping of pectin methylesterases, pectinases, and oligogalacturonides in plants during biotic stresses.,2015, 6: 331.

[32] Owens R A, Blackburn M, Ding B. Possible involvement of the phloem lectin in long-distance viroid movement., 2001, 14(7):905-909.

[33] Folimonova S Y, Achor D S. Early events of citrus greening (Huanglongbing) disease development at the ultrastructural level., 2010, 100(9):949-958.

[34] Chisholm S T, Coaker G, Day B, STASKAWICA B J. Host-microbe interactions: shaping the evolution of the plant immune response., 2006, 124(4):803-814.

[35] Jones J D, Dangl J L. The plant immune system., 2006, 444(7117): 323-329.

[36] Zou H, GOWDA S, Zhou L, HAJERI S, CHEN G, DUAN Y. The destructive citrus pathogen, ‘Liberibacter asiaticus’ encodes a functional flagellin characteristic of a pathogen-associated molecular pattern., 2012, 7(9):e46447.

[37] Hao G X, Boyle M, Zhou L J, DUAN Y P. The intracellular citrus huanglongbing bacterium, ‘Liberibacter asiaticus’ encodes two novel autotransporters., 2013, 8(7):e68921.

[38] Shi Q C, Febres V J, Zhang S J, YU F H, MCCOLLUM G, HALL D G, MOORE G A, STOVER E. Identification of gene candidates associated with huanglongbing tolerance, using ‘Liberibacter asiaticus’ flagellin 22 as a proxy to challenge citrus., 2017, 31(2): 200-211.

[39] Da GRACA J V, Douhan G W, Halbert S E, KEREMANE M L, LEE R F, VIDALAKIS G, ZHAO H. Huanglongbing: An overview of a complex pathosystem ravaging the world’s citrus., 2016, 58(4):373-387.

[40] Zhang Y, Li X. Salicylic acid: Biosynthesis, perception, and contributions to plant immunity., 2019, 50: 29-36.

[41] ZHU F, Xi D H, YUAN S, XU F, ZHANG D W, LIN H H. Salicylic acid and jasmonic acid are essential for systemic resistance againstin., 2014, 27(6):567-577.

[42] Zou X, Bai X, Wen Q, XIE Z, WU L, PENG A, he y, xu l, chen s. Comparative analysis of tolerant and susceptible citrus reveals the role of methyl salicylate signaling in the response to Huanglongbing., 2019, 38:1516-1528.

[43] Van Loon L C, Rep M, Pieterse C M J. Significance of inducible defense-related proteins in infected plants., 2006, 44(1): 135-162.

[44] Wu J, Alfrez F M, JOHNSON E G, GRAHAM J H. Up-regulation of PR 1 and less disruption of hormone and sucrose metabolism in roots is associated with lower susceptibility to ‘Liberibacter asiaticus’., 2018, 67(6): 1426-1435.

[45] Fu S, Shao J, Zhou C, HARTUNG J S. Transcriptome analysis of sweet orange trees infected with ‘Liberibacter asiaticus’ and two strains of., 2016, 17: 349.

[46] Nwugo C C, Duan Y,Lin H. Study on citrus response to Huanglongbing highlights a down-regulation of defense-related proteins in lemon plants upon ‘. Liberibacter asiaticus’ infection., 2013, 8(6):e67442.

[47] Martinelli F, Reagan R L, Dolan D, FILECCIA V, DANDEKAR A M. Proteomic analysis highlights the role of detoxification pathways in increased tolerance to Huanglongbing disease., 2016, 16(1):167.

Biologic and Transcriptomic Analysis ofResponses to ‘Liberibacter asiaticus’ at Different Infection Stages

TENG CaiLing1, ZHONG Xi1, WU HaoDi1, HU Yan2, ZHOU ChangYong1, WANG XueFeng1

(1National Citrus Engineering Research Center, Citrus Research Institute, Chinese Academy of Agricultural Sciences, Chongqing 400712;2Ganzhou Bureau of Fruit Industry, Ganzhou 341000, Jiangxi)

【】Citrus Huanglongbing (HLB), associated with phloem-colonized ‘Liberibacter asiaticus’ (CLas), severely impedes worldwide citrus production. The objective of this study is to analyze the biological symptoms, microstructures and transcriptomes ofin response to CLas infection at different stages, to reveal the tolerance mechanism of, and to provide a basis for further screening of resistance genes and HLB-tolerant/resistant citrus breeding. 【】The budwoods ofgrafted on two-year-old Carrizo citrange rootstocks (×) used in this study were graft-inoculated with budwoods from a CLas (strain GZBJT)-infected Guanximiyou pummelo maintained in a greenhouse at Citrus Research Institute. The budwoods used as inoculum were tested CLas positive and free of other potential phloem-limited pathogens, such as(CTV) or Citrus tatter leaf virus (CTLV) by PCR before grafting. Inoculated plants were kept in greenhouse along with mock-inoculated healthy control plants. Real-time quantitative PCR (qPCR) was performed every 15 days after inoculation. Four months after inoculation (the earliest establishment of CLas by qPCR) and 14 months after inoculation were defined as the early infection stage and the late infection stage, respectively.The biological symptoms and microstructures were observed to analyze the structural changes of different infection stages. combined with comparative transcriptome and RT-qPCR validation, the response mechanism ofagainst HLB was explored.【】No typical symptom was observed inat the early and late stages of infection. Light microscopy observation from the midribs of HLB-affected and uninfectedrevealed that no significant structure change was found at the early infection stage and only a few sieves in the phloem were blocked at the late stage. By comparing the RNA-seq data, 181 and 1 384 genes were found to be differentially expressed at the early stage and late stage, respectively. Differentially expressed genes (DEGs) mainly involved incell wall metabolism, host defense response, starch and sucrose metabolism, callose synthesis and signal transduction. Comparative transcriptome analysis showed that the expression of related genes in starch and sucrose metabolism and cell wall metabolism was down-regulated at the early infection stage, and the expression of related genes in salicylic acid metabolism, salicylic acid signal transduction pathway, pathogenesis-related protein and glutathione-S-transferases was up-regulated at the late infection stage.【】The early response ofto CLas infection is mainly characterized by stable physical structure, undisturbed pathways such as starch synthesis and photosynthesisSalicylic acid-mediated resistance signals, effector-triggered immunity (ETI), and glutathione-S-transferases mediated detoxification contribute to the tolerance ofagainst CLas at the late infection stage.

citrus Huanglongbing; ‘Liberibacter asiaticus’ (CLas);; disease tolerance; transcriptome

10.3864/j.issn.0578-1752.2020.07.007

2019-09-01；

2019-10-15

國家重點研發計劃（2018YFD0201500）、國家自然科學基金（31871925，31671992）

滕彩玲，E-mail：864569130@qq.com。鐘晰，E-mail：767911863@qq.com。滕彩玲和鐘晰為同等貢獻作者。通信作者王雪峰，E-mail：wangxuefeng@cric.cn

（責任編輯 岳梅）