球囊菌胁迫中华蜜蜂幼虫肠道过程中病原的转录组学研究

【目的】本研究利用RNA-seq技术对球囊菌胁迫的中华蜜蜂(中蜂)幼虫肠道进行深度测序,经趋势分析得到差异表达基因(DEGs)的显著表达模式,进而对胁迫过程中的球囊菌进行转录组学分析。【方法】利用Illumina HiSeq 2500平台对球囊菌胁迫的中蜂幼虫肠道进行深度测序,并利用相关软件进行了深入分析。最后,通过RT-qPCR对RNA-seq数据进行了验证。【结果】本研究共得到球囊菌的41133932条高质量clean reads。22865个DEGs共聚类为8个基因表达模式,其中,16769个DEGs聚类为2个显著上调趋势与2个显著下调趋势。GO富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于40与37个GO term,基因富集数最多的为细胞进程(2486 unigenes)。KEGG代谢通路(pathway)富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于119和112个pathway,基因富集数最多的分别是氨基酸生物合成(127 unigenes)与核糖体(98 unigenes)。进一步分析表明球囊菌在胁迫中蜂幼虫肠道的过程中通过提高物质合成促进其增殖,而宿主通过抑制球囊菌的蛋白合成抵御病原入侵。富集在MAPK信号通路的11个DEGs的表达水平随着胁迫时间的延长而逐渐下降,推测中蜂幼虫通过抑制该通路而阻遏球囊菌增殖。【结论】本研究不仅为揭示白垩病过程中的球囊菌-中蜂幼虫互作提供了重要信息,也为阐明不同抗性蜂种的球囊菌抗性差异奠定了基础。     

意大利蜜蜂工蜂中肠响应Nosema ceranae胁迫的高表达基因分析

东方蜜蜂微孢子虫Nosema ceranae是专性寄生蜜蜂中肠上皮细胞的单细胞真菌,对意大利蜜蜂(意蜂)具有较强的侵染性。本研究利用RNA-seq技术对正常10 d (Apis mellifera ligustica control group,Am CK)和N. ceranae胁迫10 d的意蜂工蜂中肠(Apis mellifera ligustica treatment group,Am T)进行测序,共得到160 847 237 100条原始读段,过滤后得到1 066 955 298条有效读段。主成分分析结果显示Am CK与Am T测序样品的组内重复性较好,组间的基因表达模式差异明显。GO分类结果显示,Am CK与Am T的前100位高表达基因(HEGs)分别分布于32和33个GO terms,基因富集数最多的均为代谢进程。KEGG代谢通路(pathway)富集分析结果显示,Am CK的前100位HEGs富集在21个pathways,基因富集数最多的是内吞作用、信号通路和嘌呤代谢; Am T的前100位HEGs富集在26个pathways,基因富集数最多的是内吞作用、RNA转运和蛋白质的内质网加工。前100位HEGs的Venn分析结果显示Am CK与Am T的共有HEGs为87个,特有HEGs分别均为13个。研究结果揭示了N. ceranae胁迫意蜂工蜂中肠过程中的基因表达谱信息,也为解析N. ceranae致病的分子机理提供了有益信息和线索。