瓠瓜KNOX基因家族全基因组鉴定及组织表达分析

【目的】 探明瓠瓜KNOX基因家族特征,在瓠瓜全基因组水平上对KNOX基因家族成员进行鉴定和分析,研究其组织表达模式。【方法】利用生物信息学方法,对瓠瓜KNOX家族基因成员进行鉴定、编码蛋白理化性质分析及家族成员结构域和保守基序分析,利用MEGA 5.05软件构建系统进化树,采用实时荧光定量PCR技术,分析KNOX基因在瓠瓜不同组织（根、茎、叶、花、果实）中的表达情况。【结果】鉴定到12个瓠瓜KNOX家族基因,这些基因分布在6条染色体上,且均定位于细胞核。大多数瓠瓜KNOX家族蛋白含有同源异型盒蛋白特有的4种结构域KNOX1、KNOX2、ELK 和Homeobox_KN。系统进化分析将瓠瓜KNOX基因分为KNOX Class Ⅰ和KNOX Class Ⅱ 2大类群,2个类群又可进一步分为5个分支,其中Class Ⅰ-C为瓠瓜特有的进化分支。KNOX基因启动子区有多个激素响应元件和植物生长相关元件。组织表达分析发现,瓠瓜KNOX基因具有不同的表达模式,KNOX Class Ⅰ类基因表达具有组织特异性,在根、茎和果实中表达量较高,在其余组织中表达量较低或不表达,其中Lsi04G014450在根和茎中表达量最高,Lsi11G006380在果实中表达量最高;KNOX Class Ⅱ类基因在所有组织中均有较高表达。【结论】瓠瓜KNOX基因家族有12个成员,分布于6条染色体上,在进化关系上可分为5组,具有瓠瓜特有的进化分支,在不同组织中的表达具有特异性。     

甘薯响应蔓割病病原菌侵染的IbWRKY7基因克隆与表达分析

为探究IbWRKY7基因在甘薯响应蔓割病病原菌侵染过程中的表达模式,以高抗蔓割病品种‘鄂薯11’为供试材料,克隆到1个新的WRKY家族基因IbWRKY7,进行生物信息学分析;并对‘鄂薯11’和蔓割病敏感品种‘栗子香’在蔓割病病原菌侵染后IbWRKY7基因的表达模式进行分析。结果显示,IbWRKY7的CDS序列全长为945 bp,编码314个氨基酸;推测其编码不稳定的亲水性蛋白,蛋白分子质量为33.92 kU,pI 9.82,含有1个WRKY七肽保守序列和1个C₂HC型锌指结构。IbWRKY7基因上游2 000 bp的启动子区域内存在多种类型的顺式作用元件,如植物激素应答元件Myb、胁迫响应元件MYC和MYB等。多序列比对及系统进化树分析结果表明,IbWRKY7蛋白与日本牵牛花InWRKY7和三裂叶薯ItWRKY7亲缘关系最近。亚细胞定位预测显示IbWRKY7蛋白定位于细胞核。实时荧光定量PCR结果表明,与蔓割病病原菌侵染0 h相比,侵染2、4、12、24、48、72、96和120 h后,‘鄂薯11’的IbWRKY7基因表达量均显著提高(P<0.05);而‘栗子香’的IbWRKY7基因表达量除侵染24 h外的其他7个时间点均显著高于0 h。侵染2~48 h,‘鄂薯11’的IbWRKY7基因表达量显著高于‘栗子香’。双因素方差分析结果表明,不同品种和侵染时间点的IbWRKY7基因表达量均存在显著差异。综上,甘薯IbWRKY7基因上游2 000 bp的启动子区域含有12种激素应答和胁迫应答相关的顺式作用元件。IbWRKY7具有WRKY家族的典型结构特征,属于第Ⅲ类WRKY蛋白,氨基酸序列与同源物种相似度高,进化上高度保守。IbWRKY7基因受蔓割病病原菌侵染诱导表达,且在不同蔓割病抗性的甘薯品种中表达量差异显著。     

菠萝R2R3-MYB基因家族S20亚族的鉴定

