鸭肠炎病毒感染对鸭肠道黏膜组织转录组变化分析

为了探究鸭肠炎病毒感染对鸭十二指肠黏膜组织相关免疫基因的变化影响,本试验采用鸭肠炎病毒接种35日龄的健康樱桃谷鸭,将其分为3组(2个试验组和1个对照组),在感染后24 h、48 h采集十二指肠,提取总RNA进行浓度、纯度及完整性检测后构建cDNA文库,使用BGISEQ-500平台对2个实验组和对照组黏膜组织测序,筛选差异表达基因,并用GO和KEGG数据库进行分析。结果显示,2个实验组感染鸭肠炎病毒后T24和T48分别筛选出共720个差异表达基因,其中T24表达上调58个,下调163个;T48表达上调77个,表达下调422个。GO功能注释结果显示,感染后差异基因都与代谢、免疫、炎症和调控途径等密切相关。KEGG通路富集分析发现,感染T24的差异表达基因富集到神经活性配体-受体相互作用、补体和凝血级联、氮代谢、cAMP信号通路和IL-17信号转导通路;感染T48的主要差异表达基因富集到神经活性配体-受体相互作用、细胞因子-细胞因子-受体相互作用、趋化因子信号通路、造血细胞谱系和补体和凝血级联。借助BGISEQ-500平台对鸭肠炎病毒感染后十二指肠黏膜组织相关因子的变化来探究其致病机制及其他生物学通路,弥补了非模式生物转录组研究中缺乏基因信息的不足。     

基于高通量测序技术分析番鸭呼肠孤病毒感染雏番鸭肠道转录组学的变化

番鸭呼肠孤病毒(Muscovy duck reovirus,MDRV)感染雏番鸭可导致肠道黏膜受损,造成肠黏膜免疫功能低下甚至丧失。为寻找番鸭呼肠孤病毒感染雏番鸭肠道组织免疫相关的差异表达基因,本研究建立番鸭呼肠孤病毒自然发病模型,在感染后3d、6d、21d采集同居感染组和对照组试验雏番鸭的空肠,对其进行高通量测序,对差异表达基因进行GO功能分类和KEGG数据库分析。结果显示,感染后3d、6d和21d分别筛选出2 297、5 860和3 721个差异表达基因;GO分子功能结果显示,免疫相关差异表达基因主要和免疫球蛋白的产生、T细胞活化和单核细胞趋化性等有关;KEGG通路富集结果显示,差异表达基因主要涉及在吞噬体、细胞溶质DNA传感途径、细胞因子-细胞因子受体相互作用、Toll样受体信号通路、Jak-STAT信号通路和RIG-I受体信号通路;荧光定量RT-PCR对差异表达基因进行验证,其结果与转录组学结果基本一致。本研究为进一步阐明MDRV感染导致雏番鸭肠道损伤和黏膜免疫抑制的机制奠定基础。     

基于高通量测序技术分析家兔感染猪瘟兔化弱毒株后脾脏的转录组学变化

本研究基于高通量测序技术对感染猪瘟兔化弱毒株(C株)后0h,12h,24h,30h,36h,48h的家兔脾脏组织进行测序,以家兔基因组作为参考,筛选感染后各时间点和对照组(0h)之间的差异基因。利用Uniprot和NCBI数据库查找各组中变化水平显著的前10个差异基因的生物学功能。结果发现,感染后12h,24h,30h的前10个差异基因只有少数与炎症、免疫、凋亡等相关;感染后36h,48h的前10个差异基因完全相同,其中B2M、RLA-DR-ALPHA、CD74和IGJ参与抗病毒过程。GO功能注释结果显示:感染后各时间点的差异基因都和免疫、代谢、调控途径紧密相关。KEGG通路富集分析发现,感染后24h的差异基因富集到RIG-I样受体信号通路,感染后30h的差异基因富集到黏着斑和细胞外基质-受体相互作用通路,感染后36h和48h的差异基因共同富集到的通路有20个,其中蛋白酶体,溶酶体,核糖体,趋化因子信号通路,B细胞受体信号通路,抗原加工提呈通路和FcγR介导的吞噬通路与抗病毒相关。本研究建立的家兔感染C株后脾脏组织的差异表达基因数据库,将为进一步阐述家兔适应C株的分子机制奠定基础。     

京海黄鸡柔嫩艾美耳球虫感染后盲肠转录组分析

