Dataset Information

Identification and Validation of a Four-Gene Ferroptosis Signature for Predicting Overall Survival of Lung Squamous Cell Carcinoma

ABSTRACT:

Background

Lung squamous cell carcinoma (LUSC) represents 30% of all non-small cell lung carcinoma. Targeted therapy is not sufficient for LUSC patients because of the low frequency of targeted-effective mutation in LUSC whereas immunotherapy offers more options for patients with LUSC. We explored a ferroptosis-related prognostic signature that can potentially assess the prognosis and immunotherapy efficacy of LUSC patients.

Methods

A total of 502 LUSC patients were downloaded from The Cancer Genome Atlas (TCGA). The external validation data were obtained from the Gene Expression Omnibus (GEO): GSE73403. Then, we identified the candidate genes and constructed the prognostic signature through the Cox survival regression analyses and least absolute shrinkage and selection operator (LASSO). Risk score plot, Kaplan–Meier curve, and ROC curve were used to assess the prognostic power and performance of the model. The CIBERSORT algorithm estimated the fraction of immune cell types. TIDE was utilized to predict the response to immunotherapy. IMvigor210 was used to explore the association between the risk scores and immunotherapy outcomes. A nomogram combined selected clinical characteristics, and the risk scores were constructed.

Results

We screened 132 differentially expressed ferroptosis-related genes. According to KEGG and GO pathway analyses, these genes were mainly engaged in the positive regulation of cytokine production, cytokine metabolic process, and oxidoreductase activity. We then constructed a prognostic model via LASSO regression. The proportions of CD8+ T cells, CD4+ activated T cells, and follicular helper T cells were significantly different between low-risk and high-risk groups. TIDE algorithm indicated that low-risk LUSC patients might profit more from immune checkpoint inhibitors. The predictive value of the ferroptosis gene model in immunotherapy response was further confirmed in IMvigor210. Finally, we combined the clinical characteristics with a LASSO regression model to construct a nomogram that could be easily applied in clinical practice.

Conclusion

We identified a prognostic model that provides an accurate and objective basis for guiding individualized treatment decisions for LUSC.

SUBMITTER: Wang Q

PROVIDER: S-EPMC9330609 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background: Ferroptosis is an iron-dependent programmed cell death process. Recent studies have found that ferroptosis inducers hold promising potential in the treatment of lung adenocarcinoma (LUAD). However, the comprehensive analysis about the prognostic value of ferroptosis-related genes in LUAD remains to be elucidated. Methods: The RNA sequencing data and corresponding clinical information were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. A total of 259 ferroptosis-related genes were extracted from FerrDb website. The ferroptosis-related prognostic signature was developed by least absolute shrinkage and selection operator (LASSO) Cox regression analysis in TCGA LUAD cohort, and then validated by 5 independent GEO cohorts. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA) were performed to identify the difference in biological processes and functions between different risk groups. The expression levels of core prognostic genes were then verified in LUAD samples by immunohistochemistry (IHC) and erastin-treated LUAD cell lines by real-time polymerase chain reaction (PCR). The potential roles of GPX2 and DDIT4 as ferroptosis drivers in LUAD cell line were further confirmed by in vitro experiments. Results: A total of 20 intersecting genes between 70 ferroptosis-related DEGs and 45 potential prognostic genes were obtained for LASSO Cox regression analysis. The ferroptosis-related prognostic signature was developed by 7 core prognostic DEGs, and stratified LUAD patients into two risk groups. Kaplan-Meier analysis showed that the overall survival (OS) of LUAD patients in the high-risk group was significantly worse than that of the low-risk group. External validation of 5 independent GEO cohorts further confirmed that the ferroptosis-related prognostic signature was an ideal biomarker for predicting the survival of LUAD patients. Significant enrichment of fatty acid metabolism and cell cycle-related pathways were found in different risk groups. The expression patterns of 7 core prognostic genes in LUAD and adjacent normal lung tissues were validated by IHC, which was almost consistent with the results from public database. Furthermore, the changes related to cell cycle and ferroptosis after erastin treatment were also validated in LUAD cell lines. In addition, silencing GPX2 or DDIT4 could partially reverse the erastin-induced ferroptosis. Conclusion: In summary, the ferroptosis-related prognostic signature based on 7 core prognostic DEGs indicated superior predictive performance of LUAD patients. Targeting ferroptosis holds potential to be a therapeutic alternative for LUAD.