【目的】鉴定菠萝R2R3-MYB基因家族S20亚族成员,分析其进化、结构与表达模式,为菠萝S20亚族的生物学功能研究奠定基础。【方法】以菠萝基因组数据库为基础,利用生物信息学方法对菠萝S20亚族基因进行鉴定,并对其系统进化、结构、保守基序和表达模式进行系统分析。同时,克隆并原核表达代表性基因AcMYB108a,利用Western blot进行杂交验证。【结果】在菠萝基因组上共鉴定到5个S20亚族基因,分布于5条染色体上。系统进化与保守基序分析发现,不同植物的S20蛋白具有较高的保守性且均具有motif3保守基序。顺式作用元件分析表明,菠萝S20亚族基因与植物逆境胁迫响应过程有关,且受乙烯诱导。基因表达分析结果显示,S20亚族基因的表达存在组织特异性,其中AcMYB108a在果实中的表达水平最高。利用pET原核表达系统成功表达了AcMYB108a重组蛋白,并通过了Western杂交验证。【结论】菠萝S20亚族基因与逆境胁迫响应相关,且受到乙烯诱导。     

木薯SR45亚家族基因鉴定及表达

目的 精氨酸和丝氨酸富集蛋白(Arginine/serine-rich proteins,SR)为剪接复合体的主要成员,不仅参与植物前体mRNA可变剪接过程,在植物非生物胁迫中也具有相当重要的作用。SR45基因为SR基因亚家族成员之一。本研究旨在分析木薯SR45亚家族成员蛋白结构特征与表达模式,为进一步了解该亚家族基因在木薯中的功能提供理论支持。方法 利用生物信息学技术重新构建木薯SR基因家族进化树,对木薯SR45亚家族蛋白理化性质、基因结构、保守结构域进行分析,同时利用转录组数据分析低温与干旱胁迫下SR基因表达变化,运用实时荧光定量PCR(RT-qPCR)研究各基因成员在不同组织中的特异性表达以及对低温胁迫的响应。结果 木薯SR基因家族共7个大类26个成员,SR45亚家族共有5个成员;SR45亚家族基因编码蛋白长度为135~423 aa,相对分子质量为14970~47210,等电点为5.19~12.34;预测其主要定位于细胞核和叶绿体;SR45编码RS结构域和RRM结构域,具有Motif 3、Motif 6、Motif 2和Motif 1保守基序。转录组数据和RT-qPCR分析表明,SR45基因均响应木薯低温胁迫,且MeSR45-2显著上调表达,MeSR45-4、MeSR45-5在木薯根和叶中有较高表达量。结论 木薯SR45亚家族基因显著响应低温胁迫;将MeSR45-2基因列为调控低温逆境变化的候选基因,将根和叶列为研究SR45基因亚类成员的主要组织。本研究为探索木薯SR45基因奠定了理论基础,并为进一步研究木薯SR45基因逆境胁迫应答过程指明了方向。     

巨桉GRF基因家族生物信息学分析及其在不同氮水平下的组织表达模式

【目的】从全基因组水平上鉴定巨桉生长调节因子（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蛋白均具有CX₉CX₁₀CX₂H和QX₃LX₂Q保守基序。顺式作用元件分析表明,EgrGRF基因启动子区含有茉莉酸甲酯、脱落酸、分生组织表达及低温等响应元件。实时荧光定量PCR结果显示,EgrGRF2和EgrGRF6基因在低氮处理的巨桉叶片中优势表达,EgrGRF2、EgrGRF3、EgrGRF5和EgrGRF6基因在低氮处理根中的表达量较高氮和常规氮处理显著升高,EgrGRF4基因在常规氮处理茎中的表达量最高。【结论】鉴定获得6个巨桉GRF基因家族成员,其中EgrGRF4基因可能参与茎生长发育的调控,EgrGRF2、EgrGRF3、EgrGRF5和EgrGRF6基因可能参与了巨桉对低氮胁迫的响应。     

大豆BBX基因家族的鉴定及表达

