黃啟秀 曲延英 姚正培 李夢雨 陳全家

?

海島棉枯萎病抗性與類黃酮代謝途徑基因表達量的相關性

黃啟秀 曲延英 姚正培 李夢雨 陳全家*

新疆農業(yè)大學農學院 / 新疆農業(yè)大學農業(yè)技術重點實驗室, 新疆烏魯木齊 830052

枯萎病是危害海島棉生產的重要因素之一, 研究枯萎病抗性分子機制將為培育抗病海島棉品種、解決枯萎病對海島棉的危害問題提供堅實的基礎。本研究在前期轉錄組測序的基礎上, 對海島棉枯萎病抗性差異表達基因進行分析(Differentially Expressed Gene, DEG); 以7個抗病性表現(xiàn)不同的海島棉品種為材料, 利用qRT-PCR方法研究抗病差異表達基因在接種0~40 h的表達量差異, 分析基因表達量與病情指數的相關性。結果表明, DEG分析得出類黃酮代謝通路相關基因與海島棉枯萎病抗性有關。qRT-PCR分析顯示抗病材料中類黃酮代謝通路關鍵基因的表達量顯著高于感病材料。在接菌后多個時間點, 類黃酮代謝通路中的關鍵基因、和在抗病材料中的表達量顯著或極顯著高于感病材料, 其中和基因的表達量與病情指數呈顯著負相關。綜上所述, 類黃酮代謝通路相關基因對海島棉枯萎病抗性均有影響, 且、和基因是關鍵基因。

海島棉, 枯萎病, 類黃酮, 表達分析, 轉錄組測序

海島棉(L)在棉花栽培種中纖維品質最優(yōu)[1], 枯萎病的浸染導致其產量和纖維品質降低[2]。研究表明在相同條件下對陸地棉(L.)和海島棉接種枯萎病菌, 后者的感病程度較重[3]。陸地棉枯萎病抗性遺傳與海島棉不同, 校百才等[4]研究發(fā)現(xiàn)陸地棉的枯萎病抗性以加性效應為主, 而Smith等[5]認為海島棉的抗性是由2個具有加性效應的顯性基因共同決定。隨著分子生物學技術的發(fā)展, 從分子水平上研究棉花的抗病機制已有諸多報道, Yang等[6]利用基因沉默技術證明在海島棉抗黃萎病過程中扮演重要角色。Gao等[7]研究發(fā)現(xiàn)和是影響棉花抗黃萎病性的關鍵基因。但目前關于海島棉抗枯萎病性基因研究的報道仍然較少, 挖掘海島棉抗病基因, 為海島棉抗病育種提供依據十分關鍵。

近年來, 越來越多的研究者采用轉錄組測序(RNA-seq)挖掘功能基因。Xu等[9]以海島棉抗病品種海7124為材料, 利用RNA-seq技術分析其接種黃萎病菌V991菌株后的轉錄組變化, 檢測到3442個差異表達基因, 并認為木質素代謝通路相關基因在棉花抗病防御反應中扮演重要角色。王春暉[10]對陸地棉抗黃萎病自交系1017012及其接菌后6個時間段的RNA混樣進行轉錄組測序, 基因差異表達分析發(fā)現(xiàn)接菌后上調表達基因59 856個, 下調表達基因51 976個, 128條代謝途徑相關基因富集表達。

類黃酮化合物作為重要的次生代謝物, 廣泛分布于維管植物[11]。研究表明其與植物抗病性密切相關。Malhotra等[12]和賈振華[13]發(fā)現(xiàn)類黃酮化合物槲皮素在番茄和擬南芥的抗病方面發(fā)揮積極作用。Anna等[14]研究證實, 25種類黃酮物質可抑制黃萎病菌絲的生長。關于類黃酮代謝途徑結構基因與植物抗病性的關系已有研究: 轉反義基因楊樹中抗病物質兒茶素(Catechin)含量明顯降低[15], GbCHI蛋白能顯著抑制大麗輪枝菌()孢子的萌發(fā)和菌絲的生長[16]。植物在遭受病原微生物浸染后活性會增強, 說明由它調控的類黃酮代謝途徑是植物抵抗病原菌的防衛(wèi)反應之一[17]。雖然類黃酮代謝途徑中個別基因對植物抗病的影響已有諸多研究, 但系統(tǒng)性總結并揭示海島棉類黃酮代謝途徑中各基因與海島棉枯萎病抗性的關系還未見報道。本研究以7個抗病性表現(xiàn)不同的海島棉品種為材料, 研究RNA-Seq富集到的6個類黃酮代謝途徑相關基因及另外2個關鍵基因和在0 h和其他5個接菌時間點的表達量差異, 并對基因表達量與病情指數進行相關性分析, 以期探究類黃酮代謝途徑相關基因影響海島棉枯萎病抗性的關鍵時期和關鍵基因, 為培育海島棉抗病品種提供依據。

1 材料與方法

1.1 試驗材料

