西方蜜蜂采集蜂上颚腺中高表达基因的筛选与表达分析

【目的】本研究旨在筛选西方蜜蜂Apis mellifera采集蜂上颚腺中高表达基因,为进一步筛选和研究蜜蜂采集行为相关基因提供依据。【方法】基于前期测序的西方蜜蜂5种不同职能工蜂(3日龄工蜂、10日龄哺育蜂、10日龄采集蜂、21日龄哺育蜂和21日龄采集蜂)上颚腺转录组数据,筛选采集蜂上颚腺的差异表达基因(differentially expressed genes, DEGs),并对这些DEGs进行GO和KEGG分析;qRT-PCR检测随机选取的8个DEGs在10日龄哺育蜂和10日龄采集蜂上颚腺以及两个关键DEGs(Δ-1-吡咯啉-5-羧酸合成酶基因Amp5cs和细胞色素P450 9e2基因CYP9Q3)在工蜂不同发育时期和采集蜂各组织中的表达量。【结果】筛选到22个DEGs在21日龄采集蜂上颚腺中的表达量显著高于在3日龄工蜂、10日龄哺育蜂和21日龄哺育蜂上颚腺中的表达量,同时在10日龄采集蜂上颚腺中的表达量也显著高于在10日龄哺育蜂上颚腺中的表达量。GO和KEGG富集分析显示这些DEGs主要富集在胆固醇代谢、半乳糖代谢、淀粉和蔗糖代谢、精氨酸和脯氨酸代谢、细胞凋亡-果蝇、氨基酸生物合成等方面。qRT-PCR结果表明,8个DEGs(LOC100576395, LOC411983, LOC410235, LOC725581, LOC410527, LOC406131, LOC408453和LOC410253)的表达模式与转录组数据的表达模式一致;2个关键DEGs Amp5cs和CYP9Q3在工蜂各发育阶段均有表达,且在采集蜂中表达量最高;Amp5cs在采集蜂腹、胸、上颚腺和触角中高量表达,P450 9e2在采集蜂触角 和足中表达量显著高于在其他组织中的。【结论】本研究在减小日龄因素干扰下筛选了西方蜜蜂采集蜂上颚腺中22个高表达的DEGs,这些DEGs可能主要参与采集蜂上颚腺生理发育以及能量供应、外源性物质解毒、花蜜转化等代谢通路,进而影响蜜蜂的采集行为。这些结果为西方蜜蜂上颚腺的功能研究提供理论参考,同时也为采集力强的新品种培育奠定了基础。     

亚致死剂量吡虫啉胁迫对意大利蜜蜂工蜂嗅觉学习行为的影响及其脑部转录组分析

【目的】分析吡虫啉处理对意大利蜜蜂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测序结果可靠。【结论】亚致死剂量吡虫啉胁迫意大利蜜蜂成年工蜂显著降低意大利蜜蜂的嗅觉学习能力,影响意大利蜜蜂脑部免疫解毒等相关基因的表达和酶活性及氧化还原等生物代谢过程;亚致死剂量吡虫啉短时胁迫会刺激意大利蜜蜂的嗅觉感觉过程及神经信号传导过程。     

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

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

越冬西方蜜蜂4个亚种的转录组学差异分析

【目的】为进一步探讨西方蜜蜂Apis mellifera在越冬期基因表达模式的变化规律及抗寒机理,进而更好地指导蜂群越冬管理。【方法】利用转录组学鉴定和比较了西方蜜蜂4个亚种高加索蜂A.m.caucasica、欧洲黑蜂A.m.mellifera、意大利蜜蜂A.m.ligustica和卡尼鄂拉蜂A.m.carnica在越冬期的差异表达基因(differentially expressed genes, DEGs),并进行了GO和KEGG富集分析。利用qRT-PCR验证转录组测序结果的可靠性。【结果】西方蜜蜂4个亚种在越冬期DEGs的表达量变化各具特点,其中意大利蜜蜂在越冬初期到越冬中期DEGs表达量变化活跃,卡尼鄂拉蜂和欧洲黑蜂在越冬中期到越冬末期DEGs表达量变化活跃,高加索蜂在整个越冬期DEGs表达量相对稳定。越冬初期西方蜜蜂4个亚种的DEGs表达模式存在差异,高加索蜂和欧洲黑蜂中表皮蛋白22、卵黄原蛋白、气味结合蛋白14、热休克蛋白70和细胞色素P4509e2基因的表达量比意大利蜜蜂和卡尼鄂拉蜂的高,意大利蜜蜂中几丁质酶、防御素1、视黄醇脱氢酶和细胞色素b5基因的表达量较高,而卡尼...     

