潘那利番茄GRF基因家族全基因组鉴定分析

【目的】分析研究潘那利番茄GRF基因家族全基因组鉴定分析。【方法】研究基于生物信息学方法对潘那利GRF基因家族成员进行全基因组鉴定和分析,并对其起源及进化进行了追溯。【结果】潘那利番茄中共鉴定了10个SpGRF成员,不均匀分布在7条染色体上,预测了SpGRFs蛋白的分子量、等电点、亲水性总均值等理化性质。SpGRFs基因可分为6个亚家族。所有SpGRF基因N端都含有1个QLQ和1个WRC结构域,C端则为多种保守基序。共线性结果显示基因组内有3对6个旁系同源基因,全部为片段复制,SpGRFs受自然选择压力下,有共同起源祖先,与拟南芥亲缘关系更近。【结论】鉴定了潘那利番茄中GRF基因家族的基本信息。     

黄瓜ERF基因家族鉴定及其在雌花芽分化中的表达分析

【目的】通过以黄瓜9930_V2版本基因组为参照进行生物信息学分析,对ERF基因家族在基因组中的数量、结构以及表达特征进行分析,为研究ERF转录因子在黄瓜雌花分化与发育中的作用提供数据支持。【方法】根据已报道的拟南芥ERF,利用黄瓜基因组数据库中9930_V2版本基因组进行BLAST比对,通过MEGA、MEME、TBtools、ExPASy等工具进行生物信息学分析。采用qRT-PCR方法检测不同性型黄瓜材料、雌花发育初期不同阶段中ERF基因家族成员的表达水平。采用酵母单杂交方法验证家族成员与乙烯响应元件GCC-box的互作。【结果】从黄瓜材料9930基因组中鉴定得到138个ERF基因家族成员,共分为10个亚族,编码氨基酸长度介于126—745。按照基因家族成员在染色体上的位置分布,将其命名为CsERF1-CsERF138。多序列比对和motif分析结果表明,黄瓜ERF基因家族均具有AP2/ERF结构域,其中4个成员具有B3结构域。表达分析结果显示,在不同性型材料中共有19个ERF家族成员差异表达,其中9个在FFMMAA基因型中高表达,10个在ffMMAA基因型中高表达。通过雌花芽发育初期ERF家族成员的表达趋势分析,发现31个ERF随子房发育表达上调,30个表达下调。初步证明CsERF9和CsERF31具有结合GCC-box元件的功能。【结论】从黄瓜基因组中鉴定出138个ERF基因家族成员,均拥有1个或多个AP2/ERF结构域;其中部分成员在不同性型材料中差异表达,并可能参与雌花分化初期的基因表达调控;部分成员具有结合保守元件GCC-box调控下游基因表达的功能。     

水稻WOX家族基因鉴定及其种子在逆境萌发中的表达分析

【目的】开展水稻WOX基因家族成员的鉴定及在种子逆境萌发中的表达分析,可为探究WOX家族基因的功能、作用机制及水稻遗传改良提供理论参考。【方法】对水稻WOX基因家族成员进行鉴定,利用生物信息学软件对其系统进化关系、理化性质、基因结构、保守结构域、保守基序、顺式作用元件和染色体定位等进行分析,并利用实时荧光定量PCR（q...     

大白菜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与自交不亲和性相关。     

苦荞VQ基因家族的全基因组鉴定及其在叶斑病原与激素处理下的表达谱分析

