基于生物信息学分析的非小细胞肺癌诊断预后相关基因的筛选

为了从基因层面探讨非小细胞肺癌(non-small cell lung cancer, NSCLC)发生发展的内在机制,筛选与NSCLC诊断、预后相关的基因,为NSCLC分子机制的进一步研究提供生物信息学依据,利用生物信息学方法对GEO数据库和TCGA数据库的数据集进行合并分析,筛选NSCLC组织与正常肺组织之间的差异表达基因(differentially expressed genes, DEGs),并对所取交集的DEGs进行基因集富集分析(gene set enrichment analysis, GSEA)、基因本体论(gene ontology, GO)分析、KEGG (kyoto encyclopedia of genes and genomes)通路富集分析、蛋白质相互作用(protein-protein interaction, PPI)分析、ROC曲线诊断效能分析及LASSO生存分析。文中共筛选出240个DEGs,主要涉及核分裂、染色体分离等生物学过程。GSEA分析结果显示,富集的通路主要涉及DNA修复和细胞周期。从PPI网络中筛选出20个hub基因, ROC结果显示, UBE2C (AUC=0.939)、TOP2A(AUC=0.927)、RRM2 (AUC=0.927)、CCNB1 (AUC=0.928)、MKI67 (AUC=0.930)、AURKA (AUC=0.931)、MELK(AUC=0.950)相对具有较高的诊断价值, LASSO COX回归结果则显示IL6、KIAA0101、MKI67、TPX2、AURKA、CDKN3及CDCA5与NSCLC患者的预后强相关。本研究结果表明, ZWINT、KIF2C、MELK、CDCA5可能在NSCLC中发挥着重要的作用,为阐明NSCLC的分子机制提供了新思路。     

胰腺癌中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可能与胰腺癌的发生发展最为相关,为进一步探究胰腺癌的发病机制提供了生物信息学依据。     

肝细胞癌相关差异基因的生物信息学及预后分析

肝细胞癌是全球癌症相关死亡的主要原因,目前对肝细胞癌的发病机制研究尚不完善,探索肝细胞癌发生、发展相关的分子标志物及其预后具有重要意义。从GEO数据库获得肝细胞癌组织和非癌组织的基因表达阵列数据GSE84402,利用GEO2R筛选差异表达基因;采用DAVID数据库对差异基因进行GO富集分析和KEGG通路分析;通过STRING数据库和Cytoscape软件构建差异表达基因对应的蛋白质相互作用网络,并从网络中筛选出核心基因(hub genes);结合KM plotter数据库的临床信息对hub genes进行预后分析。结果显示:共得到1 307个差异表达基因,其中上调基因741个,下调基因566个,这些差异表达基因主要涉及细胞分裂、细胞周期、DNA复制及物质代谢等生物学过程及生物通路。通过GO、KEGG及蛋白质相互作用网络筛选出BUB1、BUB1B、CCNA2、CCNB1、CCNB2、CDC20、CDK1、MAD2L1、PLK1等9个hub genes,进一步分析发现hub genes均与细胞周期的调控相关,表明细胞周期的调控失常在肝细胞癌的发生、发展过程中具有重要作用。生存分析显示9个hub genes在肝细胞癌患者中均为表达上调的基因,且与患者预后不良相关,这为寻找肝细胞癌患者预后相关生物标志物的研究提供了线索。     

基于拓扑分析方法鉴定与胃癌相关的生物标志物

胃癌(gastric cancer, GC)是最常见的恶性肿瘤之一。由于GC发病隐匿的特性,其早期检测困难。因此,研究与GC早期诊断和预后相关的生物标志物至关重要。从GEO数据库下载了3组基因表达数据集GSE79973、 GSE19826和GSE13911,通过Limma包筛选差异表达基因(differentially expressed genes, DEGs),并使用DEGs、 STRING V11数据库和Cytoscape构建了DEGs的蛋白质-蛋白质相互作用(protein-protein interaction, PPI)网络,通过4种拓扑分析方法取交集筛选hub基因,并通过单变量Cox分析、多变量Cox分析、 Lasso回归分析、生存分析、通路分析以及文献法验证hub基因。从3个数据集中分别筛选了1 599个、 333个和662个DEGs。通过拓扑分析方法筛选了4个hub基因,即CDK1、AURKA、PTTG1和UBE2C。GO和KEGG富集分析结果表明4个hub基因参与了细胞外基质-受体相互作用、糖尿病并发症中的AGE-RAGE信号通路、小细胞肺癌和蛋白质消化吸收等通路。...     

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

