中华蜜蜂6日龄幼虫响应蜜蜂球囊菌胁迫的环状RNA应答

【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)专性侵染蜜蜂幼虫而导致白垩病,危害蜜蜂健康和养蜂生产。本研究旨在探究中华蜜蜂(Apis cerana cerana,简称中蜂) 6日龄幼虫响应球囊菌胁迫的环状RNA(circular RNA,circRNA)差异表达谱及差异表达circRNA (differentially expressed circRNA,DEcircRNA)在宿主胁迫应答中的潜在功能。【方法】利用去除线性RNA的circRNA-seq技术对正常和球囊菌侵染的中蜂6日龄幼虫肠道(AcCK和AcT)进行测序。利用find_circ软件鉴定circRNA,统计circRNA的长度和环化类型。根据|log_2(Fold change)|≥1和P≤0.05的标准筛选DEcircRNA。将DEcircRNA的来源基因比对Gene ontology (GO)数据库和Kyoto Encyclopedia of Genes and Genomes (KEGG)数据库,从而获得功能及通路(pathway)注释。随机挑选3个DEcircRNA进行RT-qPCR验证。【结果】AcCK和AcT的circRNA-seq分别得到76342570和68269362条原始读段(raw reads),经严格质控得到74524108和66974392条有效读段(clean reads),Q30分别为92.75%和94%,GC含量分别为54.31%和54.90%。比对上东方蜜蜂(Apis cerana)参考基因组的短序列读段(anchor reads)共计23648400条。AcCK和AcT中分别鉴定到805和702个circRNA,长度均介于201–1000 nt,数量最多的环化类型均为已注释外显子circRNA,但分布在不同长度、不同环化类型的circRNA数量存在差异。AcCK vs AcT比较组共有494个DEcircRNA,包括257个上调circRNA和237个下调circRNA;上调和下调幅度最大的circRNA分别为novel_circ_000123和novel_circ_000726。上述DEcircRNA的来源基因可注释到11条生物学进程相关条目,9条分子功能相关条目,9条细胞组分相关条目,以及73条通路。进一步分析发现,部分DEcircRNA的来源基因注释到7条细胞免疫通路和3条体液免疫通路。【结论】中蜂6日龄幼虫响应球囊菌胁迫的过程中可能通过改变分布在不同长度和环化类型的circRNA数量,以及特异性表达一些circRNA和调节部分circRNA的表达量对病原产生应答;novel_circ_000027、novel_circ_000127、novel_circ_000312等DEcircRNA在宿主的胁迫应答过程中可能通过调控氧化磷酸化、细胞和体液免疫等通路发挥特殊作用。研究结果为深入理解中蜂幼虫对球囊菌的胁迫应答机制及二者的相互作用机制提供了新见解。     

LncRNA13164通过ace-miR-4968-y调节中华蜜蜂幼虫对蜜蜂球囊菌侵染的免疫应答

【目的】长链非编码RNA (long non-coding RNA,lncRNA)可通过竞争性内源RNA(competing endogenous RNA,ce RNA)等多种方式发挥重要的调控作用,但蜜蜂lncRNA的功能研究至今依然缺失,lncRNA在蜜蜂免疫应答中的作用尚不清楚。本研究旨在揭示中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫响应蜜蜂球囊菌(Ascosphaera apis)侵染的免疫应答中lncRNA的调控功能及作用机制。笔者团队前期通过深度测序和生物信息学分析发现,lncRNA13164与ace-miR-4968存在靶向关系且潜在参与在中蜂幼虫对蜜蜂球囊菌侵染的应答。【方法】利用实时荧光定量PCR (real-time fluorescence quantitative polymerase chain reaction,RT-qPCR)检测蜜蜂球囊菌接种后中蜂幼虫肠道内lncRNA13164的表达谱。采用ILoc-Lnc RNA软件预测lncRNA13164的亚细胞定位。联用RNAhybrid、Miranda和Target Scan软件预测lnc...     

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