【目的】对大豆（Glycine max）全基因组BBX（B-box）蛋白基因进行鉴定和生物信息学分析,探究其对非生物逆境胁迫的响应模式,为大豆BBX基因的应用提供参考。【方法】利用生物信息学方法鉴定大豆BBX基因家族成员,分析其保守结构域、系统进化树、基因复制关系、顺式作用元件、组织特异性表达和非生物逆境胁迫（干旱和100 mmol/L NaCl胁迫,均处理0,1,6,12 h）响应模式。【结果】大豆BBX基因家族（GmBBXs）有42个成员,分布在除2和16号染色体外的18条染色体上。保守结构域分析表明,GmBBXs在N端有一个或两个B-box结构域（B-box1和B-box2）,B-box1较B-box2保守;部分成员在C端存在一个CCT结构域。系统发育分析表明,GmBBXs家族成员分为5个亚家族,分别为B1+B2+CCT（Ⅰ型和Ⅱ型）、B1+CCT（Ⅲ型）、B1+B2（Ⅳ型）和B1（Ⅴ型）类型,其成员数量分别为3,8,6,18和7。共线性分析发现,在GmBBXs基因的扩展中存在串联复制和片段复制事件,且片段复制事件比例较串联复制大。GmBBXs基因启动子上含有多种顺式作用元件,如光应答元件、逆境胁迫应答元件和激素响应元件等。GmBBXs基因表达分析结果显示,GmBBXs在大豆根、茎、叶、花、荚果和种子中的表达存在差异,不同GmBBXs响应干旱和盐胁迫的方式不同,获得的36个GmBBXs基因表达数据可分为6种表达模式。【结论】在大豆中共鉴定得到42个GmBBXs基因,分布于18条染色体上,其结构域较为保守;GmBBXs的表达具有组织特异性,参与干旱和盐胁迫响应。     

石榴低温响应因子CBF基因家族鉴定及其表达分析

CBF（C-repeat binding factor）是AP2家族的一类转录激活因子,在植物响应低温胁迫和提高植物耐寒性方面具有重要作用。虽然多个石榴品种的基因组已经发布,但仍缺乏对其CBF基因家族的全面研究。为全面分析石榴（Punica granatum L.）CBF基因家族分子生物学特性,本文对鉴定到的PgCBFs成员进行了生物信息和表达分析。结果表明,石榴全基因组中有7个CBF基因家族成员,蛋白序列中除包含AP2结构域外,还具有CBF特征序列PKKPAGRxKFxETRHP和DSAWR;7个成员集中分布在染色体1和4上,编码蛋白质氨基酸长度为201~267 aa,分子质量为21.69~29.81 ku;进化分析显示PgCBFs基因家族成员分布在Group Ⅱ~Ⅳ 亚组中,与水稻CBF成员组成的Group Ⅰ存在明显区分;除 PgCBF2、 PgCBF4、 PgCBF5基因结构中含有1~2个内含子外,其余与其他物种类似,属内含子缺失型;编码的蛋白序列均含有保守基序Motif 1~7;蛋白二级结构主要为不规则卷曲和α-螺旋,三级结构较为相似。PgCBFs基因家族扩增主要来源于串联重复,并与拟南芥、苹果、桃分别存在2对、5对和4对共线性关系;GO注释显示PgCBFs多与转录调控、低温胁迫或激素诱导相关;启动子区含有多种顺势作用元件,主要为胁迫刺激、激素诱导和光响应元件;基因表达分析发现除 PgCBF6外,其余成员在根中高度表达,在低温胁迫处理下PgCBF在根和韧皮部表达量显著上调,其中 PgCBF7能够快速响应并持续应答低温诱导。结合进化树、共线性、启动子以及表达分析,推测 PgCBF7可能与石榴幼苗低温胁迫调控相关。本研究可为深入研究石榴CBF基因家族生物学功能提供理论基础。     

绿豆WRKY基因家族的全基因组鉴定及生物信息学分析