为了探究京海黄鸡感染柔嫩艾美尔球虫(E. tenella)后盲肠组织差异表达基因,以及球虫感染分子应答过程和免疫应答机制,试验采用RNA-seq技术对E. tenella感染和非感染组第7天的盲肠组织进行转录组测序,筛选差异表达基因,并进行差异基因的功能、通路富集分析。结果表明,在感染和非感染组中有显著差异的表达基因2 830个(P0.05),其中1 419个基因上调,1 411个基因下调。随机选取10个差异基因进行qRT-PCR验证,结果显示差异基因的表达倍数与RNA-seq检测结果显著相关(r=0.988,P0.000),决定系数达0.975。GO分析表明,有2 356个差异基因获得GO功能注释,显著富集的前30个GO terms主要涉及细胞交流、信号转导、血管生成、氧化还原酶活性等。KEGG分析发现差异基因显著富集的信号通路有黏着斑、细胞外基质-受体相互作用、过氧化物酶体增殖物激活受体等。这些通路中的差异基因有ANGPTL4、ACSL5、VEGFC、CD44和MAKP10等,提示这些基因在宿主柔嫩艾美耳球虫感染过程中发挥重要作用。     

鸭肠炎病毒核衣壳蛋白受体的筛选与鉴定

为进一步阐明鸭肠炎病毒(DEV)致病机理,本研究通过构建DEV感染鸭肝组织cDNA文库及其文库质粒和DEV核衣壳蛋白(NP)基因"诱饵"质粒,利用酵母双杂交系统筛选DEV NP互作蛋白(受体)基因,并采用GSTpull-down试验进行验证,结果显示:所构建鸭肝组织cDNA文库的库容量为1×106 CFU,插入片段大小集中在0.5～1kb;"诱饵"质粒pGBKT7-NP无自激活现象;筛选并初步证实DEV NP在鸭肝组织细胞的受体蛋白为蛋白激酶C抑制蛋白(PKCI)。这些结果为进一步研究DEV致病机理提供新的线索。     

基于RNA-seq研究胶质母细胞瘤中差异表达基因及蛋白相互作用网络

本研究基于RNA-seq数据分析了胶质母细胞瘤中的差异表达基因,并对差异表达基因进行了功能(GO term)和通路(KEGG)富集分析。进一步通过蛋白相互作用网络挖掘了胶质母细胞瘤的调控机理。结果表明,405个基因在肿瘤组织中差异表达(p-value≤0.05,|log FC|≥1.5),其中216个(53.3%)基因上调,189个(46.6%)基因下调。基因本体(gene ontology,GO)富集结果表明,这些差异表达基因参与了离子转运,神经冲动传递,细胞信号转导和细胞粘附等。此外,KEGG通路富集结果表明,差异基因参与了许多重要的生物学通路,包括ECM受体相互作用、黏着和钙信号等通路。进一步的蛋白相互作用网络分析鉴定了5个关键的hub基因,包括PTK2B、CDK1、FN1、CCND1和FOS。这5个关键基因对于胶质母细胞瘤的发生和发展可能起到了关键作用,可以作为潜在的调控位点和筛选的标志物。     

基因芯片筛选新疆维吾尔妇女宫颈癌差异表达基因

新疆南部维吾尔族聚居区是宫颈癌高发区. 本文旨在利用基因芯片技术筛选与维吾尔族妇女宫颈癌发生相关的基因. 首先,分别提取5例新疆维吾尔妇女宫颈癌和5例子宫肌瘤组织(对照)的mRNA,逆转录成cDNA,并用Cy3-dUTP标记子宫肌瘤组织的cDNA,用Cy5 dUTP标记宫颈癌组织的cDNA,制成芯片杂交探针.为筛选出宫颈癌组织中差异表达的基因,上述标记探针分别与含有20 000条人类基因的Affymetrix基因芯片进行杂交,杂交信号用GeneChip Scanner 3000扫描仪扫描,并用芯片图像分析软件(SAM software)分析扫描结果.筛选出的差异表达基因经GO(Gene Ontology)分析和KEGG（Kyoto Encyclopedia of Genes and Genomes）信号通路分析,确定其在宫颈癌中的作用.基因芯片筛选结果显示,在宫颈癌组织中发现2 758个差异表达基因,其中1 326个上调基因,1 432个下调基因.GO分析和KEGG信号通路分析表明,表达差异在两倍以上的基因涉及168个信号通路,包括细胞粘附分子、细胞周期以及MAPK和mTOR信号通路等.上述结果表明,基因芯片技术筛选出大量与宫颈癌发生相关的基因,其中表达差异显著的基因涉及细胞粘附分子、细胞周期和mTOR等信号通路.     