選取海島棉材料新海14 (母本, 感病), 06-146 (父本, 抗病)及其雜交后的F2:6RIL群體中的抗、感株系RIL10893 (超抗)、RIL10895 (高抗)、RIL10897 (感病)和RIL10796 (易感)共6份材料進行轉錄組測序。選取海島棉抗病材料06-146、埃棉2號(抗病)、DJ-07-136 (抗病)、5917 (抗病)、新海14、海92-3 (感病)和埃及棉424 (感病)共7份材料進行相關基因表達分析(表1和表2)。

1.2 種植及采樣方法

選取飽滿種子, 經75%乙醇沖洗3次, 雙氧水浸泡4~5 h后, 放入ddH2O里浸泡催芽24 h左右至露白, 將種子放在鋪有濕潤濾紙的發(fā)芽盒中, 28℃避光恒溫培養(yǎng)至芽長3 cm左右, 移至霍格蘭氏(Hoagland’s)營養(yǎng)液中, 在25℃, 光暗周期16 h/8 h的條件下培養(yǎng)。待第一片真葉完全展開時, 接種稀釋至1×106個mL–1的枯萎病菌(生理小種7號)。選取0 h及接菌后40 h的棉株下胚軸, 經液氮速凍處理, –80℃超低溫冰箱保存?zhèn)溆谩?/p>

1.3 轉錄組測序及差異表達基因分析

從0 h及接菌后4、10、18、28和40 h的棉株下胚軸提取RNA進行轉錄組測序, 得到的Reads用Trinity軟件進行無方向組裝, 將短的Reads拼裝成長的轉錄本, 用以鑒定轉錄本的亞型, 對unigene進行RPKM (Reads Per Kilobase per Million mapped reads, 每百萬reads中來自某基因每千堿基長度的reads數)處理, 再用DEG進行差異分析(篩選閾值為<0.005且|log2Fold Change|>1)。使用DAVID (The Database for Annotation, Visualization and Intergrated Discovery)在線平臺(https://dauid.ncifcrf.gov/)對這些基因進行功能聚類分析, 明確海島棉在抗病過程中涉及的相關代謝通路中主要基因, 以及這些基因在不同品種和不同接菌時期的表達差異。

表1 轉錄組測序材料

表2 實時熒光定量材料

1.4 實時熒光定量試驗及數據處理方法

將得到的差異表達序列, 運用NCBI中的ORF finder和BlastX找到序列的保守區(qū)域, 利用Primer 5.0設計qRT-PCR引物(表3)。使用天根生化科技有限公司的試劑盒(RNA plant Plus Reagent), 按照產品說明書提取總RNA。參照Thermo Fisher Scientific科技公司RNA純化試劑盒(Thermo Scientific Gene JET)的具體步驟對總RNA進行純化, 利用Titertek-Berthold檢測系統(tǒng)有限公司(Titertek-Berthold)的超微量分光光度計(Colibri)檢測總RNA濃度和質量, RNA濃度大于80 ngmL–1, 且A260/A280大于1.8, 采用賽默飛世爾(Thermo Fisher Scientific)公司逆轉錄試劑盒(RevertAid First strand cDNA Synthesis kit)反轉錄提取后的RNA, 用ddH2O稀釋cDNA樣品5倍后, 置–80℃?zhèn)溆谩Ｊ褂萌浇鹕锛夹g(TransGen Biotech)有限公司qRT-PCR試劑盒(TransStart Green qPCR SuperMix), 通過美國ABI公司Prism 7500 Fast System, 以為內參, 檢測7個品種(3個感病品種, 4個抗病品種), 6個接菌時期(0、4、10、18、28和40 h)的基因表達量變化。每個反應重復3次。根據Ct值, 利用公式2–DDCt計算基因的相對表達量。

2–DDCt=2–[(Ct目的基因–Ct內參基因)處理組–(Ct目的基因–Ct內參基因)對照組]

表3 熒光定量引物

1.5 病情指數的調查與計算

于2014年和2015年4月中旬將材料種植于新疆阿拉爾市新開嶺鎮(zhèn)枯萎病重病田中, 2次重復。在每年10月5日左右, 選擇連續(xù)5株棉株, 采用剖桿的方法考察海島棉的枯萎病抗性, 參照五級分類法[18]判斷棉花的發(fā)病級別并計算病情指數。

病情指數(DI) = ∑級數×每級病株數/(調查總株數×最高級值)×100

2 結果與分析

2.1 DEG分析

轉錄組測序分析DEG數據得出兩親本差異表達基因755個, 抗病親本上調表達的基因1185個。利用RIL群體中抗病材料RIL 10895和RIL 10893, 感病材料RIL 10897和RIL 10796作進一步分析(表4)。

表4 轉錄組測序中類黃酮代謝相關基因的表達

在表中出現(xiàn)同一基因的2個表達, 例如:基因為和, 實際上和序列均與擬南芥庫中的同一個基因相似, 這種現(xiàn)象是由比對到序列部位不同而導致, 但其表達趨勢均類似。因此后續(xù)基因分析過程中統(tǒng)一用代表。

The same gene has two forms, for example,gene is displayed asand, in fact,andare similar to the same sequence of TAIR, this phenomenon is caused by different parts of a sequence alignment. But its expression trends are similar. So in the following gene analysis process unified the display for.

