小麦叶锈菌休眠与萌发夏孢子的差异表达

【目的】小麦叶锈菌引起的小麦叶锈病是小麦的重要病害之一,孢子萌发和侵入气孔是保证叶锈菌正常侵染寄主的重要前提。本研究旨在研究小麦叶锈菌夏孢子萌发的差异表达特性,为揭示萌发的机制及其与致病的关系提供依据。【方法】利用小麦叶锈菌致病类型THTT的休眠夏孢子和萌发夏孢子进行RNA-seq分析。【结果】在萌发夏孢子中筛选出相比休眠夏孢子上调表达的基因为4400个,下调基因5325个。GO富集分析发现,上调差异基因主要涉及细胞进程、有机物代谢、信号传导、单个有机体过程、催化酶活性等;下调基因主要涉及单组织过程、单个有机体细胞过程、有机物代谢过程、催化活性调节等过程。利用KEGG分析发现,差异基因共参与109条通路,从中筛选出两条与孢子萌发相关的通路——丝裂原活化蛋白激酶(MAPK)信号通路和囊泡运输中的SNARE蛋白交流。10个基因的定量分析结果与基因数字表达谱分析结果一致。【结论】研究获得了叶锈菌孢子萌发及侵入气孔过程中的重要差异基因,研究结果为研究叶锈菌致病机制奠定基础。     

巴氏蘑菇子实体不同发育阶段的转录组分析

为探讨巴氏蘑菇子实体不同发育阶段基因的表达情况,本研究对巴氏蘑菇子实体不同发育时期（原基、采收期和开伞期）进行转录组测序,以本实验室已获得的巴氏蘑菇JA菌株的不育单孢菌株JA-15036基因组为参考基因组研究原基与采收期及开伞期样本间差异表达基因,并对差异表达基因进行了GO功能和Pathway富集分析。GO功能分析结果显示,差异表达基因主要富集在跨膜转运、碳水化合物代谢途径和膜组分,它们协同调控为子实体生长发育提供稳定的内环境。KEGG富集分析结果表明,原基期上调的差异表达基因主要富集在核糖体蛋白和DNA复制,表明原基期细胞代谢旺盛,其中核糖体蛋白基因上调为后期蛋白质合成提供重要场所;采收期和开伞期子实体时期差异表达基因主要富集在碳水化合物代谢、脂肪酸降解和氨基酸代谢等途径,为巴氏蘑菇子实体的生长发育与成熟提供营养与能量。     

CWI和HOG信号通路基因在斑玉蕈不同发育阶段的表达模式

为了更清楚地了解MAPK信号通路中的细胞壁完整性信号通路（cell wall integrity,CWI）和高渗透压甘油（high-osmolarity glycerol pathway,HOG）信号通路对斑玉蕈菌丝成熟、原基形成和子实体发育过程的影响及调节作用,对MAPK信号通路中的CWI和HOG信号通路基因在斑玉蕈不同菌丝培养时间（40、60、80和100d）和不同生长发育关键时期（24h、菌丝恢复期、菌丝转色期、原基期和子实体期）的表达模式进行分析,以期揭示这两条信号通路基因参与调节斑玉蕈菌丝的生长、子实体的形成和发育的作用。在斑玉蕈的CWI和HOG信号通路中经分析鉴定一共获得了15个关键基因。CWI信号通路基因表达分析表明：在菌丝培养的40-100d的过程中,大部分CWI信号通路基因在第60天时表达量最高,其中rho1、ssk1、ssk2和ste20的基因表达量上调了2-5倍,在第80-100天时出现持续下降。在HOG信号通路中的大部分基因也在菌丝培养的第60天表达量达到最高。其中sho1、ste20、ssk1和ssk2基因的表达量上调最为显著,而hog1基因的表达量在菌丝培养的第40-100天呈持续下降。子实体形成过程中两条通路的大部分基因在原基形成时期表达量最高,而在子实体时期表达量下调。其中HOG信号通路中的ssk2基因表达量上调最为显著。以上结果说明在菌丝生长过程中第60天时菌丝细胞生长增殖最为旺盛,而在第80天开始菌丝细胞基本开始停止生长,菌丝也逐渐达到成熟。同时在菌丝增殖生长过程中,斑玉蕈持续地上调CWI信号通路基因的表达来调控菌丝细胞壁的完整性,从而控制菌丝细胞壁的形成。其中bck1、mkk1和slt2基因可能对斑玉蕈菌丝细胞的分裂增殖和细胞壁的形成以及诱导子实体形成起到关键作用。     