HBV感染者肝细胞样本的关键基因和信号通路的生物信息学挖掘

为探究乙型肝炎病毒(Hepatitis B virus,HBV)感染后肝细胞基因表达和信号通路的改变情况,从Gene Expression Omnibus(GEO)数据库中下载了HBV感染者肝细胞样本制作的基因表达谱数据集GSE83148,进行质量检测、数据标准化后筛选出差异表达基因,进一步做GO和KEGG富集分析以及基因网络相互作用分析,筛选关键基因和信号通路。从HBV感染样本中筛选出fold change≥2,p-value0. 05的上调差异表达基因44个,GO分析获得关键基因BAK1和TP63,差异表达基因网络互作获得5个位于枢纽位置的基因:NDUFS1、NDUFS2、COX7B、ATP5B、OPA1。KEGG分析获得关键信号通路有:乙型肝炎信号通路、病毒癌变信号通路、Fox O信号通路、PI3K-Akt信号通路。本研究筛选出的多数基因与线粒体和氧化呼吸链有关,造成这一现象的具体机制还需进一步探究。     

未成熟树突状细胞与成熟树突状细胞的细胞骨架调控相关基因表达的差异分析

分析未成熟树突状细胞(immature dendritic cells,im DCs)向成熟的树突状细胞(mature dendritic cells,m DCs)分化的过程中,细胞骨架的调控相关基因及信号通路的表达变化,为进一步理解不同分化阶段树突状细胞(dendritic cells,DCs)的生物物理学特性、迁移能力和免疫学相关功能的改变。利用CEO数据库获得经CD14+单核细胞(monocytos,monos)诱导而成的im DCs和m DCs的m RNA的表达数据(芯片编号:GSE15076),通过R语言软件对原始数据进行处理,筛选差异表达基因(Log|FC|≥2,p0.05),利用STRING online-工具对差异表达基因进一步筛选(可信度≥0.4);Cytoscape软件构建蛋白质相互作用网络图(protein-protein interaction network,PPI),通过Cluster ONE对筛选后的差异表达基因进行聚类分析,筛选功能相关性密切的差异表达基因,并利用DAVID在线分析进一步进行GO分析和KEGG信号通路分析。m DCs相对于im DCs差异表达的基因共3 351个,上调基因1 801个,下调基因1 550个,其中C-C趋化因子受体活性、G-蛋白偶联的嘌呤核苷酸受体活性和异源三聚体G蛋白复合物等与细胞骨架调控密切相关。主要涉及的信号通路包括趋化因子/趋化因子受体信号通路、Rap1信号通路、Jak-STAT信号通路和PI3K-Akt信号通路等,其中的相关基因与细胞骨架的调控关系密切。这对进一步深入理解DCs的生物物理学特性、迁移能力和免疫学功能有一定的帮助。     

胰腺癌中CDK1的表达与预后的生物信息学分析

为探讨胰腺癌的发病机制并为胰腺癌的防治提供生物信息学依据,用GEO2R在线工具分析GSE16515中胰腺癌患者肿瘤组织和相应正常组织的差异表达基因(differentially expressed genes, DEGs),通过DAVID数据库对DEGs进行GO分析和KEGG通路富集分析,然后通过STRING数据库构建蛋白质相互作用(protein-protein interaction, PPI)网络,用Cytoscape软件进行关键基因(hub基因)筛选和功能模块分析,并在GEPIA数据库对hub基因进行验证,用CCLE数据库检测靶基因在胰腺癌组织及细胞系中的表达水平。分析结果显示胰腺癌中筛选出的376个DEGs主要涉及细胞周期、p53信号通路、蛋白质消化吸收、ECM-受体相互作用、PI3K-Akt信号通路、血小板激活信号通路。GEPIA数据库验证结果显示10个hub基因均在胰腺癌组织中高表达,其中8个hub基因与胰腺癌患者的不良预后有关。CCLE数据库检测结果显示周期蛋白依赖性激酶1 (cyclin-dependent kinase 1, CDK1)在胰腺癌组织和细胞中均有较高的表达水平。本研究结果表明CDK1可能与胰腺癌的发生发展最为相关,为进一步探究胰腺癌的发病机制提供了生物信息学依据。     

慢性鼻-鼻窦炎伴鼻息肉基因表达谱的生物信息学分析