【目的】对绿豆WRKY转录因子在全基因组上进行鉴定及生物信息学分析,为绿豆WRKY基因功能的研究以及高抗绿豆品种培育奠定理论基础。【方法】利用Hmmer软件同源搜索,并通过在线工具SMART、NCBI-CDD和Pfam数据库进行结构域再确认获得79个绿豆WRKY基因,然后通过在线软件ProtParam、ProtScale、SWISS-MODEL分析绿豆WRKY家族蛋白的理化性质、亲疏水性并构建3D立体模型,使用MEGA X、McScanX软件构建进化树和绘制共线性图,利用MEME、PlantCARE在线工具分析WRKY转录因子的保守基序并预测基因上游的顺式作用元件。【结果】鉴定得到的79个VrWRKY蛋白的氨基酸数量为64~746,分子量为7.42~80.93 ku,等电点为4.49~10.41,大多数成员为亲水性不稳定蛋白。系统进化树分析表明,绿豆WRKY基因家族的79个成员可分为Ⅰ、Ⅱ和Ⅲ 3组,其成员数量分别为16,51和12,Ⅱ组进一步划分为Ⅱ-a、Ⅱ-b、Ⅱ-c、Ⅱ-d和Ⅱ-e 5个亚组,其成员数量依次为2,17,17,7和8,其中Ⅱ-a、Ⅱ-b、Ⅱ-d和Ⅱ-e亲缘关系更近,而Ⅱ-c在进化树上的分布整体更靠近Ⅲ组。基因结构及保守基序分析显示,该家族基因均为间隔基因,且都含有保守基序WRKYGQK。基因染色体定位分析表明,所有成员在11条染色体上均有分布。顺式作用元件分析表明,绿豆WRKY基因启动子区存在丰富的激素及胁迫响应相关的元件。基因组间共线性分析结果表明,绿豆WRKY基因与拟南芥WRKY基因之间存在同源关系。基因复制结果表明,绿豆基因组有1对串联重复基因,15对大片段复制基因,片段复制在绿豆WRKY基因家族扩张中起主要作用。绿豆WRKY蛋白3D结构模拟表明,所有成员具有相似的两种蛋白结构,其中一种立体结构中心含有1个Zn²⁺。【结论】鉴定得到79个含典型WRKY保守结构域的绿豆WRKY蛋白,其生物信息学分析丰富了绿豆WRKY转录因子的分子生物学理论。     

苋菜AmSPL基因家族转录组鉴定及表达分析

为明确SPL(Squamosa Promoter Binding Protein Like)在苋菜中的功能,基于‘全红’苋菜转录组数据库筛选获得17个SPL家族成员,并进行生物信息学分析预测。结果表明,17 个AmSPL蛋白分别由186~777 个氨基酸组成,相对分子量从20.53~88.27 ku。构建AmSPL与拟南芥(Arabidopsis thaliana L.)16个AtSPL成员系统进化树,并根据AtSPLs基因家族的分类标准将AmSPL成员分为6组。实时荧光定量PCR分析表明,在黑暗条件下,AmSPLs的表达量均高于蓝光下,推测AmSPLs的表达可能通过某种方式受到蓝光的负调控。预测分析其中的8 个AmSPL可能是miR156的靶基因。     

大豆GmPIN2b调控根系响应低磷胁迫的功能研究

目的 研究生长素转运蛋白PIN基因家族在大豆Glycine max 根系适应低磷胁迫中的功能。方法 对大豆23个GmPIN家族成员进行进化树分析和表达模式分析,并进一步对GmPIN2b进行功能分析。结果 大豆23个GmPIN家族成员分散在7个不同的亚族,其中,GmPIN2a、GmPIN2b和GmPIN9a、GmPIN9d与拟南芥Arabidopsis thaliana的AtPIN2位于同一亚族。不同GmPIN家族成员在大豆中的组织表达定位及其受低磷调控的表达模式不同,其中,GmPIN2b在大豆根系中的表达水平在低磷胁迫6 d后显著上调,回补表达GmPIN2b能够部分恢复Atpin2突变体的表型。无论在低磷还是高磷条件下,回补表达GmPIN2b的转基因植株鲜质量和主根长均显著高于Atpin2突变体;而且回补表达GmPIN2b显著提高了Atpin2突变体在低磷条件下的一级侧根数,以及高磷条件下根系对重力的敏感性。结论 GmPIN2b在大豆根系形态建成响应低磷胁迫的过程中起重要的调控作用。     

