红脚艾蒿的转录组解析

[目的]利用高通量测序技术解析红脚艾(Artemisia rubripes Nakai)的转录组信息特征.[方法]通过高通量测序平台Illumina HiSeq 2500对红脚艾进行转录组测序,通过Trinity软件de novo组装获得Unigene,并基于序列同源性对Unigene进行功能注释,得到红脚艾的转录组信息.[结果]测序数据经过质控后共获得24126043条高质量的reads,通过de novo组装获得173093个转录本,对组装的转录本去冗余后共获得85991个Unigene,平均长度为616.87 bp,N50为925 bp.共有47216个Unigene在NR、KEGG、COG、KOG、GO数据库获得功能注释,40802个Unigene在NR数据库注释,显示红脚艾与向日葵(Helianthus annuus)的单基因匹配率最高,16846个Unigene被KEGG数据库注释到130条代谢途径中,26171个Unigene被注释到25个KOG功能分类中,23203个Unigene被GO注释到生物过程、细胞组成和分子功能三大类51个功能分类,12810个Unigene被注释到25个COG功能分类中.[结论]利用高通量测序技术获得了红脚艾转录组信息特征,这些数据将为后期开展功能基因鉴定、解析化合物次生代谢途径及其调控机制奠定研究基础.     

香榧转录组测序及生物信息学基础分析

香榧具有重要的经济价值,但其基因组信息相对匮乏,限制了其分子生物学和基因功能的研究。本文以不同组织的香榧作为研究对象,采用新一代高通量测序技术平台Illumina Hi Seq?2000对香榧转录组进行测序和数据分析,共得到37,349,086个reads片段,总碱基数为4.35 G。利用组装软件,对获得的高质量序列进行组装,共得到104,636个Unigene,平均长度为784 nt,N50为1,702。将Unigene序列与公共数据库进行比对,28,766个Unigenes获得了注释。其中26,856个Unigene在NR蛋白数据库中获得注释,24,003个Unigenes在NT数据库中获得注释,21,401个Unigene在Swiss-Prot蛋白数据库中获得注释,16,137个Unigene在COG数据库中获得注释,11,410个Unigene在GO数据库中获得注释。根据KEGG注释信息,18,564个Unigene被划分到256个代谢途径中。SSR位点搜索发现,在4,217个Unigene中含有4,706个SSR位点。分析所获得的转录组数据,将为香榧功能基因的克隆,基因的表达,指纹图谱构建和分子标记辅助选育奠定基础。     

荸荠优良品种‘桂蹄3号’球茎发育过程转录组测序分析

【目的】开展荸荠球茎发育过程转录组测序研究,为研究荸荠球茎发育过程相关基因表达信息提供参考。【方法】运用高通量测序技术,对荸荠球茎不同发育时期进行转录组测序研究。【结果】组装共得到223 182条转录本和90 542条Unigene,平均长度为809 bp,N50为1119。组装完整性较高,效果较好。对所得Unigene进行不同数据库注释,共有50 583条Unigene成功注释到7个数据库(NR、GO、NT、Pfam、KEGG、KOG及Swiss-prot)。对差异表达基因分析,发现球茎膨大初期的差异表达基因数目最多,为最活跃阶段。GO功能富集分析结果表明,33 205个基因获得功能注释,分为分子功能、细胞组分和生物学过程等3大类和54个亚类。COG功能分类结果表明,17 743个基因分布于25个功能区域,其中碳水化合物代谢占重要地位。KEGG代谢通路注释结果表明有20 667个基因获得功能注释,共有116条代谢途径,其中淀粉-蔗糖代谢占主要作用。【结论】利用高通量转录组测序技术首次建立了荸荠优良品种‘桂蹄3号’球茎的转录组数据库,为进一步研究荸荠球茎淀粉生物合成相关基因的功能及形成的分子机制提供了数据基础。     

黄芩高通量转录组测序数据组装和分析

采用高通量测序技术对黄芩的转录组进行测序,获得了29 099 899条reads数据,拼接后得到53 353条Unigene;将所获得的Unigene与COG,GO,KEGG,Swiss-Prot,NR这5个公共数据库进行比对,结果发现,分别有10 756,21 950,8 101,20 339,29 288条Unigene可比对到以上5个数据库中;已注释的Unigene与COG数据库比对后按功能可分为25类;根据GO功能可分为三大类57个亚类;经过与KEGG数据库比对后按照代谢通路可分为116类;利用Get ORF软件进行ORF预测,获得长度大于300 nt的ORF共20 552个;通过SSR分析,共获得5 658个SSR标记。获得的转录组信息可为今后进行黄芩分子标记的开发和关键基因的克隆及功能分析等研究提供基础数据。     

