"神蜂"牌J9311蜂胶囊抗癌作用试验研究

以J931 1蜂胶囊水溶液为样品 ,对小鼠进行了抑制S1 80 实体瘤 (小鼠肉瘤细胞系 )和H2 2 实体瘤 (小鼠肝癌细胞系 )试验以及成分合理性分析试验。试验结果显示 ,对于S1 80 实体瘤 ,J931 1蜂胶囊的有效抑瘤剂量范围是 :0 6g kg～ 4 8g kg。以0 3g kg、0 6g kg、1 2g kg为低、中、高剂量进行抑瘤试验显示 ,中、高剂量组小鼠瘤重明显降低 ,与对照组比较差异显著 (P <0 0 5 ,P <0 0 1 ) ;对于H2 2 实体瘤的有效抑瘤剂量范围是 :1 6g kg～ 6 0g kg ,重复试验差异性结果与S1 80 实体瘤相同。结果表明J931 1蜂胶囊具有抑制小鼠S1 80 和H2 2 实体瘤的作用。J931 1蜂胶囊成分合理性试验显示 ,J931 1蜂胶囊组分抑瘤效果显著优于单独成分及各成分二二组合 ,单独成分中蜂王浆和蜂胶起主要抑瘤作用 ,中药成分起辅助作用     

西方蜜蜂的基因注释信息优化及新基因鉴定

西方蜜蜂(Apis mellifera)全基因组序列早在2006年就已公布,但其基因注释信息并不完善。本研究利用前期获得的意蜂幼虫肠道RNA-seq数据对西方蜜蜂基因组的已知基因进行结构优化,对新基因进行预测、鉴定和分析。利用Cufflink软件进行转录本重构,进而与参考转录本序列进行比较,共对4 177个已公布西方蜜蜂基因的测序结果进行了优化,其中5′端延长、3′端延长、5′和3′端都延长的基因分别有2 383、2 392和2 384个。利用Cuffcompare软件将重构转录本与参考基因组进行比对,共预测出197个西方蜜蜂新基因,随机挑选20个新基因进行RT-PCR验证,有18对引物成功扩增出符合预期的目的片段。GO分类结果显示这些新基因涉及运输和代谢、能量代谢和信号分子和交互等21个GO terms。KEGG代谢通路(pathway)富集分析结果显示这些新基因富集于13个pathways,其中基因富集数最多的是细胞进程、单一有机体和结合。本研究为西方蜜蜂基因组信息提供了有益的补充和完善,也为新基因的功能研究打下了基础。     

利用转录组数据开发意大利蜜蜂的SSR分子标记

意大利蜜蜂(简称意蜂)是自然界重要的授粉昆虫,广泛用于世界各国的养蜂生产,具有很高的生态和经济价值。为了开发意蜂的SSR分子标记,利用MISA软件对基于前期获得的意蜂幼虫肠道转录组数据组装得到的12 115条unigenes进行搜索,共预测出分布于2 149条unigenes的6 312个SSR位点,其中主要的重复类型为二核苷酸重复(54.42%)和三核苷酸重复(32.49%),主要基元为AT/AT(31.2%)和AG/CT(18.6%)。进一步利用软件设计出18 444对特异性引物,随机选取24对引物对国内3个不同来源的意蜂幼虫样品进行SSR位点扩增,有22对成功扩增出目的片段。研究开发出的SSR分子标记可用于意蜂的种群遗传和分子进化等研究,结果表明利用转录组数据鉴定非模式生物SSR位点的方法可行且高效。     

蜜蜂球囊菌菌丝和孢子中微小RNA及其靶mRNA的比较分析

