侵染中华蜜蜂6日龄幼虫的蜜蜂球囊菌的微小RNA差异表达谱及调控网络

【目的】蜜蜂球囊菌(Ascosphaeraapis,简称球囊菌)是专性侵染蜜蜂幼虫的致死性真菌病原。MicroRNA(miRNA)作为一类重要的基因表达调控因子,能够广泛参与真菌及其宿主的相互作用过程。本研究通过比较分析球囊菌孢子(AaCK)和侵染中华蜜蜂(Apis cerana cerana,简称中蜂) 6日龄幼虫肠道内的球囊菌(AaT)的smallRNA(sRNA)组学数据对球囊菌的差异表达miRNA(differentiallyexpressed miRNA,DEmiRNA)、靶mRNA及二者间的调控网络进行全面解析,旨在揭示miRNA介导的球囊菌对中蜂幼虫的侵染机制。【方法】对于球囊菌侵染的中蜂6日龄幼虫肠道的small RNA-seq (sRNA-seq)数据,利用BLAST工具连续比对东方蜜蜂(Apiscerana)和球囊菌的参考基因组筛滤得到AaT的sRNA组学数据。分别将AaCK和AaT的sRNA组学数据比对miRBase数据库,对球囊菌侵染宿主前后miRNA的数量和结构特征进行分析。联用RNAhybrid+svm_light、Miranda和TargetScan软件预测AaCK vs AaT比较组中DEmiRNA的靶mRNA,进而利用相关生物信息学软件对上述靶mRNA进行GO分类和KEGG代谢通路富集分析。通过Cytoscape软件对DEmiRNA-mRNA调控网络进行可视化。利用Stem-loop RT-PCR、RT-qPCR和分子克隆验证测序结果的可靠性。【结果】在AaCK和AaT中分别鉴定到380和387个miRNA。结构特征分析结果显示,AaCK和AaT的mi RNA皆集中分布在18–25 nt,且首位碱基主要偏向于U。AaCKvsAaT比较组共有270个DEmiRNA,包含155个上调miRNA和115个下调miRNA,分别靶向结合6091和6145个mRNA。GO分类结果显示,上述靶mRNA主要涉及代谢进程、细胞进程、应激反应等15个生物学进程;细胞、细胞组分、细胞器等12个细胞组分;催化活性、结合、转运子活性等11个分子功能。KEGG代谢通路富集分析结果显示,上述靶mRNA富集在123条代谢通路,参与对氨基酸代谢、碳水化合物代谢以及核苷酸代谢等物质代谢,氧化磷酸化、硫代谢、氮代谢等能量代谢,以及MAPK和Hippo等信号通路的调控。球囊菌DEmiRNA与靶mRNA之间存在复杂的调控关系,其中miR-29-x、miR-250-x、miR-4968-y、miR-11200-x、novel-m0023-5p、novel-m0130-5p和novel-m0135-5p等DEmiRNA可靶向结合与球囊菌的半胱氨酸蛋白酶、DNA甲基化转移酶以及几丁质酶相关的mRNA;此外,miR-7-x、miR-9-z、miR-319-y和miR-5951-y等同时参与调控MAPK信号通路;进一步分析发现,miR-250-x同时参与对DNA甲基化转移酶、MAPK信号通路及其他酶类合成与代谢途径的调控,并可能参与球囊菌与中蜂6日龄幼虫之间的跨界调控。通过Stem-loopRT-PCR和RT-qPCR验证了4个DEmiRNA的差异表达,并利用分子克隆和Sanger测序证实miR-7-x的序列与测序结果一致。【结论】本研究解析了侵染中蜂6日龄幼虫的球囊菌的miRNA差异表达谱及DEmiRNA的调控网络,揭示了球囊菌DEmiRNA可能通过调控病原的物质和能量代谢、增殖、毒力、信号通路及相关mRNA参与对中蜂幼虫的侵染过程。miR-7-x、miR-250-x、novel-m0023-5p等关键DEmiRNA有望作为白垩病治疗的新型分子靶点。     

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

