Abstract
Background:
Given the age relevance of prostate cancer (PCa) and the role of mitochondrial dysfunction (MIDS) in aging, we orchestrated molecular subtypes and identified key genes for PCa from the perspective of MIDS.
Methods:
Cluster analysis, COX regression analysis, function analysis, and tumor immune environment were conducted. We performed all analyses using software R 3.6.3 and its suitable packages.
Results:
CXCL14, SFRP4, and CD38 were eventually identified to classify the PCa patients in The Cancer Genome Atlas (TCGA) database and the Gene Expression Omnibus (GEO) dataset into two distinct clusters. Patients in the cluster 2 had shorter BCR-free survival than those in the cluster 1 in terms of both TCGA database and GEO dataset. We divided the patients from the TCGA database and the GEO dataset into high- and low-risk groups according to the median of MIDS-related genetic prognostic index. For patients in the TCGA database, the biochemical recurrence (BCR) risk in high-risk group was 2.34 times higher than that in low-risk group. Similarly, for patients in the GEO dataset, the risk of BCR and metastasis in high-risk group was 2.35 and 3.04 times higher than that in low-risk group, respectively. Cluster 2 was closely associated with advanced T stage and higher Gleason score for patients undergoing radical prostatectomy or radiotherapy. For patients undergoing radical prostatectomy, the number of CD8+ T cells was significantly lower in cluster 2 than in cluster 1, while cluster 2 had significantly higher stromal score than cluster 1. For patients undergoing radical radiotherapy, cluster 2 had significantly higher level of CD8+ T cells, neutrophils, macrophages, dendritic cells, stromal score, immune score, and estimate score, but showed lower level of tumor purity than cluster 1.
Conclusions:
We proposed distinctly prognosis-related molecular subtypes at genetic level and related formula for PCa patients undergoing radical prostatectomy or radiotherapy, mainly to provide a roadmap for precision medicine.
Introduction
Prostate cancer (PCa) is the most common non-skin malignant tumor diagnosed among American men in 2021, accounting for 26% (). For localized PCa, radical radiotherapy and radical prostatectomy are the preferred treatment options. However, three-quarters of men will experience biochemical recurrence (BCR) after receiving radical treatment without evidence of overt metastatic disease (). There has been no agreement on the definition of BCR (). However, for recurrence patients, the median time to metastasis is 8 years, and the median time from metastasis to death is 5 years (). Due to the lack of prospective randomized trials with a high level of evidence, the best management for BCR has not yet been confirmed since no intervention is currently considered to extend survival, which highlights the importance of personalized therapy and deciding when to start which treatment.
Very little has been known about the cause of PCa, among which aging is the only definite risk factor (). Cellular senescence is a driver of aging and age-related diseases. The increase of age is accompanied by the accumulation of senescent cells in the tissues and the appearance of cellular senescence (). Cell senescence is a cellular stress response caused by irradiation and other macromolecular damage which was once considered to be a tumor suppressor mechanism, but recent studies have shown that senescent cells are metabolically active, and the inflammatory mediators they secrete are called senescence-associated secretory phenotype (SASP) or senescence messaging secretome (). Senescent cells exacerbate inflammation through SASP, which is called “inflammageing” (, ).
Mitochondria have been identified as one of the key regulators of the development of aging phenotypes, especially the pro-inflammatory SASP (). The role of mitochondria in PCa has gradually become clear with a large number of studies on various nuclear-encoded pathways. There is considerable crosstalk between the nucleus and mitochondria through the retrograde signal from the mitochondria to the nucleus and the anterograde signal from the nucleus to the mitochondria through the translocation of cytoplasmic translation proteins to the mitochondria (). Mitochondrial damage has been shown to be involved in the pathophysiology of PCa (), which is a highly hereditary disease (). Changes in the mitochondrial genome have been proven to be related to predictors of tumor proliferation, metastasis, and BCR (). Next-generation sequencing of mitochondrial DNA from 115 men showed a positive correlation between the total burden of acquired mitochondrial DNA variants and the elevated Gleason score at diagnosis and BCR (). Given the age relevance of PCa and the role of mitochondrial dysfunction (MIDS) in aging, we orchestrated molecular subtypes and identified key genes for PCa from the perspective of MIDS, so as to provide a roadmap for the evolution of precision medicine. In addition, we also developed an independent genetic prognosis index to quantify the recurrence risk of patients. Our study has been registered in the ISRCTN registry (No. ISRCTN11560295).
Methods
Data Preparation
For the combination of GSE46602 (), GSE32571 (), GSE62872 (), and GSE116918 () from the Gene Expression Omnibus (GEO) datasets (), R package “inSilicoMerging” () was used and “removeBatchEffect” function of the “limma (version 3.42.2)” package was used to remove the batch effects (Supplementary Figure 1). Subsequently, we extracted the differentially expressed mRNAs between tumor and normal tissues from the GSE46602 (), GSE32571 (), and GSE62872 (), and further conducted the prognosis analysis through log-rank test using the GSE116918 (). Similar methods were used to proceed the PCa data from the TCGA database in the UCSC XENA (). Differentially expressed genes (DEGs) were defined as llogFCl ≥0.4 and p.adj. <0.05. P-value of BCR-free survival or metastasis-free survival was restricted to less than 0.05. MIDS-related genes were obtained from the GeneCards (). The candidate genes were identified through the intersection of DEGs and prognosis-related genes in the GEO and TCGA databases, and the MIDS-related genes. The gene interactions and drug analysis of the candidate genes were performed through the STRING database () and GSCALite () which included drug data of the cancer therapeutics response portal (CTRP) and genomics of drug sensitivity in cancer (GDSC).
Molecular Subtypes and Genetic Prognosis Index
R packages “ConsensusClusterPlus” and “limma” were used to subtyping the patients who underwent radical prostatectomy in the TCGA database or underwent radical radiotherapy in the GSE116918 () through the three candidate genes. The consensus matrix k value denoted the number of clusters. Subsequently, we analyzed the correlations between the clinical parameters and two clusters and prognostic value of the clusters for PCa patients from the TCGA database and GSE116918 (). Gene set enrichment analysis (GSEA) of the two clusters was conducted, and p-value of <0.05 and a false discovery rate (FDR) of <0.25 were considered statistically significant (, ). Besides, we constructed a MIDS-related genetic prognostic index (MDGPI) according to the results of multivariate COX regression analysis for PC patients in the TCGA database to quantify the BCR risk of patients. The MDGPI formula was as follows: risk score = −1.601 + 0.063 ∗ CXCL14 + 0.176 − SFRP4 − 0.095 ∗ CD38. Then, we used the 248 tumor patients in the GSE116918 () to confirm the prognostic value of the MDGPI score.
Tumor Immune Microenvironment (TME) and Checkpoints
We analyzed the tumor immune microenvironment (TME) through the TIMER and ESTIMATE algorithms (, ). In addition, 54 and 47 common immune checkpoints were analyzed for PCa patients from the TCGA and GEO databases, respectively. Comparisons between TME components and immune checkpoints and the two clusters were performed through the Wilcoxon rank sum test. The Spearman analysis was used to explore the relationship between MDGPI and TME components and immune checkpoints. Immune checkpoints, which were differentially expressed between the two clusters and were significantly associated with the BCR-free survival for patients in the TCGA database and GSE116918 (), were identified as well. We presented the flowchart of this study in Figure 1.
Figure 1
Statistical Analysis
We performed all analyses using software R 3.6.3 and its suitable packages. We utilized Wilcoxon test under the circumstance of non-normal data distribution. Variables could be entered into multivariate COX regression analysis if p-value <0.1 in the univariable Cox regression analysis. Survival analysis was conducted through log-rank test and presented as Kaplan–Meier curve. Besides, the Spearman analysis was used to assess the correlations among continuous variables if they did not meet Shapiro–Wilk normality test. Statistical significance was set as two-sided p <0.05. Significant marks were as follows: ns, p ≥0.05; *, p <0.05; **, p <0.01; ***, p <0.001.
Results
Molecular Subtype and its Clinical Values
The GSE46602 (), GSE32571 (), and GSE62872 () had 209 normal and 360 tumor samples, and the GSE116918 () contained 248 PCa patients undergoing radical radiotherapy with complete data of BCR and metastasis. Besides, we also obtained 498 tumor and 52 normal samples of PCa from the TCGA database, among which 430 PCa patients undergoing radical prostatectomy had complete data of BCR. After the intersection of DEGs and prognosis-related genes in the GEO and TCGA databases, and the MIDS-related genes (Figures 2A–C), CXCL14, SFRP4, and CD38 were eventually identified to classify the PCa patients in the TCGA database into two distinct clusters (Figure 2D; consensus matrix k = 2). Moreover, patients in cluster 2 had shorter BCR-free survival than those in cluster 1 (HR: 2.18, 95% CI: 1.29–3.69, p = 0.003; Figure 2E). Similarly, we observed that these three genes could obviously distinguish cluster 2 from cluster 1 for patients undergoing radical radiotherapy in the GSE116918 () (Figure 2F; consensus matrix k = 2), and patients in cluster 2 were more prone to BCR (HR: 2.37, 95% CI: 1.39–4.04, p = 0.001; Figure 2G) and metastasis (HR: 2.94, 95% CI: 1.26–6.84, p = 0.013; Figure 2H) than their counterparts. We divided the patients from the TCGA database and the GSE116918 () into high- and low-risk groups according to the median of MDGPI score. For patients in the TCGA database, the BCR risk in high-risk group was 2.34 times higher than that in low-risk group (95% CI: 1.40–3.91; Figure 2I). Similarly, for patients in the GSE116918 (), the risk of BCR and metastasis in high-risk group was 2.35 and 3.04 times higher than that in low-risk group, respectively (Figures 2J, K). In addition, patients in cluster 2 had significantly higher levels of CXCL14, SFRP4, and MGPI score, and lower level of CD38 than those in cluster 1 for PCa patients from the TCGA database (Figure 2L) and GSE116918 () (Figure 2M).
Figure 2
For patients undergoing radical prostatectomy, we found that cluster 2 was significantly associated with older age, BCR, higher N stage, positive residual tumor, higher Gleason score, and advanced T stage (Table 1). Similarly, for patients undergoing radical radiotherapy, we observed that cluster 2 was significantly related to BCR, metastasis, higher Gleason score, and advanced T stage (Table 2). One rather interesting outcome was that cluster 2 was an independent risk factor for patients undergoing radical radiotherapy (Figure 3A).
Table 1
| Characteristic | Cluster 1 | Cluster 2 | P-value |
|---|---|---|---|
| Samples (n) | 243 | 187 | |
| Age, median (IQR) | 61 (56, 65) | 63 (57, 67) | 0.010 |
| Biochemical recurrence, n (%) | 0.018 | ||
| No | 219 (50.9%) | 153 (35.6%) | |
| Yes | 24 (5.6%) | 34 (7.9%) | |
| N stage, n (%) | <0.001 | ||
| N0 | 182 (48.5%) | 124 (33.1%) | |
| N1 | 23 (6.1%) | 46 (12.3%) | |
| Residual tumor, n (%) | 0.003 | ||
| No | 170 (40.6%) | 103 (24.6%) | |
| Yes | 68 (16.2%) | 78 (18.6%) | |
| Gleason score (GS), n (%) | <0.001 | ||
| GS = 6 | 30 (7%) | 9 (2.1%) | |
| GS = 7 | 149 (34.7%) | 57 (13.3%) | |
| GS = 8 | 29 (6.7%) | 30 (7%) | |
| GS = 9 | 35 (8.1%) | 91 (21.2%) | |
| T stage, n (%) | <0.001 | ||
| T2 | 111 (26.2%) | 44 (10.4%) | |
| T3 | 128 (30.2%) | 133 (31.4%) | |
| T4 | 0 (0%) | 8 (1.9%) |
The correlations between clinical indicators and clusters in the TCGA database.
IQR, interquartile range; GS, Gleason score.
Table 2
| Characteristic | Cluster 1 | Cluster 2 | P-value |
|---|---|---|---|
| Samples (n) | 136 | 112 | |
| Age, median (IQR) | 67 (64, 72) | 69 (62, 73) | 0.632 |
| T stage, n (%) | <0.001 | ||
| T1 | 39 (17.5%) | 12 (5.4%) | |
| T2 | 42 (18.8%) | 34 (15.2%) | |
| T3 | 41 (18.4%) | 51 (22.9%) | |
| T4 | 0 (0%) | 4 (1.8%) | |
| Gleason score (GS), n (%) | <0.001 | ||
| GS = 6 | 37 (14.9%) | 5 (2%) | |
| GS = 7 | 60 (24.2%) | 39 (15.7%) | |
| GS = 8 | 26 (10.5%) | 26 (10.5%) | |
| GS = 9 | 13 (5.2%) | 42 (16.9%) | |
| Biochemical recurrence, n (%) | 0.005 | ||
| No | 115 (46.4%) | 77 (31%) | |
| Yes | 21 (8.5%) | 35 (14.1%) | |
| Metastasis, n (%) | 0.041 | ||
| No | 129 (52%) | 97 (39.1%) | |
| Yes | 7 (2.8%) | 15 (6%) |
The correlations between clinical indicators and clusters in the GSE116918 ().
IQR, interquartile range; GS, Gleason score.
Figure 3
TME and Immune Checkpoints Analysis
For patients undergoing radical prostatectomy, the number of CD8+ T cells was significantly lower in cluster 2 than cluster 1 (p <0.001), while cluster 2 had significantly higher stromal score than cluster 1 (Figure 3B). Moreover, MDGPI score was closely associated with CD4+ T cells (r: 0.16), CD8+ T cells (r: −0.1), macrophages (r: 0.13), dendritic cells (r: 0.18), stromal score (r: 0.37), immune score (r: 0.19), and estimate score (r: 0.32) (Figure 3C). For patients undergoing radical radiotherapy, cluster 2 had significantly higher level of CD8+ T cells (p = 0.002), neutrophils (p <0.001), macrophages (p <0.001), dendritic cells (p <0.001), stromal score (p <0.001), immune score (p <0.001), and estimate score (p <0.001), but showed lower level of tumor purity than cluster 1 (p <0.001) (Figure 3D). In addition, MDGPI score showed significantly correlations with CD8+ T cells (r: 0.23), neutrophils (r: 0.35), macrophages (r: 0.31), dendritic cells (r: 0.35), stromal score (r: 0.53), immune score (r: 0.36), estimate score (r: 0.47), and tumor purity (r: −0.47) (Figure 3E).
In terms of immune checkpoints, 23 and 18 checkpoints were significantly differentially expressed between cluster 2 and cluster 1 for PCa patients from the TCGA database (Figure 3F) and the GSE116918 () (Figure 3G), respectively. For patients from the TCGA database, MDGPI was highly associated with CTLA4 (r: 0.23), CD276 (r: 0.12), KLRD1 (r: −0.11), CD44 (r: −0.19), IDO2 (r: 0.15), TNFRSF4 (r: 0.17), NRP1 (r: 0.19), TNFRSF18 (r: 0.28), TNFSF18 (r: 0.21), CD80 (r: 0.33), CD86 (r: 0.32), CD40 (r: −0.13), CD70 (r: 0.17), LAIR1 (r: 0.28), HAVCR2 (r: 0.33), and BTNL2 (r: −0.21) (Figure 3H). For patients from the GSE116918 (), MDGPI was closely related to BTNL2 (r: −0.18), CD2 (r: 0.20), CD226 (r: 0.14), CD40 (r: −0.24), CD70 (r: −0.22), CD86 (r: 0.26), CD96 (r: 0.26), HAVCR2 (r: 0.31), ICOSLG (r: −0.23), LAG3 (r: −0.18), LAIR1 (r: 0.22), LAYN (r: 0.13), LGALS9 (r: 0.22), SIGLEC15 (r: −0.20), TIGIT (r: −0.17), TNFRSF18 (r: −0.26), TNFRSF8 (r: −0.21), TNFSF14 (r: −0.21), and VTCN1 (r: −0.19) (Figure 3I). Among the above genes (Figures 3F–I), patients who had higher expression of CD96 (HR: 1.79, 95% CI: 1.06–3.03) in the GSE116918 (), and higher level of TNFSF18 (HR: 2.31, 95% CI: 1.38–3.86), CD80 (HR: 1.75, 95% CI: 1.05–2.93), and CD160 (HR: 2.22, 95% CI: 1.32–3.72) in the TCGA database were more prone to BCR than their counterparts (Figure 3J).
Function and Drug Analysis
In order of the predicted scores from the highest to the lowest in the STRING database (), the predicted functional partners of CD38 were PECAM1, CBL, NAMPT, NMNAT1, NMNAT2, NMNAT3, ENPP1, NNMT, ENPP3, and FCGR3A (Figure 4A); for CXCL14, the predicted interaction genes were CXCR4, CXCL12, CXCR3, CXCR2, CCR2, CCR1, CXCR5, CCR7, CXCR1, and CCR5 (Figure 4B); for SFRP4, the predicted interaction genes were WNT2, WNT3A, WNT7A, WNT8A, WNT1, WNT5A, WNT2B, WNT16, WNT4, and WNT10B (Figure 4C). For patients undergoing radical prostatectomy, the highly enriched pathways in cluster 2 were cell cycle, mismatch repair, spliceosome, oocyte meiosis, nucleotide excision repair, base excision repair, homologous recombination, and RNA degradation (Figure 4D). For patients undergoing radical radiotherapy, the highly enriched pathways in cluster 2 were extracellular matrix (ECM) receptor interaction, cell cycle, TGF beta signaling pathway, antigen processing and presentation, Toll like receptor signaling pathway, complement and coagulation cascades, Wnt signaling pathway, chronic myeloid leukemia, lysosome, Notch signaling pathway, circadian rhythm mammal, Fc gamma R-mediated phagocytosis, colorectal cancer, P53 signaling pathway, focal adhesion, apoptosis, small cell lung cancer, neurotrophin signaling pathway, and pancreatic cancer (Figure 4E).
Figure 4
In terms of drug analysis, CTRP drug sensitivity showed that CR-1-31B, Merck60, SB-743921, SR-II-138A, decitabine, leptomycin B, and necrosulfonamide were potentially sensitive to CXCL14, SFRP4, and CD38 (Figure 4F), while GDSC drug sensitivity showed that methotrexate and vorinostat were potentially sensitive to CXCL14, SFRP4, and CD38 (Figure 4G).
Discussion
Although the 5-year relative survival rate of PCa is as high as 98% in America, its estimated death toll is second only to lung cancer and its long-term decline in cancer mortality since 1993 has stopped (). The global population of individuals over 65 years old is growing rapidly and about 20% of the world population will be aged 65 or older by 2030 (). PCa, as an aging-related cancer with a high incidence in men over 65 (), will affirmatively attract considerable attention with an aging population worldwide. Mitochondria are highly evolved organelles that govern energy production, distribution and biosynthesis (). MIDS is one of the typical phenotypes of aging, which could lead to reactive oxygen species-(ROS) driven lipid damage, deposits, and lipofuscin accumulation (). Except for oxidation, lipid-derived aldehyde byproducts, such as 4-hydroxy-2-nonenal have been reported in senescent cells (, ). Furthermore, increased senescence has been observed in the mouse models of MIDS and elevated oxidative stress (). Decreased NAD+/NADH ratios were observed in senescent cells, the possible mechanism was the accumulated pro-inflammatory M1-like macrophages in metabolic tissues during aging which expressed high level of NAD-consuming enzyme CD38 and thereby reduced tissue NAD levels (, ). Given the important role of MIDS in cellular senescence, it is reasonable to link MIDS to cellular senescence and senescence-related diseases, the most important of which are tumors. From the perspective of cancergenesis, MIDS may be the result of mutations in oncogenes and tumor suppressor genes, since changes in the expression levels of oncogenes and tumor suppressor genes such as TP53, which is discussed in next paragraph, and also bcl-2, HIF-1α may directly affect mitochondrial respiration and metabolism, resulting in MIDS (, ). In prostate, zinc in normal prostate epithelial cells slows the tricarboxylic acid cycle (TCA) and ATP production, whereas PCa cells observe TCA and ETC activation (, ). Despite cancergenesis, there is growing evidence that MIDS is associated with cancer cell survival, proliferation, recurrence and metastasis, the strongest evidence of which attributes to mitochondrial-derived ROS (mROS) (). The binding of mitochondrial to SUMO-deficient hexokinase 2 (HK2) enhanced glucose consumption and lactate production, along with reduced mitochondrial respiration, resulted in PCa cell proliferation (). Oxidative phosphorylation-related sulfite oxidase in mitochondria was found to be associated with BCR in post-prostatectomy patients, and the elevated Ki-67 LI score suggested that the mechanism of recurrence may be related to the activation of oxidative phosphorylation and the induction of cell proliferation (). Porporato et al. proposed a hyper-invasive and hyper-metastatic tumor cell phenotype centered on multiple mitochondrial pathways, including ETC overload or partial ETC inhibition and increased succinate and superoxide production, with protein tyrosine kinases Src and Pyk2 as downstream effectors (). Furthermore, understanding the relationship between MIDS and cancer also opens up therapeutic opportunities, for example, new evidence suggests that mitochondrial remodeling is critical in apoptosis and programmed death (–), suggesting that MIDS may be a marker for cancer detection and a target for treatment, especially radiation and chemotherapy, which achieves its therapeutic purpose by inducing apoptosis (). Treatment with the mitochondria-specific superoxide scavenger mitoTEMPO can prevent metastasis in vitro (). 5-(4-methoxyphenyl)-3H-1,2-dithiole-3-thione (AOL), a member of a new class of mROS inhibitors, reduces steady-state cellular ROS levels in human lung cancer cells, expressing anticancer properties (). MIDS, while being a metabolic marker of cancer cells (), establishes a pathway for drug resistance in tumor cells. For example, biguanides, tigecycline and gamitinib inhibit tumor cell energy synthesis by reducing the mitochondrial electron transport chain (ETC), but the upregulation of glycolytic genes compensates for the lack of ATP production, resulting in drug insensitivity (–). In addition, mitochondrial metabolism-related enzymes such as SUMO-deficient HK2 and mitochondrial 2,4-dienoyl-CoA reductase (DECR1) mentioned above also associated with resistance to docetaxel () and resistance to bicalutamide, apalutamide, or enzalutamide (). Radioresistance is an adaptive response to radiation-induced damage by altering several cellular processes that sustain tumor growth. Mitochondria and metabolic reprogramming have been implicated in many cellular processes involved in radioresistance (, ). For example, enzymes important in base excision repair are localized to mitochondria or actively transported to mitochondria (). In addition, cell cycle, oncogenes, tumor suppressor genes, autophagy, cellular metabolism, and ROS are also sites for mitochondria-mediated radioresistance (). For PCa, lactate dehydrogenase A (LDHA), a major metabolic enzyme that produces lactate, is an enzyme that has been shown to be closely related to the glycolytic pathway and PCa radioresistance, and LDHA-targeted therapy combined with radiotherapy can improve the radiosensitivity of radioresistant PCa cells (). It is clear that PCa undergoes metabolic reprogramming, in this process, MIDS is indeed cross-linked with the occurrence, development, recurrence, metastasis and treatment resistance.
Pelicano et al. () found that MIDS and ROS imbalance promoted breast cell motility and the mechanism was that overexpressed CXCL14 could cell motility through elevation of cytosolic Ca (2+) by binding to the inositol 1,4,5-trisphosphate receptor on the endoplasmic reticulum. Besides, SFRP4 was an important Wnt signaling antagonist, and activation of SFRP4 could lead to Wnt signaling suppression and histone modification in PCa stem cells and thereby sensitized tumor cells to chemotherapeutic drugs, enhancing cell death (). The decrease in mitochondrial ATP could reduce calcium uptake into the endoplasmic reticulum, leading to endoplasmic reticulum stress and to impaired Wnt signaling; in turn, the recovery of the ATP level or the inhibition of endoplasmic reticulum stress restored Wnt activity (). Thus, CXCL14, CD38, and SFRP4 were closely associated with MIDS. In this study, using the above three genes, we firstly proposed distinct prognosis-related molecular subtypes from the fresh perspective of MIDS for PCa patients undergoing radical prostatectomy or radiotherapy. Moreover, the molecular subtype was highly associated with the T stage and Gleason score, both of which were closely related to the prognosis of PCa (). In addition, this classification was an independent risk factor for patients undergoing radical radiotherapy.
It is worth noting that we have observed different or even opposite results for the infiltrated immune cells under the two different treatments for PCa patients. This requires us to critically look at the role of inflammation in PCa progression. For patients undergoing radical radiotherapy, various immune cells, namely, CD8+ T cells, neutrophils, macrophages, and dendritic cells are more enriched in cluster 2, the group with a worse prognosis, and the tumor purity is lower. Besides, MDGPI had highly positive correlations with these immune cells. Radiotherapy can trigger and induce inflammation/immune response through factors such as DNA damage, cell death and senescence, immune cell response, cellular stress, hypoxia, and tumor antibodies (). Oxidative stress and DNA damage caused by radical radiotherapy are both considered the initiation events of PCa (, ). Furthermore, for PCa patients with severe inflammation, it was observed that their recurrence-free survival was shorter (). In addition, Schoenfeld et al. found that the single nucleotide polymorphisms of RNASEL, a gene implicated in inflammation, significantly reduces the risk of BCR in patients with radical radiotherapy, but no significant impact on patients with radical prostatectomy (). Research has observed that increased T cell density has been associated with PTEN loss and poorer outcome in African American men with PCa, with lower BCR-free survival, which is partly consistent with our research results (, ). This also explains why in our study, the stromal score and immune score are higher, and cluster 2 with lower tumor purity has a worse prognosis for PCa patients after radical radiotherapy. However, for patients undergoing radical prostatectomy, we did not observe any statistically significant results except for CD8+ T cells between cluster 2 and cluster 1, and a negative correlation between MDGPI and CD8+ T cells was detected. CD8+ T cells are the most powerful effectors in the anticancer immune response (), and the reduced CD8+ T cells might contribute to the worse prognosis of cluster 2 undergoing radical prostatectomy. An interesting finding was the completely opposite trends in CD8+ T cells between radical prostatectomy and radical radiotherapy groups. For radical prostatectomy, the poor prognosis group (cluster 2) showed low levels of CD8+ T cells, CD8+ T cells were higher in the radiotherapy group. CD8+ effector cells are the main effector cells for targeting antitumor immune response (). In the radical prostatectomy patient population, low levels of CD8+ T cells may predict clearance of tumor antigens and low antitumor immune responses, leading to fewer recurrences and metastases. Studies have found that the proportion of CD8+ cells transiently decreased after the first four weeks of radiotherapy, while the proliferation rate of CD8+ T cells increases at the end of radical radiotherapy for PCa and persists until three months after treatment, and the frequency and function of antigen-specific CD8+ T cells remained stable during treatment (). Radiation has proinflammatory and immunomodulatory effects and, contrary to popular belief, promotes antitumor immune responses, namely, T cell homing and tumor infiltration (, ). Lin et al. found that PD-1 expression began to increase after chemotherapy-induced increases in the number of CD8+ tumor-infiltrating lymphocytes immediately after radiotherapy, suggesting that these CD8+ T cells began to become functionally exhausted (). In addition, we observed a higher positive correlation between MDGPI and stromal score, which indicated the role of stromal components in the progression of PCa patients undergoing radical prostatectomy. Reactive stroma has been used to assess the PCa-specific mortality in diagnostic prostate needle biopsies (). TNFRSF18, CD86, CD40, CD70, LAIR1, HAVCR2, and BTNL2 were common checkpoints which were differentially expressed between the two clusters and associated with MDGPI for patients undergoing radical prostatectomy or radiotherapy. However, TNFRSF18 and CD70 showed opposite results after receiving the different treatments, both of which seemed to be involved in interactions between activated T-lymphocytes (). The paradoxical results of the two genes could partially explain the progression mechanism of PCa patients undergoing different treatment when combined the opposite results of immune-infiltrating cells. Furthermore, CD96, TNFSF18, CD80, and CD160 identified in this study might be the potential targets of PCa due to their prognostic values. CD96 may play a role in the adhesion interaction between activated T cells and NK cells in the late immune response. According to Biograph’s knowledge base, in the context of CD96, PCa ranks 4th among 6021 disease concepts, strongly indicating the potential role of this gene in the development of PCa (). The expression of T-cell immunoglobulin domain and mucin domain-containing molecule 3 (TIM-3) encoded by HAVCR2 in CD4+ and CD8+ T cells of the PCa patients was significantly increased compared with benign prostatic hyperplasia, suggesting that it may affect the development and progression of PCa (). In addition, TIM-3 expression is also related to the poor prognosis of PCa (, ). In the TCGA group, CD80 mRNA expression is correlated with BCR, suggesting that genetic variation and mRNA expression in CD80 may be a predictor and potential target of local PCa (). In fact, several clinical studies targeting CD80 (B7-1) have approved, for example, ipilimumab, an antagonistic monoclonal antibody that binds CD80 on antigen-presenting cells, providing options for PCa patients in the future (). The findings of these immune checkpoints are helpful to the choice of medications for adjuvant treatment of patients undergoing radical prostatectomy or radiotherapy.
GSEA analysis shows that the results of our research are related to a variety of cancers, such as colorectal cancer, small cell lung cancer, etc., which further proves the clinical significance of our molecular subtypes. In addition, it was also found to be related to the p53 signaling pathway. Li et al. found that p53-mediated mitochondrial dysfunction can promote PCa cell apoptosis in vitro (). In addition, abnormal activation of Notch signal has also been shown to be closely related to the occurrence and development of PCa (). The deletion of TP53, a tumor suppressor gene, can regulate mitochondrial respiration by promoting the Warburg effect in cancer cells, while increase the uptake of glucose in cancer cells through the repression of transcription of glucose transporter (GLUT) isoforms 1 and 4 and inhibition of the expression of glycolytic enzymes (–).TP53 gene improves the fidelity of DNA replication and homologous recombination through transcriptional activation of mismatch repair (MMR) genes. Abnormal MMR protein expression may be involved in the progression of PCa (). In addition, cluster 2 was found to be related to lysosome in both TCGA and GEO groups. Lysosomes are known to be involved in a variety of cancer processes. Various risk factors for PCa, such as ionizing radiation and oxidative stress, can activate the activity of lysosomal enzymes, which may cause cancer to occur by destroying proteins and other components of cells (). Meanwhile, abiraterone can inhibit the proliferation of PCa cells in vitro, thus promoting apoptosis by regulating mitochondrial autophagy (). Basic autophagy genes are ubiquitous in tumors including PCa, and autophagy defects promote tumorigenesis (). On the other hand, tumor cells are also dependent on autophagy, and the loss of autophagy gene inhibits the formation and metastasis of the primary tumor (, ). In this study, we also found some sensitive drugs to the investigated three genes, which needed to be further studied in vivo and in vitro.
As the aging of the global population continues to develop in the coming decades, PCa in elderly men will bring a huge burden of disease. At present, there is no optimal plan for the management of BCR. In this paper, we calculated that the gene prognostic index composed of CXCL14, SFRP4, and CD38, can well predict the pathogenesis of individual patients with PCa after radical prostatectomy and radiotherapy. In this way, from a clinical perspective, a timely warning can be given before the thorny problem of insufficient treatment methods and poor prognosis for BCR and metastasis patients. At the same time, the discovery of the three targeted genes also avoided tedious and expensive whole-genome sequencing. Our research integrates two high-throughput sequencing and microarray sequencing platforms, as mutual verification, the results are more reliable and have strong clinical relevance.
Conclusions
We proposed distinctly prognosis-related molecular subtypes at genetic level and related formula for PCa patients undergoing radical prostatectomy or radiotherapy, mainly to provide a roadmap for precision medicine.
Funding
This program was supported by the National Natural Science Foundation of China (Grant Nos. 81974099, 82170785, 81974098, 82170784), the programs from the Science and Technology Department of Sichuan Province (Grant No. 21GJHZ0246), the Young Investigator Award of Sichuan University 2017 (Grant No. 2017SCU04A17), the Technology Innovation Research and Development Project of Chengdu Science and Technology Bureau (2019-YF05-00296-SN), and the Sichuan University—Panzhihua Science and Technology Cooperation Special Fund (2020CDPZH-4). The funders had no role in study design, data collection or analysis, preparation of the manuscript, or the decision to publish.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
DCF proposed the project, conducted data analysis, interpreted the data, and wrote the manuscript. XS, FCZ, and QX conducted data analysis, interpreted the data. QW and LY supervised the project, and interpreted the data. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Acknowledgments
The results showed here are in whole or part based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.858479/full#supplementary-material
Abbreviations
PCa, Prostate cancer; BCR, biochemical recurrence; SASP, senescence-associated secretory phenotype; ROS, reactive oxygen species; MMR, mismatch repair; TIM-3, T-cell immunoglobulin domain and mucin domain-containing molecule 3; TME, tumor immune microenvironment; MIDS, mitochondrial dysfunction; DEGs, Differentially expressed genes; GEO, Gene Expression Omnibus; CTRP, cancer therapeutics response portal; GDSC, genomics of drug sensitivity in cancer; GSEA, Gene set enrichment analysis; FDR, false discovery rate; MDGPI, MIDS-related genetic prognostic index.
References
1
SiegelRLMillerKDFuchsHEJemalA. Cancer Statistics, 2021. CA Cancer J Clin (2021) 71(1):7–33. doi: 10.3322/caac.21654
2
ScherHISoloKValantJToddMBMehraM. Prevalence of Prostate Cancer Clinical States and Mortality in the United States: Estimates Using a Dynamic Progression Model. PloS One (2015) 10(10):e0139440. doi: 10.1371/journal.pone.0139440
3
CooksonMSAusGBurnettALCanby-HaginoEDD'AmicoAVDmochowskiRRet al. Variation in the Definition of Biochemical Recurrence in Patients Treated for Localized Prostate Cancer: The American Urological Association Prostate Guidelines for Localized Prostate Cancer Update Panel Report and Recommendations for a Standard in the Reporting of Surgical Outcomes. J Urol (2007) 177(2):540–5. doi: 10.1016/S0084-4071(08)70132-5
4
PoundCRPartinAWEisenbergerMAChanDWPearsonJDWalshPC. Natural History of Progression After PSA Elevation Following Radical Prostatectomy. JAMA (1999) 281(17):1591–7. doi: 10.1001/jama.281.17.1591
5
van DeursenJM. The Role of Senescent Cells in Ageing. Nature (2014) 509(7501):439–46. doi: 10.1038/nature13193
6
GorgoulisVAdamsPDAlimontiABennettDCBischofOBishopCet al. Cellular Senescence: Defining a Path Forward. Cell (2019) 179(4):813–27. doi: 10.1016/j.cell.2019.10.005
7
CoppeJPDesprezPYKrtolicaACampisiJ. The Senescence-Associated Secretory Phenotype: The Dark Side of Tumor Suppression. Annu Rev Pathol (2010) 5:99–118. doi: 10.1146/annurev-pathol-121808-102144
8
FranceschiCCampisiJ. Chronic Inflammation (Inflammaging) and Its Potential Contribution to Age-Associated Diseases. J Gerontol A Biol Sci Med Sci (2014) 69 Suppl 1:S4–9. doi: 10.1093/gerona/glu057
9
ChapmanJFielderEPassosJF. Mitochondrial Dysfunction and Cell Senescence: Deciphering a Complex Relationship. FEBS Lett (2019) 593(13):1566–79. doi: 10.1002/1873-3468.13498
10
XiaoJCohenPSternMCOdedinaFCarptenJReamsR. Mitochondrial Biology and Prostate Cancer Ethnic Disparity. Carcinogenesis (2018) 39(11):1311–9. doi: 10.1093/carcin/bgy133
11
ParrRLMillsJHarbottleACreedJMCrewdsonGRegulyBet al. Mitochondria, Prostate Cancer, and Biopsy Sampling Error. Discovery Med (2013) 15(83):213–20.
12
BarberLGerkeTMarktSCPeischSFWilsonKMAhearnTet al. Family History of Breast or Prostate Cancer and Prostate Cancer Risk. Clin Cancer Res (2018) 24(23):5910–7. doi: 10.1158/1078-0432.CCR-18-0370
13
KalsbeekAMChanEFGroganJPetersenDCJaratlerdsiriWGuptaRet al. Mutational Load of the Mitochondrial Genome Predicts Pathological Features and Biochemical Recurrence in Prostate Cancer. Aging (Albany NY) (2016) 8(11):2702–12. doi: 10.18632/aging.101044
14
MortensenMMHøyerSLynnerupASØrntoftTFSørensenKDBorreMet al. Expression Profiling of Prostate Cancer Tissue Delineates Genes Associated With Recurrence After Prostatectomy. Sci Rep (2015) 5:16018. doi: 10.1038/srep16018
15
KunerRFälthMPressinottiNCBraseJCPuigSBMetzgerJet al. The Maternal Embryonic Leucine Zipper Kinase (MELK) Is Upregulated in High-Grade Prostate Cancer. J Mol Med (Berl) (2013) 91(2):237–48. doi: 10.1007/s00109-012-0949-1
16
PenneyKLSinnottJATyekuchevaSGerkeTShuiIMKraftPet al. Association of Prostate Cancer Risk Variants With Gene Expression in Normal and Tumor Tissue. Cancer Epidemiol Biomarkers Prev (2015) 24(1):255–60. doi: 10.1158/1055-9965.EPI-14-0694-T
17
JainSLyonsCAWalkerSMMcQuaidSHynesSOMitchellDMet 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. doi: 10.1093/annonc/mdx637
18
EdgarRDomrachevMLashAE. Gene Expression Omnibus: NCBI Gene Expression and Hybridization Array Data Repository. Nucleic Acids Res (2002) 30(1):207–10. doi: 10.1093/nar/30.1.207
19
TaminauJMeganckSLazarCSteenhoffDColettaAMolterCet al. Unlocking the Potential of Publicly Available Microarray Data Using Insilicodb and Insilicomerging R/Bioconductor Packages. BMC Bioinf (2012) 13:335. doi: 10.1186/1471-2105-13-335
20
GoldmanMJCraftBHastieMRepečkaKMcDadeFKamathAet al. Visualizing and Interpreting Cancer Genomics Data via the Xena Platform. Nat Biotechnol (2020) 38(6):675–8. doi: 10.1038/s41587-020-0546-8
21
StelzerGRosenNPlaschkesIZimmermanSTwikMFishilevichSet al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr Protoc Bioinf (2016) 54:1.30.1–1.30.33. doi: 10.1002/cpbi.5
22
SzklarczykDGableALNastouKCLyonDKirschRPyysaloSet al. The STRING Database in 2021: Customizable Protein-Protein Networks, and Functional Characterization of User-Uploaded Gene/Measurement Sets. Nucleic Acids Res (2021) 49(D1):D605–12. doi: 10.1093/nar/gkaa1074
23
LiuCJHuFFXiaMXHanLZhangQGuoAY. GSCALite: A Web Server for Gene Set Cancer Analysis. Bioinformatics (2018) 34(21):3771–2. doi: 10.1093/bioinformatics/bty411
24
SubramanianATamayoPMoothaVKMukherjeeSEbertBLGilletteMAet 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. doi: 10.1073/pnas.0506580102
25
LiberzonASubramanianAPinchbackRThorvaldsdóttirHTamayoPMesirovJP. Molecular Signatures Database (MSigDB) 3.0. Bioinformatics (2011) 27(12):1739–40. doi: 10.1093/bioinformatics/btr260
26
LiBSeversonEPignonJCZhaoHLiTNovakJet al. Comprehensive Analyses of Tumor Immunity: Implications for Cancer Immunotherapy. Genome Biol (2016) 17(1):174. doi: 10.1186/s13059-016-1028-7
27
YoshiharaKShahmoradgoliMMartínezEVegesnaRKimHTorres-GarciaWet al. Inferring Tumour Purity and Stromal and Immune Cell Admixture From Expression Data. Nat Commun (2013) 4:2612. doi: 10.1038/ncomms3612
28
CostantinoSPaneniFCosentinoF. Ageing, Metabolism and Cardiovascular Disease. J Physiol (2016) 594(8):2061–73. doi: 10.1113/JP270538
29
AbateMFestaAFalcoMLombardiALuceAGrimaldiAet al. Mitochondria as Playmakers of Apoptosis, Autophagy and Senescence. Semin Cell Dev Biol (2020) 98:139–53. doi: 10.1016/j.semcdb.2019.05.022
30
AdemowoOSDiasHKIBurtonDGAGriffithsHR. Lipid (Per) Oxidation in Mitochondria: An Emerging Target in the Ageing Process? Biogerontology (2017) 18(6):859–79. doi: 10.1007/s10522-017-9710-z
31
WileyCDVelardeMCLecotPLiuSSarnoskiEAFreundAet al. Mitochondrial Dysfunction Induces Senescence With a Distinct Secretory Phenotype. Cell Metab (2016) 23(2):303–14. doi: 10.1016/j.cmet.2015.11.011
32
CovarrubiasAJKaleAPerroneRLopez-DominguezJAPiscoAOKaslerHGet al. Senescent Cells Promote Tissue NAD+ Decline During Ageing via the Activation of CD38+ Macrophages. Nat Metab (2020) 2(11):1265–83. doi: 10.1038/s42255-020-00305-3
33
SciacovelliMGuzzoGMorelloVFrezzaCZhengLNanniniNet al. The Mitochondrial Chaperone TRAP1 Promotes Neoplastic Growth by Inhibiting Succinate Dehydrogenase. Cell Metab (2013) 17(6):988–99. doi: 10.1016/j.cmet.2013.04.019
34
SinghRLetaiASarosiekK. Regulation of Apoptosis in Health and Disease: The Balancing Act of BCL-2 Family Proteins. Nat Rev Mol Cell Biol (2019) 20(3):175–93. doi: 10.1038/s41580-018-0089-8
35
DesoukiMMGeradtsJMilonBFranklinRBCostelloLC. Hzip2 and Hzip3 Zinc Transporters Are Down Regulated in Human Prostate Adenocarcinomatous Glands. Mol Cancer (2007) 6:37. doi: 10.1186/1476-4598-6-37
36
FranklinRBCostelloLC. Zinc as an Anti-Tumor Agent in Prostate Cancer and in Other Cancers. Arch Biochem Biophys (2007) 463(2):211–7. doi: 10.1016/j.abb.2007.02.033
37
KongHReczekCRMcElroyGSSteinertEMWangTSabatiniDMet al. Metabolic Determinants of Cellular Fitness Dependent on Mitochondrial Reactive Oxygen Species. Sci Adv (2020) 6(45):eabb7272. doi: 10.1126/sciadv.abb7272
38
ShangguanXHeJMaZZhangWJiYShenKet al. SUMOylation Controls the Binding of Hexokinase 2 to Mitochondria and Protects Against Prostate Cancer Tumorigenesis. Nat Commun (2021) 12(1):1812. doi: 10.1038/s41467-021-22163-7
39
KuroseHNaitoYAkibaJKondoROgasawaraSKusanoHet al. High Sulfite Oxidase Expression Could Predict Postoperative Biochemical Recurrence in Patients With Prostate Cancer. Med Mol Morphol (2019) 52(3):164–72. doi: 10.1007/s00795-018-00214-1
40
PorporatoPEPayenVLPerez-EscuredoJDe SaedeleerCJDanhierPCopettiTet al. A Mitochondrial Switch Promotes Tumor Metastasis. Cell Rep (2014) 8(3):754–66. doi: 10.1016/j.celrep.2014.06.043
41
KasaharaAScorranoL. Mitochondria: From Cell Death Executioners to Regulators of Cell Differentiation. Trends Cell Biol (2014) 24(12):761–70. doi: 10.1016/j.tcb.2014.08.005
42
BenardGRossignolR. Ultrastructure of the Mitochondrion and Its Bearing on Function and Bioenergetics. Antioxid Redox Signal (2008) 10(8):1313–42. doi: 10.1089/ars.2007.2000
43
OteraHMiharaK. Mitochondrial Dynamics: Functional Link With Apoptosis. Int J Cell Biol (2012) 2012:821676. doi: 10.1155/2012/821676
44
TrottaAPChipukJE. Mitochondrial Dynamics as Regulators of Cancer Biology. Cell Mol Life Sci (2017) 74(11):1999–2017. doi: 10.1007/s00018-016-2451-3
45
Dias AmoedoNDardLSarlakSMahfoufWBlanchardWRousseauBet al. Targeting Human Lung Adenocarcinoma With a Suppressor of Mitochondrial Superoxide Production. Antioxid Redox Signal (2020) 33(13):883–902. doi: 10.1089/ars.2019.7892
46
WeinbergSEChandelNS. Targeting Mitochondria Metabolism for Cancer Therapy. Nat Chem Biol (2015) 11(1):9–15. doi: 10.1038/nchembio.1712
47
ZhuZJiangWThompsonMDEcheverriaDMcGinleyJNThompsonHJ. Effects of Metformin, Buformin, and Phenformin on the Post-Initiation Stage of Chemically Induced Mammary Carcinogenesis in the Rat. Cancer Prev Res (Phila) (2015) 8(6):518–27. doi: 10.1158/1940-6207.CAPR-14-0121
48
XuZYanYLiZQianLGongZ. The Antibiotic Drug Tigecycline: A Focus on Its Promising Anticancer Properties. Front Pharmacol (2016) 7:473. doi: 10.3389/fphar.2016.00473
49
ParkHKLeeJELimJJoDEParkSASuhPGet al. Combination Treatment With Doxorubicin and Gamitrinib Synergistically Augments Anticancer Activity Through Enhanced Activation of Bim. BMC Cancer (2014) 14:431. doi: 10.1186/1471-2407-14-431
50
BlommeAFordCAMuiEPatelRNtalaCJamiesonLEet al. 2,4-Dienoyl-CoA Reductase Regulates Lipid Homeostasis in Treatment-Resistant Prostate Cancer. Nat Commun (2020) 11(1):2508. doi: 10.1038/s41467-020-16126-7
51
Lynam-LennonNMaherSGMaguireAPhelanJMuldoonCReynoldsJVet al. Altered Mitochondrial Function and Energy Metabolism Is Associated With a Radioresistant Phenotype in Oesophageal Adenocarcinoma. PloS One (2014) 9(6):e100738. doi: 10.1371/journal.pone.0100738
52
McCannEO’SullivanJMarconeS. Targeting Cancer-Cell Mitochondria and Metabolism to Improve Radiotherapy Response. Transl Oncol (2021) 14(1):100905. doi: 10.1016/j.tranon.2020.100905
53
ZinovkinaLA. Mechanisms of Mitochondrial DNA Repair in Mammals. Biochem (Mosc) (2018) 83(3):233–49. doi: 10.1134/S0006297918030045
54
HaoJGrahamPChangLNiJWasingerVBeretovJet al. Proteomic Identification of the Lactate Dehydrogenase A in a Radioresistant Prostate Cancer Xenograft Mouse Model for Improving Radiotherapy. Oncotarget (2016) 7(45):74269–85. doi: 10.18632/oncotarget.12368
55
PelicanoHLuWZhouYZhangWChenZHuYet al. Mitochondrial Dysfunction and Reactive Oxygen Species Imbalance Promote Breast Cancer Cell Motility Through a CXCL14-Mediated Mechanism. Cancer Res (2009) 69(6):2375–83. doi: 10.1158/0008-5472.CAN-08-3359
56
DeshmukhAArfusoFNewsholmePDharmarajanA. Epigenetic Demethylation of sFRPs, With Emphasis on Sfrp4 Activation, Leading to Wnt Signalling Suppression and Histone Modifications in Breast, Prostate, and Ovary Cancer Stem Cells. Int J Biochem Cell Biol (2019) 109:23–32. doi: 10.1016/j.biocel.2019.01.016
57
CostaRPeruzzoRBachmannMMontàGDVicarioMSantinonGet al. Impaired Mitochondrial ATP Production Downregulates Wnt Signaling via ER Stress Induction. Cell Rep (2019) 28(8):1949–60.e6. doi: 10.1016/j.celrep.2019.07.050
58
ChangAJAutioKARoachM3rdScherHI. High-Risk Prostate Cancer-Classification and Therapy. Nat Rev Clin Oncol (2014) 11(6):308–23. doi: 10.1038/nrclinonc.2014.68
59
McKelveyKJHudsonALBackMEadeTDiakosCI. Radiation, Inflammation and the Immune Response in Cancer. Mamm Genome (2018) 29(11-12):843–65. doi: 10.1007/s00335-018-9777-0
60
SfanosKSYegnasubramanianSNelsonWGDe MarzoAM. The Inflammatory Microenvironment and Microbiome in Prostate Cancer Development. Nat Rev Urol (2018) 15(1):11–24. doi: 10.1038/nrurol.2017.167
61
De MarzoAMPlatzEASutcliffeSXuJGrönbergHDrakeCGet al. Inflammation in Prostate Carcinogenesis. Nat Rev Cancer (2007) 7(4):256–69. doi: 10.1038/nrc2090
62
IraniJGoujonJMRagniEPeyratLHubertJSaintFet al. High-Grade Inflammation in Prostate Cancer as a Prognostic Factor for Biochemical Recurrence After Radical Prostatectomy. Pathologist Multi Center Study Group. Urology (1999) 54(3):467–72. doi: 10.1016/S0090-4295(99)00152-1
63
SchoenfeldJDMargalitDNKasperzykJLShuiIMRiderJREpsteinMMet al. A Single Nucleotide Polymorphism in Inflammatory Gene RNASEL Predicts Outcome After Radiation Therapy for Localized Prostate Cancer. Clin Cancer Res (2013) 19(6):1612–9. doi: 10.1158/1078-0432.CCR-12-2718
64
KaurHBGuedesLBLuJMaldonadoLReitzLBarberJRet al. Association of Tumor-Infiltrating T-Cell Density With Molecular Subtype, Racial Ancestry and Clinical Outcomes in Prostate Cancer. Mod Pathol (2018) 31(10):1539–52. doi: 10.1038/s41379-018-0083-x
65
KaurHBLuJGuedesLBMaldonadoLReitzLBarberJRet al. TP53 Missense Mutation Is Associated With Increased Tumor-Infiltrating T Cells in Primary Prostate Cancer. Hum Pathol (2019) 87:95–102. doi: 10.1016/j.humpath.2019.02.006
66
RaskovHOrhanAChristensenJPGögenurI. Cytotoxic CD8+ T Cells in Cancer and Cancer Immunotherapy. Br J Cancer (2021) 124(2):359–67. doi: 10.1038/s41416-020-01048-4
67
FarhoodBNajafiMMortezaeeK. CD8(+) Cytotoxic T Lymphocytes in Cancer Immunotherapy: A Review. J Cell Physiol (2019) 234(6):8509–21. doi: 10.1002/jcp.27782
68
EckertFSchaedlePZipsDSchmid-HorchBRammenseeHGGaniCet al. Impact of Curative Radiotherapy on the Immune Status of Patients With Localized Prostate Cancer. Oncoimmunology (2018) 7(11):e1496881. doi: 10.1080/2162402X.2018.1496881
69
GanssRRyschichEKlarEArnoldBHämmerlingGJ. Combination of T-Cell Therapy and Trigger of Inflammation Induces Remodeling of the Vasculature and Tumor Eradication. Cancer Res (2002) 62(5):1462–70.
70
McBrideWHChiangCSOlsonJLWangCCHongJHPajonkFet al. A Sense of Danger From Radiation. Radiat Res (2004) 162(1):1–19. doi: 10.1667/RR3196
71
LinLKaneNKobayashiNKonoEAYamashiroJMNickolsNGet al. High-Dose Per Fraction Radiotherapy Induces Both Antitumor Immunity and Immunosuppressive Responses in Prostate Tumors. Clin Cancer Res (2021) 27(5):1505–15. doi: 10.1158/1078-0432.CCR-20-2293
72
SaeterTVlatkovicLWaalerGServollENeslandJMAxcronaKet al. The Prognostic Value of Reactive Stroma on Prostate Needle Biopsy: A Population-Based Study. Prostate (2015) 75(6):662–71. doi: 10.1002/pros.22957
73
SoTIshiiN. The TNF-TNFR Family of Co-Signal Molecules. Adv Exp Med Biol (2019) 1189:53–84. doi: 10.1007/978-981-32-9717-3_3
74
LiekensAMDe KnijfJDaelemansWGoethalsBDe RijkPDel-FaveroJ. BioGraph: Unsupervised Biomedical Knowledge Discovery via Automated Hypothesis Generation. Genome Biol (2011) 12(6):R57. doi: 10.1186/gb-2011-12-6-r57
75
PiaoYJinX. Analysis of Tim-3 as a Therapeutic Target in Prostate Cancer. Tumour Biol (2017) 39(7):1010428317716628. doi: 10.1177/1010428317716628
76
WuJLinGZhuYZhangHShiGShenYet al. Low TIM3 Expression Indicates Poor Prognosis of Metastatic Prostate Cancer and Acts as an Independent Predictor of Castration Resistant Status. Sci Rep (2017) 7(1):8869. doi: 10.1038/s41598-017-09484-8
77
PiaoYRPiaoLZZhuLHJinZHDongXZ. Prognostic Value of T Cell Immunoglobulin Mucin-3 in Prostate Cancer. Asian Pac J Cancer Prev (2013) 14(6):3897–901. doi: 10.7314/APJCP.2013.14.6.3897
78
WangQYeYYuHLinSHTuHLiangDet al. Immune Checkpoint-Related Serum Proteins and Genetic Variants Predict Outcomes of Localized Prostate Cancer, a Cohort Study. Cancer Immunol Immunother (2021) 70(3):701–12. doi: 10.1007/s00262-020-02718-1
79
MadanRAMohebtashMArlenPMVergatiMRauckhorstMSteinbergSMet al. Ipilimumab and a Poxviral Vaccine Targeting Prostate-Specific Antigen in Metastatic Castration-Resistant Prostate Cancer: A Phase 1 Dose-Escalation Trial. Lancet Oncol (2012) 13(5):501–8. doi: 10.1016/S1470-2045(12)70006-2
80
KhanFSinghVKSaeedMKausarMAAnsariIA. Carvacrol Induced Program Cell Death and Cell Cycle Arrest in Androgen-Independent Human Prostate Cancer Cells via Inhibition of Notch Signaling. Anticancer Agents Med Chem (2019) 19(13):1588–608. doi: 10.2174/1871520619666190731152942
81
BansalASoniARaoPSinghLCMishraAKMohantyNKet al. Implication of DNA Repair Genes in Prostate Tumourigenesis in Indian Males. Indian J Med Res (2012) 136(4):622–32.
82
Schwartzenberg-Bar-YosephFArmoniMKarnieliE. The Tumor Suppressor P53 Down-Regulates Glucose Transporters GLUT1 and GLUT4 Gene Expression. Cancer Res (2004) 64(7):2627–33. doi: 10.1158/0008-5472.CAN-03-0846
83
KawauchiKArakiKTobiumeKTanakaN. P53 Regulates Glucose Metabolism Through an IKK-NF-kappaB Pathway and Inhibits Cell Transformation. Nat Cell Biol (2008) 10(5):611–8. doi: 10.1038/ncb1724
84
JiangPDuWMancusoAWellenKEYangX. Reciprocal Regulation of P53 and Malic Enzymes Modulates Metabolism and Senescence. Nature (2013) 493(7434):689–93. doi: 10.1038/nature11776
85
BerkersCRMaddocksODCheungECMorIVousdenKH. Metabolic Regulation by P53 Family Members. Cell Metab (2013) 18(5):617–33. doi: 10.1016/j.cmet.2013.06.019
86
TappelA. Lysosomal and Prostasomal Hydrolytic Enzymes and Redox Processes and Initiation of Prostate Cancer. Med Hypotheses (2005) 64(6):1170–2. doi: 10.1016/j.mehy.2004.11.039
87
HanJZhangJZhangWZhangDLiYZhangJet al. Abiraterone and MDV3100 Inhibits the Proliferation and Promotes the Apoptosis of Prostate Cancer Cells Through Mitophagy. Cancer Cell Int (2019) 19:332. doi: 10.1186/s12935-019-1021-9
88
AitaVMLiangXHMurtyVVPincusDLYuWCayanisEet al. Cloning and Genomic Organization of Beclin 1, a Candidate Tumor Suppressor Gene on Chromosome 17q21. Genomics (1999) 59(1):59–65. doi: 10.1006/geno.1999.5851
89
WeiHWeiSGanBPengXZouWGuanJL. Suppression of Autophagy by FIP200 Deletion Inhibits Mammary Tumorigenesis. Genes Dev (2011) 25(14):1510–27. doi: 10.1101/gad.2051011
90
KimJKunduMViolletBGuanKL. AMPK and mTOR Regulate Autophagy Through Direct Phosphorylation of Ulk1. Nat Cell Biol (2011) 13(2):132–41. doi: 10.1038/ncb2152
91
GuoJYChenHYMathewRFanJStroheckerAMKarsli-UzunbasGet al. Activated Ras Requires Autophagy to Maintain Oxidative Metabolism and Tumorigenesis. Genes Dev (2011) 25(5):460–70. doi: 10.1101/gad.2016311
Summary
Keywords
molecular subtype, prostate cancer, mitochondria dysfunction, biochemical recurrence, radical prostatectomy, radical radiotherapy
Citation
Feng D, Shi X, Zhang F, Xiong Q, Wei Q and Yang L (2022) Mitochondria Dysfunction-Mediated Molecular Subtypes and Gene Prognostic Index for Prostate Cancer Patients Undergoing Radical Prostatectomy or Radiotherapy. Front. Oncol. 12:858479. doi: 10.3389/fonc.2022.858479
Received
20 January 2022
Accepted
08 March 2022
Published
06 April 2022
Volume
12 - 2022
Edited by
Bela Ozsvari, University of Salford, United Kingdom
Reviewed by
Natasha Kyprianou, Icahn School of Medicine at Mount Sinai, United States; Zongbing You, Tulane University, United States
Updates
Copyright
© 2022 Feng, Shi, Zhang, Xiong, Wei and Yang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiang Wei, weiqiang933@126.com; Lu Yang, wycleflue@163.com
†These authors have contributed equally to this work
This article was submitted to Molecular and Cellular Oncology, a section of the journal Frontiers in Oncology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.