本研究从美国国立生物信息中心(NCBI)的子数据库基因表达数据库(GEO)中选择基因表达谱GSE36830数据集,采用GEO2R筛选正常钩突和鼻息肉组织之间的差异表达基因(DEGs),对关键通路和差异表达基因进行数据库挖掘和分析,经筛选后的差异表达基因采用戴维在线工具对其进行基因本体富集分析(GO)、京都基因和基因组百科全书(KEGG)分析,然后将DEGs导入String数据库进行蛋白质互作网络分析,绘制差异表达基因互作网络图,将其数据导入Cytoscape软件中,筛选网络中心节点和关键基因,分析关键子网络。共筛选出699个DEGs,其中475个基因为上调表达基因,224个基因为下调表达基因。在GO分析中,针对生物过程,上调的DEGs包括:炎症反应、免疫反应、细胞趋化性、炎症反应的正向调节和细胞的粘附等;下调的DEGs主要参与:唾液分泌、生物矿物组织发展、细胞氨基酸生物合成过程、视网膜内稳态及离子跨膜转运等。在KEGG分析中,上调的DEGs主要在参与造血细胞系、细胞因子-细胞因子受体相互作用、破骨细胞分化、趋化因子信号通路、癌症中的转录失调、哮喘、金黄色葡萄球菌感染等信号通路中富集,而下调的DEGs在唾液腺分泌及胆汁分泌信号通路中富集。差异表达基因互作网络图筛选出前10个关键基因:ITGAM、IL10、CD86、TLR8、ITGAX、CCL2、CCR7、SRC、EGF及ITGB2。本研究得到了一组鼻息肉差异表达基因的生物信息学分析结果,但仍需进一步用基础试验来验证。本文分析的结论为慢性鼻-鼻窦炎、鼻息肉的研究提供了新的研究方向,也为鼻息肉发病机制研究的思路提供了一定的建设性作用。     

基于基因芯片技术的转Slnac2基因拟南芥差异表达基因分析

NAC转录因子,对植物的生长发育及抵御逆境胁迫起着重要的作用。实验室前期克隆了辽宁碱蓬(Suaeda liaotungensis)Slnac2基因,转Slnac2拟南芥(Arabidopsis thaliana)提高了抵御盐胁迫的能力。本研究利用基因芯片技术筛选转SlNAC2拟南芥差异表达基因,共筛选出差异表达2倍以上的基因1 258个。GO富集度分析结果显示,分子功能相关基因占34.06%,细胞组分相关基因占33.13%,生物学过程相关基因占32.81%。KEGG分析表明,差异表达基因涉及63个信号通路,主要有植物激素信号转导、核糖体、植物病原物相互作用和抗坏血酸代谢等。通过实时荧光定量PCR对部分差异基因进行验证,所得结果与基因芯片结果一致。本研究揭示,SlNAC2可调控多个下游基因的表达,提高植物在逆境胁迫下的生存能力。     

多杀性巴氏杆菌攻毒小鼠肺脏的转录组分析

多杀性巴氏杆菌(Pasteurella multocida)可感染多种动物宿主并引发败血症和呼吸道疾病,危害动物健康,给养殖业带来严重经济损失。本研究为探究多杀性巴氏杆菌感染期间对宿主肺脏基因表达的影响,选取羊源A型多杀性巴氏杆菌对小鼠(Mus musculus)进行攻毒实验,攻毒后72 h采集小鼠肺脏组织并提取总RNA。基于转录组测序,与对照组相比,攻毒组获得2 443个差异表达基因,其中有1 437个基因上调,1 006个基因下调。对2 443个差异表达基因进行GO富集分析,结果显示,差异基因显著富集的GO terms主要有B细胞稳态、T细胞迁移的正调控、整合素复合物和胶原蛋白结合;KEGG富集分析结果显示,差异基因显著富集于11条免疫炎症信号通路,包括细胞因子与细胞因子受体互作、趋化因子信号通路等;免疫炎症相关信号通路中,差异基因经蛋白质-蛋白质互作分析,结果显示,C3、C3ar1、C5ar1、Cxcl10、Cxcr2和Gng11基因处于网络中心,说明这些基因在多杀性巴氏杆菌感染过程中发挥了重要作用;从11条免疫炎症通路的基因中挑选10个进行实时荧光定量PCR验证,发现有8个基因...     

鼻咽癌组织中基因差异表达及蛋白互作网络分析