【目的】本研究结合前期已获得的miRNA和mRNA组学数据对东方蜜蜂微孢子虫Nosema ceranae的差异表达miRNA(differentially expressed miRNA, DEmiRNA)靶向意大利蜜蜂Apis mellifera ligustica工蜂中肠的mRNA和差异表达mRNA(differentially expressed mRNA, DEmRNA)进行生物信息学预测、数据库注释和调控网络分析,旨在解析东方蜜蜂微孢子虫对意大利蜜蜂工蜂中肠的基因表达调控。【方法】比较东方蜜蜂微孢子虫感染7 d(AmT1)和10 d(AmT2)的意大利蜜蜂工蜂中肠和未感染中肠(分别为AmCK1和AmCK2)的mRNA数据,筛选意大利蜜蜂工蜂的显著性DEmRNA;通过比较侵染AmT1和AmT2的东方蜜蜂微孢子虫(分别为NcT1和NcT2)和东方蜜蜂微孢子虫纯净孢子(NcCK)的miRNA数据筛选出东方蜜蜂微孢子虫的DEmiRNA。利用TargetFinder软件预测东方蜜蜂微孢子虫DEmiRNA靶向的意大利蜜蜂工蜂中肠DEmRNA。利用相关生物信息学工具对上述意大利蜜蜂工蜂中肠靶DEmRNA进行GO和KEGG数据库注释。根据KEGG数据库注释信息筛选出意大利蜜蜂工蜂中肠免疫防御和能量代谢通路相关DEmRNA,并构建和分析上述DEmRNA与相应的东方蜜蜂微孢子虫DEmiRNA之间的调控网络。【结果】NcCK vs NcT1比较组中东方蜜蜂微孢子虫的77条显著上调miRNA和52条显著下调miRNA可分别靶向AmCK1 vs AmT1比较组中意大利蜜蜂工蜂中肠的118条显著下调mRNA和135条显著上调mRNA,这些mRNA可分别注释到31和25个GO条目,以及113和107条KEGG通路。NcCK vs NcT2比较组中东方蜜蜂微孢子虫的52条显著上调miRNA和49条显著下调miRNA可分别靶向AmCK2 vs AmT2比较组中意大利蜜蜂工蜂中肠的97条显著下调mRNA和210条显著上调mRNA,这些mRNA可分别注释到27和30个GO条目以及97和127条KEGG通路。NcCK vs NcT1和NcCK vs NcT2比较组中的11条共同显著上调miRNA和19条共同显著下调miRNA分别靶向AmCK1 vs AmT1和AmCK2 vs AmT2比较组的6条共同显著下调和14条共同显著上调mRNA,可分别注释到7和10个GO条目以及0和9条KEGG通路。NcCK vs NcT1和NcCK vs NcT2比较组中东方蜜蜂微孢子虫的DEmiRNA可靶向AmCK1 vs AmT1和AmCK2 vs AmT2比较组中意大利蜜蜂工蜂中肠的氧化磷酸化和硫代谢等能量代谢通路相关DEmRNA,以及胞吞作用、黑色素生成、溶酶体、自噬、Toll样受体信号通路、细胞凋亡、Ras信号通路、泛素介导的蛋白水解和MAPK信号通路等免疫防御通路相关DEmRNA。进一步分析发现,miR-216-x, miR-5119-y, bantam-y和miR-8-y在NcCK vs NcT1和NcCK vs NcT2比较组中皆显著上调表达,且靶向意大利蜜蜂工蜂中肠的溶酶体、黑色素生成、泛素介导的蛋白水解、MAPK信号通路及Ras信号通路等免疫防御通路相关的显著下调mRNA。【结论】在侵染过程中,东方蜜蜂微孢子虫的DEmiRNA对意大利蜜蜂工蜂的基因表达具有广泛的潜在影响,东方蜜蜂微孢子虫可能通过上调部分miRNA跨界调控意大利蜜蜂工蜂的免疫防御以促进侵染,通过下调部分miRNA跨界调控意大利蜜蜂工蜂的能量代谢以加强能量窃取并促进增殖。     

互作基因定位方法的研究和计算机软件的开发

质量性状中存在6种可能的基因互作类型, 即互补、重叠、累加、显性上位、隐性上位和抑制。在遗传研究中, 有时会遇到互作基因定位的问题, 但至今未见有关互作基因定位方法和计算机软件的系统的研究报道。文章给出了基于极大似然估计的互作基因定位方法及相应的计算机软件(IGMapping 1.0)。计算机模拟表明, 此方法可以无偏地估计一个共显性标记与一个互作基因之间的重组率或连锁距离。     

东方蜜蜂微孢子虫的基因结构优化及新基因鉴定

