睡莲转录因子bZIP家族的分子进化以及功能分析

【目的】基于睡莲基因组鉴定睡莲bZIP（basic leucine zipper）家族成员,并对其进行分析,以揭示睡莲bZIP家族的分子进化和功能。【方法】从Waterlily Pond数据库获取睡莲基因组序列,利用HMMER3.0程序识别睡莲的bZIP家族成员,并使用CDD程序进一步确认其含有的保守bZIP结构域,使用IQ-tree软件构建系统进化树。利用ExPASy和SOPMA在线网站进行蛋白质结构分析,通过MEME程序进行保守基序分析,使用MCScan和Circos软件对基因复制事件进行分析以及可视化展示。从NCBI下载睡莲转录组数据（SRA Study：SRP222853）,用R软件对睡莲bZIP家族成员表达数据的Pearson相关系数（PCC）进行计算和可视化分析,使用Cytoscape软件对NcbZIP成员之间的表达数据关系进行分析。【结果】从睡莲基因组中共鉴定出46个bZIP家族成员,按成员在染色体上的分布命名为NcbZIP01—NcbZIP46。根据系统进化分析可以将睡莲bZIP家族成员分为A、B、C、D、E、G、H、I、J和S共10个亚家族,其中A亚家族所含成员最多（11个）,相同亚家族成员具有相似的保守结构域和基因结构。理化性质分析表明,睡莲bZIP家族成员蛋白质长度介于101—1 898 aa,分子量大小介于12.04—214.64 kD。染色体定位分析发现,睡莲共有14条染色体,46个bZIP家族成员不均匀地分布在其中的10条染色体上,其中1号染色体上分布最多。睡莲bZIP基因家族发生10个复制事件,其中9个片段复制事件,1个串联复制事件,A亚家族所含基因复制事件最多（3次）。对NcbZIP成员在不同组织下的表达进行分析,根据表达情况分为Ⅰ、Ⅱ和Ⅲ组,Ⅰ组成员在所有组织中均高度表达,Ⅱ组成员几乎在所有组织中均不表达,Ⅲ组成员在不同组织中表达水平各不相同,其中C、D和G亚家族的大部分成员集中在Ⅲ组。通过睡莲bZIP成员表达量的Pearson相关系数分析,发现NcbZIP45与所有NcbZIP成员之间的相关性最高。【结论】在睡莲基因组中鉴定出46个bZIP成员,分为10个亚家族,不均匀地分布在14条染色体上,结构进化保守,组织表达模式多样。     

苹果bZIP转录因子家族生物信息学分析及其在休眠芽中的表达

【目的】鉴定苹果（Malus×domestica Borkh.）基因组上的bZIP基因（MdbZIP）,为研究苹果bZIP转录因子提供相关信息以及在芽休眠过程中的调控作用提供理论参考。【方法】通过Pfam下载bZIP隐马尔科夫模型bZIP_1（PF00170）与bZIP_2（PF07716）,利用HMMER 3.0鉴定苹果bZIP基因。使用Clustal Omega、MEGA6.0、MapInspect、DNAMAN 6.0和MEME4.10.2等软件对其蛋白序列进行生物信息学分析。采用Microarray分析与qRT-PCR技术检测苹果bZIP基因在不同处理下及其在高需冷量品种与低需冷量品种中的表达情况。【结果】鉴定得到120个苹果bZIP基因,与拟南芥的系统进化树分析将苹果bZIP分为10个亚家族（A-I和S）。染色体定位分析显示,109个苹果bZIP不均匀分布于17条染色体上,其中,11个基因无匹配的染色体定位。8号染色体上分布最多（13个）,1号染色体分布最少（1个）,一些染色体区域基因密度较高。基因结构分析表明,MdbZIP基因家族外显子数量0-23个,其中23个基因无内含子,分布于F亚家族（4）与S亚家族（19）,基因结构进化高度保守。保守元件分析表明,MdbZIP基因家族包含30个保守元件：元件1为bZIP保守结构域;在D亚家族发现的元件10与G亚家族发现的14为已知元件,另外多数元件功能未知。通过Microarray分析显示,多个MdbZIP均可能与芽休眠的解除相关。qRT-PCR结果显示在不同品种中A亚家族8个MdbZIP均呈现出ABA诱导表达,而D亚家族中随着冷处理时间延长,在需冷量不同的品种中出现多种表达模式。【结论】苹果bZIP基因家族结构高度保守,在ABA与冷处理下呈现不同表达模式,可能参与调控苹果芽休眠进程。     