大竹蛏高通量转录组测序数据组装和分析

为进一步开发大竹蛏Solen grandis的基因资源,采用2代Illumina Hi-seq测序技术对大竹蛏的鳃组织进行了转录组测序,构建了转录组数据库,获得338 483 476条Clean Reads数据;拼接组装后获得190 856条Unigene数据,平均长度为1147 bp;与NR、NT、KO、Swiss Prot、PFAM、GO、KOG等数据库进行Blast信息比对(E-value为10-5),共获得63 337个注释基因;与NR数据库比对发现,大竹蛏转录组基因序列与长牡蛎Crassostrea gigas具有较高的同源性,为53.3%;将大竹蛏转录组的Unigene的功能通过与KOG数据库进行注释比对划分为25类;GO数据库注释可分为三类,即细胞组分、生物过程和分子功能,共包括65个分支;KEGG分析发现,大竹蛏转录组数据中按照代谢通路可分为92类,利用Blast蛋白库比对和Estscan软件进行ORF预测,获得长度大于300 nt的ORF共50 681个;通过SSR分析,共获得73 089个SSR标记。本研究中获得的转录组信息可为今后进行大竹蛏分子标记的开发和关键基因的克隆及功能分析等研究提供基础数据。     

红螯螯虾转录组高通量测序及分析

【目的】为发掘和研究红螯螯虾功能基因SSR标记奠定基础。【方法】对红螯螯虾肝脏、精巢以及卵巢转录组测序选用新一代高通量测序技术Illumina Hi Seq 2000,同时借助生物信息学法展开基因表达谱探究及功能基因预测。【结果】肝脏组织获得47 629 414条clean reads,精巢组织获得47 018 968条clean reads,卵巢组织获得53 267 362条clean reads。所有reads组装后获得了67 369个Unigene。通过GO分类,16 989个Unigene被分为3个种类,即生物学过程、细胞组分、分子功能。通过COG分类,4 697个Unigene被分为25个种类。通过KEGG分类,9 842个Unigene分属331个代谢途径。【结论】通过对红螯鳌虾肝脏等3个组织转录组的测序,获得了其相关基因的表达图谱及其功能分类。     

青篱柴转录组的高通量测序及分析

青篱柴具有很好的观赏价值和药用价值,但其基因组信息相对缺乏,导致其分子生物学和基因功能的研究受到限制.为了获得和挖掘青篱柴的基因数据和功能,首次利用新一代高通量测序技术平台Illumina HiSeqTM2000对青篱柴转录组进行测序,经de novo组装后共得到66 755条单基因簇(Unigene).进一步利用7大公共数据库进行Blast同源比对,注释了50 339条Unigene.研究发现,558个基因参与了青篱柴次生代谢产物的生物合成,其中在甜菜红碱、芥子油苷、类黄酮、单菌霉素和苯丙素生物合成途径中的Unigene分别有2个、18个、44个、70个和231个,这些基因很可能参与了青篱柴药用活性物质的生物合成.SSR位点搜索发现,在18 972条Unigene中含有22 542个SSR位点.研究结果将为青篱柴功能基因的克隆、基因的表达、分子标记开发和遗传多样性研究奠定基础.同时,转录因子分析结果也为青篱柴的喀斯特环境适应性机制提供了初步线索.     

泸定百合转录组测序与特性分析

泸定百合是非常优良的种质基因资源,鳞茎研究是其科研和生产的关键环节,然而泸定百合鳞茎的分子生物学研究尚属空白。为揭示泸定百合鳞茎的分子机制,应用高通量测序技术对泸定百合的鳞茎进行转录组分析。原始数据经过生物信息学分析后,共获得4.9G有效数据,43 412条Unigene;所获得的Unigene与NR,GO,COG和KEGG等数据库进行搜索比对后,发现有30 225条和21 531条Unigene分别获得NR和GO注释,14 069条Unigene得到COG注释,GO注释中获得374条Unigene可能作为转录因子参与泸定百合鳞茎的形成与发育;共有10 822条Unigene参与了KEGG代谢途径,其中有571条Unigene注释到次生代谢合成途径,其中66条和23条Unigene分别参与了萜类和生物碱的合成。研究结果为采用生物技术手段提高百合的药用成分的含量及改善其园艺性状积累了基本资料。     

