譚西北 李鵬 孫磊 樊秀彩 劉崇懷 姜建福 張穎

摘要：【目的】葡萄白腐病是危害最嚴重的病害之一，不同的葡萄種質對白腐病的抗性存在差異。在前期的研究中篩選到了1 份抗病種質刺葡萄0941，該種質為中國野生葡萄，具有耐濕熱、抗病性強的特點。旨在通過轉錄組測序和WGCNA篩選，挖掘刺葡萄抗白腐病關鍵基因。【方法】以抗病的刺葡萄0941（Vd）和感病的歐亞種美人指（Vv）為試材，通過室內接菌，設置0 h、24 h、48 h 三個時間點，研究病原菌誘導下感病和抗病種質的基因表達差異，篩選出抗病基因，并利用qRT-PCR進行驗證?！窘Y果】根據(jù)|Log2Fold Change|≥2 和p≤0.05 的篩選條件，在0 h、24 h、48 h 抗病品種刺葡萄和感病品種美人指之間篩選到差異表達基因分別有5266、4725、4653 個。其中上調基因數(shù)分別為1737、1788、1727 個，下調基因數(shù)為3529、2937、2926 個；進一步對這些差異基因進行GO和KEGG富集分析，發(fā)現(xiàn)了具有防御反應、跨膜運輸和一些植物激素信號轉導等功能的基因或者通路。為了明確不同品種葡萄果實抗病之間的調控網(wǎng)絡，進行了WGCNA分析，得到7 個基因模塊，篩選出2 個與抗病高度相關的關鍵的模塊并且通過hub 基因分析明確模塊內基因的互作網(wǎng)絡關系，鑒定出2 個類黃酮代謝相關基因、2 個植物激素信號轉導相關基因以及2 個hub 基因，并且通過熒光定量驗證其表達情況與轉錄組數(shù)據(jù)一致?！窘Y論】獲得了刺葡萄與美人指響應白腐病侵染的差異表達基因，并顯著富集在防御反應、跨膜運輸和一些植物激素信號轉導等通路，篩選出6 個抗病候選基因，包括F35H、EIN和MYB等。

關鍵詞：葡萄；抗病；轉錄組；基因挖掘；WGCNA 分析

中圖分類號：S663.1 文獻標志碼：A? 文章編號：1009-9980（2023）04-0653-16

葡萄在我國具有悠久的栽培歷史，因其具有豐富的營養(yǎng)物質和較高的經(jīng)濟價值，深受廣大民眾的喜愛。近年來，隨著葡萄栽培面積和產(chǎn)量的持續(xù)增長，病害造成的損失也逐年增加，其中葡萄白腐病是危害最嚴重的病害之一。該病是由白腐墊殼孢[Conielladiplodiella （Speg.） Petrk & Sydow.]引起的真菌性病害[1]，19 世紀在意大利首次報道，我國最早發(fā)現(xiàn)于1899 年，目前主要發(fā)生在我國東北、華北、西北和華東北部地區(qū)，植株感病后，葉片邊緣出現(xiàn)水漬狀病斑，慢慢侵染到葉片中部，呈現(xiàn)出不太明顯的同心輪紋狀的褐色病斑；葡萄果粒逐漸變褐變軟，發(fā)病后期整個果穗腐爛，造成大面積減產(chǎn)[2-3]。目前在農(nóng)業(yè)生產(chǎn)上主要是化學藥劑防治葡萄白腐病，化學防治法雖然效果好、見效快，但是長期用藥，不僅增加成本、降低果實品質，而且還會存在農(nóng)藥殘留，嚴重危害環(huán)境及食品安全。

不同葡萄種質對白腐病的抗性存在差異。歐亞種對白腐病的抗性較差，東亞種群的葡萄抗病性較強[4]。中國是葡萄屬植物重要的起源地之一，基本包括了所有的東亞種群，野生葡萄種質資源豐富[5]，利用野生種質資源導入異質抗性基因是解決栽培品種抗性基因貧乏的有效途徑。我國的刺葡萄對包括白腐病在內的多種真菌性病害都具有很強的抗性[6-7]，是葡萄白腐病抗性研究的良好材料。因此，解析葡萄抗病機制、挖掘抗性基因、培育具有較強抗病性的葡萄品種尤為重要。