【目的】VQ基因家族在植物生长、发育以及对生物或非生物胁迫反应中发挥重要功能。在全基因组尺度上,全面鉴定苦荞（Fagopyrum tataricum L. Gaertn.）VQ（FtVQ）基因家族,分析其在苦荞叶斑病原——互格链格孢（Alternaria alternata）和黑孢霉（Nigrospora osmanthi）侵染和防御相关激素——水杨酸（SA）、茉莉酸（JA）、乙烯（ET）处理下的表达模式,为深入解析苦荞VQ基因家族在植物抗病防御中的功能及机理奠定基础,同时为优良基因资源发掘及抗病品种改良提供线索。【方法】基于VQ保守结构域的隐马尔可夫文件（PF05678）,采用HMMER 3.0对苦荞平苦一号基因组数据库进行比对搜索,鉴定VQ基因;通过DNAMAN、MapInspect、MEGA、MEME、OrthoFinder、PLACE等生物信息学工具分析基因结构、染色体分布、启动子顺式元件、蛋白质理化性质、蛋白质保守基序、蛋白质亚细胞定位和蛋白质系统进化关系;采用实时荧光定量PCR（qPCR）方法分析苦荞叶VQ基因在病原侵染或激素处理下的表达模式。【结果】从苦荞基因组中鉴定获得28个VQ基因,大小为566—1 454 bp,均无内含子,不均一地分布在8条染色体上。根据它们在染色体上的物理位置,命名为FtVQ1—FtVQ28。每一个FtVQ蛋白含有1个VQ基序——FxxxVQx（L/F/I/V/A/Y）TG（x代表任意氨基酸）。亚细胞定位预测表明,21个FtVQ蛋白定位在细胞核中,其余定位在叶绿体或细胞质中。根据蛋白质氨基酸序列与保守结构基序,FtVQ蛋白归类于5个亚家族,亚家族内基因结构和蛋白质基序相对保守。基因重复分析表明,苦荞基因组中有8对VQ旁系同源基因,均为大片段重复基因,提示大片段基因重复在FtVQ基因家族数量扩张中发挥主要作用;它们的非同义突变和同义突变的比值（Ka/Ks）均小于1,提示重复基因在进化中经历了纯化选择。启动子顺式元件预测表明,所有FtVQ基因启动子含有BIHD1OS、CGTCA、ERELEA4、W-box和类W-box等病原或SA、JA、ET反应元件,尤其在FtVQ10、FtVQ14、FtVQ15、FtVQ22、FtVQ23、FtVQ27的启动子区域密集程度更高。qPCR分析显示,在可检测的20个FtVQ基因中,有55%—70%的基因为病原或激素处理下的差异表达基因（DEGs）,其中72.7%—85.7%的DEGs的表达显著上调。【结论】苦荞基因组拥有28个VQ基因成员,部分VQ基因可能参与了苦荞对叶斑病原的抗性反应。     

葡萄SAP家族的鉴定与表达分析

