彭华^1,*, 何秀静^2,*, 高健⁴, 罗茂³, 潘光堂^2,*, 张志明^2,*
1四川旅游学院, 四川成都 610100
²四川农业大学玉米研究所, 四川温江 611130
³四川医科大学药物与功能性食品研究中心, 四川泸州 646000
⁴第三军医大学西南医院病理学研究所西南癌症中心, 重庆 400038
^*通讯作者(Corresponding authors): 潘光堂, 张志明, E-mail: panlab605@gmail.com, Tel: 028-86290917 第一作者联系方式: E-mail: 16828333@qq.com ^**同等贡献(Contributed equally to this work)
收稿日期:2015-05-12 接受日期:2015-11-20网络出版日期:2015-12-07基金:本研究由国家高技术研究发展计划(863计划)项目(2012AA10A307)和四川省科技厅青年基金项目(2015JQO021)资助

摘要 SBP基因家族是一类植物基因组特有的转录因子, 参与植物生长发育及多种生理生化过程。近来, 大量研究已在多种植物中鉴定出 SBP转录因子, 但关于玉米( Zea mays L.) SBP转录因子家族的系统分析报道尚少。本研究通过对拟南芥、水稻等植物已知的转录因子与玉米基因组数据比对, 并设置一系列严格的筛选标准从玉米基因组中挖掘 SBP转录因子, 系统发育分析显示单子叶植物玉米和水稻的 SBP基因保守性更强、亲缘关系更近; 共鉴定37个 SBP基因, 分布在9条染色体上; 通过基因分析注释以及启动子功能预测, 进一步发现 SBP家族基因参与植物生长发育、形态建成、逆境胁迫响应、花器官发育以及植物光反应等过程。并且, 玉米 SBP转录因子可通过参与赤霉素、生长素、脱落酸、水杨酸等多条激素信号调控途径来调节植物的生长发育。

关键词:玉米; SBP转录因子家族; 生物信息学
Genome-wide Identification and Function Analysis of SBP Gene Family in Maize
PENG Hua^1,**, HE Xiu-Jing^2,**, GAO Jian⁴, LUO Mao³, PAN Guang-Tang^2,*, ZHANG Zhi-Ming^2,*
1 Sichuan Tourism University, Chengdu 610100, China;

² Maize Research Institute of Sichuan Agricultural University, Wenjiang 611130, China;

³Research Center for Drug Discovery of Luzhou Medical College, Luzhou 646000, China;

⁴Institute of Pathology and Southwest Cancer Center, Southwest Hospital, Third Military Medical University, and Key Laboratory of Tumor Immunopathology, Ministry of Education, Chongqing 400038, China

Fund:This study was supported by the National High Technology Research and Development Program of China (863 Program) (2012AA10A307) and the grants for Young Scientists from Science and Technology Department of Sichuan Province (2015JQO021)
Abstract SBP gene family, as a plant special transcription factor, is involved in plant growth and development, as well as many physiological and biochemical processes. Recently, SBP transcription factor family has been identified in model plants, such as Arabidopsis and Oryza sativa; however, systematic analysis of SBP transcription factor family in maize ( Zea mays L.) is scarcely. In this study, based on homology alignment technology, we aligned all known SBP TFs from Arabidopsis and Oryza sativa with those from maize genome sequence to mine novel SBP TFs in maize. A total of 37 SBP TFs distributed in eight chromosomes were identified. Phylogenetic analysis indicated that SBP transcription factor genes have stronger homology, especially between Zea maysand Oryza sativa. Moreover, promoters- ciselements analysis of those SBP TFs demonstrated that they might be involved in plant growth and development, morphogenesis, adversity response, the development of flower organs and photosynthesis. It is probable that SBP TFs regulate plant growth and development through the multiple hormone of signaling transduction pathway, such as gibberellin, auxin, abscisic acid, and salicylic acid.

Keyword:Maize; SBP TFs gene family; Bioinformatics
Show Figures
Show Figures














