甘蔗TB1基因的克隆與生物信息學分析

李旭娟　林秀琴　劉洪博　李純佳　徐超華　劉新龍

摘 要 以負調控禾本科植物分蘗的關鍵基因TB1為研究對象，根據水稻、玉米、高粱TB1同源基因保守序列設計引物，使用RT-PCR方法和RACE技術從甘蔗莖尖處克隆到該基因的同源基因ScTB1。ScTB1 cDNA全長1 668 bp，包含274 bp的5′ UTR，1 149 bp的CDS和245 bp的3′ UTR。生物信息學分析表明該基因可編碼382個氨基酸殘基，編碼產物含有保守的SP區、TCP區和R區，屬于TCP家族轉錄因子。其分子量為40.7 ku，理論等電點為8.55，蛋白亞細胞定位預測其主要定位于細胞質。系統進化分析表明ScTB1屬于CYC/TB1亞族蛋白，與高粱TB1和玉米TB1聚為一個亞類。根據序列的保守性和預測結果可推測ScTB1很可能也參與了甘蔗分蘗性狀的調控。

關鍵詞 甘蔗；分蘗；TB1基因；克?。簧镄畔W分析

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

Abstract In this study， TB1 homolog of sugarcane was cloned from the stem apex by using RT-PCR and RACE method referring to the conserved sequence of TB1 homologs in rice， maize and sorghum. The cDNA length of ScTB1 is 1 668 bp， which contains 274 bp 5′ UTR， 1 149 bp CDS and 245 bp 3′ UTR. Bioinformatics analysis showed that ScTB1 belonged to the TCP family of transcription factors， which is composed of 382 amino acid residues and contains three conservative domains（SP， TCP and R domains）. The molecular weight and theoretical isoelectric point of ScTB1 was 40.7 ku and 8.55， respectively， and the gene was likely located in cytoplasm. Phylogenetic tree analysis indicates that ScTB1 belonged to CYC/TB1 subfamily and had a very close relationship with TB1 homologs of maize and sorghum. Based on the sequence conservation， it is suggested that ScTB1 may be required for the regulation of tillering trait. This study provides important information for further research on the functional analysis and regulating mechanisms of ScTB1 in future.

Key words Sugarcane； Tillering； TB1 gene； Cloning； Bioinformatics analysis

doi 10.3969/j.issn.1000-2561.2015.11.010

甘蔗（Saccharum spp.）是主產于熱帶和亞熱帶區域的重要糖料經濟作物，全世界近80%以上的食糖來自于甘蔗。此外，甘蔗還因蘊含較高的生物量而成為世界上最重要的能源作物之一[1-2]。甘蔗屬無性繁殖作物，蔗莖為其主要收獲產品，主要由單一的主莖和有效分蘗莖組成，其中有效分蘗莖數量的多少對于甘蔗最終產量的大小具有重要影響，由此可見，分蘗性狀表現對于甘蔗產量具有舉足輕重的作用。而從甘蔗分蘗性狀著手，研究其分子基礎，對于利用分蘗性狀改良甘蔗品種產量具有重要意義。

