意大利蜜蜂幼虫响应东方蜜蜂微孢子虫胁迫的高表达基因分析

为探究意大利蜜蜂(Apis mellifera ligustica,简称意蜂)幼虫响应东方蜜蜂微孢子虫(Nosema ceranae)胁迫的高表达基因(HEG)的作用,利用RNA-seq技术对正常的意蜂4、5和6日龄幼虫(Am4CK、Am5CK和Am6CK)及N.ceranae胁迫的意蜂4、5和6日龄幼虫(Am4T、Am5T和Am6T)进行测序。共得到268 765 544条原始读段,经过滤得到的有效读段数为266 318 442条,各样品的Q30均在92.42%及以上。根据FPKM值≥15的标准从上述6个样品中分别筛选出2 460、2 692、2 908、2 367、2 549和2 270个HEG。Venn分析结果显示,去除与对照组共有的HEG,Am4T、Am5T和Am6T的共有HEG数为3个,特有HEG数分别为373、210和79个。GO分类结果显示,Am4T、Am5T和Am6T的特有HEG分别能够注释26、27和16个GO条目,注释基因数最多的分别是代谢进程、结合和结合。KEGG分析结果显示这些特有HEG可分别富集在145、138和68条代谢通路,富集基因数最多的分别是过氧化物酶体、碳代谢和内质网蛋白加工。进一步分析结果表明N.ceranae胁迫可导致宿主的部分物质代谢、能量代谢和免疫代谢通路水平提高。从6个样品的共有HEG中随机选取10个进行RT-PCR验证,均能成功扩增出目的条带,说明本研究中的HEG真实存在。研究结果揭示了意蜂幼虫响应N.ceranae胁迫的HEG表达谱及作用,为深入解析宿主的胁迫应答机制提供了重要信息。     

意大利蜜蜂4、5和6日龄幼虫肠道发育过程中差异表达基因的趋势分析

为探讨意大利蜜蜂(简称意蜂)幼虫肠道发育过程中的基因表达谱和基因差异表达信息,研究基于前期已获得的高质量的意蜂幼虫肠道转录组数据,利用STEM软件对意蜂4、5和6日龄幼虫肠道的差异表达基因(DEGs)进行趋势分析,进一步通过相关生物学软件对显著上调趋势包含的DEGs进行GO与KEGG代谢通路富集分析及表达量聚类分析。趋势分析结果显示,5 920个DEGs共聚类在8个基因表达模式,其中935个DEGs聚类在1个显著上调趋势。GO富集分析结果显示,显著上调趋势包含DEGs富集于39个GO条目,富集基因数最多的为结合,其次是细胞进程和单组织进程。KEGG代谢通路富集分析结果显示,显著上调趋势包含DEGs富集于94个代谢通路,富集基因数最多的是嘌呤代谢、内质网蛋白质加工和内吞作用。表达量聚类分析结果显示,随着意蜂幼虫肠道发育时间的延长,代谢和免疫相关通路富集基因的表达水平逐渐增强。研究提供了意蜂幼虫肠道发育过程的基因表达谱和基因差异表达信息,揭示了基因表达规律,为深入研究意蜂幼虫肠道发育的分子机制打下了一定基础。     

蜜蜂球囊菌中长链非编码RNA的调控作用