乙醇胁迫抑制毕赤酵母表达外源蛋白的转录组学分析

在5 L发酵罐中进行毕赤酵母发酵表达猪?干扰素的实验,发现甘油培养末期乙醇的积累会抑制外源蛋白的表达。从转录组学角度系统分析不同浓度乙醇胁迫条件下,毕赤酵母甘油培养期和甲醇诱导期细胞的生理状态变化。研究结果表明,在甘油培养期,乙醇胁迫使得毕赤酵母细胞中的545个基因发生了显著差异表达(265个基因表达上调,280个基因表达下调),这些差异表达基因的功能主要涉及蛋白质合成、能量代谢、细胞周期和过氧化物酶代谢。乙醇胁迫增加了蛋白质错误折叠的情况,降低了核糖体和线粒体的结构完整性,使得甘油培养末期无法得到大量具有健全功能的酵母细胞。在甲醇诱导期,与甲醇代谢、蛋白质加工合成、氨基酸代谢等途径相关的294个基因发生了显著差异表达(171个基因表达上调,123个基因表达下调),导致内质网胁迫不能被及时解除,破坏了细胞内的氨基酸正常代谢。     

刺芹侧耳不同发育时期的转录组差异分析

为了探讨刺芹侧耳子实体生长发育时期的基因表达变化,本文利用高通量测序技术对刺芹侧耳不同发育时期（菌丝期、原基期、子实体时期）进行RNA-Seq分析,在转录水平上解析差异表达基因在刺芹侧耳生长发育过程中的作用和功能。KEGG功能富集显示,菌丝期差异表达基因主要富集在碳代谢和氨基酸代谢中,其中三羧酸循环中编码柠檬酸合酶、乌头酸水合酶、异柠檬酸脱氢酶、琥珀酰辅酶A合成酶、琥珀酸脱氢酶、苹果酸脱氢酶的基因表达量均上调,说明碳代谢和氨基酸代谢是菌丝时期的主要能量来源;原基期上调的差异表达基因主要富集在脂肪酸代谢,其中RT-PCR定量结果显示原基期编码脂肪酸合酶的基因和编码脂酰辅酶A合成酶的基因下调,编码超氧化物酶的基因和编码过氧化氢酶的基因上调,表明脂肪酸代谢和抗氧化酶对刺芹侧耳原基期维持机体的稳定和生物应激方面起着重要作用。子实体时期上调的差异表达基因主要富集在剪接体、类固醇的生物合成以及AMPK信号通路中,说明环境因子对子实体时期有一定的影响。     

茶树种子发育过程的转录组分析

该研究以‘铁观音’茶树品种的种子为试验材料,采用转录组测序技术分析种子发育的3个时期(幼果期、膨大期、成熟期)的表达差异,探究茶树种子油脂代谢的分子机制。结果表明:(1)经转录组测序、组装后共获得30 940 581个clean reads,经数据合并拼接最终得到36 951条非冗余Unigene序列,其中28 476个Unigene可得到功能注释;在转录本中能够被注释到GO分类的Unigene有11 201条(30.3%),KEGG分析发现共有17 172个基因参与了127个代谢通路。(2)经KEGG通路筛选出14条与脂肪酸代谢相关的通路,且随着茶籽的发育,大部分脂肪酸调控途径相关基因呈下调趋势,其中上调基因数最多的有α-亚麻酸代谢途径和脂肪酸降解途径(有17个基因表达量上调),下调基因数最多的是甘油磷脂代谢途径(有58个基因表达量下调);在茶籽发育幼果期α-亚麻酸代谢途径中表达量上调的基因数超过表达量下调的基因数。(3)研究发现茶籽脂肪酸合成相关的基因涉及14个脂类调控途径,共409条差异基因;随着茶树种子发育到成熟期,上调的差异表达基因数量在减少,下调的差异表达基因数量增加,其中α-亚麻酸途径中的基因PLA2G16、DAD1、pldA、FabF、FabI表达量上调显著,随后表达量下调。(4)qRT-PCR检测结果表明,7个茶树FAD和1个ACP差异表达基因的水平与转录组测序结果基本一致;随着茶籽的发育,基因CsFAD7和Δ6-CsFAD从幼果期、果实膨大期至果实成熟期都为差异下调表达,CsFAD2、CsFAD6和Δ7-CsFAD为差异上调表达,CsFAD8、Δ8-CsFAD和CsACP在幼果期至果实膨大期差异上调表达,在果实膨大期至果实成熟期差异下调表达。     

