[目的]探究MYB基因家族成员在桃雄性不育中的表达模式,以筛选出与花粉发育相关的桃MYB家族成员。[方法]选择桃可育品种红芙蓉和不育品种沙红,以其红盛花前30,24,18,15 d(依次对应S1、S2、S3、S4时期)的花芽及盛花期花蕾为材料,对花药形态进行观察与比较;对2个桃品种S1～S3时期花芽中差异表达的MYB基因家族成员进行鉴定,并根据其FPKM值绘制热图;利用MEGA6软件,对桃花芽发育过程中差异表达的MYB转录因子和拟南芥、水稻等物种的28条MYB蛋白序列进行多重比对,构建系统进化树;从桃MYB基因家族成员中选择6个(Prupe.3G046900、Prupe.1G551400、Prupe.4G192000、Prupe.1 G323400、Prupe.1G273900和Prupe.2G192100)表达趋势可能与育性相关并且与其他物种同源性较高的成员,分析其在不同育性桃品种中的表达特性。[结果]在盛花前18 d(S3时期)时,红芙蓉与沙红花药形态出现差异,红芙蓉花药最终表现为橙红色且较为饱满,而沙红则表现为淡黄色且瘪小。39个MYB基因家族成员在沙红和红芙蓉之间差异表达,其中27个属于R2R3-MYB类转录因子;根据FPKM值聚类热图,39个MYB基因家族成员被分为3簇。由系统进化树可知,所有MYB蛋白被分成16簇,其中Prupe.3G046900与AtMYB26、Prupe.3G017100与AtMYB99、Prupe.3G0006300与AtMYB108同源性较高;而Prupe.1G276300与AtMYB25和AtMYB1,Prupe.1G551400、Prupe.1G323400与AtMYB4和AtMYB32分别单独成为1簇。Prupe.3G046900在可育品种各时期花芽中的表达量高于不育品种,Prupe.1 G551400、Prupe.1G323400、Prupe.4G192000于花粉败育发生前后在不育品种花芽中的表达量高于可育品种。[结论]推测桃MYB基因家族成员Prupe.3G046900、Prupe.1G551400、Prupe.1G323400、Prupe.4G192000参与桃花粉发育的调控。

[目的 ]鉴定白桦AMT基因家族成员，分析AMT基因的表达模式。[方法 ]本研究利用生物信息学方法鉴定了该家族成员并通过实时荧光定量技术分析其基因表达模式。[结果 ]从白桦基因组中鉴定了9个AMT家族成员，并将其分为AMT1和AMT2两个亚家族，分别命名为BpAMT1.1～1.4和BpAMT2.1～2.5;BpAMTs氨基酸残基数为384～522个，等电点为4.61～8.16，均定位于质膜与细胞器膜上；BpAMT基因不均匀的分布在5条染色体上，且成员间存在串联重复现象。BpAMT基因的表达具有组织部位特异性，呈现叶>根>茎趋势；同时发现，硝酸钾、氯化铵、茉莉酸甲酯、赤霉素、脱落酸、氯化镉、干旱、4℃和日变化等均可以影响BpAMT基因表达，且不同处理下家族成员的响应程度存在差异。[结论 ]从白桦中鉴定9个BpAMT基因，聚类为两个亚家族，在调节氮素吸收转运、响应激素信号及非生物胁迫中发挥不同的作用，该研究结果为深入解析BpAMT基因在白桦生长发育和抵抗逆境中的功能奠定基础。