胁迫意大利蜜蜂幼虫肠道的球囊菌的转录组分析

【目的】本研究旨在通过趋势分析对胁迫意大利蜜蜂Apis mellifera ligustica(简称意蜂)幼虫肠道的球囊菌Ascosphaera apis的差异表达基因(DEGs)进行转录组分析。【方法】将纯化的球囊菌孢子配制为1×107孢子/m L的饲料饲喂意蜂3日龄幼虫,利用Illumina Hi Seq 2500平台对球囊菌胁迫的意蜂幼虫肠道c DNA进行测序,将过滤得到的有效读段(clean reads)映射(mapping)到核糖体数据库及意蜂参考基因组,最后将未映射上的reads映射到本课题组组装并注释的球囊菌参考转录组。利用STEM软件对DEGs进行趋势分析。利用WEGO软件对显著表达模式中的DEGs进行GO富集分析。利用Blastall对显著表达模式中的DEGs进行KEGG代谢通路富集分析。最后,通过对随机选取的6个DEGs进行RT-q PCR分析,对RNA-seq数据进行验证。【结果】球囊菌胁迫意蜂幼虫肠道样品的Illumina测序共得到球囊菌的25 454 076条原始读段(raw reads),经过滤得到24 909 820条clean reads,Q30均在93.46%以上。趋势分析结果显示,19 893个DEGs聚类为8个表达模式,其中有12 151个DEGs聚类为3个表现为显著上调趋势的表达模式。GO富集分析结果显示,表现上调趋势的DEGs富集在40个GO term,富集基因数最多的为细胞进程(cellular process)(2 601 unigenes),其次为代谢进程(metabolic process)(2 553 unigenes)和细胞(cell)(2 522unigenes)。KEGG代谢通路富集分析结果显示,上调趋势中的DEGs富集在119个代谢通路上,其中富集基因数最多的是核糖体(ribosome)(213 unigenes),其次为氨基酸生物合成(biosynthesis of amino acids)(154 unigenes)和内质网蛋白加工(protein processing in endoplasmic reticulum)(130unigenes)。共有48个DEGs富集在MAPK信号通路上,聚类分析结果显示,这些DEGs随着胁迫时间的延长表达水平逐渐增强。RT-q PCR结果显示,6个DEGs的表达水平变化趋势与RNA-seq数据一致,证明了本研究中的转录组数据真实可靠。【结论】本研究为在分子水平揭示球囊菌的致病机理提供了重要信息,也为阐释球囊菌胁迫意蜂幼虫过程中的病原-宿主互作机制奠定了基础。     

意大利蜜蜂幼虫肠道响应球囊菌早期胁迫的转录组学分析