杜仲雌雄株转录组测序数据组装及基因功能注释

以10年以上树龄的杜仲雌株当年新发枝条上的幼果、嫩芽、叶片和树皮和雄株新发枝条上嫩芽、叶片和树皮为材料,采用Illumina Hi SeqTM2000高通量测序技术进行转录组测序,获得雌株51,574,000条、雄株52,430,502条Clean Reads数据,分别包含总长度为4,641,660,000nt和4,718,745,180nt核苷酸序列数据信息;经拼接组装,获得雌株基因信息长达69,461,730nt的423,339个Contig片段,获得雄株基因信息长达94,814,201nt的542,383个Contig片段;经进一步拼接,分别获得平均长度为288nt的雌株159,434个Unigene片段和平均长度为231nt的雄株257,288个Unigene片段,共有48,761个表达序列标签(EST)。以BLAST(E-value1.0E-5)将Unigene对NR、NT、KEGG和COG数据库进行比对,获得CDS序列35,541条,再通过ESTscan分析获得CDS片段13,220条,共获得48,761条CDS片段。与NR数据库比对发现杜仲雌、雄株转录组Unigene与葡萄相似序列最多(33.8%),其次是蓖麻(11.4%)和杨树(11.2%),与拟南芥的相似序列仅2.3%;根据Unigene与COG数据库比对结果,可将有COG功能的7,571条Unigene分为24类,而根据GO数据库注释,杜仲转录组有GO功能注释的23,314条Unigene可分为生物过程、细胞组分和分子功能3大类55分支。与KEGG数据库比对,杜仲雌、雄株转录组17,468条Ungenes分属128类代谢通路,其中有2,399条属于次生物质代谢途径,314条参与萜类化合物生物合成途径。     

多油辣木转录组高通量测序及分析

为了研究多油辣木(Moringa oleifera Lam.)中黄酮类化合物生物合成的分子基础,通过Illumina Hi Seq 2000高通量测序技术对辣木茎、叶进行转录组测序,利用Trinity软件将数据组装成Unigene,基于BLAST对所有Unigene进行功能注释,共获得49 365个Unigene,平均长度为903 bp。通过GO分类,15 959个Unigene被分成生物学过程、细胞组分和分子功能3个主要类别。通过KOG分类,8 713个Unigene被分为26个种类。通过KEGG分类,8 397个Unigene分属于130个代谢途径,其中代谢所含Unigene最多,共3 620个。在代谢途径中,找到参与黄酮类化合物合成相关的Unigene 45个,其中包括查尔酮合酶、查尔酮异构酶、黄烷酮-3-羟化酶、黄酮醇合成酶和二氢黄酮醇还原酶等。研究结果为挖掘多油辣木黄酮类化合物生物合成关键基因提供了基础数据,并为下一步的资源开发和利用奠定了基础。     

茄子温敏单性结实抑制差减文库的构建与分析

 【目的】构建茄子单性结实差减文库，分离获得茄子单性结实相关基因，研究茄子单性结实形成的分子机制。【方法】以茄子单性结实品系D-10低温条件下无籽的子房和果实，以及适温条件下有籽的子房和果实为材料，采用抑制差减杂交（SSH）技术构建正向和反向差减文库。对阳性克隆测序、BLAST比对分析、GO（gene ontology）功能注释及分类和荧光定量PCR分析。【结果】正反向文库共获得了2 347个阳性克隆；测序、序列拼接后，共获得1 248个高质量的unigenes；将其与非冗余数据库BLAST比对后，有1 109个unigenes找到同源序列，139个unigenes无同源序列为新基因。1 248个unigenes中527个在GO数据库中有功能注释，其中细胞组分、分子功能和生物过程的unigenes分别为206、445和361条，分析表明，茄子单性结实果实的差异主要发生在细胞内部，单性结实性状的表达可能受到一些因子和激酶的调控。【结论】成功构建了茄子单性结实不同时期果实的抑制差减文库，获得3个可能与茄子单性结实相关的unigenes，包括MADS-box基因、雌蕊伸展相关蛋白和果实膨大相关基因。     