【目的】通过生物信息学分析明确SAP（Stress Associated Protein）家族在葡萄基因组中的数量、结构,利用qRT-PCR技术分析SAP家族的生物节律和在激素、非生物胁迫处理下的表达模式,为进一步探究SAP家族在葡萄中的作用奠定基础。【方法】在前期从山葡萄‘双优’（Vitis amurensis cv. Shuangyou）和欧洲葡萄‘红地球’（Vitis vinifera cv. Red Globe）0℃寒胁迫下转录组数据库中筛选出表达量明显上调的基因SAP15基础上,利用NCBI和葡萄基因组数据库中的BLAST功能,根据SAP家族的保守结构域,鉴定葡萄基因组中的SAP。利用生物信息学软件DNAMAN5.0、MEME、GSDS2.0、ExPASy和MEGA6等对VvSAPs序列、基因结构、蛋白质结构、理化性质、染色体定位和进化树等进行分析。采用qRT-PCR方法检测VvSAP家族的生物节律和在激素、逆境处理下的表达特征。【结果】从葡萄（Vitis vinifera）全基因组中鉴定得到15个SAP家族成员,均含有AN1保守结构域,大多含有A20保守结构域。按照其保守结构域和在染色体上的位置,依次命名为VvSAP1—VvSAP15,可分为Class Ⅰ、Class Ⅱ和Class Ⅲ三类。理化性质分析表明,VvSAP家族编码氨基酸数目介于109—293,理论等电点分布在7.99—9.68。染色体定位分析发现,该基因家族的15个基因分布于葡萄的9条染色体上,其中第8条染色体上分布最多,有3个VvSAP。亚细胞定位分析表明,VvSAP家族主要在细胞核、叶绿体和细胞骨架中进行表达。VvSAP家族的二级结构主要以无规则卷曲和α-螺旋为主。基因结构分析表明,VvSAP1—VvSAP12的DNA序列中均无内含子,VvSAP13—VvSAP15的DNA序列中有1个内含子,基因结构高度保守。qRT-PCR分析结果表明,在VvSAP家族成员中,VvSAP1和VvSAP9在无处理和各处理下（激素和非生物胁迫）均表达极低或不表达,初步鉴定为假基因。除VvSAP10在400 mmol·L ^-1 NaCl处理下呈下调表达,其余基因均在50 μmol·L ^-1 ABA、100 μmol·L ^-1SA和400 mmol·L ^-1NaCl处理下呈上调表达,其中,VvSAP10—VvSAP14显著响应50 μmol·L ^-1 ABA胁迫,处理24 h后相对表达量分别上调为0 h的37.19、36.63、21.69、58.34和267.35倍;VvSAP8和VvSAP11在NaCl处理4 h后相对表达量最高,分别为0 h的13.16和12.42倍;在4℃低温胁迫下,相对表达量上调最高的是VvSAP15,处理后8 h相对表达量为0 h的35.90倍。 【结论】从葡萄基因组中共鉴定出15个SAP家族成员,初步鉴定得出VvSAP1和VvSAP9为假基因,15个基因分别分布于9条染色体上,并且结构进化高度保守。所有成员均相响应逆境,且均有昼夜节律变化,但该基因家族成员在不同逆境胁迫下的响应程度和在不同时间点的表达模式存在一定的差异。     

芥菜全基因组中Hsf基因家族鉴定与分析

【目的】通过鉴定与分析芥菜基因组中Hsf转录因子家族成员,为芥菜Hsf基因功能研究与遗传改良提高抗逆性提供理论依据。【方法】利用生物信息学方法分析芥菜Hsf热激转录因子家族成员的功能结构、保守基序、启动子顺式作用元件、系统进化、共线性及采用RNA-seq验证Hsf低温胁迫的基因表达。【结果】芥菜基因组中鉴定出71个Hsf基因家族成员,分布于18条染色体上,归类为3个亚家族,蛋白均含有DBD和HR-A/B结构域。BjuHsf启动子区域包含与逆境胁迫、激素、生长发育等相关顺式作用元件。系统进化与共线性分析表明,芥菜Hsf家族成员与大白菜Hsf家族成员具有更近的亲缘关系。在低温胁迫下,BjuHsf基因表达分析表明8个BjuHsf基因显著上调表达。【结论】这些显著差异表达BjuHsf基因与芥菜抗寒性极其相关,可作为芥菜耐寒性遗传改良的候选基因。     

谷子NADP-ME的鉴定及其对逆境胁迫的响应