东方蜜蜂微孢子虫Nosema ceranae是一种寄生于蜜蜂中肠上皮细胞的单细胞真菌,对蜜蜂的健康危害严重,给世界各国的养蜂业造成较大损失。本研究基于前期获得的N.ceranae孢子的转录组数据对其已注释基因进行结构优化,并对未注释基因进行预测和分析。通过将测序得到的clean reads比对参考基因组和转录本重构,共对10个N.ceranae的已注释基因的5'端或3'端进行了延长。利用Cuffcompare软件将重构转录本与参考基因组进行比对,共鉴定出27个新基因,随机挑选9个新基因进行RT-PCR验证,均能扩增出符合预期的目的片段,表明预测出的新基因真实存在。有6个新基因能够注释到GO数据库和6个基因注释到KEGG数据库。进一步分析结果显示上述新基因注释到细胞等10个GO条目上,它们可能在N.ceranae的生命活动中具有重要功能。研究结果为N.ceranae的基因结构和功能注释信息的完善提供了有益补充,也为新基因的功能研究打下了基础。     

施氮对水稻土N2O释放及反硝化功能基因(narG/nosZ)丰度的影响

以紫潮泥和红黄泥两种不同质地的水稻土壤作为研究对象,通过室内培养试验,分析施用硝态氮肥对N2O释放和反硝化基因(narG/nosZ)丰度的影响,并探讨反硝化基因丰度与N2O释放之间的关系。结果表明,施用硝态氮显著增加两种水稻土的N2O释放量。在72h培养过程中,施氮改变了紫潮泥反硝化基因(narG/nosZ)的丰度,但并未明显影响红黄泥反硝化基因(narG/nosZ)丰度。通过双变量相关分析发现,除了紫潮泥narG基因外,其它的反硝化基因丰度和N2O释放之间并没有显著相关性。     

脑心通胶囊辅助治疗不稳定型心绞痛的临床疗效和安全性的Meta 分析

摘要目的：系统评价脑心通胶囊辅助治疗不稳定型心绞痛的临床疗效及其安全性。方法：计算机检索PubMed(1966～2012.3)、 EMbase(1974～2012.3)、Cochrane Library(2012 年第3 期)、CBM(1978～2012.3)、CNKI(1994～2012.3)、VIP(1989～2012.3)。纳入脑 心通胶囊联合常规药物治疗不稳定型心绞痛的随机对照试验(RCT),并补充检索纳入研究的参考文献;按Cochrane 系统评价方法 对纳入研究进行资料提取及质量评估后,采用RevMan 5.1 软件进行统计学分析。结果：共纳入11 个RCT,包括共1384 例患者, Meta 分析结果显示：与常规药物治疗相比,脑心通胶囊联合常规药物明显缓解不稳定型心绞痛患者的临床症状,差异有统计学意 义[RR=1.24,95%CI(1.18-1.31),P<0.01],提高了心电图的改善率[RR=1.35,95%CI(1.24-1.47),P<0.01]。结论：目前研究表明脑心通 胶囊辅助治疗不稳定型心绞痛的短期疗效明显优于常规治疗,但其不良反应和远期疗效尚不确定,需要更多高质量随机对照试 验进一步证实。     

意大利蜜蜂工蜂中肠响应Nosema ceranae胁迫的高表达基因分析

东方蜜蜂微孢子虫Nosema ceranae是专性寄生蜜蜂中肠上皮细胞的单细胞真菌,对意大利蜜蜂(意蜂)具有较强的侵染性。本研究利用RNA-seq技术对正常10 d (Apis mellifera ligustica control group,Am CK)和N. ceranae胁迫10 d的意蜂工蜂中肠(Apis mellifera ligustica treatment group,Am T)进行测序,共得到160 847 237 100条原始读段,过滤后得到1 066 955 298条有效读段。主成分分析结果显示Am CK与Am T测序样品的组内重复性较好,组间的基因表达模式差异明显。GO分类结果显示,Am CK与Am T的前100位高表达基因(HEGs)分别分布于32和33个GO terms,基因富集数最多的均为代谢进程。KEGG代谢通路(pathway)富集分析结果显示,Am CK的前100位HEGs富集在21个pathways,基因富集数最多的是内吞作用、信号通路和嘌呤代谢; Am T的前100位HEGs富集在26个pathways,基因富集数最多的是内吞作用、RNA转运和蛋白质的内质网加工。前100位HEGs的Venn分析结果显示Am CK与Am T的共有HEGs为87个,特有HEGs分别均为13个。研究结果揭示了N. ceranae胁迫意蜂工蜂中肠过程中的基因表达谱信息,也为解析N. ceranae致病的分子机理提供了有益信息和线索。     

