摘要: 【目的】意大利蜜蜂Apis mellifera ligustica是生产性能优良的西方蜜蜂Apis mellifera亚种。本研究旨在利用分子生物学手段揭示3种微小RNA(microRNA, miRNA) ame-miR-13b, ame-miR-100和ame-miR-bantam调控蛹期变态发育的分子机理。【方法】通过stem-loop RT-PCR验证ame-miR-13b, ame-miR-100和ame-miR-bantam在意大利蜜蜂工蜂蛹期的真实表达。通过RT-qPCR检测ame-miR-13b, ame-miR-100, ame-miR-bantam和靶mRNA在意大利蜜蜂工蜂蛹期发育过程的表达谱。利用相关生物信息学软件预测在意大利蜜蜂工蜂蛹期ame-miR-13b, ame-miR-100和ame-miR-bantam的靶mRNA,构建和分析miRNA-mRNA调控网络以及对靶mRNA进行数据库注释。【结果】ame-miR-13b, ame-miR-100和ame-miR-bantam在意大利蜜蜂工蜂蛹期发育过程真实表达且表达量总体均呈动态上升趋势;可分别靶向850, 136和506条mRNA,并与其靶mRNA之间形成较复杂的调控网络;这些靶mRNA可分别被注释到32, 24和33个GO条目,以及204, 114和171条KEGG通路。【结论】ame-miR-13b, ame-miR-100和ame-miR-bantam潜在通过调控蜕皮激素基因、yellow基因以及Hippo, FoxO, Wnt, Jak-STAT和Notch信号通路相关基因的表达共同影响意大利蜜蜂工蜂蛹期的变态发育过程。

摘要: 【目的】本研究旨在利用已获得的PacBio单分子实时(single molecule real-time, SMRT)测序数据对蜜蜂球囊菌Ascosphaera apis菌丝(AaM)和孢子(AaS)中的转录因子(TF)、融合基因和RNA编辑事件进行鉴定和分析,以期丰富蜜蜂球囊菌的相关信息,并为进一步探究它们的功能提供理论依据。【方法】利用BLASTx工具将AaM和AaS的全长转录本序列比对到Nr, Swiss-Prot和KEGG数据库以获得一致性最高的蛋白序列,再利用hmmscan软件将上述蛋白序列比对到Plant TFdb数据库以获得TF的分类及注释信息。采用TOFU软件中的fusion_finder.py程序进行融合基因的预测,进而分析融合基因的序列和位置信息。使用SAMtools预测AaM和AaS中的RNA编辑事件,再利用ANNOVAR软件对RNA编辑事件进行注释,进而采用相关生物信息学软件对RNA编辑位点基因进行GO功能和KEGG通路注释。【结果】在AaS中共鉴定到17个TF家族的213个TF,其中C2H2家族包含的TF成员最多。在AaM和AaS中分别鉴定到921和510个融合基因,二者共有的融合基因为510个,特有的融合基因分别为411和0个。在AaM和AaS中分别鉴定到547和191次RNA编辑事件,其中AaM中同义单核苷酸突变的数量最多,AaS中非同义单核苷酸突变的数量最多。此外,在AaM中鉴定到12种碱基替换类型,其中发生C->T的RNA编辑事件数量最多,达到158次;在AaS中鉴定到9种碱基替换类型,其中发生C->T和G->T的RNA编辑事件数量最多,均有42次。AaM和AaS中RNA编辑位点基因分别涉及19和24个GO功能条目;此外还能注释到11和20条KEGG通路。【结论】蜜蜂球囊菌的菌丝和孢子中含有丰富的TF、融合基因和RNA编辑位点;转录因子C2H2家族与蜜蜂球囊菌菌丝和孢子的生长发育和细胞活动具有潜在关联;RNA编辑事件的碱基替换类型在蜜蜂球囊菌和其他物种中具有物种特异性;RNA编辑可能在蜜蜂球囊菌菌丝和孢子的生长和代谢中发挥作用。

摘要: 【目的】当前全球气候变化、土地利用改变、人类活动加剧等正威胁着传粉昆虫的多样性及分布;蜜蜂是生态系统中重要的传粉昆虫类群,对气候、环境变化响应敏感。本研究以重要的传粉昆虫大蜜蜂Apis dorsata为对象,探讨全球变化格局下其潜在适生区变化以及影响其分布的关键因子。【方法】通过文献、馆藏和野外调查系统收集了全球范围内大蜜蜂的物种分布数据,使用13个环境变量通过MaxEnt模型模拟了大蜜蜂当前的潜在适生区;使用9个气候变量并结合公共地球系统模型(CCSM4)模拟了大蜜蜂过去、当前和未来的潜在适生区。【结果】AUC比率显示MaxEnt模型对大蜜蜂的潜在适生区模拟具有较高的准确性,模拟结果表明大蜜蜂的中高潜在适生区主要分布在南亚和东南亚湿润的热带雨林、热带季节性雨林和低地雨林。人类影响、温度季节性变化、等温性、最冷季均温和海拔是影响大蜜蜂潜在适生区的5个最主要因子;在人类影响下大蜜蜂的潜在适生区向山区和连片的湿润常绿森林区收缩,中高潜在适生区显著减少且呈破碎化趋势。基于9个气候变量和CCSM4气候模型对过去、当前和未来的模拟结果显示:在过去的末次冰盛期,东南亚地区可能是大蜜蜂的避难所;在未来,广布于热带地区的大蜜蜂适生区与当前的相近,且部分地区适生指数升高。【结论】基于气候的模拟结果显示大蜜蜂能积极应对未来气候变暖,但随人类活动的加剧及全球气候变化,大蜜蜂仍然面临较大的威胁,需要加强其在南亚和东南亚的中高潜在适生区的重视和保护。

