侵染中华蜜蜂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...     

温度、相对湿度和pH对蜜蜂球囊菌孢子萌发的影响

研究了温度、相对湿度和pH对蜜蜂球囊菌(Ascosphaeraapis)孢子萌发3个阶段(活化、膨大、产生萌发管)的影响．结果表明,孢子活化和膨大在15～40和(25～40)±0．5℃范围内受温度的影响不明显(P＞0．05);萌发管仅发生在25～37±0．5℃,最适温度位于(31～35)±0．5℃．相对湿度越大,越有利于孢子萌发,而相对湿度低于80％对孢子萌发极为不利．孢子萌发的3个阶段在pH为5～7．8时几乎不受pH变化的影响,而在pH值较低影响很大．可见,A．apis是一种高度专一的蜜蜂幼虫病原体．     

球囊菌侵染对意大利蜜蜂幼虫肠道免疫与解毒相关基因表达及蛋白质、脂质和糖类含量的影响

[目的]球囊菌Ascosphaera apis是一种真菌性病原体,可以侵染意大利蜜蜂Apis mellifera ligustica幼虫,导致蜂群患白垩病.本研究通过球囊菌侵染意大利蜜蜂幼虫,探究对意大利蜜蜂幼虫肠道免疫解毒相关基因表达及蛋白质、脂质和糖类含量的影响,分析球囊菌侵染胁迫下意大利蜜蜂免疫应答反应与营养物质...     

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...     

意大利蜜蜂幼虫肠道内球囊菌及其纯培养的高表达基因差异分析

【背景】球囊菌是一种典型的蜜蜂真菌性病原,特异性地侵染蜜蜂幼虫。目前,有关球囊菌在侵染过程中的基因表达规律的信息极为有限。【目的】通过深入分析胁迫意大利蜜蜂(意蜂)幼虫肠道的球囊菌及其纯培养的高表达基因(HEGs)差异,探索球囊菌在胁迫意蜂幼虫肠道后期与胁迫前的基因表达规律。【方法】利用RNA-Seq技术对球囊菌胁迫的意蜂6日龄幼虫肠道(Aam T)及球囊菌纯培养(AaCK)进行深度测序,根据FPKM值筛选得到球囊菌的HEGs,进而通过GO(Gene ontology)及KEGG(Kyoto Encyclopedia of Genes and Genomes)富集分析、Venn分析对上述HEGs进行功能预测和生物学意义挖掘。【结果】AaCK与Aam T的转录组测序共得到105 447 578条原始读段(Raw reads),经过滤得到88 466 344条有效读段(Clean reads),两端平均Q20为97.50%,平均Q30为93.81%。GO富集分析结果显示,AaCK的HEGs富集于26个GO terms,基因富集数最多的是细胞(22 Unigenes),其次是细胞组件(22 Unigenes)和代谢过程(21 Unigenes);Aam T的HEGs富集于22个GO terms,基因富集数最多的是催化活性(23 Unigenes),其次是细胞进程(18 Unigenes)和代谢过程(18 Unigenes)。KEGG代谢通路(Pathway)富集分析显示,AaCK的HEGs富集在109个Pathways上,基因富集数最多的是核糖体(179 Unigenes),其次是氨基酸生物合成(70Unigenes)和碳代谢(62 Unigenes);Aam T的HEGs富集在114个Pathways上,基因富集数量最多的是核糖体(178 Unigenes),其次是碳代谢(116 Unigenes)和氧化磷酸化(112 Unigenes)。Venn分析结果显示,AaCK与Aam T共有的HEGs有260个,二者特有的HEGs分别为2 161个和4 445个。【结论】提供了胁迫意蜂6日龄幼虫肠道的球囊菌及其纯培养的HEGs表达谱,揭示了球囊菌在胁迫前和胁迫后期的基因表达规律,也为阐明球囊菌致病的分子机理提供了有益的信息和线索。     

环状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)...     

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

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