接菌后, 類黃酮代謝途徑中6個關鍵基因(、、、、和)在抗病親本及RIL高抗, 超抗材料中均顯著上調表達(表2)?；蛟诳共∮H本及RIL超抗材料中上調表達的倍數分別為15.81和23.90, 比感病親本及RIL感病材料高出7倍和11倍。在抗病親本及RIL超抗材料中的表達量高出感病親本及RIL感病材料6倍之多。,基因在接菌前后表達量差異倍數的高低同棉花的抗性相一致,為RIL10893 (8.35) >RIL10895 (1.38) > RIL10897 (1.18) > RIL10796 (1.46),為RIL10893 (9.12) > RIL10895 (3.23) > RIL10897 (1.74) > RIL10796 (1.25);: RIL10893 (4.97) > RIL10895 (2.04) > RIL10897 (1.29) > RIL10796 (0.94)。暗示類黃酮代謝途徑相關基因可能與海島棉不同抗性表型的產生關系密切。

2.2 qRT-PCR結果分析

為進一步驗證類黃酮代謝途徑相關基因與海島棉枯萎病抗性的關系, 選取7個抗性表現(xiàn)不同的資源材料對類黃酮代謝通路中的8個基因進一步研究表明在接菌后18 h、40 h抗病材料的(花青素還原酶)表達量顯著高于感病材料, 且在這2個時間點, 感病材料的變異系數大于抗病材料, 說明該基因在感病材料之間的表達量有差異(圖1和表5)。除海92-3外, 另外2個感病材料在4、10和28 h、40 h的表達量均上調, 且4 h和10 h上升幅度明顯大于28 h和40 h, 在40 h所有材料的表達量均上調, 但是抗病材料的表達量顯著高于感病材料。推測該基因的表達量在感病材料中前期上調明顯, 在抗病材料中后期上調明顯。

(肉桂酸-4-羥化酶)(圖1和表5), 在接菌后40 h, 在抗病材料中的表達量顯著高于感病材料, 在感病材料中表達量的變異系數大于抗病材料, 說明在40 h感病材料之間的表達量也存在較大差異。感病材料中除了在4 h有上調表達外, 在其余接菌時間段均無明顯上調?？共〔牧铣鼶J-07-136外, 都在40 h擁有最大值。該基因在接菌后40 h抗病材料的瞬時表達量增加明顯。

(二氫黃酮醇還原酶)(圖1和表5), 在接菌脅迫后, 在抗病材料中的表達量均顯著或極顯著高于感病材料。除10 h外, 抗病材料在各時期表達量的變異系數均小于感病材料, 說明該基因的表達量在感病材料之間也存在差異。在感病材料中, 該基因的表達量除埃及棉424在4 h有少量上調外, 其余均為下調。而在抗病材料中, 該基因的表達量在4 h均上調表達, 并且在10 h持續(xù)上升, 在18 h全部下調, 在28 h除06-146外其余無明顯上調。該基因的表達量在抗病材料和感病材料中有顯著差異。在抗病材料中除DJ-07-136外, 其余均有上升下降在上升的趨勢, 且均在接菌后期(28 h和40 h)達到峰值, 推測該基因可能在海島棉抗枯萎病中起著重要作用。

(類黃酮3’羥化酶)(圖1和表5), 接菌脅迫后18、28和40 h在抗病材料中的表達量顯著高于感病材料, 在28 h和40 h的變異系數小于感病材料, 說明在這2個時間點該基因的表達量在感病材料之間也存在較大差異。從圖1可知, 除DJ-07-136外該基因的表達量在4~18 h均無明顯上調。在28 h, 新海14 (感病)、06-146 (抗病)、DJ-07-136 (抗病)、埃棉2號(抗病)表達量有大幅上升, 但新海14 (感病)的表達量顯著低于另外3種抗病材料, 且在這3個抗病材料中, 該基因的表達量在此時有最大值。該基因在感病材料中無明顯變化, 在抗病材料中后期明顯上調表達。

(查爾酮異構酶)(圖1和表5), 接菌脅迫后10、18、28和40 h, 在抗病材料中的表達量顯著高于感病材料, 但其變異系數小于感病材料。說明在感病材料之間該基因的表達水平有差異。抗病材料DJ-07-136和5917的表達量在18 h擁有最大值, 抗病材料06-146, 埃棉2號的表達量在40 h擁有最大值。該基因在感病材料中無明顯上調表達。推測該基因可能與海島棉抗枯萎病有關。

(查爾酮合成酶)(圖1和表5), 在接菌后18 h和40 h抗病材料的表達量顯著高于感病材料, 變異系數小于感病材料, 可能該基因的表達量在感病材料之間也存在差異。除06-146外, 該基因在抗病材料中4 h和10 h的表達量較18 h和28 h更加活躍。除埃棉2號外, 該基因在抗病材料中的表達量均有上升下降再上升的趨勢, 且第1個表達峰值出現(xiàn)在10 h, 第2個表達峰值在40 h。