猕猴桃 AcNPR1基因克隆及抗病响应表达分析

为探究猕猴桃 AcNPR1基因在抗病响应中的作用及其在不同抗性猕猴桃品种间抗性差异机制,以高感病品种‘红阳’和抗病品种‘徐香’猕猴桃叶片为材料克隆 AcNPR1基因序列;利用ExPASy-ProtParam tool、MEGA 7.0 等对‘红阳’AcNPR1蛋白的理化性质进行分析、亚细胞定位进行预测、进化关系等进行分析;通过RT-qPCR分析 AcNPR1基因在不同品种猕猴桃植株中的组织表达模式以及病原菌和水杨酸（Salcylic acid,SA）对其诱导表达模式。结果表明, AcNPR1基因开放阅读框1 770 bp（Genbank 登录号MW881148）,编码589个氨基酸,理论等电点6.27,具有典型的BTB/POZ结构域、ANK重复序列、NPR-like等NPR1蛋白保守结构域,AcNPR1蛋白预测定位于细胞核中;蛋白质的二级结构以α螺旋和无规卷曲为主;系统进化显示其与茶树的亲缘关系最为密切,相似性为83.94%。组织特异性分析表明 AcNPR1基因在‘红阳’和‘徐香’品种叶片中的表达水平最高。在丁香假单胞杆菌猕猴桃致变种（Pseudomonas syringae pv.actinidiae,Psa）处理下,抗病品种‘徐香’ AcNPR1基因的相对表达水平在处理0 h后迅速增加,为对照的3.6倍;而高感品种‘红阳’在48 h后显著增加,为对照的2.7倍。在SA+Psa处理下,‘徐香’ AcNPR1基因相对表达水平在24 h内达到最高水平,是对照的7.7倍;‘红阳’在72 h后达到最大值,仅为对照的1.6倍。 ‘徐香’ AcNPR1基因的抗病响应反应早于‘红阳’且更易受SA诱导表达。 AcNPR1基因在猕猴桃抗病胁迫方面具有一定作用,在不同抗性的猕猴桃品种抗病机制中存在差异。     

青稞籽粒花青素合成相关基因 HvnF3′M的克隆与表达分析

旨在从紫粒青稞‘达章紫’和白粒青稞‘昆仑12号’中分别克隆具有典型的P450超家族结构域的 F3′M 基因。序列分析结果表明,该基因序列全长为 1 584 bp,共编码527个氨基酸。序列比对发现,‘达章紫’与‘昆仑12号’中的 HvnF3′M序列存在2个碱基差异,一致性为99.87%;氨基酸序列也存在2个差异,一致性为99.62%。理化性质分析表明,两蛋白理论等电点分别为6.85和7.19,分子质量分别为57.01 ku和56.89 ku。蛋白质序列分析表明,HvnF3′M为亲水性稳定的碱性蛋白,‘达章紫’中该蛋白主要由α-螺旋、无规卷曲、延伸链和β转角组成;‘昆仑12号’中该蛋白主要由α-螺旋、无规卷曲和延伸链组成。系统进化分析表明,在11种禾本科植物种中,两品种的 F3′M蛋白与大麦的亲缘进化关系最近,与玉米亲缘关系最远。荧光定量PCR（RT-qPCR）分析表明,随着籽粒颜色的形成（尤其是在晚期）,‘达章紫’ HvnF3′M基因的表达量呈现极显著升高（P＜0.01）;而‘昆仑12号’ HvnF3′M基因表达量呈现逐渐下降的趋势。从研究结果推测, HvnF3′M基因与青稞籽粒颜色的形成有着密切联系,且在籽粒花青素形成过程中起正向调控作用。     

小麦OSCA基因家族全基因组鉴定及表达分析

