囊状幼虫病病毒侵染对中华蜜蜂营养和免疫反应的影响

【目的】囊状幼虫病病毒(sacbrood virus, SBV)是严重危害中华蜜蜂Apis cerana cerana蜂群健康和种群数量的病原微生物,但其对蜜蜂的致死机制不明。本研究旨在探究SBV对不同发育阶段中华蜜蜂营养代谢和免疫的影响。【方法】分别给中华蜜蜂2日龄幼虫和新羽化成虫饲喂SBV,逐日统计死亡蜜蜂数量,检测病毒对蜜蜂存活的影响;利用qPCR检测中华蜜蜂4日龄幼虫、预蛹以及10和20日龄成虫体内SBV RNA、营养代谢基因(ilp1, ilp2, hex110, hex70b, hex70c和vg)、先天性病毒免疫基因(rel, toll, apidaecin, abaecin, defensin, hymenoptaecin, jra, key和state92e)、细胞凋亡基因(atg7和LOC100577876)和抗RNA病毒基因(dis3和dicer)的表达水平。【结果】 SBV感染显著降低了中华蜜蜂幼虫的存活率,但对成虫的生存影响不明显。SBV RNA在中华蜜蜂4日龄幼虫和预蛹体内的表达量显著高于其在10和20日龄成虫体内的表达量。SBV显著降低了中华蜜蜂幼虫营养代谢基因ilp1, ilp2, hex110, hex70b和hex70c以及成虫营养代谢基因vg和hex110的表达量,但显著提高了4日龄幼虫rel, toll, apidaecin, abaecin, defensin, hymenoptaecin和jra以及成虫key和 state92e等先天性病毒免疫基因的表达量,还引起预蛹体内的细胞凋亡基因atg7和LOC100577876的表达量显著增加。【结论】SBV在中华蜜蜂幼虫和预蛹体内的感染水平远高于在成虫体内的,其对幼虫的危害也大于对成虫。SBV显著影响了中华蜜蜂正常的营养代谢,染毒中华蜜蜂能够提高自身的免疫水平来应对;预蛹期中华蜜蜂幼虫细胞凋亡水平的显著增加可能与化蛹异常及死亡有关。     

重庆中华蜜蜂囊状幼虫病病毒多聚蛋白基因分子特性分析

【目的】分析地处中国南方的重庆地区中华蜜蜂 Apis cerana cerana（简称“中蜂”）囊状幼虫病病毒（sacbrood virus, SBV）（AcSBV-S.Chn-CQ）的遗传特性。【方法】采用RT-PCR法对采自重庆23个区县的阳性中蜂幼虫进行SBV多聚蛋白基因扩增、序列测定和序列分析,采用邻接法（neighbor-joining method）基于基因组和基因组片段SB11-SB12构建系统进化树,使用RDP3检测SBV的基因重组。【结果】AcSBV-S.Chn-CQ与中国南方中蜂SBV（AcSBV-S.Chn）其他分离株、越南北部东方蜜蜂 A. cerana SBV（AcSBV-N.Vie）、越南北部西方蜜蜂 A. mellifera SBV（AmSBV-N.Vie）、韩国东方蜜蜂SBV（AcSBV-S.Kor）的同源性较高。AcSBV-S.Chn-CQ在结构蛋白编码区域存在连续51个核苷酸缺失,在非结构蛋白编码区域有3个不连续的核苷酸缺失,与AcSBV-N.Vie, AmSBV-N.Vie和AcSBV-S.Kor核苷酸缺失相同。AcSBV-S.Chn-CQ与AcSBV-S.Chn-FZ, AcSBV-N.Vie, AmSBV-N.Vie和AcSBV-S.Kor在亚洲基因型中形成一亚型,亚型内检测到重组事件。【结论】推测AcSBV-S.Chn-CQ与SBV-N.Vie和AcSBV-S.Kor可能来自一个共同祖先。     