蓝莓bZIP转录因子基因预测及其在根中表达分析

利用生物信息学方法对蓝莓bZIP转录因子基因家族进行鉴定分析,分析基因家族的基本结构、保守域以及与拟南芥bZIP基因家族之间的同源进化关系。共鉴定到131个蓝莓bZIP转录因子基因,调取blastp和HMMER检索到的序列ID,根据ID人工去冗余,后进行CDD(conserve domain database)分析,保留特征结构域筛选得到31条基因序列。根据拟南芥的分类标准将其分为9个亚族,分别是A、B、C、D、E、G、I、K、S组,每组所含有的成员数目不等。根据R和RStudio软件进行热图制作,显示不同基因在蓝莓不同组织或是不同pH条件下表达模式有差异,分析热图明确B组、C组、D组、E组bZIP转录因子家族基因成员在蓝莓根系受土壤pH胁迫下起重要抗逆调控作用。     

桑树Aux/IAA基因家族鉴定与IBA处理下表达模式分析

【目的】在全基因组范围内鉴定桑树Aux/IAA基因家族成员，并研究在IBA处理下对该家族基因表达模式的影响，为深入研究桑树Aux/IAA基因家族的功能提供依据。【方法】通过本地川桑全基因组数据对桑树Aux/IAA基因家族成员进行鉴定分析，包括蛋白质理化性质分析、系统进化树构建、基因结构分析、蛋白质三级结构及蛋白互作关联分析、启动子顺式作用元件分析。并以“强桑1号”扦插枝为试验材料，对其进行1 000 mg·L^-1 IBA处理和清水对照处理，分析处理后10、20、30、40 d的MnAux/IAA基因家族的表达模式。【结果】MnAux/IAA基因家族共有51个成员，分为8个亚家族，AtIAA基因家族分布在其中6个亚家族中；51个基因家族成员均有外显子和内含子，仅有1个成员无UTR;顺式作用元件表明，MnAux/IAA基因家族对光、激素、逆境胁迫等都有相应的调控作用；蛋白质同族进化序列的基因结构相似度较高，但族间差异较大；IBA处理后，47个基因表达量会有不同程度的提升，而基因MnAux/IAA34、MnAux/IAA33、MnAux/IAA32在表达量上有较大程度的...     

毛竹磷转运蛋白Ⅰ家族基因鉴定及表达模式

  目的  鉴定毛竹Phyllostachys edulis磷转运蛋白Ⅰ (phosphate transporter 1, PHTⅠ)家族基因,分析其表达模式。  方法  利用生物信息学方法,鉴定毛竹PHTⅠ家族成员,分析基因启动子调控元件、编码蛋白的理化性质、基因结构、氨基酸保守基序、基因在染色体的位置、组织表达特异性、基因适应性进化及系统进化等。  结果  毛竹中共鉴定出20个PHTⅠ家族基因 (PePHTs),分布在10条染色体上,均定位于细胞膜。每个基因都含有1~2个内含子,PePHTs启动子序列中包含干旱、低温等非生物胁迫以及赤霉素等激素类响应元件。毛竹PHTⅠ大部分为碱性蛋白质,分子量为48.61~76.37 kDa,理论等电点为6.84~9.30,疏水性值均大于0,都属于疏水蛋白。基因适应性进化分析显示:大多数PePHTs的选择压力值小于0,说明多数基因受到负选择压力。转录组表达图谱表明:PePHTs在不同组织中的表达存在差异性,说明该基因家族在毛竹生长发育过程中发挥着不同的作用。系统进化树表明:PePHTs都聚类在第Ⅰ亚家族,并且优先和水稻Oryza sativa聚类在同一支上。  结论  PHTⅠ家族在植物吸收和转运磷的过程中扮演了重要作用。本研究结果为深入研究毛竹PHTⅠ家族基因的功能奠定了理论基础。图5表1参45     

白花泡桐bZIP基因家族鉴定及其对丛枝植原体致病过程的响应