【目的】蜜蜂球囊菌(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的表达谱和差异表达信息,揭示了球囊菌与宿主之间在miRNA组学水平存在复杂的互作。miR-6052-x和miR-1277-x作为调控网络的核心可能通过影响细胞凋亡参与宿主的免疫防御,miR-26-x和miR-30-x可能通过调控Jak-STAT信号通路参与宿主的胁迫应答。本研究筛选出的关键DEmiRNA有望作为治疗白垩病的分子靶标。     

蜜蜂球囊菌菌丝和孢子中环状RNA的鉴定及比较分析

[目的] 本研究旨在明确蜜蜂球囊菌（Ascosphaera apis,简称球囊菌）菌丝和孢子中环状RNA（circular RNA,circRNA）的数量、种类和表达谱差异,并探讨共有circRNA、特有circRNA和差异表达circRNA（differentially expressed circRNA,DEcircRNA）在菌丝与孢子中的潜在作用。[方法] 基于前期获得的球囊菌菌丝（AaM）和孢子（AaS）的高质量RNA-seq数据,利用find_circ软件预测circRNA。通过Venn分析筛选AaM和AaS的共有circRNA和特有circRNA。根据P ≤ 0.05且|log₂ fold change|≥ 1的标准筛选AaMvs.AaS比较组的DEcircRNA。通过比对GO和KEGG数据库对circRNA的来源基因进行功能和通路注释。利用TargetFinder软件预测circRNA靶向结合的miRNA及miRNA靶向结合的mRNA。采用Cytoscape软件对竞争性内源RNA（competing endogenous RNA,ceRNA）调控网络进行可视化。通过RT-qPCR对DEcircRNA进行验证。[结果] AaM和AaS中分别含有13210156和19011000条短序列读段（anchors reads）,其中分别有6124922和11392886条能够比对上球囊菌参考基因组。在AaM和AaS中分别鉴定到1868个和2225个circRNA,二者共有的circRNA为1098个,AaM的特有circRNA为770个,AaS的特有circRNA为1127个。AaM和AaS的circRNA长度主要介于1000-2000 nt,基因间区circRNA为主要环化类型。AaMvs.AaS比较组包含456个上调circRNA和97个下调circRNA。共有circRNA的来源基因可注释到29个功能条目和14类通路;AaM的特有circRNA的来源基因可注释到31个功能和17类通路;AaS的特有circRNA的来源基因可注释到34个功能条目和16类通路;DEcircRNA的来源基因可注释到29个功能条目和40条通路。调控网络分析结果显示,36个共有circRNA靶向4个miRNA进而调控6个与内吞作用通路相关的靶mRNA;4（255）个AaM（AaS）的特有circRNA靶向2（2）个miRNA进而调控8（2）个次级代谢产物生物合成通路相关的靶mRNA;9个DEcircRNA靶向2个DEmiRNA进而调控3个MAPK信号通路相关的DEmRNA。RT-qPCR结果显示10个DEcircRNA的表达趋势与测序数据一致,证实了测序数据的可靠性。[结论] 球囊菌菌丝和孢子的共有circRNA、特有circRNA和DEcircRNA可能通过调控来源基因表达和充当ceRNA的方式调节球囊菌的物质和能量代谢、内吞作用、次级代谢产物生物合成和MAPK信号通路,进而影响球囊菌菌丝生长、孢子萌发和致病性。     

环状RNA ame_circ_000115调节蜜蜂球囊菌胁迫下意大利蜜蜂幼虫肠道基因表达

环状RNA (circular RNA, circRNA)是一类新型非编码RNA,已被证实可通过竞争性内源RNA(competing endogenous RNA,ceRNA)调控网络等多种方式参与调节昆虫基因表达和免疫应答。目前,circRNA在蜜蜂免疫应答中的作用尚不清楚。本研究对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)ame＿circ＿000115(简称ac115)的反向剪切(back splicing,BS)位点进行验证。采用RT-qPCR检测蜜蜂球囊菌(Ascosphaera apis)胁迫下意蜂幼虫肠道内ac115的表达谱,利用双荧光素酶报告基因试验验证ac115与ame-miR-13b的结合关系。通过饲喂特异性siRNA对蜜蜂球囊菌胁迫的意蜂幼虫肠道内ac115进行干扰,进而检测干扰ac115对宿主免疫应答相关的6个基因表达量的影响。结果表明,ac115的BS位点真实存在;相较于未接种组,蜜蜂球囊菌接种组4日龄幼虫肠道内ac115的表达量极显著上调(P<0.000 1),5和6日龄幼虫肠道内ac115显著上调表达(P<0.05)...     