幼虫信息素中三种酯类对中华蜜蜂和意大利蜜蜂工蜂哺育和封盖行为以及蜂王发育影响

为研究蜜蜂幼虫信息素中3种酯类成分（甲基棕榈酸酯、乙基棕榈酸酯和乙基油酸酯）对中华蜜蜂Apis cerana cerana和意大利蜜蜂Apis mellifera ligustica工蜂哺育行为、封盖行为以及蜂王发育影响, 在人造蜂蜡王台中加入1%和0.1%（w/w）的3种酯类作为实验组, 以不添加酯类（0%）的为对照组, 移入1日龄工蜂幼虫, 测定王台接受率、单个王台中幼虫和王浆重量; 另将分别添加1%和0.1%（w/w）3种酯类的石蜡假幼虫放入工蜂巢房中, 同样设对照组, 测定假幼虫的封盖率; 在新鲜王浆中以1%和0.1%（w/w）分别加入3种酯类作为实验组, 以不添加酯类（0%）的作为对照组, 再分别在1日龄、2日龄和3日龄幼虫王台中加入0.01 mL含有酯类的蜂王浆, 并测定蜂王初生重和卵巢管数量。结果表明: 0.1%甲基棕榈酸酯可以显著提高中蜂和意蜂幼虫重量; 意蜂的甲基棕榈酸酯和乙基油酸酯两个实验组(1.0%, 0.1%)假幼虫封盖率都极显著高于对照组; 乙基油酸酯两个实验组(1.0%, 0.1%)都显著降低了中蜂和意蜂蜂王初生重和卵巢管数量。这说明不同蜜蜂幼虫信息素具有不同的生物学效应。     

中华蜜蜂mrjp1 cDNA的克隆及其序列分析

构建中华蜜蜂(Apis cerana cerana)8日龄工蜂头部cDNA文库,利用中蜂基因组的mrjp3部分基因片段作为杂交探针,采用DIG标记筛选cDNA文库,获得mrjps阳性克隆120个;对阳性克隆进行PCR扩增和测序,通过NCBI的BLAST序列比对,获得12个与印度蜂(Apis cerana india)、西方蜜蜂(Apis mellifera L.)mrjp1基因同源的中蜂mtjp1 cDNA片段,并进一步对中华蜜蜂mrjp1的cDNA全序列进行测定和分析。序列比对分析表明,东方蜜蜂(Apis cerana)与西方蜜蜂mrjp1的cDNA序列相似性为93．78％,中华蜜蜂与印度蜂的相似性高达99．36％,这一结果从分子水平证实中华蜜蜂与印度蜂有较近的共同祖先,而东方蜜蜂与西方蜜蜂的亲缘关系较远。     

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

蜂王上颚腺信息素对中华蜜蜂、意大利蜜蜂雄蜂选择行为的影响

【目的】分析蜂王上颚腺信息素(Queen mandibular pheromone,QMP)对中华蜜蜂Apis cerana cerana(简称中蜂)雄蜂与意大利蜜蜂Apismelliferaligustica(简称意蜂)雄蜂行为反应的变化,探索两蜂种雄蜂之间婚飞的干扰机理。【方法】本试验通过室内Y型嗅觉仪检测QMP及其主要成分反式-9-氧代-2-癸烯酸[(E)-9-oxodec-2-enoicacid,9-ODA]对飞行、爬行状态下中蜂雄蜂与意蜂雄蜂选择行为的差异进行研究。【结果】飞行或爬行状态下的中蜂雄蜂和意蜂雄蜂,均对3.5、7.0和14.0μg/μL的9-ODA没有显著趋向性反应(P 0.05);飞行中蜂雄蜂、飞行与爬行意蜂雄蜂均对0.04、0.2、1.0及7.0μg/μL的QMP具有显著趋避反应(P 0.05),而QMP对爬行中蜂雄蜂无显著影响(P 0.05);在中蜂雄蜂和意蜂雄蜂共存时的试验中发现:飞行、爬行意蜂雄蜂均对7.0μg/μL QMP存在显著趋避反应(P 0.05),而7.0μg/μL QMP对飞行或爬行中蜂雄蜂均无显著影响(P 0.05)。【结论】在室内环境下,QMP对中蜂、意蜂雄蜂有驱避作用。     

