基于共表达网络分析鉴定与膀胱癌相关的枢纽基因和通路

Identification of hub genes and pathways associated with bladder cancer based on co-expression network analysis.

作者信息Dong-Qing Zhang, Chang-Kuo Zhou, Shou-Zhen Chen, Yue Yang, Ben-Kang Shi
PMID28693282
发布时间2017-07
DOI10.3892/ol.2017.6267

实验完整度

中

包含生物信息学分析(DEG筛选、共表达网络、通路富集)和临床样本验证(RT-qPCR和Western blot),但缺乏功能实验和动物模型验证。

主要模型

膀胱癌患者肿瘤组织样本 正常膀胱组织样本

重点核对

BC患者肿瘤组织和正常样本各10例,来源于齐鲁医院活检 筛选DEGs的cut-off标准:|logFC|≥2且P<0.05 通路富集分析使用KEGG数据库,显著性阈值为EASE score < 0.01 RT-qPCR每个样本重复3次,以β-actin作为内参 Western blot以GAPDH作为内参,使用10 μg蛋白进行SDS-PAGE

摘要

本研究的目的是利用中心性分析和通路富集分析鉴定与膀胱癌(BC)相关的枢纽基因和信号通路。从ArrayExpress数据库中筛选出正常受试者和BC患者之间的差异表达基因(DEGs)。利用差异共表达基因和链接构建BC的共表达网络,并通过BC中共表达网络的度中心性分析研究枢纽基因。基于DEGs,通过京都基因与基因组百科全书数据库分析研究富集的信号通路。使用逆转录定量聚合酶链反应(RT-qPCR)和western blotting验证BC组织中枢纽基因的表达。共筛选出329个DEGs,包括147个上调和182个下调基因。在BC和正常对照之间构建的共表达网络包含182个节点和434条边,每个基因对中的两个基因均为差异共表达基因。共表达网络的中心性分析表明,度值最高的前5个枢纽基因包括半乳糖凝集素4(LGALS4)、蛋白酪氨酸磷酸酶受体N2型(PTPRN2)、跨膜蛋白酶丝氨酸11E(TMPRSS11E)、三结构域蛋白31(TRIM31)和钾电压门控通道亚家族D成员3(KCND3)。通路分析显示,329个DEGs在5个条目中显著富集(细胞周期、DNA复制、卵母细胞减数分裂、p53信号通路和过氧化物酶体增殖物激活受体信号通路)。根据RT-qPCR和western blot分析,5个枢纽基因中有4个显著表达,包括LGALS4、PTPRN2、TMPRSS11E、TRIM31;然而,KCND3未显著表达。在本研究中,成功鉴定了5个枢纽基因(LGALS4、PTPRN2、TMPRSS11E、TRIM31和KCND3)以及5条可能作为膀胱癌早期诊断和治疗潜在生物标志物的生物学通路。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
本研究旨在通过生物信息学分析识别与膀胱癌相关的枢纽基因和信号通路,并通过实验验证其表达,以探索膀胱癌的分子机制和潜在生物标志物。
核心机制
本研究发现膀胱癌中差异表达基因显著富集于细胞周期、DNA复制、卵母细胞减数分裂、p53信号通路和PPAR信号通路,提示这些通路可能参与膀胱癌的发生发展。
主要证据
通过对两个基因表达数据集(共146个样本)的差异表达分析,共筛选出329个DEGs,构建了包含182个节点和434条边的共表达网络,并通过度中心性分析确定了前5个枢纽基因;RT-qPCR和Western blot验证显示LGALS4、TMPRSS11E和TRIM31的表达趋势与生信分析一致。
研究意义
这些基因和通路可能作为膀胱癌早期诊断和治疗的潜在生物标志物。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

数据收集与预处理

获取膀胱癌和正常对照的基因表达数据并进行标准化,为后续分析提供可靠的数据基础。

从ArrayExpress数据库获取E-MTAB-1940和E-GEOD-3167两个基因表达数据集,共包含146个样本(128例膀胱癌,18例对照)。数据经过背景校正、归一化、探针映射和合并等预处理步骤。

