徐學(xué)蕾，鞠延侖，王婉妮，吳今人，房玉林,2*

(1 西北農(nóng)林科技大學(xué) 葡萄酒學(xué)院，陜西楊陵 712100；2 西北農(nóng)林科技大學(xué)，合陽(yáng)葡萄試驗(yàn)示范站，陜西合陽(yáng) 715307)

葡萄(VitisviniferaL.)是世界上具有重要經(jīng)濟(jì)價(jià)值的果樹(shù)，然而，植物生長(zhǎng)的自然環(huán)境是由一系列復(fù)雜的生物脅迫和非生物脅迫構(gòu)成，其中干旱、高溫、低溫、高鹽等非生物脅迫輕則影響植物的生長(zhǎng)發(fā)育，重則導(dǎo)致植物的死亡，最終影響作物的質(zhì)量和產(chǎn)量[1-2]。為避免和減輕這些逆境脅迫對(duì)植物造成的影響，植物可以在分子、細(xì)胞和生理水平上進(jìn)行應(yīng)答[3]。因此，通過(guò)分子生物學(xué)技術(shù)改良植物品質(zhì)，對(duì)提高葡萄產(chǎn)量和質(zhì)量具有重要的意義。

bHLH是廣泛存在于動(dòng)物、植物及真菌中的一大類呈堿性且有螺旋-環(huán)-螺旋(basic-Helix-Loop-Helix，bHLH)結(jié)構(gòu)的轉(zhuǎn)錄因子[4]，總共含有約50～60個(gè)氨基酸[5-6]。其中，堿性氨基酸區(qū)域位于bHLH結(jié)構(gòu)域的N端，含有約15個(gè)氨基酸，包括約6個(gè)堿性殘基，其主要功能是識(shí)別并特異性結(jié)合靶基因啟動(dòng)子中的DNA序列[7-8]。螺旋環(huán)螺旋區(qū)域位于bHLH結(jié)構(gòu)域的C端，包含約40～50個(gè)氨基酸，包括兩個(gè)由疏水環(huán)連接的疏水殘基組成的兩性α螺旋，其主要功能是進(jìn)行結(jié)構(gòu)域二聚化以促進(jìn)蛋白質(zhì)和蛋白質(zhì)相互作用并允許形成同源二聚體或異源二聚體復(fù)合物以調(diào)節(jié)下游靶基因的表達(dá)[9]。

目前，bHLH93的研究集中在擬南芥、蘋(píng)果和甜菜等作物中。Sharma[10]研究發(fā)現(xiàn)，AtbHLH93在體內(nèi)和體外與開(kāi)花抑制劑MAF5啟動(dòng)子元件結(jié)合，從而調(diào)控?cái)M南芥的開(kāi)花時(shí)間，并在控制擬南芥生長(zhǎng)形態(tài)方面發(fā)揮著重要的作用；An等[11]研究表明，蘋(píng)果中AtbHLH93的同源基因MdbHLH93受MdBT2蛋白的泛素化修飾，加速ABA調(diào)節(jié)的葉片衰老，并直接激活衰老結(jié)構(gòu)基因MdSAG18的轉(zhuǎn)錄而負(fù)調(diào)節(jié)葉片的衰老；Wang等[12]在甜菜的研究中發(fā)現(xiàn)，BvbHLH93可以通過(guò)提高抗氧化活性和減少活性氧(ROS)的產(chǎn)生調(diào)節(jié)鹽脅迫耐受性。

bHLH家族是與植物抗逆性相關(guān)的轉(zhuǎn)錄因子家族的第二大家族，在有效應(yīng)對(duì)植物生物脅迫及非生物脅迫，提升作物產(chǎn)量方面起著重要的作用[13-14]。葡萄作為世界上種植最廣泛的水果作物之一，具有重要的經(jīng)濟(jì)價(jià)值[15]。盡管bHLH93在植物生長(zhǎng)發(fā)育中具有重要作用，但在葡萄中尚未見(jiàn)相關(guān)研究報(bào)道。因此，本研究在克隆葡萄VvbHLH93基因的基礎(chǔ)上進(jìn)行了生物信息學(xué)分析，并結(jié)合該基因在葡萄不同組織及果實(shí)不同生長(zhǎng)發(fā)育時(shí)期的表達(dá)模式分析，以期為深入研究bHLH家族基因的功能提供理論基礎(chǔ)。

1 材料和方法

1.1 實(shí)驗(yàn)材料