安徽省七种蜜蜂病毒的发生与流行研究

【目的】调查安徽省内7种常见蜜蜂病毒:蜜蜂畸翅病毒(Deformed wing virus,DWV)、以色列急性麻痹病毒(Israeli acute paralysis virus,IAPV)、急性蜜蜂麻痹病毒(Acute bee paralysis virus,ABPV)、慢性麻痹病毒(Chronic bee paralysis virus,CBPV)、黑蜂王台病毒(Black queen cell virus,BQCV)、囊状幼虫病病毒(Sacbrood virus,SBV)、克什米尔病毒(Kashmir bee virus,KBV)的感染发生情况,为安徽养蜂业可持续健康发展提供理论依据。【方法】运用反转录RT-PCR和序列分析比对的方法对安徽省内21个乡镇中的38个蜂场蜜蜂样品进行研究分析,以获得以上7种蜜蜂病毒的特异性发生情况。【结果】意大利蜜蜂Apis mellifera蜂场感染率:DWV(64%),IAPV(43%),CBPV(32%),ABPV(14%),BQCV(11%);中华蜜蜂Apis cerana蜂场感染率:DWV(80%),IAPV(40%),CBPV(30%),ABPV(10%),BQCV(0)。SBV和KBV在所有的蜜蜂样品中均未检测到。【结论】DWV,IAPV,CBPV,ABPV,BQCV在安徽省内大范围都存在发生流行现象,SBV和KBV对安徽蜜蜂的潜在危害可能性小。     

中华蜜蜂Orco嗅觉受体基因的克隆、表达及亚细胞定位

【目的】克隆鉴定中华蜜蜂Apis cerana cerana的Orco嗅觉受体基因, 并对其在工蜂触角上进行免疫荧光定位。【方法】利用RT-PCR技术克隆中华蜜蜂Orco基因, 并对其编码的氨基酸序列进行生物信息学分析, 使用Real-time PCR技术鉴定其在中华蜜蜂不同发育时期及不同组织的表达谱; 利用免疫荧光定位技术在中华蜜蜂工蜂触角中对Orco进行亚细胞定位。【结果】获得中华蜜蜂Orco基因的全长cDNA序列, 命名为AcerOrco (GenBank登录号： JF968610.1), 其全长为1 434 bp, 编码477个氨基酸, 预测其含7个跨膜结构以及4个位于细胞膜外的亲水区。表达谱分析显示, AcerOrco在卵、 幼虫和蛹期呈低丰度表达, 1日龄及内勤蜂时期主要在触角和足中表达, 且在1日龄的触角中表达量最高; 采集蜂时期的触角、 头(去除触角)、 胸、 腹和翅中均有较高丰度的表达。亚细胞定位结果显示, AcerOrco不仅在采集蜂触角鞭节上大量表达（尤其在触角鞭节第1亚节中表达量较高）, 而且常成对出现, 并且发现AcerOrco可能主要在触角毛形感器的外部神经元(outer dendrite, OD)以及板形感器的树突神经元中表达。【结论】成功克隆了AcerOrco基因全长, 获得了其表达谱, 且将其定位于工蜂采集蜂的触角感器神经元上, 最终推测AcerOrco与中蜂嗅觉发育和触角感器功能密切相关。     

中华蜜蜂幼虫原代细胞培养及中华蜜蜂囊状幼虫病毒在培养的细胞中的复制