摘要: 【目的】微孢子虫(Microsporidia)孢壁在孢子构成及孢子侵染宿主过程中扮演重要的角色。本研究旨在鉴定获得的东方蜜蜂微孢子虫Nosema ceranae新型孢壁蛋白,并进行基因克隆和原核表达,明确其亚细胞定位。【方法】通过在线软件对东方蜜蜂微孢子虫新型孢壁蛋白AAJ76_1400036761序列进行生物信息学分析。利用PCR法获取目的片段并将其克隆至原核表达载体pCold II中,利用IPTG诱导表达重组蛋白并通过镍柱亲和层析法纯化目的蛋白。以获得的重组蛋白为抗原免疫小鼠制备多克隆抗体,通过间接免疫荧光技术和免疫胶体金定位技术对该蛋白进行亚细胞定位分析;利用蛋白质免疫印迹法检测该蛋白与东方蜜蜂微孢子虫几丁质壳的互作。【结果】在MicrosporidiaDB数据库中获得AAJ76_1400036761基因序列,基因全长681 bp,编码226个氨基酸;预测等电点为6.84,分子量为26.19 kD。SDS-PAGE电泳和Western blot结果表明AAJ76__1400036761重组蛋白能够在大肠杆菌Eescherichia coli Rosetta中高量表达。Western blot结果表明,制备的多克隆抗体能够特异地识别东方蜜蜂微孢子虫总蛋白中的AAJ76__1400036761,说明其在成熟东方蜜蜂微孢子虫中有表达。亚细胞定位结果显示,AAJ76__1400036761定位于东方蜜蜂微孢子虫孢壁上。重组蛋白AAJ76_1400036761能够与蜜蜂微孢子虫的几丁质壳结合。【结论】AAJ76_1400036761蛋白在东方蜜蜂微孢子虫成熟孢子中有表达;该蛋白定位于东方蜜蜂微孢子虫孢壁上,为东方蜜蜂微孢子虫新的孢壁蛋白。本研究为深入研究该蛋白的生物学功能奠定了基础。

摘要: 【目的】分析吡虫啉处理对意大利蜜蜂Apis mellifera ligustica嗅觉学习行为及脑部基因转录的影响,为新烟碱类杀虫剂对蜜蜂的负面影响提供依据。【方法】实验室条件下一次性饲喂意大利蜜蜂成年工蜂含有4 ng吡虫啉的50%蔗糖溶液,以饲喂不含吡虫啉的50%蔗糖溶液为对照,通过伸吻反应(proboscis extension response, PER)行为实验测定其对意大利蜜蜂成年工蜂嗅觉学习行为的影响;收集上述测试的意大利蜜蜂工蜂脑部提取RNA,进行RNA-Seq测序和生物信息学分析,利用实时荧光定量PCR挑选6个差异表达基因(DEGs)检测其在脑部的表达量以验证RNA-Seq测序结果。【结果】一次性饲喂意大利蜜蜂成年工蜂含4 ng吡虫啉的50%蔗糖溶液后,意大利蜜蜂的嗅觉学习能力与对照组(饲喂50%蔗糖溶液)相比显著降低。RNA-Seq测序结果显示饲喂意大利蜜蜂含4 ng吡虫啉蔗糖溶液后,处理组与对照组之间共有123个DEGs[校正后的P值(padj)<0.05)],包括82个下调DEGs和41个上调DEGs。GO聚类分析发现下调DEGs主要富集在S-腺苷甲硫氨酸依赖性甲基转移酶活性、酸性磷酸酶活性、氧化还原酶活性、蛋白质异二聚化活性等,上调DEGs主要富集在跨膜受体活性、分子传感器活性、神经学系统过程等功能条目。KEGG富集分析显示下调DEGs主要富集在核糖体和溶酶体等细胞器、碳代谢和色氨酸代谢等代谢通路及Toll和Imd信号通路,上调DEGs未富集在KEGG代谢通路。荧光定量PCR结果显示测试的6个DEGs相对表达量变化趋势与RNA-Seq测序结果FPKM(fragments per kilobase million)值变化趋势一致,证明RNA-Seq测序结果可靠。【结论】亚致死剂量吡虫啉胁迫意大利蜜蜂成年工蜂显著降低意大利蜜蜂的嗅觉学习能力,影响意大利蜜蜂脑部免疫解毒等相关基因的表达和酶活性及氧化还原等生物代谢过程;亚致死剂量吡虫啉短时胁迫会刺激意大利蜜蜂的嗅觉感觉过程及神经信号传导过程。