玉米bZIP转录因子ZmbZIP26的分子特征、亚细胞定位及响应非生物胁迫的表达分析

曹丽茹, 马晨晨, 庞芸芸, 叶飞宇, 王振华, 鲁晓民

华北农学报. 2023, 38(4): 1-10

华北农学报 ›› 2023, Vol. 38 ›› Issue (4) : 1-10. DOI: 10.7668/hbnxb.20193892
作物遗传育种·种质资源·生物技术

玉米bZIP转录因子ZmbZIP26的分子特征、亚细胞定位及响应非生物胁迫的表达分析

作者信息 +

Molecular Characteristics,Subcellular Localization and Expression Analysis of Maize bZIP Transcription Factor ZmbZIP26 in Response to Abiotic Stress

Author information +
History +

摘要

bZIP转录因子广泛存在于植物中,在调节植物生长发育和非生物胁迫应答中起着重要作用。为探究bZIP转录因子在玉米干旱胁迫响应中的功能,在玉米苗期干旱胁迫处理5 d和复水3 d时,利用转录组测序技术对转录因子基因的表达变化进行分析,筛选到一个响应干旱和复水处理的bZIP转录因子(ZmbZIP26),共表达网络分析发现,ZmbZIP26处于网络调节的核心节点位置。ZmbZIP26基因含有558 bp的开放阅读框,编码185个氨基酸,属于亲水性蛋白。蛋白质进化树及保守序列分析发现,ZmbZIP26蛋白与高粱和芒草同源蛋白的同源性较高,而且在同一氨基酸位置上的保守基序相同。启动子顺式作用元件分析表明,ATG上游2 000 bp区域内含有干旱响应元件、激素响应元件和光响应元件等。qRT-PCR分析发现,ZmbZIP26是组成型表达基因,在幼茎、雌穗和根中高表达,且ZmbZIP26基因积极响应干旱、高温、高盐、氮胁迫及恢复正常的过程,可能在植物的抗逆过程中发挥着重要作用。亚细胞定位分析发现,ZmbZIP26属于核蛋白,定位在细胞核上。蛋白质互作预测发现,ZmbZIP26可能与锌指蛋白、丝氨酸蛋白、钙依赖性蛋白、谷胱甘肽转移蛋白等互作构建了一张调控网络,协同调控玉米生长发育和胁迫应答过程。

Abstract

bZIP transcription factors are widely found in plants and play an important role in regulating plant growth and development and abiotic stress response.In order to explore the functional role of bZIP transcription factor in maize drought stress response,transcriptome sequencing technology was used to analyze the expression changes of transcription factors in maize seedlings treated with drought stress for 5 days and rehydration for 3 days,and a bZIP transcription factor(ZmbZIP26)was screened from transcriptome data in response to drought and rewatering treatment.Co-expression network analysis revealed that ZmbZIP26 was at the core node of network regulation.The gene contained a 558 bp open reading frame encoding 185 amino acids,which was a hydrophilic protein.Phylogenetic tree and conserved sequence analysis showed that ZmbZIP26 protein had high homology with homologous proteins of sorghum and Miscanthus,and also had the same conserved motifs at the same amino acid positions.Cis-element analysis showed that the upstream 2 000 bp region of the ATG site contained drought response elements,hormone response elements and light response elements.qRT-PCR analysis showed that ZmbZIP26 was a constitutively expressed gene,which was highly expressed in young stems,female panicles and roots.ZmbZIP26 positively responded to drought,high temperature,high salt and nitrogen stress and the process of restoring,which might play an important role in the process of plant resistance.Subcellular localization analysis revealed that ZmbZIP26 was a nuclear protein localized in the nucleus.Protein interaction prediction showed that ZmbZIP26 might interact with zinc finger protein,serine protein,Ca-dependent protein and glutathione transfer protein to construct a regulatory network,which cooperatively regulated maize growth and development and stress response process.

关键词

玉米 / bZIP转录因子 / 干旱-复水处理 / 非生物胁迫

Key words

Maize / bZIP transcription factor / Drought-restoring treatment / Abiotic stress

