湯亞琪，蘇文文，李 瑩，王保通，李 強(qiáng)

(西北農(nóng)林科技大學(xué)植物保護(hù)學(xué)院/旱區(qū)作物逆境生物學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西楊凌 712100)

小麥條銹病是世界范圍內(nèi)小麥主產(chǎn)區(qū)最重要的病害之一[1]，其病原菌條形柄銹菌小麥?；?PucciniastriiformisWest. f. sp.triticiEriks & Henn，Pst)為專性寄生菌[2]，生理小種致病性變異迅速且頻繁，容易造成主栽小麥品種抗條銹的頻繁“喪失”[3]，導(dǎo)致小麥條銹病的爆發(fā)流行，嚴(yán)重威脅著全球小麥的安全生產(chǎn)[4]。

中4是以普通小麥(TriticumaestivumL.)(2n=6x=42，AABBDD)與中間偃麥草(Thinopyrumintermedium)(2n=6x=42，E1E1E2E2XX)經(jīng)有性雜交、選育獲得的異源八倍體小偃麥(AABBDDXX)[5]。自1983年列為中國(guó)小麥條銹菌鑒別寄主至今，對(duì)我國(guó)現(xiàn)有小麥條銹菌流行小種均表現(xiàn)為抗病[6]，其抗性穩(wěn)定且持久，具有較高的研究和利用價(jià)值。但由于中4是異源八倍體小偃麥，與六倍體普通小麥經(jīng)常規(guī)雜交后獲得穩(wěn)定遺傳后代的難度較大，導(dǎo)致其抗性研究嚴(yán)重滯后[7]。藺瑞明等[7]將中4分別與中國(guó)春phb1突變體和條銹病感病品種銘賢169雜交，對(duì)其雜交后代進(jìn)行遺傳分析，發(fā)現(xiàn)其對(duì)條銹菌小種CYR31和CYR32的抗病性由1對(duì)顯性基因控制。曹世勤等[8]利用基因推導(dǎo)法發(fā)現(xiàn)中4含有未知抗條銹病基因。甘肅省天水市農(nóng)業(yè)科學(xué)研究所以中4為親本，先后選育出中梁22號(hào)[9]、中93444[10]、中93447[11]等多個(gè)具有較好豐產(chǎn)性和抗病性的品種(系)材料。但是其抗條銹病機(jī)制目前尚未見報(bào)道。

本研究利用轉(zhuǎn)錄組測(cè)序，篩選中4接種條銹菌后的差異表達(dá)基因，以期為進(jìn)一步解析中4抗條銹性分子機(jī)制提供科學(xué)依據(jù)。

1 材料與方法

1.1 供試材料

供試小麥品種為我國(guó)小麥條銹菌鑒別寄主中4和條銹病感病品種銘賢169，小麥條銹菌為我國(guó)小麥條銹菌流行小種CYR32，均由西北農(nóng)林科技大學(xué)植物保護(hù)學(xué)院小麥病原真菌監(jiān)測(cè)及抗病遺傳實(shí)驗(yàn)室提供。

1.2 樣品采集

挑選籽粒發(fā)育良好、顏色形狀一致的中4種子，催芽后播種在7 cm × 7 cm的小花盆內(nèi)。待小麥幼苗長(zhǎng)至第一葉完全展開時(shí)，接種條銹菌CYR32，接種0、12、24和48 h后，分別采集接種葉片樣品，以未接種葉片(0 h)為對(duì)照，每個(gè)時(shí)間點(diǎn)采集3個(gè)生物學(xué)重復(fù)，用于轉(zhuǎn)錄組數(shù)據(jù)質(zhì)量測(cè)定。接種48 h后的材料在人工氣候箱(溫度15～17 ℃、濕度90%、16 h光照/8 h黑暗、光照強(qiáng)度10 000 lx)繼續(xù)培養(yǎng)16 d，用于觀察表型，以銘賢169作為感病對(duì)照。

1.3 RNA提取及文庫構(gòu)建