【目的】在全基因组水平对小麦(Triticum aestivum L.)OSCA家族成员（TaOSCAs）进行鉴定,以明确TaOSCA在小麦中的潜在生物学功能。【方法】利用BLASTP结合保守结构域筛选的方法鉴定小麦TaOSCA蛋白。用生物信息学方法分析TaOSCA基因家族成员系统发育关系,基因结构及跨膜结构域;利用公开的转录组数据分析TaOSCA在小麦不同发育阶段不同器官/组织的表达谱;以小麦品种中国春为材料,进行PEG-6000模拟的渗透胁迫处理,用实时荧光定量PCR检测叶片TaOSCAs在渗透胁迫下表达模式。【结果】从小麦中鉴定到35个TaOSCA蛋白;TaOSCAs分为TaOSCAⅠ、TaOSCAⅡ、TaOSCAⅢ和TaOSCAⅣ 4个亚家族,分别有15,16,3和1个成员。基因结构分析结果表明,除TaOSCA4.1-4B外,其余TaOSCA均包含多个外显子。跨膜结构域分析结果表明,除TaOSCA2.3-4B有5个跨膜结构域（TMs）外,其余TaOSCA成员均包含8～11个TMs。TaOSCAs在小麦不同发育时期各器官中的表达结果表明,TaOSCA1.1-3B/3D、TaOSCA1.2-1A/1B/1D、TaOSCA1.3-1A、TaOSCA1.6-1A/1D和TaOSCA3.1-2B/2D在全生育期高水平表达,TaOSCA1.4-2A/2B/2D在不同发育阶段的根组织中均高表达,TaOSCA1.5-1B/1D在生殖生长期穗中高表达,TaOSCA1.6-1B、TaOSCA2.6-3B/3D在生殖生长期籽粒中高表达,表明这些基因可能在根、穗及籽粒发育过程中发挥重要功能。渗透胁迫处理下,TaOSCA2.4和TaOSCA4.1表达下调,TaOSCA2.5的表达无显著变化,TaOSCA1.1、TaOSCA1.2等11个TaOSCA基因表达上调,说明这11个基因可能在小麦苗期响应渗透胁迫中发挥作用。【结论】TaOSCA基因可能在小麦响应渗透胁迫的过程中发挥重要作用。     

杧果C2C2型锌指蛋白MiLSDs的生物信息及胁迫响应分析

LSD（lesion simulating disease,LSD）基因编码一类特殊的C2C2型锌指蛋白,在细胞程序性死亡（Programmed cell death,PCD）和其他生物过程中发挥重要作用。通过鉴定杧果LSD基因家族成员并对其进行生物信息学分析,为探究杧果LSD转录因子在植物进化与发育中的生物功能奠定基础。根据杧果全基因组数据,对杧果LSD基因家族进行鉴定,包括理化性质、保守结构域、二级结构、三级结构、系统进化树分析等;采用实时荧光定量PCR方法,探究MiLSD基因在生物胁迫和激素处理下的表达模式。结果表明,杧果中共有3个LSD转录因子,其分子质量为10.36~17.58 ku,理论等电点（pI）为6.68~8.93;保守结构域、二级结构及三级结构分析发现MiLSDs成员具有较高保守性;系统发育树分析显示, MiLSD1、 MiLSD2、 MiLSD3分别分布在2个亚家族（ClassⅥ和ClassⅨ）。实时荧光定量PCR对胶孢炭疽菌（Colletotrichum gloeosporioides,Cg）、细菌性黑斑病菌（Xanthomonas citri pv. Mangiferaeindicae,Xcm）、水杨酸（Salicylic acid,SA）、茉莉酸甲酯（Methyl Jasmonate,MeJA）侵染下杧果叶片基因表达量的分析显示, MiLSD1在Cg侵染下表达量明显升高,说明 MiLSD1可能对胶孢炭疽菌侵染起正调控作用; MiLSD2在Xcm、SA处理下表达量也显著上调,说明 MiLSD2可能通过诱导SA来增加杧果对细菌性黑斑病的抗性; MiLSD3除在MeJA处理下表达量较高外其余3个处理均少量表达,说明 MiLSD3可能参与杧果对MeJA胁迫的响应而对于真菌和细菌侵染以及水杨酸胁迫并不表达;此外,在MeJA处理下 MiLSD1、 MiLSD2的表达量也显著上调,表明茉莉酸甲酯可能受MiLSDs的调控。     

东方百合‘索邦’组蛋白去乙酰化酶基因 LoSorHDA1的克隆与表达分析

