普通煙草SBP轉(zhuǎn)錄因子家族的全基因組鑒定及其進(jìn)化和表達(dá)分析

劉 成，李曉旭，蘇玉龍，郭永峰

（1.中國農(nóng)業(yè)科學(xué)院煙草研究所，青島 266101；2.中國農(nóng)業(yè)科學(xué)院研究生院，北京 100081）

普通煙草SBP轉(zhuǎn)錄因子家族的全基因組鑒定及其進(jìn)化和表達(dá)分析

劉 成1，2，李曉旭1，2，蘇玉龍1，2，郭永峰1*

（1.中國農(nóng)業(yè)科學(xué)院煙草研究所，青島 266101；2.中國農(nóng)業(yè)科學(xué)院研究生院，北京 100081）

SBP轉(zhuǎn)錄因子家族作為植物特異性的轉(zhuǎn)錄因子具有重要的生物學(xué)功能，廣泛參與花和果實(shí)發(fā)育等諸多生命過程。利用普通煙草（Nicotiana tabacum）TN90基因組數(shù)據(jù)，通過隱馬爾科夫模型檢索，SMART和Pfam分析，鑒定了普通煙草SBP轉(zhuǎn)錄因子家族成員；通過多序列比對鑒定并分析了普通煙草SBP結(jié)構(gòu)域的結(jié)構(gòu)特征；通過鄰接法、極大似然法進(jìn)行了系統(tǒng)發(fā)育分析；分別利用GSDS和MEME對SBP轉(zhuǎn)錄因子家族成員的基因結(jié)構(gòu)和蛋白保守結(jié)構(gòu)域進(jìn)行了分析；利用普通煙草TN90轉(zhuǎn)錄組數(shù)據(jù)分析了SBP轉(zhuǎn)錄因子家族成員在不同組織和時(shí)期的表達(dá)情況，并對鑒定出的SBP轉(zhuǎn)錄因子家族成員進(jìn)行了GO注釋分析。結(jié)果表明，普通煙草基因組編碼32個(gè)SBP轉(zhuǎn)錄因子家族成員，氨基酸序列長度差異較大；系統(tǒng)發(fā)育分析表明，所有的SBP轉(zhuǎn)錄因子家族成員可分為8個(gè)亞家族。轉(zhuǎn)錄組數(shù)據(jù)顯示，SBP基因家族成員在不同組織和時(shí)期中表達(dá)迥異，但在花器官中均有較高的表達(dá)量。GO注釋結(jié)果表明，該家族成員作為轉(zhuǎn)錄因子，在植物的生長發(fā)育、逆境響應(yīng)等過程中發(fā)揮作用。本研究為煙草及其他植物中SBP轉(zhuǎn)錄因子家族基因的鑒定和功能研究提供了基礎(chǔ)。

生物信息學(xué)分析；SBP轉(zhuǎn)錄因子家族；普通煙草

轉(zhuǎn)錄因子是參與轉(zhuǎn)錄調(diào)控的一類蛋白，通過與特定序列的DNA結(jié)合，能夠激活或抑制目標(biāo)基因的轉(zhuǎn)錄，從而調(diào)控基因的表達(dá)［1］。SQUAMOSA Promoter Binding Proteins（SBPs）轉(zhuǎn)錄因子家族是一類植物特有的重要轉(zhuǎn)錄因子，由多個(gè)成員組成，最早發(fā)現(xiàn)于金魚草中，廣泛參與植物生長、發(fā)育以及多種生理生化過程［2］。SPL（SBP-like）基因主要調(diào)控花和果實(shí)發(fā)育，有報(bào)道稱擬南芥中SPL3，SPL8，SPL14在開花、孢子形成、赤霉素信號傳導(dǎo)、抵抗毒素侵襲等生命過程中發(fā)揮著重要作用［3-6］。作為植物特異的轉(zhuǎn)錄因子，所有已知的SBP轉(zhuǎn)錄因子家族成員蛋白中均含有一段由79個(gè)氨基酸殘基組成的高度保守的DNA結(jié)合域。擬南芥SBP轉(zhuǎn)錄因子的DNA結(jié)合域由兩個(gè)分離的鋅指結(jié)構(gòu)組成，一個(gè)是C3H或C4，另一個(gè)是C2HC［7］。根據(jù)文獻(xiàn)報(bào)道，SBP轉(zhuǎn)錄因子家族很多成員都受到miR156的調(diào)控。目前在擬南芥和水稻中都各發(fā)現(xiàn)有11個(gè)SPL基因含有miR156的識別位點(diǎn)［8］，擬南芥營養(yǎng)發(fā)育過程中SPL3的時(shí)空性表達(dá)受到miR156的調(diào)控［9］，水稻中的miR156和某些OsSPL基因也會發(fā)生組織特異性互作［10］。隨著越來越多的植物基因組測序的完成，SBP轉(zhuǎn)錄因子家族在全基因組范圍內(nèi)已經(jīng)被廣泛地鑒定。