Project description:BackgroundHead and neck squamous cell carcinoma (HNSCC) is a common cancer characterized by late diagnosis and poor prognosis. The aim of this study was to identify a novel ferroptosis-related DNA methylation signature as an alternative diagnosis index for patients with HNSCC.MethodsMethylome and transcriptome data of 499 HNSCC patients, including 275 oral squamous cell carcinoma (OSCC) samples, were obtained from The Cancer Genome Atlas (TCGA). An additional independent methylation dataset of 50 OSCC patients from the NCBI Gene Expression Omnibus (GEO) database was used for validation. As an index of ferroptosis activity, the ferroptosis score (FS) of each patient was inferred from the transcriptome data using single-sample gene set enrichment analysis. Univariate, multivariate, and LASSO Cox regression analyses were used to select CpG sites for the construction of a ferroptosis-related DNA methylation signature for diagnosis of patients.ResultsWe initially inferred the FS of each TCGA HNSCC patient and divided the samples into high- and low-FS subgroups. Results showed that the high-FS subgroup displayed poor overall survival. Moreover, 378 differentially methylated CpG sites (DMCs) were identified between the two HNSCC subgroups, with 16 selected to construct a 16-DNA methylation signature for risk prediction in HNSCC patients using the LASSO and multivariate Cox regression models. Relative operating characteristic (ROC) curve analysis showed great predictive efficiency for 1-, 3-, and 5-year HNSCC survival using the 16-DNA methylation signature. Its predictive efficiency was also observed in OSCC patients from the TCGA and GEO databases. In addition, we found that the signature was associated with the fractions of immune types in the tumor immune microenvironment (TIME), suggesting potential interactions between ferroptosis and TIME in HNSCC progression.ConclusionsWe established a novel ferroptosis-related 16-DNA methylation signature that could be applied as an alternative tool to predict prognosis outcome in patients with HNSCC, including OSCC.

Project description:BackgroundThe prognosis of patients for lung adenocarcinoma (LUAD) is known to vary widely; the 5-year overall survival rate is just 63% even for the pathological IA stage. Thus, in order to identify high-risk patients and facilitate clinical decision making, it is vital that we identify new prognostic markers that can be used alongside TNM staging to facilitate risk stratification.MethodsWe used mRNA expression from The Cancer Genome Atlas (TCGA) cohort to identify a prognostic gene signature and combined this with clinical data to develop a predictive model for the prognosis of patients for lung adenocarcinoma. Kaplan-Meier curves, Lasso regression, and Cox regression, were used to identify specific prognostic genes. The model was assessed via the area under the receiver operating characteristic curve (AUC-ROC) and validated in an independent dataset (GSE50081) from the Gene Expression Omnibus (GEO).ResultsOur analyses identified a four-gene prognostic signature (CENPH, MYLIP, PITX3, and TRAF3IP3) that was associated with the overall survival of patients with T1-4N0-2M0 in the TCGA dataset. Multivariate regression suggested that the total risk score for the four genes represented an independent prognostic factor for the TCGA and GEO cohorts; the hazard ratio (HR) (high risk group vs low risk group) were 2.34 (p < 0.001) and 2.10 (p = 0.017). Immune infiltration estimations, as determined by an online tool (TIMER2.0) showed that CD4+ T cells were in relative abundance in the high risk group compared to the low risk group in both of the two cohorts (both p < 0.001). We established a composite prognostic model for predicting OS, combined with risk-grouping and clinical factors. The AUCs for 1-, 3-, 5- year OS in the training set were 0.750, 0.737, and 0.719; and were 0.645, 0.766, and 0.725 in the validation set. The calibration curves showed a good match between the predicted probabilities and the actual probabilities.ConclusionsWe identified a four-gene predictive signature which represents an independent prognostic factor and can be used to identify high-risk patients from different TNM stages of LUAD. A new prognostic model that combines a prognostic gene signature with clinical features exhibited better discriminatory ability for OS than traditional TNM staging.

