劉紅梅，胡尚連*，曹 穎，盧學(xué)琴，徐 剛，龍治堅(jiān)，任 鵬

(1.西南科技大學(xué)生命科學(xué)與工程學(xué)院，綿陽(yáng) 621010；2.四川省生物質(zhì)資源利用與改性工程技術(shù)研究中心，綿陽(yáng) 621010)

慈竹WRKY轉(zhuǎn)錄因子基因克隆及其脅迫誘導(dǎo)表達(dá)

劉紅梅1,2，胡尚連1,2*，曹 穎1,2，盧學(xué)琴1,2，徐 剛1,2，龍治堅(jiān)1,2，任 鵬1,2

(1.西南科技大學(xué)生命科學(xué)與工程學(xué)院，綿陽(yáng) 621010；2.四川省生物質(zhì)資源利用與改性工程技術(shù)研究中心，綿陽(yáng) 621010)

為了解慈竹WRKY生物學(xué)功能以及通過(guò)基因工程手段改良竹類植物。利用已獲得慈竹轉(zhuǎn)錄組數(shù)據(jù)，克隆獲得2條WRKY轉(zhuǎn)錄因子,分別命名為BeWRKY1 (GenBank：KJ462124)和BeWRKY2 (GenBank：KJ462125)，兩者分別編碼189和295個(gè)氨基酸殘基。序列分析表明，兩者編碼的蛋白均含有一個(gè)完整的WRKY superfamily的保守結(jié)構(gòu)域。BeWRKY1屬于Ⅱc亞家族；BeWRKY2則與TaWRKY16和AtWRKY39等聚為一枝，屬于Ⅱd亞家族。定量PCR分析表明，BeWRKY1和BeWRKY2分別在莖和完全展開葉中具有較高表達(dá)量。200mmol·L-1NaCl和20% PEG 6000脅迫處理激活了BeWRKY1表達(dá)，卻減少了BeWRKY2的轉(zhuǎn)錄積累。結(jié)果表明BeWRKY1可能短時(shí)間內(nèi)迅速參與逆境響應(yīng)；BeWRKY2對(duì)干旱和ABA脅迫傷害更為敏感，而可能不參與高鹽類脅迫保護(hù)反應(yīng)。

慈竹；WRKY；克?。簧镄畔W(xué)分析；脅迫誘導(dǎo)

WRKY蛋白是植物特有的轉(zhuǎn)錄因子家族[1]，含有高度保守的WRKY域的鋅指蛋白，其保守域約由60個(gè)氨基酸殘基組成[2]。根據(jù)WRKY結(jié)構(gòu)域和鋅指結(jié)構(gòu)的數(shù)量將其分為Ⅰ類、Ⅱ類、Ⅲ類[3]。Ⅰ類含有兩個(gè)WRKY結(jié)構(gòu)域，Ⅱ類和Ⅲ類雖都含有一個(gè)WRKY結(jié)構(gòu)域，但鋅指結(jié)構(gòu)不同。其中的Ⅱ類還根據(jù)不同的序列特征分為5個(gè)亞群[4]。目前已從擬南芥[5]、水稻(Oryzasativa)[6]、小麥(Triticumaestivum)[7]、(Solanumlycopersicum)[2]、毛竹(Phyllostachysheterocyclacv.Pubescens)[8]等植物中克隆得到大量的WRKY轉(zhuǎn)錄因子，但在慈竹中尚未見報(bào)道。WRKY蛋白具有快速、瞬時(shí)、具組織特異性誘導(dǎo)表達(dá)特性。其不僅與植物的生物與非生物脅迫應(yīng)答密切相關(guān)，也與植物生長(zhǎng)發(fā)育調(diào)控有關(guān)[9-11]。WRKY在非生物脅迫和非生物方面，研究表明水稻OsWRKY11[12]和OsWRKY45[ 13]、菊花WRKY3[14]、有24條WRKY基因[15]以及擬南芥AtWRKY30[16]均表現(xiàn)由抗干旱脅迫的能力，其中OsWRKY45不但具有抗干旱的能力還有抗高溫和抗病毒的能力[13]。在植物生長(zhǎng)發(fā)育以及衰老過(guò)程中，研究發(fā)現(xiàn)ScWRKY1轉(zhuǎn)錄因子在胚胎形成過(guò)程中發(fā)揮著至關(guān)重要的作用[17]。科學(xué)家發(fā)現(xiàn)在擬南芥的幼片和成熟葉片AtWRKY6轉(zhuǎn)錄因子幾乎不表達(dá)，而在其衰老葉片中檢測(cè)到很高的表達(dá)量，因此判斷其與植物衰老程序密切相關(guān)[18]。雖然目前大量的WRKY轉(zhuǎn)錄因子功能已被證實(shí)，但還有大部分WRKY轉(zhuǎn)錄因子待更深入的進(jìn)行研究。

