茶树ZF-HD转录因子基因家族的鉴定及表达分析

【目的】搜索鉴定茶树ZF-HD（Zinc finger homeodomain）转录因子基因家族成员,并对其进行生物信息学分析及在不同组织及非生物胁迫和激素处理下的表达模式,为茶树ZF-HD基因家族的功能研究及茶树遗传性状改良提供理论依据。【方法】通过隐马尔可夫模型预测和序列比对从茶树基因组中筛选鉴定出茶树ZF-HD基因家族成员,利用在线生物信息学分析软件对其基因结构、启动子元件,以及编码蛋白的理化性质、氨基酸序列、结构特征、进化关系等进行预测分析,并基于转录组测序（RNA-Seq）数据分析其在不同组织及在干旱胁迫、盐胁迫和茉莉酸甲酯（MeJA）处理下的表达模式。【结果】从茶树基因组中鉴定出17个ZF-HD基因家族成员（CsZHD1~CsZHD17）,其编码区（CDS）序列长度为369~2187 bp,编码122~728个氨基酸残基,蛋白分子量13.51~80.42 kD,理论等电点为6.09~9.19,均属于亲水性不稳定蛋白,蛋白二级结构均由β-转角、延伸链、α-螺旋和无规则卷曲构成。除CsZHD2、CsZHD5和CsZHD7蛋白定位于细胞质,CsZHD9蛋白定位于细胞外,其他13个蛋白均定位于细胞核。仅CsZHD2、CsZHD7、CsZHD9、CsZHD10和CsZHD12基因含外显子和内含子,其他12个CsZHDs基因均无内含子。除CsZHD7蛋白具有2个ZF-HD_dimer结构域外,其他16个蛋白均具有1个ZF-HD_dimer结构域。CsZHD1~CsZHD17蛋白具有2~5个保守基序（motif）,其中motif 1和motif 3为共有的保守基序。茶树ZF-HD家族蛋白可分为5个亚族（ZHDⅠ、ZHDⅡ、ZHDⅢ、ZHDⅣ和MIF）,较拟南芥少了ZHDⅤ亚族。除CsZHD2、CsZHD12和CsZHD17基因在8个组织中不表达或表达量极低外,其他14个CsZHDs基因在8个组织中呈差异性表达;除CsZHD2、CsZHD12和CsZHD17基因外,其他15个CsZHDs基因在顶芽或花中表达量较高。在干旱胁迫和盐胁迫处理下,除CsZHD1、CsZHD2、CsZHD5、CsZHD12、CsZHD14和CsZHD17基因的表达量始终处于较低水平外,其他CsZHDs基因均呈不同的变化趋势;在MeJA处理下,除CsZHD2、CsZHD5和CsZHD12基因表达量极低外,多数CsZHDs基因呈下调表达的趋势。【结论】从茶树基因组中鉴定出17个ZH-HD基因家族成员,其编码蛋白具有保守的锌指结构域和同源异型盒结构域;茶树与拟南芥ZH-HD基因家族相比缺少ZHDⅤ亚族;CsZHDs基因的表达具有组织特异性,且大多数成员的表达受非生物胁迫和MeJA的影响。     

玉米转录因子ZmbZIP77基因克隆及表达分析

碱性亮氨酸拉链蛋白(Basic Leucine Zipper,bZIP)是广泛存在于真核生物中的一类转录因子,广泛参与植物多种生理进程.从玉米(Zea mays L.)中克隆出1个碱性亮氨酸拉链蛋白家族成员ZmbZIP 77,对该基因进行生物信息学分析,并采用实时荧光定量PCR技术分析了其在玉米不同组织部位及不同胁迫处理下的表达谱.结果表明,ZmbZIP77与其他bZIP家族成员一样,均含有1个保守的GBF1(Plant G-box binding factor 1)结构域,暗示它和该家族中其他成员具有相同或相似的功能;根据bZIP序列的同源性,植物bZIP序列可以分为3类,ZmbZIP77位于第Ⅱ类,与小麦Wlip9a和拟南芥AtbZIP10亲缘关系最近;玉米ZmbZIP77基因在根中表达量最高,胚芽鞘次之,在叶中表达量最低;高温、干旱和脱落酸均能显著诱导玉米根中ZmbZIP 77的表达,提示Zm-bZIP 77可能在玉米组织发育及逆境胁迫响应中起到重要作用.     

