Development of Chromatin Regulator-related Molecular Subtypes and a Signature to Predict Prognosis and Immunotherapeutic Response in Head and Neck Squamous Cell Carcinoma

Abstract

Background: Chromatin regulators (CRs) serve as indispensable factors in tumor biological processes by influencing tumorigenesis and the immune microenvironment and have been identified in head and neck squamous cell carcinoma (HNSCC). Hence, CR-related genes (CRRGs) are considered potential biomarkers for predicting prognosis and immune infiltration in HNSCC. In this study, we established a novel signature for predicting the prognosis and immunotherapeutic response of HSNCC. Methods: A total of 870 CRRGs were obtained according to previous studies. Subsequently, patients in the TCGA-HNSC cohort were divided into different clusters based on the expression of prognostic CRRGs. Kaplan‒Meier (K‒M) survival analysis was conducted to compare the prognosis in clusters, and the CIBERSORT and ssGSEA methods assessed the immune infiltration status. In addition, the differences in immunotherapeutic responses were determined based on the TICA database. Furthermore, the differentially expressed CRRGs between clusters were identified, and the predictive signature was established according to the results of univariate Cox, least absolute shrinkage and selection operator regression analysis, and multivariate Cox. The predictive effects of the risk model were evaluated according to the area under the receiver operating characteristic (ROC) curve (AUC) in both the training and external test cohorts. A nomogram was established, and survival comparisons, functional enrichment analyses, and immune infiltration status and clinical treatment assessments were performed. In addition, the hub gene network and related analysis were conducted with the Cytohubba application. Results: Based on the expression of prognostic CRRGs, patients were divided into two clusters, in which Cluster 1 exhibited a better prognosis, more enriched immune infiltration, and a better immunotherapeutic response but exhibited chemotherapy sensitivity. The AUC values of the 1-, 3- and 5- year ROC curves for the risk model were 0.673, 0.732, and 0.692, respectively, as well as 0.645, 0.608, and 0.623 for the test set. In addition, patients in the low-risk group exhibited more immune cell enrichment and immune function activation, as well as a better immunotherapy response. The hub gene network indicated ACTN2 as the core gene differentially expressed between the two risk groups. Conclusions: We identified molecular subtypes and established a novel predictive signature based on CRRGs. This effective CRRS system can possibly provide a novel research direction for exploring the correlation between CRs and HNSCC and requires further experimental validation.