轉錄組測序已被廣泛用于解析擬南芥、番茄、葡萄和其他植物的抗病機制；研究結果能夠提供植物在病原菌侵染下應答的信號通路和相關基因表達的信息。陳哲等[8]利用轉錄組測序分析，發(fā)現(xiàn)草莓受到炭疽病病原菌侵染后能夠影響類黃酮生物合成和植物激素信號傳導等通路中關鍵基因的表達水平，并鑒定出了19 個感病相關的差異表達基因。暹羅芽孢桿菌提高了貯藏杧果果實的抗病性，Jiang 等[9]利用暹羅芽孢桿菌的誘導機制，建立了杧果果實樣品在貯藏期間的比較轉錄組分析，部分基因（JAZ、BAK1 和PR1）被暹羅芽孢桿菌處理上調，引發(fā)應激反應，并誘導酚類抗性物質的合成，從而提高杧果果實的抗病性。其次，一些基因（WRKY22、HSP90、CNGCs、SOD、PAL、4CL、CHS 和HCT）如過氧化物酶體、苯丙素、黃酮類和姜酚生物合成，在杧果果實中上調，從而增強系統(tǒng)的抗病能力，刺激免疫反應。姜長岳[10]以毛葡萄丹鳳-2 葉片為材料接種白粉菌進行轉錄組測序，從MYB家族中篩選出抗病相關轉錄因子VqMYB154，發(fā)現(xiàn)轉錄因子VqMYB154 具有正向調控3 個芪合酶靶基因表達以增強植物抗病性的功能。Jiao 等[11]通過轉錄組分析中國野生葡萄對白粉病的抗病機制，發(fā)現(xiàn)感病品種湖南一號下調基因比例明顯高于山-24，而上調基因主要與防御反應和抗病代謝物的生物合成相關，同時發(fā)現(xiàn)12 號和13 號染色體上的NLR基因對白粉菌最敏感，其中許多基因在易感種質的所有感染階段始終下調。Liu 等[12]利用轉錄組測序分析，從河岸葡萄中篩選了與抗病性相關的轉錄因子ERF、MYB、WRKY和bHLH，確定了197 個候選基因，可能參與霜霉病菌引起的防御反應，其中有5 個基因富含亮氨酸重復序列（LRR）。沈才琦等[13]通過對感炭疽病歐亞種里扎馬特和抗炭疽病刺葡萄黑珍珠的雜交后代中感炭疽病7-2-6 株系和抗炭疽病7-1-8 株系進行轉錄組測序分析，發(fā)現(xiàn)葡萄響應膠孢炭疽菌侵染的差異表達基因，顯著富集在植物與病原體相互作用、類黃酮生物合成和MAPK信號途徑等3 個通路，并篩選出8 個抗病候選基因，包括MLO蛋白及WRKY、ERF轉錄因子等。轉錄組分析已廣泛用于發(fā)掘植物響應病病原菌脅迫的調控因子，這為研究刺葡萄抗病機制提供了理論依據(jù)。筆者在本研究中對白腐病病原菌侵染前后的感病歐亞種美人指和抗病中國野生刺葡萄的果實進行轉錄組測序分析，通過轉錄組分析已經(jīng)鑒定了幾個可能在抗葡萄白腐病中起關鍵作用的基因。

1 材料和方法

1.1 試驗材料

以抗白腐病刺葡萄0941 和感白腐病歐亞種美人指轉色期的果實為試驗材料，材料采自于國家葡萄種質資源圃（鄭州）。在葡萄果實的果柄處接種直徑1 cm 的白腐病原菌菌塊，每組接種18 粒葡萄，保濕并置于28 ℃恒溫箱中培養(yǎng)，同時在果實果柄處接種無菌水作為對照。分別在接菌0 h、24 h、48 h 收集去除感染部分后的果實，每組均設置3 個生物學重復。采用多糖多酚植物總RNA提取試劑盒（華越洋，貨號0416-50），參照說明書提取樣品總RNA，置于-80 ℃超低溫冰箱保存?zhèn)溆谩?/p>

1.2 文庫構建及測序

總RNA質量檢測合格后，通過深圳華大基因科技服務有限公司進行轉錄組測序。對得到的原始數(shù)據(jù)進行質控，去除低質量的reads 獲得的高質量的質控數(shù)據(jù)，并儲存于格式為fastq 的數(shù)據(jù)文件中。以葡萄基因組12 ×（https：//urgi.versailles.inra.fr/Species/Vitis）為參考基因組進行比對。利用HISAT2 軟件將得到的clean reads 比對參考基因組進行組裝。使用RSEM 軟件計算基因表達水平，然后轉換成FPKM（每千堿基轉錄序列的預期片段數(shù)/百萬堿基對測序）值。使用對比工具Bowtie2 將clean reads 與參考基因組進行比對以統(tǒng)計基因比對率，比對率都在90%以上。利用組裝的轉錄本與GO 數(shù)據(jù)庫、KEGG數(shù)據(jù)庫、NR數(shù)據(jù)庫、Swiss-Prot 數(shù)據(jù)庫、Pfam數(shù)據(jù)庫、COG數(shù)據(jù)庫以及KOG數(shù)據(jù)庫進行基因功能注釋。

1.3 差異表達基因篩選和富集分析

比對品種間基因的表達情況，篩選差異表達的部分進一步富集分析。使用DESeq2 軟件進行差異基因的篩選，將|log2FoldChange|≥2、p≤0.05 的基因作為顯著的差異表達基因，對篩選出的差異基因進行GO功能富集和KEGG通路分析，以獲得差異表達基因的主要生物學功能及主要生化代謝途徑。

1.4 WGCNA分析通過WGCNA分析進一步探究差異表達基因之間的調控網(wǎng)絡，旨在鑒定相關基因模塊及模塊內的hub 基因。對每個基因模塊進行評估，檢驗每個模塊中基因的相關性，利用Cytoscape 3.8.0 計算每個基因的相關性，鑒定每個模塊內的hub 基因。通過計算每個模塊的特征基因值，檢驗各個基因白腐病之間的相關性。此外，為了進一步探究模塊中基因的功能，利用R包（ClusterProfiler）進行GO和KEGG富集分析。