玉米ZmEXO70s基因家族鉴定及苗期耐热性关联分析

【目的】鉴定玉米EXO70s（ZmEXO70s）基因家族成员,分析其遗传变异与玉米耐热性的关联性,为揭示其在玉米耐热方面的功能和分子机制打下基础。【方法】以拟南芥EXO70s（AtEXO70s）基因家族成员序列为参考,利用MaizeGDB和NCBI数据库从玉米全基因组中鉴定出ZmEXO70s基因家族成员,对其进行基因结构及系统进化分析,根据其在系统发育进化树上的位置对其进行系统命名,并基于转录组测序（RNA-Seq）数据分析ZmEXO70s基因家族成员在不同组织及非生物胁迫下的表达模式。最后,鉴定玉米自交群体AM368中ZmEXO70s基因SNP位点,结合玉米AM368群体苗期耐热存活率进行基因家族关联分析。【结果】从玉米全基因组序列中共鉴定出36个ZmEXO70s基因,分布在10条染色体上,长度为228～3726 bp,编码75～1241个氨基酸残基,蛋白分子量为8.6～136.6 kD,理论等电点（pI）介于4.5～9.9,大部分蛋白pI小于7.0。拟南芥、水稻和玉米的EXO70s蛋白被分为9组（A组～I组）,其中ZmEXO70s蛋白在各组中均有分布。大多数ZmEXO70s基因不含内含子,只有少数基因含有内含子,且内含子数量不同。36个ZmEXO70s蛋白共含8种基序（Motif 1～Motif 8）,主要集中在C端,说明这些蛋白端序列较保守,且Motif 1～Motif 8组成Exo70保守结构域。ZmEXO70s基因在不同组织中的表达水平均存在明显差异,其中ZmEXO70B3基因在雄穗、花药、叶片和花丝中高表达,ZmEXO70C1a基因只在雄穗和花药中高表达,ZmEXO70D2b基因在胚乳和萌发种子中低表达,ZmEXO70G1c基因在花丝和萌发种子中低表达。有23个ZmEXO70s基因至少可响应一种非生物胁迫,说明这些ZmEXO70s基因参与非生物胁迫响应过程。ZmEXO70D2a和ZmEXO70E1基因的遗传变异与玉米苗期耐热性具有显著相关性（P≤0.01）。【结论】ZmEXO70s基因家族成员在系统发育进化上较保守,其基因组复制事件可能发生在禾本科植物分化后,且大多数ZmEXO70s基因被保留,仅部分基因被丢失。ZmEXO70s基因可能在非生物胁迫响应过程中发挥重要作用,ZmEXO70D2a和ZmEXO70E1基因可作为调控玉米苗期耐热性重要候选基因。     

小麦TaBTLs基因家族全基因组鉴定与表达模式分析

【目的】鉴定和分析小麦TaBTLs基因家族成员。【方法】通过生物信息学方法，分析小麦TaBTLs家族成员的数量、生化性质、基因结构、共线性关系、系统进化关系和启动子顺式作用元件；利用小麦RNA-Seq数据，分析TaBTL基因在不同组织器官和发育时期的表达模式以及对非生物胁迫的响应水平。【结果】小麦共有35个TaBTLs基因家族成员，所有成员均含RING-H2和BZF保守结构域，相对分子量介于29.64～45.16 ku之间。TaBTLs基因分布在除2A、2B和2D染色体以外的所有染色体上，且各成员之间存在大量重复事件。不同物种间BTLs基因分为7个亚族。TaBTLs基因在启动子区域存在大量的植物生长发育和逆境响应顺式作用元件；该基因家族成员在不同的组织和器官中广泛表达，并且受盐和干旱胁迫诱导。【结论】研究结果为阐明TaBTLs基因家族的进化关系和进一步研究其家族成员的功能提供理论依据和前期工作基础。     

