王圣桃 王新庄 高 铭 韩大勇

脑胶质瘤的发病率及病死率高，即使采用综合治疗方案（手术辅以放、化疗），5年生存率仍不足35%[1，2]。肿瘤基因的表达模式和特定分子机制是目前的研究热点[3，4]。研究表明，乳腺癌易感基因相互作用蛋白1（BRCA1 interacting protein C-terminal helicase 1，BRIP1）基因异常表达与卵巢癌、乳腺癌等恶性肿瘤发生、发展密切相关[5，6]。本文采用生信分析方法探讨胶质瘤BRIP1基因表达情况及其与病人预后的关系。

1 资料与方法

1.1 数据库检索 从GEO数据库GPL570平台选取三个BRIP1基因表达谱芯片（GSE4290、GSE50161和GSE74195），其中GSE4290包括23例正常样本和157例胶质瘤样本，GSE50161包括13例正常样本和34例胶质瘤样本，GSE74195包括5例正常样本和27例胶质瘤样本。

1.2 生存预后分析 使用TCGA下载有关胶质瘤BRIP1表达、生存率和临床特征的数据，包含326例低级别胶质瘤（low-grade glioma，LGG）样本和327例高级别胶质瘤（high-grade glioma，HGG）样本的mRNA-seq数据。以胶质瘤BRIP1表达中位值为界限，分成高表达和低表达组，使用R软件Kaplan-Meier法绘制生存曲线，评价BRIP1与胶质瘤病人生存时间的相关性。

1.3 胶质瘤病人预后影响因素分析 采用多因素Cox比例回归风险模型分析BRIP1表达水平与胶质瘤病人预后的关系。

1.4 BRIP1预测胶质瘤病人生存时间的价值 利用ROC曲线分析BRIP1基因表达水平预测胶质瘤病人生存时间的效能。使用logistic回归和Kaplan-Meier法绘制ROC曲线。

1.5 BRIP1基因相关信号通路富集分析 采用GSEA4.0.2 软件对BRIP1进行基因富集分析。以胶质瘤BRIP1表达中位值为基准进行分组，进行KEGG功能富集分析，排列数参数设置为1 000，基因组数据库设置为KEGG细胞信号通路，获得BRIP1的主要富集通路。

2 结果

2.1 胶质瘤BRIP1的表达情况 从GEO数据库GPL570平台选取GSE4290、GSE50161、GSE74195基因芯片进行分析，结果显示，相对于正常样本，胶质瘤BRIP1表达明显上调（P＜0.001；……