慈竹(Bambusaemeiensis)是我國(guó)西南地區(qū)普遍栽培的大型叢生竹種[19]，廣泛用于造紙和園林綠化[20]，在所有竹種中慈竹屬于iV級(jí)紙漿竹[21]，是重要的工業(yè)原料。但慈竹易受多種不利因素的影響而大面積減產(chǎn)。既影響新竹，又影響筍的生長(zhǎng)[22]；因此改良慈竹品種迫在眉睫。近年來(lái)隨著基因組和轉(zhuǎn)錄組測(cè)序技術(shù)的廣泛應(yīng)用，毛竹、麻竹等竹種的基因組和(或)轉(zhuǎn)錄組數(shù)據(jù)已經(jīng)被公布[23-24]，但慈竹的研究相對(duì)較少。本課題組利用RNA-Seq技術(shù)進(jìn)行了慈竹筍轉(zhuǎn)錄組研究，為發(fā)現(xiàn)慈竹筍發(fā)育、特征代謝合成和分解途徑等提供了大量的候選基因，還為慈竹SSR和SNP等分子標(biāo)記的開發(fā)提供了極為豐富的序列基礎(chǔ)(數(shù)據(jù)待發(fā)表)。本研究以此轉(zhuǎn)錄組數(shù)據(jù)為基礎(chǔ)，篩選并克隆了其中2個(gè)WRKY轉(zhuǎn)錄因子，并進(jìn)行生物信息學(xué)和脅迫誘導(dǎo)表達(dá)模式分析，為了解慈竹WRKY生物學(xué)功能以及通過(guò)基因工程手段改良竹類植物奠定了一定基礎(chǔ)。

1 材料與方法

1.1 實(shí)驗(yàn)試劑和菌株

OMEG植物總RNA提取試劑盒購(gòu)自成都博瑞克生物技術(shù)有限公司；PrimerScript RT reagent Kit perferct real Time反轉(zhuǎn)錄試劑盒、TaKaRa LA Taq 聚合酶、T4-DNA連接酶、PGM-T vector、EcoRⅠ均購(gòu)至天津的TaKaRa公司；膠回收試劑盒，質(zhì)粒試劑盒，DH5α感受態(tài)細(xì)胞；IPTG、X-Gal、Amp+等購(gòu)至北京的天根公司；Nacl、PEG -6000、ABA均為國(guó)產(chǎn)分析純?cè)噭?/p>

1.2 慈竹WRKY轉(zhuǎn)錄因子的克隆

采集西南科技大學(xué)資源圃中的慈竹筍(露出地面100 cm)用于轉(zhuǎn)錄組高通量測(cè)序。以獲得的慈竹轉(zhuǎn)錄組數(shù)據(jù)為基礎(chǔ)，從100 cm高的慈竹筍克隆了2個(gè)WRKY基因：WRKY1的引物為5′ATGGCGGCTTCGTTAGGACTAA 3′(上游引物)和5′TCAGAAGAGCAGTGAGCCTGCA 3′(下游引物)；WRKY2的引物為5′ATGGAGGAAGTGGAGGAGGCCA 3′(上游引物)和5′CTAAGCTTGCGCAGACTGAGTT 3′(下游引物)。采用 LA-Taq酶(TaKaRa)進(jìn)行目的基因擴(kuò)增。將PCR回收產(chǎn)物連接到pGM-T載體，測(cè)序分析。

1.3 氨基酸多序列比對(duì)和系統(tǒng)發(fā)育樹的構(gòu)建

通過(guò)NCBI里的BLAST進(jìn)行在線分析與之相近物種的同源率。利用ClustalW2將兩者編碼的氨基酸序列與GenBank數(shù)據(jù)庫(kù)中獲得的毛竹、小麥的全長(zhǎng)WRKY編碼的氨基酸序列進(jìn)行比對(duì)分析。利用Mega5.2的Pairwise alignment和Multiple alignment程序?qū)寺～@得的慈竹BeWRKY1和BeWRKY2以及在GeneBank里搜索得到的毛竹、小麥、水稻、擬南芥、大麥(Hordeumvulgare)的151條WRKY進(jìn)行分析，再采用Construct/Test Neibhbor-Joining Tree生成系統(tǒng)發(fā)育進(jìn)化樹。