蓖麻RcFEN1基因克隆及非生物胁迫下表达模式分析

【目的】克隆蓖麻结构特异性核酸酶1基因(RcFEN1),并检测其组织表达特性及不同非生物胁迫下的时空表达模式,为揭示RcFEN1基因在蓖麻生长发育和响应低温胁迫等非生物胁迫的功能提供理论参考。【方法】根据蓖麻低温胁迫转录组注释信息,通过RT-PCR对RcFEN1基因编码区(CDS)进行克隆,利用生物信息学软件进行序列特征分析,并利用实时荧光定量PCR(q RT-PCR)检测其组织表达特性及不同非生物胁迫下的时空表达模式。【结果】克隆获得的RcFEN1基因(GenBank登录号OP205366)包含1个1154 bp的开放阅读框(ORF),编码383个氨基酸残基,含有RAD2/XPG核酸酶蛋白家族典型的结构域XPG-N和XPG-I,属于RAD2/XPG蛋白家族成员。RcFEN1蛋白含多个生物活性位点及功能区,但无跨膜区域和信号肽结构,是一个非分泌型亲水蛋白,定位在细胞核。RcFEN1蛋白与特洛花TsFEN1蛋白的氨基酸序列一致性最高,与大叶藻ZmFEN1蛋白序列的一致性最低,分别为90.28%和84.32%。RcFEN1蛋白与一串红和硬皮地星的FEN1蛋白亲缘关系较近。RcFEN1基因在子叶中的相对表达量显著高于根、茎和真叶(P<0.05,下同),在真叶中的相对表达量最低。RcFEN1基因在4种非生物胁迫下均有表达,但表达模式存在明显差异,其中,RcFEN1基因对干旱胁迫响应最敏感,处理2和4 h时相对表达量显著高于对照组(0 h)及其他处理时间;RcFEN1基因在盐胁迫和ABA胁迫下呈现相同的响应表达模式,均在处理4 h时开始表达,随后相对表达量降低;在低温胁迫下,RcFEN1基因在低温胁迫下延迟表达,且处理12 h后才被激活表达,且相对表达量持续增加。【结论】从蓖麻中克隆获得RcFEN1基因属于RAD2/XPG家族成员,是蓖麻低温胁迫下的差异表达基因,参与蓖麻低温胁迫的响应调控。     

普通小麦转录因子基因TaWRKY35的克隆及功能分析