實(shí)驗(yàn)材料為釀酒葡萄‘赤霞珠’(VitisviniferaL. cv. Cabernet Sauvignon)，采自西北農(nóng)林科技大學(xué)葡萄酒學(xué)院釀酒葡萄示范基地。試驗(yàn)于2021年6～9月采集不同時(shí)期果實(shí)；于7月份采集根、莖、葉、芽及果皮，用于不同組織及葡萄不同生長(zhǎng)發(fā)育時(shí)期的熒光定量分析。所收集材料先使用錫箔紙包好，經(jīng)液氮速凍后，于-80 ℃冰箱保存。

1.2 方 法

1.2.1 基因克隆葡萄樣品的RNA提取采用RNAprep Pure多糖多酚植物總RNA提取試劑盒(天根生化科技有限公司，北京)，使用核酸蛋白儀和電泳分析進(jìn)行質(zhì)量檢測(cè)，選擇HiScript Ⅱ Reverse Transcriptase反轉(zhuǎn)錄試劑盒(南京諾唯贊生物科技股份有限公司，南京)對(duì)檢測(cè)合格的RNA進(jìn)行反轉(zhuǎn)錄成cDNA。采用Primer Premier 5.0軟件設(shè)計(jì)引物(表1)，并進(jìn)行PCR擴(kuò)增，1%瓊脂糖凝膠檢測(cè)PCR產(chǎn)物，然后利用膠回收試劑盒(北京天漠科技開(kāi)發(fā)有限公司，北京)回收基因的片段，并利用核酸蛋白儀檢測(cè)回收產(chǎn)物的純度和濃度。引物序列見(jiàn)表1。

表1 VvbHLH93克隆及表達(dá)分析所使用的引物及序列