【目的】长链非编码RNA（long non-coding RNA,lncRNA）是一类二、三级结构高度保守,长度>200 nt且不具蛋白编码能力的RNA,在转录和转录后水平广泛参与调控剂量补偿、细胞分化和生长发育等生命活动。本研究基于前期获得的蜜蜂球囊菌（Ascosphaera apis,简称球囊菌）菌丝和孢子混合样品的高质量lncRNA组学数据进行球囊菌lncRNA的顺式（cis）作用、反义lncRNA（antisense lncRNA）作用和竞争性内源RNA（competing endogenous RNA,ceRNA）作用的分析和探讨,以期揭示lncRNA在球囊菌中的潜在功能。【方法】基于lncRNA基因在染色体上的位置,预测lncRNA上下游100 kb以内的蛋白编码基因;使用Blast软件将上下游基因比对到GO和KEGG数据库,以获得功能和通路注释。利用LncTar软件对反义lncRNA的靶mRNA进行预测,并使用Blast软件将上述靶mRNA比对到KEGG和eggNOG数据库。利用TargetFinder软件预测lncRNA靶向结合的miRNA及miRNA靶向结合的mRNA,根据靶向结合关系建立lncRNA-miRNA和lncRNA-miRNA-mRNA调控网络,进而通过Cytoscape v3.7.1软件进行调控网络的可视化。利用RT-PCR对调控网络中的lncRNA、靶miRNA和靶mRNA进行表达验证。【结果】共预测出371个lncRNA的5 852个上下游基因,可注释到细胞进程、代谢进程和应激反应等48个功能条目,以及新陈代谢途径、次生代谢产物的生物合成和抗生素的生物合成等121条通路,表明部分lncRNA可通过顺式作用调节上下游基因的表达,从而参与调控球囊菌的生长发育及物质能量代谢等基础生命活动。球囊菌的7个lncRNA与7个靶mRNA存在序列互补关系,其中5个mRNA在eggNOG数据库中仅注释为假定蛋白,仅gene3444在KEGG数据库注释为核孔复合体蛋白An-Nup120和假定蛋白,表明上述1个反义lncRNA可能参与调控核孔复合体蛋白的生物合成等生物学过程。此外,共预测出227个lncRNA与73个miRNA之间存在靶向结合关系,其中多数lncRNA(79.02%)仅能结合1—2个miRNA,部分miRNA可被多个lncRNA靶向结合;进一步构建氧化磷酸化通路和MAPK信号通路相关的lncRNA-miRNA-mRNA调控网络,分析结果显示氧化磷酸化通路相关的222个lncRNA靶向78个miRNA及50个mRNA,MAPK信号通路相关的222个lncRNA靶向76个miRNA及46个mRNA,表明部分lncRNA通过ceRNA作用调控此两条通路,从而影响球囊菌的能量合成、环境适应以及生长发育等过程。【结论】部分lncRNA可能通过顺式作用和ceRNA作用调节球囊菌的生长发育、物质能量代谢以及环境适应等生物学过程;MSTRG.5393.1可能作为反义lncRNA调控球囊菌的核孔复合体的蛋白合成。     

超声波辅助提取蜂胶总黄酮工艺的优化

[目的]研究超声波对蜂胶总黄酮提取得率的影响,优化蜂胶总黄酮提取工艺参数。[方法]选取料液比（每克粗蜂胶对应的乙醇毫升数）、乙醇浓度、超声波处理时间、超声波频率4个因素,以蜂胶总黄酮提取得率为指标进行蜂胶总黄酮提取试验。[结果]经单因素及正交试验分析,上述4个因素对蜂胶总黄酮提取得率影响的大小顺序为超声波频率〉乙醇浓度〉超声波时间〉料液比;蜂胶总黄酮提取的最佳工艺参数为：料液比1∶9 g/ml,乙醇浓度80%,超声波处理12 min,超声波频率45 kHz;经验证试验,在最佳条件下蜂胶总黄酮的提取得率为20.25%。[结论]超声波频率变化对蜂胶总黄酮提取得率影响较大,最佳取值为45 kHz,对实际生产具有理论指导意义。     

意大利蜜蜂工蜂中肠发育过程中长链非编码RNA的差异表达分析