(花色素合成酶)(圖1和表5), 在接菌后10、28和40 h抗病材料的表達量顯著高于感病材料, 且在28 h和40 h抗病材料表達量的變異系數小于感病材料, 說明該基因在接菌后28 h和40 h, 感病品種之間的表達水平也有差異。該基因在抗病材料中4 h和10 h的表達量較18 h和28 h更加活躍。在抗病材料中該基因的表達量均有上升下降再上升的趨勢, 第1個表達峰值在4 h和10 h, 第2個表達峰值在40 h。推測該基因在抗病材料中的表達趨勢可能是一個達到峰值后逐漸降低又回升的過程, 在感病材料中沒有這種趨勢。

(黃酮醇合成酶)(圖1和表5), 在接菌后10、28和40 h抗病材料的表達量顯著高于感病材料, 而其變異系數小于感病材料, 說明該基因的表達量在感病材料之間也存在差異。在接菌后該基因的表達量在感病材料中均有下調的趨勢。除DJ-07-136外, 抗病材料相對表達量的變化均有上升下降, 再上升的趨勢。

這8個基因在4個抗病材料中的表達量較感病材料上調更加明顯, 尤其是在40 h。除4 h外,基因在抗病材料中各個時期的表達量均顯著高于感病材料基因在抗病材料中各時期的相對表達量均顯著高于感病材料。的基因表達量在感病材料中無明顯變化, 在抗病材料中后期明顯上調表達。、、和在抗病材料中4 h和10 h的相對表達量高于18 h和28 h。和的表達量在抗病材料中均有上調后下調又回升的趨勢, 但在感病材料中并無這種趨勢。、和基因在感病材料接菌后4 h的表達量有明顯上調, 但是在4 h之后表達量上調不如抗病材料明顯。

2.3 抗、感材料間基因表達量方差分析

接菌40 h后, 8個基因在抗、感材料間均有極顯著差異(表6)。在除接菌后4 h以外的時期,基因的表達量在抗、感材料間均呈極顯著差異?；虻谋磉_量在抗、感材料之間, 除在4 h和28 h有顯著差異外, 在其他時期均有極顯著差異。在抗、感材料之間,基因的表達量在28 h和40 h有極顯著差異。在方差分析中發(fā)現(xiàn)抗病材料內部以及感病材料內部基因表達量的變異系數較大, 但是仍然發(fā)現(xiàn)抗、感材料之間的基因表達量有顯著差異, 說明抗、感材料之間的差異極顯著。

圖1 類黃酮代謝通路相關基因表達水平

表5 類黃酮代謝通路相關基因表達量在抗感材料之間的差異

表6 抗、感材料間基因表達量方差分析

*和**分別表示在<0.05和<0.01水平抗、感材料間差異顯著。

*and**indicate significant difference between resistant and susceptible genotypes at<0.05 and<0.01, respectively.

從表6可看出, 盡管在轉錄組測序的DEG分析中未檢測出和基因, 但是在接菌后40 h時8個基因的表達量在抗、感材料之間均差異極顯著, 這說明類黃酮途徑中的關鍵基因與海島棉枯萎病抗性關系密切。從表6還可看出, 在接菌后多個時間點和基因的表達量在抗、感材料之間均差異顯著。

2.4 海島棉病情指數與類黃酮代謝途徑相關基因表達量的相關性分析

將抗、感共7份材料2014—2015年的病情指數與類黃酮代謝通路相關基因表達量的相關性分析(表7)表明, 在接菌40 h后,,基因的表達量與病情指數呈極顯著負相關,和基因的表達量與病情指數呈顯著負相關。接菌28 h后,與基因的表達量與病情指數呈顯著負相關。接菌10 h后基因表達量與病情指數呈顯著負相關。從表7可以看出和基因分別在2個時期都被檢測出的基因表達量與病情指數之間顯著負相關, 說明這2個基因可能在海島棉枯萎病抗性中起著重要的作用。

表7 海島棉病情指數與類黃酮代謝途徑相關基因表達量的相關系數

*和**分別表示在<0.05和<0.01水平差異顯著。

*and**indicate significant difference at<0.05 and<0.01, respectively.

3 討論

因為棉花黃萎病菌孢子萌發(fā)并滲透植物表皮細胞的過程發(fā)生在真菌感染后的前12 h[19]。復雜的病菌感染反應以及信號的傳導和交換通常發(fā)生在真菌感染12~40 h[20-21]。所以本研究選取接菌后6個時間點(0、4、10、18、28和40 h)進行熒光定量分析。由于棉花在接菌40 h后會激發(fā)細胞凋亡等生理生化反應[22]。為了避免凋亡基因等影響對抗病基因的判斷, 本研究不對40 h以后的材料進行分析。

