龔道勇　胡尚連*　曹　穎　盧學(xué)琴　張慶波

(1.西南科技大學(xué)植物細(xì)胞工程實(shí)驗(yàn)室，綿陽(yáng)　621010； 2.四川省生物質(zhì)資源利用與改性工程技術(shù)研究中心，綿陽(yáng)　621010)

竹子在我國(guó)種類(lèi)多，分布廣，具有生產(chǎn)周期短和易更新的特點(diǎn)，且是一種不可或缺的可再生資源，有著廣闊的開(kāi)發(fā)利用前景[1]。慈竹(Bambusaemeiensis)為叢生竹類(lèi)，屬于禾本科(Gramineae)竹亞科(Bambusoideae)慈竹屬(Bambusaemeiensis)，是西南地區(qū)的優(yōu)勢(shì)竹種，其節(jié)間長(zhǎng)，纖維素含量多，纖維長(zhǎng)徑比較大，分子聚合度相對(duì)較高，是造紙的良好材料[2～3]。但是慈竹本身易受冷凍脅迫、土壤干旱等影響，從而造成慈竹大面積減產(chǎn)。因此，探究如何提高它對(duì)低溫或干旱脅迫的抗性十分重要。

以往對(duì)慈竹的研究主要集中在基因的克隆[4]、生理生化[5]、遺傳的多樣性[6]等方面的研究，而有關(guān)慈竹bZIP基因克隆表達(dá)方面的研究鮮有報(bào)道。有研究表明，bZIP轉(zhuǎn)錄因子廣泛參與植物的生長(zhǎng)發(fā)育：例如在小麥的187個(gè)bZIP基因中，發(fā)現(xiàn)有48個(gè)與花藥的發(fā)育有關(guān)[7]；擬南芥中C/S1亞族的bZIP轉(zhuǎn)錄因子在能量不足或者環(huán)境脅迫下參與植株生長(zhǎng)的重新分配[8]；而G亞族的GIP1則參與到植株早期的生長(zhǎng)發(fā)育[9]。光信號(hào)的轉(zhuǎn)導(dǎo)：擬南芥中的GBF1在藍(lán)光誘導(dǎo)下參與光形態(tài)的建成，如子葉擴(kuò)張的正調(diào)控[10]；而對(duì)于一些藻類(lèi)植物，在bZIP轉(zhuǎn)錄因子參與下，它可以不利用光合通路，而完成光合作用[11]。種子的萌發(fā)：如分析蓖麻種子時(shí)，發(fā)現(xiàn)了49個(gè)RcbZIP基因與種子萌發(fā)和保藏相關(guān)[12]；而水稻中OsABI5可以提高水稻種子的萌發(fā)率[13]。基因的表達(dá)：在西紅柿中，進(jìn)行了全基因組的系統(tǒng)鑒定，為蛋白質(zhì)功能分析提供理論依據(jù)[14]。病蟲(chóng)害的防御：如對(duì)四種大豆的bZIP基因進(jìn)行銹病脅迫實(shí)驗(yàn)，獲得了新的基因型可以提高對(duì)銹病的抗性[15]；在胡椒中，發(fā)現(xiàn)了一個(gè)CAbZIP1基因，它可以增強(qiáng)對(duì)胡椒病原體的抗性[16]。生物和非生物脅迫應(yīng)答：如Atkinson等從基因的角度，分析了植物生物脅迫和非生物脅迫的相互作用是通過(guò)激素信號(hào)通路使之結(jié)合在一起[17]；在小麥中bZIP轉(zhuǎn)錄因子是通過(guò)分子構(gòu)象的不同來(lái)影響其結(jié)構(gòu)與功能[18]；依賴(lài)于ABA信號(hào)的bZIP轉(zhuǎn)錄因子通過(guò)Ser/Thr激酶參與到非生物脅迫中，如水稻ABF1對(duì)冷脅迫的響應(yīng)[19]。以及ABA敏感性等各種信號(hào)的響應(yīng)[20]。鑒于此，本次研究基于慈竹轉(zhuǎn)錄組數(shù)據(jù)庫(kù)[21]，利用克隆得到兩條序列，對(duì)其進(jìn)行生物信息學(xué)分析、組織表達(dá)以及非生物脅迫表達(dá)分析，為慈竹的抗凍、耐鹽、抗旱等提供一定的理論依據(jù)。

1　材料和方法

1.1　植物材料

采自西南科技大學(xué)生命科學(xué)與工程學(xué)院竹類(lèi)資源圃的慈竹筍(高100 cm)、當(dāng)年生且完全抽枝的慈竹莖稈、展開(kāi)葉和未展開(kāi)葉作為試驗(yàn)材料，于液氮中冷凍后置于-80℃保存?zhèn)溆谩?/p>

