【目的】白垩病是困扰养蜂生产的顽疾。目前,尚无利用二代测序技术研究中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫白垩病的报道。本研究利用RNA-seq技术对健康(Ac CK)及球囊菌(Ascosphaera apis)胁迫的中蜂4日龄幼虫肠道(Ac T)进行深度测序,在转录组水平研究中蜂幼虫在球囊菌胁迫早期的胁迫应答。【方法】通过Illumina Hi Seq 2500平台对Ac CK和Ac T进行双端(PE125)测序,首先对测序数据进行质控和评估,利用edge R软件进行差异表达基因(

利用RNA-seq技术对胁迫中华蜜蜂(简称中蜂)6日龄幼虫肠道的球囊菌及其体外培养进行深度测序,根据基因的FPKM值筛选得到高表达基因(HEGs),进而对HEGs进行GO及KEGG代谢通路(pathway)富集分析.Illumina测序数据经质控和过滤得到100814558条有效读段(clean reads),平均Q30均在93.81%以上.GO富集分析结果显示处理组(Aac T)的HEGs富集于39个GO term,基因富集数最多的是细胞(1872 unigene)、细胞组件(1872 unigene)

利用前期获得的中华蜜蜂(Apis cerana cerana)幼虫肠道转录组数据对东方蜜蜂(Apis cerana)基因组的已注释基因进行结构优化、对未注释的新基因进行预测、分析和鉴定:利用Cufflink软件对reads进行转录本重构,通过与参考转录本序列进行比较,共对3 366个东方蜜蜂的基因结构进行了优化;利用Cuffcompare软件将重构转录本与参考基因组进行比对,共预测出527个东方蜜蜂的新基因,其中234个新基因在Nr数据库中存在注释信息;随机选取12个新基因进行PCR鉴定,有11个能够扩增出符合预期的目的片段,表明多数预测出的新基因真实存在;GO分类和KEGG代谢通路分析结果显示上述新基因可能在东方蜜蜂的生长发育、物质代谢、遗传信息传递等方面具有重要功能。本研究结果丰富和完善了东方蜜蜂基因组的基因结构与功能注释信息,也为新基因的功能研究打下基础。

【目的】微小RNA(microRNA,miRNA)是一类重要的基因表达调控因子,可影响宿主与病原间的互作过程。蜜蜂球囊菌(Ascosphaera apis)是一种特异性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)幼虫肠道在球囊菌胁迫前期的差异表达miRNA(differentially expressed miRNA,DEmiRNA)及其靶基因进行深入分析,在miRNA组学水平探究意蜂幼虫在球囊菌侵染前期的胁迫应答,并通过构建显著DEmiRNA的调控网络筛选出与宿主应答相关的关键miRNA。【方法】利用small RNA-seq(sRNA-seq)技术对正常及球囊菌侵染的意蜂4日龄幼虫肠道(AmCK和AmT)进行高通量测序,首先对原始数据进行质控和评估,随后将过滤后的数据与西方蜜蜂(Apis mellifera)参考基因组进行比对;将比对上的序列标签(tags)注释到miRBase数据库,得出已知miRNA的表达量;通过TPM(tags per million)算法对各样本中miRNA的表达量进行归一化处理,以|log_2 fold change|≥1且P≤0.05作为标准筛选得到显著DEmiRNA;利用TargetFinder软件预测显著DEmiRNA的靶基因,并对其进行GO和KEGG代谢通路(pathway)富集分析。利用Cytoscape软件对miRNA-mRNA调控网络进行可视化。最后,利用茎环反转录PCR(Stem-loop RT-PCR)和荧光定量PCR(qPCR)验证测序数据的可靠性。【结果】AmCK和AmT样品的测序分别得到13 553 302和10 777 534条原始读段(raw reads),经严格过滤后得到的有效读段(clean reads)数分别为13 186 921和10 480 913条。各样品的生物学重复间的Pearson相关性系数分别在0.9822和0.9508以上。共有10个显著DEmiRNA,包括4个上调miRNA和6个下调miRNA。显著DEmiRNA在AmT的整体表达水平低于AmCK。10个显著DEmiRNA可靶向结合3 788个靶基因,其中上调miRNA的1 240个靶基因可注释到GO数据库中的39个GO条目,主要富集在结合、细胞进程、代谢进程和应激反应等;下调miRNA的749个靶基因可注释到34个GO条目,主要富集在细胞进程、结合、代谢进程和应激反应等。KEGG数据库注释结果显示,上调miRNA和下调miRNA的靶基因分别注释到95和66条代谢通路,富集基因数最多的分别是Wnt信号通路、Hippo信号通路、光传导以及内吞作用、磷脂酰肌醇信号系统、嘌呤代谢。对于上调和下调miRNA,分别有31和52个靶基因注释到内吞作用,15和7个靶基因注释到泛素介导的蛋白水解,11和5个靶基因注释到Jak-STAT信号通路,1和3个靶基因注释到MAPK信号通路。显著DEmiRNA与靶mRNA之间形成复杂的调控网络,7个显著DEmiRNA靶向结合96个与Wnt信号通路相关的mRNA,8个显著DEmiRNA靶向结合55个与内吞作用相关的mRNA。Stem-loop RT-PCR和qPCR结果验证了测序数据的可靠性。【结论】对意蜂幼虫肠道在球囊菌侵染前期的DEmiRNA及其靶基因进行预测和分析,并构建和分析了DEmiRNA-mRNA调控网络,研究结果提供了宿主miRNA的表达谱和差异表达信息,揭示了DEmiRNA通过调控细胞生命活动、新陈代谢以及部分细胞和体液免疫等生物学过程参与宿主的胁迫应答。miR-4331-y、miR-4968-y、miR-8440-y、novel-m0023-5p和novel-m0025-3p共同参与了宿主的Wnt信号通路和内吞作用的调控,可作为白垩病治疗的潜在分子靶标。