【目的】全基因组范围内鉴定白花泡桐(Paulownia fortunei)基因组中的碱性亮氨酸拉链转录因子(basic leucine zipper, bZIP)家族成员，并研究该家族基因在泡桐丛枝病发生过程中的作用。【方法】基于白花泡桐基因组图谱，利用隐马尔科夫模型、拟南芥bZIP蛋白质序列和本地BlastP程序，以E-value<0.001为阈值进行PfbZIP家族成员鉴定；使用MAGE7.0软件进行蛋白质序列比对及系统发育进化树构建，使用Tbtools软件进行基因结构及启动子顺式作用元件分析，并利用泡桐丛枝病发生前后及恢复过程中的转录组数据筛选与泡桐丛枝病发生相关的基因。【结果】鉴定到PfbZIP基因家族共有89个成员，根据系统进化分析将该家族基因分为13个亚家族；共线性分析发现该家族内发生了63次片段复制事件；顺式作用元件分析表明，PfbZIP基因家族可能响应光照、植物激素调控、生物与非生物胁迫等生物学过程。对泡桐丛枝病发生前后及恢复过程中的转录组数据分析发现，30个PfbZIP基因在发病前后显著差异表达，其中PfbZIP46基因表达在恢复过程中也有明显变化，故推测Pfb...     

银杏ERF转录因子家族的全基因组学鉴定及表达模式分析

以银杏基因组数据库为基础，通过序列比对，并经过鉴定筛选，共获得112个ERF转录因子家族成员，对其进行生物信息学及表达模式分析。结果显示：银杏112个ERF成员在系统进化树上分为10个亚家族；基因结构分析表明，大部分ERF家族基因结构简单，仅少量基因含有1～4个内含子；基因表达分析显示，39个ERF基因在雄蕊中表达量较高，20个基因在幼苗中表达量较高，40个基因在胚中表达量较高；部分基因在雄蕊、胚和幼苗中具有相同的表达模式，说明其可能具有类似的功能。     

单叶蔷薇bZIP转录因子家族鉴定与表达分析

【目的】筛选单叶蔷薇bZIP转录因子家族的成员，并研究其结构特点、共线性以及在不同器官中的表达模式，为进一步揭示该家族在单叶蔷薇生长发育及抗逆响应中的作用奠定基础。【方法】以采自新疆的单叶蔷薇叶片为试验材料，采用生物信息学方法，对单叶蔷薇全基因组和转录组信息进行分析，包括理化性质、系统发育、基因结构、保守结构域、染色体位置、共线性和启动子分析，以及在根、茎、叶、花、果实中的表达模式。【结果】从单叶蔷薇全基因组中共鉴定出50个单叶蔷薇bZIP转录因子家族成员，其蛋白质长度为149～700个氨基酸，分子质量为17.14～75.47 ku，等电点4.42～10.72，其中49个成员位于细胞核上。根据拟南芥bZIP转录因子家族的分类，将单叶蔷薇bZIP转录因子家族分为12亚族（A-K和S亚族），其中包括1对串联重复和7对片段重复；单叶蔷薇bZIP多含有1～15个与非生物胁迫有关的顺式作用元件。单叶蔷薇bZIP转录因子家族基因在不同器官中均有表达，各成员的表达量存在差异，其中Rbe013649在花中特异性表达，Rbe006639、Rbe028637在根中特异性表达，Rbe004215、Rbe028400、Rbe002636、Rbe002635在果实中特异性表达，Rbe011331在叶片中特异性表达。【结论】单叶蔷薇bZIP转录因子家族基因可能广泛参与各器官生长发育，其中Rbe013649可能调控花青素的合成，Rbe006639可能在抗旱过程中起重要作用。     

桑树类甜蛋白家族鉴定与生物信息学分析

为了全面分析桑树类甜蛋白家族基因特征,从桑树全基因组中鉴定类甜蛋白家族(Thaumatin-like protein,TLP)基因,并对基因结构和组成、蛋白质结构域和系统进化以及密码子使用特性进行分析。通过研究,筛选并鉴定了25个桑树TLP家族成员。多数具有典型的索马甜家族标签和保守氨基酸残基。基因结构分析结果表明,12个TLP基因家族成员含有1个内含子,10个成员含有2个内含子,且内含子相位较保守。系统进化分析结果显示,桑树TLP分为8个聚类组,组内序列一致性较高。密码子偏性分析结果显示,桑树TLP家族基因密码子第3位的核苷酸使用偏性较弱,多数基因的进化主要是由碱基随机突变造成的,属于低表达基因,少数基因受自然选择压力影响,基因表达水平较高,在植物响应环境胁迫中发挥作用。     

甜瓜GRAS家族全基因组的鉴定与表达分析

