張玉葉等

摘 要 應用電子克隆技術，以煙草AAG59585序列為探針，獲得了甘蔗蘇氨酸脫氨酶（threonine deaminase）基因的全長cDNA，命名為ScTD。經RT-PCR擴增和驗證，該基因序列與電子克隆結果一致性高達99%。序列分析結果表明：ScTD基因全長2 120 bp，具有完整的開放閱讀框（ORF，164～1 681 bp），編碼505個氨基酸，該基因編碼蛋白定位于微體，為可溶性蛋白，不存在信號肽，二級結構原件多為α-螺旋，含有多個保守功能域，主要在氨基酸合成中發(fā)揮作用。電子表達分析結果顯示，該基因在甘蔗根尖、根頸、花序、葉片、全株、芽和莖中組成型表達，其中在花序和根中的表達量最高。熒光定量PCR結果分析表明：該基因在甘蔗各組織中均有表達，且在根中的表達量最高。此外，該基因的表達受水楊酸脅迫后表達量最高，約為對照組的43.16倍，其次為聚乙二醇和茉莉酸甲酯，分別為6.90倍和4.40倍，推測該基因的表達與甘蔗抗病蟲和抗?jié)B透脅迫有關。此結果為甘蔗ScTD基因的克隆和功能驗證以及在甘蔗基因工程中的應用奠定基礎。

關鍵詞 甘蔗；蘇氨酸脫氨酶基因；電子克??；生物信息學；表達分析

中國分類號 S566.1 文獻標識碼 A

Cloning and Expression Analysis of Threonine

Deaminase Gene in Sugarcane

ZHANG Yuye， HUANG Ning， SU Weihua， XIAO Xinhuan， LUO Jun， QUE Youxiong*

Key Laboratory of Sugarcane Biology and Genetic Breeding， Ministry of Agriculture， Fujian Agriculture and Forestry University/

Sugarcane Research & Development Center， China Agriculture Research System， Fuzhou， Fujian 350002， China

Abstract Using AAG59585 sequence from Nicotiana attenuata as the probe， the full-length cDNA sequence of one threonine deaminase gene termed as ScTD was cloned in silico from sugarcane（Saccharum Complex）. The sequence obtained from RT-PCR amplification had as high as 99% similarity with that from electronic cloning. The sequence analysis showed that the ScTD had a length of 2 120 bp containing an open reading frame（ORF， 164～1 681 bp）encoding 505 amino acids residues. The characters predicted based on the analysis of bioinformatics showed that the ScTD of sugarcane was a soluble acidic protein， which had two conserved functional domains with the main function for amino acid synthesis， and this protein was located in microbody. The main secondary structure element was α-helix. The electron expression analysis suggested that ScTD was constitutively expressed in different tissues of sugarcane， such as root tip， root collar， inflorescence， leaf， buds， whole plant and stem， especially higher in inflorescence and root. Real-time quantitative PCR（RT-qPCR）analysis revealed that the expression of ScTD in root was higher than that in other organs， while the inducible expression level of ScTD was apparent under the stresses of salicylic acid（SA）， polyethylene glycol（PEG）and methyl jasmonate（MeJA）， approximately 43.16 times， 6.90 times and 4.40 times that of control group respectively， which suggested that ScTD involves in the diseases and pest resistance and osmotic stresses in sugarcane. The results in this study will provide the basis of cloning and functional identification of other members of ScTD in sugarcane and promote the use of ScTD gene in sugarcane genetic engineering.

Key words Sugarcane； Threonine deaminase gene； In silico cloning； Bioinformatics； Expression analysis

doi 10.3969/j.issn.1000-2561.2014.01.011