为探究中华蜜蜂(Apis cerana cerana,以下简称中蜂)与意大利蜜蜂(Apis mellifera ligustica,以下简称意蜂)幼虫的球囊菌(Ascospahaera apis)抗性差异产生的原因,利用Venn分析、GO分类分析、KEGG代谢通路分析以及Ka/Ks分析对二者的转录组进行比较研究。Venn分析结果显示,中蜂与意蜂的同源基因有17 656个,归属于8 111个基因家族,二者的特有基因分别有1 158和241个,分别归属于468和86个基因家族。GO分类结果显示,中蜂和意蜂的特有基因分别富集在38和28个GO条目。KEGG代谢通路富集分析结果显示,中蜂和意蜂的特有基因分别富集于96和21个代谢通路;进一步分析发现中蜂的特有基因富集在内吞作用、溶酶体、泛素介导的蛋白水解和黑化作用等4个细胞免疫通路,以及MAPK信号通路和Toll-like受体信号通路等2个体液免疫通路,而意蜂仅有1个特有基因富集在MAPK信号通路。Ka/Ks分析结果显示受到强烈正向选择、弱正向选择、中性选择和负向选择的单拷贝同源基因分别有37、281、2 630和3 269个。对受到强烈正向选择的基因进行GO分类和KEGG代谢通路富集分析,GO结果显示中蜂和意蜂的单拷贝同源基因富集的GO条目类型相同;KEGG结果显示中蜂的单拷贝基因仅富集在氧化磷酸化。上述结果表明免疫相关基因数量的差异是二者的球囊菌抗性差异产生的重要原因之一,中蜂在与球囊菌的协同进化过程中可能通过提高能量利用从而限制球囊菌的增殖。本研究结果可为阐明中蜂及意蜂幼虫的球囊菌抗性差异产生的分子机制提供参考。

温度是影响蜜蜂发育的关键因子之一.为研究低温胁迫对西方蜜蜂幼虫基因表达的影响,利用高通量转录组测序技术分别对最适发育温度和低温处理的西方蜜蜂工蜂6日龄幼虫的基因表达情况进行分析,获得差异表达基因(DEGs),并对DEGs进行GO功能分类和KEGG富集分析.结果显示,低温处理后共检测到98个DEGs.GO功能分类显示,DEGs富集在21条二级GO功能注释上,富集最多DEGs的GO功能注释为结合、细胞进程、代谢进程、催化活性;KEGG富集分析显示,上调和下调DEGs分别富集在25和3条通路上,主要包括甘氨酸-丝氨酸-苏氨酸代谢、核黄素代谢、昆虫激素生物合成等代谢通路,以及与信号转导相关的Hippo、Notch等信号通路.qPCR结果显示,转录组测序获得的DEGs IP3K、Ex、CYP18A1、GOD、GCDH表达差异显著.本研究结果揭示西方蜜蜂工蜂6日龄幼虫主要通过提高抗氧化能力和减慢发育速度应对低温.