【目的】蜜蜂球囊菌(Ascosphaera apis,球囊菌)专性侵染蜜蜂幼虫而导致白垩病。本研究旨在通过small RNA-seq（sRNA-seq）技术和生物信息学方法对球囊菌纯化菌丝（AaM）和纯化孢子（AaS）进行深度测序和比较分析,明确球囊菌菌丝miRNA和孢子miRNA的数量、结构和表达谱差异,并揭示菌丝和孢子共有miRNA、特有miRNA和差异表达miRNA（differentially expressed miRNA,DEmiRNA）及其靶mRNA与球囊菌菌丝和孢子生长、发育和病原致病性的潜在关系。【方法】实验室条件下获得纯培养的球囊菌,利用sRNA-seq技术对AaM和AaS分别进行测序,通过对原始读段（raw reads）进行过滤和质控获得有效标签序列（clean tags）。通过Venn分析筛选菌丝和孢子共有miRNA和特有miRNA。根据P≤0.05且|log₂ fold change|≥1的标准筛选AaM vs AaS的DEmiRNA。对上述共有miRNA、特有miRNA和DEmiRNA的靶mRNA进行预测,并对靶mRNA进行GO及KEGG数据库注释。根据靶向结合关系构建DEmiRNA和靶mRNA的调控网络。利用RT-qPCR验证测序数据的可靠性。【结果】AaM和AaS中分别得到12 982 320和12 708 832条raw reads,经过滤和质控分别得到10 800 101和9 888 848条clean tags。AaM中miRNA的长度介于18—26 nt,AaS中miRNA的长度介于18—24 nt,分布miRNA数量最多的长度均为18 nt,AaM和AaS中首位碱基为U的miRNA数量最多。AaM和AaS中表达量最高的miRNA均为miR6478-x、miR10516-x和miR482-x。菌丝和孢子共有miRNA靶向结合5 946个mRNA,二者特有miRNA分别靶向结合6 141和6 346个mRNA。共有miRNA的靶mRNA主要参与代谢进程、细胞进程和催化活性等42个功能条目,以及翻译、碳水化合物代谢和能量代谢等120条通路。AaM vs AaS比较组包含93个DEmiRNA,可靶向结合6 090个mRNA,这些靶mRNA可注释到38个功能条目和120条通路。DEmiRNA与靶mRNA之间形成较为复杂的调控网络,miR-4968-y位于调控网络的中心且能够靶向结合多达118个mRNA。RT-qPCR结果显示5个DEmiRNA的表达趋势与测序数据一致,证实了本研究中测序数据的可靠性。【结论】球囊菌菌丝和孢子中的miRNA具有类似的结构特征,但表达谱表现出明显差异;菌丝和孢子可能通过特异性表达和差异表达部分miRNA对其生长、发育和生殖进行调控。     

微小RNA介导东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂的分子机制

