应用生物信息学筛选结直肠癌Hub基因及验证

系统医学 2024 年 2 月第 9 卷第 3
期
应用生物信息学筛选结直肠癌Hub 基因及验证
陈树华，温日葵，祝惠钦，谢荣章
云浮市人民医院检验科，广东云浮 527300
[摘要] 目的 通过运用生物信息学分析，筛选出与结直肠癌相关的差异表达基因（Differentially Expressed genes, DEGs ），并验证其生物学功能。

方法 云浮市人民医院检验科从基因表达综合数据库（Gene Expression Omnibus, GEO ）中下载结直肠癌芯片数据GSE 21815、GSE 31905、GSE 35279资料，应用GEO2R 语言进行处理得出结直肠癌和正常结直肠组织之间的差异表达基因，并通过生物信息学工具DAVID 、STRING 、Cytoscape 构
建差异表达基因的蛋白互作网络，筛选Hub 基因，应用基因本体论（Gene Ontology, GO ）、基因百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG ）分析筛选出的Hub 基因的生物功能，并利用MiRDB 工具找出可能调控Hub 基因的miRNA ，并于2017—2022年8月收集30例结直肠癌组织和30例正常结直肠组织样本，通过实时荧光定量PCR （Quantitative Real-time PCR, qPCR ）验证。

结果 经生物学信息分析和蛋白质相互作用网络图分析催产素受体基因、基质金属蛋白酶11基因、间质上皮转化因子基因、基质金属蛋白酶7基因、激肽释放酶8基因、激肽释放酶10基因为结直肠癌组织发生发展的关键Hub 基因。

结直肠癌组织中基质金属蛋白酶11基因（4.38±1.58）、间质上皮转化因子基因（2.69±0.29）、基质金属蛋白酶7基因（0.88±0.14）、激肽释放酶8基因（11.09±3.90）、激肽释放酶10基因mRNA （7.88±2.20）的表达，显著高于正常结直肠组织组织，差异有统计学意义（t =9.605、25.339、26.376、9.541、3.726，P 均<0.001）。

结论 结直肠癌组织中基质金属蛋白酶11基因、间质上皮转化因子基因、基质金属蛋白酶7基因、激肽释放酶8基因、激肽释放酶10基因异常表达，可能参与结直肠癌发生过程，有望为后续结直肠癌基础研究及临床诊疗提供依据。

[关键词] 结直肠癌；生物信息学；Hub 基因
[中图分类号] R735 [文献标识码] A [文章编号] 2096-1782（2024）02(a)-0042-05
Screening and Validation of Hub Gene in Colorectal Cancer by Bioinfor⁃matics
CHEN Shuhua, WEN Rikui, ZHU Huiqin, XIE Rongzhang
Departments of Laboratory, Yunfu People′s Hospital, Yunfu, Guangdong Province, 527300 China
[Abstract] Objective Through bioinformatics analysis, the differentially expressed genes (DEGs) associated with colorectal cancer were screened out and their biological functions were verified. Methods The clinical laboratory of Yunfu People′s Hospital downloaded colorectal cancer chip data GSE 21815, GSE 31905 and GSE 35279 from Gene Expression Omnibus (GEO) database. The differentially expressed genes between colorectal cancer and normal colorectal tissue were obtained by GEO2R language. The protein interaction network of differentially expressed genes was constructed by bioinformatics tools DAVID, STRING and Cytoscape, and Hub genes were screened. The biologi⁃cal functions of Hub Genes screened were analyzed using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Mirnas that may regulate hub genes were identified using the MiRDB tool, and Quantitative Real-time PCR (qPCR) was used to verify 30 colorectal cancer tissue samples and 30 normal colorectal tissue samples col⁃
lected from August 2017 to 2022. Results Oxytocin receptor gene, matrix metalloproteinase-11 gene, mesenchymal epithelial transformation factor gene, matrix metalloproteinase-7 gene, kallikrein-8 gene and kallikrein-10 gene were
analyzed as the key Hub genes in the development of colorectal cancer tissue by biological information analysis and DOI ：10.19368/ki.2096-1782.2024.03.042
[基金项目] 2021年云浮市医药卫生科技计划项目（20210310）
[作者简介] 陈树华（1983-），男，本科，副主任技师，研究方向为生化、核酸检测。

