沙葱萤叶甲卵滞育的转录组学分析

【目的】建立沙葱萤叶甲Galeruca daurica滞育卵转录组数据库,挖掘卵滞育相关的基因以及代谢和信号通路,在转录组水平探讨卵滞育的分子机制。【方法】采用Illumina NovaSeq6000高通量测序平台对沙葱萤叶甲滞育卵与解除滞育卵进行转录组测序,并进行生物信息学分析;利用DESeq软件分析沙葱萤叶甲滞育卵与解除滞育卵中的差异表达基因,对差异表达基因进行KEGG通路富集分析;利用qRT PCR技术对10个差异表达基因的表达模式进行验证。【结果】基于沙葱萤叶甲滞育卵与解除滞育卵转录组测序结果,共获得53 389个unigene,其中差异表达基因2 145个,24个差异表达基因与保幼激素信号及脂肪酸生物合成和降解相关。与解除滞育卵相比,滞育卵转录组中1 297个基因上调表达,富集于124条KEGG通路,其中核糖体通路显著富集;848个基因下调表达,富集于73条KEGG通路,其中MAPK信号通路和糖胺聚糖生物合成通路显著富集。qRT-PCR结果表明,随机选取的10个差异表达基因的表达趋势与RNA-Seq转录组测序结果完全一致。【结论】保幼激素,脂肪酸生物合成和降解,核糖体,MAPK信号及糖胺聚糖生物合成等通路可能在沙葱萤叶甲卵滞育调控中起着重要的作用。     

伞裙追寄蝇滞育关联基因的转录组学分析

【目的】本研究旨在探讨调控伞裙追寄蝇Exorista civilis滞育的重要关联基因以及代谢途径,为明确该虫在转录组水平下的滞育分子机制提供理论基础。【方法】采用新一代高通量测序平台Illumina HiSeq^TM2000对伞裙追寄蝇非滞育以及滞育的蛹进行转录组测序分析以及生物信息学分析;应用KAAS在线pathway比对分析工具对筛选出的满足padj＜0.05且fold change≥32或padj＜0.05且fold change≤1/32条件的差异表达基因进行KEGG通路富集分析。【结果】根据测序结果,共获取58 050个基因。差异倍数在32倍以上的滞育关联基因(diapause-associated genes, DAGs)有454个,其中406个表达上调,共涉及134条通路,包括氧化磷酸化、柠檬酸循环及其他重要通路;48个表达下调,涉及32条通路。KEGG通路富集分析结果表明,滞育关联基因主要集中在信号转导、内分泌系统、碳水化合物代谢等途径中。【结论】本研究获得的伞裙追寄蝇转录组数据揭示了其滞育调控的重要代谢途径与关联基因。     

沙葱萤叶甲成虫不同夏滞育阶段的转录组分析(英文)

【目的】本研究旨在揭示内蒙古草原新害虫沙葱萤叶甲Galeruca daurica专性夏滞育相关的重要基因以及代谢通路。【方法】应用RNA-Seq技术,对沙葱萤叶甲成虫不同夏滞育阶段［滞育前期(PD)、滞育期(D)及滞育后期(TD)］进行转录组测序、分析及基因功能预测,基于RNA-Seq数据筛选夏滞育不同阶段差异表达基因;利用qPCR对基于RNA-Seq数据筛选的10个差异表达基因的表达水平进行验证。【结果】从9个文库中获得202 770 198 clean reads,将12 078 060条转录本组装获得82 292 条unigene,平均长度为783.59 bp,N50为1 545 bp。沙葱萤叶甲D vs PD和TD vs D 比较组分别有2 395(2 119上调和277下调)和62(59上调和3下调)个差异基因。KEGG分析表明,D vs PD和TD vs D比较组差异表达基因分别显著富集于糖孝解/糖异生通路和脂肪酸生物合成通路;此外,许多与钙离子信号转导相关的基因在滞育期间差异表达。10个差异表达基因的qPCR分析表明,RNA-Seq与qPCR结果高度一致。【结论】糖孝解/糖异生、脂肪酸生物合成及钙离子信号通路可能在沙葱萤叶甲滞育调节中起着重要的作用。本研究为进一步研究沙葱萤叶甲成虫专性夏滞育的分子机理奠定了基础。     

茶足柄瘤蚜茧蜂碳水化合物代谢相关的滞育关联基因差异表达分析