1.4 WRKY基因編碼蛋白的三級(jí)結(jié)構(gòu)預(yù)測(cè)

利用ExPaSy(www.expasy.ch/tools/)工具中的 CPH models程序在線對(duì)慈竹、毛竹、小麥和水稻進(jìn)行同源模建，再利用RasMol預(yù)測(cè)軟件對(duì)預(yù)測(cè)的WRKY基因編碼蛋白的三維結(jié)構(gòu)進(jìn)行分析。

1.5 慈竹WRKY基因的表達(dá)模式和脅迫誘導(dǎo)表達(dá)分析

用Real-timePCR方法進(jìn)行組織表達(dá)模式和脅迫誘導(dǎo)表達(dá)分析。分別采集慈竹的竹筍(100 cm)，當(dāng)年生且完全抽枝的慈竹莖稈、完全展開葉和未完全展開葉進(jìn)行組織表達(dá)模式分析，每一個(gè)分析取三次生物學(xué)重復(fù)。取培養(yǎng)30 d且生長(zhǎng)狀態(tài)良好的慈竹實(shí)生幼苗(蛭石+1/10倍Hogland營(yíng)養(yǎng)液培養(yǎng))，參考擬南芥[25]、玉米[26]、毛竹[8,27]的脅迫處理方法，分別以NaCl(0、50、100、150、200、250 mmol·L-1)、PEG-6000(0、5、10、15、20、25%)、ABA(0、5、10、50、100、200 μmol·L-1)進(jìn)行預(yù)實(shí)驗(yàn)，每個(gè)處理重復(fù)3次(每個(gè)重復(fù)取3個(gè)獨(dú)立植株)。由預(yù)試驗(yàn)結(jié)果分析，本試驗(yàn)分別以NaCl(0、200 mmol·L-1)、PEG6000(0、20%) 、ABA(0、200 μmol·L-1)進(jìn)行脅迫處理，每個(gè)處理重復(fù)3次(每個(gè)重復(fù)取3個(gè)獨(dú)立植株)，分別在0 h(未處理的對(duì)照)、6 h、12 h、24 h、48 h及7 d(天)取樣，保存于-80 ℃冰箱備用。采用Plant RNA Kit試劑盒提取RNA，用PrimeScript?RT reagent Kit perfect real Time反轉(zhuǎn)錄試劑盒合成cDNA。分別以WRKY1F/R(5′GGCAAGAAGGCTGTCAAG3′;5′ACCCCGTCGTAGGTAGTG3′)和WRKY2F/R(5′TTTCCAGTCCTTCTTGTCG3′;5′CACTTTCCATTCCCATCC3′)為引物，以慈竹Tublin基因[19]為內(nèi)參(Tublin-F 5′GCCGTGAACTCATCCCCTT3′；Tublin-R 5′TTGTTCTTGGCATCCCACAT3′)，在Bio-Rad公司的iQ5 Multicolor上進(jìn)行Real- time PCR分析，用2-ΔΔCT計(jì)算法計(jì)算其表達(dá)量[28]。

2 結(jié)果分析

2.1 慈竹WRKY轉(zhuǎn)錄因子的克隆和氨基酸多序列比對(duì)

根據(jù)設(shè)計(jì)的引物，以100 cm高的慈竹筍cDNA為模板，進(jìn)行PCR擴(kuò)增，獲得了兩個(gè)WRKY轉(zhuǎn)錄因子。測(cè)序分析結(jié)果表明，BeWRKY1(GenBank：KJ462124)長(zhǎng)570 bp，編碼189個(gè)氨基酸殘基；BeWRKY2(GenBank：KJ462125)長(zhǎng)888bp，編碼295個(gè)氨基酸殘基。將克隆到的慈竹BeWRKY1和BeWRKY2，與毛竹的PheWRKY1(GenBank：GU944762.1)、小麥的TaWRKY16(GenBank：EU665428.1)的氨基酸序列比對(duì)。結(jié)果表明，這四個(gè)WRKY轉(zhuǎn)錄因子均含有一個(gè)WRKY保守域結(jié)構(gòu)(圖1，黑色框區(qū)域)。BLAST結(jié)果表明，慈竹BeWRKY1 和毛竹的PheWRKY1 基因的相似性最高，為92%，都屬于酸性蛋白；慈竹BeWRKY2 基因與小麥的TaWRKY16也有較高的相似性，為89%，均屬于堿性蛋白。

