鹽堿脅迫下野生大豆CBL-CIPK通路表達(dá)分析

魏志園 王 宇 司增志 魏 萊 溫曉蕾 喬亞科 李桂蘭 張 鍇

(1河北科技師范學(xué)院農(nóng)學(xué)與生物科技學(xué)院，河北秦皇島 066000；2河北科技師范學(xué)院園藝科技學(xué)院，河北秦皇島 066000)

我國是大豆的起源地，保存有豐富的野生大豆(Glycine soja)資源。野生大豆常生長在沿海灘涂、鹽堿荒地等環(huán)境中，是典型的鹽生植物。在野生大豆優(yōu)質(zhì)基因利用方面，我國學(xué)者通過野生大豆與栽培大豆(Glycine max)雜交育種技術(shù)已培育出一些優(yōu)異的種質(zhì)資源[1－6]。雖然我國擁有豐富的野生大豆種質(zhì)資源，但我國對野生大豆的育種利用還不充分。在中美貿(mào)易談判及國內(nèi)非轉(zhuǎn)基因大豆需求量增加的背景下，進(jìn)一步充分利用野生大豆種質(zhì)資源，應(yīng)用傳統(tǒng)技術(shù)培育栽培大豆新品種成為解決我國大豆供需矛盾的重要途徑。

近年來野生大豆的優(yōu)異特性尤其是耐鹽堿性已成為學(xué)者的研究重點(diǎn)。野生大豆對鹽堿脅迫的反應(yīng)比較復(fù)雜[7]。肖鑫輝等[8]研究發(fā)現(xiàn)，高耐鹽堿型野生大豆材料在鹽堿脅迫下莖、葉平均Na＋和Cl－含量顯著低于低耐鹽型，但高耐型野生大豆植株莖葉中也存在Na＋和Cl－含量高水平和低水平兩種類型，推測野生大豆可能存在高耐受性和低吸收性兩種耐鹽堿機(jī)制。此外，研究發(fā)現(xiàn)野生大豆可以通過增加超氧化物歧化酶(superoxide dismutase，SOD)、過氧化氫酶(catalase，CAT)及過氧化物酶(peroxidase，POD)等抗氧化酶活性，降低膜脂過氧化物丙二醛(malondialdehyde，MDA)含量來提高自身耐鹽性[9－10]。野生大豆具有較強(qiáng)的抗逆能力，近年來學(xué)者對其耐鹽堿性進(jìn)行了諸多探索[11－15]。但是，目前野生大豆耐鹽堿生理機(jī)制，尤其是野生大豆鹽堿脅迫下植株體內(nèi)的離子平衡機(jī)制還不清楚，耐鹽堿關(guān)鍵基因表達(dá)及調(diào)控機(jī)制尚不明確。

類鈣調(diào)磷酸酶B 蛋白(calcineurin B-like protein，CBLs)是植物中重要的Ca2＋離子感受器。擬南芥CBLs蛋白含有4 個EF 手結(jié)構(gòu)，是Ca2＋離子的結(jié)合位點(diǎn)，它們的排列順序和間距高度保守[16－17]。CBLs 蛋白N－端具有豆蔻酰化(myristoylation)和棕櫚?；?palmitoylation)兩個脂質(zhì)修飾位點(diǎn)，介導(dǎo)CBLs 蛋白在胞內(nèi)的定位。CBLs 互作蛋白激酶(CBL-interacting protein kinase，CIPKs)是一類植物特有的絲氨酸－蘇氨酸蛋白激酶，其蛋白的C 端含有天冬酰胺－丙氨酸－苯丙酰胺(NAF)結(jié)構(gòu)域，介導(dǎo)CIPKs 蛋白與CBLs 蛋白的互作；N 端含有絲氨酸/蘇氨酸蛋白激酶結(jié)構(gòu)域，具有激酶催化活性[17]。目前研究發(fā)現(xiàn)，擬南芥(Arabidopsis)中存在10 個CBLs蛋白和26 個CIPKs 蛋白，栽培大豆(G.max)中存在7個CBLs 蛋白和13 個CIPKs 蛋白[18－20]。擬南芥中CBLs蛋白定位在質(zhì)膜、液泡膜、細(xì)胞質(zhì)及細(xì)胞核中[21－24]，而由于CIPKs 蛋白缺少定位信號結(jié)構(gòu)域，在正常狀態(tài)下CIPKs 定位在細(xì)胞質(zhì)和細(xì)胞核中，當(dāng)CIPKs 與CBLs 蛋白互作形成CBL-CIPK 復(fù)合體會重新定位至質(zhì)膜或液泡膜上[21，25]。