为明确糖代谢相关途径在茶足柄瘤蚜茧蜂Lysiphlebustestaceipes蛹滞育过程中的作用,揭示滞育调控的分子机制,本试验利用转录组测序技术,对滞育组与非滞育组的茶足柄瘤蚜茧蜂蛹进行转录组测序,并结合生物信息学方法对糖代谢相关途径中的差异表达基因进行了筛选与分析.GO注释到的与碳水化合物代谢条目相关的差异基因共...     

球孢白僵菌诱导的柑橘木虱免疫应答转录组分析

【目的】筛选柑橘木虱Diaphorina citri Kuwayama应对球孢白僵菌侵染的免疫应答及其网络调控基因,以进一步探讨柑橘木虱对球孢白僵菌的免疫防御机制。【方法】采用Illumina高通量测序平台对感染球孢白僵菌24、48、72 h与健康的柑橘木虱转录组进行了测序,利用生物信息学软件对筛选到的差异表达基因进行了功能注释、分类以及参与的信号通路分析。【结果】组装得到138 313条不可延长的非冗余Unigene,其N50和N90分别为2 532 bp和413 bp,平均长度为1 191.26 bp。在CK vs.S24h、CK vs.S48h和CK vs.S72h 3个转录组里分别获得了971、1 671和752个显著差异表达基因(DEGs),GO富集分析表明分别有405、614和542个DEGs显著富集到56、83和60个GO terms中,KEGG pathway分析结果显示分别有98、333和247个DEGs显著富集到10、27和25条代谢通路里。【结论】差异表达基因主要富集在能量代谢、离子运输、转录和翻译调控、生殖和发育调控以及免疫防御反应等相关通路,大部分编码潜在的与免疫识别及调控相关的基因,并从中筛选到5个显著上调表达的免疫相关基因,为开展柑橘木虱免疫应答虫生真菌的研究奠定理论基础。下一步将进行免疫相关基因的q-PCR验证及表达量分析,研究其在柑橘木虱应对球孢白僵菌侵染过程中的作用。     

感染球孢白僵菌的烟粉虱若虫免疫应答转录组分析

【目的】筛选烟粉虱Bemisia tabaci(Gennadius)若虫应对白僵菌Beauveria bassiana侵染的应答基因,以进一步研究烟粉虱免疫反应的分子机制。【方法】采用新一代高通量测序技术对感染和非感染白僵菌的烟粉虱4龄若虫进行了测序分析,并筛选了差异表达基因;利用生物信息学工具对转录组测序得到的基因进行了功能注释、分类以及参与的信号通路展示。【结果】组装得到非冗余Unigene 232 554个,其N50和N90分别为1 153 bp和260 bp,平均长度为67 424 bp。对所有的基因进行差异表达分析发现:以P<0.05为标准筛选得到了1 166个差异表达基因,其中上调表达基因有474个,下调表达基因有692个,其中,与免疫相关的基因有405个;GO富集分析发现:有416个GO term有富集现象,包括156个生物学过程(66 402个Unigenes),89个细胞组分(27 645个Unigenes)和154个分子功能(73 417个Unigenes);KEGG代谢通路分析将烟粉虱转录组1 145个DEG匹配到309个通路上,其中76个通路得到了富集。【结论】对405个可能参与到烟粉虱若虫对白僵菌侵染的免疫识别和防御的基因进行了测序。研究结果为从分子水平上开展应用虫生真菌防治烟粉虱的研究奠定信息学基础。     

意大利蝗卵发育过程保幼激素和蜕皮激素含量的变化及其代谢相关基因的表达谱

【目的】鉴于保幼激素(juvenile hormone, JH)与蜕皮激素(ecdysone, Ecd)在昆虫生长发育中的重要作用,本研究旨在明确这两种激素及其代谢相关基因对意大利蝗Calliptamus italicus卵发育过程的调控机制。【方法】采用ELISA技术检测意大利蝗卵发育过程中JH与Ecd含量的变化;采用qRT-PCR技术检测JH与Ecd代谢通路的重要基因(JHE,JHEH,DIB和EcR)在意大利蝗卵发育过程的表达模式。【结果】意大利蝗卵滞育阶段(阶段Ⅳ-Ⅵ)的JH含量显著高于早期发育阶段(阶段Ⅰ-Ⅲ)的,滞育后发育阶段(阶段Ⅶ-Ⅸ) JH含量显著下降;Ecd含量于滞育早期(阶段Ⅳ)显著上升,而后显著下降,滞育后发育阶段再次上升。JHE的表达量在意大利蝗卵早期发育阶段和滞育后发育阶段均呈先升后降的趋势,滞育阶段JHE的表达量较低;JHEH表达量在意大利蝗卵早期发育阶段也先升后降,滞育阶段变化不显著,滞育后发育阶段显著升高;DIB表达量在意大利蝗卵滞育阶段高于早期发育阶段的,在滞育后发育阶段下降;EcR表达量在意大利蝗卵早期发育阶段及滞育阶段均无显著变化,滞育后发育阶段...     