氯虫苯甲酰胺干扰桃小食心虫交配的转录组分析

【目的】探明氯虫苯甲酰胺处理后桃小食心虫(Carposina sasakii)成虫的转录组差异,了解该药剂影响桃小食心虫交配的基因在功能分类和代谢通路等方面的生物学特征,挖掘与交配相关的功能基因。【方法】通过生物学实验观察氯虫苯甲酰胺干扰桃小食心虫成虫交配及繁殖情况。采用Illumina Hi SeqTM2500高通量测序技术对刚羽化的桃小食心虫雌雄成虫、羽化后4—6 h进入交配高峰期的雌雄成虫和经氯虫苯甲酰胺处理4—6 h的雌雄成虫进行转录组测序,利用Trinity软件对所得序列进行de novo组装及评估,之后对获得的有效序列进行功能注释,并利用q RT-PCR技术分析氯虫苯甲酰胺处理后相关基因的时空表达变化。【结果】氯虫苯甲酰胺处理后,桃小食心虫的交配率显著降低,寿命缩短,产卵量减少。通过合并组装桃小食心虫转录组有效序列共获得102 831条unigene,其中34 526个有注释信息。根据筛选标准,氯虫苯甲酰胺处理过程中,雌雄虫中分别有122个和147个基因发生变化,其中相同差异基因31个。对所存在的234个差异基因进行GO功能注释和富集分析结果显示,分子功能过程中的催化活性和结合活性及生物学过程中与代谢过程、单一生物体过程和细胞过程相关的5类基因占主导地位。KEGG分类结果显示,富集到代谢通路的最多,有25个,包括昆虫激素合成、药物代谢等。通过比对分析,鉴定c64662.graph_c0为桃小食心虫鱼尼丁受体基因,其长度为15 637 bp,与已报道Cs Ry R的一致性为99.0%。此外,在234个差异基因中,鉴别羧酸酯酶unigene 3个、细胞色素P450 unigene 4个、肌钙蛋白unigene3个、气味结合蛋白unigene 1个和生物钟unigene 1个。细胞色素P450 c40709.graph_c0参与昆虫激素生物合成。根据转录组中基因表达分析,12个基因在雌雄虫中的表达均出现不同变化趋势。q RT-PCR结果显示,氯虫苯甲酰胺诱导羧酸酯酶c51998.graph_c0基因上调表达;药剂处理后,雌雄虫的c57480.graph_c0和c53794.graph_c0的表达分别呈现显著的上调和下调变化,而c40709.graph_c0基因仅在雄虫中显著下调,处理6 h后c53281.graph_c0只在雌虫中上调表达;氯虫苯甲酰胺处理后3个肌钙蛋白基因在整个试验阶段均表现明显的下调趋势;雄虫中的生物钟unigene c60883.graph_c0和气味结合蛋白unigene c45675.graph_c0的变化趋势较为一致,进入暗期后立即上调,但均受氯虫苯甲酰胺的抑制。雌虫体内Ry R的表达量显著上调,而雄虫初次到达求偶高峰期前Ry R的表达量与对照差异不显著,到第2个求偶高峰期时,Ry R表达显著下调。除细胞色素P450c57480.graph_c0、c40709.graph_c0和c53281.graph_c0外,其余基因在雄虫中的表达均高于雌虫。且触角酯酶基因c54944.graph_c0、生物钟基因c60883.graph_c0和气味结合蛋白基因c45675.graph_c0在黑暗光照交替变化时,表达发生明显地上调或下调。【结论】通过转录组测序发现氯虫苯甲酰胺干扰桃小交配的作用机制是由靶标基因、嗅觉相关基因、代谢基因、生物钟基因等相互作用引起的。     

稻瘟病菌菌株FJ81278中新发现的一个基因注释及其表达分析

对稻瘟病菌菌株FJ81278进行基因组和转录组分析,发现一个新的特异基因MoHsbA编码的蛋白具有疏水表面结合蛋白A(HsbA)结构域和一个信号肽.进一步的生物信息学分析发现,该基因组中另外7个基因编码的蛋白也具疏水表面结合蛋白结构域.稻瘟病菌中的HsbA蛋白聚集于系统发育树的同一个大的分支上,说明它们具有较近的亲缘关...     

基于银鲳RNA－seq数据中SR标记的信息分析