煙草是重要的經(jīng)濟(jì)作物，同時(shí)也是一種重要的模式植物，經(jīng)常作為分子生物學(xué)的研究對象。在普通煙草基因組注釋信息公布前，人們往往采取同源克隆結(jié)合RACE技術(shù)等手段研究煙草中重要的功能基因，如陸瑩等［11］利用馬鈴薯的同源基因在煙草中克隆到了光周期的關(guān)鍵基因NtCO1并進(jìn)行了初步分析；劉好寶等［12］在煙草中同源克隆到4個(gè)鉀離子通道基因：NKT3、NKT4、NKT5和NKT6。隨著生物信息學(xué)的發(fā)展和煙草基因組測序、注釋的完成，煙草中如鈣依賴蛋白激酶基因家族（CDPK Gene Family）［13］等基因家族都已成為研究的熱點(diǎn)。但煙草中針對SBP轉(zhuǎn)錄因子家族的研究還未見報(bào)道。本研究在普通煙草的全基因組注釋信息的基礎(chǔ)上，利用生物信息學(xué)方法鑒定出普通煙草32個(gè)SBP轉(zhuǎn)錄因子家族成員，分析了各成員的進(jìn)化關(guān)系、基因結(jié)構(gòu)、蛋白保守結(jié)構(gòu)域、表達(dá)情況，并進(jìn)行GO注釋分析，為進(jìn)一步研究普通煙草SBP轉(zhuǎn)錄因子家族成員的功能奠定基礎(chǔ)。

1 材料與方法

1.1 普通煙草SBP轉(zhuǎn)錄因子家族成員的鑒定與序列分析

根據(jù)文獻(xiàn)報(bào)道，在TAIR數(shù)據(jù)庫（http://www. arabidopsis.org/）下載了擬南芥SBP轉(zhuǎn)錄因子家族成員的蛋白全長序列，在PlantTFDB（http:/ /planttfdb.cbi.pku.edu.cn/）數(shù)據(jù)庫［14］下載了番茄、馬鈴薯的SBP轉(zhuǎn)錄因子家族成員的蛋白全長序列。對于擬南芥、番茄、馬鈴薯的SBP蛋白全長序列，利用MAFFT分別進(jìn)行序列比對，使用HMMER 3.0分別建立隱馬爾科夫模型序列譜（profile hidden Markov models，profile HMMs）。在茄科植物數(shù)據(jù)庫（http://solgenomics.net/）下載普通煙草TN90蛋白序列數(shù)據(jù)庫，分別使用建立的3個(gè)隱馬爾科夫模型序列譜對普通煙草TN90蛋白序列數(shù)據(jù)庫進(jìn)行檢索，選取E值為0.01，手工去除冗余序列，得到候選序列。最后針對所獲得的候選序列，利用SMART和Pfam分析蛋白結(jié)構(gòu)，剔除不含完整SBP結(jié)構(gòu)域的蛋白序列。利用ExPASy Proteomics Server（http://expasy.org/）對所有SPB家族成員蛋白氨基酸序列進(jìn)行分子量、等電點(diǎn)預(yù)測分析。

1.2 多序列比對和系統(tǒng)發(fā)育分析

將擬南芥和普通煙草TN90的SBP蛋白全長序列用MAFFT［15］在默認(rèn)參數(shù)下進(jìn)行多序列比對。選取保守的SBP結(jié)構(gòu)域，利用Texshade將結(jié)果可視化。利用擬南芥、普通煙草SBP結(jié)構(gòu)域蛋白序列的比對結(jié)果，使用MGEA5構(gòu)建鄰接樹（Neighbor-Joining， NJ tree），參數(shù)設(shè)置為：Bootstrap檢驗(yàn)次數(shù)設(shè)為1000次，并選擇Poisson Model 和Pairwise Deletion進(jìn)行分析。利用同樣的比對結(jié)果，使用PhyML構(gòu)建極大似然樹［16］（Maxumum Likeihood， ML tree）。Bootstrap設(shè)為100次，并利用ProtTest 2.4［17］對擬南芥和普通煙草SBP結(jié)構(gòu)域氨基酸序列比對結(jié)果進(jìn)行分析，選擇LG+G模型構(gòu)建極大似然樹，并利用Tree View 1.4.0顯示系統(tǒng)發(fā)生樹。

1.3 基因結(jié)構(gòu)和結(jié)構(gòu)域分析

