中华蜜蜂幼虫肠道响应球囊菌早期胁迫的转录组学

【目的】白垩病是困扰养蜂生产的顽疾。目前,尚无利用二代测序技术研究中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫白垩病的报道。本研究利用RNA-seq技术对健康(Ac CK)及球囊菌(Ascosphaera apis)胁迫的中蜂4日龄幼虫肠道(Ac T)进行深度测序,在转录组水平研究中蜂幼虫在球囊菌胁迫早期的胁迫应答。【方法】通过Illumina Hi Seq 2500平台对Ac CK和Ac T进行双端(PE125)测序,首先对测序数据进行质控和评估,利用edge R软件进行差异表达基因(DEG)分析,进而对DEGs进行GO富集分析及KEGG代谢通路富集分析,最后,利用实时荧光定量PCR(q RT-PCR)验证测序数据的可靠性。【结果】Ac CK和Ac T的转录组测序共得到188 457 338条原始读段(raw reads),经过滤得到182 088 448条有效读段(clean reads),两端Q20与Q30均在97.96%和94.97%以上,说明测序数据质量良好;主成分分析(PCA)结果显示第一与第二主成分可分别解释基因表达总体差异的75.8%和10.7%;DEG分析结果显示上调基因与下调基因的数量分别为344和239个;GO富集分析结果显示DEGs共富集在36个GO分类(term)上,其中基因富集数最多的细胞(106 unigenes)、细胞组件(106 unigenes)和代谢进程(104 unigenes);KEGG代谢通路富集分析结果显示上调与下调基因分别富集在72个和45个代谢通路上,其中上调基因富集数最多的是核糖体(72 unigenes)、碳代谢(16 unigenes)和糖酵解(14 unigenes),而下调基因富集数最多的是碳代谢(9 unigenes)、二羧酸代谢(8 unigenes)和氨基酸生物合成(7 unigenes),进一步分析表明中蜂幼虫肠道的部分细胞免疫对球囊菌的胁迫产生应答,而体液免疫不产生应答,宿主的代谢相关基因受到球囊菌的显著抑制。【结论】揭示了中蜂幼虫在球囊菌入侵早期的胁迫应答,为深入解析中蜂幼虫的胁迫应答机制提供了重要信息,也为在分子水平研究关键应答基因打下了基础。     

胁迫中华蜜蜂幼虫肠道的球囊菌及其体外培养的高表达基因分析

利用RNA-seq技术对胁迫中华蜜蜂(简称中蜂)6日龄幼虫肠道的球囊菌及其体外培养进行深度测序,根据基因的FPKM值筛选得到高表达基因(HEGs),进而对HEGs进行GO及KEGG代谢通路(pathway)富集分析.Illumina测序数据经质控和过滤得到100814558条有效读段(clean reads),平均Q30均在93.81%以上.GO富集分析结果显示处理组(Aac T)的HEGs富集于39个GO term,基因富集数最多的是细胞(1872 unigene)、细胞组件(1872 unigene)和细胞进程(1748 unigenes);对照组(Aa CK)的HEGs富集于37个GO term,基因富集数最多的是细胞(785 unigenes),其次是细胞组件(785 unigenes)和代谢进程(776 unigenes).KEGG pathway富集分析显示,Aac T的HEGs富集在119个代谢通路上,基因富集数最多的是核糖体(176 unigenes)、碳代谢(147 unigenes)及氨基酸的生物合成(134 unigenes);Aa CK的HEGs富集在110个代谢通路上,基因富集数最多的是核糖体(179 unigenes)、氨基酸的生物合成(70 unigenes)以及碳代谢(62 unigenes).深入分析发现,对于富集在MAPK信号通路上的高表达基因数量,胁迫中蜂幼虫肠道的球囊菌远多于体外培养的球囊菌,说明病原的该通路在胁迫后期被显著激活.     

意大利蜜蜂幼虫响应东方蜜蜂微孢子虫胁迫的高表达基因分析