伞裙追寄蝇滞育关联基因和滞育关联代谢物的联合分析

为探索伞裙追寄蝇蛹滞育调控的分子机制,本文对伞裙追寄蝇Exorista civilis滞育蛹和非滞育蛹进行转录组测序以及代谢组检测,鉴定其关键的滞育关联基因(diapause-associated genes, DEGs)和滞育关联代谢物(diapasuse-associated metabolites, DEMs)。本研究基于高通量测序以及液质联用技术,通过筛选,在转录组中获得了差异表达基因7 513个,在代谢组中,获得差异代谢物501个,其中氨基酸占比最多。将所有差异表达基因与差异表达代谢物同时向KEGG映射,获得两者共同的pathway信息,明确差异表达基因与差异代谢物共同参与的主要生化途径和信号转导途径。本研究通过转录组和代谢组的联合分析,在正、负离子模式下,差异表达基因与差异代谢物共同富集到70条通路。在负离子模式下,滞育关联因子主要参与氨基酸代谢和神经系统;在正离子模式下,滞育关联因子主要参与消化系统和信号转导途径。本研究重点分析了柠檬酸循环、cAMP信号通路、氨酰-tRNA的生物合成,将为进一步深入研究伞裙追寄蝇滞育调控的分子机制奠定理论基础。     

七星瓢虫滞育关联基因的转录组学分析

对七星瓢虫正常发育、滞育以及滞育解除的雌成虫进行RNA测序,并对筛选出来的滞育相关基因进行KEGG通路富集分析,从分子水平解析七星瓢虫滞育发机理。本研究以正常发育产卵、滞育30 d以及滞育贮存30 d后解除产卵的七星瓢虫雌成虫为研究对象,分别抽提RNA,合成c DNA,构建c DNA文库,文库检测合格后在Illummina Hiseq 2500测序仪上进行双向测序。根据测序结果,共获取unigene 82820个。采用两两比较法对正常发育组和滞育组、滞育组和滞育解除组进行差异表达分析,分别获得差异表达基因3501个和1427个。深入分析两组比对结果,将在滞育组上调且滞育解除组下调的unigene定义为滞育关联基因,共有443个基因为滞育关联基因。应用KEGG KAAS在线pathway比对分析工具对滞育关联基因进行通路富集分析,结果发现这些基因主要集中在碳水化合物代谢、脂质代谢以及信号转导等途径中。     

茶足柄瘤蚜茧蜂滞育和非滞育蛹中与能量代谢相关的差异表达蛋白

【目的】本研究旨在从蛋白质组整体层面阐明茶足柄瘤蚜茧蜂Lysiphlebus testaceipes蛹滞育背后的多蛋白调控,重点筛选与能量代谢相关的滞育关联蛋白并分析其功能,有助于更好地理解茶足柄瘤蚜茧蜂蛹滞育的代谢机制。【方法】利用同位素标记相对和绝对定量(isobaric tags for relative and absolute quantification, iTRAQ)技术比较了茶足柄瘤蚜茧蜂滞育蛹与非滞育蛹的蛋白含量;利用GO, KEGG网络数据库等生物信息学方法分析鉴定茶足柄瘤蚜茧蜂滞育蛹与非滞育蛹中差异表达蛋白(differentially expressed proteins, DEPs)。【结果】分析得到茶足柄瘤蚜茧蜂滞蛹与非滞育蛹DEPs有135个,包括滞育蛹中上调表达蛋白有38个,下调表达蛋白有97个。GO和KEGG富集分析表明,与天冬氨酸转运、L-谷氨酸转运、胆碱脱氢酶活性和胆碱生物合成甘氨酸甜菜碱条目以及氧化磷酸化通路相关的蛋白上调表达。【结论】氧化磷酸化通路相关蛋白在茶足柄瘤蚜茧蜂滞育过程中呈显著上调表达,说明能量代谢与该蜂滞育密切相关,并推测氧化磷酸化过程除为滞育蛹提供维持生命基本活动的主要能量外,还提供热能。     

西方蜜蜂工蜂蛹响应低温胁迫的差异表达基因分析

