0 引言

多形性胶质母细胞瘤（glioblastoma multiforme, GBM）是成人最常见的恶性原发性脑肿瘤之一。目前主要治疗手段包括手术切除联合放疗，但疗效不佳，中位生存期仅12~15月^[1]。通过手术获得的肿瘤组织的组织学检查是目前确定性诊断的必要条件。尽管神经系统影像学检查可以提供GBM诊断价值，但对与GBM具有相似的影像学特征的其他脑病变诊断价值不大。因此，在组织学结果不确定或手术禁忌的情况下，筛选GBM分子标志物将为其诊断提供一定帮助。目前已经发现一些GBM中潜在的生物标志物，如PTEN^[2]、IDH1^[3]、TP53^[4]。

本研究利用生物信息学方法对基因芯片GSE7696进行分析来获得差异表达基因（differentially expressed genes, DEGs），并对差异基因进行聚类分析和功能富集分析，同时构建蛋白互作（protein-protein interaction, PPI）网络来筛选核心基因，最后通过TCGA肿瘤数据库GBM全基因组表达谱数据，对核心基因的表达水平进行验证，以提供GBM更多可能的潜在诊断和治疗的分子靶标。

1 材料与方法 1.1 芯片数据

从GEO（gene expression omnibus, http://www.ncbi.nlm.nih.gov/geo/）中下载编号为GSE7696^[5]的基因芯片，该芯片共有84例样本，包括4例正常脑组织样本（对照组）和80例GBM组织样本（实验组），比较实验组和对照组之间基因表达谱的差异情况。该芯片的平台信息：GPL570（HG-U133 Plus 2.0，Affymetrix Human Genome U133 Plus 2.0 Array）。芯片的探针注释信息来自Affymetrix公司，包含54 675条探针信息。芯片表达谱分析方法利用R软件（https://www.r-project.org/）。

1.2 样本的预处理及聚类分析

从GEO数据库下载原始CEL文件，利用R软件读取原始文件后，使用Affy包中的RMA算法标准化数据后得到基因的表达矩阵，计算样本之间的Pearson相关系数，对所有样本进行聚类分析，剔除离群样本。

1.3 差异表达基因分析

将预处理后得到的基因表达矩阵文件用R软件读入，利用R软件中limma包^[6]对80例GBM组织样本和4例正常脑组织样本进行差异表达分析，并且应用贝叶斯检验方法进行多重检验校正。差异基因筛选标准为P < 0.05，基因表达值倍数变化（fold change, FC）≥2或≤-2。

1.4 GO功能注释、KEGG信号富集以及通路网络构建

采用DAVID在线分析平台^[7]（https://david.ncifcrf.gov/）对差异基因在基因本体（gene ontology, GO）中注释其参与的生物学过程（biological process, BP），并进行京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）通路分析。采用Cytoscape软件对富集到的所有信号通路进行网络构建。

1.5 蛋白互作网络构建及核心蛋白筛选

利用String蛋白互作数据库^[8]（the Search Tool for the Retrieval of Interacting Genes, STRING, http://string-db.org/）分析GBM组织和正常脑组织差异基因之间的蛋白互作关系，构建蛋白质互作网络，以综合评分（combination score）大于0.4为阈值条件。将STRING中所得蛋白互作网络数据导入Cytoscape软件^[9]，运用其网络分析（network analyzer）插件计算节点的边（Degree, 即互作连线的数量）筛选网络中心节点(Hub Node)。中心节点对应的蛋白质为具有重要生理调节功能的核心蛋白质（基因）^[10]。

1.6 核心蛋白（基因）的验证

利用TCGA肿瘤数据库（The Cancer Genome Atlas Database, https://cancergenome.nih.gov/）获取169例GBM样本和5例正常脑组织样本全基因组表达谱数据，对PPI网络的核心基因进行进一步验证。

1.7 统计学方法

采用SPSS19.0统计软件。两组比较采用t检验，P < 0.05为差异有统计学意义。

2 结果 2.1 样本聚类分析

Pearson相关分析结果显示80例GBM组织样本和4例正常脑组织样本聚类良好，无离群样本，84例样本均可用于下一步分析，见图 1。

2.2 差异表达基因

以P < 0.05，基因表达值倍数变化≥2或≤-2为筛选条件。GBM组织和正常脑组织DEGs一共2 142个，其中上调基因968个，下调基因1 174个，见图 2。差异基因基本能够区分肿瘤组织样本和正常组织样本。

2.3 差异表达基因生物学功能注释

GO功能注释显示，GBM中差异表达基因富集到的生物学过程（biological process, BP）中富集度最高的20条，见图 3A。其中与GBM密切相关程度最高的3个生物学过程是突触传递（synaptic transmission）、有丝分裂细胞周期（mitotic cell cycle）、细胞分裂（cell division）。图 3A中每一条横柱代表一个生物学过程，横柱的长度代表富集到的差异基因数目，颜色越红表示P值越小，越具有统计学意义，下同。

2.4 差异表达基因KEGG信号通路及信号通路网络（pathway network）

利用DAVID在线富集工具将差异基因进行KEGG信号通路富集，图 3B显示富集度最高的20条通路。其中在GBM中明显富集且研究较多的3条信号通路有细胞周期（cell cycle）、MAPK信号通路（MAPK signaling pathway）以及PI3K-Akt信号通路（PI3K-Akt signaling pathway）。信号通路网络图（图 3C）显示MAPK信号通路、凋亡（apoptosis）、细胞周期以及P53信号通路（P53 signaling pathway）在通路网络中处于核心地位，说明他们在GBM的发生发展中起到极其重要的作用。图 3C中点的大小代表degree值，点越大代表通路越核心；颜色反映了通路中差异表达基因的情况，红色代表差异基因全部上调，蓝色代表下调，黄色代表既有上调又有下调。带箭头的实线表示两个信号通路之间的上下游关系，箭头起始端为上游的信号通路，箭头指向端为下游的信号通路。