擬南芥AtCBL4 基因(也稱salt-overly-sensitive 3，SOS3)可以響應(yīng)多種逆境，包括鹽堿、低鉀、干旱和低溫等脅迫。在脅迫條件下，植株胞內(nèi)Ca2＋積累，Ca2＋可介導(dǎo)AtCBL4 蛋白與AtCIPK24 蛋白(也稱SOS2)形成AtCBL4-AtCIPK24 復(fù)合物[26]，AtCBL4 蛋白N 端豆蔻?；閷?dǎo)AtCBL4-AtCIPK24 復(fù)合物移動到質(zhì)膜上。AtCBL4-AtCIPK24 復(fù)合物可磷酸化質(zhì)膜上的SOS1(salt-overly-sensitive 1)蛋白[27－28]。SOS1 是定位在質(zhì)膜上的Na＋/H＋反向轉(zhuǎn)運(yùn)蛋白，在鹽脅迫下，AtCIPK24可磷酸化SOS1 自抑制結(jié)構(gòu)域中的絲氨酸殘基，激活SOS1 的轉(zhuǎn)運(yùn)活性，導(dǎo)致Na＋轉(zhuǎn)運(yùn)出細(xì)胞，這對高鹽脅迫下保持細(xì)胞Na＋平衡，維持細(xì)胞滲透壓十分重要[29]。另外，AtCBL10 也可與AtCIPK24 互作，形成的AtCBL10-AtCIPK24 復(fù)合物移動到液泡膜上，激活A(yù)tNHX(vacuolar Na＋/H＋exchanger)蛋白活性。NHX是存在于液泡膜上的Na＋/H＋轉(zhuǎn)運(yùn)蛋白，激活后可使胞質(zhì)中的Na＋轉(zhuǎn)運(yùn)至液泡中，進(jìn)一步降低細(xì)胞內(nèi)的Na＋濃度。液泡中積累的Na＋也可以作為滲透壓降低細(xì)胞水勢，促進(jìn)植物在鹽脅迫下吸水[30]。由此可見，CBLCIPK 信號通路對植物鹽脅迫的抗性十分重要[31－32]。

本研究前期對野生大豆材料進(jìn)行耐鹽堿鑒定，篩選到高耐鹽堿野生大豆材料2010－12。為進(jìn)一步揭示野生大豆耐鹽堿機(jī)制，選擇高耐鹽堿野生大豆2010－12 和敏感型野生大豆2012－49 為試驗材料，在鹽堿脅迫處理后0(鹽堿脅迫處理后立即取樣)、3、10 h 進(jìn)行轉(zhuǎn)錄組分析。鑒于CBL-CIPK 信號通路在植物耐鹽堿脅迫中發(fā)揮重要作用，且野生大豆CBL-CIPK 信號通路對鹽堿脅迫的響應(yīng)機(jī)制尚不清楚，本研究在轉(zhuǎn)錄組分析基礎(chǔ)上，進(jìn)一步分析野生大豆中CBL-CIPK 信號通路的蛋白結(jié)構(gòu)、亞細(xì)胞定位、磷酸化位點(diǎn)及鹽堿脅迫下高耐/敏感野生大豆材料基因表達(dá)特點(diǎn)，旨在揭示CBL-CIPK 信號通路基因在野生大豆耐鹽堿脅迫中的作用。

1 材料與方法

1.1 試驗材料及方法

高耐鹽堿野生大豆材料2010－12 和敏感野生大豆材料2012－49，由河北科技師范學(xué)院野生大豆課題組提供。

