玉米中低温应答ERF基因鉴定及表达分析

ERF转录因子在植物的多种发育过程和非生物胁迫反应中发挥重要作用。在玉米中,虽然ERF家族基因对盐胁迫和干旱胁迫的响应研究已经有报道,但尚没有研究从全基因组水平鉴定ERF家族基因对低温胁迫的响应。本研究采用高通量测序技术Illumina HiSeq,分别在5℃和22℃处理下构建玉米叶组织转录组,从中筛选出57条低温响应ZmERFs,包括53条上调表达的Zm ERFs,4条下调表达的Zm ERFs。不同的Zm ERFs在不同组织部位和发育阶段的表达模式有很大变化。蛋白保守结构域分析表明ZmERFs蛋白中共有20个保守基序,EAR基序出现在组B1成员中,LWSY基序出现在A1,A5和B2成员中。启动子分析显示,多种逆境相关元件ABRE、ARE、LTR和MBS在低温响应ZmERFs的启动子中高度富集。RT-qPCR结果进一步显示GRMZM2G544539和GRMZM2G040664基因可以同时受低温、盐和干旱胁迫诱导表达,GRMZM2G434203基因可以同时受低温和干旱胁迫诱导表达,表明ERF转录因子可以参与包括低温在内的多重非生物胁迫反应。此外,启动子中含有7个ABRE元件的GRMZM2...     

利用WGCNA鉴定非生物胁迫相关基因共表达网络

加权共表达网络分析(Weighted Gene Co-expression Network Analysis, WGCNA)是用来描述不同样品之间基因关联模式的系统生物学方法,可以用来鉴定高度协同变化的基因集。本研究利用正常水稻组织共47份转录组数据,通过冷胁迫、干旱胁迫、盐胁迫不同的处理方式,使用WGCNA方法,根据已克隆基因的报道与以上3种胁迫相关的关键基因,探究不同逆境下基因之间的调控关系。通过对低表达量基因的过滤,最终利用筛选的30,339个表达的基因来构建共表达矩阵,得到15个模块。分析发现已知的水稻3种相关基因在各个模块均有存在,于是对预测到的靶基因进行GO富集分析。对3种胁迫下处理的转录组数据进行差异表达基因分析,结合已报道与胁迫相关的基因,选取各胁迫相关的2个模块进行了基因调控网络的构建。鉴于3种胁迫相关基因在green模块中大量分布,通过对green模块下各自特有的基因和共有的基因的GO功能富集分析,并对共有的基因构建调控网络,挖掘到2599个与3种胁迫都相关的基因,并预测出25个抗逆相关的关键基因,为水稻的抗逆及综合抗逆能力等研究提供了新思路。     

基于转录组及WGCNA的甘蔗干旱响应调控网络分析

干旱是制约甘蔗产业发展的重要因素之一。前人研究表明,斑茅具有良好的抗性基因,可以通过远缘杂交遗传给后代。本研究以斑茅与甘蔗杂交F1代无性系YCE96-40为材料,对苗期干旱处理0 h和24 h后的叶和根进行转录组测序分析,比较了根和叶在转录水平上响应干旱的差异,鉴定出21,885个(DR vs CR:10176, DL vs CL:7907)差异表达基因(DEGs),根中差异表达基因多于叶中,说明根对干旱胁迫响应更为剧烈。GO功能富集分析发现,根和叶中DEGs均富集到与脱水反应相关及激素信号转导过程相关的条目,比如“对渗透胁迫生物过程的响应”和“对缺水生物过程的反应”等。与叶不同的是,根中大量DEGs显著富集到与细胞膜相关的条目。在根中鉴定出多个木质素相关DEGs,表明木质素参与了根的干旱响应。通过对所有DEGs进行WGCNA分析,共鉴定出11个基因共表达模块,其中与干旱处理后的根显著相关有5个模块,与干旱处理后的叶显著相关的有2个模块。进一步筛选出了26个转录因子为甘蔗干旱响应的候选转录因子,构建了调控网络。研究结果为进一步理解甘蔗抗旱性的分子机制及甘蔗抗旱性育种提供了理论指导。     

利用WGCNA进行玉米花期基因共表达模块鉴定

