基于纳米孔全长转录组数据完善东方蜜蜂微孢子虫的基因组注释

【目的】利用已获得的纳米孔全长转录组数据对现有的东方蜜蜂微孢子虫（Nosema ceranae）参考基因组的基因序列和功能注释进行完善。【方法】采用TransDecoder软件预测东方蜜蜂微孢子虫基因的开放阅读框（open reading frame,ORF）及相应的氨基酸。利用gffcompare软件将全长转录本与参考基因组注释的转录本进行比较,对基因组注释基因的非编码区向上游或下游延伸,修正基因的边界。利用MISA软件鉴定长度在500 bp以上的全长转录本的简单重复序列（simple sequence repeat,SSR）位点,包括单核苷酸重复、双核苷酸重复、三核苷酸重复、四核苷酸重复、五核苷酸重复、六核苷酸重复、混合SSR等类型。通过Blast工具将鉴定到的新基因和新转录本比对Nr、KOG、eggNOG、GO和KEGG数据库,从而获得功能注释。【结果】共预测出2 353个完整ORF,其中长度分布在0—100个氨基酸的ORF最多,占总ORF数的72.12%。共对东方蜜蜂微孢子虫的2 340个基因进行了结构优化,其中5′端延长的基因有1 182个,3′端延长的基因有1 158个。共鉴定到1 658个SSR,其中单核苷酸重复、双核苷酸重复、三核苷酸重复、四核苷酸重复的数量分别为1 622、23、7和6个;单核苷酸重复类型的SSR密度最大,达到182.32个/Mb,其次为混合SSR、双核苷酸重复和三核苷酸重复,分别达到6.90、2.78和0.73个/Mb。共鉴定出954个新基因,其中分别有951、333、371、422和321个新基因可注释到Nr、KOG、eggNOG、GO和KEGG数据库。此外,还鉴定出6 164条新转录本,其中分别有6 141、2 808、2 932、3 196和2 585条新转录本可注释到Nr、KOG、eggNOG、GO和KEGG数据库。新基因和新转录本注释数量最多的物种均为东方蜜蜂微孢子虫,其次是蜜蜂微孢子虫（Nosema apis）。【结论】研究结果较好地完善了现有的东方蜜蜂微孢子虫参考基因组已注释基因的序列和功能注释,并补充和注释了大量参考基因组未注释的新基因和新转录本。     

营养生态条件对蜜蜂球囊菌生长及产孢的影响

从碳源、氮源、维生素、矿质元素等方面研究了营养生态条件对蜜蜂球囊菌(Ascosphaera apis Olive & Spiltoir)生长及产孢的影响.结果表明,营养生态条件的变化对蜜蜂球囊菌的影响极大.A.apis菌丝体生长最好的碳源是果糖,产孢量最高的是葡萄糖,不利用山梨糖和甘露糖;该菌优先利用有机氮源 ,不利用无机氮源中NaNO2、H2NCSNH2;维生素对A.apis 产孢量有极为显著的影响,复合维生素比单一维生素能更好地促进菌丝体生长, 尤其对产孢更为有利;矿质元素中Mg2+、P5+、K+、Zn2+促进菌丝体生长和孢子的形成,Na+对菌丝体的生长有抑制作用.     

福州地区蜜蜂养殖与蜂蜜安全生产现状调研报告

本文就福建省福州地区2013年蜜蜂养殖现状、蜂产品市场现状、消费者关注蜂产品的因素、会员对协会功能发挥的期望等进行了调查,据此在蜂农产品销售方面、蜂农组织联合方面、蜂专业协会组织引导方面提出了未来发展策略。     

蜜蜂白垩病的研究进展（连载二）

本文就蜜蜂白垩病的病原及其生长条件、症状与诊断、防治等方面作了简要的述评。     

中国养蜂学会蜜蜂饲养管理专业委员会第十二次学术研讨会暨工作会议在四川省眉山市召开

当我拜读了晋华贵老师发表在<蜜蜂杂志>2006年第9期17～18页上的文章<谈中蜂的活框饲养与传统继承>后,顿感眼前一亮.文中通篇以感性知识为背景,通过否定之否定视角,将中蜂饲养的科学方法进行精辟提炼和论述,源于实践而又高于实践,具有极强的指导意义,说出了众蜂友想说而又表达不出的心里话,是篇难得的好文章.     

蜜蜂球囊菌分泌多种胞外酶侵染蜜蜂幼虫

对健康和患病蜜蜂幼虫血淋巴进行SDS-PAGE电泳和蛋白酶、酯酶的活性染色,以研究蜜蜂球囊菌分泌的胞外酶种类及对蜜蜂侵染的作用.结果表明:健康蜜蜂幼虫血淋巴中的蛋白含量丰富,主要由4种高分子质量的蛋白组成,而患病幼虫血淋巴中的蛋白含量很少,主要蛋白组分被降解;在患病幼虫血淋巴中可以检测到健康幼虫所没有的多种蛋白酶和酯酶...     