圖1 四個(gè)WRKY轉(zhuǎn)錄因子編碼的氨基酸多序列比對(duì)Fig.1 Amino acid sequence alignmen of four WRKY transcription factors

2.2 WRKY轉(zhuǎn)錄因子系統(tǒng)發(fā)育樹的構(gòu)建

利用Mega5.2構(gòu)建進(jìn)化樹，對(duì)慈竹、毛竹、小麥、水稻、擬南芥以及大麥共計(jì)151條WRKY轉(zhuǎn)錄因子進(jìn)行分析。所有的轉(zhuǎn)錄因子聚成三大類，分別為Ⅰ類，Ⅱ類和Ⅲ類。Ⅱ類又聚成Ⅱa，Ⅱb，Ⅱc，Ⅱd及Ⅱe5個(gè)亞類。圖2A可知，BeWRKY1與PheWRKY1、AtWRKY50 和OsWRKY67 聚成一小枝，它們屬于Ⅱc亞類；BeWRKY2與TaWRKY16、AtWRKY39和AtWRKY74聚為一枝，屬于Ⅱd亞類。

用MEME軟件分析這些WRKY氨基酸序列的保守基序(圖2-B)。其中，BeWRKY2和BeWRKY2、PheWRKY1和TaWRKY16均有3個(gè)motif。其中motif 1的寬度為35，其正則表達(dá)式為[ER]GC[GP][AV][KR]K[HR]VER[CD][VGR]DDP[SCR][MY][LV]I[TV]TY[DE]G[EV]H[NK]H[AT][RT][MP][GP][FT][GQ]；motif 2的寬度為41，其正則表達(dá)式為[KR][IV][GP][AF][IR][ST][NR][KS][EI][AV][DE]I[LP][DP]D[EG][FY][KS]WRKYG[KQ]K[AP][IV]K[GN]SP[HN]PR[GN]YY[KR]CS[ST]；motif 3寬度為27，正則表達(dá)式為[AE][LMT][DS]D[HP][VHR][NP][AFS][EQ][QY][AS]P[AR][PR][RV][CF][AST]G[AR]G[DE]D[GH]N[GDE][KN][CDT]。從motif模式圖中看出，BeWRKY1和PheWRKY12、BeWRKY2和TaWRKY2結(jié)構(gòu)相似，與進(jìn)化樹分析結(jié)果一致，進(jìn)步一論證了它們同屬于WRKYⅡ類家族。

圖2 慈竹BeWRKY1和BeWRKY2基因與其他物種WRKY基因編碼氨基酸的聚類(A)和氨基酸序列保守基序分析(B)Fig.2 A：Cluster of the deduced amino acid sequences of BeWRKY1 and BeWRKY2; B：The analysis of their conserved motifs Be：Bambusa emeiensis Os：Oryza sativa At：Arabidopsis thaliana Ta：Triticum aestivum Hv：Hordeum vulgare L. Phe：Phyllostachys heterocycla cv. Pubescens

2.3 慈竹、毛竹、小麥WRKY轉(zhuǎn)錄因子編碼的蛋白三級(jí)結(jié)構(gòu)預(yù)測(cè)

蛋白質(zhì)三級(jí)結(jié)構(gòu)是在二級(jí)結(jié)構(gòu)上折疊，盤繞而成。CPH models在線預(yù)測(cè)慈竹、毛竹和小麥的WRKY轉(zhuǎn)錄因子編碼蛋白的三級(jí)結(jié)構(gòu)。如圖3所示，三個(gè)不同物種均含有豐富的β-折疊?β-轉(zhuǎn)角，無(wú)其它的結(jié)構(gòu)。慈竹BeWRKY2與小麥TaWRKY16編碼的蛋白質(zhì)的β-折疊、β-轉(zhuǎn)角數(shù)量相同，兩者結(jié)構(gòu)相似。慈竹BeWRKY1與PheWRKY1編碼的蛋白質(zhì)有相同數(shù)量的β-轉(zhuǎn)角，但慈竹的β-折疊比毛竹多3個(gè)，兩者結(jié)構(gòu)存在一定的差別。