权重基因共表达网络分析(weighted gene co-expression network analysis, WGCNA)是系统生物学的一种研究方法,在挖掘生物学数据与特定性状之间的生物学关系方面具有十分重要的作用。本研究利用玉米(Zea mays L.)自交系B73的14份不同发育阶段的转录组数据,筛选掉低表达丰度的基因,最终得到了22,426个高表达的基因用于创建基因表达矩阵;利用不同组织作为性状,创建表型矩阵。然后利用R软件中的WGCNA包建立了共表达网络,共得到20个模块。本研究将与组织相关性高于0.65的模块定义为组织特异性模块,最终鉴定到14个组织特异性模块。利用在线网站Agrigo对组织特异性模块中的基因进行GO (gene ontology)富集分析,发现14个模块中均可以得到富集种类。开花作为玉米生育周期中的一个重要生理过程,不仅代表着植物从营养生长到生殖生长的转变,也关系到产量、株高和抗逆性等农艺性状。本研究发现8个组织特异性模块中的基因可以富集到与开花调控的代谢通路。此外,有17个已经报道过的开花时间调控基因存在于共表达模块中,并且主要分布在Blue模块和Darkmagenta模块,因此本研究重点关注了这2个模块内部的基因调控网络。本研究通过计算不同组织中的基因表达丰度,并联合权重基因共表达网络分析的方法,鉴定到了具有生物学意义的共表达基因模块,挖掘到了数个开花相关的模块,有助于揭示玉米开花调控的遗传机制。     

甜玉米乳熟期果皮厚度相关基因的加权基因共表达网络分析

为深入探究甜玉米果皮厚度发育的分子机制,利用华南农业大学甜玉米课题组选育的果皮厚度差异较大的甜玉米自交系M03和M08为试验材料,对其授粉后15,19,23 d的果皮进行转录组测序。联合差异表达基因分析与加权基因共表达网络(WGCNA)分析,鉴定甜玉米乳熟期果皮厚度相关的共表达基因模块,根据模块内基因的连接度鉴定核心基因,利用qRT-PCR验证核心基因表达量的可靠性。比较M03和M08不同时期果皮转录组数据共筛选出4 748个差异表达基因(DEGs),DEGs主要富集到碳水化合物代谢过程、植物-病原体相互作用、植物MAPK信号通路。WGCNA分析共构建了18个共表达模块,通过筛选得到4个具有生物学意义且与果皮厚度显著相关(相关系数的绝对值>0.60)的特异性共表达模块,分别是Turquoise、Yellow、Magenta和Pink模块。GO和KEGG功能富集分析结果均表明,特异性模块富集到的天冬氨酸、丙氨酸和谷氨酸代谢、半胱氨酸和蛋氨酸代谢以及植物MAPK信号通路在甜玉米果皮厚度发育中具有重要的作用。选取连接度最高的前20个基因进行核心基因筛选,通过基因功能注释最终筛选到包括M...     

利用WGCNA鉴定玉米株高和穗位高基因共表达模块

株高和穗位高是玉米株型的重要影响因子,与产量性状紧密相关。加权基因共表达网络分析(weightedgene co-expression network analysis, WGCNA)是探索基因网络与特定性状间关联关系的重要方法,为株高和穗位高相关基因的挖掘提供新途径。本研究利用郑58、掖478、昌7-2和黄早四及其组配的杂交种郑单958、安玉5号、郑58/黄早四和掖478/黄早四,结合其在45,000株hm^–2和67,500株hm^–2条件下的转录组数据,采用WGCNA构建了2种密度条件下的共表达网络,分别得到24个和21个共表达模块,并鉴定到与株高和穗位高显著且高度相关(相关系数的绝对值>0.50)的共表达模块15个,其中两性状相同的模块共6个。基因功能富集分析结果表明,株高和穗位高共表达模块主要参与生长发育、光合作用、响应光刺激、植物激素、碳水化合物合成/代谢等重要活动。根据模块内基因的连接度,发现乙烯响应因子EREB14、硫胺素酶TENA2、磷酸甘油酸激酶PGK、谷胱甘肽转移酶GST2和琥珀酸脱氢酶SUDH7等是模块内的核心基因。...     

基于WGCNA的马铃薯根系抗旱相关共表达模块鉴定和核心基因发掘