黄曲霉分生孢子色素合成基因pks1的鉴定及其对生长发育和侵染性的影响

【目的】分生孢子色素是真菌细胞壁的重要成分,对真菌的生长发育极为重要,并有助于真菌抵御各种环境胁迫。本研究鉴定了黄曲霉分生孢子色素合成基因,并研究了分生孢子色素对黄曲霉生长发育及其对抗紫外照射和侵染能力的影响。【方法】通过已知真菌孢子色素合成基因蛋白序列同源比对确定了黄曲霉分生孢子色素合成基因及其所在的基因簇,利用同源重组策略对目标基因进行敲除,获得了该色素合成基因缺失的突变菌株,并研究该基因敲除后对表型、产孢、菌核形成、黄曲霉毒素产生、抗紫外照射和侵染性等影响。【结果】与野生型菌株相比,黄曲霉pks1基因缺失菌株的分生孢子颜色变为白色,生长速度、孢子产量、菌核形成和黄曲霉毒素B_1的产生均没有显著性变化,但该基因的缺失导致孢子对紫外线照射的抵御能力明显减弱,降低了黄曲霉对玉米和花生种子的侵染能力。【结论】pks1(AFLA_006170)基因是黄曲霉分生孢子色素合成的关键基因,影响黄曲霉分生孢子对紫外线照射等不利环境因子的抵抗能力和对粮食种子的侵染能力。     

草地贪夜蛾取食Cry1Ab和Cry1Fa蛋白后中肠转录组及ABC基因表达分析

【目的】为揭示草地贪夜蛾Spodoptera furgiperda幼虫取食Bt蛋白后与中肠上相关ATP结合盒转运子(ATP-binding cassette transporter, ABC)蛋白基因的表达量变化的关系。【方法】分别使用含活化晶体蛋白Cry1Ab (LC₇₀=240.2 μg/g)和Cry1Fa (LC₇₀=270.0 μg/g)蛋白的人工饲料饲喂草地贪夜蛾4龄幼虫48 h,利用高通量测序对中肠进行转录组测序并进行生物信息学分析,筛选处理后差异表达基因;利用RT-qPCR验证差异表达ABC基因的表达量。【结果】与饲喂正常人工饲料的对照相比,饲喂含240.2 μg/g Cry1Ab和270.0 μg/g Cry1Fa的人工饲料后草地贪夜蛾4龄幼虫中肠转录组中分别检测到1 305和1 202个差异表达基因。Cry1Ab和Cry1Fa处理组与对照组之间分别有994和912个差异表达基因被GO功能注释到生物学过程、分子功能和细胞组分三大类。在最终筛选到的9个差异表达的ABC家族基因中,Cry1Ab处理组与对照组之间有4个差异表达ABC基因,3个上调,1个下调; Cry1Fa处理组与对照组之间有5个差异表达ABC基因,2个上调,3个下调;Cry1Ab和Cry1Fa处理组与对照组之间有2个ABC基因(LOC118267200和LOC118267201)表达量均显著上调。RT-qPCR验证结果表明,与对照组相比,Cry1Ab处理组有3个ABC基因表达量极显著上调,2个ABC基因表达量下调;Cry1Fa处理组有5个ABC基因表达量上调,1个ABC基因表达量下调。【结论】Cry1Ab和Cry1Fa蛋白的摄入可以影响草地贪夜蛾幼虫中肠一些ABC家族基因的表达量变化,这些基因的表达量变化与昆虫抗性产生有关。经比对后发现,ABCC家族与ABCG8基因表达量变化显著。本研究为下一步明确草地贪夜蛾体内ABC转运蛋白在Bt蛋白杀虫机制中的作用,以及合理使用Bt蛋白防治草地贪夜蛾及延缓抗性提供了理论依据。     

漆酶同工酶基因在斑玉蕈生长发育过程中的表达分析