用小刀劃破野生大豆種子泥漿膜后，播種于鋪滿營養(yǎng)土的培養(yǎng)缽中，上覆薄層蛭石。培養(yǎng)缽置于溫室(25℃光16 h/15℃暗8 h)。待野生大豆長出真葉后，移栽至裝有30 mL 霍格蘭培養(yǎng)液的50 mL 離心管中，置于相同光照培養(yǎng)條件的培養(yǎng)箱中，濕度60%，5 d 后進(jìn)行鹽堿脅迫處理。處理液按NaCl ∶Na2CO3(無水)＝1 ∶3 質(zhì)量分?jǐn)?shù)配制[33]，調(diào)pH 值至8.0，濃度為100 mmol·L－1。處理0(澆灌處理液后立即取樣)、3、10 h后分別取兩種野生大豆材料的葉片1 g，迅速保存于液氮中，所有樣品設(shè)3 次生物學(xué)重復(fù)。

1.2 測序數(shù)據(jù)來源

提取野生大豆葉片RNA，構(gòu)建cDNA 文庫，送至華大基因(北京)利用BGISEQ－500 平臺進(jìn)行轉(zhuǎn)錄組測序。利用Trimmomatic(v0.36)軟件對原始數(shù)據(jù)進(jìn)行過濾，HISAT2(v2.1.0)軟件比對野生大豆參考基因序列。通過華大基因交互系統(tǒng)(https:/ /report.bgi.com/ps/login/login.html)篩選CBL 和CIPK 通路基因，并利用栽培大豆中CBL 和CIPK 通路基因序列與本地轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行BLAST 比對。

1.3 預(yù)測數(shù)據(jù)來源

通過NCBI(http:/ /www.ncbi.nlm.nih.gov/)獲得CBLs和CIPKs基因的蛋白序列，利用MEME(http:/ /alternate.meme-suite.org/tools/meme)在線分析工具進(jìn)行蛋白結(jié)構(gòu)域預(yù)測。利用BUSCA (http:/ /busca.biocomp.unibo.it/)在線分析工具進(jìn)行蛋白亞細(xì)胞定位預(yù)測。利用NetPhos 3.1 Server(http:/ /www.cbs.dtu.dk/services/NetPhos/)在線分析工具進(jìn)行蛋白磷酸化位點(diǎn)預(yù)測。

1.4 數(shù)據(jù)分析

CBLs和CIPKs基因表達(dá)量采用FPKM(fregments per kilobase per million)法進(jìn)行計算，取｜log2FoldChange ｜>1.5 為組間差異顯著。利用Heml(v 1.0.3.7)軟件對取對數(shù)值(log2)的基因表達(dá)量進(jìn)行聚類分析。利用MEGA7 軟件對野生大豆和栽培大豆的CBL-CIPK基因進(jìn)行系統(tǒng)發(fā)育分析。

2 結(jié)果與分析

2.1 野生大豆CBL-CIPK 通路基因鑒定

對高耐鹽堿及敏感野生大豆材料鹽堿脅迫處理0、3、10 h 后的轉(zhuǎn)錄組測序，結(jié)果發(fā)現(xiàn)，18 組樣品(3 次生物學(xué)重復(fù)) 共鑒定819.52 M Clean Reads，包括122.91 Gb Clean Bases。通過與已公布的野生大豆(G.soja)基因組比對，鑒定到已知轉(zhuǎn)錄本48 875 個，新轉(zhuǎn)錄本32 331 個。新轉(zhuǎn)錄本中具有蛋白編碼潛力的轉(zhuǎn)錄本25 036 個，新基因1 780 個。

鑒定到CBL-CIPK 信號通路基因共48 個，其中CBLs基因10 個，CIPKs基因38 個。對篩選到的基因進(jìn)行KO 注釋，發(fā)現(xiàn)9 個CBLs基因匹配到K06268(絲氨酸/蘇氨酸蛋白磷酸酶2B 調(diào)節(jié)亞基，serine/threonine-protein phosphatase 2B regulatory subunit)中。30 個CIPKs基因匹配到K07198(5′-AMP 激活蛋白激酶，催化α 亞單位，5′-AMP-activated protein kinase，catalytic alpha subunit)中，14 個匹配到K12761(碳代謝解壓蛋白激酶，carbon catabolite-derepressing protein kinase)中(附表1)。

附表1 轉(zhuǎn)錄組分析鑒定的野生大豆CBL-CIPK 信號通路基因表達(dá)量Enclose Table 1 Gene expression of wild soybean CBL-CIPK proteins indentified by trascriptome analysis

