馬鈴薯ERECTA 基因克隆和生物信息學分析

張 樂，陳如意，李玉雙，郭志平

（麗水學院，a.醫學院；b.生態學院；c.醫務室，浙江 麗水 323000）

ERECTA 是一種富含亮氨酸重復序列的類受體蛋白激酶，廣泛分布于各種植物的器官和組織［1］。擬南芥（Arabidopsis thaliana）ERECTA基因參與花藥發育，調控細胞增殖和分化［2］，ERECTA基因還與植株抗逆性相關，其過表達能夠提高植株對高溫干旱的耐受性［3，4］。目前，僅在少數經濟作物上有EREC?TA基因相關的研究［5-7］。

馬鈴薯（Solanum tuberosum）在世界各地廣泛種植，因其具有產量高、適應性強、營養豐富等特點，已成為世界第三大淀粉原料作物和第四大糧食作物［8］。馬鈴薯除菜用外，近年來，利用馬鈴薯生產薯片、薯條、淀粉等的產量不斷增加。為滿足不同需求，育種工作者根據不同區域的氣候特點和生產用途，在利用傳統方法積極開展馬鈴薯種質資源研究和雜交育種的同時，積極探索利用現代生物技術在分子水平開展種質資源研究和分子標記育種，目前已取得了一定的成果［9］。本研究通過RT-PCR 技術克隆馬鈴薯ERECTA基因cDNA 序列，并對其進行生物信息學分析，為進一步開展該基因的功能研究和馬鈴薯分子育種提供參考。

1 材料與方法

1.1 材料

供試材料為浙江省麗水學院生態學院實驗室保存的克新19 馬鈴薯（種質來自黑龍江省農業科學研究院克山分院）。

1.2 方法

1.2.1 馬鈴薯總RNA 提取與反轉錄 總RNA 的提取參照李飛等［10］的方法，取新鮮馬鈴薯幼苗嫩葉置于預冷的研缽中，同時加入Trizol 試劑，充分研磨并振蕩混勻，放置5 min，離心10 min 后取上清液；添加氯仿混勻，室溫下放置3 min，離心15 min 取上清液；換至新離心管中，按1 mL Trizol 加0.5 mL 異丙醇的比例配比，靜置后離心10 min，棄上清液；按1 mL Trizol 液加入1 mL 以上的75% 乙醇清洗沉淀，渦旋振蕩混勻，離心5 min，沉淀在室溫下干燥5~10 min，然后將總RNA 溶于DEPC 水。 以提取得到的總RNA 為模板，采用PrimeScriptTMRT reagent 試劑盒進行反轉錄。

1.2.2 馬鈴薯ERECTA基因cDNA 克隆 在NCBI數據庫中獲取馬鈴薯ERECTA基因cDNA 序列（NCBI登陸號：XM_006353788），設計上游引物：5′-TCAAAGGAGCTGAGTTAGTGCA-3′；下游引物：5′-TCAGCATCAGTGGTCCCAAA-3′。PCR 擴增程序為94 ℃預變性5 min；94℃變性30 s，55.7 ℃退火30 s，72 ℃延伸3 min，共35個循環，最后72 ℃延伸10 min。瓊脂糖凝膠電泳檢查PCR 擴增產物，割膠回收后送生工生物工程（上海）股份有限公司檢測。

1.2.3 生物信息學分析 利用ORF finder（https：//www.ncbi.nlm.nih.gov/orffinder/）將馬鈴薯ERECTA基因的cDNA 翻譯成氨基酸序列。使用ProtParam tool（https：//web.expasy.org/protparam/）預測等電點、分子質量等理化性質，并通過ProP 1.0 Server（http：//www.cbs.dtu.dk/services/ProP-1.0/）分析馬鈴薯ERECTA蛋白中有無信號肽、前導肽。利用MEGA6.0 構建系統進化樹，1 000 次Bootstrap。 用SWISS-MODLE（https：//swissmodel.expasy.org/interactive）對馬鈴薯ERECTA 蛋白進行三維建模。

2 結果與分析

2.1 馬鈴薯ERECTA 基因cDNA 克隆

利用Trizol 試劑從馬鈴薯幼苗嫩葉中提取到了總RNA，經Nanodrop 和瓊脂糖凝膠電泳確定總RNA濃度和質量（圖1A）。隨后反轉錄獲得cDNA 第一鏈，并利用Taq酶體外擴增ERECTA基因cDNA，電泳檢測條帶大小為3 281 bp，與預期一致（圖1B）。