1.2　主要試劑和菌種

寶生物(大連)TaKaRa公司：DL2000 DNA Marker、LA-Taq聚合酶、pMD19-T載體、IPTG、PrimeScript? RT reagent Kit(Perfect DNA)反轉(zhuǎn)錄試劑盒、ApaⅠ、XbaⅠ、KpnⅠ等;天根公司：質(zhì)粒小提試劑盒、X-Gal和膠回收試劑盒;OMEGA BIO-TEK公司：Plant RNA Kit試劑盒;BioBRK公司：Green View染料;GeneTech公司：酵母提取物、瓊脂糖和胰蛋白胨;大腸桿菌(E.coli)DH5α為本試驗(yàn)室制備保存。

1.3　試驗(yàn)方法

1.3.1慈竹總RNA的提取和cDNA的合成

以-80℃超低溫冰箱保存的慈竹筍(高100 cm)、莖稈、展開(kāi)葉和未展開(kāi)葉為材料，按照OMEGA公司Plant RNA Kit試劑盒說(shuō)明書(shū)提取總RNA。按照寶生物公司的PrimeScript?RT reagent Kit(Perfect DNA)反轉(zhuǎn)錄試劑盒將RNA反轉(zhuǎn)錄合成cDNA。

1.3.2慈竹BebZIP2和BebZIP6的克隆

從慈竹筍轉(zhuǎn)錄組數(shù)據(jù)庫(kù)中篩選出2條序列，采用Premier 5.10軟件設(shè)計(jì)兩對(duì)引物(BebZIP2-F和BebZIP2-R；BebZIP6-F和BebZIP6-R，表1)，對(duì)其進(jìn)行克隆。按照擴(kuò)增反應(yīng)程序：95℃預(yù)變性3 min，95℃變性30 s，55℃*BebZIP2基因退火溫度為55℃，BebZIP6為58℃。退火30 s，72℃延伸1 min，30個(gè)循環(huán)，72℃延伸10 min，4℃保持。1.0%的瓊脂糖凝膠電泳，回收后用大連寶生物T4ligase試劑盒進(jìn)行pMD19-T載體連接，42℃熱激轉(zhuǎn)化E.coliDH5α,篩選陽(yáng)性克隆酶切驗(yàn)證后送至上海英濰捷基公司進(jìn)行測(cè)序，保留測(cè)序正確的質(zhì)粒進(jìn)行下一步實(shí)驗(yàn)。

1.4　生物信息學(xué)分析

從NCBI數(shù)據(jù)庫(kù)(http://www.ncbi.nlm.nih.gov/guide/sequence-analysis/)中獲得其它物種bZIP基因的功能性氨基酸和核苷酸序列，用DAMBE檢測(cè)飽和度[22]；為不飽和狀態(tài)，適合建樹(shù)。用DNAMAN軟件和在線(xiàn)程序Clustal Omega(http://www.ebi.ac.uk/Tools/msa/clustalo)對(duì)基因序列及其編碼的氨基酸序列進(jìn)行多序列比對(duì)分析；用Mega7.0軟件采用最大似然法(Maximum Likehood,ML)構(gòu)建系統(tǒng)發(fā)育樹(shù)。本次研究選擇Models(Find Best DNA/Protein Models)為JTT+G+F，Bootstrap參數(shù)設(shè)為1 000 replicates；用MEME(http://meme-suite.org/index.html)、Phyre2.0(http://www.sbg.bio.ic.ac.uk/～phyre2/html/)等在線(xiàn)軟件工具分別對(duì)bZIP蛋白序列進(jìn)行保守基序和三級(jí)結(jié)構(gòu)分析；用Yloc[23～24](http://abi.inf.uni-tuebingen.de/Services/YLoc/webloc.cgi)進(jìn)行亞細(xì)胞定位預(yù)測(cè)分析(選擇模型：Plants YLoc-HighRes)。將已克隆得到的慈竹2條BebZIP基因的編碼框CDS序列提交至Selecton(http://selecton.tau.ac.il/)在線(xiàn)程序[25～26]，設(shè)置模型為力學(xué)—經(jīng)驗(yàn)型模型(MEC)，對(duì)這2個(gè)基因進(jìn)行選擇性壓力分析。

1.5　慈竹BebZIP基因組織表達(dá)模式分析

根據(jù)測(cè)序的慈竹BebZIP2和BebZIP6基因的全長(zhǎng)序列，設(shè)計(jì)real-time PCR引物(BebZIP2-RT-F和BebZIP2-RT-R；BebZIP6-RT-F和BebZIP6-RT-R，表1)。分別提取慈竹各組織的總RNA并反轉(zhuǎn)錄合成cDNA；在iQ Multicolor Real-Time PCR自動(dòng)擴(kuò)增儀進(jìn)行反應(yīng)，每個(gè)樣品重復(fù)3次，反應(yīng)體系20 μL。以?xún)?nèi)參基因Tubulin為對(duì)照，用2-ΔΔCT計(jì)算法[27]計(jì)算其相對(duì)表達(dá)量。