目的:利用生物信息学方法从鼻咽癌组织基因芯片中筛选差异表达基因,并分析它们之间的互相作用。方法:利用R语言程序包等相关软件分析鼻咽癌组织和正常组织中差异表达的基因,并对其进行在线GO分析、KEGG富集分析及蛋白互作分析。结果:从25例鼻咽癌组织样本和3例正常组织对照样本中共分析出1103个差异表达基因,包括600个上调基因和503个下调基因(P0.05)。GO分析结果显示,差异表达基因在生物过程中显著富集,包括细胞分裂、DNA复制、G1/S有丝分裂细胞周期的转变和有丝分裂核分裂;在分子功能中,包括与蛋白质结合、与ATP结合、蛋白同二聚化活性、与DNA复制起点结合以及与受损的DNA结合;在细胞组分中,包括细胞质、核质、胞外体和船体中部。KEGG通路分析显示,差异表达基因在细胞周期、DNA复制、癌症途径、p53信号通路、错配修复、小细胞肺癌、卵母细胞减数分裂、氨基糖和核苷酸糖代谢、基部切除修复和人T淋巴细胞病毒(HTLV-Ⅰ)感染等信号通路中显著富集。在蛋白互作分析网络中筛选出10个节点度最高核心基因CDC45、MCM10、MCM3、MCM5、MCM2、CDC7、CDT1、CDC6、SMC2和NCAPG。结论:通过鼻咽癌组织和正常组织的对比筛选,发现多个基因表达发生变化,为相关临床研究提供了重要的信息基础。     

黄芩-黄连药对治疗幽门螺杆菌相关性消化性溃疡的核心基因与关键microRNA筛选研究

本研究旨在应用网络药理学和生物信息学方法探讨黄芩-黄连药对(SRCR)治疗幽门螺杆菌(Hp)相关性消化性溃疡(PU)的分子网络调控机制,并筛选出其作用的核心基因和关键microRNA。首先通过中药系统药理学数据库和分析平台(TCMSP)检索SRCR的主要活性成分和可能的作用靶点;从GEO数据库获取GSE70394基因表达谱数据,采用R软件对Hp感染人胃黏膜上皮细胞(GES)涉及差异表达基因(DEG)的芯片数据进行均一化处理,并绘制DEG火山图和聚类热图;再通过STRING数据库构建SRCR成分靶点与Hp感染人GES细胞DEG的蛋白交互网络;运用DAVID数据库对上述特征性基因进行GO功能和KEGG通路富集分析;然后利用Cytohubba筛选出SRCR治疗Hp相关性PU涉及的核心基因;采用TargetScan和miRBase数据库对调控核心基因的关键microRNA进行预测。最终在TCMSP中共检索到SRCR的41个主要活性成分和216个可能的作用靶点。GEO芯片数据质量良好,以P0.05和︱logFC︱≥2为筛选条件,共获得128个Hp感染人GES细胞的DEG,其中上调基因77个,下调基因51个。取交集网络得到SRCR治疗HpPU的特征性基因靶点127个。对这些靶点进行GO分析发现,生物学过程涉及炎症反应、免疫反应等55个条目;细胞组分包含细胞外空间、外泌体、细胞核等11个条目;分子功能包含生长因子活性、表皮生长因子受体结合等15个条目。KEGG结果显示细胞因子受体相互作用通路、NOD样受体信号通路、NF-κB信号通路等33条信号通路与SRCR治疗Hp相关性PU密切相关。经Cytohubba筛选得到CXCL8、IL10、IL1B等十个SRCR调控的核心基因。进一步筛选得到miR-93-5p、miR-374b-5p、miR-3937等十个SRCR调控的最具功能性的microRNA。上述结果表明,SRCR治疗HpPU具有多成分、多靶点、多通路协同作用的特点,SRCR可通过调控编码RNA(基因)和非编码RNA(microRNA)发挥其抗HpPU的效应。     

烟草脉斑驳病毒侵染本氏烟的转录组研究