【目的】非生物逆境是制约小麦生产的主要因素。WRKY转录因子是植物特有的转录因子家族,参与对多种非生物逆境胁迫的应答。普通小麦基因组中至少含有200个WRKY,目前仅对少数成员的功能进行了鉴定。通过克隆小麦WRKY转录因子TaWRKY35并揭示其功能,为改良小麦的抗逆性提供基因资源。【方法】利用小麦全长cDNA数据信息,从强抗旱小麦品种旱选10号中克隆逆境胁迫应答转录因子基因TaWRKY35;采用实时定量PCR技术,分析目标基因在不同发育时期、不同组织器官中的表达水平,以及在PEG-6000、NaCl、ABA和低温等逆境胁迫下的表达模式;构建目标基因与GFP融合的表达载体,转化小麦原生质体,以确定目标蛋白在细胞中的作用位置;构建TaWRKY35过表达载体,通过农杆菌介导的方法转化拟南芥,揭示目标基因的功能。【结果】TaWRKY35的开放阅读框全长1 134 bp,编码377个氨基酸的蛋白。TaWRKY35的N端有一个典型的WRKY结构域,C端有一个C₂HC型的锌指结构域,属于WRKY家族第Ⅲ亚家族。测序发现TaWRKY35编码区序列非常保守,在32份多态性较高的六倍体小麦材料中未检测到DNA多态性。TaWRKY35在小麦的不同发育时期、不同组织中均有表达,其中在幼苗根基部表达量最高,约为苗期叶片中表达量的26倍;其次为幼芽根基,约为苗期叶片中表达量的21倍;接下来依次为孕穗期茎节、孕穗期穗、芽期根、孕穗期根、苗期根、胚芽、孕穗期叶片,而在幼苗叶片中表达水平最低;同时在ABA、NaCl、PEG和低温胁迫条件下,小麦幼苗TaWRKY35的表达量均显著上升,表明TaWRKY35参与对4种逆境胁迫的应答,但在不同胁迫条件下表达模式差异显著,其中对盐胁迫最敏感。亚细胞定位发现,TaWRKY35特异定位在细胞核上,是典型的核定位蛋白。在高盐胁迫条件下,TaWRKY35过表达转基因拟南芥子叶白化率显著低于对照,TaWRKY35过表达转基因拟南芥的细胞膜稳定性和幼苗存活率均显著高于野生型对照,表明TaWRKY35能增强植物的耐盐性。【结论】TaWRKY35编码蛋白含有WRKY和C₂HC结构域,为WRKY家族第Ⅲ亚家族成员。TaWRKY35在小麦发育的不同时期、不同组织中均有表达,且受ABA、PEG、NaCl及低温等逆境胁迫诱导。过表达TaWRKY35能显著提高转基因拟南芥的耐盐性。     

白桦BpERF1和BpERF2基因的克隆与表达模式分析

【目的】对白桦Bp ERF1和Bp ERF2基因进行生物信息学分析及非生物胁迫下的表达模式分析,为白桦抗逆育种提供重要基因源和理论依据。【方法】从白桦转录组数据库中找到2条ERF基因,通过生物学网站及软件对其进行生物信息学分析;通过RT-PCR分析盐和干旱胁迫下2个ERF基因的表达模式。【结果】Bp ERF1蛋白含有2个核定位信号,Bp ERF2蛋白含有1个核定位信号;两个ERF蛋白的结构域均含有3个反向平行的β折叠和一个α螺旋;Bp ERF1和Bp ERF2都属于ERF亚家族成员;在Na Cl胁迫下,Bp ERF1基因在根和叶中的相对表达量与对照相比呈上调表达的趋势;Bp ERF2基因在根茎叶中均呈下调表达的趋势;在PEG胁迫下,两个基因的表达模式与Na Cl相似。【结论】生物信息学分析结果推测Bp ERF1和Bp ERF2基因可能参与植物对生物和非生物胁迫的应答;2个ERF基因均可以不同程度地被高盐和干旱胁迫所诱导,但不同ERF基因在胁迫下的表达模式存在差异,ERF基因的表达具有组织特异性,说明2个ERF基因可能通过不同的调控机制来应答逆境胁迫。     

小麦蛋白磷酸酶2A调节亚基基因TaBβ-1的克隆及其在非生物胁迫下的表达特性

【目的】克隆小麦蛋白磷酸酶2A(PP2A)调节亚基(PR55)基因TaBβ-1,分析其在非生物胁迫下的表达特性,为小麦抗逆育种提供候选基因。【方法】以小麦品种旱选10号为材料,通过电子克隆和RT-PCR获得TaBβ-1的全长cDNA序列,采用生物信息学软件分析TaBβ-1及其编码蛋白TaBβ-1的序列特征,预测其功能,利用实时荧光定量PCR(real-time quantitative PCR)技术分析该基因在小麦孕穗期不同组织中、不同生育时期新叶中的表达情况,以及在PEG、NaCl、低温及外源激素ABA等非生物胁迫下的表达模式,检测不同水分条件下成株期转基因拟南芥的叶片细胞膜稳定性。【结果】获得TaBβ-1的全长cDNA序列1931bp,其开放阅读框(ORF)为1539bp。该基因编码512个氨基酸,预测TaBβ-1蛋白分子量为57.1kD,等电点为5.87,含有PP2A调节亚基的1个CDC55保守结构域、1个alpha/beta结构域、2个PR55保守结构域和6个WD重复子。TaBβ-1在小麦孕穗期的根、穗、叶中均有表达,表达量依次为根穗叶;在不同生育时期的新叶中,苗期叶片的表达量最高。TaBβ-1的表达明显受PEG、NaCl、低温以及外源ABA的诱导。【结论】小麦蛋白磷酸酶2A调节亚基家族基因TaBβ-1在小麦苗期叶片中的表达量明显高于根、穗以及其它时期的叶片;TaBβ-1参与对高渗、高盐、低温等多种胁迫以及ABA处理的应答反应,但表达模式不同;在水分胁迫条件下,转基因拟南芥比野生型具有较高的细胞膜稳定性。     