表1　慈竹BebZIP2和BebZIP6基因克隆和Real-time PCR分析所用引物

1.6　慈竹BebZIP基因非生物脅迫表達(dá)分析

采用1/10 Hoagland營(yíng)養(yǎng)液,將慈竹幼苗培養(yǎng)30 d左右，選擇生長(zhǎng)情況幾乎一致且生長(zhǎng)良好的幼苗作為試驗(yàn)材料。每個(gè)處理分別設(shè)置6個(gè)試驗(yàn)梯度：NaCl(0、50、100、150、200、250 mmol·L-1)、ABA(0、5、10、50、100、200 μmol·L-1)、PEG6000(0%、5%、10%、15%、20%、25%)進(jìn)行預(yù)實(shí)驗(yàn)，每個(gè)處理重復(fù)3次(每個(gè)重復(fù)取3個(gè)獨(dú)立植株)。

預(yù)實(shí)驗(yàn)結(jié)果顯示：250 mmol·L-1NaCl處理3 d后，慈竹幼苗枯萎死亡，200 mmol·L-1NaCl處理的幼苗葉尖部位變黑，其他濃度的NaCl處理后，整株幼苗無(wú)明顯變化；200 μmol·L-1ABA處理2 d后，幼苗的葉邊緣的綠色明顯變?yōu)辄S色，而其他濃度的ABA處理的幼苗均無(wú)明顯變化；20% PEG6000處理2 d后，可觀察到幼苗葉呈卷曲狀，25% PEG6000處理的幼苗葉卷曲比較嚴(yán)重且干癟，其他濃度的PEG6000處理后，整株幼苗無(wú)明顯變化。根據(jù)預(yù)實(shí)驗(yàn)結(jié)果最終確定以NaCl(0、200 mmol·L-1)、ABA(0、200 μmol·L-1)、PEG6000(0%、20%)進(jìn)行后續(xù)的脅迫試驗(yàn)，每個(gè)處理重復(fù)3次(每個(gè)重復(fù)取3個(gè)獨(dú)立植株)，分別處理0(未處理的對(duì)照組)、6、12、24、48 h及7 d后，選取脅迫處理后的慈竹葉作為試驗(yàn)材料，液氮處理后于-80℃冰箱保存?zhèn)溆谩Ｌ崛NA;合成cDNA;以?xún)?nèi)參基因Tubulin為對(duì)照，用2-ΔΔCT計(jì)算法計(jì)算其相對(duì)表達(dá)量。

2　結(jié)果與分析

2.1　慈竹BebZIP2和BebZIP6基因的克隆

用幼嫩的慈竹筍cDNA為模版，通過(guò)設(shè)計(jì)基因特異引物，經(jīng)PCR擴(kuò)增獲得了2條DNA序列(圖1)，將其克隆到pMD19-T載體上，獲得陽(yáng)性克隆并測(cè)序。用NCBI在線(xiàn)工具Conserved分別對(duì)其進(jìn)行保守結(jié)構(gòu)預(yù)測(cè)并將這兩條序列提交至NCBI數(shù)據(jù)庫(kù)注冊(cè)；保守結(jié)構(gòu)預(yù)測(cè)屬于bZIP超家族，命名為BebZIP2和BebZIP6(GenBank注冊(cè)號(hào)分別為KU554693和KU560573)。其中BebZIP2和BebZIP6編碼區(qū)長(zhǎng)度分別為504 bp(編碼167個(gè)氨基酸)和720 bp(編碼239個(gè)氨基酸)。

2.2　慈竹BebZIP2和BebZIP6蛋白進(jìn)化樹(shù)分析及多序列比對(duì)分析

