Immunization Combined with Ferroptosis Related Genes to Construct a New Prognostic Model for Head and Neck Squamous Cell Carcinoma

SIMPLE SUMMARY: Immunity combined with ferroptosis is being considered as a new tumor treatment modality, and its regulation in head and neck squamous cell carcinoma is still unknown. The purpose of this study was to look into the potential molecular biological roles of immune ferroptosis genes in head and neck squamous cell carcinoma. The 12-IFRM signatures were successfully constructed and classified into high- and low-risk groups using the TCGA database and related data resources. In patients with head and neck squamous cell carcinoma, feature-based risk scores were more predictive of survival than traditional clinicopathological features. Furthermore, the expression of CD8(+)T cells and macrophage M0 differed significantly between the two groups. The expression of TNFSF9 and CD44 in the high-risk groups was significantly increased compared with the low-risk groups. Next, we found a higher proportion of high-risk mutations than in the low-risk group. In addition, the high-risk group was more sensitive to some chemotherapy drugs. Finally, we performed correlation analysis on the model genes. In this paper, the 12-IFRM signatures was developed with promising application prospects for predicting the clinical outcomes and treatment outcomes in head and neck squamous cell carcinoma. ABSTRACT: Ferroptosis is a new type of programmed cell death that plays a pivotal role in a variety of tumors. Moreover, immunity is closely related to ferroptosis. However, immune-ferroptosis-related mRNAs (IFRMs) are still not fully understood in the regulation of head and neck squamous cell carcinoma (HNSC). The purpose of this paper was to investigate the IFRMs prediction of HNSC and its possible molecular biological role. RNA-Seq and related clinical data were mined from the TCGA database, ImmPort database, GeneCards database, FerrDb database, and previous data. In R software, the “DESeq2” package was used to analyze the differential expression of IFRMs. We used univariate Cox analysis to judge the prognosis of the IFRMs. Using the least absolute shrinkage and selection operator (LASSO) and Cox regression, a prediction model for 12 IFRMs was established. In this study, the Kaplan–Meier survival curve and receiver operating characteristic (ROC) curve analysis were used to evaluate the prediction results. Moreover, factors such as immune landscape, somatic mutations, and drug susceptibility are also discussed. We successfully constructed the signature of 12-IFRMs. The two risk groups were classified according to the risk score obtained by this signature. Compared with conventional clinicopathological features, the characteristic-based risk score was more predictive of survival in patients with HNSC. Furthermore, the expression of CD8(+)T cells and macrophage M0 differed significantly between the two groups. Moreover, the expression of TNFSF9 and CD44 in high-risk groups was significantly increased compared with the low-risk groups. Then, we found a higher proportion of high-risk mutations than in the low-risk group. Next, the high-risk group was more sensitive to chemotherapy drugs such as bosutinib, docetaxel, erlotinib, gefitinib, imatinib, lapatinib, and sorafenib. Finally, an in-depth analysis of the association and potential value of the 12 genes was performed. In summary, the 12-IFRM signatures established in this paper had good application prospects and could be effectively used to predict the clinical outcome and treatment response of head and neck squamous cell carcinoma.