1.5 實時熒光定量PCR（qRT-PCR）分析

為驗證RNA-seq 測序檢測到DEGs，本研究基于分析挑取9 個差異表達基因進行qRT-PCR 分析。使用Primer3plus 在線網(wǎng)站（https：//www.primer3plus.com/index.html）設計引物，并由上海生工生物公司進行合成，引物具體信息見表1。然后用Fasting-King RT Kit（With gDNA）反轉錄試劑盒（天根）將RNA 反轉錄為cDNA。每個樣本3 個重復，使用FastStart Essential DNA Green Master 試劑盒按照其說明在ABI-7300 系統(tǒng)進行qRT-PCR 分析，反應程序：94 ℃ 30 s；94 ℃ 5 s，60 ℃ 30 s，40 個循環(huán)。以VvEF1-α 為內參基因[14]，用2- ΔΔCT 公式計算相對表達量。

2 結果與分析

2.1 轉錄組測序質量分析

轉錄組測序后對原始數(shù)據(jù)進行過濾，一共得到798.77 百萬條高質量的過濾序列，共產(chǎn)生119.82 Gb過濾堿基，每個樣本的Q30 均大于92%，清潔讀數(shù)率都在96%以上。轉錄組測序結果與葡萄參考基因組的比對率高，測序質量好，可用于后續(xù)分析（表2）。對各組重復之間進行相關性分析，每組樣品之間的相關系數(shù)均高于0.8（圖1），以確保試驗結果的可靠性。

2.2 差異表達基因篩選

篩選|log2FoldChange|≥2、p≤0.05 的基因作為顯著差異表達基因，在0 h、24 h、48 h 抗病品種刺葡萄（Vd）和感病品種美人指（Vv）之間的差異表達基因分別有5266、4725、4653 個（圖2-A）。其中上調基因數(shù)分別為1737、1788、1727 個，下調基因數(shù)分別為3529、2937、2926 個。在抗病品種和感病品種之間的差異表達基因均是下調基因數(shù)目遠大于上調基因數(shù)目，且24 h 和48 h 較0 h 下調基因數(shù)目大幅減少，3 個時期對比組共有差異表達基因有2684 個（圖2-B）。

2.3 GO、KEGG功能注釋及富集分析

將篩選出來的所有時期的差異表達基因進行GO富集分析，結果分為3 類，分別為生物學過程、細胞組分和分子功能。結果僅展示出富集程度最為顯著的前20 條的代謝途徑（圖3）。在膜、細胞、細胞部分等組分中具有較高程度的富集，在生物學過程中防御反應、跨膜運輸、代謝過程、細胞過程等具有較高程度的富集，在分子功能中轉錄調節(jié)活性和催化活性等具有較高程度的富集，結果表明這些功能在刺葡萄、美人指與白腐病的防御反應中發(fā)揮著重要作用。

對刺葡萄和美人指之間的差異表達基因進行KEGG代謝通路富集分析，取富集最為顯著的前20個通路進行分析（圖4）。通過KEGG分析，表明上調差異基因主要富集在類黃酮生物合成、黃酮和黃酮醇生物合成、苯丙烷生物合成、芪類和二芳基庚烷類及姜酚生物合成、α-亞麻酸代謝、亞油酸代謝等代謝途徑中，而下調差異基因主要富集在植物與病原菌的相互作用、類黃酮生物合成、苯丙烷生物合成、植物激素信號和絲裂原活化蛋白激酶（MAPK）信號轉導途徑、ABC運輸、苯丙氨酸代謝、氨基糖和核苷酸糖代謝、二萜生物合成、DNA復制等通路中。

2.4 刺葡萄抗白腐病的相關差異表達基因

刺葡萄和美人指接種白腐病病原菌0 h、24 h、48 h 后進行轉錄組分析，從KEGG富集的通路中選擇白腐病病原菌誘導后表達差異較大的基因，作為重要的差異表達基因。差異表達基因主要來自次生代謝和植物激素信號轉導兩個通路。

葡萄果實接種白腐病病原菌后，與苯丙烷、類黃酮、黃酮和黃酮醇生物合成等次生代謝途徑的相關基因在24 h 和48 h 均受到誘導而產(chǎn)生表達差異。圖5展示了0 h、24 h和48 h中苯丙烷、類黃酮、黃酮和黃酮醇生物合成途徑相關的差異表達基因（differentiallyexpressed genes，DEGs）及其轉錄水平。其中苯丙氨酸解氨酶PAL（Vitvi00g01431、Vitvi16g01503）、四香豆酸酯CoA 連接酶4CL（Vitvi16g00139）、查耳酮合成酶CHS（Vitvi16g00990、Vitvi16g01480）、莽草酸/奎寧酸羥基肉桂酰轉移酶HCT（Vitvi11g00735）、黃酮醇合酶FLS（Vitvi08g01575）這7個基因在受到白腐病原菌誘導48 h 后，在抗病品種刺葡萄中的表達量遠高于感病品種美人指。芪合酶STS（Vitvi16g01465、Vitvi16g01466）的表達量在受到病原菌誘導后，在刺葡萄和美人指中均呈現(xiàn)上升趨勢，且美人指中的表達量高于刺葡萄。