用慈竹轉(zhuǎn)錄組數(shù)據(jù)庫(kù)中篩選出的bZIP(BebZIP1-BebZIP6)蛋白序列，結(jié)合單子葉植物：水稻(Oryzasativa：OsZIP-1a；OsZIP-2b；TRAB1；RF2a；RF2b；RISBZ1；RISBZ4；RISBZ5；RITA-1；OsBZ8)；小麥(Triticumaestivum：EmBP-1；HBP-1a(1)；HBP-1b(c1)；HBP-1b(c38);HALF1；SPA)；玉米(Zeamays：CPRF-2；DFL1；LG2；OHP1；OHP2；Maize O2；ZmGBF1；ZmBP-1a；OBF3.1；OBF3.2)和雙子葉植物擬南芥(Arabidopsisthaliana：BZO2H1；BZO2H2；BZO2H3；BZO2H4)；大豆(Glycinemax：HBF1)等已知功能的bZIP蛋白構(gòu)建進(jìn)化樹(shù)，分析得知BebZIP2與玉米的CPRF-2置信度達(dá)到了65%；BebZIP2和BebZIP6與水稻C亞族的RISBZ5蛋白聚在一枝，置信度為90%(圖2)。CPRF-2即bZIP17，它與光誘導(dǎo)相關(guān)[28]。RISBZ5即OsbZIP52，它參與水稻對(duì)脅迫的響應(yīng)，在低溫脅迫下其表達(dá)量可迅速提升，說(shuō)明它可能是調(diào)控冷和滲透壓脅迫途徑中的一部分，同時(shí)相關(guān)研究表明，OsbZIP52基因的過(guò)表達(dá)可以提高水稻幼苗對(duì)冷和干旱脅迫的敏感程度[29]。因而，推測(cè)BebZIP2和BebZIP6基因的功能參與到慈竹對(duì)冷和干旱等非生物脅迫的響應(yīng)。

圖1　慈竹bZIP基因PCR擴(kuò)增產(chǎn)物　M.DL2000 Maker；1.BebZIP6 PCR擴(kuò)增產(chǎn)物；2.BebZIP2 PCR擴(kuò)增產(chǎn)物Fig.1　Electrophoresis of PCR products of bZIP genes in B.emeiensis　M. DL2000 Maker; 1.Amplification products of BebZIP6 PCR; 2.Amplification products of BebZIP2 PCR

選取與BebZIP2和BebZIP6同源性較高的氨基酸序列，用DNAMAN軟件和在線(xiàn)程序Clustal Omega進(jìn)行多序列比對(duì)分析(圖3)。可以發(fā)現(xiàn)它們存在典型的堿性亮氨酸拉鏈區(qū)域[N—X(7)—R/K—X(9)—L—X(6)—L—X(6)—L]，圖中黃色區(qū)域且該區(qū)域非常的保守。

圖2　不同物種bZIP蛋白的進(jìn)化樹(shù)分析Fig.2　Phylogenetic analysis of deduced amino acids of bZIP proteins from different plants

圖3　BebZIP氨基酸序列與其他物種bZIP蛋白氨基酸序列多重比對(duì)結(jié)果　“*”表示保守氨基酸；“：”表示保守替換；“.”表示非保守替換Fig.3　Multiple alignment analysis of bZIP conservative structure domain in different plants　“*” represents a conserved amino acid; “:” represents a conservative replacement; “.” represents non conservative replacement

圖4　慈竹BebZIP2和BebZIP6蛋白保守基序分析Fig.4　The conserved motif analysis of BebZIP2 and BebZIP6 protein

圖5　bZIP蛋白質(zhì)的三級(jí)結(jié)構(gòu)預(yù)測(cè)　A. c2wt7B的三級(jí)結(jié)構(gòu)；B. RISBZ5的三級(jí)結(jié)構(gòu)；C. BebZIP2的三級(jí)結(jié)構(gòu)；D. BebZIP6的三級(jí)結(jié)構(gòu)Fig.5　bZIP protein tertiary structure prediction　A. Tertiary structure of model c2wt7B; B. Tertiary structure of RISBZ5; C. Tertiary structure of BebZIP2; D. Tertiary structure of BebZIP6

2.3　慈竹BebZIP氨基酸保守基序分析

用MEME軟件分析慈竹2條bZIP氨基酸序列的保守基序，設(shè)置基序數(shù)量為6個(gè)，寬度不限。結(jié)果顯示，這些氨基酸序列均含有高度保守的motif1-motif5核心保守基序(圖4)。其中，BebZIP6的結(jié)構(gòu)與水稻RISBZ5蛋白的保守基序幾乎完全相同，與進(jìn)化分析結(jié)果一致。BebZIP2的motif2、motif4和motif5則組合在了一起。

2.4　慈竹BebZIP基因編碼蛋白三級(jí)結(jié)構(gòu)預(yù)測(cè)分析

