- Open Access
A gene prognostic index from cellular senescence predicting metastasis and radioresistance for prostate cancer
Journal of Translational Medicine volume 20, Article number: 252 (2022)
Senescent cells have been identified in the aging prostate, and the senescence-associated secretory phenotype might be linked to prostate cancer (PCa). Thus, we established a cellular senescence-related gene prognostic index (CSGPI) to predict metastasis and radioresistance in PCa.
We used Lasso and Cox regression analysis to establish the CSGPI. Clinical correlation, external validation, functional enrichment analysis, drug and cell line analysis, and tumor immune environment analysis were conducted. All analyses were conducted with R version 3.6.3 and its suitable packages.
We used ALCAM and ALDH2 to establish the CSGPI risk score. High-risk patients experienced a higher risk of metastasis than their counterparts (HR: 10.37, 95% CI 4.50–23.93, p < 0.001), consistent with the results in the TCGA database (HR: 1.60, 95% CI 1.03–2.47, p = 0.038). Furthermore, CSGPI had high diagnostic accuracy distinguishing radioresistance from no radioresistance (AUC: 0.938, 95% CI 0.834–1.000). GSEA showed that high-risk patients were highly associated with apoptosis, cell cycle, ribosome, base excision repair, aminoacyl-tRNA biosynthesis, and mismatch repair. For immune checkpoint analysis, we found that PDCD1LG2 and CD226 were expressed at significantly higher levels in patients with metastasis than in those without metastasis. In addition, higher expression of CD226 significantly increased the risk of metastasis (HR: 3.65, 95% CI 1.58–8.42, p = 0.006). We observed that AZD7762, PHA-793887, PI-103, and SNX-2112 might be sensitive to ALDH2 and ALCAM, and PC3 could be the potential cell line used to investigate the interaction among ALDH2, ALCAM, and the above drugs.
We found that CSGPI might serve as an effective biomarker predicting metastasis probability and radioresistance for PCa and proposed that immune evasion was involved in the process of PCa metastasis.
The world is now in an era of an aging population, which contributed 16% of cancer cases between 2005 and 2015 . By 2030, approximately 20% of the world’s population will be aged 65 or older, with an exponential augmentation in the prevalence of prostate cancer (PCa) because this disease is most common in men 65 and older [2, 3]. Moreover, the prevalence of metabolic disorders apparently increases and further facilitates the morbidity and mortality of PCa [4, 5]. PCa has the highest morbidity and mortality among urothelial malignancies, with an estimated 1.4 million new male cases and 375,000 deaths worldwide in 2020 . Thus, with population aging and increased life expectancy globally, the improvement of prognosis is an increasingly important area for PCa patients.
PCa could be a pernicious disease for patients with intermediate- or high-risk localized and locally advanced cancer who ask for local curative treatment in most cases . The curative treatments of PCa include radical prostatectomy, intensity-modulated radiotherapy, and proton beam therapy with or without hormone therapy [8, 9]. Unfortunately, approximately one-third of such patients develop this disease due to metastasis and radiotherapy resistance . Cellular senescence is a predominant trait of aged organisms, and excessive accumulation of senescent cells in tissues can contribute to the onset and progression of various age-related diseases, including cancer [11, 12]. The features of the senescence phenotype usually consist of the activation of a chronic DNA damage response, the involvement of various cyclin-dependent kinase inhibitors, increased secretion of proinflammatory and tissue-remodeling factors, induction of antiapoptotic genes, altered metabolic rates, and endoplasmic reticulum stress . To date, the extent to which cellular senescence contributes to PCa remains elusive. In view of this, we developed and validated a cellular senescence-related gene prognostic index (CSGPI) to predict metastasis and tumor radioresistance and explored the related changes in the tumor immune microenvironment (TME) for PCa patients undergoing radical radiotherapy.
Data sources and clinical analysis
Our study has been registered in the ISRCTN registry (No. ISRCTN11560295). In total, we obtained GSE32571 , GSE62872 , GSE79021 , and GSE116918  from the Gene Expression Omnibus (http://www.ncibi.nlm.nih.gov/geo/)  to develop CSGPIs related to metastasis. The specific process of combing GEO datasets could be seen in our previous study . The R package “inSilicoMerging” and the “removeBatchEffect” function of R pachkage “limma” were adopted to merge the four datasets and to further remove the batch effect, respectively . Prostate adenocarcinoma data in the TCGA database and GSE21034  were used for external validation of the prognostic value of the CSGPI. Seventy percent of patients in GSE116918  were extracted randomly to internally validate the prognostic value of the CSGPI. Moreover, radiotherapy resistance was tested using GSE53902 . GSE32571 , GSE62872 , and GSE79021  were used to identify differentially expressed genes (DEGs) and tumor-related genes using weighted gene coexpression network analysis (WGCNA). DEGs were defined by llogFCl ≥ 0.4 and p.adj. < 0.05. Tumor-related genes were defined by a coefficient ≥ 0.3 and p.adj. < 0.001. Cellular senescence-related genes were obtained from the GeneCards database . Subsequently, we performed an intersection of tumor-related genes, DEGs, and cellular senescence-related genes to determine the candidate genes (Fig. 1). A total of 248 samples with complete clinical data in GSE116918  were used to determine the definitive genes through Lasso and Cox regression analysis (Fig. 1). Thereafter, we established a formula for risk stratification: CSGPI risk score = 0.97428*ALCAM-0.85073*ALDH2. We divided the patients into high- and low-risk groups based on the median CSGPI risk score and tested its clinical correlation and prognostic value for metastasis-free survival (MFS) in GSE116918 . External validation using the TCGA database, GSE21034  and GSE53902 , and internal validation were also conducted (Fig. 1). In addition, to further determine the prognostic role of the CSGPI risk score for PCa patients, we stratified the 248 patients according to the latest European Association of Urology (EAU) and National Comprehensive Cancer Network (NCCN) guidelines [24, 25] and compared these factors using Cox regression analysis in terms of MFS.
Protein–protein interaction and competing endogenous RNA (ceRNA) network
We used the GeneMANIA  database to explore the possible links of gene interactions. In addition, we determined the long noncoding RNAs (lncRNAs) that were differentially expressed and associated with MFS, and we further established the ceRNA network according to LncBase  and the miWalk database .
Functional enrichment analysis
Gene Ontology (GO) referred to three aspects as follows: biological process, cell composition, and molecular function. We conducted GO and Kyoto Encyclopedia of Genes and Genome (KEGG) analyses to explore the potential biological functions and signaling pathways of the candidate genes using the package R “clusterProfiler”. We divided the patients in GSE116918  into high- and low-risk groups based on the median CSGPI risk score. Subsequently, gene set enrichment analysis (GSEA) was performed using GSEA software (version 3.0) (http://www.gsea-msigdb.org) . “h.all.v7.4.symbols.gmt” and “c2.cp.kegg.v7.4.symbols.gmt” from the molecular signature database  were used to detect pathways and molecular mechanisms. Considering the gene expression profile and risk groups, the minimum gene set was 5, and the maximum was 5000. p < 0.05 and false discovery rate < 0.25 were considered statistically significant.
DNA mismatch repair (MMR) gene mutation and DNA methylation analysis
The MMR genes and methyltransferases were obtained from a previous study . We analyzed the correlation between these genes and the CSGPI score through Spearman analysis.
TME, drug, and cell line analysis
Seventeen common immune checkpoint genes were used for the analysis in this study. In addition, we used the Xcell  algorithm to analyze the TME through the package R “IOBR” . Differential expression and Spearman analyses of these parameters were performed. We analyzed the drug sensitivity of ALCAM and ALDH2 through GSCALite . Subsequently, we analyzed PCa-related cell lines of ALCAM and ALDH2 and the possible sensitive drugs through the canSAR database .
All analyses were conducted with R version 3.6.3 and its suitable packages. Cytoscape 3.8.2  was used to establish the ceRNA network. Normality tests were performed using the Shapiro–Wilk method, and when the sample did not conform to a normal distribution, Spearman’s correlation analysis was conducted to describe the correlation between quantitative variables. The significance of two groups of samples was tested by the Wilcoxon test and Kruskal–Wallis test for three or more groups. We also carried out survival analysis using the log rank test. Only variables that were statistically significant in the univariable Cox regression analysis were included in the multivariable Cox regression models. ROC curves were generated using the R packages “timeROC” and “pROC”. Each outcome was regarded as statistically significant with a two-sided p value of < 0.05. Significant mark: ns, p ≥ 0.05; *, p < 0.05; **, p < 0.01; ***, p < 0.001.
Development of CSGPI and its clinical value
We detected 64 candidate genes through the intersection of tumor-related genes, DEGs, and cellular senescence-related genes (Fig. 2A–C). A total of 248 tumor samples in GSE116918  were used to identify prognostic genes. Subsequently, we identified ALCAM and ALDH2 as independent prognostic genes through Lasso and Cox regression analyses (Fig. 2D, E) and further established the CSGPI score based on the following formula: CSGPI risk score = 0.97428*ALCAM-0.85073*ALDH2. The Sanky plot showed clinical indicators and CSGPI scores. We divided the patients into high- and low-risk groups according to the median CSGPI score. We confirmed that the CSGPI score could be used as an independent risk factor for metastasis-free survival (MFS) (HR: 9.787, 95% CI 2.213–43.286, p = 0.003; Fig. 2F) and MFS after biochemical recurrence (BCR) (HR: 6.334, 95% CI 1.479–27.139, p = 0.013; Fig. 2G). Furthermore, the prognostic role of the CSGPI was confirmed again when compared to the EAU and NCCN risk groups. In addition, we observed an increasing trend in the CSGPI score with increasing Gleason score (Fig. 2H), T stage (Fig. 2I), and the presence of BCR (Fig. 2J). Furthermore, patients in the high-risk group had a higher risk of BCR than those in the low-risk group (HR: 2.20, 95% CI 1.30–3.72, p = 0.004; Fig. 2K).
We observed that high-risk patients experienced a higher risk of metastasis than their counterparts regardless of MFS (HR: 10.37, 95% CI 4.50–23.93, p < 0.001; Fig. 3A) or MFS after BCR (HR: 6.26, 95% CI 2.65–14.77, p = 0.004; Fig. 3B). The CSGPI score had better diagnostic accuracy for MFS (Fig. 3C, D) and MFS after BCR (Fig. 3E, F). In terms of internal validation, we observed a similar result for MFS (HR: 16.61, 95% CI 6.42–43.00, p < 0.001; Fig. 3G). For external validation of diagnosis, we also detected a consistent outcome using GSE21034  (AUC: 0.823, 95% CI 0.698–0.947; Fig. 3H). For prostate adenocarcinoma in the TCGA database, high-risk patients had a higher risk of metastasis than low-risk patients (HR: 1.60, 95% CI 1.03–2.47, p = 0.038; Fig. 3I). Furthermore, CSGPI had high diagnostic accuracy distinguishing radioresistance from no radioresistance (AUC: 0.938, 95% CI 0.834–1.000; Fig. 3J). We detected that the long noncoding RNA (lncRNA) PART1 was closely related to MFS in the high- and low-risk groups (HR: 0.28, 95% CI 0.12–0.65%; Fig. 3K). Subsequently, we found that lncRNA PART1 might modulate the expression of ALCAM and ALDH2 through interaction with 75 possible miRNAs (Fig. 3L) and further established the ceRNA network (Fig. 3M). Moreover, ALDH2 and ALCAM might interact through coexpression and genetic interactions (Fig. 3N).
GO analysis showed that the 64 candidate genes might engage in the process of negative regulation of epithelial cell proliferation, reproductive system development, tissue remodeling, focal adhesion, collagen-containing extracellular matrix, endoplasmic reticulum chaperone complex, actin binding, scaffold protein binding, and cell adhesion molecule binding (Fig. 4A). KEGG analysis indicated that these genes might be involved in glutathione metabolism, focal adhesion, the MAPK signaling pathway, vascular smooth muscle contraction, drug metabolism, proteoglycans in cancer, pyruvate metabolism, the oxytocin signaling pathway, and PCa (Fig. 4B). GSEA showed that high-risk patients were highly associated with apoptosis, cell cycle, ribosome, base excision repair, aminoacyl-tRNA biosynthesis, and mismatch repair, while low-risk patients were closely associated with glutathione metabolism, arginine and proline metabolism, drug metabolism cytochrome P450, focal adhesion, fatty acid metabolism, arachidonic acid metabolism, regulation of actin cytoskeleton, and so on (Fig. 4C). For the hallmarks, high-risk patients were closely associated with E2F targets (genes encoding cell cycle-related targets of E2F transcription factors), mitotic spindle (genes important for mitotic spindle assembly), MYC targets V1, G2 M checkpoint (genes involved in the G2/M checkpoint, as in progression through the cell division cycle), and MYC targets V2, whereas the low-risk patients were enriched in myogenesis and apical junction (genes encoding components of apical junction complex.), xenobiotic metabolism (genes encoding proteins involved in the processing of drugs and other xenobiotics), estrogen response late, and estrogen response early (Fig. 4D).
TME, drug, and cell line analysis
The MMR gene and methyltransferase analyses showed that the CSGPI score was positively related to MSH2 (r: 0.18), EPCAM (r: 0.27), and DNMT3B (r: 0.13) (Fig. 5A). For immune checkpoint analysis, we found that PDCD1LG2 and CD226 were expressed at significantly higher levels in patients with metastasis than in those without metastasis (Fig. 5B). In addition, higher expression of CD226 significantly increased the risk of metastasis (HR: 3.65, 95% CI 1.58–8.42, p = 0.006; Fig. 5C). Spearman analysis showed a negative correlation between the CSGPI score and CD274 (r: − 0.13), CD47 (r: − 0.16) and CD200 (− 0.21) (Fig. 5D). Compared to patients in the no metastasis group, patients in the metastasis group scored significantly higher for natural killer T cells (p = 0.046) and plasmacytoid dendritic cells (p = 0.041) but scored lower for mast cells (p = 0.03) (Fig. 5E). Spearman analysis showed that CSGPI was positively associated with immature dendritic cells (r: 0.13, p = 0.045), gamma delta T cells (Tgd) (r: 0.17, p = 0.008), and macrophages (r: 0.14, p = 0.030) but negatively related to endothelial cells (r: − 0.14, p = 0.032), preadipocytes (r: − 0.14, p = 0.030), skeletal muscle cells (r: − 0.15, p = 0.012), stromal score (r: − 0.24, p < 0.001), and microenvironment score (r: − 0.15, p = 0.019) (Fig. 5F). In addition, patients aged and over 65 years scored significantly higher in terms of immature dendritic cells (p = 0.01) and mast cells (p = 0.019) than their counterparts (Fig. 5G). Spearman analysis showed that age was significantly associated with immature dendritic cells (r: 0.14, p = 0.02), mast cells (r: 0.13, p = 0.04), and pericytes (r: 0.16, p = 0.11) (Fig. 5H). We observed that AZD7762, PHA-793887, PI-103, and SNX-2112 might be sensitive to ALDH2 and ALCAM (Fig. 5I), and PC3 could be the potential cell line used to investigate the interaction among ALDH2, ALCAM, and the above drugs (Fig. 5J).
PCa has long been a question of great interest in the field of urology. It is well known that metastatic castration-resistant PCa is the leading cause of death, and the prevalence of metastasis is increasing . With an aging population worldwide, the problem will only grow worse. Prior studies have noted that a doubling time of prostate-specific antigen (PSA) ≤ 7.5 months or PSA ≥ 0.5 ng/mL are independent risk factors for MFS ; however, BCR is not a specific indicator of overall survival and PCa-related mortality because a subset of patients only undergo rising PSA levels and will not progress .
PCa is an age-related disease, and senescent cells accumulate with age in all tissues. Although senescent cells cannot replicate, these cells are metabolically active and form an inflammatory microenvironment through the senescence-associated secretory phenotype (SASP) . SASP consists of various proinflammatory mediators, including cytokines, chemokines, growth factors and proteases, and enables senescence of adjacent nonsenescent cells through paracrine pathways and further contributes to the inflammatory microenvironment . Many age-related diseases, such as osteoarthritis and atherosclerosis, have been demonstrated to be correlated with cellular senescence and the accumulation of senescent cells . Senescent cells have also been found in the aging prostate, and previous studies indicate that SASP plays an important role in tumorigenesis despite the tumor suppression of senescence [40, 43, 44]. However, few previous studies have developed a simple and practical genetic biomarker to predict metastasis for patients undergoing radical radiotherapy. In this study, we established and validated a CSGPI score with two genes to predict MFS for PCa patients undergoing radical prostatectomy or radiotherapy. The origins of senescence-triggering mechanisms consist of therapy-induced, oncogene-induced and age-induced senescence . In this study, the high diagnostic accuracy of the CSGPI indirectly demonstrated the link between cellular senescence and radioresistance. The proinflammatory process of senescence derives from DNA damage caused by various stimuli and amplifying effects of SASP [40, 45], which was consistent with our findings that high-risk patients were highly enriched in cell cycle and MMR. MSH2 and EPCAM might be involved in the process of metastasis and radioresistance in patients undergoing radiotherapy. Previous studies have reported that several miRNAs, such as miR-30 families, are dysregulated in PCa and interact with the p16INK4A/Rb pathway, which is associated with cellular senescence [40, 46]. Our ceRNA network indicated that hsa-miR-30a-3p, hsa-miR-30d-3p, and hsa-miR-30e-3p might contribute to the process of cellular senescence.
PDCD1LG2 and CD226 showed significantly higher expression in patients with metastasis, and patients with higher expression levels of CD226 had susceptibility to metastasis. In addition, natural killer T cells and plasmacytoid dendritic cells scored significantly higher in the metastasis group than in their counterparts. In the human body, CD226 is highly expressed on the surface of NK cells and CD8 + T cells and can activate the function of these cells . Thus, we proposed that immune evasion might be involved in the process of metastasis, and the positive correlation between macrophages and CSGPI score supported this opinion as well. Immunosenescence, defined as the changes in the immune system associated with age, has been sought to produce a progressive deterioration in the ability to respond to infections and to develop immunity after vaccination, both of which are associated with a higher mortality rate in the elderly . However, our findings indicated that immunosenescence was not necessarily correlated with age in patients with metastatic PCa. Moreover, we must admit the following limitations. First, gene expression signatures are subject to sampling bias caused by intratumor genetic heterogeneity. In addition, the microenvironment features might be distinct in different tumor regions, such as the tumor core and invasive margin. More importantly, all findings, such as the ceRNA network in this study, still warrant further confirmation.
We found that CSGPI might serve as an effective biomarker predicting metastasis probability and radioresistance for PCa and proposed that immune evasion was involved in the process of PCa metastasis.
Availability of data and materials
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article.
Global Burden of Disease Cancer Collaboration, Fitzmaurice C, Allen C, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 32 cancer groups, 1990 to 2015: a systematic analysis for the global burden of disease study. JAMA Oncol. 2017;3(4):524–48.
Costantino S, Paneni F, Cosentino F. Ageing, metabolism and cardiovascular disease. J Physiol. 2016;594(8):2061–73.
Daniyal M, Siddiqui ZA, Akram M, et al. Epidemiology, etiology, diagnosis and treatment of prostate cancer. Asian Pac J Cancer Prev. 2014;15(22):9575–8.
Fadini GP, Ceolotto G, Pagnin E, et al. At the crossroads of longevity and metabolism: the metabolic syndrome and lifespan determinant pathways. Aging Cell. 2011;10(1):10–7.
Feng D, Yang Y, Wang Z, et al. Inflammatory bowel disease and risk of urinary cancers: a systematic review and pooled analysis of population-based studies. Transl Androl Urol. 2021;10(3):1332–41.
Sung H, Ferlay J, Siegel RL, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.
Teo MY, Rathkopf DE, Kantoff P. Treatment of advanced prostate cancer. Annu Rev Med. 2019;70:479–99.
Schroeck FR, Jacobs BL, Bhayani SB, et al. Cost of new technologies in prostate cancer treatment: systematic review of costs and cost effectiveness of robotic-assisted laparoscopic prostatectomy, intensity-modulated radiotherapy, and proton beam therapy. Eur Urol. 2017;72(5):712–35.
Feng D, Liu S, Li D, et al. Analysis of conventional versus advanced pelvic floor muscle training in the management of urinary incontinence after radical prostatectomy: a systematic review and meta-analysis of randomized controlled trials. Transl Androl Urol. 2020;9(5):2031–45.
Liu CM, Kao CL, Tseng YT, et al. Ginger phytochemicals inhibit cell growth and modulate drug resistance factors in docetaxel resistant prostate cancer cell. Molecules. 2017;22(9):1477.
López-Otín C, Blasco MA, Partridge L, et al. The hallmarks of aging. Cell. 2013;153(6):1194–217.
Lecot P, Alimirah F, Desprez PY, et al. Context-dependent effects of cellular senescence in cancer development. Br J Cancer. 2016;114(11):1180–4.
Hernandez-Segura A, Nehme J, Demaria M. Hallmarks of cellular senescence. Trends Cell Biol. 2018;28(6):436–53.
Kuner R, Fälth M, Pressinotti NC, et al. The maternal embryonic leucine zipper kinase (MELK) is upregulated in high-grade prostate cancer. J Mol Med. 2013;91(2):237–48.
Penney KL, Sinnott JA, Tyekucheva S, et al. Association of prostate cancer risk variants with gene expression in normal and tumor tissue. Cancer Epidemiol Biomarkers Prev. 2015;24(1):255–60.
Sinnott JA, Peisch SF, Tyekucheva S, et al. Prognostic utility of a new mRNA expression signature of gleason score. Clin Cancer Res. 2017;23(1):81–7.
Jain S, Lyons CA, Walker SM, et al. Validation of a Metastatic Assay using biopsies to improve risk stratification in patients with prostate cancer treated with radical radiation therapy. Ann Oncol. 2018;29(1):215–22.
Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.
Feng D, Shi X, Xiong Q, Zhang F, Li D, Yang L. A gene prognostic index associated with epithelial-mesenchymal transition predicting biochemical recurrence and tumor chemoresistance for prostate cancer. Front Oncol. 2022;11:805571. https://doi.org/10.3389/fonc.2021.805571.
Taminau J, Meganck S, Lazar C, et al. Unlocking the potential of publicly available microarray data using inSilicoDb and inSilicoMerging R/Bioconductor packages. BMC Bioinform. 2012;13:335.
Taylor BS, Schultz N, Hieronymus H, et al. Integrative genomic profiling of human prostate cancer. Cancer Cell. 2010;18(1):11–22.
Cojoc M, Peitzsch C, Kurth I, et al. Aldehyde dehydrogenase is regulated by β-Catenin/TCF and promotes radioresistance in prostate cancer progenitor cells. Cancer Res. 2015;75(7):1482–94.
Stelzer G, Rosen N, Plaschkes I, et al. The GeneCards Suite: from gene data mining to disease genome sequence analyses. Curr Protoc Bioinform. 2016;54:1.30.1-1.30.33.
EAU Guidelines. Edn. presented at the EAU Annual Congress Amsterdam 2022. Arnhem: EAU Guidelines Office; 2022.
Schaeffer E, Srinivas S, Antonarakis ES, et al. NCCN guidelines insights: prostate cancer, version 1.2021. J Natl Compr Canc Netw. 2021;19(2):134–43.
Warde-Farley D, Donaldson SL, Comes O, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38:W214-20.
Paraskevopoulou MD, Vlachos IS, Karagkouni D, et al. DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts. Nucleic Acids Res. 2016;44(D1):D231–8.
Sticht C, De La Torre C, Parveen A, et al. miRWalk: an online resource for prediction of microRNA binding sites. PLoS ONE. 2018;13(10): e0206239.
Subramanian A, Tamayo P, Mootha VK, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.
Liberzon A, Subramanian A, Pinchback R, et al. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27(12):1739–40.
Zhang Z, Zhang X, Huang A. Aggresome-autophagy associated gene HDAC6 is a potential biomarker in pan-cancer, especially in colon adenocarcinoma. Front Oncol. 2021;11: 718589.
Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017;18(1):220.
Zeng D, Ye Z, Shen R, et al. IOBR: multi-omics immuno-oncology biological research to decode tumor microenvironment and signatures. Front Immunol. 2021;12: 687975.
Liu CJ, Hu FF, Xia MX, et al. GSCALite: a web server for gene set cancer analysis. Bioinformatics. 2018;34(21):3771–2.
Mitsopoulos C, Di Micco P, Fernandez EV, et al. canSAR: update to the cancer translational research and drug discovery knowledgebase. Nucleic Acids Res. 2021;49(D1):D1074–82.
Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Weiner AB, Matulewicz RS, Eggener SE, et al. Increasing incidence of metastatic prostate cancer in the United States (2004–2013). Prostate Cancer Prostatic Dis. 2016;19(4):395–7.
Markowski MC, Chen Y, Feng Z, et al. PSA doubling time and absolute PSA predict metastasis-free survival in men with biochemically recurrent prostate cancer after radical prostatectomy. Clin Genitourin Cancer. 2019;17(6):470-475.e1.
Artibani W, Porcaro AB, De Marco V, et al. Management of biochemical recurrence after primary curative treatment for prostate cancer: a review. Urol Int. 2018;100(3):251–62.
Fiard G, Stavrinides V, Chambers ES, et al. Cellular senescence as a possible link between prostate diseases of the ageing male. Nat Rev Urol. 2021;18(10):597–610.
Huda N, Liu G, Hong H, et al. Hepatic senescence, the good and the bad. World J Gastroenterol. 2019;25(34):5069–81.
van Deursen JM. The role of senescent cells in ageing. Nature. 2014;509(7501):439–46.
Castro P, Giri D, Lamb D, et al. Cellular senescence in the pathogenesis of benign prostatic hyperplasia. Prostate. 2003;55(1):30–8.
Coppé JP, Desprez PY, Krtolica A, et al. The senescence-associated secretory phenotype: the dark side of tumor suppression. Annu Rev Pathol. 2010;5:99–118.
de Bono JS, Guo C, Gurel B, et al. Prostate carcinogenesis: inflammatory storms. Nat Rev Cancer. 2020;20(8):455–69.
Valentino A, Reclusa P, Sirera R, et al. Exosomal microRNAs in liquid biopsies: future biomarkers for prostate cancer. Clin Transl Oncol. 2017;19(6):651–7.
Huang Z, Qi G, Miller JS, et al. CD226: an emerging role in immunologic diseases. Front Cell Dev Biol. 2020;8:564.
Aw D, Silva AB, Palmer DB. Immunosenescence: emerging challenges for an ageing population. Immunology. 2007;120(4):435–46.
The results showed here are in whole or part based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga.
This program was supported by the National Natural Science Foundation of China (Grant Nos. 81974099, 82170785, 81974098, 82170784), programs from Science and Technology Department of Sichuan Province (Grant Nos. 21GJHZ0246), Young Investigator Award of Sichuan University 2017 (Grant No. 2017SCU04A17), Technology Innovation Research and Development Project of Chengdu Science and Technology Bureau (2019-YF05-00296-SN), Sichuan University–Panzhihua science and technology cooperation special fund (2020CDPZH-4). The funders had no role in the study design, data collection or analysis, preparation of the manuscript, or the decision to publish.
Ethics approval and consent to participate
The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Consent for publication
The authors have no competing interests to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Dengxiong Li and Xu Shi are co-first authors.
About this article
Cite this article
Feng, D., Li, D., Shi, X. et al. A gene prognostic index from cellular senescence predicting metastasis and radioresistance for prostate cancer. J Transl Med 20, 252 (2022). https://doi.org/10.1186/s12967-022-03459-8
- Cellular senescence
- Prognostic index
- Prostate cancer
- Tumor immune microenvironment
- Metastasis-free survival
- Immune checkpoint