Weichen Wang
Zhichao Lu
Maoyu Wang
Zongheng Liu
Bing Wu
Chengkai Yang
Peipei Gong- Department of Neurosurgery, Affiliated Hospital of Nantong University, Medical School of Nantong University, Nantong, China
Background: Copper ions are essential for cellular physiology. Cuproptosis is a novel method of copper-dependent cell death, and the cuproptosis-based signature for glioma remains less studied.
Methods: Several glioma datasets with clinicopathological information were collected from TCGA, GEO and CGGA. Robust Multichip Average (RMA) algorithm was used for background correction and normalization, cuproptosis-related genes (CRGs) were then collected. The TCGA-glioma cohort was clustered using ConsensusClusterPlus. Univariate Cox regression analysis and the Random Survival Forest model were performed on the differentially expressed genes to identify prognostic genes. The cuproptosis-signature was constructed by calculating CuproptosisScore using Multivariate Cox regression analysis. Differences in terms of genomic mutation, tumor microenvironment, and enrichment pathways were evaluated between high- or low-CuproptosisScore. Furthermore, drug response prediction was carried out utilizing pRRophetic.
Results: Two subclusters based on CRGs were identified. Patients in cluster2 had better clinical outcomes. The cuproptosis-signature was constructed based on CuproptosisScore. Patients with higher CuproptosisScore had higher WHO grades and worse prognosis, while patients with lower grades were more likely to develop IDH mutations or MGMT methylation. Univariate and Multivariate Cox regression analysis demonstrated CuproptosisScore was an independent prognostic factor. The accuracy of the signature in prognostic prediction was further confirmed in 11 external validation datasets. In groups with high-CuproptosisScore, PIK3CA, MUC16, NF1, TTN, TP53, PTEN, and EGFR showed high mutation frequency. IDH1, TP53, ATRX, CIC, and FUBP1 demonstrated high mutation frequency in low-CuproptosisScore group. The level of immune infiltration increased as CuproptosisScore increased. SubMap analysis revealed patients with high-CuproptosisScore may respond to anti-PD-1 therapy. The IC50 values of Bexarotene, Bicalutamide, Bortezomib, and Cytarabine were lower in the high-CuproptosisScore group than those in the low-CuproptosisScore group. Finally, the importance of IGFBP2 in TCGA-glioma cohort was confirmed.
Conclusion: The current study revealed the novel cuproptosis-based signature might help predict the prognosis, biological features, and appropriate treatment for patients with glioma.
Introduction
Glioma is the most common primary brain tumor, accounting for about 40% of all brain tumors (1), among which glioblastoma (GBM) is the most malignant brain tumor. According to the classification of World Health Organization (WHO), gliomas are classified into four different histopathological grades: Grade I, II, III and IV, of which WHO II and III are considered low-grade gliomas (LGG). Glioblastoma (GBM, WHO IV) characterized by new angiogenesis is the most aggressive molecular subtype of glioma (2, 3). The median survival of LGG can be achieved from five to ten years through the administration of surgery, radiotherapy, and chemotherapy combination treatments, whereas the median survival of GBM is normally around one or two years (4, 5). The prognosis of glioma patients is divergent, which may be related to different tumor grades, mutation of isocitrate dehydrogenase (IDH) (6), amplification of epidermal growth factor receptor (EGFR) (7) and other factors. The current glioma prognostic evaluation model is mainly based on clinical factors, which has limited predictive ability (8–11). Therefore, a better prognostic evaluation system is needed.
Gene-regulated cell death known as “programmed cell death” is crucial for tissue homeostasis and growth, it also takes part in several pathological processes (12). At present, various types of cell death, such as necroptosis, ferroptosis, and pyroptosis, have been found to belong to necrotic programmed cell death (12). Researchers have found that cell death is closely related to tumorigenesis and prognosis. In the process of tumor development, cell death often occurs in the intratumoral area due to metabolic stress, such as hypoxia or glucose deprivation (13). Consequently, triggering programmed cell death could be a potential strategy for novel tumor therapy. Currently, knowledge of programmed cell death in cancer is continuously updated as more types of programmed cell death are discovered and recognized. Biology has long recognized copper as a vital component in all living things, from bacteria and fungi to mammals and plants, where it is a must for survival (14, 15). In humans, it binds to enzymes that assist in blood clotting, hormone maturation, and cell processing of energy, however, excessive copper can cause cell death (14). Cuproptosis as a new type of cell death is modulated and regulated by copper in cells. Copper ion directly binds to the lipoacylated components in the tricarboxylic acid cycle, leading to abnormal aggregation of lipoacylated proteins and loss of iron-sulfur cluster proteins, which leads to protein toxic stress and ultimately leads to cell death (15). Cuproptosis in glioma, however, has not yet been studied in depth.
The establishment of a glioma prognosis prediction model based on transcriptome data combined with clinical data can improve the prediction ability to a certain extent, which has highly significant clinical significance. In this study, glioma gene expression data and clinical data were collected from open databases, combined with reported cuproptosis-related genes (CRGs), was used to establish and verify cuproptosis-related clusters and signature. Subsequently, the prognosis, immune status, and treatment response of patients was also explored based on the cuproptosis-related clusters and signature.
Materials and methods
Collection and preprocessing of data for glioma
Transcripts and clinical data of glioma samples, including survival status, IDH status, grade, gender, and age, were collected from TCGA database based on UCSC Xena platform (16, 17). A total of 656 glioma patients with corresponding data were enrolled. Meanwhile, the gene expression profile of the control (non-tumoral samples) were also obtained from Genome Tissue Expression (GTEx) project (https://www.gtexportal.org) (18, 19). In addition, 11 glioma-cohorts (CGGA311, CGGA668, GSE108474, GSE13041, GSE16011, GSE43289, GSE43378, GSE4412, GSE4412, GSE68838, and GSE83300) were collected from Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/) or Chinese Glioma Genome Atlas (CGGA, http://www.cgga.org.cn/). Robust Multichip Average (RMA) algorithm was used for background correction and normalization (20). Data in the form of fragments per kilobase million (FPKM) was transformed into transcripts per kilobase million (TPM). The list of cuproptosis-related genes (CRGs) refers to the previous literature (14). In the end, ten CRGs were included in our study: FDX1, LIAS, LIPT1, DLD, DLAT, PDHA1, POHB, MTF1, GLS, and CDKN2A.
Establishment of cuproptosis-clusters and cuproptosis-signature
Based on the collected ten CRGs, the TCGA-glioma cohort was clustered using ConsensusClusterPlus package (21). Next, principal component analysis (PCA) (22) was further carried out to assess patterns associated with cuproptosis. Limma package was used to identify the differentially expressed genes (DEGs) in cuproptosis-clusters (logFC>1, P<0.05) (23). Subsequently, Univariate Cox regression analysis was performed to identify prognostic DEGs preliminarily (24). Subsequently, more valuable prognostic genes were screened out based on Random Survival Forest model (variable importance>0.25) (25, 26). To construct a cuproptosis-signature, Multivariate Cox regression analysis was used to estimate and weight the regression coefficients of the prognostic genes, and the CuproptosisScore for each glioma sample was calculated. According to the best optimal cutoff, the patients were divided into high- or low-CuproptosisScore subgroups. The association between overall survival (OS) and CuproptosisScore was analyzed using Kaplan-Meier curves. ROC curves were further utilized to validate the efficiency and accuracy of CuproptosisScore in predicting outcomes at one-, two-, and three-year. In addition, Univariate or Multivariate Cox regression analyses of CuproptosisScore and several clinical factors were performed to verify the independence of CuproptosisScore in predicting prognosis.
Genomic mutation analysis for cuproptosis-signature
The data of somatic mutations (16, 27) or copy number variation (CNV) (16) was acquired from TCGA. Genomic Identification of Significant Targets in Cancer (GISTIC) (28) algorithm was used to assess genomic characterization and CNV landscape.
Analysis of immune infiltration
Immune cell abundance (immune score), stromal cell infiltrating level (stromal score), and tumor purity (ESTIMATE score) were estimated via ESTIMATE (The Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression) algorithm (29). Using Tumor Immune Estimation Resource 2.0 (TIMER 2.0, http://timer.cistrome.org/) (30), a comprehensive analysis of immune infiltration in the tumor microenvironment of glioma was carried out. MCPcounter algorithm was used to estimate the relative proportions of ten immune cells in glioma (31). The infiltration of 28 immune cells was indicated by enrichment scores, which were calculated by single sample gene set enrichment analysis (ssGSEA) using Gene Set Variation Analysis (GSVA) R package (32, 33). Immunomodulators associated with seven different immune processes (Antigen presentation, Cell adhesion, Co-inhibitor, Co-stimulator, Ligand, Receptor and Other) were obtained from the previous literature (34). The response of glioma to anti-PD1 and anti-CTLA4 therapy was evaluated by Submap algorithm (35–37).
Enrichment pathway analysis
All gene sets from Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were downloaded from the MSigDB database (38). Gene Set Enrichment Analysis (GSEA) (39) and Gene Set Variation Analysis (GSVA) (33) were carried out according to clusterProfiler and GSVA packages, respectively.
Drug response prediction
Pharmacogenomic data from Genomics of Drug Sensitivity in Cancer (GDSC) (40) database was used to predict drug sensitivity in the enrolled glioma cases. The half maximal inhibitory concentration (IC50) value was calculated by pRRophetic package to reflect the drug response (41).
Immunohistochemistry (IHC) staining
The tissue sections through deparaffinization and dehydration were incubated with polyclonal rabbit anti-human IGFBP2 antibodies (1:50, Proteintech, 11065-3-AP) overnight at 4°C after epitope retrieval, H2O2 treatment, and non-specific antigens blocking. Next, sections were incubated with secondary antibodies (1:1000, Proteintech, SA00001-2) for two hours at room temperature, and then the signal was detected by an enhanced DAB staining kit (Proteintech, China).
Western blot
Tumor tissues as well as normal tissues were lysed in RIPA buffer (Solarbio, Beijing, China), protease and phosphatase inhibitors were added, and then denatured at 100°C for 15 min. The protein samples were then separated by 10% SDS-PAGE and transferred to polyvinylidene fluoride (PVDF) membranes. Next, PVDF membranes were blocked with 5% skim milk powder solution for 1 hour, and incubated with primary antibodies, including anti-IGFBP2 antibody (1:1000, Proteintech, 11065-3-AP), anti-GAPDH antibody (1:5000, Abcam, ab9485) overnight, followed by secondary antibodies (1:2000, Proteintech, SA00001-2) for 2 hours at room temperature, observed with the ECL kit chemiluminescence reagent (Billerica Millipore, MA, USA). Protein band signals were detected by the Chemidoc detection system (Bio-Rad, Hercules, CA, USA) and quantified by ImageJ software (National Institutes of Health, USA).
Statistical analysis
The Wilcoxon test was used to compare non-normally distributed data. The T-test was used to compare normally distributed variables between two groups. The R package survminer was used to estimate OS between two groups using Kaplan-Meier survival plots. Cox regression of survival analysis was also performed by survival. Time-dependent receiver operating characteristics (ROC) curves were plotted using R package timeROC. All heatmaps were performed through pheatmap package. The data were visualized by ggplot2 (V4.1.2). P<0.05 was considered statistically significant.
Results
Characteristics of cuproptosis-clusters for TCGA-glioma
The clinical information of patients from TCGA was listed in the Supplemental Table S1. The correlations among the ten CRGs were mostly positive, and the most strongly associated variables were DLD and DLAT (Figure 1A). Based on ConsensusClusterPlus, the optimal number of clusters was determined, k=2 (Figure 1B). Furthermore, PCA analysis was further used to validate that patients in the two subclusters clustered separately, which confirmed the reliability of the clustering results (Figure 1C). Patients with glioma in cluster2 had significantly better clinical outcomes than those in cluster1 (Figure 1D). Compared with the other subcluster, the expressions of FDX1, DLD, DLAT, PDHB, GLS were relatively high in cluster1, while the expression of CDKN2A was relatively high in cluster2, which indicates that these CRGs may be genetic markers for identifying different clusters (Figure 1E). Interestingly, we discovered that cluster2 had a greater percentage of IDH mutation status, which may be one of the factors contributing to a better prognosis of this subcluster (Figure 1E).
Figure 1 Characteristics of CuproptosisCluster in TCGA-glioma cohort. (A) The correlations among the ten cuproptosis-related genes (CRGs). The color represents the correlation coefficient. (B) Cluster diagram for subtype analysis of glioma samples. The intragroup correlations were the highest and the inter-group correlations were low when k=2. (C) PCA analysis for the two subclusters. (D) Kaplan-Meier survival curve showing survival probability of cluster1 and cluster2. (E) Heatmap showing the expression levels of the ten cuproptosis-related genes (CRGs) in different clinical features and clusters. *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001.
Establishment of cuproptosis-signature
In the two subclusters, a total of 27 differentially expressed genes (DEGs) were identified (logFC>1, P<0.05), and the volcano map accurately reflected the gene expression differences between the two subclusters (Figure 2A). After Univariate Cox regression analysis, 16 potential pro-oncogenes (HR>1; CYTOR, EMP3, OCIAD2, PLA2G5, FABP5P7, IGFBP2, TSTD1, TIMP1, RBP1, METTL7B, POSTN, CHI3L1, H19, CXCL14, LTF and ENC1) and 9 potential suppressor genes were identified (HR<1; CAMK2A, LINC01088, CDKN2B, LINC00689, TPTEP1, C5orf38, KLRC2, VIPR2, and SMOC1) (Figure 2B). Figure 2C showed the distribution of error rates in Random Survival Forest model, after which the relative importance of seven genes (H19, CYTOR, IGFBP2, EMP3, KLRC2, C5orf38, and CHI3L1) was established (variable importance>0.25, Figure 2D). Multivariate Cox regression analysis was used to develop the cuproptosis-signature, and the CuproptosisScore for each glioma sample was calculated according to the following formula: 0.0621*ExprH19+0.0196*ExprCYTOR+0.2739*ExprIGFBP2+0.0183*ExprKLRC2+0.0036*ExprC5orf38+0.1406*ExprCHI3L1. Heatmap displayed the distribution of six genes in cuproptosis-signature, CuproptosisScore and the clinical characteristics (Figure 2E). It was clear that a higher CuproptosisScore was linked to higher expressions of H19, CHI3L1, CYTOR, and IGFBP2 and, in contrast, was associated with lower expressions of KLRC2 and C5orf38 (Figure 2E). In the meantime, IDH mutation status was more likely to be present in glioma patients with lower CuproptosisScore. (Figure 2E).
Figure 2 Establishment of cuproptosis-signature. (A) The volcano map reflects the differentially expressed genes identified (logFC > 1, P < 0.05). (B) The forest figure for Univariate Cox regression analysis of the differentially expressed genes. (C) The distribution of error rates in Random Survival Forest model. (D) The variable relative importance of the seven genes. (E) Heat map showing the relationship between six genes in the cuproptosis-signature and CuproptosisScore distribution and its clinical characteristics. ****P < 0.0001.
Prognostic potential of cuproptosis-signature
Next, we analyzed the CuproptosisScore of TCGA patients among WHO grades, mutation status and MGMT methylation status (Figure 3A). Patients with higher CuproptosisScore had higher WHO grades, while patients with lower grades were more likely to develop IDH mutations or MGMT methylation (Figure 3A), all of which may explain the significantly better clinical outcomes of patients with lower CuproptosisScore (P<0.001, Figure 3B). In addition, Univariate and Multivariate Cox regression analysis of CuproptosisScore and clinicopathologic features demonstrated that both CuproptosisScore and Grade were independent prognostic factors for patients with glioma (Figure 3C). The survival ROC curves predicted by the cuproptosis-signature showed that the AUCs were all greater than 0.8, indicating the effectiveness of the cuproptosis-signature in predicting prognosis for glioma at the 1-year (AUC=0.898), 2-year (AUC=0.922), 3-year (AUC=0.918), 4-year (AUC=0.867), and 5-year (AUC=0.828) time points (Figure 3D). Furthermore, we conducted Univariate Cox regression analysis on the OS (overall survival) of glioma patients based on the external validation data sets, and the results showed that HR was greater than 1 in all of the 11 data sets, which further validated the accuracy of our constructed cuproptosis-signature in prognostic prediction (Figure 3E).
Figure 3 Prognostic potential of cuproptosis-signature. (A) The violin figures for comparing the CuproptosisScore of TCGA patients among WHO grades, mutation status and MGMT methylation status. (B) Kaplan-Meier survival curve showing survival probability of high-CuproptosisScore or low-CuproptosisScore subgroups. (C) The forest figure for Univariate or Multivariate Cox regression analysis of CuproptosisScore and clinicopathologic features. (D) The 1-year, 2-year, 3-year, 4-year, and 5-year survival ROC curves are predicted by the cuproptosis-signature. (E) Univariate Cox regression analysis of the cuproptosis-signature in 11 external validation data sets. ***P < 0.001.
Genomic mutation analysis for cuproptosis-signature
GISTIC algorithm was used to assess the genomic characterization landscape between high- CuproptosisScore and low-CuproptosisScore subgroups, which was shown in Figure 4A. In patients with high-CuproptosisScore, PIK3CA, MUC16, NF1, TTN, TP53, PTEN, and EGFR had high mutation frequency (over 10%, Figure 4B), while in those with low-CuproptosisScore, IDH1, TP53, ATRX, CIC, and FUBP1 had high mutation frequency (over 10%, Figure 4C). TP53 had high mutation rates in both groups (26% and 51%, respectively). In agreement with the findings above, the mutation rate of IDH1 was particularly high in the low-CuproptosisScore group, reaching as high as 89% (Figure 4C).
Figure 4 Genomic mutation analysis for cuproptosis-signature. (A) Genomic characterization landscape of high-CuproptosisScore or low-CuproptosisScore subgroups. (B) Gene mutation frequency in high-CuproptosisScore. (C) Gene mutation frequency in low-CuproptosisScore.
Immune status for cuproptosis-signature
Based on ESTIMATE, MCPcounter, ssGSEA, and TIMER algorithms mentioned in the Methods section, the heatmap showed the abundance of infiltrating immune cell populations at different CuproptosisScores (Figure 5A). In general, the level of immune infiltration increased as the CuproptosisScore increased (Figure 5A). However, it was observed that patients with lower CuproptosisScores had more tumor purity (Figure 5A). In addition, our results showed that glioma patients with high CuproptosisScores also had higher levels of TMB (Figure 5B), GEP (Figure 5C), and CYT (Figure 5D). GSVA analysis also suggested that patients with high CuproptosisScores were enriched in immune-related pathways, such as negative regulation of macrophage apoptotic process, macrophage fusion, B cell receptor signaling pathway, T cell receptor signaling pathway, and primary immunodeficiency (Figure 5E).
Figure 5 Immune status for cuproptosis-signature. (A) The heatmap shows the abundance of infiltrating immune cell populations at different CuproptosisScores. (B–D) Glioma patients with high CuproptosisScores had higher levels of TMB (B), GEP (C), and CYT (D). (E) The heatmap shows CuproptosisScores, clinical features, and immune-related pathways based on GSVA analysis. ****P < 0.0001.
Immunotherapy and chemotherapy of cuproptosis-signature
Immunomodulators (IMs) are closely related to the immunotherapy of malignant tumors, as well, agonists and antagonists for immunomodulators are also being studied (42). The expression of IM-related genes varied across high-CuproptosisScore or low-CuproptosisScore subgroups (Figure 6A). As the heatmap demonstrated that the level of TNFSF9, IL13, and TIGIT, showed no difference between the two groups, VTCN1, TNF, CX3CL1, IL12A, HMGB1, EDNRB, and TLR4 were highly expressed in the low-CuproptosisScore group, and the remaining genes were highly expressed in the high-CuproptosisScore group (Figure 6A). In addition, SubMap analysis revealed patients with high-CuproptosisScore may respond to anti-PD-1 therapy (Figure 6B). This may be due to the high expression of IMs in this group of patients. We also investigated the IC50 values of four chemotherapeutics (Bexarotene, Bicalutamide, Bortezomib, and Cytarabine) between the high- and low-CuproptosisScore groups. Results showed that IC50 values of patients in the high- CuproptosisScore group were lower than those in the low-CuproptosisScore group, suggesting that patients in the high-CuproptosisScore group were more likely to benefit from these four drugs (Figure 6C).
Figure 6 Immunotherapy and chemotherapy of cuproptosis-signature. (A) Correlation of CuproptosisScore with seven immunomodulators in gliomas. (B) SubMap analysis for cuproptosis-signature in gliomas. (C) Box plots of estimated IC50 for several chemotherapeutic agents in the high- or low-CuproptosisScore groups. *P < 0.05; **P < 0.01; ***P < 0.001, ****P < 0.0001.
The importance of IGFBP2 in TCGA-glioma cohort
The aforementioned results showed that the cuproptosis-signature we created has substantial clinical significance. Next, we randomly selected one gene from this signature, IGFBP2, and explored its important value in gliomas. As is evident from WB and IHC results that IGFBP2 was significantly higher in the six tumor tissues than in the paired adjacent tumor tissues (Figures 7A–C). The patients with IGFBP2 expression values were listed in the Supplemental Table S2. The expression level of IGFBP2 was further compared between the glioma sample and the healthy control sample, it was found that IGFBP2 was significantly overexpressed in the cancer tissue (Figure 7D). In terms of the survival curve, glioma patients with low expression of IGFBP2 had better survival, indicating that IGFBP2 may be a promoter of the malignant progression of glioma. The 1-year (AUC=0.877), 2-year (AUC=0.92), 3-year (AUC=0.91), 4-year (AUC=0.858), and 5-year (AUC=0.822) survival ROC curves predicted by IGFBP2 revealed that the AUCs were all higher than 0.8, indicating the efficiency of IGFBP2 in predicting prognosis for glioma (Figures 7E, F). Therefore, IGFBP2 is highly likely to be the oncogenic gene of glioma. The heatmap (Figure 7G) showed that the expression value of IGFBP2 was positively correlated with the expression value of eight immune checkpoints (LAG3, CD274, PDCD1LG2, TNFRSF9, PDCD1, CTLA4, CD247, and TNFRSF4). Moreover, we carried out GSEA analysis to explore cancer and immune-related signaling pathways positively modulated by IGFBP2. We found six signaling pathways (Figure 7H): immune response, T cell receptor signaling pathway, regulation of immune response, pathways in cancer, p53 signaling pathway, and TF signaling pathway.
Figure 7 (A) WB for IGFBP2 in 3 pairs patients from Nantong cohort. (B) Represented IHC for IGFBP2 in three parents with different WHO stage from Nantong cohort. (C) Boxplot of IHC for IGFBP2 in six pairs parents from Nantong cohort. (D) The expression level of IGFBP2 in glioma sample and the control normal sample. (E) Kaplan-Meier survival curve showing survival probability of high- or low-expression IGFBP2. (F) The 1-year, 2-year, 3-year, 4-year, and 5-year survival ROC curves are predicted by the expression of IGFBP2. (G) The heat map shows the correlation between IGFBP2 and eight immune checkpoints in TCGA. (H) GSEA maps of cancer and immune-related signaling pathways positively modulated by IGFBP2. ∗∗P < 0.01; ∗∗∗P < 0.001; ∗∗∗∗P < 0.0001.
Discussion
Gliomas, especially glioblastoma (GBM), are the most destructive brain tumors within the human nervous system (43). Despite improvements in glioma diagnosis and treatment in recent years, gliomas are still difficult to treat with surgery alone due to their invasive and quickly proliferating nature. Patients with postoperative recurrence have a poor prognosis, with the median survival time only being extended by a few months (44, 45). One of the crucial characteristics of tumor cells is their resistance to cell death. Unrestricted proliferation is typical for tumor cells, and they overcome growth inhibition by resisting death and avoiding being killed by immune cells. However, due to metabolic stress, such as hypoxia and glucose deprivation, necrotizing cell death often occurs in the interior of solid tumors, which affects the occurrence and development of tumors by reshaping the tumor microenvironment. With the discovery of ever-more programmed death modes and the elucidation of associated molecular mechanisms, our understanding of the role of cell death in tumor is constantly updated. Since multiple forms of cell death occur simultaneously in tumors, an in-depth study of cell death on the occurrence and development of tumors can help us better understand their pathogenesis and pave the way for the creation of effective anti-tumor medications. For example, abnormalities of apoptosis pathways play critical roles in tumorigenesis, and tumor cell avoidance of apoptosis has long been thought to lead to primary or acquired therapeutic resistance (46). Necroptosis has both pro-tumor and anti-tumor effects in different types of cancer (13). Inducing necroptosis of tumor cells is an important way to overcome chemotherapy resistance of tumor cells. Finding a novel way to precisely regulate necroptosis might be an essential research target in the field of tumor therapy in the future (47). Numerous pieces of evidence suggest that pyrotopia plays an important role in tumor progression, and inducing pyrotopia has become one of the focuses of cancer immunotherapy (48–51). Ferroptosis is a type of cell death induced by oxidative stress. Cancer cells metabolize more efficiently than normal cells, with a higher ROS load (52) and require large amounts of iron, thus they are more sensitive to ferroptosis than normal cells (53). However, cancer cells also employ additional genetic or epigenetic mechanisms to combat elevated ROS levels, thereby reducing their sensitivity to ferroptosis (54). Therefore, ferroptosis is closely related to the occurrence and development of tumors.
Copper ions can be combined with a variety of proteins or enzymes, as cofactors or structural components, involved in the regulation of energy metabolism, mitochondrial respiration, antioxidant, and other physiological processes (55, 56). The content of copper ions maintains a dynamic balance, which can lead to oxidative stress (55) and abnormal autophagy (56), and thus induce a variety of copper or copper ion-related diseases. Tsvetkov et al. proposed for the first time that a new method of cell death with copper dependence, which was called cuprotosis (14). Several studies had shown that copper metabolism was associated with tumorigenesis, and cancer cells have a higher demand for copper than normal cells (57–60). Wang et al. found that blocking Cu2+ transport can cause oxidative stress and decrease cellular ATP levels, which in turn activates AMP-activated protein kinase (AMPK), leading to reduced adipogenesis and inhibiting tumor cell proliferation (61). Studies have confirmed that copper is closely related to the expression level of hypoxia-inducible factor 1α (HIF-1α) (62). The use of copper chelating agent tetrathiomolybdate can significantly reduce the content of Cu2+ in vivo, and dramatically reduce tumor angiogenesis, restrain tumor growth, and reduce the invasion of breast cancer cells (63). In conclusion, cuprotosis is a novel kind of cellular regulatory death that impacts copper metabolism. The identification of cuprotosis molecular pathways has implications for the mechanism of cuprotosis, cancer drug discovery, and a deeper understanding of copper metabolic diseases. In this study, we included ten cuproptosis-related genes (CRGs): FDX1, LIAS, LIPT1, DLD, DLAT, PDHA1, POHB, MTF1, GLS, and CDKN2A. The correlations among the ten CRGs were primarily positive in TCGA-glioma cohort. In addition, a prognostic signature based on CuproptosisScore was established to explore its prognostic and clinical value in glioma. The current research provides a reference for exploring the mechanism of cuprotosis in the development of glioma.
As for the ten CRGs, the phosphorylation and dephosphorylation of PDHA1 (Pyruvate Dehydrogenase E1 alpha Subunit) are key modulators of deactivation and activation of PDC(Pyruvate Dehydrogenase Complex) (64). It was reported that the increasing level of PDHA1 was observed in the higher grade of glioma and PDHA1 could regulate the migration of glioma cells (65). LINC00665 promoted MTF1 degradation, and MTF1 bound to the promoter region of GTSE1 and transcription promoted GTSE1 expression, which proved that LINC00665/MTF1/GTSE1 axis played an important role in regulating the biological behavior of glioma cells (66). GLS are oncogenic genes of glioma (67, 68) and Qiangzhen Huang et al. found that GLS could regulate the effect of SNAP25 in glioma (66, 69). CDKN2A homozygous deletion was reported to serve as an adverse prognostic factor for IDH-mutant gliomas (70–72). However, to our knowledge, the role of the remaining CRGs in glioma has not been reported in the literature.
The complexity of gliomas is mainly reflected by their molecular heterogeneity. Molecular subtypes can well predict the occurrence and development of glioma polymorphism, which can assist us in developing better treatments (73). Mesenchymal subtypes are particularly malignant compared to other subtypes (neurogenic, canonical, and preneurotic) according to TCGA classification, with recurrent GBM always fatal and often presenting as a mesenchymal phenotype (74–76). In addition, mesenchymal subtypes of gliomas expressed higher levels of angiogenic markers in addition to higher levels of necrosis (74, 77). It has been reported that the transition from the former neural subtype to the mesenchymal subtype is closely associated with treatment resistance and poor prognosis (78). Currently, no fully verified and feasible classification system has been applied to clinical practice, and the glioma classification system needs to be continuously explored and improved. In this study, we determined the optimal number (k=2) of clusters based on R package (Figure 1B). Furthermore, the reliability of the clustering results was confirmed by PCA analysis (Figure 1C). Patients with glioma in cluster2 had significantly better clinical outcomes than those in cluster1 (Figure 1D). Interestingly, it was observed that cluster2 had a higher proportion of IDH mutation status, which may be one of the reasons for the better prognosis of this subcluster (Figure 1E). In the two subclusters, a total of 27 differentially expressed DEGs were identified (Figure 2A). When combined with the aforementioned findings, our study sheds light on the need for a new glioma typing system.
Although the WHO classification system has been used for many years to predict the prognosis of patients with glioma, it is occasionally inaccurate due to the heterogeneity of the tumor. In addition to identifying potential biomarkers, new advances in bioinformatics and genome sequencing can help predict cancer patient outcomes and treatment strategies (79, 80). Studies have shown that the prognostic value of a single biomarker is limited, and it is better to integrate multiple biomarkers into a single model (81). For example, three IncRNAs can predict the prognosis of colorectal cancer (CRC) based on a network of metastasis-related competing endogenous RNAs (ceRNA) (82). By extracting TCGA-related data, four IncRNA signals can effectively predict the survival time of lung adenocarcinoma (LUAD) (83). Importantly, recent studies have confirmed the predictive power of some IncRNA prognostic signatures in gliomas, such as immune-associated IncRNAs and autophagy-associated IncRNAs, which have strong prognostic potential for glioma patients (44, 45). In this study, we also used bioinformatics methods to identify a signature containing multiple genes based on CuproptosisScore. Patients with higher CuproptosisScore had higher WHO grades, while patients with lower grades were more likely to develop IDH mutations or MGMT methylation (Figure 3A) and patients with lower CuproptosisScore had the significantly better clinical outcomes (Figure 3B). In addition, Univariate and Multivariate Cox regression analysis of the signature demonstrated that the cuproptosis-signature was an independent prognostic factor for patients with glioma (Figure 3C). The survival ROC curves indicated the efficiency of the cuproptosis-signature in predicting prognosis for glioma (Figure 3D). Furthermore, we conducted Univariate Cox regression analysis on the OS of glioma patients based on the external validation data sets, and the results demonstrated the accuracy of our constructed cuproptosis-signature in prognostic prediction. In conclusion, the signature we have identified has excellent prognostic power (Figure 3E).
We randomly selected one gene from this signature, IGFBP2, and explored its important value in gliomas. In glioma, IGFBP2 is often involved in the activation of PTEN, AKT and other related pathways, leading to enhanced invasiveness and malignancy (84, 85). Studies have shown that overexpression of IGFBP2 can increase the malignant degree of glioma and up-regulate the expression of invasion protein MMP2, thereby enhancing the invasion ability of glioma cells (86). Previous studies have confirmed that the expression levels of IGFBP2 transcripts and proteins are positively correlated with the malignant degree of glioma, suggesting that IGFBP2 plays an important role in malignant transformation, tumor necrosis and metastasis of glioma (87, 88). In our study, we also found that IGFBP2 is highly likely to be the oncogenic gene of glioma.
However, our study still has several shortcomings. Due to the lack of data of clinical samples collected by us, the prognostic cuproptosis-signature constructed in this study was based on the public database, which may have a certain bias in the source of samples. The hypothesis obtained in this study was not verified by experimental results, and the next step is to be confirmed by various in vivo and in vitro experiments and larger multicenter studies. In conclusion, this prognostic cuproptosis-signature still needs to be further tested, evaluated and applied in a wide range of clinical settings.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
WW, ZCL, and MW carried out experiments and analysis. WW, PG, ZCL, MW, ZHL, BW, HH, and CY wrote the manuscript. WW and PG conceived the study. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Medical Scientific Research Project of Jiangsu Provincial Health Commission (H2019058).
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.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.998236/full#supplementary-material
References
1. Gusyatiner O, Hegi ME. Glioma epigenetics: From subclassification to novel treatment options. Semin Cancer Biol (2018) 51:50–8. doi: 10.1016/j.semcancer.2017.11.010
2. Samioti SE, Benos LT, Sarris IE. Effect of fractal-shaped outer boundary of glioblastoma multiforme on drug delivery. Comput Methods Programs BioMed (2019) 178:191–9. doi: 10.1016/j.cmpb.2019.06.031
3. Polivka J Jr., Polivka J, Holubec L, Kubikova T, Priban V, Hes O, et al. Advances in experimental targeted therapy and immunotherapy for patients with glioblastoma multiforme. Anticancer Res (2017) 37:21–33. doi: 10.21873/anticanres.11285
4. Lapointe S, Perry A, Butowski NA. Primary brain tumours in adults. Lancet (2018) 10145:432–46. doi: 10.1016/s0140-6736(18)30990-5
5. Garcia CR, Slone SA, Dolecek TA, Huang B, Neltner JH, Villano JL. Primary central nervous system tumor treatment and survival in the united states, 2004-2015. J Neurooncol (2019) 144:179–91. doi: 10.1007/s11060-019-03218-8
6. Eckel-Passow JE, Lachance DH, Molinaro AM, Walsh KM, Decker PA, Sicotte H, et al. Glioma groups based on 1p/19q, IDH, and TERT promoter mutations in tumors. N Engl J Med (2015) 372:2499–508. doi: 10.1056/NEJMoa1407279
7. Eskilsson E, Røsland GV, Solecki G, Wang Q, Harter PN, Graziani G, et al. EGFR heterogeneity and implications for therapeutic intervention in glioblastoma. Neuro Oncol (2018) 20:743–52. doi: 10.1093/neuonc/nox191
8. Zhao YY, Chen SH, Hao Z, Zhu HX, Xing ZL, Li MH. A nomogram for predicting individual prognosis of patients with low-grade glioma. World Neurosurg (2019) 130:e605–12. doi: 10.1016/j.wneu.2019.06.169
9. Kiran M, Chatrath A, Tang X, Keenan DM, Dutta A. & dutta, a. a prognostic signature for lower grade gliomas based on expression of long non-coding RNAs. Mol Neurobiol (2019) 56:4786–98. doi: 10.1007/s12035-018-1416-y
10. Razzak M. Genetics: new molecular classification of gastric adenocarcinoma proposed by the cancer genome atlas. Nat Rev Clin Oncol (2014) 11:499. doi: 10.1038/nrclinonc.2014.138
11. Deng S, Li Y, Yi G, Lei B, Guo M, Xiang W, et al. Overexpression of COX7A2 is associated with a good prognosis in patients with glioma. J Neurooncol (2018) 136:41–50. doi: 10.1007/s11060-017-2637-z
12. Vanden Berghe T, Linkermann A, Jouan-Lanhouet S, Walczak H, Vandenabeele P. Regulated necrosis: the expanding network of non-apoptotic cell death pathways. Nat Rev Mol Cell Biol (2014) 15:135–47. doi: 10.1038/nrm3737
13. Qin X, Ma D, Tan YX, Wang HY, Cai Z. The role of necroptosis in cancer: A double-edged sword? Biochim Biophys Acta Rev Cancer (2019) 1871:259–66. doi: 10.1016/j.bbcan.2019.01.006
14. Tsvetkov P, Coy S, Petrova B, Dreishpoon M, Verma A, Abdusamad M, et al. Copper induces cell death by targeting lipoylated TCA cycle proteins. Science (2022) 375:1254–61. doi: 10.1126/science.abf0529
15. Kim BE, Nevitt T, Thiele DJ. Mechanisms for copper acquisition, distribution and regulation. Nat Chem Biol (2008) 4:176–85. doi: 10.1038/nchembio.72
16. Wang Z, Jensen MA, Zenklusen JC.A practical guide to the cancer genome atlas (TCGA).Methods Mol Biol (2016) 1418:111–41. doi: 10.1007/978-1-4939-3578-9_6
17. Goldman MJ, Craft B, Hastie M, Repečka K, McDade F, Kamath A, et al. Visualizing and interpreting cancer genomics data via Xena Platform. Nat Biotechnol (2020) 38:675–8. doi: 10.1038/s41587-020-0546-8
18. Human genomics. the genotype-tissue expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science (2015) 348:648–60. doi: 10.1126/science.1262110
19. Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med (2015) 21:938–45. doi: 10.1038/nm.3909
20. Katz S, Irizarry RA, Lin X, Tripputi M, Porter MW. A summarization approach for affymetrix GeneChip data using a reference training set from a large, biologically diverse database. BMC Bioinf (2006) 7:464. doi: 10.1186/1471-2105-7-464
21. Seiler M, Huang CC, Szalma S, Bhanot G. ConsensusCluster: a software tool for unsupervised cluster discovery in numerical data. Omics (2010) 14:109–13. doi: 10.1089/omi.2009.0083
22. Ringnér M. What is principal component analysis? Nat Biotechnol (2008) 26:303–4. doi: 10.1038/nbt0308-303
23. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43:e47. doi: 10.1093/nar/gkv007
24. van Dijk PC, Jager KJ, Zwinderman AH, Zoccali C, Dekker FW. The analysis of survival data in nephrology: basic concepts and methods of cox regression. Kidney Int (2008) 74:705–9. doi: 10.1038/ki.2008.294
25. Taylor JM. Random survival forests. J Thorac Oncol (2011) 6:1974–5. doi: 10.1097/JTO.0b013e318233d835
26. Ishwaran H, Lu M. Standard errors and confidence intervals for variable importance in random forest regression, classification, and survival. Stat Med (2019) 38:558–82. doi: 10.1002/sim.7803
27. Ng PK, Li J, Jeong KJ, Shao S, Chen H, Tsang YH, et al. Systematic functional annotation of somatic mutations in cancer. Cancer Cell (2018) 33:450–62.e410. doi: 10.1016/j.ccell.2018.01.021
28. Mermel CH, Schumacher SE, Hill B, Meyerson ML, Beroukhim R, Getz G. GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol (2011) 12:R41. doi: 10.1186/gb-2011-12-4-r41
29. Yoshihara K, Shahmoradgoli M, Martinez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun (2013) 4:2612. doi: 10.1038/ncomms3612
30. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res (2017) 77:e108–10. doi: 10.1158/0008-5472.Can-17-0307
31. Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol (2016) 17:218. doi: 10.1186/s13059-016-1070-5
32. Xiao B, Liu L, Li A, Xiang C, Wang P, Li H, et al. Identification and verification of immune-related gene prognostic signature based on ssGSEA for osteosarcoma. Front Oncol (2020) 10:607622. doi: 10.3389/fonc.2020.607622
33. Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinf (2013) 14:7. doi: 10.1186/1471-2105-14-7
34. Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang T-H, et al. The immune landscape of cancer. Immunity (2019) 51:411–2. doi: 10.1016/j.immuni.2019.08.004
35. Hoshida Y, Brunet JP, Tamayo P, Golub TR, Mesirov JP. Subclass mapping: identifying common subtypes in independent disease data sets. PloS One (2007) 2:e1195. doi: 10.1371/journal.pone.0001195
36. McGranahan N, Furness AJ, Rosenthal R, Ramskov S, Lyngaa R, Saini SK, et al. Clonal neoantigens elicit T cell immunoreactivity and sensitivity to immune checkpoint blockade. Science (2016) 351:1463–9. doi: 10.1126/science.aaf1490
37. Lu X, Jiang L, Zhang L, Zhu Y, Hu W, Wang J, et al. Immune signature-based subtypes of cervical squamous cell carcinoma tightly associated with human papillomavirus type 16 expression, molecular features, and clinical outcome. Neoplasia (2019) 21:591–601. doi: 10.1016/j.neo.2019.04.003
38. Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst (2015) 1:417–25. doi: 10.1016/j.cels.2015.12.004
39. Subramanian A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
40. Yang W, Soares J, Greninger P, Edelman EJ, Lightfoot H, Forbes S, et al. Genomics of drug sensitivity in cancer (GDSC): A resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res (2013) 41:D955–961. doi: 10.1093/nar/gks1111
41. Geeleher P, Cox N, Huang RS. pRRophetic: an r package for prediction of clinical chemotherapeutic response from tumor gene expression levels. PloS One (2014) 9:e107468. doi: 10.1371/journal.pone.0107468
42. Tang J, Shalabi A, Hubbard-Lucey VM. Comprehensive analysis of the clinical immuno-oncology landscape. Ann Oncol (2018) 29:84–91. doi: 10.1093/annonc/mdx755
43. Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol (2015) 16:e173–180. doi: 10.1016/s1470-2045(14)71116-7
44. Luan F, Chen W, Chen M, Yan J, Chen H, Yu H, et al. An autophagy-related long non-coding RNA signature for glioma. FEBS Open Bio (2019) 9:653–67. doi: 10.1002/2211-5463.12601
45. Xia P, Li Q, Wu G, Huang Y. An immune-related lncRNA signature to predict survival in glioma patients. Cell Mol Neurobiol (2021) 41:365–75. doi: 10.1007/s10571-020-00857-8
46. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell (2011) 144:646–74. doi: 10.1016/j.cell.2011.02.013
47. Chefetz I, Grimley E, Yang K, Hong L, Vinogradova EV, Suciu R, et al. A pan-ALDH1A inhibitor induces necroptosis in ovarian cancer stem-like cells. Cell Rep (2019) 26:3061–75.e3066. doi: 10.1016/j.celrep.2019.02.032
48. Wang Q, Wang Y, Ding J, Wang C, Zhou X, Gao W, et al. A bioorthogonal system reveals antitumour immune function of pyroptosis. Nature (2020) 579:421–6. doi: 10.1038/s41586-020-2079-1
49. Kaplanov I, Carmi Y, Kornetsky R, Shemesh A, Shurin GV, Shurin MR, et al. Blocking IL-1β reverses the immunosuppression in mouse breast cancer and synergizes with anti-PD-1 for tumor abrogation. Proc Natl Acad Sci U.S.A. (2019) 116:1361–9. doi: 10.1073/pnas.1812266115
50. Kaplanski G. Interleukin-18: Biological properties and role in disease pathogenesis. Immunol Rev (2018) 281:138–53. doi: 10.1111/imr.12616
51. Wang Y, Gao W, Shi X, Ding J, Liu W, He H, et al. Chemotherapy drugs induce pyroptosis through caspase-3 cleavage of a gasdermin. Nature (2017) 547:99–103. doi: 10.1038/nature22393
52. Tsoi J, Robert L, Paraiso K, Galvan C, Sheu KM, Lay J, et al. Multi-stage differentiation defines melanoma subtypes with differential vulnerability to drug-induced iron-dependent oxidative stress. Cancer Cell (2018) 33:890–904.e895. doi: 10.1016/j.ccell.2018.03.017
53. Carlisle AE, Lee N, Matthew-Onabanjo AN, Spears ME, Park SJ, Youkana D, et al. Selenium detoxification is required for cancer-cell survival. Nat Metab (2020) 2:603–11. doi: 10.1038/s42255-020-0224-7
54. Rojo de la Vega M, Chapman E, Zhang DD. NRF2 and the hallmarks of cancer. Cancer Cell (2018) 34:21–43. doi: 10.1016/j.ccell.2018.03.022
55. Guo H, Wang Y, Cui H, Ouyang Y, Yang T, Liu C, et al. Copper induces spleen damage through modulation of oxidative stress, apoptosis, DNA damage, and inflammation. Biol Trace Elem Res (2022) 200:669–77. doi: 10.1007/s12011-021-02672-8
56. Jian Z, Guo H, Liu H, Cui H, Fang J, Zuo Z, et al. Oxidative stress, apoptosis and inflammatory responses involved in copper-induced pulmonary toxicity in mice. Aging (Albany NY) (2020) 12:16867–86. doi: 10.18632/aging.103585
57. Atakul T, Altinkaya SO, Abas BI, Yenisey C. Serum copper and zinc levels in patients with endometrial cancer. Biol Trace Elem Res (2020) 195:46–54. doi: 10.1007/s12011-019-01844-x
58. Feng Y, Zeng JW, Ma Q, Zhang S, Tang J, Feng JF, et al. Serum copper and zinc levels in breast cancer: A meta-analysis. J Trace Elem Med Biol (2020) 62:126629. doi: 10.1016/j.jtemb.2020.126629
59. Pavithra V, Sathisha TG, Kasturi K, Mallika DS, Amos SJ, Ragunatha S, et al. Serum levels of metal ions in female patients with breast cancer. J Clin Diagn Res (2015) 9:BC25–c27. doi: 10.7860/jcdr/2015/11627.5476
60. Zhang M, Shi M, Zhao Y. Association between serum copper levels and cervical cancer risk: A meta-analysis. Biosci Rep (2018) 38:4. doi: 10.1042/bsr20180161
61. Wang J, Luo C, Shan C, You Q, Lu J, Elf S, et al. Inhibition of human copper trafficking by a small molecule significantly attenuates cancer cell proliferation. Nat Chem (2015) 7:968–79. doi: 10.1038/nchem.2381
62. Zimna A, Kurpisz M. Hypoxia-inducible factor-1 in physiological and pathophysiological angiogenesis: Applications and therapies. BioMed Res Int (2015) 2015:549412. doi: 10.1155/2015/549412
63. MacDonald G, Nalvarte I, Smirnova T, Vecchi M, Aceto N, Dolemeyer A, et al. Memo is a copper-dependent redox protein with an essential role in migration and metastasis. Sci Signal (2014) 7:ra56. doi: 10.1126/scisignal.2004870
64. Holness MJ, Sugden MC. Regulation of pyruvate dehydrogenase complex activity by reversible phosphorylation. Biochem Soc Trans (2003) 31:1143–51. doi: 10.1042/bst0311143
65. Yang Z, Wang Y, Zhang L, Zhao C, Wang D. Phosphorylated form of pyruvate dehydrogenase α1 mediates tumor necrosis factor α-induced glioma cell migration. Oncol Lett (2021) 21:176. doi: 10.3892/ol.2021.12437
66. Ruan X, Zheng J, Liu X, Liu Y, Liu L, Ma J, et al. lncRNA LINC00665 stabilized by TAF15 impeded the malignant biological behaviors of glioma cells via STAU1-mediated mRNA degradation. Mol Ther Nucleic Acids (2020) 20:823–40. doi: 10.1016/j.omtn.2020.05.003
67. Martín-Rufián M, Nascimento-Gomes R, Higuero A, Crisma AR, Campos-Sandoval JA, Gómez-García MC, et al. Both GLS silencing and GLS2 overexpression synergize with oxidative stress against proliferation of glioma cells. J Mol Med (Berl) (2014) 92:277–90. doi: 10.1007/s00109-013-1105-2
68. Szeliga M, Bogacińska-Karaś M, Różycka A, Hilgier W, Marquez J, Albrecht J. Silencing of GLS and overexpression of GLS2 genes cooperate in decreasing the proliferation and viability of glioblastoma cells. Tumour Biol (2014) 35:1855–62. doi: 10.1007/s13277-013-1247-4
69. Huang Q, Lian C, Dong Y, Zeng H, Liu B, Xu N, et al. SNAP25 inhibits glioma progression by regulating synapse plasticity via GLS-Mediated. Glutaminolysis Front Oncol (2021) 11:698835. doi: 10.3389/fonc.2021.698835
70. Appay R, Dehais C, Maurage CA, Alentorn A, Carpentier C, Colin C, et al. CDKN2A homozygous deletion is a strong adverse prognosis factor in diffuse malignant IDH-mutant gliomas. Neuro Oncol (2019) 21:1519–28. doi: 10.1093/neuonc/noz124
71. Pessôa IA, Amorim CK, Ferreira WAS, Sagica F, Brito JR, Othman M, et al. Detection and correlation of single and concomitant TP53, PTEN, and CDKN2A alterations in gliomas. Int J Mol Sci (2019) 20:11. doi: 10.3390/ijms20112658
72. Lu VM, O'Connor KP, Shah AH, Eichberg DG, Luther EM, Komotar RJ, et al. The prognostic significance of CDKN2A homozygous deletion in IDH-mutant lower-grade glioma and glioblastoma: A systematic review of the contemporary literature. J Neurooncol (2020) 148:221–9. doi: 10.1007/s11060-020-03528-2
73. Wirsching HG, Galanis E, Weller M. Glioblastoma. Handb Clin Neurol (2016) 134:381–97. doi: 10.1016/b978-0-12-802997-8.00023-2
74. Phillips HS, Kharbanda S, Chen R, Forrest WF, Soriano RH, Wu TD, et al. Molecular subclasses of high-grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. Cancer Cell (2006) 9:157–73. doi: 10.1016/j.ccr.2006.02.019
75. Kwon SM, Kang SH, Park CK, Jung S, Park ES, Lee JS, et al. Recurrent glioblastomas reveal molecular subtypes associated with mechanistic implications of drug-resistance. PloS One (2015) 10:e0140528. doi: 10.1371/journal.pone.0140528
76. Wood MD, Reis GF, Reuss DE, Phillips JJ. Protein analysis of glioblastoma primary and posttreatment pairs suggests a mesenchymal shift at recurrence. J Neuropathol Exp Neurol (2016) 75:925–35. doi: 10.1093/jnen/nlw068
77. Verhaak RG, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson MD, et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell (2010) 17:98–110. doi: 10.1016/j.ccr.2009.12.020
78. Yoon SJ, Shim JK, Chang JH, Moon JH, Roh TH, Sung KS, et al. Tumor mesenchymal stem-like cell as a prognostic marker in primary glioblastoma. Stem Cells Int (2016) 2016:6756983. doi: 10.1155/2016/6756983
79. Desany B, Zhang Z. Bioinformatics and cancer target discovery. Drug Discovery Today (2004) 9:795–802. doi: 10.1016/s1359-6446(04)03224-6
80. Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res (2009) 37:1–13. doi: 10.1093/nar/gkn923
81. Tian X, Zhu X, Yan T, Yu C, Shen C, Hong J, et al. Differentially expressed lncRNAs in gastric cancer patients: A potential biomarker for gastric cancer prognosis. J Cancer (2017) 8:2575–86. doi: 10.7150/jca.19980
82. Liu S, Cao Q, An G, Yan B, Lei L. Identification of the 3-lncRNA signature as a prognostic biomarker for colorectal cancer. Int J Mol Sci (2020) 21:24. doi: 10.3390/ijms21249359
83. Sui J, Yang S, Liu T, Wu W, Xu S, Yin L, et al. Molecular characterization of lung adenocarcinoma: A potential four-long noncoding RNA prognostic signature. J Cell Biochem (2019) 120:705–14. doi: 10.1002/jcb.27428
84. Moore LM, Holmes KM, Smith SM, Wu Y, Tchougounova E, Uhrbom L, et al. IGFBP2 is a candidate biomarker for Ink4a-arf status and a therapeutic target for high-grade gliomas. Proc Natl Acad Sci USA (2009) 106:16675–9. doi: 10.1073/pnas.0900807106
85. Azar WJ, Zivkovic S, Werther GA, Russo VC. IGFBP-2 nuclear translocation is mediated by a functional NLS sequence and is essential for its pro-tumorigenic actions in cancer cells. Oncogene (2014) 33:578–88. doi: 10.1038/onc.2012.630
86. Zhou YH, Hess KR, Liu L, Linskey ME, Yung WK. Modeling prognosis for patients with malignant astrocytic gliomas: quantifying the expression of multiple genetic markers and clinical variables. Neuro Oncol (2005) 7:485–94. doi: 10.1215/s1152851704000730
87. Zheng S, Houseman EA, Morrison Z, Wrensch MR, Patoka JS, Ramos C, et al. DNA Hypermethylation profiles associated with glioma subtypes and EZH2 and IGFBP2 mRNA expression. Neuro Oncol (2011) 13:280–9. doi: 10.1093/neuonc/noq190
Keywords: glioma, cuproptosis, signature, clusters, bioinformatics
Citation: Wang W, Lu Z, Wang M, Liu Z, Wu B, Yang C, Huan H and Gong P (2022) The cuproptosis-related signature associated with the tumor environment and prognosis of patients with glioma. Front. Immunol. 13:998236. doi: 10.3389/fimmu.2022.998236
Received: 19 July 2022; Accepted: 03 August 2022;
Published: 30 August 2022.
Edited by:
Nan Zhang, Harbin Medical University, ChinaReviewed by:
Xu Qi-Tong, Nanjing Medical University, ChinaZewei Tu, The Second Affiliated Hospital of Nanchang University, China
Copyright © 2022 Wang, Lu, Wang, Liu, Wu, Yang, Huan and Gong. 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: Peipei Gong, bnRuc2dwcEAxNjMuY29t
†These authors have contributed equally to this work