白腐病病原菌侵染葡萄果實時，植物激素信號轉導途徑中生長素、細胞分裂素、赤霉素、脫落酸、油菜素甾醇、茉莉酸和水楊酸等重要防御相關的信號通路均發(fā)生富集，82 個差異表達基因顯著富集于植物激素信號轉導途徑（圖6），39 個DEGs（TCH4，4個；AUXIAA，1 個；CH3，2 個；TF，1 個；JAZ，2 個；SAUR，16 個；PR1，2 個；MYC2，1 個；ETR；2 個；TIR，1 個）在刺葡萄中的表達水平明顯高于美人指。其中，水楊酸信號途徑下游基因PR1（Vitvi03g00860）和茉莉酸信號途徑下游基因MYC2（Vitvi02g00231）在刺葡萄接菌處理后表達上調，且PR1（Vitvi03g00860）在接菌48 h 表達水平很高，最高上調1000 倍（0.27＜FPKM＜1100），MYC2（Vitvi02g00231）在接菌0 h 和24 h 表達水平無明顯變化，但在48 h 上調10 倍（0＜FPKM＜13）；茉莉酸信號途徑下游基因JAZ（Vitvi01g02293）在刺葡萄接菌24 h 下調表達3 倍，而后上調100 倍；乙烯信號途徑上游基因ETR（Vitvi05g00684）在刺葡萄和美人指接菌后都上調。

2.5 接菌處理后轉錄因子的表達模式

DEGs 中轉錄因子有344 個，屬于39 個轉錄因子家族，主要包含40 個AP2/ERF、66 個MYB、22 個WRKY、35 個bHLH、24 個NAC、16 個LBD、16 個C2H2、11 個HD-ZIP 等（圖7）。在接菌處理后，其中10 個轉錄因子在美人指中的表達量明顯高于刺葡萄，例如AP2/ERF（Vitvi05g01073）在接菌24 h 美人指中的表達量是刺葡萄的60 倍，在接菌48 h 美人指的表達量是刺葡萄的35 倍。16 個轉錄因子在刺葡萄中的表達量顯著高于美人指，例如參與植物苯丙烷類次生代謝途徑的MYB（Vitvi02g01019）在接菌24 h 刺葡萄中的表達量是美人指的100 倍，在接菌48 h 刺葡萄的表達量是美人指的40 倍。8 個轉錄因子在受到白腐病原菌的誘導后的刺葡萄和美人指中表達模式相反，例如響應激素信號轉導的AP2/ERF（Vitvi16g01438）在接菌24 h 刺葡萄中表達上調12倍，而在美人指中下調表達10 倍；在接菌48 h 刺葡萄中表達下調表達10 倍，而美人指上調3 倍。

2.6 WGCNA分析鑒定刺葡萄抗病的關鍵基因

經(jīng)過濾FPKM＜1 的基因，4085 個差異基因用于構建加權基因共表達網(wǎng)絡。通過計算18 個樣本表達水平的相關系數(shù)聚類分析，樣本間聚類良好，未出現(xiàn)離群樣本。根據(jù)差異基因的FPKM 值進行相關度分析并聚類，相關度較高的基因被分配到同一個模塊中，圖中聚類樹的不同分支代表不同的基因共表達模塊，不同顏色代表不同的模塊，共劃分為7 個模塊，并繪制基因共表達網(wǎng)絡熱圖（圖8）。根據(jù)基因表達模塊的趨勢與抗病性狀的相關性分析，以皮爾遜相關系數(shù)（r＞0.5）和顯著性p 值的大?。?em>p＜0.05）為篩選條件MEyellow、MEblack、MEgreen模塊與抗病性狀呈顯著正相關，其中MEblack模塊在Vd0h、Vd24h、Vd48h 與抗病性狀呈顯著正相關，MEyellow 模塊在Vd48h 與抗病性狀呈顯著正相關。對這兩個抗病相關模塊中的差異表達基因進行KEGG富集分析，其富集的通路主要集中在類黃酮生物合成、MAPK信號途徑、苯丙烷生物合成、植物相互作用（圖9）。

選取Meyellow、MEblack 模塊內連通度最高的前150 個基因作為模塊的hub 基因，其中候選轉錄因子共有13 個，并通過WGCNA 預測與它們可能存在互作關系的DEGs，利用Cytoscape 軟件進行可視化分析（圖10）。在MEblcak 模塊中有4 個轉錄因子（MYB，2 個；EIN 和TCP，各1 個）與植物防御反應有關，其中與MYB4-Like（Vitvi17g00232）相關DEGs 的數(shù)量最多，共有104 個；這些基因主要與植物的代謝反應有關，如苯丙氨酸生物合成、色氨酸代謝、糖酵解、泛素介導的蛋白水解等。在MEyellow 模塊中有9 個轉錄因子（MYB，3 個；ERF，3 個；bHLH，2 個；bZIP，1 個）與植物防御反應有關，其中MYB1R1（Vitvi04g00076）有相關性的DEGs 數(shù)量為134 個，這些基因的功能注釋發(fā)現(xiàn)與多種植物激素有關，如生長素、乙烯、水楊酸和茉莉酸，還和氧化還原反應、類黃酮生物合成等有關。