为了探明漆酶在斑玉蕈生长发育过程中的功能,对斑玉蕈转录测序预测的13个漆酶基因序列进行分析、鉴定和构建分子系统发育树;检测了不同生长发育时期漆酶的活性和漆酶基因表达水平。研究结果显示：13个基因片段中有10个是漆酶基因。不同的漆酶同工酶之间进化关系存在明显差异,大多数漆酶与木腐菌（金针菇Flammulina filiformis和侧耳属Pleurotus）进化关系较近。对斑玉蕈不同生长发育时期的酶活检测结果显示,从斑玉蕈的菌丝恢复期到钉头期,漆酶活性逐渐升高,而在子实体形成后期酶活逐渐降低。对培养40d、60d和80d的菌丝样品以及不同生长发育时期的样品进RT-qPCR检测,结果显示在菌丝营养生长时期,大多数漆酶基因在第40-60天表达量持续增加1-3倍,而在第60-80天时表达量出现降低的情况。而在生殖生长时期,大多数漆酶基因在转色期或者原基期相对表达量达到最大值,并在子实体期出现降低,这与漆酶活性的检测具有一致性。lcc3、lcc7、lcc8和lcc9在斑玉蕈生殖生长过程中相对表达量出现了10-100倍的上调。这说明从菌丝培养到菌丝扭结形成子实体和子实体发育的过程中,不同的漆酶可能发挥着不同的作用,表达量较高的漆酶基因可能对基质降解和子实体形成起主要作用。     

猪苓菌丝形成菌核差异蛋白质组学研究

药用真菌猪苓菌丝形成菌核分子机制一直是人们关注的热点问题。本研究首次对猪苓进行了系统的蛋白质组学研究。基于液质联用技术并在数据依赖性采集模式下,成功鉴定了猪苓菌丝形成菌核的初始期、菌核生长期和成熟期的蛋白质组,共鉴定蛋白1 391个(global FDR 1%)。采用先进的SWATH MSALL的非标记定量方式对初始期菌丝和菌核进行了差异蛋白质组学研究,从中定量蛋白1 234个,占鉴定蛋白数量的88.7%。猪苓菌丝形成菌核过程中有378个蛋白差异表达。差异蛋白GO注释结果表明,在猪苓菌丝形成菌核过程中,差异蛋白参与的蛋白功能和生物学过程多样,以催化和结合为主,还参与感知环境刺激、信号转导、电子转移、过氧化活性、菌丝生长及胞壁合成等功能。进一步推测分析表明,猪苓菌丝形成菌核可能是由于缺氧或温差等多种环境胁迫诱导所致,生成活性氧并形成氧化应激状态。信号通过丝裂原活化蛋白激酶通路(MAPK通路)途径向下游传递,可能调控了胞壁蛋白糖基化修饰,进而促进猪苓菌丝极性生长和菌丝形态学变化而形成菌核。KEGG代谢通路分析表明,差异蛋白还主要参与了糖、脂质和氨基酸等初级产物以及甾醇、多聚乙酰等次级代谢物的合成与代谢。因此,对差异蛋白质组的深入解析和研究有利于揭示猪苓菌丝形成菌核的蛋白分子机制,具有较强的理论和现实意义。     

果生刺盘孢CfHAC1调控应答二硫苏糖醇胁迫的转录组分析

果生刺盘孢Colletotrichum fructicola是油茶炭疽病优势病原菌。前期研究发现bZIP转录因子CfHac1参与调控该菌的生长发育和致病性。为了揭示转录因子CfHac1调控果生刺盘孢响应内质网压力和致病机理,本研究测定了ΔCfhac1突变体对内质网压力胁迫剂的敏感性,发现突变体对二硫苏糖醇（dithiothreitol,DTT）的耐受性下降,说明CfHAC1基因可能参与调控果生刺盘孢响应内质网压力胁迫过程。进一步利用高通量RNA-seq技术对该病菌野生型菌株和CfHAC1敲除突变体菌株在DTT胁迫下的转录组进行了比较分析,结果表明差异表达基因共有2 680个,其中上调表达基因有1 181个,下调表达基因有1 499个。Gene Ontology 功能分析结果显示,差异表达基因主要参与催化活性、结合、代谢过程、细胞过程、细胞成分合成、生物过程调控和应激反应等生物学过程。KEGG功能富集分析表明,上调表达基因主要被富集到核糖体、真核细胞的核糖体生物合成、RNA转运和氰基氨基酸代谢通路中;下调表达基因显著富集在内质网蛋白质加工、N-聚糖生物合成、类固醇合成和蛋白质分泌等通路中。分析发现转录因子CfHac1调控内质网胁迫应答和致病相关基因的表达。本研究提供了在全基因组水平上对CfHAC1基因与内质网压力胁迫应答之间关联的新认识,为阐明果生刺盘孢响应内质网压力胁迫和致病机制奠定了基础。     