为了研究甜瓜(Cucumis melo L.)GRAS基因家族的结构特征、进化模式和生长发育过程中的表达模式，对甜瓜GRAS全基因组的基本信息、染色体位置、共线性、系统发育树、基因结构、表达模式、顺式元件进行了全面分析。结果表明，从甜瓜中鉴定出37个GRAS基因；这些基因在甜瓜的11条染色体上呈不均匀分布，并在6、7、8号染色体上出现串联重复；甜瓜GRAS家族可分为HAM、LS、SCR、SCL3、DELLA、LISCL、SHR、PAT1 8个亚族；基因结构分析显示，含有1～3个内含子的基因有8个，没有内含子的基因有29个；启动子分析，鉴定到生长素、脱落酸、赤霉素、低温、干旱、光信号响应作用元件以及参与分生组织和胚乳表达的元件；表达模式分析结果表明，HAM、DELLA亚族的基因几乎在所有组织中都有表达，LS亚族在根和茎中表达含量较高，PAT1亚族在甜瓜的整个生长发育过程中发挥着重要作用。得出结论，甜瓜全基因组共获得GRAS基因37个，共可分为8个亚族，不同的亚族具有不同的基因结构、蛋白质结构和表达模式；不同亚族在甜瓜生长发育的不同阶段具有差异性表达，它们共同参与调控甜瓜的发育过程。     

苹果属S-RNase基因的进化研究

【目的】通过分析苹果属植物自交不亲合性位点S-RNase基因的序列,研究S-RNase基因在苹果属的进化历史,序列分歧特点和遗传多态性。【方法】利用栽培苹果的S-RNase基因序列在GenBank数据库中检索和鉴定所有已知的苹果属植物的S-RNase基因,同时利用S-RNase基因的保守引物获得陇东海棠和变叶海棠的S-RNase基因序列。通过系统发育分析研究S-RNase基因的进化历史,进而计算dN/dS比值揭示S-RNase基因序列分歧的特点,最后估计并比较苹果属不同类群S-RNase基因的遗传多态性及其遗传分化。【结果】苹果属植物的S-RNase基因可以分为16个亚类,同一亚类S-RNase基因的序列分歧较小,亚类之间的分歧很大。S-RNase基因的成对dN/dS比值中有50%的大于1,并且滑动窗分析显示S-RNase基因具有多个dN/dS比值显著大于1的区域。栽培苹果和野生苹果在S-RNase基因的遗传多态性上没有明显差异。在所涉及的苹果属野生类群中,栽培苹果与塞威士苹果在S-RNase基因上的遗传分化最小。【结论】适应性氨基酸替换在苹果属植物S-RNase基因的序列分歧上发挥了重要作用。现有的S-RNase基因数据显示栽培驯化没有导致栽培苹果S-RNase基因遗传多态性的降低,基于S-RNase基因的遗传分化支持栽培苹果是由塞威士苹果驯化而来的观点。     

植物TGA转录因子研究进展

TGA基因家族是bZIP转录因子家族中十分重要的一组,其能够与靶基因启动子上的as-1区域相结合,激活或抑制下游靶标基因的转录,从而调控植株抗性或花器官发育。自烟草中首个TGA基因被鉴定以来,从拟南芥、水稻和苹果等多个物种中均有该家族基因被分离和鉴定出来。拟南芥基因组中共有10个TGA转录因子,依据其序列相似性可将其分为5组（Ⅰ组包含TGA1与TGA4;TGA2、TGA5与TGA6组成第Ⅱ组;TGA3与TGA7组成第Ⅲ组;TGA9和TGA10构成第Ⅳ组;PAN为第Ⅴ组）。其中Ⅰ、Ⅱ与Ⅲ组成员广泛参与植株的抗病反应。酵母双杂交和pull-down等结果显示,TGA1-TGA7均可与水杨酸信号途径关键调节因子NPR1相互作用。EMSA结果表明,这种相互作用能够促进TGA对下游PR1启动子的结合并上调其表达,提高植株抗病性。但这3组基因在植物基础抗性与系统获得抗性中的作用有所不同：tga3突变体对病菌的基础抗性存在缺陷,但植株诱导抗性却并不受影响;TGA1与TGA4对基础抗性和系统抗性均有一定影响;TGA2、TGA5与TGA6之间存在功能冗余,仅tga2/5/6三突变体才表现出与npr1突变体类似的系统获得性抗性缺乏的表型。酵母双杂交筛选发现：Ⅱ组TGA能够与谷氧还蛋白GRX480相互作用并介导SA对JA途径标记基因的抑制作用,同时,该组蛋白还能够与GRAS家族蛋白SCL14相互作用,增强下游CYP81D11与GSTU7等解毒相关基因的表达,以一种不依赖于NPR1的信号途径提高植物对外源化学物质毒害的耐性。此外,tga1/4双突变体与nrt2.1/2.2双突变体的初生根和侧生根生长在低氮条件下较野生型显著下降,ChIP和酵母单杂交结果显示,TGA1能够与硝酸盐转运蛋白基因NRT2.1及NRT2.2的启动子相互结合,通过调节这两个基因的表达来调节植物的氮响应。而TGA3在镉长距离运输中发挥重要作用。Ⅳ与Ⅴ组成员在调控花器官发育中具有重要作用。tga9/10表现出与roxy1/2双突变体类似的花药发育缺陷的表型;PAN能够与NPR1类蛋白BOP1和BOP2相互作用,并且pan与bop1/2双突变体均可表现出5枚萼片的表型,暗示花发育与抗病可能具有类似的信号调节机制。在文章最后介绍了翻译后修饰对TGA功能的影响,并对TGA未来的研究方向进行了探讨,以期为该领域研究者提供参考。     

