束紅梅,郭書巧,鞏元勇,蔣 璐,朱靜雯,倪萬潮
(江蘇省農(nóng)業(yè)科學院 經(jīng)濟作物研究所,農(nóng)業(yè)部長江下游棉花與油菜重點實驗室,江蘇 南京 210014)
陸地棉油菜素內(nèi)酯信號基因GhBES1家族的鑒定及表達分析
束紅梅,郭書巧,鞏元勇,蔣 璐,朱靜雯,倪萬潮
(江蘇省農(nóng)業(yè)科學院 經(jīng)濟作物研究所,農(nóng)業(yè)部長江下游棉花與油菜重點實驗室,江蘇 南京 210014)
BES1基因是植物激素油菜素內(nèi)酯信號轉(zhuǎn)導(dǎo)過程中的重要轉(zhuǎn)錄因子,為了解陸地棉油菜素內(nèi)酯信號基因GhBESI家族 ,通過對陸地棉異源四倍體標準系TM-1基因組數(shù)據(jù)庫的全面分析和鑒定,獲得21個GhBES1基因,這些基因分布在不同的亞基因組,有9對基因在A亞組和D亞組上存在對應(yīng)關(guān)系,AD亞組對應(yīng)基因序列一致性高于95%。根據(jù)進化分析,21個GhBES1基因分為4個亞家族,除了亞家族Ⅲ,其他3個亞家族在擬南芥中均有對應(yīng)基因;除了亞家族Ⅳ基因和亞家族Ⅱ中GhBES1-1基因外,其他GhBES1基因均有2個外顯子。亞細胞定位結(jié)果表明,21個GhBES1蛋白主要定位到細胞核、細胞質(zhì)等部位,其中7對基因的蛋白亞細胞定位完全一致。GhBES1基因在根、莖、葉、頂端分生組織中均有表達,但不同亞家族基因在不同組織中的表達水平存在差異。研究結(jié)果為棉花GhBES1基因家族功能的深入研究奠定了基礎(chǔ)。
陸地棉;GhBES1基因家族;系統(tǒng)進化分析;基因表達
油菜素內(nèi)酯(Brassinosteroid,BR)是植物正常生長發(fā)育所必需的重要植物激素,同時也參與了植物的逆境響應(yīng)[1-3]。近年來,關(guān)于油菜素內(nèi)酯信號轉(zhuǎn)導(dǎo)研究已取得了突破性的進展[4-5],在模式植物擬南芥中已經(jīng)明確油菜素內(nèi)酯信號通路[6-9]。油菜素內(nèi)酯從受體BRI1傳遞信號給核靶基因依賴于BES1(BRI1 EMS SUPPRESSOR 1)家族的轉(zhuǎn)錄因子。
BES1是在功能增強型突變體BES1-D中發(fā)現(xiàn)的[10-11],油菜素內(nèi)酯信號的傳遞最終使去磷酸化狀態(tài)的BES1/BZR1在細胞內(nèi)累積,激活油菜素內(nèi)酯信號通路下游的轉(zhuǎn)錄調(diào)控[10]。在植物轉(zhuǎn)錄因子數(shù)據(jù)庫中近80個物種存在BES1轉(zhuǎn)錄因子家族,但是每個物種BES1家族的基因數(shù)目存在差異,在模式植物擬南芥中BES1轉(zhuǎn)錄因子家族有8個成員[11-13]。研究人員通過ChIP-chip鑒定了BES1轉(zhuǎn)錄因子的靶基因,Yu等[14]鑒定的1 609個BES1靶基因中至少有250個是受油菜素內(nèi)酯調(diào)控的,這些基因參與調(diào)控生長發(fā)育、細胞信號轉(zhuǎn)導(dǎo)、環(huán)境脅迫響應(yīng)以及與其他激素信號通路間的交互等過程。而且,元件序列分析證明BES1有非典型的螺旋-環(huán)-螺旋(bHLH)DNA結(jié)合結(jié)構(gòu)域可以綁定到E-盒(CANNTG)或BRRE(BR response element,CGTGT/CG)以調(diào)節(jié)靶基因的表達[15-16],而BRRE和E-盒在油菜素內(nèi)酯誘導(dǎo)或抑制基因的啟動子中都存在。因此,BES1轉(zhuǎn)錄因子在油菜素內(nèi)酯調(diào)控植物生長發(fā)育和逆境響應(yīng)過程中起著重要作用。
陸地棉(Gossypiumhirsutum)是當今最重要的棉花栽培品種。二倍體D亞組雷蒙德氏棉和A亞組亞洲棉測序工作的完成[17-18],以及陸地棉遺傳標準系TM-1全基因組序列的公布[19-20],為棉花分子生物學研究奠定了良好的基礎(chǔ)。然而,目前在基因組水平上關(guān)于棉花BES1基因家族的研究較少。因此,本研究利用生物信息學方法,分析鑒定陸地棉全基因組序列中的BES1家族基因,從基因組水平上分析BES1基因在陸地棉中的數(shù)目、序列特征、進化與組織表達特征,為進一步研究BES1家族成員的生物學功能奠定理論基礎(chǔ)。
1.1 試驗材料
陸地棉異源四倍體標準系TM-1基因組數(shù)據(jù)來自南京農(nóng)業(yè)大學作物遺傳與種質(zhì)作物創(chuàng)新重點實驗室(http://mascotton.njau.edu.cn/html/Data/Genomefhsequence/2015/05/05/16ab0945-19e9-49f7-a09e-8e956866bf.html)。
1.2GhBES1家族成員鑒定及序列分析
用HMMER及Blast程序搜索鑒定陸地棉全基因組中GhBES1家族的所有成員。用Perl程序從基因組數(shù)據(jù)庫文件Gossypium_hirsutum_v1.1.cds.fa和Gossypium_hirsutum_v1.1.pep.fa中提取GhBES1基因的CDS序列和氨基酸序列。
利用ProtParam工具(http://web.expasy.org/protparam/)在線分析獲得陸地棉GhBES1蛋白氨基酸序列的分子量、等電點等信息。
1.3BES1基因的進化和結(jié)構(gòu)分析
用GSDS(Gene structure display server,http://gsds.cbi.pku.edu.cn/index.php)繪制BES1基因結(jié)構(gòu)圖。
用DNAMAN軟件,對陸地棉BES1基因家族成員進行序列同源性的比對,在氨基酸水平上分析他們的遺傳多樣性;使用ClustalX(1.83)軟件對氨基酸序列進行多序列聯(lián)配,將多序列聯(lián)配的結(jié)果通過Mega 5.0軟件進行分析,構(gòu)建了系統(tǒng)進化樹,采用鄰近法(Neighbor-joining,NJ)并利用Bootstrap方法進行評估,數(shù)值設(shè)置為1 000。
1.4 亞細胞定位預(yù)測
利用WoLF PSORT軟件(http://www.genscript.com/wolf-psort.html)分析氨基酸序列,進行亞細胞定位預(yù)測。
1.5GhBES1基因家族成員的表達分析
用OMIGA公司植物RNA提取試劑盒提取陸地棉根、莖、葉、頂端分生組織的總RNA。以梯度總RNA量及梯度循環(huán)數(shù)確定可區(qū)分表達強度差異的最佳PCR循環(huán)數(shù)。以等量總RNA進行RT-PCR,所用引物見表1,擴增產(chǎn)物用1.2%瓊脂糖凝膠電泳分離。反轉(zhuǎn)錄所用生化酶等購自Promega公司,PCR所用酶等購自TaKaRa公司,RNA提取過程中所需槍頭、離心管等購自Axygen公司。
表1 引物Tab.1 PCR primers
2.1 陸地棉GhBES1基因家族鑒定及序列分析
根據(jù)已公布的BES1基因的結(jié)構(gòu)域DUF822,從陸地棉異源四倍體標準系TM-1基因組數(shù)據(jù)中共鑒定出21個GhBES1基因,根據(jù)這些基因在染色體出現(xiàn)的先后順序分別命名為GhBES1-1~GhBES1-21 (表1)。對鑒定出來的基因進行分析,在A亞組共有10個GhBES1基因,在D亞組共有11個GhBES1基因,其中A亞組中9個基因在D亞組有對應(yīng)基因。并且在A05和D05、A11和D11這2對亞基因組中有2對GhBES1基因,其他成對的亞基因組只有1對平行同源進化基因(表2、圖1)?;蛐蛄幸恢滦苑治霭l(fā)現(xiàn)(圖1),A、D亞組相對應(yīng)的9對基因的序列一致性高于95%。但其他GhBES1基因序列的一致性差異較大,有些基因的序列一致性只有38%。
表2 陸地棉GhBES1基因家族基本信息Tab.2 Basic information of GhBES1 gene family of upland cotton
這21個GhBES1基因的編碼區(qū)和CDS最長的是GhBES1-19,最短的是GhBES1-13。除了GhBES1-3和GhBES1-13、GhBES1-9和GhBES1-19這2對基因外,其余8對基因的CDS序列長度完全一致。蛋白質(zhì)理化性質(zhì)分析表明,21個GhBES1蛋白的氨基酸序列個數(shù)從176(GhBES1-13)~737(GhBES1-19) 不等,其中有15個蛋白的氨基酸序列集中在320個氨基酸左右;分子量最大的是83.036 kDa (GhBES1-19),最小的是19.740 kDa(GhBES1-13);理論等電點最大的是9.26(GhBES1-13),最小的是5.39(GhBES1-17)(表2)。
2.2 陸地棉GhBES1家族基因結(jié)構(gòu)及系統(tǒng)進化分析
基于陸地棉GhBES1基因的蛋白全長序列構(gòu)建了系統(tǒng)進化樹(圖2)。對陸地棉中21個GhBES1基因和擬南芥中8個AtBES1基因進行遺傳多樣性分析,結(jié)果表明:BES1基因家族成員之間同源性變異幅度較大。按照進化情況陸地棉中21個GhBES1基因分為4個亞家族,除了亞家族Ⅲ外,每個陸地棉GhBES1亞家族在擬南芥中均有對應(yīng)基因。亞家族Ⅰ中有6個GhBES1基因,對應(yīng)擬南芥中的AtBES1、AtBZR1、AtBEH1、AtBEH2;亞家族Ⅱ中有8個GhBES1基因,對應(yīng)擬南芥中的AtBEH3、AtBEH4;亞家族Ⅳ中有5個GhBES1基因,對應(yīng)擬南芥中的AtBMY2、AtBMY7。
對陸地棉GhBES1家族成員的基因結(jié)構(gòu)分析表明(圖2),GhBES1基因家族大多數(shù)成員的結(jié)構(gòu)比較簡單,其中亞家族Ⅰ、Ⅱ、Ⅲ中基因除GhBES1-1外均含有2個外顯子,亞家族Ⅳ的5個GhBES1基因含有5~12個外顯子。從基因的外顯子長度(圖2、表3)看出:在10對基因中,有7對A、D亞基因組對應(yīng)的基因外顯子個數(shù)及片段大小完全一致。GhBES1-1和GhBES1-11外顯子個數(shù)相差1個,但長度相同;GhBES1-3和GhBES1-13外顯子個數(shù)相同,但長度存在差異;GhBES1-9和GhBES1-19的差異較大,分別含有5,12個外顯子,后者的外顯子序列長度明顯長于前者,而且后者的內(nèi)含子序列長度達到6 842 bp。
圖1 GhBES1基因序列的一致性分析Fig.1 Analysis of sequence identities among GhBES1 genes
2.3 陸地棉GhBES1蛋白的亞細胞定位
陸地棉TM-1基因組中21個GhBES1蛋白成員的亞細胞定位預(yù)測結(jié)果表明(表2),亞家族Ⅱ中8個GhBES1蛋白均只有1個亞細胞定位,定位到細胞核;亞家族Ⅰ中的6個蛋白和GhBES1-3、GhBES1-7、GhBES1-17共9個成員具有2個亞細胞定位,定位到細胞核、線粒體或細胞質(zhì);亞家族Ⅳ中GhBES1-9、GhBES1-19、GhBES1-21和亞家族Ⅲ中GhBES1-13 4個成員具有3個或3個以上亞細胞定位,定位到細胞核、葉綠體、細胞質(zhì)等。9對A、D亞組對應(yīng)基因中,有7對基因具有相同的亞細胞定位。
2.4 陸地棉GhBES1基因的組織表達分析
從陸地棉GhBES1基因的4個亞家族選擇代表基因(表3中加粗的基因)進行組織表達分析(圖3)。所選的5個代表基因在根、莖、葉、頂端分生組織中均有表達,其中GhBES1-1(亞家族Ⅱ)在根中表達量較高,GhBES1-2(亞家族Ⅰ)在葉中的表達量略低,GhBES1-3(亞家族Ⅲ)在根中的表達量較低,GhBES1-7和GhBES1-21同屬于亞家族Ⅳ,兩者在各個組織中的表達水平差異不明顯。可以看出不同亞家族GhBES1基因在各個器官中的表達有差異,推測這些GhBES1基因在各個器官中的功能可能也存在差異。
圖2 BES1基因家族的系統(tǒng)進化樹和基因結(jié)構(gòu)Fig.2 The phylogenetic tree and gene structure of BES1 genes
表3 陸地棉21個GhBES1基因外顯子的比較Tab.3 Comparison of exon of 21 GhBES1 genes of cotton bp
D.頂端分生組織;Y.葉片;J.莖;G.根。D.Apical meristem.Y.Leaf;J.Stem;G.Root.
BES1基因是油菜素內(nèi)酯信號轉(zhuǎn)導(dǎo)過程中的重要基因[14-16],但關(guān)于其在陸地棉中的情況研究較少。本研究在陸地棉異源四倍體標準系TM-1基因組數(shù)據(jù)中發(fā)現(xiàn)有21條GhBES1基因的家族成員,其中9對為A亞組和D亞組對應(yīng)基因,這與前人相關(guān)研究認為二倍體A 組亞洲棉和D組雷蒙德氏棉基因組之間存在很強的共線性的結(jié)論相吻合[20-21];9對共線性基因中,7對共線性基因的氨基酸長度、結(jié)構(gòu)完全相同,并且這7對基因具有完全相同的亞細胞定位。但是除了共線性基因,其他GhBES1基因的序列相似度較低。說明陸地棉21個GhBES1基因中有共線性基因的存在,但也有部分基因序列差異明顯。
擬南芥BES1基因家族有8個成員,其中部分成員的功能都已研究得較為清楚[7,11,22-24]。對陸地棉和擬南芥中的所有BES1進行進化樹分析發(fā)現(xiàn),在陸地棉標準系TM-1基因組中鑒定的21個GhBES1基因家族成員分為4個亞家族,除了亞家族Ⅲ,其余每個亞家族在擬南芥中有同源基因,由此推測,棉花和擬南芥BES1基因家族成員在功能上可能具有一定的相似性。因此,在研究棉花GhBES1基因家族成員功能時,在一定程度上可以參考擬南芥中BES1基因的研究結(jié)果。
組織表達分析表明,TM-1全基因組中的GhBES1基因家族的代表基因在根、莖、葉和頂端分生組織中均有表達,這在一定程度上暗示棉花GhBES1基因家族多數(shù)成員可能在多個組織的生長發(fā)育過程中起作用,但是不同亞家族的GhBES1基因在不同組織中的表達存在差異,說明這些BES1基因的功能可能也存在差別。
前人研究表明,油菜素內(nèi)酯信號基因BES1通過介導(dǎo)調(diào)控下游基因的表達進而調(diào)控植物的生長發(fā)育和逆境響應(yīng)[7,11,22-24]。在本研究中陸地棉21個GhBES1基因中部分基因的序列結(jié)構(gòu)、組織表達存在差異,那么它們在植物生長發(fā)育、逆境響應(yīng)中哪些途徑起作用、功能存在哪些差異,還有待深入研究。
[1] Bajguz A.Brassinosteroid enhanced the level of abscisic acid inChlorellavulgarissubjected to short-term heat stress [J].Journal of Plant Physiology,2009,166(8):882-886.
[2] 束紅梅,郭書巧,鞏元勇,等.油菜素內(nèi)酯對NaCl脅迫下棉花葉片生理特征和基因表達譜的影響 [J].應(yīng)用生態(tài)學報,2016,37(1):150-156.
[3] Divi U K,Rahman T,Krishna P.Brassinosteroid-mediated stress tolerance inArabidopsisshows interactions with abscisic acid,ethylene and salicylic acid pathways [J].BMC Plant Biology,2010,10(1):151.
[4] Clouse S D. Plant Hormones: Brassinosteroid signal transduction and action[M].Netherlands: Springer,2010.
[5] Maharjan P M,Schulz B,Choe S.BIN2/DWF12 antagonistically transduces brassinosteroid and auxin signals in the roots ofArabidopsis[J].J Plant Biology,2011,54(2):126-134.
[6] Clouse S D.Brassinosteroid signal transduction:clarifying the pathway from ligand perception to gene expression [J].Molecular Cell,2002,10(5):973-982.
[7] Kang S,Yang F,Li L,et al.TheArabidopsistranscription factor BES1 is a direct substrate of MPK6 and regulates immunity[J].Plant Physiology,2015,167(3):1076-1086
[8] Belkhadir Y,and Chory J.Brassinosteroid signaling:a paradigm for steroid hormone signaling from the cell surface [J].Science,2006,314(5804):1410-1411.
[9] Li J,Jin H.Regulation of brassinosteroid signaling [J].Trends in Plant Science,2007,12(12):37-41.
[10] 衛(wèi)卓赟,黎 家.受體激酶介導(dǎo)的油菜素內(nèi)酯信號轉(zhuǎn)導(dǎo)途徑[J].生命科學,2011,23(11):1106-1113.
[11] Jiang J,Zhang C,Wang X.A recently evoloed isoform of the transcription factor BES1 promotes brassinosteroid signaling and development inArabidopsisthaliana[J].Plant Cell,2015,27(2):361-374
[12] Wang Z Y,Nakano T,Gendron J,et al.Nuclear-localized BZR1 mediates brassinosteroid-induced growth and feedback suppression of brassinosteroid biosynthesis [J].Developmental Cell,2002,2(4):505-513.
[13] Zhao J,Peng P,Schmitz R J,et al.Two putative BIN2 substrates are nuclear components of brassinosteroid signaling [J].Plant Physiology,2002,130(3):1221-1229.
[14] Yu X,Li L,Zola J,et al.A brassinosteroid transcriptional network revealed by genome-wide identification of BES1 target genes inArabidopsisthaliana[J].Plant Journal,2011,65(4):634-646.
[15] Yin Y,Vafeados D,Tao Y,et al.A new class of transcription factors mediate brassinosteroid-regulated gene expression inArabidopsis[J].Cell,2005,120(2):249-259.
[16] He J X,Gendron J M,Sun Y,et al.BZR1 is a transcriptional repressor with dual roles in brassinosteroid homeostasis and growth responses [J].Science,2005,307(5715):1634-1638.
[17] Wang K,Wang Z,Li F,et al.The draft genome of a diploid cottonGossypiumraimondii[J].Nature Genetics,2012,44(10):1098-1103.
[18] Li F,F(xiàn)an G,Wang K,et al.Genome sequence of the cultivated cottonGossypiumarboreum[J].Nature Genetics,2014,46(6):562-574.
[19] Li F,F(xiàn)an Z,Lu C,et al.Genome sequence of cultivated Upland cotton (GossypiumhirsutumTM-1) provides insights into genome evolution [J].Nature Biotechnology,2015,33(5):524-530.
[20] Zhang T,Hu Y,Jiang W,et al.Sequencing of allotetraploid cotton (GossypiumhirsutumL.acc.TM-1) provides a resource for fiber improvement [J].Nature Biotechnology,2015,33(5):531-537.
[21] 倪萬潮,鞏元勇,徐珍珍,等.陸地棉GhZIP基因家族全基因組分析[J].華北農(nóng)學報,2015,30(6):8-16.
[22] Li L,Yu X,Thompson A,et al.ArabidopsisMYB30 is a direct target of BES1 and cooperates with BES1 to regulate brassinosteroid-induced gene expression [J].Plant Journal,2009,58(2):275-286.
[23] Guo R,Qian H,Shen W,et al. BZR1 and BES1 participate in regulation of glucosinolate biosynthesis by brassinosteroids inArabidopsis[J].Journal of Experimental Botany,2013,64(8):2401-2412.
[24] Li L,Ye H X,Guo H Q,et al.ArabidopsisIWS1 interacts with transcription factor BES1 and is involved in plant steroid hormone brassinosteroid regulated gene expression [J].Proceeding of the National Academy of Sciences of the United States of America,2010,107(8):3918-3923.
Identification and Expression Analysis of the Brassinosteroid Signal GeneGhBES1 Family in Upland Cotton
SHU Hongmei,GUO Shuqiao,GONG Yuanyong,JIANG Lu,ZHU Jingwen,NI Wanchao
(Institute of Industrial Crops,Jiangsu Academy of Agricultural Sciences,Key Laboratory of Cotton and Rapeseed (Nanjing),Ministry of Agriculture,Nanjing 210014,China)
BES1 gene is a major transcription factor responsible for brassinosteroid-regulated gene expression and plant growth.In order to know the family ofGhBES1 gene in upland cotton,21GhBES1 genes were identified from the genome sequence database of upland cotton (GossypiumhirsutumL.acc.TM-1).TheseGhBES1 genes were found to be distributed on different subgenomes.There were 9 pairs ofGhBES1 genes in the A and D subgenomes with one-to-one correspondence relationship,and the gene sequence identity of the correspondence genes were higher than 95%.Based on the results of phylogenetic tree,the 21GhBES1 genes could be divided into four subfamilies,and there were homologous genes ofArabidopsisthalianaL.in each subfamily,except for subfamily Ⅲ;GhBES1 genes had two exons, except for genes in subfamily Ⅳ andGhBES1-1 in subfamily Ⅱ.The result of subcellular localization indicated that 21 GhBES1 proteins were primarily located in nuclear and cytoplasmic,and 7 pairs ofGhBES1 genes had the same location.The tissue expression analysis showed thatGhBES1 genes were expressed in the leaf,root,stem and shoot apical meristem,but the expression level ofGhBES1 genes in different subfamily were different.These results provide a theoretical basis on cloning and functional research ofGhBES1 genes.
Upland cotton;GhBES1 gene family; Phylogenetic analysis;Gene expression
2017-04-16
國家自然科學基金項目(31201139);國家轉(zhuǎn)基因生物新品種培育項目(2016ZX08005-001)
束紅梅(1982-),女,江蘇射陽人,副研究員,博士,主要從事農(nóng)業(yè)生物技術(shù)研究。
倪萬潮(1962-),男,江蘇盱眙人,研究員,碩士,主要從事農(nóng)業(yè)生物技術(shù)研究。
S562;Q78
A
1000-7091(2017)04-0007-06
10.7668/hbnxb.2017.04.002