TB1（TEOSINTE BRANCHED 1）基因，是一個作用于獨腳金內酯（Strigolactones， SLs）下游的基因，屬植物所特有的TCP家族轉錄因子，編碼含有不典型的螺旋-環-螺旋結構（non-canonical basic-Helix-Loop-Helix strueture，bHLH）的DNA結合蛋白，主要在分生組織中起作用，與細胞分化和生長有很大的關系[3]。研究發現，TB1對禾本科植物側枝的發生有著重要的作用，是多條信號調節通路的一個節點，它和上下游基因互作共同參與禾本科植物分蘗性狀的調控[4-6]。TB1最先從玉米中克隆出來[4]，后又在水稻[7]、高粱[8]、擬南芥[9]等作物中克隆出其同源基因并進行了功能分析，結果表明TB1在單子葉植物和雙子葉植物中結構和功能方面都表現出高度的保守性。目前有關甘蔗TB1基因的研究報道還十分少，Pribil等[10]曾在研究報道中簡單提到轉入水稻或甘蔗TB1基因的轉基因甘蔗株高有顯著增加，但分蘗數有所降低；Aitken[11]等檢測了一個從甘蔗TB1基因EST序列開發的EST-SSR標記，發現該標記對甘蔗分蘗性狀影響的效應較??；而關于甘蔗TB1基因的結構和功能尚不清楚，限制了利用甘蔗TB1基因調控甘蔗品種分蘗性狀表現、形態建成和產量改良。鑒于TB1基因在植物分枝（蘗）調控中的重要作用和在甘蔗研究中的滯后，本研究采用RT-PCR和RACE技術克隆到甘蔗TB1基因（ScTB1），并對其開展生物信息學分析，以期初步探索ScTB1的結構和功能，為后續該基因的功能解析、分子調控機制研究和分子輔助育種利用奠定良好的基礎。

1 材料與方法

1.1 材料

以國內甘蔗主栽品種新臺糖22號（ROC22）為研究材料，選分蘗盛期的分蘗苗莖尖組織部位，用于后續RNA的提取，所有研究材料由國家甘蔗種質資源圃提供。

1.2 方法

1.2.1 總RNA的提取及反轉錄 取甘蔗品種ROC22分蘗苗莖尖新鮮組織于液氮中迅速研磨成粉末，使用全式金TransZolTM Plant RNA提取試劑盒提取總RNA，并用1.0%的瓊脂糖凝膠和Bio Drop Lite PC超微量可見紫外分光光度計檢測RNA質量和濃度，然后使用反轉錄試劑盒Trans Script One-step gDNA Removal and cDNA Synthesis Super Mix反轉錄成cDNA，于-20 ℃保存備用。

1.2.2 基因克隆 （1）基因片段擴增。使用Primer Premier 5.0軟件，根據水稻、玉米、高粱等作物TB1基因保守區序列設計RT-PCR引物：ScTB1-113-F和ScTB1-815-R（見表1）。以上述反轉錄cDNA為模板，使用Trans Taq DNA Polymerase High Fidelity（HiFi）PCR擴增獲得相應片段，并用pEASY-T5 Zero Cloning Kit進行連接后轉化大腸桿菌DH5α感受態細胞，選取5個克隆交由深圳華大基因研究院測序，測序結果使用DNAMAN軟件進行比對分析。

（2）RACE擴增ScTB1 3′和5′端。以比對正確的克隆序列為模板，設計3′和5′ RACE基因特異引物（GSP）（見表1），使用SMARTERTM RACE cDNA Amplification Kit試劑盒獲得5′ RACE cDNA和3′ RACE cDNA模板，按照試劑盒說明書進行5′ RACE和3′ RACE產物PCR擴增，PCR產物回收、連接轉化、測序和序列比對同上。使用Vector NTI 11.5軟件將基因片段與比對正確的RACE產物序列拼接得到甘蔗ScTB1 cDNA全長序列，根據所拼接的cDNA序列的編碼區外側設計引物ScTB1-131-F和ScTB1-1609-R擴增真實的ScTB1編碼區序列，測序分析進一步驗證拼接序列的正確性。

1.2.3 生物信息學分析 首先利用在線工具ORF Finder查找所得序列的開放閱讀框（ORF）；然后使用ExPASy服務器中的ProtScale和ProtParam軟件、SOP-MA在線軟件、PSORT軟件、ProtFun軟件、TMHMM軟件和SIGNA-IP軟件分別預測ScTB1的理化性質、二級結構、亞細胞定位、功能、跨膜結構和信號肽；再通過DNAMAN和MEGA等生物學軟件對TB1同源蛋白進行系統進化分析。

2 結果與分析

2.1 總RNA提取及ScTB1基因片段擴增