改良海島棉抗枯萎病性是提高新疆海島棉生產的根本途徑。通過揭示海島棉的抗病機制可以快速培育抗病品種。體外施加類黃酮化合物可顯著提高植物的抗病性[23-24]。研究者們已通過RNA-seq的方法證實類黃酮代謝相關基因與馬鈴薯、油棕和向日葵抗病有關[25-27]。前人對于棉花類黃酮的研究也主要集中在棉纖維顏色方面[28]。而對于類黃酮與海島棉枯萎病抗性的關系研究甚少。Sun等[29]在對接種黃萎病菌前后陸地棉及海島棉的轉錄組數據分析后指出, 在海島棉中苯丙烷代謝途徑相關基因呈上調表達趨勢。苯丙烷代謝通路在棉花抵抗病原菌方面的積極作用已被許多研究所證實[11,30-32]。類黃酮代謝途徑(圖2)[9]是苯丙烷代謝途徑通過查爾酮合成酶()延伸的另一個重要分支[13]。本研究RNA-seq分析表明, 在接菌后40 h, 類黃酮代謝途徑6個關鍵基因(和)的表達量在抗感材料之間差異顯著。其中和基因在抗病親本接菌前后表達量的差異倍數高達15.81、12.99、13.30和11.54。說明類黃酮代謝途徑與海島棉抗病性存在密切聯(lián)系。

雖然研究者們通過多種方法證明類黃酮代謝途徑與植物抗病有關, 且關鍵的轉錄因子以及調控基因也在多種植物中被克隆, 但由于類黃酮代謝途徑比較復雜, 相關基因大多數為多基因家族。本研究在RNA-seq的基礎上進行DEG分析, 檢測到類黃酮代謝通路相關基因與海島棉枯萎病抗性密切相關, 利用熒光定量的方法對不同抗性資源材料進一步分析表明, 接菌40 h后這些基因在抗病材料中的表達量顯著高于感病材料。進一步證實類黃酮代謝通路對海島棉抵抗枯萎病起積極的作用。

在除40 h外的其他時間點,和基因的表達量在抗、感材料間也存在顯著或極顯著差異,和基因在接菌后2個時間點的表達量與病情指數有顯著或極顯著負相關。其中基因的表達量在接菌后的4個時間點均有極顯著差異?；蚴穷慄S酮代謝途徑中的一個關鍵酶, 在轉基因煙草中活性與總黃酮含量呈正相關[33]。推測在抗病材料中基因受病菌誘導后表達量上調, 使黃酮含量增加, 而類黃酮化合物與相關酚類物質能夠保護植物免受微生物浸染[34], 所以黃酮含量的增加可以導致棉花獲得抗性。在接菌后各個時間點基因的表達量在抗感材料之間均有顯著或者極顯著差異。根據馬春雷等[35]的研究結果可知,和基因的表達量與茶樹兒茶素質量分數呈一定的正相關?；虻母弑磉_可能會導致兒茶素含量增加。前人研究發(fā)現(xiàn), 在枯萎病菌浸染后棉苗組織中兒茶素質量分數明顯升高, 且在抗病品種棉苗組織中的兒茶素質量分數高于感病品種, 另外還發(fā)現(xiàn)棉花體內的兒茶素可抑制產孢及孢子的萌發(fā)[36]。兒茶素的增加可能會抑制枯萎病菌的菌絲生長。以上兩基因CHI后續(xù)可通過沉默該基因后測定黃酮含量來探究該基因的抗病機制。后續(xù)可以結合兒茶素含量測定或者將正反義基因轉入棉花進一步研究基因的功能。

圖2 棉花中的類黃酮途徑

Phenylalanine: 苯丙氨酸; Cinnamate: 肉桂酸;-Cinnamate: 對香豆酸;-Cinnamate-CoA: 香豆酸輔酶A; Chalcone: 查爾酮; Flavanone: 黃烷酮; Isoflavones: 異黃酮; Isoflavanones: 二氫異黃酮類; Dihydro flavonols: 二氫黃酮醇; Flavonols: 黃酮醇; Anthocyanins: 花青素; Anthocyanins: 花青素類; Procyanidins: 原花青素; Flavan-4-ols: 黃烷-4-醇; Flavones: 黃酮類; PAL: phenylalanine ammonia-lyase, 苯丙氨酸解氨酶; C4H: cinnamate 4-hydroxylase, 肉桂酸羥化酶; 4CL: 4-coumarate:CoA ligase, 4-香豆酸: 輔酶A連接酶; CHS: chalcone synthase, 查爾酮合酶; CHI: chalcone isomerase, 查爾酮異構酶; IFS: isoflavone synthase, 異黃酮合酶; IFR: isoflavone reductase, 異黃酮還原酶; F3’H: lavanone 3’-hydroxylase, 類黃酮3’羥化酶; FLS: flavonol synthase, 黃酮醇合酶; DFR: dihydroflavonol-4-reductase, 二氫黃酮醇還原酶; ANS: anthocyanidin synthase, 花青素合酶; ANR: anthocyanidin reductase, 花青素還原酶; UFGT, UDP-flavonoid glucosyltransferase, UDP類黃酮糖基轉移酶。粗體字示本文所研究的基因。Bold is the gene studied in this article.

