意大利蜜蜂工蜂中肠piRNA的鉴定与分析

【目的】Piwi蛋白互作RNA(PIWI-interacting RNA, piRNA)在昆虫的发育和免疫等重要生物学过程发挥重要的调控作用。本研究旨在丰富西方蜜蜂Apis mellifera的piRNA信息,并为进一步探究差异表达piRNA(differentially expressed piRNA, DEpiRNA)调控意大利蜜蜂Apis mellifera ligustica工蜂中肠发育的分子机理提供基础。【方法】基于已获得的意大利蜜蜂7日龄(Am7)和10日龄(Am10)工蜂中肠的small RNA(sRNA)组学数据,对piRNA进行预测和分析。将质控后的sRNA数据比对西方蜜蜂参考基因组,再将比对上的序列标签(tags)进一步比对数据库以滤除rRNA和tRNA等小非编码RNA(non-coding RNA,ncRNA),进而根据piRNA的长度特征鉴定piRNA。采用TPM(tags per million)算法对piRNA表达量进行计算和归一化处理。根据|log₂ fold change|≥1且P≤0.05的标准筛选Am7 vs Am10 比较组的DEpiRNA。通过相关软件预测DEpiRNA的靶mRNA并进行GO和KEGG数据库注释。利用Cytoscape软件对DEpiRNA-mRNA调控网络进行可视化。通过Stem-loop RT-PCR对随机选取的Am7与Am10两组共有的6个piRNA的表达进行验证。利用RT-qPCR对随机挑选的6个DEpiRNA的表达趋势进行验证。【结果】从意大利蜜蜂工蜂中肠中共鉴定到596个piRNA,长度介于24~33 nt,且Am7和Am10组中不同长度分布范围的piRNA的首位碱基偏向性具有明显差异。在Am7 vs Am10比较组共筛选出41个DEpiRNA,其中piR-ame-11093, piR-ame-1111451,piR-ame-190949和piR-ame-932156可分别靶向1 195, 1 018, 4 040和1 063条mRNA。上述靶mRNA可分别注释到45个功能条目和45条通路。Stem-loop RT-PCR验证结果显示6个piRNA(piR-ame-1084826, piR-ame-11093, piR-ame-14476, piR-ame-24995, piR-ame-39500和piR-ame-774987)均真实表达。RT-qPCR结果显示共有6个DEpiRNA (piR-ame-1084826,piR-ame-11093, piR-ame-14476, piR-ame-24995, piR-ame-39500和piR-ame-774987)的表达趋势与测序数据中的表达趋势一致,证实了本研究中sRNA-seq数据的可靠性和piRNA差异表达趋势的真实性。【结论】本研究从意大利蜜蜂工蜂中肠中鉴定到596个piRN,长度介于24~33 nt。不同长度分布范围的piRNA的首位碱基偏向性具有明显差异;piR-ame-11093, piR-ame-1111451,piR-ame-190949和piR-ame-932156具有靶向调控基因表达参与工蜂中肠发育过程的潜力。     

意大利蜜蜂工蜂唾腺和中肠消化酶活性比较研究

对意大利蜜蜂 Apis mellifera ligustica L工蜂唾腺和中肠中淀粉酶、蛋白酶、蔗糖酶、果胶酶和海藻糖酶5种消化酶进行活性测定,发现巾肠的蛋门酶和蔗糖酶活性较唾腺高,有显著性差异;唾腺的果胶酶活性较中肠高,有显著性差异;唾腺与中肠的淀粉酶和海藻糖酶没有显著性差异.表明意大利蜜蜂工蜂对蔗糖和蛋白的消化主要由中肠分泌的消化酶来完成,唾腺起辅助作用;对果胶的消化主要在唾腺,中肠起辅助作用;而在淀粉和海藻糖的消化过程中唾腺和中肠都有很重要的作用.     

