侵染中华蜜蜂6日龄幼虫的蜜蜂球囊菌的微小RNA差异表达谱及调控网络

【目的】蜜蜂球囊菌(Ascosphaeraapis,简称球囊菌)是专性侵染蜜蜂幼虫的致死性真菌病原。MicroRNA(miRNA)作为一类重要的基因表达调控因子,能够广泛参与真菌及其宿主的相互作用过程。本研究通过比较分析球囊菌孢子(AaCK)和侵染中华蜜蜂(Apis cerana cerana,简称中蜂) 6日龄幼虫肠道内的球囊菌(AaT)的smallRNA(sRNA)组学数据对球囊菌的差异表达miRNA(differentiallyexpressed miRNA,DEmiRNA)、靶mRNA及二者间的调控网络进行全面解析,旨在揭示miRNA介导的球囊菌对中蜂幼虫的侵染机制。【方法】对于球囊菌侵染的中蜂6日龄幼虫肠道的small RNA-seq (sRNA-seq)数据,利用BLAST工具连续比对东方蜜蜂(Apiscerana)和球囊菌的参考基因组筛滤得到AaT的sRNA组学数据。分别将AaCK和AaT的sRNA组学数据比对miRBase数据库,对球囊菌侵染宿主前后miRNA的数量和结构特征进行分析。联用RNAhybrid+svm_light、Miranda和Tar...     

参与调控意大利蜜蜂工蜂中肠基因表达的东方蜜蜂微孢子虫miRNA的组学解析及其调控网络

[目的]本研究结合前期已获得的miRNA和mRNA组学数据对东方蜜蜂微孢子虫Nosema ceranae 的差异表达 miRNA(differentially expressed miRNA,DEmiRNA)靶向意大利蜜蜂 Apis mellifera ligustica 工蜂中肠的 mRNA 和差异表达 mRNA(d...     

中华蜜蜂6日龄幼虫响应蜜蜂球囊菌侵染的长链非编码RNA应答研究

[目的]本研究旨在探究长链非编码RNA (long non-coding RNA,lncRNA)在中华蜜蜂(Apis cerana cerana,简称中蜂)6日龄幼虫应答蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)侵染过程中的差异表达谱及调控作用.[方法]利用链特异性cDNA建库的RNA-seq技术对未被...     

蜜蜂球囊菌菌丝和孢子中毒力因子相关全长转录本的鉴定及分析

[目的]蜜蜂球囊菌Ascosphaera apis是一种专性侵染蜜蜂幼虫而导致白垩病的致死性真菌病原.基于前期已获得的高质量纳米孔(Nanopore)长读段测序数据,对蜜蜂球囊菌菌丝(AaM)和孢子(AaS)中的毒力因子相关全长转录本进行鉴定和分析,为毒力因子相关剪接异构体的功能研究提供参考信息和基础.[方法]利用BLAST工具,将AaM和AaS中的所有全长转录本比对Nr数据库,以鉴定蜜蜂球囊菌的毒力因子(几丁质酶、脂肪酶、水解酶和蛋白酶)相关的全长转录本.使用minimap2软件,将AaM和AaS中的全长转录本序列与蜜蜂球囊菌参考基因组注释的已知转录本序列进行比对,将比对到参考基因组的全长转录本进行归一化处理,再通过每百万计数法(Counts per million,CPM)计算毒力因子相关全长转录本的表达量.利用百迈克云平台的相关工具绘制转录本的表达量聚类热图.通过IGV浏览器对部分毒力因子相关全长转录本结构进行可视化.[结果]在AaM鉴定到毒力因子相关的367个基因及407个全长转录本,包括12条几丁质酶相关全长转录本,48条脂肪酶相关全长转录本,289条水解酶相关全长转录本,58条蛋白酶相关全长转录本.在AaS鉴定到毒力因子相关367个基因及400个全长转录本,包括14条几丁质酶相关全长转录本,63条脂肪酶相关全长转录本,267条水解酶相关全长转录本,56条蛋白酶相关全长转录本.另外,AaM和AaS特有的毒力因子(几丁质酶、脂肪酶)相关全长转录本分别有0条和17条,共有的毒力因子相关全长转录本有60条.进一步分析发现,蜜蜂球囊菌的部分毒力因子基因可通过可变剪接形成多条剪接异构体.[结论]共鉴定到蜜蜂球囊菌毒力因子(几丁质酶、脂肪酶、水解酶和蛋白酶)相关的367个基因和486条全长转录本;相比于蜜蜂球囊菌参考基因组注释的转录本,绝大多数毒力因子基因对应的全长转录本数量更多且结构更为复杂.研究结果丰富了蜜蜂球囊菌毒力因子相关基因和转录本的注释信息,为毒力因子相关剪接异构体的功能研究提供了基础,也为白垩病防控提供了潜在靶点.     

