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: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: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.