系统医学 2024 年 2 月第 9 卷第 3
期
protein interaction network diagram. The mRNA expression of matrix metalloproteinase 11 gene (4.38±1.58), stromal epithelial transformation factor gene (2.69±0.29), matrix metalloproteinase 7 gene (0.88±0.14), kallikrein 8 gene (11.09±3.90), and kallikrein 10 gene (7.88±2.20) in colorectal cancer tissue were significantly higher than those in normal colorectal tissue, and the differences were statistically significant (t=9.605, 25.339, 26.376, 9.541, 3.726, all P<0.001). Conclusion The abnormal expression of matrix metalloproteinase 11 gene, stromal epithelial transformation factor gene, matrix metalloproteinase 7 gene, kallikrein 8 gene, and kallikrein 10 gene in colorectal cancer tissue may be involved in the development process of colorectal cancer, and is expected to provide a basis for subsequent basic research and clinical diagnosis and treatment of colorectal cancer.
[Key words] Colorectal cancer; Bioinformatics; Hub gene
结直肠癌为临床中常见的消化道恶性肿瘤之一，其病死率在癌症中排第3位，严重威胁人类健康。

因此，早期识别结直肠癌诊断和预后相关生物标志物至关重要[1]。

生物信息学为生物学、信息学相结合的内容，在阐述疾病分子机制方面起着积极的作用。

基因芯片主要是通过微阵列技术将高密度DNA片段阵列依附着玻璃、尼龙等材料上，筛选出有价值的基因进一步研究分析，目前常用于收集疾病表达谱数据[2]。

Hub基因是经生物信息学筛选出的核心基因，也是当前临床疾病治疗的潜在靶点[3]。

云浮市人民医院检验科运用生物信息学方法，筛选出结直肠癌相关的差异表达基因予以生物通路富集分析、制作蛋白-蛋白互作网络，筛选出Hub基因，并于2017—2022年8月收集30例结直肠癌组织和30例正常结直肠组织样本进一步验证Hub基因表达，为今后临床诊治结直肠癌、评估预后转归提供新型的辅助工具。

现报道如下。

1 对象与方法
1.1研究对象及芯片数据资料
收集2017—2022年8月本院就诊的30例结直肠癌患者组织标本作为异常组，另收集30例非癌结肠组织作为正常组。

从美国国立生物技术信息中心的基因表达综合数据库（Gene Expression Omni⁃bus, GEO）中下载结直肠癌芯片数据GSE 21815、GSE 31905、GSE 35279资料，其中GSE 21815包括
12例结直肠癌组织和9例正常结直肠对照组织；GSE 31905包括55例结直肠癌组织和7例正常结直肠组织；GSE 35279包括74例癌组织和5例正常组织。

1.2方法
1.2.1数据预处理运用tidyverse软件对获取的芯片数据资料予以数据处理，对各样本的基因文件进行命名，并把基因的表达量置于同一表达矩阵内。

在metadata元数据库中把样本注释内容提取出来，使其与表达矩阵相符。

最后通过R语言的biomaRt 包对标准基因符号进行标注，进而取得标准化的基因表达矩阵。

1.2.2筛选差异基因应用edgeR包、limma包将数据预处理获得的基因表达矩阵进行差异表达基因分析，再应用GEO中的分析工具analyze with geo2r分为结直肠癌组织组和正常组织组，以∣logFC∣>2，P< 0.01作为有效基因的纳入标准[4]，筛选出GSE 21815、31905、35279数据集中的差异基因，并绘制韦恩图，找出GSE 21815、31905、35279数据集中的共有基因。

1.2.3提取lncRNA在Gencode数据库中收集和整理经过筛选的差异基因lncRNA相关信息，通过R 软件lncRNA内ensembl ID及其表达量等信息。

1.2.4生物功能富集分析运用生物信息学工具DAVID，将与正常结直肠组织有差异的结直肠癌相关差异表达基因（Differentially Expressed Genes, DEGs）筛选出来，通过基因本体（Gene Ontology, GO）开展功能富集分析（包括细胞组分、分子功能、生物过程3个部分），并进行京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）通路富集分析。