为探究意大利蜜蜂(Apis mellifera ligustica,简称意蜂)幼虫响应东方蜜蜂微孢子虫(Nosema ceranae)胁迫的高表达基因(HEG)的作用,利用RNA-seq技术对正常的意蜂4、5和6日龄幼虫(Am4CK、Am5CK和Am6CK)及N.ceranae胁迫的意蜂4、5和6日龄幼虫(Am4T、Am5T和Am6T)进行测序。共得到268 765 544条原始读段,经过滤得到的有效读段数为266 318 442条,各样品的Q30均在92.42%及以上。根据FPKM值≥15的标准从上述6个样品中分别筛选出2 460、2 692、2 908、2 367、2 549和2 270个HEG。Venn分析结果显示,去除与对照组共有的HEG,Am4T、Am5T和Am6T的共有HEG数为3个,特有HEG数分别为373、210和79个。GO分类结果显示,Am4T、Am5T和Am6T的特有HEG分别能够注释26、27和16个GO条目,注释基因数最多的分别是代谢进程、结合和结合。KEGG分析结果显示这些特有HEG可分别富集在145、138和68条代谢通路,富集基因数最多的分别是过氧化物酶体、碳代谢和内质网蛋白加工。进一步分析结果表明N.ceranae胁迫可导致宿主的部分物质代谢、能量代谢和免疫代谢通路水平提高。从6个样品的共有HEG中随机选取10个进行RT-PCR验证,均能成功扩增出目的条带,说明本研究中的HEG真实存在。研究结果揭示了意蜂幼虫响应N.ceranae胁迫的HEG表达谱及作用,为深入解析宿主的胁迫应答机制提供了重要信息。     

基于比较转录组学分析揭示中华蜜蜂及意大利蜜蜂幼虫的球囊菌抗性差异机制

为探究中华蜜蜂(Apis cerana cerana,以下简称中蜂)与意大利蜜蜂(Apis mellifera ligustica,以下简称意蜂)幼虫的球囊菌(Ascospahaera apis)抗性差异产生的原因,利用Venn分析、GO分类分析、KEGG代谢通路分析以及Ka/Ks分析对二者的转录组进行比较研究。Venn分析结果显示,中蜂与意蜂的同源基因有17 656个,归属于8 111个基因家族,二者的特有基因分别有1 158和241个,分别归属于468和86个基因家族。GO分类结果显示,中蜂和意蜂的特有基因分别富集在38和28个GO条目。KEGG代谢通路富集分析结果显示,中蜂和意蜂的特有基因分别富集于96和21个代谢通路;进一步分析发现中蜂的特有基因富集在内吞作用、溶酶体、泛素介导的蛋白水解和黑化作用等4个细胞免疫通路,以及MAPK信号通路和Toll-like受体信号通路等2个体液免疫通路,而意蜂仅有1个特有基因富集在MAPK信号通路。Ka/Ks分析结果显示受到强烈正向选择、弱正向选择、中性选择和负向选择的单拷贝同源基因分别有37、281、2 630和3 269个。对受到强烈正向选择的基因进行GO分类和KEGG代谢通路富集分析,GO结果显示中蜂和意蜂的单拷贝同源基因富集的GO条目类型相同;KEGG结果显示中蜂的单拷贝基因仅富集在氧化磷酸化。上述结果表明免疫相关基因数量的差异是二者的球囊菌抗性差异产生的重要原因之一,中蜂在与球囊菌的协同进化过程中可能通过提高能量利用从而限制球囊菌的增殖。本研究结果可为阐明中蜂及意蜂幼虫的球囊菌抗性差异产生的分子机制提供参考。     

球囊菌侵染对意大利蜜蜂成年工蜂的影响

为了探究球囊菌侵染对意大利蜜蜂成年工蜂存活率、免疫基因、发育基因及肠道菌群的影响,将纯化培养的球囊菌孢子接种6日龄意大利蜜蜂成年工蜂,利用荧光定量PCR(qRT-PCR)检测球囊菌侵染对意大利蜜蜂成年工蜂免疫基因、发育基因及肠道菌群的影响。结果表明,连续3 d接种球囊菌孢子对意大利蜜蜂成年工蜂存活率无显著影响(χ²=0.133,P=0.715)。意大利蜜蜂工蜂中肠与后肠球囊菌孢子含量在接种球囊菌孢子4 d时分别显著高于接种8 d时(P<0.05)。与未接种组相比,接种球囊菌孢子4 d,意大利蜜蜂成年工蜂免疫相关基因abaecin、apidaecin、defensin-1、defensin-2、hymenoptaecin、GNBP-1和IMD均显著上调表达(P<0.05);接种球囊菌孢子8 d,意大利蜜蜂成年工蜂免疫相关基因相对表达量无显著差异(P>0.05);接种球囊菌孢子4 d和8 d,发育相关基因VGMC、usp、kr-h1、AMHex10869相对表达量无显著差异(P>0.05)。与未接种组相比,接种球囊菌孢子4 d,意大利蜜蜂成年工...     

