陈大福，郭睿，熊翠玲，梁勤，郑燕珍，徐细建，张曌楠，黄枳腱，张璐，王鸿权，解彦玲，童新宇



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

陈大福，郭睿，熊翠玲，梁勤，郑燕珍，徐细建，张曌楠，黄枳腱，张璐，王鸿权，解彦玲，童新宇

（福建农林大学蜂学学院，福州 350002）

【目的】白垩病是困扰养蜂生产的顽疾。目前，尚无利用二代测序技术研究中华蜜蜂（,简称中蜂）幼虫白垩病的报道。本研究利用RNA-seq技术对健康（AcCK）及球囊菌()胁迫的中蜂4日龄幼虫肠道（AcT）进行深度测序，在转录组水平研究中蜂幼虫在球囊菌胁迫早期的胁迫应答。【方法】通过Illumina HiSeq 2500平台对AcCK和AcT进行双端（PE125）测序，首先对测序数据进行质控和评估，利用edgeR软件进行差异表达基因（DEG）分析，进而对DEGs进行GO富集分析及KEGG代谢通路富集分析，最后，利用实时荧光定量PCR（qRT-PCR）验证测序数据的可靠性。【结果】 AcCK和AcT的转录组测序共得到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），进一步分析表明中蜂幼虫肠道的部分细胞免疫对球囊菌的胁迫产生应答，而体液免疫不产生应答，宿主的代谢相关基因受到球囊菌的显著抑制。……