2.2 野生大豆CBL-CIPK 蛋白結(jié)構(gòu)域預(yù)測

對野生大豆CBLs 蛋白結(jié)構(gòu)域預(yù)測發(fā)現(xiàn)，10 個CBLs 蛋白均含有3 個EF 手結(jié)構(gòu)(圖1-A)。雖然EF手結(jié)構(gòu)在10 個CBLs 蛋白氨基酸序列的位置不同，但3 個EF 手結(jié)構(gòu)之間的排列及它們之間的距離均一致。分析EF1、EF2 和EF3 蛋白基序發(fā)現(xiàn)，10 個EF 手結(jié)構(gòu)氨基酸序列保守程度較高，說明該結(jié)構(gòu)域?qū)BLs 蛋白發(fā)揮功能較重要。

對野生大豆38 個CIPKs 蛋白結(jié)構(gòu)域預(yù)測發(fā)現(xiàn)，32個CIPKs 蛋白含有ATP 結(jié)合結(jié)構(gòu)域(ATP binding domain)，36 個含有絲氨酸/蘇氨酸蛋白激酶活性結(jié)構(gòu)域(serine/threonine protein kinase active domain，ST)，17 個含有天冬酰胺－丙氨酸－苯丙氨酸保守結(jié)構(gòu)域(Asn-Ala-Phe 結(jié)構(gòu)域，NAF 結(jié)構(gòu)域)。從3 個結(jié)構(gòu)域基序的分析結(jié)果可知(圖1－B)，ATP 結(jié)合結(jié)構(gòu)域和絲氨酸/蘇氨酸蛋白激酶活性結(jié)構(gòu)域高度保守，說明這兩個結(jié)構(gòu)域?qū)IPKs 蛋白較重要。NAF 結(jié)構(gòu)域兩端不保守，只在Asn-Ala-Phe 氨基酸位點(diǎn)保守程度較高。除這3 個結(jié)構(gòu)域外，一些CIPKs 蛋白還含有一些保守結(jié)構(gòu)域，絲氨酸/蘇氨酸蛋白激酶活性結(jié)構(gòu)域兩端的結(jié)構(gòu)域，它們共同介導(dǎo)了CIPKs 的蛋白激酶活性。

圖1 野生大豆CBL-CIPK 蛋白結(jié)構(gòu)域及重要結(jié)構(gòu)域基序預(yù)測結(jié)果Fig.1 Prediction results of CBL-CIPK proteins domain and domain motif of wild soybean

2.3 野生大豆CBL-CIPK 蛋白亞細(xì)胞定位預(yù)測

對10 個CBLs 蛋白亞細(xì)胞定位預(yù)測發(fā)現(xiàn)，7 個CBLs蛋白位于細(xì)胞核中，1 個CBLs 蛋白位于細(xì)胞質(zhì)中，1 個位于葉綠體中，1 個位于內(nèi)膜系統(tǒng)上。對38 個CIPKs 蛋白亞細(xì)胞定位預(yù)測發(fā)現(xiàn)，20 個CIPKs 蛋白定位在細(xì)胞質(zhì)中，占52.6%，9 個定位在葉綠體中，占23.6%，6 個定位在細(xì)胞核中，占15.7%，另有3 個分別位于細(xì)胞器膜、線粒體和內(nèi)膜系統(tǒng)上，分別均占2.6%(表1)。

表1 CBL-CIPK 蛋白亞細(xì)胞定位及磷酸化位點(diǎn)預(yù)測結(jié)果Table 1 Subcellular localization and phosphorylation site prediction of CBL-CIPK proteins

表1(續(xù))

2.4 野生大豆CBL-CIPK 蛋白磷酸化位點(diǎn)預(yù)測

對CBL-CIPK 蛋白磷酸化位點(diǎn)預(yù)測發(fā)現(xiàn)，CBLs 蛋白磷酸化修飾以絲氨酸(Ser)位點(diǎn)最多，占總數(shù)的52.3%～73.1%，蘇氨酸(Thr)位點(diǎn)次之，占19.2%～45.8%，酪氨酸(Tyr)位點(diǎn)最少，占0～14.3%。CIPKs蛋白磷酸化修飾也以絲氨酸(Ser)位點(diǎn)最多，占40.6%～77.8%，蘇氨酸(Thr)位點(diǎn)次之，占13.3%～46.7%，酪氨酸(Tyr)位點(diǎn)最少，占0～22.7%(表1)。