权重基因共表达网络分析(Weighted Gene Co-expression Network Analysis, WGCNA)是系统生物学皀一种研究方法,在多样本转录组数据中挖掘与目标性状相兲皀基因模块有较广泛皀应用。为了深入探究马铃薯应对干旱胁迫皀分子机制,本研究以国际马铃薯中心引迚栻培种C16 (CIP 397077.16)和C119 (CIP 398098.119)为试验材料,将其无菌组培苗利用甘露醇模拟干旱胁迫,处理0 h、2 h、6 h、12 h和24 h,取其根系迚行转录组测序,每样设3个生物学重复,共30个样本。基于以上转录组数据,利用WGCNA构建与抗逆生理性状相兲联皀权重基因共表达网络,得到15个与根系抗旱密切相兲皀基因共表达模块,幵从4个与目标性状兲联度最高皀模块中収掘到数个核心基因,功能注释表明其中大部分参与干旱胁迫调控通路。这些结果为迚一步研究马铃薯根系抗旱皀分子遗传机制提供了线索。     

玉米自交系响应花粒期高温胁迫差异表达基因的分析

为了探明玉米重要自交系花粒期高温胁迫下转录组基因的表达差异,发掘玉米响应高温胁迫的关键基因和蛋白质,明确高温胁迫引起玉米减产的机理,从而利用分子生物学技术手段提高玉米耐高温胁迫性。首先,利用Ampha~(TM)Z30花粉活性分析仪检测高温胁迫和正常生长条件下玉米自交系昌7-2、郑58的花粉活性。然后收集花粉提取总RNA,利用Illumina Hiseq2000~(TM)高通量测序技术分别对昌7-2、郑58花粒期高温胁迫材料和正常生长条件下的对照材料进行全转录组测序,序列结果分析后获得差异显著表达基因。通过差异基因维恩图(VENN图)重叠区域分析、GO功能分类和富集统计、KEGG pathway通路分类和富集分析以及转录因子分析,获得玉米花粒期高温胁迫相关的重要差异表达基因和蛋白质。花粉活性检测结果显示,高温胁迫后和正常生长条件下,昌7-2花粉活性分别为24. 37%,42. 89%,郑58花粉活性分别为35. 57%,64. 83%。高温胁迫后在昌7-2花粉转录组中共检测到4 176个显著差异表达基因,郑58花粉转录组中共检测到5 487个显著差异表达基因。KEGG pathway通路分类和富集分析结果显示,高温胁迫后郑58的花粉转录组中有2 062个差异表达基因富集到399个相关通路上,昌7-2的花粉转录组中有1 943个差异表达基因富集到352个相关通路上。转录因子分析结果表明,共获得55个转录因子家族,其中,有45个转录因子包含10个以上差异表达基因。研究表明,通过对玉米自交系花粒期高温胁迫后花粉转录组测序获得了大量的差异显著表达基因信息,昌7-2、郑58对高温胁迫响应机制存在明显差异,热激蛋白(Hsps)、热激转录因子(Hsfs)、生长素(IAA)基因和磷脂酰基醇-4,5-二磷酸基因家族(PIP2)在响应玉米花粒期高温胁迫过程中起着重要作用。     

玉米ZmWRKY53基因克隆及诱导表达分析

WRKY是植物细胞中进化较为保守的一类转录因子,参与了生物和非生物逆境胁迫应答反应的调控。本研究基于干旱诱导的玉米(Zea mays L.)B73转录组数据库,成功克隆到1个WRKY基因,命名为Zm WRKY53。序列比对发现,该基因含有904 bp的完整开放读码框,编码302个氨基酸,在C端含有1个保守的WRKYGQK结构域,同源比对ZmWRKY53属于IIc类WRKY转录因子。生物信息学预测,Zm WRKY53为碱性不稳定亲水性蛋白,无信号肽和跨膜结构域。洋葱内表皮亚细胞定位显示,Zm WRKY53蛋白位于细胞核中。实时荧光定量PCR(qRT-PCR)分析表明,Zm WRKY53基因在玉米的根和叶中有较高的表达水平,而在茎、玉米须、雄穗和胚中的表达量较低。不同的非生物胁迫如低温(8℃)、NaCl、H2O2和PEG均能显著提高Zm WRKY53的表达水平,尤其是PEG诱导使其表达水平上调9.5倍。接种玉米纹枯病的叶片中,Zm WRKY53的表达量提高了13倍。对比不同的激素处理情况,脱落酸(abscisic acid,ABA)和水杨酸(salicylic acid,SA)处理后,Zm WRKY53的表达水平持续上调;而在赤霉素(gibberellin,GA)和乙烯处理的玉米根中,Zm WRKY53呈现无规律的变化特征。蛋白互作预测,十个候选基因可作为后续Zm WRKY53调节机制的靶目标研究。上述实验表明,Zm WRKY53基因通过对不同激素信号的调节,不仅参与玉米对纹枯病的抗性反应,还参与对非生物胁迫信号的转导。     