蛋白質(zhì)的三級(jí)結(jié)構(gòu)是在二級(jí)結(jié)構(gòu)的基礎(chǔ)上進(jìn)一步折疊盤(pán)繞而成。利用已經(jīng)被X-Ray解析出晶體模型的c2wt7B(NDB:NA0777)為模型[30],c2wt7B具有典型的亮氨酸拉鏈結(jié)構(gòu)。用Phyre 2.0對(duì)BebZIP2和BebZIP6進(jìn)行三級(jí)結(jié)構(gòu)進(jìn)行預(yù)測(cè)。結(jié)果顯示這兩個(gè)蛋白與c2wt7B的三級(jí)構(gòu)型高度相似，猶如一條拉鏈。如圖5可以看出，BebZIP2和BebZIP6的三級(jí)結(jié)構(gòu)主要由α-螺旋組成。Yloc亞細(xì)胞定位預(yù)測(cè)結(jié)果表明，BebZIP2主要集中于細(xì)胞核的概率為99.99%(置信度為1)，BebZIP6主要集中于細(xì)胞核的概率為100.00%(置信度為1)；實(shí)際情況需進(jìn)一步實(shí)驗(yàn)去驗(yàn)證。

2.5　慈竹BebZIP基因的選擇性進(jìn)化分析

運(yùn)用Selecton程序中的MEC模型在慈竹bZIP氨基酸序列上探測(cè)正選擇位點(diǎn)。蛋白質(zhì)一級(jí)結(jié)構(gòu)位點(diǎn)的標(biāo)注顏色大部分為白色到深紫色之間，表明這些位點(diǎn)經(jīng)歷了強(qiáng)的凈化選擇位點(diǎn)，而標(biāo)注為黃色的位點(diǎn)，說(shuō)明這些位點(diǎn)受到陽(yáng)性選擇壓力。由圖6可見(jiàn)在MEC模型下，慈竹bZIP基因中凈化選擇壓力和陽(yáng)性選擇壓力共同主導(dǎo)了它們的進(jìn)化。

圖6　慈竹bZIP基因選擇性壓力檢測(cè)Fig.6　bZIP genes selective pressure detection in B.emeiensis

2.6　慈竹BebZIP基因的組織表達(dá)分析

組織表達(dá)分析結(jié)果顯示，這2條BebZIP基因在慈竹筍、展開(kāi)葉和未展開(kāi)葉以及莖中均有表達(dá)，但其在各個(gè)部位的表達(dá)量有差異(圖7)?？傮w上看，它們?cè)诓煌M織中的表達(dá)趨勢(shì)不盡相同，但均在未展開(kāi)葉中表達(dá)量較高，在筍和莖中表達(dá)量的差異不大。同一基因在不同組織中的表達(dá)量差異較大，表達(dá)量的大小為未展開(kāi)葉>展開(kāi)葉>莖>筍。

2.7　慈竹BebZIP基因的非生物脅迫分析

經(jīng)過(guò)脅迫處理后，觀察預(yù)實(shí)驗(yàn)幼苗(圖8)，與CK比較，200 mmol·L-1NaCl處理2 d后,幼苗葉尖部位變黑；20% PEG-6000處理2天后，可觀察到幼苗葉呈卷曲狀；200 μmol·L-1ABA處理2 d后，幼苗的葉邊緣的綠色明顯變?yōu)辄S色。

圖7　實(shí)時(shí)定量PCR分析2個(gè)BebZIP基因在不同組織中的表達(dá)Fig.7　Quantitative real-time PCR analysis of the expression pattern of 6 BebZIP genes in different tissues

圖8　脅迫處理后植株的變化Fig.8　The phenotype of the plant after treatments

在NaCl(200 mmol·L-1)、PEG6000(20%)和ABA(200 μmol·L-1)3種脅迫條件的處理下(圖9)，BebZIP2基因的表達(dá)情況基本一致，表現(xiàn)為先下降，后上升，最后再下降。在NaCl、PEG6000和ABA處理6 h后其表達(dá)量顯著下降，但隨著脅迫時(shí)間的延長(zhǎng)，它的表達(dá)量緩慢上升，脅迫7天后BebZIP2的表達(dá)量又開(kāi)始下降。

然而，BebZIP6在3種脅迫環(huán)境下的表達(dá)情況與BebZIP2明顯不同，在NaCl和PEG處理?xiàng)l件下，大多數(shù)脅迫時(shí)間點(diǎn)下其表達(dá)水平與對(duì)照初始水平差異不大，且脅迫7天后其表達(dá)量均與對(duì)照基本持平。而在ABA處理?xiàng)l件下，BebZIP6基因的表達(dá)量呈現(xiàn)波浪線(xiàn)形狀變化，依次表現(xiàn)為上升，下降，上升，下降，最后上升到與對(duì)照組水平基本持平。