通過KEGG富集分析、轉錄因子家族分析以及WGCNA中與抗病性相關聯(lián)的模塊，結合基因注釋，篩選出6 個與抗病相關的候選基因，分別是Vitvi01g02263（FAOMT）、Vitvi06g01183（F35H）、Vitvi06g01036（EIN3）、Vitvi04g00076（MYB1R1）、Vitvi17g00232（MYB4）和Vitvi07g02070（ERF098）。

2.7 DEGs的RT-qPCR

為了驗證RNA-Seq 數(shù)據(jù)的可靠性，從DEGs 中選擇9 個基因，包括2 個類黃酮代謝相關的DEGs、2個植物激素信號轉導相關的DEGs、2 個共表達模塊篩選到的關鍵基因和3 個隨機挑選的DEGs，通過RT-qPCR方法分析了這些基因在刺葡萄和美人指受白腐病病原菌侵染后的表達量變化，熒光定量顯示以上基因在接種病菌前后的相對表達量與轉錄組數(shù)據(jù)一致，這說明轉錄組數(shù)據(jù)是可靠的（圖11）。其中FAOMT、F35H、EIN3、MYB4 的表達量在所有時期中刺葡萄均明顯高于美人指；MYB1R1、ERF098 在接菌48 h 后，在刺葡萄中的表達量較其他時期和美人指中的均有顯著提升。

3 討論

本研究中，利用RNA-Seq 技術對接種白腐病原菌后不同時間點的抗病和感病葡萄葉片進行轉錄組測序，在0 h、24 h 和48 h 三個時期兩者之間分別篩得差異基因5266、4725和4653 個。當白腐病病原菌侵染葡萄果實后，多個代謝通路的基因表達發(fā)生改變，其中次生代謝物可能積極參與葡萄抗病反應。苯丙烷代謝物對植物的發(fā)育和生存至關重要，苯丙烷代謝物的生物合成和多樣性可通過加氧酶、氧化還原酶、連接酶和轉移酶等一系列酶來實現(xiàn)，這些酶通過?；?、甲基化、糖基化和羥基化對代謝物的基本骨架進行化學修飾[15]。苯丙氨酸由莽草酸途徑產(chǎn)生，苯丙類化合物是類黃酮、異黃酮和孜然芹素等多種酚類化合物的前體，研究發(fā)現(xiàn)苯丙氨酸在植物體內積累，能顯著降低植物對病原菌的敏感性[16]。苯丙氨酸解氨酶（PAL）是苯丙烷代謝途徑中的關鍵酶，能夠將苯丙氨酸轉為肉桂酸，是催化苯丙烷類次級代謝的第一步反應，起到樞紐的作用[17]。PAL參與花色苷、植物保護素、酚類代謝等物質的合成，該基因在大多數(shù)植物中是以多基因家族存在，PAL基因家族成員調控的代謝途徑也各有差異。有研究表明，植物中PAL活性降低時，抵抗外界干擾的能力下降[18]。肉桂酸4-羥基化酶（C4H）為通用的苯丙氨酸途徑限速酶，參與苯丙烷代謝途徑第二步催化反應，起到承上啟下的作用。該酶能提高自身的活性，來抵抗外界的環(huán)境脅迫。C4H將反式肉桂酸催化為β-香豆酸，通過4-香豆酰輔酶A連接酶轉化為香豆酞輔酶A，為下游類黃酮、木質素等次級代謝途徑提供前體物質[19]。PAL催化苯丙烷代謝合成途徑中的第一步反應產(chǎn)生的物質與C4H催化生成4-香豆酸。PAL 與C4H這兩個酶通常情況下均以協(xié)同的方式進行表達，來控制黃酮物質的合成，但通常情況下，C4H的活性比PAL的活性低很多[20]。在本試驗中PAL（Vitvi16g01503）、4CL（Vitvi16g00139）基因在受到病原菌誘導48 h 后，其在抗病品種刺葡萄中的表達量遠高于美人指。在植物體內PAL 作為苯丙烷類代謝定速酶，與植物的抗病性有重要關系。已報道在金針菇感染異型葡枝霉菌之后PAL活性明顯增強[21]，β-1，3-葡寡糖10 倍誘抗劑處理后的葡萄葉片苯丙氨酸解氨酶（PAL）活性增強，說明處理之后提高了體內防御酶和抗病相關蛋白活性，進而誘導植株產(chǎn)生系統(tǒng)抗病性[22]，以上結果表明PAL與植物的抗病性呈現(xiàn)正相關關系。