1.2.5蛋白互作网络（Protein Protein Interaction, PPI）分析在STRING数据库上参与富集分析的DEGs基因进行PPI分析，应用Cytoscape构建PPI网络。

具体操作为：借助网络分析算出网络中心性参数—度中心性，用MCODE插件提取关键子网，获得相对集中的表达子集，识别出结直肠癌的潜在调控基因。

最后，用ClueGO、CluePedia把KEGG通路可视化，绘制生物通路之间的基因互作网络，提取与通路之间相互作用程度排名前10的差异基因作为Hub基因。

系统医学 2024 年 2 月第 9 卷第 3
期
1.2.6 结直肠癌组织样本Hub 基因表达验证 应用
实时荧光定量PCR （Quantitative Real-Time PCR, qPCR ）测量癌组织、正常结直肠组织中hub 基因表达量没具体方法。

①提取总RNA ：取0.3~0.5 g 组织样本，将其研磨成粉末状，置于1.5 mL 离心管内，根据RNA 提取试剂盒说明书流程提取总RNA 。

②逆转录反应：根据逆转录试剂盒的操作流程进行逆转录反应，取200 ng RNA 量，滴加1 μl 特异性逆转录引物+12 μL 无菌水，孵育5 min 后，依次滴加Ribol⁃lock RNase 抑制剂、RevertAid M-MuLV 反转录酶。

③RT-PCR 反应检测：以SYBR green 法，设置反应体系：95℃ 10 min ，95℃ 15 s ，60℃ 1 min ，40个循环，以U6作为内参基因，以2-△△Ct 作为基因表达量，以正常结直肠组织为对照，得出结直肠癌组织中hub 表达量。

△△Ct=（Ct 目的基因-Ct U6）结直肠癌组织-
（Ct 目的基因-Ct U6）正常结直肠组织。

1.3 统计方法
应用SPSS 26.0统计学软件进行数据统计分析，计量资料（mRNA 表达量）符合正态分布，以（x ˉ±s ）描
述，行t 检验，P <O.05为差异有统计学意义。

2 结果
2.1 差异基因筛选结果
GSE 21815（937）、31905（1470）、35279（1458）
数据集共有141例结直肠癌组织和21例正常结
直肠组织，3个数据集分别检索到937个、1 470个、1 458个DEGs 基因，在软件上绘制韦恩图后获得共有上调基因105个、共有下调基因140个。

2.2 核心基因验证
经PPI 网络图分析可知，245个共有差异基因导入到STRING 数据库生成PPI 网络图，并对PPI 网络图进行可视化分析得出与PPI 网络相互作用程度最高的10个Hub 基因，分别为催产素受体基因（Oxyto⁃cin Receptor, OXTR ）、基质金属蛋白酶11基因（Ma⁃trix Metalloproteinase-11, MMP11）、酪氨酸蛋白激酶受体B2基因（Erythropoietin-Producing Hepatoma-
B2, EPHB2）、间质上皮转化因子基因（Mesenchymal-Epithelial Transition Factor, MET ）、转化生长因子β1
基因（transforming Growth Factor-β1, TGF-β1）、抑制斯钙素-2基因（Stanniocalcin-2, STC2）、基质金属蛋白酶7基因（Matrix Metalloproteinase-7, MMP7）、激肽
释放酶8基因（Kallikrein-8, KLK8）、激肽释放酶10基因（Kallikrein-10, KLK10）、角蛋白23基因（Kera⁃tin-23, KRT23）。

结直肠癌组织的MMP11、c-MET 、MMP7、KLK8、KLK10 mRNA 表达量明细高于正常结直肠组织，差异有统计学意义（P 均<0.05）。

而结直肠癌组织与正常结直肠组织的OXTR 、EPHB2、TGF-β1、STC2、KRT23 mRNA 表达比较，差异无统计学意义（P 均>0.05），见表1。

3 讨论
目前已有大量基础研究、临床研究揭示了结直肠癌发生发展的病因病机[5-6]，但其患病率、病死率近年仍居高不下，有学者认为其原因主要是大多数研究是针对单个遗传学事件进行队列研究，造成结果偏倚[7]。

