意大利蜜蜂幼虫肠道在球囊菌胁迫后期的差异表达微小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可参与到意蜂幼虫肠道与球囊菌间的互作。     

中华蜜蜂6日龄幼虫响应蜜蜂球囊菌胁迫的环状RNA应答

【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)专性侵染蜜蜂幼虫而导致白垩病,危害蜜蜂健康和养蜂生产。本研究旨在探究中华蜜蜂(Apis cerana cerana,简称中蜂) 6日龄幼虫响应球囊菌胁迫的环状RNA(circular RNA,circRNA)差异表达谱及差异表达circRNA (differentially expressed circRNA,DEcircRNA)在宿主胁迫应答中的潜在功能。【方法】利用去除线性RNA的circRNA-seq技术对正常和球囊菌侵染的中蜂6日龄幼虫肠道(AcCK和AcT)进行测序。利用find_circ软件鉴定circRNA,统计circRNA的长度和环化类型。根据|log_2(Fold change)|≥1和P≤0.05的标准筛选DEcircRNA。将DEcircRNA的来源基因比对Gene ontology (GO)数据库和Kyoto Encyclopedia of Genes and Genomes (KEGG)数据库,从而获得功能及通路(pathway)注释。随机挑选3个DEcircRNA进行RT-qPCR验证。【结果】AcCK和AcT的circRNA-seq分别得到76342570和68269362条原始读段(raw reads),经严格质控得到74524108和66974392条有效读段(clean reads),Q30分别为92.75%和94%,GC含量分别为54.31%和54.90%。比对上东方蜜蜂(Apis cerana)参考基因组的短序列读段(anchor reads)共计23648400条。AcCK和AcT中分别鉴定到805和702个circRNA,长度均介于201–1000 nt,数量最多的环化类型均为已注释外显子circRNA,但分布在不同长度、不同环化类型的circRNA数量存在差异。AcCK vs AcT比较组共有494个DEcircRNA,包括257个上调circRNA和237个下调circRNA;上调和下调幅度最大的circRNA分别为novel_circ_000123和novel_circ_000726。上述DEcircRNA的来源基因可注释到11条生物学进程相关条目,9条分子功能相关条目,9条细胞组分相关条目,以及73条通路。进一步分析发现,部分DEcircRNA的来源基因注释到7条细胞免疫通路和3条体液免疫通路。【结论】中蜂6日龄幼虫响应球囊菌胁迫的过程中可能通过改变分布在不同长度和环化类型的circRNA数量,以及特异性表达一些circRNA和调节部分circRNA的表达量对病原产生应答;novel_circ_000027、novel_circ_000127、novel_circ_000312等DEcircRNA在宿主的胁迫应答过程中可能通过调控氧化磷酸化、细胞和体液免疫等通路发挥特殊作用。研究结果为深入理解中蜂幼虫对球囊菌的胁迫应答机制及二者的相互作用机制提供了新见解。     

LncRNA13164通过ace-miR-4968-y调节中华蜜蜂幼虫对蜜蜂球囊菌侵染的免疫应答

【目的】长链非编码RNA (long non-coding RNA,lncRNA)可通过竞争性内源RNA(competing endogenous RNA,ce RNA)等多种方式发挥重要的调控作用,但蜜蜂lncRNA的功能研究至今依然缺失,lncRNA在蜜蜂免疫应答中的作用尚不清楚。本研究旨在揭示中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫响应蜜蜂球囊菌(Ascosphaera apis)侵染的免疫应答中lncRNA的调控功能及作用机制。笔者团队前期通过深度测序和生物信息学分析发现,lncRNA13164与ace-miR-4968存在靶向关系且潜在参与在中蜂幼虫对蜜蜂球囊菌侵染的应答。【方法】利用实时荧光定量PCR (real-time fluorescence quantitative polymerase chain reaction,RT-qPCR)检测蜜蜂球囊菌接种后中蜂幼虫肠道内lncRNA13164的表达谱。采用ILoc-Lnc RNA软件预测lncRNA13164的亚细胞定位。联用RNAhybrid、Miranda和Target Scan软件预测lnc...     

中华蜜蜂6日龄幼虫响应蜜蜂球囊菌侵染的长链非编码RNA应答研究