高抗旱性二倍体马铃薯材料转录组分析及抗旱基因的初步挖掘

为探究不同干旱胁迫时间处理下,二倍体马铃薯材料对干旱胁迫响应的分子机制,挖掘抗旱相关基因。以高抗旱二倍体资源A90为试验材料,利用转录组测序技术对PEG-6000胁迫下不同时间的差异表达基因进行分析,通过GO富集、KEGG通路分析与转录因子预测参与马铃薯干旱胁迫响应的差异表达基因,初步挖掘抗旱调控关键基因;并通过RT-qPCR对3个抗旱候选基因进行逆境胁迫表达分析。结果表明：A90在20%PEG-6000胁迫处理3,6,24 h与对照相比,共有2 519个差异表达基因,这些基因主要富集在植物激素信号转导、谷胱甘肽代谢、苯丙烷生物合成、戊糖和葡萄糖醛酸相互转化等抗旱相关过程中。并且显著富集在通路上的340个共差异表达基因中有53个基因注释到了RLKs、AP2/ERF和Tify等15个转录因子家族中,其中StST/K1、StERF1和StTify1的表达量较高,有可能是抗旱调控关键基因。基因StST/K1和StERF1受到干旱和低温胁迫主要在材料的根和叶中上调表达,基因StERF1在盐胁迫下在材料的根、茎和叶中均上调表达;基因StTify1在受到干旱在材料的茎和叶中上调表达,受到盐胁迫在材...     

甜菜ERF转录因子基因克隆及非生物胁迫响应分析

ERF(Ethylene-responsive element binding factors)是调控植物生长发育和响应胁迫的重要转录因子，在植物应对生物及非生物胁迫反应、调控胁迫相关功能基因的表达、提高植物的抗逆性中起着重要的作用。分析甜菜ERF转录因子基因在非生物胁迫下的表达规律，旨在为深入研究ERF转录因子在甜菜逆境胁迫应答中的作用机制提供理论依据。本研究以高糖型甜菜品系‘BS02’为试材，利用RT-PCR方法克隆得到Bv＿ammr基因，并对其进行生物信息学分析。该基因CDS区长906 bp，编码301个氨基酸。蛋白质分子量为33.19 kD，理论等电点为8.61，其二级结构以无规则卷曲和α-螺旋为主，具有一个AP2保守结构域；氨基酸序列的同源性和进化树分析显示其编码蛋白与马兜铃菌毛所编码的蛋白亲缘关系较近。采用实时荧光定量PCR方法观测非生物胁迫下该基因表达模式，结果表明，Bv＿ammr基因对低温、高温、干旱、盐、ABA等非生物胁迫有不同程度响应。     

利用WGCNA筛选马铃薯块茎发育候选基因

块茎是马铃薯的主要经济器官,解析其形成和发育机制对于马铃薯高产育种具有重要意义。虽然结薯相关基因已有报道,但仍有大量参与结薯的基因有待进一步挖掘和鉴定。随着测序技术的不断革新和发展,转录组数据愈加丰富和繁杂,利用加权共表达网络分析(weighted gene co-expression network analysis, WGCNA)筛选特定性状、组织或发育阶段相关基因的方法被广泛应用。本研究利用WGCNA分别对马铃薯DM1-3516R44(DM)材料以及RH89-039-16(RH)材料各15个组织器官的转录组数据进行分析,构建共表达网络,筛选与块茎发育相关的共表达模块。对筛选表达模块中的基因进行GO (Gene Ontology)和KEGG (Kyoto Encyclopedia of Genes and Genomes)通路富集分析。以模块中连通度前10的基因为核心基因,用PGSC数据库和NCBI数据库进行注释,并利用qRT-PCR进行验证。通过WGCNA分析,在DM和RH材料中均得到13个共表达模块,其中DM材料的Cyan模块和RH材料的Black模块与块茎显著相关。分析表明,...     