[目的]开发银鲳分子标记技术。[方法]通过对银鲳（Pampus argenteus）进行高通量转录组测序（RNA－seq ）获得银鲳转录组原始试验数据，经过拼接后，获得3715603条unigene序列。采用生物信息学分析软件MISA对所有unigene进行简单重复序列（ simple se-quence repeat，SSR）位点鉴定。同时，利用软件对银鲳转录组SSR的多态性进行评价。[结果]银鲳转录组水平上，共鉴定出107007个SSR位点，分布在97289条unigene中，发生频率为2．62％，SSR平均密度为476个／Mbp。在银鲳转录组SSRs中，单核苷酸与二核苷酸重复序列为主要重复类型，分别占总SSRs的48．14％和34．10％。银鲳转录组数据中SSR序列共包括424种重复基元类型，单核苷酸重复基元A占较高比例，占同一重复类型SSRs的49．57％，二核苷酸重复基元TG／AC和三核苷酸重复基元GAG／AAC是优势重复基元，分别占同一重复类型SSRs的42．43％和9．61％。重复序列长度在12 bp以上的SSR标记位点数占总SSR的76．95％，具有丰富的多态性。[结论]银鲳SSRs位点具有极大的可开发性，可为银鲳遗传多样性研究、遗传图谱构建及分子辅助育种提供有效工具。     

黄瓜抗霜霉病相关EST-SSR的信息分析

对霜霉病菌诱导的黄瓜cDNA文库测序获得的3091条ESTs进行分析,经过除前处理和拼接后得到2 507条Unigenes.应用生物信息学方法,对黄瓜抗霜霉病相关EST-SSRs的分布频率和碱基重复特点进行了归纳及分析.结果表明,2 507条Unigenes中共检索到231个EST-SSR,出现频率为7.47％,平均分...     

基于高通量测序的杂交鳢转录组数据分析

运用高通量技术对杂交鳢进行转录本测序和数据分析,经拼接与组装,最终获得52 065条unigenes,长度范围201～12 407 bp,序列平均长度为710 bp。从长度分布、GC含量和基因表达量等方面对unigenes进行评估,数据显示测序质量较好。进一步的数据分析,共鉴定出20 535个CDS,37 324个SNP位点和7 830个微卫星。用6大数据库Swiss–Prot、Nr、Pfam、KEGG、KOG和GO注释杂交鳢转录组unigenes,分别对应有20 955、28 938、19 339、14 052、19 358和18 635条unigenes获得注释。根据GO数据库的注释,可将这些unigenes分为生物过程、细胞组分和分子功能3大类共50亚类。KEGG代谢通路分析表明,这些unigenes参与了5大类30亚类共268个代谢通路,其中,参与信号转导通路的unigenes数量最多,有1 716条。     

中华蜜蜂幼虫肠道参考转录组的de novo组装及SSR分子标记鉴定

【目的】利用RNA seq技术对中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫肠道参考转录组进行de novo组装,并进行功能及代谢通路注释,进而利用该转录组数据进行中蜂幼虫的SSR分子标记鉴定。【方法】实验室条件下饲养中蜂幼虫,将纯化的蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)孢子饲喂3日龄幼虫,剖取4、5和6日龄幼虫肠道,液氮速冻。将健康幼虫肠道与感染球囊菌的幼虫肠道同时进行Illumina测序。通过对raw reads的过滤得到clean reads,利用Trinity软件组装得到unigenes。通过BLASTx(E-value10-5)比对NCBI Nr、Swiss-Prot、KOG和KEGG数据库,对unigenes进行功能和代谢通路注释。利用MISA软件对所有unigenes进行SSR搜索,并利用Primer Premier 5软件设计特异性SSR引物,通过常规PCR对来源于北京、辽宁兴城和四川成都的中蜂幼虫肠道样本进行SSR位点鉴定。【结果】中蜂幼虫肠道的RNA seq共得到35 670 000条reads,de novo组装得到43 557个unigenes,平均长度为898 nt。共有18 225个unigenes可被注释到上述公共蛋白数据库,单独注释到NCBI Nr、Swiss-Prot、KOG和KEGG数据库的unigenes数分别为3 899、443、37和10个。KOG注释结果显示,11 442条unigenes分布于25个基因家族,其中注释到RNA加工和修饰家族的基因数最多,达1 249个。9 679个unigenes的GO分类结果显示,在生物学进程分类中,注释到细胞进程的基因最多,达4 201个,在分子功能和细胞组分类中,注释到结合与细胞的基因数最多,分别为4 935和2 900个。4 517个unigenes可注释到KEGG数据库中的216个代谢通路,注释到核糖体的基因数最多,达385个。利用MISA软件,可在7 763个unigenes搜索到13 448个SSR位点,随机选取20对SSR引物对国内3个不同来源的中蜂幼虫肠道样本的SSR位点进行扩增,有6对引物可鉴定出SSR分子标记。【结论】成功组装并注释了中蜂幼虫肠道参考转录组,可为中蜂及其幼虫的分子生物学及组学研究提供重要的参考信息,也可用于补充、丰富和检验东方蜜蜂的参考基因组,基于此转录组数据开发出6个中蜂的SSR分子标记,可应用于中蜂的基因图谱构建、基因多样性分析、基因定位等研究,也说明利用转录组数据开发非模式生物SSRs的方法可行。     