花生下胚轴曝光与不曝光差异表达基因的数字表达谱分析

为研究花生播种后通过覆土延长幼苗曝光,可促使下胚轴伸长,引伸子叶节出土,控制早花下针,有利于提高产量的机理,本研究采用数字基因表达谱的方法分析了花生幼苗曝光(处理)和不曝光(对照)下胚轴中的基因表达差异。结果表明,对照中检测到40 275个基因,处理中检测到40 554个基因;处理与对照相比,共检测出2 722个表达差异基因,其中1 887个基因上调表达,835个基因下调表达;根据GO分析标准,2 722个表达差异基因可以分为3个本体类别。生物过程本体有22类功能的基因,其中代谢过程类有关的基因最多,其次是细胞过程类有关的基因;细胞组成本体有14类功能的基因,其中细胞和细胞组分类有关的基因最多,其次是细胞器类有关的基因;分子功能本体有12类功能的基因,其中催化类基因最多,其次是结合类基因。进一步对差异表达基因富积的植物激素代谢路径分析表明,生长素代谢路径中有6个差异表达基因,其中3个基因上调表达,3个基因下调表达;细胞分裂素代谢路径中有3个差异表达基因,全部为下调表达;赤霉素代谢路径中有3个差异表达基因,其中2个上调表达,1个下调表达;脱落酸代谢路径中有4个差异表达基因,其中3个表达上调,1个表达下调;乙烯代谢路径中有2个差异表达基因,1个上调表达,1个下调表达;茉莉酸代谢路径有9个差异表达基因,2个表达上调,7个表达下调;另外,油菜素甾醇代谢路径中有一个表达上调差异基因。本研究结果对了解延长幼苗曝光,促进下胚轴生长,控制下针有利于花生高产的机理具有重要意义。     

基于iTRAQ技术分析模拟空间诱变酵母菌株的差异蛋白表达

目的：探究模拟空间诱变处理后,生长速率提高的酿酒酵母细胞的诱变机制。方法：采用同位素标记相对和绝对定量（iTRAQ）技术对模拟空间诱变处理前后的酿酒酵母进行定量蛋白质组学研究,探讨对照组与诱变组在蛋白质组表达水平上的差异,以及受影响的生物功能。结果：共鉴定到764个差异蛋白,占检测蛋白总数的15.54%。其中,上调差异表达蛋白378个、下调差异表达蛋白386个,主要涉及细胞核、线粒体、核糖体和质膜等重要细胞器和细胞结构。结论：提高重要代谢的效率、降低DNA修复功能与环境防御能力,可能是模拟空间诱变提高酿酒酵母细胞生长速率的重要机制之一。     

UVA对白念珠菌生物膜作用的转录组学研究

目的应用转录组技术,对UVA照射后的白念珠菌菌丝及生物膜进行基因表达分析。方法体外构建白念珠菌(ATCC90028)生物膜,采用UVA-LED抗菌设备(复旦大学光电所制,波长365nm,强度13mW/cm2,频率100 Hz),照射白念珠菌生物膜5min、10min、15min,对照组不光照。用转录组技术分析UVA实验组和对照组之间差异表达的基因,进行GO、KEGG遗传差异富集分析。结果 UVA光照组相对于对照组有大量差异表达基因,核糖体合成RIO2基因、氨基糖和核苷糖代谢相关基因MCR1、CBR1、内质网蛋白加工途径相关基因SHP1显著上调;RNA降解途径相关基因MPP6和细胞内吞作用相关基因CHMP6显著下调。结论 UVA直接或间接调节白念珠菌菌丝及生物膜相关生理代谢途径,参与诱发细胞凋亡过程。     

感染球孢白僵菌的家蚕免疫应答差异表达基因分析