蘇氨酸脫氨酶（threonine deaminase，TD）也叫蘇氨酸脫氫酶（threonine dehydratase）。在植物和細菌中生物合成的蘇氨酸脫氨酶包含一個N端磷酸吡哆醛結合催化域和C端調控域，該酶在原核細胞及真核細胞中均可能與甘氨酸C-乙酰基轉移酶（glycine C-acetyltransferase）結合成為復合物協同作用，將蘇氨酸轉化生成α-酮丁酸和氨，然后以α-酮丁酸為反應物合成異亮氨酸，此途徑即為異亮氨酸合成的關鍵步驟[1]。L-異亮氨酸是合成人體激素、酶類的原料，具有促進蛋白質合成和抑制其分解的效果，在人體生命活動中起著重要作用，在食品和醫(yī)藥行業(yè)中具有廣泛的應用及商業(yè)價值[2]。蘇氨酸脫氨酶受異亮氨酸的負反饋調控[3]。目前，該酶在大腸桿菌、鼠傷寒桿菌、谷氨酸棒桿菌等[4-5]細菌中已有報道。在植物中，蘇氨酸脫氨酶與植物抗蟲有關，分別在番茄等茄科植物[6]，煙草[7]中有相關報道，但在甘蔗中尚未見有關蘇氨酸脫氨酶的報道。

電子克?。╥n silico cloning）是通過EST或基因組的序列組裝和拼接，快速地獲得目標基因部分或全長cDNA序列的技術。該技術在人、鼠和水稻等[8-9]EST豐富的模式生物中首先得到運用，隨后在植物中運用此技術獲得新基因的報道日益增加，此技術在甘蔗中也得到了有效利用[10-11]。電子表達分析是通過某物種特定基因所有相關表達序列標簽（Expressed Sequence Tag，EST）的整合分析，從而獲得該基因表達相關信息的一種新型基因表達分析技術。目前，電子表達分析與電子克隆技術結合起來挖掘新基因和分析基因功能已被成功應用于包括甘蔗等多種植物目標基因的表達分析[12]。本研究以甘蔗為材料，電子克隆與實驗驗證結合獲得甘蔗蘇氨酸脫氨酶基因，并對其功能進行探討，為異亮氨酸的生產以及研究該基因在甘蔗抗逆性中的作用提供基礎數據。

1 材料與方法

1.1 材料

1.1.1 植物材料 甘蔗品種新臺糖22號（ROC22）由福建農林大學農業(yè)部福建甘蔗生物學與遺傳育種重點實驗室提供。

1.1.2 試劑 PrimeScriptTM RT-PCR Kit反轉錄試劑盒、SYBR R Green PCR Master Mix Kit （Roche）均購自大連寶生物（TaKaRa）公司。

1.2 方法

1.2.1 植物材料處理 從田間選取3株健康、長勢一致的ROC22植株，整株取回。將根部的泥土洗凈，選取白色的幼根作為材料。用于組織表達的材料為甘蔗第七或第八節(jié)，分別取蔗芽，蔗皮，蔗髓，第一片葉的葉片和葉鞘。研究目的基因受不同外源脅迫后的表達特性時，將ROC22蔗莖作為種莖，砍成單芽莖段后，采用沙培法，在塑料大棚中育苗至其長出第四或第五葉。選取長勢一致的蔗苗水培1周，分為對照和4種脅迫組，每組每個時間點均設3株蔗苗，作為生物學重復。脅迫組處理條件為水楊酸（salicylic acid，SA）（5 mmol/L）、茉莉酸甲酯（methyl jasmonate，MeJA）（100 μmol/L）、脫落酸（abscisic acid，ABA）（100 μmol/L）、模擬干旱（PEG8000）（25.0%），每1 h噴撒葉面1次，未處理的蔗葉作為對照。分別取6 h后各小組內3棵蔗苗的蔗葉。以上所有植物材料均立即投入-80 ℃液氮中保存，備用。

1.2.2 電子克隆獲得基因序列 以煙草TD基因（AAG59585）蛋白序列作為探針，運行NCBI的tblastn檢索甘蔗EST數據庫，獲得與之相似的甘蔗EST序列群，用BioEdit中的contig assembly program（CAP）（http：//pbil. univ-lyon1.fr/cap3.php）軟件拼接組裝為重疊群（contig），獲得組裝全長cDNA序列的EST候選清單。然后以最終拼好的contig重疊群為新探針，再次進行搜索，直到沒有新的甘蔗EST可供拼接為止。最后利用NCBI的ORF finder（http：//www.ncbi.nlm.nih.gov/gorf/gorf.html）服務器，對最終所獲得的cDNA序列進行開放閱讀框分析并進行翻譯。