大白菜ACA基因家族的全基因组鉴定与表达分析

【目的】通过对大白菜ACA（Ca²⁺-ATPase）基因家族鉴定与表达分析,研究其家族基因间的共性与特性,为进一步揭示ACA家族进化关系提供数据支撑,为深入解析BraACAs在低温胁迫、盐胁迫以及自交不亲和方面的功能研究奠定基础。【方法】根据已报道的拟南芥ACA基因家族,同源比对出大白菜ACA基因家族,利用在线软件Expasy预测其分子量、理论等电点等理化性质;采用MEGA 5.0软件构建系统进化树;运用在线软件GSDS 2.0绘制基因结构图谱;TBtools对其染色体定位;McscanX软件进行拟南芥与大白菜ACA家族基因共线性分析;利用在线软件PlantCARE预测大白菜ACA基因家族启动子元件;通过在线工具Pfam和MEME进行蛋白保守结构域分析;利用qRT-PCR技术检测BraACAs在不同组织、非生物胁迫和自交异交授粉后的表达量。【结果】大白菜ACA基因家族有18个基因成员,分布在大白菜10条染色体上;根据进化树关系分成4组,分别包含3、4、4和7个成员;蛋白结构域分析显示,有13个成员包含N端自抑结构域。qRT-PCR结果表明,BraACAs主要在花与果荚中高表达;低温胁迫下,Bra002762与Bra035649表达量总体上调;盐胁迫下,Bra031701表达量显著上调;自交和异交授粉中,Bra003276和Bra024117差异性表达。对Bra002762、Bra035649、Bra031701、Bra003276和Bra024117亚细胞定位分析,发现这些基因均定位在细胞质膜上。【结论】大白菜ACA家族基因蛋白结构均含有4个ACA基因特有的高度保守结构域。该家族在大白菜不同组织中表达模式不同,5个ACA家族基因成员编码蛋白定位于细胞膜上,其中Bra002762、Bra035649、Bra031701与低温和盐胁迫响应有关;Bra003276和Bra024117与自交不亲和性相关。     

睡莲转录因子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条染色体上,结构进化保守,组织表达模式多样。     

苹果LRR-RLK基因家族鉴定和表达分析