圖1 馬鈴薯ERECTA 基因cDNA 克隆

2.2 馬鈴薯ERECTA 蛋白分子特征

馬鈴薯ERECTA基因cDNA 序列包括一個2 973 bp 的開放閱讀框，可編碼990 個氨基酸。預測分析馬鈴薯ERECTA 蛋白理化性質，結果顯示該蛋白的分子質量為108.95 kDa，理論等電點為5.85。馬鈴薯ERECTA 蛋白具有一個信號肽，信號肽切割位點位于32 位和33 位氨基酸之間；馬鈴薯ERECTA不含前導肽（圖2）。

圖2 馬鈴薯ERECTA 蛋白信號肽和前導肽預測

2.3 ERECTA 氨基酸序列對比

多序列比對結果表明，馬鈴薯ERECTA 蛋白的氨基酸序列中含有10 個保守基序，其中8 個為富含亮氨酸的LRR 基序，它們由長度為8~45 的氨基酸殘基組成，通常折疊成弧形或馬蹄形；ERECTA 有1 個Pkinase 結構域，該結構域是具有催化功能的蛋白激酶結構域（圖3）。

圖3 馬鈴薯ERECTA 及其同源序列的多重比對

2.4 同源性分析

序列同源性分析表明，馬鈴薯ERECTA與辣椒（Capsicum annuum）ERECTA的 一 致 性 最 高（94.04%），其次是與煙草（Nicotiana tabacum）ERECTA的一致性（93.65%）較高，與黃花風鈴木（Handroanthus impetiginosus）ERECTA 的一致性最低（85.03%）。這與多序列比對的結果一致（表1）。

表1 馬鈴薯ERECTA 與其他植物ERECTA 一致性分析

2.5 ERECTA 蛋白的系統進化樹分析

為了分析馬鈴薯ERECTA 與其他植物ERECTA之間的親緣關系，利用MEGA 軟件構建了系統進化樹。結果表明，馬鈴薯ERECTA 與智利番茄、芝麻ERECTA 聚為一簇，進而和小粒咖啡聚為一簇，說明馬鈴薯ERECTA 在進化上的保守性（圖4）。

圖4 馬鈴薯ERECTA 與其他植物ERECTA 的系統進化樹分析

2.6 馬鈴薯ERECTA 蛋白質結構分析

馬鈴薯ERECTA 的二級結構預測結果表明，該蛋白質含34.9% 的α-螺旋，50.9% 的無規則卷曲，14.2% 的β-折疊。三級結構預測結果表明，馬鈴薯ERECTA 包含8 個富含亮氨酸的LRR 基序（圖5）。

3 討論

ERECTA基因對植株生長發育有多個方面作用，如控制植株的矮化過程［11］。對擬南芥的研究表明，ERETCA 通過調控分生組織的細胞分裂周期實現植株矮化［12，13］。ERECTA基因家族還可以促進細胞增殖和抑制氣孔分化，從而促進植物葉片的伸長［14］。相關研究表明，ERETCA基因與植物抗逆性也密切相關，能夠提高植物對高溫干旱的耐受性［15］。例如高粱（Sorghum bicolor）ERECTA1和ERECTA2基因均在高粱莖和葉中表達，而且其表達水平均隨著干旱脅迫程度的加深而逐漸提高［15］。因此，進一步研究ERECTA基因在馬鈴薯生長發育中的作用，對馬鈴薯種質研究和育種有指導意義。

馬鈴薯ERECTA蛋白分子的質量為108.95 kDa，這與Shen 等［3］研究的擬南芥ERECTA 分子質量有所不同，但與辣椒ERECTA 的分子質量相近。通過多序列比對和系統發育分析表明，馬鈴薯與辣椒等的親緣關系較近，但與擬南芥的關系較遠。研究同時表明，馬鈴薯ERECTA 與擬南芥ERECTA 在結構和功能上有所差異［7］。ERECTA 為編碼富含亮氨酸重復序列的類受體蛋白激酶（LRR-RLK），LRR-RLKs是普遍存在于植物中的信號受體亞家族，它們調控多種信號通路［15］。馬鈴薯ERECTA 蛋白同樣擁有8個LRR 保守基序以及1 個激酶構域，這些保守度較高的基序在其他植物ERECTA 蛋白中也存在［15，16］。馬鈴薯ERECTA基因的克隆和生物信息學分析，為后續的功能研究和分子育種奠定了基礎。