圖9　BebZIP2和BebZIP6對(duì)脅迫處理下的表達(dá)變化A. NaCl處理；B. PEG6000處理；C. ABA處理Fig.9　Differential expression of BebZIP2 and BebZIP6 under the different stress treatments　A. NaCl treatment; B. PEG6000 treatment; C. ABA treatment

3　討論

bZIP轉(zhuǎn)錄因子家族在整個(gè)真核生物界中扮演著重要角色。在植物中，bZIP可以調(diào)控包括光形態(tài)建成、葉和種子的發(fā)育、能量平衡、以及生物和非生物脅迫響應(yīng)等一系列重要的生物學(xué)過(guò)程[31]。由于慈竹的這兩個(gè)基因?qū)儆贑亞族。有研究表明大豆C亞族bZIP轉(zhuǎn)錄因子GmbZIP62在ABA信號(hào)傳導(dǎo)中通過(guò)上調(diào)ABI1和ABI2等逆境相關(guān)基因的表達(dá)，從而可能在逆境脅迫中起作用[32]；玉米C亞族的Opaque2通過(guò)和PBF蛋白的相互作用調(diào)控種子儲(chǔ)存蛋白的產(chǎn)量[33]；歐芹C亞族的CPFR-2和G/HBF-1可能還與對(duì)環(huán)境和病害侵襲的防御反應(yīng)有關(guān)[34～35]；在水稻中OsbZIP23的過(guò)表達(dá)能顯著提高水稻的耐鹽性、抗旱性和ABA敏感性；在對(duì)照植株中，水稻的耐鹽性、抗旱性和ABA敏感性也隨之降低比較明顯，這表明OsbZIP23基因在水稻的耐鹽和抗旱過(guò)程中起到了非常重要的作用[36]。

慈竹這兩個(gè)基因與水稻C亞族的OsbZIP52同源性較高，而亞細(xì)胞定位表明OsbZIP52定位于細(xì)胞核[28]；推測(cè)BebZIP2和BebZIP6也定位于為細(xì)胞核，這與Yloc亞細(xì)胞定位預(yù)測(cè)結(jié)果相一致，但需要進(jìn)一步實(shí)驗(yàn)去驗(yàn)證。水稻中，OsbZIP52作為反式激活因子能與順式作用元件G-Box特異性結(jié)合，提高對(duì)寒冷和干旱的敏感性[29]。在水稻中，通過(guò)Realtime PCR分析OsbZIP52的過(guò)表達(dá)，其非生物脅迫相關(guān)基因，如OsLEA3,OsTPP1,Rab25,gp1precursor,β-gal,LOC_Os05g11910和LOC_Os05g39250等基因都下調(diào)了；這暗示OsbZIP52/RISBZ5可能在寒冷和干旱中作為負(fù)調(diào)控因子[29]。同樣推測(cè)慈竹的這兩個(gè)基因也可能有類(lèi)似功能。

在脅迫處理下，BebZIP2基因相對(duì)表達(dá)量下降的很快，很可能是由于它的啟動(dòng)更快，使得它對(duì)脅迫非常敏感能夠迅速反饋給信號(hào)通路，從而減少它的表達(dá)量；而B(niǎo)ebZIP6則相對(duì)滯后，使得它在初始條件下相對(duì)表達(dá)量較低，在脅迫處理下，對(duì)脅迫處理很敏感，為了適應(yīng)環(huán)境的改變其相對(duì)表達(dá)量上升，待適應(yīng)環(huán)境后又開(kāi)始下降，為了穩(wěn)定持續(xù)地適應(yīng)環(huán)境，它需要多次的調(diào)控其表達(dá)，來(lái)回的上下波動(dòng)最終趨于恒定，這有助于提高慈竹對(duì)多種非生物脅迫的適應(yīng)能力。這些推斷是否正確均需進(jìn)一步實(shí)驗(yàn)去驗(yàn)證。

1.夏登云.慈竹組培快繁體系的研究[D].合肥:安徽農(nóng)業(yè)大學(xué),2009.

Xia D Y.Study on tissure cultule and rapid propagation system ofNeosinocalmusaffinis(Rendle) Keng F.[D].Hefei:Anhui Agricultural University,2009.

2.火樹(shù)華.樹(shù)木學(xué):2版[M].北京:中國(guó)林業(yè)出版社,1992:456-459.

Huo S H.Dendrology:2nd ed[M].Beijing:China Forestry Publishing House,1992:456-459.