總RNA提取后電泳結果顯示3條完整的條帶（如圖1-A）， OD260/280在1.98～2.00之間，OD260/230在2.00～2.20間，表明所提RNA質量完整，鹽離子去除干凈，符合實驗要求。RT-PCR產物電泳后在750 bp上下處有一條帶（圖1-B），回收克隆測序并比對分析證實其為ScTB1基因的部分片段，長706 bp。

2.2 RACE擴增ScTB1基因3'端和5'端

ScTB1 RACE結果如圖2所示，3′ RACE經兩輪PCR反應在500 bp附近出現2條條帶（圖2-B），大片段回收測序比對后為ScTB1基因的3′端序列，長度為479 bp。5′ RACE經兩輪PCR反應在700 bp左右出現一條帶（圖2-D），回收測序比對后為ScTB1基因的5′端序列，長度為670 bp。RT-PCR獲得片段和兩段RACE產物序列拼接后得到ScTB1基因的cDNA全長序列，該序列總長為1 668 bp，包含274 bp的5′ UTR，1 149 bp的CDS和245 bp的3′ UTR。利用cDNA編碼區全長引物RT-PCR擴增得到1條長度為1 478 bp的片段（圖2-E），回收測序比對后為正確序列，證實所拼接獲得的cDNA序列正確，同時將序列提交至GenBank數據庫，登錄號：KP137673.1。

2.3 ScTB1基因序列的生物信息學分析

2.3.1 ScTB1編碼氨基酸序列結構和理化性質分析

通過在線工具找到ScTB1中最大的ORF并翻譯成氨基酸序列，如圖3所示。該ORF長1 149 bp，共編碼382個氨基酸，其中包括SP、TCP、R 3個保守結構域（圖3中加下劃線部分即為相應結構域），其中R域可能形成一個親水性α螺旋，這與Lukens等[12]的研究報道一致。

ScTB1蛋白質的理化性質預測結果表明該蛋白不穩定，不穩定系數為50.51；平均疏水性為-0.621，推測其可能是親水性氨基酸[13]，其親/疏水信號如圖4，圖4中峰值分布在-0.5以下的比分布在0.5以上的多，再次表明該蛋白為親水性蛋白。

2.3.2 ScTB1蛋白亞細胞定位、跨膜結構域和信號肽預測 ScTB1蛋白亞細胞定位預測結果顯示該蛋白可能定位于細胞質、線粒體基質和溶酶體，其概率分別為0.650、0.100和0.100，分布在細胞質的可能性較大。跨膜結構和信號肽預測結果（圖5、6）表明ScTB1不存在跨膜結構和信號肽，可初步判斷ScTB1不屬于膜蛋白也不屬于分泌蛋白，因而可能定位于細胞質或細胞器基質中，這與上述亞細胞定位結果相吻合。

2.3.3 ScTB1的二級結構和功能預測 ScTB1的二級結構預測結果如圖7所示，該蛋白主要由 α 螺旋（h）、延伸鏈（e）和無規則卷曲（c）構成，其中無規則卷曲最多（53.93%），α螺旋次之（37.17%），延伸鏈最少（8.90%）。ScTB1功能預測結果表明，該基因行使復制和轉錄作用概率為0.373，調控作用概率為0.257，參與中心代謝的概率為0.153，參與翻譯的概率為0.077，推測其可能為轉錄因子。

2.3.4 ScTB1同源性分析 為了進一步明確ScTB1的功能，了解它與其他植物TB1同源基因的系統進化關系，筆者選擇了高粱（GenBank： XP_002466597.1、AAK37490.1、AAK37494.1和AAK37483.1）、玉米（AAL17103.1）、Chionachne koenigii（GenBank：AAK37504.1）、Tripsacum cundinamarcae（GenBank： AAK37480.1）、薏苡（Coix lacryma-jobi）（AAK37484.1）等幾種植物TB1蛋白序列與ScTB1蛋白序列進行序列比對分析，結果見圖8，發現它們在保守結構域上氨基酸相似度很高，非保守結構域上存在一定變化。