本文二维码

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用
曹丽茹 , 马晨晨 , 庞芸芸 , 叶飞宇 , 王振华 , 鲁晓民. 玉米bZIP转录因子ZmbZIP26的分子特征、亚细胞定位及响应非生物胁迫的表达分析. 华北农学报. 2023, 38(4): 1-10 https://doi.org/10.7668/hbnxb.20193892
Liru CAO , Chenchen MA , Yunyun PANG , Feiyu YE , Zhenhua WANG , Xiaomin LU. Molecular Characteristics,Subcellular Localization and Expression Analysis of Maize bZIP Transcription Factor ZmbZIP26 in Response to Abiotic Stress. Acta Agriculturae Boreali-Sinica. 2023, 38(4): 1-10 https://doi.org/10.7668/hbnxb.20193892
随着全球气候变暖,加剧了极端干旱天气发生的程度和频率,导致植物在生长过程中受到干旱的影响愈加严重,使粮食生产受到严重威胁。玉米(Zea mays L.)是世界上主要的粮食作物之一,在全球的种植范围较为广泛,其产量在粮食作物中居于首位[1]。我国是农业大国,玉米不仅是我国重要的粮食作物和饲料作物,工业生产中也是必不可少的原料[2]。随着我国人口数量的增加以及对能源的需求逐渐增长,在有限的耕地面积下,保证玉米的高产、稳产是保障粮食安全的重要途径[3]。玉米在生长过程中需水量大,对水分胁迫十分敏感,干旱会导致玉米产量降低,甚至严重时导致玉米绝收[4],应对干旱胁迫最有效的途径是培育抗旱玉米新品种,而利用分子育种技术改良创制抗旱新种质是最有效的途径之一,因此,对玉米抗旱基因的挖掘和遗传解析是当前重中之重的工作。
植物体内都存在着自我调节机制和自我恢复能力,当环境发生改变时,作物会出现恢复原生长的趋势。干旱胁迫条件下,作物为适应环境会减慢生长,在干旱后复水又会加快生长,在一定程度上弥补干旱胁迫造成的损失,产生了补偿效应或者超补偿效应[5]。补偿效应在自然界普遍存在,是植物应对逆境胁迫的自我调节机制[6]。Wang等[7]研究发现,籼稻幼苗在干旱胁迫程度未超过其耐受阈值时,对其补水后可以产生补偿或超补偿效应,表现出其恢复力和适应性。刘婷婷等[8]通过对高粱幼苗干旱复水处理,发现干旱时,高粱的相对含水量、叶绿素含量等生理指标均有所下降,复水后各项生理指标恢复正常,高粱干旱适应能力与复水后的恢复性呈正相关。饲料大麦的株高、穗长、单株产量等指标随着干旱胁迫时间的增加显著下降,在复水后这些指标都得到了不同程度的恢复,且恢复程度随着复水量的增加而增加[9]。
植物面临干旱胁迫时体内除了生理生化发生变化,在分子水平和蛋白质水平也会发生一系列的变化,分子水平上的改变起着决定性作用[10]。响应干旱胁迫的基因按照其作用的方式可以分为2种,即功能基因和调控基因,其中调控基因包括蛋白激酶和转录因子[11]。转录因子(Transcription factors,TF) 是一种能和基因启动子区中的顺式作用元件特异性结合的蛋白质分子,调控植物的生长发育过程,并在抵御逆境胁迫过程中发挥着重要作用[12]。bZIP转录因子在转录因子中家族成员最多,且广泛存在于植物中,在调节植物响应非生物胁迫的生理过程中起着十分重要的作用[13]。拟南芥中存在多个与抗逆性相关的AtbZIP 基因。干旱胁迫下,AtbZIP62基因被强烈诱导,该基因过表达后拟南芥的抗旱能力明显增强[14];AtbZIP1过表达使拟南芥的抗盐、抗旱能力都提高[15];在拟南芥中过表达玉米ZmbZIP72基因,发现其激活了Rd29B、Rab18、His1-3 等脱落酸(ABA)诱导基因的高表达,增强了转基因植株的抗旱性[16]。沉默OsbZIP71基因后,发现水稻转基因植株对盐和干旱胁迫十分敏感[17]。 SlbZIP1基因通过调节ABA的合成从而影响番茄对盐和旱胁迫的抗性[18]。
目前,有关bZIP转录因子与植物生长发育关系的研究较多,其在植物抗逆中的作用研究也逐渐增多, 但关于其响应玉米干旱-复水处理的研究较少。本研究利用高通量测序技术对玉米正常、干旱和恢复浇水条件下的叶片和根系进行转录组测序,筛选到一个bZIP转录因子, 暂时命名为ZmbZIP26。对ZmbZIP26进行生物信息学分析,同时运用qRT-PCR技术分析其在玉米不同组织器官中的表达情况,以及在干旱、高温、高盐、氮胁迫处理下的表达特性,为深入研究玉米ZmbZIP26基因的生物学功能奠定基础,也为挖掘抗逆相关分子资源、培育抗逆品种提供依据。

1 材料和方法

1.1 试验材料及处理方法

以玉米优良自交系郑36作为试验材料,当玉米生长到三叶一心时进行干旱胁迫(土壤含水量45%~50%)和复水浇水(土壤含水量大于95%)处理,取胁迫前、干旱胁迫5 d以及复水3 d后的叶片和根系,分别命名为T0dY、T5dY、TR3dY和T0dG、T5dG、TR3dG,5株为一个样品,每个样品3个生物学重复,用于转录组测序和后期差异基因表达量的验证。取在田间正常生长到三叶一心的玉米幼叶、幼茎、幼根,以及抽雄后的成熟叶、成熟茎、成熟根、雄穗、雌穗组织样品进行基因时空表达分析。温室中采用光培养16 h、培养温度为30 ℃,暗培养8 h、培养温度为26 ℃,相对湿度33%~55%。当郑36生长到三叶一心时,取生长一致的玉米苗进行处理,干旱组用20% PEG-6000(Hoagland营养液+20% PEG)处理,盐胁迫组用浓度200 mmol/L的NaCl处理,氮胁迫组将玉米幼苗移至缺乏铵态氮的 Hoagland 营养液中,高温组42 ℃培养,在处理0,4,8,12,24,32,36,48,52,60,72,84 h以及解除胁迫后的12,24 h取玉米根系,3株为一个重复,取3个生物学重复,用于分析基因响应非生物胁迫的表达量变化。取完样品迅速放入液氮中,并存放在-80 ℃的超低温冰箱中。

1.2 转录组测序与差异表达转录因子基因筛选