乙烯诱导下柑橘果实离层SSH文库的构建及分析

为发掘柑橘果实脱落过程中的关键基因,构建了受乙烯诱导的柑橘果实离层SSH文库,经蓝、白斑筛选,菌液PCR检测和测序验证,共获得385条有效EST.经聚类拼接得到单基因簇30个,其中片段重叠群28个、单一EST序列2个.经BLASTX比对NCBI的非冗余蛋白数据库,在这些ESTs中,28个找到同源序列,2个无匹配;BLAST2GO功能注释表明:这些单一基因主要涉及衰老、抗逆防御、细胞壁水解、信号转导等功能.     

青花菜花蕾转录组测序分析及蜡粉合成相关基因挖掘

【目的】对青花菜花蕾进行转录组测序分析，并挖掘与蜡粉合成相关基因，为探明青花菜花球表面蜡粉形成的分子机制提供理论参考。【方法】分别提取野生型和蜡粉缺失型青花菜花球总RNA，采用Illumina HiSeqTM 2500平台进行转录组测序，获得高质量Clean reads，采用Trinity进行序列组装后获得青花菜Unigene库，将获得的Unigene序列与Nr、Nt、KEGG、Pfam、KOG/COG、Swiss-Prot和GO数据库比对，获得基因功能注释信息；使用DESeq2进行差异表达分析。【结果】共获得44.68 Gb Clean data，De novo组装得到41244条Unigenes，N50长度为1847 bp。从所获得的Unigenes中筛选出8685个差异表达基因（DEGs）（上调基因5747个，下调基因2938个），共有8038个基因被注释到不同数据库，其中，5220个基因注释到Pfam数据库； 2066个基因注释到COG数据库，3866个基因注释到KOG数据库； 2580个差异表达基因被注释到75个转录因子家族中，注释最多的是MYB家族（235个）；GO数据库中6095个差异表达基因注释到细胞组分、分子功能和生物学过程三大类的52个功能分类； KEGG数据库中，1671个差异表达基因富集到138条代谢通路，其中13个差异表达基因与脂肪酸合成有关，7个差异表达基因与蜡粉生物合成途径有关。【结论】转录因子MYB家族在调控青花菜蜡粉合成中发挥重要作用。蜡粉合成过程中相关酶基因的差异表达是调控青花菜蜡粉合成的关键，尤其是野生型和蜡粉缺失突变体中特异性表达的差异表达基因，可作为后续研究青花菜花球表面蜡粉形成分子机制的对象。     

胁迫中华蜜蜂幼虫肠道的球囊菌及其体外培养的高表达基因分析

利用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)和细胞进程(1748 unigenes);对照组(Aa CK)的HEGs富集于37个GO term,基因富集数最多的是细胞(785 unigenes),其次是细胞组件(785 unigenes)和代谢进程(776 unigenes).KEGG pathway富集分析显示,Aac T的HEGs富集在119个代谢通路上,基因富集数最多的是核糖体(176 unigenes)、碳代谢(147 unigenes)及氨基酸的生物合成(134 unigenes);Aa CK的HEGs富集在110个代谢通路上,基因富集数最多的是核糖体(179 unigenes)、氨基酸的生物合成(70 unigenes)以及碳代谢(62 unigenes).深入分析发现,对于富集在MAPK信号通路上的高表达基因数量,胁迫中蜂幼虫肠道的球囊菌远多于体外培养的球囊菌,说明病原的该通路在胁迫后期被显著激活.