10 個CBLs 蛋白平均每個蛋白發(fā)生12.0 次絲氨酸(Ser)位點(diǎn)磷酸化修飾，7.0 次蘇氨酸(Thr)位點(diǎn)磷酸化修飾，0.8 次酪氨酸(Tyr)位點(diǎn)磷酸化修飾；38 個CIPKs 蛋白平均每個蛋白發(fā)生20.6 次絲氨酸(Ser)位點(diǎn)磷酸化修飾，9.9 次蘇氨酸(Thr)位點(diǎn)磷酸化修飾，3.9 次酪氨酸(Tyr)位點(diǎn)磷酸化修飾(表1)。

圖2 鹽堿脅迫下野生大豆CBLs 蛋白(A)和CIPKs 蛋白(B)基因表達(dá)聚類分析Fig.2 The cluster analysis of gene expression of CBLs protein (A) and CIPKs protein (B) in wild soybean under saline-alkali stress

2.5 野生大豆鹽堿脅迫下CBL-CIPK 基因表達(dá)模式分析

對10 個CBLs基因表達(dá)模式進(jìn)行分析發(fā)現(xiàn)，高耐鹽堿野生大豆材料GsCBL3(D0Y65_047783)基因在鹽堿脅迫處理后3、10 h 的表達(dá)量與0 h 相比下調(diào)。高耐鹽堿野生大豆材料GsCBL1 (D0Y65 _045867) 和GsCBL4(D0Y65_014543)基因在鹽堿脅迫處理后10 h的表達(dá)量與0、3 h 相比均上調(diào)。兩個耐性不同的野生大豆材料比較，在鹽堿脅迫處理0、10 h 后，GsCBL10(D0Y65_046727)基因在高耐鹽堿材料中的表達(dá)量比敏感材料上調(diào)，其余基因各時間點(diǎn)相比均無差別。

38 個CIPKs基因表達(dá)模式分析發(fā)現(xiàn)，兩個耐性不同的野生大豆材料對比，在鹽堿脅迫處理0 h 后，GsCIPK10_like(D0Y65_004857)基因在高耐鹽堿野生大豆材料中的表達(dá)量比敏感材料上調(diào)，其余基因各時間點(diǎn)相比均無差別。

2.6 野生大豆與栽培大豆CBL-CIPK 基因系統(tǒng)發(fā)育分析

通過NCBI 檢索到栽培大豆(G.max)CBLs基因6個，CIPKs基因15 個。利用MEGA7 軟件分別分析野生大豆與栽培大豆CBL家族基因和CIPK家族基因的進(jìn)化關(guān)系。結(jié)果發(fā)現(xiàn)，野生大豆與栽培大豆CBLs基因親緣關(guān)系較近，野生大豆與栽培大豆的CBL1、CBL3、CBL4、CBL7、CBL10 分別聚為一類，證明其同源性很高。野生大豆與栽培大豆CIPKs基因也高度同源，野生大豆與栽培大豆的CIPK1、CIPK2、CIPK3、CIPK4、CIPK5、CIPK6、CIPK8、CIPK9、CIPK10、CIPK11、CIPK12、CIPK14和CIPK23 分別聚為一類，證明這些基因同源關(guān)系很高。但野生大豆GsCIPK5_like 基因與栽培大豆GmCIPK5 親緣關(guān)系較遠(yuǎn)，而與栽培大豆GmCIPK25 基因親緣關(guān)系較近。栽培大豆GmCIPK21 基因與野生大豆GsCIPK21 親緣關(guān)系較遠(yuǎn)，而與野生大豆GsCIPK1_like 基因親緣關(guān)系較近(圖3)。

3 討論