茶树CsADH基因家族的鉴定和表达模式分析

ADH转录因子在植物的生长发育和非生物胁迫响应中起重要的调控作用。本研究基于全基因组数据,从茶树基因组中鉴定到51个ADH基因,对其基因结构、进化关系、保守域、染色体定位进行分析,同时分析它们在PEG诱导的干旱胁迫、盐胁迫和冷处理中的转录组数据。结果表明:51个茶树CSADH基因在茶树染色体上分布不均;根据进化关系将茶树ADH基因分为Ⅱ类;基因结构和保守基序分析发现相同亚家族的基因结构和保守结构域基本一致;基因表达分析显示,CSADH基因在花和果实组织中具有较高的表达水平,部分基因随着叶片成熟度的增加,表达水平升高;不同的CSADH基因在PEG诱导的干旱胁迫和冷处理下存在差异表达,说明CSADH基因广泛参与茶树生长发育和响应非生物胁迫中发挥重要作用。茶树ADH基因家族的鉴定与表达分析,为深入解析ADH基因响应茶树生长发育和非生物胁迫中的功能及分子机制提供理论依据,进一步为茶树抗逆育种提供更多基因资源。     

miRNA在调控植物种子发育及响应非生物胁迫中的研究进展

microRNA(miRNA)是一类内源性小分子非编码RNA,它通过引导mRNA的裂解或抑制翻译调控靶基因在植物种子发育和响应非生物胁迫过程中起着关键作用。为了进一步鉴定和明确与种子发育及响应非生物胁迫相关的miRNAs功能和调控机制,归纳了植物中参与种子胚和胚乳发育调控及响应低温、盐、干旱等非生物胁迫的miRNAs类型、靶基因及功能。miRNA在进化上高度保守,其表达在生物发育过程中具有明显的组织特异性和时间特异性,但在不同植物之间又有着相似性。然而,目前miRNAs生物发生和功能的调控因子是如何在转录或转录后被调控的以及miRNAs是如何利用转录裂解和翻译抑制机制来调控其靶点的还有待进一步阐明。未来对这些问题的研究不仅能为植物种子发育和植物响应非生物胁迫机制提供新的见解,而且能为基因的转录后调控研究提供更多的思路。     

盐、低磷以及干旱胁迫下乌拉尔甘草根转录组分析

对乌拉尔甘草根进行转录组测序分析,挖掘响应盐、低磷和干旱胁迫的相关基因,探讨甘草抵御逆境胁迫的分子机制。以150 mmol/L NaCl模拟盐胁迫、10μmol/L KH₂PO₄模拟低磷胁迫和15%聚乙二醇(PEG-6000)模拟干旱胁迫处理下的乌拉尔甘草根作为研究材料,应用Illumina HiSeq 2000测序平台对根进行转录组测序,对测序结果进行基因功能注释分析和差异表达基因(DEGs)筛选。分别获得盐、低磷和干旱胁迫处理下甘草根中4 294、3 201和4 719个DEGs。GO富集结果表明,3种逆境胁迫下甘草根中的DEGs均显著富集到细胞过程、代谢过程、刺激响应、细胞、细胞器、催化活性等功能过程中;KEGG途径富集分析表明,3种胁迫处理下的DEGs显著富集的代谢途径主要涉及基因表达调控、生物合成及代谢、信号转导等。转录因子鉴定结果显示3种胁迫处理下差异表达数量最多的转录因子主要来自ERF、bHLH、MYB-related、WRKY、Dof、NAC等转录因子家族;基因表达水平分析显示,参与植物激素代谢及信号转导相关的DEGs在盐胁迫...     

红麻海藻糖生物合成关键酶基因HcTPPJ的克隆及响应逆境的表达分析