【目的】GRAS基因家族成员在调节植物生长发育中发挥着关键作用。通过生物信息学分析GRAS在桃基因组中的分布、结构及进化,研究家族成员在不同组织中的表达特异性及其对UV-B的响应,解析GRAS基因家族的生物学功能。【方法】对设施油桃‘中油5号’(Prunus persica var. nectarina cv. Zhongyou5)补充适量UV-B(Ultraviolet-B)剂量,利用Plant TFDB数据库鉴定桃GRAS基因家族成员。采用Clustal W、MEGA6.0、ProtParam tool、MCScanX、Circos、SMART、NCBI-CDD、ExPASy、GSDS和MEME等软件构建系统进化树,绘制染色体定位图,预测蛋白的相对分子质量与等电点等理化性质等,分析GRAS基因家族成员在不同组织中的表达模式,利用qRT-PCR技术检测桃GRAS在UV-B处理下的表达情况。【结果】从桃全基因组中鉴定出48个GRAS转录因子家族基因,构建系统进化树将这48个成员分为9个亚家族,PpGRAS在桃的8条染色体上呈不均匀分布。对GRAS基因家族进行理化性质分析发现,其蛋白平均长度为590.52 aa,等电点在4.36—7.56。GRAS家族基因结构分析表明,有40个基因不含内含子,8个基因含有1个内含子。保守元件分析显示,GRAS家族包含20个保守元件,其中Motif 2和4在GRAS的家族中高度保守,同一个亚家族成员含有相同的保守元件,其可能具备相似的功能。然而,有些亚家族成员的表达模式不同,这可能与其保守基序之外的序列有关。PpGRAS在不同组织中具有不同的表达模式。叶片中PpGRAS5经UV-B处理后上调表达最显著,而有多达15个基因表达量呈现下调。果实中PpGRAS13经UV-B处理后上调表达,但有9个基因表达下调。韧皮部中,UV-B处理后有14个基因上调表达,而PpGRAS16经处理后在韧皮部中表达下调最明显。【结论】从桃基因组中共鉴定出48个GRAS基因家族成员,分布于8条染色体上;多数PpGRAS能响应UV-B处理,但在不同组织中的表达不尽相同。

【目的】从全基因组水平上鉴定巨桉生长调节因子(growth-regulating factor,GRF)基因家族成员,分析其在不同氮素水平下的组织表达模式,为巨桉GRF基因功能研究及氮高效利用桉树品种的培育奠定基础。【方法】在NCBI网站中选择巨桉基因组,用拟南芥和毛果杨的GRF蛋白序列与巨桉基因组数据库进行BLAST蛋白序列同源比对,并通过InterPro和SMART数据库进行保守结构域分析,获得巨桉GRF基因。利用在线网站Ex-pasy protparam、SignalP-5.0、WoLF PSORT和SOPMA,对巨桉GRF蛋白进行基本理化性质分析、信号肽预测、亚细胞定位及二级结构预测。利用MEGA7.0软件构建系统进化树,使用MEME在线平台分析EgrGRF蛋白的保守结构域和基序,用Plant CARE预测基因上游的顺式作用元件,并使用TBtools软件将结果可视化。以2年生巨桉苗为供试材料,浇灌氮素浓度分别为45(高氮)、15(常规氮)和1.5(低氮)mmol/L的营养液,4 d后取样,采用实时荧光定量PCR技术检测EgrGRF基因在3种氮素水平下的叶片、茎和根中的表达情况。【结果】鉴定得到6个巨桉GRF基因家族成员(EgrGRF1～EgrGRF6),分布于4条染色体上;6个EgrGRF蛋白的氨基酸数目为296～605个,分子质量为33 415.86～64 630.89 u,理论等电点为7.29～8.85,脂溶指数为39.93～64.79,均为亲水性蛋白;无规则卷曲和α-螺旋为主要二级结构元件,均定位于细胞核上,未检测到信号肽存在。系统进化树分析表明,巨桉、毛果杨和拟南芥的GRF家族成员可分为4组,各组中EgrGRF蛋白数量分别为0,3,1和2个,多数EgrGRF成员与毛果杨亲缘关系较近。基因结构及蛋白基序分析结果显示,巨桉GRF基因家族成员具有2～4个内含子;6个EgrGRF蛋白均具有CX9CX10CX2H和QX3LX2Q保守基序。顺式作用元件分析表明,EgrGRF基因启动子区含有茉莉酸甲酯、脱落酸、分生组织表达及低温等响应元件。实时荧光定量PCR结果显示,EgrGRF2和EgrGRF6基因在低氮处理的巨桉叶片中优势表达,EgrGRF2、EgrGRF3、EgrGRF5和EgrGRF6基因在低氮处理根中的表达量较高氮和常规氮处理显著升高,EgrGRF4基因在常规氮处理茎中的表达量最高。【结论】鉴定获得6个巨桉GRF基因家族成员,其中EgrGRF4基因可能参与茎生长发育的调控,EgrGRF2、EgrGRF3、EgrGRF5和EgrGRF6基因可能参与了巨桉对低氮胁迫的响应。

