康桂娟 黎瑜 曾日中 聶智毅

（中國熱帶農(nóng)業(yè)科學院橡膠研究所 農(nóng)業(yè)部橡膠樹生物學與遺傳資源利用重點實驗室和省部共建國家重點實驗室培育基地—海南省熱帶作物栽培生理學重點實驗室，儋州 571737）

植物特有的NAC（NAM，ATAF1/2和CUC2）轉(zhuǎn)錄因子組成了植物基因組中最大的轉(zhuǎn)錄因子家族之一［1]。膜結(jié)合NAC轉(zhuǎn)錄因子是植物NAC轉(zhuǎn)錄因子家族中一類C端具有跨膜結(jié)構(gòu)域的轉(zhuǎn)錄調(diào)控因子，該類膜結(jié)合NAC蛋白依靠此結(jié)構(gòu)錨著在細胞膜或內(nèi)質(zhì)網(wǎng)膜上［3-8]。到目前為止，已經(jīng)在擬南芥（Arabidopsis thaliana）、大豆（Glycine maxL.）、玉米（Zea mays）、葡萄（Vitis viniferaL.）、水稻（Oryza sative）和毛果楊（Populus trichocarpa）等多種陸生植物中發(fā)現(xiàn)具有多個NTLs［1，4，6，8-11]。

NAC家族轉(zhuǎn)錄因子廣泛參與調(diào)控植物體生長發(fā)育、生物和非生物脅迫應答以及激素信號轉(zhuǎn)導過程［1，12-13]。越來越多的研究發(fā)現(xiàn)膜結(jié)合NAC轉(zhuǎn)錄因子的表達受非生物脅迫的誘導，在植物的生長發(fā)育和逆境脅迫應答中發(fā)揮重要作用［1，6，8]。Kim等［5，8]研究發(fā)現(xiàn)擬南芥NTLs基因受多種生物和非生物脅迫條件誘導，并且可以調(diào)節(jié)擬南芥的開花時間。研究結(jié)果顯示，水稻OsNTL基因與擬南芥相似，都可以被多種非生物脅迫誘導表達［8]。Puranik等［1]研究發(fā)現(xiàn)谷子膜結(jié)合SiNAC基因可能與脅迫和生長發(fā)育調(diào)控相關(guān)，可以被滲透脅迫、鹽脅迫、乙烯和茉莉酸甲酯誘導表達。本實驗室已經(jīng)從橡膠樹中克隆了HbNTL1基因的全長cDNA序列，在此基礎上，利用基因組步移（Genome Walker）方法克隆了該基因的啟動子序列，并通過生物信息學軟件分析其核心啟動子結(jié)構(gòu)、順式作用元件及其生物學功能，旨在為進一步研究HbNTL1轉(zhuǎn)錄因子的表達調(diào)控以及其在橡膠樹抵御逆境脅迫中的作用提供信息和依據(jù)。

1 材料與方法

1.1 材料

巴西橡膠樹（Hevea brasiliensis）無性系熱研7-33-97種植在中國熱帶農(nóng)業(yè)科學院實驗農(nóng)場，幼葉液氮研磨后凍存在-80℃，用于基因組DNA的提取。

1.2 方法

1.2.1 巴西橡膠樹基因組DNA的提取 橡膠樹基因組DNA的提取采用改進的CTAB法［16]。

1.2.2HbNTL1基因啟動子的克隆 啟動子克隆采用基因組步移的方法，Genome Walker文庫由本研究組夏可燦等按Universal Genome WalkerTMKit（Clontech）使用說明書構(gòu)建［17]。根據(jù)HbNTL1轉(zhuǎn)錄因子的cDNA序列，設計一對巢式PCR特異性引物：5'-CTCTGTGGCCCGGTTTGTCTTGGAACCA-3'（NTLGSP1）和5'-GTGAAGCCGCAGAATCACCACGTGCCA-3'（NTLGSP2）。上游引物為Universal Genome WalkerTMKit自帶的接頭引物AP1和AP2，以PvuII，EcoRV，StuI和DraI酶切連接的產(chǎn)物為模板，反應程序為：94℃ 25 s，72℃ 3 min，7個循環(huán)；94℃ 25 s，67℃3 min，32個循環(huán)。