【目的】鉴定谷子NADP-ME家族成员,研究不同成员对非生物逆境胁迫的响应,为揭示SiNADP-ME在谷子逆境应答信号途径中的作用奠定基础。【方法】 利用生物信息学方法鉴定谷子基因组中的NADP-ME家族成员。采用GSDS2.0、plantCARE、Clustalx、MEGA6.0等软件及网站ExPASy对鉴定成员蛋白和基因序列进行生物信息学分析。采用qRT-PCR方法检测SiNADP-ME在苗期不同逆境、不同生育期干旱胁迫及不同光照强度下的表达情况。【结果】 谷子NADP-ME家族由7个成员组成,它们在谷子的第2、3、5、7染色体上呈不均匀分布。保守功能域分析显示7个基因都含有NADP-ME特征保守功能域。序列比对发现谷子NADP-ME成员之间序列非常保守,相似性较高,7个谷子成员序列一致性为77.30%,而不同物种NADP-ME序列之间相似性为56.52%。序列分析显示SiNADP-ME1、SiNADP-ME4、SiNADP-ME5、SiNADP-ME6序列较长,分别编码576、639、652和636个氨基酸,而SiNADP-ME2、SiNADP-ME3、SiNADP-ME7序列较短,分别编码213、265和149个氨基酸。基因结构分析显示SiNADP-ME1有2个可变剪切,SiNADP-ME5有3个可变剪切,其他基因无可变剪切。SiNADP-ME1、SiNADP-ME2、SiNADP-ME3、SiNADP-ME7含内含子较少,而SiNADP-ME4、SiNADP-ME5、SiNADP-ME6含内含子较多。蛋白参数预测显示谷子NADP-ME成员间分子量跨度较大,在161.94—725.43 kD,等电点为5.32—8.05,不稳定指数为23.01—45.01,脂肪系数介于89.19—107.77,平均疏水指数介于-0.218—0.004。亚细胞定位预测显示SiNADP-ME成员主要被定位在叶绿体、线粒体和细胞质中。顺式元件分析显示SiNADP-ME成员启动子区域主要包括激素类应答、逆境应答、光应答以及其他类生长调控相关的顺式元件。聚类分析发现谷子SiNADP-ME基因在单、双子叶植物分离之前就已存在。不同物种同源基因对在进化树中广泛存在揭示它们在进化上可能存在共同祖先,也暗示它们在某些信号通路中可能具有相似的功能。苗期逆境表达分析表明所有谷子SiNADP-ME家族基因表达量在本文应用的4种逆境胁迫下都被明显诱导。SiNADP-ME1在ABA、低温、NaCl处理后被诱导的最高相对表达量分别为对照的460.53、411.50和15.24倍;SiNADP-ME6在ABA、低温、PEG、NaCl处理后被诱导的最高相对表达量分别为对照的211.13、15.21、772.41和643.99倍。进一步分析表明SiNADP-ME1和SiNADP-ME6在拔节期、抽穗期和灌浆期干旱胁迫下表达量上调。【结论】 从谷子基因组中鉴定了7个NADP-ME基因家族成员;7个成员间序列非常保守并且都含有NADP-ME基因典型特征结构域;7个谷子NADP-ME家族基因参与了植物非生物逆境应答,特别是SiNADP-ME1和SiNADP-ME6可能在ABA、盐、干旱、低温等逆境应答信号途径中起重要作用。     

油棕脱落酸受体PYL基因家族的全基因组鉴定及表达分析

【目的】鉴定油棕（Elaeis guineensis）脱落酸（ABA）受体PYR/PYL/RCARs （PYL）基因家族成员,分析其表达特性,为探究ABA信号通路在油棕果肉成熟过程中的功能研究提供理论依据。【方法】以拟南芥和水稻的PYL蛋白氨基酸序列作为参考序列,通过BLASTp比对及保守结构域预测分析从油棕基因组中鉴定出PYL基因家族成员,利用生物信息学软件对其染色体定位、基因结构、启动子顺式作用元件及编码蛋白的理化性质、保守结构域、进化关系进行分析,并采用实时荧光定量PCR对PYL家族基因在不同组织、不同发育期果实及外源ABA处理下的表达特性进行检测。【结果】从油棕基因组中共鉴定出12个油棕PYL基因家族成员（EgPYL1~EgPYL112）,分布在8条染色体和1个Scaffolds上,含有1~3个外显子,开放阅读框（ORF）为564~765 bp,编码187~254个氨基酸,蛋白分子量为20.95~28.33 kD,等电点（pI）为5.26~7.95,不稳定指数为32.67~52.87,脂溶指数为73.87~87.60,总平均亲水性为-0.68~-0.17。12个PYL家族蛋白均含有特征结构域PYR/PYL/RCAR,分为3个亚族。EgPYL1和EgPYL6基因具有共线性,EgPYL4、EgPYL5、EgPYL9和EgPYL11基因具有共线性。EgPYLs基因的启动子上含有大量植物激素响应元件、逆境胁迫响应元件和光响应元件。EgPYLs基因在根、茎尖、叶、花和果肉中均有表达,但表达量差异较明显。EgPYL7、EgPYL8和EgPYL9基因的表达量随果肉成熟度增加逐渐升高,在23周达峰值。11个Eg PYLs基因均受外源ABA诱导表达。【结论】大多数PYL基因家族成员参与油棕对ABA的响应,且部分成员（如EgPYL7、EgPYL8和EgPYL9）在油棕果实发育中发挥重要的调控作用。     