【目的】探究水曲柳光敏色素互作因子(FmPIFs)在激素调控与非生物胁迫应答过程中的重要作用,为揭示水曲柳抗逆分子机制和制定林木遗传育种策略提供理论依据。【方法】在水曲柳中克隆FmPIFs基因,并对其基因结构、蛋白质理化性质、保守基序、系统进化关系等进行生物信息学分析。采用qRT-PCR方法分析水曲柳中FmPIFs基因在不同组织及不同激素与胁迫条件下的表达模式。【结果】获得5个水曲柳FmPIFs基因家族成员,分别命名为FmPIF1、FmPIF3、FmPIF4、FmPIF7和FmPIF8,其对应的蛋白均为亲水性不稳定蛋白,全部定位在细胞核内。多序列比对结果表明FmPIFs均存在APB保守结构域,成员FmPIF1与FmPIF3存在特有的APA结构域。组织特异性分析显示FmPIFs均在叶中表达量最高,其中成员FmPIF8表达量最高,为对照的3.96倍;但在茎中少量表达,茎中表达量最高的是FmPIF3,仅是对照的0.21倍;而在根中表达量均极低。胁迫响应分析表明,FmPIFs正调控水曲柳植株盐、碱和干旱胁迫抗性,而对植株抗寒性起负调控作用,其中成员FmPIF3对寒冷及盐胁迫明显响应,FmPIF8在碱胁迫下表达量明显上调,FmPIF1在干旱胁迫下表达量明显上调。在激素响应结果中,FmPIFs对脱落酸(ABA)、水杨酸(SA)和赤霉素(GA_3)的响应较为一致,而对生长素(IAA)和茉莉酸甲酯(MeJA)的响应存在差异。FmPIF1在施加MeJA后剧烈应答且表达量显著上调,FmPIF7在SA处理后表达量明显上调,FmPIF3和FmPIF4在GA3处理后表达量明显上调。【结论】FmPIFs各成员在基因及蛋白结构上表现出较高的一致性。RT-qPCR结果表明,FmPIFs在水曲柳叶片中表达量最高。在盐、碱、干旱和寒冷胁迫下,FmPIFs被诱导表达,且大部分表达模式相似。FmPIFs也在水曲柳响应IAA、ABA、MeJA、SA、GA3激素调控过程中发挥重要作用。