1.2.3 生物信息學分析 利用Softberry（http：//linu x1.softberry.com/berry.phtml）網(wǎng)站的TSSP（植物啟動子預測）數(shù)據(jù)庫分析了所克隆得到的可能為巴西橡膠樹HbNTL1轉(zhuǎn)錄因子基因啟動子的序列；利用Berkeley Drosophila Genome Project（BDGP）的在線軟件Neural Network Promoter Prediction Promoter（http：//www.fruitfly.org/seq_tools/promoter.html）預測核心啟動子區(qū)域和轉(zhuǎn)錄起始位點；利用在線軟件PlantCARE（http：//bioinformatics.psb.ugent.be/webtools/plantcare/html/）和搜索植物順式作用元件數(shù)據(jù)庫PLACE（http：//www.dna.affrc.go.jp/PLACE/）對啟動子序列中可能存在的順式作用元件進行分析［18]。

2 結(jié)果

2.1 橡膠樹HbNTL1基因啟動子的克隆

根據(jù)橡膠樹HbNTL1基因的cDNA序列設計一對巢式PCR特異性引物NTLGSP1和NTLGSP2，進行PCR擴增，兩輪PCR擴增后得到約2.0 kb的條帶（圖1），將PCR產(chǎn)物連接pMD18-T vector，轉(zhuǎn)化大腸桿菌DH5α，篩選陽性克隆進行測序。將測序結(jié)果與已知序列比對分析，去除載體、接頭引物和與讀碼框重疊序列，得到橡膠樹HbNTL1基因上游非編碼區(qū)1 718 bp的序列。

2.2 HbNTL1基因啟動子的生物信息學分析

在線軟件TSSP預測結(jié)果表明，HbNTL1基因5'非編碼區(qū)序列存在一個可能的啟動子區(qū)，位置在1 513 bp處，一個對應的TATA框，位置在1 484 bp處，一個增強子CAAT-box，在1 436 bp處（圖2）。BDGP在線軟件Neural Network Promoter Prediction Promoter預測HbNTL1基因的核心啟動子序列和轉(zhuǎn)錄起始位點，結(jié)果（圖2）顯示，HbNTL1基因5'非編碼區(qū)序列中含有真核生物典型的核心啟動子區(qū)域（1 473-1 522 bp），其中+1位置的堿基A是轉(zhuǎn)錄起始位點，位于翻譯起始位點（ATG）上游206 bp處。HbNTL1基因的轉(zhuǎn)錄起始位點堿基為A，TATA-box在轉(zhuǎn)錄起始位點的-30 bp處，符合高等植物啟動子的基本結(jié)構(gòu)特征。

利用在線軟件PlantCARE和PLACE對HbNTL1基因啟動子序列中可能含有的順式作用元件進行分析，結(jié)果發(fā)現(xiàn)，除具有12個拷貝的TATA-box和CAAT-box外，還具有9個赤霉素、脫落酸和茉莉酸等激素響應元件（ABRE、GAREAT、 P-box、TATC-box和TGACG-motif）、56個逆境脅迫響應元件（ACGTATERD1、ARE、CCAATBOX1、CURECORECR、DOFCOREZM、HSE、MYBCORE、TC-rich repeats、W-box）、6個光順式反應元件（GAG-motif、GATA-box和G-box）、28個組織特異性表達元件（CACTFTPPCA1、GTGANTG10、POLLEN1LELAT52、ROOTMOTIFTAPOX1和RY-repeat）以及2個拷貝的晝夜節(jié)律控制元件（Circadian）等（表1）。

表1 HbNTL1啟動子區(qū)域主要順式作用元件分析

3 討論

NAC轉(zhuǎn)錄因子是植物特有的一類轉(zhuǎn)錄調(diào)控因子，在植物的生長發(fā)育、器官建成、激素調(diào)節(jié)以及生物和非生物脅迫應答等方面具有重要作用［1，2]。但是對于NAC類轉(zhuǎn)錄因子的研究還主要集中在基因的克隆與轉(zhuǎn)基因表達分析上，而對他們的上游調(diào)控因子和下游靶標基因都不是很清楚。本研究在已知巴西橡膠樹膜結(jié)合NAC轉(zhuǎn)錄因子HbNTL1基因cDNA序列的基礎上，利用Genome Walker的方法獲得了其5'調(diào)控序列1 718 bp，通過在線軟件的預測與分析該序列具有真核生物典型的核心啟動子區(qū)域和大量的順式作用元件。HbNTL1基因啟動子序列中存在多個激素響應元件、光作用順式元件和組織特異性表達元件，因此，HbNTL1的表達可能受到光信號和激素的調(diào)控，并且存在組織特異性。