Shih等[37]在研究高粱的抗病機制時認為,基因在接菌后24 h表達量開始上調, 在接菌后72 h表達量達到最大值。本研究與其類似的是基因的表達量均在接菌后期上調明顯。另外研究結果還指出在接菌后期,基因在抗病材料中的表達量顯著高于感病材料。這個表達量差異說明其與抗病密切相關。轉煙草與正常煙草相比苯丙氨酸解氨酶(phenylalanineammonia-lyase, PAL)和過氧化物酶(peroxidase, POD)活性增加, 可減輕病原微生物對煙草自身產生的傷害[38]。另外,基因的過量表達可以導致具有抗氧化和抗菌活性的類黃酮化合物槲皮素(Quercetin)的含量增加[23-24]。因此推測在抗病材料中的高表達可能導致PAL和 POD活性增加或者槲皮素含量增加, 從而提高海島棉枯萎病抗性。

關于類黃酮代謝途徑相關基因與植物抗病性的關系, 前人的研究主要集中在陸地棉中。Deshika等[39]研究證明棉花在接種枯萎病菌后24 h基因的表達量與對照并無差異。本研究證實基因在接菌后40 h表達量上調明顯, 而在28 h大多數品種并無明顯上調表達。該基因在海島棉和陸地棉中的差異表達暗示海島棉抗枯萎病機制不同。Xu等[9]以接種黃萎病菌前后的海7124為材料進行轉錄組測序, 并利用qPCR的方法對差異表達基因進行分析, 發(fā)現(xiàn)抗病材料中基因和基因在接菌后各個時期(4、12、24和48 h)的表達量均高于感病材料, 且基因的表達量要明顯高于基因的表達量。另外Xu等[9]還發(fā)現(xiàn)和基因的表達量均在在接菌后24 h達峰值, 在48 h表達量有少量下降。而本研究中這2個基因的表達量都在40 h有最大值。在接種黃萎病菌和枯萎病病菌后和基因的差異表達說明枯萎病與黃萎病的抗性機制也存在差異。

4 結論

抗病材料中類黃酮代謝通路相關基因的表達量顯著高于感病材料。在接菌后40 h, 類黃酮代謝途徑8個基因在抗病材料中的表達量均顯著高于感病材料。其中、在抗病材料中其他時間點的表達量也顯著高于感病材料。在接菌后的兩個時間點,和基因的表達量與病情指數呈顯著負相關。在接菌后40 h,基因的表達量與病情指數呈極顯著負相關。推測、和可能與海島棉抗枯萎病有關。

[1] 孔慶平. 我國海島棉生產概況及比較優(yōu)勢分析. 中國棉花, 2002, 12: 19–23 Kong Q P. Production situation and comparative advantage analysis ofin China., 2002, 12: 19–23 (in Chinese with English abstract)

[2] 徐秋華, 張獻龍, 聶以春, 馮純大. 我國棉花抗枯萎病品種的遺傳多樣性分析. 中國農業(yè)科學, 2002, 35: 272–276 Xu Q H, Zhang X L, Nie Y C, Feng C D. Genetic diversity evaluation of cultivars (L.) resistant to fusarium wilt by RAPD markers., 2002, 35: 272–276 (in Chinese with English abstract)

[3] 彭姜龍, 曲延英, 王莉萍, 王聰, 郝維維, 汪鋮亮. 不同溫度條件下接種棉花枯萎病菌海島棉和陸地棉抗病性研究. 新疆農業(yè)科學, 2013, 50: 89–93 Peng J L, Qu Y Y, Wang L P, Wang C, Hao W W, Wang C L. A study on properties of resistance to Fusarium wilt in leaves of sea island cotton () and upland cotton (L.) under different temperature conditions., 2013, 50: 89–93 (in Chinese with English abstract)

[4] 校百才. 陸地棉抗枯、黃萎病性狀配合力、遺傳力的初步研究. 作物學報, 1985, 11: 267–273 Xiao B C. Preliminary study on cotton Fusarium and Verticillium wilt resistance traits, combining ability and heritability., 1985, 11: 267–273 (in Chinese with English abstract)

[5] Smith A L, Dick J B. Inheritance of resistance to Fusarium wilt in upland and sea island cottons as complicated by nematodes under field conditions., 1960, 50: 44–48

[6] Yang J, Ma Q, Zhang Y, Wang X F, Zhang G Y, Ma Z Y. Molecular cloning and functional analysis of GbRVd, a gene inthat plays an important role in conferring resistance to Verticillium wilt., 2016, 575: 687–694

[7] Gao X Q, Wheeler T, Li Z H, Charles M K, He P, Shan L B. Silencing GhNDR1 and GhMKK2 compromises cotton resistance to Verticillium wilt., 2011, 66: 293–375

[8] Snyder M. RNA-seq: a revolutionary tool for transcriptomics., 2009, 10: 57–63

[9] Xu L, Zhu L F, Tu L L, Liu L L, Yuan D J, Jin L, Long L, Zhang X L. Lignin metabolism has a central role in the resistance of cotton to the wilt fungusas revealed by RNA-Seq-dependent transcriptional analysis and histochemistry., 2011, 62: 5607–5621

[10] 王春暉. 陸地棉抗病自交系在黃萎病菌誘導下的轉錄組測序研究. 中國農業(yè)科學院碩士學位論文, 北京, 2013 Wang C H. Study on Transcriptome Sequencing of Cotton (L.) Resistant Inbred Line Induced by Verticillium. MS Thesis of Chinese Academy of Agricultural Sciences, Beijing, China, 2013 (in Chinese with English abstract)