【目的】深入挖掘巨桉脱水素基因EgrDHN家族成员并探讨其应对逆境胁迫的表达模式，从中筛选具有抗逆功能的EgrDHN基因对于培育桉树抗逆新品种具有重要意义。【方法】以巨桉全基因组数据为参考，利用生物信息学筛选鉴定得到巨桉EgrDHN基因家族成员，对其进行生物信息学分析，具体包括进化树分析、蛋白序列比对、启动子顺式作用元件预测等，并基于生信分析结果，对EgrDHN基因家族成员进行低温、高盐和ABA胁迫处理，采用RT-qPCR测定其表达模式，同时对其组织表达特异性进行了测定。【结果】从巨桉基因组中共鉴定出EgrDHN家族9个成员，将其命名为EgrDHN1～EgrDHN9，分布于5条染色体上，根据蛋白质保守序列的特点可分为K_nS、SK_n和Y_nSK_n 3种类型。系统进化树显示植物DHN家族成员可分为4个进化枝，其中第Ⅳ进化枝包含巨桉EgrDHN7、EgrDHN8、EgrDHN9，但不包含任何毛果杨Populus trichocarpa DHN同源基因（PtrDHNs），两者存在明显差异。启动子元件预测分析发现巨桉EgrDHN启动子中含有低温响应（low-temperature response,LTR）等胁迫响应类元件，MYB、MYC等与抗逆性相关以及植物激素响应类顺式作用元件。组织特异性表达分析发现EgrDHN家族成员在根、成熟叶片、子叶中表达量较高，其中EgrDHN2在成熟叶片中表达量最高。非生物逆境胁迫下经RT-qPCR检测发现EgrDHN家族成员中，EgrDHN1在高盐胁迫下表达量下调，在低温和ABA胁迫时上调，而其他成员在低温、高盐、ABA胁迫下均存在上调，其中低温条件下以EgrDHN2、EgrDHN3、EgrDHN4响应最为明显。【结论】大部分EgrDHNs基因参与逆境胁迫响应，且除EgrDHN1在盐胁迫下的表达量下降外其它成员在逆境胁迫下的表达量均上升。

[目的]探究bHLH亚家族成员MYC基因在杨树生长发育及对环境响应中的表达模式,为解析茉莉酸信号通路关键基因MYC调控杨树生长发育及抗逆响应提供参考。[方法]利用生物信息学方法鉴定杨树基因组中的MYC基因家族成员并对各成员的基因结构、保守基序进行了系统分析;在此基础上,采用实时荧光定量PCR技术检测了MYC成员在不同组织、不同植物激素响应及逆境胁迫下的表达模式。[结果]毛果杨基因组中含有10个MYC成员,这些成员在进化上相对保守,均具有典型的bHLH结构域,根据进化关系将其分为3个分支。组织表达分析发现,大部分MYC家族成员在根中均具有较高表达量,分枝Ⅱ两对基因在茎中呈相反表达模式。不同激素处理及非生物胁迫条件下MYC基因存在显著表达差异,但基因对之间多存在类似的表达模式。[结论]杨树MYC基因家族成员可能参与不同的生物学过程,研究结果为深入解析杨树MYC基因的功能奠定了基础。

【目的】对毛竹Phyllostachys edulis ICE基因家族进行鉴定及分析，找出响应毛竹抗寒关键家族成员，研究毛竹ICE基因的生物学功能、响应低温胁迫的分子机制及遗传转化，为提高毛竹抗寒性奠定理论基础。【方法】利用生物信息学方法分析毛竹ICE基因家族成员，并对4、0、-2℃低温处理0(对照)、0.5、1.0、24.0、48.0 h的毛竹生理指标和ICE基因的表达模式进行分析。【结果】共鉴定了4个毛竹ICE基因。保守结构域和多重序列比对分析表明：PeICE基因结构高度相似。系统发育关系及启动子顺式作用元件分析显示：PeICE基因与水稻Oryza sativa亲缘关系更近，同时存在大量与非生物胁迫相关的顺式作用元件。活性氧自由基(ROS)染色发现随着处理时间增长，ROS染色逐渐加深，但是其0℃处理24.0 h、-2℃处理1.0 h后染色逐渐减弱。脯氨酸(Pro)质量摩尔浓度、超氧化物歧化酶(SOD)活性显示：4和0℃条件下，Pro质量摩尔浓度和SOD活性整体增加，但-2℃时低于对照。过氧化物酶(POD)活性显示：在3个低温处理下均增加。ICE基因表达模式分析发现：4、0℃处理时PeICE表达量整体增加，且都以PeICE3增量最明显；而-2℃处理下PeICE整体表达量水平低于对照。【结论】随着温度降低和处理时间增强，毛竹受到的损伤不断增强，其内酶活系统以及ICE基因积极响应低温胁迫，其中，PeICE3对低温胁迫最为敏感，但在-2℃时，ICE基因表达量并未增加，推测该基因家族响应了寒冷胁迫而非冷冻胁迫。图7表1参42