苦荞NAC基因家族的生物信息学分析

以苦荞基因组数据库为基础,利用比对程序,经鉴定和筛选,共获得72个苦荞NAC家族基因的序列,对其进行生物信息学分析。基因结构分析表明,除FtPinG0002967400.01.T01基因外,苦荞NAC家族其余71个基因均含有内含子。蛋白质理化性质预测分析表明,72条蛋白序列具有NAC结构域。72个NAC基因共有7个保守基序,在苦荞8条染色体上均有分布。有33个NAC基因在种子发育过程中差异表达。系统进化分析表明,苦荞NAC蛋白序列与拟南芥的可以分为9个亚家族。     

龙眼乙烯合成途径基因鉴定及响应ACC处理的分析

目的S-腺苷甲硫氨酸合成酶（SAMS）、1-氨基环丙烷-1-羧酸合成酶（ACS）和1-氨基环丙烷-1-羧酸氧化酶（ACO）是植物合成乙烯的3个关键酶,从全基因组水平鉴定龙眼SAMS、ACS和ACO（DlSAMS、DlACS和DlACO）基因家族,并进行生物信息学、体细胞胚早期不同阶段及不同浓度1-氨基环丙烷-1-羧酸（ACC）和不同处理时间下龙眼胚性愈伤组织的表达模式分析,为进一步研究和利用DlSAMS、DlACS和DlACO基因家族奠定基础。方法 从拟南芥数据库（TAIR）下载SAMS、ACS和ACO家族蛋白质序列作为参考序列,TBtools、NCBI Blast等工具搜索龙眼基因组数据库,鉴定DlSAMS、DlACS和DlACO基因家族。使用ExPASy、PrediSi、TMHMM Server 2.0、NetPhos 3.1 Server、Plant-PLoc、MEME、PlantCARE、STRING、TBtools等软件和在线工具预测DlSAMS、DlACS和DlACO家族成员蛋白质基本理化性质、保守基序（Motif）与互作关系,定位染色体,分析基因共线性与选择压力、基因结构及启动子顺式作用元件;Clustal W和MEGA-X软件对龙眼、拟南芥、番茄、水稻和玉米的SAMS、ACS和ACO家族成员蛋白质序列分别进行多序列比对、构建系统进化树;根据龙眼转录组数据库的FPKM值,应用TBtools软件绘制DlSAMS、DlACS和DlACO家族成员在体细胞胚早期不同阶段的表达热图;采用qRT-PCR法分析不同浓度ACC和不同处理时间下龙眼胚性愈伤组织DlSAMS、DlACS和DlACO家族成员的表达情况。结果 龙眼第三代基因组中分别鉴定获得SAMS、ACS和ACO家族成员4个、11个和4个,在家族成员数量上与拟南芥、番茄等物种差异不大。DlSAMS家族成员最保守,所有成员的Motif相同;DlACS和DlACO家族成员也含有较多保守的Motif。3个家族共19个成员分布在11条染色体上,有6对共线性基因,与拟南芥有30对共线性基因,均受纯化选择作用。3个家族的成员均含有大量光、激素和逆境响应元件。家族成员内部的蛋白互作关系较弱,但家族成员之间互作关系非常紧密。多物种系统进化关系表明,SAMS、ACS和ACO家族均可分为3个亚家族,DlSAMS家族成员仅分布在SAMS-Ⅰ和SAMS-Ⅱ中,SAMS-Ⅲ可能为单子叶植物所特有;DlACS和DlACO家族成员在3个亚族中均有分布,且分布较为均匀;DlSAMS、DlACS和DlACO均与番茄和拟南芥的SAMS、ACS和ACO亲缘关系较近。对转录组FPKM值的分析表明,DlSAMS家族的所有成员以及DlACO4B在体细胞胚早期的3个阶段均高表达,可能在体细胞胚发生过程发挥着重要作用;DlACS1和DlACS6B在球形胚阶段高表达,可能与球形胚的形成密切相关。在胚性愈伤组织继代培养基中添加不同浓度的ACC能显著影响龙眼胚性愈伤组织的增殖量及DlSAMS、DlACS和DlACO家族成员的表达。培养20 d时,主要表现为基因的上调表达,且ACC浓度越高,上调越明显;在25—35 d则主要表现为基因的下调表达。但是在20 d时,0.01 mmol∙L^-1 ACC处理的基因主要表现为下调表达,这可能是导致其增殖量显著大于对照的原因。结论 龙眼第三代基因组中分别有SAMS、ACS和ACO家族成员4个、11个和4个,在进化过程中的保守性均较高,且包含大量激素和逆境响应元件;3个家族成员对龙眼体细胞胚胎的发生起重要调控作用,0.01 mmol∙L^-1 ACC处理可能通过调控DlSAMS、DlACS和DlACO家族成员的表达促进龙眼胚性愈伤组织增殖。     