蜜蜂球囊菌分泌多种胞外酶侵染蜜蜂幼虫

对健康和患病蜜蜂幼虫血淋巴进行SDS-PAGE电泳和蛋白酶、酯酶的活性染色,以研究蜜蜂球囊菌分泌的胞外酶种类及对蜜蜂侵染的作用.结果表明:健康蜜蜂幼虫血淋巴中的蛋白含量丰富,主要由4种高分子质量的蛋白组成,而患病幼虫血淋巴中的蛋白含量很少,主要蛋白组分被降解;在患病幼虫血淋巴中可以检测到健康幼虫所没有的多种蛋白酶和酯酶...     

基于RNA-seq技术分析意大利蜜蜂幼虫肠道响应球囊菌胁迫的高表达基因

对意蜂幼虫肠道组织基因进行深度RNA测序,并对高表达基因进行GO注释和KEGG分析,旨在揭示意蜂幼虫免疫抗病机制。分别取3、4、5、6日龄意蜂幼虫中肠组织抽提RNA,合成cDNA。得到的总读段数都在26 715 638条以上,有效数据达到97%以上,且测序饱和度良好,说明测序所得数据真实可信。将测序有效数据和意蜂公布DNA数据库进行比对,发现意蜂幼虫肠道中有13 789个基因有不同程度表达。选取各样本FPKM排序最前且表达水平最高的100个基因进行GO注释,这里主要涉及细胞组分、分子功能和生物学进程功能分析。利用KEGG分析显示这100个高表达基因分布在能量代谢、翻译、组装和降解以及信号转导等多条通路上。初步分析相应共有基因和特有基因在相关通路上的作用,并对意蜂幼虫肠道受球囊菌胁迫后的相关免疫机制进行一定分析。     

基于RNA-seq数据分析玉米抗虫响应基因的可变剪接事件

采用Illumina Hiseq~(TM)分别对黏虫取食的玉米叶片及对照组材料进行转录组测序,鉴定和分析响应黏虫取食基因的可变剪接事件。结果表明,对照组中鉴定出15 701个基因对应79 363个可变剪接事件,黏虫取食组中鉴定到11 791个基因的39 385个可变剪接事件。2组玉米基因组在不同的可变剪接类型中,均是以第1个外显子可变剪切、最后1个外显子可变剪切、单内含子滞留和可变5′或3′端剪切4种类型为主。对2组发生可变剪接事件的基因进行比对发现,黏虫处理组新增加了1 121个可变剪接基因,差异表达分析鉴定出177个表达显著差异的可变剪接基因。GO功能富集分析表明,发生可变剪接的差异基因主要富集于氧化还原酶活性、转移酶活性、DNA结合、ATP结合、代谢过程、以DNA为模板的转录及调控相关的功能,表明这些可变剪接差异表达基因参与了玉米的抗虫响应过程。     

意大利蜜蜂幼虫肠道在球囊菌侵染前期的 差异表达microRNA及其调控网络