意大利蜜蜂幼虫肠道在球囊菌胁迫后期的差异表达微小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可参与到意蜂幼虫肠道与球囊菌间的互作。     

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

【目的】本研究利用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的表达水平随着胁迫时间的延长而逐渐下降,推测中蜂幼虫通过抑制该通路而阻遏球囊菌增殖。【结论】本研究不仅为揭示白垩病过程中的球囊菌-中蜂幼虫互作提供了重要信息,也为阐明不同抗性蜂种的球囊菌抗性差异奠定了基础。     

中华蜜蜂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在宿主的胁迫应答过程中可能通过调控氧化磷酸化、细胞和体液免疫等通路发挥特殊作用。研究结果为深入理解中蜂幼虫对球囊菌的胁迫应答机制及二者的相互作用机制提供了新见解。     

吡虫啉对成年意大利蜜蜂脑神经细胞致凋亡作用

【目的】蜜蜂是重要的授粉昆虫, 但一直受到新烟碱类杀虫剂（如吡虫啉）的危害, 该类杀虫剂主要作用于蜜蜂脑神经细胞。本研究旨在明确吡虫啉对意大利蜜蜂Apis mellifera ligustica脑神经细胞致凋亡作用。【方法】利用亚致死剂量吡虫啉饲喂成年意大利蜜蜂, 通过原位末端标记法（TUNEL）检测脑神经细胞的凋亡情况, 利用免疫荧光法检测Caspase-1激活情况, 通过透射电镜观察脑神经细胞的超微结构。【结果】在对选取的脑部7个部位进行总体观测后, 摄入亚致死剂量的吡虫啉（9.90 ng/蜜蜂）诱导成年意蜂工蜂脑神经细胞凋亡率随时间递增而增加, 在处理9 d和12 d时, 处理组与空白对照组差异极显著（P<0.01）; Caspase-1阳性细胞率也随时间递增而增加, 在处理3, 6, 9和12 d时, 处理组与空白对照组差异极显著（P<0.01）。凋亡率和Caspase-1阳性细胞率均表现出时间效应, 且两者成正相关性。超微结构表明凋亡的神经细胞呈现凋亡和自噬双重特征, 包括： 细胞固缩、 染色质浓缩、 凋亡小体和自噬体出现; 线粒体肿胀, 有的被自噬泡包裹后进行线粒体自噬。【结论】本研究从细胞水平证实了亚致死剂量吡虫啉对成年意蜂工蜂脑神经细胞具有致凋亡作用, 并且其凋亡途径与Caspase-1和细胞自噬有关, 这为利用细胞凋亡法评价杀虫剂对非靶标生物的神经毒性提供一定的理论基础。     

意大利蜜蜂和小峰熊蜂在温室桃园的传粉行为及其影响因素

意大利蜜蜂Apis mellifera ligustica和小峰熊蜂Bombus hypocrita是我国设施农业中的两种主要授粉昆虫, 二者的传粉方式和传粉效率因不同作物而各有不同。为了选用最理想的蜂种为温室作物授粉, 提高作物的授粉效率, 我们于2008-2010年连续3年在北京的温室桃园进行了意大利蜜蜂和小峰熊蜂的传粉行为及其影响因素的研究。结果表明, 两种蜂的日活动规律不同, 和意大利蜜蜂相比, 小峰熊蜂的活动起点温度低、 日工作时间长、 单花访问时间长, 采集蜂在温室内距离蜂巢不同距离之间扩散比较均匀。而意大利蜜蜂采集蜂的数量随授粉半径的增大而明显减少。温室内环境因子对两种蜂传粉活动的影响基本一致, 温度对两种蜂的传粉活动影响最大, 其次是湿度、 单花花蜜浓度和光照强度, 而单花花蜜体积对蜂活动无明显影响。本研究认为对于温室桃授粉应优先选用授粉性能稳定的小峰熊蜂, 并且适当调节温室内环境条件, 以提高其授粉效率。     

意蜂和中蜂四种同工酶的研究