在茄科植物數(shù)據(jù)庫（http://solgenomics.net/），下載普通煙草TN90基因組信息，通過Perl程序分析并提取SBP基因家族成員的基因結(jié)構(gòu)信息，使用GSDS（http://gsds.cbi. pku.edu.cn/）工具對分析結(jié)果進(jìn)行基因結(jié)構(gòu)的可視化展示。利用Multiple Expectation Maximization for Motif Elicitation（MEME， http://meme.nbcr.net/meme/cgi-bin/mem e.cgi）［18］對SBP轉(zhuǎn)錄因子家族成員的蛋白序列基序（motif）進(jìn)行分析。

1.4 普通煙草SBP轉(zhuǎn)錄因子家族基因表達(dá)分析及GO注釋分析

在NCBI SRA（http://www.ncbi.nlm.nih.gov/sra）數(shù)據(jù)庫檢索并下載TN90轉(zhuǎn)錄組測序（RNA-seq）數(shù)據(jù)，提取表達(dá)數(shù)據(jù)，利用R語言對轉(zhuǎn)錄組表達(dá)數(shù)據(jù)進(jìn)行分析并提取SBP基因家族成員的基因表達(dá)數(shù)據(jù)，將表達(dá)數(shù)據(jù)標(biāo)準(zhǔn)化后，使用R語言pheatmap程序繪制熱圖（heatmap），將表達(dá)數(shù)據(jù)可視化。使用普通煙草SBP轉(zhuǎn)錄因子家族蛋白序列，利用Blast2GO進(jìn)行GO（Gene Ontology， GO）注釋［19］，WEGO工具［20］分析注釋結(jié)果并進(jìn)行可視化展示。

2 結(jié) 果

2.1 普通煙草SBP轉(zhuǎn)錄因子家族成員鑒定及理化特征分析

利用文獻(xiàn)報(bào)道的擬南芥SBP轉(zhuǎn)錄因子家族成員蛋白序列以及PlantTFDB中番茄和馬鈴薯的SBP轉(zhuǎn)錄因子家族成員蛋白序列，分別構(gòu)建了隱馬爾可夫模型序列譜。利用構(gòu)建的3個(gè)隱馬爾科夫模型序列譜分別檢索普通煙草TN90蛋白數(shù)據(jù)庫，綜合3個(gè)檢索結(jié)果，人工去除冗余序列，獲得候選蛋白序列。通過SMART和Pfam對候選蛋白序列進(jìn)行結(jié)構(gòu)域分析，剔除不含完整SBP功能域的蛋白序列，最終預(yù)測得到32個(gè)煙草SBP基因家族成員，命名為NtSPL1-NtSPL32（表1）。

表1 普通煙草SBP轉(zhuǎn)錄因子家族成員特征Table 1 Characteristics of SBP transcription factor family in N. tabacum

對普通煙草SBP轉(zhuǎn)錄因子家族成員進(jìn)行蛋白質(zhì)序列及理化性質(zhì)分析發(fā)現(xiàn)，這些氨基酸的長度變化差異較大，在133～999 aa，蛋白分子量則介于15.27～110.97 kD。家族成員蛋白的理論等電點(diǎn)變化較小，在5.15～9.76，其中第I-V亞族（7.64～9.43）理論等電點(diǎn)都在堿性范圍內(nèi)，表明這些亞家族蛋白富含堿性氨基酸，而第VIII亞族（NtSPL6和NtSPL32）理論等電點(diǎn)較低（5.15～5.29），其蛋白含有較多的酸性氨基酸。所有SBP轉(zhuǎn)錄因子家族蛋白疏水性均為負(fù)值，表明此家族均屬于親水性蛋白，其中第Ⅶ亞族蛋白疏水性普遍較低，從-1.313到-1.084，說明第VII亞族蛋白均具有較強(qiáng)的親水性。

2.2 普通煙草SBP結(jié)構(gòu)域的鑒定和分析

圖1 擬南芥、普通煙草SBP轉(zhuǎn)錄因子家族蛋白SBP結(jié)構(gòu)域分析Fig. 1 Conserved SBP domain analysis of SBP family proteins in N. tabacum and A. thaliana

