苹果WRKY基因家族生物信息学及表达分析

【目的】鉴定苹果（Malus domesticaBorkh.）基因组上132个WRKY基因,为研究苹果WRKY转录因子在非生物和生物胁迫以及生长和发育过程中的调控作用奠定相关理论基础,也为进一步分析苹果WRKY基因提供信息。【方法】利用HMMER 3.0软件,通过WRKY保守域全蛋白序列PF03106用于鉴定苹果WRKY基因。采用WebLogo 3、DNAMAN 5.0、MapInspect、MEME和MEGA5.1等软件对其蛋白序列进行生物信息学分析。采用RT-PCR技术检测苹果WRKY基因的组织表达情况。【结果】鉴定得到132个苹果WRKY基因。分组鉴定和进化树分析结果显示,苹果WRKY蛋白分为I、II和III类型,I组共有24个成员可进一步分为I-C和I-N亚组,其锌指结构是C₂H₂类型（CX₄CX_22-23HXH）。II组含有1个WRKY区域共有79个成员,可进一步分为II-a、II-b、II-c、II-d和II-e亚组,分别有8、12、31、14和14个成员,其锌指结构为C₂H₂类型（CX_4-5CX₂₃HXH）。III组共有29个成员,其锌指结构为C₂HC类型（CX₇CX_23-24HXC）;WRKY结构域分析显示,其高度保守,绝大多数都含有WRKYGQK七肽和锌指结构;染色体定位分析显示,苹果WRKY分布于苹果17条染色体中,呈不均匀分布。染色体1和9上分布最多,为13个;其次是染色体12,分布12个;染色体2、5和14分布最少,为4个;基因结构分析表明,MdWRKY基因家族多数由2-5个外显子组成,基因结构进化高度保守;保守元件分析表明,MdWRKY基因家族包含10个保守元件：元件1-6为WRKY盒;元件7-10为未知盒。MdWRKY基因家族都包含有WRKY盒,I组中含有2个WRKY盒,II-a和II-b亚组中含有未知元件8,III组中含有未知元件7和9。半定量结果显示,12个MdWRKY均在根、茎、叶、花和果中表达,且呈现出多种相对表达模式。【结论】苹果WRKY基因家族结构高度保守,可能参与调控苹果生长和发育等过程。     

冬瓜WRKY基因家族鉴定及表达分析

【目的】WRKY转录因子家族在植物生长发育、抵御胁迫等过程起着至关重要的作用。挖掘参与冬瓜（Benincasa hispida Cogn.）非生物胁迫的WRKY基因家族成员,探究其生物学功能。【方法】基于冬瓜基因组数据库鉴定冬瓜WRKY成员,利用生物信息学的方法系统分析其WRKY基因家族成员的蛋白理化性质、系统发育、保守基序、顺式作用元件,通过qPCR分析冬瓜WRKY的表达情况。【结果】冬瓜WRKY基因家族有57个成员,氨基酸数目在126～650之间,相对分子质量在13.92～71.68 kD之间;蛋白等电点在4.61～9.69之间。根据系统进化树将该家族分为3个类群,第二类群又分为5个亚组。冬瓜WRKY外显子有2～6个。染色体定位分析发现,有56个基因定位在12条染色体上,1个基因未定位在染色体上。保守基序分析显示,Motif 1和Motif 3存在于56个基因,且同一类群具有类似的保守基序结构。同时,冬瓜WRKY基因的启动子上均含有应激反应相关元件、激素响应元件。WRKY在冬瓜不同组织和果实发育时期特异表达。RT-PCR结果表明非生物胁迫和激素处理冬瓜可诱导部分WRKY基因的表达。...     

绿豆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转录因子的分子生物学理论。     

甘薯全基因组WRKY转录因子的基因鉴定与逆境胁迫表达分析