【目的】 微小RNA（microRNA,miRNA）是一类重要的基因表达调控因子,可影响宿主与病原间的互作过程。蜜蜂球囊菌（Ascosphaera apis）是一种特异性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在对意大利蜜蜂（Apis mellifera ligustica,简称意蜂）幼虫肠道在球囊菌胁迫前期的差异表达miRNA（differentially expressed miRNA,DEmiRNA）及其靶基因进行深入分析,在miRNA组学水平探究意蜂幼虫在球囊菌侵染前期的胁迫应答,并通过构建显著DEmiRNA的调控网络筛选出与宿主应答相关的关键miRNA。【方法】 利用small RNA-seq（sRNA-seq）技术对正常及球囊菌侵染的意蜂4日龄幼虫肠道（AmCK和AmT）进行高通量测序,首先对原始数据进行质控和评估,随后将过滤后的数据与西方蜜蜂（Apis mellifera）参考基因组进行比对;将比对上的序列标签（tags）注释到miRBase数据库,得出已知miRNA的表达量;通过TPM（tags per million）算法对各样本中miRNA的表达量进行归一化处理,以|log₂ fold change|≥1且P≤0.05作为标准筛选得到显著DEmiRNA;利用TargetFinder 软件预测显著DEmiRNA的靶基因,并对其进行GO和KEGG代谢通路（pathway）富集分析。利用Cytoscape软件对miRNA-mRNA调控网络进行可视化。最后,利用茎环反转录PCR（Stem-loop RT-PCR）和荧光定量PCR（qPCR）验证测序数据的可靠性。【结果】 AmCK和AmT样品的测序分别得到13 553 302和10 777 534条原始读段（raw reads）,经严格过滤后得到的有效读段（clean reads）数分别为13 186 921和10 480 913条。各样品的生物学重复间的Pearson相关性系数分别在0.9822和0.9508以上。共有10个显著DEmiRNA,包括4个上调miRNA和6个下调miRNA。显著DEmiRNA在AmT的整体表达水平低于AmCK。10个显著DEmiRNA可靶向结合3 788个靶基因,其中上调miRNA的1 240个靶基因可注释到GO数据库中的39个GO条目,主要富集在结合、细胞进程、代谢进程和应激反应等;下调miRNA的749个靶基因可注释到34个GO条目,主要富集在细胞进程、结合、代谢进程和应激反应等。KEGG数据库注释结果显示,上调miRNA和下调miRNA的靶基因分别注释到95和66条代谢通路,富集基因数最多的分别是Wnt信号通路、Hippo信号通路、光传导以及内吞作用、磷脂酰肌醇信号系统、嘌呤代谢。对于上调和下调miRNA,分别有31和52个靶基因注释到内吞作用,15和7个靶基因注释到泛素介导的蛋白水解,11和5个靶基因注释到Jak-STAT信号通路,1和3个靶基因注释到MAPK信号通路。显著DEmiRNA与靶mRNA之间形成复杂的调控网络,7个显著DEmiRNA靶向结合96个与Wnt信号通路相关的mRNA,8个显著DEmiRNA靶向结合55个与内吞作用相关的mRNA。Stem-loop RT-PCR和qPCR结果验证了测序数据的可靠性。【结论】 对意蜂幼虫肠道在球囊菌侵染前期的DEmiRNA及其靶基因进行预测和分析,并构建和分析了DEmiRNA-mRNA调控网络,研究结果提供了宿主miRNA的表达谱和差异表达信息,揭示了DEmiRNA通过调控细胞生命活动、新陈代谢以及部分细胞和体液免疫等生物学过程参与宿主的胁迫应答。miR-4331-y、miR-4968-y、miR-8440-y、novel-m0023-5p和novel-m0025-3p共同参与了宿主的Wnt信号通路和内吞作用的调控,可作为白垩病治疗的潜在分子靶标。     

玉米ZmcpSRP54基因可变剪接分析

cpSRP54基因在植物叶绿体发育中发挥重要作用。本研究用生物信息学方法对玉米ZmcpSRP54蛋白进行了进化分析,利用RT-PCR方法对玉米ZmcpSRP54基因全长cDNA进行了克隆,进而对其进行可变剪接分析。结果表明,玉米ZmcpSRP54蛋白与其它7个物种的cpSRP54蛋白有很高的相似性。鉴定到玉米ZmcpSRP54基因25个不同转录本,1个为标准剪接转录本,其它24个为可变剪接转录本。24个可变剪接转录本共使用了17种非标准剪接位点,4种可变剪接方式,主要以可变的3′端位点、可变的5′端位点和外显子跳跃3种可变剪接方式为主。预测出18种不同长度的ORF,编码18种不同长度的多肽。11个多肽其保守结构域氨基酸序列发生部分缺失,6个多肽其保守结构域氨基酸序列发生全部缺失,1个编码全新的多肽。这些结果将为玉米ZmcpSRP54基因的功能研究提供重要信息。     

蜜蜂球囊菌基因结构优化及新基因鉴定