意大利蜜蜂工蜂幼虫nkd基因的生物信息学分析及功能研究

【目的】意大利蜜蜂(Apis mellifera ligustica,简称意蜂)是西方蜜蜂(Apis mellifera)的亚种之一。蜜蜂球囊菌(Ascosphaera apis)侵染意蜂幼虫导致白垩病。本研究对西方蜜蜂裸表皮蛋白(naked cuticle, Nkd)进行保守基序预测和系统进化分析,并通过RNAi明确nkd基因对意蜂工蜂幼虫体重及宿主响应蜜蜂球囊菌胁迫的免疫应答的影响,以期丰富西方蜜蜂基因nkd的信息,并揭示意蜂幼虫nkd的功能。【方法】通过MEME软件预测西方蜜蜂和其他9个物种Nkd蛋白的保守基序。采用MEGA X软件对西方蜜蜂及其他9个物种的Nkd蛋白进行系统进化分析。通过饲喂dsRNA对意蜂幼虫肠道内的nkd进行RNAi。使用电子天平对幼虫进行称重。利用RT-qPCR检测nkd基因的干扰效率及免疫基因的相对表达量。【结果】西方蜜蜂与东方蜜蜂、柑橘凤蝶、家蚕和金凤蝶的Nkd蛋白均含有3个保守基序(motif 1、motif 2 和motif 3),说明上述5个昆虫物种的Nkd具有较高的保守性。西方蜜蜂与东方蜜蜂的Nkd蛋白聚为一支,说明二者的亲缘关系近。与dsRNA-egfp组相比,dsRNA-nkd组5日龄和6日龄幼虫肠道内nkd的表达量均极显著下调(P<0.001),干扰效率分别为49.60%和56.40%。另外,dsRNA-nkd组幼虫体重较dsRNA-egfp组显著下降,说明nkd显著影响幼虫体重。RT-qPCR结果显示,4日龄幼虫肠道内abaecin、apidaecin、birc5、defensin-1和PGRP-S2均被激活表达;5日龄幼虫肠道内abaecin、apidaecin、birc5和defensin-1均被激活表达,PGRP-S2的表达受到抑制;6日龄幼虫肠道内abaecin被激活表达,而apidaecin、birc5、defensin-1和PGRP-S2的表达均受到抑制,说明上述5个免疫基因在宿主响应胁迫的过程中呈不同的表达趋势,均参与宿主的免疫应答,nkd与abaecin和apidaecin的表达存在负向调控关系。【结论】西方蜜蜂的Nkd蛋白含有3个保守基序(motif 1、motif 2和motif 3),西方蜜蜂与东方蜜蜂的Nkd蛋白亲缘关系最近,通过饲喂dsRNA能有效干扰意蜂工蜂幼虫肠道内nkd表达,nkd影响意蜂工蜂幼虫体重及宿主对蜜蜂球囊菌胁迫的免疫应答。     