類黃酮代謝途徑是苯丙烷代謝途徑通過查爾酮合成酶（CHS）延伸的另一個重要分支，類黃酮代謝途徑廣泛存在于各種植物體內，類黃酮化合物是類黃酮代謝途徑中產(chǎn)生的大量次生代謝物質，大自然中類黃酮化合物一般分為六大類：黃酮、黃酮醇、黃烷酮、黃烷醇、花色素、異黃酮[23]。在植物中，葉片和花等器官的類黃酮化合物含量比較高，類黃酮化合物作為次生代謝物質在植物中可以作為植物和微生物之間的關鍵信號物質發(fā)揮植物防御素的功能，在植物中有很強的抗病抑菌作用。CHS 可以通過催化形成不同的類黃酮化合物來提高植物的抗病能力，也可通過改變植保素、黃酮素、抗氧化酶SOD和POD 等的積累以及防御酶的活性來提高防御能力[24-27]。研究表明，在番茄和蘋果中過表達CHS 能夠增強對病原菌等環(huán)境脅迫的抗性[28-29]?？箍菸〉拿藁ㄆ贩N中CHS表達量明顯上調后，類黃酮的含量提高從而增強其對枯萎病的抗性[30]。此外，郭澤西[31]已經(jīng)證明抗性品種中查爾酮合成酶基因CHSI的表達量比感病品種高，過表達CHSl能夠提高感病葡萄品種對灰霉病和霜霉病的抗性。在本研究中發(fā)現(xiàn)接種白腐病病原菌誘導48 h 后，在抗病品種刺葡萄中的表達量遠高于感病品種美人指，這與前人的研究一致。同時筆者通過加權共表達網(wǎng)絡分析發(fā)現(xiàn)的兩個候選基因：F35H和FAOMT，均是類黃酮代謝途徑中的關鍵基因，且FAOMT是F35H的下游基因[32]，二者在刺葡萄中的表達量均明顯高于美人指。類黃酮代謝通路由結構基因和調節(jié)基因共同控制，調節(jié)基因包括MYB、bHLH 等轉錄因子[33]，此前發(fā)現(xiàn)MYB4 能夠控制肉桂酸4-羥化酶（C4H）的負表達。筆者發(fā)現(xiàn)類黃酮代謝途徑中的關鍵基因和轉錄因子在抗病品種刺葡萄和感病品種美人指中表現(xiàn)出明顯的差異，因此，在刺葡萄響應白腐病病原菌侵染時，類黃酮物質可能起到重要作用。

此外，EIN3 和ERF098 是模塊中篩選出刺葡萄響應白腐病侵染的關鍵基因，轉錄因子EIN3 是乙烯信號轉導途徑中的調控因子，可激活乙烯反應因子（ERF）[34]。本研究中乙烯信號途徑下游基因EIN3的表達水平在處理不同時間點中刺葡萄均高于美人指，而EIN3 下游基因ERF098 在48 h 抗病品種刺葡萄中的表達量明顯高于其他時間和感病品種美人指，而在刺葡萄中48 h 的表達量明顯上升這可能與反應時間有關。由此可知，乙烯可能在刺葡萄抗白腐病反應的激素信號轉導中發(fā)揮了重要作用。

參考文獻References：

[1] AKKURT M，WELTER L，MAUL E，TOPFER R，ZYPRIAN E.Development of scar markers linked to powdery mildew （Uncinulanecator） resistance in grapevine （Vitis vinifera L.and Vitissp.）[J]. Molecular Breeding，2007，19（2）：103-111.

[2] ALI K，MALTESE F，ZYPRIAN E，REX M，CHOI Y H，VERPOORTER. NMR metabolic fingerprinting based identificationof grapevine metabolites associated with downy mildew resistance[J]. Journal of Agricultural and Food Chemistry，2009，57（20）：9599-9606.

[3] 王躍進，賀普超. 中國葡萄屬野生種抗白腐病機制研究[J]. 北方園藝，1988（1）：5-8.

WANG Yuejin，HE Puchao. Study of white rot disease resistantmechanism in ChineseVitis wild species[J]. Northern Horticulture，1988（1）：5-8.

[4] 劉崇懷. 中國葡萄屬（Vitis L.）植物分類與地理分布研究[D].鄭州：河南農(nóng)業(yè)大學，2012.

LIU Chonghuai. Studies on taxonomy and geographical distributionof Chinese wild grape species[D]. Zhengzhou：Henan AgriculturalUniversity，2012.

[5] 董陽輝，徐佩娟，王藝平，錢劍銳，何鐵海，朱朝磊，陳寶忠，王洪興，張建萍. 葡萄白腐病研究[J]. 江西農(nóng)業(yè)學報，2011，23（2）：107-110.

DONG Yanghui，XU Peijuan，WANG Yiping，QIAN Jianrui，HETiehai，ZHU Chaolei，CHEN Baozhong，WANG Hongxing，ZHANG Jianping. Study on white rot of grape[J]. Acta AgriculturaeJiangxi，2011，23（2）：107-110.

[6] 張穎，樊秀彩，孫海生，姜建福，劉崇懷. 葡萄種質資源對葉片白腐病的抗性鑒定及評價[J]. 果樹學報，2017，34（9）：1095-1105.

ZHANG Ying，F(xiàn)AN Xiucai，SUN Haisheng，JIANG Jianfu，LIUChonghuai. Identification and evaluation of resistance to whiterot in grape resources[J]. Journal of Fruit Science，2017，34（9）：1095-1105.

