Yikang He
Yangfan Ye2,3
Wei Tian
Huaide Qiu- 1Department of Rehabilitation, Zhongda Hospital, School of Medicine, Southeast University, Nanjing, China
- 2Department of Neurosurgery, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China
- 3The First Clinical Medical College, Nanjing Medical University, Nanjing, China
- 4Center of Rehabilitation Medicine, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China
Objective: To establish a lncRNA panel related to ferroptosis, tumor progression, and microenvironment for prognostic estimation in patients with glioma.
Methods: LncRNAs associated with tumor progression and microenvironment were screened via the weighted gene co-expression network analysis (WGCNA). Overlapped lncRNAs highlighted in WGCNA, related to ferroptosis, and incorporated in Chinese Glioma Genome Atlas (CGGA) were identified as hub lncRNAs. With expression profiles of the hub lncRNA, we conducted the least absolute shrinkage and selection operator (LASSO) regression and built a ferroptosis-related lncRNA signature to separate glioma patients with distinct survival outcomes. The lncRNA signature was validated in TCGA, the CGGA_693, and CGGA_325 cohorts using Kaplan-Meier survival analysis and ROC curves. The ferroptosis-related lncRNA panel was validated with 15 glioma samples using quantitative real-time PCR (qRT-PCR). Multivariate Cox regression was performed, and a nomogram was mapped and validated. Immune infiltration correlated to the signature was explored using TIMER and CIBERSORT algorithms.
Results: The present study identified 30 hub lncRNAs related to ferroptosis, tumor progression, and microenvironment. With the 30 hub lncRNAs, we developed a lncRNA signature with distinct stratification of survival chance in patients with glioma in two independent cohorts (HRs>1, p < 0.05). The lncRNA signature revealed a panel of 14 lncRNAs, i.e., APCDD1L-AS1, H19, LINC00205, LINC00346, LINC00475, LINC00484, LINC00601, LINC00664, LINC00886, LUCAT1, MIR155HG, NEAT1, PVT1, and SNHG18. These lncRNA expressions were validated in clinical specimens using qRT-PCR. Robust predictive accuracies of the signature were present across different datasets at multiple timepoints. With univariate and multivariate regressions, we demonstrated that the risk score based on the lncRNA signature is an independent prognostic indicator after clinical factors were adjusted. A nomogram was constructed with these prognostic factors, and it has demonstrated decent classification and accuracy. Additionally, the signature-based classification was observed to be correlated with multiple clinical characteristics and molecular subtypes. Further, extensive immune cells were upregulated in the high-risk group, such as CD8+ T cell, neutrophil, macrophage, and myeloid dendritic cell, indicating increased immune infiltrations.
Conclusion: We established a novel ferroptosis-related lncRNA signature that could effectively stratify the prognosis of glioma patients with adequate predictive performance.
Introduction
Glioma is the most common intracranial malignancy, accounting for 80% of brain malignancies and 30% of brain tumors (de Robles et al., 2015). Adult diffused gliomas were classified into lower-grade glioma (LGG) and glioblastoma multiforme (GBM), as indicated in the 2016 World Health Organization (WHO) classification (Louis et al., 2016; Wen and Huse, 2017). While GBM presents an inferior prognosis with a median survival of 14.6 months (Lara-Velazquez et al., 2017), most LGG progresses to GBM, leading to an unfavorable overall survival rate (Ostrom et al., 2014; Claus et al., 2015; Lim et al., 2018; Kirby and Finnerty, 2020). Interestingly, patients with LGG could decease within months post-surgery while those with GBM lived for years in some cases, which indicates that the current classification of gliomas does not reflect the distinct survival outcomes (Chen et al., 2017). To evaluate prognosis and develop therapeutic strategies, molecular biomarkers including O6-methylguanine-DNA methyltransferase (MGMT) methylation, codeletion of chromosome 1 and 19 (1p/19q), mutations in isocitrate dehydrogenase (IDH) are under active investigation in patients with glioma (Molinaro et al., 2019). Although these molecular biomarkers are implicated in the regulation of cancer cell proliferation and death (Ceccarelli et al., 2016), it remains a challenge to inform individual prognosis and treatment. Therefore, novel biomarkers are warranted for prognostic stratification and therapeutic strategies.
Ferroptosis, a novel form of programmed cell death, is characterized by the aberrant accumulation of lipid peroxide related to iron metabolism (Dixon et al., 2012). Various types of cancer cells were observed to be sensitive to ferroptosis (Eling et al., 2015; Belavgeni et al., 2019; Carbone and Melino, 2019), which may alter the process of cancer development and therapies (Gan, 2019; Liang et al., 2019; Stockwell and Jiang, 2019). Additionally, genes with regulatory effects on ferroptosis, such as P53 (Junttila and Evan, 2009), DPP4 (Enz et al., 2019), and HSPB1 (Arrigo and Gibert, 2012), are associated with the development and progression of cancers. Besides, ferroptosis-related genes have shown potentials to be prognostic biomarkers in glioma (Liu et al., 2020; Zhuo et al., 2020), indicating the significance of ferroptosis in glioma.
Long non-coding RNAs (lncRNAs) are a type of transcripts with more than 200 nucleotides which does not code for proteins (Clark et al., 2012). lncRNAs are the major component of human transcriptome with a proportion up to 80% (Hong et al., 2020). Owing to the ubiquitous regulatory effect of lncRNAs, it is involved in tumorigenesis in numerous cancers. Moreover, dysregulation of lncRNAs has been reported to be associated with the pathogenicity of gliomas (Bhan et al., 2017). Further, there has been emerging reports of lncRNA signatures in predicting survival of glioma patients (Wang et al., 2018; Luan et al., 2019; Xia et al., 2020). A ferroptosis-related lncRNAs signature was observed to be associated with prognosis, tumor microenvironment, and response to radiotherapy in glioma (Zheng et al., 2021). However, the previous study used median risk scores for each dataset, leading to half-by-half classifications and potential false-positive results. In machine learning models, a consistent cutoff value for different datasets enhances generalizability, and thereby increasing applicability in real world.
In the present study, we sought to establish a ferroptosis-related lncRNA signature with transcriptomic and clinical data from the Cancer Genome Atlas (TCGA) and two cohorts from the Chinese Glioma Genome Atlas (CGGA). We analyzed the expression profiles of ferroptosis-related lncRNAs and fitted key lncRNAs into a prediction model, which was rigorously validated in three datasets with the same cutoff value. Our data indicate that the ferroptosis-related lncRNA signature is ready to serve as a novel prognostic panel, and therapies targeting these lncRNAs might yield effective results.
Materials and Methods
Data Collection
RNA sequencing (RNA-seq) data and corresponding clinical information of 698 glioma samples were taken from TCGA (https://portal.gdc.cancer.gov), while the test datasets with 693 glioma samples (CGGA_693) as well as with 325 samples (CGGA_325) was retrieved from the CGGA (http://www.cgga.org.cn). Then ferroptosis-related genes were retrieved from previous literature (Liu et al., 2020).
Coexpression Network Identified Clinically Relevant Long Non-Coding RNAs
The lncRNAs associated with tumor progression and microenvironment were screened using the weighted gene co-expression network analysis (WGCNA) package (Langfelder and Horvath, 2008). Samples and lncRNAs were filtered with statistical tests, while outliers were detected with hierarchical clustering. Then we selected the soft threshold and engendered the hierarchical clustering tree. Module eigengene (ME) representing the overall expression level of the module was calculated, and the Pearson correlation coefficients between lncRNA modules and clinical traits were calculated to identify modules of interest. Modules with correlation coefficients >0.5 with regarding to tumor grade and ESTIMATE scores were identified, and lncRNAs in these modules were identified as clinically relevant lncRNAs.
Subsequently, we calculated the correlation coefficients between each gene expression in the specific module and clinical trait, and the value of coefficients were defined as gene significance (GS). Module memberships (MMs) were defined as the coefficients between gene expression and MEs. Correlations between GS and MMs were analyzed to verify the significance of lncRNAs in the specific module.
Identification of Hub Long Non-Coding RNAs and Development of a Prognostic Signature
Based on the list of clinically relevant lncRNAs, we identified specific ferroptosis-related lncRNAs whose expression was correlated to ferroptosis-related genes (Pearson Cor>0.5 and p < 0.05). Overlapped lncRNAs indexed in the ferroptosis-related lncRNA list, highlighted in WGCNA, and incorporated in CGGA were identified as hub lncRNAs. With expression profiles of the hub lncRNA, we conducted the least absolute shrinkage and selection operator (LASSO) regression analysis (Meier et al., 2008) and selected key lncRNAs in TCGA dataset. Subsequently, we calculated the individualized risk score with coefficients and expression levels of the key lncRNAs, which stands for a ferroptosis-related prognostic signature that separates the high-risk glioma patients from the low-risk ones. The formulae for risk score were as follows.
(Xi represents the identified hub lncRNA expression, and βi represents the corresponding coefficient.)
TCGA-glioma patients were divided into high-risk and low-risk groups based on the median risk score as the cutoff. Then the clinical relevance was validated using Kaplan-Meier survival analysis between risk groups. The prognostic signature was presented as a risk plot that includes the distribution of risk score and classification of survival status for patients in different groups. The time-dependent receiver operating characteristic (ROC) analysis was performed (Heagerty et al., 2000), and the area under the curve (AUC) was calculated to evaluate the prognostic accuracy in the TCGA-glioma patients at 1-year, 3-years, and 5-years post-operation. AUC>0.7 was considered adequate accuracy. The prognostic signature with the same risk score formula and cutoff value was then validated in the CGGA-glioma patients. It was then validated in CGGA using Kaplan-Meier survival analysis, ROC curves, as well as a risk plot.
The Independently Prognostic Value of the Signature
The univariate and multivariate Cox regressions were conducted to determine whether the prognostic signature was independent of clinical characteristics as well as to identify other independent prognostic factors. To determine whether the signature was independent of additional treatment, we also conducted subgroup survival analysis with variables of chemotherapy and radiotherapy. A nomogram was constructed with these prognostic factors using the rms R package, and it was validated using calibration plot (Qiu et al., 2020). Additionally, the correlation between the prognostic signature and clinical characteristics were explored by chi-square test using TCGA.
Quantitative Real-Time PCR Analysis of Long Non-Coding RNAs in Glioma Tissues
A total of 15 glioma specimens (nine LGG vs six high-grade gliomas [HGG]) were obtained from the Department of Neurosurgery, First Affiliated Hospital of Nanjing Medical University. All specimens were stored in liquid nitrogen, and all patients provided informed consent. The collecting of specimens was approved by the Ethics Committee of First Affiliated Hospital of Nanjing Medical University (No: 2020-SRFA-167).
Total RNA extraction was performed by standard procedures. quantitative real-time PCR (qRT-PCR) was conducted using Oligo primers (Tsingke Biotechnology, for indicated lncRNAs with AceQ qPCR SYBR Master Mix (Vazyme, Q141-02). In all assays, GAPDH served as normalization control. Data were analyzed using the −ΔΔCt method with each test performed in triplicate.
Immune Infiltration Correlated to the Long Non-Coding RNAs Signature
To explore the immune infiltration in the glioma samples stratified by the signature, immune cell abundances were calculated using TIMER (Li et al., 2020) and CIBERSORT (Newman et al., 2015). Between-group comparisons were conducted using Wilcoxon test, and immune cell types with p < 0.05 were kept for visualization.
Statistical Analysis
Cox regression analysis and Kaplan-Meier method were performed for survival analysis, where Cox p-values and log-rank p-values were calculated. Pearson correlation was applied to determine significant correlations between continuous variables, while chi-square test was adopted to delineate correlations between categorical variables. PCR data were analyzed with independent T tests between groups. Multiple T tests were corrected using the Two-stage linear step-up procedure of Benjamini, Krieger and Yekutieli, with q = 1% being statistically significant (Benjamini et al., 2006). Data were analyzed using R software 4.1.1 and GraphPad Prism 8.3.0, where p-value<0.05 was regarded as statistically significant.
Results
Identification of Hub Long Non-Coding RNAs and Enrichment Analysis
The flow chart of the present study is presented in Figure 1. Using the WGCNA package, we performed the hierarchical clustering and detected one outlier (Figure 2A). With the cutoff value of 0.95 in Scale-free R2, the soft-thresholding power was set to be five (Figure 2B). The clustering dendrograms identified 14 modules (Figure 2C), out of which the turquoise module was significantly correlated with tumor grade as well as tumor microenvironment (Cor>0.5, Figure 2D). To be specific, features including tumor grade, ESTIMATE combined score, stromal score, and immune score were highlighted. Scatterplots of GS for the aforementioned features versus MM in the turquoise module exhibit a marked correlation, suggesting the functional significance of these lncRNAs in the turquoise module (Figures 2E–H). Therefore, a total of 1223 clinically relevant lncRNAs embedded in the turquoise module were selected.
FIGURE 1. Work flow of the present study.
FIGURE 2. WGCNA using lncRNA expressions (A) Sample clustering to detect outliers (B) Soft threshold selection (C) Cluster dendrogram to distinguish modules (D) Correlation matrix between modules and clinical traits (E–H) Scatter plot with GS and MM in the turquoise module with regard to tumor grade (E), ESTIMATE combined score (F), ESTIMATE stromal score (G), and ESTIMATE immune score (H).
Development and Validation of the Prognostic Signature
Out of the 1223 lncRNAs in the turquoise module, 303 lncRNAs were related to ferroptosis. After intersecting the 303 ferroptosis-related lncRNAs with the gene list of CGGA, we retrieved 30 hub lncRNAs (Figure 3A). The correlation between the expressions of 30 hub lncRNAs (columns) and ferroptosis-related genes (rows) was presented in Figure 3B. The forest plot showing prognostic value of the 30 hub lncRNAs was presented in Figure 3C. With the 30 hub lncRNAs, we selected 14 key lncRNAs with Lasso regression using TCGA dataset (Figures 4A,B). As shown in Figure 4C, the 14 lncRNAs were as follows (with absolute value of model coefficients increasing): APCDD1L-AS1, H19, LINC00205, LINC00346, LINC00475, LINC00484, LINC00601, LINC00664, LINC00886, LUCAT1, MIR155HG, NEAT1, PVT1, and SNHG18. Except for LINC00205 exhibits a negative coefficient on risk score, which indicates its protective effects on prognosis, the other lncRNAs were observed to be risk factors with positive coefficients (Figure 4C). Subsequently, we calculated individualized risk scores with coefficient-weighted expression levels of 14 lncRNAs after extracting the coefficients. TCGA-glioma patients were divided into high-score and low-score groups based on the median risk score of 3.04 as the cutoff. Significant difference was shown between these two groups in survival analysis (Figure 4D, HR = 3.1, 95% CI = 2.7–3.6, p < 0.001). ROC curve analysis of the signature in TCGA-glioma patients achieved AUC of 0.853, 0.856, and 0.825 at 1-, 3- and 5-years (Figure 4E). The risk plot presented a clear separation of survival status between risk groups with red dots being ceased cases and blue ones alive (Figure 4F). Likewise, we calculated individualized risk scores in CCGA_693 and CGGA_325 glioma patients and divided them into high-risk and low-risk groups with the same cutoff value (3.04). The patients in high-risk group had worse prognosis compared to those in low-risk group (Figures 4G,J). And ROC curve analysis of the signature in CGGA_693 achieved AUCs of 0.719, 0.783, and 0.782 at 1-, 3- and 5-years respectively, indicating a moderately high sensitivity and specificity of the signature in predicting patients’ survival (Figure 4H). Similar result was observed in ROC curve analysis of the signature based on the CGGA_325 (Figure 4K). Further, we revealed the distribution patterns of risk scores and survival status, and found that the survival time of patients with high-risk scores was markedly lower than that of patients with low-risk scores (Figures 4I, 4L).
FIGURE 3. Identification of hub lncRNAs (A) Venn Gram for identification of hub lncRNAs (B) Correlation matrix of hub lncRNAs and genes implicated in ferroptosis (C) Forest plot showing prognostic value of the 30 hub lncRNAs.
FIGURE 4. Construction and validation of the ferroptosis-related lncRNA signature (A–B) Visualization of LASSO regression (C) Coefficients of the 14 lncRNAs fitted in the signature (D) Kaplan-Meier curve presenting survival of high-risk and low-risk groups in TCGA (E) Time-dependent ROC for the signature in TCGA (F) Discriminative plot of the signature using TCGA (G, J) Kaplan-Meier curve presenting survival of high-risk and low-risk groups in CGGA_693 or CGGA_325 respectively (H, K) Time-dependent ROC for the signature in CGGA_693 or CGGA_325 respectively (I, L) Discriminative plot of the signature using CGGA_693 or CGGA_325 separately.
Assessment of the Independence of the Prognostic Signature
We next sought to determine whether the risk score based on the lncRNA signature is an independent prognostic indicator of OS in glioma patients. The risk score and several clinicopathological variables such as gender, age, grade, IDH mutation status, 1p19q-codeletion status, MGMT promoter methylation status, radiotherapy and chemotherapy were included in univariate and multivariate Cox regression analyses. Univariate Cox regression analysis showed that risk score, age, grade, radiotherapy, IDH mutation status, 1p19q-codeletion status and MGMT promoter methylation status were significantly correlated with the OS of glioma patients, whereas gender and chemotherapy did not show any correlation (Figure 5A, left panel). Four variables were further identified as independent prognostic factors of OS by multivariate analysis, including risk score (HR = 2.257, 95% CI = 1.437–3.544, p < 0.001), age (HR = 1.040, 95% CI = 1.022–1.059, p < 0.001), grade (HR = 2.130, 95% CI = 1.365–3.324, p < 0.001) and radiotherapy (HR = 0.507, 95% CI = 0.279–0.920, p = 0.026) (Figure 5A, right panel). Further, subgroup survival analysis showed the signature predicts overall survival in groups with chemotherapy (Figure 5B), without chemotherapy (Figure 5C), with radiotherapy (Figure 5D), and without radiotherapy (Figure 5E).
FIGURE 5. Evaluation of the prognostic signature and establishment of a nomogram (A) Univariate and Multivariate Cox regression analysis of the ferroptosis-related lncRNA signature and other clinicopathological features in gliomas (B–E) Subgroup survival analysis using groups with chemotherapy (B), without chemotherapy (C), with radiotherapy (D), without radiotherapy (E) (F) The nomogram incorporating independent prognostic factors (G) ROC curves for the nomogram-based classification at multiple timepoints (H) Calibration plots for 1-, 3-, and 5-years survival predictions.
With the aforementioned independent prognostic factors, we established a nomogram for predicting individual OS. As presented in Figure 5F, the clinicians can predict survival outcome for individuals with known age, tumor grade, and risk score in an intuitive manner. The classification accuracy was 0.941, 0.964, and 0.897 for 1-, 3-, and 5-years survival prediction (Figure 5G). Further, the prediction accuracy was also evaluated with calibration plots, which showed that the nomogram-based prediction was highly consistent with the actual survival at 1-, 3-, and 5-years post-operation (Figure 5H).
Clinical Correlations of the Prognostic Classifier
As shown in Figure 6A, the prognostic classifier was significantly correlated with multiple clinical features, including age, grade, chemotherapy, radiotherapy and MGMT promoter methylation status, 1p19q-codeletion status, and IDH mutation status (p < 0.001). Notably, the correlation between the lncRNA signature and grade was prominent. Further, we analyzed the expressions of 14 lncRNA signatures between high-grade and lower-grade samples, and the between-group results were distinct (p < 0.001) (Figure 6B).
FIGURE 6. Clinical correlation of the prognostic signature (A) Heatmap presenting correlations between the signature and clinical features (B) Boxplot showing lncRNA expressions between high-grade and lower-grade samples. ***p < 0.001, ****p < 0.0001 (C) Validation for the lncRNA signature in glioma samples using qRT-PCR. ****q < 0.0001. (D) Valcano plot showing expressions of the 14 lncRNAs across LGG and HGG. Red dots represent upregulated lncRNAs in HGG and the green dot represents downregulated lncRNA.
Validations for Long Non-Coding RNAs Expressions With Quantitative Real-Time PCR
Out of the 14 lncRNAs, high expressions of PVT1 and SNHG18 in glioblastoma samples as opposed to LGG was validated in the previous study (Zheng et al., 2021). Likewise, lncRNA MIR155HG was reported in high-grade glioma as compared to LGG (He et al., 2021). As such, we hereby validated the other 11 lncRNAs in clinical specimens. As shown in Figure 6C, we successfully validated 11 lncRNAs, i.e., APCDD1L-AS1, H19, LINC00205, LINC00346, LINC00475, LINC00484, LINC00601, LINC00664, LINC00886, LUCAT1, NEAT1 (q < 0.0001). Apart from LINC00205, all the other lncRNAs were upregulated in HGG (Figures 6C,D).
Immune Infiltration in the High-Risk Group
Using the TIMER algorithm, we identified four types of immune cells upregulated in the high-risk group, i.e., CD8+ T cell, neutrophil, macrophage, and myeloid dendritic cell (Figure 7). With CIBERSORT algorithm, more extensive immune cell components were identified, such as B cell memory, T cell CD4+ memory (resting and activated), macrophage (M0, M1, and M2), and mast cell activated. These results indicate there was increased immune infiltration in the high-risk group defined by the ferroptosis-related lncRNA signature (Figure 7).
FIGURE 7. Heatmap presenting immune infiltration correlated to the lncRNA signature.
Discussion
Due to the inherent heterogeneity of glioma, the current WHO classification of gliomas does not sufficiently characterize the distinct survival outcomes. Even with well-established molecular markers, the prognosis stratification and treatment of glioma patients remain challenges. Ferroptosis is involved in the development and progression of cancers and its pivotal role in glioma has been recently unveiled (Liu et al., 2020). In this context, the present study aimed to establish a novel ferroptosis-related lncRNA signature for glioma.
Our study identified 30 ferroptosis-related lncRNAs in association with tumor grade and microenvironment via WGCNA. With the 30 hub lncRNAs, we developed a lncRNA signature for prognosis prediction. The ferroptosis-related lncRNA signature was clinically relevant with distinct stratification of survival chance in patients with glioma in two independent cohorts. Robust predictive accuracies of the signature were present across different datasets at multiple timepoints. With univariate and multivariate regressions, we demonstrated that the risk score based on the lncRNA signature is an independent prognostic indicator after clinical factors were adjusted. The nomogram incorporating age, tumor grade, radiotherapy and risk score was observed to be accurate in both classification and prediction performance. What’s more, the signature-based stratification was observed to be correlated with extensive clinical characteristics, including tumor grade, MGMT promoter methylation status, 1p19q-codeletion status, and IDH mutation status. Further, increased immune infiltration was observed in the high-risk group as defined by the signature, which contributed to detrimental survival outcomes.
The ferroptosis-related lncRNA signature highlighted a total of 14 lncRNAs, i.e., APCDD1L-AS1, H19, LINC00205, LINC00346, LINC00475, LINC00484, LINC00601, LINC00664, LINC00886, LUCAT1, MIR155HG, NEAT1, PVT1, and SNHG18. Out of the 14 lncRNAs, eight of them were previously reported. H19 was reported to be upregulated in high-grade glioma and facilitated the proliferation of diffuse intrinsic pontine glioma (Roig-Carles et al., 2021). LINC00346 promotes cell migration, proliferation, and apoptosis by targeting ROCK1(Chen et al., 2020) and miR-128-3p/SZRD1 axis (Geng et al., 2020) in glioma. Overexpression of LINC00475 was observed in hypoxic gliomas and silencing LINC00475 results in suppression of tumor proliferation, migration, as well as invasion (Yu et al., 2020a). By regulating the miR-141-3p/YAP1 axis, LINC00475 facilitates tumor progression of glioma (Yu et al., 2020b). LUCAT1 is an oncogenic molecule in glioma, and its knockdown induced inhibition of cell viability and invasion by regulating miR-375 in glioma (Gao et al., 2018). LncRNA miR155HG contributes to tumor growth and progression in glioblastoma (Wu et al., 2019), and its prognostic value was also identified in previous bioinformatic analysis (Zheng et al., 2021). NEAT1 demonstrated to be a contributor to glioma cell migration, invasion, and tumor progression (Chen et al., 2018; Zhou et al., 2018). PVT1 promotes tumorigenesis and cancer progression of glioma via regulation of miR-128-3p/GREM1 Axis and BMP signaling pathway (Fu et al., 2018). Upregulation of SNHG18 promotes resistance to radiotherapy in glioma by repressing Semaphorin 5A (Zheng et al., 2016). High expressions of PVT1 and SNHG18 in glioblastoma samples as opposed to LGG was validated in the previous study (Zheng et al., 2021). Notably, the role of lncRNAs APCDD1L-AS1, LINC00205, LINC00484, LINC00601, LINC00664, and LINC00886 was first reported and validated in the present study.
Zheng et al. (Zheng et al., 2021) established a prognostic ferroptosis-Related lncRNAs signature correlated with immune landscape and radiotherapy response in Glioma. However, they applied median risk scores for three cohorts with a half-by-half classification in each cohort, thereby leading to potential false-positive results. Notably, there was no other report of ferroptosis-related lncRNA signature in gliomas. By contrast, the present study reports a ferroptosis-related lncRNA signature, which exhibits a distinct separation of survival rates in glioma patients with the same cutoff value across different cohorts. Further, expression levels of the ferroptosis-related lncRNAs signature were also validated between high-grade and lower-grade samples using qPCR. It could be translated to clinical settings as a test panel for informing individualized prognosis. Therapies targeting these lncRNAs could hold promise for enhanced efficacy. However, a lack of functional validation is considered as a major limitation in the present study. The functional implications of the lncRNAs require further validation with empirical data.
Conclusion
The present study established a ferroptosis-related lncRNA signature that could effectively stratify the prognosis of glioma patients, with superior predictive performance to all clinical index. These findings may pave the way for developing novel biomarkers for prognosis and treatments in gliomas.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article.
Ethics Statement
All specimens were stored in liquid nitrogen, and all patients provided informed consent. The collecting of specimens was approved by the Ethics Committee of First Affiliated Hospital of Nanjing Medical University (No: 2020-SRFA-167).
Author Contributions
HQ conceived the project and designed the study. YH and HQ analyzed and interpreted the data. WT and YY collected clinical samples and conducted lab investigation. HQ and WT contributed to the revision of this manuscript. All authors have read and approved the final manuscript.
Funding
This study was funded by the Postgraduate Research; Practice Innovation Program of Jiangsu Province (SJCX20_0478 & SJCX20_0485).
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.
References
Arrigo, A.-P., and Gibert, B. (2012). HspB1 Dynamic Phospho-Oligomeric Structure Dependent Interactome as Cancer Therapeutic Target. Curr. Mol. Med. 12 (9), 1151–1163. doi:10.2174/156652412803306693
Belavgeni, A., Bornstein, S. R., von Mässenhausen, A., Tonnus, W., Stumpf, J., Meyer, C., et al. (2019). Exquisite Sensitivity of Adrenocortical Carcinomas to Induction of Ferroptosis. Proc. Natl. Acad. Sci. USA 116 (44), 22269–22274. doi:10.1073/pnas.1912700116
Benjamini, Y., Krieger, A. M., and Yekutieli, D. (2006). Adaptive Linear Step-Up Procedures that Control the False Discovery Rate. Biometrika 93 (3), 491–507. doi:10.1093/biomet/93.3.491
Bhan, A., Soleimani, M., and Mandal, S. S. (2017). Long Noncoding RNA and Cancer: a New Paradigm. Cancer Res. 77 (15), 3965–3981. doi:10.1158/0008-5472.can-16-2634
Carbone, M., and Melino, G. (2019). Stearoyl CoA Desaturase Regulates Ferroptosis in Ovarian Cancer Offering New Therapeutic Perspectives. Cancer Res. 79 (20), 5149–5150. doi:10.1158/0008-5472.can-19-2453
Ceccarelli, M., Barthel, F. P., Malta, T. M., Sabedot, T. S., Salama, S. R., Murray, B. A., et al. (2016). Molecular Profiling Reveals Biologically Discrete Subsets and Pathways of Progression in Diffuse Glioma. Cell 164 (3), 550–563. doi:10.1016/j.cell.2015.12.028
Chen, Q., Cai, J., Wang, Q., Wang, Y., Liu, M., Yang, J., et al. (2018). Long Noncoding RNA NEAT1, Regulated by the EGFR Pathway, Contributes to Glioblastoma Progression through the WNT/β-Catenin Pathway by Scaffolding EZH2. Clin. Cancer Res. 24 (3), 684–695. doi:10.1158/1078-0432.ccr-17-0605
Chen, R., Smith-Cohn, M., Cohen, A. L., and Colman, H. (2017). Glioma Subclassifications and Their Clinical Significance. Neurotherapeutics 14 (2), 284–297. doi:10.1007/s13311-017-0519-x
Chen, X., Li, D., Chen, L., Hao, B., Gao, Y., Li, L., et al. (2020). Long Noncoding RNA LINC00346 Promotes Glioma Cell Migration, Invasion and Proliferation by Up‐regulating ROCK1. J. Cel. Mol. Med. 24 (22), 13010–13019. doi:10.1111/jcmm.15899
Clark, M. B., Johnston, R. L., Inostroza-Ponta, M., Fox, A. H., Fortini, E., Moscato, P., et al. (2012). Genome-wide Analysis of Long Noncoding RNA Stability. Genome Res. 22 (5), 885–898. doi:10.1101/gr.131037.111
Claus, E. B., Walsh, K. M., Wiencke, J. K., Molinaro, A. M., Wiemels, J. L., Schildkraut, J. M., et al. (2015). Survival and Low-Grade Glioma: the Emergence of Genetic Information. Neurosurg Focus 38 (1), E6. doi:10.3171/2014.10.focus12367
de Robles, P., Fiest, K. M., Frolkis, A. D., Pringsheim, T., Atta, C., St. Germaine-Smith, C., et al. (2015). The Worldwide Incidence and Prevalence of Primary Brain Tumors: a Systematic Review and Meta-Analysis. Neuro Oncol. 17 (6), 776–783. doi:10.1093/neuonc/nou283
Dixon, S. J., Lemberg, K. M., Lamprecht, M. R., Skouta, R., Zaitsev, E. M., Gleason, C. E., et al. (2012). Ferroptosis: an Iron-dependent Form of Nonapoptotic Cell Death. Cell 149 (5), 1060–1072. doi:10.1016/j.cell.2012.03.042
Eling, N., Reuter, L., Hazin, J., Hamacher-Brady, A., and Brady, N. R. (2015). Identification of Artesunate as a Specific Activator of Ferroptosis in Pancreatic Cancer Cells. Oncoscience 2 (5), 517–532. doi:10.18632/oncoscience.160
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
Fu, C., Li, D., Zhang, X., Liu, N., Chi, G., and Jin, X. (2018). LncRNA PVT1 Facilitates Tumorigenesis and Progression of Glioma via Regulation of MiR-128-3p/GREM1 Axis and BMP Signaling Pathway. Neurotherapeutics 15 (4), 1139–1157. doi:10.1007/s13311-018-0649-9
Gan, B. (2019). DUBbing Ferroptosis in Cancer Cells. Cancer Res. 79 (8), 1749–1750. doi:10.1158/0008-5472.can-19-0487
Gao, Y.-S., Liu, X.-Z., Zhang, Y.-G., Liu, X.-J., and Li, L.-Z. (2018). Knockdown of Long Noncoding RNA LUCAT1 Inhibits Cell Viability and Invasion by Regulating miR-375 in Glioma. Oncol. Res. 26 (2), 307–313. doi:10.3727/096504017x15088061795756
Geng, Y. B., Pan, C. C., Xu, C., Zuo, P. C., Wang, Y., Li, X. O., et al. (2020). Long Non-coding RNA LINC00346 Regulates Proliferation and Apoptosis by Targeting miR-128-3p/SZRD1 axis in Glioma. Eur. Rev. Med. Pharmacol. Sci. 24 (18), 9581–9590. doi:10.26355/eurrev_202009_23046
He, X., Sheng, J., Yu, W., Wang, K., Zhu, S., and Liu, Q. (2021). LncRNA MIR155HG Promotes Temozolomide Resistance by Activating the Wnt/β-Catenin Pathway via Binding to PTBP1 in Glioma. Cell Mol Neurobiol 41 (6), 1271–1284. doi:10.1007/s10571-020-00898-z
Heagerty, P. J., Lumley, T., and Pepe, M. S. (2000). Time-Dependent ROC Curves for Censored Survival Data and a Diagnostic Marker. Biometrics 56 (2), 337–344. doi:10.1111/j.0006-341x.2000.00337.x
Hong, W., Liang, L., Gu, Y., Qi, Z., Qiu, H., Yang, X., et al. (2020). Immune-related lncRNA to Construct Novel Signature and Predict the Immune Landscape of Human Hepatocellular Carcinoma. Mol. Ther. - Nucleic Acids 22, 937–947. doi:10.1016/j.omtn.2020.10.002
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
Kirby, A. J., and Finnerty, G. T. (2020). New Strategies for Managing Adult Gliomas. J. Neurol. 268, 3666–3674. doi:10.1007/s00415-020-09884-3
Langfelder, P., and Horvath, S. (2008). WGCNA: an R Package for Weighted Correlation Network Analysis. BMC Bioinformatics 9, 559. doi:10.1186/1471-2105-9-559
Lara-Velazquez, M., Al-Kharboosh, R., Jeanneret, S., Vazquez-Ramos, C., Mahato, D., Tavanaiepour, D., et al. (2017). Advances in Brain Tumor Surgery for Glioblastoma in Adults. Brain Sci. 7 (12), 166. doi:10.3390/brainsci7120166
Li, T., Fu, J., Zeng, Z., Cohen, D., Li, J., Chen, Q., et al. (2020). TIMER2.0 for Analysis of Tumor-Infiltrating Immune Cells. Nucleic Acids Res. 48 (W1), W509–w514. doi:10.1093/nar/gkaa407
Liang, C., Zhang, X., Yang, M., and Dong, X. (2019). Recent Progress in Ferroptosis Inducers for Cancer Therapy. Adv. Mater. 31 (51), 1904197. doi:10.1002/adma.201904197
Lim, M., Xia, Y., Bettegowda, C., and Weller, M. (2018). Current State of Immunotherapy for Glioblastoma. Nat. Rev. Clin. Oncol. 15 (7), 422–442. doi:10.1038/s41571-018-0003-5
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
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
Luan, F., Chen, W., Chen, M., Yan, J., Chen, H., Yu, H., et al. (2019). An Autophagy‐related Long Non‐coding RNA Signature for Glioma. FEBS open bio 9 (4), 653–667. doi:10.1002/2211-5463.12601
Meier, L., Van De Geer, S., and Bühlmann, P. (2008). The Group Lasso for Logistic Regression. J. R. Stat. Soc. Ser. B (Stat Methodol) 70 (1), 53–71. doi:10.1111/j.1467-9868.2007.00627.x
Molinaro, A. M., Taylor, J. W., Wiencke, J. K., and Wrensch, M. R. (2019). Genetic and Molecular Epidemiology of Adult Diffuse Glioma. Nat. Rev. Neurol. 15 (7), 405–417. doi:10.1038/s41582-019-0220-2
Newman, A. M., Liu, C. L., Green, M. R., Gentles, A. J., Feng, W., Xu, Y., et al. (2015). Robust Enumeration of Cell Subsets from Tissue Expression Profiles. Nat. Methods 12 (5), 453–457. doi:10.1038/nmeth.3337
Ostrom, Q. T., Bauchet, L., Davis, F. G., Deltour, I., Fisher, J. L., Langer, C. E., et al. (2014). The Epidemiology of Glioma in Adults: a "state of the Science" Review. Neuro-Oncology 16 (7), 896–913. doi:10.1093/neuonc/nou087
Qiu, H., Hu, X., He, C., Yu, B., Li, Y., and Li, J. (2020). Identification and Validation of an Individualized Prognostic Signature of Bladder Cancer Based on Seven Immune Related Genes. Front. Genet. 11, 12. doi:10.3389/fgene.2020.00012
Roig-Carles, D., Jackson, H., Loveson, K. F., Mackay, A., Mather, R. L., Waters, E., et al. (2021). The Long Non-coding RNA H19 Drives the Proliferation of Diffuse Intrinsic Pontine Glioma with H3K27 Mutation. Int. J. Mol. Sci. 22 (17). doi:10.3390/ijms22179165
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
Wang, W., Zhao, Z., Yang, F., Wang, H., Wu, F., Liang, T., et al. (2018). An Immune-Related lncRNA Signature for Patients with Anaplastic Gliomas. J. Neurooncol. 136 (2), 263–271. doi:10.1007/s11060-017-2667-6
Wen, P. Y., and Huse, J. T. (2017). 2016 World Health Organization Classification of Central Nervous System Tumors. Continuum (Minneap Minn) 23 (6), 1531–1547. doi:10.1212/CON.0000000000000536
Wu, W., Yu, T., Wu, Y., Tian, W., Zhang, J., and Wang, Y. (2019). The miR155HG/miR-185/ANXA2 Loop Contributes to Glioblastoma Growth and Progression. J. Exp. Clin. Cancer Res. 38 (1), 133. doi:10.1186/s13046-019-1132-0
Xia, P., Li, Q., Wu, G., and Huang, Y. (2020). An Immune-Related lncRNA Signature to Predict Survival in Glioma Patients. Cell Mol. Neurobiol. 41 (2), 365–375. doi:10.1007/s10571-020-00857-8
Yu, L., Gui, S., Liu, Y., Qiu, X., Qiu, B., Zhang, X., et al. (2020). Long Intergenic Non-protein Coding RNA 00475 Silencing Acts as a Tumor Suppressor in Glioma under Hypoxic Condition by Impairing microRNA-449b-5p-dependent AGAP2 Up-Regulation. Ther. Adv. Med. Oncol. 12, 1758835920940936. doi:10.1177/1758835920940936
Yu, M., Yi, B., Zhou, W., Gong, W., Li, G., and Yu, S. (2020). Linc00475 Promotes the Progression of Glioma by Regulating the miR-141-3p/YAP1 axis. J. Cel Mol Med. 25 (1), 463–472. doi:10.1111/jcmm.16100
Zheng, J., Zhou, Z., Qiu, Y., Wang, M., Yu, H., Wu, Z., et al. (2021). A Prognostic Ferroptosis-Related lncRNAs Signature Associated with Immune Landscape and Radiotherapy Response in Glioma. Front. Cel Dev. Biol. 9, 675555. doi:10.3389/fcell.2021.675555
Zheng, R., Yao, Q., Ren, C., Liu, Y., Yang, H., Xie, G., et al. (2016). Upregulation of Long Noncoding RNA Small Nucleolar RNA Host Gene 18 Promotes Radioresistance of Glioma by Repressing Semaphorin 5A. Int. J. Radiat. Oncology*Biology*Physics 96 (4), 877–887. doi:10.1016/j.ijrobp.2016.07.036
Zhou, K., Zhang, C., Yao, H., Zhang, X., Zhou, Y., Che, Y., et al. (2018). Knockdown of Long Non-coding RNA NEAT1 Inhibits Glioma Cell Migration and Invasion via Modulation of SOX2 Targeted by miR-132. Mol. Cancer 17 (1), 105. doi:10.1186/s12943-018-0849-2
Keywords: long non-coding RNAs, prognostic signature, glioma, the cancer genome atlas, chinese glioma genome atlas
Citation: He Y, Ye Y, Tian W and Qiu H (2021) A Novel lncRNA Panel Related to Ferroptosis, Tumor Progression, and Microenvironment is a Robust Prognostic Indicator for Glioma Patients. Front. Cell Dev. Biol. 9:788451. doi: 10.3389/fcell.2021.788451
Received: 02 October 2021; Accepted: 16 November 2021;
Published: 07 December 2021.
Edited by:
Maryam Mehrpour, Université Paris Descartes, FranceReviewed by:
Ian Kerr, University of Nottingham, United KingdomElla L. Kim, Johannes Gutenberg University Mainz, Germany
Copyright © 2021 He, Ye, Tian and Qiu. 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: Wei Tian, YTk4MzMxNDUxOEAxMjYuY29t; Huaide Qiu, SmFtZXM1NUBuam11LmVkdS5jbg==