【目的】本研究结合前期已获得的miRNA和mRNA组学数据对东方蜜蜂微孢子虫Nosema ceranae的差异表达miRNA(differentially expressed miRNA,DEmiRNA)靶向意大利蜜蜂Apis mellifera ligustica工蜂中肠的mRNA和差异表达mRNA(differentially...【目的】本研究结合前期已获得的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跨界调控意大利蜜蜂工蜂的能量代谢以加强能量窃取并促进增殖。展开更多
【背景】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种专性侵染蜜蜂幼虫的真菌病原,可引起成年蜜蜂数量和蜂群群势的急剧下降。长链非编码RNA(long non-coding RNA,lncRNA)是一类新近发现的非编码RNA,在表观遗传、细胞周期、剂量补...【背景】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种专性侵染蜜蜂幼虫的真菌病原,可引起成年蜜蜂数量和蜂群群势的急剧下降。长链非编码RNA(long non-coding RNA,lncRNA)是一类新近发现的非编码RNA,在表观遗传、细胞周期、剂量补偿等众多生命活动中发挥重要生物学功能。【目的】明确球囊菌菌丝和孢子中lncRNA的数量、种类和表达谱差异,并探究共有lncRNA、特有lncRNA和差异表达lncRNA(differentially expressed lncRNA,DElncRNA)在菌丝与孢子中的潜在功能。【方法】利用基于链特异性建库的lncRNA-seq技术对球囊菌的纯化菌丝(AaM)和纯化孢子(AaS)分别进行测序。根据FPKM(Fragment Per Kilobase of per Million mapped reads)法计算lncRNA在AaM和AaS中的表达水平。通过Venn分析筛选AaM与AaS的共有lncRNA和特有lncRNA。按照P≤0.05且|log2 fold change|≥1的标准筛选AaM vs AaS比较组中的DElncRNA。通过Blast工具将共有lncRNA、特有lncRNA和DElncRNA的上下游基因比对GO和KEGG数据库,以进行功能及通路注释。根据靶向结合关系构建共有lncRNA、特有lncRNA和DElncRNA的竞争性内源RNA(competing endogenous RNA,ceRNA)调控网络并利用Cytoscape软件进行可视化。利用RT-qPCR验证测序数据的可靠性。【结果】AaM和AaS分别测得108614646和105675408条原始读段(raw reads),经严格过滤得到107780032和104621402条有效读段(clean reads),Q20分别为98.76%和98.72%,Q30分别达到95.84%和95.78%。共鉴定到850个lncRNA。AaM和AaS的共有lncRNA为701个,二者的特有lncRNA分别为39和110个。上述共有lncRNA通过顺式作用调控3992个上下游基因,它们涉及细胞进程、代谢进程和催化活性等42个功能条目,以及代谢途径、次生代谢物的生物合成和抗生素的生物合成等117条通路;AaM的特有lncRNA和AaS的特有lncRNA通过顺式作用分别调控243和672个上下游基因。AaM vs AaS比较组包含的255个DElncRNA通过顺式作用调控1479个上下游基因,它们涉及代谢进程、细胞进程和催化活性等41个功能条目,以及代谢途径、次生代谢产物的生物合成和抗生素的生物合成等107条通路。从共有lncRNA、孢子特有lncRNA和DElncRNA中分别预测到41、5和13个微小RNA(microRNA,miRNA)的前体序列。调控网络分析结果显示,菌丝lncRNA、孢子lncRNA形成较为复杂的ceRNA调控网络;菌丝lncRNA可靶向结合8个miRNA,进而调控77个mRNA;孢子lncRNA可靶向结合7个miRNA,进而调控87个mRNA;2个DElncRNA(TCONS00008630与TCONS00009302)可靶向结合miR-4968-y,进而调控10个mRNA。RT-qPCR验证结果显示4个DElncRNA的差异表达趋势与测序结果一致,表明测序数据真实可靠。【结论】共有lncRNA、特有lncRNA和DElncRNA可能通过调控上下游基因的表达,作为miRNA的前体,以及充当ceRNA影响菌丝和孢子的物质和能量代谢、自噬、转录、MAPK信号通路、泛素介导的蛋白水解、蛋白酶体以及次生代谢产物的生物合成等生物学过程,从而调节球囊菌的生长、发育、生殖和致病性。展开更多
【目的】东方蜜蜂微孢子虫(Nosema ceranae)感染意大利蜜蜂(Apis mellifera ligustica,简称意蜂)导致蜜蜂微孢子虫病。本研究结合前期已获得的miRNA和mRNA组学数据,通过生物信息学方法对意蜂工蜂中肠的差异表达miRNA(differentially exp...【目的】东方蜜蜂微孢子虫(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分别靶向病原的1345和1046条mRNA;进一步分析发现,宿主的47条显著上调miRNA和34条显著下调miRNA可分别靶向NcCK vs NcT1比较组中病原的584条显著下调mRNA和265条显著上调mRNA,它们可分别注释到19和22个功能条目以及66和64条通路。AmCK2 vs AmT2比较组中宿主的56条显著上调miRNA和51条显著下调miRNA分别靶向病原的1260和1317条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影响病原的侵染和增殖。展开更多
【目的】本研究旨在通过差异表达基因(differentially expressed gene,DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于...【目的】本研究旨在通过差异表达基因(differentially expressed gene,DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log 2(Fold change)|≥1的标准,通过比较分析筛选出NcCK vs NcT1,NcCK vs NcT2和NcT1 vs NcT2比较组的DEG。通过相关生物信息学软件对上述DEG进行Venn分析、GO分类和KEGG代谢通路富集分析。根据Nr和KEGG数据库注释信息和相关文献进行对东方蜜蜂微孢子虫的毒力因子和侵染相关因子的统计和分析。通过RT-qPCR验证转录组数据及DEG表达趋势。【结果】从NcCK vs NcT1,NcCK vs NcT2和NcT1 vs NcT2比较组分别鉴定出1397,1497和52个DEG。Venn分析结果显示各比较组共有的上调和下调基因分别为10和1个。GO分类结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG富集数最多的功能条目为代谢进程、细胞进程、单组织进程、细胞、细胞组件、细胞器、催化活性和结合,而NcT1 vs NcT2中DEG富集数最多的是代谢进程、细胞进程、单组织进程、催化活性和结合。KEGG代谢通路富集分析结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG分别富集到80和79条通路;富集在糖酵解/糖异生和MAPK信号通路的上调基因数量多于下调基因。毒力因子分析结果显示,孢壁蛋白9基因和孢壁蛋白12基因在NcCK vs NcT1和NcCK vs NcT2中均下调表达,孢壁蛋白8基因仅在NcCK vs NcT1中表达量下调;此外孢壁蛋白前体基因、孢壁和锚定盘复合蛋白基因、几丁质合酶基因、极管蛋白基因、蓖麻毒素B凝集素基因的表达水平在NcCK vs NcT1和NcCK vs NcT2中表现为上调。侵染相关因子分析结果表明,糖酵解途径的3个关键酶基因在NcCK vs NcT1和NcCK vs NcT2中上调表达;3个涉及ATP/ADP移位酶的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,但有1个表达量下调;2个涉及ABC转运蛋白的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,另有4个下调表达。RT-qPCR结果证实了本研究中转录组数据及DEG表达趋势的真实可靠性。【结论】本研究通过比较分析解析东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程的转录组动态,揭示了孢壁蛋白、孢壁和锚定盘复合蛋白、几丁质酶、极管蛋白和蓖麻毒素B凝集素等毒力因子编码基因,以及己糖激酶、丙酮酸激酶、6-磷酸果糖激酶、ATP/ADP移位酶和ABC转运蛋白等侵染相关因子编码基因在病原增殖中扮演重要角色,为阐明东方蜜蜂微孢子虫的侵染机制提供了基础。展开更多
【目的】通过对东方蜜蜂微孢子虫(Nosema ceranae)纯化孢子与侵染意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂的东方蜜蜂微孢子虫的差异表达miRNA(DEmiRNA)及其靶mRNA进行系统分析,筛选、分析和探讨病原毒力因子和侵染因子相关...【目的】通过对东方蜜蜂微孢子虫(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,|log2 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分别预测出1885、1733和1524个靶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等侵染因子基因表达进行调控,从而适应宿主细胞内的环境并促进自身的增殖与侵染。展开更多
文摘【目的】本研究结合前期已获得的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跨界调控意大利蜜蜂工蜂的能量代谢以加强能量窃取并促进增殖。
文摘【背景】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种专性侵染蜜蜂幼虫的真菌病原,可引起成年蜜蜂数量和蜂群群势的急剧下降。长链非编码RNA(long non-coding RNA,lncRNA)是一类新近发现的非编码RNA,在表观遗传、细胞周期、剂量补偿等众多生命活动中发挥重要生物学功能。【目的】明确球囊菌菌丝和孢子中lncRNA的数量、种类和表达谱差异,并探究共有lncRNA、特有lncRNA和差异表达lncRNA(differentially expressed lncRNA,DElncRNA)在菌丝与孢子中的潜在功能。【方法】利用基于链特异性建库的lncRNA-seq技术对球囊菌的纯化菌丝(AaM)和纯化孢子(AaS)分别进行测序。根据FPKM(Fragment Per Kilobase of per Million mapped reads)法计算lncRNA在AaM和AaS中的表达水平。通过Venn分析筛选AaM与AaS的共有lncRNA和特有lncRNA。按照P≤0.05且|log2 fold change|≥1的标准筛选AaM vs AaS比较组中的DElncRNA。通过Blast工具将共有lncRNA、特有lncRNA和DElncRNA的上下游基因比对GO和KEGG数据库,以进行功能及通路注释。根据靶向结合关系构建共有lncRNA、特有lncRNA和DElncRNA的竞争性内源RNA(competing endogenous RNA,ceRNA)调控网络并利用Cytoscape软件进行可视化。利用RT-qPCR验证测序数据的可靠性。【结果】AaM和AaS分别测得108614646和105675408条原始读段(raw reads),经严格过滤得到107780032和104621402条有效读段(clean reads),Q20分别为98.76%和98.72%,Q30分别达到95.84%和95.78%。共鉴定到850个lncRNA。AaM和AaS的共有lncRNA为701个,二者的特有lncRNA分别为39和110个。上述共有lncRNA通过顺式作用调控3992个上下游基因,它们涉及细胞进程、代谢进程和催化活性等42个功能条目,以及代谢途径、次生代谢物的生物合成和抗生素的生物合成等117条通路;AaM的特有lncRNA和AaS的特有lncRNA通过顺式作用分别调控243和672个上下游基因。AaM vs AaS比较组包含的255个DElncRNA通过顺式作用调控1479个上下游基因,它们涉及代谢进程、细胞进程和催化活性等41个功能条目,以及代谢途径、次生代谢产物的生物合成和抗生素的生物合成等107条通路。从共有lncRNA、孢子特有lncRNA和DElncRNA中分别预测到41、5和13个微小RNA(microRNA,miRNA)的前体序列。调控网络分析结果显示,菌丝lncRNA、孢子lncRNA形成较为复杂的ceRNA调控网络;菌丝lncRNA可靶向结合8个miRNA,进而调控77个mRNA;孢子lncRNA可靶向结合7个miRNA,进而调控87个mRNA;2个DElncRNA(TCONS00008630与TCONS00009302)可靶向结合miR-4968-y,进而调控10个mRNA。RT-qPCR验证结果显示4个DElncRNA的差异表达趋势与测序结果一致,表明测序数据真实可靠。【结论】共有lncRNA、特有lncRNA和DElncRNA可能通过调控上下游基因的表达,作为miRNA的前体,以及充当ceRNA影响菌丝和孢子的物质和能量代谢、自噬、转录、MAPK信号通路、泛素介导的蛋白水解、蛋白酶体以及次生代谢产物的生物合成等生物学过程,从而调节球囊菌的生长、发育、生殖和致病性。
文摘【目的】东方蜜蜂微孢子虫(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分别靶向病原的1345和1046条mRNA;进一步分析发现,宿主的47条显著上调miRNA和34条显著下调miRNA可分别靶向NcCK vs NcT1比较组中病原的584条显著下调mRNA和265条显著上调mRNA,它们可分别注释到19和22个功能条目以及66和64条通路。AmCK2 vs AmT2比较组中宿主的56条显著上调miRNA和51条显著下调miRNA分别靶向病原的1260和1317条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影响病原的侵染和增殖。
文摘【目的】本研究旨在通过差异表达基因(differentially expressed gene,DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log 2(Fold change)|≥1的标准,通过比较分析筛选出NcCK vs NcT1,NcCK vs NcT2和NcT1 vs NcT2比较组的DEG。通过相关生物信息学软件对上述DEG进行Venn分析、GO分类和KEGG代谢通路富集分析。根据Nr和KEGG数据库注释信息和相关文献进行对东方蜜蜂微孢子虫的毒力因子和侵染相关因子的统计和分析。通过RT-qPCR验证转录组数据及DEG表达趋势。【结果】从NcCK vs NcT1,NcCK vs NcT2和NcT1 vs NcT2比较组分别鉴定出1397,1497和52个DEG。Venn分析结果显示各比较组共有的上调和下调基因分别为10和1个。GO分类结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG富集数最多的功能条目为代谢进程、细胞进程、单组织进程、细胞、细胞组件、细胞器、催化活性和结合,而NcT1 vs NcT2中DEG富集数最多的是代谢进程、细胞进程、单组织进程、催化活性和结合。KEGG代谢通路富集分析结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG分别富集到80和79条通路;富集在糖酵解/糖异生和MAPK信号通路的上调基因数量多于下调基因。毒力因子分析结果显示,孢壁蛋白9基因和孢壁蛋白12基因在NcCK vs NcT1和NcCK vs NcT2中均下调表达,孢壁蛋白8基因仅在NcCK vs NcT1中表达量下调;此外孢壁蛋白前体基因、孢壁和锚定盘复合蛋白基因、几丁质合酶基因、极管蛋白基因、蓖麻毒素B凝集素基因的表达水平在NcCK vs NcT1和NcCK vs NcT2中表现为上调。侵染相关因子分析结果表明,糖酵解途径的3个关键酶基因在NcCK vs NcT1和NcCK vs NcT2中上调表达;3个涉及ATP/ADP移位酶的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,但有1个表达量下调;2个涉及ABC转运蛋白的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,另有4个下调表达。RT-qPCR结果证实了本研究中转录组数据及DEG表达趋势的真实可靠性。【结论】本研究通过比较分析解析东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程的转录组动态,揭示了孢壁蛋白、孢壁和锚定盘复合蛋白、几丁质酶、极管蛋白和蓖麻毒素B凝集素等毒力因子编码基因,以及己糖激酶、丙酮酸激酶、6-磷酸果糖激酶、ATP/ADP移位酶和ABC转运蛋白等侵染相关因子编码基因在病原增殖中扮演重要角色,为阐明东方蜜蜂微孢子虫的侵染机制提供了基础。
文摘【目的】通过对东方蜜蜂微孢子虫(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,|log2 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分别预测出1885、1733和1524个靶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等侵染因子基因表达进行调控,从而适应宿主细胞内的环境并促进自身的增殖与侵染。