【目的】球囊菌特异性地侵染蜜蜂幼虫而导致白垩病,该病造成成年蜜蜂数量的大幅下降,从而严重影响蜂蜜等产品的产量。目前,尚无利用二代测序技术研究白垩病的报道。本研究利用RNA-seq技术对正常及球囊菌胁迫后的意大利蜜蜂4日龄幼虫肠道进行深度测序,为解析宿主响应球囊菌胁迫的应答机制提供了重要的参考信息。【方法】利用Illumina Hiseq2500平台对对照组(AmCK)和处理组(AmT)幼虫肠道进行双端(PE125)测序;利用Perl脚本进行下机数据过滤;利用R软件进行测序饱和度分析、计算各样品之间的相关性系数;利用SOAP aligner/soap2软件将未比对上核糖体的读段(Reads)比对到意大利蜜蜂参考基因组;基于GO数据库和KEGG数据库进行GO分类和KEGG代谢通路(Pathway)富集分析。【结果】RNA-seq共得到181 096 194条原始读段(Raw reads),经过滤得到179078764条有效读段(Clean reads)。差异表达基因(DEGs)分析结果显示上调与下调基因的数量分别为4和20个。DEG的基因本体论(Gene ontology,GO)分析结果显示,这些DEGs分布于10个GO分类,除归入结合(Binding)的基因和催化活性(Catalytic activity)的部分基因外,绝大多数GO分类上的基因均表现为下调。DEGs的KEGG代谢通路分析结果显示各有1个DEG分别富集在核糖体和氧化磷酸化且均下调表达。【结论】本研究揭示了意大利蜜蜂幼虫肠道在球囊菌入侵早期的胁迫应答,为解析宿主响应球囊菌胁迫的应答机制奠定了基础。     

不同行为表现的蜜蜂脑部多巴胺及其受体基因的检测分析

【目的】对不同行为表现的意大利蜜蜂Apis mellifera ligustica工蜂脑部多巴胺及其受体基因进行检测。【方法】从蜜蜂观察箱中采集不同行为表现的工蜂,解剖其脑部,用高效液相色谱-电化学检测器和实时荧光PCR仪分别检测多巴胺含量和受体基因相对表达量,并进行计算分析研究。【结果】建立了不同行为表现的蜜蜂脑部多巴胺及受体基因的研究方法。意大利蜜蜂哺育蜂脑部的多巴胺含量与新出房工蜂、采集蜂、休息状态下的工蜂存在显著差异。多巴胺受体1基因及多巴胺转运体基因在哺育蜂脑部处于高表达状态,4种不同分工的意大利蜜蜂脑部多巴胺受体2基因、受体3基因相对表达量差异不显著。【结论】多巴胺与蜜蜂不同行为表现密切相关,多巴胺受体1基因与转运体基因可能参与蜜蜂哺育行为。     

亚致死剂量吡虫啉胁迫对意大利蜜蜂内勤蜂与外勤蜂免疫解毒相关基因表达及免疫解毒酶系活力的影响

摘要: 【目的】吡虫啉(imidaclorprid)是广泛使用的新烟碱类杀虫剂之一。大量研究表明亚致死剂量吡虫啉影响意大利蜜蜂Apis mellifera ligustica(简称“意蜂”)幼虫的发育和成年蜜蜂的采集、学习等行为。本实验旨在探究亚致死剂量吡虫啉对意大利蜜蜂内勤蜂(1日龄成年工蜂)与外勤蜂(21日龄成年工蜂)免疫解毒相关基因表达及免疫解毒酶系活力的影响,进而为蜜蜂健康的维护提供科学依据。【方法】测定饲喂含0.1 ng/μL吡虫啉的50%蔗糖溶液不同时间后意蜂成年工蜂的存活率;利用荧光定量PCR检测饲喂含0.1 ng/μL吡虫啉的50%蔗糖溶液6 d后其体内免疫基因多酚氧化酶基因(PPOA3, GenBank登录号： GB43738), Abaecin类抗菌肽基因(ABA, GenBank登录号： GB18323),葡萄糖脱氢酶基因(GLD, GenBank登录号： GB43007)和解毒基因细胞色素P450基因(CYP450 6a2, GenBank登录号： GB49876)的表达,并采用双抗体一步夹心法酶联免疫吸附试验测定其体内细胞色素P450酶(cytochrome P450, CYP450)含量和多酚氧化酶(polyphenol oxidase,PPO)活力。【结果】1日龄和21日龄意蜂成年工蜂连续饲喂6 d含0.1 ng/μL吡虫啉的50%蔗糖溶液后,其存活率与对照组(饲喂含0.1 ng/μL丙酮的50%蔗糖溶液)无显著差异;连续饲喂9 d含0.1 ng/μL吡虫啉的50%蔗糖溶液后,1日龄意蜂成年工蜂存活率与对照组无显著差异,而21日龄意蜂成年工蜂存活率与对照组有显著差异。1日龄意蜂成年工蜂自由取食含0.1 ng/μL吡虫啉的蔗糖溶液6 d后, PPOA3, CYP450 6a2, ABA和GLD表达水平,细胞色素P450含量以及多酚氧化酶活力与对照组相比均有显著下调趋势;而21日龄意蜂成年工蜂取食该药液6 d后,CYP450 6a2, ABA和GLD表达水平及多酚氧化酶活力与对照组相比均有显著下调趋势,PPOA3表达水平和细胞色素P450含量有显著上调趋势。【结论】亚致死剂量吡虫啉影响意大利蜜蜂内勤蜂与外勤蜂免疫解毒相关基因的表达及免疫解毒酶系活力;吡虫啉短期胁迫对意大利蜜蜂内勤蜂与外勤蜂的存活无显著影响,长期胁迫则会影响意大利蜜蜂内勤蜂与外勤蜂的存活。     