黄瓜DIR家族基因的全基因组鉴定及其表达分析

【目的】基于黄瓜基因组信息和转录组测序大数据,利用生物信息学手段,对黄瓜中DIR基因家族进行鉴定,并分析其在不同组织器官和胁迫响应过程中的表达模式,为后续深入研究黄瓜DIR基因的生物学功能奠定重要基础。【方法】基于已报道的DIR基因HMM模型文件,利用HMMER软件包的hmmsearch程序从黄瓜蛋白数据库中筛选出可能的DIR基因ID,并利用在线工具Pfam和SMART进行验证,最终确定黄瓜DIR家族基因。利用ExPASy、TBtools、GSDS、MEME、MEGA、MCScanX和Circos等工具分析黄瓜DIR基因家族成员的理化特征、染色体定位、基因结构、系统进化树和共线性。基于黄瓜在不同组织和胁迫响应下的转录组测序大数据,利用黄瓜V3版本基因组信息进行转录组分析,检索黄瓜DIR基因在不同转录组测序分析中的表达情况,利用TBtools软件绘制表达热图,分析黄瓜DIR基因在不同组织和胁迫响应过程中的表达模式。【结果】在黄瓜中鉴定到23个DIR家族基因,分布于7条染色体上,编码氨基酸个数在78—684,分子量8.70—73.82 kD;系统进化分析将黄瓜DIR基因家族划分为3个亚族,...     

龙眼SPL基因家族全基因组鉴定及表达分析