【目的】建立一种适用于蜜蜂源的细胞培养方法,为蜜蜂源细胞培养和蜂病毒的研究奠定基础。【方法】比较在Grace和WH2两种昆虫细胞培养基中培养的中华蜜蜂 Apis cerana cerana 幼虫原代细胞状况,筛选出适用于中华蜜蜂幼虫原代细胞培养的最佳培养基,并通过细胞活力比较,确定用于中华蜜蜂幼虫原代细胞培养的适宜胎牛血清(fetal bovine serum, FBS)浓度。在此基础上,用该原代细胞接种中华蜜蜂囊状幼虫病毒（Chinese sacbrood bee disease, CSBV）,并通过实时定量RT-PCR方法对病毒复制情况进行检测。【结果】相对于WH2培养基,在Grace培养基中生长的细胞大而圆,透明,边缘整齐,无颗粒物,活力明显高于WH2培养,且含15% FBS的Grace培养基更适合于中华蜜蜂幼虫原代细胞培养。CSBV接种在该原代细胞,能够复制增殖,同时伴随宿主细胞快速分裂。【结论】中华蜜蜂幼虫原代细胞培养基在含15%FBS的Grace中能够良好生长,并且CSBV可以在该原代细胞中进行复制。     

意大利蜜蜂幼虫肠道的miRNAs的生物信息学预测及分析

【目的】蜜蜂肠道是食物消化和营养吸收的主要部位,同时也是抵御病原侵染的重要场所。本研究旨在对意大利蜜蜂Apismelliferaligustica(简称意蜂)幼虫肠道的microRNAs(miRNAs)及其靶基因进行深入分析,进而解析miRNAs在幼虫肠道发育和生长过程中的作用。【方法】利用small RNA-seq(sRNA-seq)技术对意蜂4、5和6日龄幼虫肠道进行测序,通过相关生物信息学软件对意蜂幼虫肠道的miRNAs进行预测及分析。利用TargetFinder软件预测miRNAs的靶基因,然后将靶基因通过BLAST软件进行GO和KEGG数据库的功能注释。利用Cytoscape软件构建miRNAs与其靶向结合的mRNAs的调控网络。【结果】意蜂幼虫肠道样品的测序共获得96 329 456条有效标签序列(Clean tags),预测出560个miRNAs,包括45个novelmiRNAs。上述miRNAs的长度主要分布在17-27nt之间,且不同长度的miRNAs的首位碱基偏向性具有明显差异。表达量聚类分析结果显示,331个miRNAs为4、5和6日龄幼虫肠道所共有且表达量较为稳定,随着发育时间延长,特有miRNAs的表达水平呈总体增加的趋势。利用TargetFinder软件共预测出16479个意蜂幼虫肠道的靶基因,其中有8132和3361个可分别注释到GO和KEGG数据库,进一步分析发现有224个靶基因注释在Wnt信号通路,分别有2个和27个靶基因与保幼激素和蜕皮激素的调控密切相关。【结论】本研究在全基因组水平对意蜂幼虫肠道的miRNAs进行预测和分析,研究结果不仅丰富了意蜂的miRNAs信息,也为深入研究意蜂幼虫肠道的生长发育机理打下了初步基础。     

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

【目的】本研究利用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应答分析

