论文部分内容阅读
目的:胶质母细胞瘤(GBM)是一种高度恶化的弥漫性神经胶质瘤,具有极强的浸润性。GBM的主要治疗方法为手术切除联合放化疗,但术后易复发,且预后不良,其中位生存期仅有14个月左右。因此,本文基于生物信息学分析,寻找与GBM不良预后相关的分子标志物,将有助于筛选可靠生物分子标志物,并对其进行风险评估,并具有关键意义。方法:从癌症基因组图谱(TCGA)中下载GBM相关数据集,利用R语言进行差异基因分析,FDR<0.05,fold change>1,并绘制相关火山图和热图。以P<0.05作为显著性检验判别标准,对差异基因进行生存分析。利用DAVID8.4对影响总生存率的差异基因进行基因本体论(GO)和京都基因组百科全书(KEGG)富集分析,探究相关靶基因参与的生物过程、细胞组成、分子功能和信号通路;利用String数据库获取对总生存有影响的差异基因的蛋白-蛋白相互作用网络,联合Cytoscape软件中MCODE插件对前述基因进行模块分析。然后,利用Cytoscape中ClueGO插件将与GBM相关的基因簇进行GO和KEGG基因富集分析,观察通路之间的关系和基因在通路中的富集情况。同时,对该基因簇及其临床资料进行Cox单因素回归分析,将可能影响患者生存期的变量纳入到Cox多变量回归分析中。根据风险评分中位数将患者分为高风险组和低风险组进行预后生存分析,并检验模型的灵敏度和特异性。结果:1、本研究从TCGA中下载的数据包括5例正常组织和156例GBM组织样本。通过差异分析得到6690个基因,其中上调基因为3460个和下调基因为3230个。2、对差异表达的6690个基因进行预后分析,得到对预后有影响的323个上调基因和177个下调基因。GO注释显示靶基因参与的生物过程主要集中在细胞外基质组织、整合素介导的信号通路等;靶基因参与的细胞组成主要有胶原三聚体、粘着斑等,靶基因参与的分子功能有胶原结合、整联蛋白结合等。KEGG富集分析显示靶基因参与ECM-受体相互作用、粘着斑等信号通路。3、利用String数据库和Cytoscape软件对差异表达的500个基因进行蛋白-蛋白相互作用网络可视化分析后,并通过特定条件筛选,得出与GBM的发生、发展生物学进程与信号通路密切相关的涵盖25个基因的基因簇。4、对基因簇所包含的基因进行Cox单因素回归分析得出23个具有显著预后的基因。对上述基因进行Cox多因素回归分析,得出涵盖8个基因的线性预测模型。风险评分=0.373×ExpRCN1+0.349×ExpFN1+0.523×ExpLAMB1-0.727×ExpLAMC1+0.747×ExpITGB8-0.448×ExpTNC+0.227×ExpHSPG2+0.325×ExpITGB5。5、根据中位风险评分,将156位GBM患者分为高风险组(n=78)和低风险组(n=78)进行预后分析,结果显示高风险组预后不良。ROC曲线面积AUC=0.822,表明8基因模型对GBM患者的预后具有良好的灵敏度和特异性。结论:1、CKAP4、COL1A1、COL5A1、COL6A1、COL6A2、COL6A3、DMP1、FAM20C、FN1、HSPG2、IGFBP1、ITGA1、ITGA3、ITGA5、ITGB5、ITGB8、LAMA4、LAMB1、LAMC1、MXRA8、PTK2、PXN、RCN1、SCG3、TNC基因簇参与胶原分解代谢、整联素所介导的信号通路、细胞外基质结构和组织、基底层、血小板衍生生长因子结合、粘着斑等GBM发生、发展生物学相关进程与信号通路。2、筛选得到RCN1、FN1、LAMB1、LAMC1、ITGB8、TNC、HSPG2、ITGB5可以作为GBM患者预后的联合标志物,相对于单个独立预后因子具有良好的灵敏度和特异性。