【目的】鉴定龙眼SPL（DlSPL）基因家族所有成员,并对其成员进行表达分析,为探究SPL在龙眼生长发育中的功能研究提供参考。【方法】采用生物信息学方法对龙眼SPL基因家族成员进行鉴定,并对其基本理化性质、基因结构、系统进化关系、启动子顺式作用元件、FPKM值进行分析;利用qRT-PCR技术检测其在非胚性愈伤组织（nonembryonic callus,NEC）、胚性培养物及不同激素处理下的胚性愈伤组织（embryonic callus,EC）中的表达情况。【结果】DlSPL基因家族共有14个成员,其基因结构和蛋白结构均高度保守且只有一个SBP结构域,DlSPL启动子存在大量的光反应元件、组织特异性调控元件、激素应答元件、胁迫响应元件和植物生长发育有关的顺式调控元件。RNA-seq表达量分析（FPKM值）表明,11个DlSPL成员在不同光质处理的EC中表达,只有DlSPL8在蓝光和白光处理下的EC中呈下调表达趋势;11个DlSPL成员在激素2,4-D和KT中表达,其中DlSPL3和DlSPL13在2,4-D和激动素（kinetin,KT）共同处理下的表达量高于2,4-D和KT分开处理;有13个DlSPL成员在非胚性及胚性培养物中表达,其中7个成员在NEC阶段表达量最高。qRT-PCR结果表明,在不同胚性培养物中,DlSPL1和DlSPL14在EC阶段表达最高,DlSPL5、DlSPL7和DlSPL13在不完全胚性紧实结构（incomplete embryonic compact structure,ICpEC）阶段表达量最高;DlSPL1、DlSPL5、DlSPL7和DlSPL13成员响应脱落酸（abscisic acid,ABA）信号并显著下调;DlSPL5响应茉莉酸甲酯（methyl jasmonate,MeJA）信号且呈上调表达趋势,其他3个成员呈下调表达趋势。【结论】共鉴定出龙眼SPL基因家族成员14个,均含有一个高度保守的SBP结构域;DlSPL可能参与龙眼体胚的形态建成,并响应ABA和MeJA的应答。     

基于盐胁迫转录组信息的蚕豆F-box基因家族分析

【目的】通过生物信息学方法分析蚕豆F-box基因家族成员的分布、结构及进化,研究家族成员在不同处理时间条件下的表达模式及对盐胁迫的响应,为该类基因生物学功能和盐胁迫机制的研究提供参考。【方法】基于蚕豆盐胁迫转录组测序（RNA-seq）数据,利用NR、Swiss-prot和PFAM 3个数据库和NCBI网站,对蚕豆F-box基因进行筛选注释;利用Web Logo 3、Prot Comp 9.0、MEGA-X和MEME等软件进行保守结构域、亚细胞定位、系统进化树和Motif等生物信息学分析。基于盐胁迫转录组数据分析蚕豆（yz17134耐盐和yz17078不耐盐）F-box基因家族在盐胁迫下的差异表达模式,并采用实时荧光定量PCR技术（qRT-PCR）检测部分家族成员在16和24 h的表达情况。【结果】基于盐胁迫转录组测序数据,注释得到161个蚕豆F-box基因,均含有F-box保守结构域。根据C端结构域的不同,将其分成11个亚族：FBX、FBXFBA、FBXLRR、FBXPP2、FBXKelch、FBXTUB、FBXFBD、FBXDUF、FBXACTIN、FBXWD40和FBO。保守结构域分析表明,F-box保守基序中包含1个极度保守的色氨酸残基。比较分析蚕豆F-box家族和拟南芥F-box家族共同构建的进化树,发现同一C端结构域的基因大多聚集在一起。亚细胞定位预测结果显示,124个F-box基因定位于细胞外,37个定位于细胞核中。基因结构分析表明,蚕豆F-box家族基因的DNA序列中均无内含子,且均由UTR区和CDS区组成。基于盐胁迫转录组数据的F-box差异表达模式分析表明,蚕豆F-box基因在2个不同处理时间点上的表达各不相同,在盐处理16 h的表达较为明显。qRT-PCR分析结果表明,在F-box家族成员中,共存在5个差异基因。其中Vf056266.1、Vf062764.1和Vf024236.1在盐处理16 h的表达量均上调,Vf060904.1和Vf045761.1在盐处理16 h的表达量均下调。【结论】蚕豆F-box基因家族注释得到161个蚕豆F-box基因,分为11个亚族。其中5个重要的F-box基因在不同盐处理时间的表达量存在差异。     

番茄DIR基因家族鉴定及其对非生物胁迫响应的分析