为补充和完善蜜蜂球囊菌Ascosphaera apis基因组序列和功能注释信息,本研究利用RNA-seq技术对球囊菌菌丝和孢子进行深度测序,并基于高质量的转录组数据对已注释基因进行结构优化,对未注释基因进行预测、鉴定和分析。将测序得到的有效读段进行参考基因组比对和转录本重构;利用Cuffcompare软件将重构转录本与参考基因组进行比对。结果表明:共对101个已注释基因的5′端或3′端进行了延长;共鉴定出373个新基因,随机挑选10个新基因进行RT-PCR验证,其中8个能扩增出符合预期的目的片段,表明预测出的多数新基因真实存在;共计147个球囊菌新基因可注释到Nr和eggNOG数据库中,85个新基因注释到29个GO条目,66个新基因注释到33条代谢通路。本研究结果补充了球囊菌的基因结构和功能注释信息,也为新基因的功能研究打下了初步基础。     

意大利蜜蜂响应东方蜜蜂微孢子虫胁迫的免疫应答

【目的】通过对意大利蜜蜂（Apis mellifera ligustica,简称意蜂）工蜂中肠响应东方蜜蜂微孢子虫（Nosema ceranae）胁迫的差异表达基因（differentially expressed gene,DEG）及免疫通路进行深入细致的分析,揭示宿主的细胞和体液免疫应答,为关键免疫应答基因的筛选及功能研究打下基础。【方法】基于前期获得的正常及东方蜜蜂微孢子虫胁迫的意蜂7日龄和10日龄工蜂中肠（Am7CK、Am7T、Am10CK、Am10T）转录组数据,按照FDR≤1,P≤0.05和|log₂ fold change|≥1的标准筛选出各组的DEG,利用相关生物信息学软件对DEG进行Pearson相关性分析、Venn分析、GO分类和KEGG代谢通路富集分析,进而对免疫通路富集基因进行统计和分析,利用实时荧光定量PCR（real-time quantitative PCR,RT-qPCR）验证转录组数据的可靠性。【结果】 差异分析结果显示,Am7CK vs Am7T比较组包含472个上调基因和385个下调基因;Am10CK vs Am10T比较组包含611个上调基因和360个下调基因。Venn分析结果显示,两个比较组特有的DEG分别为739和853个,共有的DEG为118个。GO分类结果显示,Am7CK vs Am7T比较组中上调和下调基因分别涉及23和29个功能条目,其中富集上调基因数最多的前5位分别为结合、催化活性、代谢进程、细胞进程和单组织进程;富集下调基因数最多的前5位分别为代谢进程、单组织进程、催化活性、细胞进程和结合。Am10CK vs Am10T比较组中上调和下调基因分别涉及36和26个功能条目,其中富集上调基因数最多的前5位分别为单组织进程、结合、细胞进程、催化活性和代谢活性;富集下调基因数最多的前5位分别为结合、细胞进程、催化活性、代谢进程和单组织进程。KEGG代谢通路富集分析结果显示,Am7CK vs Am7T比较组中上调和下调基因分别富集在38和33条代谢通路,富集上调基因数最多的前5条分别为胆汁分泌、内质网蛋白加工、泛素介导的蛋白水解、PI3K-Akt信号通路和神经营养因子信号通路;富集下调基因数最多的前5条分别为胞质DNA传感通路、嘌呤代谢、嘧啶代谢、RNA聚合酶和核糖体;涉及泛素介导的蛋白水解等3条细胞免疫通路,以及PI3K-Akt信号通路等7条体液免疫通路。Am10CK vs Am10T比较组中上调和下调基因分别富集在54和43条代谢通路,富集上调基因数最多的前5条分别为Hippo信号通路、药物代谢-细胞色素P450、P450对外源物质的代谢、泛素介导的蛋白水解和鞘脂类代谢;富集下调基因数最多的前5条分别为mRNA监测、鞘脂类信号通路、果糖和甘露糖代谢、半乳糖代谢和鞘脂类代谢;涉及泛素介导的蛋白水解等7条细胞免疫通路,以及NF-κB信号通路等2条体液免疫通路。RT-qPCR验证结果显示6个随机挑选的DEG的表达量变化趋势与测序结果一致,证实了本研究中测序数据的可靠性。进一步分析发现意蜂7日龄和10日龄工蜂中肠的NF-κB信号通路均被东方蜜蜂微孢子虫激活,随即启动了3种抗菌肽apidaecin、defensin-1和hymenoptaecin的合成,体现了它们在宿主抵御东方蜜蜂微孢子虫入侵中的重要性。【结论】在转录组水平解析了意蜂工蜂中肠对东方蜜蜂微孢子虫的免疫应答,揭示宿主在胁迫前期同时作出细胞和体液免疫应答,前者可能在抵御病原入侵方面发挥主要作用;宿主的细胞免疫在胁迫后期持续增强,但体液免疫较大程度地减弱;泛素介导的蛋白水解通路及富集DEG,NF-κB信号通路及富集DEG,以及抗菌肽编码基因apidaecin、defensin-1和hymenoptaecin可能在意蜂工蜂对东方蜜蜂微孢子虫的免疫应答中起到关键作用。     