HbNTL1基因啟動子序列中含有大量的逆境脅迫響應元件，占到了總數(shù)的50%左右。MYBCORE和MYCCONSENSUSAT元件在HbNTL1基因啟動子序列中共發(fā)現(xiàn)有15個拷貝，這些元件分別是MYB和MYC轉(zhuǎn)錄因子的識別位點，與其相應轉(zhuǎn)錄因子在ABA信號通路和非生物脅迫應答中具有重要作用［19，20]。ABRE元件是脫落酸響應元件，在干旱、脫水、鹽和低溫等非生物脅迫時誘導植物體內(nèi)積累脫落酸，從而調(diào)控相關(guān)基因的表達［21]。啟動子分析發(fā)現(xiàn)HbNTL1基因啟動子序列中包含了14個拷貝的DOFCOREZM元件，該元件是Dof 轉(zhuǎn)錄因子的識別位點，調(diào)控保衛(wèi)細胞特異基因的表達，在植物防御過程中具有重要作用［22]。在HbNTL1基因啟動子中還發(fā)現(xiàn)有14個拷貝的W-box元件，該順式作用元件是WRKY轉(zhuǎn)錄因子的識別位點。植物特有的WRKY基因通過與順式作用元件W-box特異結(jié)合調(diào)控下游目標基因表達，與植物的防御反應、逆境應答、衰老信號和生長發(fā)育等密切相關(guān)［23，24]。ABRE、DOFCOREZM、MYBCORE、MYCCONSENSUSATHSE和W-box等反應元件在其他植物脅迫逆境相關(guān)NAC轉(zhuǎn)錄因子基因的啟動子區(qū)域都有報導［1，2，10，25-28]。HbNTL1基因啟動子序列中這些反應元件的大量存在表明，HbNTL1可能是一個逆境脅迫相關(guān)NAC轉(zhuǎn)錄因子（stress-responsive NAC），在橡膠樹抵抗逆境脅迫的生理過程中具有重要功能。本實驗室將利用巴西橡膠樹HbNTL1基因的5'調(diào)控序列，通過酵母單雜交系統(tǒng)篩選調(diào)控HbNTL1轉(zhuǎn)錄因子的上游調(diào)控因子，從而進一步研究HbNTL1的表達調(diào)控及其在橡膠樹逆境脅迫中的作用。

4 結(jié)論

本研究利用Genome Walker 的方法克隆了巴西橡膠樹膜結(jié)合NAC轉(zhuǎn)錄因子HbNTL15'調(diào)控序列1 718 bp，軟件預測與分析該序列具有真核生物典型的核心啟動子區(qū)域和大量的順式作用元件，其中逆境脅迫相關(guān)元件占到50%，HbNTL1可能是一個逆境脅迫相關(guān)NAC轉(zhuǎn)錄因子，在橡膠樹抵抗逆境脅迫的生理過程中具有重要功能。

［1] Puranik S, Sahu PP, Srivastava PS, et al. NAC proteins：regulation and role in stress tolerance［J] . Trends Plant Sci, 2012, 176（6）：369-381.

［2] Nakashima K, Takasaki H, Mizoi J, et al. NAC transcription factors in plant abiotic stress responses［J] . BBA-Gene Regul Mech,2012, 1819（2）：97-103.

［3] Nuruzzaman M, Manimekalai R, Sharoni AM, et al. Genome-wide analysis of NAC transcription factor family in rice［J] . Gene, 2010,465（1-2）：30-44.

［4] Le DT, Nishiyama R, Watanabe Y, et al. Genome-wide survey and expression analysis of the plant-specific NAC transcription factor family in soybean during development and dehydration stress［J] .DNA Res, 2011, 18（4）：263-276.

［5] Kim SG, Park CM. Membrane-mediated salt stress signaling in flowering time control［J] . Plant Signal Behav, 2007, 2（6）：517-518.

［6] Seo PJ, Kim SG, Park CM. Membrane-bound transcription factors in plants［J] . Trends Plant Sci, 2008, 13（10）：550-556.

［7] Tran LSP, Nishiyama R, Yamaguchi-Shinozaki K, et al. Potential utilization of NAC transcription factors to enhance abiotic stress tolerance in plants by biotechnological approach［J] . GM Crops,2010, 1（1）：32-39.

［8] Kim SG, Lee SM, Seo PJ, et al. Genome-scale screening and molecular characterization of membrane-bound transcription factors inArabidopsisand rice［J] . Genomics, 2010, 95（1）：56-65.

