Novel Immune-Related Ferroptosis Signature in Esophageal Cancer: An Informatics Exploration of Biological Processes Related to the TMEM161B-AS1/hsa-miR-27a-3p/GCH1 Regulatory Network

Background: Considering the role of immunity and ferroptosis in the invasion, proliferation and treatment of cancer, it is of interest to construct a model of prognostic-related differential expressed immune-related ferroptosis genes (PR-DE-IRFeGs), and explore the ferroptosis-related biological processes in esophageal cancer (ESCA). Methods: Four ESCA datasets were used to identify three PR-DE-IRFeGs for constructing the prognostic model. Validation of our model was based on analyses of internal and external data sets, and comparisons with past models. With the biological-based enrichment analysis as a guide, exploration for ESCA-related biological processes was undertaken with respect to the immune microenvironment, mutations, competing endogenous RNAs (ceRNA), and copy number variation (CNV). The model’s clinical applicability was measured by nomogram and correlation analysis between risk score and gene expression, and also immune-based and chemotherapeutic sensitivity. Results: Three PR-DE-IRFeGs (DDIT3, SLC2A3, and GCH1), risk factors for prognosis of ESCA patients, were the basis for constructing the prognostic model. Validation of our model shows a meaningful capability for prognosis prediction. Furthermore, many biological functions and pathways related to immunity and ferroptosis were enriched in the high-risk group, and the role of the TMEM161B-AS1/hsa-miR-27a-3p/GCH1 network in ESCA is supported. Also, the KMT2D mutation is associated with our risk score and SLC2A3 expression. Overall, the prognostic model was associated with treatment sensitivity and levels of gene expression. Conclusion: A novel, prognostic model was shown to have high predictive value. Biological processes related to immune functions, KMT2D mutation, CNV and the TMEM161B-AS1/hsa-miR-27a-3p/GCH1 network were involved in ESCA progression.