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

摘要

利用RNA-seq技术对胁迫中华蜜蜂(简称中蜂)6日龄幼虫肠道的球囊菌及其体外培养进行深度测序,根据基因的FPKM值筛选得到高表达基因(HEGs),进而对HEGs进行GO及KEGG代谢通路(pathway)富集分析.Illumina测序数据经质控和过滤得到100814558条有效读段(clean reads),平均Q30均在93.81%以上.GO富集分析结果显示处理组(AacT)的HEGs富集于39个GO term,基因富集数最多的是细胞(1872 unigene)、细胞组件(1872 unigene)和细胞进程(1748 uni-genes);对照组(AaCK)的HEGs富集于37个GO term,基因富集数最多的是细胞(785 unigenes),其次是细胞组件(785 uni-genes)和代谢进程(776 unigenes).KEGG pathway富集分析显示,AacT的HEGs富集在119个代谢通路上,基因富集数最多的是核糖体(176 unigenes)、碳代谢(147 unigenes)及氨基酸的生物合成(134 unigenes);AaCK的HEGs富集在110个代谢通路上,基因富集数最多的是核糖体(179 unigenes)、氨基酸的生物合成(70 unigenes)以及碳代谢(62 unigenes).深入分析发现,对于富集在MAPK信号通路上的高表达基因数量,胁迫中蜂幼虫肠道的球囊菌远多于体外培养的球囊菌,说明病原的该通路在胁迫后期被显著激活.%In this study, Ascosphaera apis spores (AacT), from the gut of 6-day-old Apis cerana cerana larvae, and A.c.cerana were sequenced utilizing RNA-seq technology. Pure A. apis spore was used as the control ( AaCK) . A total of 100814558 clean reads with a mean Q30 of 93.81% were obtained after quality control and filtration. GO enrichment analysis suggested that the highly ex-pressed genes ( HEGs) in AacT were involved in 39 terms, which was mostly enriched in cell ( 1872 unigenes) , followed by cellular component ( 1872 unigenes) and cellular process ( 1748 unigenes) . And the HEGs in AaCK were engaged in 37 terms, with the lar-gest group in cell ( 785 unigenes) , followed by cellular component ( 785 unigenes) and metabolic process ( 776 unigenes) . KEGG enrichment analysis showed that the HEGs in AacT were enriched in 119 pathways, dominating in ribosme (176 unigenes), fol-lowed by carbon metabolism (147 unigenes) and biosynthesis of amino acids (134 unigenes). While the HEGs in AaCK were in-volved in 119 pathways, and the largest groups were ribosome ( 179 unigenes) , biosynthesis of amino acids ( 70 unigenes) and car-bon metabolism ( 62 unigenes) .