目前, SBP转录因子已在大量植物中被鉴定, 具有多种功能, 在植物的生长发育过程中发挥重要调控作用。例如, 拟南芥SBP8和SBP14基因分别与调控花粉发育和真菌毒素FB1诱导的程序性死亡的抗性相关^[4], 拟南芥SPL3基因的转录受发育调控, 主要集中在花序顶端分生组织、花分生组织和花器官原基等, 可能与控制开花时间有关^[5]。水稻基因组SBP转录因子主要在花和愈伤组织表达^[6]。番茄SBP基因(LeSPL-CNR)是控制果实发育的关键基因^[7]。玉米SBP转录因子lg1与玉米舌叶和叶耳的发育密切相关, 若该基因发生突变, 植株将不能形成正常形态的舌叶和叶耳组织^{[8, 9]}。Andrea等^[10]发现玉米SBP转录因子tsh4与玉米雄花分枝形态构成有密切关系。Eveland等研究表明玉米SBP转录因子tga1调控玉米果穗的颖片发育, 是玉米进化过程中的一个关键基因, 同时Eveland等和Lä nnenpä ä 等^{[11, 12]}研究揭示玉米SBP转录因子家族成员果实构造基因tga1控制有麸玉米的种皮形态特征。Chuck等^[13]研究显示, 玉米SBP转录因子ub2和ub3基因可以通过控制玉米雌雄花序的分支和分化影响玉米的产量。研究进一步提示, ub2和ub3双突变玉米植株的雄花分支数减少, 花药数量减少, 散粉量严重减少, 同时玉米雌穗穗长变短、穗行数增加, 产量下降^[13]。另外, 当tsh4基因也缺失时, 玉米产量下降更严重。Chuck等^[14]研究显示, 玉米Cg1基因受Zma-miR156b/c的调控, 与玉米开花期SBP-box基因表达密切相关, 提示Cg1基因能够通过抑制玉米SBP-box基因表达, 进而调控玉米幼嫩叶片发育和雄花分枝。
目前, 大量研究陆续从在拟南芥(Arabidopsis thaliana)、水稻(Oryza sativa)、葡萄(Vitis amurensis)、苹果(Malus pumilaMill)、大豆[Glycine max(Linn.) Merr]、草莓(Fragaria ananassa)、番茄(Lycopersicon esculentumMiller)等植物基因组中鉴定出SBP基因^{[3, 6, 15, 16, 17, 18, 19]}。随着玉米基因组测序的完成为加快玉米的遗传育种和分子生物学研究提供了机遇, 然而, 关于玉米SBP基因的信息仍然较少。本研究利用生物信息学方法, 在玉米全基因水平鉴定SBP转录因子, 同时结合拟南芥、水稻SBP家族基因成员, 分析植物中SBP基因的进化关系及其表达模式, 研究结果不仅有助于鉴定玉米SBP基因家族的功能, 还可为以后玉米育种过程提供重要的理论依据。
1 材料与方法1.1 材料水稻和拟南芥SBP蛋白序列来源于植物转录因子数据PlnTFDB网站(http://plntfdb.bio.uni-potsdam. de/v3.0/)。
1.2 玉米SBP转录因子的鉴定 从蛋白质家族数据库Pfam网站(http://pfam. sanger.ac.uk/)下载SBP基因家族的隐马氏模型序列谱(PF03110)文件, 应用HMMER 3.0程序包中的hmmbuild 命令将其转化成隐马氏模型序列模式。首先应用GENSCAN将在NCBI下载的玉米基因组B73 V5b+参考序列基因组(http://www.maizesequence. org/v5b+)数据预测得到蛋白序列, 并利用HMMER 3.0软件中的hmmsearch对这些蛋白进行转录因子预测(参数为默认设置)^[20]。然后, 利用在线数据库SMART (http://smart.embl-heidelberg.de/)对候选SBP转录因子结构域进行鉴定, 确定真实的玉米SBP转录因子基因。利用玉米基因组mRNA数据库(http:// www.maizesequence.org/index.html)进行玉米SBP转录因子对应mRNA序列的比对分析(参数设置为BLASTp, 数据库基因集B73_RefGen_v3, E=1e-10, 其他参数默认); 收集所有阈值小于1.0的非冗余的数据集, 与PlnTFDB (http://plntfdb.bio.uni-potsdam. de/v3.0/)和PlantTFDB (http://planttfdb.cbi.pku.edu.cn) 中的SBP家族成员比对; 重新注释的序列通过其SMART (http://smart.embl-heidelberg.de/), CDD (http:// www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi)和Inter- ProScan (http://www.ebi.ac.uk/Tools/InterProScan/)来确认SBP家族结构域的存在, 进而查明玉米SBP转录因子mRNA序列。
1.3 玉米SBP转录因子的分类 利用ClustalX 2.1软件^[21]对玉米SBP转录因子的SBP保守域蛋白序列进行多序列比对, 结合比对结果与水稻、拟南芥SBP转录因子的分类依据, 进行玉米SBP转录因子的分类。
1.4 玉米SBP转录因子的系统发育关系 首先, 从玉米、拟南芥和水稻的SBP转录因子中提取完整SBP结构域(包括I组成员N端和C端的SBP结构域, 发生缺失的SBP结构域被舍弃)的蛋白序列以及完整的SBP结构域。然后, 利用ClustalX 2.1软件进行蛋白质多序列比对。利用Mega 5.0软件Neighbor-Joining算法构建系统进化树^{[22, 23]}, 设置JTT模型, 位点替换速率Gamma分布, Bootstrap值为1000。
1.5 玉米SBP 转录因子蛋白序列的氨基酸保守域分析 利用在线软件MEME 4.8.1 (http://meme.sdsc. edu/meme/cgi-bin/meme.cgi)对玉米SBP转录因子蛋白序列进行氨基酸保守域的预测分析。搜索得到的保守域数目上限为30; 保守域氨基酸残基的数量为10~300; 保守域可重复出现。
1.6 玉米SBP转录因子的表达模式分析 利用玉米转录组数据库qTeller (http://qteller. com/)中公布的RNAseq在各种组织中的表达结果分析玉米SBP转录因子(ZmSBP)表达模式, 并用MEV软件对所检测到的ZmSBP家族成员进行聚类分析。进一步利用PlantCARE (htp://bioinformatics.psb.ugent. be/webtols/plantcare/html/)网站预测候选ZmSBP可能包含的转录因子顺式作用元件, 预测ZmSBP功能。

2 结果与分析2.1 玉米SBP转录因子的鉴定与染色体定位 利用Pfam工具搜索和在线数据库SMART、CDD和Inter-ProScan比对鉴定, 共获得37个典型结构域的玉米SBP基因, 分别命名为ZmSBP1~ ZmSBP37 (表1)。如图1所示, 除第9染色体外, 玉米SBP基因被锚定于所有染色体, 几乎呈均匀分布状态; 第4和第5染色体上SBP基因数目分布最多, 分别有8个和7个家族成员; 第3、第6、第7染色体上SBP基因数目最少, 为2个成员; 其余染色体上包含SBP基因3~5个。
图1
Fig. 1

	Figure OptionViewDownloadNew Window
	图1 玉米SBP 基因的染色体定位Fig. 1 Chromosome mapping of SBP genes in maize

表1
Table 1
表1(Table 1)

表1 玉米SBP基因家族信息 Table 1 Information of SBP-domain gene family in maize 蛋白 Protein 编码序列 Coding sequence 名称 Name 重命名 Rename 染色体 Chromosome AC225342.3_FGP001 AC225342.3_FGT001 AC225342.3_FG001 ZmSBP1 5 AC233751.1_FGP002 AC233751.1_FGT002 AC233751.1_FG002 ZmSBP2 4 GRMZM2G017091_P02 GRMZM2G017091_T02 GRMZM2G017091 ZmSBP3 10 GRMZM2G020906_P02 GRMZM2G020906_T02 GRMZM2G020906 ZmSBP4 6 GRMZM2G036297_P01 GRMZM2G036297_T01 GRMZM2G036297 ZmSBP5 2 GRMZM2G058588_P01 GRMZM2G058588_T01 GRMZM2G058588 ZmSBP6 10 GRMZM2G061734_P01 GRMZM2G061734_T01 GRMZM2G061734 ZmSBP7 2 GRMZM2G065451_P02 GRMZM2G065451_T02 GRMZM2G065451 ZmSBP8 4 GRMZM2G067624_P02 GRMZM2G067624_T02 GRMZM2G067624 ZmSBP9 7 GRMZM2G081127_P02 GRMZM2G081127_T02 GRMZM2G081127 ZmSBP10 8 GRMZM2G097275_P04 GRMZM2G097275_T04 GRMZM2G097275 ZmSBP11 5 GRMZM2G098557_P01 GRMZM2G098557_T01 GRMZM2G098557 ZmSBP12 4 GRMZM2G101499_P01 GRMZM2G101499_T01 GRMZM2G101499 ZmSBP13 6 GRMZM2G101511_P01 GRMZM2G101511_T01 GRMZM2G101511 ZmSBP14 4 GRMZM2G106798_P03 GRMZM2G106798_T03 GRMZM2G106798 ZmSBP15 4 GRMZM2G108522_P01 GRMZM2G108522_T01 GRMZM2G108522 ZmSBP16 1 GRMZM2G109354_P01 GRMZM2G109354_T01 GRMZM2G109354 ZmSBP17 8 GRMZM2G111136_P02 GRMZM2G111136_T02 GRMZM2G111136 ZmSBP18 5 GRMZM2G113779_P01 GRMZM2G113779_T01 GRMZM2G113779 ZmSBP19 2 GRMZM2G126018_P02 GRMZM2G126018_T02 GRMZM2G126018 ZmSBP20 2 GRMZM2G126827_P01 GRMZM2G126827_T01 GRMZM2G126827 ZmSBP21 5 GRMZM2G133646_P01 GRMZM2G133646_T01 GRMZM2G133646 ZmSBP22 3 GRMZM2G138421_P01 GRMZM2G138421_T01 GRMZM2G138421 ZmSBP23 5 GRMZM2G148467_P01 GRMZM2G148467_T01 GRMZM2G148467 ZmSBP24 10 GRMZM2G156621_P01 GRMZM2G156621_T01 GRMZM2G156621 ZmSBP25 5 GRMZM2G156756_P01 GRMZM2G156756_T01 GRMZM2G156756 ZmSBP26 1 GRMZM2G160917_P01 GRMZM2G160917_T01 GRMZM2G160917 ZmSBP27 1 GRMZM2G163813_P03 GRMZM2G163813_T03 GRMZM2G163813 ZmSBP28 4 GRMZM2G168229_P01 GRMZM2G168229_T01 GRMZM2G168229 ZmSBP29 4 GRMZM2G169270_P02 GRMZM2G169270_T02 GRMZM2G169270 ZmSBP30 1 GRMZM2G179666_P01 GRMZM2G179666_T01 GRMZM2G179666 ZmSBP31 1 GRMZM2G307588_P01 GRMZM2G307588_T01 GRMZM2G307588 ZmSBP32 7 GRMZM2G371033_P01 GRMZM2G371033_T01 GRMZM2G371033 ZmSBP33 8 GRMZM2G414805_P04 GRMZM2G414805_T04 GRMZM2G414805 ZmSBP34 5 GRMZM2G418011_P01 GRMZM2G418011_T01 GRMZM2G418011 ZmSBP35 10 GRMZM2G460544_P01 GRMZM2G460544_T01 GRMZM2G460544 ZmSBP36 4 GRMZM5G878561_P01 GRMZM5G878561_T01 GRMZM5G878561 ZmSBP37 3

表1 玉米SBP基因家族信息 Table 1 Information of SBP-domain gene family in maize

2.2 玉米SBP转录因子的结构域比对及分类 利用Clustal X程序对玉米SBP蛋白家族氨基酸序列进行多序列比对表明(图2), 除ZmSBP16、ZmSBP31、ZmSBP1外, 各玉米SBP蛋白的SBP结构域包含约79个氨基酸残基并具备2个锌指结构, 分别为C3 (C-C-C-H)和C2HC (C-C-H-C)类型, 命名为Zn1和Zn2; 其中, ZmSBP16、ZmSBP31和ZmSBP1中只含有Zn1结构而Zn2缺失结构; 另外, 除ZmSBP35、ZmSBP3、ZmSBP4外, 所有转录因子在SBP保守结构域的C端都有一个NLS (nuclear localization signal)核定信号位点, 其他ZmSBP均具有完全保守的SBP结构域。
图2
Fig. 2

Figure OptionViewDownloadNew Window

图2 玉米SBP蛋白在不同亚家族中的比对分析 利用Clustal软件对玉米SBP蛋白家族氨基酸序列进行序列比对。各玉米SBP蛋白的SBP结构域包含2个锌指结构, 分别为Zn1和Zn2。蛋白序列中不同背景色代表各个位点的保守程度, 不同字母代表相应的氨基酸。Fig. 2 Alignment of multiple ZmSBP protein in different groups Multiple alignment of the SBP domains of the maize SBP-box proteins was obtained with Clustal software. The two conserved zinc finger structures (Zn1, Zn2) are indicated. In addition, sequence logo of the SBP domain of ZmSBPs was showed. The overall height of each stack represents the degree of conservation at this position, while the height of the letters within each stack indicates the relative frequency of the corresponding amino acids.

进一步依据SBP转录因子氨基酸序列的ClustalW多序列比对结果, 结合Meighbor-joining (J)和Bootstrap方法, 分析玉米ZmSBP转录因子蛋白的同源性和系统进化关系。结果显示(图3), 玉米SBP基因家族成员可以分为5种类型(I~V), 其中I类只有1个家族成员ZmSBP15, II和III类包含10个ZmSBP家族成员, IV和V类包含8个ZmSBP家族成员。
图3
Fig. 3

	Figure OptionViewDownloadNew Window
	图3 SBP转录因子家族基因的进化关系Fig. 3 Phylogenetic relationships of maize SBP family genes

2.3 玉米SBP转录因子的蛋白保守序列分析 利用MEME软件分析结果, 构建ZmSBP转录因子蛋白保守序列的结构特征图, 即motif分布图。结果表明(图4), ZmSBP转录因子蛋白保守序列主要包括4个motif结构, 分别为motif 1~4, 其中所有氨基酸位点都包括C2HC结构域(motif 1)和C3H结构域(motif 2)。进一步分析显示, ZmSBP1和ZmSBP31在进化过程中缺失motif 1结构, ZmSBP19缺失motif 2结构, ZmSBP17及ZmSBP22在进化过程中motif 3结构发生了2~3次复制; ZmSBP12和ZmSBP30在进化过程中motif 4结构发生了2次复制, 推测这些ZmSBP蛋白在进化过程中, 通过对功能区序列的复制产生了新的功能。
图4
Fig. 4

	Figure OptionViewDownloadNew Window
	图4 玉米SBP转录因子家族基因的基因结构及基序元件组成 基序结构图是利用MEME软件分析玉米SBP蛋白的保守位点获得, 图中黑线部分代表非保守序列, 各种颜色的盒子代表不同的基序元件。Fig. 4 Gene structure and motif compositons of ZmSBP genes Schematic representation of the conserved motifs in the SBP proteins from maize elucidated by MEME. A number in the colored box represents each motif. The black lines represent the non-conserved sequences.

2.4 玉米SBP转录因子的表达模式及功能分析利用玉米转录组数据库qTeller中公布的RNA-seq数据, 对玉米不同组织SBP转录因子的表达模式及功能分析。图5表明, ZmSBP10、ZmSBP15、ZmSBP22、ZmSBP28、ZmSBP23、ZmSBP26、ZmSBP12和ZmSBP30在成熟玉米根茎、叶、果穗、胚及胚乳中大量表达。大部分ZmSBP基因在胚珠、幼穗和未受精果穗中大量表达, 但在花粉、维管束鞘组织和叶肉细胞中表达量很低。其中, ZmSBP10是唯一一个在维管束鞘和叶肉细胞中大量表达的转录因子, 推测ZmSBP10在进化过程中具有不同于其他玉米SBP家族成员的功能。
图5
Fig. 5

	Figure OptionViewDownloadNew Window
	图5 玉米SBP转录因子家族基因表达模式分析 图中用数值(0~14.025)代表不同表达水平的基因, 数值越大表示基因表达水平越高。Fig. 5 Expression pattern of ZmSBP transcription factors gene in different organs Different expression levels are identified using number indexes (0-14.025), the larger numbers represent the higher expression levels.

利用在线软件PlantCARE (http://bioinformatics. psb.ugent.be/webtools/plantcare/html/)预测ZmSBP转录因子启动子包含的顺式作用元件, 进而预测ZmSBP转录因子成员的功能显示(表2), 37个转录因子在启动子区的各种功能元件有所不同, 除了最基本的TATA-box、CAAT-box外, 个别转录因子还包括如W-box、MBS、AuxRR-core等其他转录因子的结合位点。其中, W-box是AP2/EREBP类转录因子结合位点, 涉及植物生长发育, 如糖酵解和脂肪酸循环途径中关键酶的合成。同时W-box也是真菌诱导反应元件, 能够参与植物生物胁迫响应过程; MBS是MYB转录因子家族特有的结合位点, 推测其可能参与诱导抗性基因表达, 进而参与调控植物逆境胁迫响应; AuxRR-core一类是生长素响应因子, 在植物形态建成, 如胚乳、子叶和侧根发育中起着重要作用。另外, 还发现大量光刺激应答元件, 推测其可能与植物光合作用相关, 参与调控植物开花等光反应。
表2
Table 2
表2(Table 2)

表2 玉米SBP转录因子启动子区预测结果 Table 2 Forecasting result of promoter area of SBP transcription factors of maize

表2 玉米SBP转录因子启动子区预测结果 Table 2 Forecasting result of promoter area of SBP transcription factors of maize

2.5 植物SBP基因的进化发育关系 玉米、水稻和拟南芥的SBP基因系统发育树, 表明(图6) 3个物种的基因家族成员可分为5个类群(I~V), 其中第II类群仅包括水稻和玉米SBP基因家族成员, 其余4个类群均包括来自3个物种的SBP家族成员。另外, 进化发育树末端同一分支的两个外部结点很可能就是亲缘关系比较近的同源基因对。进一步对进化树末端分析显示, 水稻和玉米中存在大量直系同源基因, 如LOC_Os05g33810.1和ZmSBP17、LOC_Os05g44860.1和ZmSBP13、LOC_ Os01g69830.1和ZmSBP37等, 但是在玉米和拟南芥、水稻和拟南芥中未能发现直系同源基因, 提示单子叶植物的SBP转录因子家族进化模式存在相似性, 推测双子叶植物(如拟南芥) SBP转录因子家族的基本特征形成时间早于单子叶植物。此外, 在物种内还发现了大量SBP基因旁系同源基因, 如玉米ZmSBP29和ZmSBP18、ZmSBP5和ZmSBP6、ZmSBP36和ZmSBP27、ZmSBP32和ZmSBP20、ZmSBP2和ZmSBP14等均为旁系同源基因, 提示SBP基因家族在物种内以大量同源基因形式存在, 这些基因家族在物种中按照各自物种特异性方式扩展。
玉米、水稻和拟南芥的SBP基因系统发育树, 表明(图6) 3个物种的基因家族成员可分为5个类群(I~V), 其中第II类群仅包括水稻和玉米SBP基因家族成员, 其余4个类群均包括来自3个物种的SBP家族成员。另外, 进化发育树末端同一分支的两个外部结点很可能就是亲缘关系比较近的同源基因对。进一步对进化树末端分析显示, 水稻和玉米中存在大量直系同源基因, 如LOC_Os05g33810.1和ZmSBP17、LOC_Os05g44860.1和ZmSBP13、LOC_ Os01g69830.1和ZmSBP37等, 但是在玉米和拟南芥、水稻和拟南芥中未能发现直系同源基因, 提示单子叶植物的SBP转录因子家族进化模式存在相似性, 推测双子叶植物(如拟南芥) SBP转录因子家族的基本特征形成时间早于单子叶植物。此外, 在物种内还发现了大量SBP基因旁系同源基因, 如玉米ZmSBP29和ZmSBP18、ZmSBP5和ZmSBP6、ZmSBP36和ZmSBP27、ZmSBP32和ZmSBP20、ZmSBP2和ZmSBP14等均为旁系同源基因, 提示SBP基因家族在物种内以大量同源基因形式存在, 这些基因家族在物种中按照各自物种特异性方式扩展。
图6
Fig. 6

	Figure OptionViewDownloadNew Window
	图6 玉米、拟南芥和水稻SBP转录因子基因的系统进化树Fig. 6 Phylogenetic trees constructed with SBP domain protein sequences from maize (ZmSBP), rice (OsSPL), and Arabidopsis (AtSPL)

利用所有SBP基因的核心序列用于构建基因的系统进化树, 在图上标注每个SBP基因的序号、ID位点和数据来源。
SBP domain sequences of all genes were used for phylogenetic tree construction, the accession numbers or locus IDs and data sources of SBP-box genes are mached in the figure.

3 讨论随着植物基因组学研究的深入尤其是测序技术的迅猛发展, 已陆续完成大量植物全基因组序列鉴定, 转录因子全基因组的鉴定及表达模型研究逐渐成为当前植物基因功能研究的热点。近年来, 不同植物种类的SBP基因家族成员及大小已得到分析, 这些研究结果提示植物SBP转录因子家族成员在调节植物生长发育及多种生理生化过程中发挥重要作用。同时, 玉米基因组测序的完成为我们在全基因
组水平分析SBP转录因子奠定了基础。本研究基于玉米全基因组测序数据库, 在玉米全基因水平鉴定SBP转录因子并分析其表达功能, 共鉴定出37个SBP基因, 包括1个I型基因, 10个II型基因, 10个III型基因, 8个IV型基因, 8个V型基因。通过与水稻、拟南芥转录因子系统发育关系比对分析, 在拟南芥、水稻中分别发现了25个和27个SBP转录因子, 并且水稻和玉米中大量SBP基因具有直系同源关系, 而在玉米和拟南芥、水稻和拟南芥中未能检出, 提示与双子叶植物如拟南芥比较, 单子叶植物玉米和水稻中SBP基因保守性更强、亲缘关系更近。利用系统发育树将3类物种SBP基因分为5个类群(I~V), 单子叶植物的SBP转录因子家族进化模式存在相似起源方式, 推测双子叶植物(如拟南芥) SBP转录因子家族的基本特征形成时间早于单子叶植物。
SBP转录因子保守结构域的C端是核定位信号区域, 能够引导蛋白进入细胞核行使功能^[3]。本研究在分析SBP基序时发现, 玉米蛋白的结构域包含约79个氨基酸残基并具备2个锌指结构, 分别为C3 (C-C-C-H)和C2HC (C-C-H-C)类型, 推测其可能特异识别SQUAMOSA (SQ-UA)启动子。另外, ZmSBP转录因子蛋白保守序列分析结果显示, ZmSBP蛋白在进化过程中可能经历片段重复事件, 通过功能区序列复制而产生新的功能。
本研究进一步利用玉米转录组数据库所公布RNA-seq数据, 分析玉米SBP转录因子表达模式及功能, 结果显示, 大部分ZmSBP基因在玉米胚珠、幼穗和未受精果穗中大量表达, 但在玉米花粉、维管束鞘组织和叶肉细胞中表达量很低, 其中ZmSBP10是唯一一个在维管束鞘和叶肉细胞中大量表达的转录因子。研究表明, 植物花序发育与其产量密切相关, 如拟南芥同源基因SPL3、SPL4和SPL5能够调控拟南芥花发育, 拟南芥SBP8和SBP14基因能够调控其花粉发育^{[4, 5]}, 进而影响拟南芥产量。随着玉米SBP基因的揭示及功能研究迅猛发展, 大量研究显示玉米控制花序发育的部分SBP基因与玉米产量密切相关。例如, 玉米Cg1基因能够通过抑制玉米SBP基因表达, 调控玉米幼嫩叶片发育和雄花分枝^[14]; 玉米ub2和ub3基因能够通过控制玉米雌雄花序的分支和分化^[13], 调控玉米产量, 提示增强ub2和ub3基因表达可以增加玉米雄花分支和花药数, 进而增长玉米雌穗穗长, 提高玉米产量。部分研究揭示, 玉米tsh4能够调控玉米雄花分枝形态构成^[10], 提示增强tsh4基因表达, 能够有效增加玉米产量。
Chen等^[24]研究显示AtSPL8 (SPB转录因子)功能缺失会导致拟南芥花丝变短、萼片和表皮毛数量减少等半不育表型, 而外施赤霉素(GA)药剂, 其育性得到恢复, 这表明AtSPL8参与了GA信号转导途径, 推测SPB基因能通过参与调控赤霉素(GA)信号转导途径来影响植物的生长发育; Martin等^[25]研究表明SBP转录因子可以通过参与信号传导途径调控种子萌发、根系生长、下胚轴及茎的伸长、植物发育阶段转变、成花诱导以及花器官发育过程。本研究进一步预测ZmSBP转录因子启动子顺式作用元件并分析ZmSBP转录因子成员的功能结果显示, SBP家族基因与花发育等植物生长发育、形态建成、逆境胁迫以及植物光反应等相关。并且, 玉米SBP转录因子可通过参与赤霉素、生长素、脱落酸、水杨酸等多条激素信号调控途径来调节植物的生长发育。
近几年研究表明, SBP转录因子不仅参与玉米根、茎、叶、果穗、胚、胚乳、胚珠等器官的生长发育, 而且涉及多条激素信号调控途径, 提示SBP转录因子在植物生长发育各进程具有重要的作用, 然而SBP转录因子上游调控机制尚不清楚。本文基于玉米全基因组测序的完成, 鉴定分析玉米SBP的基因结构、进化关系以及表达模型与功能, 为未来深入解析玉米SBP转录因子在玉米生长发育和调节激素信号传递途径中作用提供基因资源和理论依据。
4 结论共鉴定分布在9条染色体上的37个SBP基因, 单子叶植物玉米和水稻的SBP基因保守性更强、亲缘关系更近; 玉米SBP家族基因参与植物生长发育、形态建成、逆境胁迫响应、花器官发育以及植物光反应等过程, 可通过参与多条激素信号调控途径来调节植物的生长发育。
The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.


参考文献View Option
原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1]	Huijser P, Klein J, Lönnig W E, Meijer H, Saedler H, Sommer H. Bracteomania, an inflorescence anomaly, is caused by the loss of function of the MADS-box gene squamosa in Antirrhinum majus. EMBO J, 1992, 11: 1239-1249[本文引用:1]
[2]	Yamasaki K, Kigawa T, Inoue M, Tateno M, Yamasaki T, Yabuki T, Aoki M, Seki E, Matsuda T, Nunokawa E, Ishizuka Y, Terada T, Shirouzu M, Osanai T, Tanaka A, Seki M, Shinozaki K, Yokoyama S. A novel zinc-binding motif revealed by solution structures of DNA-binding domains of Arabidopsis SBP-family transcription factors. J Mol Biol, 2004, 337: 49-63[本文引用:1]
[3]	Birkenbihl R P, Jach G, Saedler H, Huijser P. Functional dissection of the plant-specific SBP-domain: overlap of the DNA- binding and nuclear localization domains. J Mol Biol, 2005, 352: 585-596[本文引用:2]
[4]	Schmid M, Uhlenhaut N H, Godard F, Demar M, Bressan R, Weigel D, Lohmann J U. Dissection of floral induction pathways using global expression analysis. Development, 2003, 130: 6001-6012[本文引用:2]
[5]	Cardon G H, Höhmann S, Nettesheim K, Saedler H, Huijser P. Functional analysis of the Arabidopsis thaliana SBP-box gene SPL3: a novel gene involved in the floral transition. Plant J, 1997, 12: 367-377[本文引用:2]
[6]	Yang Z, Wang X, Gu S, Hu Z, Xu H, Xu C. Comparative study of SBP-box gene family in Arabidopsis and rice. Gene, 2008, 407: 1-11[本文引用:2]
[7]	Manning K, Tör M, Poole M, Hong Y, Thompson A J, King G J, Giovannoni J J, Seymour G B. A naturally occurring epigenetic mutation in a gene encoding an SBP-box transcription factor inhibits tomato fruit ripening. Nat Genet, 2006, 38: 948-952[本文引用:1]
[8]	Liu H, Yang X, Liao X, Zuo T, Qin C, Cao S, Dong L, Zhou H, Zhang Y, Liu S, Shen Y, Lin H, Lübberstedt T, Zhang Z, Pan G. Genome-wide comparative analysis of digital gene expression tag profiles during maize ear development. Genomics, 2015, 106: 52-60[本文引用:1]
[9]	Moreno M A, Harper L C, Krueger R W, Dellaporta S L, Freeling M. liguleless1 encodes a nuclear-localized protein required for induction of ligules and auricles during maize leaf organogenesis. Genes Dev, 1997, 11: 616-628[本文引用:1]
[10]	Chuck G, Whipple C, Jackson D, Hake S. The maize SBP-box transcription factor encoded by tasselsheath4 regulates bract development and the establishment of meristem boundaries. Development, 2010, 137: 1243-1250[本文引用:2]
[11]	Eveland A L, Goldshmidt A, Pautler M, Morohashi K, Liseron-Monfils C, Lewis M W, Kumari S, Hiraga S, Yang F, Unger-Wallace E, Olson A, Hake S, Vollbrecht E, Grotewold E, Ware D, Jackson D. Regulatory modules controlling maize inflorescence architecture. Genome Res, 2014, 24: 431-443[本文引用:1]
[12]	Lännenpää M, Jänönen I, Hölttä-Vuori M, Gardemeister M, Porali I, Sopanen T. A new SBP-box gene BpSPL1 in silver birch (Betula pendula). Physiol Plant, 2004, 120: 491-500[本文引用:1]
[13]	Chuck G S, Brown P J, Meeley R, Hake S. Maize SBP-box transcription factors unbranched2 and unbranched3 affect yield traits by regulating the rate of lateral primordia initiation. Proc Natl Acad Sci USA, 2014, 111: 18775-18780[本文引用:3]
[14]	Chuck G, Cigan A M, Saeteurn K, Hake S. The heterochronic maize mutant Corngrass1 results from overexpression of a tand em microRNA. Nat Genet, 2007, 39: 544-549[本文引用:2]
[15]	曹雪, 上官凌飞, 于华平, 杨光, 王晨, 谭洪花, 房经贵. 葡萄SBP基因家族生物信息学分析. 基因组学与应用生物学, 2010, 29: 791-798 Cao X, Shang-Guan L F, Yu H P, Yang G, Wang C, Tan H H, Fang J G. Bioinformatics analysis of the SBP gene family in grapevine. Genomics Appl Biol, 2010, 29: 791-798 (in Chinese with English abstract)[本文引用:1]
[16]	刘更森, 慕茜, 戴洪义, 上官凌飞, 张玉刚. 苹果SBP基因家族生物信息学分析. 江西农业学报, 2011, 23(12): 23-27 Liu G S, Mu Q, Dai H Y, Shang-Guan L F, Zhang Y G. Bioinformatics analysis of SBP gene family in apple. Acta Agric Jiangxi, 2011, 23(12): 23-27 (in Chinese with English abstract)[本文引用:1]
[17]	朱命喜, 刘洋, 吴琼, 刘春燕, 徐晶, 陈庆山, 胡国华. 大豆SBP转录因子家族的预测分析. 大豆科学, 2011, 30: 177-183 Zhu M X, Liu Y, Wu Q, Liu C Y, Xu J, Chen Q S, Hu G H. Forecasting analysis of SBP transcription factor families in soybean. Soybean Sci, 2011, 30: 177-183 (in Chinese with English abstract)[本文引用:1]
[18]	葛安静, 张春华, 董清华, 赵密珍, 宋长年, 张希. 草莓SBP基因家族生物信息学初步分析. 中国农学通报, 2012, 28(13): 215-220 Ge A J, Zhang C H, Dong Q H, Zhao M Z, Song C N, Zhang X. Primary bioinformatics analysis of the SBP gene family in strawberry. Chin Agric Sci Bull, 2012, 28(13): 215-220 (in Chinese with English abstract)[本文引用:1]
[19]	万红建, 袁伟, 俞锞, 刘云飞, 李志邈, 叶青静, 王荣青, 阮美颖, 周国治, 姚祝平, 杨悦俭. 番茄SBP基因家族的全基因组鉴定、结构特征及表达分析. 分子植物育种, 2013, 11: 299-306 Wan H J, Yuan W, Yu K, Liu Y F, Li Z M, Ye Q J, Wang R Q, Ruan M Y, Zhou G Z, Yao Z P, Yang Y J. Genome-wide identification, structure characterization and expression analysis of SBP gene family in tomato. Mol Plant Breed, 2013, 11: 299-306 (in Chinese with English abstract)[本文引用:1]
[20]	Wang L, Cao C, Ma Q, Zeng Q, Wang H, Cheng Z, Zhu G, Qi J, Ma H, Nian H, Wang Y. RNA-seq analyses of multiple meristems of soybean: novel and alternative transcripts, evolutionary and functional implications. BMC Plant Biol, 2014, 14: 169[本文引用:1]
[21]	Huang L, Zhao X, Yu Q, Cui W Z, Liu Q X. Evidence for the coevolution of axon guidance molecule Netrin and its receptor Frazzled. Gene, 2014, 544: 25-31[本文引用:1]
[22]	Li P S, Yu T F, He G H, Chen M, Zhou Y B, Chai S C, Xu Z S, Ma Y Z, He G H. Genome-wide analysis of the Hsf family in soybean and functional identification of GmHsf-34 involvement in drought and heat stresses. BMC Genomics, 2014, 15: 1009[本文引用:1]
[23]	Haag J R, Brower-Toland B, Krieger E K, Sidorenko L, Nicora C D, Norbeck A D, Irsigler A, LaRue H, Brzeski J, McGinnis K, Ivashuta S, Pasa-Tolic L, Chand ler V L, Pikaard C S. Functional diversification of maize RNA polymerase IV and V subtypes via alternative catalytic subunits. Cell Rep, 2014, 9: 378-390[本文引用:1]
[24]	Chen X, Zhang Z, Liu D, Zhang K, Li A, Mao L. SQUAMOSA promoter-binding protein-like transcription factors: star players for plant growth and development. J Integr Plant Biol, 2010, 52: 946-951[本文引用:1]
[25]	Martin R C, Asahina M, Liu P P, Kristof J R, Coppersmith J L, Pluskota W E, Bassel G W, Goloviznina N A, Nguyen T T, Pupel P, Nonogaki H. The microRNA156 and microRNA172 gene regulation cascades at post-germinative stages in Arabidopsis. Seed Sci Res, 2010, 20: 79-87[本文引用:1]