圖3 慈竹、毛竹和小麥WRKY轉(zhuǎn)錄因子編碼蛋白質(zhì)的三級(jí)結(jié)構(gòu)Fig.3 Model of the tertiary structure of protein coded by WRKY transcription factors from Bambusa emeiensis，Phyllostachys heterocycla cv.Pubescens and Triticum aestivum

圖4 實(shí)時(shí)定量PCR分析BeWRKY1，BeWRKY2在不同組織的表達(dá)Fig.4 Quantitative real-time PCR analysis of the expression patterns of BeWRKY1 and BeWRKY2 in different tissues

2.4 慈竹BeWRKY1、BeWRKY2基因組織表達(dá)分析

分析BeWRKY1和BeWRKY2基因在慈竹中組織特異性表達(dá)。結(jié)果表明，這兩個(gè)轉(zhuǎn)錄因子具有一定的組織表達(dá)特異性(圖4)。BeWRKY1相對(duì)表達(dá)量依次為莖>筍>未展開葉>完全展開葉，各個(gè)部位的相對(duì)表達(dá)量的差異較大，在莖中表達(dá)量最高。而BeWRKY2的相對(duì)表達(dá)量為完全展開葉>莖>未展開葉>筍，展開葉的相對(duì)表達(dá)量高于其他各部位的相對(duì)表達(dá)量。

2.5 脅迫處理對(duì)葉片中BeWRKY1和BeWRKY2表達(dá)影響

用實(shí)時(shí)定量PCR分析了在不同脅迫處理下慈竹幼苗葉片中BeWRKY1和BeWRKY2的表達(dá)情況，并進(jìn)行顯著性分析(Tublin基因?yàn)閮?nèi)參)。如圖5A所示，在NaCl(200 mmol·L-1)和PEG- 6000(20%)脅迫6 h，BeWRKY1表達(dá)水平在6 h時(shí)顯著性上調(diào)，分別是未處理的3.3倍、1.5倍，但隨著脅迫時(shí)間的延長(zhǎng)，它的表達(dá)量緩慢下降，說(shuō)明此基因可能短時(shí)間內(nèi)迅速參與逆境響應(yīng)。在NaCl脅迫處理48 h和7 d時(shí)，BeWRKY2的表達(dá)水平并沒(méi)有明顯變化，可見鹽脅迫對(duì)BeWRKY2轉(zhuǎn)錄積累的影響較?。欢赑EG和ABA脅迫處理的慈竹幼苗葉片中，BeWRKY2的表達(dá)在6 h后顯著下調(diào)：脅迫至12 h時(shí)，BeWRKY2的表達(dá)水平降至最低，隨后逐漸回升，7 d時(shí)與對(duì)照基本持平(圖5B),說(shuō)明BeWRKY2對(duì)干旱和ABA脅迫更為敏感。

圖5 脅迫處理下BeWRKY1和BeWRKY2的表達(dá)量變化Fig.5 Differential expression of BeWRKY1 and BeWRKY2 under stress treatments

3 討論與結(jié)論

WRKY轉(zhuǎn)錄因子廣泛存在于植物中，是近年來(lái)發(fā)現(xiàn)的植物體內(nèi)特有的一類轉(zhuǎn)錄因子，其與植物的非生物及生物脅迫應(yīng)答[12-13]、抗病[8]、生長(zhǎng)發(fā)育[29]、以及衰老息息相關(guān)[15]。近年來(lái)WRKY 轉(zhuǎn)錄因子在植物抗病反應(yīng)中的研究已經(jīng)成為植物分子生物學(xué)領(lǐng)域的前沿課題。

本研究以慈竹轉(zhuǎn)錄組數(shù)據(jù)為基礎(chǔ)，克隆了兩條慈竹WRKY轉(zhuǎn)錄因子序列BeWRKY1和BeWRKY2。其中BeWRKY1與已經(jīng)報(bào)道的毛竹PheWRKY1轉(zhuǎn)錄因子同源率高達(dá)92%，都屬于酸性蛋白；BeWRKY2與小麥TaWRKY16轉(zhuǎn)錄因子同源最高為89%，均屬于堿性蛋白。系統(tǒng)進(jìn)化樹分析表明，BeWRKY1與PheWRKY1、AtWRKY50 和OsWRKY67 聚成一小枝，屬于Ⅱc亞家族類。這一家族中的某些成員，如與PheWRKY1等，在抗病方面以及抗非生物脅迫具有一定的功能[8]。而BeWRKY2與AtWRKY39、AtWRKY74和TaWRKY2聚為一枝，屬于Ⅱd亞家族類。研究表明AtWRKY39突變體種子在高溫處理下萌發(fā)率降低，存活率也大大降低。反之，過(guò)表達(dá)植株對(duì)熱高溫脅迫耐受力增強(qiáng)。在高溫脅迫下擬南芥AtWRKY39突變體植物中水楊酸調(diào)控的PR1和水楊酸相關(guān)的MBF1c基因均表現(xiàn)為下調(diào)，在過(guò)表達(dá)植株中兩基因的表達(dá)增強(qiáng)，表明AtWRKY39正調(diào)控水楊酸和茉莉酮酸通路相互作用應(yīng)答高溫脅迫[30]。推斷這一亞家族中的成員可能與抗高溫干旱有關(guān)。