本研究选择GSE 21815、GSE 31905、GSE 352793个数据库运用生物信息学深入研究，筛选获得245个结直肠癌相关差异表达基因，而且对上调、下调245个共有差异基因进行GO 生物功能分析结果显示，①上调DEGs 基因的细胞组分（Cellu⁃
lar Component, CC ）涉及细胞外基质、细胞外区域、细胞质膜等；下调基因细胞组分涉及细胞外区域、细胞外谜题、细胞间膜蛋白等。

②上调基因的分子功能主要是序列特异性DNA 结合蛋白活性、RNA 聚合酶Ⅱ、生长因子激活等；下调基因分子功能主要为激素激活体系、结合锌离子的蛋白质。

③上调基因的生物过程（Biological Process, BP ）包括细胞增殖、细胞凋亡、药物反应、蛋白质水解等；下调基因生物过程包括蛋白质水解、小分子物质跨膜转运等。

经KEGG 富集分析可知，上调DEGs 基因主要
表1　结直肠癌组织与正常结直肠组织10个Hub 基因mRNA 表达量比较(x
ˉ±s )组别
结直肠癌组织（n =30）正常结直肠组织（n =30）t 值P 值
OXTR
0.87±0.090.79±0.121.6610.102
MMP114.38±1.581.57±0.279.605
<0.001EPHB22.99±0.512.85±0.381.2130.230
MET 2.69±0.291.01±0.2125.339
<0.001TGF-β11.40±0.311.24±0.361.8560.068
STC2
1.43±0.231.37±0.211.0240.310
MMP7
0.88±0.140.17±0.0426.376
<0.001KLK811.09±3.904.00±1.189.541
<0.001KLK10
7.88±2.205.96±1.763.726
<0.001KRT23
0.96±0.150.98±0.150.5450.588
注：OXTR ：催产素受体基因，MMP11：基质金属蛋白酶11基因，EPHB2：酪氨酸蛋白激酶受体B2基因，MET ：间质上皮转化因子基因，TGF-β1：转化生长因子β1基因，STC2；抑制斯钙素-2基因，MMP7：基质金属蛋白酶7基因，KLK8：激肽释放酶8基因，KLK10:激肽释放酶10基因，KRT23：角蛋白23基因。

富集分布转化生长因子-β（Transforming Growth Factor-β, TGF-β)信号通路、Wnt信号通路；下调DEGs基因富集于PPAR信号通路、氮代谢信号通路
等。

由此可见，上调基因主要分布在细胞核、细胞膜、细胞外，通过调节DNA复制、参与细胞周期等过程以参与结直肠癌细胞的细胞迁移、增殖和凋亡。

而下调基因主要分布在细胞外区域，参与机体代谢，进而对癌症发生发展起作用。

本研究制作差异基因蛋白互作PPI网络获得10个Hub基因，并经临床组织样本验证发现，结直肠癌组织中基质金属蛋白酶11基因（4.38±1.58）、间质上皮转化因子基因（2.69±0.29）、基质金属蛋白酶7基因（0.878±0.143）、激肽释放酶8基因（11.09±3.90）、激肽释放酶10基因mRNA（7.88±2.20）mRNA表达显著高于正常结直肠组织组织（P均<0.05）。

基质金属蛋白酶（MMPs）对细胞外基质、基底膜蛋白酶有溶解作用，在正常情况下该物质在体内的表达水平较低，但若发生病理改变，特别是出现癌细胞增殖分化、迁移等，其表达水平会显著升高[8]。

MMP11就是MMPs 中的一员，可通过重构细胞外基质诱导肿瘤进展，并可通过抑制细胞凋亡使肿瘤细胞存活[9]。

MMP7是MMPs中分子量最少的分泌蛋白，其能够溶解细胞外基质，并可水解机体蛋白多糖、胶原底物，促进肿瘤细胞增殖分化[10]。

目前已有研究报道，MMP11（4.25±1.19）、MMP7（0.886±0.125）在多种癌组织中呈高表达[11]。

李军彦等[12]的研究显示，MMP7表达能够提示结直肠癌病情发生及其进展。

温凌等[13]研究也发现，检测MMP-11能够辅助诊断结直肠癌发生发展情况。

C-Met为原癌基因的一种，能够使活化信号传导通路发生磷酸化，促使细胞运动，进而引起上皮细胞分散、内皮细胞迁移，肿瘤细胞发生浸润迁移。