[7] 賈云飛. 刺葡萄0943VdWRKY70 轉錄因子抗病功能分析[D].洛陽：河南科技大學，2018.

JIA Yunfei. Function analysis ofVdWRKY70 transcription factordisease resistance fromVitis davidii 0943[D]. Luoyang：HenanUniversity of Science and Technology，2018.

[8] 陳哲，黃靜，趙佳，梁宏. 草莓應答炭疽菌侵染的轉錄組分析[J].植物保護，2020，46（3）：138-146.

CHEN Zhe，HUANG Jing，ZHAO Jia，LIANG Hong. Transcriptomicsanalysis of strawberry response toColletotrichum theobromicolainfection[J]. Plant Protection，2020，46（3），138-146.

[9] JIANG Z C，LI R，TANG Y，CHENG Z Y，QIAN M J，LI W，SHAO Y Z. Transcriptome analysis reveals the inducing effectofBacillus siamensis on disease resistance in postharvest mango[J]. Foods，2022，11（1）：107.

[10] 姜長岳. 中國野生毛葡萄VqMYB154 轉錄因子調控抗白粉病機理研究[D]. 楊凌：西北農(nóng)林科技大學，2021.

JIANG Changyue. Mechanism analysis of transcription factorVqMYB154 regulating the resistance to powdery mildew in ChinesewildVitis quinquangularis[D]. Yangling：Northwest A&FUniversity，2021.

[11] JIAO C，SUN X P，YAN X X，XU X Z，YAN Q，GAO M，F(xiàn)EI ZJ，WANG X P. Comparative transcriptome profiling of Chinesewild grapes provides insights into powdery mildew resistance[J].Phytopathology，2021，111（11）：2041-2051.

[12] LIU L，ZHANG B，WANG H，YU S Y，GUAN T S，HUANG YF，LIU C Y. Candidate resistance genes selection and transcriptomeanalysis for the early responses toPlasmopara viticola infectionin grape cultivars[J]. Journal of Plant Pathology，2020，102（3）：857-869.

[13] 沈才琦，孫磊，張穎，姜建福，劉崇懷，樊秀彩. 膠孢炭疽菌侵染不同抗性葡萄葉片的轉錄組學分析[J]. 果樹學報，2022，39（5）：730-742.

SHEN Caiqi，SUN Lei，ZHANGYing，JIANG Jianfu，LIU Chonghuai，F(xiàn)AN Xiucai. Transcriptomic analysis of different resistantleaves infected byColletotrichum gloeosporioides in grape[J].Journal of Fruit Science，2022，39（5）：730-742.

[14] 查倩，奚曉軍，蔣愛麗，田益華，王世平. 葡萄實時定量PCR 中穩(wěn)定內參基因的篩選[J]. 果樹學報，2016，33（3）：268-274.

ZHAQian，XI Xiaojun，JIANG Aili，TIANYihua，WANG Shiping.Identification of the appropriate reference gene through using areal- time quantitative PCR in grapes[J]. Journal of Fruit Science，2016，33（3）：268-274.

[15] WANG S C，ALSEEKH S，F(xiàn)ERNIE A R，LUO J. The structureand function of major plant metabolite modifications[J]. MolecularPlant，2019，12（7）：899-919.

[16] KUMAR V，HATAN E，BAR E，DAVIDOVICH-RIKANATI R，DORON- FAIGENBOIM A，SPITZER- RIMON B，ELAD Y，ALKAN N，LEWINSOHN E，OREN-SHAMIR M. Phenylalanineincreases chrysanthemum flower immunity againstBotrytiscinerea attack[J]. The Plant Journal，2020，104（1）：226-240.

[17] 梁冬. 大豆PAL2-3 基因的功能研究[D]. 哈爾濱：東北農(nóng)業(yè)大學，2017.

LIANG Dong. Functional study ofPAL2-3 gene in soybean[D].Harbin：Northeast Agricultural University，2017.

[18] 董艷珍. 植物苯丙氨酸解氨酶基因的研究進展[J]. 生物技術通報，2006（S1）：31-33.

DONG Yanzhen. Research advances on plant phenylalanine ammonia-lyase gene[J]. Biotechnology Bulletin，2006（S1）：31-33.

[19] 鄒廣平，李雅麗，馮夢薇，楊修，許智偉. 甘草懸浮細胞肉桂酸-4-羥基化酶（C4H）基因的克隆及表達分析[J]. 生物技術通報，2017，33（11）：96-100.

ZOU Guangping，LI Yali，F(xiàn)ENG Mengwei，YANG Xiu，XU Zhiwei.Cloning and expression analysis of cinnamate 4-hydroxylasegene fromGlycyrrhiza uralensis suspension cultured cells[J].Biotechnology Bulletin，2017，33（11）：96-100.

[20] 陳蒙. 山葡萄苯丙烷代謝相關基因的克隆表達與功能驗證[D].延吉：延邊大學，2019.

CHEN Meng. Cloning and expression and functional verificationof phenypropane pathway gene inVitis amurensis[D]. Yanji：Yanbian University，2019.