1.2.3 甘蔗TD基因的RT-PCR及測序 采用Trizol法提取甘蔗ROC22葉片的總RNA。用反轉錄試劑盒合成cDNA，作為PCR的模板。以上述電子克隆獲得的cDNA 序列為模板，設計一對特異性RT-PCR擴增引物，5′端引物（5′-TGCCCGTGACAA

CGCCAGAA-3′）和3′端引物（5′-GGTGCCGCAAACC

GATAAA-3′）。PCR反應體系為25 μL：10×GC Buffer 2.5 μL，10 mmol/L dNTPs 2 μL，引物（20 μmol/L）各1.0 μL，Ex Taq酶0.2 μL，cDNA模板1 μL，ddH2O 17.3 μL。PCR反應程序：94 ℃ 4 min；94 ℃ 1 min，60 ℃ 45 s，72 ℃ 1 min，35個循環(huán)；72 ℃ 10 min。用Gel Extraction Kit對擴增產物進行純化回收，并與pMD-19T載體連接轉化至DH5α感受態(tài)細胞中，在抗性Ampr的LB平板上篩選陽性克隆，隨機挑取單菌落進行菌落PCR鑒定，將陽性單菌落菌液送至上海生工生物工程技術服務有限公司測序。

1.2.4 甘蔗ScTD基因序列的生物信息學分析

利用ExPASy中ProtParam pI/Mw（http：//web.expasy.org/compute_pi/）在線工具對序列的理化性質、氨基酸組成等一級結構進行預測，利用SignalP 4.1 Server（http：//www.cbs.dtu.dk/services/SignalP/）進行信號肽分析，ProtScale（http：//web.expasy.org/protscale/）進行疏水性/親水性預測，Profun 2.2 Server（http：//www.cbs.dtu.dk/services/ProtFun/）進行功能預測，編碼蛋白亞細胞定位用Psort（http：//www.psort.org/）在線工具預測，GOR IV（http：//npsa-pbil.ibcp.fr/cgi-bin/npsa_automat.pl？page=npsa_gor4.html）二級結構預測，SWISSMODEL（http：//swissmodel.expasy.org/）進行蛋白質三級結構預測，Blast上查找氨基酸同源性序列。

1.2.5 甘蔗ScTD基因的電子表達分析 根據植物EST序列組織來源或脅迫材料的來源不同，進行電子表達分析，分析該基因組織特異性表達或生物、非生物脅迫后表達特性[12]。

1.2.6 甘蔗ScTD基因的熒光定量PCR分析 采用Trizol 法提取ROC22葉片、葉鞘、側芽、蔗皮、蔗髓、根及經SA、MeJA、ABA、PEG脅迫處理和對照組ROC22葉片的總RNA，用反轉錄試劑盒合成cDNA，作為熒光定量PCR模板，以25 S rRNA作為內參基。合成熒光定量引物5′端引物（5′-TGC

AATGGCATTATCCTTGT-3′）和3′端引物（5′-TGACA

GCTACACCATCAGCA-3′），25S rRNA的5′端引物（5′-TGCTTTGCTCGGTTTATAGC-3′）和3'端引物（5′-GCAGCCAAGCGTTCATAGC-3′）[13]。按照SYBR R Green PCR Master Mix Kit（Roche）說明書配置定量反應體系。熒光定量PCR擴增程序為：50 ℃ 2 min；95 ℃ 10 min；95 ℃ 15 s，60 ℃ 1 min，45個循環(huán)，設3次重復。

1.3 數據處理

使用DNAMAN軟件，將重組質粒序列與電子克隆序列進行同源性比對。并用DNAMAN軟件分析氨基酸同源性及多序列比對，完成并構建系統(tǒng)進化樹。在Excel工作表中采用 2-△△CT算法[14-15]分析熒光定量PCR實驗結果，畫圖并計算3次重復數據的標準誤。

2 結果與分析

2.1 甘蔗ScTD基因的獲得

電子克隆得到一個全長為2 120 bp的基因序列（圖1），該cDNA序列具有完整的開放閱讀框（ORF，164～1 681 bp），編碼505個氨基酸，說明此序列為一條完整的甘蔗蘇氨酸脫氨酶基因的cDNA序列，命名為ScTD。

經過RT-PCR 擴增獲得了約1 400 bp的單一條帶（圖2），經克隆與測序，證實其序列與電子克隆所獲得的序列一致性達到99%，僅有14個堿基不同，經推測可能是在PCR擴增時堿基錯配造成的。因此，推測電子克隆所獲得的甘蔗ScTD基因序列是正確的。