侵染中华蜜蜂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和TargetScan软件预测AaCK vs AaT比较组中DEmiRNA的靶mRNA,进而利用相关生物信息学软件对上述靶mRNA进行GO分类和KEGG代谢通路富集分析。通过Cytoscape软件对DEmiRNA-mRNA调控网络进行可视化。利用Stem-loop RT-PCR、RT-qPCR和分子克隆验证测序结果的可靠性。【结果】在AaCK和AaT中分别鉴定到380和387个miRNA。结构特征分析结果显示,AaCK和AaT的mi RNA皆集中分布在18–25 nt,且首位碱基主要偏向于U。AaCKvsAaT比较组共有270个DEmiRNA,包含155个上调miRNA和115个下调miRNA,分别靶向结合6091和6145个mRNA。GO分类结果显示,上述靶mRNA主要涉及代谢进程、细胞进程、应激反应等15个生物学进程;细胞、细胞组分、细胞器等12个细胞组分;催化活性、结合、转运子活性等11个分子功能。KEGG代谢通路富集分析结果显示,上述靶mRNA富集在123条代谢通路,参与对氨基酸代谢、碳水化合物代谢以及核苷酸代谢等物质代谢,氧化磷酸化、硫代谢、氮代谢等能量代谢,以及MAPK和Hippo等信号通路的调控。球囊菌DEmiRNA与靶mRNA之间存在复杂的调控关系,其中miR-29-x、miR-250-x、miR-4968-y、miR-11200-x、novel-m0023-5p、novel-m0130-5p和novel-m0135-5p等DEmiRNA可靶向结合与球囊菌的半胱氨酸蛋白酶、DNA甲基化转移酶以及几丁质酶相关的mRNA;此外,miR-7-x、miR-9-z、miR-319-y和miR-5951-y等同时参与调控MAPK信号通路;进一步分析发现,miR-250-x同时参与对DNA甲基化转移酶、MAPK信号通路及其他酶类合成与代谢途径的调控,并可能参与球囊菌与中蜂6日龄幼虫之间的跨界调控。通过Stem-loopRT-PCR和RT-qPCR验证了4个DEmiRNA的差异表达,并利用分子克隆和Sanger测序证实miR-7-x的序列与测序结果一致。【结论】本研究解析了侵染中蜂6日龄幼虫的球囊菌的miRNA差异表达谱及DEmiRNA的调控网络,揭示了球囊菌DEmiRNA可能通过调控病原的物质和能量代谢、增殖、毒力、信号通路及相关mRNA参与对中蜂幼虫的侵染过程。miR-7-x、miR-250-x、novel-m0023-5p等关键DEmiRNA有望作为白垩病治疗的新型分子靶点。     

意大利蜜蜂幼虫肠道在球囊菌胁迫后期的差异表达微小RNA及其靶基因分析

【目的】蜜蜂球囊菌Ascosphaera apis(简称球囊菌)是一种特异性侵染蜜蜂幼虫肠道的致死性真菌病原。微小RNA(microRNA,mi RNA)是一种在转录后水平对m RNA进行负调控的关键调控因子。本研究旨在对意大利蜜蜂Apis mellifera ligustica(简称意蜂)幼虫肠道在球囊菌胁迫后期的差异表达mi RNAs(DEmiRNAs)及其靶基因进行深入分析,为揭示DEmiRNAs在胁迫应答中的作用提供重要信息。【方法】利用small RNA-seq(s RNA-seq)技术对正常及球囊菌侵染的意蜂6日龄幼虫肠道(分别表示为Am CK和Am T)进行测序,通过相关生物信息学软件对DEmiRNAs进行预测和分析。利用TargetFinder软件预测DEmiRNAs的靶基因,然后利用Blast软件对靶基因进行GO和KEGG数据库的功能注释。通过Cytoscape软件构建DEmiRNAs与其靶m RNAs的调控网络。通过茎环实时荧光定量PCR(stem-loop RT-qPCR)验证测序数据的可靠性。【结果】Am CK和Am T的测序分别得到9 230 496和10 823 667条有效序列标签; Am CK vs Am T比较组中包含15个上调和6个下调mi RNAs,分别结合3 503和3 252个靶基因,它们可分别注释到40和39个GO terms以及104和99条代谢通路;进一步分析发现调控网络中的17个DEmiRNAs靶向结合116个与丝氨酸蛋白酶相关的m RNAs,14个DEmiRNAs靶向结合54个与泛素介导的蛋白水解相关的m RNAs。Stem-loop RT-qPCR验证结果显示随机选择的4个DEmiRNAs (mi R-251-x,mi R-9277-y,mi R-1672-x和mi R-4968-y)的表达量变化趋势与测序结果一致,表明测序数据真实可信。【结论】本研究率先对意蜂幼虫肠道在球囊菌胁迫后期的DEmiRNAs及其靶基因进行预测与分析,提供了mi RNAs的表达谱和差异表达信息。ame-miR-927b,mi R-429-y和mi R-8440-y等DEmiRNAs可能参与对丝氨酸蛋白酶的调控和对泛素介导的蛋白水解的调控,DEmiRNAs与靶基因之间存在复杂的调控网络关系。DEmiRNAs可参与到意蜂幼虫肠道与球囊菌间的互作。     

影响蜜蜂球囊菌侵染蜜蜂幼虫的因素及侵染过程观察