12-氧-植物二烯酸还原酶(OPR)为黄素单核苷酸(FMN)依赖的氧化还原酶，是合成茉莉酸的关键酶，其对植物生长发育及防御调控具有重要作用。为研究OPR基因在玉米中的抗病作用，利用生物信息学方法，在31个玉米不同自交系中对OPR家族成员进行鉴定，并分析受玉米大斑病菌侵染后的表达规律。结果表明，在玉米B73品系中鉴定出了8个OPR基因，这些基因不均匀地分布于7条染色体上，且其表达蛋白富含酸性氨基酸。进一步分析发现，玉米OPR家族只有1种结构域Oxidored＿FMN,并且所有成员均包含已鉴定的10个蛋白质保守基序。利用MEGA软件对玉米、小麦、水稻和拟南芥的OPR家族成员进行系统发育关系分析，发现这些植物OPR基因家族进化关系具有明显差异，其中玉米和水稻的进化关系最为接近。运用OrthoFinder软件对其同源组分析发现，玉米OPR基因家族具有高度保守性，所有OPR基因均为核心基因，但其功能略有差别。进一步利用河北省植物生理与分子病理学重点实验室前期的RNA-seq数据，对玉米B73品系OPR基因响应大斑病菌侵染的表达模式进行分析，并通过qRT-PCR对基因表达模式进行验证，发现这些基因在病菌侵染玉米过程中呈现3种不同的表达模式。本研究系统的鉴定了玉米泛基因组OPR家族基因，以及确定了其在应对玉米大斑病菌侵染后的表达模式。

【目的】探究赤苍藤WRKY转录因子在不同组织及其生长发育过程中的功能作用，为深入地解析赤苍藤生长发育的分子调控机制提供理论依据。【方法】利用生物信息学方法对赤苍藤WRKY基因家族进行全基因组鉴定，分析赤苍藤WRKY基因家族成员的结构和功能，并通过转录组测序分析WRKYs基因在不同组织的表达模式。【结果】共鉴定66个EsWRKYs基因，在11条染色体中均有分布，亚细胞定位预测有62个EsWRKY蛋白定位在细胞核，系统发育分析将赤苍藤WRKY蛋白家族成员分为三大类（Ⅰ～III），分别有18、44和4个成员，其中Ⅱ类可分为5个亚类（Ⅱa～Ⅱe），分别有3、7、18、8和8个成员，66个EsWRKYs基因的外显子有2～11个，保守基序分析结果显示Motif1和Motif3包含WRKY结构域，Motif1是EsWRKYs蛋白最主要的保守基序，EsWRKYs基因家族成员启动子区有与生长发育、逆境响应、激素响应等相关的19种顺式作用元件。表达模式分析显示EsWRKY14、EsWRKY22、EsWRKY29、EsWRKY30、EsWRKY47和EsWRKY61可能参与赤苍藤茎叶的生长发育过程，多数EsWRKYs基因在种仁发育过程中显著下调，表明这些EsWRKYs基因可能在种仁发育前期发挥主要作用。【结论】66个赤苍藤WRKY基因家族成员的理化性质具有较大差异，均具有WRKY保守结构域，EsWRKYs基因在不同组织中表达模式不同，在生长发育、逆境胁迫及次生代谢过程等方面可能发挥重要作用。