意大利蜜蜂幼虫肠道在球囊菌侵染前期的 差异表达microRNA及其调控网络

【目的】 微小RNA（microRNA,miRNA）是一类重要的基因表达调控因子,可影响宿主与病原间的互作过程。蜜蜂球囊菌（Ascosphaera apis）是一种特异性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在对意大利蜜蜂（Apis mellifera ligustica,简称意蜂）幼虫肠道在球囊菌胁迫前期的差异表达miRNA（differentially expressed miRNA,DEmiRNA）及其靶基因进行深入分析,在miRNA组学水平探究意蜂幼虫在球囊菌侵染前期的胁迫应答,并通过构建显著DEmiRNA的调控网络筛选出与宿主应答相关的关键miRNA。【方法】 利用small RNA-seq（sRNA-seq）技术对正常及球囊菌侵染的意蜂4日龄幼虫肠道（AmCK和AmT）进行高通量测序,首先对原始数据进行质控和评估,随后将过滤后的数据与西方蜜蜂（Apis mellifera）参考基因组进行比对;将比对上的序列标签（tags）注释到miRBase数据库,得出已知miRNA的表达量;通过TPM（tags per million）算法对各样本中miRNA的表达量进行归一化处理,以|log₂ fold change|≥1且P≤0.05作为标准筛选得到显著DEmiRNA;利用TargetFinder 软件预测显著DEmiRNA的靶基因,并对其进行GO和KEGG代谢通路（pathway）富集分析。利用Cytoscape软件对miRNA-mRNA调控网络进行可视化。最后,利用茎环反转录PCR（Stem-loop RT-PCR）和荧光定量PCR（qPCR）验证测序数据的可靠性。【结果】 AmCK和AmT样品的测序分别得到13 553 302和10 777 534条原始读段（raw reads）,经严格过滤后得到的有效读段（clean reads）数分别为13 186 921和10 480 913条。各样品的生物学重复间的Pearson相关性系数分别在0.9822和0.9508以上。共有10个显著DEmiRNA,包括4个上调miRNA和6个下调miRNA。显著DEmiRNA在AmT的整体表达水平低于AmCK。10个显著DEmiRNA可靶向结合3 788个靶基因,其中上调miRNA的1 240个靶基因可注释到GO数据库中的39个GO条目,主要富集在结合、细胞进程、代谢进程和应激反应等;下调miRNA的749个靶基因可注释到34个GO条目,主要富集在细胞进程、结合、代谢进程和应激反应等。KEGG数据库注释结果显示,上调miRNA和下调miRNA的靶基因分别注释到95和66条代谢通路,富集基因数最多的分别是Wnt信号通路、Hippo信号通路、光传导以及内吞作用、磷脂酰肌醇信号系统、嘌呤代谢。对于上调和下调miRNA,分别有31和52个靶基因注释到内吞作用,15和7个靶基因注释到泛素介导的蛋白水解,11和5个靶基因注释到Jak-STAT信号通路,1和3个靶基因注释到MAPK信号通路。显著DEmiRNA与靶mRNA之间形成复杂的调控网络,7个显著DEmiRNA靶向结合96个与Wnt信号通路相关的mRNA,8个显著DEmiRNA靶向结合55个与内吞作用相关的mRNA。Stem-loop RT-PCR和qPCR结果验证了测序数据的可靠性。【结论】 对意蜂幼虫肠道在球囊菌侵染前期的DEmiRNA及其靶基因进行预测和分析,并构建和分析了DEmiRNA-mRNA调控网络,研究结果提供了宿主miRNA的表达谱和差异表达信息,揭示了DEmiRNA通过调控细胞生命活动、新陈代谢以及部分细胞和体液免疫等生物学过程参与宿主的胁迫应答。miR-4331-y、miR-4968-y、miR-8440-y、novel-m0023-5p和novel-m0025-3p共同参与了宿主的Wnt信号通路和内吞作用的调控,可作为白垩病治疗的潜在分子靶标。     

感染意大利蜜蜂工蜂的东方蜜蜂微孢子虫及其纯化孢子的高表达基因分析

