基于癌症基因图谱挖掘前列腺癌不同Gleason分级癌组织相关基因分析

doi:10.16150/j.1671-2870.2018.06.013

摘要/Abstract

摘要： 目的：运用生物信息学探究不同Gleason分级（Gleason Pattern,GP）前列腺癌组织间的分子表达差异及其基因间的相互影响。方法：下载癌症基因图谱（the cancer genome atlas,TCGA）前列腺癌患者的转录组数据及GP资料,运用生物信息学分析GP 4级与GP 3级前列腺癌组织间的基因表达差异,进而通过基因本体（gene ontology, GO）、京都基因和基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）信号通路富集及蛋白相互作用网络（protein protein interaction network, PPI）,探究不同GP的前列腺癌组织间差异表达基因的相互作用及影响的生物学过程。结果：生信分析共鉴定出312个差异表达基因,且相较于GP 3级癌组织,GP 4级前列腺癌组织中存在157个基因表达上调,155个基因表达下调。生信分析显示,在基因相互作用网络中富集到的22个显著相关基因主要参与的生物学过程为细胞周期及有丝分裂。结论：GP 4级前列腺癌组织内的肿瘤细胞有丝分裂更为活跃,临床可通过检测细胞周期相关蛋白,协助对前列腺癌患者作出更合理的疾病风险及预后评估;同时,针对Gleason积分（Gleason score,GS）≥7分的患者,细胞周期相关蛋白有望成为其有效的治疗靶点。

关键词: 前列腺癌, Gleason评分, 癌症基因图谱, 生物信息学, 细胞周期

