姜 騁 張 曦 田 晴 李 莉
(林木遺傳育種國家重點實驗室,東北林業大學,哈爾濱 150040)
堿性螺旋—環—螺旋蛋白(Basic helix-loop-helix protein,bHLH)是一類含有眾多成員的轉錄因子家族,廣泛存在于動物、植物以及微生物中[1]。bHLH基因家族因包含有bHLH結構域而得名,該結構域由60個左右的保守氨基酸序列組成,包含兩個保守但功能不同的結構域,一個是位于N端的堿性氨基酸區,其作用是識別并特異性結合靶基因啟動子中的DNA基序,另一個是位于C端的螺旋—環—螺旋區,其主要功能是進行域二聚,促進蛋白—蛋白相互作用,并允許形成同型二聚體或異型二聚體來控制基因轉錄。bHLH蛋白結構中除了高度保守的bHLH結構域外,其他序列的保守性較差[2]。
近年來的研究發現,植物bHLH基因參與了廣泛的生物學過程中,包括光信號調節[3];激素信號轉導[4];響應高鹽、干旱、機械損傷、氧化脅迫[5]、低溫脅迫[6]、缺鐵反應[7]和類黃酮類物質合成[8]等。目前通過大量的研究已經發現了許多響應逆境脅迫的bHLH轉錄因子。例如Liu等研究發現,過表達bHLH122增強了轉基因擬南芥對干旱、高鹽以及滲透壓的抵抗力,同時bHLH122可以直接調控參與抗逆響應的基因的轉錄表達,表明bHLH122是一個干旱及其他逆境信號通路的轉錄激活子[9]。水稻OsbHLH148和RERJ1(OsbHLH006)基因能夠通過茉莉酸信號途徑參與干旱脅迫和損傷反應的應答[10~11]。小麥的TabHLH39在擬南芥中過表達后,增強了轉基因植株對鹽害、干旱及氧化作用的耐受性[12]。Jin等報道,秋子梨PubHLH1基因能夠增強轉基因植株對冷害的抵抗性,并且激活了逆境應答基因的轉錄表達[6]。
白樺(BetulaplatyphyllaSuk.)是我國東北地區具有重要應用價值的樹種之一。本課題組前期對高鹽、干旱脅迫下白樺基因的轉錄組進行了分析,獲得了逆境脅迫誘導下基因的表達信息,根據高鹽、干旱不同脅迫時間下基因表達量比較分析的結果,選出一個差異表達顯著的bHLH基因,命名為BpbHLH112。為了進一步研究該基因調控白樺應答逆境脅迫響應的功能,本研究克隆了BpbHLH112的編碼區序列,對其序列進行了生物信息學分析,并對鹽和干旱脅迫下BpbHLH112的表達模式進行了分析,進一步克隆獲得BpbHLH112基因啟動子序列。順式作用元件預測分析結果表明,該區域包含多個逆境響應和激素響應等順式作用元件,因而構建了BpbHLH112啟動子片段融合GUS報告基因的表達載體,采用煙草中瞬時表達的方法對該啟動子在正常條件和脅迫條件下的GUS基因表達的特性進行了比較分析,本研究結果為深入研究BpbHLH112基因調控白樺應對非生物脅迫作用的分子機理提供了依據,具有重要的理論意義。
白樺無性系組培苗由本課題組保存。將組培生根苗移栽到黑土∶蛭石∶珍珠巖=5∶3∶2的基質中培養8周,選取株高和長勢一致的幼苗,分別進行200 mmol·L-1NaCl、18%(w/v)PEG6000澆灌處理,處理時間為0、6、12、24、48、72 h,每個時間點均設3個生物學重復。分別剪取處理后植株的根組織,液氮速凍后保存于-80℃,用于總RNA的提取。
參照“植物RNA提取試劑盒”(BioTeke)的操作流程,提取白樺葉片和根組織的總RNA,采用PrimeScriptTMRT reagent Kit with gDNA Eraser(TakaRa)反轉錄試劑盒將總RNA反轉錄成cDNA。在BpbHLH112基因開放閱讀框兩端設計基因上游引物(5′-ATGGCAGAGGAATTTCAGGC-3′)和下游引物(5′-CCTGAATGTTCCCCCAAATG-3′),以白樺cDNA為模板,擴增白樺BpbHLH112基因全長編碼序列(ORF),使用DNA凝膠回收試劑盒、PCR產物,然后克隆到pMD18-T載體上,并轉入到大腸桿菌DH5α感受態細胞中,挑取單克隆擴大培養后進行測序驗證。
根據“新型快速植物基因組DNA提取試劑盒”(BioTeke)說明書,以白樺無性系葉片為材料提取基因組DNA,根據序列比對分析結果,在BpbHLH112基因啟動子區的兩端設計上游引物bHPro-F(5′-CGAGCTCCTGATGGCGATTTTTGTGACT-3′)和下游引物bHPro-R(5′-AGATCTATGTAGCCATCCACCAGTTT-3′),以白樺DNA為模板,擴增BpbHLH112基因啟動子片段,PCR產物膠回收純化后與克隆載體pMD18-T連接,命名為pMD18-T-BpbHLH112pro,連接產物轉化大腸桿菌DH5α中,PCR驗證陽性克隆進行測序。
利用在線軟件ExPASy的ProtParam tool(http://web.expasy.org/cgi-bin/protparam/)程序分析BpbHLH112編碼蛋白的相對分子量、等電點等理化性質,使用在線網站http://smart.emblheidelberg.de進行保守結構域的預測。通過NCBI中的BLASTP(https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastp)獲取其同源蛋白的氨基酸序列,利用DNAMAN進行多序列比對,采用MEGA5.1軟件中鄰位相連法(Neighbor-joining)構建系統進化樹。
運用在線分析工具PSORT Prediction(http://psort1.hgc.jp/form.html)對BpbHLH112基因進行亞細胞定位預測分析。
使用New PLACE[13](http://www.dna.affrc.go.jp/PLACE/)和Plant CARE[14](http://bioinformatics.psb.ugent.be/webtools/plantcare/html/)數據庫預測BpbHLH112基因啟動子的主要順式作用元件。
采用軟件Primer 5設計實時熒光定量PCR上游引物(5′-CGGTCTCCGATCAGAACTCC-3′)和下游引物(5′-TTGCTCTCACCTCTCCCACT-3′)。實時熒光定量PCR反應體系按照TaKaRa公司的SYBR Premix Ex TaqTMⅡ試劑盒使用說明。PCR反應條件為95℃預變性30 s;95℃變性5 s,60℃退火30 s,40個循環。以白樺Tubulin基因(GenBank number:FG067376)作為內參,其擴增引物為5′-TCAACCGCCTTGTCTCTCTCAGG-3′和5′-TGGCTCGAATGCACTGTTGG-3′,每個樣品重復3次,用2-ΔΔCt法進行基因的相對表達量分析。
根據pCAMBIA1301載體和BpbHLH112基因啟動子序列的融合特性,結合ClonExpress Ⅱ One Step Cloning Kit(Vazyme)的使用說明,設計帶限制性內切酶XbaⅠ和BglⅡ酶切位點的啟動子上下游引物(見表1),其中p1301-GUS-F和p1301-GUS-R為pCAMBIA1301載體序列引物,用于驗證目的片段是否連接到載體上。
表1BpbHLH112啟動子片段與GUS基因融合所用引物序列
Table 1 A list of primers used to construct theBpbHLH112 promoter:GUS expression vector
注:‘—’為線性化載體兩末端同源序列;‘=’為用于載體線性化的酶切位點
Note:‘—’the labeled nucleotides were the homologous sequences to vector;‘=’ nucleotides underlined were used for producingXbaⅠ andBglⅡ restriction site
以攜帶有BpbHLH112啟動子片段的pMD18-T-BpbHLH112pro質粒為模板進行PCR擴增,產物經膠回收和雙酶切驗證,同時以pCAMBIA1301空載體為對象,用XbaⅠ和BglⅡ進行雙酶切,經電泳檢測后回收較大載體片段。將啟動子片段BpbHLH112pro與去除了CaMV35S啟動子的pCAMBIA1301載體片段進行無縫連接,使BpbHLH112啟動子片段代替pCAMBIA1301載體中的CaMV35S啟動子,并與GUS報告基因融合。
農桿菌注射滲透法瞬時轉化煙草過程參考Yang等[15],分別將含有陽性對照、陰性對照和重組載體的質粒轉化大腸桿菌DH5α中,把經PCR驗證和測序比對正確的轉化子導入農桿菌GV3101中,經菌液PCR驗證的單菌落擴大培養,用滲透液[10 mmol·L-1MgCl2,10 mmol·L-1MES,150 mmol·L-1AS,pH=5.7]重懸后,注射到6周齡的煙草植株葉片中,浸染后的植株繼續培養2 d,用直徑為1 cm的葉片打孔器制備煙草葉盤,并進行脅迫處理。
將浸染2 d后的葉盤分為4組,分別用H2O(對照)、200 NaCl mmol·L-1、18% PEG6000及100 mmol·L-1ABA脅迫處理24 h,然后進行GUS染色分析和GUS酶活定量分析。
按照Jefferson[16]等的方法對轉基因煙草葉盤進行GUS組織化學染色分析。侵染后的葉盤洗凈后置于GUS染色液中,37℃保溫過夜,用脫色液脫色,觀察染色結果并拍照記錄。
稱取100 mg脅迫處理后的葉盤加入3倍體積的GUS提取緩沖液[50 mmol·L-1磷酸鈉緩沖液,10 mmol·L-1EDTA,0.1% Triton X-100,0.1%(w/v) SDS和0.1% β-巰基乙醇],研磨成勻漿,離心收集上清。利用TAKARA Bradford protein assay kit測定總蛋白含量,使用Infinie200 PRO多功能酶標儀測定樣品的熒光強度。
利用BpbHLH112基因特異性引物,以白樺根和葉片cDNA為模板,通過PCR擴增獲得BpbHLH112基因開放閱讀框(ORF)擴增片段(見圖1)。根據PCR產物測序結果以及對ORF序列的生物信息學分析表明,白樺BpbHLH112基因ORF序列全長為1 431 bp,編碼1個由476個氨基酸殘基組成的蛋白質(見圖1)。預測分子量為52.05 kD,理論等電點(Theoretical pI)為5.78。該蛋白主要由絲氨酸(Ser)、蘇氨酸(Thr)、谷氨酰胺(Gln)等22種氨基酸構成,其中絲氨酸(Ser)的含量最高,占14.7%。脂溶系數(Aliphatic index)為55.92,總平均親水性平均指數(Grand average of hydropathicity)為-0.593,說明BpbHLH112蛋白是一個親水性蛋白。不穩定系數(Instability index)分析顯示該蛋白不穩定系數為57.80,可能是一個不穩定蛋白。

圖2 白樺BpbHLH112與其他物種同源蛋白多序列比對Fig.2 Amino acid sequence alignment of BpbHLH112 and bHLHs from the other plants

圖3 白樺BpbHLH112與其他物種同源蛋白系統進化樹 分支上的數值表示Bootstrap驗證中基于1 000次重復該節點的可信度Fig.3 Phylogentic tree analysis of BpbHLH112 and bHLHs from the other plants The numerical value on the branch indicates the credibility of the node based on 1 000 repetitions in bootstrap verification
利用NCBI上的BLASTP工具,對BpbHLH112進行氨基酸序列同源性的檢索,找出與其序列一致性相對較高的來自9個物種的bHLH蛋白序列。多重序列比對分析結果表明(見圖2),白樺BpbHLH112轉錄因子與其他bHLH轉錄因子都含有一個高度保守的DNA結構域,對應于白樺bHLH112的349~409位氨基酸殘基處,這一保守結構域由60個氨基酸構成,分為兩個功能不同的區域:basic區域(堿性氨基酸區)含16個氨基酸,HLH區域含44個氨基酸殘基,該區域是具有螺旋環螺旋區域基因家族所特有的bHLH保守結構域。不同植物的bHLH蛋白存在幾段長度不等的序列高度保守區域,可以推測這些序列可能是該蛋白的重要功能位點。白樺bHLH112蛋白與櫟樹(Quercusruber)、歐洲葡萄(Vitisvinifera)、銀白楊(Populusalba)的氨基酸序列一致性較高,分別為75%、64%和63%。
利用MEGA5.1將白樺BpbHLH112氨基酸序列與其他物種bHLH序列進行聚類分析,采用最大似然法繪制系統進化樹,結果表明白樺BpbHLH112與櫟樹bHLH同源蛋白的親緣關系較近(見圖3)。
運用在線分析工具PredictNLS預測BpbHLH112蛋白的亞細胞定位,結果表明,BpbHLH112可能定位在細胞核上,預測可信度為92。

圖4 鹽和干旱脅迫下白樺根中BpbHLH112基因的表達分析Fig.4 Expression analysis of BpbHLH112 under salt and PEG6000 stress
為了研究白樺BpbHLH112基因在鹽和干旱脅迫處理下的表達模式,分別用200 mmol·L-1NaCl或18% PEG6000脅迫處理2個月大的白樺實生苗,在0、6、12、24、48、72 h分別對各處理的根組織取材,提取RNA反轉錄成cDNA作模板進行實時RT-PCR,結果如圖4所示。
經NaCl處理6~48 h的時間內,白樺BpbHLH112基因的表達呈上調趨勢,在處理48 h時達到最高,其表達量高于0 h對照(即未脅迫脅迫)的9倍以上,隨著脅迫時間的延長,在72 h時又回落到6 h的水平。在PEG6000處理6~12 h的時間內,白樺BpbHLH112基因的相對表達量沒有發生顯著變化,在處理24 h后,與0 h對照相比,其表達量呈顯著上升趨勢,在處理72 h時達到最高,為對照的7倍以上。以上分析表明,白樺BpbHLH112基因能夠對鹽或干旱脅迫產生應答反應。
以白樺DNA為模板,采用上游引物bHPro-F和下游引物bHPro-R,擴增BpbHLH112基因啟動子片段,使用1%的瓊脂糖凝膠電泳分別進行檢測,結果顯示其條帶大小與目的片段實際大小相一致(見圖5A)。將啟動子片段與載體pMD18-T的連接片段pMD-18-T-pBpbHLH112pro轉化DH5α感受態細胞中,經菌落PCR驗證的陽性克隆菌液送測序。測序結果顯示,獲得的啟動子片段與白樺全基因組測序結果中該基因啟動子序列的比對完全一致。
使用New PLACE和PlantCARE在線軟件分析預測BpbHLH112啟動子所含有的順式作用元件,分析結果顯示:BpbHLH112基因啟動子區域除了含有TATA-box和CAAT-box等核心元件之外,還包含了一些激素相關元件,如脫落酸應答相關元件ABRELATERD1和ABRERATCAL;赤霉素應答元件GAREAT、MYBGAHV和WRKY71OS;生長素、水楊酸和光響應元件ASF1MOTIFCAMV;水楊酸響應元件WBOXATNPR1等;逆境脅迫響應元件如干旱響應元件MYB1AT和MYBCORE;干旱及低溫響應相關元件MYCCONSENSUSAT;鹽及病原物響應相關元件GT1GMSCAM4;鹽響應、水楊酸及光響應元件GT1CONSENSUS;損傷響應元件WBOXNTERF3;抗病響應相關元件BIHD1OS等,光響應元件如GATABOX、IBOXCORE、INRNTPSADB等。除此之外,BpbHLH112基因啟動子還含有多個參與保衛細胞中K+滲入通道基因表達的TAAAGSTKST1元件、Ca2+響應相關順式元件ABRERATCAL、參與Ca2+/calmodulin信號響應的“CGCG box”、光合作用中碳代謝相關元件DofCOREZM、參與轉錄激活ARR1AT元件、花粉特異表達的順式作用元件POLLEN1LELAT52等其他控制植物生長發育的順式作用元件(見表2),這表明BpbHLH112基因啟動子的活性可能受多種因素的調節。
以pMD-18-T-pBpbHLH112為模板,并以bHPro-F/bHPro-R為引物進行PCR擴增,使用1%的瓊脂糖凝膠電泳進行檢測,結果如圖4A所示,其條帶大小與BpbHLH112基因啟動子片段實際大小相一致,膠回收后連接去除35S啟動子的pCAMBIA1301載體的骨架片段,轉化大腸桿菌DH5α后,提取質粒并用pCAMBIA1301載體引物p1301-GUS-F/1301-GUS-R和啟動子擴增片段引物交叉進行PCR驗證(見圖5B),證明BpbHLH112基因啟動子片段已經連接到pCAMBIA1301載體上,采用液氮凍融法導入農桿菌GV3101感受態細胞中,在含有50 mg·L-1Kana和50 mg·L-1Rif的LB培養基上篩選,經菌落PCR鑒定獲得陽性克隆,可用于后續的試驗。

圖5 BpbHLH112基因啟動子融合GUS載體構建 A.BpbHLH112啟動子克隆(M. DL2000 DNA Marker;1.BpbHLH112基因啟動子全長);B.BpbHLH112啟動子重組質粒檢測(M. DL2000 DNA Marker; 1. BpbHLH112啟動子引物和pCAMBIA1301載體引物交叉PCR產物)Fig.5 Construction of GUS expression vector fused with BpbHLH112 promoter fragment A.Cloning of the promoter of BpbHLH112 gene(M. DL2000 DNA Marker; 1. Full length of BpbHLH112 promoter); B.Detection of recombinant plasmid of BpbHLH112 promoter(M. DL2000 DNA Marker; 1. PCR product using combined primers derived from vector and BpbHLH112 promoter)
利用注射法將攜帶有BpbHLH112啟動子片段重組質粒的農桿菌浸染煙草葉片,陽性對照為帶有CaMV 35S啟動子的pCAMBIA1301載體質粒,而去除CaMV 35S啟動子的pCAMBIA1301載體質粒作為陰性對照。
GUS組織化學染色分析結果如圖6A所示,陰性對照未被染上藍色,而陽性對照被染成了明顯的深藍色。在實驗組中, 轉入BpbHLH112啟動子片段載體質粒的葉盤在H2O對照處理下染色較淺,而經過NaCl、PEG6000和ABA處理后,浸染后的葉片染色效果明顯加深,表明BpbHLH112啟動子片段能夠響應干旱、鹽和ABA脅迫,驅動GUS基因在煙草葉片中的表達。進一步對BpbHLH112啟動子片段的轉基因煙草葉片在脅迫處理后的GUS酶活性進行測定,結果表明(見圖6B),在鹽、干旱和ABA脅迫誘導下,轉入BpbHLH112啟動子片段的葉盤中GUS酶活性均有顯著的變化,分別是H2O對照處理葉盤中GUS酶活性的9.82、12.14和8.09倍。由此可進一步表明,BpbHLH112啟動子受干旱、鹽及ABA脅迫誘導,屬于逆境脅迫誘導型啟動子。

表2 白樺BpBHLH112啟動子的部分順式作用元件
注:堿基序列:R=G/A,W=A/T,N=A/T/C/G
Note:Base sequence:R=G/A,W=A/T,N=A/T/C/G

圖6 BpbHLH112啟動子片段對非生物脅迫響應分析 A.NaCl、PEG6000和ABA脅迫下轉基因煙草葉盤GUS組織化學染色分析;B. NaCl、PEG6000和ABA脅迫下轉基因煙草GUS酶活分析,所有數據為3次重復Fig.6 Analysis of BpbHLH112 promoter construct in tobacco plants under salt,drought and ABA treatments A. GUS histochemical staining of BpbHLH112 promoter construct in tobacco plants under NaCl,PEG6000 and ABA treatments; B. GUS activity of BpbHLH112 promoter construct in tobacco plants under NaCl,PEG6000 and ABA treatments Values represent the means±SD from three repeats
植物在生長發育的過程中會受到各種不利環境的影響,為了在變化多端的環境中更好地生存,植物會在環境發生變化時調控特定抗逆基因的表達[17]。逆境相關基因的表達在轉錄水平上的調控是通過啟動子中逆境誘導相關的順式作用元件與反式作用因子的相互作用來實現的[18~19]。本研究基于鹽、干旱脅迫下白樺轉錄組表達譜分析結果,選出白樺BpbHLH112基因,對其表達模式分析表明,BpbHLH112基因能夠響應鹽和干旱脅迫反應,進一步分析發現,BpbHLH112啟動子能夠被鹽、干旱及ABA脅迫因子激活,驅動下游基因的表達,屬于多逆境誘導型啟動子。
BpbHLH112啟動子活性受干旱誘導,這可能與啟動子中含有3個干旱應答相關順式元件有關。在BpbHLH112啟動子中存在有1個MYB1AT基序、6個MYCCONSENSUSAT基序和2個MYBCORE基序。Abe等人分別于1997[20]和2003年[21]證明,擬南芥的干旱響應基因RD22啟動子中MYB1AT和MYCCONSENSUSAT元件可以與MYB或bHLH轉錄因子相結合,在干旱和ABA脅迫下顯著增強靶基因的表達。Urao等發現,擬南芥MYB轉錄因子可以結合MYBCORE元件,調節脫水誘導基因的表達[22],我們推測BpbHLH112啟動子中的干旱應答相關順式元件對啟動子在干旱處理下活性的升高起重要作用。
BpbHLH112啟動子活性也受NaCl一定程度的誘導。前人研究發現,許多植物基因的啟動子中含有序列保守的GT元件[23],擬南芥調控因子TFIIA-TBP-TATA可以與GT-1元件結合形成復雜的復合體,調控植物基因的轉錄過程[24~25]。Park等報道了GT-1-like轉錄因子可以結合GT-1 box(GT1GMSCAM4)元件,在大豆SCaM-4基因啟動子響應病原菌和NaCl誘導的表達調控中起正調控作用[26]。從序列分析結果可知,BpbHLH112啟動子中存在有5個GT1GMSCAM4元件,它們可能對BpbHLH112啟動子在高鹽脅迫下的表達活性提高起到調控作用。這與早前對豌豆PsSEOF1基因[27]和甜瓜CmLOX08[28]基因啟動子功能研究的結果相一致。
ABA是重要的植物激素之一,不僅參與多個植物生長發育過程,而且作為脅迫激素參與植物對外界脅迫條件的適應[29]。植物轉錄因子參與調節干旱、高鹽、低溫等脅迫信號應答的途徑有兩種:一種是依賴于ABA的信號轉導途徑,另一種是不依賴于ABA的信號轉導途徑[30]。研究表明,轉錄因子通過結合ABA脅迫相關的順式作用元件如ABREs、DRE/CRTs或者含有MYB和MYC蛋白結合位點(MYBR和MYCR)等激活植物對逆境的應激反應[31]。本研究發現,在BpbHLH112啟動子區域分布有1個ABRELATERD1元件、1個ACGTATERD1和3個ABRERATCAL元件。Nakashima等報道,擬南芥脫水誘導基因RD29B和RD29A啟動子中ABA響應元件ABREs與亮氨酸拉鏈類(bZIP)轉錄因子AREBs結合,激活靶基因在特定環境條件下表達[32]。因此,BpbHLH112啟動子中的ABREs元件對該啟動子在ABA處理下的活性升高可能是一種正調控作用,但是BpbHLH112基因是通過依賴或不依賴ABA信號途徑對高鹽和干旱脅迫產生應答反應尚需進一步的研究。
鈣離子(Ca2+)作為第二信使參與各種細胞內以及細胞間信號轉導,以鈣離子及其靶蛋白鈣調素為核心的鈣信使系統在植物應對逆境脅迫的信號轉導過程中起著中心作用[33~34]。Yang報道,6個擬南芥信號應答基因AtSR1-6能夠被多種環境脅迫如高溫、高鹽、UV、損傷以及多種信號分子如茉莉酸甲脂、H2O2和SA處理誘導表達,AtSR1-6基因編碼鈣調素結合蛋白,它們與CGCGBOX基序結合促進抗逆基因的激活與表達[35]。在BpbHLH112啟動子區域存在有4個CGCG box元件,對其可能結合的鈣調素結合蛋白的鑒定及其功能可以進一步深入研究。此外,BpbHLH112啟動子還包含有3個TAAAG元件,該種元件參與了馬鈴薯保衛細胞中K+滲入通道基因KST1的表達[36],因此,對BpbHLH112啟動子TAAAG元件功能的進一步研究可能了解TAAAG元件結合的調控因子以及控制保護細胞特異性基因表達的機制。
本研究初步分析了BpbHLH112基因及其啟動子響應高鹽、干旱和外源ABA脅迫下的表達模式,發現了一些與非生物脅迫以及逆境信號傳導相關的元件,為后續研究該基因啟動子轉錄應答反應機制提供了理論基礎,后續研究可構建啟動子5′缺失突變體和主要調節區域定點突變,從而找到核心啟動子區域和功能作用元件,進一步應用于林木遺傳育種改良研究。