慈竹兩個(gè)WRKY轉(zhuǎn)錄因子均有WRKY superfamily的保守結(jié)構(gòu)域，但編碼的蛋白理化參數(shù)差異較大。實(shí)時(shí)定量PCR結(jié)果表明，BeWRKY1和BeWRKY2分別在莖和展開葉的表達(dá)量較高，表明這兩個(gè)轉(zhuǎn)錄因子在抗逆方面可能具有不同的功能。逆境誘導(dǎo)表達(dá)的結(jié)果也證實(shí)了這一推論。在200 mmol·L-1NaCl和20% PEG 6 000處理下，BeWRKY1在脅迫處理6 h時(shí)表達(dá)量均明顯升高，說(shuō)明此基因可能短時(shí)間內(nèi)迅速參與逆境響應(yīng)；脅迫后期其表達(dá)量又緩慢下降至與0 h基本持平，可能是植物體避免過(guò)度應(yīng)答造成傷害的一種自我保護(hù)特性[2]。研究報(bào)道多個(gè)WRKY轉(zhuǎn)錄因子參與植物抗旱和抗鹽脅迫[31]。如與BeWRKY1同屬于Ⅱc家族的TaWRKY10(GenBank：HQ700327)基因在小麥?zhǔn)芨啕}、干旱脅迫下均表現(xiàn)為先急劇上調(diào)后緩慢下降與0h持平。在過(guò)表達(dá)TaWRKY10的煙草植株中，TaWRKY10通過(guò)降低ROS積累和激活脅迫防御相關(guān)基因的表達(dá)來(lái)提高抗旱和抗鹽的能力[32]。另外Babitha等證明Ⅱc 家族的AtWRKY28過(guò)表達(dá)的擬南芥植物在高鹽、干旱和氧化處理下抗脅迫能力明顯增強(qiáng)[33]。在干旱和ABA脅迫下，BeWRKY2表達(dá)明顯降低，表明BeWRKY2對(duì)干旱和ABA脅迫傷害更為敏感，而可能不參與高鹽類脅迫保護(hù)反應(yīng)。據(jù)研究報(bào)道屬于Ⅱd家族的水稻OsWRKY11基因增強(qiáng)了水稻耐高溫和干旱的能力[12]，因此BeWRKY2可能與抗干旱等有關(guān)。BeWRKY1和BeWRKY2的抗非生物脅迫功能仍需進(jìn)一步驗(yàn)證。

慈竹WRKY轉(zhuǎn)錄因子的研究處于初級(jí)階段，其功能還需要進(jìn)一步深入研究。BeWRKY1和BeWRKY2轉(zhuǎn)錄因子的克隆、信息學(xué)分析、表達(dá)模式的分析和功能初步探究可為以慈竹品種改良提供理論依據(jù)?，F(xiàn)階段本研究室已經(jīng)成功構(gòu)建了BeWRKY1和BeWRKY2的過(guò)表達(dá)載體和RNAi表達(dá)載體，并轉(zhuǎn)化煙草和楊樹，為深入研究其功能奠定了基礎(chǔ)。

[1] Eulgem T,Somssich IE.Networks of WRKY transcription factors in defense signaling[J].Current Opinion in Plant Biology,2007,(10)：366-371.

[2] 金慧,欒雨時(shí).番茄WRKY基因的克隆與分析[J].西北農(nóng)業(yè)報(bào),2011,20(4)：96-101.

[3] Park CJ,Shin YC,Lee BJ,et al.A hot pepper gene encoding WRKY transcription factor is induced during hypersensitive response to Tobacco mosaic virus and Xanthom-onas campestris[J]. Planta,2006,223(2)：168-179.