【目的】本研究旨在从转录组水平筛选和分析西方蜜蜂Apis mellifera工蜂不同时期蛹对低温胁迫反应的差异表达基因(DEGs)。【方法】将西方蜜蜂工蜂封盖后3 d预蛹以及封盖后6和9 d蛹分别置于低温环境(20℃)和最适发育温度(35℃)中4 h,分别作为处理组和对照组,通过转录组学技术筛选低温处理组与对应的对照组之间的DEGs,并进行GO功能分类和KEGG通路分析。利用RT qPCR分别对封盖后3 d预蛹以及封盖后6和9 d蛹的8, 6和5个DEGs的表达谱进行验证。【结果】与对照组相比,封盖后3 d预蛹以及封盖后6和9 d蛹受低温胁迫后的DEGs分别有220, 50和26个;GO功能分类发现,DEGs富集数最多的条目为代谢进程、细胞进程、催化活性和结合,封盖后3 d预蛹的DEGs在生物学进程调控、细胞部分和细胞器等有较多富集。KEGG通路分析显示,处理组和对照组间西方蜜蜂各日龄DEGs在整体概述图、氨基酸代谢、信号转导、运输和分解代谢有富集。封盖后3 d预蛹以及封盖后6和9 d蛹受低温胁迫后共有的DEG 3-磷酸肌醇依赖性蛋白激酶1基因PDK1上调表达,同时富集在自噬-动物、mTOR和FoxO信号通路。低温胁迫后,封盖后3 d预蛹的胰岛素受体底物1-B基因IRS1-B和Kruppel同源物1基因Kr-h1上调表达,而核激素受体FTZ-F1基因Ftz-F1与蜕皮启动激素基因Eth显著下调表达,说明封盖后3 d预蛹响应低温细胞自噬和蜕皮受到抑制程度更大。在低温胁迫后封盖后6 d蛹中与昆虫角质层着色与免疫相关的酪氨酸羟化酶基因TyHyd下调表达;与对照组相比,在低温胁迫后封盖后9 d蛹中DEGs最少,说明在3个蛹期中,其受低温的影响较小。【结论】本研究测定了西方蜜蜂3个不同发育阶段蛹响应低温的DEGs,结果显示大部分DEGs为阶段特有,说明西方蜜蜂不同发育阶段响应低温的机制不同。一些共有DEGs以及阶段特有DEGs的功能研究和其作用机制是进一步研究蜜蜂对低温响应机制的重点内容,为探究蜜蜂蛹期响应低温胁迫的分子机制提供了基础数据。     

低温胁迫下意大利蜜蜂预蛹的差异表达基因趋势分析

【目的】蜜蜂是典型的具有发育狭温性的全变态昆虫。本研究以对低温最敏感的意大利蜜蜂Apis mellifera ligustica预蛹为研究对象,通过低温胁迫不同时间的差异表达基因(differentially expressed genes, DEGs)趋势分析,探讨低温胁迫对蜜蜂发育影响的关键基因。【方法】对3日龄意大利蜜蜂封盖子预蛹进行20℃低温胁迫18 h(T18)和36 h(T36),以未经低温胁迫的预蛹为对照(CK),通过Illumina HiSeq~(TM)平台进行转录组学测定。利用Short Time-series Expression Miner(STEM)软件对2个处理组与对照组比较共有的差异表达基因进行趋势分析,再进一步对显著富集趋势模式中富集的差异表达基因进行GO分类和KEGG pathway分析。利用RT-qPCR对随机挑选的5个DEGs的表达模式进行验证。【结果】对检测到1 062个T18 vs CK和T36 vs CK共有的DEGs进行趋势分析,发现3个显著的基因表达模式,包括2个上调表达模式(Profile 6,有539个基因;Profile 7,有271个基因),1个下调表达模式(Profile 1,有183个基因)。对3个显著富集趋势模式DEGs分别进行GO富集分析和KEGG pathway分析,在Profile 6中找到同时富集在FoxO信号通路和寿命调节信号通路上的胰岛素样肽基因ILP和叉头蛋白O基因FoxO持续上调,说明低温胁迫会影响蜜蜂预蛹蜕皮激素信号传递;在Profile 7中CREB结合蛋白基因CBP持续上调,蜜蜂预蛹受到低温胁迫会影响细胞发育;在Profile 1中细胞色素P450基因CYP450 306a1 (phm)和WNT1显著下调,说明化蛹时受到低温胁迫会影响蜕皮激素的合成。通过RT-qPCR分析挑选的5个基因的表达结果和高通量测序结果一致。【结论】通过对蜜蜂响应低温胁迫的差异表达基因趋势分析发现,蜜蜂胰岛素和蜕皮激素共同调控的FoxO可能是低温胁迫抑制蜜蜂化蛹的一个关键基因。本研究为探索低温胁迫影响蜜蜂发育变态的分子调控机制奠定了基础。     