为探究东方蜜蜂微孢子虫(Nosema ceranae)在侵染意大利蜜蜂(Aips mellifera ligustica)过程的基因表达谱及高表基因(highly expressed gene,HEG)的功能,本研究利用RNA-seq技术对N.ceranae感染的意蜂工蜂7和10 d工蜂中肠(NcT1和NcT2)及其纯化孢子(NcCK)进行高通量测序。根据FPKM值≧15的标准筛选出各组的HEG并进行Venn分析,进而通过GO分类和KEGG代谢通路富集分析对各组的共有HEG进行功能注释。测序得到的原始读段经严格质控后分别得到5 816 465、28 501 225和210 824 312条高质量的有效读段,分别从NcT1、NcT2和NcCK中筛选出1 263、1 233和1 511个HEG,其中有1 079个基因在各组均高量表达,分别有15和7个基因仅在NcT1和NcT2中高量表达。GO分类结果显示,共有HEG涉及9个生物学进程相关条目,10个细胞组分相关条目,以及5个分子功能相关条目。KEGG代谢通路富集分析结果表明,共有HEG富集在218条代谢通路,富集基因数最多的是核糖体、Epstein-Barr病毒感染、细胞周期、内质网蛋白加工和真菌核糖体生物合成。深入分析结果显示分别有5和5个共有HEG富集在MAPK和cAMP信号通路,表明此2条信号通路及富集的HEG可能在N.ceranae的环境应激和侵染过程中发挥重要作用。研究结果提供了N.ceranae感染意蜂工蜂过程中的基因表达谱,揭示了HEG的作用,为深入解析N.ceranae的侵染机制打下了初步基础。     

基于RNA-seq技术分析意大利蜜蜂幼虫肠道响应球囊菌胁迫的高表达基因

对意蜂幼虫肠道组织基因进行深度RNA测序,并对高表达基因进行GO注释和KEGG分析,旨在揭示意蜂幼虫免疫抗病机制。分别取3、4、5、6日龄意蜂幼虫中肠组织抽提RNA,合成cDNA。得到的总读段数都在26 715 638条以上,有效数据达到97%以上,且测序饱和度良好,说明测序所得数据真实可信。将测序有效数据和意蜂公布DNA数据库进行比对,发现意蜂幼虫肠道中有13 789个基因有不同程度表达。选取各样本FPKM排序最前且表达水平最高的100个基因进行GO注释,这里主要涉及细胞组分、分子功能和生物学进程功能分析。利用KEGG分析显示这100个高表达基因分布在能量代谢、翻译、组装和降解以及信号转导等多条通路上。初步分析相应共有基因和特有基因在相关通路上的作用,并对意蜂幼虫肠道受球囊菌胁迫后的相关免疫机制进行一定分析。     

中华蜜蜂幼虫肠道参考转录组的de novo组装及SSR分子标记鉴定

【目的】利用RNA seq技术对中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫肠道参考转录组进行de novo组装,并进行功能及代谢通路注释,进而利用该转录组数据进行中蜂幼虫的SSR分子标记鉴定。【方法】实验室条件下饲养中蜂幼虫,将纯化的蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)孢子饲喂3日龄幼虫,剖取4、5和6日龄幼虫肠道,液氮速冻。将健康幼虫肠道与感染球囊菌的幼虫肠道同时进行Illumina测序。通过对raw reads的过滤得到clean reads,利用Trinity软件组装得到unigenes。通过BLASTx(E-value10-5)比对NCBI Nr、Swiss-Prot、KOG和KEGG数据库,对unigenes进行功能和代谢通路注释。利用MISA软件对所有unigenes进行SSR搜索,并利用Primer Premier 5软件设计特异性SSR引物,通过常规PCR对来源于北京、辽宁兴城和四川成都的中蜂幼虫肠道样本进行SSR位点鉴定。【结果】中蜂幼虫肠道的RNA seq共得到35 670 000条reads,de novo组装得到43 557个unigenes,平均长度为898 nt。共有18 225个unigenes可被注释到上述公共蛋白数据库,单独注释到NCBI Nr、Swiss-Prot、KOG和KEGG数据库的unigenes数分别为3 899、443、37和10个。KOG注释结果显示,11 442条unigenes分布于25个基因家族,其中注释到RNA加工和修饰家族的基因数最多,达1 249个。9 679个unigenes的GO分类结果显示,在生物学进程分类中,注释到细胞进程的基因最多,达4 201个,在分子功能和细胞组分类中,注释到结合与细胞的基因数最多,分别为4 935和2 900个。4 517个unigenes可注释到KEGG数据库中的216个代谢通路,注释到核糖体的基因数最多,达385个。利用MISA软件,可在7 763个unigenes搜索到13 448个SSR位点,随机选取20对SSR引物对国内3个不同来源的中蜂幼虫肠道样本的SSR位点进行扩增,有6对引物可鉴定出SSR分子标记。【结论】成功组装并注释了中蜂幼虫肠道参考转录组,可为中蜂及其幼虫的分子生物学及组学研究提供重要的参考信息,也可用于补充、丰富和检验东方蜜蜂的参考基因组,基于此转录组数据开发出6个中蜂的SSR分子标记,可应用于中蜂的基因图谱构建、基因多样性分析、基因定位等研究,也说明利用转录组数据开发非模式生物SSRs的方法可行。