2.2 甘蔗ScTD基因的生物信息學分析

2.2.1 甘蔗ScTD基因編碼氨基酸的一級結構預測

甘蔗ScTD基因編碼氨基酸進行一級結構預測，結果見表1。

2.2.2 甘蔗ScTD蛋白二級結構預測和分析 結果見表2。表2顯示，α-螺旋所占的比例最高為45.64%，其次是無規(guī)則卷曲占37.15%和最低的延伸鏈占比例17.19%。

2.2.3 甘蔗ScTD蛋白信號肽預測和分析 結果見表3所示。第48位苯丙氨酸殘基具有最高的原始剪切位點分值0.126和最高的信號肽分值0.109，第24 位丙氨酸殘基具有最高的綜合剪切位點分值0.131。由于最后算得氨基酸殘基的加權平均值較小，為0.103（遠遠<0.5），則推測ScTD基因所編碼的蛋白不存在信號肽，而為非分泌蛋白，說明該蛋白在細胞質中合成后，不能進行蛋白轉運。

2.2.4 甘蔗ScTD蛋白疏水性/親水性的預測和分析

由圖3可知，第185位具有最高分值，為2.300，疏水性最強；第361位和第499位具有最低分值，為-2.322，親水性最強。整條多肽鏈表現為親水性，沒有明顯的疏水區(qū)，故推測甘蔗ScTD蛋白是一種親水蛋白。

2.2.5 甘蔗ScTD蛋白質三級結構預測 ScTD蛋白的空間結構以無規(guī)則卷曲和螺旋為主（圖4）。比較甘蔗、高粱、玉米和水稻中TD蛋白的三級空間結構虛擬圖，結果顯示，甘蔗的三級空間結構與高粱、玉米和水稻蛋白的三級空間結構高度相似。 2.2.6 甘蔗ScTD蛋白的功能預測 該蛋白的最主要功能為氨基酸生物合成，也有一定可能性參與了輔酶因子生物合成、翻譯能量代謝及脂肪酸代謝（表4）。由于蘇氨酸脫氨酶是異亮氨酸合成的關鍵酶之一，則預測結果與實際相符。

2.2.7 甘蔗ScTD蛋白的亞細胞定位 甘蔗ScTD蛋白可能定位于微體的過氧化酶體中（表5）。

2.2.8 甘蔗ScTD蛋白的保守結構域分析 如圖5所示，甘蔗ScTD蛋白包含2個保守結構域，分別為Trp-synth-bete_ⅡSuperfamily和ACT Superfamily。這和前文所提到此基因在植物和細菌中N端和C端的結構特征相符[1]。

2.2.9 甘蔗ScTD蛋白的氨基酸同源性分析 對甘蔗ScTD蛋白的氨基酸進行同源性分析。在結果中選擇高粱（Sorghum bicolor_gb|XP_002466574.1|）、玉米（Zea mays_gb|DAA51003.1|）、水稻（Oryza sativa Japonica Group_gb|AAL 58211.1|）、二穗短柄草（Brachypodium distachyon_gb|XP_003559531.1|）、番茄（Solanun lycopersicum_gb|NP_ 001234234.1|）、毛楊果（Populus trichocarpa_gb|XP_002326230.1|）和黃瓜（Cucumis sativus_gb|XP_004143291.1|）進行分析，蛋白的氨基酸序列相似性分別為98%、96%、94%、88%、80%、80%和76%。

從圖6可以看出，同屬單子葉植物的甘蔗和高粱同源性最高。系統(tǒng)進化樹顯示（圖7）：單子葉植物甘蔗、玉米、水稻和二穗短柄草為一分支，同屬雙子葉植物的番茄、黃瓜和毛楊果為另一分支。因此，該基因在不同物種間具有很強的保守性，它們對應的蛋白屬于同一個家族。

2.3 甘蔗ScTD基因的表達分析

2.3.1 甘蔗ScTD基因的電子表達分析 電子表達分析結果見圖8。圖8顯示，甘蔗ScTD基因在甘蔗根尖、根頸、花序、葉片、莖、全株和芽中組成型表達，其中在花序中的表達量比其他類型組織表達量高。此外，該基因的表達受葡萄糖桿菌調控。