【目的】筛选、鉴定并分析金鱼草萜烯合酶（Terpenoid synthase， TPS）基因家族成员，确定金鱼草主要花香挥发性萜烯类物质合成的TPS基因，为进一步解析该基因功能奠定基础。【方法】基于生物信息学分析，鉴定AmTPS基因家族成员，并对其基因结构、蛋白保守基序、顺式作用元件、系统进化、染色体定位和共线性、蛋白质理化性质、亚细胞定位预测等进行分析；通过实时荧光定量（qRT-PCR）分析TPS候选基因在不同组织中的表达水平；最后利用气相色谱-质谱联用仪（GC-MS）对金鱼草花香挥发性成分进行分析鉴定。【结果】本研究共鉴定出30个AmTPS候选基因，基因结构及保守基序显示AmTPS基因外显子数量为3～13个，保守基序DDxxD是主要基序，出现31次。AmTPS启动子的顺式作用元件主要集中在生长与发育类元件中，抗逆、胁迫响应元件较少。系统发育分析显示，AmTPS共分为4个亚家族TPS-a，TPS-b，TPS-e/f，TPS-g，且TPS-a、TPS-g亚家族成员较多，TPS-e/f成员仅1个。蛋白理化性质结果表明，预测氨基酸数量（AA）为291～1124 aa，分子量（MW）为34.392 ～131.778 kDa，理论等电点（pI）为4.98 ～ 6.98，疏水平均（GRAVY）值为-0.464～-0.24。qRT-PCR表达分析显示，AmTPS在花器官中的转录水平普遍高于其它部位，特别是TPS-g亚家族基因。利用GC-MS分析金鱼草花香成分，发现TPS-g亚家族基因合成的无环单萜类物质相对含量占比80.51%，与该家族基因表达模式吻合。【结论】金鱼草AmTPS基因家族成员在进化上比较保守，TPS-g亚家族的表达具有组织特异性，在金鱼草花香挥发性无环单萜的大量释放中发挥重要作用。本研究为进一步探究金鱼草萜烯类花香挥发性化合物提供参考。

【目的】鉴定与分析甘薯KUP/HAK/KT基因家族，为进一步深入对IbHAK基因功能研究提供理论基础。【方法】基于生物信息学手段，利用泰中6号基因组数据，在全基因组水平对甘薯KUP/HAK/KT基因家族进行鉴定，同时进一步对家族成员的基本特征、系统进化、基因复制和表达模式进行分析。【结果】在甘薯中共鉴定到22个IbHAK 基因，蛋白质长度为227aa-1252aa，分子量为18 297.95～140 640.46 u，等电点为5.54～9.31，全部定位于质膜上，蛋白完整的二级结构有4～13个跨膜区。依据系统进化关系，将其分为ClusterⅠ、ClusterⅡ、ClusterⅢ和ClusterⅣ4个进化簇。染色体定位结果显示22个基因分布在11条染色体上，其中Chr4染色体上分布最多，有4个IbHAK基因，22个IbHAK基因中共产生8个复制基因。22个IbHAK基因都含有TATA和CAAT控植物生长发育相关元件以及光响应元件，且所有的IbHAK基因至少含有1种逆境响应元件和1种激素响应元件。IbHAK基因呈现出组织表达特异性，其中IbHAK6、IbHAK8、IbHAK9、IbHAK22在根中表达量最高，逆境胁迫下多个IbHAK基因出现表达量变化。【结论】从甘薯全基因组鉴定22个KUP/HAK/KT基因家族成员，甘薯KUP/HAK/KT基因家族成员在进化上比较保守，IbHAK6、IbHAK8、IbHAK9、IbHAK22在甘薯根形成过程中发挥作用，大部分IbHAK基因能够响应逆境胁迫。

