李津津，劉凱，王芳，王紅霞，張志華

(1河北農(nóng)業(yè)大學(xué) 山區(qū)研究所，河北省山區(qū)農(nóng)業(yè)工程技術(shù)研究中心，國(guó)家北方山區(qū)農(nóng)業(yè)工程技術(shù)研究中心，河北 保定 071001；2 河北農(nóng)業(yè)大學(xué) 園藝學(xué)院，河北 保定 071001)

低溫脅迫對(duì)植物的地理分布和生長(zhǎng)發(fā)育有著重要影響。當(dāng)植物受到低溫脅迫時(shí)，轉(zhuǎn)錄因子通過(guò)結(jié)合基因啟動(dòng)子上的順式作用元件激活低溫響應(yīng)基因COR(Cold-Related)的表達(dá)[1]。bHLH(Basic/helix-loop-helix)因含有大約60個(gè)氨基酸組成高度保守的bHLH結(jié)構(gòu)域而得名，包含basic和HLH區(qū)域，而且是轉(zhuǎn)錄因子大家族成員之一，廣泛存在于植物中，在植物生長(zhǎng)、脅迫應(yīng)答、信號(hào)轉(zhuǎn)導(dǎo)和次生代謝調(diào)控中有重要作用[2-3]。ICE屬于一種類(lèi)似MYC的bHLH轉(zhuǎn)錄因子，可與CBF3啟動(dòng)子中MYC作用元件特異性結(jié)合，并誘導(dǎo)CBF/DREB1下游基因的轉(zhuǎn)錄表達(dá)[4]。目前，在模式植物擬南芥中ICE1(Inducer of CBF Expression1)-CBF(C-Repeat-Binding-Factor)-COR(Cold-Related)低溫信號(hào)途徑已經(jīng)基本清晰，為其他植物研究抗寒分子機(jī)制的發(fā)展提供了理論依據(jù)[5]。自2003年Chinnudamy等[6]首次發(fā)現(xiàn)擬南芥突變體ice1中CBF3及CBFs下游基因的表達(dá)出現(xiàn)了下降，而且抗寒能力下降。研究人員陸續(xù)在黑枸杞[7]、棗[8]桑樹(shù)[9]、桃[10]、苦蕎[11]鑒定出了bHLH基因家族轉(zhuǎn)錄因子，并且在忍冬[12]、獼猴桃[13]、油菜[14]、茶樹(shù)[15]、赤桉[4]等植物中鑒定了ICE1同源基因，但未發(fā)現(xiàn)在核桃中有報(bào)道。ICE是ICE-CBF-COR冷響途徑最上游調(diào)控因子[5]。在響應(yīng)低溫應(yīng)答過(guò)程中，ICE1轉(zhuǎn)錄因子被激活，其下游基因CBFs及COR被誘導(dǎo)表達(dá)，以提高自身抵御低溫脅迫的能力[16]。已有研究表明，將油菜ICE1基因轉(zhuǎn)化煙草后，在低溫脅迫下轉(zhuǎn)基因煙草明顯比野生型對(duì)照有較高的承受能力[14]。過(guò)表達(dá)ICE1基因水稻在4 ℃低溫脅迫后，轉(zhuǎn)基因株系存活率和脯氨酸含量明顯增加，相對(duì)電導(dǎo)率和丙二醛含量積累速率明顯下降，提高了抗低溫脅迫能力[17]。已有多項(xiàng)研究證實(shí)ICE1基因通過(guò)不同作用方式參與植物自身提高耐寒性能力的過(guò)程[18，21]。因此，抗寒機(jī)制研究對(duì)促進(jìn)培育抗寒核桃新品種有重要意義。

