胁迫意大利蜜蜂幼虫肠道的球囊菌的转录组分析  

Transcriptomic analysis of Ascosphaera apis stressing larval gut of Apis mellifera ligustica(Hyemenoptera:Apidae)

作　　者：陈大福[1] 郭睿[1] 熊翠玲[1] 梁勤[1] 郑燕珍[1] 徐细建[1] 黄枳腱 张曌楠 张璐[1] 李汶东 童新宇 席伟军 

机构地区：[1]福建农林大学蜂学学院,福州350002 [2]金华农业科学研究院,浙江金华321000

出　　处：《昆虫学报》2017年第4期401-411,共11页Acta Entomologica Sinica

基　　金：国家自然科学基金项目(30800806);现代农业产业技术体系建设专项资金(CARS-45-KXJ7);福建农林大学科技发展资金(KF2015123)

摘　　要：【目的】本研究旨在通过趋势分析对胁迫意大利蜜蜂Apis mellifera ligustica(简称意蜂)幼虫肠道的球囊菌Ascosphaera apis的差异表达基因(DEGs)进行转录组分析。【方法】将纯化的球囊菌孢子配制为1×107孢子/m L的饲料饲喂意蜂3日龄幼虫,利用Illumina Hi Seq 2500平台对球囊菌胁迫的意蜂幼虫肠道c DNA进行测序,将过滤得到的有效读段(clean reads)映射(mapping)到核糖体数据库及意蜂参考基因组,最后将未映射上的reads映射到本课题组组装并注释的球囊菌参考转录组。利用STEM软件对DEGs进行趋势分析。利用WEGO软件对显著表达模式中的DEGs进行GO富集分析。利用Blastall对显著表达模式中的DEGs进行KEGG代谢通路富集分析。最后,通过对随机选取的6个DEGs进行RT-q PCR分析,对RNA-seq数据进行验证。【结果】球囊菌胁迫意蜂幼虫肠道样品的Illumina测序共得到球囊菌的25 454 076条原始读段(raw reads),经过滤得到24 909 820条clean reads,Q30均在93.46%以上。趋势分析结果显示,19 893个DEGs聚类为8个表达模式,其中有12 151个DEGs聚类为3个表现为显著上调趋势的表达模式。GO富集分析结果显示,表现上调趋势的DEGs富集在40个GO term,富集基因数最多的为细胞进程(cellular process)(2 601 unigenes),其次为代谢进程(metabolic process)(2 553 unigenes)和细胞(cell)(2 522unigenes)。KEGG代谢通路富集分析结果显示,上调趋势中的DEGs富集在119个代谢通路上,其中富集基因数最多的是核糖体(ribosome)(213 unigenes),其次为氨基酸生物合成(biosynthesis of amino acids)(154 unigenes)和内质网蛋白加工(protein processing in endoplasmic reticulum)(130unigenes)。共有48个DEGs富集在MAPK信号通路上,聚类分析结果显示,这些DEGs随着胁迫时间的延长表达水平逐渐增强。RT-q PCR结果显示,6个DEGs的表达水平变化趋势与RNA-seq数据一致,证明了本研究中的转录组数据真实可靠。【结论】本研�[ Aim ] This research is designed to conduct transcriptomic analysis of the differentially expressed genes （DEGs） of Ascosphaera apis stressing larval guts of Apis mellifera ligustica via trend analysis. [ Methods] The purified A. apis spores at a concentration of 1 x 107 spores/mL was used to feed 3-day-old larvae of A. m. ligustica, and then the eDNA of stressed larval guts was sequenced at Illumina HiSeq 2500 platform. After filtration, the clean reads were used to mapping the ribosome database and the reference genome of A. m. ligustica, and the unmapped reads were used to mapping the reference transcriptome of A. apis assembled previously. The STEM software was used to analyze the gene expression patterns. Gene ontology （GO） enrichment analysis for DEGs involved in significant expression patterns was performed using WEGO software. KEGG enrichment analysis for DEGs associated with significant expression patterns was carried out by using Blastall. Finally, RT-qPCR analysis of six randomly selected DEGs was performed to validate the RNA-seq data. [ Results ] The RNA-seq of A. apis produced 25 454 076 raw reads, and after filtration, 24 909 820 clean reads with Q30 above 93.46% were obtained. Trend analysis results showed that 19 893 DEGs were grouped into eight gene expression patterns, among which 12 151 DEGs were assigned to three expression patterns with significantly up- regulated expression trend. GO enrichment analysis results indicated that all DEGs within significantly up-regulated expression patterns were enriched in 40 GO terms, and the mostly enriched one was cellular process （2 601 unigenes） , followed by metabolic process （2 553 unigenes） and cell （2 522 unigenes）. KEGG enrichment analysis result displayed that all DEGs within significantly up-regulated expression patterns were enriched in 119 metabolism pathways, and the mostly enriched one was ribosome （213 unigenes）, followed by biosynthesis of amino acids （154 unigenes） and protein processing in endoplasmic reticulum

关 键 词：意大利蜜蜂 球囊菌 幼虫肠道 RNA-SEQ 转录组 差异表达基因 

分 类 号：Q966[生物学—昆虫学]