Meador 等[34]和Sánchez-Barrena 等[35]指出，擬南芥中CBLs 蛋白含有4 個EF 手結(jié)構(gòu)，這些EF 手結(jié)構(gòu)是CBLs 蛋白結(jié)合Ca2＋的位置，對其發(fā)揮正常功能十分重要。蛋白結(jié)構(gòu)域預(yù)測結(jié)果表明，野生大豆10 個CBLs 通路蛋白均含有3 個EF 手結(jié)構(gòu)。通過比較擬南芥與野生大豆的EF 手結(jié)構(gòu)，發(fā)現(xiàn)擬南芥和野生大豆的EF1 結(jié)構(gòu)相同，野生大豆EF2 序列較長，包括了擬南芥EF2 和EF3 的序列，野生大豆EF3 與擬南芥EF4結(jié)構(gòu)一致。曹齊衛(wèi)等[36]對黃瓜(Cucumis sativus)CBLs蛋白結(jié)構(gòu)分析也發(fā)現(xiàn)其含有3 個EF 手結(jié)構(gòu)。說明野生大豆可能與黃瓜CBLs 蛋白結(jié)構(gòu)相似。CIPKs 蛋白催化結(jié)構(gòu)域有許多保守的區(qū)域，其中最重要的是ATP結(jié)合結(jié)構(gòu)域和絲氨酸/蘇氨酸蛋白激酶活性結(jié)構(gòu)域[37]。這兩個結(jié)構(gòu)域分別具有絲氨酸/蘇氨酸激酶特異性和酪氨酸激酶的特異性。除此之外，CIPKs 蛋白的NAF 結(jié)構(gòu)域由25 個氨基酸構(gòu)成，由其保守的天冬酰胺－丙氨酸－苯丙氨酸(Asn-Ala-Phe)殘基命名，其功能是與CBLs 蛋白相互作用[21]。本研究通過對野生大豆的CIPKs 蛋白結(jié)構(gòu)域基序分析發(fā)現(xiàn)，CIPKs 蛋白保守度最高的結(jié)構(gòu)域為絲氨酸/蘇氨酸蛋白激酶活性結(jié)構(gòu)域，其次為ATP 結(jié)合結(jié)構(gòu)域，NAF 結(jié)構(gòu)域保守度最低，這與擬南芥中CIPKs 蛋白結(jié)構(gòu)基本一致。

CBLs 蛋白在細(xì)胞內(nèi)的位置受蛋白N－端影響。擬南芥中AtCBL1、AtCBL4、AtCBL5 和AtCBL9 蛋白N－端均含有一小段定位于質(zhì)膜的N－豆蔻?；Y(jié)構(gòu)域，AtCBL2、AtCBL3、AtCBL6 和AtCBL10 蛋白含有一段延長的N－端結(jié)構(gòu)域，定位在液泡膜上。AtCBL7 和AtCBL8 蛋白由于其N－端結(jié)構(gòu)域突變而定位在細(xì)胞質(zhì)和細(xì)胞核中[38]。另外，Oh 等[23]研究指出AtCBL3 定位于細(xì)胞質(zhì)中。由于CIPKs 不含有胞內(nèi)定位信號，當(dāng)其獨(dú)立存在時可能定位在細(xì)胞核、細(xì)胞質(zhì)和質(zhì)膜上。當(dāng)CIPKs 與CBLs 互作形成CBL-CIPK 復(fù)合體后，CIPKs 的胞內(nèi)定位則由CBLs 蛋白決定[21，25]。如AtCBL4-AtCIPK24 復(fù)合物定位于質(zhì)膜，AtCBL2-AtCIPK1 定位于液泡膜[32]。本研究通過對野生大豆CBL-CIPK 蛋白的亞細(xì)胞定位預(yù)測發(fā)現(xiàn)，GsCBL1－2、GsCBL3 isoform X1、GsCBL4、GsCBL7、GsCBL9、GsCBL10 和GsCBL10－3 定位于細(xì)胞核中，這是野生大豆與擬南芥CBLs 蛋白胞內(nèi)定位的主要差異，代表了野生大豆的特異性。GsCBL3 定位于細(xì)胞質(zhì)中，與擬南芥中AtCBL3 的胞內(nèi)定位一致。但這僅是蛋白亞細(xì)胞定位的預(yù)測結(jié)果，野生大豆CBL-CIPK 蛋白在胞內(nèi)的定位還需試驗證實。