2.5 通过Cytoscape软件构建的蛋白互作网络

网络分析显示，根据每个基因的节点数目排序，节点数目最多的基因，代表在整个PPI网络中所起的作用越大，即最相关的核心基因，有CDK1、KIF2C、BIRC5、CDC20、KIF11、CENPA、KIF20A、TOP2A、NDC80等9个核心基因，见图 4。这9个核心基因中全部在GBM中表达升高，差异有统计学意义，见表 1。

2.6 核心基因TCGA肿瘤数据库验证

从TCGA肿瘤数据库获取169例GBM样本和5例正常脑组织样本全基因组表达谱数据，比较这9个核心基因在GBM和正常组织中的转录水平，相对于正常组织，9个基因在GBM样本中表达明显升高，差异有统计学意义（P < 0.001），见图 5（略）。

3 讨论

目前基因芯片在疾病的诊断，特别是肿瘤和正常组织对比、肿瘤组织来源的鉴别等方面具有重要作用。它可用于肿瘤基因表达检测，寻找新基因，为寻找肿瘤的分子靶标提供了一个重要手段。

本研究从GEO芯片数据库获取了美国Affymatrix公司全基因组表达谱芯片，包括80例GBM组织和4例正常的脑组织样本，进行深入分析，以Fold Change绝对值≥2和差异显著性指标P < 0.05为条件，我们获得差异表达基因的数目为2 142个，其中上调968个，下调1 174个。这一结果提示在GBM发生发展的过程中，多种肿瘤相关的基因表达异常或者肿瘤抑制基因表达失活，说明GBM是由多个分子生物学调节异常共同导致。

为了探究这些调节异常基因在GBM中所起的生物学作用，我们将获得的2 142个差异基因进行GO功能注释，结果发现这些基因群主要参与细胞周期、细胞增殖、信号转导以及小分子物质代谢过程，由此说明这些生物学过程的调节异常是导致GBM发生发展的重要因素。同时，为了探究这些基因群所参与的信号通路过程，我们将差异基因进行了KEGG信号通路富集分析，同样发现细胞周期的异常调节在GBM中占有重要地位，其他明显富集到的信号通路有PI3K-AKT信号通路，肿瘤相关信号通路、钙离子信号通路以及MAPK信号通路。对富集到的信号通路进行网络分析后发现，MAPK信号通路处于通路网络的核心，由此说明GBM的发生与该信号通路密切相关。MAPK又称丝裂原活化蛋白激酶，是信号从细胞表面转导到细胞核内部的重要传递者，许多肿瘤的发生和发展都与MAPK信号通路的异常调节有关^[11]。以上结果对于研究GBM的生物学发生过程具有很好的提示作用，有待进一步验证。

在9个核心基因中，与细胞周期调节有关的有CDK1和CDC20。CDK1又称细胞周期蛋白依赖性激酶1，是一种高度保守的蛋白质，起到丝氨酸/苏氨酸激酶的作用，是细胞周期调控的关键参与者。CDK1的过度表达与多种肿瘤细胞的增殖密切相关^[12]。CDC20又称细胞分裂周期蛋白20，通过与其他几种蛋白质相互作用从而在细胞周期的多个点上起到调节作用^[13]。

与细胞分裂有关的基因有KIF2C、KIF11、CENPA、KIF20A、TOP2A、NDC80。KIF2C属于驱动蛋白样蛋白家族的成员，该家族的大多数蛋白质是微管依赖的分子马达，在细胞分裂期间运输细胞器中的细胞器并移动染色体。最新的研究表明，沉默KIF11可以导致染色体的不稳定从而促进肿瘤的发生^[14]。CENPA是修饰的核小体或核小体样结构的组成部分，其含有靶向着丝粒所需的组蛋白H3相关蛋白结构域，在细胞的有丝分裂中发挥重要作用，研究表明CENPA是乳腺癌潜在的预后生物标志物，其表达增加与乳腺癌的不良预后有关^[15]。KIF20A过表达与宫颈癌患者进展和不良预后有关^[16]，并且，研究发现KIF20A在神经胶质瘤组织中高表达并且提示预后不良^[17]。TOP2A基因编码DNA拓扑异构酶，是一种控制和改变转录过程中DNA拓扑状态的酶，在乳腺癌^[18]、神经母细胞瘤^[19]、肺癌^[20]等肿瘤中高表达，并且促进肿瘤细胞的增殖和迁移。NDC80在细胞有丝分裂的过程中主要参与染色体的正确分离的调控，有研究发现过表达NDC80后能够明显促进结肠癌细胞的增殖和迁移^[21]。

BIRC5又称survivin蛋白或者存活蛋白，凋亡抑制剂家族的成员。survivin蛋白抑制半胱天冬酶活性的作用，从而负性调控导致凋亡。抑制survivin蛋白导致细胞凋亡增加和肿瘤生长减少^[22]。数据表明survivin蛋白可能作为癌症治疗的新靶标。最近研究发现survivin蛋白的核表达是GBM患者预后不良的一个因素，survivin蛋白亚细胞定位可以帮助预测用标准方案治疗的GBM患者的总生存期^[23]。

综上，应用基因表达谱芯片可以筛选GBM相关基因及信号通路，是否是GBM特异的敏感基因需要更多病例验证以及细胞实验的进一步研究，希望从中能筛选出有助于GBM早期诊断和治疗的分子标志物或基因治疗靶点。