【目的】从苹果全基因组中鉴定LRR-RLK家族蛋白成员,并进行生物信息学和表达模式分析,为研究苹果LRR-RLK的潜在功能提供理论基础。【方法】利用BLASTp基于GDR数据库鉴定苹果LRR-RLK家族成员,通过ExPASy Proteomics Server、Cell-PLoc、CD-Search Tool、MEGAX、MG2C等软件分析LRR-RLK蛋白序列基本信息、亚细胞定位情况、结构域组成、系统进化关系以及染色体定位情况。利用实时荧光定量PCR技术检测苹果12个LRR-RLK的组织表达和诱导表达特性。【结果】苹果LRR-RLK基因家族包含378个成员,这些LRR-RLK蛋白包括318—1 827个不等的氨基酸残基,等电点分布在5.16—9.75。亚细胞定位结果显示LRR-RLK蛋白均定位在细胞膜。系统进化分析可将其分为15类,各亚家族基因数量分布在1—111。染色体定位结果显示,LRR-RLK在苹果17条染色体上均有分布,其中第7条染色体数量最多,为40个。LRR-RLK家族基因具有2个特定的保守结构域,分别是LRR结构和蛋白激酶结构。蛋白二级结构以无规则卷曲为主,其次是α-螺旋,β-转角所占比例最小。通过定量检测发现筛选的12个家族成员在各组织中均有表达（除MD00G1105400外）,且多数基因在茎中表达水平相对较高。低温条件下,7个基因上调表达,其中MD09G1153800上调最明显,最高为对照的6.8倍,而MD06G1170200和MD05G1061600均下调表达;在干旱条件下,8个基因上调表达,其中MD00G1105400上调最明显,最高为对照的9.6倍;在盐胁迫条件下,MD04G1150400、MD13G1108000和MD02G1071800始终处于上调表达状态,其中MD02G1071800上调最明显,最高为对照的14.9倍。苹果新梢接种轮纹病菌后,12个LRR-RLK家族基因表达基本上呈先上升后下降的趋势。并且在‘望山红’中,1 d时表达水平较高,而在‘鸡冠’中,3 d时表达水平较高。MD09G1153800和 MD05G1065800在‘鸡冠’响应轮纹病菌侵染过程中显著上调表达,而在‘望山红’中无响应,可作为进一步开展抗病研究和功能分析的候选基因。【结论】苹果LRR-RLK基因家族包含378个成员,进化上可分为15组,在17条染色体上均有分布,多数基因具有在茎中高表达的组织表达特征,多数基因受逆境和轮纹病菌调控。     

茶树Whirly基因家族的鉴定与表达分析

【目的】对茶树Whirly基因家族成员进行鉴定及表达分析,为进一步探究Whirly基因家族在茶树生长发育和逆境响应的作用机理提供参考。【方法】以茶树品种北茶36、龙井43和大面白的叶片为材料,基于Whirly基因家族的保守结构特征,从茶树基因组鉴定CsWhys基因家族成员,利用生物信息学软件对其进行可视化分析,并对组织特异性和时空表达模式分析。【结果】从舒茶早茶树基因组中鉴定获得2个CsWhy基因,命名为CsWhy1(XP_028085682.1)和CsWhy2(XP_028102803.1)。系统进化树分析结果显示CsWhy1和CsWhy2蛋白分别与已知的定位于叶绿体和线粒体的Whirly蛋白聚为一类。CsWhys基因启动子区域富含激素、胁迫和MYB转录因子结合位点等顺式作用元件,种类呈多样性。组织表达特性分析显示,CsWhy1在茶树叶片和芽中具有较高的表达丰度;CsWhy2则在茶树不同组织中均具有较高的表达丰度。时空表达模式分析结果发现,CsWhys的表达在不同时间点受低温、干旱、ABA、GA和H₂O₂处理诱导上调。冷驯化条件下,CsWhy1在抗寒品种龙井43和不抗寒品种大面白中呈现相反的表达模式,而CsWhy2在2个品种中均下调表达。【结论】茶树基因组中2个CsWhys基因具有Whirly基因家族的典型特征,茶树CsWhys基因的表达具有明显的组织特异性,推测CsWhys基因家族在茶树芽叶生长和抗逆过程中起发挥了重要调控作用。     

芒果VOZ转录因子基因家族成员鉴定及生物信息学分析