[4] 劉戈宇,胡鳶雷,祝建波.植物WRKY蛋白家族的結(jié)構(gòu)及其功能[J].生命的化學(xué),2006,26(3)：231-233.

[5] Johnson CS，Kolevski B，Smyth DR. a trichome and seed coat development gene of Arabidopsis,encodes a WRKY transcription factor[J].Plant Cell,2002,14(6)：1359-1375.

[6] Qiu YP，Jing SJ，F(xiàn)u J，et al.Cloning and analysis of expression profile of 13 WRKY genes in rice[J]. Chinese Science Bulletin,2004,49(20)：1860-1869.

[7] 劉磊.小麥WRKY基因的克隆、表達(dá)分析及其遺傳轉(zhuǎn)化[D].華中科技大學(xué),2011.1.

[8] 崔曉偉.PheWRKY1基因的克隆與功能分析[D].中國(guó)林業(yè)科學(xué)研究院,2011,7.

[9] Tripathi P，Rabara RC，Rushton PJ. A systems biology perspective on the role of WRKY transcription factors in drought responses in plants[J].Planta,2014,239(2)：255-266.

[10] Rushton DL,Tripathi P,Rabara RC,et al.WRKY transcription factors：key components in abscisic acid signalling[J].Plant Biotechnol J，2012,10(1)：2-11.

[11] Thao NP，Tran L-SP.Potentials toward genetic engineering of drought-tolerant soybean[J].Crit Rev Biotechnol,2012,32(4)：349-362.

[12] Wu X，Shiroto Y，Kishitani S，et al.Enhanced heat and drought tolerance in transgenic rice seedlings overexpressingOsWRKY11 under the control ofHSP101 promoter[J].Plant Cell Rep，2009,28(1)：21-30.

[13] Qiu Y，Yu D.Over-expression of the stress-inducedOsWRKY45 enhances disease resistance and drought tolerance inArabidopsis[J]. Environ Exp Bot,2009,65(1)：35-47.

[14] Liu QL,Zhong M，Li S,et al. Overexpression of achrysanthemum transcription factor gene,DgWRKY3，in tobacco enhances tolerance to salt stress[J]. Plant Physiology and Biochemistry,2013(69)：27-33.

[15] Zhou QY，Tian AG，Zou HF,et al. Soybean WRKY-type transcription factor genes，GmWRKY13，GmWRKY21，andGmWRKY54，confer differential tolerance to abiotic stresses in transgenicArabidopsisplants[J].Plant Biotechnol J,2008，6(5)：486-503.

[16] Scarpeci TE，Zanor MI，Mueller-Roeber B，et al.Overexpression ofAtWRKY30 enhances abiotic stress tolerance during early growth stages inArabidopsisthaliana[J]. Plant Mol Biol,2013,83(3)：265-277.

[17] Lagace M，Matton DP. Characterization of a WRKY transcription factor expressed in late torpedo-stage embryos ofSolanumchacoense[J]. Planta，2004，219(1)：185-189.

[18] Robatzek S，Somssich IE. A new member of the Arabidopsis WRKY transcription factor family，AtWRKY6，is associated with both senescence- and defence-related processes[J]. Plant J，2001,28(2)：123-133.

[19] 吳曉宇,胡尚連,曹穎,等.慈竹CCoAOMT基因的克隆及生物信息學(xué)分析[J].南京林業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版,2012,36(3)：17-22.

[20] 段春香,董文淵,劉時(shí)才,等.慈竹無(wú)性系種群生長(zhǎng)與立地條件關(guān)系[J].林業(yè)科技開發(fā),2008,22(3)：42-44.

[21] 陳其兵,蔣瑤,盧學(xué)琴,等.四川不同地區(qū)慈竹的遺傳多樣性研究[J].西北農(nóng)林科技大學(xué)學(xué)報(bào),2009,37(6)：187-193.

[22] 劉懷,趙志模,鄧永學(xué),等.竹裂爪螨毛竹種群與慈竹種群對(duì)不同寄主植物的適應(yīng)及其生殖隔離[J]. 應(yīng)用生態(tài)學(xué)報(bào),2004,15(2)：299-302.

[23] Liu M,Qiao G,Jiang J,et al.Transcriptome Sequencing and De Novo Analysis for Ma Bamboo (DendrocalamuslatiflorusMunro)Using the Illumina Platform [J]. PloS One,2012,7(10)：e46766.