为探究蜜蜂球囊菌Ascosphaera apis只侵染封盖前后蜜蜂幼虫的原因及相关侵染机制, 本研究利用实验室饲养的意大利蜜蜂Apis mellifera ligustica幼虫, 给其接种球囊菌孢子, 探究不同接种量（0, 1.0×10², 1.0×10³, 1.0×10⁴, 1.0×10⁵ 和1.0×10⁶ 孢子/mL）、 接种时期（3, 4, 5和6龄幼虫）以及28℃低温处理6 h对蜜蜂球囊菌侵染的影响。同时对处于不同侵染阶段的蜜蜂幼虫做病理学切片, 探究球囊菌侵染的过程。结果显示： 球囊菌孢子接种量与蜜蜂的发病率密切相关（r=0.9883）, 蜜蜂幼虫对低于1.0×10³孢子/mL的侵染有抗性, 与不接孢子的对照组差异不显著（P>0.05）。蜜蜂不同龄期接种发病率的差异是因不同龄幼虫食量不同导致的摄入孢子剂量的不同引起的, 28℃低温处理能够显著提高处于幼虫到蛹转化期蜜蜂的发病率（P<0.05）, 而对取食阶段的蜜蜂幼虫没有影响。病理学研究表明, 在整个幼虫期, 摄入的孢子因中肠没有氧气不生长, 对幼虫没有致病性, 幼虫的取食和发育过程正常, 至幼虫期结束进入蛹期后, 蜜蜂的中后肠接通, 摄入的孢子伴随蜜蜂的蛹便进入后肠并在此迅速萌发生长, 在1~2 d内菌丝即突破体表, 导致蜜蜂死亡。蜜蜂球囊菌选择营养物质储存最多而防御能力较低的幼虫到蛹转化期侵染, 降低了侵染成本, 提高成功率。该研究阐明了蜜蜂球囊菌侵染蜜蜂的机制, 丰富了昆虫与病原菌之间相互作用的内容。     

结合三代测序与二代测序技术揭示蜜蜂球囊菌孢子转录组的复杂性

【目的】蜜蜂球囊菌(Ascosphaera apis)是一种专性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在利用PacBio单分子实时(singlemoleculereal-time,SMRT)测序技术对蜜蜂球囊菌孢子(AaS)中基因的可变剪切(alternative splicing,AS)和可变多聚腺苷酸化(alternative polyadenylation,APA)以及长链非编码RNA (long non-coding RNA,lncRNA)进行鉴定和分析,进而揭示蜜蜂球囊菌孢子中转录组的复杂性。【方法】采用Suppa软件对蜜蜂球囊菌孢子中基因的AS事件进行鉴定。通过RT-PCR对不同类型的AS事件进行验证。采用TAPIS pipeline对蜜蜂球囊菌孢子基因的APA位点进行鉴定。利用MEME软件分析孢子全长转录本的poly(A)剪接位点上游50bp的序列特征并鉴定motif。联用CPC和CNCI软件和比对Swiss-prot数据库的方法预测lncRNA,取三者的交集作为高可信度的lncRNA集合。进一步比较lncRNA和mRNA的转录本长度,外显子数量与长度,内含子长度,GC含...     

球囊菌胁迫中华蜜蜂幼虫肠道过程中病原的转录组学研究

【目的】本研究利用RNA-seq技术对球囊菌胁迫的中华蜜蜂(中蜂)幼虫肠道进行深度测序,经趋势分析得到差异表达基因(DEGs)的显著表达模式,进而对胁迫过程中的球囊菌进行转录组学分析。【方法】利用Illumina HiSeq 2500平台对球囊菌胁迫的中蜂幼虫肠道进行深度测序,并利用相关软件进行了深入分析。最后,通过RT-qPCR对RNA-seq数据进行了验证。【结果】本研究共得到球囊菌的41133932条高质量clean reads。22865个DEGs共聚类为8个基因表达模式,其中,16769个DEGs聚类为2个显著上调趋势与2个显著下调趋势。GO富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于40与37个GO term,基因富集数最多的为细胞进程(2486 unigenes)。KEGG代谢通路(pathway)富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于119和112个pathway,基因富集数最多的分别是氨基酸生物合成(127 unigenes)与核糖体(98 unigenes)。进一步分析表明球囊菌在胁迫中蜂幼虫肠道的过程中通过提高物质合成促进其增殖,而宿主通过抑制球囊菌的蛋白合成抵御病原入侵。富集在MAPK信号通路的11个DEGs的表达水平随着胁迫时间的延长而逐渐下降,推测中蜂幼虫通过抑制该通路而阻遏球囊菌增殖。【结论】本研究不仅为揭示白垩病过程中的球囊菌-中蜂幼虫互作提供了重要信息,也为阐明不同抗性蜂种的球囊菌抗性差异奠定了基础。     