2.3.2 甘蔗ScTD基因的熒光定量PCR表達分析

分析甘蔗ScTD基因在不同組織中的表達，結果見圖9所示。由圖9可知，該基因在不同組織中均有表達，在蔗髓中表達量最少；其次是芽中；根中表達量最高，是蔗髓中的51.02倍，是芽中表達量的10.87倍。分析甘蔗ScTD基因在不同脅迫下的應答反應，結果見圖10。由圖10可知，該基因在不同脅迫下的相對表達量呈上調趨勢，其中在水楊酸脅迫下的相對表達量最高，為對照組的43.16倍，其次是在聚二乙醇和茉莉酸甲酯脅迫下的表達量，分別為6.90和4.40倍。推測該基因表達與甘蔗的抗病性和滲透有關，與前人的報道相符[6-7]。

3 討論與結論

本研究采用電子克隆與實驗克隆相結合的方法，成功地克隆了甘蔗ScTD基因，重組子序列與電子克隆序列同源性高達99%，推測本研究電子克隆所獲得的序列是正確的，進一步證明電子克隆與實驗克隆相結合是發(fā)現新基因的有效途徑。而后，采用多種生物信息學手段分析了該基因的特性，甘蔗ScTD基因在不同物種間具有較強的保守性，不同物種相應的蛋白同源序列均達到80%左右。所編碼的蛋白位于微體，為親水、非分泌酸性蛋白，在調控異亮氨酸的生物合成中發(fā)揮關鍵作用。

本研究利用Unigene數據庫分析了甘蔗蘇氨酸脫氨酶基因的電子表達，檢索到22條同源性較高的EST序列，甘蔗EST數據庫中包含了283 951條EST，已經能夠非常好地覆蓋了甘蔗的全基因組，此結果能夠如實反映該基因在甘蔗基因組中的表達情況。電子表達分析結果顯示，該基因在花序和根中都有很高的表達量，這與熒光定量PCR結果大體一致，都是在根中相對表達量高，說明電子表達分析技術可以用于甘蔗中目的基因的表達分析，這與本課題組以往的的研究結果一致[12]。

水楊酸（SA）在超敏反應（hypersensitive response，HR）中能介導導致細胞死亡的氧化迸發(fā)（oxidative burst），誘導病程相關蛋白（pathogenesis-related protein，PR）基因的表達，引發(fā)系統(tǒng)獲得性抗性（systemic acquired resistance，SAR）的產生，并且能誘導一些防衛(wèi)基因（如PR）的激活或表達[16]。PEG為一種理想的模擬滲透脅迫物質[17]。茉莉酸和茉莉酸甲酯能夠誘導植物產生對害蟲有毒、抗營養(yǎng)和抗消化作用的化合物，對害蟲生理活動產生不利影響[18]。本文通過熒光定量PCR實驗對甘蔗ScTD基因在不同組織和不同脅迫下的相對表達量行進分析。不同組織中的表達分析結果顯示，該基因在蔗髓、蔗皮、側芽、葉片、葉鞘和根中都有表達，其中在根中的相對表達量最高，芽次之。在SA、ABA、MeJA和PEG等脅迫下的表達分析結果顯示，該基因在本實驗所有脅迫下均上調表達，其中SA脅迫下的相對表達量最高，PEG和MeJA脅迫下的表達量次之。有研究結果表明，氨基酸合成所需的氮素主要來源于大氣和土壤，根吸收土壤中的硝酸鹽，以此來進行氮代謝[19]，而蘇氨酸脫氨酶是合成異亮氨酸關鍵酶之一，主要作用為氨基酸的合成，這可能是該基因在根中表達量會相對其他組織要高的原因。芽的生長發(fā)育需要較多氨基酸，推測這是蘇氨酸脫氨酶在芽中相對含量較高的原因。因此，本研究推測蘇氨酸脫氨酶基因與甘蔗的抗病蟲和抗?jié)B透脅迫存在一定的關聯性。

綜上所述，本研究采用電子克隆與實驗驗證的方法成功克隆到甘蔗ScTD基因，并進行了該基因的生物信息學分析、電子表達分析和熒光定量PCR分析。研究結果為甘蔗中ScTD基因功能的深入解析及其應用積累了基礎資料。

