中华蜜蜂幼虫肠道响应球囊菌胁迫的microRNA应答分析  

MicroRNA responses in the larval gut of Apis cerana cerana to Ascosphaera apis stress

作　　者：杜宇[1] 童新宇 周丁丁 陈大福[1] 熊翠玲[1] 郑燕珍[1] 徐国钧[1] 王海朋 陈华枝 郭意龙 隆琦 郭睿[1] Yu Du;Xinyu Tong;Dingding Zhou;Dafu Chen;Cuiling Xiong;Yanzhen Zheng;Guojun Xu;Haipeng Wang;Huazhi Chen;Yilong Guo;Qi Long;Rui Guo(College of Bee Science,Fujian Agriculture and Forestry University,Fuzhou 350002,Fujian Province,China)

机构地区：[1]福建农林大学蜂学学院

出　　处：《微生物学报》2019年第9期1747-1764,共18页Acta Microbiologica Sinica

基　　金：国家自然科学基金(31702190);现代农业产业技术体系建设专项资金(CARS-44-KXJ7);福建省教育厅中青年教师教育科研项目(JAT170158);福建农林大学科技创新专项基金(CXZX2017342,CXZX2017343);福建省大学生创新创业训练计划(201610389053,201810389082,201810389029)~~

摘　　要：【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种能够侵染中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫的致死性真菌病原。微小RNA(microRNA,miRNA)可通过在转录后水平靶向抑制或降解mRNA而参与宿主与病原互作过程。本研究旨在对球囊菌胁迫的中蜂6日龄幼虫肠道的差异表达miRNA(DEmiRNA)及其靶基因进行深入分析,进而揭示DEmiRNA在中蜂响应球囊菌胁迫应答过程中的作用。【方法】利用Illumina MiSeq平台对正常及球囊菌胁迫的中蜂6日龄幼虫肠道(AcCK和AcT)进行测序,通过相关生物信息学软件预测DEmiRNA及其靶基因。通过Blast将靶基因注释到GO和KEGG数据库。利用Cytoscape软件构建DEmiRNA与其靶mRNA的调控网络。通过Stem-loop RT-PCR和qPCR验证测序数据的可靠性。【结果】本研究共预测出537个miRNA,其长度分布介于16–35 nt之间,且不同长度的miRNA首位碱基偏向性差异明显。通过Stem-loop RT-PCR证实了10个novel miRNA的表达。AcCK vs AcT比较组共有54个DEmiRNA,包含31个上调和23个下调miRNA,可分别靶向结合6170和8199个靶基因。GO分类结果显示上调和下调miRNA的靶基因分别涉及47和47个条目,富集基因数最多的皆为结合细胞进程和催化活性。KEGG代谢通路(pathway)富集分析结果表明上调和下调miRNA的靶基因分别富集在134和126条pathway,富集基因数最多的均为内吞作用和内质网中的蛋白质加工。调控网络分析结果表明,DEmiRNA及其靶mRNA形成十分复杂的调控关系;31个DEmiRNA可靶向结合51个与泛素介导的蛋白水解相关的mRNA,18个DEmiRNA可靶向结合14个与Jak-STAT信号通路相关的mRNA;miR-1277-x、miR-26-x、miR-27-y、miR-30-x、miR-6052-x等16个miRNA共同参与了上述两条免疫通路的调控。最后,随机挑选3个DEmiRNA进行qPCR验证,结果证明了测序数据的可靠性。【结论】本研究提供了中蜂幼虫肠道在球囊菌胁迫后期的miRNA的表达谱和差异表�[Objective]Ascosphaera apis is a lethal fungal pathogen for Apis cerana cerana larvae.microRNA can participate in host-pathogen interaction processes by inhibition or degradation of mRNA via targeting at post-transcriptional level.The aim of this study was to analyze the differentially expressed miRNAs(DEmiRNAs)and their target genes in the 6-day-old larval gut of A.c.cerana under A.apis stress and reveal DEmiRNAs’roles in the stress response process.[Methods]Normal and A.apis-challenged 6-day-old larval guts of A.c.cerana(AcCK and AcT)were sequenced using Illumina MiSeq platform,followed by prediction and analysis of DEmiRNAs and their target genes using related softwares.Target genes of DEmiRNAs were annotated to Gene Ontology(GO)and Kyoto Encyclopedia of Genes and Genomes(KEGG)databases with Blast.Regulation networks between DEmiRNAs and target mRNAs were constructed using Cytoscape.Stem-loop RT-PCR and qPCR were conducted to verify the reliability of sequencing data.[Results]Deep sequencing of larval gut samples generated 537 miRNAs,the length of which was distributed between 16 nt and 35 nt.The first base bias of miRNAs with various length had apparent difference.The expression of 10 novel miRNAs were validated using Stem-loop RT-PCR.There were 54 DEmiRNAs in AcCK vs.AcT comparison group,including 31 up-regulated and 23 down-regulated miRNAs,which can respectively link 6170 and 8199 target genes.GO classification suggested that target genes of up-and down-regulated miRNAs were respectively involved in 47 and 47 terms,and the largest ones were binding,cellular process,and catalytic activity.KEGG pathway enrichment analysis demonstrated that target genes of up-and down-regulated miRNAs were respectively engaged in 134 and 126 pathways,and the mostly enriched ones were endocytosis and protein processing in endoplasmic reticulum.Analysis of regulation networks revealed that very complex regulation relationships existed between DEmiRNAs and corresponding target mRNAs;31 miRNAs could bind 51 mRNAs associated wit

关 键 词：中华蜜蜂 幼虫 球囊菌 MICRORNA 靶基因 调控网络 

分 类 号：R73[医药卫生—肿瘤]