Abstract: Objective: To investigate the differentially expressed genes involved in distinguished Gleason patterns of prostate cancer (PC) and to reveal the potential molecular mechanisms by applying integrated bioinformatics. Methods: The expression profiles and information for Gleason score were downloaded from the cancer genome atlas (TCGA) database. The differentially expressed genes (DEGs) were obtained and were further analyzed for their interaction and biological processes involved by bioinformatics methods. The gene ontology(GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichments of DEGs were used for processing by DAVID online analyses. The protein-protein interaction (PPI) networks of the DEGs were constructed from the STRING database. Results: A total of 312 DEGs were identified in the TCGA datasets, of which 157 genes were upregulated and 155 genes were downregulated, in the Gleason pattern 4 PC. The 22 most closely related genes among DEGs were identified from the PPI network and were focused primarily on cell cycle and mitotic division. Conclusions: This study indicated that cancer cells in Gleason pattern 4 PC are more dynamic in the process of mitotic cell cycle. In clinical practice, detection of protein associated with cell cycle may allow the more rational prediction for risk and prognosis of PC. Meanwhile, for patients with Gleason score ≥7 PC, genes associated with cell cycle may provide the effective targets for treatment of prostate cancer.

Key words: Prostate cancer, Gleason score, The cancer genome atlas, Bioinformatics, Cell cycle

中图分类号:

R737.25

王涛, 邓玉, 赵萍, 于宝华, 王翔, 王朝夫. 基于癌症基因图谱挖掘前列腺癌不同Gleason分级癌组织相关基因分析[J]. 诊断学理论与实践, 2018, 17(06): 694-700.

WANG Tao, DENG Yu, ZHAO Ping, YÜ Baohua, WANG Xiang, WANG Chaofu. Identification of genes associated with distinguished Gleason patterns of prostate cancer by analyzing TCGA database[J]. Journal of Diagnostics Concepts & Practice, 2018, 17(06): 694-700.

参考文献

[1] 韩苏军, 张思维, 陈万青, 等. 中国前列腺癌发病现状和流行趋势分析[J]. 临床肿瘤学杂志,2013,18(4):330-334.
[2] Nwosu V, Carpten J, Trent J M, et al.Heterogeneity of genetic alterations in prostate cancer: evidence of the complex nature of the disease[J]. Hum Mol Genet,2001, 10(20):2313-2318.
[3] Allsbrook WC Jr, Mangold KA, Johnson MH, et al.Inte-robserver reproducibility of Gleason grading of prostatic carcinoma: urologic pathologists[J]. Hum Pathol,2001,32(1):74-80.
[4] Shen M M, Abate-Shen C.Molecular genetics of prostate cancer: new prospects for old challenges[J]. Genes Dev,2010,24(18):1967-2000.
[5] Tomczak K, Czerwińska P, Wiznerowicz M.The Cancer Genome Atlas(TCGA): an immeasurable source of know-ledge[J]. Contemp Oncol(Pozn),2015,19(1A):A68-A77.
[6] Cancer Genome Atlas Network. Comprehensive molecular portraits of human breast tumours[J]. Nature,2012,490(7418):61-70.
[7] Cancer Genome Atlas Research Network. Comprehensive genomic characterization of squamous cell lung cancers[J]. Nature,2012,489(7417):519-525.
[8] von Mering C, Huynen M, Jaeggi D, et al. STRING: a database of predicted functional associations between proteins[J]. Nucleic Acids Res,2003,31(1):258-261.
[9] Eggener SE, Scardino PT, Walsh PC, et al.Predicting 15-year prostate cancer specific mortality after radical prostatectomy[J]. J Urol,2011,185(3):869-875.
[10] Wang H, Zhang C, Rorick A, et al.CCI-779 inhibits cell-cycle G2-M progression and invasion of castration-resistant prostate cancer via attenuation of UBE2C transcription and mRNA stability[J]. Cancer Res,2011,71(14):4866-4876.
[11] Horning A M, Wang Y, Lin C K, et al.Single-Cell RNA-seq Reveals a Subpopulation of Prostate Cancer Cells with Enhanced Cell-Cycle-Related Transcription and Attenuated Androgen Response[J]. Cancer Res,2018,78(4):853-864.
[12] Djamgoz MB, Coombes RC, Schwab A.Ion transport and cancer: from initiation to metastasis[J]. Philos Trans R Soc Lond B Biol Sci, 2014,369(1638): 20130092.
[13] Cho E A, Kim E J, Kwak S J, et al.cAMP signaling inhibits radiation-induced ATM phosphorylation leading to the augmentation of apoptosis in human lung cancer cells[J]. Mol Cancer, 2014, 13:36.
[14] Follin-Arbelet V, Torgersen ML, Naderi EH, et al.Death of multiple myeloma cells induced by cAMP-signaling involves downregulation of Mcl-1 via the JAK/STAT pathway[J]. Cancer Lett, 2013,335(2): 323-331.
[15] Willder JM, Heng SJ, McCall P, et al. Androgen receptor phosphorylation at serine 515 by Cdk1 predicts bioche-mical relapse in prostate cancer patients[J]. Br J Cancer,2013,108(1):139-148.
[16] Tsaur I, Makarević J, Hudak L, et al.The cdk1-cyclin B complex is involved in everolimus triggered resistance in the PC3 prostate cancer cell line[J]. Cancer Lett,2011, 313(1):84-90.
[17] Chang WL, Yu CC, Chen CS, et al.Tubulin-binding agents down-regulate matrix metalloproteinase-2 and -9 in human hormone-refractory prostate cancer cells - a critical role of Cdk1 in mitotic entry[J]. Biochem Pharmacol,2015,94(1):12-21.
[18] Chen Z, Zhang C, Wu D, et al.Phospho-MED1-enhanced UBE2C locus looping drives castration-resistant prostate cancer growth[J]. EMBO J,2011,30(12):2405-2419.
[19] de Resende MF, Vieira S, Chinen LT,et al. Prognostication of prostate cancer based on TOP2A protein and gene assessment: TOP2A in prostate cancer[J]. J Transl Med,2013,11:36.
[20] Zhou W, Yang Y, Xia J, et al.NEK2 induces drug resistance mainly through activation of efflux drug pumps and is associated with poor prognosis in myeloma and other cancers[J]. Cancer Cell,2013,23(1):48-62.
[21] Sugimasa H, Taniue K, Kurimoto A, et al.Heterogeneous nuclear ribonucleoprotein K upregulates the kinetochore complex component NUF2 and promotes the tumorigenicity of colon cancer cells[J]. Biochem Biophys Res Commun,2015,459(1):29-35.
[22] Subhash VV, Tan SH, Tan WL, et al.GTSE1 expression represses apoptotic signaling and confers cisplatin resistance in gastric cancer cells[J]. BMC Cancer,2015, 15: 550.
[23] Fristrup N, Birkenkamp-Demtr?der K, Reinert T, et al. Multicenter validation of cyclin D1, MCM7, TRIM29, and UBE2C as prognostic protein markers in non-muscle-invasive bladder cancer[J]. Am J Pathol,2013,182(2):339-349.
[24] Psyrri A, Kalogeras KT, Kronenwett R, et al.Prognostic significance of UBE2C mRNA expression in high-risk early breast cancer. A Hellenic Cooperative Oncology Group (HeCOG) Study[J]. Ann Oncol,2012,23(6):1422-1427.
[25] Zhang Z, Liu P, Wang J, et al.Ubiquitin-conjugating enzyme E2C regulates apoptosis-dependent tumor progression of non-small cell lung cancer via ERK pathway[J]. Med Oncol,2015,32(5):149.
[26] Guo G, Sun X, Chen C, et al.Whole-genome and whole-exome sequencing of bladder cancer identifies frequent alterations in genes involved in sister chromatid cohesion and segregation[J]. Nat Genet,2013,45(12):1459-1463.
[27] de Voer RM, Geurts van Kessel A, Weren RD, et al. Germline mutations in the spindle assembly checkpoint genes BUB1 and BUB3 are risk factors for colorectal cancer[J]. Gastroenterology,2013,145(3):544-547.
[28] Nyati S, Schinske-Sebolt K, Pitchiaya S, et al. The kinase activity of the Ser/Thr kinase BUB1 promotes TGF-β signaling[J]. Sci Signal,2015,8(358):ra1.
[29] Sun X, Clermont PL, Jiao W, et al.Elevated expression of the centromere protein-A(CENP-A)-encoding gene as a prognostic and predictive biomarker in human cancers[J]. Int J Cancer,2016,139(4):899-907.
[30] Liang Y, Ahmed M, Guo H, et al.LSD1-Mediated Epigenetic Reprogramming Drives CENPE Expression and Prostate Cancer Progression[J]. Cancer Res,2017,77(20):5479-5490.
[31] Aytes A, Mitrofanova A, Lefebvre C, et al.Cross-species regulatory network analysis identifies a synergistic interaction between FOXM1 and CENPF that drives prostate cancer malignancy[J]. Cancer Cell,2014,25(5):638-651.
[32] Imai K, Hirata S, Irie A, et al.Identification of HLA-A2-restricted CTL epitopes of a novel tumour-associated antigen, KIF20A, overexpressed in pancreatic cancer[J]. Br J Cancer,2011,104(2):300-307.
[33] Tomita Y, Yuno A, Tsukamoto H, et al.Identification of promiscuous KIF20A long peptides bearing both CD4+ and CD8+ T-cell epitopes: KIF20A-specific CD4+ T-cell immunity in patients with malignant tumor[J]. Clin Cancer Res,2013,19(16):4508-4520.
[34] Wang J, Guo X, Xie C, et al.KIF15 promotes pancreatic cancer proliferation via the MEK-ERK signalling pathway[J]. Br J Cancer,2017,117(2):245-255.
[35] Taniwaki M, Takano A, Ishikawa N, et al.Activation of KIF4A as a prognostic biomarker and therapeutic target for lung cancer[J]. Clin Cancer Res,2007,13(22 Pt 1): 6624-6631.
[36] Speers C, Zhao S G, Kothari V, et al.Maternal Embryo-nic Leucine Zipper Kinase (MELK) as a Novel Mediator and Biomarker of Radioresistance in Human Breast Cancer[J]. Clin Cancer Res,2016,22(23):5864-5875.
[37] Xia H, Kong S N, Chen J, et al.MELK is an oncogenic kinase essential for early hepatocellular carcinoma recurrence[J]. Cancer Lett,2016,383(1):85-93.
[38] Wang Y, Begley M, Li Q, et al.Mitotic MELK-eIF4B signaling controls protein synthesis and tumor cell survival[J]. Proc Natl Acad Sci U S A,2016,113(35):9810-9815.
[39] Goto Y, Kurozumi A, Arai T, et al.Impact of novel miR-145-3p regulatory networks on survival in patients with castration-resistant prostate cancer[J]. Br J Cancer, 2017, 117(3):409-420.
[40] Wang W Y, Hsu C C, Wang T Y, et al.A gene expression signature of epithelial tubulogenesis and a role for ASPM in pancreatic tumor progression[J]. Gastroenterology,2013,145(5):1110-1120.
[41] Jeffery J, Sinha D, Srihari S, et al.Beyond cytokinesis: the emerging roles of CEP55 in tumorigenesis[J]. Oncogene,2016,35(6):683-690.
[42] Gordon C A, Gulzar Z G, Brooks J D.NUSAP1 expression is upregulated by loss of RB1 in prostate cancer cells[J]. Prostate,2015,75(5):517-526.
[43] Yamada Y, Arai T, Kojima S, et al.Regulation of antitumor miR-144-5p targets oncogenes: Direct regulation of syndecan-3 and its clinical significance[J]. Cancer Sci,2018,109(9):2919-2936.