TCP家族轉錄因子根據功能、進化差異等可分為兩個亞家族（CYC/TB1和PCFs），CYC/TB1亞族與植物頂端優勢的維持和分蘗發生密切相關，PCFs亞族主要參與植物細胞周期調控[14]。本研究使用MEGA 6.06軟件，選擇鄰接法構建了竹類、水稻、小麥、擬南芥等植物TCP家族蛋白的系統進化樹（見圖9）。根據系統進化樹，可將25個TCP家族蛋白聚為明顯的兩個大類群（clusterⅠ和clusterⅡ），大類群下又可聚為不同的亞類群，而且多數類群bootstrap值超過90%。ClusterⅠ主要由CYC/TB1亞家族蛋白組成，內部又可分為兩個亞類，第1亞類主要由竹類、高粱、玉米、甘蔗、小麥和水稻等禾本科植物TB1家族蛋白組成，第2亞類主要由擬南芥、菊類的部分TCP家族基因組成。clusterⅡ主要由以水稻PCF1、PCF2和擬南芥部分TCP轉錄因子為代表的PCFs亞家族蛋白組成。ScTB1蛋白被歸到clusterⅠ類群，與高粱TB1和玉米TB1聚為一個亞類，并得到99%的bootstrap，因此屬于CYC/TB1亞族，由此推測ScTB1可能參與了植物頂端優勢的維持和分蘗性狀的調控。

3 討論與結論

分蘗性狀是甘蔗品種重要的農藝性狀，對于最終產量的形成具有重要作用，因此挖掘調控甘蔗分蘗性狀表現的基因，對于利用基因工程和分子育種手段精準調控甘蔗有效分蘗莖數和理想株型的塑造都具有重要意義，本研究旨在初步預測ScTB1的結構和功能，以期最終闡明其功能后通過基因沉默等分子技術將其運用于甘蔗生產。TB1屬于TCP家族蛋白，是調控植物分枝（蘗）的關鍵調控基因[4-9]。TCP家族蛋白屬植物特有轉錄因子，含有不典型bHLH DNA結合蛋白，根據其結構和功能可將TCP家族分為兩個亞家族（CYC/TB1和PCF），CYC/TB1亞家族主要參與植物器官形態建成，PCF亞家族參與細胞周期調控[14]。TCP家族具SP、TCP和R結構域，其中TCP和SP結構域為兩個亞家族共有，而R結構域僅存在于CYC/TB1亞家族的某些基因當中[3]。本研究克隆到的ScTB1編碼蛋白含有R結構域，顯然屬于CYC/TB1亞族，在系統進化上，其與玉米、水稻、高粱等TB1蛋白聚為一類，也初步證實筆者克隆的ScTB1屬于CYC/TB1亞族。

文中對ScTB1進行亞細胞定位預測結果顯示其主要可能定位于細胞質，這與其屬于轉錄因子這一說法似乎矛盾，但軟件預測結果只能作為初步參考，而且李超[15]碩士論文中也明確指出TCP家族蛋白在細胞質中的分布較細胞核中多，說明預測結果也具備一定參考性。此外，李超[15]論文中也指出TCP結構域的堿性區域包含一個核定位信號，所以，TB1蛋白也可能定位于細胞核中，與其屬于轉錄因子一說吻合。生物信息學分析結果只是一個初步的預測，具體亞細胞定位還有待實驗進一步證實。

前人研究表明TB1基因與分蘗的發生密切相關[16]，Doebley等[17]和Hubbard等[18]研究表明TB1基因在玉米葉腋分生組織表達，抑制側芽的生長。水稻OsTB1基因過表達植株表現出分蘗芽萌發受抑制和分蘗莖明顯減少現象[7]，玉米、小麥TB1基因在水稻中過表達水稻也表現出分蘗減少的現象[6]。Kebrom等[8]針對高粱TB1基因的功能分析也表明SbTB1對高粱分蘗芽的萌發具有明顯的抑制作用。本研究所克隆的ScTB1基因與高粱、玉米和水稻的TB1在基因序列和蛋白結構上表現出較高的保守性，系統進化樹上聚為較為緊密的一類，加之TB1蛋白理化性質、結構和功能預測等可初步推測ScTB1可能參與了甘蔗分蘗性狀的調控，該作用機制將在后續的功能分析中得到驗證。