【目的】对甘薯(Ipomoea batatas)全基因组WRKY转录因子进行基因鉴定与逆境胁迫表达分析,为甘薯WRKY基因的应用提供参考。【方法】在甘薯全基因组序列的基础上,应用生物信息学方法对WRKY基因进行挖掘,分析基因潜在复制关系、保守结构域、亚家族系统进化树、保守基序和抗逆表达模式,并对甘薯SPF1基因与IbWRKY基因进行了比较和分析。【结果】甘薯全基因组上共有82个WRKY基因,可分为Ⅰ、Ⅱ和Ⅲ3个类型,其中类型Ⅰ和类型Ⅲ的成员数量分别为18和5,类型Ⅱ可进一步分为Ⅱ-a、Ⅱ-b、Ⅱ-c、Ⅱ-d和Ⅱ-e 5个亚类,其成员数量分别为4,14,20,8和13。甘薯WRKY基因不均匀分布于15条染色体上,其中25对基因存在潜在复制关系。保守结构域分析发现,大部分WRKY的结构域高度保守,但个别基因存在保守结构域缺失现象。拟南芥与甘薯的WRKY转录因子在系统进化树上呈现按物种少量聚集的现象。在甘薯WRKY转录因子家族中共发现10个保守基序,包括含有WRKY七肽结构域的基序1和基序7、含有锌指结构域的基序2和基序3,其中基序3为Ⅰ类型WRKY转录因子N端WRKY结构域所特有。对逆境胁迫下的甘薯转录组数据进行分析发现,甘薯苗期在蔓割病菌(Fusarium oxysporum f. sp.batatas)侵染后共有19个WRKY基因差异表达,甘薯块根贮藏期在低温胁迫下有34个WRKY基因差异表达。蛋白序列比较和分析结果表明,甘薯SPF1的氨基酸序列与IbWRKY32的氨基酸序列一致度最高,为85.4%。【结论】甘薯全基因组中鉴定得到82个WRKY基因,其结构域较为保守,在蔓割病胁迫与低温胁迫条件下均存在差异表达。     

番木瓜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是有潜力的抗炭疽病候选基因。     

盐芥GRP7基因的克隆和逆境表达分析

【目的】从抗逆植物盐芥中分离GRP7(glycine-rich RNA-binding protein7,GRP7)基因,为进一步揭示GRPs基因在盐芥逆境应答中的作用机制及生物学功能奠定基础。【方法】利用盐芥GRP7基因的EST序列设计特异引物,克隆GRP7基因的全长序列,对其保守结构域和系统进化树进行分析,并对TsGRP7基因的启动子区进行预测;利用qRT-PCR技术,检测GRP7基因在不同逆境胁迫、不同组织中的表达情况。【结果】TsGRP7基因编码区全长522bp,编码173个氨基酸;TsGRP7蛋白的N端第9~82位氨基酸之间有1个RNA识别基序(RRM),RRM中含有保守的RNP-1和RNP-2亚结构域;多重序列比对和系统进化树分析表明,盐芥和拟南芥亲缘关系较近;启动子序列分析表明,盐芥GRP7启动子区含有HSE热激响应元件、LTR低温应答元件等多个逆境相关的顺式作用元件,表明盐芥GRP7基因参与逆境响应。qRT-PCR分析表明,盐芥GRP7基因在不同逆境胁迫条件及不同组织中表达不同。【结论】盐芥GRP7基因参与低温、高盐、PEG等逆境响应。     

瓠瓜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组,具有瓠瓜特有的进化分支,在不同组织中的表达具有特异性。     

大豆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的表达具有组织特异性,参与干旱和盐胁迫响应。     

谷子转录因子SiNAC18通过ABA信号途径正向调控干旱条件下的种子萌发

【目的】谷子(Setaria italica L.)具有显著耐旱性。研究旨在通过反向遗传学方法分析并鉴定在干旱条件下影响植物萌发过程的重要调控因子,为研究作物干旱条件下种子萌发的调控机制创造条件。【方法】使用Clustal X 2.0和MEGA 5.05软件对谷子SiNAC18蛋白序列及其同源序列进行多序列比对,并构建系统进化树;利用real-time PCR方法检测SiNAC18在不同胁迫条件下的表达模式;通过瞬时转化的方法分析SiNAC18蛋白亚细胞定位;在拟南芥中过表达SiNAC18,分析SiNAC18的生物学功能;分析SiNAC18在转基因拟南芥中可能控制的下游基因。【结果】SiNAC18全长1 074 bp,编码由357个氨基酸组成的亲水性蛋白,分子量约为38.8 k D;系统进化树分析表明SiNAC18属于NAC转录因子家族第Ⅰ组的NAP亚组,与拟南芥基因At NAC29同源性最高;氨基酸序列比对结果显示,SiNAC18与其他物种包括水稻、拟南芥、大豆和玉米中同源性最高的NAC类转录因子蛋白的N端都具有A、B、C、D和E这5个保守结构域,蛋白C端具有高度多态性,证明SiNAC18的N端序列与其结合下游基因启动子元件相关;real-time PCR结果显示,SiNAC18在干旱(PEG)、ABA、高盐(Na Cl)及过氧化氢(H2O2)处理条件下的表达量明显上升;亚细胞定位结果表明SiNAC18蛋白定位于细胞核中;基因功能分析结果显示,在ABA和PEG胁迫处理下,SiNAC18转基因拟南芥与野生型种子的萌发率存在明显差异:在正常生长条件下,野生型拟南芥WT和SiNAC18转基因拟南芥的萌发率基本一致,在PEG浓度为10%和15%的MS培养基上,SiNAC18转基因拟南芥的萌发率显著高于WT。在2和5μmol·L-1 ABA处理条件下,转基因拟南芥的萌发率显著低于WT;下游基因表达分析结果显示,ABA信号途径相关基因At RD29A,脯氨酸合成相关基因At P5CR和At PRODH以及过氧化物酶基因At PRX34在SiNAC18转基因株系中的表达量高于WT中的表达量,表明SiNAC18通过调控这些下游基因影响转基因植物在干旱条件下的萌发率。【结论】谷子NAC类转录因子基因SiNAC18可能通过ABA信号途径、氧化胁迫调控等途径正向调控植物在干旱条件下的萌发过程。     