【目的】利用small RNA-seq(sRNA-seq)技术和生物信息学方法对中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫肠道的微小RNA(microRNA,miRNA)进行全转录组鉴定和分析,旨在丰富中蜂的miRNA信息,并为深入研究miRNA调控中蜂幼虫肠道发育的分子机理提供依据。【方法】利用sRNA-seq技术对中蜂4、5和6日龄幼虫肠道样品(Ac1、Ac2和Ac3)进行测序,通过数据质控获得有效标签序列(clean tags)。采用Blast工具将clean tags连续比对东方蜜蜂(Apis cerana)基因组和miRBase数据库,以鉴定保守miRNA和新miRNA。采用TPM法对miRNA的表达量进行归一化处理。通过GraphPad Prism7软件统计各组肠道样品中sRNA占比、miRNA长度分布及首位碱基偏向性。利用相关软件预测上述miRNA靶向的mRNA并进行GO和KEGG数据库注释。进一步根据靶向结合关系构建和分析注释到发育和免疫相关通路的基因及其靶向miRNA的调控网络,并利用Cytoscape软件进行可视化。利用茎环反转录PCR(Stem-loop RT-PCR)、分子克隆和Sanger测序验证miRNA的表达和序列的真实性。【结果】共鉴定到中蜂的371个保守miRNA和64个新miRNA;这些miRNA的长度介于18—25 nt且首位碱基主要偏向于U;上述miRNA共靶向14 750条mRNA,涉及离子结合、金属离子结合、细胞膜、细胞膜组件和单一有机体进程等2 270个GO条目,以及内吞作用、细胞凋亡、mTOR信号通路、RNA转运和昆虫激素的生物合成等332条KEGG通路。进一步分析结果显示156个miRNA与注释到Wnt、Hippo、Notch和mTOR等生长发育相关通路的67个靶基因存在调控关系,145个miRNA与注释到Toll、Imd/JNK、Jak-STAT和抗菌效应因子等免疫相关途径的21个靶基因存在调控关系。Stem-loop RT-PCR结果显示miR-8-y、miR-9-z、miR-14-y、miR-281-y、miR-283-x和miR-306-x均能扩增出预期的特异性片段;Sanger测序结果显示上述6个miRNA的序列与深度测序结果一致。【结论】提供了中蜂miRNA的数量、结构特征和表达谱;揭示中蜂幼虫肠道miRNA潜在调控诸多生命进程与细胞活动;中蜂幼虫肠道的部分miRNA可通过靶向结合相应的mRNA参与调节发育和免疫相关途径。

对健康和患病蜜蜂幼虫血淋巴进行SDS-PAGE电泳和蛋白酶、酯酶的活性染色,以研究蜜蜂球囊菌分泌的胞外酶种类及对蜜蜂侵染的作用.结果表明:健康蜜蜂幼虫血淋巴中的蛋白含量丰富,主要由4种高分子质量的蛋白组成,而患病幼虫血淋巴中的蛋白含量很少,主要蛋白组分被降解;在患病幼虫血淋巴中可以检测到健康幼虫所没有的多种蛋白酶和酯酶的活性.可见,蜜蜂球囊菌在侵染蜜蜂幼虫的过程中涉及到胞外酶的作用,这些胞外酶参与降解蜜蜂的组织,为蜜蜂球囊菌的生长提供营养.