Project description:BackgroundLung squamous cell carcinoma (LUSC), one of the main pathological types of lung cancer, has led to consequential socioeconomic burden. Ferroptosis is an iron-dependent form of cell death process with potentials for therapeutic target in various kinds of tumors. However, whether ferroptosis-related genes (FRGs) are associated with the prognosis of LUSC patients is still unclear. The aim of this study was to establish a FRGs-based signature which could stratify patients with LUSC.MethodsThe RNA sequencing profiles and corresponding clinical data of LUSC patients were retrieved from The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO) dataset. A FRG-based signature was developed using the TCGA-LUSC cohort and validated in the GEO cohort. Gene set enrichment analysis (GSEA) and analysis of immune cell characteristics were conducted to assess the relationship between FRGs and biological function or immune status. A nomogram based on selected clinical factors and the risk scores which were generated from the FRG-based signature was developed using the TCGA cohort and validated in the GEO cohort.ResultsA set of 16 FRGs, significantly associated with overall survival (OS) in the TCGA cohort, was identified and could classify LUSC patients into two risk groups. Kaplan-Meier analysis illustrated that the survival rate of the high-risk group was significantly lower than the low-risk group. Assessment and external validation of the signature showed that the survival predictive performance of this signature was adequate. Additionally, multiple pathways and functions were enriched through GSEA and the analysis of immune cell characteristics showed significantly different abundances of immune cells among the two risk groups. Finally, a nomogram integrating the FRG-based signature and selected clinical factors was also developed and assessed in both the TCGA and GEO cohort.ConclusionThis study indicated the association between the FRGs and prognosis of patients with LUSC. Targeting ferroptosis may serve as a novel potential therapeutic alternative for LUSC.

Project description:Cumulative studies have shown that RNA binding proteins (RBPs) play an important role in numerous malignant tumors and are related to the occurrence and progression of tumors. However, the role of RBPs in kidney renal clear cell carcinoma (KIRC) is not fully understood. In this study, we first downloaded gene expression data and corresponding clinical information of KIRC from the Cancer Genome Atlas (TCGA) database, International Cancer Genome Consortium (ICGC), and Gene Expression Omnibus (GEO) database, respectively. A total of 137 differentially expressed RBPs (DERBPs) were then identified between normal and tumor tissue, including 38 downregulated and 99 upregulated RBPs. Nine RBPs (EIF4A1, RPL36A, EXOSC5, RPL28, RPL13, RPS19, RPS2, EEF1A2, and OASL) were served as prognostic genes and exploited to construct a prognostic model through survival analysis. Kaplan-Meier curves analysis showed that the low-risk group had a better survival outcome when compared with the high-risk group. The area under the curve (AUC) value of the prognostic model was 0.713 in the TCGA data set (training data set), 0.706 in the ICGC data set, and 0.687 in the GSE29609 data set, respectively, confirming a good prognostic model. The prognostic model was also identified as an independent prognostic factor for KIRC survival by performing cox regression analysis. In addition, we also built a nomogram relying on age and the prognostic model and internal validation in the TCGA data set. The clinical benefit of the prognostic model was revealed by decision curve analysis (DCA). Gene set enrichment analysis revealed several crucial pathways (ERBB signaling pathway, pathways in cancer, MTOR signaling pathway, WNT signaling pathway, and TGF BETA signaling pathway) that may explain the underlying mechanisms of KIRC. Furthermore, potential drugs for KIRC treatment were predicted by the Connectivity Map (Cmap) database based on DERBPs, including several important drugs, such as depudecin and vorinostat, that could reverse KIRC gene expression, which may provide reference for the treatment of KIRC. In summary, we developed and validated a robust nine-RBP signature for KIRC prognosis prediction. A nomogram with risk score and age can be applied to promote the individualized prediction of overall survival in patients with KIRC. Moreover, the two drugs depudecin and vorinostat may contribute to KIRC treatment.