谷子SiARGOS1的克隆、表达分析和功能标记开发

【目的】从谷子中分离受激素诱导表达、参与器官大小控制的拟南芥ARGOS(Auxin-regulated gene involved in organ size)基因家族的同源基因,进行生物信息学分析,明确其在不同组织器官及其受植物激素诱导的表达模式,分析基因编码区及其启动子序列差异,开发功能标记,为谷子产量性状相关基因的改良提供依据。【方法】通过对已有ARGOS蛋白保守结构域进行BLAST,明确谷子ARGOS家族成员数目并进行蛋白序列分析,采用同源克隆方法获得谷子ARGOS家族成员之一——SiARGOS1编码区及其启动子序列,用生物信息学方法分析SiARGOS1启动子的顺式作用元件,通过实时荧光定量PCR分析该基因在谷子各器官中以及不同植物激素条件下的诱导表达模式,利用基因编码区及启动子序列的SNP和插入缺失序列开发分子标记,同时利用85份谷子品种的穗重(panicle weight,PW)、穗粒重(grain weight,GW)和千粒重(thousand-grain weight,TGW)等产量性状数据进行基因型间的差异显著性分析,挖掘用于检测该基因与谷子产量性状相关优异等位变异的功能标记。【结果】获得6个谷子ARGOS家族成员,均具有典型的保守OSR(organ size related)结构域,包含2个跨膜螺旋结构和1个高度保守富含亮氨酸区域,克隆了与拟南芥AtARGOS同源的家族成员之一——SiARGOS1编码区及其启动子序列,该基因位于谷子第8染色体上,开放阅读框为342 bp,无内含子,编码113个氨基酸,启动子区域为2 109 bp,含有与生长素、乙烯、茉莉酸和赤霉素等多种植物激素调控有关的元件。表达分析发现,SiARGOS1在谷子根、茎、叶和穗等器官中均有表达,在根中表达量最高,其次为茎和叶,穗中表达量最低。SiARGOS1对生长素吲哚乙酸(indole-3-acetic acid,IAA)不敏感,但受乙烯利(ethephon,ETH)上调表达。不同基因型谷子SiARGOS1序列分析发现,SiARGOS1编码区151 bp(起始密码子83 bp)处存在1个SNP(C/G),导致该基因第28个氨基酸发生突变(Ala/Gly),据此设计一个CAPS-AccⅡ标记;另外,启动子区存在19个SNP和2个In Del,根据-1 652—-1 651处(TA)_(2/3)和-1 165—-1 163处(TCA)_(1/2)的序列差异分别设计SSR引物AP-1和AP-2。同时,用这些标记对85份谷子品种进行检测,CPAS-AccⅡ和AP-1检测到不同基因型的穗重、穗粒重和千粒重的差异均不显著,而AP-2检测的2种基因型间除千粒重差异不显著外,穗重和穗粒重在2015和2016两年间差异均达到显著水平。【结论】在谷子中发现6个ARGOS家族成员,均具有保守OSR结构域,其中,谷子SiARGOS1开放阅读框为342 bp,无内含子,与拟南芥AtARGOS同源,该基因对生长素不敏感,但受乙烯利上调表达。在该基因启动子区开发的SSR标记AP-2可作为功能标记,用于谷子穗重和穗粒重等产量性状相关优异等位变异的鉴定和筛选。     

茶树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的影响。     

基于盐胁迫转录组信息的蚕豆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基因在不同盐处理时间的表达量存在差异。     

大麦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基因存在一定的功能分化,具有明显的组织和时空特异性,推测其在大麦逆境响应和激素调节中具有重要调控作用。     