基于前期已获得的意大利蜜蜂(Apis mellifera ligustica)工蜂中肠转录组数据,在全转录组水平上对意大利蜜蜂的可变剪接基因(alternatively spliced gene, ASG)进行鉴定和分析,旨在揭示ASG在宿主响应东方蜜蜂微孢子虫(Nosema ceranae)胁迫中的作用.在正常意大利蜜蜂7和10日龄工蜂中肠(Am7CK、Am10CK)和东方蜜蜂微孢子虫胁迫的意大利蜜蜂7和10日龄工蜂中肠(Am7T、Am10T)样品中检测到可变剪接事件数分别为73 383、85 618、79 813和74 693次,对应的基因数分别为9 586、10 063、9 829和9 622个,平均每个ASG上发生的可变剪接事件数分别为7.66、8.50、8.12和7.76次.各样品的可变剪接事件中外显子跨越的数量最多,分别为4 411、5 247、4 993和4 629个.Venn分析结果显示,Am7CK、Am10CK、Am7T和Am10T中的共有ASG数为9 056个,特有ASG数分别为100、153、275和95个.GO分类结果显示,这些共有ASG分布在50个功能条目,特有ASG分别分布于26、30、33和24个功能条目.KEGG代谢通路富集分析结果显示,上述共有ASG富集在325个代谢通路,包括内质网蛋白加工、核糖体和RNA转运等;特有ASG分别富集在31、104、126和22个代谢通路.研究结果提供了意大利蜜蜂响应东方蜜蜂微孢子虫胁迫过程中的ASG数量和种类信息,揭示宿主ASG在胁迫响应中参与了新陈代谢和细胞免疫,为关键ASG的筛选和功能研究建立了数据基础.