以东方百合‘索邦’（Sorbonne）为材料,结合前期转录组数据,克隆组蛋白去乙酰化酶基因HDA1(Histone deacetylases),命名为LoSorHDA1。LoSorHDA1序列长1 518 bp,编码505个氨基酸,生物信息学预测含有一个高度保守的组蛋白去乙酰化结构域。系统进化分析表明,LoSorHDA1与拟南芥等物种的HDA1同源蛋白聚为一支,而与HDAC家族其他成员关系较远;LoSorHDA1编码的氨基酸序列与单子叶植物小果野焦（Musa acuminata）的相似度最高,而双子叶植物HDA1则聚为另一支。RT-PCR结果显示：LoSorHDA1在‘索邦’不同组织中普遍表达,其中在茎生根、嫩茎、花被片、花丝、柱头、花柱、子房中表达量较高,在下部叶和花药中表达稍弱。亚细胞定位表达结果表明,LoSorHDA1定位在细胞核和细胞质中。荧光定量PCR结果显示：LoSorHDA1在不同花发育阶段中的各个花部组织中均有不同程度的表达,在雌雄蕊中表达趋势相反,而在内外花被片中的表达趋势较为一致,说明LoSorHDA1可能参与东方百合花发育进程,这为后续LoSorHDA1的去乙酰化修饰在百合花发育中的功能研究奠定基础。     

番木瓜WRKY转录因子家族基因鉴定及其响应短孢炭疽菌侵染的表达分析

【目的】在全基因组水平鉴定番木瓜(Carica papaya L.) WRKY（CpWRKY）转录因子家族基因,并对其在短孢炭疽菌侵染下的表达量进行分析,为CpWRKYs 家族的功能研究和利用奠定基础。【方法】采用生物信息学方法,鉴定和筛选CpWRKYs转录因子家族基因成员,并对其进行系统进化分析、基因结构分析、蛋白保守基序预测和启动子顺式作用元件分析;同时以番木瓜炭疽病耐性品种和敏感性品种为材料,利用短孢炭疽菌侵染采后果实,于0,24,48 h 取样,对其转录组学数据进行分析,筛选2个品种中差异表达的CpWRKY基因,并利用实时荧光定量PCR(RT-qPCR)试验对筛选结果进行验证。【结果】经系统分析从番木瓜中共鉴定出48个CpWRKYs成员,分为Ⅰ、Ⅱ和Ⅲ 3类;其中第Ⅱ类成员最多,可分为Ⅱa、Ⅱb、Ⅱc、Ⅱd和Ⅱe 5个亚类。CpWRKYs家族成员的氨基酸残基数目为97~747 个,等电点为4.83~9.71,为亲水性蛋白,大部分CpWRKYs的结构域高度保守,但有个别蛋白保守域缺失。在 CpWRKYs家族中共预测到10个保守基序,其中包括含有WRKY七肽结构域的基序1、含锌指结构的基序2和同时含有七肽序列和锌指结构的基序3,其中基序3为大多数I类CpWRKY转录因子N端WRKY结构域所特有。CpWRKYs含有1~6个外显子,同一家族或亚家族成员在基因结构上表现出一定的相似性,同时在CpWRKYs启动子中检测到大量与生物胁迫和非生物胁迫相关的元件。对短孢炭疽菌侵染番木瓜的转录组数据进行分析,在番木瓜耐性品种中筛选出了特异表达的CpWRKY22、CpWRKY11、CpWRKY2、CpWRKY33和CpWRKY25,RT-qPCR检测结果显示,这些基因在耐性品种中特异上调表达,在敏感性品种中的表达量无明显变化。【结论】番木瓜WRKY家族共包括48个成员,该家族基因结构和蛋白基序具有组间差异性和组内保守性。CpWRKYs家族成员强烈响应炭疽菌侵染,其中CpWRKY22、CpWRKY11、CpWRKY2、CpWRKY33和CpWRKY25是有潜力的抗炭疽病候选基因。     

汉麻Alfin-like转录因子家族鉴定及表达分析

