【目的】本研究旨在通过差异表达基因(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等侵染因子基因表达进行调控,从而适应宿主细胞内的环境并促进自身的增殖与侵染。展开更多
文摘【目的】本研究旨在通过差异表达基因(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等侵染因子基因表达进行调控,从而适应宿主细胞内的环境并促进自身的增殖与侵染。
