结合三代测序与二代测序技术揭示蜜蜂球囊菌孢子转录组的复杂性

【目的】蜜蜂球囊菌(Ascosphaera apis)是一种专性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在利用PacBio单分子实时(singlemoleculereal-time,SMRT)测序技术对蜜蜂球囊菌孢子(AaS)中基因的可变剪切(alternative splicing,AS)和可变多聚腺苷酸化(alternative polyadenylation,APA)以及长链非编码RNA (long non-coding RNA,lncRNA)进行鉴定和分析,进而揭示蜜蜂球囊菌孢子中转录组的复杂性。【方法】采用Suppa软件对蜜蜂球囊菌孢子中基因的AS事件进行鉴定。通过RT-PCR对不同类型的AS事件进行验证。采用TAPIS pipeline对蜜蜂球囊菌孢子基因的APA位点进行鉴定。利用MEME软件分析孢子全长转录本的poly(A)剪接位点上游50bp的序列特征并鉴定motif。联用CPC和CNCI软件和比对Swiss-prot数据库的方法预测lncRNA,取三者的交集作为高可信度的lncRNA集合。进一步比较lncRNA和mRNA的转录本长度,外显子数量与长度,内含子长度,GC含...     

西方蜜蜂AmAGO1蛋白的分子特性、时空表达谱及抗体制备

【目的】Argonaute (AGO)家族是一类在进化上高度保守的蛋白家族,其在真核生物中主要参与形成RNA诱导沉默复合体(RNA-induced silencing complex, RISC),进而通过沉默基因表达参与诸多生物学过程。蜜蜂AGO蛋白的相关研究迄今仍然缺失。本研究拟通过预测西方蜜蜂Apis mellifera的AGO1(AmAGO1)理化性质和分子特性,解析AmAGO1在西方蜜蜂中的时空表达谱,并制备AmAGO1的多克隆抗体,为深入开展AmAGO1的功能与机制研究提供参考和基础。【方法】PCR扩增西方蜜蜂AmAGO1的编码序列(coding sequence, CDS);通过生物信息学预测AmAGO1蛋白的理化性质和分子特性;利用RT-qPCR检测AmAGO1在西方蜜蜂工蜂卵、3日龄幼虫、7日龄预蛹、8日龄预蛹、12日龄蛹以及1, 2, 6, 12, 15和18日龄成虫以及工蜂成虫触角、咽下腺、脑、表皮、中肠、脂肪体和毒腺中的表达量;构建原核表达质粒后诱导表达AmAGO1融合蛋白,并鉴定其表达形式;制备AmAGO1多克隆抗体,利用ELISA, Western blot...     

东方蜜蜂微孢子虫孢子中微小RNA的鉴定与分析

【目的】丰富东方蜜蜂微孢子虫Nosema ceranae的微小RNA(microRNA, miRNA)信息,并为深入探究miRNA在病原孢子和病原侵染中的功能提供理论和实验依据。【方法】基于已获得的small RNA-seq数据,利用生物信息学软件对东方蜜蜂微孢子虫的纯净孢子中的miRNA进行鉴定和分析。采用茎环反转录PCR(stem-loop RT-PCR)检测已鉴定的miRNA的表达;通过分子克隆与Sanger测序验证miRNA的序列。使用TargetFinder软件预测这些miRNA的靶基因,并对靶基因进行数据库注释。根据miRNA与靶基因的靶向结合关系构建调控网络,再利用Cytoscape软件进行可视化。【结果】在东方蜜蜂微孢子虫孢子中共鉴定到10个miRNA;这些miRNA的长度分布介于21~25 nt,首位碱基表现出U偏向性,每一位碱基的偏向性差异明显。Stem-loop RT-PCR检测结果表明这10个miRNA均真实表达;Sanger测序结果证实了随机选取的其中2个miRNA的序列真实性。共预测出249个靶基因,其中分别有249, 118, 136和3个靶基因可注释到Nr,Swiss-Prot, KOG和eggNOG数据库。此外,分别有134和71个靶基因可分别注释到GO数据库的30个功能条目和KEGG数据库的54条通路。【结论】本研究揭示了东方蜜蜂微孢子虫孢子中miRNA的存在和表达;这些miRNA通过调控潜在靶基因的表达参与孢子的生命活动。     

东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程中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在病原侵染中的功能及调控机制提供了理论和实验依据。     

环状RNA ame_circ_000115调节蜜蜂球囊菌胁迫下意大利蜜蜂幼虫肠道基因表达

环状RNA (circular RNA, circRNA)是一类新型非编码RNA,已被证实可通过竞争性内源RNA(competing endogenous RNA,ceRNA)调控网络等多种方式参与调节昆虫基因表达和免疫应答。目前,circRNA在蜜蜂免疫应答中的作用尚不清楚。本研究对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)ame＿circ＿000115(简称ac115)的反向剪切(back splicing,BS)位点进行验证。采用RT-qPCR检测蜜蜂球囊菌(Ascosphaera apis)胁迫下意蜂幼虫肠道内ac115的表达谱,利用双荧光素酶报告基因试验验证ac115与ame-miR-13b的结合关系。通过饲喂特异性siRNA对蜜蜂球囊菌胁迫的意蜂幼虫肠道内ac115进行干扰,进而检测干扰ac115对宿主免疫应答相关的6个基因表达量的影响。结果表明,ac115的BS位点真实存在;相较于未接种组,蜜蜂球囊菌接种组4日龄幼虫肠道内ac115的表达量极显著上调(P<0.000 1),5和6日龄幼虫肠道内ac115显著上调表达(P<0.05)...