海藻糖生物合成关键酶基因TPP在植物响应各种非生物胁迫过程中具有重要作用。为了明确红麻TPP基因在应对非生物胁迫过程中的作用,本研究根据转录组CL541.Contig2Unigene序列设计特异性引物,通过PCR获得红麻TPP基因的cDNA全长序列。生物信息学分析表明,该基因最大开放读码框为1128 bp,编码一个含有375个氨基酸的蛋白质。序列一致性分析发现,该蛋白质氨基酸序列与其他物种TPPJ氨基酸序列的一致性为71.18%,故将该基因命名HcTPPJ。表达模式分析表明,该基因在根、茎、叶中均有表达,在盐、干旱胁迫下,随着胁迫时间的延长,HcTPPJ表达显著上调,表明该基因参与红麻的盐、干旱胁迫响应过程。在盐、干旱胁迫下, HcNCED3显著上调表达,而外源喷施脱落酸时,HcTPPJ表达变化不明显。由此推测,在红麻响应盐、干旱胁迫过程中,该基因的表达可能不受脱落酸信号分子的调控。这将为进一步阐明该基因在红麻响应盐、干旱胁迫过程的作用奠定基础。     

玉米分子伴侣基因ZmBiP2在逆境下的功能分析

BiP(binding protein)基因编码的蛋白是一类重要的分子伴侣,在动植物的逆境胁迫响应过程中具有重要的作用。从玉米抗旱自交系旱21中分离到分子伴侣基因ZmBiP2,序列分析结果显示,该基因开放阅读框长1989 bp,编码含663个氨基酸的蛋白,包含热激蛋白家族特有的TVIGIDLGTTYSC保守结构域和ATP结合位点。实时荧光定量PCR结果显示,Zm Bi P2基因在玉米的雄穗和子房中表达量最高;在盐、甘露醇胁迫条件下,该基因在植株的地上部分上调表达。过量表达Zm Bi P2基因的拟南芥转基因株系在种子萌发时期对盐和甘露醇胁迫的耐受能力减弱,在苗期对盐胁迫敏感。推测在植物响应非生物胁迫的过程中,过量表达的ZmBiP2蛋白可能发挥负调节蛋白的功能。     

山羊草谷胱甘肽S-转移酶基因家族鉴定及表达分析

谷胱甘肽S-转移酶是一种多功能蛋白酶,在植物体内参与干旱、盐、低温、重金属等多种非生物胁迫的调节;山羊草是普通小麦D染色体组的供体物种,深入挖掘山羊草中GST基因,对进一步分析六倍体小麦GST基因的功能具有重要意义。本研究利用信息生物学手段,在山羊草中共发现114条GST基因序列,分属于6个亚族;基因复制分析发现共4对基因发生了基因复制,且均为纯化选择;采用荧光定量PCR对部分GST基因在非生物胁迫下的表达分析发现,8个GST基因在响应干旱和盐胁迫时,主要在根部显著上调表达,3个GST基因在响应低温胁迫时,在根和叶中均显著上调表达,说明山羊草中的GST基因在应答非生物胁迫时,在不同组织中的表达存在着差异。     

百合响应非生物胁迫的分子机制研究

非生物胁迫严重影响了百合的生产、生长发育过程,百合有应对高温、寒冷、干旱、高盐或激素ABA等非生物胁迫的能力,百合可以通过胁迫信号感知、信号激活、信号转录和转导,并通过一系列相应胁迫的基因的表达以及生理反应来响应非生物胁迫。百合对非生物胁迫的响应过程比较复杂,响应过程分为胁迫信号的接受与感受,信号的转导,转录调控基因的表达。本研究分别从转录组测序和蛋白质组学方面,对百合响应非生物胁迫的基因表达机理及百合在非生物胁迫分子机制成果进展进行总结阐述,并讨论了培育抗逆性强的百合品种的前景,为百合分子育种提供理论参考。     

玉米发芽期响应盐胁迫的转录组分析

玉米种子发芽期对盐分较为敏感,为了进一步明确玉米发芽期对盐胁迫的应答机制,本研究利用高通量转录组测序技术,对玉米耐盐自交系昌7-2在150 mmol/L NaCl胁迫下的胚芽进行转录组测序和数据分析。结果表明：共有40 290个Unigenes获得功能注释,其中差异表达基因615个。在差异表达基因的GO富集分析中,对富集程度明显的前20个类别做出功能分析。有21个DEGs富集到苯丙烷类的生物合成途径,是KEGG富集分析最显著的通路。对差异表达基因进行COG分类,共得到25个不同的COG功能注释。通过对目标基因的挖掘共得到13个较重要的转录因子,涉及到MYB、WRKY、AP2、bZIP、bHLH和NAC等6个家族。本研究为挖掘玉米耐盐相关基因、解析玉米响应盐胁迫的分子调控机制提供科学依据。