意大利蜜蜂工蜂中肠响应东方蜜蜂微孢子虫胁迫的可变剪接基因分析

基于前期已获得的意大利蜜蜂(Apis mellifera ligustica)工蜂中肠转录组数据,在全转录组水平上对意大利蜜蜂的可变剪接基因(alternatively spliced gene, ASG)进行鉴定和分析,旨在揭示ASG在宿主响应东方蜜蜂微孢子虫(Nosema ceranae)胁迫中的作用.在正常意大利蜜蜂7和10日龄工蜂中肠(Am7CK、Am10CK)和东方蜜蜂微孢子虫胁迫的意大利蜜蜂7和10日龄工蜂中肠(Am7T、Am10T)样品中检测到可变剪接事件数分别为73 383、85 618、79 813和74 693次,对应的基因数分别为9 586、10 063、9 829和9 622个,平均每个ASG上发生的可变剪接事件数分别为7.66、8.50、8.12和7.76次.各样品的可变剪接事件中外显子跨越的数量最多,分别为4 411、5 247、4 993和4 629个.Venn分析结果显示,Am7CK、Am10CK、Am7T和Am10T中的共有ASG数为9 056个,特有ASG数分别为100、153、275和95个.GO分类结果显示,这些共有ASG分布在50个功能条目,特有ASG分别分布于26、30、33和24个功能条目.KEGG代谢通路富集分析结果显示,上述共有ASG富集在325个代谢通路,包括内质网蛋白加工、核糖体和RNA转运等;特有ASG分别富集在31、104、126和22个代谢通路.研究结果提供了意大利蜜蜂响应东方蜜蜂微孢子虫胁迫过程中的ASG数量和种类信息,揭示宿主ASG在胁迫响应中参与了新陈代谢和细胞免疫,为关键ASG的筛选和功能研究建立了数据基础.     

金属离子对蜜蜂球囊菌几丁质酶活力的影响

探讨了不同金属离子对蜜蜂球囊菌几丁质酶活力的影响.结果表明:在0-100 mmol.L-1范围内,一价金属离子Na+、K+随离子浓度的增大,对几丁质酶活力的影响表现为先促进后抑制;在1-10 mmol.L-1范围内,二价金属离子Mg2+、Mn2+起抑制作用,Zn2+基本没有影响,Ca2+、Fe2+起促进作用;在1-10 mmol.L-1范围内,三价金属离子A l3+起促进作用,Fe3+基本没有影响.     

人类1号染色体可变剪接与普通剪接基因同义密码子的使用分析I.同义密码子偏爱使用分析

人类1号染色体可变剪接(选择性剪接)基因344非冗余蛋白质编码序列(188183密码子)和普通剪接(非可变剪接)基因的386蛋白质编码序列(223116密码子)被用于研究人类密码子使用偏爱模式.全部密码子使用数据分析表明,人类可变剪接基因密码子的偏爱水平显著高于普通剪接基因.在人类1号染色体基因中,密码子第三位置的G C含量有很大的异质性(0.24～0.95),并且可变剪接基因密码子第三位置平均G C含量(64.66%)大于普通剪接基因(59.97%).Nc值对GC3s图显示密码子偏爱使用除了受核苷酸组成制约外,其它的因子可能也影响密码子的使用变化.此外,可变剪接基因中以G 或C结尾的密码子比普通剪接基因出现的频率高.密码子使用的差异可能是由可变剪接基因pre-mRNA特有的结构特征和多种剪接模式决定的.     