对汉麻Alfin-like转录因子家族进行鉴定和生物信息学分析,为Alfin-like转录因子调控汉麻生长发育及抗性奠定基础。通过利用拟南芥Alfin-like转录因子蛋白序列作为模板,利用TBtools、MEGA等生物信息学工具分析汉麻全基因组数据。在汉麻基因中共鉴定出5个CsALs基因,对其进行了理化性质分析、亚细胞定位、保守基序、基因结构、染色体定位、顺式元件、共线性分析、进化分析及聚类分析等基础分析。结果表明：鉴定的5个汉麻AL转录因子均含有DUF3594结构域（PAL结构域）和PHD手指状结构域;经理化性质分析得出,CsALs蛋白质长度在239～256 aa之间,等电点在4.97～5.31之间,亲水性为-0.809～-0.623,属于酸性亲水蛋白;CsALs基因分布在4条染色体上,亚细胞定位发现AL转录因子均定位在细胞核;上游2 000 bp启动子元件分析与光响应、胁迫响应和激素相关;根据与拟南芥、大豆、水稻、玉米、烟草等系统发育分析共分为B、E、F 3个亚族;热图聚类分析表明CsALs基因在汉麻中具有较高表达。     

结球甘蓝SET基因家族鉴定与表达分析

植物SET基因是一类含有SET结构域的高度保守的基因家族,参与调控植物多种生命过程。为探讨结球甘蓝SET基因家族及其表达情况,采用生物信息学方法对SET基因家族进行鉴定,并利用实时荧光定量PCR（qRT-PCR）和半定量PCR（RT-PCR）对SET基因在不同组织中的表达模式进行分析。结果鉴定出28个结球甘蓝SET基因,它们不均匀分布在除1号染色体外的其他8条染色体上,编码蛋白的相对分子质量在29.20 ~ 185.16 kDa之间,预测等电点在5.02 ~ 9.06之间。基因结构分析表明结球甘蓝SET基因有0 ~ 23个内含子。系统进化分析将该家族成员分为7个亚族,其中第Ⅱ和第Ⅴ亚家族成员最多。在7个亚家族中,每个家族挑选1个基因,qRT-PCR结果显示选取的7个SET基因在结球甘蓝的叶片中不表达,而在根、茎、花蕾和花中均有表达,其中花蕾中表达量较高,进一步RT-PCR检测发现其在花粉母细胞时期高度表达。对结球甘蓝SET基因家族进行了鉴定,分析其进化关系和表达模式,可为深入研究该家族基因功能奠定一定的理论基础。     

水稻lncRNA SVR及邻近SAUR基因在种子低温萌发中的表达

目的 长链非编码RNA(Long non-coding RNA，lncRNA)长度大于200 nt，一般不具有编码蛋白质功能。探究低温下lncRNA对邻近SAUR基因表达的影响，以期为研究水稻种子低温萌发的能力提供理论依据。方法 lncRNA SVR由华南农业大学国家植物航天育种工程技术研究中心前期研究筛选，可以响应低温胁迫，通过生物信息学方法分析lncRNA SVR的二级结构并寻找lncRNA SVR内部的冷胁迫基序，通过qRT-PCR分析lncRNA SVR和SAUR基因的表达特性。结果 lncRNA SVR序列中存在高度相似的冷胁迫响应基序，且在茎环连接处。表达特性分析表明，在种子萌发过程中低温胁迫会持续降低lncRNA SVR的表达，邻近的多个SAUR基因在低温胁迫下的表达量明显高于在常温下的表达量，表明lncRNA SVR的邻近SAUR基因一定程度上能够和lncRNA SVR一样响应低温胁迫。表达相关性分析表明，在低温萌发中lncRNA SVR与这些SAUR基因的表达均呈负相关关系，lncRNA SVR与OsSAUR55的表达呈显著负相关关系。进一步分析种子低温萌发中SAUR基因在lncRNA SVR敲除系中的表达情况，结果表明，在敲除系中，lncRNA SVR的下降幅度低于其在野生型中的下降幅度，OsSAUR41、OsSAUR53、OsSAUR54、OsSAUR55表达量的上升幅度均显著低于其在野生型中的上升幅度。结论 lncRNA SVR可能负调控OsSAUR55的表达，进而响应低温胁迫。