參考文獻

[1] Hartshorne D， Greenberg D M. Studies on liver threonine dehydrogenase[J]. Archives of Biochemistry and Biophysics， 1964， 105： 173-178.

[2] Blomstrand E. A role for branched-chain amino acids in reducing central fatigue[J]. Journal of Nutrition， 2006， 136（2）： 544-547.

[3] Gallagher D T， Gilliland G L， Xiao G Y， et al. Structure and control of pyridoxal phosphate dependent allosteric threonine deaminase[J]. Structure， 1998， 6： 465-475.

[4] Eisenstein E. Cloning， expression， purification， and characterization of biosynthetic threonine deaminase from Escherichia coli[J]. Journal of biological chemistry， 1991， 266（9）： 5 801-5 807.

[5] 謝希賢， 杜 軍， 徐慶陽， 等. 谷氨酸棒桿菌蘇氨酸脫氨酶基因敲除及對L-亮氨酸發(fā)酵的影響[J]. 中國食品學報， 2010， 10（2）： 40-45.

[6] Gonzales-Vigil E， Bianchetti C M， Phillips G N， et al. Adaptive evolution of threonine deaminase in plant defense against insect herbivores[J]. Proceedings of the National Academy of Science， 2011， 108（14）： 5 897-5 902.

[7] Kang J H， Wang L， Giri A， et al. Silencing threonine deaminase and JAR4 in Nicotiana attenuata impairs jasmonic acid-isoleucine-mediated defenses against Manduca sexta[J]. Plant Cell， 2006， 18： 3 303-3 320.

[8] 王海濤， 暢繼武， 馬小兵， 等. 人和大鼠自噬相關新基因WIPI-3的電子克隆和序列分析[J]. 現代生物醫(yī)學進展， 2006， 6（4）： 16-21.

[9] 黃 驥， 王建飛， 張紅生， 等. 水稻葡萄糖-6-磷酸脫氫酶cDNA的電子克隆[J]. 遺傳學報， 2002， 29（11）： 1 012-1 016.

[10] 黃 寧， 張玉葉， 蘇亞春， 等. 甘蔗天冬氨酰半醛脫氫酶基因的電子克隆與分析[J]. 生物信息學， 2013， 11（2）： 91-99.

[11] 李國印， 闕友雄， 許莉萍， 等. 甘蔗MYB2轉錄因子的電子克隆和生物信息學分析[J]. 生物信息學， 2011， 09（1）： 24-27.

[12] 黃 寧， 張玉葉， 凌 輝， 等. 甘蔗二氨基庚二酸異構酶基因的克隆與表達分析[J]. 熱帶作物學報， 2013， 34（11）： 2 200-2 209.

[13] 闕友雄， 許莉萍， 徐景升， 等. 甘蔗基因表達定量PCR分析中內參基因的選擇[J]. 熱帶作物學報， 2009， 30（3）： 274-277.

[14] Livak K J， Schmittgen T D. Analysis of relative gene expression data using Real-Time quantitative PCR and the 2- ΔΔCT method[J]. Methods， 2001， 25（4）： 402-408.

[15] 闕友雄， 許莉萍， 林劍偉， 等. 甘蔗NBS-LRR類抗病基因同源序列的分離與鑒定[J]. 作物學報， 2009， 35（4）： 631-639.

[16] Rasin I. Role of salicylic acid in plants[J]. Annual Review of Plant Physiology and Plant Molecular Biology， 1992， 43： 439-463.

[17] Jia G L， Dai H P， Feng B L， et al. Biochemical characteristics in Broomcorn Millet（Panicum miliaceum L）seedlings under PEG simulated drought stress[J]. Acta Botanica Boreali-Occidentalia Sinica， 2008， 28： 2 073-2 079.

[18] 桂連友， 劉樹生， 陳宗懋. 外源茉莉酸和茉莉酸甲酯誘導植物抗蟲作用及其機理[J]. 昆蟲學報， 2004， 47（4）： 507-514.

[19] 許振柱， 周廣勝. 植物氮代謝及其環(huán)境調節(jié)研究進展[J]. 應用生態(tài)學報， 2004， 15（3）： 511-516.

責任編輯：趙軍明