蜜蜂球囊菌的microRNA鉴定及其调控网络分析

【目的】本研究利用small RNA-seq技术对球囊菌的纯培养进行测序,对球囊菌的micro RNAs miRNAs)进行预测、鉴定和分析,进而构建miRNAs-mRNAs的调控网络。【方法】利用Illumina Hiseq Xten平台对球囊菌菌丝与孢子进行测序,通过相关生物信息学软件对球囊菌的miRNAs进行预测和分析,通过茎环(Stem-loop)PCR对部分miRNAs进行鉴定,利用Cytoskype软件构建miRNAs-mRNAs的调控网络。【结果】本研究共获得48268696条clean reads,预测出118个球囊菌的miRNAs,它们的长度分布介于18–25 nt之间,不同长度的mi RNA的首位碱基偏好性差异明显。Stem-loop PCR验证结果显示共有10个miRNAs能够扩增出符合预期的目的片段,说明多数miRNAs可能真实存在。共预测出6529个球囊菌miRNAs的靶基因,其中5725个能够注释到Nr、Swissprot、KOG、GO和KEGG数据库。进一步分析结果显示有24个靶基因注释在MAPK信号通路。Cytoskype软件分析结果显示球囊菌的miRNAs与mRNAs之间存在复杂的调控网络,绝大多数的miRNAs处于调控网络的内部且同时结合多个mRNAs。【结论】本研究率先对球囊菌的miRNAs及miRNAs-mRNAs调控网络进行全面分析,研究结果丰富了对球囊菌miRNAs的认识,为其基础生物学信息提供了有益补充,也为阐明球囊菌致病的分子机理打下了一定基础。     

基于第三代长读段测序数据解析蜜蜂球囊菌基因的可变剪切与可变腺苷酸化

【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种专性侵染蜜蜂幼虫的真菌病原,导致的白垩病是严重影响养蜂生产的顽疾,每年给养蜂业造成较大损失。本研究旨在基于已获得的第三代长读段测序数据对球囊菌菌丝(Aam)和孢子(Aas)中基因的可变剪切(alternative splicing,AS)和可变多聚腺苷酸化(alternative polyadenylation,APA)进行深入分析。【方法】利用Astalavista软件鉴定Aam和Aas中基因的AS事件类型。利用IGV浏览器对部分剪切异构体(isoform)的结构进行可视化。通过TAPIS pipeline对Aam和Aas中基因的APA位点进行鉴定。通过MEME软件对Aam和Aas中全长转录本的APA位点上游50bp的序列特征进行分析并对motif进行鉴定。【结果】在Aam中共鉴定到286次AS事件,包括162次RI(Retained intron),87次A3(Alternative 3'splice-site)、32次A5(Alternative 5'splice-site)和5次SE (Skipping exon);在Aas中共鉴定到559次AS事件,包括305次RI、155次A3、85次A5、13次SE和1次MEE (Mutually exclusive exon)。进一步分析发现,现有参考基因组上的多数注释基因结构并不完整;部分注释基因在菌丝和孢子中转录形成的isoform在数量和结构方面均存在差异;部分isoform在参考基因组中没有对应的注释基因。Aam中共鉴定到2748个基因含有1个及以上的APA位点,其中含有1个APA位点的基因数量最多(726,26.42%);Aas中共鉴定到2768个基因含有1个及以上的APA位点,其中含有5个以上APA位点的基因数量最多(1180,42.63%)。部分基因在菌丝和孢子中含有不同的APA位点数。序列特征分析结果显示,球囊菌全长转录本的3'UTR的上下游表现出明显的碱基倾向性,U和A分别富集在3'UTR的上游和下游。此外,在球囊菌全长转录本的APA位点上游鉴定到4个motif,分别是UCUCCU、UCUUCU、CCCACC和CCCCCU。【结论】本研究通过对球囊菌菌丝和孢子中基因的AS和APA进行深入分析,揭示了球囊菌转录组的复杂性,为完善现有的基因组和转录组注释提供了宝贵信息,也为探究AS和APA在球囊菌的基因表达调控中的作用提供了关键基础。     

