- Research
- Open access
- Published:
Identification of the cytoplasmic DNA-Sensing cGAS-STING pathway-mediated gene signatures and molecular subtypes in prostate cancer
BMC Cancer volume 24, Article number: 732 (2024)
Abstract
Background
Considering the age relevance of prostate cancer (PCa) and the involvement of the cGAS-STING pathway in aging and cancer, we aim to classify PCa into distinct molecular subtypes and identify key genes from the novel perspective of the cGAS-STING pathway. It is of significance to guide personalized intervention of cancer-targeting therapy based on genetic evidence.
Methods
The 430 patients with PCa from the TCGA database were included. We integrated 29 key genes involved in cGAS-STING pathway and analyzed differentially expressed genes and biochemical recurrence (BCR)-free survival-related genes. The assessments of tumor stemness and heterogeneity and tumor microenvironment (TME) were conducted to reveal potential mechanisms.
Results
PCa patients were classified into two distinct subtypes using AURKB, TREX1, and STAT6, and subtype 1 had a worse prognosis than subtype 2 (HR: 21.19, p < 0.001). The findings were validated in the MSKCC2010 cohort. Among subtype 1 and subtype 2, the top ten mutation genes were MUC5B, DNAH9, SLC5A10, ZNF462, USP31, SIPA1L3, PLEC, HRAS, MYOM1, and ITGB6. Gene set variation analysis revealed a high enrichment of the E2F target in subtype 1, and gene set enrichment analysis showed significant enrichment of base excision repair, cell cycle, and DNA replication in subtype 1. TME evaluation indicated that subtype 1 had a significantly higher level of T cells follicular helper and a lower level of plasma cells than subtype 2.
Conclusions
The molecular subtypes mediated by the cGAS-STING pathway and the genetic risk score may aid in identifying potentially high-risk PCa patients who may benefit from pharmacologic therapies targeting the cGAS-STING pathway.
Introduction
Prostate cancer (PCa) is the most common malignancy among males in about 60% country worldwide, which is the second most frequent cancer and the fifth leading cause of cancer death among males [1]. The precise pathogenesis of PCa remains uncertain, however, it is acknowledged that PCa is a heterogeneous age-related ailment, and its occurrence exhibits a quadratic rise with advancing age [2,3,4]. This observation implies that aging might exert a conceivably pivotal function in the initiation and/or progression of PCa. Aging denotes the outcome of time-dependent aggregation of cellular impairments accompanied by a gradual decline in physiological integrity, resulting in compromised functionality. Furthermore, the cumulative cellular damage during aging may sporadically confer anomalous benefits upon specific cells, ultimately instigating carcinogenesis [5]. Deoxyribonucleic acid (DNA) damage accumulation represents a shared feature of the aging process and assumes a causal function in various premature aging disorders [6, 7]. Despite the existence of an intricate network of DNA repair mechanisms aimed at addressing a considerable portion of nuclear DNA damage induced by external and internal stresses [8], the repair capacity is limited and the unrepaired DNA damage will persist and accumulate with age, eventually leading to aging and death. In addition, the integrity of mitochondrial DNA (mtDNA) is an important part of the genomic stability systems, most mtDNA mutations caused by replication errors occur in the origin of life and increase by polyclonal expansion in aging. The accumulated mtDNA mutations damage respiratory chain dysfunction, leading to mitochondrial dysfunction in different tissues[9]. The progressive deterioration of mitochondrial function during the aging process leads to apoptosis, resulting in the release of mtDNA into the cytoplasm via BAX/BAK macropores or voltage-dependent anion channel (VDAC) pores. Additionally, the mitochondrial permeability transition pore (MPTP) may serve as an alternative pathway for inner mitochondrial membrane traversal of mtDNA, enabling its release from mitochondria to the cytoplasm [10]. Cancer cells frequently exhibit an abundance of cytoplasmic double-stranded DNA (dsDNA) derived from various sources, including the genome, mitochondria, and exogenous origins. Nevertheless, chromosomal instability (CIN) constitutes a major contributor to the presence of cytoplasmic DNA [11].
The occurrence of anomalous cytoplasmic dsDNA serves as damage-related molecular patterns (DAMP) that can be detected by the DNA sensor known as cyclic GMP-AMP synthase (cGAS) within the cytoplasm. The binding of cGAS to dsDNA induces an allosteric effect, facilitating the catalysis of cyclic-GMP-AMP (cGAMP) synthesis. cGAMP functions as a second messenger and robust agonist of stimulator of interferon genes (STING). The interaction between cGAMP and STING stimulates the production of type-I interferons (IFN) and initiates the immune responses of the host [12]. Inflammation exerts significant influence as a potent promoter of tumorigenesis [13,14,15]. The cGAS-STING pathway serves as a mediator of the transcriptional activity encompassing a wide range of molecular processes, including inflammation and the oncogenic progression facilitated by cytokines, chemokines, and growth factors, which promote cellular proliferation, survival, and angiogenesis [16]. The DNA damage response (DDR) acts as a protective mechanism against the replication of DNA lesions by activating multiple cellular events, thereby playing a crucial role in maintaining genomic integrity [17]. CIN stands as a hallmark of human cancer and is closely associated with tumor advancement, therapeutic resistance, and distant metastasis. Cancer cells harboring unstable genomes are susceptible to chromosome missegregation during mitosis, resulting in the formation of micronuclei. These micronuclei possess fragile nuclear envelopes, leading to the exposure of genomic contents to the cytoplasm and subsequently triggering the cGAS-STING pathway, culminating in the production of IFN and inflammatory cytokines [18,19,20]. According to reports, CIN displays a strong correlation with the progression and metastatic burden observed in PCa [21]. Metastatic and castration-resistant PCa cases featuring defects in DDR showcase characteristics of genomic instability, which significantly impact the prognosis of PCa patients [22].
Considering the age-related nature of PCa and the involvement of the cGAS-STING pathway in both aging and cancer, as indicated by our previous study [23], we have conducted an analysis of molecular subtypes and identified key genes for PCa through the lens of the cGAS-STING pathway. This effort aims to provide a roadmap for the advancement of precision medicine in PCa. Furthermore, we have developed an independent genetic prognosis index to quantitatively assess the risk of recurrence in PCa patients.
Methods
Data preparation
We integrated 29 key genes involved in cGAS-STING pathway from previous literatures [11, 24,25,26,27,28,29]. The PCa gene matrix were obtained from our previous study [30]. Differentially expressed genes and biochemical recurrence (BCR)-free survival-related genes were analyzed. Differential expression was defined based on an absolute log fold change value greater than 0.4 and a Padj value smaller than 0.05. Subsequently, we performed an intersection analysis of differentially expressed genes with BCR-related and cGAS-STING-related genes. This allowed us to identify a specific set of genes for grouping 430 PCa patients who underwent radical prostatectomy in the TCGA database using the nonnegative matrix factorization (NMF) algorithm. Furthermore, we established a related risk score based on these genes. MSKCC2010 cohort [31] containing 140 PCa patients undergoing radical prostatectomy were used to externally validated the molecular subtypes and risk score. The prognosis and clinical traits of molecular subtypes were analyzed.
Mutation landscape and functional differences between two subtypes
RNA-sequencing profiles, genetic mutation data, and corresponding clinical information for PCa were obtained from the TCGA database (https://portal.gdc.com). The maftools package in R software was utilized to download and visualize mutation data. Furthermore, a comparison of mutation frequencies between two subtypes was conducted. The baseline features of two TCGA subtypes in prostate cancer patients was shown in Table 1.
Functional analysis involved gene set enrichment analysis utilizing the "c2.cp.kegg.v7.4.symbols.gmt" and "h.all.v7.4.symbols.gmt" datasets from the molecular signatures database [32,33,34]. Gene set enrichment analysis was performed based on gene expression and subtypes, with a minimum gene set defined as 5 and a maximum gene set defined as 5000. Resampling was conducted 1000 times. A p-value of < 0.05 and a false discovery rate of < 0.10 were considered statistically significant. Gene set variation analysis employed a minimum and maximum gene set of 5 and 5000, respectively. The "wilcox.test" function was used to assess the differences in each pathway between the two clusters. A log fold change of 0.4 was considered, and statistical significance was defined as a Padj value < 0.01 and a false discovery rate < 0.01.
Tumor stemness and heterogeneity analyses
The tumor stemness indexes comprised various scores, namely differentially methylated probes-based stemness scores (DMPss), DNA methylation-based stemness scores (DNAss), enhancer elements/DNA methylation-based stemness scores (ENHss), epigenetically regulated DNA methylation-based stemness scores (EREG-METHss), epigenetically regulated RNA expression-based stemness scores (EREG.EXPss), RNA expression-based stemness scores (RNAss), and mRNAsi [35, 36]. Tumor heterogeneity was assessed based on factors such as homologous recombination deficiency (HRD), loss of heterozygosity (LOH), neoantigen (NEO) burden, tumor ploidy, tumor purity, mutant-allele tumor heterogeneity (MATH), tumor mutation burden (TMB), and microsatellite instability (MSI) [37, 38]. The results of the aforementioned indicators were obtained from our previous study [39]. To compare the differences between the two subtypes, the Wilcoxon rank sum test was employed.
Tumor microenvironment assessment
The comprehensive evaluation of the tumor microenvironment and immune components was conducted using the MCPCounter and ESTIMATE algorithms [40]. The TIDE algorithm was employed to predict the potential response to immune checkpoint blockade (ICB) therapy [41]. A high TIDE score indicates a lower efficacy of ICB. The differences in the expression levels of 54 immune checkpoints and the tumor microenvironment scores between the two subtypes were analyzed using the Wilcoxon rank sum test. Figure 1 illustrates the flowchart outlining the methodology employed in this study.
Statistical analysis
Analyses were performed using R software version 3.6.3 and relevant packages. The Wilcoxon test was employed for comparisons involving abnormal distributions. Survival analysis was conducted using the log-rank test, and the results were presented as Kaplan–Meier curves. Statistical significance was defined as a two-sided p-value of less than 0.05. The notation for indicating significance levels was as follows: not significance (ns), p ≥ 0.05; *, p < 0.05; **, p < 0.01; ***, p < 0.001.
Results
Identification of cGAS-STING-mediated molecular subtypes and key genes
Through gene differential expression analysis, we found that three genes including AURKB, TREX1, and STAT6 were with statistical difference (Fig. 2A). Prognosis analyses were performed, which showed different clinical implication among these genes. To be specific, STAT6 (HR:0.53, 95%CI:0.37–0.76), XYLT2 (HR:0.43, 95%CI:0.21–0.90), and TREX1 (HR:0.60, 95%CI:0.38–0.94) were linked to favorable prognosis, but IKBKB (HR:2.04, 95%CI:1.37–3.06) and AURKB (HR:1.37, 95%CI:1.11–1.69) were associated with poor prognosis (Fig. 2B). The differentially expressed genes in TCGA database (TCGA-DEGs) had intersection with cGAS-STING genes and BCR-free survival-related genes. There were 11 genes involving in both cGAS-STING pathway and TCGA-DEGs. There were 5 genes having relationship with both cGAS-STING pathway and BCR-free survival. Thereinto, 3 genes were associated with cGAS-STING pathway, TCGA-DEGs, and BCR-free survival (Fig. 2C). Thus, we conducted subsequent classification analysis using AURKB, TREX1, and STAT6. These genes demonstrated the ability to distinctly classify 430 PCa patients in the TCGA database into two subtypes (Fig. 2D), where subtype 1 had worse prognosis than subtype 2 (HR: 21.19, p < 0.001; Fig. 2E). We further validated our findings using MSKCC2010 cohort and similar results were observed (Fig. 2F-G). Using multivariate Cox analysis, we established cGAS-STING related gene prognostic index (CSRGPI). [CSRGPI = -0.38327171969959*STAT6 + 0.341526882562623*AURKB-0.641776317397796*TREX1.] Based on the interquartile range, PCa patients were stratified into high- and low-risk groups. Analysis of both cohorts revealed a significantly elevated risk of BCR in the high-risk group compared to the low-risk group (Fig. 2H-I).
Mutation landscape and functional analysis
Between subtype 1 and subtype 2, the top ten mutation genes were MUC5B, DNAH9, SLC5A10, ZNF462, USP31, SIPA1L3, PLEC, HRAS, MYOM1 and ITGB6 (Fig. 3A). Gene set variation analysis (GSVA) revealed a significant enrichment of E2F target genes in subtype 1 (Fig. 3B) and it exhibited a substantial enrichment of base excision repair, cell cycle, and DNA replication, as revealed by gene set enrichment analysis (GSEA) (Fig. 3C).
Tumor heterogeneity and stemness and tumor immune microenvironment (TME)
In terms of tumor heterogeneity, subtype 1 had significantly higher levels of HRD and TMB than subtype 2. However, there was no significant difference on levels of LOH, MATH, MSI, Neoantigen, Ploidy, Purity among these two subtypes (Fig. 3D). Likewise, we did not identify obvious distinctions in the indicators of tumor stemness between subtype 1 and subtype 2, such as mRNAsi, DMPss, DNAss, ENHss, et al. (Fig. 3D). For TME evaluation, there was no significant difference on TIDE score and ESTIMATE algorithm (including stromal score, immune score, and ESTIMATE score) (Fig. 3E). However, for the levels of significant immune cells, subtype 1 demonstrated a notable elevation in T cells follicular helper and a decrease in plasma cells compared to subtype 2 (Fig. 3E). Regarding the levels of significant immune checkpoints, subtype 1 exhibited significantly higher levels of ICOS, CTLA4, CD80, LAG3, TNFRSF8, and CD86 compared to subtype 2 (Fig. 3E).
Discussion
There is an exponential increase in the prevalence of PCa due to the fact that approximately additional 1.4 million people had PCa and 375,000 people died for the disease worldwide in 2020 [1]. Although the major determining factors of PCa have not been described, the fact that age and family history are the major determinants has been identified so far [42].
Aging encompasses a temporal progression of functional decline characterized by the gradual deterioration of physiological integrity. This process gives rise to various abnormal human pathologies, including age-related diseases such as cancer, neurodegenerative disorders, diabetes, and cardiovascular ailments. Under physiological circumstances, DNA remains confined within the nucleus and mitochondria, while cytosolic and endolysosomal compartments harbor nucleases capable of DNA degradation. Genomic Instability and mitochondrial dysfunction are two phenotypes of aging, which are the major causes of aberrant release of dsDNA to cytoplasm with age [5].
Endogenous cytoplasmic dsDNA is a major driver of sterile inflammation, also called ‘‘inflamm-aging’’ that regarded as a hallmark of the aging process[43]. Genomic instability is associated with the formation of micronuclei. Micronuclei are small membrane-bounded compartments encapsulated by a nuclear envelope with chromosomes or chromosome fragments. Nevertheless, the nuclear envelope of micronuclei exhibits structural deficiencies, rendering it susceptible to fragmentation upon exposure to both internal and external stressors. Consequently, the DNA content within these micronuclei can be released directly into the cytoplasm. Mitochondria, on the other hand, consist of two distinct membranes, namely the outer and inner membranes, which play vital roles in facilitating intracellular metabolism and energy production. Within the inner mitochondrial matrix, discrete nucleoids organize the mtDNA [44, 45].
Age-related mitochondrial dysfunction is marked by reduced respiratory capacity per individual mitochondrion, accompanied by diminished mitochondrial membrane potential and an associated increase in the generation of reactive oxygen species (ROS). This impairment of mitochondrial function acts as both a cause and a consequence of cellular senescence, prompting the permeabilization of mitochondrial membranes and initiating either apoptosis or necrosis [46, 47]. The proapoptotic proteins BAX/BAK are expressed on the outer mitochondrial membrane, where they assemble into macropores, leading to the induction of mitochondrial outer membrane permeabilization (MOMP) upon sensing apoptotic signals [48]. Following the activation of BAK/BAX and subsequent loss of cytochrome c, the integrity of the mitochondrial network becomes compromised. Additionally, the pores formed by BAX/BAK on the outer mitochondrial membrane gradually expand, leading to the protrusion of the inner mitochondrial membrane, which contains the mtDNA, into the cytoplasm [49]. cGAS is a protein consisting of 520 amino acids that typically remains in an inactive state. It comprises an unstructured and highly basic N-terminus, spanning 160 amino acids, and a globular structural domain composed of 360 amino acids [12].
Recent investigations have revealed that cGAS is not solely confined to the cytoplasm, but also exhibits constitutive and predominant localization within the nucleus. Despite its role as a DNA sensor and interaction with nuclear DNA, the catalytic function of cGAS is impeded due to its strong tethering to the nucleus. This tethering occurs as cGAS binds to a negatively charged acidic region created by histones, utilizing its second DNA-binding site. The high-affinity binding to nucleosomes obstructs the binding of double-stranded DNA and consequently induces the inactivation of cGAS [50]. Nonetheless, sequestering cGAS in an inactive state serves to maintain a low level of soluble cGAS. Upon stimulation by agonistic double-stranded DNA (dsDNA), the high-affinity nucleosome-bound cGAS swiftly replenishes the soluble cGAS pool without requiring de novo protein synthesis. This mechanism plays a crucial role in achieving an optimal dynamic range and facilitating a rapid response in scenarios where soluble cGAS becomes depleted due to the excessive presence of aberrant dsDNA leakage in the cytoplasm [51]. Upon encountering dsDNA, cytoplasmic cGAS undergoes a conformational change, leading to its activation and the subsequent generation of the second messenger molecule, cGAMP. Subsequently, cGAMP binds to STING, which is situated within the endoplasmic reticulum (ER). This binding event triggers a reconfiguration of STING's conformation, facilitating its translocation from the ER to the Golgi apparatus via the ER-Golgi intermediate compartment (ERGIC). Within the Golgi, STING undergoes polymerization and activates TANK-binding kinase 1 (TBK1) [24, 52]. Following activation, TBK1 phosphorylates STING, leading to the recruitment of interferon regulatory factor 3 (IRF3). Subsequently, TBK1 further phosphorylates IRF3, prompting its dimerization and subsequent translocation into the nucleus. This cascade of events initiates the production of interferons (IFNs) and other target genes, including those encoding inflammatory cytokines such as IL6 and IL12. The transcriptional response induced by IRF3 varies across different cell types and depends on the specific triggering of STING as well as the precise in vivo environment [10].
The cGAS-STING pathway has emerged as a significant mechanism involved in inflammation-driven tumorigenesis. The pathway exerts its effects through downstream effector molecules, including TBK1, which has been associated with sustained inflammation and the progression of cancer [53, 54]. Considering the significant involvement of the cGAS-STING pathway in the processes of aging and inflammation, it is plausible to establish a connection between cGAS-STING and age-related diseases as well as inflammation-associated tumors. PCa stands as a notable example of such diseases [55, 56].
Our findings provide evidence that AURKB, STAT6, and TREX1 are among the genes associated with the cGAS-STING pathway, contributing to the BCR-free survival of PCa patients. AURKB, belonging to the Aurora kinase gene family, is a highly conserved gene involved in threonine/serine regulation in mammals. It forms the chromosomal passenger complex (CPC) along with proteins such as Survivin, Borealin, and the Inner Centromere Protein. Throughout mitosis, the CPC dynamically localizes to different regions of the centromere, including the prophase and metaphase centromeres, the cortex and spindle midzone during anaphase, and the mid-body during telophase. Its crucial functions encompass chromosome segregation and cytokinesis processes [57, 58].
A protein complex consisting of AURKB and P53 has been identified, wherein AURKB phosphorylates serine and threonine residues at position 284 of P53. Perturbations in AURKB expression result in inhibiting the transcriptional activity of p53, thereby compromising its role as a tumor suppressor [59]. AURKB exhibits high expression levels in primary human prostate cancer (PCa) and its cell lines. Nuclear localization of AURKB in prostatic intraepithelial neoplasia lesions is associated with clinical staging [60, 61]. Additionally, elevated levels of transforming growth factor β (TGFβ) have been observed in certain advanced cancers, contributing to tumor progression. AURKB can form a complex with TGFβ receptor (TβR), known as the AURKB-TβRI complex, which correlates with the malignancy of PCa and serves as a potential prognostic biomarker for patients at risk of developing aggressive PCa [62]. Addepalli et al. demonstrated that the knockdown of AURKB through RNAi inhibited PCa growth in athymic nude mice [63]. STAT6 is a transcription factor involved in tumorigenesis and the regulation of the tumor microenvironment. Its activation depends on the phosphorylation of conserved tyrosine residues on the receptor by cytokines. Phosphorylated STAT6 forms homodimers and translocates to the nucleus, where its DNA-binding domain recognizes and binds to specific DNA sequences, activating the transcription and translation of target genes [64]. Das et al. examined the expression of STAT6 in clinical PCa tissue specimens and found its presence in prostate intraepithelial neoplasia, malignant epithelial layers, and particularly high expression in the fibromuscular stroma. Furthermore, the expression level of STAT6 was significantly associated with tumor size and high histological grades of PCa [65]. TREX1 is a 3'-5' DNA exonuclease associated with the ER. Following the rupture of the micronuclear envelope, TREX1 translocates into the nucleus and degrades micronuclear DNA. This process leads to a decrease in cytoplasmic DNA levels and the inhibition of activation of the cGAS-STING pathway [66]. Moreover, TREX1 is capable of degrading tumor-derived DNA in the cytoplasm, thereby compromising the IFN-dependent antitumor immunity induced by the cGAS-STING pathway. Consequently, inhibition of TREX1 exonuclease activity is considered a potential immunotherapeutic strategy to restore antitumor immunity by reducing the degradation of tumor-derived DNA [67]. Interestingly, cancer cells exhibit cytoplasmic accumulation of DNA, which is associated with the limited activity of TREX1 induced by various mechanisms in cancers. This accumulation leads to chronic inflammation and the suppression of an effective immune response [68]. Additionally, the downregulation of TREX1 expression in senescent cells may contribute to the aberrant activation of the cGAS-STING pathway and the accumulation of cytoplasmic DNA observed in cancers [69].
Considering the strong association between aging and the cGAS-STING pathway in PCa, we identified key genes involved in this pathway and computed a gene prognostic index comprising AURKB, STAT6, and TREX1. Using the NMF algorithm as described earlier, we stratified PCa patients into two subtypes. Subtype 1 exhibited a significantly higher risk of BCR compared to subtype 2. Additionally, we developed a risk score to quantify the BCR risk in PCa patients, and the high-risk group showed a significantly greater BCR risk than the low-risk group. Analysis of mutation genes between subtype 1 and subtype 2 revealed several noteworthy findings. MUC5B was associated with the acquisition of hormone independence and preferentially expressed in hormone-independent PCa [70]. PLEC played a crucial role as a regulator in PCa development and metastasis [71]. ITGB6 encoded the β6 integrin subunit, with significantly higher expression observed in metastatic castrate-resistant androgen receptor-negative prostate tumors compared to androgen receptor-positive tumors [72]. Approximately 10% of CRPC patients exhibited HRD, and patients with M1 stage, high Gleason score, and IDC-P pathology demonstrated higher HRD scores, which correlated with poor prognosis in PCa [73, 74]. Consistently, we observed that subtype 1, associated with worse prognosis, exhibited a significantly higher HRD score than subtype 2. Furthermore, patients in subtype 1 might display sensitivity to immune therapy, as they exhibited a significantly higher TMB than subtype 2.
The limitations of our study must be mentioned. First, we did not conduct experimental validations following our analyses, which may hamper the clinical transformation of our findings. Further in-depth researches are promptly warranted. Second, besides BCR-free survival, other outcomes such as overall survival or cancer-specific survival deserve to be carefully considered as well. Third, due to the complexity of molecular pathways in tumor biology, the interaction between cGAS-STING and others may cause undetected influence on PCa, making it difficult to elucidate the detailed impacts induced by a single mechanism.
The global population is experiencing rapid aging, presenting numerous challenges. Among these challenges is the increased morbidity of PCa, which disproportionately affects older males and leads to disability, diminished quality of life, escalated costs, and even mortality. Consequently, it becomes crucial to carefully select PCa patients and identify those who are at a high risk of BCR and disease progression. Our study results have provided sufficient evidence that the prognostic index is a reliable predictor of BCR-free survival in individual PCa patients.
Conclusions
Our study revealed that the utilization of molecular subtypes and genetic risk score that associated with the cGAS-STING pathway could assist clinicians in identifying PCa patients who harbor unfavorable tumor feature. This identification enables the selection of appropriate pharmacologic therapies that specifically target the cGAS-STING pathway, thereby maximizing the potential benefits for these patients.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available in the TCGA (https://www.cancer.gov/tcga) and GEO (https://www.ncbi.nlm.nih.gov/geo/, accession number: GSE116918 and GSE46602) repositories. All data can be accessible from corresponding authors with reasonable requests.
Abbreviations
- PCa:
-
Prostate cancer
- cGAS:
-
Cyclic GMP-AMP synthase
- cGAMP:
-
Cyclic-GMP-AMP
- STING:
-
Stimulator of interferon genes
- TCGA:
-
The cancer genome atlas
- CSRGPI:
-
CGAS-STING related gene prognostic index
- DNA:
-
Deoxyribonucleic acid
- mtDNA:
-
Mitochondrial DNA
- VDAC:
-
Voltage- dependent anion channel
- MPTP:
-
Mitochondrial permeability transition pore
- dsDNA:
-
Double-stranded DNA
- CIN:
-
Chromosomal instability
- DAMP:
-
Damage-related molecular patterns
- IFN:
-
Type-I interferons
- DDR:
-
DNA damage response
- BCR:
-
Biochemical recurrence
- NMF:
-
Nonnegative matrix factorization
- DMPss:
-
Differentially methylated probes-based stemness scores
- DNAss:
-
DNA methylation-based stemness scores
- ENHss:
-
Enhancer elements/DNA methylation-based stemness scores
- EREG-METHss:
-
Epigenetically regulated DNA methylation-based stemness scores
- EREG.EXPss:
-
Epigenetically regulated RNA expression-based stemness scores
- RNAss:
-
RNA expression-based stemness scores
- HRD:
-
Homologous recombination deficiency
- LOH:
-
Loss of heterozygosity
- NEO:
-
Neoantigen
- MATH:
-
Mutant-allele tumor heterogeneity
- TMB:
-
Tumor mutation burden
- MSI:
-
Microsatellite instability
- ICB:
-
Immune checkpoint blockade
- ROS:
-
Reactive oxygen species
- MOMP:
-
Mitochondrial outer-membrane permeabilization
- ER:
-
Endoplasmic reticulum
- ERGIC:
-
ER-Golgi intermediate compartment
- IRF3:
-
Interferon regulatory factor 3
- CPC:
-
Chromosomal passenger complex
- TGFβ:
-
Transforming growth factor β
- GSVA:
-
Gene set variation analysis
- GSEA:
-
Gene set enrichment analysis
References
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:209–49.
Ghayour-Mobarhan M, Ferns GA, Moghbeli M. Genetic and molecular determinants of prostate cancer among Iranian patients: an update. Crit Rev Clin Lab Sci. 2020;57:37–53.
Chang AY, Skirbekk VF, Tyrovolas S, et al. Measuring population ageing: an analysis of the Global Burden of Disease Study 2017. Lancet Public Health. 2019;4:e159–67.
Feng D, Li D, Wu R, Han P. Scientific advancements in drug development and trials for urothelial carcinoma: insights from the 2023 ASCO-GU cancers symposium. Aging Dis. 2023;14(6):1953–7.
López-OtÃn C, Blasco MA, Partridge L, et al. The hallmarks of aging. Cell. 2013;153:1194–217.
Moskalev AA, Shaposhnikov MV, Plyusnina EN, et al. The role of DNA damage and repair in aging through the prism of Koch-like criteria. Ageing Res Rev. 2013;12:661–84.
Burtner CR, Kennedy BK. Progeria syndromes and ageing: what is the connection? Nat Rev Mol Cell Biol. 2010;11:567–78.
Lord CJ, Ashworth A. The DNA damage response and cancer therapy. Nature. 2012;481:287–94.
Ameur A, Stewart JB, Freyer C, et al. Ultra-deep sequencing of mouse mitochondrial DNA: mutational patterns and their origins. PLoS Genet. 2011;7:e1002028.
Decout A, Katz JD, Venkatraman S, Ablasser A. The cGAS-STING pathway as a therapeutic target in inflammatory diseases. Nat Rev Immunol. 2021;21:548–69.
Kwon J, Bakhoum SF. The cytosolic DNA-Sensing cGAS-STING pathway in cancer. Cancer Discov. 2020;10:26–39.
Sun L, Wu J, Du F, et al. Cyclic GMP-AMP synthase is a cytosolic DNA sensor that activates the type I interferon pathway. Science. 2013;339:786–91.
Colotta F, Allavena P, Sica A, et al. Cancer-related inflammation, the seventh hallmark of cancer: links to genetic instability. Carcinogenesis. 2009;30:1073–81.
Coussens LM, Werb Z. Inflammation and cancer. Nature. 2002;420:860–7.
Grivennikov SI, Greten FR, Karin M. Immunity, inflammation, and cancer. Cell. 2010;140:883–99.
Nowarski R, Gagliani N, Huber S, Flavell RA. Innate immune cells in inflammation and cancer. Cancer Immunol Res. 2013;1:77–84.
Karanika S, Karantanos T, Li L, et al. DNA damage response and prostate cancer: defects, regulation and therapeutic implications. Oncogene. 2015;34:2815–22.
Crasta K, Ganem NJ, Dagher R, et al. DNA breaks and chromosome pulverization from errors in mitosis. Nature. 2012;482:53–8.
Harding SM, Benci JL, Irianto J, et al. Mitotic progression following DNA damage enables pattern recognition within micronuclei. Nature. 2017;548:466–70.
Mackenzie KJ, Carroll P, Martin CA, et al. cGAS surveillance of micronuclei links genome instability to innate immunity. Nature. 2017;548:461–5.
Nguyen B, Fong C, Luthra A, et al. Genomic characterization of metastatic patterns from prospective clinical sequencing of 25,000 patients. Cell. 2022;185:563-575.e511.
Miller ET, You S, Cadaneanu RM, et al. Chromosomal instability in untreated primary prostate cancer as an indicator of metastatic potential. BMC Cancer. 2020;20:398.
Zheng W, Feng D, Xiong X, Liao X, Wang S, Xu H, Le W, Wei Q, Yang L. The role of cGAS-STING in age-related diseases from mechanisms to therapies. Aging Dis. 2023;14(4):1145–65.
Fang R, Jiang Q, Guan Y, et al. Golgi apparatus-synthesized sulfated glycosaminoglycans mediate polymerization and activation of the cGAMP sensor STING. Immunity. 2021;54:962-975 e968.
Kato K, Nishimasu H, Oikawa D, et al. Structural insights into cGAMP degradation by Ecto-nucleotide pyrophosphatase phosphodiesterase 1. Nat Commun. 2018;9:4424.
Li T, Huang T, Du M, et al. Phosphorylation and chromatin tethering prevent cGAS activation during mitosis. Science. 2021;371:eabc5386.
Qi Z, Yan F, Chen D, et al. Identification of prognostic biomarkers and correlations with immune infiltrates among cGAS-STING in hepatocellular carcinoma. Biosci Rep. 2020;40:BSR20202603.
Senegas A, Gautheron J, Maurin AG, Courtois G. IKK-related genetic diseases: probing NF-kappaB functions in humans and other matters. Cell Mol Life Sci. 2015;72:1275–87.
Wu Z, Lin Y, Liu LM, et al. Identification of cytosolic DNA sensor cGAS-STING as immune-related risk factor in renal carcinoma following pan-cancer analysis. J Immunol Res. 2022;2022:7978042.
Feng D, Shi X, Zhang F, et al. Energy metabolism-related gene prognostic index predicts biochemical recurrence for patients with prostate cancer undergoing radical prostatectomy. Front Immunol. 2022;13: 839362.
Cerami E, Gao J, Dogrusoz U, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–4.
Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14: 7.
Liberzon A, Subramanian A, Pinchback R, et al. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27:1739–40.
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 U S A. 2005;102:15545–50.
Malta TM, Sokolov A, Gentles AJ, et al. Machine learning identifies stemness features associated with oncogenic dedifferentiation. Cell. 2018;173:338-354 e315.
Lian H, Han YP, Zhang YC, et al. Integrative analysis of gene expression and DNA methylation through one-class logistic regression machine learning identifies stemness features in medulloblastoma. Mol Oncol. 2019;13:2227–45.
Bonneville R, Krook MA, Kautto EA, et al. Landscape of microsatellite instability across 39 cancer types. JCO Precis Oncol. 2017;2017:PO.17.00073.
Thorsson V, Gibbs DL, Brown SD, et al. The immune landscape of cancer. Immunity. 2018;48: 812-830 e814.
Feng D, Shi X, Zhu W, Zhang F, Li D, Han P, Wei Q, Yang L. A pan-cancer analysis of the oncogenic role of leucine zipper protein 2 in human cancer. Exp Hematol Oncol. 2022;11(1):55.
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.
Jiang P, Gu S, Pan D, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med. 2018;24:1550–8.
Dean M, Lou H. Genetics and genomics of prostate cancer. Asian J Androl. 2013;15:309–13.
Franceschi C, Capri M, Monti D, et al. Inflammaging and anti-inflammaging: a systemic perspective on aging and longevity emerged from studies in humans. Mech Ageing Dev. 2007;128:92–105.
Dadsena S, Zollo C, GarcÃa-Sáez AJ. Mechanisms of mitochondrial cell death. Biochem Soc Trans. 2021;49:663–74.
Taanman JW. The mitochondrial genome: structure, transcription, translation and replication. Biochim Biophys Acta. 1999;1410:103–23.
Miwa S, Kashyap S, Chini E, von Zglinicki T. Mitochondrial dysfunction in cell senescence and aging. J Clin Invest. 2022;132:e158447.
Green DR, Galluzzi L, Kroemer G. Mitochondria and the autophagy-inflammation-cell death axis in organismal aging. Science. 2011;333:1109–12.
Tait SW, Green DR. Mitochondrial regulation of cell death. Cold Spring Harb Perspect Biol. 2013;5:a008706.
McArthur K, Whitehead LW, Heddleston JM, et al. BAK/BAX macropores facilitate mitochondrial herniation and mtDNA efflux during apoptosis. Science. 2018;359:eaao6047.
Zhao B, Xu P, Rowlett CM, et al. The molecular basis of tight nuclear tethering and inactivation of cGAS. Nature. 2020;587:673–7.
Hopfner KP, Hornung V. Molecular mechanisms and cellular functions of cGAS-STING signalling. Nat Rev Mol Cell Biol. 2020;21:501–21.
Dobbs N, Burnaevskiy N, Chen D, et al. STING activation by translocation from the ER is associated with infection and autoinflammatory disease. Cell Host Microbe. 2015;18:157–68.
Ahn J, Xia T, Konno H, et al. Inflammation-driven carcinogenesis is mediated through STING. Nat Commun. 2014;5:5166.
Barbie DA, Tamayo P, Boehm JS, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009;462:108–12.
Göbel A, Dell’Endice S, Jaschke N, et al. The role of inflammation in breast and prostate cancer metastasis to bone. Int J Mol Sci. 2021;22:5078.
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:597–610.
Sasai K, Treekitkarnmongkol W, Kai K, et al. Functional significance of aurora kinases-p53 protein family interactions in cancer. Front Oncol. 2016;6:247.
Carmena M, Wheelock M, Funabiki H, Earnshaw WC. The chromosomal passenger complex (CPC): from easy rider to the godfather of mitosis. Nat Rev Mol Cell Biol. 2012;13:789–803.
Wu L, Ma CA, Zhao Y, Jain A. Aurora B interacts with NIR-p53, leading to p53 phosphorylation in its DNA-binding domain and subsequent functional suppression. J Biol Chem. 2011;286:2236–44.
Lee EC, Frolov A, Li R, et al. Targeting Aurora kinases for the treatment of prostate cancer. Cancer Res. 2006;66:4996–5002.
Chieffi P, Cozzolino L, Kisslinger A, et al. Aurora B expression directly correlates with prostate cancer malignancy and influence prostate cell proliferation. Prostate. 2006;66:326–33.
Song J, Zhou Y, Yakymovych I, et al. The ubiquitin-ligase TRAF6 and TGFβ type I receptor form a complex with Aurora kinase B contributing to mitotic progression and cytokinesis in cancer cells. EBioMedicine. 2022;82:104155.
Addepalli MK, Ray KB, Kumar B, et al. RNAi-mediated knockdown of AURKB and EGFR shows enhanced therapeutic efficacy in prostate tumor regression. Gene Ther. 2010;17:352–9.
Karpathiou G, Papoudou-Bai A, Ferrand E, et al. STAT6: a review of a signaling pathway implicated in various diseases with a special emphasis in its usefulness in pathology. Pathol Res Pract. 2021;223:153477.
Das S, Roth CP, Wasson LM, Vishwanatha JK. Signal transducer and activator of transcription-6 (STAT6) is a constitutively expressed survival factor in human prostate cancer. Prostate. 2007;67:1550–64.
Mohr L, Toufektchan E, von Morgen P, et al. ER-directed TREX1 limits cGAS activation at micronuclei. Mol Cell. 2021;81:724-738.e729.
Hemphill WO, Simpson SR, Liu M, et al. TREX1 as a novel immunotherapeutic target. Front Immunol. 2021;12:660184.
Anindya R. Cytoplasmic DNA in cancer cells: several pathways that potentially limit DNase2 and TREX1 activities. Biochim Biophys Acta Mol Cell Res. 2022;1869: 119278.
Takahashi A, Loo TM, Okada R, et al. Downregulation of cytoplasmic DNases is implicated in cytoplasmic DNA accumulation and SASP in senescent cells. Nat Commun. 2018;9:1249.
Legrier ME, de Pinieux G, Boyé K, et al. Mucinous differentiation features associated with hormonal escape in a human prostate cancer xenograft. Br J Cancer. 2004;90:720–7.
Buckup M, Rice MA, Hsu EC, et al. Plectin is a regulator of prostate cancer growth and metastasis. Oncogene. 2021;40:663–76.
Krishn SR, Garcia V, Naranjo NM, et al. Small extracellular vesicle-mediated ITGB6 siRNA delivery downregulates the αVβ6 integrin and inhibits adhesion and migration of recipient prostate cancer cells. Cancer Biol Ther. 2022;23:173–85.
Zhu S, Zhao J, Nie L, et al. Homologous recombination deficiency (HRD) score in aggressive prostatic adenocarcinoma with or without intraductal carcinoma of the prostate (IDC-P). BMC Med. 2022;20:237.
van Wilpe S, Simnica D, Slootbeek P, et al. Homologous recombination repair deficient prostate cancer represents an immunologically distinct subtype. Oncoimmunology. 2022;11: 2094133.
Acknowledgements
We appreciated the Chengdu Basebiotech Co,Ltd for their assistance in drawing and data process.
Funding
This program was supported by the National Natural Science Foundation of China (Grant Nos. 81974099, 82170785, 81974098, 82170784), National Key Research and Development Program of China (2021YFC2009303), programs from Science and Technology Department of Sichuan Province (Grant Nos. 2021YFH0172), 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 APC was funded by School-level Science Development fundation of North Sichuan Medical College for Young Scholars (CBY21-QB01). The funders had no role in the study design, data collection or analysis, preparation of the manuscript, or the decision to publish.
Author information
Authors and Affiliations
Contributions
QW and LY: conceptualization and methodology; ZHX and JY: data acquisition and analysis; YFL and ZHX: data interpretation; JY, ZHX, and WTZ: writing-original draft. QW and LY: writing-review and editing, supervision and funding acquisition.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Yang, J., Xu, Z., Zheng, W. et al. Identification of the cytoplasmic DNA-Sensing cGAS-STING pathway-mediated gene signatures and molecular subtypes in prostate cancer. BMC Cancer 24, 732 (2024). https://doi.org/10.1186/s12885-024-12492-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12885-024-12492-3