3.劉曉光,劉海英,張東遠(yuǎn),等.竹子生物技術(shù)育種研究進(jìn)展[J].河北農(nóng)業(yè)大學(xué)學(xué)報(bào),2002,25(S):151-154.

Liu X G,Liu H Y,Zhang D Y,et al.Advances of biotechnological approach on bamboo breeding[J].Journal of Agricultural University of Hebei,2002,25(S):151-154.

4.胡尚連,曹穎,黃勝雄,等.慈竹4CL基因的克隆及其生物信息學(xué)分析[J].西北農(nóng)林科技大學(xué)學(xué)報(bào):自然科學(xué)版,2009,37(8):204-210.

Hu S L,Cao Y,Huang S X,et al.Cloning and bioinformation analysis of 4CL gene inNeosinocalamusaffinis[J].Journal of Northwest A&F University:Natural Science Edition,2009,37(8):204-210.

5.胡尚連,蔣瑤,陳其兵,等.四川2種叢生竹理化特性及纖維形態(tài)研究[J].植物研究,2010,30(6):708-712.

Hu S L,Jiang Y,Chen Q B,et al.Physical and chemical properties of 2 species in bamboos from the different regions in sichuan province[J].Bulletin of Botanical Research,2010,30(6):708-712.

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

Chen Q B,Jiang Y,Lu X Q,et al.Assessment of genetic diversity inNeosinocalamusaffinisfrom different regions of Sichuan province[J].Journal of Northwest A&F University:Natural Science Edition,2009,37(6):187-193.

7.Li X Y,Gao S Q,Tang Y M,et al.Genome-wide identification and evolutionary analyses of bZIP transcription factors in wheat and its relatives and expression profiles of anther development related TabZIP,genes[J].BMC Genomics,2015,16:976.

8.Weltmeier F,Rahmani F,Ehlert A,et al.Expression patterns within theArabidopsisC/S1 bZIP transcription factor network:availability of heterodimerization partners controls gene expression during stress response and development[J].Plant Molecular Biology,2009,69(1-2):107-119.

9.Shaikhali J.GIP1 protein is a novel cofactor that regulates DNA-binding affinity of redox-regulated members of bZIP transcription factors involved in the early stages ofArabidopsisdevelopment[J].Protoplasma,2014,252(3):867-883.

10.Mallappa C,Yadav V,Negi P,et al.A basic leucine zipper transcription factor,G-box-binding factor 1,regulates blue light-mediated photomorphogenic growth inArabidopsis[J].Journal of Biological Chemistry,2006,281(31):22190-22199.

11.Takahashi F.Blue-light-regulated transcription factor,Aureochrome,in photosynthetic stramenopiles[J].Journal of Plant Research,2016,129(2):189-197.

12.Jin Z W,Xu W,Liu A Z.Genomic surveys and expression analysis ofbZIPgene family in castor bean(RicinuscommunisL.)[J].Planta,2014,239(2):299-312.

13.Zou M J,Guan Y C,Ren H B,et al.A bZIP transcription factor,OsABI5,is involved in rice fertility and stress tolerance[J].Plant Molecular Biology,2008,66(6):675-683.

14.Li D Y,Fu F Y,Zhang H J,et al.Genome-wide systematic characterization of the bZIP transcriptional factor family in tomato(SolanumlycopersicumL.)[J].BMC Genomics,2015,16:771.

15.Alves M S,Soares Z G,Vidigal P M P,et al.Differential expression of four soybean bZIP genes duringPhakopsorapachyrhiziinfection[J].Functional & Integrative Genomics,2015,15(6):685-696.

16.Lee S C,Choi H W,Hwang I S,et al.Functional roles of the pepper pathogen-induced bZIP transcription factor,CAbZIP1,in enhanced resistance to pathogen infection and environmental stresses[J].Planta,2006,224(5):1209-1225.

17.Atkinson N J,Urwin P E.The interaction of plant biotic and abiotic stresses:from genes to the field[J].Journal of Experimental Botany,2012,63(10):3523-3544.

18.Sornaraj P,Luang S,Lopato S,et al.Basic leucine zipper(bZIP) transcription factors involved in abiotic stresses:a molecular model of a wheat bZIP factor and implications of its structure in function[J].Biochimica et Biophysica Acta(BBA)-General Subjects,2016,1860(1):46-56.

19.Banerjee A,Roychoudhury A.Abscisic-acid-dependent basic leucine zipper(bZIP) transcription factors in plant abiotic stress[J].Protoplasma,2017,254(1):3-16.