辣椒HD-Zip基因家族鉴定、系统进化及表达分析

【目的】鉴定辣椒HD-Zip基因家族,并利用生物信息学方法系统分析其在基因组中的分布、基因结构、进化分化特征及在不同组织中的时空表达特异性,解析该家族的进化特征及生物学功能。【方法】根据已报道及PlantTFDB数据库中的拟南芥HD-Zip序列,利用本地BLAST工具在我国辣椒测序品种‘遵辣1号’基因组中比对,并利用Pfam、SMART工具进一步验证。采用EMBOSS Programs、MEGA、GSDS、MEME、MCScanX、OrthoMCL、Circos等软件预测辣椒HD-Zip基因家族成员蛋白理化性质,构建系统进化树,定位染色体,分析基因结构、基因复制类型及直系、旁系同源基因。基于GEO数据库,运用R软件、本地perl语言及Cytoscape分析辣椒HD-Zip组织表达差异并绘制共表达网络。【结果】本研究在‘遵辣1号’基因组中鉴定获得42条辣椒HD-Zip,命名为CaHDZ01—CaHDZ42。CaHDZs长度跨度较大,70% CaHDZ蛋白的pI小于7.0。除CaHDZ42,其余基因不均匀地分布在12条染色体上,部分基因为片段复制。该基因家族可分为4个亚族,分别含有18、9、5、10个HD-Zip,基因结构及蛋白结构域差别显著。辣椒、番茄和拟南芥3个物种中的直系同源基因对数目大体相同,但同为茄科的辣椒和番茄之间的稍多;辣椒中的旁系同源基因少于番茄和拟南芥,说明辣椒基因组的倍增事件并没有使CaHDZs明显扩增。对无油樟、水稻、玉米、番茄、马铃薯、辣椒‘CM334’、辣椒‘Zunla-1’、毛果杨、葡萄以及拟南芥9个代表物种的HD-Zip进化特征分析结果表明,从被子植物开始,HD-Zip基因家族就稳定存在4个亚族。推测在形成4个亚族前,HD-Zip分为两组,其中一组分化成I和II亚族,而另一组则分化成为III和IV亚族。CaHDZs在根、茎、叶、花芽、花和果实不同发育时期的表达模式分析结果显示,4个亚族具有不同程度的表达趋势。其中I亚族基因在辣椒不同组织中的表达量均较高,且不同成员间表达模式不同,CaHDZ22在茎中的表达最高,表明该基因可能对辣椒茎的生长有重要作用。II、III和IV亚族基因在不同组织中的表达量相对较低,但部分基因在特定组织中具有较大的表达量。如CaHDZ34在辣椒果实成熟后期具有较大高的表达量,CaHDZ02和CaHDZ28在果实膨大时表达较高,CaHDZ04在果实成熟前期具有较高的表达量。CaHDZs表达网络中有33对基因表达趋势的相关系数（PCC）大于0.8,6对大于0.9,表明CaHDZs协同调控了辣椒的生长发育,不同亚族之间也具有协同性。【结论】在‘遵辣1号’基因组中鉴定获得42条CaHDZs,可分为4个亚族,不同亚族的基因结构、蛋白保守结构域及表达模式不同。在进化过程中,辣椒HD-Zip保守性高,数目没有明显扩增,I和II亚族、III和IV亚族关系更近。CaHDZs具有组织表达差异性,协同调控了辣椒的生长发育。     