筛选髓母细胞瘤(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的增殖和分裂起着关键性的作用,相关的分子机制值得进一步深入研究。     

人肝细胞癌预后不良相关基因的生物信息学分析及其临床意义

目的:筛选肝细胞癌(HCC)预后不良相关基因,并探讨其临床意义。方法:在基因表达综合数据库(GEO)中获取符合分析条件的肝细胞癌全基因组表达谱数据并分析得到差异表达基因(DEGs),再运用生物学信息注释及可视化数据库 (DAVID) 和蛋白相互作用数据库 (String) 分别进行功能富集分析和蛋白质互作用网络的构建。利用癌症基因组图谱数据库(TCGA)和Cox比例风险回归模型对相关差异基因进行预后分析。结果:找到一个符合条件的人类HCC数据库 (GSE84402),共筛选出1141个差异表达基因(DEGs),其中上调基因720个,下调基因421个。基因功能富集分析和蛋白质互作用分析结果显示CDK1、CDC6、CCNA2、CHEK1、CENPE 、PIK3R1、RACGAP1、BIRC5、KIF11和CYP2B6为HCC预后的关键基因。TCGA数据库和Cox回归模型分析显示CDC6、PIK3R1、RACGAP1和KIF11的表达升高,CENPE的表达降低与HCC预后不良密切相关。结论:CDC6、CENPE、PIK3R1、RACGAP1和KIF11可能和HCC的预后不良相关,可作为未来HCC预后研究的参考标志物。     

肺腺癌诊断标志物筛选及免疫细胞浸润分析

用生物信息学方法筛选肺腺癌(Lung adenocarcinoma,LUAD)的诊断生物标志物,并分析肺腺癌中免疫细胞浸润情况。从GEO和TCGA数据库下载肺腺癌的表达数据集,利用R软件筛选肺腺癌与正常肺组织间的差异表达基因(DEGs),使用DAVID网站对DEGs进行GO及KEGG富集分析,使用STRING及Cytoscape等工具对DEGs构建蛋白相互作用网络并筛选hub基因;利用Kaplan-Meier法对DEGs进行生存分析,并对hub基因进行ROC分析筛选诊断生物标志物,利用GSEA预测有预后价值的基因参与的信号通路;并用Cibersort软件反卷积算法分析肺腺癌中免疫细胞浸润情况。共得到肺腺癌的234个DEGs,这些基因主要参与信号转导、物质代谢、免疫反应等相关信号通路;构建PPI网络筛选出的20个hub基因中8个存在预后价值(CCNA2、DLGAP5、HMMR、MMP1、MMP9、MMP13、SPP1、TOP2A),ROC分析中DLGAP5、SPP1值分别是0.703、0.706;DLGAP5、SPP1基因表达水平与肺腺癌组织浆细胞、未活化的CD4⁺记忆细胞、调节T细胞、巨噬细胞M0、M1、M2及中性粒细胞浸润密切相关(P＜0.05)。肺腺癌中DLGAP5、SPP1具有较高诊断价值且参与肺腺癌组织免疫细胞浸润;DLGAP5、SPP1基因可作为肺腺癌诊断的生物标志物,可为肺腺癌的靶向治疗研究提供新思路。     

小儿脓毒症休克相关基因的筛选及网络构建

为探究脓毒症休克与SIRS的差异表达基因及网络的构建,筛选潜在的核心基因,从GEO数据库下载相关基因表达谱GSE26378,数据分为脓毒症休克与SIRS各29个样本,通过在线软件GCBI对其进行标准化及差异基因筛选;对差异基因进行GO分析;基于KEGG进行功能通路分析以及基因信号网络分析;差异基因共表达网络分析。结果表明:两组中总共有1 456个基因被识别为差异基因(P0.05),与SIRS组相比,脓毒症休克组中有条859条下调基因,597条上调基因。GO功能富集分析显示差异基因主要参与了细胞周期、细胞免疫、细胞代谢。KEGG功能通路分析显示差异基因主要参与了MAPK信号通路、P53信号通路、wnt信号通路、细胞凋亡信号通路,细胞周期受体信号通路等。共表达分析发现基因CCNB1、NUSAP1、OIP5、SHCBP1、ZWINT、TOP2A、DLGAP5等位于网络中央部位,而基因信号网络分析发现基因PLCB1、PIK3CA、STAT3、CAMK2D、PRKCB、CREB1位于网络核心。基因芯片分析有助于发现脓毒症休克与SIRS患儿外周血单核细胞在转录组学上的改变,而生物信息学网络分析有助于发现潜在的靶点。     

生长激素对感染J亚群禽白血病病毒的鸡胚成纤维细胞系的影响

生长激素(growth hormone, GH)不仅能调节动物的代谢、生长、发育和生殖,也能调节动物的免疫系统。本研究旨在通过转录组测序及生物信息学技术分析GH影响J亚群禽白血病病毒(subgroup J avian leukosis viruses, ALV-J)感染的鸡胚成纤维细胞系(chicken fibroblast cells, DF-1)中的关键基因和信号通路。将空载和GH的过表达质粒分别转染至DF-1,接着均接种ALV-J毒株SCAU-HN06,采用转录组测序技术检测相关的基因表达谱变化。结果显示,与对照组相比,实验组中有上调差异表达基因(differentially expressed genes, DEGs)122个(49%),下调DEGs 127个(51%)。通过GO功能富集分析发现,DEGs主要富集于新陈代谢过程、生物调节、细胞过程、细胞部分、催化活性和结合等过程。经KEGG富集分析发现DEGs主要在免疫系统、信号转导、信号分子和相互作用、细胞生长和死亡等通路上富集。进一步构建蛋白质-蛋白质相互作用(protein-protein interaction, PPI...     

IgA 肾病患者与膜性肾病患者外周血单核细胞mRNA分析

为了寻找诊断、鉴别IgA肾病(IgAN)和膜性肾病(MN)的血液特异性标记物,利用公共数据库中的IgAN和MN患者的外周血单核细胞(PBMCs)的转录组表达谱数据集识别特异性生物标记物,为诊断和鉴别提供简便、可靠的依据补充。从公共基因表达数据库(GEO)下载IgAN患者组(n=15)和MN患者组(n=8)芯片数据集,筛选前250个差异表达基因(DEGs)。通过分析筛选关键基因和途径,进行基因本体(GO)富集分析、京都基因与基因组百科全书(KEGG)通路分析和蛋白质与蛋白质相互作用关系(PPI)分析等进一步了解DEGs。通过分析共发现75个显著DEGs,其中73个上调基因,2个下调基因。GO富集分析的生物学过程(BP)主要包括蛋白质转运、内溶酶体到溶酶体转运、趋化因子介导的信号通路作用等。显著富集差异表达基因KEGG通路分析包括Endocytosis和Hepatitis B的相关信号通路。PPI筛选出EPS15、STAT4、CCL2、SUN2、SEC24C、SEC31A、GOLGB1、F2R,RAB12和PTK2B等关键基因。成功筛选出核心差异表达基因,为IgAN和MN的诊断和鉴别提供简便、可靠的依据补充,甚至提供治疗的新靶点。     

FGF13敲低的肺癌A549细胞株构建及转录组学分析

本研究的目的是探究成纤维细胞生长因子13(fibroblast growth factor13, FGF13)影响非小细胞肺癌发生发展的分子机制。通过构建干扰FGF13的真核表达载体并进行慢病毒包装,获得FGF13敲低的肺癌A549细胞株。利用转录组测序(RNA Sequencing, RNA-seq)技术检测基因表达谱并筛选差异表达基因(differentially expressed genes, DEGs),并对差异基因进行GO和KEGG富集分析,对部分DEGs的表达水平进行验证。采用STRING在线数据库筛选DEGs编码的蛋白相互作用网络(protein-protein interaction, PPI)中的关键基因,并通过UALCAN数据库对关键基因进行预后分析。与对照组细胞相比,FGF13敲低的肺癌A549细胞株共筛选出1 114个DEGs(P≤0.05,差异倍数≥2),其中表达上调和下调的DEGs分别有672和442个。实时荧光定量PCR(real-time quantitative PCR, RT-qPCR)验证CXCL10和SELE等6个基因与测序结果一致。GO富集分...     

肾纤维化关键基因的生信筛选及分析

[目的]利用生物信息学技术筛选与肾纤维化相关基因，并预测与基因相关通路，疾病、细胞类型和有关药物。[方法]在公共基因芯片数据库(GEO)中获取与肾纤维化相关三个基因数据集，利用Venn图鉴定出共表达的差异基因，Metascape工具对差异基因进行功能(GO)和通路富集(KEGG)分析，还利用STRING构建蛋白相互作用网络(PPI)网络以及可视化工具Cytoscape筛选关键基因，最后用及Enrichr和Comparative Toxicogenomics Database(CTD)数据分析有关疾病、细胞类型和药物。[结果]与肾纤维化相关的数据集GSE148420、GSE38117、GSE54441经Venn图共得到83个DEGs。经Cytoscape对STRING工具得到的PPI网络可视化后，筛选出10个与肾纤维化相关的关键基因，分别是F13B、ALDH8A1、A1CF、PAH、KMO、ALDH6A1、SPP2、ACAT1、ABAT、CAT。GO和KEGG富集显示这些基因与氧化还原酶活性，血小板致密颗粒，色氨酸、结氨酸、亮氨酸和异亮氨酸等氨基酸代谢通路有关。利用Enrichr和CTD...     

家兔前体脂肪细胞分化不同时期基因表达谱分析

脂肪的过度积累严重危害人类健康。前体脂肪细胞分化是脂肪发育的关键过程,研究前体脂肪细胞分化相关基因的表达有助于认识脂肪沉积的机理。尽管家兔是一种理想的研究脂肪发育的动物模型,但是针对其前体脂肪细胞分化不同时期基因表达谱的研究鲜见报道。本研究通过诱导家兔前体脂肪细胞分化,在分化第0 d、3 d和9 d收集脂肪细胞,利用转录组测序(RNA-seq),在分化第3 d样本与第0 d样本的比较中筛选出1352个差异表达基因(differentially expressed genes, DEGs),在分化第9 d样本与第3 d样本的比较中筛选出888个DEGs。GO (gene ontology)功能富集和KEGG (kyoto encyclopedia of genes and genomes)通路分析发现,0~3 d分化期上调的DEGs显著富集在PPAR信号通路和PI3K-Akt信号通路上,3~9d分化期上调的DEGs显著富集到与细胞周期调控有关的GO条目和KEGG信号通路,0~3d和3~9d阶段特异上调的DEGs可能分别作用于细胞质和细胞核。通过DEGs的蛋白-蛋白互作(protein-protein interaction, PPI)网络分析发现,筛选出的核心节点(hub node)基因可能通过调控细胞周期而影响家兔前体脂肪细胞分化。     

脊髓损伤机制研究的生物信息学分析

用生物信息学方法筛选参与脊髓损伤(spinal cord injury, SCI)发展过程的关键分子和通路,可为脊髓损伤发展机制的研究提供指导。从GEO数据库下载基因芯片数据,并将数据集中的样本分为脊髓损伤组(SCI组)和正常组(normal组)。应用R语言处理来自不同数据集样本间的批次效应,同时对基因芯片的表达数据进行标准化处理,并通过PCA分析监测标准化处理后数据的质量。应用R语言中的limma包分析标准化后的基因表达矩阵,以得到差异基因。将差异基因导入DAVID数据库进行GO (gene ontology)分析,并通过KEGG数据库进行通路分析。然后应用STRING数据库构建PPI网络,并通过Cytoscape中的cytoHubba插件分析得到10个hub基因。最后应用箱式图监测hub基因在不同样本中的表达,并用GeneCards数据库查询hub基因的功能。此外,为了补充差异基因筛选的不足,通过R语言对基因表达矩阵进行了GSEA (gene set enrichment analysis)分析。结果显示:TYROBP、ITGB2、PTPRC和FCER1G等基因在脊髓损伤发展过程中发挥重要的作用;细胞外基质的炎症反应、葡糖醛酸基转移酶活性的变化和星形胶质细胞的迁移等与脊髓损伤的发展机制关系密切; TNF信号通路、NF-κB信号通路和p53信号通路在脊髓损伤的发展机制中发挥重要的作用。这些关键的分子和通路在脊髓损伤中的作用值得我们进行更深入的探讨。     

肾透明细胞癌Caki-1细胞系差异表达基因的生物信息学分析

目的比较肾透明细胞癌Caki-1细胞系与正常肾上皮细胞系ASE-5063中的差异表达基因（DEGs）,寻找潜在的肾透明细胞癌特异性分子标志物。 方法利用GEO数据库自带的GEO2R在线分析工具分析基因芯片GSE78179,将筛选出的DEGs分别导入Metascape、STRING以及Cytoscape进行综合分析并筛选出核心基因。最后使用FunRich等软件对筛选出的核心基因进行GO和KEGG富集分析。 结果共筛选出562个DEGs,其中上调基因345个,下调基因217个。进一步使用MCODE筛选出36个关键基因,GO功能分析发现这些基因与细胞粘附分子活性、趋化因子活性、细胞通讯和信号转导等密切相关;KEGG通路富集结果则表明差异基因主要集中在趋化因子信号通路、TNF信号通路以及NF-κB信号通路等多种与肿瘤相关的通路上。 结论运用生物信息学方法筛选出肾透明细胞癌Caki-1细胞系中DEGs,其中数个核心基因广泛参与多种肿瘤的病理进程,但尚未在肾透明细胞癌有相关研究报道,提示其可能是治疗肾透明细胞癌的潜在靶点。     

基于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个关键基因对于胶质母细胞瘤的发生和发展可能起到了关键作用,可以作为潜在的调控位点和筛选的标志物。     

GWGS芯片整合分析识别外周血细胞中与骨质疏松症相关的枢纽基因

为了对骨质疏松症基因芯片数据集进行整合分析并识别出外周血细胞中与骨质疏松症相关的枢纽基因,通过检索GEO和ArrayExpress数据库获得骨质疏松症相关的表达谱芯片数据集;运用GWGS (genome-wide global significance)方法对纳入的数据集进行整合分析,筛选出差异表达基因(differentially expressed genes, DEGs);然后,运用GO (gene ontology)富集分析和KEGG (kyoto encyclopedia of genes and genomes)通路富集分析对差异表达基因进行功能注释,并建立蛋白质相互作用(protein-protein interaction, PPI)网络,筛选出骨质疏松症相关的枢纽基因。公共数据库检索得到3个符合纳入排除标准的研究集, GWGS整合分析筛选出排序前200的DEGs,这些基因主要富集的GO条目为脂多糖的细胞反应、凋亡过程和炎症反应,与骨质疏松症相关的KEGG富集通路为破骨细胞分化等。PPI分析进一步检测到与骨质疏松症相关的10个枢纽基因,其中9个基因已有研究报道和骨质疏松症的发生发展相关,而ELANE基因还未有研究报道与骨质疏松症有关。ELANE基因同时在人的骨髓组织、小鼠骨髓和骨组织中高表达,这个基因很可能与骨质疏松症有潜在的联系。本研究的结果将有助于进一步理解骨质疏松症的分子致病机理。     

强直性脊柱炎差异表达基因的生物信息学分析

目的探讨强直性脊柱炎(AS)患者差异表达基因,并基于差异基因探讨强直性脊柱炎发病相关的可能生物学过程和信号通路。方法检索基因表达谱数据库(GEO)并筛选AS相关基因表达谱数据集。应用GEO在线分析功能GEO2R分析AS组和正常对照组的差异表达基因,用Cytoscape软件clueGO插件进行基因本体论和京都基因与基因组百科全书分析,采用String蛋白-蛋白相互作用(PPI)数据库分析差异表达基因编码蛋白间的相互作用;应用Cytoscape绘制蛋白相互作用网络图,并软件筛选信号通路关键基因分析。结果选取AS患者全血表达数据集GSE25101为研究对象,分析获得差异表达基因72个。72个差异表达基因分子功能主要为参与高迁移率族盒染色体蛋白1(HMGB1)转导机制;生物学过程主要富集于巨噬细胞迁移、骨髓细胞凋亡过程、线粒体呼吸链复合体装配、ATP合成偶联电子传输、线粒体ATP合成耦合电子输运等;细胞成分主要富集于呼吸链复合体、线粒体呼吸体等。信号通路富集于氧化磷酸化信号通路和帕金森综合征相关信号通路。PPI网络经过cytohubba插件筛选,ATP5J、NDUFS4、UQCRB、UQCRH、NDUFB3、COX7B、LSM3、ATP5EP2、ENY2、PSMA4被筛选为网络中的核心基因。结论通过生物信息学方法进行预测了AS的潜在机制,并筛选出10个潜在的与AS相关的重要分子,其中氧化磷酸化可能在AS发病机制中发挥了重要的作用。     

关于食管癌亚型食管基底细胞样鳞癌的基因芯片数据分析

分析食管基底细胞样鳞癌(basaloid squamous cell carcinoma of the esophagus,BSCCE)的组织样本基因表达情况,探讨其病理学过程和肿瘤免疫逃逸机制,同时对BSCCE独特的肿瘤微环境及潜在的临床指标在对其远期预后的病理基础进行讨论,以期为临床工作中更好的治疗BSCCE提供一定的理论基础。利用临床术中获取的组织样本及其基因芯片数据,进一步获取BSCCE组织样本与同体正常组织样本的基因差异表达情况(log|FC|≥2,p0.05),其中BSCCE组织相对于正常组织的上调表达差异表达基因共489个,下调差异表达基因922个,利用DAVID在线分析系统对差异表达基因进行京都基因与基因组百科全书(kyoto encyclopedia of genes and genomes,KEGG)信号通路分析获取对BSCCE组织病理情况相关的分析结果,进一步利用STRING-OL工具对差异表达基因的可信度(可信度≥0.4)进行筛选,获取差异表达基因的蛋白质相互作用分析(protein-protein interaction,PPI)结果,进一步利用Cytoscape软件对筛选后的PPI结果制作相应的基因相互作用网状图,并进一步筛选核心表达的差异基因。在上调差异基因的KEGG分析中获得了88个结果,如p53signaling pathway、Cytokine-cytokine receptor interaction pathway、Toll-like receptor signaling pathway等,在下调的差异基因KEGG分析中获得了125个结果,如Drug metabolism-cytochrome P450 pathway、MAPK signaling pathway、Calcium signaling pathway、Cytokine-cytokine receptor interaction pathway等,这些相关的信号通路对进一步理解BSCCE的肿瘤微环境变化及肿瘤病理相关机制有重要意义,通过PPI分析,我们筛选出了一些与肿瘤远期生存率(overall survival,OS)相关的基因:CDC6、CDC20、TOP2A、NEK2、CDC25A等,这为我们进一步研究BSCCE及其肿瘤病理标志物和远期生存的关系提供了理论基础。