對高耐鹽堿和敏感野生大豆材料在鹽堿脅迫下CBL-CIPK基因轉(zhuǎn)錄水平進(jìn)行檢測，重點(diǎn)分析脅迫后同一時間點(diǎn)兩材料間轉(zhuǎn)錄水平差異的基因。結(jié)果發(fā)現(xiàn)，兩個野生大豆材料差異表達(dá)基因全部集中在脅迫處理后0 h，且高耐鹽堿野生大豆材料中GsCBL10 和GsCIPK10 基因轉(zhuǎn)錄水平在脅迫處理后0 h 與敏感材料相比均增加。由于0 h 是脅迫處理后立即取樣檢測，故推測耐性野生大豆在受到鹽堿脅迫后CBL-CIPK基因響應(yīng)十分迅速。余義和等[39]研究了葡萄(Vitis vinifera)在干旱、低溫和鹽脅迫下VvCIPK10 基因的表達(dá)量，發(fā)現(xiàn)VvCIPK10 受到干旱、低溫和鹽脅迫誘導(dǎo)表達(dá)，并且干旱和鹽脅迫后2 h 基因表達(dá)量即達(dá)到峰值。Kim 等[40]研究發(fā)現(xiàn)，擬南芥的 AtCBL10 可與AtCIPK24 直接互作并轉(zhuǎn)移至液泡膜上，激活液泡膜上AtNHX1(液泡Na＋/H＋逆轉(zhuǎn)運(yùn)子)的活性，可介導(dǎo)胞質(zhì)中的Na＋轉(zhuǎn)運(yùn)至液泡中，這對增強(qiáng)擬南芥在鹽脅迫下的耐性有重要作用。高耐鹽堿野生大豆材料中GsCBL10 和GsCIPK10 的上調(diào)表達(dá)可能受到鹽堿脅迫的誘導(dǎo)，這可能與野生大豆對鹽堿脅迫的抗性有關(guān)，但需要進(jìn)一步試驗證明。

通過對野生大豆與栽培大豆CBL-CIPK基因系統(tǒng)發(fā)育分析，發(fā)現(xiàn)野生大豆與栽培大豆的CBL-CIPK基因親緣關(guān)系較近。但是，野生大豆與栽培大豆在進(jìn)化上也存在差異，如野生大豆GsCIPK5_like 基因和栽培大豆GmCIPK5 親緣關(guān)系較遠(yuǎn)，栽培大豆GmCIPK21 基因與野生大豆GsCIPK21 親緣關(guān)系較遠(yuǎn)。值得注意的是，轉(zhuǎn)錄組分析表明，高耐鹽堿野生大豆中GsCIPK5_like 基因在鹽堿脅迫后0 h 表達(dá)量高達(dá)109.88，而敏感材料在同時間點(diǎn)的表達(dá)量只有37.12，說明野生大豆GsCIPK5_like 基因受鹽堿脅迫的誘導(dǎo)，其在脅迫后表達(dá)量迅速上調(diào)可能對野生大豆耐鹽堿脅迫響應(yīng)較重要。

4 結(jié)論

本研究在野生大豆鹽堿脅迫處理轉(zhuǎn)錄組分析基礎(chǔ)上，鑒定到10 個野生大豆CBLs基因和38 個CIPKs基因。野生大豆CBLs 蛋白均含有3 個EF 手結(jié)構(gòu)，CIPKs 蛋白含有ATP 結(jié)合結(jié)構(gòu)域、絲氨酸/蘇氨酸蛋白激酶活性結(jié)構(gòu)域和NAF 結(jié)構(gòu)域。CBLs 蛋白主要定位于細(xì)胞核中，CIPKs 蛋白定位在細(xì)胞質(zhì)、葉綠體、細(xì)胞核、細(xì)胞器膜上等。耐性野生大豆CBL-CIPK 通路基因響應(yīng)鹽堿脅迫主要發(fā)生在脅迫處理后0 h，推測野生大豆CBL-CIPK基因響應(yīng)鹽堿脅迫十分迅速。野生大豆與栽培大豆CBL-CIPK 通路基因親緣關(guān)系很近，但存在個別親緣關(guān)系較遠(yuǎn)的基因。

附表1(續(xù))