【目的】鉴定番茄DIR基因家族所有成员,并对其基因结构、编码蛋白的理化性质、系统进化、染色体定位、共线性、启动子元件、表达模式、互作转录因子及内源竞争RNA预测等进行了解析,为探究DIR在番茄生长发育和逆境胁迫中的作用提供参考。【方法】采用生物信息学方法,基于全基因组数据对番茄DIR基因家族成员进行鉴定,运用Phytozome、MEME、PlantCARE、opsRNATarget和plantcircnet等在线网站获取染色体位置、保守基序、顺式作用元件、存在互作的转录因子、miRNA和circRNA等信息。利用Mapteace、TBtools、Cytoscape、OmicShare等作图软件及工具绘制染色体定位图、进化树、DIR与转录因子关系图、ceRNA网络等。采用NCBI的基因数据库结合转录组测序及qRT-PCR试验研究DIR基因家族在逆境下的表达情况。【结果】共鉴定到番茄中27个DIR基因,将其命名为SlDIR1—SlDIR27,分别位于12条染色体上,且大部分基因位于染色体末端,其基因结构、基序和结构域相对保守,其中22个SlDIR具有一个外显子的经典结构。番茄DIR基因家族同拟南芥的共线性关系远高于水稻和大豆。基于系统进化关系将27个番茄DIR成员分为3个不同的亚家族。转录组数据表明大部分DIR基因在番茄根部具有较高的表达量。此外,DIR基因的启动子区域含有多个与干旱、低温等非生物胁迫,以及MeJA、ABA、SA等激素诱导相关的顺式作用元件,且预测到与激素、生长发育、非生物胁迫相关的ERF、E2F/DP、MYB等互作转录因子。结合转录组数据分析,分别有5、10、10和13个SlDIR在农药、干旱、盐和冷胁迫后显著上调表达,其中,SlDIR23受到以上4种胁迫的诱导表达,而SlDIR8、SlDIR13和SlDIR20特异性响应冷胁迫的诱导,SlDIR17特异性响应盐胁迫。番茄DIR的ceRNA调控表明,miR-156与靶基因SlDIR8可能共同作用调控番茄的逆境胁迫。【结论】共鉴定出番茄DIR家族基因27个,不均匀地分布在12条染色体上,在根部有较高的表达量。SlDIR1、SlDIR13和SlDIR14等具有MeJA、ABA、SA等激素响应元件,其中,SlDIR6仅具有MeJA元件,SlDIR27仅具有SA响应元件。另外,SlDIR2、SlDIR14、SlDIR23等参与干旱、盐、低温等多种逆境胁迫,其中SlDIR23在不同胁迫处理下均可被激活。此外,DIR基因和转录因子、非编码RNA相互作用,共同参与调控番茄植株的逆境胁迫。     

辣椒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具有组织表达差异性,协同调控了辣椒的生长发育。     

苦荞蔗糖转运体家族FtSUCs的鉴定与生物信息学分析

植物蔗糖转运体家族蛋白( sucrose transporter, SUC)是介导蔗糖跨膜运输的重要载体。本文基于苦荞幼苗发育过程中转录组测序获得的数据库,筛选鉴定到10个FtSUCs基因,序列长度在1 436~7 109 bp,编码氨基酸148~605个,其中8个FtSUCs可能定位于内质网,其序列一致性为34.30%,保守位点达到20个。苦荞FtSUCs与拟南芥AtSUCs基因被聚为3类,FtSUCs中4个归入类型Ⅰ,2个归入类型Ⅱ,4个归入类型Ⅲ。苦荞FtSUCs蛋白主要由α-螺旋、不规则卷曲、延长链与β-折叠组成,其中9个FtSUCs包含跨膜结构域。幼苗发育过程中,FtSUCs基因间表达模式存在差异,其中,FtPinG0002608200.01、FtPinG0001943900.01、FtPinG0001944100.01在子叶期高表达,FtPinG0000342600.01表达量逐渐提高,FtPinG0001943500.01和FtPinG0009584000.01在各时期稳定表达,推测它们可能在不同发育阶段发挥特定功能。苦荞中多个FtSUCs基因存在显著正相关或负相关关系,暗示FtSUCs基因间可能存在协同效应或功能互补,进而共同调控蔗糖的跨膜运输。