微小RNA ame-miR-31a和ame-miR-13b在意大利蜜蜂哺育蜂和采集蜂脑部的表达差异

【目的】研究微小RNA ame-miR-31a和ame-miR-13b在意大利蜜蜂Apis mellifera ligustica哺育蜂和采集蜂脑部的表达情况,探索其在蜜蜂哺育蜂和采集蜂行为转变中的调节作用。【方法】利用qPCR分别检测ame-miR-31a和ame-miR-13b在意大利蜜蜂工蜂正常蜂群和同龄蜂群中哺育蜂和采集蜂脑部的表达量。【结果】不论在意大利蜜蜂工蜂正常蜂群中还是在同龄蜂群中,ame-miR-31a在哺育蜂脑部的表达量总是显著高于其在采集蜂脑部的表达量;在正常蜂群中,ame-miR-13b在采集蜂脑部的表达量显著高于其在哺育蜂脑部中的表达量;在同龄蜂群中,ame-miR-13b在年老哺育蜂脑部的表达量极显著高于其在同龄采集蜂脑部的表达量,在年轻哺育蜂脑部与同龄年轻采集蜂脑部的表达量差异不显著。【结论】ame-miR-13b在意大利蜜蜂哺育蜂和采集蜂脑部的表达没有明显规律,而ame-miR-31a总是在哺育蜂脑部高表达,提示ame-miR-31a在蜜蜂哺育蜂和采集蜂行为转变中可能发挥重要的调控作用。     

RNAi沉默王浆主蛋白1(Mrjp1)基因对意大利蜜蜂工蜂学习和记忆的影响

【目的】前期研究发现经吡虫啉处理的意大利蜜蜂Apis mellifera ligustica(简称“意蜂”)工蜂学习能力下降,转录组学分析表明王浆主蛋白1(major royal jelly protein 1, MRJP1)基因在吡虫啉处理的蜜蜂脑中显著下调,MRJP1可能参与调控蜜蜂学习能力。本研究旨在采用RNA干扰(RNAinference, RNAi)技术将Mrjp1特异性沉默,验证MRJP1在意蜂工蜂嗅觉学习中的关键作用。【方法】通过克隆技术获得Mrjp1基因cDNA 序列,经测序验证后,设计引物,合成用于RNAi干扰Mrjp1基因表达的dsRNA。注射dsMrjp1的意蜂工蜂作为处理组(dsMrjp1注射组),注射dsEGFP的意蜂工蜂作为对照组(dsEGFP注射组),随后通过伸吻反应(proboscis extension response, PER)实验比较两组的嗅觉学习与记忆能力差异。最后采用实时荧光定量PCR(quantitative real-time PCR, qRT-PCR)检测注射dsMrjp1后意大利蜜蜂工蜂脑中Mrjp1的相对表达量。【结果】dsMrjp1注射组与dsEGFP注射组意蜂工蜂学习能力差异显著,dsMrjp1注射组意蜂工蜂的学习能力显著降低。学习后2 h,两组意蜂工蜂的记忆力无显著差异。qRT-PCR结果显示Mrjp1的表达水平在dsMrjp1注射组意蜂工蜂脑中显著低于dsEGFP注射组,表明学习能力降低的处理组意蜂脑内对应的Mrjp1表达水平也降低。【结论】通过RNAi抑制意蜂工蜂Mrjp1基因的表达后,其嗅觉学习能力受到显著性抑制,但记忆力未受到显著影响,提示Mrjp1可能是调控意蜂学习的重要基因之一。本研究结果有助于后续进一步研究蜜蜂嗅觉学习相关的分子机制。     