將文獻(xiàn)報(bào)道的16條擬南芥SBP轉(zhuǎn)錄因子家族成員蛋白序列和預(yù)測的32條普通煙草SBP轉(zhuǎn)錄因子家族成員蛋白序列，使用MAFFT在默認(rèn)參數(shù)下進(jìn)行多序列比對，發(fā)現(xiàn)所有序列均含保守的SBP結(jié)構(gòu)域。選取保守的SBP結(jié)構(gòu)域，利用Texshade將結(jié)果可視化（圖1），擬南芥和普通煙草的SBP結(jié)構(gòu)域高度相似，均含79個(gè)氨基酸殘基。利用核磁共振方法，擬南芥中兩個(gè)SBP結(jié)構(gòu)域的三維溶液構(gòu)象已經(jīng)被報(bào)道［21］，兩者都含2個(gè)鋅指結(jié)構(gòu)（Zinc finger）。AtSPL3的N端鋅指結(jié)構(gòu)（Zn-1）為Cys3His（C3H），C端鋅指結(jié)構(gòu)（Zn-2）為Cys2HisCys（C2HC）。AtSPL7的N端鋅指結(jié)構(gòu)為Cys4（C4），C端鋅指結(jié)構(gòu)與同樣為C2HC。經(jīng)過對擬南芥和普通煙草SBP結(jié)構(gòu)域的分析，發(fā)現(xiàn)所有的AtSPL和NtSPL轉(zhuǎn)錄因子均包含兩個(gè)鋅指結(jié)構(gòu)：AtSPL7、NtSPL6和NtSPL32的N段鋅指結(jié)構(gòu)為C4，其余成員的N端鋅指結(jié)構(gòu)均為C3H；所有的C端鋅指結(jié)構(gòu)均為C2HC。與擬南芥類似，普通煙草SBP轉(zhuǎn)錄因子家族成員蛋白SBP結(jié)構(gòu)域C端存在由大量的堿性氨基酸殘基組成的高度保守的核定位信號（nuclear localization signal， NLS）［21］。

2.3 擬南芥和普通煙草SBP轉(zhuǎn)錄因子家族的系統(tǒng)發(fā)育分析

為了揭示普通煙草SBP轉(zhuǎn)錄因子家族的進(jìn)化關(guān)系，參照Guo等［2］的方法，使用擬南芥和普通煙草SBP轉(zhuǎn)錄因子家族成員SPB結(jié)構(gòu)域氨基酸序列，利用MEGA和PhyML分別構(gòu)建了鄰接樹（圖2）和極大似然樹（圖3），兩者具有高度相似的拓?fù)浣Y(jié)構(gòu)。構(gòu)建的鄰接樹顯示，擬南芥和普通煙草SBP轉(zhuǎn)錄因子家族成員被分為8個(gè)亞家族每一個(gè)亞家族都包含了兩個(gè)物種中至少一個(gè)的SBP轉(zhuǎn)錄因子家族成員，這說明了SBP轉(zhuǎn)錄因子家族成員的分化時(shí)間可能早于兩個(gè)物種的分化時(shí)間。其中亞家族I-VIII含有擬南芥SBP轉(zhuǎn)錄因子家族成員的個(gè)數(shù)分別為：1、2、2、3、1、3、3和1；含有普通煙草SBP轉(zhuǎn)錄因子家族成員的個(gè)數(shù)分別為：7、2、2、4、3、4、8和2。普通煙草SBP轉(zhuǎn)錄因子家族成員的數(shù)量約為擬南芥的兩倍，但是擬南芥和普通煙草在不同亞家族中SBP轉(zhuǎn)錄因子家族成員的倍數(shù)關(guān)系差異較大，這說明在擬南芥和普通煙草分化后，SBP各個(gè)亞家族可能發(fā)生了不同的進(jìn)化事件。此外，在亞家族I中普通煙草有7個(gè)SBP轉(zhuǎn)錄因子家族成員，而擬南芥只有一個(gè)SBP轉(zhuǎn)錄因子家族成員位于該分支中，說明在普通煙草物種形成后，該亞家族內(nèi)SBP轉(zhuǎn)錄因子家族成員可能經(jīng)歷了多次復(fù)制事件。在8個(gè)亞家族中，亞家族VIII與其他亞家族的親緣關(guān)系較遠(yuǎn)，形成一個(gè)較為獨(dú)立的分支，這與以前的報(bào)道是一致的［22］。進(jìn)一步分析發(fā)現(xiàn)，亞家族8中3個(gè)成員AtSPL7、NtSPL6和NtSPL32在SBP結(jié)構(gòu)域中的N段鋅指結(jié)構(gòu)為C4，而其他亞家族所有成員的N端鋅指結(jié)構(gòu)均為C3H。

圖2 普通煙草、擬南芥SBP轉(zhuǎn)錄因子家族成員鄰接系統(tǒng)發(fā)育樹Fig. 2 Neighbor-joining phylogenetic tree of SBP family in N. tabacum and A. thaliana

2.4 擬南芥和普通煙草SBP基因家族的基因結(jié)構(gòu)分析