基于蜜蜂球囊菌纳米孔测序数据的基因非翻译区延长、SSR位点发掘及未注释基因和转录本鉴定

【目的】利用已获得的纳米孔长读段测序数据完善现有的蜜蜂球囊菌Ascosphaera apis参考基因组注释信息,并对未注释的新基因和新转录本进行鉴定和功能注释。【方法】基于已获得的纳米孔长读段测序数据,采用gffcompare软件将蜜蜂球囊菌全长转录本与参考基因组注释的转录本进行比较,进而对参考基因组注释基因的非翻译区(untranslated region, UTR)进行延长。利用TransDecoder软件对蜜蜂球囊菌基因的开放阅读框(open reading frame, ORF)及相应的氨基酸序列进行预测。通过MISA软件发掘长度在500 bp以上的全长转录本的SSR位点。通过Blast工具将鉴定到的新基因和新转录本比对Nr, KOG, eggNOG, Swiss-Prot, Pfam, GO和KEGG数据库进行功能注释。【结果】共对蜜蜂球囊菌的9 481个基因进行了UTR延长,其中5′UTR和3′UTR延长的基因分别有4 744和4 737个。共预测出10 492个完整ORF,其中编码长度分布在0～100和100～200个氨基酸的ORF最多,分别占ORF总数的38.96%和3...     

蜜蜂球囊菌菌丝和孢子中全长转录本的差异表达分析

[目的]本研究旨在探究蜜蜂球囊菌Ascosphaera apis(简称“球囊菌”)菌丝(Aam)和孢子(Aas)中的全长转录本与菌丝生长、孢子萌发以及有性生殖的相关性.[方法]将前期获得的Aam和Aas的Nanopore长读段测序数据中的有效读段与球囊菌参考基因组注释的已知转录本进行序列比对,获取全长转录本与已知转录本...     

过表达和敲减ame-miR-13b对意大利蜜蜂幼虫肠道内基因表达的影响

【目的】微小RNA(microRNAs, miRNAs)在昆虫的繁殖、发育和免疫等重要生命活动的调控过程中扮演关键角色。本研究旨在探究在意大利蜜蜂Apis mellifera ligustica幼虫个体水平过表达与敲减ame-miR-13b对其靶基因表达的影响,为深入探究ame-miR-13b调控意大利蜜蜂幼虫肠道发育的分子机理提供理论和实验依据。【方法】根据ame-miR-13b的核苷酸序列设计合成相应的模拟物mimic-ame-miR-13b和抑制物inhibitor-ame-miR-13b及其对照mimic-NC和inhibitor-NC,混入饲料后饲喂意大利蜜蜂3日龄幼虫,每12 h更换一次饲料,连续饲喂6次,以分别过表达和敲减ame-miR-13b;通过RT-qPCR检测过表达和敲减ame-miR-13b后意大利蜜蜂幼虫肠道中ame-miR-13b及其靶基因Ecr, Egfr和P450 18a1的表达量。【结果】与mimic-NC饲喂组相比, mimic-ame-miR-13b饲喂组中ame-miR-13b在意大利蜜蜂4和6日龄幼虫肠道中均显著上调表达;与饲喂inhibitor-NC组相比,inhibitor-ame-miR-13b饲喂组中ame-miR-13b在4日龄幼虫肠道中下调表达,在5和6日龄幼虫肠道中显著下调表达。与mimic-NC饲喂组相比,过表达ame-miR-13b组中意大利蜜蜂6日龄幼虫肠道中Egfr的表达量显著降低,而Ecr和P450 18a1的表达量均显著升高。与inhibitor-NC饲喂组相比,敲减ame-miR-13b组中意大利蜜蜂6日龄幼虫肠道中Egfr上调表达但不显著,Ecr显著上调表达,而P450 18a1显著下调表达。【结论】通过饲喂法在意大利蜜蜂幼虫个体中成功实现了miRNA的过表达和敲减;ame-miR-13b与Egfr之间存在潜在的负调控关系。