【目的】蜜蜂球囊菌(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有望作为治疗白垩病的分子靶标。     

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

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

[目的] 本研究旨在探究长链非编码RNA（long non-coding RNA,lncRNA）在中华蜜蜂（Apis cerana cerana,简称中蜂）6日龄幼虫应答蜜蜂球囊菌（Ascosphaera apis,简称球囊菌）侵染过程中的差异表达谱及调控作用。[方法] 利用链特异性cDNA建库的RNA-seq技术对未被侵染及球囊菌侵染的中蜂6日龄幼虫肠道（AcCK和AcT）进行深度测序。通过相关生物信息学软件分析lncRNA的结构特征和表达谱。筛选并分析差异表达lncRNA（differentially expressed lncRNA,DElncRNA）的顺式（cis）作用及竞争性内源RNA（competing endogenous RNA,ceRNA）调控网络。采用RT-qPCR验证测序数据及DElncRNA差异变化趋势的可靠性。[结果] AcCK和AcT共预测出642个已知lncRNA和487个新lncRNA。与蛋白编码基因相比,上述中蜂lncRNA外显子数更少、长度更短且表达量更低。43个antisense lncRNA与40个正义链mRNA之间互补配对。AcCK和AcT比较组包含367个上调lncRNA和268个下调lncRNA。有194个DElncRNA潜在调控461个上下游基因,并涉及细胞进程、代谢进程和催化活性等38个功能条目以及氨基酸代谢、内吞作用和MAPK等191条通路。此外,有180个DElncRNA可靶向结合50个DEmiRNA,进而调控6365个mRNA;三者之间形成较为复杂的ceRNA调控网络。[结论] 中蜂6日龄幼虫肠道的部分lncRNA可作为antisense lncRNA参与应答球囊菌侵染;部分DElncRNA可通过cis作用调节物质代谢和免疫途径相关的上下游基因,从而介导宿主的侵染应答;TCONS_00010661和TCONS_00003104等DElncRNA可通过ceRNA网络调控Jak-STAT和氧化磷酸化等通路及富集基因,进而参与宿主的侵染应答。     

侵染中华蜜蜂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有望作为白垩病治疗的新型分子靶点。     

意大利蜜蜂工蜂幼虫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影响意蜂工蜂幼虫体重及宿主对蜜蜂球囊菌胁迫的免疫应答。     

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

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

茶树花香气及茶叶气味对中华蜜蜂的引诱效应

我国茶园面积约占全球的60%,茶树花已成为我国主要植物蜜源之一,尤其在秋冬季。茶树花招引蜜蜂的机理尚不清晰,遂以中华蜜蜂Apis cerana cerana为试虫,以8个国家级良种茶树的鲜花、茶树花主要香气组分和茶鲜叶重要挥发性组分为味源,用Y形嗅觉仪进行行为测定,结果表明:①在0.25—5.00 g质量范围内,不仅8个良种茶树的鲜花明显引诱该蜂,而且每种茶树花质量为某一确定值时,其引诱的蜂数与CK(洁净空气)引诱的蜂数差异显著(P0.05);②在茶树花27个主要香气组分中,发现苯乙酮(10~(-6) g/mL)、芳樟醇(10~(-6)、10~(-2) g/mL)、莰烯(10~(-4) g/mL)、顺-3-己烯-1-醇(10~(-10) g/mL)、α-法尼烯(10~(-6)、10~(-4)、10~(-2) g/mL)、癸醛(10~(-6) g/mL)、β-紫罗兰酮(10~(-6) g/mL)、亚油酸(10~(-4)、10~(-2) g/mL)、顺-2-戊烯-1-醇(10~(-4) g/mL)和苯甲醛(10~(-2) g/mL)显著引诱该蜂(P0.05),而水杨酸甲酯(10~(-2) g/mL)、橙花醇(10~(-2) g/mL)和辛醛(10~(-6) g/mL)排斥该蜂(P0.05);③在茶鲜叶24个重要挥发性组分中,发现十八醇(10~(-4) g/mL)、正辛醇(10~(-2) g/mL)、吲哚(10~(-6) g/mL)、柠檬醛(10~(-2) g/mL)、薄荷醇(10~(-6) g/mL)、β-石竹烯(10~(-4) g/mL)、薄荷酮(10~(-4)、 10~(-2) g/mL)显著引诱该蜂(P0.05),而邻甲苯酚(10~(-2) g/mL)、香茅醇(10~(-6) g/mL)、1,3-丙二醇(10~(-2) g/mL)、α-松油烯(10~(-6) g/mL)和正戊醇(10~(-4) g/mL)排斥该蜂(P0.05)。认为:茶树花气味及其主要香气组分和茶鲜叶主要挥发性组分显著引诱中华蜜蜂,尤其是茶树花香气组分中含量分别超过20%的苯乙酮和芳樟醇强烈引诱该蜂,导致花季时期广袤的茶园及其中的茶树花对该蜂具有强大的引诱效应。     

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