意大利蜜蜂对两种油菜花朵气味的偏爱性

【目的】明确意大利蜜蜂Apis mellifera ligustica（简称“意蜂”）对两种油菜（甘蓝型胜利油菜 Brassica napus cv. Shengli和芥菜型马尾油菜 B. juncea cv. Mawei）的气味偏爱性及选择行为特征。【方法】本研究开展大田访花偏爱性试验、室外群体气味偏爱性试验以及室内“Y”型嗅觉仪行为反应和学习记忆试验,测定意蜂蜂群及青年采集蜂对甘蓝型胜利油菜和芥菜型马尾油菜两种油菜花朵的选择次数,并通过采集蜂训练时及训练后的喙伸反应率分析两种油菜花朵气味对采集蜂学习记忆能力的影响,最终评估意蜂蜂群及个体对这两种油茶花朵气味的选择性。【结果】大田试验12:00-13:00时间段,访问胜利油菜的意蜂数量显著高于访问马尾油菜的数量(P<0.05）。群体气味偏爱性试验14:00-15:00时间段意蜂访问胜利油菜的数量显著高于访问马尾油菜的数量（P<0.05）。行为反应试验结果显示意蜂选择胜利油菜的次数（3.86±2.83）显著高于访问马尾油菜的次数（2.28±1.87）（P<0.05）。在学习记忆试验中,随着训练次数的增加,意蜂对两种油菜花朵气味的喙伸反应率逐渐升高,相同的训练次数,意蜂对两种油菜气味的学习能力差异不显著（P>0.05）;训练结束后,随着时间的增加,意蜂对胜利油菜及马尾油菜的气味记忆逐渐下降,在24 h对胜利油菜的气味记忆明显高于对马尾油菜的气味记忆（P<0.05）。【结论】本研究证明意蜂较偏爱胜利油菜的花朵气味,花朵气味是影响其采集偏爱的重要因素,气味响应性试验可定量分析蜜蜂对花朵气味的嗅觉敏感性。     

中华蜜蜂气味结合蛋白AcerOBP14的表达及配体结合特性分析

[目的]气味结合蛋白(odorant binding proteins,OBPs)在昆虫寄主定位、产卵地选择等行为中发挥着重要作用,明确中华蜜蜂Apis cerana cerana AcerOBP14与配体的结合特性有助于阐明中华蜜蜂嗅觉识别的分子机制.[方法]通过qRT-PCR测定OBP14在20日龄中华蜜蜂成年工蜂...     

转录组分析揭示东方蜜蜂微孢子虫侵染意大利蜜蜂的分子机制

【目的】本研究旨在通过差异表达基因(differentially expressed gene, DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log_2(Fold change)|≥1的标准,通过比较分析筛选出NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组的DEG。通过相关生物信息学软件对上述DEG进行Venn分析、GO分类和KEGG代谢通路富集分析。根据Nr和KEGG数据库注释信息和相关文献进行对东方蜜蜂微孢子虫的毒力因子和侵染相关因子的统计和分析。通过RT-qPCR验证转录组数据及DEG表达趋势。【结果】从NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组分别鉴定出1 397, 1 497和52个DEG。Venn分析结果显示各比较组共有的上调和下调基因分别为10和1个。GO分类结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG富集数最多的功能条目为代谢进程、细胞进程、单组织进程、细胞、细胞组件、细胞器、催化活性和结合,而NcT1 vs NcT2中DEG富集数最多的是代谢进程、细胞进程、单组织进程、催化活性和结合。KEGG代谢通路富集分析结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG分别富集到80和79条通路;富集在糖酵解/糖异生和MAPK信号通路的上调基因数量多于下调基因。毒力因子分析结果显示,孢壁蛋白9基因和孢壁蛋白12基因在NcCK vs NcT1和NcCK vs NcT2中均下调表达,孢壁蛋白8基因仅在NcCK vs NcT1中表达量下调;此外孢壁蛋白前体基因、孢壁和锚定盘复合蛋白基因、几丁质合酶基因、极管蛋白基因、蓖麻毒素B凝集素基因的表达水平在NcCK vs NcT1和NcCK vs NcT2中表现为上调。侵染相关因子分析结果表明,糖酵解途径的3个关键酶基因在NcCK vs NcT1和NcCK vs NcT2中上调表达;3个涉及ATP/ADP移位酶的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,但有1个表达量下调;2个涉及ABC转运蛋白的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,另有4个下调表达。RT-qPCR结果证实了本研究中转录组数据及DEG表达趋势的真实可靠性。【结论】本研究通过比较分析解析东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程的转录组动态,揭示了孢壁蛋白、孢壁和锚定盘复合蛋白、几丁质酶、极管蛋白和蓖麻毒素B凝集素等毒力因子编码基因,以及己糖激酶、丙酮酸激酶、6-磷酸果糖激酶、ATP/ADP移位酶和ABC转运蛋白等侵染相关因子编码基因在病原增殖中扮演重要角色,为阐明东方蜜蜂微孢子虫的侵染机制提供了基础。     

低温胁迫下意大利蜜蜂预蛹的差异表达基因趋势分析

【目的】蜜蜂是典型的具有发育狭温性的全变态昆虫。本研究以对低温最敏感的意大利蜜蜂Apis mellifera ligustica预蛹为研究对象,通过低温胁迫不同时间的差异表达基因(differentially expressed genes, DEGs)趋势分析,探讨低温胁迫对蜜蜂发育影响的关键基因。【方法】对3日龄意大利蜜蜂封盖子预蛹进行20℃低温胁迫18 h(T18)和36 h(T36),以未经低温胁迫的预蛹为对照(CK),通过Illumina HiSeq~(TM)平台进行转录组学测定。利用Short Time-series Expression Miner(STEM)软件对2个处理组与对照组比较共有的差异表达基因进行趋势分析,再进一步对显著富集趋势模式中富集的差异表达基因进行GO分类和KEGG pathway分析。利用RT-qPCR对随机挑选的5个DEGs的表达模式进行验证。【结果】对检测到1 062个T18 vs CK和T36 vs CK共有的DEGs进行趋势分析,发现3个显著的基因表达模式,包括2个上调表达模式(Profile 6,有539个基因;Profile 7,有271个基因),1个下调表达模式(Profile 1,有183个基因)。对3个显著富集趋势模式DEGs分别进行GO富集分析和KEGG pathway分析,在Profile 6中找到同时富集在FoxO信号通路和寿命调节信号通路上的胰岛素样肽基因ILP和叉头蛋白O基因FoxO持续上调,说明低温胁迫会影响蜜蜂预蛹蜕皮激素信号传递;在Profile 7中CREB结合蛋白基因CBP持续上调,蜜蜂预蛹受到低温胁迫会影响细胞发育;在Profile 1中细胞色素P450基因CYP450 306a1 (phm)和WNT1显著下调,说明化蛹时受到低温胁迫会影响蜕皮激素的合成。通过RT-qPCR分析挑选的5个基因的表达结果和高通量测序结果一致。【结论】通过对蜜蜂响应低温胁迫的差异表达基因趋势分析发现,蜜蜂胰岛素和蜕皮激素共同调控的FoxO可能是低温胁迫抑制蜜蜂化蛹的一个关键基因。本研究为探索低温胁迫影响蜜蜂发育变态的分子调控机制奠定了基础。     

ame-miR-79负调控意大利蜜蜂工蜂幼虫肠道中靶基因CYP450和FG的表达

【目的】本研究旨在通过饲喂ame-miR-79的模拟物(mimic)和抑制物(inhibitor)对意大利蜜蜂Apis mellifera ligustica工蜂幼虫肠道内的ame-miR-79分别进行过表达和敲减,探究ame-miR-79对幼虫肠道内靶基因表达的调控作用。【方法】通过分子克隆与Sanger测序验证意大利蜜蜂工蜂幼虫肠道内ame-miR-79的序列真实性。通过饲喂mimic-miR-79和inhibitor-miR-79对意大利蜜蜂工蜂4-6日龄幼虫肠道内的ame-miR-79分别进行过表达和敲减。采用相关生物信息学软件进行ame-miR-79的靶基因预测与分析。通过RT-qPCR检测ame-miR-79的过表达和敲减效果及过表达和敲减ame-miR-79后靶基因的相对表达量。【结果】ame-miR-79在意大利蜜蜂工蜂幼虫肠道内真实存在。与无义模拟物(nonsense mimic, mimic-NC)组相比,mimic-miR-79组的4-6日龄幼虫肠道内ame-miR-79的表达量皆极显著上调;与无义抑制物(nonsense inhibitor, inhibitor-NC)组相比,inhibitor-miR-79组的4和5日龄幼虫肠道内ame-miR-79的表达量显著下调,6日龄幼虫肠道内ame-miR-79的表达量极显著下调。ame-miR-79共靶向303个基因,涉及27个GO条目和179条KEGG通路。相较于mimic-NC组,靶基因细胞色素P450基因CYP450在mimic-miR-79组的4和5日龄幼虫肠道内均极显著下调表达,而在6日龄幼虫肠道内上调表达;靶基因fringe糖基化转移酶(fringe glycosyltransferase, FG)基因在4日龄幼虫肠道内下调表达,在5日龄幼虫肠道内显著下调表达,而在6日龄幼虫肠道内上调表达。相较于inhibitor-NC组,CYP450在inhibitor-miR-79组的4日龄幼虫肠道内极显著上调表达,在6日龄幼虫肠道内上调表达,而在5日龄幼虫肠道内下调表达;FG的表达水平在4日龄幼虫肠道内显著上调,在5和6日龄幼虫肠道内极显著上调。【结论】ame-miR-79在意大利蜜蜂工蜂幼虫肠道内真实存在;通过饲喂模拟物和抑制物能分别实现意大利蜜蜂工蜂幼虫肠道内ame-miR-79的有效过表达和敲减;ame-miR-79负调控意大利蜜蜂工蜂幼虫肠道内CYP450和FG的表达。     

ame-miR-14在意大利蜜蜂工蜂幼虫肠道发育过程的调控作用

【目的】本研究旨在揭示ame-miR-14在意大利蜜蜂Apis mellifera ligustica工蜂幼虫肠道发育过程的调控作用。【方法】通过Stem-loop RT-PCR和Sanger测序分别验证ame-miR-14在意大利蜜蜂工蜂6日龄幼虫肠道中的表达和序列真实性。饲喂ame-miR-14的模拟物(mimic-ame-miR-14)和抑制物(inhibitor-ame-miR-14)及其相应的阴性对照mimic-NC和inhibitor-NC对ame-miR-14分别进行过表达和敲降,利用RT-qPCR检测意大利蜜蜂工蜂4-6日龄幼虫肠道中ame-miR-14的表达量。利用生物信息学软件预测ame-miR-14的靶基因并进行相关分析。利用RT-qPCR检测ame-miR-14的过表达和敲降后其靶基因FoxO和Hedgehog在4-6日龄幼虫肠道中的相对表达量。【结果】ame-miR-14在意大利蜜蜂工蜂6日龄幼虫肠道中真实存在和表达。相较于饲喂mimic-NC,饲喂mimic-ame-miR-14后ame-miR-14表达量在意大利蜜蜂工蜂4-6日龄幼虫肠道中均为显著上调;相...