[11] Hichri I, Barrieu F, Bogs J, Kappel C, Delrot S, Lauvergeat V. Recent advances in the transcriptional regulation of the flavonoid biosynthetic pathway., 2011, 62: 2465–2483

[12] Malhotra B, Onyilagha J C, Bohm B A, Towers G H N, James D, Harborne J B, French C J. Inhibition ofby flavonoids., 1996, 43: 1271–1276

[13] 賈振華. 植物黃酮類化合物槲皮素與轉錄因子誘導和調控植物防衛(wèi)反應的研究. 南京農業(yè)大學博士學位論文, 江蘇南京, 2010 Jia Z H. Study on the Effects of Quercetin and Transcription Factoron Plant Defense Response. PhD Dissertation of Nanjing Agricultural University, Nanjing, China, 2010 (in Chinese with English abstract)

[14] Anna K P, Edward F. S, Jaroslav P. Effect of flavonoids on mycelial growth of., 1995, 23: 683–693

[15] 左濤, 趙樹堂, 盧孟柱, 孫愛東, 王延偉, 賀偉. 楊樹二氫黃酮醇-4-還原酶基因(DFR)的克隆及反義表達對兒茶素合成的影響. 東北林業(yè)大學學報, 2016, 10: 49–55 Zuo T, Zhao S T, Lu M Z, Sun A D, Wang Y W, He W. Cloning dihydroflavonol-4-reductase gene (DFR) of poplar and its antisense expression effects on catechin synthesis., 2016, 10: 49–55 (in Chinese with English abstract)

[16] 馬銀平, 王付欣, 楊淳淋, 沈法富, 夏桂先. 海島棉幾丁質酶基因的克隆與功能分析. 遺傳, 2012, 34: 240–247 Ma Y P, Wang F X, Yang C L, Shen F F, Xia G X. Cloning and functional analysis of chitinase genefrom sea-island cotton ().(Beijing), 2012, 34: 240–247 (in Chinese with English abstract)

[17] Leslie A W, Li G Q, Doreen W, Imre E S, Keith R D. The phenylalanine ammonialyase gene family in., 1995, 27: 327–338

[18] Wajad N, Abdul L T, Saghir A, Khalid M, Abid M, Zhou B L. Evaluation of cotton leaf curl virus resistance in BC1, BC2, and BC3progenies from an interspecific cross betweenand., 2014, 9: e111861

[19] Emilie F F, Bart P H J T. Physiology and molecular aspects of Verticillium wilt diseases caused byand., 2006, 7: 71–86

[20] Zhang S, Klessig D F. MAPK cascades in plant defense signaling., 2001, 6: 520–527

[21] Kunkel B N, Brooks D M. Cross talk between signaling pathways in pathogen defense., 2002, 5: 325–331

[22] Patil M A, Pierce M L, Phillips A L, Venters B J, Essenberg M. Identification of genes up-regulated in bacterial-blight-resistant upland cotton in response to inoculation withpv.., 2005, 67: 319–335

[23] Atsushi N, Chikara M, Mineo S, Hideyuki M, Atsushi K, Hong J S, Keisuke K, Jun A, Akira K. Functional analysis of soybean genes involved in flavonoid biosynthesis by virus-induced gene silencing., 2007, 5: 778–790

[24] Boots A W, Haenen G R M M, Bast A. Health effects of quercetin: from antioxidant to nutraceutical., 2008, 585: 325–337

[25] 包改紅, 畢陽, 李永才, 王毅, 王婷, 唐瑛, 馬朝玲, 白小東. 硫色鐮刀菌接種對抗病/易感品種馬鈴薯塊莖苯丙烷代謝的影響比較. 食品科學, 2015, 36: 251–256 Bao G H, Bi Y, Li Y C, Wang Y, Wang T, Tang Y, Ma Z L, Bai X D. Comparison of phenylpropanoid pathway metabolism in slices of susceptible and resistant potato cultivars inoculated with., 2015, 36: 251–256 (in Chinese with English abstract)

[26] Guo S, Zuo Y, Zhang Y, Wu C, Su W, Jin W, Yu H, An Y, Li Q. Large-scale transcriptome comparison of sunflower genes responsive to., 2017, 18: 42

[27] Tan B A, Daim L D J, Ithnin N, Ooi T E K, Md-Noh N, Mohamed M, Mohd-Yusof H, Appleton D R, Kulaveerasingam H. Expression of phenylpropanoid and flavonoid pathway genes in oil palm roots during infection by., 2016, 7: 11–20

[28] Sun W G, Xu Y. Study on the composition of the pigments in the nature colored cotton., 2009, 23: 119–124

[29] Sun Q, Jiang H Z, Zhu X Y, Wang W N, He X H, Shi Y Z,Yuan Y L, Du X M, Cai Y F. Analysis of sea-island cotton and upland cotton in response toinfection by RNA sequencing., 2013, 14: 852