【目的】筛选家蚕Bombyx mori应对白僵菌Beauveria bassiana侵染的应答基因, 以进一步研究家蚕抵御真菌侵染的分子机制。【方法】采用新一代Solexa高通量测序技术对感染白僵菌及未感染白僵菌的对照组家蚕进行了测序分析, 筛选差异表达基因; 结合生物信息学工具分析差异表达基因的功能注释、 分类及涉及的信号通路等; 应用荧光定量PCR技术验证10个基因的差异表达。【结果】通过测序和生物信息学分析共获得377个差异表达基因, 其中表达上调基因236个, 下调基因141个; KEGG通路分析表明, 各通路中既有表达上调的基因, 也有下调基因; 12个上调基因、 26个下调基因参与3个显著性富集的KEGG通路, 即核糖体、 氨酰tRNA生物合成和剪接体通路。定量PCR与测序结果显示, 溶菌酶、 热激蛋白、 谷胱甘肽S-转移酶、 肽聚糖识别蛋白等与免疫应激相关的蛋白基因均呈现表达上调。【结论】本研究筛选获得的差异表达基因, 特别是上调表达的基因可能与家蚕应对白僵菌侵染的应答机制有关, 其中与免疫应激相关的蛋白基因如溶菌酶、 热激蛋白、 谷胱甘肽S 转移酶、 肽聚糖识别蛋白基因等可能直接参与了家蚕对白僵菌的免疫识别和防御, 研究结果为从分子水平阐明家蚕抵御真菌侵染的防御机制和白僵菌对家蚕的致病机理提供新的依据。     

信息素信号通路基因在斑玉蕈生长发育过程中的差异表达

为了更清楚地了解斑玉蕈菌丝成熟、原基形成和子实体发育的过程,本研究对不同菌丝培养时期的栽培瓶进行出菇实验,并对其不同培养时期和生长发育关键时期的信息素通路基因进行差异表达分析,以期揭示信息素信号通路基因参与调节斑玉蕈菌丝的生长、子实体形成和发育的作用。研究结果表明：斑玉蕈菌丝培养40-80d过程中,子实体产量呈上升的趋势,说明菌丝的成熟程度对产量会产生重要影响。对斑玉蕈基因组中的信息素信号通路基因进行分析鉴定共获得了8个关键基因。信息素通路基因差异表达分析表明：在菌丝培养40-80d过程中,大部分信息素信号通路基因在第60天时表达量最高,其中ste20、cdc24和ste12上调了4-20倍,而在第80天出现下降。从菌丝恢复到扭结形成原基和子实体发育的过程中,大多数基因在原基时期表达量最高,其中ste20、cdc24、ste11和ste12表达量上调最为显著,在子实体成熟期这些基因表达量下降。因此,这说明在菌丝营养生长过程中,在第60天菌丝细胞增殖生长最为旺盛,而在第80天菌丝细胞基本停止生长,菌丝也逐渐达到成熟。同时,在菌丝生殖生长过程中,斑玉蕈持续地上调信息素通路基因表达使菌丝细胞不断地分裂增殖,从而使新生的菌丝扭结形成原基,其中ste3、ste20、cdc24、ste11和ste12基因可能对斑玉蕈菌丝细胞的分裂增殖和诱导子实体形成起到关键的作用。     

苜蓿根瘤菌自生与共生状态下核糖体蛋白的比较

苜蓿根瘤菌在与宿主植物建立共生关系的过程中,以自生状态进入宿主植物细胞,经过分化发育转变为共生状态的类菌体(Bacteroid)。由于生存环境发生了变化,类菌体在形态、结构和功能方面都产生了很大的改变,其中最为明显的改变是类菌体获得了共生固氮的能力。此时,类菌体中许多与共生相关的基因被激活,蛋白的表达量显著增加。为了探明这种改变是否与合成蛋白质的细胞器-核糖体有关,比较分析了苜蓿根瘤菌在自生和共生状态下核糖体蛋白的表达谱。蛋白质双向电泳结果显示二者之间没有明显的差别,说明类菌体的分化发育过程中核糖体蛋白的形成没有改变。     

转录组分析揭示东方蜜蜂微孢子虫侵染意大利蜜蜂的分子机制