[24] Peng Z，Zhang C，Zhang Y，et al. Transcriptome Sequencing and Analysis of the Fast Growing Shoots of Moso Bamboo (Phyllostachysedulis)[J]. PloS one，2013，8(11)：e78944.

[25] Shan H，Chen SM，Jiang JF，et al. Heterologous Expression of the ChrysanthemumR2R3-MYBTranscription FactorCmMYB2 Enhances Drought and Salinity Tolerance，Increases Hypersensitivity to ABA and Delays Flowering inArabidopsisthaliana[J]. Molecular Biotechnology，2012，Volume 51(2),160-173.

[26] Luan M，Xu M，Lu Y，et al. Family-wide survey of miR169s and NF-YAs and their expression profiles response to abiotic stress in maize roots [J]. PLoS One，2014，9(3)：e91369.

[27] Liu L，Cao XL，Bai R，et al. Isolation and Characterization of the Cold-Induced Phyllostachys edulis AP2/ERF Family Transcription Factor，peDREB1 [J]. Plant Mol Biol Rep，2012,30：679-689.

[28] Livak KJ，Schmittqen TD. Analysis of relative gene expression data using real time quantitative PCR and the 2-△△Ctmethod [J]. Methods,2001,25(4)：402-408.

[29] Mochida K，Yoshida T，Sakurai T,et al.In silico analysis of transcription factor repertoire and prediction of stress responsive transcription factors in soybean[J]. DNA Res，2009,16(6)：353-369.

[30] Li SJ，Zhou X，Chen LG，et al. Functional characterization ofAradopsisbidopsisthalianaWRKY39 in heat stress[J].Mol.Cells,2010(29)：475-483.

[31] Golldack D，Luking I，Yang O. Plant tolerance to drought and salinity：stress regulating transcription factors and their functional significance in the cellular transcriptional network[J].Plant Cell Rep, 2011,30：1383-1391.

[32] 王晨.小麥基因組中WRKY轉(zhuǎn)錄因子家族基因鑒定、克隆和TaWRKY10基因的功能分析[D].華中科技大學(xué),2013,5.

[33] Babitha K，Ramu S，Pruthvi V，et al.Co-expression ofAtbHLH17 andAtWRKY28 confers resistance to abiotic stress inArabidopsis[J]. Transgenic Res,2013,22(2)：327-341.

Cloning of WRKY Transcription Factors inBambusaemeiensisand Stress-induced Expressions

LIU Hong-mei1,2，HU Shang-lian1,2*,CAO Ying1,2，LU Xue-qin1,2， XU Gang1,2，LONG Zhi-jian1,2，REN Peng1,2

(1.College of Life Science and Engineering，Southwest University of Science and Technology,Mianyang 621010,Sichuan,China；2.Engineering Research Center for Biomass Resource Utilization and Modification of Sichuan Province， Mianyang 621010，Sichuan，China)

In order to explore the biological function of WRKY inBambusaemeiensisthen modify bamboo plants through genetic engineering. Two WRKY genes were cloned based on RNA transcriptome data ofBambusaemeiensis，which were named asBeWRKY1 (GenBank：KJ462124) andBeWRKY2(GenBank：KJ462125)，and encoded 189 and 295 amino acids，respectively. Blastp analysis indicated that both of them contained a WRKY conserved domain. BeWRKY1 belonged to Ⅱc subfamily，whereas BeWRKY2，aligned withTaWRKY16 andAtWRKY39，belonged to Ⅱd subfamily. Real-time PCR analysis showed thatBeWRKY1 exhibited higher expression in stems，while the expression level ofBeWRKY2 was higher in the full expanded leaves. 200 mmol·L-1NaCl and 20% PEG 6000 stresses triggeredBeWRKY1 expression，but decreased the transcript accumulation ofBeWRKY2.The results showed thatBeWRKY1 may response to stress in short time.BeWRKY2 was more sensitive to drought and ABA stresses，but did not participate in the protection of high salt stress response.

Bambusaemeiensis;WRKY;Cloning;Bioinformatics analysis;Stress-induced expression

2016-01-25

國(guó)家自然科學(xué)基金青年基金項(xiàng)目(31400257，31400333)；四川省“十三五”重點(diǎn)攻關(guān)資助項(xiàng)目(2016NYZ0038)

劉紅梅(1988-)，女，漢，碩士研究生,從事植物遺傳與品種改良研究。通信作者：胡尚連(1966-)，女，漢，教授，博士，從事植物生理與生物技術(shù)研究。E-mail：hushanglian@126.com