超表达苹果细胞分裂素响应基因MdCRF6影响花青苷积累和盐胁迫抗性

【目的】克隆苹果细胞分裂素响应因子基因Md CRF6(cytokinin response factor 6),鉴定其在调节花青苷积累和盐胁迫抗性中的作用,揭示Md CRF6的功能,为研究细胞分裂素信号途径和果树生长发育调控提供理论依据。【方法】以‘嘎啦’苹果(Malus domestica‘Gala’)为研究材料,利用同源序列比对和PCR技术,分离细胞分裂素响应基因Md CRF6。使用MEGA 5.0软件构建苹果Md CRF6与拟南芥CRFs间系统进化树;通过SMART软件和DNAMAN软件分析Md CRF6蛋白的保守结构域。利用实时荧光定量PCR方法检测该基因对细胞分裂素和盐胁迫的响应。通过电泳迁移率试验(EMSA),验证Md CRF6原核表达蛋白对DRE作用元件的绑定。构建Md CRF6植物超表达载体,并通过农杆菌介导的遗传转化获得转Md CRF6苹果愈伤组织。比较野生型和转基因苹果愈伤组织在花青苷积累和盐抗性方面的差异,结合基因表达分析,初步鉴定Md CRF6在调节花青苷积累和盐胁迫方面的生物学功能。【结果】分离得到了苹果细胞分裂素响应因子基因Md CRF6(基因序列号:MDP0000783818),该基因开放阅读框(ORF)为1 047 bp,编码含有348个氨基酸的蛋白。进化树分析和氨基酸序列比对结果表明,苹果Md CRF6蛋白在N端包含保守的CRF结构域,在C端包含保守的AP2/ERF结构域,并且与拟南芥At CRF6同源性最高。基因表达分析显示该基因具有细胞分裂素和盐胁迫响应,分别在10μmol·L-1 BA和100 mmol·L-1 Na Cl处理3 h和6 h时表达量最高。EMSA试验结果验证Md CRF6原核表达蛋白能够绑定DRE序列。在苹果愈伤组织中超表达Md CRF6,发现Md CRF6转基因苹果愈伤组织花青苷积累受到抑制,同时苹果愈伤组织的抗盐性降低。基因表达分析显示,Md CRF6显著抑制花青苷合成基因和盐响应相关基因的表达。【结论】苹果Md CRF6与拟南芥At CRF6具有高度的同源性,其参与植物对细胞分裂素和盐胁迫的响应。在苹果愈伤组织中超量表达Md CRF6抑制其花青苷积累,降低植物抗盐性。推测苹果Md CRF6可能通过结合花青苷合成基因以及抗盐相关基因的启动子,抑制基因的表达,从而负调节花青苷积累和抗盐性。     

陆地棉CRK基因家族的鉴定及其表达分析

【目的】富含半胱氨酸类受体激酶（CRK）是植物中最大的类受体激酶家族之一,在植物生长发育、激素信号传导和抗逆境胁迫中发挥重要作用。从全基因组水平鉴定陆地棉CRK基因家族并进行生物信息学和表达模式分析,为研究和利用陆地棉CRK基因家族奠定基础。【方法】从Pfam数据库下载stress-antifung结构域氨基酸序列,应用BLASTp程序搜索棉花基因组数据库,鉴定棉花CRK基因家族;利用Compute pI/Mw tool、SignalP、TMHMM Server V2.0、WoLF POSRT等在线工具预测陆地棉CRK家族蛋白的分子量、信号肽、跨膜结构域和亚细胞定位等;用ClustalX1.8软件对棉花和拟南芥CRK蛋白质进行氨基酸序列比对,MEGA5.0分析棉花和拟南芥CRK蛋白的系统进化关系;使用TBtools制作陆地棉CRK基因家族的染色体定位、基因结构和蛋白质结构域示意图;应用植物顺式调控元件数据库PlantCARE分析棉花启动子序列;通过植物磷酸化位点数据库PlantPhos预测陆地棉CRK家族蛋白的磷酸化位点;从NCBI数据库下载RNA-Seq数据,利用转录组定量工具Kallisto计算TPM值,通过在线工具Morpheus绘制陆地棉CRK家族基因表达热图。【结果】陆地棉基因组中有70个CRK基因,分布于14条染色体,其中52个基因（74.3%）集中串联成簇分布于A6/D6、A9/D9、A10/D10染色体,且在A/D染色体组之间呈现高度共线性关系。编码302-901个氨基酸,58个蛋白质（82.9%）具有跨膜结构域,主要定位于叶绿体、质膜和胞外。磷酸化位点预测结果表明,陆地棉和拟南芥CRK有5个相同的磷酸化位点基序,包括3种丝氨酸磷酸化位点基序和2种苏氨酸磷酸化位点基序。65个陆地棉CRK基因的启动子区（92.9%）至少含有一种逆境激素响应元件,69个基因启动子区（98.6%）至少含有一种生物或非生物胁迫响应元件。根据RNA-Seq数据分析结果,陆地棉CRK基因可分为3种不同的组织表达特征类型;盐、干旱、冷、热胁迫以及接种大丽轮枝菌均可以导致部分陆地棉CRK基因表达水平的改变。GhCRK25在根、茎、叶和胚珠中优势表达,在纤维中几乎不表达,ABA、GA3、SA、PEG-6000、氯化钠和大丽轮枝菌Vd991处理均能刺激GhCRK25迅速上调表达。应用病毒诱导的基因沉默技术（VIGS）沉默GhCRK25可导致棉花对大丽轮枝菌Vd991更为敏感。【结论】陆地棉CRK基因家族有70个成员,具有保守的基因结构和功能结构域,多样化的组织表达特征,大多数基因受激素和逆境调控。     