樣品總RNA的提取按照Plant RNA Midi Kit(OMEGA，美國(guó))試劑盒說明書進(jìn)行，獲取總RNA后，檢測(cè)RNA樣品的純度、濃度和完整性，用帶有Oligo(dT)的磁珠富集真核生物mRNA，加入Fragmentation buffer將mRNA隨機(jī)打斷，以First Strand cDNA Synthesis Kit(OMEGA，美國(guó))合成第一條cDNA鏈，然后加入緩沖液、dNTPs、RNase H和DNA polymerase I合成第二條cDNA鏈，利用AMPure XP beads純化雙鏈cDNA，純化后的雙鏈cDNA須先進(jìn)行末端修復(fù)、加A尾并連接測(cè)序接頭，再用AMPure XP beads進(jìn)行片段大小選擇，最后通過PCR富集得到 cDNA文庫。

1.4 轉(zhuǎn)錄組測(cè)序及功能注釋

通過Illumina高通量測(cè)序獲得12個(gè)cDNA文庫，用FastQC對(duì)cDNA文庫的堿基質(zhì)量值(Q30)、堿基含量分布(GC含量)和重復(fù)相關(guān)性(r2)等指標(biāo)進(jìn)行檢測(cè)，以判斷文庫數(shù)據(jù)質(zhì)量并評(píng)估其是否滿足后續(xù)分析。用HISAT2對(duì)reads進(jìn)行高效比對(duì)，用StringTie進(jìn)行從頭組裝，基于中國(guó)春全基因組(https://wheat-urgi.versailles.inra.fr/Seq-Repository)進(jìn)行有參組裝。

將完成組裝的轉(zhuǎn)錄本數(shù)據(jù)用非冗余(non-redundant，NR)數(shù)據(jù)庫和基因本體(gene ontology，GO)數(shù)據(jù)庫進(jìn)行比對(duì)注釋，用真核生物同源蛋白簇(eukaryotic orthologous groups，KOG)、蛋白質(zhì)直系同源簇(cluster of orthologous groups of proteins，COG)以及基因的進(jìn)化譜系：非監(jiān)督直系群(evolutionary genealogy of genes:non-supervised orthologous groups，eggNOG)等數(shù)據(jù)庫進(jìn)行基因富集聚類，用京都基因與基因組百科全書(kyoto encyclopedia of genes and genomes，KEGG)數(shù)據(jù)庫進(jìn)行代謝通路的富集分析。

1.5 實(shí)時(shí)熒光定量PCR

選取轉(zhuǎn)錄組測(cè)序得到的差異表達(dá)基因，根據(jù)其全長(zhǎng)cDNA序列，使用Primer 5設(shè)計(jì)特異性引物(表1)，以小麥TaEF為內(nèi)參基因。提取中4各個(gè)樣品的總RNA并反轉(zhuǎn)錄為cDNA，按照ChamQTMSYBR○RqPCR Master Mix(Vazyme，中國(guó))說明書的反應(yīng)體系和條件，使用QuantStudio 5(ABI，美國(guó))進(jìn)行qRT-PCR，按照2-△△Ct計(jì)算基因相對(duì)表達(dá)量，重復(fù)3次。

表1 熒光定量PCR所用引物Table 1 Primers used for fluorescence quantitative PCR

2 結(jié)果與分析

2.1 中4抗條銹菌鑒定結(jié)果

從圖1可以看出，中4葉片只有少許壞死斑且不連片，表現(xiàn)為高抗，而銘賢169產(chǎn)生大量夏孢子堆。

2.2 轉(zhuǎn)錄組數(shù)據(jù)質(zhì)量

取接種條銹菌0、12、24和48 h后的中4葉片樣品，每個(gè)處理3個(gè)重復(fù)，共計(jì)12個(gè)樣品進(jìn)行轉(zhuǎn)錄組測(cè)序分析，共獲得130.51 Gb Clean Data，各樣品Clean Data均達(dá)到8.42 Gb；堿基質(zhì)量值Q30的堿基百分比為92.83%～94.14%，表明堿基識(shí)別可靠，準(zhǔn)確度高；G和C、A和T的含量為58.33%～60.83%。分別將各樣品的Clean Reads與中國(guó)春全基因組IWGSC RefSeq v1.0進(jìn)行序列比對(duì)，比對(duì)效率為60.29%～69.03%，由于中4為異源八倍體小偃麥(AABBDDXX)[6]，其X染色體組與普通小麥親緣關(guān)系較遠(yuǎn)，所以轉(zhuǎn)錄組的比對(duì)效率沒有出現(xiàn)高度匹配。基于比對(duì)結(jié)果進(jìn)行基因表達(dá)量分析、功能注釋、富集分析以及新基因的發(fā)掘[8]，發(fā)掘新基因20 788個(gè)，其中 15 069個(gè)得到功能注釋。關(guān)于樣品相關(guān)性，r2值越接近1，說明兩個(gè)重復(fù)樣品相關(guān)性越強(qiáng)且樣品重復(fù)性良好，測(cè)序結(jié)果顯示，中4接種條銹菌CYR32 0 h后，3個(gè)樣品的r2值為0.960～ 0.983；接種條銹菌CYR32 12 h后，3個(gè)樣品的r2值為0.958～0.975；接種條銹菌CYR32 24 h后，3個(gè)樣品的r2值為0.972～ 0.982；而接種條銹菌CYR32 48 h后，3個(gè)樣品的r2值為0.824～0.932，表明樣品間整體相關(guān)性較強(qiáng)，重復(fù)性良好。

2.3 差異表達(dá)基因數(shù)目及染色體定位

以Fold Change≥2且FDR<0.01為篩選標(biāo)準(zhǔn)[12]，對(duì)不同時(shí)間點(diǎn)差異表達(dá)基因的分布進(jìn)行調(diào)查，發(fā)現(xiàn)0 h VS 12 h(表示兩個(gè)侵染時(shí)間點(diǎn)差異表達(dá)基因的比較，下同)、0 h VS 24 h、0 h VS 48 h分別篩選出7 816、10 412、9 563個(gè)上調(diào)表達(dá)差異基因和7 209、10 159、9 612個(gè)下調(diào)表達(dá)差異基因；12 h VS 24 h、12 h VS 48 h、24 h VS 48 h分別篩選出937、1 150、57個(gè)上調(diào)表達(dá)差異基因和 2 866、 1 643、229個(gè)下調(diào)表達(dá)差異基因。結(jié)果表明，以0 h為對(duì)照，24 h獲得的差異表達(dá)基因最多，48 h次之，12 h最少；以12 h為對(duì)照，24 h獲得的差異表達(dá)基因多于48 h；48 h較24 h新增加了286個(gè)差異表達(dá)基因。

為獲得潛在的抗病相關(guān)基因，本研究重點(diǎn)關(guān)注上調(diào)表達(dá)差異基因，從圖2可以看出，共獲得13 674個(gè)上調(diào)表達(dá)差異基因，12 h與24 h、48 h共有的上調(diào)表達(dá)差異基因分別有6 241、5 611個(gè)，24 h與48 h共有的上調(diào)表達(dá)差異基因7 416個(gè)，0～48 h篩選到持續(xù)上調(diào)表達(dá)的差異基因5 151個(gè)，這些基因可能與中4的的抗條銹性有關(guān)。對(duì)5 151個(gè)上調(diào)差異基因進(jìn)行染色體定位，其中 4 482個(gè)獲得定位信息，其余669個(gè)無法確定其染色體位置，為新基因。進(jìn)一步分析發(fā)現(xiàn)，差異表達(dá)基因在B染色體組上分布最多，為1 581個(gè)；D染色體組次之，為1 523個(gè)；A染色體組最少，為 1 378個(gè)。上調(diào)表達(dá)差異基因在2B(7.12%)、5B (6.38%)、3B(6.34%)、2D(5.85%)、7D(5.42%)染色體上所占比例較高。

2.4 共同上調(diào)表達(dá)差異基因的分類注釋及富集

2.4.1 共同上調(diào)表達(dá)差異基因的GO功能注釋

中4轉(zhuǎn)錄組測(cè)序數(shù)據(jù)中有58 059個(gè)轉(zhuǎn)錄本在GO數(shù)據(jù)庫中被注釋，5 151個(gè)共同上調(diào)表達(dá)差異基因中3 411個(gè)基因獲得GO注釋。獲得注釋的轉(zhuǎn)錄本可分為分子功能、細(xì)胞組分和生物學(xué)過程3大類，46個(gè)亞類(圖3)。分子功能包含13個(gè)亞類，其中蛋白結(jié)合、催化活性通路富集到的上調(diào)表達(dá)差異基因數(shù)目較多，此外，在轉(zhuǎn)運(yùn)活性、受體活性、酶調(diào)節(jié)活性、抗氧化活性、蛋白結(jié)合轉(zhuǎn)錄因子活性等通路也富集到上調(diào)表達(dá)差異基因；細(xì)胞組分包含13個(gè)亞類，其中細(xì)胞部件、細(xì)胞、細(xì)胞器等通路富集到的上調(diào)表達(dá)差異基因數(shù)目較多，在大分子復(fù)合物、細(xì)胞連接、超分子復(fù)合物等通路也富集到上調(diào)表達(dá)差異基因；生物學(xué)過程包含20個(gè)亞類，其中細(xì)胞過程、代謝過程、單有機(jī)體過程等通路富集到的上調(diào)表達(dá)差異基因數(shù)目較多，此外，在生殖過程、多有機(jī)體過程、信號(hào)轉(zhuǎn)導(dǎo)等通路也富集到上調(diào)表達(dá)差異基因。

2.4.2 共同上調(diào)表達(dá)差異基因的功能富集分析

利用COG、KOG、eggNOG數(shù)據(jù)庫對(duì)5 151個(gè)上調(diào)表達(dá)差異基因分別進(jìn)行直系同源分類，以期獲得更準(zhǔn)確的分類及功能注釋。COG數(shù)據(jù)庫中有2 638個(gè)上調(diào)表達(dá)差異基因獲得注釋分類，最大類別為R類(僅一般功能預(yù)測(cè))，包含555個(gè)基因，占21.04%；其次是K類(轉(zhuǎn)錄)，包含306個(gè)基因，占11.60%(圖4A)；KOG數(shù)據(jù)庫中有 2 486個(gè)上調(diào)表達(dá)差異基因獲得注釋分類，最大類別為R類，包含436個(gè)基因，占17.54%，其次是O類(翻譯后修飾，蛋白質(zhì)轉(zhuǎn)換，分子伴侶)，包含271個(gè)基因，占10.90%(圖4B)；eggNOG數(shù)據(jù)庫中有3 906個(gè)上調(diào)表達(dá)差異基因獲得注釋分類，最大類別為S類(未知功能)，包含859個(gè)基因，占21.99%，其次是R類，包含666個(gè)基因，占 17.05%(圖4C)。此外，COG、KOG、eggNOG三個(gè)數(shù)據(jù)庫分別注釋到與防御相關(guān)的V類(防御機(jī)制)上調(diào)表達(dá)差異基因，分別為47、28、25個(gè)。雖然eggNOG數(shù)據(jù)庫獲得的注釋基因最多，但是注釋到與防御相關(guān)的基因最少，在未知功能和一般功能預(yù)測(cè)富集到的基因最多，為1 525個(gè)，占 39.04%；KOG數(shù)據(jù)庫在所占比例為 5.83%～ 5.39%處出現(xiàn)分類富集，且數(shù)據(jù)庫25個(gè)分類中均注釋到上調(diào)表達(dá)差異基因；COG和KOG數(shù)據(jù)庫雖然在最大類別注釋相同，但是其余類別差異很大。三個(gè)數(shù)據(jù)庫注釋分類結(jié)果有較大差異，因此，應(yīng)綜合分析三個(gè)數(shù)據(jù)庫彌補(bǔ)單一數(shù)據(jù)庫的不足。

2.4.3 共同上調(diào)表達(dá)差異基因的KEGG通路富集分析

5 151個(gè)差異表達(dá)基因在KEGG數(shù)據(jù)庫中注釋到431個(gè)上調(diào)表達(dá)差異基因，主要分為5大類20個(gè)小類，共涉及109個(gè)代謝通路。其中，細(xì)胞過程、環(huán)境信息處理、遺傳信息加工、代謝和有機(jī)系統(tǒng)5大類注釋到的上調(diào)表達(dá)差異基因數(shù)目分別為39、24、146、203和19。有機(jī)系統(tǒng)類別注釋到14個(gè)參與植物與病原菌互作通路的差異表達(dá)基因。將顯著性設(shè)為qvalue<0.05，對(duì)上調(diào)差異基因進(jìn)行富集分析，發(fā)現(xiàn)差異基因在剪接體、苯丙氨酸代謝、過氧化物酶體、糖和氨基酸的代謝等通路顯著富集 (圖5)。

2.5 R基因和新基因的發(fā)掘

基于所選參考基因組序列的比對(duì)結(jié)果，共注釋到包括7種類型的33個(gè)R(resistance)基因，其中，12個(gè)與富含半胱氨酸類受體蛋白激酶有關(guān)，10個(gè)與抗病性蛋白有關(guān)(表2)。此外，還篩選到β-半乳糖苷酶、纖維素合酶、過氧化物酶等病程相關(guān)蛋白。

表2 中4抗條銹性上調(diào)差異基因中R基因的功能分析Table 2 Functional analysis of R gene in up-regulated differential genes of Zhong 4 resistance to stripe rust

將篩選到的669個(gè)上調(diào)表達(dá)的差異新基因在7種數(shù)據(jù)庫注釋分析，結(jié)果顯示，有335個(gè)獲得GO注釋，123個(gè)獲得COG注釋，79個(gè)獲得KOG注釋，511個(gè)獲得eggNOG注釋；其中，eggNOG數(shù)據(jù)庫最大的兩類注釋分別為功能未知和一般功能預(yù)測(cè)，占48.57%，KOG數(shù)據(jù)庫中細(xì)胞周期控制、細(xì)胞分裂和染色體分裂通路中獲得注釋信息的基因比例較高。151個(gè)上調(diào)表達(dá)差異新基因獲得KEGG注釋，且其中85個(gè)差異新基因共涉及89個(gè)信號(hào)通路；562個(gè)獲得NR注釋。其中，苯丙素的生物合成和苯丙氨酸代謝通路分別富集到7個(gè)上調(diào)表達(dá)差異基因，谷胱甘肽代謝通路分別富集到6個(gè)，內(nèi)質(zhì)網(wǎng)蛋白加工和碳代謝通路分別富集到5個(gè)，植物與病原菌互作、剪接體、糖酵解和糖異生以及淀粉和蔗糖代謝通路分別富集到4個(gè)，說明顯著富集的差異新基因主要與糖的合成與代謝、氨基酸的代謝以及植物與病原菌互作有關(guān)。上調(diào)表達(dá)差異新基因還注釋到R基因、WRKY 和NAC轉(zhuǎn)錄因子等應(yīng)答非生物脅迫的基因。

2.6 qRT-PCR對(duì)RNA-seq數(shù)據(jù)的驗(yàn)證

為驗(yàn)證轉(zhuǎn)錄組測(cè)序RNA-seq數(shù)據(jù)的準(zhǔn)確度，在差異基因中隨機(jī)選取與NB-ARC、鉀離子轉(zhuǎn)運(yùn)蛋白和熱激蛋白相關(guān)的3個(gè)基因進(jìn)行qRT-PCR，其中，轉(zhuǎn)錄本數(shù)據(jù)經(jīng)NCBI Blast確定其特異序列位置，用Premier 6設(shè)計(jì)與內(nèi)參引物TaEF-1α 退火溫度值相近的特異引物。以0 h為對(duì)照，計(jì)算基因相對(duì)表達(dá)水平與轉(zhuǎn)錄組數(shù)據(jù)表達(dá)量并進(jìn)行對(duì)比。結(jié)果(表3)表明，轉(zhuǎn)錄組測(cè)序RNA-seq數(shù)據(jù)與qRT-PCR數(shù)據(jù)整體趨勢(shì)一致，轉(zhuǎn)錄組測(cè)序數(shù)據(jù)準(zhǔn)確可信。

表3 熒光定量PCR結(jié)果與轉(zhuǎn)錄組測(cè)序結(jié)果對(duì)比Table 3 Comparison of real-time PCR results with RNA-seq results

3 討 論

條銹病的流行嚴(yán)重威脅小麥的產(chǎn)量和品質(zhì)，培育抗病品種是防治該病最有效、經(jīng)濟(jì)和環(huán)保的措施[1-4]。中4自育成至今，對(duì)我國(guó)現(xiàn)有條銹菌小種均表現(xiàn)為抗病，然而其抗條銹病機(jī)制尚不明確。

本研究利用轉(zhuǎn)錄組測(cè)序，分析了條銹菌CYR32侵染中412、24和48 h后的差異表達(dá)基因，獲得 5 151個(gè)持續(xù)上調(diào)表達(dá)差異基因。經(jīng)多個(gè)數(shù)據(jù)庫注釋及富集分析，篩選得到R基因、PR基因、蛋白酶體、WRKY和NAC轉(zhuǎn)錄因子等與調(diào)控小麥條銹菌有關(guān)的基因。其中33個(gè)上調(diào)表達(dá)差異基因參與植物與病原菌互作過程；大量上調(diào)表達(dá)差異新基因獲得注釋，在植物與病原菌互作過程、糖的合成與代謝以及氨基酸的代謝通路顯著富集。

Chen等[13]通過轉(zhuǎn)錄組學(xué)分析小種?；钥剐?、全生育期抗性和高溫成株期抗性，發(fā)現(xiàn)高溫成株期抗性中基因多樣性最高，但是基因誘導(dǎo)表達(dá)緩慢，解釋了其抗性持久性的原因。Coram等[14]分析Yr39介導(dǎo)的與條銹菌親和與不親和的兩個(gè)F7重組自交系的轉(zhuǎn)錄本，發(fā)現(xiàn)條銹菌侵染48 h后時(shí)，轉(zhuǎn)錄本數(shù)量均達(dá)到峰值；對(duì)轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行注釋后發(fā)現(xiàn)，50%以上的轉(zhuǎn)錄本參與防御和信號(hào)轉(zhuǎn)導(dǎo)過程，例如R基因同源物、PR蛋白、苯丙氨酸的生物合成和蛋白激酶等。Coram等[15]分析Yr5介導(dǎo)的與條銹菌親和與不親和的相關(guān)轉(zhuǎn)錄本，發(fā)現(xiàn)在侵染24 h時(shí)轉(zhuǎn)錄本積累量出現(xiàn)峰值，注釋后發(fā)現(xiàn)被誘導(dǎo)的轉(zhuǎn)錄本均參與基礎(chǔ)防御相關(guān)過程，Yr5快速響應(yīng)的信號(hào)通路和防御反應(yīng)有R基因參與，包括蛋白激酶信號(hào)轉(zhuǎn)導(dǎo)和活性氧導(dǎo)致的過敏性壞死反應(yīng)。本研究結(jié)果發(fā)現(xiàn)，條銹菌侵染24 h時(shí)，差異表達(dá)基因數(shù)目出現(xiàn)峰值，表明小麥在應(yīng)對(duì)條銹菌的誘導(dǎo)下，大量與寄主抗性及相關(guān)抗性基因轉(zhuǎn)錄表達(dá)來調(diào)控免疫反應(yīng)。Hao等[16]研究發(fā)現(xiàn)，在條銹菌侵染24 h時(shí)，KEGG顯著富集的通路有苯丙氨酸的生物合成和代謝，該途徑可以產(chǎn)生水楊酸、類黃酮和木質(zhì)素等化合物，并且苯丙氨酸代謝途徑在條銹菌侵染初期參與植物防御反應(yīng)；在條銹菌侵染120 h時(shí)，KEGG富集到與淀粉、蔗糖和碳水化合物代謝相關(guān)的基因，推測(cè)小麥光合作用的強(qiáng)度影響其對(duì)條銹菌的抗性。本研究KEGG通路分析發(fā)現(xiàn)，上調(diào)表達(dá)差異基因顯著富集在蛋白質(zhì)加工、苯丙氨酸代謝、過氧化物酶體、糖和氨基酸代謝等通路，與前人研究結(jié)果一致。Garnica等[17]研究發(fā)現(xiàn)，ABC轉(zhuǎn)運(yùn)蛋白基因在孢子萌發(fā)中差異表達(dá)，并在附著期高度富集；MFS轉(zhuǎn)運(yùn)蛋白優(yōu)先在孢子萌發(fā)期表達(dá)，主要影響其生長(zhǎng)過程；己糖轉(zhuǎn)運(yùn)蛋白和氨基酸轉(zhuǎn)運(yùn)體在吸器中高度表達(dá)，條銹菌的生活方式對(duì)轉(zhuǎn)運(yùn)體類型起決定作用。葡萄糖轉(zhuǎn)運(yùn)體在發(fā)芽夏孢子中高表達(dá)，說明條銹菌侵染小麥后消耗寄主碳水化合物吸收營(yíng)養(yǎng)。

PR蛋白在系統(tǒng)獲得抗性和植物防御中尤為重要。研究表明，PR1蛋白能夠抑制小麥條銹菌的生長(zhǎng)[18]。另外，TaWRKY49負(fù)調(diào)控小偃6號(hào)對(duì)小麥條銹菌的抗性，并誘導(dǎo)水楊酸、茉莉酸和活性氧響應(yīng)相關(guān)基因的表達(dá)，同時(shí)抑制乙烯相關(guān)基因的表達(dá)；TaWRKY45和aWRKY62正調(diào)控小偃6號(hào)對(duì)小麥條銹菌的抗性，且水楊酸、乙烯、茉莉酸甲酯、過氧化氫和脫落酸均能夠誘導(dǎo)TaWRKY62基因表達(dá)[19-20]。NAC基因家族編碼的蛋白質(zhì)是目前發(fā)現(xiàn)的最大的植物特異性轉(zhuǎn)錄因子之一[21-23]，在植物應(yīng)答非生物和生物脅迫以及調(diào)控植物發(fā)育等方面發(fā)揮著重要的作用，TaNAC4[21]、TaNAC21/22[22]和TaNAC30[23]負(fù)調(diào)控小麥對(duì)條銹菌的抗性，且茉莉酸甲酯、乙烯和脫落酸能夠誘導(dǎo)TaNAC4基因的表達(dá)，而水楊酸對(duì)TaNAC4基因的表達(dá)無影響[21]；TaNAC21/22作為tae-miR164的靶標(biāo)在小麥對(duì)條銹菌的抗性方面發(fā)揮著重要作用[22]；沉默TaNAC30可顯著增加H2O2的積累量，說明TaNAC30可增加小麥對(duì)條銹菌的抗性[23]。

本研究分析了條銹菌CYR32侵染中4的轉(zhuǎn)錄組數(shù)據(jù)，通過基因組序列比對(duì)獲得7種類型R基因并得到大量的新基因，新基因中也預(yù)測(cè)到R基因和抗條銹性相關(guān)基因，表明小麥品種中4中蘊(yùn)含豐富的抗性基因。