Abstract
PURPOSE: To identify molecular subtypes of oxidative stress-related genes in head and neck squamous cell carcinoma (HNSCC) and to construct a scoring model of oxidative stress-related genes. METHODS: R language based scRNA-seq and bulk RNA-seq analyses were used to identify molecular isoforms of oxidative stress-related genes in HNSCC. An oxidative stress-related gene scoring (OSRS) model was constructed, which were verified through online data and immunohistochemical staining of clinical samples. RESULTS: Using TCGA-HNSCC datasets, nine predictive genes for overall patient survival, rarely reported in previous similar studies, were screened. AREG and CES1 were identified as prognostic risk factors. CSTA, FDCSP, JCHAIN, IFFO2, PGLYRP4, SPOCK2 and SPINK6 were identified as prognostic factors. Collectively, all genes formed a prognostic risk signature model for oxidative stress in HNSCC, which were validated in GSE41613, GSE103322 and PRJEB23709 datasets. Immunohistochemical staining of SPINK6 in nasopharyngeal cancer samples validated the gene panel. Subsequent analysis indicated that subgroups of the oxidative stress prognostic signature played important roles during cellular communication, the immune microenvironment, the differential activation of transcription factors, oxidative stress and immunotherapeutic responses. CONCLUSIONS: The risk model might predict HNSCC prognosis and immunotherapeutic responses.