意大利蜜蜂工蜂中肠的环状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在意大利蜜蜂中肠发育及胁迫响应过程中的功能奠定了基础。     

蜜蜂球囊菌的microRNA鉴定及其调控网络分析

【目的】本研究利用small RNA-seq技术对球囊菌的纯培养进行测序,对球囊菌的micro RNAs miRNAs)进行预测、鉴定和分析,进而构建miRNAs-mRNAs的调控网络。【方法】利用Illumina Hiseq Xten平台对球囊菌菌丝与孢子进行测序,通过相关生物信息学软件对球囊菌的miRNAs进行预测和分析,通过茎环(Stem-loop)PCR对部分miRNAs进行鉴定,利用Cytoskype软件构建miRNAs-mRNAs的调控网络。【结果】本研究共获得48268696条clean reads,预测出118个球囊菌的miRNAs,它们的长度分布介于18–25 nt之间,不同长度的mi RNA的首位碱基偏好性差异明显。Stem-loop PCR验证结果显示共有10个miRNAs能够扩增出符合预期的目的片段,说明多数miRNAs可能真实存在。共预测出6529个球囊菌miRNAs的靶基因,其中5725个能够注释到Nr、Swissprot、KOG、GO和KEGG数据库。进一步分析结果显示有24个靶基因注释在MAPK信号通路。Cytoskype软件分析结果显示球囊菌的miRNAs与mRNAs之间存在复杂的调控网络,绝大多数的miRNAs处于调控网络的内部且同时结合多个mRNAs。【结论】本研究率先对球囊菌的miRNAs及miRNAs-mRNAs调控网络进行全面分析,研究结果丰富了对球囊菌miRNAs的认识,为其基础生物学信息提供了有益补充,也为阐明球囊菌致病的分子机理打下了一定基础。     

意大利蜜蜂幼虫肠道在球囊菌胁迫后期的差异表达微小RNA及其靶基因分析

【目的】蜜蜂球囊菌Ascosphaera apis(简称球囊菌)是一种特异性侵染蜜蜂幼虫肠道的致死性真菌病原。微小RNA(microRNA,mi RNA)是一种在转录后水平对m RNA进行负调控的关键调控因子。本研究旨在对意大利蜜蜂Apis mellifera ligustica(简称意蜂)幼虫肠道在球囊菌胁迫后期的差异表达mi RNAs(DEmiRNAs)及其靶基因进行深入分析,为揭示DEmiRNAs在胁迫应答中的作用提供重要信息。【方法】利用small RNA-seq(s RNA-seq)技术对正常及球囊菌侵染的意蜂6日龄幼虫肠道(分别表示为Am CK和Am T)进行测序,通过相关生物信息学软件对DEmiRNAs进行预测和分析。利用TargetFinder软件预测DEmiRNAs的靶基因,然后利用Blast软件对靶基因进行GO和KEGG数据库的功能注释。通过Cytoscape软件构建DEmiRNAs与其靶m RNAs的调控网络。通过茎环实时荧光定量PCR(stem-loop RT-qPCR)验证测序数据的可靠性。【结果】Am CK和Am T的测序分别得到9 230 496和10 823 667条有效序列标签; Am CK vs Am T比较组中包含15个上调和6个下调mi RNAs,分别结合3 503和3 252个靶基因,它们可分别注释到40和39个GO terms以及104和99条代谢通路;进一步分析发现调控网络中的17个DEmiRNAs靶向结合116个与丝氨酸蛋白酶相关的m RNAs,14个DEmiRNAs靶向结合54个与泛素介导的蛋白水解相关的m RNAs。Stem-loop RT-qPCR验证结果显示随机选择的4个DEmiRNAs (mi R-251-x,mi R-9277-y,mi R-1672-x和mi R-4968-y)的表达量变化趋势与测序结果一致,表明测序数据真实可信。【结论】本研究率先对意蜂幼虫肠道在球囊菌胁迫后期的DEmiRNAs及其靶基因进行预测与分析,提供了mi RNAs的表达谱和差异表达信息。ame-miR-927b,mi R-429-y和mi R-8440-y等DEmiRNAs可能参与对丝氨酸蛋白酶的调控和对泛素介导的蛋白水解的调控,DEmiRNAs与靶基因之间存在复杂的调控网络关系。DEmiRNAs可参与到意蜂幼虫肠道与球囊菌间的互作。