澜沧江流域北部中华蜜蜂有毒蜂蜜孢粉学和营养生态位分析

为分析澜沧江流域北部人食用蜂蜜中毒的原因,于2013年6—9月份,对该区域蜜蜂和蜜源植物的分布情况进行了调查,观察蜜蜂采集有毒蜜源植物的行为,并进一步调查了蜜蜂巢内蜂蜜、蜂花粉的储存情况,采集了中华蜜蜂蜂蜜样品,进行蜂蜜孢粉学与营养生态位分析。该区域有大量采用传统方式进行人工饲养的中华蜜蜂(Apis cerana cerana),及少量野生中华蜜蜂、黑色小蜜蜂(Apis audreniformis)、黑色大蜜蜂(Apis laboriosa smith)群体分布;人工饲养的中华蜜蜂蜂巢内部结构与野生中华蜜蜂蜂巢相似,为自然蜂巢,内有充足的蜜粉储存,部分蜂群蜂巢内虫害严重。该区域内主要蜜源植物为荞麦(Fagopyrum esculentum Moench),其他零星辅助蜜源较多,部分地点南烛(Vaccinium bracteatum Thunb)、昆明山海棠(Tripterygium hypoglaucum(Levl.)Hutch)连片集中分布。对中华蜜蜂蜂蜜进行孢粉学和营养生态位分析,结果表明:中华蜜蜂蜂蜜标本中含有有毒蜜源植物南烛、昆明山海棠花粉,部分样品中南烛、昆明山海棠的花粉含所占比例较高;中华蜜蜂的营养生态位宽度为0.22,比其他地区中华蜜蜂生态位指数小,推测澜沧江水电枢纽的修建等人为原因已对蜜蜂种类、蜜源植物的物种组成、群落结构造成了较大影响。     

澜沧江流域北部中华蜜蜂食源和营养生态位随海拔梯度的变化特征

为了解澜沧江流域北部中华蜜蜂(Apis cerana cerana)的分布,探究中华蜜蜂的食源和营养生态位沿海拔梯度的变化特征,调查了澜沧江流域北部各海拔区域中华蜜蜂的种群分布,运用蜂蜜孢粉学(melissopalynology)分析了各海拔区域中华蜜蜂蜂蜜中花粉的组成特征和变化规律,并综合分析了海拔、中华蜜蜂营养生态位和蜂蜜中花粉的相关性,探讨了自然环境与中华蜜蜂的分布,海拔梯度与蜜粉源植物、中华蜜蜂的食源和营养生态位的关系。结果表明:澜沧江流域北部中华蜜蜂的食源种类丰富,中华蜜蜂种群主要分布于2200—2800 m海拔区域。不同海拔区域中华蜜蜂采食花粉的种类和数量不同,中华蜜蜂种群分布多的海拔区域,蜂蜜中花粉的种类较多,但花粉数量相对少。随着海拔梯度升高,中华蜜蜂蜂蜜中花粉的数量表现为先降后增,而花粉的种类则表现为先增后降。不同海拔区域的中华蜜蜂营养生态位存在差异,推测各海拔区域蜜粉源植物分布、中华蜜蜂种内和种间授粉昆虫及食草动物等竞争因素不同,但各海拔梯度间的变化差异不显著。海拔与中华蜜蜂营养生态位呈正相关(r=0.051),相关性不显著;海拔与蜂蜜中花粉数量呈正相关(r=0.047),与蜂蜜中花粉种类呈正相关(r=0.144),相关性都不显著;中华蜜蜂营养生态位与蜂蜜中花粉的种类呈负相关(r=-0.305),相关性显著(P!0.05);与花粉的数量呈负相关(r=-0.064),相关性不显著。蜂蜜中花粉的数量与种类呈正相关(r=-0.303),且相关性显著(P!0.05)。