【目的】通过对东方蜜蜂微孢子虫（Nosema ceranae）纯化孢子与侵染意大利蜜蜂（Apis mellifera ligustica,简称意蜂）工蜂的东方蜜蜂微孢子虫的差异表达miRNA（DEmiRNA）及其靶mRNA进行系统分析,筛选、分析和探讨病原毒力因子和侵染因子相关的DEmiRNA及调控网络,在miRNA组学层面揭示东方蜜蜂微孢子虫对意蜂的侵染机制。【方法】利用small RNA-seq（sRNA-seq）技术对东方蜜蜂微孢子虫感染7 d和10 d的意蜂工蜂中肠和东方蜜蜂微孢子虫纯化孢子（NcCK）进行深度测序,通过连续比对rRNA数据库、西方蜜蜂（Apis mellifera）基因组和东方蜜蜂微孢子虫基因组筛滤出处于侵染过程的东方蜜蜂微孢子虫（NcT1和NcT2）数据和东方蜜蜂微孢子虫孢子的测序数据。根据P≤0.05,|log₂ fold change|≥1的标准,通过比较分析筛选出各比较组中的差异表达miRNA（differentially expressed miRNA,DEmiRNA）。通过相关生物信息学软件对DEmiRNA进行表达谱分析,靶mRNA预测及功能和代谢通路注释,以及调控网络的构建与分析。通过Stem-loop RT-qPCR验证DEmiRNA的差异表达趋势及测序数据的可靠性。【结果】NcCK vs NcT1、NcCK vs NcT2和NcT1 vs NcT2比较组分别包含164、122和60个DEmiRNA。Venn分析结果显示,3个比较组共有的上调和下调miRNA分别为5和6个。上述DEmiRNA分别预测出1 885、1 733和1 524个靶mRNA。这些靶mRNA分别注释到27、25和26个功能条目,其中注释数量最多的是新陈代谢进程、催化活性、细胞进程、结合和细胞。上述靶mRNA可分别注释到84、84和84条代谢通路,其中注释数量最多的是代谢途径、核糖体和次级代谢产物生物合成。此外,对于NcCK vs NcT1、NcCK vs NcT2和NcT1 vs NcT2中的DEmiRNA,分别有35、26和12个靶向结合MAPK信号通路相关靶mRNA,分别有49、40和17个DEmiRNA靶向结合糖酵解/糖异生通路相关靶mRNA。进一步分析发现,东方蜜蜂微孢子虫的DEmiRNA参与调控蓖麻毒素B凝集素、细胞凋亡抑制因子、极管蛋白和孢壁蛋白等病原毒力因子的基因表达,以及己糖激酶、ATP/ADP移位酶、ABC转运蛋白和转录因子ste12等侵染因子的基因表达。【结论】通过对东方蜜蜂微孢子虫纯化孢子与侵染意蜂工蜂的东方蜜蜂微孢子虫进行深入细致的miRNA组学分析和探讨,解析了病原侵染过程的miRNA差异表达谱,揭示了东方蜜蜂微孢子虫可能通过调节相应miRNA的表达水平对蓖麻毒素B凝集素、细胞凋亡抑制因子、极管蛋白和孢壁蛋白等毒力因子及己糖激酶、ATP/ADP移位酶、ABC转运蛋白和转录因子ste12等侵染因子基因表达进行调控,从而适应宿主细胞内的环境并促进自身的增殖与侵染。     

中华蜜蜂幼虫肠道中微小RNA的鉴定及分析

【目的】利用small RNA-seq（sRNA-seq）技术和生物信息学方法对中华蜜蜂（Apis cerana cerana,简称中蜂）幼虫肠道的微小RNA（microRNA,miRNA）进行全转录组鉴定和分析,旨在丰富中蜂的miRNA信息,并为深入研究miRNA调控中蜂幼虫肠道发育的分子机理提供依据。【方法】利用sRNA-seq技术对中蜂4、5和6日龄幼虫肠道样品（Ac1、Ac2和Ac3）进行测序,通过数据质控获得有效标签序列（clean tags）。采用Blast工具将clean tags连续比对东方蜜蜂（Apis cerana）基因组和miRBase数据库,以鉴定保守miRNA和新miRNA。采用TPM法对miRNA的表达量进行归一化处理。通过GraphPad Prism 7软件统计各组肠道样品中sRNA占比、miRNA长度分布及首位碱基偏向性。利用相关软件预测上述miRNA靶向的mRNA并进行GO和KEGG数据库注释。进一步根据靶向结合关系构建和分析注释到发育和免疫相关通路的基因及其靶向miRNA的调控网络,并利用Cytoscape软件进行可视化。利用茎环反转录PCR（Stem-loop RT-PCR）、分子克隆和Sanger测序验证miRNA的表达和序列的真实性。【结果】共鉴定到中蜂的371个保守miRNA和64个新miRNA;这些miRNA的长度介于18—25 nt且首位碱基主要偏向于U;上述miRNA共靶向14 750条mRNA,涉及离子结合、金属离子结合、细胞膜、细胞膜组件和单一有机体进程等2 270个GO条目,以及内吞作用、细胞凋亡、mTOR信号通路、RNA转运和昆虫激素的生物合成等332条KEGG通路。进一步分析结果显示156个miRNA与注释到Wnt、Hippo、Notch和mTOR等生长发育相关通路的67个靶基因存在调控关系,145个miRNA与注释到Toll、Imd/JNK、Jak-STAT和抗菌效应因子等免疫相关途径的21个靶基因存在调控关系。Stem-loop RT-PCR 结果显示miR-8-y、miR-9-z、miR-14-y、miR-281-y、miR-283-x和miR-306-x均能扩增出预期的特异性片段;Sanger测序结果显示上述6个miRNA的序列与深度测序结果一致。【结论】提供了中蜂miRNA的数量、结构特征和表达谱;揭示中蜂幼虫肠道miRNA潜在调控诸多生命进程与细胞活动;中蜂幼虫肠道的部分miRNA可通过靶向结合相应的mRNA参与调节发育和免疫相关途径。     