圖3 擬南芥、番茄SBP轉(zhuǎn)錄因子家族極大似然系統(tǒng)發(fā)育樹及基因結(jié)構(gòu)Fig. 3 The Maximum-Likelihood phylogenetic tree and gene structure of SBP family members in N. tabacum and A. thaliana

為了分析普通煙草SBP基因的結(jié)構(gòu)特征，利用GSDS構(gòu)建了擬南芥和普通煙草SBP基因家族成員的基因結(jié)構(gòu)圖（圖3）。內(nèi)含子-外顯子結(jié)構(gòu)圖表明，不同亞家族的SPL基因之間內(nèi)含子數(shù)目變異比較大，最多的有9個(gè)內(nèi)含子，最少的有1個(gè)內(nèi)含子。在亞家族I、II、III和V中，大部分成員含有2個(gè)內(nèi)含子；亞家族IV中，大部分成員含有3個(gè)內(nèi)含子；亞家族VII中，大部分成員含有1個(gè)內(nèi)含子；亞家族VI和VIII中，大部分成員含有9個(gè)內(nèi)含子。同一亞家族內(nèi)，擬南芥與普通煙草直系同源的SPL基因，其外顯子、內(nèi)含子的組織結(jié)構(gòu)高度相似并區(qū)別于其他亞家族內(nèi)的SPL基因，這同樣說明了SBP基因家族成員的分化時(shí)間可能要早于擬南芥和普通煙草的物種分化時(shí)間。進(jìn)一步分析發(fā)現(xiàn)，同一亞家族內(nèi)，擬南芥和普通煙草的旁系同源SPL基因的內(nèi)含子的組織結(jié)構(gòu)也相對一致，這也從側(cè)面反映了系統(tǒng)發(fā)育分析的可靠性。

2.5 擬南芥和普通煙草SBP轉(zhuǎn)錄因子家族的蛋白保守結(jié)構(gòu)域分析

利用序列保守結(jié)構(gòu)域識別工具M(jìn)EME，對擬南芥SBP轉(zhuǎn)錄因子和預(yù)測的普通煙草SBP轉(zhuǎn)錄因子的蛋白序列的保守結(jié)構(gòu)域進(jìn)行分析（圖4）。分析結(jié)果表明，所有的擬南芥和普通煙草SBP轉(zhuǎn)錄因子家族成員均含有79個(gè)氨基酸殘基組成的SBP結(jié)構(gòu)域（Motif 1）。除了SBP結(jié)構(gòu)域，如Motif 2等一些新的保守結(jié)構(gòu)域被鑒定出來，同一亞家族內(nèi)的SBP轉(zhuǎn)錄因子家族成員的保守結(jié)構(gòu)域的種類、數(shù)目和組織形式具有很強(qiáng)的一致性，同樣從側(cè)面反映了系統(tǒng)發(fā)育分析的可靠性。在亞家族VII內(nèi)，擬南芥和普通煙草的SBP轉(zhuǎn)錄因子家族成員除了SBP結(jié)構(gòu)域外沒有其他的保守結(jié)構(gòu)域，而亞家族VI內(nèi)存在Motif2、Motif3、Motif5和Motif7等其他亞家族中沒有或不典型的保守結(jié)構(gòu)域，這說明亞家族VI的成員可能具有其他亞家族成員所不具有的功能。此外，亞家族I、II、III和IV中擬南芥有8個(gè)成員、普通煙草有13個(gè)成員有一個(gè)位于Motif 8中的保守基序ALSLLS，據(jù)文獻(xiàn)報(bào)道，其對應(yīng)的mRNA序列可能是mRNA156的靶序列［8］。

2.6 普通煙草SBP基因家族的基因表達(dá)分析

使用TN90轉(zhuǎn)錄組測序（RNA-seq）數(shù)據(jù)進(jìn)行分析并提取SBP基因家族成員的基因表達(dá)數(shù)據(jù)，使用R語言將表達(dá)數(shù)據(jù)可視化。從圖5中可以看出，SBP基因家族成員在根、莖、幼葉、成熟葉、衰老葉、幼花等6個(gè)組織中均有表達(dá)。進(jìn)一步的分析發(fā)現(xiàn)，NtSPL基因的表達(dá)具有一定的組織特異性。所有NtSPL基因在幼花中都有表達(dá)，而且絕大多數(shù)成員還具有較高的表達(dá)量，這與擬南芥AtSPL基因在花中高度表達(dá)的報(bào)道相一致，表明SPL基因在植物的花發(fā)育中起重要作用［23］。進(jìn)一步的分析發(fā)現(xiàn)在根、成熟葉、老葉等組織中，有些NtSPL基因表達(dá)量很低，如根中的第VII亞家族基因，成熟葉和老葉中的第I、III、V亞家族基因表達(dá)量普遍較低。同一亞家族的NtSPL基因在同一組織中的表達(dá)具有差異，成員較多的亞家族如第Ⅰ、第Ⅶ亞家族，在同一組織如莖、衰老葉中表達(dá)差異較大。同一亞家族的NtSPL基因在不同組織中的表達(dá)差異較為明顯，只有第Ⅵ亞家族的NtSPL基因在所檢測的6個(gè)組織中都有較高的表達(dá)量。此外，有些NtSPL基因隨著組織的發(fā)育而呈規(guī)律性的表達(dá)，如NtSPL2隨著葉片的衰老表達(dá)量逐漸升高，NtSPL10在衰老葉中表達(dá)量顯著增加，初步推測它們可能與植物的衰老進(jìn)程相關(guān)。

