- 1State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, National Clinical Research Center for Infectious Diseases, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 2Division of Hepatobiliary and Pancreatic Surgery, Department of Surgery, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 3Department of Pathology, Hangzhou Tongchuang Medical Laboratory, Hangzhou, China
- 4Department of Infectious Diseases, ShuLan (Hangzhou) Hospital Affiliated to Zhejiang Shuren University, Shulan International Medical College, Hangzhou, China
Background: Low-grade glioma (LGG) is considered a fatal disease for young adults, with overall survival widely ranging from 1 to 15 years depending on histopathologic and molecular subtypes. As a novel type of programmed cell death, ferroptosis was reported to be involved in tumorigenesis and development, which has been intensively studied in recent years.
Methods: For the discovery cohort, data from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) were used to identify the differentially expressed and prognostic ferroptosis-related genes (FRGs). The least absolute shrinkage and selection operator (LASSO) and multivariate Cox were used to establish a prognostic signature with the above-selected FRGs. Then, the signature was developed and validated in TCGA and Chinese Glioma Genome Atlas (CGGA) databases. By combining clinicopathological features and the FRG signature, a nomogram was established to predict individuals’ one-, three-, and five-year survival probability, and its predictive performance was evaluated by Harrell’s concordance index (C-index) and calibration curves. Enrichment analysis was performed to explore the signaling pathways regulated by the signature.
Results: A novel risk signature contains seven FRGs that were constructed and were used to divide patients into two groups. Kaplan–Meier (K−M) survival curve and receiver-operating characteristic (ROC) curve analyses confirmed the prognostic performance of the risk model, followed by external validation based on data from the CGGA. The nomogram based on the risk signature and clinical traits was validated to perform well for predicting the survival rate of LGG. Finally, functional analysis revealed that the immune statuses were different between the two risk groups, which might help explain the underlying mechanisms of ferroptosis in LGG.
Conclusion: In conclusion, this study constructed a novel and robust seven-FRG signature and established a prognostic nomogram for LGG survival prediction.
Background
Lower-grade glioma (LGG) comprises of World Health Organization (WHO) grade II (diffuse lower-grade) and III (intermediate-grade) glioma tumors. Despite relatively more favorable clinical outcomes for LGG than the grade IV tumors, the survival of LGG patients widely varies from 1 to 15 years, and 70% of patients from the lower grade progress to high grade within 10 years (van den Bent, 2014). In the past decades, the histological classification has been the standard clinical approach for LGG, but it lacks the prediction of grades due to the wide range of clinical heterogeneities of LGG (Louis et al., 2016). Although some molecular markers, such as the presence or absence of mutation in isocitrate dehydrogenase (IDH) genes, have been integrated in the WHO classification for LGG, to identify novel prognostic biomarkers remains imperative.
Ferroptosis is a kind of programmed non-apoptotic cell death caused by an iron-dependent lethal accumulation of lipid peroxidation. As a promising therapeutic alternative, ferroptosis was reported to be closely related to multiple diseases such as neurodegenerative diseases (Do Van et al., 2016) and cancer (Shen et al., 2018). Aberrant expressions of ferroptosis-related genes (FRGs), such as tumor protein p53 (TP53) (Junttila and Evan, 2009), Fanconi anemia complementation group D2 (FANCD2) (Han et al., 2017), glutathione peroxidase 4 (GPX4) (Liu H. et al., 2018), heat shock protein beta 1 (HSPB1) (Arrigo and Gibert, 2012), and dipeptidyl-peptidase-4 (DPP4) (Enz et al., 2019), were reported to be correlated with tumor genesis and progression. A ferroptosis inducer, named erastin, could inhibit the tumor growth and enhance the sensitivity of chemotherapeutic drugs (Chen et al., 2015; Yu et al., 2015). Meanwhile, the induction of ferroptosis synergistically enhanced the antitumor activity of immune checkpoint inhibitors (ICIs), even in ICI-resistant tumors (Badgley et al., 2020). Pseudolaric acid B (PAB) triggered the ferroptosis and inhibited the viabilities of glioma cells (Wang Z. et al., 2018). These studies indicated the activation or inhibition of ferroptosis has potential clinical value, and better understanding of ferroptosis might provide prognostic value and therapeutic candidates for management of LGG.
In the present study, we analyzed the aberrant expression of FRGs by analyzing the LGG samples from The Cancer Genome Atlas (TCGA) and normal brain tissue from the Genotype-Tissue Expression (GTEx) database. The overall process is shown in Figure 1. To be specific, a prognostic signature was established based on seven FRGs’ expression levels which dissected the patients with LGG into two risk groups. Meanwhile, its predictive performance was elucidated by Kaplan–Meier (K–M) survival and receiver-operating characteristic (ROC) curves and further validated by two external cohorts consisting of LGG patients from the Chinese Glioma Genome Atlas (CGGA) database (Zhao et al., 2017). Furthermore, a survival nomogram was developed by integrating the risk signature with other four clinicopathological factors, and its predicted accuracy in LGG patients was assessed. Finally, enrichment analyses were performed to investigate the underlying mechanisms between the signature and tumor immunology. Overall, our data suggested that FRGs play pivotal roles in glioma progression and could be prognostic markers and therapeutic targets for glioma. The FRG signature constructed in this study might enhance the ability to predict the prognosis of LGG patients and provide potential explanations on the underlying mechanisms of ferroptosis in LGG progression.
Methods
Dataset and Source
MRNA expression data, DNA methylation 450K data, and corresponding clinical information of 529 LGG patients were obtained from TCGA data portal (https://gdc-portal.nci.nih.gov/). For differentially expressed gene analysis, gene expression data in normal tissue of Genotype-Tissue Expression (GTEx) were downloaded from the UCSC Xena browser (Goldman et al., 2014; Gentles et al., 2015). RNA-seq data and clinical information of other LGG sets were obtained from the Chinese Glioma Genome Atlas (CGGA) (http://www.cgga.org.cn/) (Bao et al., 2014; Wang et al., 2015; Zhao et al., 2017; Liu X. et al., 2018). Gene expression values of fragments per kilobase of exon per million reads mapped (FPKM) were normalized and subjected to the following analyses. A comprehensive list of ferroptosis-related genes (FRGs) was chosen and is provided in Supplementary Table S1.
Development and Evaluation of the Ferroptosis-Related Prognostic Signature
To identify differentially expressed FRGs (DE-FRGs) between TCGA and GTEx databases, the absolute value of the log2-transformed fold change (log2FC) > 2 and the false discovery rate (FDR) < 0.05 was considered to be statistically significant (Zhang et al., 2018). To investigate the mechanisms of dysregulation for FRGs in LGG, Pearson’s correlation between FRGs’ transcriptional expression levels and promotor methylation levels was calculated, and those with |R| > 0.3 and p < 0.05 were considered significantly correlated (Xu et al., 2019).
To screen genes for constructing a risk signature, univariate Cox regression models (p < 0.05) were performed to select genes that are associated with the overall survival of LGG patients in TCGA. The overlapped candidate prognostic DE-FRGs were visualized in a Venn diagram. A protein–protein interaction (PPI) network for the prognostic DE-FRGs was generated by the STRING database (version 11.0) (Szklarczyk et al., 2010).
Subsequently, the least absolute shrinkage and selection operator (LASSO) Cox regression was performed by using the “glmnet” R package to shrink the scope of gene screening. The multivariate Cox regression analysis was carried out to identify highly correlated genes and construct the survival risk score model. The risk score was calculated as follows:
where n, expi, and βi represent the number of genes, the expression level, and the coefficient of the gene i, respectively. According to the median value of the risk scores, the LGG patients were classified into low- and high-risk groups. Principal component analysis (PCA) was carried out to explore the distribution of different groups using the “prcomp” function in R. The Kaplan–Meier (K–M) curves were used to evaluate the statistical significance of the survival rates between different risk groups. The receiver-operating characteristic (ROC) curve analyses were executed to calculate the area under the curve (AUC) and to determine the predictive ability of the risk score signature. We calculated the false negative rate (FNR) of the survival prognosis by our FRG model.
In order to study whether the FRG-based prognostic signature could be used as an independent predictor for LGG patients, univariate and multivariate Cox regression analyses were conducted in TCGA and CGGA datasets. The risk score, age, gender, grade, and IDH status were used as covariates.
Construction and Evaluation of a Predictive Nomogram
To establish a clinically applicable method for predicting the prognosis of LGG patients, we formulated a prognostic nomogram to predict the survival probability at 1, 3, and 5 years based on TCGA dataset. Five independent prognostic parameters, including age, gender, grade, IDH1 status, and risk score, were used to establish the nomogram for predicting one-, three-, and five-year survival rates of LGG patients in the “rms,” “foreign,” and “survival” R packages. Bootstraps with 1,000 resamples were used for these activities. The performance of the prognostic nomogram was assessed by calculating Harrell’s concordance index (C-index) (Harrell et al., 1996), and calibration curves of the nomogram for one-, three-, and five-year overall survival were carried out to estimate the accuracy of actual observed rates with the predicted survival probability. During the external validation of the nomogram, the total points of each patient in the CGGA cohort were calculated according to the established nomogram, and then the C-index was calculated and calibration curves were derived.
Functional Annotation and Enrichment Analysis
The “clusterProfiler,” “org.Hs.eg.db,” “enrichplot,” and “ggplot2” packages in R were used to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation enrichment analyses based on the differentially expressed genes (|log2FC| ≥ 2, FDR <0.05) between the low- and high-risk groups. The top 30 GO terms and KEGG pathways were identified with a cutoff of p < 0.05 and FDR < 0.05.
The single-sample gene set enrichment analysis (ssGSEA) was conducted using the “GSVA” package in R to evaluate the enrichment score (ES) of immune parameters, including 15 immune cells and 13 immune-related pathways (Goldman et al., 2014). The annotated gene set file is provided in Supplementary Table S2. Pearson’s correlation coefficients were calculated to test the correlation between the ES and the risk score.
Statistical Analysis
All statistical analyses were performed using R (version 4.0.3). Differences between two sample groups were tested by the Wilcoxon test. Univariate and multivariate Cox regression analyses were performed to determine risk factors for low and high risk scores of the signature in LGG. The Venn diagrams, heatmaps, boxplots, forest plots, nomograms, and calibration plots were drawn using R language. The K–M curve used the log-rank test to evaluate the statistical difference of survival rates between different risk groups. The significance of survival time differences was calculated using the log-rank test with a threshold of p < 0.05.
Results
Characteristics of LGG Patients
The workflow of this study is shown in Figure 1. A total of 508 LGG patients from TCGA-LGG cohort and 171 LGG patients from the CGGA cohort were finally enrolled. The clinical characteristics of patients from two databases are summarized in Table 1.
Aberrant Expression of FRGs in LGG
32 FRG expressions were identified to have a significant difference between LGG samples from TCGA and normal samples from GTEx (Figure 2A). Fifteen genes were upregulated in LGG, whereas another seventeen genes were downregulated in tumor. Since methylation of promotor plays a critical role in regulating gene expression, the effect of promotor DNA methylation on the gene expression of DE-FRGs was evaluated. As shown in Figures 2A,B, a significant negative correlation was found between gene expression and DNA methylation in five genes, including FTH1, ACSL3, SLC1A5, HSPB1, and STEAP3, whereas significant positive correlations were found in CRYAB and ZEB1. These results demonstrated that promoter DNA methylation might regulate the expression of FRGs in LGG with complexly regulative patterns.
FIGURE 2. Dysregulation of ferroptosis-related genes (FRGs) and constructing a risk signature. (A) Heatmap of the differentially expressed FRGs between LGG and normal tissue. (B) Pearson’s correlation of FRGs between transcriptional expression and promoter methylation. (C) Venn diagram showing the overlapped gene number of differentially expressed FRGs and prognostic FRGs that correlated with OS. (D) Forest plot of hazard ratios demonstrating the prognostic value of overlapped FRGs. (E) PPI network indicating the interactions among the overlapped FRGs.
Construction of a Prognostic FRG Signature
A total of 37 FRGs were found for potent prognostic value (p < 0.05) through univariate Cox proportional hazard regression. 20 of them were differentially expressed between LGG and normal tissue (Figures 2C,D). The PPI network among these genes was constructed using the STRING online database, with parameters including a minimum required interaction score >0.4 (medium confidence). As shown in Figure 2E, TP53, SLC7A11, and HMOX1 were hub genes for ferroptosis in LGG (Figure 2E). LASSO regression analysis was performed on the 20 overlapped genes, and nine genes were retained according to the optimal lambda value (Supplementary Figure S1). By using multivariate Cox regression analysis, a prognostic model based on seven genes was developed. The prognostic risk score was imputed as follows: Risk score = (0.0086 × expression level of ACSL3) + (−0.0190 × expression level of CBS) + (0.0013 × expression level of CD44) + (−0.0027 × expression level of FADS2) + (0.0017 × expression level of HSPB1) + (0.0104 × expression level of PGD) + (0.0271 × expression level of STEAP3). The 509 patients were classified into high- and low-risk groups according to the median value as cutoff (Figure 3A and Supplementary Figures S2A,C). The PCA plot showed the patients in two risk groups were distributed in discrete directions (Figure 3B). Kaplan–Meier (K–M) plots indicated patients with high risk scores had a significantly worse OS probability (p = 9.178e-9) (Figure 3C). To verify the diagnostic competence of the ferroptosis-related risk signature, the AUC for OS prediction was calculated. The AUC of the ROC curve reached 0.855 at 1 year, 0.834 at 3 years, and 0.753 at 5 years, which indicated that the risk score presented excellent performance on prognostic prediction of TCGA-LGG (Figure 3D). The false negative rate (FNR) of prognosis for OS at 1 , 3 , and 5 years was 0.077, 0.118, and 0.169, respectively (Supplementary Table S3).
FIGURE 3. Prognostic analysis of the seven-FRG signature model in TCGA cohort. (A) Distribution and median value of the risk scores. (B) PCA plot of TCGA cohort. (C) Kaplan–Meier curves for the OS of patients in the high-risk group and low-risk group. (D) The AUC of time-dependent ROC curves verified the prognostic performance of the risk score.
Validation of the Molecular Signature in CGGA-LGG Dataset
The prognostic capability of the FRG risk signature was externally validated by the CGGA dataset. The patients from the CGGA cohort were also divided into low- or high-risk groups by the same cutoff value as that in TCGA cohort (Figure 4A, Supplementary Figures S2B,D, and Supplementary Figure S3). Consistent with the above findings, PCA found that the patients in different risk groups were distributed in two directions (Figure 4B). Additionally, the K–M survival curves showed a significantly worse survival rate with high-risk patients (Figure 4C and Supplementary Figure S3A). The AUC for one-, three-, and five-year survival was 0.752, 0.795, and 0.794, respectively (Figure 4D). These results confirmed that the seven-FRG prognostic model presented good performance in prognostic prediction of LGG.
FIGURE 4. External validation of the seven-FRG signature model in the CGGA cohort (n = 170). (A) Distribution of the risk scores. (B) PCA plot of the CGGA cohort. (C) Kaplan–Meier curves for the OS of patients in the high-risk group and low-risk group. (D) The AUC of time-dependent ROC curves verified the prognostic performance of the risk score.
Identification of the FRG Signature as an Independent Prognostic Factor
To evaluate the independent prognostic force of the signature, both the univariate and multivariate Cox proportion hazard regression models were applied. The risk score in univariate analysis was significantly correlated with OS in both TCGA and CGGA cohorts (HR = 3.814, 95% CI = 2.021–7.196, p < 0.001; HR = 4.250, 95% CI = 2.719–6.644, p < 0.001, respectively) (Figures 5A,B). The results from multivariate analysis found this signature was an independent prognostic factor for OS in TCGA-LGG (HR = 2.352, 95% CI = 1.179–4.690, p = 0.015) and CGGA-LGG (HR = 2.758, 95% CI = 1.572–4.840, p < 0.001) cohorts. These results confirmed the independent prognostic value of the seven-FRG–based prognostic signature for patients with LGG (Figures 5C,D).
FIGURE 5. Cox proportion hazard regression for the OS of LGG. (A, B) Univariate Cox regression analysis in TCGA derivation cohort and the CGGA validation cohort. (C, D) Multivariate Cox regression analysis in TCGA derivation cohort and the CGGA validation cohort.
Construction and Verification of Nomogram
Next, by integrating the above risk signature with clinicopathological factors, including age, gender, grade, and IDH status, a nomogram was cautiously constructed for prediction of one-, three-, and five-year survival rates of LGG. As shown in Figure 6A, each level of factors was assigned one point, and a total point was calculated by adding up the points in each individual. Subsequently, the C-indexes of the nomogram for predicting OS were 0.87 (95% CI, 0.83–0.91) and 0.68 (95% CI, 0.63–0.73) in TCGA and CGGA cohorts, respectively. Calibration curves also showed good agreement between predictive and observational values at the probabilities of one-, three-, and five-year survival in TCGA and CGGA cohorts (Figures 6B,C). These results revealed that the nomogram signified good accuracy in predicting the one-, three-, and five-year survival of patients with LGG.
FIGURE 6. The nomogram can predict the survival probability in LGG. (A) A nomogram of TCGA cohort used to predict the OS. (B, C) Calibration curves used to predict the one-, three-, and five-year survival in TCGA derivation cohort and the CGGA validation cohort. The x-axis and y-axis represent the predicted and actual survival rates of the nomogram, respectively.
Functional Annotation and Enrichment Analysis of the FRG Prognostic Signature
To further elucidate the biological functions and pathways associated with the risk signature, we performed GO and KEGG pathway enrichment analyses based on differentially expressed genes (DEGs) between the high- and low-risk groups. As shown in Figures 7A,C, GO analysis showed 13 clusters sets were enriched in both TCGA and CGGA cohorts. For biological processes (BPs), DEGs were mainly enriched into immune-related terms, including neutrophil degranulation and neutrophil activation involved in immune response. For cellular components (CCs), cell-substrate junction, focal adhesion, cell leading edge, and secretory granule lumen were found to be enriched. With respect to molecular function (MF), differentially expressed genes were mainly enriched in actin binding, integrin binding, cadherin binding, actin filament binding, and growth factor binding. For the KEGG pathway, 10 out of top 30 pathways were enriched in both cohorts, and most of these enriched pathways were related to infections of bacterial or viral (Figures 7B,D). Interestingly, the term of focal adhesion, which was also reported closely related to the tumor immune microenvironment (TME) (Murphy et al., 2020), was enriched in both GO and KEGG analyses. These results revealed that the risk signature was correlated with the infection and immune-related pathway.
FIGURE 7. Altered functional characteristics related to the seven-FRG signature. (A) Representative results of GO in TCGA cohort. (B) KEGG analyses in TCGA cohort. (C) Representative results of GO in the CGGA cohort. (D) KEGG analyses in the CGGA cohort. The blue rectangles indicate the overlapped terms between the two cohorts.
Furthermore, we quantified the enrichment scores (ESs) of diverse immune cell subtypes and related functions or pathways with ssGSEA to confirm the correlation between the risk signature and immune status. As shown by the boxplot in Figure 8A, in TCGA-LGG cohort, 25 out of 28 immune-related terms were significantly higher in the high-risk group (p < 0.05, Figure 8A), except DCs (dendritic cells), neutrophils, and NK cells (natural killer cells). Consistently, twenty-four immune-related terms were validated by the ICGC cohort, except Tfh cells (T follicular helper cells), while NK cells were significantly higher in the high-risk group (p < 0.01, Figure 8C).
FIGURE 8. Correlation of the risk score with infiltrative immune cells. (A) Comparison of the ssGSEA scores between different risk groups in TCGA cohort. (B) Correlations between ssGSEA scores and risk scores in TCGA cohort. (C) Comparison of the ssGSEA scores between different risk groups in the CGGA cohort. (D) Correlations between ssGSEA scores and risk scores in the CGGA cohort. P values are shown as follows: ns, not significant; *p < 0.05; **p < 0.01; ***p < 0.001.
We further investigated the association between expressions of the seven-gene signature and the immune status ES in detail by calculating Pearson’s correlation. The results showed that most of the seven genes were positively correlated with these immune cells and functions, but CBS (cystathionine beta-synthase) and FADS2 (fatty acid desaturase 2) were negatively correlated with these immune terms in both cohorts (p < 0.05, Figures 8B,D). However, inconsistent correlations of acyl-CoA synthetase long-chain family member 3 (ACSL3) with immune characters were observed between TCGA and CGGA cohorts (Figures 8B,D).
Discussion
LGG, generally acknowledged as benign tumor and universally prone to progressing into high-grade glioma, is the major cause of mortality for young adults (Hottinger et al., 2016; Picca et al., 2018). In recent years, some molecular characteristics were found to be useful for LGG classification, prognosis prediction, and treatment guidance (Takano et al., 2016). However, demand for a novel prognostic biomarker of LGG is still urgent.
As a newly defined form of programmed cell death, ferroptosis and its role in oncology have attracted increasing levels of attention from researchers. It is well known that ferroptosis is involved in cancer pathology as an inhibitor or promotor. But its effect is exerted depending on the tumor type, stage, and microenvironment (Daher et al., 2019; Dai et al., 2020). It was also demonstrated that FRGs play roles in glioma cell death and drug resistance (Wang Z. et al., 2018; Hu et al., 2020). A 19-FRG signature also showed high correlation with clinical outcomes of glioma (Liu et al., 2020). These studies indicated the potential role of FRGs in glioma.
In this study, we identified that there were 29 differentially expressed FRGs in LGG, and seven of these genes were significantly correlated with promotor methylation. We used univariate Cox analysis to screen the relations between 29 DE-FRGs and prognosis, and 20 genes were found to be related to the survival of LGG. Followed by LASSO and multivariate Cox regression analysis, we developed an efficient risk model which consisted of seven FRGs and divided LGG patients into low- and high-risk groups by the risk score. The predictive performance of this risk model was validated by K–M and ROC curve analyses in both TCGA and CGGA cohorts.
Some clinical characters, including age (≤40 years vs. >40 years), tumor grade (II vs. III), and IDH status (wild-type vs. mutation), have been well-practiced in prognostic analysis for LGG patients (Nabors et al., 2020). Multivariate Cox analysis showed that the risk group was an independent prognostic indicator as well. Furthermore, we used the risk score and other clinical characteristics to construct a nomogram prognosis model. The nomogram performed well in prediction of OS probability, which was validated by time-dependent calibration plots and showed satisfactory performance in both TCGA and CGGA cohorts. Therefore, what we established is a clinically applicable prediction model to assist individualized survival prediction of patients with LGG.
Among the seven prognostic FRGs, only HSPB1 and CBS were downregulated in LGG. Previous reports identified that they are both negative regulators for ferroptosis (Sun et al., 2015; Wang L. et al., 2018). Higher HSPB1 was found to contribute to glioma development (Ye et al., 2016), while CBS expression, consistent with our results, was reduced in glioma progression (Takano et al., 2014). ACSL3, a key regulator of lipid metabolism, is required for inducing ferroptosis resistance (Magtanong et al., 2019). The previous studies indicate that ACSL3 overexpression was associated with worse clinical outcomes in patients with NSCLC and melanoma (Chen et al., 2016; Padanad et al., 2016); in contrast, high ACSL3 expression predicted a better prognosis in ovarian cancer (Chen et al., 2016). CD44 expression suppressed ferroptosis in cancer cells and might exert its function as a glioma promoter by increasing tumor cell invasion and proliferation, which therefore could be a promising therapeutic target (Mooney et al., 2016; Liu et al., 2019). FADS2 is upregulated in glioblastoma, and its deficiency inhibited glioblastoma cell proliferation, which offers a potential novel therapeutic target (Affronti and Wellen, 2019). PGD has been reported to overexpress in a number of cancer types and to be associated with poor prognosis. In addition, the upregulated expression of PGD also performs a key role in the development of radiochemotherapeutic resistance in cancer (Yang et al., 2018; Sarfraz et al., 2020). Six-transmembrane epithelial antigen of prostate 3 (STEAP3) is a ferrireductase, which is vital for cellular iron uptake and homeostasis (Ohgami et al., 2005). A recent investigation suggested that STEAP3 might exert its function as an oncogenic mediator in glioma progression (Han et al., 2018). In our results, five of these seven FRGs (ACSL3, CD44, HSPB1, PGD, and STEAP3) were associated with poor prognosis, while the remaining two genes (CBS and FADS2) were the opposite.
Based on the DEGs between the two risk groups, we performed GO and KEGG annotation and surprisingly found that the term of focal adhesion was enriched in both GO and KEGG analyses. Focal adhesion–related terms, such as “integrin and growth factor binding,” were also enriched, which were considered closely related to cancer cell metastasis and tumor microenvironment (Murphy et al., 2020). A previous report that the expression of focal adhesion kinase was affected by high intracellular iron oxide nanoparticle concentrations (Soenen et al., 2010) may partially explain the high correlation between focal adhesion and ferroptosis catalyzed by iron. Additionally, some immune-related biological processes and pathways such as neutrophil-related terms in GO and infection pathways in KEGG were also enriched. The previous study found the ferroptosis of immune cells during infection is advantageous for infectious agents (Matsushita et al., 2015), but the regulatory mechanism between ferroptosis and infectious pathways in cancer is still unknown. Emerging evidence suggests that ferroptosis is associated with tumor immunity (Stockwell and Jiang, 2019; Liu et al., 2020; Tang et al., 2020). Consistently, our results on ssGSEA also indicated that the ferroptosis-related signature might positively regulate immune signaling pathways. High- and low- risk score groups defined by our risk model exhibited distinct immune landscapes. Specifically, we found that the ESs of universal immune features were significantly higher in the high-risk group. Pearson’s correlation test further revealed the expressions of CD44, HSPB1, PGD, and STEAP3 were positively correlated with the general 38 immune ESs, while the expression of CBS and FADS2 had an inverse correlation. The coefficients of CBS and FADS2 in the risk model were negative, and they, together with other genes’ expression, contributed to the positive correlations of risk scores with immune ESs. In the CGGA cohort, it should be pointed out that ACSL3 expression did not significantly correlate with immune feature ESs except the B cell. Recent studies confirmed that CD8+ T cells are required for glioma growth, and patients with higher abundance of CD8+ T cells in the TME were more likely to benefit from immune checkpoint inhibitor treatment (Tumeh et al., 2014; Kane et al., 2020). Taken together, a positive correlation between the risk score of our FRG model and immune characteristics indicated that the risk score could be a potential predictive biomarker for immune checkpoint blockade and worthy of further investigation. Understanding the role of immune checkpoint inhibitors in combination with ferroptosis inducers for gliomas is urgent to improve the efficacy of immunotherapy.
However, there were some limitations in our study. Firstly, the data came from two datasets with limited sample size, and additional robust analyses (Wu and Ma, 2015; Ren et al., 2019) in large-scale independent cohorts are needed in the future. Secondly, the prognostic factors in the nomogram cannot include some important prognostic factors in surgical resection, radiation treatment, or chemotherapy. The collection of detailed clinical information from samples might help to overcome these shortcomings and assess the diagnostic efficiency of the nomogram in future studies. Finally, the biological mechanisms and association with immunotherapeutic effect of the prognostic signature are still unknown. Further study is required to extend the understanding of ferroptosis and promote novel therapeutic strategies for LGG patients.
Conclusion
In summary, we developed a seven-FRG prognostic signature and nomogram which have good performance in the prediction of survival of LGG patients from TCGA and CGGA datasets. Our preliminary observation on correlations between FRGs and tumor immune features also suggested that further investigation of underlying mechanisms between ferroptotic regulation and tumor immunity in LGG is warranted.
Data Availability Statement
The data of lower-grade glioma in this study were downloaded from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). The data of the normal brain tissues were downloaded from Genotype-Tissue Expression Database (GTEx).
Author Contributions
JZ and XZ contributed to literature search and data collection. JZ and ZL conducted the data analysis. JZ and HG wrote the original draft and reviewed and edited the manuscript. LL supervised the work. All authors approved the final manuscript.
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.
Acknowledgments
The authors would like to thank The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Chinese Glioma Genome Atlas (CGGA) databases for the availability of the data.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.753680/full#supplementary-material
Abbreviations
AUC, area under the ROC curve; CGGA, Chinese Glioma Genome Atlas; FDR, false discovery rate; FRGs, ferroptosis-related genes; GO, Gene Ontology; GSEA, gene set enrichment analysis; GTEx, Genotype-Tissue Expression; C-index, Harrell’s concordance index; KEGG, Kyoto Encyclopedia of Genes and Genomes; K–M, Kaplan–Meier; LASSO, least absolute shrinkage and selection operator; LGG, lower-grade glioma; OS, overall survival; ROC, receiver-operating characteristic; TCGA, The Cancer Genome Atlas
References
Affronti, H. C., and Wellen, K. E. (2019). Epigenetic Control of Fatty-Acid Metabolism Sustains Glioma Stem Cells. Cancer Discov. 9 (9), 1161–1163. doi:10.1158/2159-8290.CD-19-0733
Arrigo, A.-P., and Gibert, B. (2012). HspB1 Dynamic Phospho-Oligomeric Structure Dependent Interactome as Cancer Therapeutic Target. Cmm 12 (9), 1151–1163. doi:10.2174/156652412803306693
Badgley, M. A., Kremer, D. M., Maurer, H. C., DelGiorno, K. E., Lee, H.-J., Purohit, V., et al. (2020). Cysteine Depletion Induces Pancreatic Tumor Ferroptosis in Mice. Science 368 (6486), 85–89. doi:10.1126/science.aaw9872
Bao, Z.-S., Chen, H.-M., Yang, M.-Y., Zhang, C.-B., Yu, K., Ye, W.-L., et al. (2014). RNA-seq of 272 Gliomas Revealed a Novel, recurrentPTPRZ1-METfusion Transcript in Secondary Glioblastomas. Genome Res. 24 (11), 1765–1773. doi:10.1101/gr.165126.113
Chen, L., Li, X., Liu, L., Yu, B., Xue, Y., and Liu, Y. (2015). Erastin Sensitizes Glioblastoma Cells to Temozolomide by Restraining xCT and Cystathionine-γ-Lyase Function. Oncol. Rep. 33 (3), 1465–1474. doi:10.3892/or.2015.3712
Chen, W.-C., Wang, C.-Y., Hung, Y.-H., Weng, T.-Y., Yen, M.-C., and Lai, M.-D. (2016). Systematic Analysis of Gene Expression Alterations and Clinical Outcomes for Long-Chain Acyl-Coenzyme A Synthetase Family in Cancer. PloS one 11 (5), e0155660. doi:10.1371/journal.pone.0155660
Daher, B., Parks, S. K., Durivault, J., Cormerais, Y., Baidarjad, H., Tambutte, E., et al. (2019). Genetic Ablation of the Cystine Transporter xCT in PDAC Cells Inhibits mTORC1, Growth, Survival, and Tumor Formation via Nutrient and Oxidative Stresses. Cancer Res. 79 (15), 3877–3890. doi:10.1158/0008-5472.CAN-18-3855
Dai, E., Han, L., Liu, J., Xie, Y., Kroemer, G., Klionsky, D. J., et al. (2020). Autophagy-dependent Ferroptosis Drives Tumor-Associated Macrophage Polarization via Release and Uptake of Oncogenic KRAS Protein. Autophagy 16 (11), 2069–2083. doi:10.1080/15548627.2020.1714209
Do Van, B., Gouel, F., Jonneaux, A., Timmerman, K., Gelé, P., Pétrault, M., et al. (2016). Ferroptosis, a Newly Characterized Form of Cell Death in Parkinson's Disease that Is Regulated by PKC. Neurobiol. Dis. 94, 169–178. doi:10.1016/j.nbd.2016.05.011
Enz, N., Vliegen, G., De Meester, I., and Jungraithmayr, W. (2019). CD26/DPP4 - a Potential Biomarker and Target for Cancer Therapy. Pharmacol. Ther. 198, 135–159. doi:10.1016/j.pharmthera.2019.02.015
Gentles, A. J., Newman, A. M., Liu, C. L., Bratman, S. V., Feng, W., Kim, D., et al. (2015). The Prognostic Landscape of Genes and Infiltrating Immune Cells across Human Cancers. Nat. Med. 21 (8), 938–945. doi:10.1038/nm.3909
Goldman, M., Craft, B., Swatloski, T., Cline, M., Morozova, O., Diekhans, M., et al. (2014). The UCSC Cancer Genomics Browser: Update 2015. Nucleic Acids Res. 43 (D1), D812–D817. doi:10.1093/nar/gku1073
Han, B., Shen, Y., Zhang, P., Jayabal, P., Che, R., Zhang, J., et al. (2017). Overlooked FANCD2 Variant Encodes a Promising, Portent Tumor Suppressor, and Alternative Polyadenylation Contributes to its Expression. Oncotarget 8 (14), 22490–22500. doi:10.18632/oncotarget.14989
Han, M., Xu, R., Wang, S., Yang, N., Ni, S., Zhang, Q., et al. (2018). Six-Transmembrane Epithelial Antigen of Prostate 3 Predicts Poor Prognosis and Promotes Glioblastoma Growth and Invasion. Neoplasia 20 (6), 543–554. doi:10.1016/j.neo.2018.04.002
Harrell, F. E., Lee, K. L., and Mark, D. B. (1996). Multivariable Prognostic Models: Issues in Developing Models, Evaluating Assumptions and Adequacy, and Measuring and Reducing Errors. Statist. Med. 15 (4), 361–387. doi:10.1002/(sici)1097-0258(19960229)15:4<361::aid-sim168>3.0.co;2-4
Hottinger, A. F., Hegi, M. E., and Baumert, B. G. (2016). Current Management of Low-Grade Gliomas. Curr. Opin. Neurol. 29 (6), 782–788. doi:10.1097/wco.0000000000000390
Hu, Z., Mi, Y., Qian, H., Guo, N., Yan, A., Zhang, Y., et al. (2020). A Potential Mechanism of Temozolomide Resistance in Glioma-Ferroptosis. Front. Oncol. 10 (897). doi:10.3389/fonc.2020.00897
Junttila, M. R., and Evan, G. I. (2009). p53 - a Jack of All Trades but Master of None. Nat. Rev. Cancer 9 (11), 821–829. doi:10.1038/nrc2728
Kane, J. R., Zhao, J., Tsujiuchi, T., Laffleur, B., Arrieta, V. A., Mahajan, A., et al. (2020). CD8+ T-Cell-Mediated Immunoediting Influences Genomic Evolution and Immune Evasion in Murine Gliomas. Clin. Cancer Res. 26 (16), 4390–4401. doi:10.1158/1078-0432.CCR-19-3104
Liu, H.-J., Hu, H.-M., Li, G.-Z., Zhang, Y., Wu, F., Liu, X., et al. (2020). Ferroptosis-Related Gene Signature Predicts Glioma Cell Death and Glioma Patient Progression. Front. Cel Dev. Biol. 8, 538. doi:10.3389/fcell.2020.00538
Liu, H., Schreiber, S. L., and Stockwell, B. R. (2018a). Targeting Dependency on the GPX4 Lipid Peroxide Repair Pathway for Cancer Therapy. Biochemistry 57 (14), 2059–2060. doi:10.1021/acs.biochem.8b00307
Liu, T., Jiang, L., Tavana, O., and Gu, W. (2019). The Deubiquitylase OTUB1 Mediates Ferroptosis via Stabilization of SLC7A11. Cancer Res. 79 (8), 1913–1924. doi:10.1158/0008-5472.CAN-18-3037
Liu, X., Li, Y., Qian, Z., Sun, Z., Xu, K., Wang, K., et al. (2018b). A Radiomic Signature as a Non-invasive Predictor of Progression-free Survival in Patients with Lower-Grade Gliomas. NeuroImage: Clin. 20, 1070–1077. doi:10.1016/j.nicl.2018.10.014
Louis, D. N., Perry, A., Reifenberger, G., von Deimling, A., Figarella-Branger, D., Cavenee, W. K., et al. (2016). The 2016 World Health Organization Classification of Tumors of the Central Nervous System: a Summary. Acta Neuropathol. 131 (6), 803–820. doi:10.1007/s00401-016-1545-1
Magtanong, L., Ko, P.-J., To, M., Cao, J. Y., Forcina, G. C., Tarangelo, A., et al. (2019). Exogenous Monounsaturated Fatty Acids Promote a Ferroptosis-Resistant Cell State. Cel Chem. Biol. 26 (3), 420–432.e429. doi:10.1016/j.chembiol.2018.11.016
Matsushita, M., Freigang, S., Schneider, C., Conrad, M., Bornkamm, G. W., and Kopf, M. (2015). T Cell Lipid Peroxidation Induces Ferroptosis and Prevents Immunity to Infection. J. Exp. Med. 212 (4), 555–568. doi:10.1084/jem.20140857
Mooney, K. L., Choy, W., Sidhu, S., Pelargos, P., Bui, T. T., Voth, B., et al. (2016). The Role of CD44 in Glioblastoma Multiforme. J. Clin. Neurosci. 34, 1–5. doi:10.1016/j.jocn.2016.05.012
Murphy, J. M., Rodriguez, Y. A. R., Jeong, K., Ahn, E.-Y. E., and Lim, S.-T. S. (2020). Targeting Focal Adhesion Kinase in Cancer Cells and the Tumor Microenvironment. Exp. Mol. Med. 52 (6), 877–886. doi:10.1038/s12276-020-0447-4
Nabors, L. B., Portnow, J., Ahluwalia, M., Baehring, J., Brem, H., Brem, S., et al. (2020). Central Nervous System Cancers, Version 3.2020, NCCN Clinical Practice Guidelines in Oncology. NCCN Clin. Pract. Guidel. Oncol. 18 (11), 1537–1570. doi:10.6004/jnccn.2020.0052
Ohgami, R. S., Campagna, D. R., Greer, E. L., Antiochos, B., McDonald, A., Chen, J., et al. (2005). Identification of a Ferrireductase Required for Efficient Transferrin-dependent Iron Uptake in Erythroid Cells. Nat. Genet. 37 (11), 1264–1269. doi:10.1038/ng1658
Padanad, M. S., Konstantinidou, G., Venkateswaran, N., Melegari, M., Rindhe, S., Mitsche, M., et al. (2016). Fatty Acid Oxidation Mediated by Acyl-CoA Synthetase Long Chain 3 Is Required for Mutant KRAS Lung Tumorigenesis. Cel Rep. 16 (6), 1614–1628. doi:10.1016/j.celrep.2016.07.009
Picca, A., Berzero, G., and Sanson, M. (2018). Current Therapeutic Approaches to Diffuse Grade II and III Gliomas. Ther. Adv. Neurol. Disord. 11, 175628561775203. doi:10.1177/1756285617752039
Ren, J., Du, Y., Li, S., Ma, S., Jiang, Y., and Wu, C. (2019). Robust Network-Based Regularization and Variable Selection for High-Dimensional Genomic Data in Cancer Prognosis. Genet. Epidemiol. 43 (3), 276–291. doi:10.1002/gepi.22194
Sarfraz, I., Rasul, A., Hussain, G., Shah, M. A., Zahoor, A. F., Asrar, M., et al. (2020). 6‐Phosphogluconate Dehydrogenase Fuels Multiple Aspects of Cancer Cells: From Cancer Initiation to Metastasis and Chemoresistance. BioFactors 46 (4), 550–562. doi:10.1002/biof.1624
Shen, Z., Liu, T., Li, Y., Lau, J., Yang, Z., Fan, W., et al. (2018). Fenton-Reaction-Acceleratable Magnetic Nanoparticles for Ferroptosis Therapy of Orthotopic Brain Tumors. ACS Nano 12 (11), 11355–11365. doi:10.1021/acsnano.8b06201
Soenen, S. J. H., Nuytten, N., De Meyer, S. F., De Smedt, S. C., and De Cuyper, M. (2010). High Intracellular Iron Oxide Nanoparticle Concentrations Affect Cellular Cytoskeleton and Focal Adhesion Kinase-Mediated Signaling. Small 6 (7), 832–842. doi:10.1002/smll.200902084
Stockwell, B. R., and Jiang, X. (2019). A Physiological Function for Ferroptosis in Tumor Suppression by the Immune System. Cel Metab. 30 (1), 14–15. doi:10.1016/j.cmet.2019.06.012
Sun, X., Ou, Z., Xie, M., Kang, R., Fan, Y., Niu, X., et al. (2015). HSPB1 as a Novel Regulator of Ferroptotic Cancer Cell Death. Oncogene 34 (45), 5617–5625. doi:10.1038/onc.2015.32
Szklarczyk, D., Franceschini, A., Kuhn, M., Simonovic, M., Roth, A., Minguez, P., et al. (2010). The String Database in 2011: Functional Interaction Networks of Proteins, Globally Integrated and Scored. Nucleic Acids Res. 39 (Suppl. l_1), D561–D568. doi:10.1093/nar/gkq973
Takano, N., Sarfraz, Y., Gilkes, D. M., Chaturvedi, P., Xiang, L., Suematsu, M., et al. (2014). Decreased Expression of Cystathionine β-Synthase Promotes Glioma Tumorigenesis. Mol. Cancer Res. 12 (10), 1398–1406. doi:10.1158/1541-7786.MCR-14-0184
Takano, S., Ishikawa, E., Sakamoto, N., Matsuda, M., Akutsu, H., Noguchi, M., et al. (2016). Immunohistochemistry on IDH 1/2, Atrx, P53 and Ki-67 Substitute Molecular Genetic Testing and Predict Patient Prognosis in Grade III Adult Diffuse Gliomas. Brain Tumor Pathol. 33 (2), 107–116. doi:10.1007/s10014-016-0260-x
Tang, R., Xu, J., Zhang, B., Liu, J., Liang, C., Hua, J., et al. (2020). Ferroptosis, Necroptosis, and Pyroptosis in Anticancer Immunity. J. Hematol. Oncol. 13 (1), 110. doi:10.1186/s13045-020-00946-7
Tumeh, P. C., Harview, C. L., Yearley, J. H., Shintaku, I. P., Taylor, E. J. M., Robert, L., et al. (2014). PD-1 Blockade Induces Responses by Inhibiting Adaptive Immune Resistance. Nature 515 (7528), 568–571. doi:10.1038/nature13954
van den Bent, M. J. (2014). Practice Changing Mature Results of RTOG Study 9802: Another Positive PCV Trial Makes Adjuvant Chemotherapy Part of Standard of Care in Low-Grade Glioma. Neuro-Oncology 16 (12), 1570–1574. doi:10.1093/neuonc/nou297
Wang, L., Cai, H., Hu, Y., Liu, F., Huang, S., Zhou, Y., et al. (2018a). A Pharmacological Probe Identifies Cystathionine β-synthase as a New Negative Regulator for Ferroptosis. Cell Death Dis 9 (10), 1005. doi:10.1038/s41419-018-1063-2
Wang, Y., Qian, T., You, G., Peng, X., Chen, C., You, Y., et al. (2015). Localizing Seizure-Susceptible Brain Regions Associated with Low-Grade Gliomas Using Voxel-Based Lesion-Symptom Mapping. Neuro-Oncology 17 (2), 282–288. doi:10.1093/neuonc/nou130
Wang, Z., Ding, Y., Wang, X., Lu, S., Wang, C., He, C., et al. (2018b). Pseudolaric Acid B Triggers Ferroptosis in Glioma Cells via Activation of Nox4 and Inhibition of xCT. Cancer Lett. 428, 21–33. doi:10.1016/j.canlet.2018.04.021
Wu, C., and Ma, S. (2015). A Selective Review of Robust Variable Selection with Applications in Bioinformatics. Brief Bioinform 16 (5), 873–883. doi:10.1093/bib/bbu046
Xu, W., Xu, M., Wang, L., Zhou, W., Xiang, R., Shi, Y., et al. (2019). Integrative Analysis of DNA Methylation and Gene Expression Identified Cervical Cancer-specific Diagnostic Biomarkers. Sig Transduct Target. Ther. 4 (1), 55. doi:10.1038/s41392-019-0081-6
Yang, X., Peng, X., and Huang, J. (2018). Inhibiting 6-phosphogluconate Dehydrogenase Selectively Targets Breast Cancer through AMPK Activation. Clin. Transl Oncol. 20 (9), 1145–1152. doi:10.1007/s12094-018-1833-4
Ye, H., Huang, H., Cao, F., Chen, M., Zheng, X., and Zhan, R. (2016). HSPB1 Enhances SIRT2-Mediated G6PD Activation and Promotes Glioma Cell Proliferation. PloS one 11 (10), e0164285. doi:10.1371/journal.pone.0164285
Yu, Y., Xie, Y., Cao, L., Yang, L., Yang, M., Lotze, M. T., et al. (2015). The Ferroptosis Inducer Erastin Enhances Sensitivity of Acute Myeloid Leukemia Cells to Chemotherapeutic Agents. Mol. Cell Oncol. 2 (4), e1054549. doi:10.1080/23723556.2015.1054549
Zhang, L., He, L., Lugano, R., Roodakker, K., Bergqvist, M., Smits, A., et al. (2018). IDH Mutation Status Is Associated with Distinct Vascular Gene Expression Signatures in Lower-Grade Gliomas. Neuro-Oncology 20 (11), 1505–1516. doi:10.1093/neuonc/noy088
Keywords: ferroptosis, lower-grade glioma, prognostic signature, TCGA, CGGA, immune statuses
Citation: Zhao J, Liu Z, Zheng X, Gao H and Li L (2021) Prognostic Model and Nomogram Construction Based on a Novel Ferroptosis-Related Gene Signature in Lower-Grade Glioma. Front. Genet. 12:753680. doi: 10.3389/fgene.2021.753680
Received: 11 August 2021; Accepted: 04 October 2021;
Published: 08 November 2021.
Edited by:
Jijun Tang, University of South Carolina, United StatesCopyright © 2021 Zhao, Liu, Zheng, Gao and Li. 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: Lanjuan Li, bGpsaUB6anUuZWR1LmNu; Hainv Gao, Z2FvaGFpbnZAMTYzLmNvbQ==