水稻WRKY80转录调节蛋白基因的分离与表达模式

【目的】分离水稻WRKY80,分析其编码蛋白的序列结构特征,了解其在各器官中以及病原接种、激素处理下的表达模式,为阐明其生物学功能奠定基础。【方法】基于水稻基因组数据库注释的Loc_Os03g63810基因的编码序列设计特异引物,用RT-PCR法从茉莉酸甲酯（MeJA）处理的水稻叶片中克隆WRKY80 cDNA 序列,运用生物信息学手段对推定的蛋白和启动子序列进行分析,用Northern杂交或实时荧光定量PCR方法研究基因的表达模式。【结果】 获得了WRKY80 cDNA序列,长1 392 bp,包含1个长1 164 bp 的完整开放读码框,编码387个氨基酸。WRKY80具有1个典型的WRKY 保守结构域,锌指类型为C2H2,归类于WRKY 第Ⅱ组;C-端为酸性结构区,且含连续的6个谷氨酰胺和8个丝氨酸,提示具有转录激活活性;预测其定位于细胞核。WRKY80与玉米和高粱等单子叶植物WRKY 的氨基酸序列同一性较高。WRKY80在各器官中组成性表达,在叶、根和穗中表达丰度较高,花中次之,在茎和颖果中丰度较低,且表达具有发育阶段相关性,在成熟叶、根中的表达分别高于幼叶和幼根;受稻瘟病菌、纹枯病菌、MeJA和乙烯利诱导表达,而水杨酸对其表达无影响。其表达模式与启动子顺式元件预测分析的结果基本一致。【结论】WRKY80具有转录因子的结构特征,推测其可能通过茉莉酸/乙烯介导的信号途径参与水稻的防御反应,也可能参与发育调控。     

苹果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条染色体上均有分布,多数基因具有在茎中高表达的组织表达特征,多数基因受逆境和轮纹病菌调控。     

白菜硫代葡萄糖甙合成中MAM 的进化分析

【目的】探索硫代葡萄糖甙合成中关键基因-甲基硫代烷烃基苹果酸合成酶（methylthioalkylmalate,MAM）基因在白菜基因组中的进化情况,为研究白菜MAM的功能提供理论指导。【方法】根据白菜数据库BARD信息分析拟南芥、盐芥和白菜MAM的共线性关系,利用MEME在线软件预测MAM序列的高度保守基序（motif）,通过R软件分析93份白菜材料MAM的表达模式,并利用MEGA5.05软件构建MAM的系统进化树。【结果】白菜的3个亚基因组均保留着MAM的同源基因。其中与拟南芥具有共线性的5个同源基因可能源自基因组古三倍化复制,其它两个同源基因可能是由转座扩增产生的;基序分析结果表明,与拟南芥的MAM相比,白菜的MAM同源基因存在基序的缺失,并且大多发生在序列的两端;在93份白菜材料中,白菜MAM同源基因之间的表达量有很大差异,Bra021947与Bra018524只在少数材料中存在极低的表达,并且Bra029356的表达量低于可检测水平;进化树分析发现,除Bra018524外,其他的白菜MAM同源基因与拟南芥的MAM聚集在不同的分支上。【结论】白菜的MAM可能在拟南芥和白菜物种分化之后独立起源并进化出不同的功能。