【目的】对芒果维管束锌指蛋白（VOZ）基因家族成员进行鉴定并对其生物信息学进行分析,为深入探究芒果VOZ转录因子在免疫反应中的生物学功能提供理论参考。【方法】基于芒果全基因组数据,利用生物信息学方法鉴定芒果VOZ转录因子基因家族成员,并分析其理化性质、保守基序及系统发育进化。通过实时荧光定量PCR检测胶孢炭疽菌（Colletotrichum gloeosporioides）与细菌性黑斑病菌（Xanthomonas citri pv.mangiferaeindicae）侵染下的相对表达量。【结果】从芒果全基因组中间鉴定出4个VOZ转录因子家族成员,分别为MiVOZ1、MiVOZ2、MiVOZ3和MiVOZ4基因,开放阅读框（ORF）长度为1290~1482 bp,编码氨基酸数量为429~493,蛋白分子量为47.26~54.96 kD,等电点（pI）为5.47~6.00,亲/疏水性指数为-0.663~-0.552,不稳定指数为37.66~50.86,二级结构均以无规则卷曲和α-螺旋为主要元件。芒果和其他9个物种的49个VOZ蛋白聚为十个分支,在ClassⅠ中3个MiVOZs蛋白（MiVOZ2、MiVOZ3和MiVOZ4）与苹果、木薯和毛果杨聚类在一起,在ClassⅤ中MiVOZ1与拟南芥、木薯和毛果杨聚类在一起。在胶孢炭疽菌侵染下,仅MiVOZ1和MiVOZ2基因的相对表达量较对照显著升高（P<0.05,下同）;在细菌性黑斑病菌侵染下,仅MiVOZ2基因的相对表达量较对照显著升高。【结论】芒果VOZ转录因子在抵御不同病原菌侵染的生物学功能方面存在差异,其中MiVOZ2基因在抵御胶孢炭疽病和细菌性黑斑病侵染的免疫反应中具有相似生物学功能,属于正调节因子。     

小麦PIN基因家族的鉴定及表达分析

【目的】PIN（puroindoline）是植物所特有的一类蛋白家族,对控制小麦籽粒硬度有重要功能。分析小麦PIN家族成员在全基因组的分布、结构及进化,研究其在不同组织的表达特异性以及在不同硬度种子的表达模式,为阐明小麦PIN基因家族的生物学功能奠定基础。【方法】根据已报道的小麦PIN基因和大麦HIN基因,利用BLASTP和HMM方法,在最新发布的小麦中国春参考序列中鉴定小麦PIN基因家族成员。利用UniProt、URGI、PFAM、CDD、expVIP等数据库,采用Clustal X、MEGA 7.0、ExPASy、MEME、GSDS、TB tools、GraphPad Prism5等软件进行生物信息学分析。采用qRT-PCR方法检测TaPIN基因家族在不同籽粒硬度小麦样品中的表达情况。【结果】共鉴定出19个小麦PIN基因,集中成簇分布于第1、5和7染色体同源群,编码148—327个氨基酸,编码蛋白相对分子量为16.39—37.19 kD,等电点为6.35—9.34。通过结构域和系统发育分析,可将19个TaPIN基因分为A和B两大类。大部分TaPINs基因仅有1个外显子,没有内含子,顺式作用元件分析发现其上游序列包含大量抗逆和种子发育相关的调控元件。转录组分析表明该基因家族在小麦籽粒中相对表达量很高,而在根茎叶等其他组织几乎不表达。实时荧光定量PCR表明,各基因间相对表达量差异显著,TaPIN9和TaPIN10表达量较高。随着小麦籽粒硬度降低,TaPIN9和TaPIN10表达上调,且表达比例增加,而TaPIN16和TaPIN6则呈现相反的趋势。【结论】小麦籽粒硬度的调节以Pina和Pinb为主,该基因家族其他成员也具有相同结构域,推测也具有相似功能,但受表达量低的限制,对籽粒硬度影响较小。从该基因的进化关系看,粗山羊草与小麦亲缘关系最近,其次是燕麦、黑麦和大麦。     

大麦TIFY基因家族成员鉴定及表达分析