2.7 普通煙草SBP轉(zhuǎn)錄因子家族GO注釋分析

通過對普通煙草SBP蛋白序列進(jìn)行GO注釋發(fā)現(xiàn)，經(jīng)鑒定普通煙草SBP轉(zhuǎn)錄因子家族成員的細(xì)胞組成（Cellular Component）、分子功能（Molecular Function）和所參與的生物學(xué)過程（Biological Process）與擬南芥SBP轉(zhuǎn)錄因子家族成員相似。GO注釋分析結(jié)果（圖6）表明，擬南芥和普通煙草的SBP蛋白均在細(xì)胞核中，發(fā)揮DNA結(jié)合和轉(zhuǎn)錄調(diào)節(jié)相關(guān)的功能，這與SBP轉(zhuǎn)錄因子家族成員作為轉(zhuǎn)錄因子的功能相一致。在細(xì)胞生物學(xué)過程方面，SBP轉(zhuǎn)錄因子家族成員的注釋呈多樣化。所有的擬南芥和普通煙草SBP蛋白

被注釋為參與到花器官發(fā)育，這與SBP成員在花器官中有較高的表達(dá)量是相吻合的。大部分蛋白還參與到葉片、莖和花粉管等植物器官的生長發(fā)育過程中，一些蛋白參與植物的防御響應(yīng)等過程中，這與已有的文獻(xiàn)報(bào)道是相符的［2］。

圖4 擬南芥、番茄SBP轉(zhuǎn)錄因子家族蛋白保守結(jié)構(gòu)域分析Fig. 4 Conserved domains within different subgroups of SBP family members in N. tabacum and A. thaliana

圖5 普通煙草SBP基因家族成員表達(dá)模式分析Fig. 5 Expression patterns of the SBP gene family members in N. tabacum

圖6 擬南芥、普通煙草SBP蛋白GO注釋Fig. 6 Gene ontology analysis of SBP proteins in N. tabacum and A. thaliana

3 討 論

SBP基因家族所編碼的具有植物特異性的轉(zhuǎn)錄因子具有許多重要的生物學(xué)功能，廣泛參與植物生長發(fā)育調(diào)節(jié)、逆境脅迫響應(yīng)及孢子、花和果實(shí)發(fā)育等諸多生命過程。本研究采用生物信息學(xué)手段，首次從普通煙草基因組中，鑒定出32個(gè)SBP轉(zhuǎn)錄因子家族成員，所有成員均含有保守的SBP結(jié)構(gòu)域。通過對SBP結(jié)構(gòu)域進(jìn)行分析發(fā)現(xiàn)，所有的NtSPL均具有兩個(gè)典型的鋅指結(jié)構(gòu)和一個(gè)位于C端的保守核定位信號。針對擬南芥和普通煙草SBP轉(zhuǎn)錄因子家族成員的蛋白序列，通過鄰接法和極大似然法構(gòu)建的系統(tǒng)發(fā)育樹具有高度相似的拓?fù)浣Y(jié)構(gòu)。系統(tǒng)發(fā)育分析表明擬南芥和普通煙草的SBP轉(zhuǎn)錄因子家族成員被分成了8個(gè)亞家族，每一個(gè)亞家族都包含了至少一個(gè)擬南芥和普通煙草SBP轉(zhuǎn)錄因子家族成員，說明SBP轉(zhuǎn)錄因子家族成員的分化形成可能早于擬南芥和普通煙草的分化時(shí)間。通過對擬南芥和普通煙草SBP基因的結(jié)構(gòu)特征分析發(fā)現(xiàn)，不同亞家族中擬南芥和普通煙草的SPL基因之間內(nèi)含子數(shù)目變異較大。然而，通過比較擬南芥與普通煙草的直系同源SPL基因，發(fā)現(xiàn)其外顯子、內(nèi)含子的組織結(jié)構(gòu)高度相似，這同樣印證了SBP基因家族成員的分化形成時(shí)間可能要早于擬南芥和普通煙草的物種分化時(shí)間。通過對SBP轉(zhuǎn)錄因子家族成員結(jié)構(gòu)域的分析發(fā)現(xiàn)，相較其他亞家族而言，亞家族VI內(nèi)存在許多較為復(fù)雜的保守結(jié)構(gòu)域。此外，普通煙草中有13個(gè)成員存在一個(gè)位于C端的保守基序ALSLLS，其對應(yīng)的mRNA序列可能是mRNA156的靶序列。SBP基因家族表達(dá)分析發(fā)現(xiàn)，NtSPLs基因的表達(dá)具有一定的組織特異性，所有NtSPL基因在幼花中都有表達(dá)，而且絕大多數(shù)成員還具有較高的表達(dá)量。GO注釋表明，普通煙草SBP轉(zhuǎn)錄因子家族成員注釋為參與到花器官發(fā)育，這與煙草SBP基因家族成員在花器官中有較高的表達(dá)量是相吻合的。大部分SBP轉(zhuǎn)錄因子家族成員被注釋為參與到葉片、莖和花粉管等植物器官的生長發(fā)育過程中，一些成員被注釋為參與植物的防御響應(yīng)等過程中，但具體功能有待于利用分子生物學(xué)等手段進(jìn)行進(jìn)一步的基因功能驗(yàn)證。