参与调控意大利蜜蜂工蜂中肠基因表达的东方蜜蜂微孢子虫miRNA的组学解析及其调控网络

[目的]本研究结合前期已获得的miRNA和mRNA组学数据对东方蜜蜂微孢子虫Nosema ceranae 的差异表达 miRNA(differentially expressed miRNA,DEmiRNA)靶向意大利蜜蜂 Apis mellifera ligustica 工蜂中肠的 mRNA 和差异表达 mRNA(d...     

枣花病对意大利蜜蜂消化酶和解毒酶系活性的影响

[目的]枣花期蜜蜂易发生枣花中毒(又称枣花病),致使大量采集蜂死亡.本研究初步探讨了枣花病对蜜蜂的毒性作用机理.[方法]以正常的意大利蜜蜂Apis mellifera ligustica工蜂为对照,采用分光光度法测定山西临猗、太谷和临县患枣花病的意大利蜜蜂工蜂体内淀粉酶、蛋白酶、蔗糖酶、果胶酶、海藻糖酶、酯酶(este...     

意大利蜜蜂工蜂中肠的环状RNA及其调控网络分析

【目的】环状RNA(circRNA)在可变剪接、转录调控和来源基因的表达调控等方面具有重要功能。本研究旨在分析意大利蜜蜂Apis mellifera ligustica工蜂中肠circRNA的数量、种类、结构特征和作用,并通过构建和分析circRNA的调控网络探索circRNA的调控功能。【方法】在实验室条件下人工饲养意大利蜜蜂工蜂,利用circRNA-seq技术对意大利蜜蜂7和10日龄成年工蜂中肠样品进行深度测序。利用find_circ软件从质控后的数据中预测circRNA。通过BLAST比对GO和KEGG数据库,对circRNA的来源基因进行功能和代谢通路注释。利用TargetFinder软件预测circRNA靶向结合的miRNA及miRNA靶向结合的mRNA,通过Cytoscape v.3.2.1软件对circRNA-miRNA和circRNA-miRNA-mRNA调控网络进行构建及可视化。通过设计背靠背引物和线性扩增引物RT-PCR对预测出的circRNA进行验证。【结果】意大利蜜蜂工蜂中肠样品的测序平均得到136 463 071条clean reads,去除rRNA后各样品的anchor reads均在136 779 122条及以上。共预测出10 833个circRNA,长度主要介于15~1 000 nt;上述circRNA的类型丰富,其中已注释的外显子circRNA数量最多,分布在西方蜜蜂1号染色体的circRNA数量最多,其次为8号染色体。CircRNA的来源基因可注释到包括结合、细胞进程和细胞在内的45个GO条目,以及包括内吞作用、内质网蛋白加工及核糖体在内的121条KEGG代谢通路,表明circRNA在意大利蜜蜂工蜂中肠的生长、发育、新陈代谢和细胞生命活动等生物学过程中发挥重要作用。进一步构建circRNA-miRNA和circRNA-miRNA-mRNA调控网络,分析结果显示部分circRNA可能作为竞争性内源RNA吸附结合microRNA,从而调控基因的表达水平。最后,对随机选择的3个circRNA的RT-PCR结果验证了其真实存在。【结论】本研究对意大利蜜蜂工蜂中肠中的circRNA进行预测、分析及鉴定。研究结果提供了中肠circRNA的数量、种类、结构特征、作用和调控网络的信息,揭示了circRNA可能通过作用于来源基因和作为竞争性内源RNA在意大利蜜蜂工蜂中肠的生长发育和免疫防御中发挥作用,为深入研究circRNA在意大利蜜蜂中肠发育及胁迫响应过程中的功能奠定了基础。     

意大利蜜蜂工蜂中肠的长链非编码RNA的预测、分析及鉴定

【目的】预测、分析和鉴定意大利蜜蜂Apismelliferaligustica工蜂中肠的长链非编码RNA(lncRNA),探究lncRNA的作用。【方法】结合lncRNA-seq技术和链特异性cDNA建库方法对意大利蜜蜂工蜂中肠进行深度测序;利用Perl脚本对原始数据进行过滤;联用CPC和CNCI软件对测序数据进行lncRNA预测,并比较lncRNA基因与其邻近的蛋白编码基因的结构特征;通过RT-PCR对随机选取的lncRNA进行鉴定;利用相关生物信息学软件对lncRNA、的上下游基因进行GO分类和KEGG代谢通路富集分析。【结果】意大利蜜蜂工蜂中肠样品的lncRNA-seq共得到1 956 118 858原始读段,经过滤得到1 946 489 304有效读段,共预测出6 353个lncRNA,它们较蛋白编码基因含有较少的外显子数、较短的转录本长度。通过RT-PCR验证了9个lncRNA的表达。lncRNA的上下游基因富集在42个GO条目和256条代谢通路,进一步分析结果表明这些lncRNA参与调控意大利蜜蜂工蜂中肠的新陈代谢和细胞生命活动等生物学过程。【结论】研究结果丰富了蜜蜂的lncRNA信息,也为阐明lncRNA在意大利蜜蜂工蜂中肠发育及胁迫响应过程中的作用打下了基础。     

意大利蜜蜂工蜂卵巢发育的研究

蜜蜂Apis mellifera L.蜂群中的工蜂卵巢发育和工蜂产卵现象受多种因素控制,了解其影响因素对养蜂生产具有重要意义。本研究将意大利蜜蜂蜂群设置为囚王群(tg1)、无王有子群(tg2)、无王无子群(tg3)以及正常有王蜂群(CK)4个试验组,通过对工蜂卵巢管的显微观察,确定不同处理组工蜂在不同时间段内卵巢的发育情况。结果表明:随着时间的延长,与CK处理组相比,tg3处理组中的工蜂卵巢发育水平最高,tg2次之,tg1最低;在31d时,4个处理组两两之间差异均达到显著水平(P<0.05)。tg1、tg2和tg3处理组中工蜂产卵前期时间分别为35、22和17d,而CK蜂群在试验期内未出现工蜂产卵现象;tg2和tg3处理组的工蜂产卵的封盖前期时间分别为8和6d,而tg1和CK组在试验期内未出现子房封盖现象。蜂群失王时间过长会刺激工蜂卵巢发育,并导致其产卵;蜂群的短期失王和蜂王老化也会刺激工蜂卵巢发育,但是刺激程度较低;蜂群中的蜂子能抑制工蜂卵巢管的发育,因此在蜂群短时间失王时可以适当地补充子脾延缓工蜂卵巢发育。     

意大利蜜蜂工蜂中肠SNP和InDel突变位点的挖掘及分析

本研究利用已获得的意大利蜜蜂Apis mellifera ligustica工蜂中肠的转录组数据对意蜂的单核苷酸多态性（Single nucleotide polymorphism, SNP）和插入缺失（Insertion-Deletion, InDel）突变位点进行挖掘和分析,共鉴定到232 678个SNP位点,其中发生转换和颠换的SNP位点数分别为196 087和36 591个;最丰富的突变类型是G/A,最少的突变类型为T/G;分布在内含子的SNP位点最多,其次为外显子和基因间区;密码子突变类型为同义突变的SNP位点数最多,其次是非同义突变、终止子增加和终止子减少;此外,SNP位点所在基因可注释到 50个 GO条目和351条KEGG通路。共鉴定到38 715个InDel位点,最丰富的突变类型为移码插入,其次是移码缺失;分布InDel位点数较多的基因组区域为内含子和基因间区;另外,InDel位点所在基因可注释到50个功能条目和340条通路。研究结果丰富了西方蜜蜂Apis mellifera的SNP和InDel位点信息,并为开发和利用意蜂的新型分子标记提供基础。     

狄斯瓦螨对意大利蜜蜂工蜂呼吸代谢的影响

【目的】意大利蜜蜂Apis mellifera ligustica(简称"意蜂")是世界重要的授粉昆虫,对维持生态平衡、提高农产品产量和保障粮食安全至关重要。但受狄斯瓦螨Varroa destructor侵害,全世界的意蜂蜂群每年都受到巨大损失。本研究主要探讨狄斯瓦螨对意蜂工蜂呼吸代谢的影响,为今后深入研究狄斯瓦螨与意蜂之间的相互作用奠定基础。【方法】应用SSI多通道昆虫呼吸仪测定感染狄斯瓦螨意蜂与健康意蜂不同发育阶段(4日龄幼虫、5日龄蛹、新出房成蜂和6日龄成蜂)工蜂的耗氧量、CO₂释放量和呼吸商,并分析其变化特征;同时采用实时荧光定量PCR(RT-qPCR)方法分析狄斯瓦螨感染意蜂与健康意蜂成年工蜂中呼吸代谢相关基因(Coq7,COXⅠ,CytB,CytC和IF-2mt)的表达量变化。【结果】受狄斯瓦螨侵染后,意蜂工蜂耗氧量、CO₂释放量呈现出先上升后下降的趋势,其中耗氧量从5日龄蛹期的0.0244 mL/min下降至6日龄成蜂期的0.0093 mL/min; CO₂释放量从5日龄蛹期的0.0174 mL/min下降...     

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

意大利蜜蜂胚后发育过程中中肠上皮组织细胞的更替

中肠是昆虫消化、 吸收营养物质的主要部位。本研究通过形态解剖、 BrdU免疫组织化学和原位末端转移酶标记(TUNEL)细胞凋亡检测等技术, 对意大利蜜蜂Apis mellifera ligustica中肠胚后发育过程中细胞的增殖和凋亡模式进行了比较研究。结果表明：意大利蜜蜂幼虫发育早期, 中肠的增加主要来自于上皮细胞的分裂以及再生细胞的增殖。在变态发育期间, 中肠上皮经历了广泛的重组, 由再生细胞重新形成的蛹上皮替代了幼虫上皮。再生细胞在蜜蜂中肠的整个发育阶段始终存在, 为中肠的生长和更替提供了主要的细胞来源。本研究为昆虫组织细胞自噬和凋亡机制的研究提供一定的证据。     

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

[目的]分析吡虫啉处理对意大利蜜蜂Apis mellifera ligustica嗅觉学习行为及脑部基因转录的影响,为新烟碱类杀虫剂对蜜蜂的负面影响提供依据.[方法]实验室条件下一次性饲喂意大利蜜蜂成年工蜂含有4 ng吡虫啉的50％蔗糖溶液,以饲喂不合吡虫啉的50％蔗糖溶液为对照,通过伸吻反应(proboscis ex...     

基于small RNA组学分析揭示意大利蜜蜂响应东方蜜蜂微孢子虫胁迫的免疫应答机制

【目的】通过深度测序和组学分析在small RNA组学层面揭示意大利蜜蜂(Apis mellifera ligustica,简称意蜂)响应东方蜜蜂微孢子虫(Nosema ceranae)的免疫应答机制。【方法】利用small RNA-seq技术对正常及N. ceranae胁迫7 d和10 d的意蜂工蜂中肠(Am7CK、Am7T、Am10CK和Am10T)进行深度测序。利用相关生物信息学软件对测序数据进行质控、已知micro RNA (mi RNA)鉴定、新mi RNA预测及mi RNA的结构特征分析。通过Stem-loop RT-PCR验证新mi RNA的表达。按照|log2(Fold change)|≥1和P≤0.05的标准筛选出Am7CKvs.Am7T和Am10CKvs.Am10T比较组的差异表达mi RNA(differentially expressed mi RNA,DEmi RNA)。利用软件预测DEmi RNA靶向结合的m RNA并进行GO和KEGG数据库注释,根据注释信息对细胞和体液免疫相关通路及富集靶m RNA进行统计和分析。根据靶向结合关系构建DEmi RNA和免疫通路相关差异表达m RNA (differentially expressed m RNA,DEm RNA)的调控网络。采用RT-q PCR对数据的可靠性和DEmi RNA的差异表达进行验证。【结果】共获得165895574条原始读段和132028990条有效序列标签,各组的组内Pearson相关性平均在87.92%及以上。共鉴定到928个已知mi RNA和56个新mi RNA。这些mi RNA的长度介于18–28nt,多数的长度为18nt和22nt且首位碱基主要偏向U。验证了12个新mi RNA的真实表达。Am7CKvs.Am7T比较组包含48个上调mi RNA和36个下调mi RNA;Am10CK vs. Am10T比较组包含56个上调mi RNA和51个下调mi RNA。两个比较组的DEmi RNA可分别靶向结合9827个和10720个m RNA。这些靶m RNA可分别注释到50和47条功能条目,以及138和135条KEGG通路。DEmi RNA与免疫通路相关靶m RNA的调控网络分析结果显示,Am7CK vs. Am7T中有26个DEmi RNA靶向与内吞作用等免疫通路相关的10个DEm RNA;Am10CK vs. Am10T中有15个DEmi RNA靶向与MAPK信号通路等免疫通路相关的10个DEm RNA。验证了测序数据和4个DEmi RNA差异表达的可靠性。【结论】研究结果揭示了宿主DEmi RNA可能通过调控物质和能量代谢、细胞和体液免疫对N. ceranae产生应答,但DEmi RNA不参与抗菌肽基因的表达调控;mi R-1-z可能参与宿主的细胞增殖、细胞凋亡和免疫进程;氧化磷酸化通路可能在宿主免疫应答及宿主-病原互作中发挥特殊作用。     

基于代谢组方法筛选蜜蜂副伤寒相关差异代谢物的研究

[目的]蜜蜂副伤寒是西方蜜蜂Apis mellifera越冬期的主要易感疾病.本研究旨在探究蜜蜂副伤寒疾病相关的代谢物和代谢通路.[方法]以越冬末期的西方蜜蜂为研究对象,利用基于高效液相质谱联用的非靶向代谢组学方法检测副伤寒患病群和正常越冬群工蜂肠道中代谢物的变化.[结果]蜜蜂副伤寒患病个体和正常蜜蜂个体在正、负离子模式下分别获得626个和518个差异代谢物,在正、负离子模式下最显著差异的10种代谢物中筛选到镰刀菌氧萘满酮、奎尼丁、褪黑素、去氧紫草素等与蜜蜂物质代谢障碍、细胞凋亡和抗氧化应激相关的差异代谢物.此外,差异代谢物显著富集于与氨基酸代谢密切相关的氨酰-tRNA合成、蛋白质消化吸收、谷胱甘肽代谢、矿物质代谢、多种氨基酸代谢等7个代谢通路.[结论]蜜蜂患副伤寒后肠道中代谢物水平发生显著变化,其中氨基酸和蛋白质的代谢异常可能是蜜蜂副伤寒发病的主要原因.     

东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程中nce-miR-10660及其靶基因表达谱

【目的】本研究旨在为探究nce-miR-10660调控东方蜜蜂微孢子虫Nosema ceranae侵染的作用机制提供理论和实验依据。【方法】采用Stem-loop RT-PCR对前期鉴定到的东方蜜蜂微孢子虫nce-miR-10660进行表达验证,再通过Sanger测序验证nce-miR-10660的序列。利用相关生物信息学软件预测和分析nce-miR-10660的靶基因。通过RT-qPCR检测nce-miR-10660及其靶基因在东方蜜蜂微孢子虫侵染意大利蜜蜂Apis mellifera ligustica工蜂过程中中肠中的表达谱。【结果】Stem-loop RT-PCR和Sanger测序结果分别证实了nce-miR-10660在东方蜜蜂微孢子虫孢子中的表达和真实存在。靶向预测结果显示nce-miR-10660共靶向RRDRP和RCDP 42等9个基因;分别有2和6个靶基因可被分别注释到KEGG数据库中的4条通路和GO数据库中的23个条目。RT-qPCR结果显示,相较于东方蜜蜂微孢子虫侵染后1 d时意大利蜜蜂工蜂中肠中nce-miR-10660的表达量,东方蜜蜂微孢子虫侵染后2 d时意...     

基于Illumina HiSeq测序技术研究吡虫啉对意大利蜜蜂雄蜂的影响

[目的]雄蜂对蜂群繁衍有着非常重要的作用.本研究旨在探究吡虫啉对意大利蜜蜂Apis mellifera Ligustica雄蜂生长发育和基因表达产生的影响.[方法]以意大利蜜蜂雄蜂为研究对象,分别以0.00001、0.0001和0.001 μg/μL浓度的吡虫啉对雄蜂幼虫进行连续饲喂处理.每天观察并记录幼虫的发育形态及死亡率,在雄蜂幼虫后期(移虫后6d)测量幼虫体重.利用Illumina HiSeq测序技术对经吡虫啉处理的雄蜂进行转录组测序,进而对差异表达基因进行深入分析.[结果]取食吡虫啉后的雄蜂幼虫,体重低于正常雄蜂,当浓度高于0.0001μg/μL时差异显著;雄蜂幼虫取食吡虫啉后出现死亡现象,且死亡率随吡虫啉浓度的升高而增大;差异表达基因分析结果上调与下调基因数量分别为390个和130个.GO富集分析结果上调基因共分布于55个GO条目,富集基因数量最多的是细胞进程、细胞、细胞组件、细胞膜、细胞膜组件、结合,下调基因共分布于48个GO条目,富集基因数量最多的是细胞进程、细胞、细胞组件.富集在有关生殖功能的差异表达基因中,上调基因数量为21个,下调基因数量为5个.KEGG代谢通路富集分析结果上调基因富集在159个通路上,其中富集基因数最多的是蛋白质消化吸收和神经活性配体-受体相互作用通路.下调基因富集在71个通路上,其中富集基因数最多的是溶酶体、胰液分泌、神经活性配体-受体相互作用通路.[结论]吡虫啉能抑制意大利蜜蜂雄蜂的生长发育,甚至造成幼虫死亡,同时,可以影响雄蜂的神经系统、代谢系统和生殖系统等.本研究结果为蜜蜂资源保护提供理论依据.     

基于Illumina HiSeq测序技术研究吡虫啉对意大利蜜蜂雄蜂的影响

[目的]雄蜂对蜂群繁衍有着非常重要的作用.本研究旨在探究吡虫啉对意大利蜜蜂Apis mellifera Ligustica雄蜂生长发育和基因表达产生的影响.[方法]以意大利蜜蜂雄蜂为研究对象,分别以0.00001、0.0001和0.001 μg/μL浓度的吡虫啉对雄蜂幼虫进行连续饲喂处理.每天观察并记录幼虫的发育形态及死亡率,在雄蜂幼虫后期(移虫后6d)测量幼虫体重.利用Illumina HiSeq测序技术对经吡虫啉处理的雄蜂进行转录组测序,进而对差异表达基因进行深入分析.[结果]取食吡虫啉后的雄蜂幼虫,体重低于正常雄蜂,当浓度高于0.0001μg/μL时差异显著;雄蜂幼虫取食吡虫啉后出现死亡现象,且死亡率随吡虫啉浓度的升高而增大;差异表达基因分析结果上调与下调基因数量分别为390个和130个.GO富集分析结果上调基因共分布于55个GO条目,富集基因数量最多的是细胞进程、细胞、细胞组件、细胞膜、细胞膜组件、结合,下调基因共分布于48个GO条目,富集基因数量最多的是细胞进程、细胞、细胞组件.富集在有关生殖功能的差异表达基因中,上调基因数量为21个,下调基因数量为5个.KEGG代谢通路富集分析结果上调基因富集在159个通路上,其中富集基因数最多的是蛋白质消化吸收和神经活性配体-受体相互作用通路.下调基因富集在71个通路上,其中富集基因数最多的是溶酶体、胰液分泌、神经活性配体-受体相互作用通路.[结论]吡虫啉能抑制意大利蜜蜂雄蜂的生长发育,甚至造成幼虫死亡,同时,可以影响雄蜂的神经系统、代谢系统和生殖系统等.本研究结果为蜜蜂资源保护提供理论依据.     

东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程的nce-miR-23928及其靶基因表达谱

【背景】东方蜜蜂微孢子虫(Nosema ceranae)专性侵染成年蜜蜂中肠上皮细胞而导致的微孢子虫病给养蜂业造成严重损失。【目的】检测东方蜜蜂微孢子虫nce-miR-23928及其靶基因在侵染意大利蜜蜂(Apis mellifera ligustica)工蜂过程的表达谱,为深入探究nce-miR-23928在东方蜜蜂微孢子虫侵染中的功能及调控机制提供依据。【方法】通过RNAhybrid、miRanda和TargetScan软件预测nce-miR-23928的靶基因。使用BLAST工具将上述靶基因比对到基因本体论(geneontology,GO)、京都基因和基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)、Nr和Swiss-Prot数据库以获得相应注释。采用实时荧光定量PCR(realtimequantitativePCR,RT-qPCR)技术检测nce-miR-23928及其靶基因在东方蜜蜂微孢子虫侵染意蜂工蜂过程中的相对表达量。【结果】相较于接种后1 d (1 day post infection, 1 dpi),nce-...     

东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程中nce-miR-12220及其靶基因的表达谱

东方蜜蜂微孢子虫Nosema ceranae侵染成年蜜蜂导致蜜蜂微孢子虫病。本研究旨在验证东方蜜蜂微孢子虫nce-miR-12220的存在和表达,并检测nce-miR-12220及其靶基因在病原侵染意大利蜜蜂Apis mellifera ligustica (简称意蜂)工蜂过程的表达谱。Stem-loop RT-PCR和Sanger测序结果显示nce-miR-12220真实存在和表达。靶向预测结果显示nce-miR-12220共靶向KRAB-A和γ tubulin等15个基因。上述靶基因可注释到19个GO条目和3条KEGG通路。RT-qPCR结果显示,相较于接种后1 d (1 day post infection,1 dpi),nce-miR-12220在2 dpi上调表达,而在3-12 dpi阶段总体表现出显著下调表达的趋势。类似地,与1 dpi相比,靶基因KRAB-A在2 dpi上调表达,而在3-12 dpi阶段总体呈下调表达的趋势。另外,与1 dpi相比,靶基因γ tubulin在2-12 dpi阶段总体表现出显著下调表达的趋势。上述结果表明nce-miR-12220与KRAB-A和γ tubulin之间存在潜在的靶向结合和正向调控关系;东方蜜蜂微孢子虫通过下调表达nce-miR-12220抑制KRAB-A的表达进而促进增殖;意蜂工蜂可能通过抑制东方蜜蜂微孢子虫的γ tubulin表达抵御病原侵染。研究结果明确了nce-miR-12220及其靶基因KRAB-A和γ tubulin在东方蜜蜂微孢子虫侵染意蜂工蜂过程中的动态表达规律,为深入探究nce-miR-12220在病原侵染中的功能及调控机制提供了理论和实验依据。     

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

【目的】本研究旨在通过差异表达基因(differentially expressed gene, DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log₂(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 ...