【目的】长链非编码RNA(lncRNA)在真核生物的基因表达、表观遗传和细胞周期调控等方面发挥重要功能。本研究旨在探究意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂中肠发育过程中lncRNA的表达谱及其作用。【方法】利用RNA-seq技术和链特异性建库方法对意蜂7和10日龄工蜂中肠(Am7、Am10)进行深度测序,下机的原始数据经过Perl脚本过滤得到高质量有效读段。利用bowtie工具将有效读段比对核糖体数据库,进一步利用Top Hat2软件将未比对到核糖体数据库上的数据比对到参考基因组。利用CPC和CNCI软件对转录本的编码能力进行预测。通过RT-PCR对部分lncRNA进行鉴定。利用edgeR软件进行差异表达lncRNA(DElncRNA)分析,进而预测lncRNA的上下游基因,并对上下游基因进行GO及KEGG代谢通路富集分析。联用RNAhybrid、Miranda和Target Scan软件预测DElncRNA靶向结合的mi RNA及mi RNA靶向结合的靶基因,并通过Cytoscape软件构建DElncRNAs-mi RNAs-m RNAs的调控网络。最后,通过RT-qPCR验证测序数据的可靠性。【结果】Am7和Am10的深度测序分别获得134 802 058和147 051 470条原始读段,经严格过滤分别得到134 166 157和146 293 288条有效读段;共得到3 890个DElncRNA,包括2 005个上调lncRNA与1 885个下调lncRNA。RT-PCR验证结果显示共有8个lncRNA能扩增出符合预期的目的片段,表明预测出的lncRNA真实存在。DElncRNA的上下游基因数为1 793个,它们涉及42个GO条目,包括代谢进程、发育进程、细胞进程、应激反应和免疫系统进程等;这些上下游基因还涉及251条代谢通路,包括碳代谢、嘌呤代谢和脂肪酸的生物合成等物质代谢通路,硫代谢、甲烷代谢和氧化磷酸化等能量代谢通路,Hippo信号通路、Wnt信号通路和Notch信号通路等信号通路,溶酶体、内吞作用和泛素介导的蛋白水解等细胞免疫通路,以及MAPK信号通路、Jak-STAT信号通路和NF-kappa B信号通路等体液免疫通路,上述结果表明DElncRNA在意蜂中肠发育过程中参与物质和能量代谢、细胞生命活动和免疫调控。进一步分析发现TCONS_00020918可通过调控西方蜜蜂王浆主蛋白1编码基因在意蜂工蜂中肠的营养吸收、级型分化中发挥功能。DElncRNA的调控网络分析结果显示DElncRNA与mi RNA、m RNA间存在复杂的调控关系,部分DElncRNA处于调控网络的中心位置且能靶向结合较多的mi RNA,也有部分mi RNA可被多个DElncRNA共同靶向,表明这些DElncRNA可能在中肠发育中发挥重要作用。随机挑取5个DElncRNA进行RT-qPCR验证,结果显示它们的表达量变化趋势与测序结果一致,证实了本研究测序数据的可靠性。【结论】差异表达长链非编码RNA(DElncRNA)广泛参与意蜂工蜂中肠的新陈代谢、细胞活动和免疫调控并作为竞争性内源RNA(ce RNA)发挥作用,研究结果为关键lncRNA的筛选和功能研究提供了必要的数据支持。     

利用第三代纳米孔长读段测序技术构建和注释蜜蜂球囊菌的全长转录组

[目的]利用第三代纳米孔(nanopore)长读段测序技术对蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)的纯化菌丝(Aam)和孢子(Aas)进行测序,构建和注释球囊菌的高质量全长转录组.[方法]通过Oxford Nanopore PromethION平台对Aam和Aas进行测序.利用Guppy软件对原始读...     

蜂群群势对蜂王浆生产的影响

蜂王浆是意大利蜜蜂生产的主要蜂产品之一,也是养蜂者收入的主要经济来源之一。测试了群势大小对意大利蜂群生产蜂王浆的影响,结果表明:王台的接受率、单台产量、酸度、10-HDA含量均随着蜂群群势的增加呈上升趋势,蜂王浆中水分含量刚好相反,呈下降趋势。当蜂群群势达10框时,王台接受率较9框有极显著提升,达72.98%,11框较10框也有极显著提升,达82.91%,12框较11框增加不显著,达83.20%;当蜂群群势达10框时,单台产浆量较9框有极显著增加,达336.49 mg,11框较10框显著增加,达371.76mg,12框较11框略有降低,但差异不显著,为368.92 mg;蜂王浆酸度11框较10框、10框较9框均显著提高,12框较11框提高不显著;10-HDA含量11框较10框显著提高,其他虽整体呈增加趋势,但相邻间差异不显著;蜂王浆中水分含量10框较9框显著下降,其他相邻间差异不显著。     

油菜蜂花粉软胶囊复配乳化剂的研制

采用均匀设计的方法,筛选油菜蜂花粉软胶囊复配乳化稳定剂配方.采用SAS统计软件对试验结果进行多元回归分析,确定了最优配方为:蜂蜡6.5 g;单甘脂5.5 g;吐温0.35 g;司盘80 1.10 g;硬脂酸0.35 g.     

"神蜂"牌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,其中基因富集数最多的是细胞进程、单一有机体和结合。本研究为西方蜜蜂基因组信息提供了有益的补充和完善,也为新基因的功能研究打下了基础。