聶佩顯 王來平 韓雪平 路 超 薛曉敏

(山東省果樹研究所，山東 泰安 271000)

山梨醇不僅是薔薇科果樹葉片光合作用的主要產(chǎn)物，而且還是韌皮部碳水化合物的主要運輸物質(zhì)[1-2]，與植物的生長、果實品質(zhì)與產(chǎn)量等密切相關[3]，在響應生物脅迫和非生物脅迫過程中具有重要作用[4]。山梨醇在薔薇科植物源葉中生成，經(jīng)韌皮部運輸卸載至庫器官[1]。盡管韌皮部運輸?shù)纳嚼娲妓急壤芨撸嵌鄶?shù)薔薇科植物的果實并不積累山梨醇[4-5]。很大一部分山梨醇通過韌皮部卸載進入果實后，在山梨醇脫氫酶(sorbitol dehydrogenase，SDH)和山梨醇氧化酶(sorbitol oxidase，SOX)的催化下轉(zhuǎn)化為果糖和葡萄糖[3]。其中SDH 已被證明是薔薇科植物果實轉(zhuǎn)化山梨醇的最主要酶[5]。SDH 主要有兩種:一種是依賴NAD＋的山梨醇脫氫酶(NAD＋-dependent sorbitol dehydrogenase，NAD＋-SDH)，其主要作用是催化山梨醇與果糖的相互轉(zhuǎn)化；另一種是依賴NADP＋的山梨醇脫氫酶(NADP＋-dependent sorbitol dehydrogenase，NADP＋-SDH)，其主要作用是催化山梨醇與葡萄糖的相互轉(zhuǎn)化[3]。NAD＋-SDH 在催化山梨醇和果糖相互轉(zhuǎn)化時，伴隨著酶活性增強，果實中碳水化合物含量逐漸積累，并促進果實發(fā)育[3]。因此，NAD＋-SDH 在調(diào)控“庫”強過程中起著關鍵作用[5]。此外，基于擬南芥突變體的研究發(fā)現(xiàn)，SDH 能夠調(diào)節(jié)植物體內(nèi)多元醇的水平，進而影響植物對干旱脅迫的響應[6]。因此，研究SDH基因表達和酶活性變化的調(diào)控機制具有重要的意義。

Yamada 等[7]首先從蘋果果實中克隆得到NAD＋-SDH基因，通過原核表達研究證明該基因編碼的蛋白具有NAD＋-SDH 活性。目前已有在蘋果[8]、梨[9]、桃[10]及枇杷[11]等薔薇科植物中克隆該基因的報道。杏(Prunus armeniacaL.)作為薔薇科李屬的重要經(jīng)濟樹種，起源于中國[12]，現(xiàn)研究多集中在種質(zhì)資源創(chuàng)新利用[13-15]、經(jīng)濟性狀分析[16]、抗寒生理[17-19]等方面，鮮見有關杏NAD＋-SDH 的相關研究。本研究選用同源克隆的方法從金太陽杏中分離得到NAD＋-SDH基因，并結(jié)合生物信息學方法，初步分析其基因序列特征及編碼蛋白結(jié)構(gòu)特征，并研究其在杏果實發(fā)育過程中的表達特性及酶活性變化，旨在為進一步驗證和解析杏NAD＋-SDH基因的功能提供理論依據(jù)。

1 材料與方法

1.1 試驗材料

供試材料金太陽杏取自山東省果樹研究所苗圃(山東省泰安市，北緯36°11′，東經(jīng)117°6′，海拔153 m)。選取生長相對一致的3 株8年生成年樹，每棵樹在4 個方位隨機采取。果實自盛花后7 d 開始取樣，每隔2 周取一次，直至商品成熟期，即在花后7、21、35、49 和56 d 采樣，每次采樣20 個果實，取果肉液氮速凍，-80℃保存?zhèn)溆谩?/p>

1.2 RNA 提取和反轉(zhuǎn)錄cDNA

利用HiPure Plant RNA Maxi Kit (美基生物Magen，廣州)試劑盒提取供試材料總RNA。cDNA 第一鏈的合成按5X All-In-One RT MasterMix(艾可莘生物科技有限公司，北京)反轉(zhuǎn)錄試劑盒說明書步驟進行，反轉(zhuǎn)錄得到的cDNA 于-20℃保存?zhèn)溆谩?/p>

1.3 NAD+-SDH 基因的克隆

根據(jù)NCBI 登記的蘋果、梨、桃等薔薇科果樹的NAD＋-SDH基因序列，設計全長引物。