[目的]本研究旨在探究长链非编码RNA (long non-coding RNA,lncRNA)在中华蜜蜂(Apis cerana cerana,简称中蜂)6日龄幼虫应答蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)侵染过程中的差异表达谱及调控作用.[方法]利用链特异性cDNA建库的RNA-seq技术对未被...     

过表达和敲减ame-miR-13b对意大利蜜蜂幼虫肠道内基因表达的影响

【目的】微小RNA(microRNAs, miRNAs)在昆虫的繁殖、发育和免疫等重要生命活动的调控过程中扮演关键角色。本研究旨在探究在意大利蜜蜂Apis mellifera ligustica幼虫个体水平过表达与敲减ame-miR-13b对其靶基因表达的影响,为深入探究ame-miR-13b调控意大利蜜蜂幼虫肠道发育的分子机理提供理论和实验依据。【方法】根据ame-miR-13b的核苷酸序列设计合成相应的模拟物mimic-ame-miR-13b和抑制物inhibitor-ame-miR-13b及其对照mimic-NC和inhibitor-NC,混入饲料后饲喂意大利蜜蜂3日龄幼虫,每12 h更换一次饲料,连续饲喂6次,以分别过表达和敲减ame-miR-13b;通过RT-qPCR检测过表达和敲减ame-miR-13b后意大利蜜蜂幼虫肠道中ame-miR-13b及其靶基因Ecr, Egfr和P450 18a1的表达量。【结果】与mimic-NC饲喂组相比, mimic-ame-miR-13b饲喂组中ame-miR-13b在意大利蜜蜂4和6日龄幼虫肠道中均显著上调表达;与饲喂inhibitor-NC组相比,inhibitor-ame-miR-13b饲喂组中ame-miR-13b在4日龄幼虫肠道中下调表达,在5和6日龄幼虫肠道中显著下调表达。与mimic-NC饲喂组相比,过表达ame-miR-13b组中意大利蜜蜂6日龄幼虫肠道中Egfr的表达量显著降低,而Ecr和P450 18a1的表达量均显著升高。与inhibitor-NC饲喂组相比,敲减ame-miR-13b组中意大利蜜蜂6日龄幼虫肠道中Egfr上调表达但不显著,Ecr显著上调表达,而P450 18a1显著下调表达。【结论】通过饲喂法在意大利蜜蜂幼虫个体中成功实现了miRNA的过表达和敲减;ame-miR-13b与Egfr之间存在潜在的负调控关系。     

蜜蜂球囊菌菌丝和孢子中环状RNA的鉴定及比较分析