4 結(jié) 論

本研究利用普通煙草TN90基因組數(shù)據(jù)，通過生物信息學(xué)方法鑒定了普通煙草基因組編碼的32個(gè)SBP轉(zhuǎn)錄因子家族成員，并將這些成員分為8個(gè)亞家族。通過分析轉(zhuǎn)錄組數(shù)據(jù)，發(fā)現(xiàn)SBP基因家族成員在花器官中的表達(dá)水平較高。GO注釋結(jié)果表明，SBP家族在植物的器官發(fā)育、逆境防御等生物過程中發(fā)揮重要作用。

本研究通過生物信息學(xué)手段對SBP轉(zhuǎn)錄因子家族進(jìn)行了鑒定和一些功能預(yù)測，但其具體的生物學(xué)功能以及調(diào)控機(jī)制，還需要通過基因克隆、表達(dá)分析及轉(zhuǎn)基因等實(shí)驗(yàn)進(jìn)行分析驗(yàn)證。

［1］ Zhang S D， Ling L Z. Genome-wide identification and evolutionary analysis of the SBP-box gene family in castor bean［J］. PLoS One， 2014， 9（1）: e86688.

［2］ Guo A Y， Zhu Q H， Gu X， et al. Genome-wide identification and evolutionary analysis of the plant specific SBP-box transcription factor family［J］. Gene，2008， 418（1-2）: 1-8.

［3］ Cardon G H， H?hmann S， Nettesheim K， et al. Functional analysis of the Arabidopsis thaliana SBP-box gene SPL3: A novel gene involved in the floral transition［J］. The Plant Journal， 1997， 12（2）: 367-377.

［4］ Zhang Y， Schwarz S， Saedler H， et al. SPL8， a local regulator in a subset of gibberellin-mediated developmental processes in Arabidopsis［J］. Plant Molecular Biology， 2007， 63（3）: 429-439.

［5］ Unte U S. SPL8， an SBP-Box Gene That Affects Pollen Sac Development in Arabidopsis［J］. The Plant Cell Online， 2003， 15（4）: 1009-1019.

［6］ Stone J M， Liang X， Nekl E R， et al. Arabidopsis AtSPL14， a plant-specific SBP-domain transcription factor， participates in plant development and sensitivity to fumonisin B1［J］. The Plant Journal， 2005， 41（5）: 744-754.

［7］ Yamasaki K， Kigawa T， Inoue M， et al. A novel zincbinding motif revealed by solution structures of DNA-binding domains of Arabidopsis SBP-family transcription factors［J］. Journal of Molecular Biology，2004， 337（1）: 49-63.

［8］ Rhoades M W， Reinhart B J， Lim L P， et al. Prediction of plant microRNA targets［J］. Cell， 2002， 110（4）: 513-520.

［9］ Schwab R， Palatnik J F， Riester M， et al. Specific effects of microRNAs on the plant transcriptome［J］. Development Cell， 2005， 8（4）: 517-527.

［10］ Xie K， Wu C， Xiong L. Genomic organization，differential expression， and interaction of SQUAMOSA promoter-binding-like transcription factors and microRNA156 in rice［J］. Plant Physiology， 2006，142（1）: 280-293.