另有研究报道，C-Met能够诱导肿瘤血管新生[14]。

曹明等[15]研究证实，结直肠癌患者癌组织C-
Met表达水平与癌细胞侵袭转移呈正相关，对于判断预后有一定指导意义。

已有研究证实，KLKs与恶性肿瘤发生发展密切相关，其作用机理是通过编码hK 蛋白溶解肿瘤细胞外基质蛋白，促使癌细胞迁移和肿瘤血管新生[16]。

相关文献报道，KLK8能够促进结直肠癌细胞增殖、转移，其过表达能够抑制肿瘤细胞凋亡，敲低KLK8后会促使癌细胞凋亡[17]。

另有研究发现，KLK10在结直肠癌组织中异常表达，与肿瘤病理分期、肝转移发生呈正相关[18]。

由此可见，MMP11、c-MET、MMP7、KLK8、KLK10高表达，能够促进结直肠癌细胞增殖分化、浸润迁移、
综上所述，本文利用生物学信息方法进行筛选和验证了一组与结直肠癌发生发展密切相关的基因，有望成为作为早期预测结直肠癌的标志物，为后续结直肠癌基础研究及临床诊疗提供依据。

[参考文献]
[1] 王晓娜,刘忠于,龚慧,等.结直肠癌相关生物标志物
研究进展[J].结直肠肛门外科,2022,28(1):34-38. [2] 中国抗癌协会肿瘤靶向治疗专业委员会.结直肠癌
分子检测高通量测序中国专家共识[J].临床肿瘤学杂志,2021,26(3):253-264.
[3] 范馨月,邵春莹,崔雷.基于PPI网络预测和验证结直
肠癌相关Hub基因[J].中华医学图书情报杂志,2019, 28(6):17-25.
[4] 殷晓晓,薛小峰,赵华.基于生信分析筛选胆囊癌关键
基因及相关通路的研究[J].中国肿瘤外科杂志,2021, 13(1):12-19.
[5] 陈素华,马天江,张智慧,等.长链非编码RNA LINC01419调控结直肠癌细胞增殖、凋亡、迁移侵袭的机制[J].安徽医药,2022,26(3):572-577,封3. [6] 赵莉,韩涛涛,李可敏,等.结肠癌与直肠癌癌变相关
基因及蛋白表达的差异分析[J].北京医学,2023,45(4): 284-288.
[7] 龚慧,黄燹,陈晖娟,等.结直肠癌表观遗传标志物研
究进展[J].结直肠肛门外科,2022,28(1):14-18. [8] 易默,吕飒美,都庆国,等.基质金属蛋白酶11、巨噬
细胞抑制细胞因子-1联合检测对结直肠癌诊断及预后评估的价值[J].癌症进展,2022,20(5):505-507,512. [9] 史春桃,程杰,刘慧,等.SPOCK1、MMP11与MMP9在
结直肠癌中的表达及临床价值[J].重庆医学,2023,52 (1):102-107.
[10] 陈伟锋.免疫组织化学法检测MMP-7表达判断直
肠癌安全手术切缘的参考价值[J].黑龙江医药,2021,
34(3):708-709.
[11] 骆苏彦,李婷婷.血清基质金属蛋白酶-11、血管生
（下转第53页）
[3] 党红民,崔丹丹.复方聚乙二醇电解质散剂用于结肠
镜检查前肠道准备的效果[J].中国肛肠病杂志,2020, 40(9):54-55.
[4] 赵明星,朱玉侠,孔祥云,等.小剂量复方聚乙二醇电
解质散剂联合二甲硅油散在便秘患者肠道准备中的应用[J].医学综述,2019,25(21):4364-4368. [5] 傅才招,林美英,林香花,等.复方聚乙二醇电解质散
剂联合二甲硅油散的肠道准备效果[J].中国老年保健医学,2019,17(5):67-68,71.
[6] 李华萍,纪桂岚.复方聚乙二醇电解质散(Ⅱ)在妇科HIFU治疗前肠道准备中的应用[J].护理实践与研究, 2018,15(22):75-77.
[7] 张淑芬,薛挥,李建宏,等.聚乙二醇电解质散剂联合
二甲硅油在结肠镜检查肠道准备中的应用效果观察[J].陕西医学杂志,2020,49(2):227-230. [8] 唐学萍.聚乙二醇电解质散剂不同口服方法在结直
肠癌术前肠道准备中的效果[J].中国乡村医药,2017, 24(21):12-13.
[9] 郭松,李军文,谢林娟,等.聚乙二醇电解质散剂联合
乳果糖对肠道清洁有效性和安全性的Meta分析[J].
包头医学院学报,2022,38(7):92-96.
[10] 周娟,薛晓玲,陈亚莉,等.妇科术前聚乙二醇电解
质散剂不同使用方法进行肠道准备的临床效果比较[J].临床医学研究与实践,2017,2(25):83-84. [11] 林五连,张晓兰,何小建,等.聚乙二醇电解质散剂联
合低剂量硫酸镁与西甲硅油在结肠镜检查前肠道准备中的应用[J].现代肿瘤医学,2017,25(4):587-590.
[12] 陆娟,曹勤,王晓春,等.乳果糖口服液联合聚乙二醇
电解质散剂在便秘患者行结肠镜检查前肠道准备
中的应用研究[J].现代生物医学进展,2015,15(20):
3920-3922.
[13] 孙付坤,章敏.不同时间服用复方聚乙二醇电解质散
在接台妇科手术术前肠道准备中的效果比较[J].护
士进修杂志,2019,34(23):2196-2198. [14] 邹伟聪.针灸联合加味小承气汤灌肠治疗术后早期
炎性肠梗阻86例分析[J].中医临床研究,2016,8(12):
111-112.
[15] 张艳芹.肠梗阻导管联合小承气汤治疗妇科术后早
期炎性肠梗阻临床观察[J].光明中医,2019,34(13):
2050-2052.
[16] 劳玲娟.复方聚乙二醇电解质散加小承气汤用于结
肠镜检前肠道准备的效果观察[J].中国肛肠病杂志,
2012,32(11):54-56.
[17] 王慈英,周江妍.聚乙二醇电解质散剂用于妇科腹腔
镜术前肠道准备的效果评价[J].中国妇幼保健,2012,
27(21):3352-3354.
[18] 李福民,蒋自乐,丁杰,等.小承气汤在胸腰椎骨折术
后便秘治疗中的疗效分析[J].中医临床研究,2022,14
(24):97-99.
[19] 李记彬,胡秀芬,刘放,等.大肠癌根治术腹腔镜手术
与开腹手术长期疗效和安全性的比较[J].河北医学,
2015,21(2):210-212.
[20] 郭桂英.小承气汤结合护理干预对妇科腹部术后肠
道功能恢复的影响[J].光明中医,2019,34(9):1441-
1443.
（收稿日期：2023-11-06）
成素-2在结直肠癌诊断中的应用价值[J].肿瘤基础
与临床,2022,35(3):198-200.
[12] 李军彦,陈文超,李中辉,等.结直肠癌患者血清
MMP7、CD44v7、Survivin、PPARγ表达及其临床
意义[J].实用癌症杂志,2022,37(7):1068-1071. [13] 温凌,王巧瑜,麦绮彬,等.5种基因mRNA在结直肠癌
患者粪便中的表达及临床意义[J].临床消化病杂志,
2023,35(4):301-306.
[14] 樊桂玲,梁晓霞,高宝红,等.结直肠癌c-MET、
CXCR4蛋白和微血管密度与肝转移的关系[J].肿瘤
研究与临床,2021,33(11):830-833.
[15] 曹明,怀建国,杨小苗.MACC1、c-met及ERCC1的表
达在结直肠癌病理评估中的应用价值研究[J].当代
医药论丛,2023,21(9):5-9.
[16] 李晓松,郑军,姚汝铖,等.人组织激肽释放酶对结直
肠癌患者的预后判断价值研究进展[J].山东医药,
2016,56(34):99-102.
[17] 花晴,申雪芳,许平波.KLK8通过调节EGF对结直肠
癌细胞凋亡的影响及其机制[J].中国癌症杂志,2021,
31(5):383-389.
[18] 郭佩,何宁宁,夏玉军.基于数据挖掘分析人组织激
肽释放酶家族在结直肠癌中的表达和预后研究[J].
癌症进展,2021,19(2):149-153,177.
（收稿日期：2023-11-05）
（上接第45页）。