【目的】分析鉴定中华蜜蜂(Apis cerana cerana)幼虫是否含有与西方蜜蜂(Apis mellifera)幼虫相同的幼虫利它素和信息素及其在不同日龄阶段幼虫的分布情况。【方法】以中华蜜蜂为试验材料,取2日龄、6日龄(即将封盖)和7日龄(封盖后)的工蜂幼虫及2日龄、7日龄(即将封盖)和8日龄(封盖后)的雄蜂幼虫,经过正己烷浸泡碾碎,取上清液并进行层析过柱分离,应用气相色谱来研究中华蜜蜂幼虫利它素和信息素。【结果】本研究首次发现中华蜜蜂幼虫信息素含有甲基棕榈酸酯(MP)、甲基油酸酯((MO)、甲基硬脂酸酯(MS)、甲基亚油酸酯(ML)、甲基亚麻酸酯(MN)、乙基棕榈酸酯(EP)、乙基油...

【目的】通过对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂中肠响应东方蜜蜂微孢子虫(Nosema ceranae)胁迫的差异表达基因(differentially expressed gene,DEG)及免疫通路进行深入细致的分析,揭示宿主的细胞和体液免疫应答,为关键免疫应答基因的筛选及功能研究打下基础。【方法】基于前期获得的正常及东方蜜蜂微孢子虫胁迫的意蜂7日龄和10日龄工蜂中肠(Am7CK、Am7T、Am10CK、Am10T)转录组数据,按照FDR≤1,P≤0.05和|log2 fold change|≥1的标准筛选出各组的DEG,利用相关生物信息学软件对DEG进行Pearson相关性分析、Venn分析、GO分类和KEGG代谢通路富集分析,进而对免疫通路富集基因进行统计和分析,利用实时荧光定量PCR(real-time quantitative PCR,RT-q PCR)验证转录组数据的可靠性。【结果】差异分析结果显示,Am7CK vs Am7T比较组包含472个上调基因和385个下调基因;Am10CK vs Am10T比较组包含611个上调基因和360个下调基因。Venn分析结果显示,两个比较组特有的DEG分别为739和853个,共有的DEG为118个。GO分类结果显示,Am7CK vs Am7T比较组中上调和下调基因分别涉及23和29个功能条目,其中富集上调基因数最多的前5位分别为结合、催化活性、代谢进程、细胞进程和单组织进程;富集下调基因数最多的前5位分别为代谢进程、单组织进程、催化活性、细胞进程和结合。Am10CK vs Am10T比较组中上调和下调基因分别涉及36和26个功能条目,其中富集上调基因数最多的前5位分别为单组织进程、结合、细胞进程、催化活性和代谢活性;富集下调基因数最多的前5位分别为结合、细胞进程、催化活性、代谢进程和单组织进程。KEGG代谢通路富集分析结果显示,Am7CK vs Am7T比较组中上调和下调基因分别富集在38和33条代谢通路,富集上调基因数最多的前5条分别为胆汁分泌、内质网蛋白加工、泛素介导的蛋白水解、PI3K-Akt信号通路和神经营养因子信号通路;富集下调基因数最多的前5条分别为胞质DNA传感通路、嘌呤代谢、嘧啶代谢、RNA聚合酶和核糖体;涉及泛素介导的蛋白水解等3条细胞免疫通路,以及PI3K-Akt信号通路等7条体液免疫通路。Am10CK vs Am10T比较组中上调和下调基因分别富集在54和43条代谢通路,富集上调基因数最多的前5条分别为Hippo信号通路、药物代谢-细胞色素P450、P450对外源物质的代谢、泛素介导的蛋白水解和鞘脂类代谢;富集下调基因数最多的前5条分别为m RNA监测、鞘脂类信号通路、果糖和甘露糖代谢、半乳糖代谢和鞘脂类代谢;涉及泛素介导的蛋白水解等7条细胞免疫通路,以及NF-κB信号通路等2条体液免疫通路。RT-q PCR验证结果显示6个随机挑选的DEG的表达量变化趋势与测序结果一致,证实了本研究中测序数据的可靠性。进一步分析发现意蜂7日龄和10日龄工蜂中肠的NF-κB信号通路均被东方蜜蜂微孢子虫激活,随即启动了3种抗菌肽apidaecin、defensin-1和hymenoptaecin的合成,体现了它们在宿主抵御东方蜜蜂微孢子虫入侵中的重要性。【结论】在转录组水平解析了意蜂工蜂中肠对东方蜜蜂微孢子虫的免疫应答,揭示宿主在胁迫前期同时作出细胞和体液免疫应答,前者可能在抵御病原入侵方面发挥主要作用;宿主的细胞免疫在胁迫后期持续增强,但体液免疫较大程度地减弱;泛素介导的蛋白水解通路及富集DEG,NF-κB信号通路及富集DEG,以及抗菌肽编码基因apidaecin、defensin-1和hymenoptaecin可能在意蜂工蜂对东方蜜蜂微孢子虫的免疫应答中起到关键作用。

【目的】微小RNA(microRNA,miRNA)是一类在转录后水平对mRNA进行负调控的关键调控因子。本研究旨在通过分析意大利蜜蜂(Apis mellifera ligustica,简称意蜂)幼虫肠道发育过程中差异表达miRNA(differentially expressed miRNA,DEmiRNA)及其调控网络,提供miRNA的表达谱和差异表达信息,揭示DEmiRNA在幼虫肠道发育中的作用。【方法】利用small RNA-seq(sRNA-seq)技术对意蜂4、5和6日龄幼虫肠道样品(Am4、Am5和Am6)进行测序,将质控后的数据与西方蜜蜂(Apis mellifera)参考基因组进行比对,然后将比对上的序列标签(tags)注释到miRBase数据库,并利用TPM(tags per million)算法归一化处理所得miRNA的表达量,再通过相关生物信息学软件对miRNA进行表达量聚类、前体二级结构预测及差异表达分析。利用TargetFinder软件预测DEmiRNA的靶基因,使用Blast软件对靶基因进行GO和KEGG数据库注释,进而通过Cytoscape软件构建mi RNA-m RNA的调控网络。采用茎环实时荧光定量PCR(Stem-loop RT-qPCR)验证测序数据的可靠性。【结果】意蜂幼虫肠道样品的测序分别得到10 841 644、12 037 678和9 230 496条有效序列标签;Am4 vs Am5比较组包含16个上调和10个下调mi RNA,Am5 vs Am6比较组包含5个上调和7个下调mi RNA。其中,novel-m0031-3p为两个比较组所共有,并结合5个与蜕皮激素诱导蛋白相关的靶基因;二者的特有DEmi RNA数分别为25和11个。Am4 vs Am5的26个DEmi RNA结合5 742个靶基因,其中2 725个靶基因可注释到GO数据库中的46个GO term,并主要富集在结合、细胞进程、代谢进程和单组织进程等;Am5 vs Am6的12个DEmi RNA结合3 733个靶基因,且其中2 725个靶基因富集在结合、细胞进程、单组织进程和代谢进程等41个GO term;此外,两个比较组中的DEmi RNA分别有1 046和676个靶基因可注释到116和92条KEGG代谢通路,且Am4 vs Am5比较组的DEmi RNA的靶基因富集在Wnt信号通路、Hippo信号通路、嘌呤代谢和内吞作用等通路上的数量均高于Am5 vs Am6比较组。进一步分析结果显示Am4 vs Am5中的上调和下调mi RNA可分别结合611和85个靶基因,其中ame-mi R-6052结合的靶基因数最多,可通过结合5个靶基因,参与对细胞色素P450的调控;mi R-281-x结合49个靶基因,并间接调控组氨酸代谢、TGF-β信号通路以及Hippo信号通路;Am5 vs Am6中的上调和下调mi RNA可分别结合43和431个靶基因,其中mi R-iab-4-x结合靶基因数量最多,并广泛参与调控背腹轴的形成、Hippo信号通路、Wnt信号通路、Fox O信号通路、Notch信号通路以及m TOR信号通路等与生长发育相关的通路。调控网络分析结果表明,DEmi RNA与靶基因间形成较为复杂的调控网络,DEmi RNA居于调控网络的中心位置,而m RNA位于调控网络的外周。最后,通过茎环实时荧光定量PCR对随机选取的3个DEmiRNA进行验证,结果证实了测序数据的可靠性。【结论】在全基因组水平对意蜂幼虫肠道的DEmi RNA及其靶基因进行预测和分析,并对DEmi RNA的调控网络进行构建及分析,发现意蜂幼虫通过调节包括ame-mi R-6052、miR-iab-4-x、mi R-281-x和novel-m0031-3p在内的多个mi RNA的表达水平对肠道生长和发育进行调控,研究结果不仅提供了意蜂肠道发育过程的mi RNA表达谱及差异表达信息,也为阐明意蜂幼虫肠道的发育机理打下基础。

应用Gateway重组技术构建CSBV非结构蛋白,RNA依赖的RNA聚合酶(RNA-dependent RNA polymerase,RdRp)的原核表达载体,将重组表达载体转化大肠杆菌Rosetta,经异丙基硫代-β-D-半乳糖苷(Isopropylβ-D-1-Thiogalactopyranoside,IPTG)诱导获得RdRp基因原核表达蛋白.以重组蛋白为抗原免疫新西兰兔制备相应的多克隆抗体,Western blot检测表明,该抗体能与RdRp重组蛋白和感病中华蜜蜂幼虫蛋白特异性结合,说明获得的抗体特异性高.利用制备的抗体进行免疫荧光标记,发现其能特异地定位在感病中华蜜蜂血淋巴细胞内.

从越夏期意大利蜜蜂(Apis mellifera ligustica)工蜂成虫肠道中分离、纯化、培养,获得9株细菌.通过形态学观察和生理生化特性研究,从9株细菌中鉴定出芽孢杆菌属(Bacillus)、柠檬酸杆菌属(Citrobacter)、沙雷氏菌属(Serratia)、泛生菌属(Pantoea)、肠杆菌属(Enterobacter).

【目的】研究眼斑芫菁幼虫肠道共生菌的多样性，并鉴定其相关功能，以探寻眼斑芫菁肠道共生菌在其生长发育中的作用。【方法】分别以人工饲料、蝗卵为食物饲喂眼斑芫菁幼虫，饲养至不同龄期的眼斑芫菁为材料，按相应龄期收集试验所需幼虫，提取眼斑芫菁幼虫肠道细菌基因组DNA，通过16S rDNA全长扩增子测序，分析其肠道细菌的多样性和相对丰度，并探究其肠道共生菌发挥的潜在功能。【结果】获得最终有效序列81 759条，聚类共获得230个OTUs，注释到17个门，27个纲，41个目，71个科，109个属。变形菌门（Proteobacteria）为3个组中的最优势门，各分组样品的优势菌属不相同，劳尔氏菌属为3个组的共同优势菌属。α多样性分析发现，蝗卵饲养的眼斑芫菁二龄幼虫（L.H2）肠道菌群多样性和丰度较四龄幼虫（L.H4）高，且均高于转接入人工饲料饲养的四龄幼虫（L.S4）；α与β多样性分析结果均显示眼斑芫菁肠道内共生菌多样性各组内差异大，组间有一定差异，但不同组之间差异不显著（P＞0.05）。功能分析结果发现，L.S4在氨基酸代谢、核苷酸代谢、环境适应、免疫系统等方面明显高于L.H4和L.H2，且L.S4的氨基酸代谢比L.H4显著增强（P=0.035）。【结论】明确了眼斑芫菁幼虫肠道菌多样性和群落结构受食物影响，其肠道细菌在代谢功能上能协助宿主适应相应的食物。研究眼斑芫菁幼虫肠道细菌多样性与食物之间的关系，为研究眼斑芫菁肠道微生物变化奠定了基础。