［9] Shen H, Yin YB, Chen F, et al. A bioinformatic analysis of NAC genes for plant cell wall development in relation to lignocellulosic bioenergy production［J] . Bioenerg Res, 2009, 2：217-232.

［10] Puranik S, Bahadur RP, Srivastava PS, et al. Molecular cloning and characterization of a membrane associated NAC family gene,SiNAC from foxtail millet［Setaria italica（L.）P. Beauv.] ［J] .Mol Biotechnol, 2011, 49：138-150.

［11] Tang YM, Liu MY, Gao SQ, et al. Molecular characterization of novel TaNAC genes in wheat and overexpression of TaNAC2a confers drought tolerance in tobacco［J] . Physiol Plant, 2012,144（3）：210-224.

［12] Hao YJ, Wei W, Song QX, et al. Soybean NAC transcription factors promote abiotic stress tolerance and lateral root formation in transgenic plants［J] .Plant J, 2011, 68（2）：302-313.

［13] Distelfeld A, Pearce SP, Avni R, et al. Divergent functions of orthologous NAC transcription factors in wheat and rice［J] . Plant Mol Biol, 2012, 78：515-524.

［14] Yang SD, Seo PJ, Yoon HK, et al. TheArabidopsisNAC transcription factor VNI2 integrates abscisic acid signals into leaf senescence via the COR/RD genes［J] . Plant Cell, 2011, 23（6）：2155-2168.

［15] Zhong RQ, Lee CH, Ye ZH. Global analysis of direct targets of secondary wall NAC master switches inArabidopsis［J] . Mol Plant, 2010, 3（6）：1087-1103.

［16] An ZW, Wang QT, Hu YS, et al. Co-extraction of high-quality RNA and DNA from rubber tree（Hevea brasiliensis）［J] . Afr J Biotechnol, 2012, 11：9308-9314.

［17] 夏可燦, 康桂娟, 黎瑜, 等.橡膠樹膠乳HbPLDα1基因機器啟動子的克隆與生物信息學分析［J] .熱帶作物學報, 2012, 1（33）：63-69.

［18] Higo K, Ugawa Y, Iwamoto M, et al. Plant cis-acting regulatory DNA elements（PLACE）database：1999 Nucleic［J] . Acids Res, 1999, 27（1）：297-300.

［19] Yamaguchi-Shinozaki K, Shinozaki K. Organization ofcisacting regulatory elements in osmotic-and cold-stress-responsive promoters［J] . Trends Plant Sci, 2005, 10（2）：88-94.

［20] Abe H, Yamaguchi-Shinozaki K, Urao T, et al. Role ofArabidopsisMYC and MYB homologs in drought-and abscisic acid-regulated gene expression［J] . The Plant Cell Online, 1997, 9（10）：1859-1868.

［21] Narusaka Y, Nakashima K, Shinwari ZK, et al. Interaction between twocis-acting elements, ABRE and DRE, in ABA-dependent expression ofArabidopsisrd29A gene in response to dehydration and high-salinity stresses［J] . Plant J, 2003, 34（2）：137-148.

［22] Plesch G, Ehrhardt T, Mueller-Roeber B. Involvement of TAAAG elements suggests a role for Dof transcription factors in guard cellspecific gene expression［J] . Plant J, 2001, 28（4）：455-464.

［23] Eulgem T, Rushton PJ, Robatzek S, et al. The WRKY superfamily of plant transcription factors［J] . Trends Plant Sci, 2000, 5（5）：199-205.

［24] Rushton PJ, Somssich IE, Ringler P, et al. WRKY transcription factors［J] . Trends Plant Sci, 2010, 15（5）：247-258.

［25] Han QQ, Zhang JH, Li HX, et al. Identification and expression pattern of one stress-responsive NAC gene fromSolanum lycopersicum［J].Mol Biol Rep, 2011, 39：1713-1720.

［26] Nakashima K, Tran LSP, Van Nguyen D, et al. Functional analysis of a NAC-type transcription factor OsNAC6 involved in abiotic and biotic stress-responsive gene expression in rice［J] . Plant J, 2007,51（4）：617-630.

［27] Fang YJ, You J, Xie KB, et al. Systematic sequence analysis and identification of tissue-specific or stress-responsive genes of NAC transcription factor family in rice［J] . Mol Genet Genomics, 2008,280：547-563.

［28] Zheng XN, Chen B, Lu GJ, et al. Overexpression of a NAC transcription factor enhances rice drought and salt tolerance［J] . Biochem Biophys Res Commun, 2009, 379（4）：985-989.