微小RNA介导意大利蜜蜂工蜂对东方蜜蜂微孢子虫的跨界调控

【目的】东方蜜蜂微孢子虫（Nosema ceranae）感染意大利蜜蜂（Apis mellifera ligustica,简称意蜂）导致蜜蜂微孢子虫病。本研究结合前期已获得的miRNA和mRNA组学数据,通过生物信息学方法对意蜂工蜂中肠的差异表达miRNA（differentially expressed miRNA,DEmiRNA）靶向结合的东方蜜蜂微孢子虫的mRNA和差异表达mRNA（DEmRNA）进行预测、数据库注释和调控网络分析,以期在组学水平解析miRNA介导意蜂工蜂对东方蜜蜂微孢子虫的跨界调控机制。【方法】通过比较东方蜜蜂微孢子虫侵染7 d和10 d的意蜂工蜂中肠（AmT1、AmT2）和未受侵染的工蜂中肠（AmCK1、AmCK2）的miRNA组学数据筛选出宿主的显著性DEmiRNA,通过比较侵染意蜂工蜂中肠的东方蜜蜂微孢子虫（NcT1、NcT2）和东方蜜蜂微孢子虫纯净孢子（NcCK）的mRNA数据筛选出病原的DEmRNA。利用TargetFinder软件预测宿主显著性DEmiRNA靶向结合的病原mRNA和DEmRNA。利用相关生物信息学工具对上述靶DEmRNA进行GO和KEGG数据库注释。结合前期研究结果筛选出孢壁蛋白、极管蛋白、蓖麻毒素B凝集素、ABC转运蛋白、ATP/ADP移位酶和糖酵解/糖异生途径等毒力因子和能量代谢通路相关的病原DEmRNA及与其存在靶向结合关系的宿主显著性DEmiRNA,并构建和分析二者的调控网络。【结果】AmCK1 vs AmT1比较组中宿主的48条显著上调miRNA和36条显著下调miRNA分别靶向病原的1 345和1 046条mRNA;进一步分析发现,宿主的47条显著上调miRNA和34条显著下调miRNA可分别靶向NcCK vs NcT1比较组中病原的584条显著下调mRNA和265条显著上调mRNA,它们可分别注释到19和22个功能条目以及66和64条通路。AmCK2 vs AmT2比较组中宿主的56条显著上调miRNA和51条显著下调miRNA分别靶向病原的1 260和1 317条mRNA;进一步分析发现,宿主的52条显著上调miRNA和49条显著下调miRNA可分别靶向NcCK vs NcT2比较组中病原的587条显著下调mRNA和336条显著上调mRNA,它们可分别注释到20和23个功能条目以及64和65条通路。AmCK1 vs AmT1和AmCK2 vs AmT2比较组的8条共同显著上调miRNA和1条共同显著下调miRNA分别靶向NcCK vs NcT1和NcCK vs NcT2比较组中的144条共同显著下调和10条共同显著上调mRNA,可分别注释到18和13个功能条目以及38和7条通路。此外,AmCK1 vs AmT1和AmCK2 vs AmT2比较组中宿主的显著上调miRNA可靶向结合NcCK vs NcT1和NcCK vs NcT2比较组中与RNAi途径,孢壁蛋白和蓖麻毒素B凝集素等毒力因子,糖酵解/糖异生途径以及MAPK信号通路相关的病原下调表达mRNA。【结论】在东方蜜蜂微孢子虫的侵染过程中,意蜂工蜂中肠的DEmiRNA与病原的DEmRNA之间存在复杂的靶向结合关系以及潜在的跨界调控关系;宿主的DEmiRNA可能通过抑制或降解病原的RNAi途径、毒力因子、糖酵解/糖异生通路、ATP/ADP移位酶、ABC转运蛋白及MAPK信号通路相关靶DEmRNA影响病原的侵染和增殖。     