梨小食心虫幼虫中肠中高表达消化酶和解毒酶基因的筛选

【目的】挖掘梨小食心虫Grapholita molesta幼虫中肠中高表达消化酶和解毒酶基因,为今后研究以肠道为靶标的新型农药和转基因作物提供理论依据。【方法】基于梨小食心虫4龄幼虫中肠转录组高通量测序数据的FPKM值,筛选高表达基因,进行GO功能注释和KEGG通路富集分析,并使用BLAST软件进行比对筛选高表达的消化酶和解毒酶基因,利用MEGA对这些高表达的消化酶和解毒酶及其他鳞翅目昆虫的同源蛋白进行系统发育分析。利用qRT-PCR技术对梨小食心虫幼虫不同龄期中肠中的高表达代表性消化酶和解毒酶基因表达量进行定量分析和验证。【结果】在GO数据库中注释了103 677个在梨小食心虫4龄幼虫中肠中高表达基因,包括细胞组分、分子功能和生物学进程三大类功能共41个分支。KEGG通路分析表明,10 846个高表达基因参与了5类生化代谢通路。筛选到具有完整开放阅读框的消化酶基因17个[5个胰蛋白酶(trypsin, TRY)基因、3个氨肽酶(aminopeptidase, APN)基因和9个羧肽酶(carboxypeptidase, CP)基因]和解毒酶基因32个[11个谷胱甘肽S-转移酶(glu...     

意大利蝗卵发育过程中血蓝蛋白基因表达分析

【目的】意大利蝗Calliptamus italicus是新疆草原的主要优势危害种,以卵在土壤中越冬。呼吸代谢可反映蝗卵的生理状态,呼吸蛋白对于呼吸系统不完善的蝗卵尤为重要。本研究旨在明确意大利蝗卵发育过程中血蓝蛋白基因的表达情况。【方法】采用实时荧光定量PCR方法检测不同发育阶段的蝗卵以及1龄蝗蝻的血蓝蛋白2个亚基基因Hc1和Hc2的表达量。【结果】根据解剖形态观察,将意大利蝗越冬卵的整个发育过程分为10个阶段,包括9个卵发育阶段(C-Ⅰ-C-Ⅹ)和1龄蝗蝻阶段(C-Ⅹ)。Hc1和Hc2在越冬蝗卵各发育阶段以及1龄蝗蝻中均有表达。其中,在蝗卵早期发育阶段(C-Ⅰ, C-Ⅱ和C-Ⅲ),Hc1表达量逐渐增加,C-Ⅲ阶段表达量显著高于C-Ⅰ和C-Ⅱ阶段;滞育阶段(C-Ⅳ, C-Ⅴ和C-Ⅵ),胚胎发育停滞,Hc1表达量较C-Ⅲ,C-Ⅶ和C-Ⅷ阶段低;滞育后发育阶段(C-Ⅶ和C-Ⅷ),蝗卵解除滞育,快速发育,Hc1表达量较早期发育阶段和滞育阶段高,其中,C-Ⅷ阶段Hc1表达量最高(212.3156±10.5470),显著高于其他所有阶段;1龄蝗蝻(C-Ⅹ)的Hc1表达量最低,为0.4017±0.1010。Hc2表达量在C-Ⅴ阶段最高(679.7511±54.5719),显著高于其他所有阶段;除C-Ⅴ阶段外,其他各阶段之间Hc2表达量差异均不显著。Hc1在蝗卵滞育后阶段高表达,而Hc2在蝗卵滞育阶段高表达。【结论】血蓝蛋白亚基基因Hc1和Hc2在整个意大利蝗卵发育过程均有表达,且具有阶段特异性。Hc1与Hc2协同作用为蝗卵发育供氧,其中,Hc1主要负责蝗卵滞育后发育期间的氧气运载,而Hc2主要维持滞育期间的氧气运载,且载氧效率较低。研究结果可为进一步探讨意大利蝗卵的抗逆机制提供科学依据。     

转录组分析揭示东方蜜蜂微孢子虫侵染意大利蜜蜂的分子机制

【目的】本研究旨在通过差异表达基因(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比较组分别鉴定出1 397, 1 497和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转运蛋白等侵染相关因子编码基因在病原增殖中扮演重要角色,为阐明东方蜜蜂微孢子虫的侵染机制提供了基础。