核桃(JuglansregiaL.)，為胡桃科胡桃屬的落葉喬木，分布和栽培遍及世界多個(gè)國(guó)家和地區(qū)，是重要的木材樹(shù)種、生態(tài)樹(shù)種、經(jīng)濟(jì)樹(shù)種和生物能源樹(shù)種。核桃是喜溫樹(shù)種，適宜生長(zhǎng)的溫度范圍是年均溫 9～16 ℃，極端最低氣溫-25 ℃[22]。由于核桃自身特性，低溫脅迫對(duì)核桃的生長(zhǎng)發(fā)育是一個(gè)限制因子。特別是晚霜，使核桃的組織器官受凍，影響產(chǎn)量，甚至使核桃絕收。目前為止，關(guān)于核桃的研究主要集中在栽培、物理防寒等方面，用分子水平來(lái)提高核桃抗寒性還未見(jiàn)報(bào)道。本研究通過(guò)核桃低溫脅迫下高通量轉(zhuǎn)錄組測(cè)序及RT-PCR技術(shù)分析并克隆得到核桃中ICE同源基因，并對(duì)其進(jìn)行了理化性質(zhì)、蛋白性質(zhì)、種間同源基因比對(duì)及系統(tǒng)進(jìn)化等生物信息學(xué)分析及轉(zhuǎn)錄表達(dá)分析，旨在為培育核桃耐低溫新品種，為研究ICE基因在核桃中的低溫響應(yīng)機(jī)制奠定基礎(chǔ)。

1 材料與方法

1.1 試驗(yàn)材料

供試材料是種植于河北農(nóng)業(yè)大學(xué)標(biāo)本園的2 a生盆栽核桃(JuglansregiaL.)，為同一株砧木上嫁接“遼寧8號(hào)”(耐低溫)和“清香”(不耐低溫)[23]的共計(jì)20株嫁接苗。試驗(yàn)材料分為正常組(冷處理0 h)和0.5 ℃(冷處理48 h)，正常組和冷處理組各處理10株嫁接苗。低溫處理48 h后采集每株第1片復(fù)葉上第2片單葉送至北京諾禾致源科技股份有限公司用于RNA-Seq測(cè)序，3次重復(fù)。其中，研究JrICE1在不同組織表達(dá)的供試材料為河北農(nóng)業(yè)大學(xué)標(biāo)本園“清香”核桃，分別采集其根、嫩莖(4月26日)、嫩葉(4月26日)、雄花、雌花、老葉(7月15日)組織，采集后液氮速凍，送至北京諾禾致源科技股份有限公司用于RNA-Seq測(cè)序，3次重復(fù)。

1.2 核桃ICE1基因的分子克隆

該研究通過(guò)轉(zhuǎn)錄組測(cè)序[padj<0.05 and log2(golden change>0)]及利用RT-PCR技術(shù)結(jié)合分析并克隆得到核桃中1個(gè)ICE1同源基因。核桃JrICE1全長(zhǎng)克隆引物序列,見(jiàn)表1。

由表1可知，找到其全長(zhǎng)序列并設(shè)計(jì)引物，以核桃葉片cDNA為模板擴(kuò)增出基因的編碼序列。PCR擴(kuò)增產(chǎn)物經(jīng)1.2 %瓊脂糖凝膠電泳檢測(cè)，將符合預(yù)期片段大小的片段回收連接pMD18-T載體上由上海生工公司完成DNA測(cè)序。

1.3 核桃ICE基因的生物信息學(xué)分析