[21] 黃靜，李長田，賈傳文，王曉敏，呂瑞娜. 異型葡枝霉感染不同金針菇后防御酶系的活性變化[J]. 中國食用菌，2021，40（1），103-108.

HUANG Jing，LI Changtian，JIA Chuanwen，WANG Xiaomin，L? Ruina. Changes of defensive enzyme activities in differentFlammulina filiformis varieties after inoculation withCladobotryumvarium[J]. Edible Fungi of China，2021，40（1）：103-108.

[22] 王曉琳，吳琴燕，彭燕瓊，鄔劼，黃潔雪，毛妮妮，吉沐祥. 2 種生物誘抗劑對葡萄霜霉病的誘導抗病作用[J]. 中國農(nóng)學通報，2021，37（32）：127-131.

WANG Xiaolin，WU Qinyan，PENG Yanqiong，WU Jie，HUANG Jiexue，MAO Nini，JI Muxiang. Effects of two inducerson induced resistance to grapevine downy mildew[J]. ChineseAgricultural Science Bulletin，2021，37（32）：127-131.

[23] TAYLOR L P，GROTEWOLD E. Flavonoids as developmentalregulators[J]. Current Opinion in Plant Biology，2005，8（3）：317-323.

[24] QIAO Y，MENG X J，JIN X X，DING G H. Genes expression ofkey enzymes in phenylpropanes metabolism pathway in cucumberwith RT-PCR[J]. Advanced Materials Research，2013，746：53-57.[25] CHEN L J，GUO H M，LIN Y，CHENG H M. Chalcone synthaseEaCHS1 fromEupatorium adenophorum functions in saltstress tolerance in tobacco[J]. Plant Cell Reports，2015，34（5）：885-894.

[26] 陳帥. 煙草類黃酮代謝途徑中關鍵酶CHS 基因與R2R3MYB 類轉錄抑制因子功能研究[D]. 雅安：四川農(nóng)業(yè)大學，2017.

CHEN Shuai. Function analysis ofCHS genes and R2R3 MYBrepressors related to flavonoid biosynthesis pathway inNicotianatabacum[D]. Yaan：Sichuan Agricultural University，2017.

[27] RAHMAN F U，KHAN I A，ASLAM A，LIU R T，SUN L，WUY D，ASLAM M M，KHAN A U，LI P，JIANG J F，F(xiàn)AN X C，LIU C H，ZHANG Y. Transcriptome analysis reveals pathogenesis-related gene 1 pathway against salicylic acid treatment ingrapevine （Vitis vinifera L.）[J]. Front Genet，2022，13，1033288.

[28] KROL P，IGIELSKI R，POLLMANN S，KEPCZYNSKA E.Priming of seeds with methyl jasmonate induced resistance tohemi-biotrophFusarium oxysporum f. sp.lycopersici in tomatovia 12-oxo-phytodienoic acid，salicylic acid，and flavonol accumulation[J]. Journal of Plant Physiology，2015，179：122-132.

[29] JIAO C，SUN X P，YAN X X，XU X Z，YAN Q，GAO M，F(xiàn)EI ZJ，WANG X P. Grape transcriptome response to powdery mildewinfection：comparative transcriptome profiling of Chinese wildgrapes provides insights into powdery mildew resistance[J]. Phytopathology，2022，111（11）：2041-2051.

[30] 黃啟秀. 海島棉類黃酮代謝途徑中與枯萎病抗性相關基因的克隆及功能驗證[D]. 烏魯木齊：新疆農(nóng)業(yè)大學，2017.

HUANG Qixiu. Cloning and functional verification of genes relatedtoFusarium wilt resistance in flavonoid mtabolism pathwayof island cotton （Gossypium barbadense L.） [D]. Urumqi：Xinjiang Agricultural University，2017.

[31] 郭澤西，孫大運，曲俊杰，潘鳳英，劉露露，尹玲. 查爾酮合成酶基因在葡萄抗灰霉病和霜霉病中的作用[J]. 中國農(nóng)業(yè)科學，2022，55（6）：1139-1148.

GUO Zexi，SUN Dayun，QU Junjie，PAN Fengying，LIU Lulu，YIN Ling. The role of chalcone synthase gene in grape resistanceto gray mold and downy mildew[J]. Scientia AgriculturaSinica，2022，55（6）：1139-1148.

[32] LI P Q，RUAN Z，F(xiàn)EI Z X，YAN J J，TANG G H. Integratedtranscriptome and metabolome analysis revealed that flavonoidbiosynthesis may dominate the resistance ofZanthoxylum bungeanumagainst stem canker[J]. Journal of Agricultural andFood Chemistry，2021，69（22）：6360-6378.

[33] SHVACHKO N A，LOSKUTOV I G，SEMILET T V，POPOV VS，KOVALEVA O N，KONAREV A V. Bioactive components inoat and barley grain as a promising breeding trend for functionalfood production[J]. Molecules，2021，26（8）：2260.

[34] GU S Y，WANG L C，CHEUH C M，LO W S.CHITINASELIKE1 regulates root development of dark- grown seedlings bymodulating ethylene biosynthesis inArabidopsis thaliana[J].Frontiers in Plant Science，2019，10：600.