蜜蜂球囊菌几种胞外蛋白酶的特性

 【目的】研究由蜜蜂球囊菌引起的真菌性疾病白垩病的侵染机制,以及蜜蜂球囊菌胞外蛋白酶的种类及其特性。【方法】利用人工合成的培养基,在体外诱导蜜蜂球囊菌分泌胞外蛋白酶,用硫酸铵沉淀和透析,获得初酶液,利用相应的底物对酪蛋白酶、偶氮酪蛋白酶、胶原蛋白酶、弹性蛋白酶和胰蛋白酶的活力进行测定。【结果】蜜蜂球囊菌分泌的胞外酶主要有酪蛋白酶、偶氮酪蛋白酶、胶原蛋白酶3种,无弹性蛋白酶和胰蛋白酶。对这3种酶的最适温度、pH值以及热稳定性和酸碱稳定性进行了测定。这3种胞外酶的最适反应温度都在40℃左右,最适pH为7左右,均不耐高温,热稳定性差,在中性的pH环境下比较稳定。【结论】蜜蜂球囊菌可以分泌多种胞外蛋白酶,它们的最适pH和温度与蜜蜂消化道内的环境一致。     

观赏海棠'火焰'果皮中可变剪接基因分析

[目的]研究观赏海棠'火焰'果皮中可变剪接基因.[方法]利用已有苹果属观赏海棠品种'火焰'果皮转录组测序数据,分析可变剪接类型和数量,并对发生可变剪接的基因进行KEGG分析,研究观赏海棠可变剪接基因与其发育过程的关系.[结果]在观赏海棠果皮发育过程中,主要发生的可变剪接是第一个外显子可变剪接和最后一个外显子可变剪接两种类型.可变剪接在果皮发育前期发生较多,同时KEGG通路分析表明可变剪接主要发生在嘌呤代谢、氧化磷酸化、糖酵解/糖异生等通路过程中.[结论]可变剪接在观赏海棠果皮的整个发育过程中普遍存在.观赏海棠果皮发育过程中共有的可变剪接基因主要与磷酸化和糖代谢有关.     

基于RNA-seq数据大规模挖掘蜜蜂球囊菌的SSR分子标记

基于已获得的球囊菌转录组数据预测微卫星标记(SSR),并进行SSR位点的信息分析和SSR引物的挖掘.利用软件MISA对球囊菌转录组中42610条unigenes进行搜索,共预测出7968个SSRs,分布于5233条unigenes中,其中最主要的重复类型为三核苷酸重复(53.15%),其次为二核苷酸重复(32.32%)和四核苷酸重复(8.46%).二核苷酸重复中的基序主要是AG/CT(占总量的15.8%).针对所有的SSR位点,利用Primer Premier 5软件设计出6956对SSR特异性引物,随机选取20对引物对两个不同来源的球囊菌样品进行SSR位点扩增,共有6对引物成功扩增出符合预期的目的片段.研究结果表明,利用球囊菌转录组数据开发SSR引物是可行的.     

意大利蜜蜂6日龄幼虫肠道内球囊菌的差异表达基因分析

为检测蜜蜂球囊菌(Ascosphaera apis)在胁迫意大利蜜蜂(Apis mellifera ligustica,简称意蜂)幼虫后期的基因表达谱,利用RNA-seq技术对纯培养的球囊菌孢子及球囊菌胁迫的意蜂6日龄幼虫肠道进行深度测序。通过比较处理组和对照组得到21 361个差异表达基因(DEGs),包括15 306个上调基因和6 055个下调基因。GO富集分析结果显示,上调基因富集于39个GO条目(term),基因富集数最多的为细胞进程(2 416unigenes)、细胞(2 408 unigenes)和细胞组件(2 408 unigenes);下调基因富集于38个GO term,基因富集数最多的为代谢进程(1 227 unigenes)、细胞进程(1 185 unigenes)和细胞(1 131 unigenes)。KEGG代谢通路富集分析结果显示,上调基因富集在119个通路,其中基因富集数最多的是核糖体(221 unigenes),其次为内质网中蛋白加工(190 unigenes)和内吞作用(177 unigenes);下调基因富集在113个通路上,基因富集数最多的是核糖体(144 unigenes),其次为RNA转运(113 unigenes)和氨基酸生物合成(108 unigenes)。进一步分析发现,富集在MAPK信号通路上的64个基因上调表达,说明该通路在球囊菌胁迫后期被激活。研究结果不仅为揭示球囊菌在胁迫意蜂幼虫后期的病原-宿主互作提供了重要信息,也为阐明球囊菌致病的分子机制奠定了基础。