采用NCBI中Blast進(jìn)行序列的查找與分析，利用ORF Finder (https://www.ncbi.nlm.nih.gov /orffnder/)查找ICE基因的ORF并完成核苷酸的翻譯。利用NCBI中的Conserved Domain(https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi)預(yù)測(cè)該基因的保守區(qū)域；利用在線(xiàn)分析工具Gene Structure Display Server預(yù)測(cè)基因結(jié)構(gòu)；利用蛋白質(zhì)的理化性質(zhì)使用ExPASy (https://web.expasy.org/ protparam/)在線(xiàn)網(wǎng)站分析；利用PSIPRED(http://bioinf.cs. ucl.ac.uk /psipred)分析蛋白的二級(jí)結(jié)構(gòu)，利用NetPhos 3.1 Server和Net NGlyc 1.0 Server分別預(yù)測(cè)蛋白的糖基化和磷酸化位點(diǎn)；利用 SignalP 3.0 Server(http://www.cbs.dtu.dk/services/SignalP)在線(xiàn)工具對(duì)信號(hào)肽預(yù)測(cè)；利用TMHMM 法分析蛋白質(zhì)的跨膜區(qū)：基于 HMM 方法(http://www.cbs.dtu.dk/services/TMHMM/)的蛋白質(zhì)跨膜區(qū)預(yù)測(cè)工具；利用ngLOC與Nuc-PLoc在線(xiàn)軟件預(yù)測(cè)亞細(xì)胞定位。利用DNAMAN8.0對(duì)其氨基酸進(jìn)行多序列同源比對(duì)；利用MEGA 7.0構(gòu)建ML系統(tǒng)發(fā)育樹(shù)。

2 結(jié)果與分析

2.1 核桃JrICE1基因克隆

根據(jù)設(shè)計(jì)的引物，以核桃葉片總RNA反轉(zhuǎn)錄的cDNA為模板，進(jìn)行PCR擴(kuò)增后，產(chǎn)物經(jīng)1.2 %瓊脂糖凝膠電泳檢測(cè)顯示1條1 000～2 000 bp 的特異性條帶(圖1)，與預(yù)期一致。測(cè)序結(jié)果見(jiàn)圖2，經(jīng) NCBI 數(shù)據(jù)庫(kù)檢索分析(表2)，比對(duì)得到1條全長(zhǎng)基因序列。經(jīng)NCBI的ORF進(jìn)行開(kāi)放閱讀框分析1 659 bp 全長(zhǎng)序列，編碼552個(gè)氨基酸，命名為JrICE1。

表2 核桃JrICE1基本信息Table 2 Basic information of the JrICE1 in walnut

圖1 核桃JrICE1基因的PCR擴(kuò)增Figure 1 PCR amplification of the JrICE1 gene in walnut注: 泳道1是DL 2 000 DNA Marker; 泳道2是JrICE1 全長(zhǎng)。

將上述回收基因片段送至生工生物工程股份有限公司測(cè)序，測(cè)序結(jié)果見(jiàn)圖2。

2.2 核桃JrICE1的生物信息分析

利用ExPasy對(duì)基因的等電點(diǎn)(pI)、蛋白大小(kDa)等蛋白理化性質(zhì)進(jìn)行預(yù)測(cè)，見(jiàn)表3。

注: *終止密碼子。圖2 核桃JrICE1基因的 cDNA序列及其預(yù)測(cè)的氨基酸序列Figure 2 The cDNA sequence and its predicted amino acid sequence of walnut JrICE1 gene

表3 核桃JrICE1蛋白質(zhì)理化性質(zhì)和一級(jí)結(jié)構(gòu)分析Table 3 Analysis of physicochemical properties and primary structure of protein of the walnut JrICE1

由表3可知，JrICE1等電點(diǎn)為5.73，偏酸性；蛋白的分子量為60.51 kDa；不穩(wěn)定系數(shù)為57.26，為不穩(wěn)定蛋白；脂肪系數(shù)為72.30，親水系數(shù)為-0.530，為親水性蛋白。

根據(jù)JrICE1測(cè)序獲得實(shí)際序列，利用NCBI在線(xiàn)分析軟件CD-SEARCH對(duì)核桃JrICE1蛋白質(zhì)結(jié)構(gòu)域預(yù)測(cè)分析見(jiàn)圖3。

圖3 核桃JrICE1保守結(jié)構(gòu)域及位置Figure 3 Conserved domains and locations of the JrICE1

由圖3可知，發(fā)現(xiàn)該蛋白均包含bHLH結(jié)構(gòu)，并且包含3個(gè)功能位點(diǎn)：DNA結(jié)合區(qū)域、E-box的特異性結(jié)合位點(diǎn)以及ACT結(jié)構(gòu)位點(diǎn)，可以看出符合bHLH家族類(lèi)型特點(diǎn)。此外，該蛋白還有1個(gè)與磷酸絲氨酸磷酸化酶相關(guān)基因結(jié)合的位點(diǎn)。

利用Net NGlyc 1.01.0預(yù)測(cè)蛋白的糖基化見(jiàn)圖4。

由圖4可知，發(fā)現(xiàn)JrICE1蛋白有2個(gè)糖基化位點(diǎn)(位置分別為：131和335)；在JrICE1蛋白發(fā)現(xiàn)60個(gè)磷酸化位點(diǎn)(絲氨酸45個(gè)、蘇氨酸12個(gè)、酪氨酸3個(gè))；該蛋白存在多個(gè)蛋白激酶結(jié)合位點(diǎn)(細(xì)胞周期蛋白依賴(lài)性激酶I(CKI)、細(xì)胞周期蛋白依賴(lài)性激酶II(CKII)、胰島素受體(INSR)、蛋白激酶A (PKA)、蛋白激酶C (PKC)、蛋白激酶G (PKG)、ATM激酶、細(xì)胞周期素依賴(lài)蛋白激酶5 (CDK5)、細(xì)胞分裂周期激酶2 (CDC2)、DNA依賴(lài)性蛋白激酶(DNAPK)、促分裂素原活化蛋白激酶(6p38 MAPK)、糖原合成酶激酶3 (GSK3)等12個(gè)。

利用SignalP 3.0 Server在線(xiàn)工具對(duì)JrICE1蛋白的信號(hào)肽進(jìn)行分析，發(fā)現(xiàn)JrICE1蛋白沒(méi)有信號(hào)肽(圖4-c)，推測(cè)該蛋白為非分泌蛋白；基于HMM方法的蛋白質(zhì)跨膜區(qū)預(yù)測(cè)(參見(jiàn)圖4-d)，預(yù)測(cè)JrICE1蛋白均不存在跨膜結(jié)構(gòu)域；利用蛋白質(zhì)二級(jí)結(jié)構(gòu)預(yù)測(cè)在線(xiàn)工具PSIPRED對(duì)JrICE1蛋白二級(jí)結(jié)構(gòu)預(yù)測(cè)(圖4-e)，發(fā)現(xiàn)JrICE1有9處螺旋結(jié)構(gòu)。利用Nuc-PLoc在線(xiàn)軟件分析得出，JrICE1蛋白位于細(xì)胞核。分析JrICE1蛋白的三級(jí)結(jié)構(gòu)，結(jié)果顯示存在螺旋、折疊等結(jié)構(gòu)(圖4-e)。

2.3 核桃JrICE1蛋白系統(tǒng)發(fā)育分析

利用DNAMAN 8軟件對(duì)核桃JrICE1蛋白序列與其他15種植物進(jìn)行比對(duì)，見(jiàn)圖5。

如圖5可知，核桃和其他植物一樣具有bHLH基因家族的典型功能區(qū)域，包括富S區(qū)(A)、NLS區(qū)(B)、bHLH結(jié)構(gòu)域(B+C)，下游含轉(zhuǎn)膜區(qū)(E)，且在bHLH結(jié)構(gòu)域高度保守，其他區(qū)域還有部分變異，但是僅有JrICE1含有類(lèi)泛素化(SUMO)結(jié)合位點(diǎn)(D)。

利用MEGA7.0軟件，建立16種植物ICE蛋白序列的系統(tǒng)發(fā)育樹(shù)，并通過(guò)MEME在線(xiàn)分析軟件分析15種植物的保守基序，結(jié)果見(jiàn)圖6。

由圖6可知，核桃JrICE1與毛果楊(Populustrichocarpa)、甜楊(Populussuaveolens)關(guān)系最近。

注: PtICE1(毛果楊，Populus trichocarpa，ABN58427.1)；HvICE2(大麥，Hordeum vulgare，ABA25897.1)；AtiICE1(節(jié)節(jié)麥，Aegilops tauschii，EMT15292.1)；AtICE1(擬南芥，Arabidopsis thaliana，AAP14668.1)；TaICE87(小麥，Triticum aestivum，ACB69502.1)，TuICE1(烏拉圖小麥，Triticum urartu，EMS66007.1)；OSICE1(水稻，Oryza sativa，NP_001045272.1)；BrsICE1(大白菜，Brassica rapa subsp. Chinensis，ACB70963.1)ZmICE2(玉米，Zea mays，ACG46593.1)；GmICE1(大豆，Glycine max，NP_001238560.1)；SlICE1(番茄，Solanum lycopersicum，AGG38826.1)；MdbHLH1(蘋(píng)果，Malus domestica, ABS50251.1)。

注: A是系統(tǒng)進(jìn)化樹(shù)和基序；B是10個(gè)保守基序。

2.4 基于RNA-Seq核桃JrICE1基因轉(zhuǎn)錄表達(dá)分析

基于RNA-Seq核桃JrICE1基因轉(zhuǎn)錄表達(dá)分析，見(jiàn)圖7。

注: QCK是清香對(duì)照；Q48是清香48 h低溫；LCK是遼寧8號(hào)對(duì)照；L48是遼寧8號(hào)48 h低溫。

由圖7可知，利用轉(zhuǎn)錄組分析JrICE1在“遼寧8號(hào)”(耐低溫)和“清香”(不耐低溫)低溫脅迫48 h后中的表達(dá)模式表明，JrICE1在“遼寧8號(hào)”葉片呈現(xiàn)顯著上升，比對(duì)照上升0.53倍；JrICE1在“清香”葉片呈現(xiàn)顯著上升，比對(duì)照上升4.5倍。

利用轉(zhuǎn)錄組分析JrICE1在“清香”核桃不同組織中的表達(dá)情況，結(jié)果見(jiàn)圖8。

注: Y是芽；NY是嫩葉；NJ是嫩莖；XH是雄花；CH是雌花；G是根；LY是老葉。

由圖8可知，JrICE1在“清香”核桃的芽、嫩葉、嫩莖、雄花、雌花、根、老葉中都有表達(dá)，其中表達(dá)量最高的是雌花，表達(dá)量最低的是雄花，其次是老葉、芽、嫩莖、嫩葉和根，暗示了JrICE1可能在雄花器官中起到更加重要的作用。

3 討論

植物在生長(zhǎng)的過(guò)程中，往往會(huì)受到生物脅迫或者非生物脅迫，嚴(yán)重影響植物的生長(zhǎng)發(fā)育，植物需要自身抗逆的調(diào)節(jié)來(lái)應(yīng)對(duì)脅迫維持自身正常生長(zhǎng)。有研究表明，ICE1基因可以調(diào)控冷反應(yīng)基因，使植物抵御寒冷的逆境環(huán)境[6]。本研究通過(guò)同源克隆與RNA-Seq結(jié)合的方式，在核桃葉片中分離出ICE基因，通過(guò)進(jìn)化分析表明核桃JrICE1與毛果楊(Populustrichocarpa)、甜楊(Populussuaveolens)關(guān)系最近。JrICE1在蛋白序列比對(duì)表明含有GAQPTLFQKRA這11個(gè)保守的氨基酸序列(圖5)，這與Mohamed[24]等人研究一致。JrICE1具有1 659 bp的完整開(kāi)放閱讀框，且蛋白是一種偏酸性的不穩(wěn)定親水蛋白，亞細(xì)胞定位于細(xì)胞核，這和已研究的其他植物相似[12-13,24]。利用NCBI中的Conserved Domian-Search預(yù)測(cè)保守結(jié)構(gòu)域，發(fā)現(xiàn)該基因包含bHLH結(jié)構(gòu)，并且包含3個(gè)功能位點(diǎn)：DNA結(jié)合區(qū)域、E-box的特異性結(jié)合位點(diǎn)以及ACT結(jié)構(gòu)位點(diǎn)，結(jié)果表明符合bHLH家族類(lèi)型特點(diǎn)，并屬于核桃ICE基因，證明了基因克隆的準(zhǔn)確性。核桃JrICE1與已報(bào)道其他典型植物的ICE序列比對(duì)，發(fā)現(xiàn)JrICE1上游含有富含二硫鍵的S區(qū)(圖5-A)，主要是維持基因的穩(wěn)定性[25]。通過(guò)NLS區(qū)(圖5-B)和轉(zhuǎn)膜區(qū)(圖5-E)保證蛋白轉(zhuǎn)運(yùn)到細(xì)胞核，而且JrICE1含有一個(gè)類(lèi)泛素化(SUMO)結(jié)合位點(diǎn)，這與已報(bào)道的甜楊[3]、忍冬[12]、火龍果[24]結(jié)構(gòu)相似，研究表明符合核桃ICE1作為轉(zhuǎn)錄因子特性。

已有研究報(bào)道，在擬南芥[6]、煙草[26]、水稻[27]、番茄[28]中證實(shí)了ICE可以正調(diào)控植物的抗寒性。在本研究中，將“遼寧8號(hào)”(耐低溫)和“清香”(不耐低溫)低溫(0.5 ℃)脅迫48 h后通過(guò)葉片轉(zhuǎn)錄組測(cè)序表明，JrICE1表達(dá)呈現(xiàn)顯著上升，說(shuō)明JrICE1在核桃中主要起正調(diào)控作用；不同組織轉(zhuǎn)錄表達(dá)分析表明，JrICE1在葉與花器官中轉(zhuǎn)錄表達(dá)較高，而且在不同組織中均有表達(dá)，說(shuō)明JrICE1在核桃生長(zhǎng)發(fā)育中具有重要功能，但具體的調(diào)控機(jī)制還有待于進(jìn)一步研究。

4 結(jié)論

本研究利用同源克隆與轉(zhuǎn)錄組測(cè)序技術(shù)結(jié)合的方法，克隆了核桃中ICE同源基因，并進(jìn)一步對(duì)其進(jìn)行生物信息分析及轉(zhuǎn)錄表達(dá)分析。結(jié)果表明：JrICE1的 CDS序列全長(zhǎng)為1 659 bp，最大開(kāi)放閱讀框1 659 bp，編碼552個(gè)氨基酸。JrICE1蛋白為不穩(wěn)定親水性蛋白，理論等電點(diǎn)5.73。同源比對(duì)發(fā)現(xiàn)JrICE1蛋白具有bHLH基因家族的典型功能區(qū)域，且在bHLH結(jié)構(gòu)域高度保守。通過(guò)進(jìn)化分析表明核桃JrICE1與毛果楊(Populustrichocarpa)、甜楊(Populussuaveolens)關(guān)系最近?！斑|寧8號(hào)”(耐低溫)和“清香”(不耐低溫)低溫脅迫48 h后葉片轉(zhuǎn)錄組測(cè)序結(jié)果分析表明，JrICE1在“遼寧8號(hào)”和“清香”葉片中表達(dá)均顯著上升，說(shuō)明JrICE1在核桃響應(yīng)低溫脅迫的過(guò)程中主要發(fā)揮正調(diào)控作用。不同組織轉(zhuǎn)錄表達(dá)分析表明，JrICE1在不同組織中均有表達(dá)，在葉與花器官中轉(zhuǎn)錄表達(dá)水平較高，說(shuō)明JrICE1在核桃生長(zhǎng)發(fā)育過(guò)程具有重要功能。以上研究說(shuō)明ICE基因可能在核桃的生長(zhǎng)發(fā)育特別是應(yīng)對(duì)低溫脅迫中起著重要作用。本研究為揭示JrICE1在核桃低溫脅迫響應(yīng)過(guò)程中的功能提供參考，為解析核桃低溫響應(yīng)機(jī)制，培育核桃耐低溫品種奠定基礎(chǔ)。