玉米大斑病菌Homeobox转录因子家族全基因组鉴定及其表达规律分析

【目的】从全基因组水平上鉴定玉米大斑病菌(Setosphaeria turcica)Homeobox转录因子家族及其分布,阐明该家族的序列及进化特征,分析该家族基因在病菌不同生长发育时期的表达规律。【方法】利用生物信息学手段搜索玉米大斑病菌全基因组数据库,鉴定Homeobox转录因子家族;采用MEGA 5.0软件进行系统进化树分析;利用在线工具GSDS(gene structure display server)(http://gsds1.cbi.pku.edu.cn/index.php)绘制基因结构图;利用Clustal X 1.83软件分析Homeobox保守结构域(HOX保守结构域)的氨基酸序列特征;利用SOPMA(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.plpage=npsa_sopma.html)对Homeobox蛋白的二级结构进行在线预测;利用实时荧光定量PCR(q RT-PCR)技术分析Homeobox转录因子家族在病菌不同发育时期的表达模式。【结果】在玉米大斑病菌中鉴定了8个Homeobox转录因子家族成员(St HTF1-8),根据基因结构及系统进化特征将其分为4类;亚细胞定位预测分析表明,这8个蛋白全部定位在细胞核中;该家族成员均含有HOX保守结构域,其二级结构具有特征性的"螺旋-转角-螺旋"(helix-turn-helix)结构;利用q RT-PCR技术对该家族成员在菌丝、分生孢子形成、芽管形成、附着胞及侵入丝形成等5个时期的表达规律分析,发现不同基因在病菌不同发育时期具有不同的表达水平,其中St HTF1在菌丝发育、分生孢子及附着胞形成等3个时期的表达水平相对较高,St HTF3、St HTF4在分生孢子形成时期表达水平最高,St HTF6在芽管形成时期的表达水平最高,St HTF2、St HTF5、St HTF7和St HTF8在附着胞形成时期的表达水平均较高。【结论】玉米大斑病菌包括Homeobox转录因子家族包含8个成员,在进化上分为4大类,全部成员均分布在细胞核内,其编码蛋白质均含有保守的HOX结构域及"螺旋-转角-螺旋"空间结构;该基因家族成员在病菌不同发育时期呈现不同的表达规律。     

烟草bZIP转录因子BZI－4的生物信息分析

[目的]探讨烟草在采收前的生理期和烘烤过程中的变黄阶段抗逆机制。[方法]利用生物信息学方法对烟草b ZIP转录因子BZI-4的理化特性、信号肽、跨膜、系统进化和亚细胞定位进行分析。[结果]烟草b ZIP转录因子BZI-4为碱性、水溶性、耐高温和非分泌型的游离蛋白,该蛋白为叶绿体中的转运肽,仅剪切修复后才具有活性,对环境适应能力强。[结论]通过烟草b ZIP转录因子BZI-4研究,为烟草抗逆境胁迫研究提供理论基础。     

桃NAC家族基因生物信息学分析及其响应低温胁迫的表达特征

NAC(NAM-ATAF1/2-CUC2)是植物特有的具有多种生物功能的一类重要转录因子,在植物应对非生物胁迫、激素信号应答与器官形成中发挥重要作用。为了解桃NAC蛋白(PpNAC)及其在桃低温胁迫过程中的响应,采用生物信息学方法和qRT-PCR系统分析了桃NAC家族基因的成员、结构、功能和低温胁迫下的表达。结果表明,桃NAC家族共有119个成员,命名为PpNAC001~PpNAC119。PpNACs的氨基酸长度在68~862 aa,平均氨基酸长度为352 aa,分子量为7.2~96.1ku,等电点为4.3~9.5。除PpNAC119定位在染色体骨架上之外,其余基因均定位在桃的8条染色体上。基因结构和保守基序分析表明,位于同一亚族的成员基因结构相似,外显子数量和长度基本相同,并且同一亚族中的成员所含motif的种类和排列顺序基本相似,推测同一亚族的成员可能具有相同的功能,而不同亚族之间的差异可能与他们的功能特异性有关。对PpNACs蛋白启动子区2 000 bp的序列进行顺式作用元件分析,结果表明,脱落酸响应原件、光响应原件、茉莉酸甲酯响应原件、赤霉素响应原件、低温响应原件和生长素响应原件在筛选到的所有原件中占比较大。qRT-PCR结果显示,在4 ℃低温胁迫2、6 h,部分基因表达量上调不明显;低温胁迫12 h,16个NAC基因的表达量均显著上调。