Abstract
Cuproptosis is a new type of cell death that is associated with mitochondrial respiration of the tricarboxylic acid cycle. Previous studies showed that long non-coding RNAs (lncRNAs) regulated low-grade glioma (LGG) progression. However, the potential applications of cuproptosis-related lncRNAs (CRLs) in LGG were not explored. A comprehensive analysis was performed in The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) cohorts. We first screened two distinct cuproptosis subtypes based on prognostic CRLs using consensus clustering. To facilitate individualized survival prediction in LGG, we constructed a prognostic signature (including CRNDE, HAR1A, and FAM181A-AS1) in the TCGA dataset. The prognostic signature exhibited excellent predictive ability and reliability, which was validated in the CGGA_325 and CGGA_693 datasets. Notably, patients in the high-risk group had increased immune cell infiltration and expression of immune checkpoints, which indicated that they may benefit more from immune checkpoint blockade (ICB) therapy. Finally, the prognostic signature screened the population with sensitivity to chemotherapy and ICB therapy. In summary, this study initially explored the mechanism of CRLs in LGG and provides some insights into chemotherapy and ICB therapy of LGG.
Introduction
Gliomas are the most common primary intracranial tumor in the central nervous system and are derived from the neuroglial stem or progenitor cells (Weller et al., 2015). Low-grade glioma (LGG, WHO grades II-III) accounts for 13%–16% of gliomas (). Most patients with LGG inevitably progress to high-grade glioma (HGG) (). Various treatment methods, such as surgery, radiotherapy, chemotherapy, immunotherapy, electric field therapy, and neutron therapy are used, but tumor recurrence appears to be unavoidable (Lu et al., 2021; Li et al., 2022). Therefore, it is urgent to identify novel potential therapeutic targets for LGG.
Copper is an essential metallic element in the human body, and it is linked to a series of biological processes, such as energy metabolism, autophagy, and oxidative stress (). Since its pivotal role in the genesis, severity, and progression of cancer was recognized, copper has attracted more attention in targeted therapy (). A previous study found that copper commonly exhibited anti-tumor effects via the formation of metal chelators and ionophores (Steinbrueck et al., 2020). For the first time, Tsvetkov et al. defined Cu-dependent cytotoxicity, which leads to regulated cell death (RCD), as “cuproptosis” (Tsvetkov et al., 2019; Tsvetkov et al., 2022). Previous research suggested that copper-dependent cell death was correlated with the progression and treatment of glioma (). Therefore, this new form of cell death raises hope for successful future therapies for glioma patients.
Long noncoding RNAs (lncRNAs) are a group of transcripts longer than 200 nucleotides that lack protein-coding capacity (). LncRNAs act as master regulators of gene expression and play vital roles in various biological processes (Peng et al., 2017). LncRNAs are implicated in the progression of stemness, proliferation, angiogenesis, and drug resistance of gliomas (Peng et al., 2018). Up to now, the underlying mechanisms of cuproptosis-related lncRNAs (CRLs) in LGG have not been established.
The present study used genome sequencing technology and bioinformatics analysis to explore the function of CRLs in LGG patients. According to the prognostic CRLs screened from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) cohorts, LGG patients were stratified into two cuproptosis subtypes with distinct over survival, clinicopathological features, immune cell infiltration, and biological processes. To quantify the difference between individuals, we established a new prognostic signature based on three CRLs (including CRNDE, HAR1A, and FAM181A-AS1). This prognostic signature had a powerful ability to predict the prognosis of LGG patients. We initially demonstrated the potential of our prognostic signature in predicting immune checkpoint blockade (ICB) therapy and chemotherapy sensitivity.
Materials and methods
Data collection and pre-processing
The transcriptional data and clinical information of LGG patients were sourced from the TCGA (https://portal.gdc.cancer.gov) and CGGA (http://www.cgga.org.cn/) datasets. The gene expression profiles of 105 normal brain tissues from the Genome Tissue Expression (GTEx) project was obtained from the UCSC Xena website (https://xena.ucsc.edu/). LGG samples were excluded when they had missing survival information and definitive histopathological diagnosis. Finally, a total of 509 LGG samples from the TCGA dataset and 591 samples from the CGGA dataset (172 from the CGGA_325 cohort and 419 from the CGGA_693 cohort) were filtered for further study. The transcriptional data from the TCGA and CGGA cohorts were TPM normalized. The characteristics of the LGG patients are listed in Supplementary Table S1.
Identification of the CRLs in the TCGA and CGGA_325 datasets
Seventeen cuproptosis-related genes (CRGs, Supplementary Table S2) were retrieved based on a published study (Tsvetkov et al., 2022). To screen the CRLs, Pearson correlation analysis was performed between the CRGs and lncRNAs in the TCGA and CGGA_325 cohorts. An absolute value of correlation coefficients >0.5 corresponding to a p-value < 0.05 was considered eligible.
Consistent clustering to determine the cuproptosis subtypes
Univariate Cox regression analysis was used to screen the prognostic CRLs. The prognostic CRLs shared by the TCGA and CGGA_325 cohorts were selected for unsupervised clustering according to the R package “ConsensusClusterPlus”. The candidate cluster number ranged from two to six. To produce the most stable consensus matrix, the procedure was repeated 100 times with an 80% sample selected in each iteration (Wilkerson and Hayes, 2010).
Evaluation of the immune landscape
Infiltrating immune and stromal cells are mainly structural components of the tumor microenvironment (TME). We used the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm to calculate the abundance of immune and stromal cells in each LGG patient (Yoshihara et al., 2013).
To better understand immune cell infiltration in the TME of LGG samples, the Tumor Immune Estimation Resource (TIMER) algorithm (Li et al., 2017) was used in the “IOBR” R package (Zeng et al., 2021). For validation, the single-sample gene set enrichment analysis (ssGSEA) algorithm was also used to evaluate the relative abundance of infiltrating immune cells. Jia et al. (2018)summarized the characteristic gene panels for TME-infiltrating immune cell types (Supplementary Table S3).
Tumor Immunophenotype Profiling (TIP, http://biocc.hrbmu.edu.cn/TIP/index.jsp) is a web-based tool that makes the process of anticancer immunity easy to visualize (Xu et al., 2018). According to the TIP, the anti-tumor immune response may be simplified as a seven-step cycle event (Supplementary Table S4): the release of cancer cell antigens (Step 1), cancer antigen presentation (Step 2), priming and activation (Step 3), trafficking of immune cells to tumors (Step 4), infiltration of immune cells into tumors (Step 5), recognition of cancer cells by T cells (Step 6), and killing of cancer cells (Step 7). In this study, we assessed the anti-tumor activity score of each LGG sample in the TCGA database.
Functional enrichment analyses
Gene set variation analysis (GSVA) is a functional enrichment analysis that estimates the difference in pathway activity of the samples using an unsupervised method (). To investigate the difference in biological processes between distinct clusters, we performed GSVA using the “GSVA” package in R. Kyoto Encyclopedia of Genes and Genomes (KEGG) gene sets were obtained from the Molecular Signatures Database (MSigDB, http://www.gsea-msigdb.org/gsea/index.jsp, v7.5.1) (Liberzon et al., 2011). KEGG pathways between different subtypes with a false discovery rate (FDR) < 0.05 were considered significant. Sixteen gene sets that represent classical biological processes were also analyzed for validation (Supplementary Table S5) (Mariathasan et al., 2018).
Chemotherapeutic and immune checkpoint blockade therapy response prediction
LGG patients greatly benefit from chemotherapy (Shaw et al., 2012; Viaccoz et al., 2012). According to the Genomics of Drug Sensitivity in Cancer (GDSC, https://www.cancerrxgene.org/) database, the “pRRophetic” package in R was used to evaluate the chemotherapy drug sensitivity of each LGG sample (). Four commonly used chemotherapeutics (temozolomide, tamoxifen, bleomycin, and vinblastine) were selected in the present study.
The subclass mapping method from GenePattern (SubMap, https://www.genepattern.org/) was used to predict the response to ICB therapy (). The SubMap method is an unsupervised algorithm that reveals common subtypes between independent datasets. The gene expression profiles of 47 melanoma patients who received ICB therapy were sourced from a published dataset (Roh et al., 2017).
Identification of differentially expressed cuproptosis-related lncRNAs between cuproptosis subtypes
According to the “limma” package in R (Ritchie et al., 2015), we further analyzed the differential expression of 37 prognostic CRLs between distinct cuproptosis subgroups. Differentially expressed CRLs (DE-CRLs) with an FDR <0.05 and | log2fold change (FC) | values >1 were considered significant.
Construction and validation of cuproptosis-related lncRNAs signature
The overlapping DE-CRLs from the TCGA and CGGA_325 cohorts were incorporated into the least absolute shrinkage and selection operator (LASSO) regression analysis using the “glmnet” R package (). The risk score for each sample was calculated by the expression level and regression coefficient of each CRL. The formula is described below:
Coef (Xi) is the regression coefficient of the CRLs, and Exp (Xi) represents the expression levels of CRLs. LGG samples were split into high- and low-risk score groups by the median value.
The Kaplan-Meier (K-M) survival curve was employed to compare the overall survival between the high- and low-risk groups. The receiver operating characteristic (ROC) curves were drawn to measure the predictive power of cuproptosis-related lncRNAs signature through the “survival-ROC” R package ().
Nomogram construction
Univariate and multivariate Cox regression analyses were used to evaluate the independent prognostic value of cuproptosis-related lncRNA signature. To evaluate the Cox regression model, the PH hypothesis test was performed. Nomogram was constructed to predict the survival probability of 1-, 3-, and 5-years overall survival (). The calibration plot was also performed to verify the accuracy.
Real-time quantitative PCR
Four tumor tissues and corresponding peritumoral brain tissues from LGG patients were collected in the Second Affiliated Hospital of Harbin Medical Universit. This research was approved by all the patients and the Ethics Committee of the hospital. Total RNA was isolated from LGG tissues using TRIzol reagent (Invitrogen, United States) according to the manufacturer’s protocol. According to the manufacturer’s instructions of the Nanodrop ND-2000 spectrophotometer (Thermo Scientific, United States), 2 μg of the total RNA was transcribed into cDNA. RT-qPCR was performed with the SYBR Green PCR kit (Takara, Japan). Independent experiments were conducted in triplicate, and ACTB served as an internal control. The primers (Tsingke Biotechnology Co., Ltd, Beijing, China) were used are displayed in Supplementary Table S6.
Statistical analysis
All statistical analyses and visualizations were executed in R 4.1.2 and GraphPad Prism 8.0.2. Survival analysis was performed using the “survival” package in R. The nomogram and the calibration curves were generated by the “rms” R package (Zhang and Kattan, 2017). The R package “maftools” was used to process and present the mutation data (Mayakonda et al., 2018). Immune cell infiltration, ssGSEA score, mRNA expression, and TIP score were compared between the two groups using the Wilcoxon test. The student’s t-test was used to analyze differences between different risk groups. P < 0.05 was considered statistically significant.
Results
Identification of prognostic CRLs from TCGA and CGGA_325 datasets
Supplementary Figure S1 provides a flow chart of this study. Pearson correlation analysis was performed between CRGs and CRLs. Combined with univariate Cox regression analyses, a total of 517 and 137 CRLs were obtained from the TCGA and CGGA_325 datasets, respectively, using | cor | > 0.5 and p < 0.05 as cutoff values (Supplementary Tables S7, S8). Thirty-seven prognostic CRLs that were shared by the two datasets were screened (Figure 1A).
FIGURE 1
Determine the cuproptosis subtypes in LGG patients
To better examine the role of cuproptosis in LGG patients, unsupervised consensus clustering of the 37 prognostic CRLs was performed using the R package “Consensus ClusterPlus”. As shown in the cumulative distribution function (CDF) plots and consensus matrix, k = 2 was the most suitable choice (Figures 1B,C). We obtained two cuproptosis subtypes from the TCGA dataset, which were labeled cluster-1 and cluster-2. A total of 263 samples were classified into cluster-1, and 246 samples were classified into cluster-2 (Supplementary Table S9). To confirm the robustness of unsupervised clustering, the same algorithm was used in the CGGA_325 cohort (Supplementary Figures S2A,B). We also obtained two distinct subtypes. Eighty samples were classified into cluster-1, and 92 samples were classified into cluster-2 (Supplementary Table S10). These results demonstrated the effectiveness of our grouping.
K-M survival curves showed that patients in cluster-2 had a significant survival advantage (Figure 1D; Supplementary Figure S2C). The distribution of various clinical factors and 37 CRL expression levels between different subgroups are intuitively shown in Figure 1E. Compared to cluster-2, patients in cluster-1 were more related to the clinicopathological features of WHO III, isocitrate dehydrogenase (IDH) wild-type, unmethylated O6-methylguanine-DNA methyltransferase promoter (MGMTp), and 1p19q non-codeletion. A similar result was observed in the CGGA_325 cohort (Supplementary Figure S2D).
Tumor microenvironment immune cell infiltration in different cuproptosis subtypes
We examined the difference in immune cell infiltration between different cuproptosis subgroups. The immune and stromal scores positively correlated with the number of immune or stromal components in the TME. The ESTIMATE score represents the comprehensive proportion of the immune and stromal scores in the TME (Yoshihara et al., 2013). The immune score, stromal score, and ESTIMATE score were significantly higher in cluster-1 in this study, which indicates an increased infiltration of immune and stromal cells (Figures 2A–C; Supplementary Figures S3A–C). Based on the TIMER algorithm, we calculated the abundance of six immune cells, including B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and myeloid dendritic cells, in the TME. The infiltration levels of six immune cells in cluster-1 were significantly increased (Figure 2D; Supplementary Figure S3D). Patients with higher immune cell infiltration corresponded with poor outcomes in LGG (Figure 2E; Supplementary Figure S3E). We also evaluated the infiltration levels of 28 immune cells using the ssGSEA score to validate these results (Figure 2F; Supplementary Figure S3F).
FIGURE 2
To better understand the anti-tumor immune response in LGG, the process was visualized using the TIP website. Patients in cluster-1 had increased immune activity scores in Step 1, Step 4, and Step 5. However, the immune activity scores in Step 3, Step 6, and Step 7 of cluster-2 were higher (Figure 2G). Although the recruitment and infiltration of immune cells in cluster-1 were higher, the abilities to recognize and kill cancer cells were lower than in cluster-2. This difference also explained why the increasing level of immune infiltration was associated with shorter overall survival in LGG.
Characteristics of the biological process in distinct cuproptosis subtypes
To further examine the differences in biological processes between distinct cuproptosis subtypes, we performed the GSVA enrichment analysis. Cluster-1 was markedly enriched in stromal, immune activation, oncogenic, and DNA damage repair (DDR) pathways, such as the ECM receptor interaction, antigen processing and presentation, TGF-β signaling pathway, and DNA replication. Patients in cluster-2 were enriched in pathways related to metabolism, including taurine and hypotaurine metabolism and the biosynthesis of unsaturated fatty acids (Figure 2H; Supplementary Table S11). Differences in 16 typical biological processes between cluster-1 and cluster-2 were also identified in this study. The stromal, immune activation, and DDR pathways were enriched in cluster-1 (Figure 2I). These results were verified in the CGGA_325 cohort (Supplementary Figures S3G,H, Supplementary Table S12).
Chemotherapy sensitivity and the immune checkpoint blockade treatment responsiveness in different cuproptosis subtypes
Chemotherapy is the main adjuvant therapy for LGG, and it provided great benefit. The present study assessed the half-maximal inhibitory concentration (IC50) of four commonly used chemotherapeutics (temozolomide, tamoxifen, bleomycin, and vinblastine) in each LGG sample. Patients in cluster-1 had a higher sensitivity for four chemotherapeutics (Figures 3A,B).
FIGURE 3
The abundance of infiltrating immune cells and enrichment of immune checkpoint pathways in cluster-1 indicated a positive response to ICB therapy. Considering the importance of ICB therapy in LGG, we next examined the response to ICB therapy in LGG. The results demonstrated that patients in cluster-1 were more responsive to anti-PD-1 treatment (Figures 3C,D).
Construction of the cuproptosis-related lncRNAs signature
With an FDR <0.05 and | logFC | > 1, we obtained eleven and thirteen DE-CRLs from the TCGA and CGGA_325 datasets, respectively (Supplementary Table S13). Eleven overlapped prognostic CRLs were selected from the two cohorts (Figure 4A). Considering the individual heterogeneity, we established a CRL-based prognostic model to quantify the difference between the individuals in LGG (Figures 4B,C). The regression coefficients of the three CRLs are shown in Supplementary Table S14. We divided the LGG patients into low- and high-risk score groups at the median cut-off. Patients in cluster-1 showed an increased risk score compared to cluster-2 (Supplementary Figures S4A, S5A). The clinicopathological features of age ≤ 40 years, WHO II, 1p19q codeletion, IDH mutation, and MGMTp methylation were associated with a decreased risk score (Supplementary Figures S4B–E, S5B–E, S6B–E).
FIGURE 4
K-M curves indicated that patients with a high-risk score had a poor prognosis (Figure 4E). The distribution plot of the risk score and survival status showed that patients in the low-risk group had a higher survival rate (Figure 4F). Principal component analysis (PCA) showed that the LGG patients were easily distinguished according to the different risk groups (Figure 4G). The AUC values of the signature for predicting 1-, 3-, and 5-years survival rates in the TCGA dataset were 0.841, 0.839, and 0.768, respectively. The AUC values were 0.773 (1-year), 0.834 (3-years), and 0.849 (5-years) in the CGGA_325 cohort, and the AUC values in the CGGA_693 cohort were 0.699, 0.725, and 0.740 for 1, 3, and 5 years, respectively. (Figure 4H).
We further evaluated the predictive ability of the prognostic signature between different clinical subgroups. The results showed that patients in the low-risk group always had a better outcome than patients in the high-risk group (Figure 5A–L; Supplementary Figures S7A–L, S8A–L). These findings demonstrated the favorable predictive ability of the prognostic signature.
FIGURE 5
Building a predictive nomogram for overall survival prediction
The results of univariate and multivariate Cox regression analyses showed that the risk score was always an independent prognostic factor in LGG (Figures 6A,B; Supplementary Figures S9A,B; Supplementary Figures S10A,B). After the PH hypothesis test (Figure 6C; Supplementary Figures S9C, S10C), a nomogram prediction model was constructed (Figure 6D; Supplementary Figures S9D, S10D). We found that the calibration curves of the nomogram were close to the standard curves (Figure 6E; Supplementary Figures S9E, S10E). These results also indicated the clinical applicability of the nomogram.
FIGURE 6
Tumor microenvironment immune cell infiltration in different risk groups
The difference in immune landscape between different risk groups was also assessed in this research. Patients in the high-risk group had a higher immune score, stromal score, and ESTIMATE score (Figures 7A–C; Supplementary Figures S11A–C, S12A–C). The abundance of immune cell infiltration in the high-risk group was also increased compared to the low-risk group (Figures 7D,E; Supplementary Figures S11D,E, Supplementary Figures S12D,E). We examined the expression levels of immune checkpoint (ICP) regulators including CD80, CD86, CD274 (PD-L1), IDO1, CTLA4, HAVCR2 (TIM-3), LAG3, PDCD1 (PD-1), and PDCD1LG2 (PD-L2) in the different risk groups. We found that most of the ICPs were upregulated in the high-risk group, which indicated that more benefits may be gained from ICB therapy (Figure 7F; Supplementary Figures S11F, S12F).
FIGURE 7
Relationship between prognostic signature and mutational status
Tumor mutation burden (TMB) has been used to predict prognosis and ICB efficacy across many cancer types (Snyder et al., 2014; Rizvi et al., 2015; ). To better understand the role of TMB in LGG, we calculated the TMB level of each sample in LGG. Patients in the high-risk group had a higher TMB level (Figure 8A). The TMB levels increased with the risk score (Figure 8B). We also found that the TMB level negatively correlated with overall survival (Figure 8C). Patients in the low-risk group showed an apparent survival advantage in the high and low TMB groups (Figure 8D).
FIGURE 8
The top 20 mutated genes in different risk groups were visualized using waterfall plots (Figures 8E,F). IDH mutation is an important factor that is associated with favorable outcomes in LGG patients (). Glioblastoma with wild-type IDH may directly transform from LGG with wild-type IDH (). Among the top 20 mutated genes, patients in the low-risk group (IDH1-89%, IDH2-7%) had a higher IDH mutation rate than patients in the high-risk group (IDH1-66%). This result supports the central role of IDH mutation in LGG. CIC mutation was also highly expressed in the low-risk group, which is associated with better survival of glioma (). The mutation rates of potential targets, including TP53, ATRX, EGFR, and TTN, were higher in the high-risk group (; ), which indicated a positive response to ICB therapy.
Chemotherapy sensitivity and the immune checkpoint blockade treatment responsiveness in different risk groups
We assessed the predictive value of the risk score for chemotherapy and ICB therapy to help construct individualized treatment plans. The IC50 values of the four chemotherapy drugs showed that patients in the high-risk group were more sensitive to chemotherapy (Figures 9A–C). As shown in Figure 9D, patients in the high-risk group may have greater benefits from ICB therapy.
FIGURE 9
Validation of the 3 CRLs expressions in LGG tissue samples
We further characterized these 3 CRLs based on gene expression and survival prognosis using the TCGA and GTEx data. As is shown in Figures 10A–C, CRNDE, and FAM181A-AS1 were expressed at high levels in LGG tissues, while HAR1A was increased in normal brain tissues. K-M curves proved the clinical outcomes between different expression levels of 3 CRLs (Figures 10D–F). We further verified the expression levels of 3 CRLs between LGG and peritumoral brain tissues (PBT) with RT-qPCR (Figures 10G–I).
FIGURE 10
Discussion
Compared to the most malignant glioblastoma, the median survival of LGG is longer. However, LGG generally converts into high-grade glioma approximately 4–5 years after diagnosis despite receiving standard care (Kumthekar et al., 2015). Cuproptosis is a form of copper-dependent cell death that is different from known cell death pathways. Dysregulated copper metabolism has been found in many types of cancer, which suggests an irreplaceable characteristic in the development of cancer (Kucharzewski et al., 2003; Jouybari et al., 2020). LncRNAs participate important roles in tumor incidence, development, and metastasis (Tang et al., 2022). The underlying mechanism of CRLs in LGG is not clear. To the best of our knowledge, this report is the first study to assess the correlation between CRLs and biological and clinical features of LGG using a bioinformatics method.
The present study first screened 37 prognostic CRLs from the TCGA and CGGA cohorts using the CRGs. According to the consistent clustering, samples in the two datasets were divided into two subtypes with different expression levels of CRLs. Patients in different cuproptosis subgroups had distinct prognostic and clinicopathological features. Cluster-1 was enriched in the stromal, immune activation, oncogenic, and DDR pathways. In contrast, patients in cluster-2 were more active in pathways related to metabolism. Therefore, combined with a previous study on cuproptosis in glioma, patients in cluster-1 were regarded as the cuproptosis-resistant type. Cluster-2 patients were characterized as a cuproptosis-sensitive type.
We also found that patients in cluster-1 had increased immune cell infiltration. The abundance of immune cells was negative for the overall survival of LGG. One of the most important functions of dendritic cells (DCs) is antigen uptake and processing. Using the antigen-processing machinery (APM), tumor-derived epitopes are cross-presented to T cells by DCs (Yamanaka, 2009). CD8+ cytotoxic T cells and NK cells directly target and kill tumor cells (Kim and Cantor, 2014). However, the TIP analysis showed that patients in cluster-2 had a decreased immune activity score for the infiltration of immune cells but an increased score in the steps of cancer cell recognition and killing. This result corresponds to the better outcomes for cluster-2 patients. Glioma is characterized by a “clod” tumor, and immune cells are trapped in the stroma rather than penetrating the tumor parenchyma (; ). Combined with the TIP findings, we found that the effectiveness of anti-tumor immune response may depend on the “effective immune cells”. Therefore, methods to increase “effective immune cells” in LGG will be a new direction of immunotherapy.
Based on the DE-CRLs between different cuproptosis subgroups, we then constructed a 3 CRL prognostic signature, including CRNDE, HAR1A, and FAM181A-AS1, to predict the individualized treatment through the LASSO regression. CRNDE is involved in cancer progression, neuronal differentiation, gametogenesis, and other developmental processes (). CRNDE also plays a regulatory role in temozolomide chemoresistance to glioma (Zhao et al., 2021). HAR1A acts as a tumor suppressor in many cancer types. demonstrated that lower HAR1A expression may result in a worse outcome for glioma patients (Waters et al., 2021). The expression level of lncRNA FAM181A-AS1 correlates with advanced tumor stage and survival of glioma (Jiang and Chen, 2020).
We further demonstrated that patients in the high- and low-risk groups exhibited significantly different survival outcomes. The AUC values of the ROC plots demonstrated the reliability and accuracy of the prognostic signature. The effectiveness and stability of the prognostic signature were validated in the TCGA and CGGA cohorts. Univariate and multivariate Cox regression analyses showed that the risk score was an independent prognostic factor in LGG. Nomogram and calibration plots also showed that the prognostic signature had excellent predictive power.
As an important component of LGG treatment, chemotherapy has attracted more attention. The present research evaluated four chemotherapeutic drug sensitivities in each LGG sample. Temozolomide is a classic chemotherapeutic drug that is used for the treatment of glioma. It is characterized as easy to administer and better tolerated (van den Bent, 2015). Tamoxifen is an estrogen receptor modulator that induces cell death in glioma (). The combined administration of tamoxifen and temozolomide was well tolerated (). Bleomycin functions by inducing single- and double-stranded DNA breaks, which is similar to ionizing radiation (Mathews et al., 2012). Research showed that the effect of bleomycin on F98 glioma cells was stronger than temozolomide (). One large study revealed that vinblastine showed low toxicity and maintained the quality of life in pediatric LGG (Lassaletta et al., 2016). The predicted IC50 values of the above four chemotherapeutic drugs are lower in cluster-1 or high-risk groups, which indicated more sensitivity to the four chemotherapeutic drugs. According to the RTOG 9802 results, LGG patients had increased overall survival after adjuvant chemotherapy and radiation. The chemotherapy drugs procarbazine, lomustine, and vincristine (PCV) were the mainstream treatment strategy for improving LGG patients’ survival rates (). However, we could not assess the sensitivity to PCV in LGG patients due to the limited candidate drugs in the predictive algorithm.
Although ICB therapy has become a promising treatment strategy against a variety of tumors, only a minority of patients obtain favorable benefits from it (). PD-1, PD-L1, and CTLA-4 are the main immune checkpoint molecules in glioma immunotherapy (Topalian et al., 2015). However, response to anti-PD-1/PD-L1 therapy commonly ranges from 10 to 40% (Zou et al., 2016). Due to the characterization of glioma as a “cold tumor”, the response frequency may be lower (). A clinical trial of recurrent glioma showed that only 8% of patients exhibited dramatic responses to anti-PD-1 therapy (). Therefore, preliminary screening for patients who are sensitive to ICB therapy using a bioinformatics method may provide optimal clinical treatment. The present study found that LGG patients with higher immune cell infiltration, and TMB levels, may produce more clinical responses to ICB therapy. PD-1 and its ligands PD-L1/PD-L2 are the most comprehensively studied immune checkpoint molecules (Qi et al., 2020). CTLA-4 is the first ICP molecule that was used in ICB therapy, its also effect on the CD80 and CD86 expressed by dendritic cells (Ohue and Nishikawa, 2019). Research found that combining select therapies with IDO1, LAG-3, and TIM-3 blockade tend to benefit against tumor growth (; Kim et al., 2017; Zhai et al., 2018). In this study, the expression of these ICP regulators was increased in the high-risk group. The above findings may provide a means for the effective application of ICB therapy for LGG.
Here are some limitations in the current work. First, the sample size was relatively small as a validation group. Second, we selected CRGs from recently published research. With further study of cuproptosis, an increasing number of CRGs may be identified. Finally, this study was based on bioinformatics analysis, and further experimental studies in vitro and in vivo are needed.
Conclusion
In summary, we identified two cuproptosis subtypes in LGG with different outcomes, clinicopathological features, and immune landscapes. We also constructed and validated a cuproptosis-related signature that exhibited robust capacity in predicting the survival outcomes of LGG patients. Notably, we also evaluated chemotherapy sensitivity and ICB treatment responsiveness in LGG patients. This study provides a new reference for the chemotherapy and ICB treatment of LGG and may be beneficial in individualized treatment strategies.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Committee of the Second Affiliated Hospital of Harbin Medical University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
XY, YL, ZL, and SH conceived and designed the study and drafted the manuscript. HZ, JD, and ZL provided analytical technical support. FW, NW, and JZ participated in the production of charts and pictures. All authors have read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (No. 61575058) and the Talent Introduction Project of Zhejiang Provincial People’s Hospital (No. C-2021-QDJJ03-01).
Acknowledgments
The authors gratefully acknowledge databases like TCGA, CGGA, and MSigDB for offering convenient access to datasets. In addition, we thank Siwen Wang for the help in data analysis and encouragement.
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/fmolb.2022.966843/full#supplementary-material
References
1
CarrabbaG.LanfranchiG.MenghettiC.RampiniP.CaroliM. (2013). Continuous tamoxifen and dose-dense temozolomide in recurrent glioblastoma. Anticancer Res.33, 3383–3389. PubMed Abstract | Google Scholar
2
BellE. H.ZhangP.ShawE. G.BucknerJ. C.BargerG. R.BullardD. E.et al (2020). Comprehensive genomic analysis in nrg oncology/RTOG 9802: A phase III trial of radiation versus radiation plus procarbazine, lomustine (CCNU), and vincristine in high-risk low-grade glioma. J. Clin. Oncol.38, 3407–3417. 10.1200/jco.19.02983PubMed Abstract | CrossRef Full Text | Google Scholar
3
BratD. J.VerhaakR. G.AldapeK. D.YungW. K.SalamaS. R.CooperL. A.et al (2015). Comprehensive, integrative genomic analysis of diffuse lower-grade gliomas. N. Engl. J. Med.372, 2481–2498. 10.1056/NEJMoa1402121PubMed Abstract | CrossRef Full Text | Google Scholar
4
BuccarelliM.D'AlessandrisQ. G.MatarreseP.MollinariC.SignoreM.CappanniniA.et al (2021). Elesclomol-induced increase of mitochondrial reactive oxygen species impairs glioblastoma stem-like cell survival and tumor growth. J. Exp. Clin. Cancer Res.40, 228. 10.1186/s13046-021-02031-4PubMed Abstract | CrossRef Full Text | Google Scholar
5
ChenD. S.MellmanI. (2017). Elements of cancer immunity and the cancer-immune set point. Nature541, 321–330. 10.1038/nature21349PubMed Abstract | CrossRef Full Text | Google Scholar
6
ChenR.CohenA. L.ColmanH. (2016). Targeted therapeutics in patients with high-grade gliomas: Past, present, and future. Curr. Treat. Options Oncol.17, 42. 10.1007/s11864-016-0418-0PubMed Abstract | CrossRef Full Text | Google Scholar
7
ChenR.Smith-CohnM.CohenA. L.ColmanH. (2017). Glioma subclassifications and their clinical significance. Neurotherapeutics14, 284–297. 10.1007/s13311-017-0519-xPubMed Abstract | CrossRef Full Text | Google Scholar
8
ChenY.GuoY.ChenH.MaF. (2020). Long non-coding RNA expression profiling identifies a four-long non-coding RNA prognostic signature for isocitrate dehydrogenase mutant glioma. Front. Neurol.11, 573264. 10.3389/fneur.2020.573264PubMed Abstract | CrossRef Full Text | Google Scholar
9
ChoiY. S.BaeS.ChangJ. H.KangS. G.KimS. H.KimJ.et al (2021). Fully automated hybrid approach to predict the IDH mutation status of gliomas via deep learning and radiomics. Neuro. Oncol.23, 304–313. 10.1093/neuonc/noaa177PubMed Abstract | CrossRef Full Text | Google Scholar
10
DenoyerD.ClatworthyS. A. S.CaterM. A. (2018). 16. Copper complexes in cancer therapy. Mater. Ions Life Sci., 469–506. 10.1515/9783110470734-022CrossRef Full Text | Google Scholar
11
DuffauH. (2018). Diffuse low-grade glioma, oncological outcome and quality of life: A surgical perspective. Curr. Opin. Oncol.30, 383–389. 10.1097/cco.0000000000000483PubMed Abstract | CrossRef Full Text | Google Scholar
12
FilleyA. C.HenriquezM.DeyM. (2017). Recurrent glioma clinical trial, CheckMate-143: The game is not over yet. Oncotarget8, 91779–91794. 10.18632/oncotarget.21586PubMed Abstract | CrossRef Full Text | Google Scholar
13
FriedmanJ.HastieT.TibshiraniR. (2010). Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw.33, 1–22. 10.18637/jss.v033.i01PubMed Abstract | CrossRef Full Text | Google Scholar
14
GeE. J.BushA. I.CasiniA.CobineP. A.CrossJ. R.DeNicolaG. M.et al (2022). Connecting copper and cancer: From transition metal signalling to metalloplasia. Nat. Rev. Cancer22, 102–113. 10.1038/s41568-021-00417-2PubMed Abstract | CrossRef Full Text | Google Scholar
15
GederaasO. A.HaugeA.EllingsenP. G.BergK.AltinD.BardalT.et al (2015). Photochemical internalization of bleomycin and temozolomide--in vitro studies on the glioma cell line F98. Photochem. Photobiol. Sci.14, 1357–1366. 10.1039/c5pp00144gPubMed Abstract | CrossRef Full Text | Google Scholar
16
GeeleherP.CoxN.HuangR. S. (2014). pRRophetic: an R package for prediction of clinical chemotherapeutic response from tumor gene expression levels. PloS one9, e107468. 10.1371/journal.pone.0107468PubMed Abstract | CrossRef Full Text | Google Scholar
17
GibbE. A.BrownC. J.LamW. L. (2011). The functional role of long non-coding RNA in human carcinomas. Mol. Cancer10, 38. 10.1186/1476-4598-10-38PubMed Abstract | CrossRef Full Text | Google Scholar
18
HanP.LiJ. W.ZhangB. M.LvJ. C.LiY. M.GuX. Y.et al (2017). The lncRNA CRNDE promotes colorectal cancer cell proliferation and chemoresistance via miR-181a-5p-mediated regulation of Wnt/β-catenin signaling. Mol. Cancer16, 9. 10.1186/s12943-017-0583-1PubMed Abstract | CrossRef Full Text | Google Scholar
19
HänzelmannS.CasteloR.GuinneyJ. (2013). Gsva: Gene set variation analysis for microarray and RNA-seq data. BMC Bioinforma.14, 7. 10.1186/1471-2105-14-7CrossRef Full Text | Google Scholar
20
HarmalkarM.UpraityS.KaziS.ShirsatN. V. (2015). Tamoxifen-Induced cell death of malignant glioma cells is brought about by oxidative-stress-mediated alterations in the expression of BCL2 family members and is enhanced on miR-21 inhibition. J. Mol. Neurosci.57, 197–202. 10.1007/s12031-015-0602-xPubMed Abstract | CrossRef Full Text | Google Scholar
21
HeagertyP. J.LumleyT.PepeM. S. (2000). Time-dependent ROC curves for censored survival data and a diagnostic marker. Biometrics56, 337–344. 10.1111/j.0006-341x.2000.00337.xPubMed Abstract | CrossRef Full Text | Google Scholar
22
HoshidaY.BrunetJ. P.TamayoP.GolubT. R.MesirovJ. P. (2007). Subclass mapping: Identifying common subtypes in independent disease data sets. PloS one2, e1195. 10.1371/journal.pone.0001195PubMed Abstract | CrossRef Full Text | Google Scholar
23
HsuE. J.CaoX.MoonB.BaeJ.SunZ.LiuZ.et al (2021). A cytokine receptor-masked IL2 prodrug selectively activates tumor-infiltrating lymphocytes for potent antitumor therapy. Nat. Commun.12, 2768. 10.1038/s41467-021-22980-wPubMed Abstract | CrossRef Full Text | Google Scholar
24
HuangR. Y.EppolitoC.LeleS.ShrikantP.MatsuzakiJ.OdunsiK. (2015). LAG3 and PD1 co-inhibitory molecules collaborate to limit CD8+ T cell signaling and dampen antitumor immunity in a murine ovarian cancer model. Oncotarget6, 27359–27377. 10.18632/oncotarget.4751PubMed Abstract | CrossRef Full Text | Google Scholar
25
HugoW.ZaretskyJ. M.SunL.SongC.MorenoB. H.Hu-LieskovanS.et al (2016). Genomic and transcriptomic features of response to anti-PD-1 therapy in metastatic melanoma. Cell.165, 35–44. 10.1016/j.cell.2016.02.065PubMed Abstract | CrossRef Full Text | Google Scholar
26
HwangI.PanH.YaoJ.ElementoO.ZhengH.PaikJ. (2020). CIC is a critical regulator of neuronal differentiation. JCI insight, 135826. 10.1172/jci.insight.135826PubMed Abstract | CrossRef Full Text | Google Scholar
27
IasonosA.SchragD.RajG. V.PanageasK. S. (2008). How to build and interpret a nomogram for cancer prognosis. J. Clin. Oncol.26, 1364–1370. 10.1200/jco.2007.12.9791PubMed Abstract | CrossRef Full Text | Google Scholar
28
JacksonC. M.ChoiJ.LimM. (2019). Mechanisms of immunotherapy resistance: Lessons from glioblastoma. Nat. Immunol.20, 1100–1109. 10.1038/s41590-019-0433-yPubMed Abstract | CrossRef Full Text | Google Scholar
29
JiaQ.WangJ.HeN.HeJ.ZhuB. (2019). Titin mutation associated with responsiveness to checkpoint blockades in solid tumors. JCI insight, 127901. 10.1172/jci.insight.127901PubMed Abstract | CrossRef Full Text | Google Scholar
30
JiaQ.WuW.WangY.AlexanderP. B.SunC.GongZ.et al (2018). Local mutational diversity drives intratumoral immune heterogeneity in non-small cell lung cancer. Nat. Commun.9, 5361. 10.1038/s41467-018-07767-wPubMed Abstract | CrossRef Full Text | Google Scholar
31
JiangX.ChenD. (2020). LncRNA FAM181A-AS1 promotes gliomagenesis by sponging miR-129-5p and upregulating ZRANB2. Aging12, 20069–20084. 10.18632/aging.103391PubMed Abstract | CrossRef Full Text | Google Scholar
32
JouybariL.KianiF.IslamiF.SanagooA.SayehmiriF.HosnedlovaB.et al (2020). Copper concentrations in breast cancer: A systematic review and meta-analysis. Curr. Med. Chem.27, 6373–6383. 10.2174/0929867326666190918120209PubMed Abstract | CrossRef Full Text | Google Scholar
33
KimH. J.CantorH. (2014). CD4 T-Cell subsets and tumor immunity: The helpful and the not-so-helpful. Cancer Immunol. Res.2, 91–98. 10.1158/2326-6066.Cir-13-0216PubMed Abstract | CrossRef Full Text | Google Scholar
34
KimJ. E.PatelM. A.MangravitiA.KimE. S.TheodrosD.VelardeE.et al (2017). Combination therapy with anti-PD-1, anti-TIM-3, and focal radiation results in regression of murine gliomas. Clin. Cancer Res.23, 124–136. 10.1158/1078-0432.Ccr-15-1535PubMed Abstract | CrossRef Full Text | Google Scholar
35
KucharzewskiM.BraziewiczJ.MajewskaU.GózdzS. (2003). Selenium, copper, and zinc concentrations in intestinal cancer tissue and in colon and rectum polyps. Biol. Trace Elem. Res.92, 1–10. 10.1385/bter:92:1:1PubMed Abstract | CrossRef Full Text | Google Scholar
36
KumthekarP.RaizerJ.SinghS. (2015). Low-grade glioma. Cancer Treat. Res.163, 75–87. 10.1007/978-3-319-12048-5_5PubMed Abstract | CrossRef Full Text | Google Scholar
37
LassalettaA.ScheinemannK.ZelcerS. M.HukinJ.WilsonB. A.JabadoN.et al (2016). Phase II weekly vinblastine for chemotherapy-naïve children with progressive low-grade glioma: A Canadian pediatric brain tumor consortium study. J. Clin. Oncol.34, 3537–3543. 10.1200/jco.2016.68.1585PubMed Abstract | CrossRef Full Text | Google Scholar
38
LiS.GaoP.DaiX.YeL.WangZ.ChengH. (2022). New prognostic biomarker CMTM3 in low grade glioma and its immune infiltration. Ann. Transl. Med.10, 206. 10.21037/atm-22-526PubMed Abstract | CrossRef Full Text | Google Scholar
39
LiT.FanJ.WangB.TraughN.ChenQ.LiuJ. S.et al (2017). Timer: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res.77, e108–e10. 10.1158/0008-5472.Can-17-0307PubMed Abstract | CrossRef Full Text | Google Scholar
40
LiberzonA.SubramanianA.PinchbackR.ThorvaldsdóttirH.TamayoP.MesirovJ. P. (2011)., 27. Oxford, England, 1739–1740. 10.1093/bioinformatics/btr260Molecular signatures database (MSigDB) 3.0BioinformaticsPubMed Abstract | CrossRef Full Text | Google Scholar
41
LuL.HuY.WangC.JiangF.WuC. (2021). Methylation and expression of the exercise-related TLR1 gene is associated with low grade glioma prognosis and outcome. Front. Mol. Biosci.8, 747933. 10.3389/fmolb.2021.747933PubMed Abstract | CrossRef Full Text | Google Scholar
42
MariathasanS.TurleyS. J.NicklesD.CastiglioniA.YuenK.WangY.et al (2018). TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature554, 544–548. 10.1038/nature25501PubMed Abstract | CrossRef Full Text | Google Scholar
43
MathewsM. S.BlickenstaffJ. W.ShihE. C.ZamoraG.VoV.SunC. H.et al (2012). Photochemical internalization of bleomycin for glioma treatment. J. Biomed. Opt.17, 058001. 10.1117/1.Jbo.17.5.058001PubMed Abstract | CrossRef Full Text | Google Scholar
44
MayakondaA.LinD. C.AssenovY.PlassC.KoefflerH. P. (2018). Maftools: Efficient and comprehensive analysis of somatic variants in cancer. Genome Res.28, 1747–1756. 10.1101/gr.239244.118PubMed Abstract | CrossRef Full Text | Google Scholar
45
OhueY.NishikawaH. (2019). Regulatory T (Treg) cells in cancer: Can Treg cells be a new therapeutic target?Cancer Sci.110, 2080–2089. 10.1111/cas.14069PubMed Abstract | CrossRef Full Text | Google Scholar
46
PengW. X.KoiralaP.MoY. Y. (2017). LncRNA-mediated regulation of cell signaling in cancer. Oncogene36, 5661–5667. 10.1038/onc.2017.184PubMed Abstract | CrossRef Full Text | Google Scholar
47
PengZ.LiuC.WuM. (2018). New insights into long noncoding RNAs and their roles in glioma. Mol. Cancer17, 61. 10.1186/s12943-018-0812-2PubMed Abstract | CrossRef Full Text | Google Scholar
48
QiY.LiuB.SunQ.XiongX.ChenQ. (2020). Immune checkpoint targeted therapy in glioma: Status and hopes. Front. Immunol.11, 578877. 10.3389/fimmu.2020.578877PubMed Abstract | CrossRef Full Text | Google Scholar
49
RitchieM. E.PhipsonB.WuD.HuY.LawC. W.ShiW.et al (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res.43, e47. 10.1093/nar/gkv007PubMed Abstract | CrossRef Full Text | Google Scholar
50
RizviN. A.HellmannM. D.SnyderA.KvistborgP.MakarovV.HavelJ. J.et al (2015). Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Sci. (New York, N.Y.)348, 124–128. 10.1126/science.aaa1348CrossRef Full Text | Google Scholar
51
RohW.ChenP. L.ReubenA.SpencerC. N.PrietoP. A.MillerJ. P.et al (2017). Integrated molecular analysis of tumor biopsies on sequential CTLA-4 and PD-1 blockade reveals markers of response and resistance. Sci. Transl. Med., eaah3560. 10.1126/scitranslmed.aah3560PubMed Abstract | CrossRef Full Text | Google Scholar
52
ShawE. G.WangM.CoonsS. W.BrachmanD. G.BucknerJ. C.StelzerK. J.et al (2012). Randomized trial of radiation therapy plus procarbazine, lomustine, and vincristine chemotherapy for supratentorial adult low-grade glioma: Initial results of RTOG 9802. J. Clin. Oncol.30, 3065–3070. 10.1200/jco.2011.35.8598PubMed Abstract | CrossRef Full Text | Google Scholar
53
SnyderA.MakarovV.MerghoubT.YuanJ.ZaretskyJ. M.DesrichardA.et al (2014). Genetic basis for clinical response to CTLA-4 blockade in melanoma. N. Engl. J. Med.371, 2189–2199. 10.1056/NEJMoa1406498PubMed Abstract | CrossRef Full Text | Google Scholar
54
SteinbrueckA.SedgwickA. C.BrewsterJ. T.2ndYanK. C.ShangY.KnollD. M.et al (2020). Transition metal chelators, pro-chelators, and ionophores as small molecule cancer chemotherapeutic agents. Chem. Soc. Rev.49, 3726–3747. 10.1039/c9cs00373hPubMed Abstract | CrossRef Full Text | Google Scholar
55
TangX.JiangF.WangX.XiaY.MaoY.ChenY. (2022). Identification of the ferroptosis-related long non-coding RNAs signature to improve the prognosis prediction in papillary renal cell carcinoma. Front. Surg.9, 741726. 10.3389/fsurg.2022.741726PubMed Abstract | CrossRef Full Text | Google Scholar
56
TopalianS. L.DrakeC. G.PardollD. M. (2015). Immune checkpoint blockade: A common denominator approach to cancer therapy. Cancer Cell.27, 450–461. 10.1016/j.ccell.2015.03.001PubMed Abstract | CrossRef Full Text | Google Scholar
57
TsvetkovP.CoyS.PetrovaB.DreishpoonM.VermaA.AbdusamadM.et al (2022). Copper induces cell death by targeting lipoylated TCA cycle proteins. Sci. (New York, N.Y.)375, 1254–1261. 10.1126/science.abf0529CrossRef Full Text | Google Scholar
58
TsvetkovP.DetappeA.CaiK.KeysH. R.BruneZ.YingW.et al (2019). Mitochondrial metabolism promotes adaptation to proteotoxic stress. Nat. Chem. Biol.15, 681–689. 10.1038/s41589-019-0291-9PubMed Abstract | CrossRef Full Text | Google Scholar
59
van den BentM. J. (2015). Chemotherapy for low-grade glioma: When, for whom, which regimen?Curr. Opin. Neurol.28, 633–938. 10.1097/wco.0000000000000257PubMed Abstract | CrossRef Full Text | Google Scholar
60
ViaccozA.LekoubouA.DucrayF. (2012). Chemotherapy in low-grade gliomas. Curr. Opin. Oncol.24, 694–701. 10.1097/CCO.0b013e328357f503PubMed Abstract | CrossRef Full Text | Google Scholar
61
WatersE.PucciP.HirstM.ChapmanS.WangY.CreaF.et al (2021). HAR1: An insight into lncRNA genetic evolution. Epigenomics13, 1831–1843. 10.2217/epi-2021-0069PubMed Abstract | CrossRef Full Text | Google Scholar
62
WellerM.WickW.AldapeK.BradaM.BergerM.PfisterS. M.et al (2015). Glioma. Nat. Rev. Dis. Prim.1, 15017. 10.1038/nrdp.2015.17PubMed Abstract | CrossRef Full Text | Google Scholar
63
WilkersonM. D.HayesD. N. (2010)., 26. Oxford, England, 1572–1573. 10.1093/bioinformatics/btq170ConsensusClusterPlus: A class discovery tool with confidence assessments and item trackingBioinformaticsPubMed Abstract | CrossRef Full Text | Google Scholar
64
XuL.DengC.PangB.ZhangX.LiuW.LiaoG.et al (2018). Tip: A web server for resolving tumor Immunophenotype profiling. Cancer Res.78, 6575–6580. 10.1158/0008-5472.Can-18-0689PubMed Abstract | CrossRef Full Text | Google Scholar
65
YamanakaR. (2009). Dendritic-cell- and peptide-based vaccination strategies for glioma. Neurosurg. Rev.32, 265–273. 10.1007/s10143-009-0189-1PubMed Abstract | CrossRef Full Text | Google Scholar
66
YoshiharaK.ShahmoradgoliM.MartínezE.VegesnaR.KimH.Torres-GarciaW.et al (2013). Inferring tumour purity and stromal and immune cell admixture from expression data. Nat. Commun.4, 2612. 10.1038/ncomms3612PubMed Abstract | CrossRef Full Text | Google Scholar
67
ZengD.YeZ.ShenR.YuG.WuJ.XiongY.et al (2021). Iobr: Multi-Omics immuno-oncology biological research to decode tumor microenvironment and signatures. Front. Immunol.12, 687975. 10.3389/fimmu.2021.687975PubMed Abstract | CrossRef Full Text | Google Scholar
68
ZhaiL.LadomerskyE.LenzenA.NguyenB.PatelR.LauingK. L.et al (2018). Ido1 in cancer: A gemini of immune checkpoints. Cell.. Mol. Immunol.15, 447–457. 10.1038/cmi.2017.143PubMed Abstract | CrossRef Full Text | Google Scholar
69
ZhangZ.KattanM. W. (2017). Drawing nomograms with R: Applications to categorical outcome and survival data. Ann. Transl. Med.5, 211. 10.21037/atm.2017.04.01PubMed Abstract | CrossRef Full Text | Google Scholar
70
ZhaoZ.LiuM.LongW.YuanJ.LiH.ZhangC.et al (2021). Knockdown lncRNA CRNDE enhances temozolomide chemosensitivity by regulating autophagy in glioblastoma. Cancer Cell. Int.21, 456. 10.1186/s12935-021-02153-xPubMed Abstract | CrossRef Full Text | Google Scholar
71
ZouW.WolchokJ. D.ChenL. (2016). PD-L1 (B7-H1) and PD-1 pathway blockade for cancer therapy: Mechanisms, response biomarkers, and combinations. Sci. Transl. Med.8, 328rv4. 10.1126/scitranslmed.aad7118PubMed Abstract | CrossRef Full Text | Google Scholar
Summary
Keywords
cuproptosis, lncRNAs, prognosis, low-grade glioma, immune checkpoint blockade, chemotherapy
Citation
Yan X, Wang N, Dong J, Wang F, Zhang J, Hu X, Zhao H, Gao X, Liu Z, Li Y and Hu S (2022) A cuproptosis-related lncRNAs signature for prognosis, chemotherapy, and immune checkpoint blockade therapy of low-grade glioma. Front. Mol. Biosci. 9:966843. doi: 10.3389/fmolb.2022.966843
Received
11 June 2022
Accepted
25 July 2022
Published
17 August 2022
Volume
9 - 2022
Edited by
Yu Xiao, Wuhan University, China
Updates
Copyright
© 2022 Yan, Wang, Dong, Wang, Zhang, Hu, Zhao, Gao, Liu, Li and Hu.
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: Zhihui Liu, ZhihuiLiu1995@163.com; Yongzhe Li, liyongzhe18@163.com; Shaoshan Hu, shaoshanhu421@163.com
This article was submitted to RNA Networks and Biology, a section of the journal Frontiers in Molecular Biosciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.