NAD＋-SDH-F(5′-CACGGGGGACTCTAGAATGGGC AAGGGAGGGATGT-3′)；

NAD＋-SDH-R(5′-GACCACCCGGGGATCCCAAATT AAACATGACTTTTATGGCA-3′)。

以杏果實cDNA 為模板進行PCR 擴增。反應體系50 μL:模板1 μL，上下游引物(10 μmol·L-1)各1.5 μL，dNTP(2.5 mmol·L-1)4 μL，酶(2.5 U·μL-1)0.5 μL，5×buffer 10 μL，ddH2O 補足至50 μL。PCR 反應條件:94℃預變性2 min；98℃變性10 s，55℃退火5 s，72℃延伸100 s，35 個循環(huán)；72℃延伸5 min。利用1%瓊脂糖凝膠電泳對擴增產(chǎn)物進行檢測，使用美基生物公司(廣州)凝膠回收試劑盒回收目的基因，回收目的產(chǎn)物連接到克隆載體pTOPO-Blunt 上，轉(zhuǎn)化大腸桿菌DH52 感受態(tài)細胞，挑選陽性克隆，提取質(zhì)粒進行測序(上海生工生物工程有限公司北京測序部)。

1.4 生物信息學分析

使用DNAMAN 軟件對該基因編碼蛋白進行氨基酸序列的分析，并將氨基酸序列提交NCBI，利用BLAST 進行序列比對，分析與其他物種的相似性及其功能保守結(jié)構(gòu)域，利用MEGA-X 鄰近法構(gòu)建進化樹，將自展循環(huán)次數(shù)(No.of Bootstrap Replications)設為1 000； 將目的氨基酸序列提交ExPASy 網(wǎng)站(https:/ /web.expasy.org/protparam)進行在線分析預測該蛋白的理化性質(zhì)；在SOPMA 和Phyre2 網(wǎng)站對目的蛋白的二級結(jié)構(gòu)和三級結(jié)構(gòu)進行在線分析[20]。

1.5 NAD+-SDH 基因在杏果實發(fā)育過程中的表達分析

分別提取金太陽杏果實5 個發(fā)育時期的總RNA，反轉(zhuǎn)錄合成cDNA 備用。根據(jù)從金太陽杏中擴增的NAD＋-SDH的開放閱讀框架(open read frame，ORF)序列設計特異性實時熒光定量PCR(quantitative realtime PCR，RT-qPCR)引物，如下所示:

PaNAD＋-SDH-F(5′-ACATTGTTCTTAGTGGTGGGT C-3′)；

PaNAD＋-SDH-R(5′-GCTCACATCTACTCCACCTTT C-3′)。

以Actin作為內(nèi)參基因[21]，如下所示:

ACTIN-F(5′-ACATTGTTCTTAGTGGTGGGTC-3′)；

ACTIN-R(5′-AGATTCGTCATACTCTGCCTTT-3′)。

RT-qPCR 采用Trans Start Top Green qPCR Super Mix 熒光定量試劑盒(全式金生物技術有限公司，北京)，使用20 μL 體系于7500 型PCR 儀(Applied Biosystems，美國)檢測NAD＋-SDH基因表達水平。反應程序:94℃預變性30 s；94℃變性5 s，58℃退火15 s，35 個循環(huán)。每樣品重復3 次，采用2-ΔΔCt的方法分析相對表達量。

1.6 SDH 酶活性測定

NAD＋-SDH 酶液的提取及活性測定參考梁東等[22]的方法，所有操作均在4℃下進行。

1.7 數(shù)據(jù)統(tǒng)計與分析

數(shù)據(jù)采用Microsoft Office Excel 2010 和SPSS 18.0軟件進行統(tǒng)計分析。

2 結(jié)果與分析

2.1 PaNAD+-SDH 基因的克隆及序列分析

以金太陽杏果實cDNA 為模板，通過PCR 擴增，電泳后得到一條約1 100 bp 的單一條帶(圖1)，與預期長度一致。測序結(jié)果顯示，該擴增片段包含一個長度為1 104 bp 的ORF，編碼367 個氨基酸(圖2)。比對結(jié)果顯示，該基因堿基序列與薔薇科李屬梅、桃和李的NAD＋-SDH基因cDNA 序列相似性達98%以上，表明該基因為杏的NAD＋-SDH基因，命名為PaSDH。

2.2 PaSDH 編碼蛋白理化性質(zhì)與結(jié)構(gòu)的預測分析

圖1 NAD+-SDH 基因ORF 擴增Fig.1 The ORF amplification of NAD+-SDH gene