2

差异表达基因筛选

筛选膀胱癌和正常对照之间显著差异表达的基因。

使用LIMMA包对预处理后的数据进行差异表达分析,通过t检验和F检验,并采用经验贝叶斯和错误发现率校正P值,最终根据|logFC|≥2且P<0.05的标准筛选出329个DEGs(147个上调,182个下调)。

3

共表达网络构建与枢纽基因识别

构建膀胱癌的差异共表达网络,并通过中心性分析识别关键枢纽基因。

使用DCGL 2.0包中的DCp和DCe方法识别差异共表达基因和链接,构建包含182个节点和434条边的共表达网络。通过度中心性分析,定义度值≥9的基因作为枢纽基因,共获得18个,其中前5个为LGALS4、PTPRN2、TMPRSS11E、TRIM31和KCND3。

4

通路富集分析

识别差异表达基因显著富集的信号通路。

基于KEGG数据库,使用DAVID工具对329个DEGs进行通路富集分析,以EASE score < 0.01为阈值,发现5条显著富集的通路:细胞周期、DNA复制、卵母细胞减数分裂、p53信号通路和PPAR信号通路。

5

枢纽基因表达验证

在膀胱癌患者组织中验证5个枢纽基因的mRNA和蛋白表达水平。

收集10例膀胱癌患者的肿瘤组织和距肿瘤2 cm的正常组织,使用RT-qPCR和Western blotting检测5个枢纽基因的表达。结果发现LGALS4、TMPRSS11E和TRIM31的表达趋势与生信分析一致,PTPRN2趋势相反,KCND3无显著差异。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
BCA 蛋白定量试剂盒CW Biotech--
TRIzol 试剂Beyotime Institute of Biotechnology--
AMV 逆转录酶Invitrogen--
Taq DNA 聚合酶Invitrogen--
硝酸纤维素膜----
Amersham ECL 蛋白质印迹检测试剂盒GE Healthcare Life Sciences--
兔抗人 LGALS4 抗体Abgent Biotech (SuZhou) Co., Ltd.AP12391a
兔抗人 PTPRN2 抗体BIOSSbs-19591R
兔抗人 TMPRSS11E 抗体Abgent Biotech (SuZhou) Co., Ltd.AP16520C
兔抗人 TRIM31 抗体BIOSSbs-6220R
兔抗人 KCND3 抗体BIOSSbs-20219R
兔抗人 GAPDH 抗体Abgent Biotech (SuZhou) Co., Ltd.AP50811
辣根过氧化物酶标记的羊抗兔免疫球蛋白G二抗Beijing Zhongshan Jinqiao Biotechnology Co., Ltd.SC-2048

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
数据收集与预处理
数据集E-MTAB-1940和E-GEOD-3167的来源和样本量(共146个样本:128例BC,18例对照),预处理方法(RMA背景校正、分位数归一化、MAS 5.0算法等),数据集合并方法(Batch Mean-centering和inSilicoMerging包)
阅读提示:Materials and methods: Data collection and preprocessing
差异表达分析
DEG筛选标准(|logFC|≥2,P<0.05),使用的软件包(LIMMA)和统计方法(t检验、F检验、经验贝叶斯、FDR校正)
阅读提示:Materials and methods: Detection of DEGs
共表达网络分析
DCGL 2.0包的具体参数(DCp和DCe方法,置换检验次数N),枢纽基因定义(度值≥9)
阅读提示:Materials and methods: Construction of differential co-expression networks by DCGL 和 Centrality analysis
通路富集分析
使用的数据库(KEGG),工具(DAVID),显著性阈值(EASE score < 0.01)
阅读提示:Materials and methods: Pathway enrichment analysis
枢纽基因表达验证
患者样本来源(齐鲁医院活检,10例BC患者,正常组织取自距肿瘤2cm处),RT-qPCR条件(表II中的退火温度和循环数),Western blot条件(抗体稀释比、蛋白上样量10 μg等)
阅读提示:Materials and methods: Patients 和 Validation of hub genes by RT-qPCR and western blotting;Table I和Table II;Figure 3和Figure 4图注