1.2.2 基因生物信息學(xué)分析使用NCBI中ORF finder(https://www.ncbi.nlm.nih.gov/orffinder/)對(duì)測(cè)序獲得的cDNA基因進(jìn)行開(kāi)放閱讀框預(yù)測(cè)；利用在線ProtParam軟件(https://web.expasy. org/protparam/)對(duì)氨基酸序列的基本理化性質(zhì)進(jìn)行預(yù)測(cè)和分析；蛋白質(zhì)的跨膜結(jié)構(gòu)利用TMHMM工具(http://www.cbs.dtu.dk/services/TMHMM/)進(jìn)行分析；信號(hào)肽利用SignalP 4.1(http://www.cbs. dtu. dk/services/SignalP-4.1/)進(jìn)行分析；利用在線NetPhos 3.1(https://services.healthtech.dtu.dk/service.php?NetPhos-3.1)進(jìn)行磷酸化位點(diǎn)的預(yù)測(cè)；疏水性區(qū)域預(yù)測(cè)利用ProtScale(https://web.expasy.org/ cgi-bin/protscale/protscale.pl)進(jìn)行；蛋白質(zhì)二級(jí)結(jié)構(gòu)利用SCOPMA工具(https://npsa-prabi.ibcp.fr/ cgi-bin/npsa_automat.pl?page=/NPSA/npsa_sopma.html)進(jìn)行預(yù)測(cè)；利用SWISS-MODEL在線軟件預(yù)測(cè)VvbHLH93蛋白三級(jí)結(jié)構(gòu)；利用DNAMAN軟件對(duì)其編碼的氨基酸序列進(jìn)行比對(duì)；構(gòu)建系統(tǒng)進(jìn)化樹(shù)利用MEGA 6.06軟件中鄰接法法進(jìn)行；使用PlantCARE(http://bioinformatics.psb.ugent. be/webtools/plantcare/html/)分析啟動(dòng)子上的順式作用元素。

1.2.3 基因表達(dá)分析分別提取‘赤霞珠’的根、莖、葉、芽及果皮的RNA，各取1 μg為模板，使用HiScript Ⅱ Reverse Transcriptase反轉(zhuǎn)錄試劑盒反轉(zhuǎn)錄合成cDNA第一鏈，于-20 ℃長(zhǎng)期保存。以葡萄a(bǔ)ctin7為內(nèi)參基因并設(shè)計(jì)引物(表1)，用RT-qPCR法分析bHLH93在不同組織和果實(shí)不同發(fā)育時(shí)期的表達(dá)量。將cDNA稀釋10倍，參考SYBR qPCR Master Mix(南京諾唯贊生物科技股份有限公司，南京)試劑盒的反應(yīng)體系進(jìn)行加樣。用BIO-RAD Cycler IQ6熒光定量PCR儀檢測(cè)基因表達(dá)量，重復(fù)測(cè)定3次。采用2-ΔΔCT法和相關(guān)軟件進(jìn)行數(shù)據(jù)處理與分析。

1.3 數(shù)據(jù)處理與分析利用Excel 2016和Origin 2021軟件進(jìn)行數(shù)據(jù)處理與分析，采用DPS 6.55軟件進(jìn)行顯著性差異分析。

2 結(jié)果與分析

2.1 VvbHLH93基因克隆及編碼蛋白的理化特性

采用PCR擴(kuò)增基因產(chǎn)物檢測(cè)得到一條單一明亮的條帶(圖1)。經(jīng)測(cè)序并NCBI blast比對(duì)，克隆得到的基因編碼核酸序列全長(zhǎng)1 319 bp，最大開(kāi)放閱讀框939 bp，共編碼312個(gè)氨基酸(圖2)。ProtParam在線預(yù)測(cè)分析顯示，VvbHLH93分子量為35.18 kD，理論等電點(diǎn)為4.68，表明為酸性蛋白；蛋白質(zhì)不穩(wěn)定性指數(shù)56.22，表現(xiàn)為不穩(wěn)定蛋白；通過(guò)親/疏水性預(yù)測(cè)發(fā)現(xiàn)，屬于親水性蛋白；未發(fā)現(xiàn)跨膜結(jié)構(gòu)，無(wú)信號(hào)肽，不屬于分泌蛋白。二級(jí)結(jié)構(gòu)預(yù)測(cè)發(fā)現(xiàn)，在VvbHLH93氨基酸組成中，以無(wú)規(guī)則卷曲(random coil)為主，占50.96%，其次為α-螺旋(α-helix)，占38.78%，第三位的是β-螺旋(β-helix)，占2.24%。采用SWISS-MODEL在線軟件預(yù)測(cè)VvbHLH93蛋白三級(jí)結(jié)構(gòu)顯示，該蛋白主要由無(wú)規(guī)則卷曲、α-螺旋和β-螺旋組成，與二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果一致。NetPhos-3.1預(yù)測(cè)分析發(fā)現(xiàn)，VvbHLH93蛋白具有多個(gè)磷酸化位點(diǎn)。其中，絲氨酸磷酸化位點(diǎn)最多，有14個(gè)；其次為蘇氨酸磷酸化位點(diǎn)，有10個(gè)；而酪氨酸磷酸化位點(diǎn)只有3個(gè)。

圖1 VvbHLH93的cDNA擴(kuò)增Fig.1 Amplification of VvbHLH93 cDNA

圖2 VvbHLH93編碼氨基酸預(yù)測(cè)序列Fig.2 Deduced amino acid sequence of VvbHLH93

2.2 VvbHLH93進(jìn)化樹(shù)及氨基酸序列同源分析

為了解VvbHLH93與同源序列之間的進(jìn)化關(guān)系，基于Clustal W對(duì)同源序列進(jìn)行序列多重比對(duì)，在此基礎(chǔ)上利用MEGA 6.06構(gòu)建系統(tǒng)發(fā)育樹(shù)。結(jié)果(圖3)表明，VvbHLH93與荷花親緣關(guān)系較近，聚為一個(gè)亞類，在進(jìn)化上與蘋(píng)果、擬南芥親緣關(guān)系較遠(yuǎn)。

圖3 VvbHLH93與其他物種bHLH93蛋白的系統(tǒng)進(jìn)化樹(shù)分析Fig.3 Phylogenetic tree analysis of VvbHLH93 and bHLH93 in various species

將VvbHLH93氨基酸序列和擬南芥、蘋(píng)果、白梨、核桃等17種植物bHLH93氨基酸序列進(jìn)行同源分析。結(jié)果(圖4)表明，VvbHLH93與這些植物bHLH93氨基酸序列都擁有保守的bHLH域。經(jīng)氨基酸序列比對(duì)分析，VvbHLH93和其他植物相似性由高到低排列分別為：荷花69.43%，哥倫比亞錦葵62.30%，栓皮櫧62.30%，核桃62.21%，白梨61.93%，蘋(píng)果61.76%，可可61.64%，梅花60.98%，櫻桃60.67%，桃60.67%，草莓60.37%，月季59.76%，黃麻59.12%，毛果楊58.39%，胡楊58.07%，異色山黃麻56.80%，擬南芥55%。

紅色方框標(biāo)記bHLH結(jié)合保守基序圖4 VvbHLH93與其他物種bHLH93蛋白序列比對(duì)分析Red box marks bHLH binding conserved motifFig.4 Protein sequence alignment analysis of VvbHLH93 and bHLH93 in various species

2.3 VvbHLH93啟動(dòng)子順式元件的預(yù)測(cè)與分析

啟動(dòng)子順式元件預(yù)測(cè)表明，VvbHLH93基因的啟動(dòng)子含有激素信號(hào)和脅迫響應(yīng)的相關(guān)元件，包括光響應(yīng)元件(AE-box、GT1-motif、TCT-motif等)和低溫、干旱等應(yīng)答元件(表2)。表明該基因可能參與植物生理調(diào)控及對(duì)逆境的響應(yīng)過(guò)程。

表2 VvbHLH93啟動(dòng)子順式元件預(yù)測(cè)

不同小寫(xiě)字母表示不同組織基因的表達(dá)差異顯著(P<0.05)，下同圖5 VvbHLH93在赤霞珠不同組織中的相對(duì)表達(dá)量Different normal letters indicate differences between different tissues for the same gene at 0.05 level. The same as belowFig.5 The relative expression levels of VvbHLH93 in different tissues

2.4 VvbHLH93的表達(dá)分析

2.4.1 在葡萄不同組織中的表達(dá)如圖5所示，VvbHLH93在‘赤霞珠’不同組織中表達(dá)水平存在顯著差異，在根、莖、葉、芽、果皮不同組織中均有表達(dá)。從相對(duì)表達(dá)量來(lái)看，VvbHLH93在莖中的表達(dá)量最高，顯著高于VvbHLH93在其他組織中的表達(dá)水平；其次是芽和果皮，但VvbHLH93在芽和果皮中的表達(dá)量無(wú)顯著差異；VvbHLH93在根中的表達(dá)量較低，在葉片中的表達(dá)量顯著低于其他組織。說(shuō)明VvbHLH93在‘赤霞珠’不同組織中的表達(dá)具有組織特異性。

2.4.2 在葡萄果實(shí)不同發(fā)育時(shí)期的表達(dá)如圖6所示，隨著赤霞珠葡萄果實(shí)的發(fā)育成熟，VvbHLH93相對(duì)表達(dá)量不斷降低，第5周后相對(duì)表達(dá)量基本為0，總體呈現(xiàn)出降低的變化趨勢(shì)。

0.幼果期(直徑> 2 mm)；1～11依次為幼果期后1～11周采樣圖6 VvbHLH93在不同發(fā)育時(shí)期中的相對(duì)表達(dá)量0. Young fruit stage (diameter > 2 mm); 1—11 samples were collected at 1—11 weeks after young fruit stage Fig.6 The relative expression levels of VvbHLH93 in different stages

2.4.3 在不同脅迫處理下的表達(dá)通過(guò)順式作用元件分析，在VvbHLH93啟動(dòng)子序列中預(yù)測(cè)到了多個(gè)與植物激素和逆境脅迫相關(guān)的響應(yīng)元件，表明該基因可能在響應(yīng)非生物脅迫的轉(zhuǎn)錄調(diào)控機(jī)制發(fā)揮作用。進(jìn)一步通過(guò)qRT-PCR分析發(fā)現(xiàn)，VvbHLH93基因在高溫脅迫、鹽脅迫、低溫脅迫和充分灌溉脅迫下具有表達(dá)特異性(圖7)。脅迫試驗(yàn)數(shù)據(jù)表明，低溫、鹽、高溫和充分灌溉處理‘赤霞珠’盆栽苗6 h后，VvbHLH93的表達(dá)量顯著下調(diào)(P<0.05)，分別為對(duì)照的52%、10%、13%和7%。

CK. 對(duì)照；LW. 4 ℃低溫；NaCl. 鹽；HW. 40 ℃高溫；FI. 充分灌溉圖7 VvbHLH93在不同脅迫下的相對(duì)表達(dá)量CK. 28 ℃ control group； LW. 4 ℃ low temperature； NaCl. Salt；HW. 40 ℃ high temperature； FI. Full irrigation Fig.7 The relative expression levels of VvbHLH93 under different stresses

3 討 論

bHLH轉(zhuǎn)錄因子是植物細(xì)胞生長(zhǎng)發(fā)育、應(yīng)對(duì)各種脅迫的關(guān)鍵調(diào)控因子，具有促進(jìn)植物開(kāi)花[16]、調(diào)節(jié)植物抗旱性[17-18]和耐鹽性[19-20]等重要作用。目前葡萄bHLH轉(zhuǎn)錄因子的功能研究相對(duì)較少，還有大量成員的功能未知。

本研究以釀酒葡萄‘赤霞珠’為研究對(duì)象，克隆獲得VvbHLH93的全長(zhǎng)cDNA序列，并對(duì)其進(jìn)行了初步的生物信息學(xué)分析和表達(dá)特性分析。該基因編碼312個(gè)氨基酸，具有高度保守的bHLH結(jié)構(gòu)域，屬于bHLH轉(zhuǎn)錄因子家族成員。氨基酸序列進(jìn)行同源性分析和系統(tǒng)進(jìn)化樹(shù)分析表明VvbHLH93與荷花bHLH93同源關(guān)系最近。

qRT-PCR分析結(jié)果表明VvbHLH93在‘赤霞珠’不同組織中表達(dá)水平存在顯著差異，從相對(duì)表達(dá)量來(lái)看，VvbHLH93在莖中的表達(dá)量最高，在葉片中的表達(dá)量最低。表明VvbHLH93在‘赤霞珠’不同組織中的表達(dá)具有組織特異性。與煙草的研究存在不同，NtbHLH93在根中表達(dá)量最高[21]。VvbHLH93在葡萄果實(shí)不同發(fā)育時(shí)期中的表達(dá)量差異比較顯著，發(fā)育后期基本不表達(dá)，總體呈現(xiàn)降低的變化趨勢(shì)。

AtbHLH93已被證實(shí)是參與調(diào)節(jié)擬南芥生長(zhǎng)發(fā)育和開(kāi)花時(shí)間的重要因子之一[10,22]，NtbHLH93基因在煙草甾醇代謝中發(fā)揮著重要作用[23]，在蘋(píng)果中，MdbHLH93和MdBT2蛋白相互作用，并激活下游結(jié)構(gòu)基因MdSAG18的轉(zhuǎn)錄而負(fù)調(diào)節(jié)葉片的衰老[11]。然而，關(guān)于VvbHLH93功能的相關(guān)研究報(bào)道較少，Muoz-Espinoza等[15]通過(guò)qPCR實(shí)驗(yàn)研究發(fā)現(xiàn)VvbHLH93在大小不同的漿果中存在差異表達(dá)，可能參與調(diào)控漿果重量，改善漿果品質(zhì)；Song等[24]研究發(fā)現(xiàn)VvbHLH93與MYB激活劑和抑制劑相關(guān)，并且受低溫脅迫的誘導(dǎo)。此外，轉(zhuǎn)錄組學(xué)分析發(fā)現(xiàn)，外源碳處理和光照處理后VvbHLH93的表達(dá)量發(fā)生變化，說(shuō)明其可能參與調(diào)節(jié)葡萄的生長(zhǎng)發(fā)育及品質(zhì)形成[25-26]。在VvbHLH93啟動(dòng)子順式元件預(yù)測(cè)中，VvbHLH93基因的啟動(dòng)子含有多種激素信號(hào)和脅迫響應(yīng)的相關(guān)元件，包括光響應(yīng)元件和低溫、干旱等應(yīng)答元件。在煙草的研究發(fā)現(xiàn)，bHLH93的表達(dá)會(huì)受光的影響而上調(diào)[27]。以上研究為進(jìn)一步明確VvbHLH93的功能提供了思路。在本研究中，我們發(fā)現(xiàn)VvbHLH93響應(yīng)高溫、鹽、低溫和充分灌溉脅迫，與對(duì)照相比，脅迫處理下VvbHLH93基因表達(dá)量顯著下調(diào)，尤其是充分灌溉脅迫下VvbHLH93基因表達(dá)量下降最多，僅為對(duì)照的7%，這與已有的研究一致[28]。

綜上所述，本研究以釀酒葡萄‘赤霞珠’為試驗(yàn)材料，克隆獲得了轉(zhuǎn)錄因子bHLH轉(zhuǎn)錄因子家族成員之一VvbHLH93，并發(fā)現(xiàn)其在不同組織和果實(shí)不同發(fā)育階段具有表達(dá)特異性。通過(guò)對(duì)VvbHLH93在多種非生物脅迫下的表達(dá)模式分析發(fā)現(xiàn)，VvbHLH93不同程度地響應(yīng)高溫、鹽、低溫和充分灌溉脅迫，尤其是在鹽、低溫和充分灌溉脅迫下，該基因的表達(dá)量顯著下調(diào)。因此推測(cè)VvbHLH93是負(fù)調(diào)控葡萄抗逆脅迫的關(guān)鍵轉(zhuǎn)錄因子。本研究為了解VvbHLH93基因在葡萄生長(zhǎng)發(fā)育和抗逆境機(jī)制提供了一定的理論參考，以便為葡萄的育種和抗逆性增強(qiáng)工作奠定基礎(chǔ)。