[目的]本研究旨在明确蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)菌丝和孢子中环状RNA(circular RNA,circRNA)的数量、种类和表达谱差异,并探讨共有circRNA、特有circRNA和差异表达circRNA (differentially expressed circRNA,DEcir...     

球囊菌侵染对意大利蜜蜂幼虫肠道免疫与解毒相关基因表达及蛋白质、脂质和糖类含量的影响

[目的]球囊菌Ascosphaera apis是一种真菌性病原体,可以侵染意大利蜜蜂Apis mellifera ligustica幼虫,导致蜂群患白垩病.本研究通过球囊菌侵染意大利蜜蜂幼虫,探究对意大利蜜蜂幼虫肠道免疫解毒相关基因表达及蛋白质、脂质和糖类含量的影响,分析球囊菌侵染胁迫下意大利蜜蜂免疫应答反应与营养物质...     

意大利蜜蜂工蜂幼虫nkd基因的生物信息学分析及功能研究

【目的】意大利蜜蜂(Apis mellifera ligustica,简称意蜂)是西方蜜蜂(Apis mellifera)的亚种之一。蜜蜂球囊菌(Ascosphaera apis)侵染意蜂幼虫导致白垩病。本研究对西方蜜蜂裸表皮蛋白(naked cuticle, Nkd)进行保守基序预测和系统进化分析,并通过RNAi明确nkd基因对意蜂工蜂幼虫体重及宿主响应蜜蜂球囊菌胁迫的免疫应答的影响,以期丰富西方蜜蜂基因nkd的信息,并揭示意蜂幼虫nkd的功能。【方法】通过MEME软件预测西方蜜蜂和其他9个物种Nkd蛋白的保守基序。采用MEGA X软件对西方蜜蜂及其他9个物种的Nkd蛋白进行系统进化分析。通过饲喂dsRNA对意蜂幼虫肠道内的nkd进行RNAi。使用电子天平对幼虫进行称重。利用RT-qPCR检测nkd基因的干扰效率及免疫基因的相对表达量。【结果】西方蜜蜂与东方蜜蜂、柑橘凤蝶、家蚕和金凤蝶的Nkd蛋白均含有3个保守基序(motif 1、motif 2 和motif 3),说明上述5个昆虫物种的Nkd具有较高的保守性。西方蜜蜂与东方蜜蜂的Nkd蛋白聚为一支,说明二者的亲缘关系近。与dsRNA-egfp组相比,dsRNA-nkd组5日龄和6日龄幼虫肠道内nkd的表达量均极显著下调(P<0.001),干扰效率分别为49.60%和56.40%。另外,dsRNA-nkd组幼虫体重较dsRNA-egfp组显著下降,说明nkd显著影响幼虫体重。RT-qPCR结果显示,4日龄幼虫肠道内abaecin、apidaecin、birc5、defensin-1和PGRP-S2均被激活表达;5日龄幼虫肠道内abaecin、apidaecin、birc5和defensin-1均被激活表达,PGRP-S2的表达受到抑制;6日龄幼虫肠道内abaecin被激活表达,而apidaecin、birc5、defensin-1和PGRP-S2的表达均受到抑制,说明上述5个免疫基因在宿主响应胁迫的过程中呈不同的表达趋势,均参与宿主的免疫应答,nkd与abaecin和apidaecin的表达存在负向调控关系。【结论】西方蜜蜂的Nkd蛋白含有3个保守基序(motif 1、motif 2和motif 3),西方蜜蜂与东方蜜蜂的Nkd蛋白亲缘关系最近,通过饲喂dsRNA能有效干扰意蜂工蜂幼虫肠道内nkd表达,nkd影响意蜂工蜂幼虫体重及宿主对蜜蜂球囊菌胁迫的免疫应答。     

意大利蜜蜂幼虫肠道内球囊菌及其纯培养的高表达基因差异分析

【背景】球囊菌是一种典型的蜜蜂真菌性病原,特异性地侵染蜜蜂幼虫。目前,有关球囊菌在侵染过程中的基因表达规律的信息极为有限。【目的】通过深入分析胁迫意大利蜜蜂(意蜂)幼虫肠道的球囊菌及其纯培养的高表达基因(HEGs)差异,探索球囊菌在胁迫意蜂幼虫肠道后期与胁迫前的基因表达规律。【方法】利用RNA-Seq技术对球囊菌胁迫的意蜂6日龄幼虫肠道(Aam T)及球囊菌纯培养(AaCK)进行深度测序,根据FPKM值筛选得到球囊菌的HEGs,进而通过GO(Gene ontology)及KEGG(Kyoto Encyclopedia of Genes and Genomes)富集分析、Venn分析对上述HEGs进行功能预测和生物学意义挖掘。【结果】AaCK与Aam T的转录组测序共得到105 447 578条原始读段(Raw reads),经过滤得到88 466 344条有效读段(Clean reads),两端平均Q20为97.50%,平均Q30为93.81%。GO富集分析结果显示,AaCK的HEGs富集于26个GO terms,基因富集数最多的是细胞(22 Unigenes),其次是细胞组件(22 Unigenes)和代谢过程(21 Unigenes);Aam T的HEGs富集于22个GO terms,基因富集数最多的是催化活性(23 Unigenes),其次是细胞进程(18 Unigenes)和代谢过程(18 Unigenes)。KEGG代谢通路(Pathway)富集分析显示,AaCK的HEGs富集在109个Pathways上,基因富集数最多的是核糖体(179 Unigenes),其次是氨基酸生物合成(70Unigenes)和碳代谢(62 Unigenes);Aam T的HEGs富集在114个Pathways上,基因富集数量最多的是核糖体(178 Unigenes),其次是碳代谢(116 Unigenes)和氧化磷酸化(112 Unigenes)。Venn分析结果显示,AaCK与Aam T共有的HEGs有260个,二者特有的HEGs分别为2 161个和4 445个。【结论】提供了胁迫意蜂6日龄幼虫肠道的球囊菌及其纯培养的HEGs表达谱,揭示了球囊菌在胁迫前和胁迫后期的基因表达规律,也为阐明球囊菌致病的分子机理提供了有益的信息和线索。     

中华蜜蜂幼虫肠道响应球囊菌胁迫的microRNA应答分析

【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种能够侵染中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫的致死性真菌病原。微小RNA(microRNA,miRNA)可通过在转录后水平靶向抑制或降解mRNA而参与宿主与病原互作过程。本研究旨在对球囊菌胁迫的中蜂6日龄幼虫肠道的差异表达miRNA(DEmiRNA)及其靶基因进行深入分析,进而揭示DEmiRNA在中蜂响应球囊菌胁迫应答过程中的作用。【方法】利用Illumina MiSeq平台对正常及球囊菌胁迫的中蜂6日龄幼虫肠道(AcCK和AcT)进行测序,通过相关生物信息学软件预测DEmiRNA及其靶基因。通过Blast将靶基因注释到GO和KEGG数据库。利用Cytoscape软件构建DEmiRNA与其靶mRNA的调控网络。通过Stem-loop RT-PCR和qPCR验证测序数据的可靠性。【结果】本研究共预测出537个miRNA,其长度分布介于16–35 nt之间,且不同长度的miRNA首位碱基偏向性差异明显。通过Stem-loop RT-PCR证实了10个novel miRNA的表达。AcCK vs AcT比较组共有54个DEmiRNA,包含31个上调和23个下调miRNA,可分别靶向结合6170和8199个靶基因。GO分类结果显示上调和下调miRNA的靶基因分别涉及47和47个条目,富集基因数最多的皆为结合细胞进程和催化活性。KEGG代谢通路(pathway)富集分析结果表明上调和下调miRNA的靶基因分别富集在134和126条pathway,富集基因数最多的均为内吞作用和内质网中的蛋白质加工。调控网络分析结果表明,DEmiRNA及其靶mRNA形成十分复杂的调控关系;31个DEmiRNA可靶向结合51个与泛素介导的蛋白水解相关的mRNA,18个DEmiRNA可靶向结合14个与Jak-STAT信号通路相关的mRNA;miR-1277-x、miR-26-x、miR-27-y、miR-30-x、miR-6052-x等16个miRNA共同参与了上述两条免疫通路的调控。最后,随机挑选3个DEmiRNA进行qPCR验证,结果证明了测序数据的可靠性。【结论】本研究提供了中蜂幼虫肠道在球囊菌胁迫后期的miRNA的表达谱和差异表达信息,揭示了球囊菌与宿主之间在miRNA组学水平存在复杂的互作。miR-6052-x和miR-1277-x作为调控网络的核心可能通过影响细胞凋亡参与宿主的免疫防御,miR-26-x和miR-30-x可能通过调控Jak-STAT信号通路参与宿主的胁迫应答。本研究筛选出的关键DEmiRNA有望作为治疗白垩病的分子靶标。     

球囊菌胁迫中华蜜蜂幼虫肠道过程中病原的转录组学研究

【目的】本研究利用RNA-seq技术对球囊菌胁迫的中华蜜蜂(中蜂)幼虫肠道进行深度测序,经趋势分析得到差异表达基因(DEGs)的显著表达模式,进而对胁迫过程中的球囊菌进行转录组学分析。【方法】利用Illumina HiSeq 2500平台对球囊菌胁迫的中蜂幼虫肠道进行深度测序,并利用相关软件进行了深入分析。最后,通过RT-qPCR对RNA-seq数据进行了验证。【结果】本研究共得到球囊菌的41133932条高质量clean reads。22865个DEGs共聚类为8个基因表达模式,其中,16769个DEGs聚类为2个显著上调趋势与2个显著下调趋势。GO富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于40与37个GO term,基因富集数最多的为细胞进程(2486 unigenes)。KEGG代谢通路(pathway)富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于119和112个pathway,基因富集数最多的分别是氨基酸生物合成(127 unigenes)与核糖体(98 unigenes)。进一步分析表明球囊菌在胁迫中蜂幼虫肠道的过程中通过提高物质合成促进其增殖,而宿主通过抑制球囊菌的蛋白合成抵御病原入侵。富集在MAPK信号通路的11个DEGs的表达水平随着胁迫时间的延长而逐渐下降,推测中蜂幼虫通过抑制该通路而阻遏球囊菌增殖。【结论】本研究不仅为揭示白垩病过程中的球囊菌-中蜂幼虫互作提供了重要信息,也为阐明不同抗性蜂种的球囊菌抗性差异奠定了基础。     

基于small RNA组学分析揭示意大利蜜蜂响应东方蜜蜂微孢子虫胁迫的免疫应答机制

【目的】通过深度测序和组学分析在small RNA组学层面揭示意大利蜜蜂(Apis mellifera ligustica,简称意蜂)响应东方蜜蜂微孢子虫(Nosema ceranae)的免疫应答机制。【方法】利用small RNA-seq技术对正常及N. ceranae胁迫7 d和10 d的意蜂工蜂中肠(Am7CK、Am7T、Am10CK和Am10T)进行深度测序。利用相关生物信息学软件对测序数据进行质控、已知micro RNA (mi RNA)鉴定、新mi RNA预测及mi RNA的结构特征分析。通过Stem-loop RT-PCR验证新mi RNA的表达。按照|log2(Fold change)|≥1和P≤0.05的标准筛选出Am7CKvs.Am7T和Am10CKvs.Am10T比较组的差异表达mi RNA(differentially expressed mi RNA,DEmi RNA)。利用软件预测DEmi RNA靶向结合的m RNA并进行GO和KEGG数据库注释,根据注释信息对细胞和体液免疫相关通路及富集靶m RNA进行统计和分析。根据靶向结合关系构建DEmi RNA和免疫通路相关差异表达m RNA (differentially expressed m RNA,DEm RNA)的调控网络。采用RT-q PCR对数据的可靠性和DEmi RNA的差异表达进行验证。【结果】共获得165895574条原始读段和132028990条有效序列标签,各组的组内Pearson相关性平均在87.92%及以上。共鉴定到928个已知mi RNA和56个新mi RNA。这些mi RNA的长度介于18–28nt,多数的长度为18nt和22nt且首位碱基主要偏向U。验证了12个新mi RNA的真实表达。Am7CKvs.Am7T比较组包含48个上调mi RNA和36个下调mi RNA;Am10CK vs. Am10T比较组包含56个上调mi RNA和51个下调mi RNA。两个比较组的DEmi RNA可分别靶向结合9827个和10720个m RNA。这些靶m RNA可分别注释到50和47条功能条目,以及138和135条KEGG通路。DEmi RNA与免疫通路相关靶m RNA的调控网络分析结果显示,Am7CK vs. Am7T中有26个DEmi RNA靶向与内吞作用等免疫通路相关的10个DEm RNA;Am10CK vs. Am10T中有15个DEmi RNA靶向与MAPK信号通路等免疫通路相关的10个DEm RNA。验证了测序数据和4个DEmi RNA差异表达的可靠性。【结论】研究结果揭示了宿主DEmi RNA可能通过调控物质和能量代谢、细胞和体液免疫对N. ceranae产生应答,但DEmi RNA不参与抗菌肽基因的表达调控;mi R-1-z可能参与宿主的细胞增殖、细胞凋亡和免疫进程;氧化磷酸化通路可能在宿主免疫应答及宿主-病原互作中发挥特殊作用。     

侵染中华蜜蜂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和Tar...     

东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程的nce-miR-23928及其靶基因表达谱

【背景】东方蜜蜂微孢子虫(Nosema ceranae)专性侵染成年蜜蜂中肠上皮细胞而导致的微孢子虫病给养蜂业造成严重损失。【目的】检测东方蜜蜂微孢子虫nce-miR-23928及其靶基因在侵染意大利蜜蜂(Apis mellifera ligustica)工蜂过程的表达谱,为深入探究nce-miR-23928在东方蜜蜂微孢子虫侵染中的功能及调控机制提供依据。【方法】通过RNAhybrid、miRanda和TargetScan软件预测nce-miR-23928的靶基因。使用BLAST工具将上述靶基因比对到基因本体论(geneontology,GO)、京都基因和基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)、Nr和Swiss-Prot数据库以获得相应注释。采用实时荧光定量PCR(realtimequantitativePCR,RT-qPCR)技术检测nce-miR-23928及其靶基因在东方蜜蜂微孢子虫侵染意蜂工蜂过程中的相对表达量。【结果】相较于接种后1 d (1 day post infection, 1 dpi),nce-...