用聚丙烯酰胺凝胶等电聚焦电泳分析了意蜂和中蜂的酯酶(Est)、异柠檬酸脱氢酶(Idh)、苹果酸酶(Me)和苹果酸脱氢酶(Mdh)同工酶.两个蜂种的四种同工酶谱有不同程度的差别、意蜂酯酶Ⅳ和苹果酸脱氢酶Ⅲ是多态性的;中蜂的四种同工酶没有多态现象.     

中华蜜蜂dynactin p62基因的克隆及不同发育阶段表达分析

为探究中华蜜蜂Apis cerana cerana dynactin p62基因的表达特性, 本研究克隆了中华蜜蜂dynactin p62的基因组DNA序列(GenBank登录号： JX101463) 和mRNA序列（GenBank登录号： JX101464）; 采用荧光定量PCR检测了中华蜜蜂dynactin p62在不同发育时期（3日龄和6日龄幼虫、 刚羽化出房蜜蜂）三型蜂中mRNA的表达量。结果表明： 该基因基因组DNA序列全长为2 403 bp, mRNA序列全长为1 491 bp, 编码496个氨基酸残基, 预测的蛋白分子量为56.49 kD, 等电点为8.31。系统发育分析表明中华蜜蜂dynactin p62与西方蜜蜂Apis mellifera dynactin p62聚成一支。该基因在不同发育时期均有表达, 在雌性蜜蜂（蜂王和工蜂）中, 刚羽化成虫期的表达量显著高于幼虫期(P<0.05), 并且同一发育时期相比, 工蜂的表达量显著高于蜂王(P<0.05); 而该基因在雄蜂中表达量没有明显的规律性。这些结果提示该基因可能与中华蜜蜂级型分化有关。     

意大利蜜蜂刚出房与性成熟雄蜂的蛋白质组差异分析

【目的】对意大利蜜蜂Apis mellifera ligustica刚出房雄蜂与性成熟雄蜂的蛋白质组进行比较,探讨雄蜂在性成熟过程中蛋白质表达变化,为进一步研究雄蜂发育生物学获得差异表达蛋白质方面的依据。【方法】采用双向电泳法建立意大利蜜蜂雄蜂发育过程中刚出房时与性成熟期的蛋白质表达谱,通过质谱分析与数据库检索,鉴定部分差异蛋白。【结果】在意大利蜜蜂刚出房雄蜂和性成熟雄蜂中分别检测到2 490和2 317个蛋白点,其中差异表达蛋白点有157个。在刚出房雄蜂中高度表达的蛋白点有102个;在性成熟雄蜂中高度表达的蛋白点有55个。对部分差异蛋白进行质谱分析,共鉴定了18个蛋白点,其中在刚出房雄蜂中上调表达的蛋白有肌钙蛋白、SEC13蛋白、DJ蛋白等,在性成熟雄蜂中上调表达的蛋白有副肌球蛋白、精氨酸激酶、肌动蛋白解聚因子等。【结论】意大利蜜蜂雄蜂在性成熟发育过程中,其体内大量蛋白表达发生了变化,其差异表达的蛋白质可能与骨骼、飞行肌以及精子发育等机能有关。     

中华蜜蜂咽下神经节的结构和胚后发育

咽下神经节是昆虫腹神经索的第一个复合神经节, 主要调节口器附肢的活动。本研究通过形态解剖、 BrdU免疫组织化学等技术, 对中华蜜蜂Apis cerana cerana咽下神经节的组织结构和胚后发育过程进行了比较研究。结果表明： 中华蜜蜂的咽下神经节由上颚、 下颚和下唇3个神经节组成。在胚后发育过程中, 细胞增殖的活跃期主要集中在预蛹和蛹发育的第1天, 增殖活动一直持续到蛹发育的第4天结束。根据神经胶质细胞的位置和形态, 咽下神经节中的神经胶质可分为3种类型--表面神经胶质、 皮层神经胶质和神经纤维网神经胶质。本研究为蜜蜂神经系统的发育和功能研究提供了理论基础。