腺苷三磷酸结合盒转运蛋白(ATP binding cassette transporter， ABC transporter)家族是最古老的膜蛋白家族之一，广泛存在于原核生物和真核生物体内，利用ATP水解释放的能量对氨基酸、脂质、抗生素等物质进行跨膜运输，参与生物体内多种生理活动。目前，对软体动物ABC转运蛋白家族的鉴定，仅在虾夷扇贝(Patinopecten yessoensis)、栉孔扇贝(Chlamys farreri)和长牡蛎(Crassostrea gigas) 3种双壳类中有系统研究，对缢蛏(Sinonovacula constricta) ABC转运蛋白家族的鉴定和表达模式分析尚未见报道。利用缢蛏基因组和转录组数据，运用本地及NCBI在线BLAST程序、FGENESH+、SMART、ExPASy、MEGA X、Mev 4.90和MapInspect等生物信息学分析工具，在全基因组水平系统地鉴定出52个缢蛏ABC转运蛋白，并对转运蛋白基因外显子数目、染色体定位等信息进行分析；通过系统进化分析，将缢蛏ABC转运蛋白家族分为8个亚家族，即ABCA～ABCH；通过比较分析不同物种ABC转运蛋白家族，推测基因串联复制事件对软体动物ABC转运蛋白家族成员数目的增加有影响。ABC转运蛋白基因在缢蛏不同发育时期和组织的表达分析显示，ABCC和ABCG亚家族多个基因在缢蛏稚贝时期表达量达到最高值，ABC转运蛋白基因在鳃、肝胰腺与性腺中表达种类较多。软体动物作为第二大动物类群，目前仅有3个物种的ABC转运蛋白得到了系统分类，缢蛏ABC转运蛋白基因家族的系统鉴定与表达模式分析为研究软体动物ABC转运蛋白基因的进化和缢蛏ABC转运蛋白功能提供了基础资料。

【目的】探究桑树NAC基因家族的生物信息学特征，并研究1 000 mg·L^-1生根粉（ABT）处理对该基因产生的影响，为深入研究桑树NAC转录因子的功能提供依据。【方法】利用生物信息学方法对川桑NAC基因家族进行基因鉴定，并分析其理化性质、保守基序、顺式作用元件、系统进化树、蛋白三级结构及互作关系，并选用‘果桑大十’作为试验材料，经ABT处理后，分析其4个时期（10、20、30、40 d）NAC基因家族的表达模式。【结果】桑树NAC基因家族共有92个成员，分为22个亚家族，亚家族成员广泛分布于细胞核，且具有相似的保守结构。顺式作用元件分析表明，桑树NAC基因家族对激素具有较强的响应能力。蛋白三级结构显示同一亚家族的蛋白空间结构相似，且功能相同，这可能与其进化同源性有关。ABT处理后，NAC基因家族的表达水平整体呈上升趋势，而MnNAC91、MnNAC67、MnNAC28在对照组（CK）中呈现较高的表达水平。【结论】NAC基因家族功能相对稳定，在整个进化过程中无明显变化；ABT可促进NAC基因家族表达，尤其是MnNAC75、MnNAC104、MnNAC7、MnNAC30、MnNAC101、MnNAC83、MnNAC102具有较高的同源性，并与拟南芥的结构功能相似，推测对抗干旱胁迫有一定作用。

采用生物信息学方法从生防菌短短芽孢杆菌(Brevibacillus brevis)X23全基因组水平鉴定出13个TetR家族基因，将其命名为BbTetR1、…、BbTetR13，对其基本性质、基因结构和高级结构进行分析。结果表明：短短芽孢杆菌TetR家族成员在进化上具有保守性，基因序列长度为477～681 bp，编码的蛋白由168～226 个氨基酸组成，相对分子质量为1.9×10~(4)～2.57×10~(4)，等电点为4.84～7.77；其基因家族成员可划分为3个亚族，亚族Ⅰ包含5个成员，亚族Ⅱ、亚族Ⅲ各包含4个成员；所有BbTetR蛋白均含有大量α–螺旋，且这些蛋白的空间结构呈现多样化；转录组分析结果表明，不同的BbTetR在短短芽孢杆菌不同生长时期的表达模式和表达量差异显著，BbTetR1、BbTetR5、BbTetR8和BbTetR10在短短芽孢杆菌不同生长阶段均有较高表达，推测其可能发挥更重要的调控作用。