使用在線工具ExPASy 分析顯示該基因編碼蛋白的相對分子量為39.19 kDa，理論等電點為6.46，不穩(wěn)定系數(shù)為25.34，屬于穩(wěn)定蛋白，其總平均親水性為0.108，屬于疏水性蛋白。InterProScan 分析表明，PaSDH 蛋白具有含鋅乙醇脫氫酶特征，在序列中具有催化鋅結(jié)合位點、結(jié)構(gòu)鋅結(jié)合位點和NAD 結(jié)合位點等功能位點(圖2)。Signal IP3.0 預測PaSDH 蛋白氨基酸序列沒有信號肽。利用SOPMA 對PaSDH 蛋白進行二級結(jié)構(gòu)分析，其含有30.52%的α-螺旋、6.54%的β-轉(zhuǎn)角、23.16%的延伸鏈和39.78%的無規(guī)則卷曲(圖3)。利用Phyre2 對其三級結(jié)構(gòu)預測，以c1pl6A 山梨醇脫氫酶的三維結(jié)構(gòu)為模板，對349 個殘基(95%的序列)建模，置信度為100%，結(jié)果如圖4 所示。

圖2 NAD+-SDH 基因的核苷酸序列及其編碼氨基酸序列Fig.2 The nucleotide acid sequence and amino acid sequence of NAD+-SDH gene

圖3 PaSDH 基因編碼蛋白的二級結(jié)構(gòu)Fig.3 The secondary structure prediction of PaSDH gene encoding protein

圖4 PaSDH 基因編碼蛋白的三級結(jié)構(gòu)預測Fig.4 Tertiary structure prediction of the PaSDH

2.3 PaSDH 氨基酸多序列比對及系統(tǒng)進化樹構(gòu)建

將PaSDH 氨基酸序列與GeneBank 中登錄的薔薇科植物NAD＋-SDH 氨基酸序列進行同源性比較，采用MEGA-X 軟件中的鄰接法構(gòu)建進化樹(圖5)。從系統(tǒng)進化關系來看，杏PaSDH 氨基酸序列與薔薇科李屬植物梅(Prunus mume，XP ＿008234305.1)、桃(Prunus persica，XP ＿ 007218159.1)、 李(Prunus salicina，ACL18054.1 )、 甜 櫻 桃 (Prunus avium，XP ＿021826929.1)、酸櫻桃(Prunus cerasus，AAK71492.1)和日本櫻花(Prunus yedoensis，PQM34387.1)的NAD＋-SDH 氨基酸序列相似性較高，分別為99.46%、98.91%、98.64%、97.82%、96.20%和94.55%，首先聚為一類。

2.4 杏果實發(fā)育過程中的SDH 酶活性變化和PaSDH 基因的表達分析

分別在杏果實花后7、21、35、49 和56 d 測定SDH酶活性(圖6)，結(jié)果表明，在杏果實發(fā)育過程中，SDH酶活性呈現(xiàn)“高-低-高-低”的變化趨勢，花后7 d，果實中SDH 酶活性最高，花后21 d，SDH 酶活性顯著降低，在花后35 d，SDH 酶活性出現(xiàn)第二個小高峰，但明顯低于花后7 d，之后SDH 酶活逐漸降低。

利用RT-qPCR 對PaSDH在杏果實花后7、21、35、49 和56 d 的相對表達量進行分析。結(jié)果表明，花后7 d，PaSDH表達量最高，之后隨著果實的發(fā)育成熟，PaSH相對表達量逐漸降低，至花后56 d 時，PaSDH相對表達量最低。

3 討論

