摘要：【目的】微小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的表达量进行归一化处理,以|log2 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条目,

关键词：意大利蜜蜂 幼虫肠道 发育 差异表达微小rna 调控网络 

单位：福建农林大学蜂学学院; 福州350002