20.Tu M X,Wang X H,Huang L,et al.Expression of a grape bZIP transcription factor,VqbZIP39,in transgenicArabidopsisthalianaconfers tolerance of multiple abiotic stresses[J].Plant Cell,Tissue and Organ Culture(PCTOC),2016,125(3):537-551.

21.楊傳鳳,曹穎,胡尚連,等.基于慈竹轉(zhuǎn)錄組MYB基因的克隆及脅迫誘導(dǎo)表達(dá)[J].森林與環(huán)境學(xué)報(bào),2015,35(1):60-66.

Yang C F,Cao Y,Hu S L,et al.Cloning of MYB gene based onBambusaemeiensistranscriptome and their induced expression under stresses[J].Journal of Forest and Environment,2015,35(1):60-66.

22.Xia X H,Lemey P.Assessing substitution saturation with DAMBE[M].//Lemey P,Salemi M,Vandamme A M.The phylogenetic handbook.2nd ed.Cambridge,UK:Cambridge University Press,2009:611-626.

23.Briesemeister S,Rahnenführer J,Kohlbacher O.Going from where to why-interpretable prediction of protein subcellular localization[J].Bioinformatics,2010,26(9):1232-1238.

24.Briesemeister S,Rahnenführer J,Kohlbacher O.YLoc-an interpretable web server for predicting subcellular localization[J].Nucleic Acids Research,2010,38(S2):W497-W502.

25.Stern A,Doron-Faigenboim A,Erez E,et al.Selecton 2007:advanced models for detecting positive and purifying selection using a Bayesian inference approach[J].Nucleic Acids Research,2007,35(S2):W506-W511.

26.張慶波,胡尚連,徐剛,等.毛竹bZIP轉(zhuǎn)錄因子的基因結(jié)構(gòu)與進(jìn)化分析[J].森林與環(huán)境學(xué)報(bào),2016,36(1):54-61.

Zhang Q B,Hu S L,Xu G,et al.Analysis on gene structure and evolution of bZIP transcription factor ofPhyllostachysedulis[J].Journal of Forest and Environment,2016,36(1):54-61.

27.Pfaffl M W.A new mathematical model for relative quantification in real-time RT-PCR[J].Nucleic Acids Research,2001,29(9):e45.

28.Schnable P S,Ware D,Fulton R S,et al.The B73 maize genome:complexity,diversity,and dynamics[J].Science,2009,326(5956):1112-1115.

29.Liu C T,Wu Y B,Wang X P.bZIP transcription factorOsbZIP52/RISBZ5:a potential negative regulator of cold and drought stress response in rice[J].Planta,2012,235(6):1157-1169.

30.Pogenberg V,Textor L C,Vanhille L,et al.Design of a bZip transcription factor with homo/heterodimer-induced DNA-binding preference[J].Structure,2014,22(3):466-477.

31.Corrêa L G,Riao-Pachón D M,Schrago C G,et al.The role of bZIP transcription factors in green plant evolution:adaptive features emerging from four founder genes[J].PLoS One,2008,3(8):e2944.

32.Liao Y,Zou H F,Wei W,et al.SoybeanGmbZIP44,GmbZIP62 andGmbZIP78 genes function as negative regulator of ABA signaling and confer salt and freezing tolerance in transgenicArabidopsis[J].Planta,2008,228(2):225-240.

33.Vicente-Carbajosa J,Moose S P,Parsons R L,et al.A maize zinc-finger protein binds the prolamin box in zein gene promoters and interacts with the basic leucine zipper transcriptional activator Opaque2[J].Proceedings of the National Academy of Sciences of the United States of America,1997,94(14):7685-7690.

34.Dr?ge-Laser W,Kaiser A,Lindsay W P,et al.Rapid stimulation of a soybean protein-serine kinase that phosphorylates a novel bZIP DNA-binding protein,G/HBF-1,during the induction of early transcription-dependent defenses[J].The EMBO Journal,1997,16(4):726-738.

35.Weisshaar B,Armstrong G A,Block A,et al.Light-inducible and constitutively expressed DNA-binding proteins recognizing a plant promoter element with functional relevance in light responsiveness[J].The EMBO Journal,1991,10(7):1777-1786.

36.陳祥福,潘陽(yáng)露,魏娟娟,等.植物bZIP轉(zhuǎn)錄因子及其耐鹽性調(diào)控研究進(jìn)展[J].長(zhǎng)江蔬菜,2014,(20):4-7.

Chen X F,Pan Y L,Wei J J,et al.Research progress of plant bZIP transcription factors and their regulation on salt resistance[J].Journal of Changjiang Vegetables,2014,(20):4-7.