參考文獻

[1] Henry R J， Kole C. Genetics， genomics and breeding of sugarcane[M]. Science Publishers， Inc， 2010.

[2] 陳如凱. 現代甘蔗遺傳育種[M]. 北京： 中國農業出版社， 2011： 1-19.

[3] Cubas P， Lauter N， Doebley J， et al. The TCP domain： a motif found in proteins regulating plant growth and development[J]. The Plant Journal， 1999， 18（2）： 215-222.

[4] Doebley J， Stec A， Hubbard L. The evolution of apical dominance in maize[J]. Nature， 1997， 386： 485-488.

[5] Finlayson S A. Arabidopsis TEOSINTE BRANCHED1-LIKE 1 regulates axillary bud outgrowth and is homologous to monocot TEOSINTE BRANCHED1[J]. Plant and Cell Physiology， 2007， 48（5）： 667-677.

[6] Choi M S， Woo M O， Koh E B， et al. Teosinte Branched 1 modulates tillering in rice plants[J]. Plant Cell Reports， 2012， 31（1）： 57-65.

[7] Takeda T， Suwa Y， Suzuki M， et al. The OsTB1 gene negatively regulates lateral branching in rice[J]. The Plant Journal， 2003， 33（3）： 513-520.

[8] Kebrom T H， Burson B L， Finlayson S A. Phytochrome B represses Teosinte Branched1 expression and induces sorghum axillary bud outgrowth in response to light signals[J]. Plant Physiology， 2006， 140（3）： 1 109-1 117.

[9] Aguilar-Martínez J A， Poza-Carrión C， Cubas P. Arabidopsis BRANCHED1 acts as an integrator of branching signals within axillary buds[J]. The Plant Cell， 2007， 19（2）： 458-472.

[10] Pribil M， Hermann S， Dun G， et al. Altering sugarcane shoot architecture through genetic engineering： prospects for increasing cane and sugar yield[J]. Proceedings of the Conference of the Australian Society of Sugar Cane Technologists， 2007， 29： 251-257.

[11] Aitken K S， Hermann S， Karno K， et al. Genetic control of yield related stalk traits in sugarcane[J]. Theoretical & Applied Genetics， 2008， 117（7）： 1 191-1 203.

[12] Lukens L， Doebley J. Molecular evolution of the teosinte branched gene among maize and related grasses[J]. Molecular Biology and Evolution， 2001， 18（4）： 627-638.

[13] 付海輝， 辛培堯， 許玉蘭， 等. 幾種經濟植物UFGT基因的生物信息學分析[J]. 基因組學與應用生物學， 2010， 30（1）： 92-102.

[14] 王玉紅， 徐立華， 李 軍. 玉米TCP家族基因生物信息學鑒定與分析[J]. 山東農業科學， 2014， 46（4）： 1-6.

[15] 李 超. 農桿菌介導的豌豆花瓣瞬間表達系統的建立及利用該系統對百脈根TCP基因亞細胞定位的研究[D]. 上海： 上海交通大學， 2009.

[16] Hu W， Zhang S， Zhao Z， et al. The Analysis of the Structure and Expression of OsTBl Gene in rice[J]. Journal of Plant Physiology and Molecular Biology， 2002， 29（6）： 507-514.

[17] Doebley J， Stec A， Gustus C. teosinte branched1 and the origin of maize： evidence for epistasis and the evolution of dominance[J]. Genetics， 1995， 141（1）： 333.

[18] Hubbard L， McSteen P， Doebley J， et al. Expression patterns and mutant phenotype of teosinte branched1 correlate with growth suppression in maize and teosinte[J]. Genetics， 2002， 162（4）： 1 927-1 935.