蜜蜂病毒的传播方式

病毒是一类专性的细胞内寄生物,几乎从所有的生命形式中都发现了病毒。到目前为止,已经在蜜蜂上至少分离鉴定出18种病毒[1]。近十年里,蜜蜂病毒传播方式的研究一直是发展较为快速的领     

基于比较转录组学分析揭示中华蜜蜂及意大利蜜蜂幼虫的球囊菌抗性差异机制

为探究中华蜜蜂(Apis cerana cerana,以下简称中蜂)与意大利蜜蜂(Apis mellifera ligustica,以下简称意蜂)幼虫的球囊菌(Ascospahaera apis)抗性差异产生的原因,利用Venn分析、GO分类分析、KEGG代谢通路分析以及Ka/Ks分析对二者的转录组进行比较研究。Venn分析结果显示,中蜂与意蜂的同源基因有17 656个,归属于8 111个基因家族,二者的特有基因分别有1 158和241个,分别归属于468和86个基因家族。GO分类结果显示,中蜂和意蜂的特有基因分别富集在38和28个GO条目。KEGG代谢通路富集分析结果显示,中蜂和意蜂的特有基因分别富集于96和21个代谢通路;进一步分析发现中蜂的特有基因富集在内吞作用、溶酶体、泛素介导的蛋白水解和黑化作用等4个细胞免疫通路,以及MAPK信号通路和Toll-like受体信号通路等2个体液免疫通路,而意蜂仅有1个特有基因富集在MAPK信号通路。Ka/Ks分析结果显示受到强烈正向选择、弱正向选择、中性选择和负向选择的单拷贝同源基因分别有37、281、2 630和3 269个。对受到强烈正向选择的基因进行GO分类和KEGG代谢通路富集分析,GO结果显示中蜂和意蜂的单拷贝同源基因富集的GO条目类型相同;KEGG结果显示中蜂的单拷贝基因仅富集在氧化磷酸化。上述结果表明免疫相关基因数量的差异是二者的球囊菌抗性差异产生的重要原因之一,中蜂在与球囊菌的协同进化过程中可能通过提高能量利用从而限制球囊菌的增殖。本研究结果可为阐明中蜂及意蜂幼虫的球囊菌抗性差异产生的分子机制提供参考。     

蜜蜂球囊菌基因结构优化及新基因鉴定

为补充和完善蜜蜂球囊菌Ascosphaera apis基因组序列和功能注释信息,本研究利用RNA-seq技术对球囊菌菌丝和孢子进行深度测序,并基于高质量的转录组数据对已注释基因进行结构优化,对未注释基因进行预测、鉴定和分析。将测序得到的有效读段进行参考基因组比对和转录本重构;利用Cuffcompare软件将重构转录本与参考基因组进行比对。结果表明:共对101个已注释基因的5′端或3′端进行了延长;共鉴定出373个新基因,随机挑选10个新基因进行RT-PCR验证,其中8个能扩增出符合预期的目的片段,表明预测出的多数新基因真实存在;共计147个球囊菌新基因可注释到Nr和eggNOG数据库中,85个新基因注释到29个GO条目,66个新基因注释到33条代谢通路。本研究结果补充了球囊菌的基因结构和功能注释信息,也为新基因的功能研究打下了初步基础。