将样品委托广州基迪奥科技有限公司做高通量测序分析,基因表达量的计算采用FPKM(Fragments per kilobase of exon model per million mapped fragments)算法,根据FPKM值计算基因在处理组和对照组之间的差异倍数,差异表达转录因子基因筛选标准为∣log2(fold change)∣≥1,筛选玉米干旱-复水过程中的差异表达bZIP转录因子,参考Maize GDB(https://www.maizegdb.org/)数据库进行命名。

1.3 ZmbZIP26基因的克隆

对玉米郑36总RNA提取利用RNAprep Pure 多糖多酚植物总 RNA 提取试剂盒(天根生化科技有限公司)的方法,根据反转录试剂盒PrimeScriptTM 1st Strand cDNA Synthesis Kit(TaKaRa)提供的反应程序(25 ℃,5 min;55 ℃,15 min;85 ℃,5 min)合成cDNA。
在MaizeGDB网站上搜索ZmbZIP26(Zm00001d034199)基因的序列,设计特异性引物(ZmbZIP26-F:ATGGACTGGGCAGCGGCGT和ZmbZIP26-R:CTACTGCCAAGAGAATGGC),以cDNA为模板对ZmbZIP26基因的开放阅读框(ORF)进行扩增,将扩增产物连接到pMDTM 18-T载体(TaKaRa)上,转化大肠杆菌DH5α感受态细胞,挑选出阳性单克隆(至少5个样品)送到华大公司进行测序,测序结果与MaizeGDB数据库进行比对。

1.4 ZmbZIP26基因的生物信息学分析

利用ProtParam对ZmbZIP26基因编码蛋白质的分子质量、疏水性及稳定性等进行预测;采用NetPhoS对蛋白质磷酸化位点进行预测;利用SOPMA对ZmbZIP26蛋白的二级结构进行分析;利用软件MEGA 中的PHYLOGENY构建ZmbZIP26蛋白与不同物种同源蛋白的进化树;利用MEME对蛋白质的保守基序(Motif)进行分析,将基序的最大数值设置为6;利用STRING在线工具预测ZmbZIP26的互作蛋白;从数据库 Ensembl Plants(http://plants.ensembl.org/index.html)获取 ZmbZIP26基因上游 2 000 bp序列,用Plantcare(http://bioinformatics.psb.ugent.be/webtools/plantcare/html/)对 ZmbZIP26 基因启动子区的顺式作用元件进行分析。

1.5 实时荧光定量PCR(qRT-PCR)分析

根据ZmbZIP26基因序列设计特异的荧光定量引物(DL-Zm00001d034199-F:CCTTTTCCCCGATTCTCTACTC和DL-Zm00001d034199-R:ACTACGATCTACGACTCTACGA),以玉米18S rDNA作为内参(引物18S-L:CCTGCGGCTTAATTGACTC和18S-R:GTTAGCAGGCTGAGGTCTGG),以上面不同组织、不同胁迫环境的样品RNA反转录cDNA作为模板,参考SYBR Premix Ex TaqTM试剂盒上的程序,采用两步法(95 ℃ 5 min;95 ℃ 10 s,60 ℃ 30 s,40个循环)进行qRT-PCR,采用2-ΔΔCt法对数据进行分析。

1.6 ZmbZIP26蛋白的亚细胞定位分析

亚细胞定位使用表达载体pMDC83-GFP,结合酶切位点SpeⅠ和AscⅠ设计引物(ZmbZIP26-83-GFP-F:CTAGACTAGTATGGACTGGGCAGCGGCGT和ZmbZIP26-83-GFP-R:TTGGCGCGCCCCTACTGCCAAGAGAATGGC)对ZmbZIP26基因进行扩增,经SpeⅠ和AscⅠ双酶切后,用T4-DNA连接酶连接表达载体并转化。将测序正确的重组质粒转到农杆菌EHA105感受态中,利用注射法将农杆菌菌液浸染七叶期生长强壮的烟草叶片,注射48 h后对浸染的叶片通过共聚焦显微镜观察烟草表皮组织细胞内绿色荧光蛋白的表达情况。

2 结果与分析

2.1 ZmbZIP26基因的筛选

对三叶期玉米幼苗进行干旱和复水处理,取根系和叶片部位进行转录组测序,对基因的FPKM值和功能注释进行分析,发现了27个bZIP转录因子成员响应干旱和复水处理,其中ZmbZIP26基因的表达量在正常、干旱和复水处理下达到差异极显著水平(P<0.01)。使用FPKM值计算27个转录因子之间的相关系数(R>0.5),并将这些基因导入软件Cytoscape来构建基因共表达网络,鉴定处于核心节点的关键基因。通过分析,发现ZmbZIP26基因在共表达网络上与其他基因有较高的关联度(图1),说明ZmbZIP26基因是共表达网络的核心关键基因。对玉米叶片和根系中ZmbZIP26 进行qRT-PCR分析,发现在干旱胁迫后,该基因在叶片和根系中相对表达量均显著升高,复水后逐渐恢复正常(图2),且定量结果与转录组分析结果一致,表明ZmbZIP26基因在玉米干旱-复水过程中发挥着重要的作用。
图1 差异表达bZIP转录因子基因的共表达网络

Fig.1 Co-expression network map of differentially expressed transcription factors bZIP-related genes

Full size|PPT slide

图2 ZmbZIP26基因表达的qRT-PCR验证
不同小写字母表示差异显著(P<0.05)。图5—6同。

Fig.2 qRT-PCR verification of ZmbZIP26 gene expression

Different lowercase letters indicate significant difference (P<0.05).The same as Fig.5—6.

Full size|PPT slide

2.2 ZmbZIP26基因的克隆和编码蛋白质性质分析

利用RT-PCR技术对ZmbZIP26基因进行扩增,利用1.2%的琼脂糖凝胶电泳,得到了一条大小约为600 bp的条带(图3),切胶回收后连接克隆载体,筛选出阳性单克隆。测序结果表明,该段序列是一个完整的558 bp的开放阅读框,编码了185个氨基酸,与MaizeGDB数据库中B73的参考序列是一致的。
图3 ZmbZIP26基因PCR扩增结果

Fig.3 PCR amplification of ZmbZIP26 gene

M.2 000 bp DNA Marker;1.ZmbZIP26。

Full size|PPT slide

采用ProtParam软件对蛋白质性质进行分析,结果表明,ZmbZIP26的相对分子质量为20.188 07 ku,分子式为C856H1 439N285O263S8,理论等电点为10.50;蛋白质的不稳定系数为71.82,属于不稳定蛋白;ProtScale分析发现,亲水性值(GRAVY值)为-0.443,属于亲水性蛋白;预测PHD跨膜螺旋区显示,该蛋白质无跨膜结构域;利用软件SOPMA进行预测,发现该蛋白质α螺旋占比73.51%,β转角占比0.54%,无规则卷曲占比20.54%,表明该蛋白质以α 螺旋结构为主(表1)。NetPhoS分析显示,该蛋白质含有14个丝氨酸和5个苏氨酸磷酸化位点,表明ZmbZIP26蛋白可能被丝氨酸和苏氨酸激酶激活,调控下游响应基因,从而参与植物的生长发育、信号传导及逆境胁迫过程等。
表1 ZmbZIP26蛋白性质分析

Tab.1 Property analysis of ZmbZIP26 protein

蛋白质性质
Protein property		 数值
Value
相对分子质量/ku Relative molecular mass 20.188 07
理论等电点 Theory PI value 10.50
不稳定系数 Instability index 71.82
亲水性值 Hydrophilic value -0.443
α螺旋/% α-helix 73.51
β转角/% β-turn 0.54
无规则卷曲/% Random coil 20.54

2.3 ZmbZIP26蛋白的系统进化以及保守序列分析

为了揭示ZmbZIP26的进化关系,利用软件MEGA构建了多物种同源蛋白的进化树,发现ZmbZIP26蛋白与高粱(Sorghum bicolor)和芒草(Miscanthus lutarioriparius)同源蛋白的关系比较紧密,且与水稻(Oryza sativa)同源蛋白的进化关系比较远(图4)。为进一步分析ZmbZIP26蛋白与其他物种同源蛋白的保守性,利用软件MEME鉴定6个保守基序,发现玉米不仅与高粱和芒草同源性较高,而且在同一氨基酸位置上的保守基序相同,与水稻的保守基序相差较大(图4),说明了不同物种间的同源性越高,保守基序相似度越高,其亲缘关系也就越近。
图4 ZmbZIP26蛋白与其他物种同源蛋白进化树及保守基序分析

Fig.4 Phylogenetic tree and conserved motifs analysis of ZmbZIP26 protein and homologous proteins from other species

Full size|PPT slide

2.4 ZmbZIP26基因启动子顺式作用元件分析

为分析ZmbZIP26基因的潜在功能,利用在线软件Plantcare对ZmbZIP26基因ATG上游2 000 bp的顺式作用元件进行了预测,结果显示,该转录因子基因除了启动子本身含有的最基本元件CAAT-box、TATA-box等,还含有ABRE、MBS、TGACG-motif、CGTCA-motif 、TCA-element等多个结合位点(表2)。ABRE参与脱落酸信号传导途径,响应干旱胁迫;MBS是MYB结合位点,参与干旱诱发过程;TGACG-motif、CGTCA-motif参与茉莉酸信号传导途径;TCA-element是水杨酸响应元件,诱导水杨酸途径相关基因的表达。此外,还有大量的响应光刺激的元件Sp1、G-Box、Box 4、I-box等,因此,推测ZmbZIP26基因还可能调控植物开花、参与光合作用等。
表2 ZmbZIP26基因启动子顺式作用元件分析

Tab.2 Cis-element analysis of ZmbZIP26 gene promoter

元件
Element		 序列
Sequence		 位点
Site		 来源
Source		 功能注释
Function annotation
GTGGC-motif GATTCTGTGGC 1 177 菠菜 光响应元件的一部分
CAAT-box CAAAT 1 028 豌豆 启动子和增强子区域常见的顺式作用元件
Sp1 GGGCGG 985 水稻 光敏元件
GATA-motif GATAGGA 315 拟南芥 光响应元件的一部分
CGTCA-motif CGTCA 620 青稞 参与MeJA响应的顺式作用调节元件
G-Box CACGTG 744 豌豆 参与光反应的顺式作用调节元件
Box 4 ATTAAT 444 欧芹 参与光响应的保守DNA模块的一部分
P-box CCTTTTG 867 水稻 赤霉素响应元件
I-box CGATAAGGCG 1 562 玉米 光响应元件的一部分
GC-motif CCCCCG 1 835 玉米 参与缺氧特异性诱导的增强剂样元件
O2-site GATGATGTGG 1 178 玉米 顺式作用调节元件参与玉米醇溶蛋白代谢调节
ABRE ACGTG 618 拟南芥 参与脱落酸反应的顺式作用元件
G-box CACGTG 744 拟南芥 参与光反应的顺式作用调节元件
A-box CCGTCC 1 336 欧芹 顺式调节元件
TATA-box TATA 130 拟南芥 转录起始点-30左右的核心启动子元件
MBS CAACTG 553 拟南芥 MYB结合位点参与干旱诱导
GA-motif ATAGATAA 626 拟南芥 光响应元件的一部分
TGACG-motif TGACG 620 青稞 参与MeJA响应的顺式作用调节元件
TCA-element CCATCTTTTT 1 742 烟草 参与水杨酸反应的顺式作用元件

2.5 ZmbZIP26基因组织表达模式分析

为了解ZmbZIP26基因的表达模式,采用qRT-PCR技术对该基因在玉米组织中表达情况进行了分析,发现其在幼叶、幼茎、幼根、成熟叶、成熟茎、成熟根、雄穗、雌穗中均有表达,其中在幼茎、雌穗、幼根、成熟根中的表达量较高(图5)。ZmbZIP26是一个组成型表达基因,在不同组织、不同时期的表达水平存在差异,表明ZmbZIP26基因在玉米生长发育的过程中起着重要的作用。
图5 ZmbZIP26基因在不同组织中的相对表达量

Fig.5 Relative expression levels of ZmbZIP26 gene in different tissues

Full size|PPT slide

2.6 ZmbZIP26基因在干旱、高温、高盐、氮胁迫下的表达模式分析

通过对ZmbZIP26基因启动子顺式调控元件分析,发现其中含有多个胁迫响应元件,因此对干旱、高温、高盐、氮胁迫处理条件下ZmbZIP26基因在根系中的表达量进行分析,结果显示,干旱、高温、高盐、氮胁迫条件下,ZmbZIP26表达量显著上调,其中在高盐、氮胁迫下表达量上调15~20倍,干旱、高温胁迫下表达量上调5~10倍(图6)。在干旱胁迫24 h,表达量达到峰值,复水后表达量逐步下降;高温胁迫84 h,表达量达到峰值,解除胁迫24 h后表达量又显著降低;高盐胁迫48 h,其表达量达到最高,解除胁迫后表达量显著降低,但解除胁迫24 h后,表达量仍高于处理前,推测高盐胁迫严重伤害了玉米生长;氮处理36 h,其表达量达到最高,解除胁迫后表达量显著降低。以上结果表明,ZmbZIP26基因积极响应干旱、高温、高盐、氮胁迫等非生物胁迫,并且在恢复生长过程中也起着重要作用。
图6 ZmbZIP26基因响应非生物胁迫的表达模式
R.解除胁迫。

Fig.6 Expression pattern of ZmbZIP26 gene in response to abiotic stress

R.Restoring.

Full size|PPT slide

2.7 ZmbZIP26蛋白的亚细胞定位分析

构建载体ZmbZIP26-pMDC83-GFP,利用农杆菌法瞬时转化生长较壮的烟草叶片,通过激光共聚焦显微镜观察发现,pMDC83-GFP处理在细胞膜和细胞核中都检测到了绿色荧光信号,但ZmbZIP26-pMDC83-GFP处理中,仅在细胞核中观察到了绿色荧光信号(图7),表明ZmbZIP26定位在细胞核上,属于核蛋白。
图7 ZmbZIP26蛋白的亚细胞定位

Fig.7 Cellular localization of ZmbZIP26 protein

Full size|PPT slide

2.8 ZmbZIP26的互作蛋白网络预测

为探索ZmbZIP26基因可能的网络作用机制,采用STRING在线预测了ZmbZIP26的互作蛋白,结果显示,与ZmbZIP26互作的蛋白质主要有6个功能结构域(表3)。GRMZM2G159402编码的蛋白质包含锌指蛋白结构,是对基因调控起着十分重要作用的一类转录因子,主要参与细胞分化过程[19];GRMZM2G008109编码的蛋白质包含枯草杆菌蛋白酶sbt3.5结构,是一类丝氨酸蛋白酶,影响细胞壁的形成;GRMZM2G018943编码的蛋白质包含翻译起始因子eIF-2B亚基结构,eIF-2B是与细胞生存紧密相关的真核起始因子,是多亚基复合物,参与调节信号转导途径;GRMZM2G056612编码的蛋白质具有钙依赖性蛋白激酶(Calcium-dependent protein kinase,CDPK)活性,是植物钙信号系统中重要的蛋白激酶;GRMZM2G471517编码的是NPR1(Nonexpressor of pathogenesis-related genes 1)互作子,能激活细胞的免疫反应,保护植物的生长;gst21是谷胱甘肽转移酶,作为谷胱甘肽结合反应的关键酶,催化谷胱甘肽结合反应的起始步骤,在生物体的防御系统中发挥着非常重要的作用。推测ZmbZIP26蛋白与互作蛋白共同构建了一张调控网络,调控着植物的生长发育和逆境胁迫应答过程。
表3 ZmbZIP26的互作蛋白预测

Tab.3 Interaction protein prediction of ZmbZIP26

蛋白质编号
Protein number		 蛋白质描述
Protein description		 氨基酸数目/个
Amino acid number		 互作系数
Interaction coefficient
GRMZM2G024159 无注释 258 0.679
GRMZM2G159402 锌指蛋白8 231 0.679
GRMZM2G404316 无特性蛋白loc100277730 290 0.679
GRMZM2G008109 枯草杆菌样蛋白酶sbt3.5 568 0.678
GRMZM2G018943 翻译起始因子eIF-2B δ亚基 727 0.678
GRMZM2G109805 无特性蛋白loc100501360前体 538 0.678
GRMZM2G056612 CDPK相关激酶5 598 0.674
GRMZM2G471517 NPR1互作子 150 0.616
GRMZM2G409309 无注释 370 0.614
GRMZM2G428168 谷胱甘肽s-转移酶21 235 0.614

3 结论与讨论

目前,在分析非生物胁迫下植物的代谢通路中,转录组测序技术已经被广泛应用[20],前人对茶树、棉花等植物响应非生物胁迫的调控机制进行了转录组测序分析,从分子层面解析了非生物胁迫下功能基因和调节基因的调控网络[21-22]。转录因子在植物非生物胁迫应答生理过程中起着重要的调节作用,近年来,在不同物种中bZIP基因家族的功能得到了广泛的研究,尤其是在拟南芥和水稻中[23-24], 在玉米中也有170多个bZIP蛋白被鉴定[25],但迄今为止,只有少数玉米bZIP基因的功能被报道[26]。本研究通过对玉米幼苗在干旱-复水过程的根系和叶片进行转录组测序,筛选到一个响应干旱和复水处理的bZIP转录因子(ZmbZIP26),基因共表达网络分析发现,ZmbZIP26处于网络调节的核心节点位置。进化树和保守基序分析显示,玉米中该蛋白质与高粱和芒草的同源性较高,且在相同氨基酸位置具有相同的保守基序;对ZmbZIP26基因ATG上游2 000 bp启动子序列进行顺式调控元件分析,结果显示,启动子序列中包含着与植物的生长发育、光反应、逆境胁迫和信号转导等相关的结合位点,且参与玉米脱落酸、茉莉酸、水杨酸等多种激素的信号传导途径。
本研究发现,ZmbZIP26在玉米的多个组织部位均有表达,表明ZmbZIP26是组成型表达基因,该结果与ZmbZIP9、 ZmbZIP17、ZmbZIP19等大部分ZmbZIP的研究结果是一致的[27],但ZmbZIP26在幼茎、雌穗和根中表达量较高,与其他ZmbZIP基因的表达模式存在着明显的差异,可见ZmbZIP基因的表达模式存在保守性和多样性。ZmbZIP26基因在干旱、高温、高盐、氮胁迫处理条件下根系中的表达量显著上调5~20倍,这与拟南芥和苹果中的表达趋势是一致的[28],解除胁迫后,ZmbZIP26基因表达量显著下降,说明ZmbZIP26基因有较强的抗逆性和恢复补偿能力。亚细胞定位分析发现,ZmbZIP26定位在细胞核上,这与大多数的ZmbZIP定位结果是一致的。转录因子进入细胞核大多数借助核定位信号,使靶基因表达,杨艳歌[29]证实玉米中ZmbZIP确实发挥着转录因子功能,但也有报道显示ZmbZIP60定位在内质网上[30],拟南芥中AtbZIP60也定位在内质网上,在内质网应激信号转导中发挥作用[31],说明同一家族的转录因子调控作用机制也会存在不同。
为进一步研究ZmbZIP26基因的互作网络机制,本研究对ZmbZIP26的互作蛋白进行预测,发现互作蛋白中包含锌指蛋白结构、枯草杆菌蛋白酶sbt3.5结构、翻译起始因子eIF-2B δ亚基、钙依赖性蛋白激酶、NPR1互作子、谷胱甘肽s-转移酶21等。锌指蛋白在植物生长发育和逆境应答中发挥重要作用,在烟草中CsZFP2基因能够积极响应干旱胁迫,过表达CsZFP2的转基因植株比野生型植株更具抗旱性[32];OsC3H10的过表达通过调控胁迫相关基因提高了水稻的抗旱性[33]。枯草杆菌蛋白酶sbt3.5结构在植物果胶甲基酯酶PME17的处理过程中起作用,在根系发育过程中高度共表达[34]。钙依赖性蛋白激酶是一类主要的钙信号受体,在植物的生长发育和非生物应激反应中发挥着重要作用[35]。低温胁迫下CDPK20基因在山药中具有调控功能,提高山药的耐低温能力[36];水稻基因OsCDPK7/13过表达后,转基因植株提高了对盐胁迫和干旱胁迫的耐受性[37]。NPR1在植物抗病信号转导过程中发挥着重要作用,甘蔗ShNPR1基因在烟草对青枯菌(Ralstonia solanacearum)和茄病镰刀菌蓝色变种(Fusarium solani var.coeruleum)的防御反应中起积极作用[38]。因此推测,ZmbZIP26可能与锌指蛋白、丝氨酸蛋白、钙依赖性蛋白、谷胱甘肽转移蛋白等互作构建了一张调控网络,协同调节玉米生长发育和胁迫应答过程。
综上,从转录组数据中筛选到一个响应干旱和复水处理的bZIP转录因子ZmbZIP26,基因共表达网络分析发现,ZmbZIP26处于网络调节的核心节点位置;qRT-PCR分析发现,该基因属于组成型表达基因,在幼叶、幼茎、幼根、成熟叶、成熟茎、成熟根、雄穗、雌穗中均有表达,且在幼茎、雌穗和根中的表达量较高;通过非生物胁迫处理发现,ZmbZIP26基因积极响应干旱、高温、高盐、氮胁迫和恢复正常过程,可能在植物的抗逆过程中发挥着重要作用;亚细胞定位分析发现,ZmbZIP26属于核蛋白,定位在细胞核上;蛋白质互作预测发现,ZmbZIP26可能与锌指蛋白、丝氨酸蛋白、钙依赖性蛋白、谷胱甘肽转移蛋白等互作构建了一张调控网络,协同调控玉米的生长发育过程和胁迫应答过程。研究结果为揭示ZmbZIP26基因在玉米生长发育过程和非生物胁迫应答反应中的作用机制提供了参考,也为创制抗逆种质提供了分子资源。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序
[1]
陆伟婷, 于欢, 曹胜男, 陈长青. 近20年黄淮海地区气候变暖对夏玉米生育进程及产量的影响[J]. 中国农业科学, 2015, 48(16):3132-3145.doi:10.3864/j.issn.0578-1752.2015.16.004.
Lu W T, Yu H, Cao S N, Chen C Q. Effects of climate warming on growth process and yield of summer maize in Huang-Huai-Hai plain in last 20 years[J]. Scientia Agricultura Sinica, 2015, 48(16):3132-3145.
本文引用 [1]
[2]
肖祖栋, 陈先敏, 李斌彬, 申思, 邓涛, 李凤元, 周顺利. 不同品种夏玉米生长发育和产量形成对播期和密度的响应特征[J]. 华北农学报, 2023, 38(1):110-116.doi:10.7668/hbnxb.20193312.
Xiao Z D, Chen X M, Li B B, Shen S, Deng T, Li F Y, Zhou S L. Response characteristics of growth-development and yield-formation of different summer maize cultivars to sowing date and planting density[J]. Acta Agriculturae Boreali-Sinica, 2023, 38(1):110-116.
本文引用 [1]
[3]
张镇涛, 杨晓光, 高继卿, 王晓煜, 白帆, 孙爽, 刘志娟, 明博, 谢瑞芝, 王克如, 李少昆. 气候变化背景下华北平原夏玉米适宜播期分析[J]. 中国农业科学, 2018, 51(17):3258-3274.doi:10.3864/j.issn.0578-1752.2018.17.003.
Zhang Z T, Yang X G, Gao J Q, Wang X Y, Bai F, Sun S, Liu Z J, Ming B, Xie R Z, Wang K R, Li S K. Analysis of suitable sowing date for summer maize in North China plain under climate change[J]. Scientia Agricultura Sinica, 2018, 51(17):3258-3274.
本文引用 [1]
[4]
张晓旭, 孙忠富, 郑飞翔, 刘江, 李崇瑞, 王一昊. 基于作物水分亏缺指数的黄淮海平原夏玉米全生育期干旱分布特征[J]. 中国农业气象, 2021, 42(6):495-506.doi:10.3969/j.issn.1000-6362.2021.06.005.
Zhang X X, Sun Z F, Zheng F X, Liu J, Li C R, Wang Y H. Characteristics of drought distribution for summer maize over whole growth period in Huang-Huai-Hai plain based on crop water deficit index[J]. Chinese Journal of Agrometeorology, 2021, 42(6):495-506.
本文引用 [1]
[5]
曹丹, 陈道钳, 吴茜, 殷俐娜, 邓西平, 王仕稳. 复水对旱后不同玉米品种植株生长恢复能力及其生理响应特性的影响[J]. 西北植物学报, 2015, 35(6):1222-1228.doi:10.7606/j.issn.1000-4025.2015.06.1222.
Cao D, Chen D Q, Wu Q, Yin L N, Deng X P, Wang S W. Research on rewatering post-drought growth recovery capacity and physiological characteristics of different maize varieties[J]. Acta Botanica Boreali-Occidentalia Sinica, 2015, 35(6):1222-1228.
本文引用 [1]
[6]
厉广辉. 花生抗旱性状鉴定及不同品种抗旱的生理机制研究[D]. 泰安: 山东农业大学, 2014.doi:10.7666/d.Y2587813.
Li G H. Identification of drought resistance traits of peanut and study on physiological mechanism of drought resistance of different varieties[D]. Taian: Shandong Agricultural University, 2014.
本文引用 [1]
[7]
Wang Y S, Ni F, Yin D H, Chen L J, Li Y H, He L X, Zhang Y L. Physiological response of Lagerstroemia indica (L.)Pers.seedlings to drought and rewatering[J]. Tropical Plant Biology, 2021, 14(4):360-370.doi:10.1007/s12042-021-09294-3.
本文引用 [1]
[8]
刘婷婷, 陈道钳, 王仕稳, 殷俐娜, 邓西平. 不同品种高粱幼苗在干旱复水过程中的生理生态响应[J]. 草业学报, 2018, 27(6):100-110.doi:10.11686/cyxb2017410.
Liu T T, Chen D Q, Wang S W, Yin L N, Deng X P. Physio-ecological responses to drought and subsequent re-watering in sorghum seedlings[J]. Acta Prataculturae Sinica, 2018, 27(6):100-110.
本文引用 [1]
[9]
徐银萍, 潘永东, 刘强德, 任诚, 姚元虎, 贾延春, 陈文庆, 火克仓, 包奇军, 赵锋, 张华瑜. 生长后期干旱复水对饲草大麦产量、品质及叶绿素含量的影响[J]. 中国土壤与肥料, 2020(2):192-197.doi:10.11838/sfsc.1673-6257.19175.
Xu Y P, Pan Y D, Liu Q D, Ren C, Yao Y H, Jia Y C, Chen W Q, Huo K C, Bao Q J, Zhao F, Zhang H Y. Effects of rewatering after drought stress on yield quality and chlorophyll of forage barley in late growth period[J]. Soils and Fertilizers Sciences in China, 2020(2):192-197.
本文引用 [1]
[10]
Chinnusamy V, Schumaker K, Zhu J K. Molecular genetic perspectives on cross-talk and specificity in abiotic stress signalling in plants[J]. Journal of Experimental Botany, 2004, 55(395):225-236.doi:10.1093/jxb/erh005.
本文引用 [1]
[11]
Avramova V, AbdElgawad H, Zhang Z F, Fotschki B, Casadevall R, Vergauwen L, Knapen D, Taleisnik E, Guisez Y, Asard H, Beemster G T S. Drought induces distinct growth response,protection,and recovery mechanisms in the maize leaf growth zone[J]. Plant Physiology, 2015, 169(2):1382-1396.doi:10.1104/pp.15.00276.
本文引用 [1]
[12]
Zhao Y, Cheng X Y, Liu X D, Wu H F, Bi H H, Xu H X. The wheat MYB transcription factor TaMYB31 is involved in drought stress responses in Arabidopsis[J]. Frontiers in Plant Science, 2018, 9:1426.doi:10.3389/fpls.2018.01426.
本文引用 [1]
[13]
Baloglu M C, Eldem V, Hajyzadeh M, Unver T. Genome-wide analysis of the bZIP transcription factors in cucumber[J]. PLoS One, 2014, 9(4):e96014.doi:10.1371/journal.pone.0096014.
本文引用 [1]
[14]
Rolly N K, Imran Q M, Shahid M, Imran M, Khan M, Lee S U, Hussain A, Lee I J, Yun B W. Drought-induced AtbZIP62 transcription factor regulates drought stress response in Arabidopsis[J]. Plant Physiology and Biochemistry, 2020, 156:384-395.doi:10.1016/j.plaphy.2020.09.013.
本文引用 [1]
[15]
Sun X L, Li Y, Cai H, Bai X, Ji W, Ding X D, Zhu Y M. The Arabidopsis AtbZIP1 transcription factor is a positive regulator of plant tolerance to salt,osmotic and drought stresses[J]. Journal of Plant Research, 2012, 125(3):429-438.doi:10.1007/s10265-011-0448-4.
本文引用 [1]
[16]
Ying S, Zhang D F, Fu J, Shi Y S, Song Y C, Wang T Y, Li Y. Cloning and characterization of a maize bZIP transcription factor,ZmbZIP72,confers drought and salt tolerance in transgenic Arabidopsis[J]. Planta, 2012, 235(2):253-266.doi:10.1007/s00425-011-1496-7.
本文引用 [1]
[17]
Liu C T, Mao B G, Ou S J, Wang W, Liu L C, Wu Y B, Chu C C, Wang X P. OsbZIP71,a bZIP transcription factor,confers salinity and drought tolerance in rice[J]. Plant Molecular Biology, 2014, 84(1/2):19-36.doi:10.1007/s11103-013-0115-3.
本文引用 [1]
[18]
Zhu M K, Meng X Q, Cai J, Li G, Dong T T, Li Z Y. Basic leucine zipper transcription factor SlbZIP1 mediates salt and drought stress tolerance in tomato[J]. BMC Plant Biol, 2018, 18(1):83.doi:10.1186/s12870-018-1299-0.
本文引用 [1]
[19]
赵楠, 赵飞, 李玉花. 锌指蛋白结构及功能研究进展[J]. 生物技术通讯, 2009, 20(1):131-134.doi:10.3969/j.issn.1009-0002.2009.01.037.
Zhao N, Zhao F, Li Y H. Advances in research on zinc finger protein[J]. Letters in Biotechnology, 2009, 20(1):131-134.
本文引用 [1]
[20]
Grassmann F. Conduct and quality control of differential gene expression analysis using high-throughput transcriptome sequencing(RNASeq)[M]// Methods in Molecular Biology. New York: Springer New York, 2019, 1834:29-43.doi:10.1007/978-1-4939-8669-9_2.
本文引用 [1]
[21]
Xie T, Chen C J, Li C H, Liu J R, Liu C Y, He Y H. Genome-wide investigation of WRKY gene family in pineapple:Evolution and expression profiles during development and stress[J]. BMC Genomics, 2018, 19(1):490.doi:10.1186/s12864-018-4880-x.
本文引用 [1]
[22]
Bowman M J, Park W, Bauer P J, Udall J A, Page J T, Raney J, Scheffler B E, Jones D C, Campbell B T. RNA-Seq transcriptome profiling of upland cotton(Gossypium hirsutum L.)root tissue under water-deficit stress[J]. PLoS One, 2013, 8(12):e82634.doi:10.1371/journal.pone.0082634.
本文引用 [1]
[23]
Yang O, Popova O V, Süthoff U, Lüking I, Dietz K J, Golldack D. The Arabidopsis basic leucine zipper transcription factor AtbZIP24 regulates complex transcriptional networks involved in abiotic stress resistance[J]. Gene, 2009, 436(1/2):45-55.doi:10.1016/j.gene.2009.02.010.
本文引用 [1]
[24]
Lakra N, Nutan K K, Das P, Anwar K, Singla-Pareek S L, Pareek A. A nuclear-localized histone-gene binding protein from rice(OsHBP1b)functions in salinity and drought stress tolerance by maintaining chlorophyll content and improving the antioxidant machinery[J]. Journal of Plant Physiology, 2015, 176:36-46.doi:10.1016/j.jplph.2014.11.005.
本文引用 [1]
[25]
Wei K F, Chen J, Wang Y M, Chen Y H, Chen S X, Lin Y N, Pan S, Zhong X J, Xie D X. Genome-wide analysis of bZIP-encoding genes in maize[J]. DNA Research, 2012, 19(6):463-476.doi:10.1093/dnares/dss026.
本文引用 [1]
[26]
Ma H Z, Liu C, Li Z X, Ran Q J, Xie G N, Wang B M, Fang S, Chu J F, Zhang J R. ZmbZIP4 contributes to stress resistance in maize by regulating ABA synthesis and root development[J]. Plant Physiology, 2018, 178(2):753-770.doi:10.1104/pp.18.00436.
本文引用 [1]
[27]
贾利强, 赵秋芳, 陈曙. 11个玉米bZIP基因应答逆境胁迫的表达分析[J]. 安徽农业大学学报, 2022, 49(2):197-204.doi:10.13610/j.cnki.1672-352x.20220520.002.
Jia L Q, Zhao Q F, Chen S. The expression profiling of 11 ZmbZIP genes under abiotic stress in maize[J]. Journal of Anhui Agricultural University, 2022, 49(2):197-204.
本文引用 [1]
[28]
Zhao J, Guo R R, Guo C L, Hou H M, Wang X P, Gao H. Evolutionary and expression analyses of the apple basic leucine zipper transcription factor family[J]. Frontiers in Plant Science, 2016, 7:376.doi:10.3389/fpls.2016.00376.
本文引用 [1]
[29]
杨艳歌. 玉米ZmbZIPs基因家族的克隆及抗逆基因的筛选与功能研究[D]. 大庆: 黑龙江八一农垦大学, 2013.
Yang Y G. Cloning of Zea mays ZmbZIPs gene family in maize,screening and functional study of stress resistance genes[D]. Daqing: Heilongjiang Bayi Agricultural University, 2013.
本文引用 [1]
[30]
Wang B, Zheng J, Liu Y J, Wang J H, Wang G Y. Cloning and characterization of the stress-induced bZIP gene Zmb ZIP60 from maize[J]. Molecular Biology Reports, 2012, 39(5):6319-6327.doi:10.1007/s11033-012-1453-y.
本文引用 [1]
[31]
Iwata Y, Koizumi N. An Arabidopsis transcription factor,AtbZIP60,regulates the endoplasmic reticulum stress response in a manner unique to plants[J]. Proceedings of the National Academy of Sciences of the United States of America, 2005, 102(14):5280-5285.doi:10.1073/pnas.0408941102.
本文引用 [1]
[32]
李婉雪, 姚新转, 张宝会, 刘洋. CsZFP2基因的克隆及其在烟草中的抗旱功能分析[J]. 福建农业学报, 2022, 37(4):453-459.doi:10.19303/j.issn.1008-0384.2022.004.005.
Li W X, Yao X Z, Zhang B H, Liu Y. Connection between tea ZFP2 and drought resistance of tobacco plants[J]. Fujian Journal of Agricultural Sciences, 2022, 37(4):453-459.
本文引用 [1]
[33]
Seong S Y, Shim J S, Bang S W, Kim J K. Overexpression of OsC3H10,a CCCH-zinc finger,improves drought tolerance in rice by regulating stress-related genes[J]. Plants (Basel,Switzerland), 2020, 9(10):E1298.doi:10.3390/plants9101298.
本文引用 [1]
[34]
Sénéchal F, Graff L, Surcouf O, Marcelo P, Rayon C, Bouton S, Mareck A, Mouille G, Stintzi A, Höfte H, Lerouge P, Schaller A, Pelloux J. Arabidopsis PECTIN METHYLESTERASE17 is co-expressed with and processed by SBT3.5,a subtilisin-like serine protease[J]. Annals of Botany, 2014, 114(6):1161-1175.doi:10.1093/aob/mcu035.
本文引用 [1]
[35]
麻浩, 王爽, 周亚丽. 植物中钙依赖蛋白激酶的研究进展[J]. 南京农业大学学报, 2017, 40(4):565-572. doi:10.7685/jnau.201701036.
Ma H, Wang S, Zhou Y L. Research progress of calcium-dependent protein kinases in plants[J]. Journal of Nanjing Agricultural University, 2017, 40(4):565-572.
本文引用 [1]
[36]
王莹, 王晓宇, 孙德慧, 霍红雁, 刘海臣, 徐惠, 张继星. 蓖麻钙依赖蛋白激酶29基因克隆与表达分析[J]. 华北农学报, 2023, 38(2):85-92. doi:10.7668/hbnxb.20192958.
Wang Y, Wang X Y, Sun D H, Huo H Y, Liu H C, Xu H, Zhang J X. Cloning and expression analysis of castor calcium dependent protein kinase 29 gene[J]. Acta Agriculturae Boreali-Sinica, 2023, 38(2):85-92.
本文引用 [1]
[37]
周子馨, 兰海燕. 植物CDPK信号转导途径及其互作组分的功能研究进展[J]. 新疆大学学报(自然科学版(中英文)), 2021(6):705-714.doi:10.13568/j.cnki.651094.651316.2020.12.16.0002.
Zhou Z X, Lan H Y. CDPK signaling pathway and its interaction components in plants[J]. Journal of Xinjiang University (Natural Science Edition in Chinese and English), 2021(6):705-714.
本文引用 [1]
[38]
Zang S J, Qin L Q, Zhao Z N, Zhang J, Zou W H, Wang D J, Feng A Y, Yang S L, Que Y X, Su Y C. Characterization and functional implications of the nonexpressor of pathogenesis-related genes 1(NPR1)in Saccharum[J]. Int J Mol Sci, 2022, 23(14):7984.doi:10.3390/ijms23147984.
本文引用 [1]

基金

中央引导地方科技发展资金(Z20221343040)
河南省农业科学院杰出青年科技基金(2022JQ02)
河南省青年人才托举工程项目(2022HYTP034)
国家自然科学基金青年基金(32201708)

文章所在专题

玉米

108

Accesses

0

Citation

Detail
段落导航
相关文章

/

〈 〉