[30] Federico P, Marta N, María A B, Fuencisla M, Barceló A R. Changes in stem lignins (monomer composition and crosslinking) and peroxidase are related with the maintenance of leaf photosynthetic integrity during Verticillium wilt in., 2004, 163: 111–123

[31] Gayoso C, Pomar F, Novo-Uzal E, Merino F, Martínez D I O. The Ve-mediated resistance response of the tomato toinvolves H2O2, peroxidase and lignins and drivesgene expression., 2010, 10: 1–19

[32] Irina O V, Robert T R, Jane E G. David J R, Kent A M, Nicholi V. Characterization of flavonols inpowder., 2003, 52: 188–195

[33] Li F, Jin Z, Qu W, Zhao D, Ma F. Cloning of a cDNA encoding thechalcone isomerase and its expression in transgenic tobacco., 2006, 44: 455–461

[34] Fawe A, Abou Z M, Menzies J G, Belanger R R. Silicon-mediated accumulation of flavonoid phytoalexins in cucumber., 1998, 88: 396–401

[35] 馬春雷, 喬小燕, 陳亮. 茶樹無色花色素還原酶基因克隆及表達分析. 茶葉科學, 2010, 30: 27–36 Ma C L, Qiao X Y, Chen L. Cloning and expression analysis ofreducase gene of tea plant ()., 2010, 30: 27–36 (in Chinese with English abstract)

[36] 宋鳳鳴, 鄭重, 童賢明, 葛秀春. 兒茶素對棉枯萎病菌胞壁降解酶的抑制及在棉花抗病性中的作用. 真菌學報, 1996, 4: 297–303 Song F M, Zheng Z, Tong X M, Ge X C. Inhibition of tea catechins on cell wall degrading enzymes of cottonwilt and its role in disease resistance of cotton., 1996, 4: 297–303 (in Chinese with English abstract)

[37] Shih C H, ChuI K, Yip W K, Lo C. Differential expression of two flavonoid 3’-hydroxylase cDNAs involved in biosynthesis of anthocyanin pigments and 3’-deoxyanthocyanidin phytoalexins in sorghum., 2006, 47: 1412–1419

[38] 張松煥, 李春奇, 郭惠明, 裴熙祥, 程紅梅. 過量表達紫莖澤蘭類黃酮3’-羥化酶基因對轉基因煙草POD、PAL活性的影響. 中國農業(yè)科技導報, 2009, 11(3): 98–101 Zhang S H, Li C Q, Guo H M, Pei X X, Cheng H M. Effects of eupatorium adenophorum flavonoid 3'-hydroxylase overexpression on pod and pal activity in transgenic tobacco., 2009, 11(3): 98–101 (in Chinese with English abstract)

[39] Kohli D, Joshi G, Deokar A A, Bhardwaj A R, Agarwal M, Katiyar-Agarwal S, Srinivasan R, Jain P K. Identification and characterization of wilt and salt stress-responsive microRNAs in chickpea through high-throughput sequencing., 2014; 9: e108851

Correlation between Fusarium wilt Resistance and Expression Levels of Genes Involved in Flavonoid Metabolism Pathway inL.

HUANG Qi-Xiu, QU Yan-Ying, YAO Zheng-Pei, LI Meng-Yu, and CHEN Quan-Jia*

College of Agronomy / Key Laboratory of Agriculture Biological Technology, Xinjiang Agricultural University, Urumqi 830052, China

Fusarium wilt is one of the important factors threatening the production of island cotton (L.). Understand the molecular mechanism of Fusarium wilt resistance will facilitate cotton breeding against Fusarium wilt and therefore alleviate the problem in cotton production. On the basis of RNA-seq results in previous studies, differentially expressed genes (DEG) were further analyzed in this study. After inoculating the pathogenic fungus for 0–40 hours, the DEGs expression levels of seven varieties in different wilt resistance were studied using qRT-PCR method. The correlation between gene expression level and disease index was analyzed. The flavonoid biosynthetic pathway genes were found to be related to Fusarium wilt resistance. The qRT-PCR result indicated the significantly higher expression levels of flavonoid metabolic pathway genes in resistant genotypes than those susceptible genotypes. The expression levels of,, and, the key genes in flavonoid metabolic pathways, were significantly higher in resistant genotypes than in susceptible genotypes (< 0.05 or< 0.01) at several time points after inoculation. Particularly, the expression levels ofandwere negatively correlated with disease index. Therefore, we conclude that the flavonoid metabolic pathway related genes might be involved in the resistance to Fusarium wilt in cotton, among which,, andare the key genes.

; Fusarium wilt; Flavonoids; Expression analysis; RNA-seq

10.3724/SP.J.1006.2017.01791

本研究由國家自然科學基金項目(31560409)資助。

This study was supported by the National Natural Science Foundation of China (31560409).

陳全家, E-mail: chqjia@126.com第一作者聯(lián)系方式: E-mail: 13899525177@163.com, Tel: 13899525177

2017-03-09; Accepted(接受日期): 2017-07-23; Published online(網絡出版日期): 2017-08-10.

http://kns.cnki.net/kcms/detail/11.1809.S.20170810.1616.010.html