为了研究烟草脉斑驳病毒（TVMV）侵染本氏烟的抗病分子机制,本文建立了TVMV转录组数据库,挖掘抗病相关基因以及代谢和信号通路,采用Illumina Novaseq 6000高通量测序平台对侵染TVMV的本氏烟进行转录组测序,并进行生物信息学分析;利用R package:edge R等软件分析了有关抗病的差异表达基因,对差异表达基因进行基因本体论（GO）及京都基因和基因组百科全书（KEGG）基于途径的通路富集分析。基于TVMV侵染本氏烟转录组的测序结果中共获得4 593个差异表达基因,其中3 564个基因上调表达,1 029个基因下调表达。共筛选出10个抗病相关的差异表达显著基因,6个上调,4个下调。GO数据库中注释到生物学过程、细胞组分及分子功能等三大类共50个功能组。其中,在第三大类分子功能中,蛋白质结合功能以及ATP结合功能所涉及的差异表达基因最显著,分别为501个和453个差异基因。KEGG共富集20条通路,注释到核糖体途径的差异基因富集程度最高,基因最显著,差异表达基因数为307个。蛋白质结合、ATP结合、核糖体、真核生物核糖体的生物反应、DNA复制（DNA replicat...     

髓母细胞瘤基因表达谱芯片关键基因的生物信息学分析

筛选髓母细胞瘤(medulloblastoma, MB)发生发展的关键基因,可为MB分子机制的进一步研究提供生物信息学依据。本文通过下载GEO (Gene Expression Omnibus)数据库GSE50161原始数据,利用R语言对正常脑组织与髓母细胞瘤组织中差异表达的基因进行分析;通过生物信息学分析工具(DAVID、STRING和Cytoscape)对差异基因进行生物学功能和蛋白质相互作用(protein-protein interaction, PPI)分析,并通过PPI筛选互作调控的关键基因。结果显示,总共筛选出999个差异表达的基因,鉴定了CCNB1、AURKB、MAD2L1、CENPE、KIF2C、BUB1、BUB1B、NDC80、CENPF、CDC20十个关键基因。差异基因生物学功能主要富集于有丝分裂的核分裂、染色体分离、微管蛋白结合、RAGE受体结合等生物过程。KEGG信号通路分析结果显示差异基因主要富集于细胞周期、NF-κB、IL-17和T细胞受体等信号通路。10个关键基因的生物学功能和信号通路主要富集于细胞有丝分裂和细胞周期通路。因此,细胞周期通路对MB的增殖和分裂起着关键性的作用,相关的分子机制值得进一步深入研究。     

基于RNA-Seq技术分析江西豆豉对肠上皮细胞基因差异表达的影响

目的：基于RNA-Seq研究江西传统特色豆豉与肠上皮细胞的相互作用。方法：Caco-2细胞与豆豉水提样共孵育后,采用RNA-Seq筛选差异表达基因,进行KEGG、GO和PPI网络分析。结果：筛选差异表达基因17 892个,其中显著性240个,含上调表达125个、下调表达115个。KEGG功能注释215个信号通路,主要富集在新陈代谢、免疫系统、细胞调控、细胞粘附等类别。GO分析主要涉及免疫应答、细胞结构及功能调控、营养素代谢、神经传导及激素分泌调节等生物学过程。PPI网络分析在STRING数据库内有效注释的204个基因编码的蛋白互作,其中IL8、VCL、NFKBIA等节点度较高。结论：江西传统特色豆豉可引起肠上皮细胞多类基因差异表达,为阐明豆豉对肠道健康作用的分子机制提供科学依据。     

小菜蛾应对球孢白僵菌感染的免疫应答基因分析

通过室内生物测定,筛选对小菜蛾Plutella xylostella具有高毒力的球孢白僵菌Beauveria bassiana菌株,分析小菜蛾应对球孢白僵菌浸染的免疫应答基因。采用新一代Illumina高通量测序技术对接种/感染48h与未接种的小菜蛾样品分别进行转录组测序,筛选差异表达基因;结合生物信息学分析差异表达基因的功能、分类及涉及的信号通路等。结果表明接种小菜蛾48h后诱导2 434个基因发生差异表达,包括1 080个上调基因和1 354个下调基因。GO分类结果表明,这些差异表达基因（DEGs）被注释到45个GO term中,包括23个生物学过程,12个细胞组分和10个分子功能。KEGG分析表明,1 308个DEGs被富集到25个功能途径,这些DEGs包括497个上调基因和811个下调基因。分析发现,DEGs编码蛋白包含肽聚糖识别蛋白、丝氨酸蛋白酶55、丝氨酸蛋白酶抑制剂以及细胞色素P450 6K1类的免疫相关蛋白。另外,组蛋白、磷酸泛酰半胱氨酸脱羧酶、双链RNA结合蛋白、黑芥子酶等基因在球孢白僵菌侵染过程中也高表达,推测这些基因可能参与小菜蛾应对球孢白僵菌侵染的免疫反应。研究结果为进一步研究小菜蛾免疫相关基因的功能奠定了基础。