NAD＋-SDH 酶催化山梨醇和果糖之間的相互轉(zhuǎn)化，其對果糖的米氏常數(shù)Km 值高于對山梨醇的Km值，表明該酶更利于向形成果糖的方向起催化作用，是薔薇科植物中山梨醇代謝的關鍵酶[5]。植物通過山梨醇脫氫酶活性的變化調(diào)節(jié)山梨醇含量以抵御非生物脅迫(干旱脅迫、鹽脅迫)[23]。本研究成功克隆了杏PaSDH基因的全長cDNA，該基因編碼區(qū)長為1 104 bp，可編碼367 個氨基酸。生物信息學分析發(fā)現(xiàn)PaSDH 蛋白具有含鋅醇脫氫酶特征性氨基酸序列、結(jié)構(gòu)鋅結(jié)合位點和催化鋅結(jié)合位點(圖2)，表明PaSDH編碼的蛋白屬于含鋅乙醇脫氫酶家族。這與前人在蘋果[24]和草莓[25]中的研究結(jié)果相似。本研究中，NAD＋-SDH 進化分析表明，PaSDH編碼氨基酸與梅NAD＋-SDH 的氨基酸序列相似性最高，為99.46%，其與薔薇科李屬桃、李、甜櫻桃、酸櫻桃和日本櫻花的SDH 氨基酸序列相似性較高，達到98.91%～94.55%，說明在薔薇科植物中NAD＋-SDH 的氨基酸序列是高度保守的。蘋果和梨中NAD＋-SDH基因以基因家族的形式存在，在蘋果中已經(jīng)分離得到超過10 個編碼NAD＋-SDH基因的全長cDNA[7-8，26]，在梨中也克隆得到10 個編碼NAD＋-SDH基因片段[9，27]。本研究得到1個編碼NAD＋-SDH基因的全長cDNA，隨著最新杏基因組的公布[12]，還將對杏NAD＋-SDH 家族其他成員進行深入探究。

圖5 杏與其他薔薇科植物NAD+-SDH 氨基酸序列的系統(tǒng)進化樹Fig.5 Phylogenetic tree of NAD+-SDH from apricot and other Rosaceae species

通過RT-qPCR 研究杏果實發(fā)育過程中NAD＋-SDH基因的表達，并檢測了酶活性的變化模式。本研究發(fā)現(xiàn)在花后7 d 杏果實中NAD＋-SDH 的酶活性及PaSDH基因的相對表達量最高。因此，認為可能在杏坐果和果實早期的發(fā)育階段，PaSDH起著重要作用，這與梁東等[22]在蘋果上的研究結(jié)果一致。前人在多種薔薇科植物中對SDH 酶活性、蛋白與轉(zhuǎn)錄之間的關系進行了研究。Yamada 等[10，28]發(fā)現(xiàn)在蘋果和桃果實生長過程中NAD＋-SDH基因相對表達量和蛋白含量變化趨勢與酶活性變化很相似，Li 等[29]發(fā)現(xiàn)枇杷果實中NAD＋-SDH 活性、蛋白含量與轉(zhuǎn)錄量的積累密切相關，Kim等[30]在李果實中也發(fā)現(xiàn)相同現(xiàn)象。以上研究結(jié)果均說明NAD＋-SDH基因轉(zhuǎn)錄水平的調(diào)控可能是調(diào)節(jié)其酶活性的關鍵步驟。與之不同的是，梨葉片中的NAD＋-SDH 酶活性在幼葉中較高，而在成熟葉中較低，基因表達量與酶活性變化呈相反趨勢[31]；在蘋果葉片和果實以及梨莖尖和花器官中，SDH 酶活性變化與基因表達模式也存在不一致的情況，說明SDH基因表達可能存在轉(zhuǎn)錄后調(diào)控機制。此外，這種不一致的情況也可能是基因家族成員共同作用的結(jié)果[32]。本研究中，隨著杏果實的生長，NAD＋-SDH 酶活性與基因相對表達量總體呈現(xiàn)下降趨勢，這與杏果實成熟過程中果糖含量逐漸下降的變化趨勢一致[16]。但在花后35 d，酶活性出現(xiàn)小幅度升高，此時基因的相對表達量卻依然呈下降趨勢，同樣觀察到酶活性與基因的轉(zhuǎn)錄水平的變化趨勢并不完全一致的現(xiàn)象，很可能也是由PaSDH 酶活性變化與SDH 基因家族成員的表達模式不一致造成的，說明在杏果實發(fā)育早期NAD＋-SDH 酶活性有可能受PaSDH基因轉(zhuǎn)錄的調(diào)控。此外，酶活性與基因轉(zhuǎn)錄水平變化趨勢不完全一致說明可能還存在其他層次的調(diào)控機制。

4 結(jié)論

本研究克隆的杏PaSDH基因cDNA 全長1 104 bp，編碼367 個氨基酸，蛋白質(zhì)分子量為39.19 kDa。PaSDH基因相對表達水平和PaSDH 酶活性在杏果實發(fā)育早期最高，在果實成熟過程中呈總體下降趨勢。但在花后35 d 存在酶活性和基因表達水平變化不一致的現(xiàn)象，說明可能存在SDH基因家族成員表達不一致的模式，或是存在除轉(zhuǎn)錄水平調(diào)控外的其他層次的調(diào)控機制。本研究為解析杏SDH基因表達模式及調(diào)控機制奠定了分子生物學基礎。