【目的】本研究旨在通过差异表达基因(differentially expressed gene, DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log_2(Fold change)|≥1的标准,通过比较分析筛选出NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组的DEG。通过相关生物信息学软件对上述DEG进行Venn分析、GO分类和KEGG代谢通路富集分析。根据Nr和KEGG数据库注释信息和相关文献进行对东方蜜蜂微孢子虫的毒力因子和侵染相关因子的统计和分析。通过RT-qPCR验证转录组数据及DEG表达趋势。【结果】从NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组分别鉴定出1 397, 1 497和52个DEG。Venn分析结果显示各比较组共有的上调和下调基因分别为10和1个。GO分类结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG富集数最多的功能条目为代谢进程、细胞进程、单组织进程、细胞、细胞组件、细胞器、催化活性和结合,而NcT1 vs NcT2中DEG富集数最多的是代谢进程、细胞进程、单组织进程、催化活性和结合。KEGG代谢通路富集分析结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG分别富集到80和79条通路;富集在糖酵解/糖异生和MAPK信号通路的上调基因数量多于下调基因。毒力因子分析结果显示,孢壁蛋白9基因和孢壁蛋白12基因在NcCK vs NcT1和NcCK vs NcT2中均下调表达,孢壁蛋白8基因仅在NcCK vs NcT1中表达量下调;此外孢壁蛋白前体基因、孢壁和锚定盘复合蛋白基因、几丁质合酶基因、极管蛋白基因、蓖麻毒素B凝集素基因的表达水平在NcCK vs NcT1和NcCK vs NcT2中表现为上调。侵染相关因子分析结果表明,糖酵解途径的3个关键酶基因在NcCK vs NcT1和NcCK vs NcT2中上调表达;3个涉及ATP/ADP移位酶的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,但有1个表达量下调;2个涉及ABC转运蛋白的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,另有4个下调表达。RT-qPCR结果证实了本研究中转录组数据及DEG表达趋势的真实可靠性。【结论】本研究通过比较分析解析东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程的转录组动态,揭示了孢壁蛋白、孢壁和锚定盘复合蛋白、几丁质酶、极管蛋白和蓖麻毒素B凝集素等毒力因子编码基因,以及己糖激酶、丙酮酸激酶、6-磷酸果糖激酶、ATP/ADP移位酶和ABC转运蛋白等侵染相关因子编码基因在病原增殖中扮演重要角色,为阐明东方蜜蜂微孢子虫的侵染机制提供了基础。     

小麦条锈菌PsCdc2基因的克隆及在条锈菌侵染小麦后的转录表达分析

【目的】克隆小麦条锈菌细胞分裂基因PsCdc2,分析该基因在条锈菌接种小麦后不同时间点的表达特征。【方法】利用PCR和RT-PCR技术克隆PsCdc2的cDNA序列和基因组序列,采用生物信息学技术预测分析该基因编码蛋白的保守结构域及基本特性,对该蛋白进行系统发育分析,构建进化树;运用实时荧光定量RT-PCR技术,以PsCdc2在夏孢子时期的表达情况为对照,分析该基因在亲和及非亲和互作中不同时间点的表达特征。【结果】PsCdc2基因组序列长2279 bp,由11个外显子和10个内含子构成,开放阅读框为885 bp,编码294个氨基酸,分子量为33.14 kDa,等电点为6.26。编码蛋白含两个保守的激酶特征位点,一个跨膜螺旋区域。PsCdc2基因编码蛋白与小麦秆锈菌、新型隐球菌、玉米瘤黑粉菌等多种真菌的Cdc2高度相似,其中与小麦秆锈菌的Cdc2亲缘关系最近,序列相似性达73.1%。实时荧光定量RT-PCR结果表明,在亲和组合中,该基因在条锈菌接种小麦的前期上调表达,其中接种后12 h时表达量最高,约为夏孢子中表达量的1.62倍,接种后24-268 h,基因表达基本呈下调趋势,其中96 h基因表达量最低,仅为夏孢子时期的0.07倍。在非亲和组合中,该基因表达基本呈下调趋势,在接种后各个时间点的表达量均低于在夏孢子中的表达量,其中接种后12 h时表达量最高,但仅为夏孢子中表达量的0.34倍;接种后96 h表达量最低,为夏孢子中表达量的0.02倍。【结论】PsCdc2可能通过调控条锈菌的细胞周期循环参与了侵染前期初生菌丝生长和吸器母细胞的形成,与条锈菌的致病性相关。本文首次报道了小麦条锈菌的Cdc2基因,为进一步揭示条锈菌细胞周期调控的本质及研究开发靶向Cdc2的新型农药,以及实现对小麦条锈病的新型药剂防治提供了理论基础。