【目的】对大麦TIFY基因家族成员进行鉴定及表达分析,为进一步探究TIFY基因家族在大麦生长发育与胁迫响应中的作用机理打下基础。【方法】基于TIFY家族蛋白的保守域特征,利用HMMER从大麦中鉴定TIFY基因家族成员,利用采用生物信息学软件对其理化性质、保守基序、特征结构域、顺式作用元件、基因结构、系统进化及表达模式进行预测分析。【结果】从大麦中鉴定出15个HvTIFYs基因（HvTIFY1～HvTIFY15）,分布于5条染色体上,且大多数基因在染色体上成簇分布。15个HvTIFYs蛋白均具有TIFY家族蛋白的特征结构域（TIFY）,根据所含保守结构域的不同,可分为ZML（4个）和JAZ亚族（11个）,且亲水性蛋白（14个）和偏碱性蛋白（11个）居多,但均定位于细胞核;二级结构相似度较高,均由α-螺旋、β-转角和无规则卷曲组成,除HvTIFY7蛋白外,其余蛋白二级结构所占比排序:无规则卷曲>α-螺旋>β-转角。HvTIFYs基因结构存在明显差异,其中,JAZ亚族11个基因的内含子数为0~6; ZML亚族4个基因的内含子数为6～7个,系统发育进化树上相邻分支的基因具有较相似的基因结构。HvTIFYs基因启动子区域富含光、激素和胁迫等顺式作用元件,种类及分布均呈多样性。5个物种的79条TIFY蛋白分为4个组,恰好与TIFY家族的4个亚族对应,其中,ZML、TIFY和JAZ亚族包含单、双子叶植物的TIFY蛋白,而PPD亚族仅含有双子叶植物的TIFY蛋白。15个HvTIFYs基因在不同组织器官中的表达量存在明显差异,其中HvTIFY1、HvTIFY2和HvTIFY8基因在8个组织中的表达量均较高,HvTIFY10和HvTIFY15基因表达量中等,HvTIFY6基因表达量较低; HvTIFY11基因不表达。15个基因在根的不同组织中对盐胁迫的敏感程度不同。【结论】从大麦中鉴定出的15个HvTIFYs基因存在一定的功能分化,具有明显的组织和时空特异性,推测其在大麦逆境响应和激素调节中具有重要调控作用。     

油棕SAUR基因家族的全基因组鉴定及生物信息学分析

【目的】挖掘并鉴定油棕生长素上调小RNA(SAUR)基因家族成员,并对其进行生物信息学分析,为油棕生长素信号转导机制的研究提供参考。【方法】根据拟南芥SAUR基因序列,在油棕全基因组数据库中同源序列搜索,并利用CDD和Pfam数据库进行SAUR蛋白结构域分析,剔除无保守结构域的序列,最后利用生物信息学软件对油棕SAUR基因家族成员进行可视化分析。【结果】从油棕全基因组中共鉴定出40个SAUR基因家族成员(EgSAUR1~EgSAUR40),有31个基因在油棕14条染色体上呈不均匀分布,其中4号染色体上分布的基因最多为7个;除EgSAUR14和EgSAUR24基因含有2个外显子外,其余38个基因均仅含1个外显子,无内含子。40个EgSAURs蛋白亚细胞定位于细胞核,其二级结构均由α-螺旋、β-转角、无规卷曲和延伸链组成,以无规卷曲和α-螺旋所占比例较高,其中有11个EgSAURs蛋白呈酸性,29个EgSAURs蛋白呈碱性。40个油棕SAUR蛋白被划分为四大类,其中第Ⅰ(除EgSAUR27蛋白外)、Ⅱ(除EgSAUR3和EgSAUR22外)、Ⅲ、Ⅳ类蛋白均有Auxin_inducible保守结构域。40个EgSAURs蛋白共含10种保守基序(Motif),其中Motif 2存在于所有基因成员中。在花中特异表达的EgSAURs基因(EgSAUR4、EgSAUR9、EgSAUR10、EgSAUR24、EgSAUR29、EgSAUR30、EgSAUR32和EgSAUR33)数量最多,其次为根、茎和叶,在中果皮中特异表达的基因数最少。【结论】基因重复和组织表达特异性是油棕SAUR家族基因的两大特点,推测该家族基因对油棕花的生长发育及授粉受精过程发挥调控作用。