［11］ 陸瑩，劉艷華，任民，等. 煙草CONSTANS同源基因的克隆與分析［J］. 中國煙草科學(xué)，2013，34（3）：60-64

［12］ 王倩，劉好寶. 煙草重要基因篇：2. 煙草鉀吸收與轉(zhuǎn)運(yùn)相關(guān)基因［J］. 中國煙草科學(xué)，2014，35（2）:139-142

［13］ 劉貫山. 煙草重要基因篇：4. 煙草鈣依賴蛋白激酶基因［J］. 中國煙草科學(xué)，2014，35（4）：109-111.

［14］ Jin J， Zhang H， Kong L， et al. PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors［J］. Nucleic acids research， 2014， 42: D1182-D1187.

［15］ Katoh K， Kuma K， Toh H， et al. MAFFT version 5: improvement in accuracy of multiple sequence alignment［J］. Nucleic Acids Research， 2005， 33（2）: 511-518.

［16］ Guindon S， Lethiec F， Duroux P， et al. PHYML Online—a web server for fast maximum likelihood-based phylogenetic inference［J］. Nucleic Acids Research，2005， 33（suppl 2）: W557-W559.

［17］ Abascal F， Zardoya R， Posada D. ProtTest: selection of best-fit models of protein evolution［J］. Bioinformatics，2005， 21（9）: 2104-2105.

［18］ Bailey T L， Boden M， Buske F A， et al. MEME SUITE: tools for motif discovery and searching［J］. Nucleic acids research， 2009， 37: W202-W208.

［19］ Conesa A， G?tz S， García-Gómez J M， et al. Blast2GO: a universal tool for annotation， visualization and analysis in functional genomics research［J］. Bioinformatics， 2005， 21（18）: 3674-3676.

［20］ Ye J， Fang L， Zheng H， et al. WEGO: a web tool for plotting GO annotations［J］. Nucleic Acids Research，2006， 34（suppl 2）: W293-W297.

［21］ Birkenbihl R P， Jach G， Saedler H， et al. Functional dissection of the plant-specific SBP-domain: overlap of the DNA-binding and nuclear localization domains［J］. Journal of Molecular Biology， 2005， 352（3）: 585-596.

［22］ Riese M， H?hmann S， Saedler H， et al. Comparative analysis of the SBP-box gene families in P. patens and seed plants［J］. Gene， 2007， 401（1）: 28-37.

［23］ 朱紅霞，胡利宗，鄧小莉，等. 大豆SBP基因家族的序列特征，表達(dá)及進(jìn)化分析［J］. 東北農(nóng)業(yè)大學(xué)學(xué)報(bào)，2012，43（7）：26-33.

Genome-Wide Identification， Phylogenetic Analysis and Expression Profiling of the SBP Transcription Factor Family in Nicotiana tabacum

LIU Cheng1，2， LI Xiaoxu1，2， SU Yulong1，2， GUO Yongfeng1*
（1. Tobacco Research Institute of Chinese Academy of Agricultural Sciences， Qingdao 266101， China； 2. Graduate School of Chinese Academy of Agricultural Sciences， Beijing 100081， China）

The members of the SBP transcription factor family have been reported to play significant roles in regulating flower and fruit development as well as other biological processes. In this study， 32 SBP transcription factor family members were identified in Nicotiana tabacum with HMMER. Both Maximum-Likelihood phylogenetic tree and neighbor-joining tree were constructed and were found to possess similar topologies using the protein sequences of the SBP-domain. Results of phylogenetic analysis revealed that the SBP transcription factor family members could be classified into 8 subfamilies. The patterns of exon-intron structure and conserved domains in Arabidopsis and tobacco were consistent with the phylogenetic results. Transcriptome analysis showed that the expression patterns of NtSPLs were different in different tissue types and the expression level of NtSPLs in flower was found to be significantly higher than the other tissues. GO analysis suggested that as transcription factors， the SBP family members could be involved in a series of biological processes such as developmental regulation and defense response. The results of this study provide insight into the evolution of the SBP transcription factor family in Nicotiana tabacum and provide useful information for future research on these genes.

bioinformatics analysis； SBP transcription factor family； Nicotiana tabacum

S572.01

1007-5119（2015）04-0001-11

10.13496/j.issn.1007-5119.2015.04.001

農(nóng)業(yè)部948項(xiàng)目“高通量植物衰老基因克隆技術(shù)（TASSEL-tagging系統(tǒng)）的引進(jìn)、創(chuàng)新及利用”（2013-Z4）

劉 成（1988-），男，碩士研究生，研究方向?yàn)闊煵莘肿佑N。E-mail：liuch1221@163.com。*通信作者，E-mail：guoyongfeng@caas.cn

2015-03-20

2015-07-20