付家旭,閆亞麗,謝小文,張心悅,溫鵬飛,關(guān)小康,王同朝,衛(wèi) 麗
(河南農(nóng)業(yè)大學(xué) 農(nóng)學(xué)院,河南 鄭州 450002)
玉米(Zea mays)是我國(guó)主要的糧食作物,也是重要的飼料和工業(yè)原料。近年來,隨著全球氣候變暖,干旱發(fā)生的頻率和強(qiáng)度逐年增加,對(duì)玉米生產(chǎn)造成嚴(yán)重影響[1]。尤其在播種-出苗階段,若遭遇干旱,會(huì)嚴(yán)重影響玉米出苗和生長(zhǎng)。干旱脅迫條件下,水分虧缺會(huì)引起氣孔關(guān)閉,導(dǎo)致光合作用下降,造成減產(chǎn),嚴(yán)重時(shí)甚至絕收[2]。
轉(zhuǎn)錄因子是一種具有特殊結(jié)構(gòu)、行使調(diào)控基因表達(dá)功能的蛋白質(zhì),通過與靶基因啟動(dòng)子區(qū)的順式作用元件結(jié)合,激活或抑制基因的表達(dá),在植物逆境脅迫中發(fā)揮著重要作用[3]。目前,已發(fā)現(xiàn)WRKY、bZIP(Basic leucine zipper)、NAC(NAM、ATAF1/ATAF2 和CUC2)、MYB(v-myb avian myeloblastosis viral oncogene homolog )和bHLH(Basic helix-loophelix)等轉(zhuǎn)錄因子家族與逆境脅迫密切相關(guān)[4-8]。WRKY 是植物中較龐大的一個(gè)轉(zhuǎn)錄因子家族,自ISHIGURO 等[9]從甘薯中克隆出第一個(gè)WRKY 基因SPF1(Sweet potato factor 1)后,在 擬 南 芥(Arabidopsis thaliana)[10]、水稻(Oryza sativaL.)[11]、高粱(Sorghum bicolor)[12]、小 麥(Triticum aestivumL.)[13]、大麥(Hordeum vulgareL.)[14]、玉米[15]和陸地棉(Gossypium hirsutum)[16]等多個(gè)物種中鑒定到WRKY基因。WRKY 轉(zhuǎn)錄因子通過特異性結(jié)合靶基因啟動(dòng)子的順式作用元件W-box[(T)TGAC(C/T)]而調(diào)控基因表達(dá),在植物生長(zhǎng)發(fā)育及逆境脅迫中具有重要作用[17]。過表達(dá)小麥TaWRKY2基因提高了干旱脅迫條件下轉(zhuǎn)基因小麥葉片中的游離脯氨酸、可溶性糖和葉綠素含量,進(jìn)而提高了對(duì)干旱的耐受性[18]。干旱和鹽脅迫誘導(dǎo)大豆(Glycine max)GmWRKY12表達(dá),過表達(dá)GmWRKY12基因提高了大豆葉片中脯氨酸含量,降低了丙二醛(MDA)含量,最終提高了抗旱性和耐鹽性[19]。過表達(dá)陸地棉GhWRKY33基因可以提高擬南芥對(duì)脫落酸(ABA)和干旱的耐受能力[20]。前期,河南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院旱作節(jié)水課題組以豫882為試驗(yàn)材料,在苗期進(jìn)行干旱-復(fù)水處理轉(zhuǎn)錄組測(cè)序分析,從中得到11 002 個(gè)差異表達(dá)的基因,包括56 個(gè)轉(zhuǎn)錄因子家族,2 332 個(gè)轉(zhuǎn)錄因子[5,21-23]。在前期[23]研究基礎(chǔ)上,對(duì)干旱-復(fù)水處理下差異表達(dá)的玉米WRKY 基因進(jìn)行鑒定,并對(duì)其蛋白質(zhì)理化性質(zhì)、系統(tǒng)進(jìn)化、染色體分布和復(fù)制關(guān)系、基因結(jié)構(gòu)和蛋白質(zhì)保守基序、啟動(dòng)子區(qū)順式作用元件及干旱-復(fù)水條件下的表達(dá)量等進(jìn)行分析,為WRKY基因的進(jìn)一步利用提供理論依據(jù)。
供試玉米材料為豫882,選取籽粒飽滿、均勻的種子,滅菌后播種在土壤∶蛭石為3∶1 的混合土樣中,放置于光照培養(yǎng)箱中[溫度28 ℃/22 ℃,光照周期14 h/10 h,光照強(qiáng)度120 μmol/(m2·s),相對(duì)濕度60%]。待玉米生長(zhǎng)至有3 片展開葉時(shí),將其移栽到霍格蘭(Hoagland)營(yíng)養(yǎng)液中,營(yíng)養(yǎng)液每2 d 更換一次。處理組在營(yíng)養(yǎng)液中加入20%聚乙二醇(PEG)6000 模擬干旱,分別處理60 h 和96 h,復(fù)水3 d,記為T60、T96、TR3d;不進(jìn)行干旱處理的對(duì)照組(霍格蘭營(yíng)養(yǎng)液)分別記為CK60、CK96、CK3d。取葉片立即放入液氮中速凍,然后放置于-80 ℃超低溫冰箱中保存。
使 用Pfam 數(shù) 據(jù) 庫(kù)(http://pfam.xfam.org/)和SMART(http://smart.embl-heidelberg.de/)在線網(wǎng)站進(jìn)一步分析玉米WRKY 基因相對(duì)應(yīng)的蛋白質(zhì),去除冗余和重復(fù)的基因并刪除缺少完整結(jié)構(gòu)域的蛋白質(zhì)序列,并將其結(jié)構(gòu)域與WRKY 結(jié)構(gòu)域進(jìn)行比對(duì),保留含有WRKY 結(jié)構(gòu)域的基因。采用RPKM 值計(jì)算基因在對(duì)照組與處理組之間的表達(dá)量差異倍數(shù),以|log2(差異倍數(shù))|≥1、錯(cuò)誤率(FDR)<0.05 為標(biāo)準(zhǔn),篩選玉米干旱-復(fù)水處理中差異表達(dá)的WRKY 基因,并 參 考Ensembl Plants(http://plants. ensembl. org/index.html)和MaizeGDB(https://www.maizegdb.org/)數(shù)據(jù)庫(kù)進(jìn)行命名。
1.3.1 蛋白質(zhì)分子質(zhì)量和等電點(diǎn) 利用在線網(wǎng)站ExPASy(https://web.expasy.org/protparam/)預(yù)測(cè)和分析玉米WRKY 基因編碼蛋白質(zhì)的分子質(zhì)量和等電點(diǎn)。
1.3.2 系統(tǒng)進(jìn)化 從玉米數(shù)據(jù)庫(kù)MaizeGDB(https://www.maizegdb.org/)中下載51 個(gè)WRKY 蛋白序列,從TAIR數(shù)據(jù)庫(kù)(https://www.arabidopsis.org/)中下載擬南芥的WRKY蛋白序列,采用DNAMAN軟件對(duì)玉米和擬南芥WRKY 蛋白的氨基酸序列進(jìn)行比對(duì),篩選出同源序列。采用MEGA 7.0軟件NJ(Neighbor-joining)法構(gòu)建系統(tǒng)進(jìn)化樹,步長(zhǎng)設(shè)置為1 000。
1.3.3 基因在染色體上的分布和復(fù)制關(guān)系 利用Ensembl Plants 下載玉米基因組注釋信息,使用TBtools分析基因在染色體上的分布和復(fù)制關(guān)系。
1.3.4 基因結(jié)構(gòu)和蛋白質(zhì)保守基序(Motif) 利用GSDS(https://gsds.cbi.pku.edu.cn/)在線分析WRKY 基因外顯子-內(nèi)含子結(jié)構(gòu)。利用MEME Suite 5.4.1(https://meme-suite.org/meme/)在線預(yù)測(cè)玉米WRKY蛋白的保守基序,Motif 個(gè)數(shù)設(shè)置為10 個(gè),基序長(zhǎng)度為6~70 aa,其他參數(shù)均為默認(rèn)值。
1.3.5 啟動(dòng)子區(qū)順式作用元件 從數(shù)據(jù)庫(kù)Ensembl Plants(http://plants.ensembl.org/index.html)獲取WRKY基 因 上 游2 000 bp 序 列,用Plant CARE(http://bioinformatics.psb.ugent.be/webtools/plantcare/html/)對(duì)WRKY基因啟動(dòng)子區(qū)的順式作用元件進(jìn)行分析。
根據(jù)Ensembl Plants 中檢索到的候選WRKY 基因的cDNA 序列,使用Primer Premier 5.0 進(jìn)行qRTPCR 引物設(shè)計(jì)(表1)。利用Trizol 試劑提取葉片總RNA,使用PrimeScriptTMRT reagent Kit with gDNA Eraser(TaKaRa)反轉(zhuǎn)錄試劑盒合成cDNA,按照TB GreenPremix ExTaqTMⅡ(TaKaRa)試劑盒說明書進(jìn)行qRT-PCR,采 用 玉 米18S RNA(GenBank No.AF168884.1)作為內(nèi)參基因。qRT-PCR 反應(yīng)體系為20 μL:cDNA 1 μL,TB Green Premix ExTaqTMⅡ10 μL,上、下游引物(10 μmol/L)各0.8 μL,ddH2O 7.4 μL。qRT-PCR 反應(yīng)程序:95 ℃30 s;95 ℃5 s、60 ℃30 s,40 個(gè)循環(huán)。試驗(yàn)進(jìn)行3 次生物學(xué)重復(fù),采用2-ΔΔCt方法計(jì)算基因的相對(duì)表達(dá)量[24]。
表1 玉米WRKY基因qRT-PCR引物序列Tab.1 qRT-PCR primer sequences for WRKY genes in maize
本研究在前期研究結(jié)果[23]基礎(chǔ)上鑒定出51 個(gè)干旱-復(fù)水處理下差異表達(dá)的玉米WRKY 基因(表2),開放閱讀框(ORF)為690~2 248 bp;編碼的氨基酸數(shù)目最多為729 個(gè),最少為99 個(gè);編碼蛋白質(zhì)分子質(zhì)量為11.22~78.73 ku。編碼蛋白質(zhì)等電點(diǎn)最小為4.58,最大為12.26,平均為7.6。其中,有20個(gè)小于7,占總數(shù)的39.22%;有31 個(gè)大于7,占總數(shù)的60.78%,說明這51 個(gè)差異表達(dá)的玉米WRKY基因編碼的蛋白質(zhì)偏向于堿性。
表2 干旱-復(fù)水處理下差異表達(dá)玉米WRKY基因基本信息Tab.2 The basic information of differentially expressed WRKY genes in maize under drought-rewatering treatment
續(xù)表2 干旱-復(fù)水處理下差異表達(dá)玉米WRKY基因基本信息Tab.2(Continued) The basic information of differentially expressed WRKY genes in maize under drought-rewatering treatment
本研究篩選出31個(gè)干旱-復(fù)水處理下差異表達(dá)玉米WRKY 基因的擬南芥同源基因,構(gòu)建玉米、擬南芥WRKY 蛋白系統(tǒng)進(jìn)化樹(圖1)。82 個(gè)WRKY蛋白可劃分為3 組,其中,Ⅰ組含有9 個(gè)玉米WRKY蛋白,占玉米總數(shù)的17.65%;Ⅱ組有28 個(gè)玉米WRKY 蛋白,占玉米總數(shù)的54.90%,占比最高,說明Ⅱ組可能在玉米進(jìn)化中占主體地位;Ⅲ組有14個(gè)玉米WRKY蛋白,占玉米總數(shù)的27.45%。Ⅱ組又可分為5 個(gè)亞組(Ⅱa—Ⅱe),其中,Ⅱc 亞組中含有12 個(gè)玉米WRKY蛋白,占玉米總數(shù)的27.45%。每個(gè)組都有擬南芥WRKY 蛋白,說明玉米和擬南芥的WRKY基因存在平行進(jìn)化。
圖1 干旱-復(fù)水處理下差異表達(dá)玉米WRKY蛋白與擬南芥同源蛋白系統(tǒng)進(jìn)化分析Fig.1 Phylogenetic evolution analysis between differentially expressed WRKY proteins in maize under drought-rewatering treatment and part of homologous WRKY proteins in Arabidopsis thaliana
51個(gè)干旱-復(fù)水處理下差異表達(dá)的玉米WRKY基因不均勻地分布在10 條染色體上(圖2)。其中,第3 和第8 染色體上分布的WRKY 基因較多,分別有10、11個(gè);分布最少的是第2和第10染色體,各分布2 個(gè)基因。 發(fā)生了2 對(duì)串聯(lián)復(fù)制事件(ZmWRKY26/ZmWRKY80、ZmWRKY27/ZmWRKY56),分別發(fā)生在第8 染色體、第3 染色體上。發(fā)生了16對(duì)片段復(fù)制事件(ZmWRKY102/ZmWRKY109、ZmWRKY73/ZmWRKY82、ZmWRKY106/ZmWRKY86、ZmWRKY87/ZmWRKY52、ZmWRKY43/ZmWRKY57、ZmWRKY27/ZmWRKY80、ZmWRKY74/ZmWRKY86、ZmWRKY48/ZmWRKY124、ZmWRKY25/ZmWRKY38、ZmWRKY29/ZmWRKY111、ZmWRKY43/ZmWRKY115、ZmWRKY56/ZmWRKY26、ZmWRKY103/ZmWRKY111、ZmWRKY82/ZmWRKY68、ZmWRKY57/ZmWRKY115、ZmWRKY74/ZmWRKY106)。其中,7 對(duì)片段復(fù)制事件(ZmWRKY27/ZmWRKY80、ZmWRKY74/ZmWRKY86、ZmWRKY29/ZmWRKY111、ZmWRKY43/ZmWRKY115、ZmWRKY56/ZmWRKY26、ZmWRKY103/ZmWRKY111、ZmWRKY74/ZmWRKY106))發(fā)生在第3 染色體和第8染色體上,表明第3染色體和第8染色體可能與玉米WRKY家族擴(kuò)張密切相關(guān)。
圖2 干旱-復(fù)水處理下差異表達(dá)玉米WRKY基因染色體分布及復(fù)制關(guān)系Fig.2 Chromosome distribution and duplication of differentially expressed WRKY genes in maize under droughtrewatering treatment
干旱-復(fù)水處理下51 個(gè)差異表達(dá)玉米WRKY基因結(jié)構(gòu)如圖3 所示,玉米WRKY 基因的外顯子數(shù)是1~12 個(gè),含有3 個(gè)外顯子的占62.74%。除ZmWRKY25基因外,其余WRKY 基因均含有內(nèi)含子。大部分同一組或亞組的WRKY 基因具有相對(duì)保守的外顯子數(shù),如Ⅱc 中除了ZmWRKY30和ZmWRKY108基因外,均含有3 個(gè)外顯子,Ⅱd 亦是如此。而Ⅰ中外顯子數(shù)量變化較大,介于2~8個(gè),其中 有 3 個(gè) 基 因(ZmWRKY115、ZmWRKY43和ZmWRKY1)含有5 個(gè)外顯子,2 個(gè)基因(ZmWRKY83和ZmWRKY62)含 有3 個(gè) 外 顯 子,ZmWRKY92、ZmWRKY87、ZmWRKY57、ZmWRKY52分別含有2、4、6、8 個(gè)外顯子;Ⅲ中ZmWRKY25只含有1 個(gè)外顯子,ZmWRKY65含有12個(gè)外顯子。
圖3 干旱-復(fù)水處理下差異表達(dá)玉米WRKY基因結(jié)構(gòu)和蛋白質(zhì)保守基序分析Fig.3 Analysis of differentially expressed WRKY gene structure and encoded protein conserved motif in maize under drought-rewatering treatment
41個(gè)玉米WRKY 蛋白含有2~4個(gè)保守基序,占80.39%。同一組中的WRKY 成員具有相似的基序組成,例如,Ⅰ中都含有Motif1、Motif2、Motif4 和Motif9;Ⅱc 中 除ZmWRKY100(含 有Motif1 與Motif4)外 其 余 成 員 都 含 有Motif1、Motif2 和Motif4。
為進(jìn)一步預(yù)測(cè)玉米WRKY 基因的功能,對(duì)51個(gè)玉米WRKY 基因上游2 000 bp 啟動(dòng)子區(qū)域進(jìn)行順式作用元件分析。由圖4 可知,51 個(gè)玉米WRKY 基因啟動(dòng)子區(qū)域含有ABA 響應(yīng)元件(ABRE)、生長(zhǎng)素響應(yīng)元件(AuxRR-core、TGA-element)、水楊酸響應(yīng)元件(TCA-element)、赤霉素響應(yīng)元件(GAREmotif、P-box、TATC-box)、茉 莉 酸 甲 酯 響 應(yīng) 元 件(CGTCA-motif、TGACG-motif)、防衛(wèi)和脅迫響應(yīng)元件(TC-rich repeats)、低溫響應(yīng)元件(LTR)、干旱響應(yīng)元件(MBS)、缺氧誘導(dǎo)有關(guān)元件(GC-motif)等植物激素及非生物脅迫相關(guān)的順式作用元件。其中,AREB、LTR、MBS、TCA-element 和TGACG-motif 分布較多,而GARE-motif、P-box、TC-rich 和TATCbox 分 布 較 少。 除 了ZmWRKY1、ZmWRKY52、ZmWRKY96和ZmWRKY106外,其余基因均含有ABRE。 含 有 CGTCA-motif 較 多 的 基 因 有ZmWRKY15、ZmWRKY29、ZmWRKY40、ZmWRKY53、ZmWRKY92和ZmWRKY104;ZmWRKY28、ZmWRKY52、ZmWRKY108和ZmWRKY115等24個(gè)基因含有MBS;88.24%的ZmWRKY基因含有TGACGmotif,其 中ZmWRKY38、ZmWRKY43、ZmWRKY53、ZmWRKY92、ZmWRKY96和ZmWRKY99含 有4 個(gè) 及以上。以上結(jié)果說明玉米WRKY 基因在激素信號(hào)轉(zhuǎn)導(dǎo)和非生物脅迫響應(yīng)中可能起著重要的作用。
圖4 干旱-復(fù)水處理下差異表達(dá)玉米WRKY基因啟動(dòng)子區(qū)順式作用元件分析Fig.4 Analysis of cis-acting elements in the promoter region of differentially expressed WRKY genes in maize under drought-rewatering treatment
依據(jù)轉(zhuǎn)錄組結(jié)果,在干旱-復(fù)水處理下,51個(gè)玉米WRKY基因呈現(xiàn)不同的表達(dá)模式(圖5)。在干旱處理60 h、96 h 和復(fù)水3 d 3 個(gè)處理時(shí)間點(diǎn),上調(diào)表達(dá)的基因數(shù)分別為25、19、3 個(gè),下調(diào)表達(dá)的基因數(shù)分 別 為26、32、48 個(gè)。ZmWRKY1、ZmWRKY10、ZmWRKY16、ZmWRKY28、ZmWRKY30、ZmWRKY33、ZmWRKY42、ZmWRKY65、ZmWRKY68、ZmWRKY78、ZmWRKY96、ZmWRKY99、ZmWRKY100、ZmWRKY102和ZmWRKY111在干旱處理下,其表達(dá)量高于對(duì)照組,在復(fù)水處理后,其表達(dá)量低于對(duì)照組,說明這15個(gè)基因正向響應(yīng)干旱脅迫。ZmWRKY87基因在干旱處理下下調(diào)表達(dá),復(fù)水處理后,表達(dá)量高于對(duì)照組,是負(fù)向響應(yīng)干旱脅迫基因。其他基因的表達(dá)模式無(wú)規(guī)律,如ZmWRKY15、ZmWRKY52、ZmWRKY103等。
圖5 干旱-復(fù)水處理下差異表達(dá)玉米WRKY基因的表達(dá)分析Fig.5 Expression analysis of differentially expressed WRKY genes in maize under drought-rewatering treatment
為驗(yàn)證轉(zhuǎn)錄組數(shù)據(jù)的可靠性,隨機(jī)挑選了8 個(gè)WRKY 基因(ZmWRKY1、ZmWRKY16、ZmWRKY30、ZmWRKY34、ZmWRKY48、ZmWRKY68、ZmWRKY82、ZmWRKY102)進(jìn)行qRT-PCR 分析。由圖6 可知,ZmWRKY1、ZmWRKY16、ZmWRKY30、ZmWRKY68、ZmWRKY102在干旱處理60 h 和96 h 條件下表達(dá)量上升,復(fù)水3 d 后表達(dá)量下降,正向響應(yīng)干旱脅迫。其中,ZmWRKY1、ZmWRKY30、ZmWRKY68在3 個(gè)處理時(shí)間點(diǎn)下均達(dá)到顯著差異水平。ZmWRKY1、ZmWRKY16、ZmWRKY30、ZmWRKY34、ZmWRKY48、ZmWRKY68、ZmWRKY82、ZmWRKY102這8 個(gè)基因在3 個(gè)處理時(shí)間點(diǎn)的qRT-PCR 結(jié)果和轉(zhuǎn)錄組結(jié)果趨勢(shì)基本一致,進(jìn)一步驗(yàn)證了轉(zhuǎn)錄組結(jié)果的準(zhǔn)確性。
圖6 干旱-復(fù)水處理下部分差異表達(dá)玉米WRKY基因的表達(dá)驗(yàn)證Fig.6 Expression validation of partial differentially expressed WRKY genes in maize under drought-rewatering treatment
玉米是高需水、需肥作物,而我國(guó)玉米主產(chǎn)區(qū)大多都存在水資源不足問題,尤其是近年來全球氣候變暖,干旱導(dǎo)致玉米產(chǎn)量降低的問題日益突出。玉米生長(zhǎng)發(fā)育進(jìn)程中,經(jīng)常處于干旱-濕潤(rùn)-干旱的動(dòng)態(tài)環(huán)境中,尤其是河南省夏季,雨量分配不均。夏玉米播種期或苗期,經(jīng)常遭遇干旱,而生育中期或后期,降雨量偏多。適宜播期遇旱導(dǎo)致出苗較差,缺苗斷壟,苗期干旱導(dǎo)致發(fā)育受阻。因此,發(fā)掘抗旱基因,培育抗旱品種,提高玉米自身的抗旱能力對(duì)玉米生產(chǎn)的發(fā)展尤為重要。本研究在轉(zhuǎn)錄組結(jié)果分析的基礎(chǔ)上,鑒定篩選出響應(yīng)干旱-復(fù)水處理的51 個(gè)差異表達(dá)的WRKY 基因。通過生物信息學(xué)分析發(fā)現(xiàn),51 個(gè)WRKY 基因不均勻地分布在10條染色體上,部分基因存在片段復(fù)制關(guān)系,暗示著這些基因在玉米WRKY 基因家族進(jìn)化中發(fā)揮重要作用[25]。進(jìn)化樹分析結(jié)果顯示,51 個(gè)WRKY 基因分成了3 個(gè)組,與小麥、擬南芥[26]分類結(jié)果一致,進(jìn)一步說明WRKY 轉(zhuǎn)錄因子家族具有平行進(jìn)化關(guān)系。51 個(gè)玉米WRKY 基因的啟動(dòng)子區(qū)域含有多個(gè)與植物激素和逆境脅迫相關(guān)的順式作用元件,如ABRE、TGACG-motif、TGA-element、LTR、MBS 等,進(jìn)而更加證明WRKY 轉(zhuǎn)錄因子在植物逆境脅迫響應(yīng)中起重要的調(diào)控作用。
逆境脅迫下,植物自身通過激活或抑制基因的表達(dá)來適應(yīng)環(huán)境條件的變化,進(jìn)而達(dá)到一種平衡。HU 等[27]研究發(fā)現(xiàn),在超表達(dá)陸地棉GhWRKY1-like基因的擬南芥植株中,GhWRKY1-like 與ABA 生物合 成 所 必 需 的AtNCED2、AtNCED5、AtNCED6和AtNCED9基因啟動(dòng)子的W-box 順式作用元件結(jié)合,促進(jìn)這些靶基因的表達(dá),從而提高植株的抗旱性。XIONG 等[3]研究發(fā)現(xiàn),超表達(dá)文冠果XsWRKY20基因通過ABA 途徑提高煙草的抗旱性。本研究發(fā)現(xiàn),在干旱-復(fù)水處理下,51 個(gè)玉米WRKY 基因呈現(xiàn)出不 同 的 表 達(dá) 模 式。ZmWRKY1、ZmWRKY10、ZmWRKY16、ZmWRKY28、ZmWRKY30、ZmWRKY33、ZmWRKY42、ZmWRKY65、ZmWRKY68、ZmWRKY78、ZmWRKY96、ZmWRKY99、ZmWRKY100、ZmWRKY102和ZmWRKY111這15 個(gè)正向響應(yīng)干旱脅迫基因與1個(gè)負(fù)向響應(yīng)干旱脅迫基因(ZmWRKY87)是今后深入研究ZmWRKY家族干旱脅迫響應(yīng)的候選基因,可為玉米抗旱分子育種以及培育玉米新品種奠定基礎(chǔ)。