CD8 + T-cell marker genes reveal different immune subtypes of oral lichen planus by integrating single-cell RNA-seq and bulk RNA-sequencing.

BACKGROUND: Oral lichen planus (OLP) is a local autoimmune disease induced by T-cell dysfunction that frequently affects middle-aged or elderly people, with a higher prevalence in women. CD8 + T cells, also known as killer T cells, play an important role in the progression and persistence of OLP. In order to identify different OLP subtypes associated with CD8 + T cell pathogenesis, consensus clustering was used. METHODS: In this study, we preprocessed and downscaled the OLP single-cell dataset GSE211630 cohort downloaded from Gene Expression Omnibus (GEO) to finally obtain the marker genes of CD8 + T cells. Based on the expression of marker genes, we classified OLP patients into CMGs subtypes using unsupervised clustering analysis. The gene expression profiles were analyzed by WGCNA using the "WGCNA" R package based on the clinical disease traits and typing results, and 108 CD8 + T-cell related OLP pathogenicity-related genes were obtained from the intersection. Patients were once again classified into gene subtypes based on intersection gene expression using unsupervised clustering analysis. RESULTS: After obtaining the intersecting genes of CD8 + T cells related to pathogenesis, OLP patients can be precisely classified into two different subtypes based on unsupervised clustering analysis, and subtype B has better immune infiltration results, providing clinicians with a reference for personalized treatment. CONCLUSIONS: Classification of OLP into different subtypes improve our current understanding of the underlying pathogenesis of OLP and provides new insights for future studies.

Deciphering a Prognostic Signature Based on Soluble Mediators Defines the Immune Landscape and Predicts Prognosis in HNSCC.

BACKGROUND: The study on Head and Neck Squamous Cell Carcinoma (HNSCC), a prevalent and aggressive form of head and neck cancer, focuses on the often-overlooked role of soluble mediators. The objective is to leverage a transcriptome-based risk analysis utilizing soluble mediator-related genes (SMRGs) to provide novel insights into prognosis and immunotherapy efficacy in HNSCC patients. METHODS: We analyzed the expression and prognostic significance of 10,859 SMRGs using 502 HNSCC and 44 normal samples from the TCGA-HNSC cohort in The Cancer Genome Atlas (TCGA). The samples were divided into training and test sets in a 7:3 ratio, with an additional external validation using 40 tumor samples from the International Cancer Genome Consortium (ICGC). Key differentially expressed genes (DEGs) with prognostic significance were identified through univariate and Lasso-Cox regression analyses. A prognostic model based on 20 SMRGs was developed using Lasso and multivariate Cox regression. We assessed the clinical outcomes and immune status in high-risk (HR) and low-risk (LR) HNSCC patients utilizing the BEST databases and single-sample Gene Set Enrichment Analysis (ssGSEA). RESULTS: The 20 SMRGs were crucial in predicting the prognosis of HNSCC, with the SMRG signature emerging as an independent prognostic indicator. Patients classified in the HR group exhibited poorer outcomes compared to those in the LR group. A nomogram, integrating clinical characteristics and risk scores, demonstrated substantial prognostic value. Immunotherapy appeared to be more effective in the LR group, possibly attributed to enhanced immune infiltration and expression of immune checkpoints. CONCLUSIONS: The model based on soluble mediator-associated genes offers a fresh perspective for assessing the pre-immune efficacy and showcases robust predictive capabilities. This innovative approach holds significant promise in advancing the field of precision immuno-oncology research, providing valuable insights for personalized treatment strategies in HNSCC.

Corrigendum: Machine learning to construct sphingolipid metabolism genes signature to characterize the immune landscape and prognosis of patients with uveal melanoma.

[This corrects the article DOI: 10.3389/fendo.2022.1056310.].

An Aggrephagy-Related LncRNA Signature for the Prognosis of Pancreatic Adenocarcinoma.

Pancreatic adenocarcinoma (PAAD) is a common, highly malignant, and aggressive gastrointestinal tumor. The conventional treatment of PAAD shows poor results, and patients have poor prognosis. The synthesis and degradation of proteins are essential for the occurrence and development of tumors. Aggrephagy is a type of autophagy that selectively degrades aggregated proteins. It decreases the formation of aggregates by degrading proteins, thus reducing the harm to cells. By breaking down proteins, it decreases the formation of aggregates; thus, minimizing damage to cells. For evaluating the response to immunotherapy and prognosis in PAAD patients, in this study, we developed a reliable signature based on aggrephagy-related genes (ARGs). We obtained 298 AGGLncRNAs. Based on the results of one-way Cox and LASSO analyses, the lncRNA signature was constructed. In the risk model, the prognosis of patients in the low-risk group was noticeably better than that of the patients in the high-risk group. Additionally, the ROC curves and nomograms validated the capacity of the risk model to predict the prognosis of PAAD. The patients in the low-risk and high-risk groups showed considerable variations in functional enrichment and immunological analysis. Regarding drug sensitivity, the low-risk and high-risk groups had different half-maximal inhibitory concentrations (IC50).

T-cell exhaustion signatures characterize the immune landscape and predict HCC prognosis via integrating single-cell RNA-seq and bulk RNA-sequencing.

Background: Hepatocellular carcinoma (HCC), the third most prevalent cause of cancer-related death, is a frequent primary liver cancer with a high rate of morbidity and mortality. T-cell depletion (TEX) is a progressive decline in T-cell function due to continuous stimulation of the TCR in the presence of sustained antigen exposure. Numerous studies have shown that TEX plays an essential role in the antitumor immune process and is significantly associated with patient prognosis. Hence, it is important to gain insight into the potential role of T cell depletion in the tumor microenvironment. The purpose of this study was to develop a trustworthy TEX-based signature using single-cell RNA-seq (scRNA-seq) and high-throughput RNA sequencing, opening up new avenues for evaluating the prognosis and immunotherapeutic response of HCC patients. Methods: The International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA) databases were used to download RNA-seq information for HCC patients. The 10x scRNA-seq. data of HCC were downloaded from GSE166635, and UMAP was used for clustering descending, and subgroup identification. TEX-related genes were identified by gene set variance analysis (GSVA) and weighted gene correlation network analysis (WGCNA). Afterward, we established a prognostic TEX signature using LASSO-Cox analysis. External validation was performed in the ICGC cohort. Immunotherapy response was assessed by the IMvigor210, GSE78220, GSE79671, and GSE91061cohorts. In addition, differences in mutational landscape and chemotherapy sensitivity between different risk groups were investigated. Finally, the differential expression of TEX genes was verified by qRT-PCR. Result: 11 TEX genes were thought to be highly predictive of the prognosis of HCC and substantially related to HCC prognosis. Patients in the low-risk group had a greater overall survival rate than those in the high-risk group, according to multivariate analysis, which also revealed that the model was an independent predictor of HCC. The predictive efficacy of columnar maps created from clinical features and risk scores was strong. Conclusion: TEX signature and column line plots showed good predictive performance, providing a new perspective for assessing pre-immune efficacy, which will be useful for future precision immuno-oncology studies.

Circadian rhythm-related genes index: A predictor for HNSCC prognosis, immunotherapy efficacy, and chemosensitivity.

Background: Head and neck squamous cell carcinoma (HNSCC) is the most common head and neck cancer and is highly aggressive and heterogeneous, leading to variable prognosis and immunotherapy outcomes. Circadian rhythm alterations in tumourigenesis are of equal importance to genetic factors and several biologic clock genes are considered to be prognostic biomarkers for various cancers. The aim of this study was to establish reliable markers based on biologic clock genes, thus providing a new perspective for assessing immunotherapy response and prognosis in patients with HNSCC. Methods: We used 502 HNSCC samples and 44 normal samples from the TCGA-HNSCC dataset as the training set. 97 samples from GSE41613 were used as an external validation set. Prognostic characteristics of circadian rhythm-related genes (CRRGs) were established by Lasso, random forest and stepwise multifactorial Cox. Multivariate analysis revealed that CRRGs characteristics were independent predictors of HNSCC, with patients in the high-risk group having a worse prognosis than those in the low-risk group. The relevance of CRRGs to the immune microenvironment and immunotherapy was assessed by an integrated algorithm. Results: 6-CRRGs were considered to be strongly associated with HNSCC prognosis and a good predictor of HNSCC. The riskscore established by the 6-CRRG was found to be an independent prognostic factor for HNSCC in multifactorial analysis, with patients in the low-risk group having a higher overall survival (OS) than the high-risk group. Nomogram prediction maps constructed from clinical characteristics and riskscore had good prognostic power. Patients in the low-risk group had higher levels of immune infiltration and immune checkpoint expression and were more likely to benefit from immunotherapy. Conclusion: 6-CRRGs play a key predictive role for the prognosis of HNSCC patients and can guide physicians in selecting potential responders to prioritise immunotherapy, which could facilitate further research in precision immuno-oncology.

Uncovering the potential role of oxidative stress in the development of periodontitis and establishing a stable diagnostic model via combining single-cell and machine learning analysis.

Background: The primary pathogenic cause of tooth loss in adults is periodontitis, although few reliable diagnostic methods are available in the early stages. One pathological factor that defines periodontitis pathology has previously been believed to be the equilibrium between inflammatory defense mechanisms and oxidative stress. Therefore, it is necessary to construct a model of oxidative stress-related periodontitis diagnostic markers through machine learning and bioinformatic analysis. Methods: We used LASSO, SVM-RFE, and Random Forest techniques to screen for periodontitis-related oxidative stress variables and construct a diagnostic model by logistic regression, followed by a biological approach to build a Protein-Protein interaction network (PPI) based on modelled genes while using modelled genes. Unsupervised clustering analysis was performed to screen for oxidative stress subtypes of periodontitis. we used WGCNA to explore the pathways correlated with oxidative stress in periodontitis patients. Networks. Finally, we used single-cell data to screen the cellular subpopulations with the highest correlation by scoring oxidative stress genes and performed a proposed temporal analysis of the subpopulations. Results: We discovered 3 periodontitis-associated genes (CASP3, IL-1ß, and TXN). A characteristic line graph based on these genes can be helpful for patients. The primary hub gene screened by the PPI was constructed, where immune-related and cellular metabolism-related pathways were significantly enriched. Consistent clustering analysis found two oxidative stress categories, with the C2 subtype showing higher immune cell infiltration and immune function ratings. Therefore, we hypothesized that the high expression of oxidative stress genes was correlated with the formation of the immune environment in patients with periodontitis. Using the WGCNA approach, we examined the co-expressed gene modules related to the various subtypes of oxidative stress. Finally, we selected monocytes for mimetic time series analysis and analyzed the expression changes of oxidative stress genes with the mimetic time series axis, in which the expression of JUN, TXN, and IL-1ß differed with the change of cell status. Conclusion: This study identifies a diagnostic model of 3-OSRGs from which patients can benefit and explores the importance of oxidative stress genes in building an immune environment in patients with periodontitis.

Identification and validation of neurotrophic factor-related gene signatures in glioblastoma and Parkinson's disease.

Background: Glioblastoma multiforme (GBM) is the most common cancer of the central nervous system, while Parkinson's disease (PD) is a degenerative neurological condition frequently affecting the elderly. Neurotrophic factors are key factors associated with the progression of degenerative neuropathies and gliomas. Methods: The 2601 neurotrophic factor-related genes (NFRGs) available in the Genecards portal were analyzed and 12 NFRGs with potential roles in the pathogenesis of Parkinson's disease and the prognosis of GBM were identified. LASSO regression and random forest algorithms were then used to screen the key NFRGs. The correlation of the key NFRGs with immune pathways was verified using GSEA (Gene Set Enrichment Analysis). A prognostic risk scoring system was constructed using LASSO (Least absolute shrinkage and selection operator) and multivariate Cox risk regression based on the expression of the 12 NFRGs in the GBM cohort from The Cancer Genome Atlas (TCGA) database. We also investigated differences in clinical characteristics, mutational landscape, immune cell infiltration, and predicted efficacy of immunotherapy between risk groups. Finally, the accuracy of the model genes was validated using multi-omics mutation analysis, single-cell sequencing, QT-PCR, and HPA. Results: We found that 4 NFRGs were more reliable for the diagnosis of Parkinson's disease through the use of machine learning techniques. These results were validated using two external cohorts. We also identified 7 NFRGs that were highly associated with the prognosis and diagnosis of GBM. Patients in the low-risk group had a greater overall survival (OS) than those in the high-risk group. The nomogram generated based on clinical characteristics and risk scores showed strong prognostic prediction ability. The NFRG signature was an independent prognostic predictor for GBM. The low-risk group was more likely to benefit from immunotherapy based on the degree of immune cell infiltration, expression of immune checkpoints (ICs), and predicted response to immunotherapy. In the end, 2 NFRGs (EN1 and LOXL1) were identified as crucial for the development of Parkinson's disease and the outcome of GBM. Conclusions: Our study revealed that 4 NFRGs are involved in the progression of PD. The 7-NFRGs risk score model can predict the prognosis of GBM patients and help clinicians to classify the GBM patients into high and low risk groups. EN1, and LOXL1 can be used as therapeutic targets for personalized immunotherapy for patients with PD and GBM.

FAM family gene prediction model reveals heterogeneity, stemness and immune microenvironment of UCEC.

Background: Endometrial cancer (UCEC) is a highly heterogeneous gynecologic malignancy that exhibits variable prognostic outcomes and responses to immunotherapy. The Familial sequence similarity (FAM) gene family is known to contribute to the pathogenesis of various malignancies, but the extent of their involvement in UCEC has not been systematically studied. This investigation aimed to develop a robust risk profile based on FAM family genes (FFGs) to predict the prognosis and suitability for immunotherapy in UCEC patients. Methods: Using the TCGA-UCEC cohort from The Cancer Genome Atlas (TCGA) database, we obtained expression profiles of FFGs from 552 UCEC and 35 normal samples, and analyzed the expression patterns and prognostic relevance of 363 FAM family genes. The UCEC samples were randomly divided into training and test sets (1:1), and univariate Cox regression analysis and Lasso Cox regression analysis were conducted to identify the differentially expressed genes (FAM13C, FAM110B, and FAM72A) that were significantly associated with prognosis. A prognostic risk scoring system was constructed based on these three gene characteristics using multivariate Cox proportional risk regression. The clinical potential and immune status of FFGs were analyzed using CiberSort, SSGSEA, and tumor immune dysfunction and rejection (TIDE) algorithms. qRT-PCR and IHC for detecting the expression levels of 3-FFGs. Results: Three FFGs, namely, FAM13C, FAM110B, and FAM72A, were identified as strongly associated with the prognosis of UCEC and effective predictors of UCEC prognosis. Multivariate analysis demonstrated that the developed model was an independent predictor of UCEC, and that patients in the low-risk group had better overall survival than those in the high-risk group. The nomogram constructed from clinical characteristics and risk scores exhibited good prognostic power. Patients in the low-risk group exhibited a higher tumor mutational load (TMB) and were more likely to benefit from immunotherapy. Conclusion: This study successfully developed and validated novel biomarkers based on FFGs for predicting the prognosis and immune status of UCEC patients. The identified FFGs can accurately assess the prognosis of UCEC patients and facilitate the identification of specific subgroups of patients who may benefit from personalized treatment with immunotherapy and chemotherapy.

A Bioinformatics-Based Analysis of an Anoikis-Related Gene Signature Predicts the Prognosis of Patients with Low-Grade Gliomas.

Low-grade glioma (LGG) is a highly aggressive disease in the skull. On the other hand, anoikis, a specific form of cell death induced by the loss of cell contact with the extracellular matrix, plays a key role in cancer metastasis. In this study, anoikis-related genes (ANRGs) were used to identify LGG subtypes and to construct a prognostic model for LGG patients. In addition, we explored the immune microenvironment and enrichment pathways between different subtypes. We constructed an anoikis-related gene signature using the TCGA (The Cancer Genome Atlas) cohort and investigated the differences between different risk groups in clinical features, mutational landscape, immune cell infiltration (ICI), etc. Kaplan-Meier analysis showed that the characteristics of ANRGs in the high-risk group were associated with poor prognosis in LGG patients. The risk score was identified as an independent prognostic factor. The high-risk group had higher ICI, tumor mutation load (TMB), immune checkpoint gene expression, and therapeutic response to immune checkpoint blockers (ICB). Functional analysis showed that these high-risk and low-risk groups had different immune statuses and drug sensitivity. Risk scores were used together with LGG clinicopathological features to construct a nomogram, and Decision Curve Analysis (DCA) showed that the model could enable patients to benefit from clinical treatment strategies.

Natural killer cell-related prognosis signature characterizes immune landscape and predicts prognosis of HNSCC.

Background: Head and neck squamous cell carcinoma (HNSCC), the most common head and neck cancer, is highly aggressive and heterogeneous, resulting in variable prognoses and immunotherapeutic outcomes. Natural killer (NK) cells play essential roles in malignancies' development, diagnosis, and prognosis. The purpose of this study was to establish a reliable signature based on genes related to NK cells (NRGs), thus providing a new perspective for assessing immunotherapy response and prognosis of HNSCC patients. Methods: In this study, NRGs were used to classify HNSCC from the TCGA-HNSCC and GEO cohorts. The genes were evaluated using univariate cox regression analysis based on the differential analysis of normal and tumor samples in TCGA-HNSCC conducted using the "limma" R package. Thereafter, we built prognostic gene signatures using LASSO-COX analysis. External validation was carried out in the GSE41613 cohort. Immunity analysis based on NRGs was performed via several methods, such as CIBERSORT, and immunotherapy response was evaluated by TIP portal website. Results: With the TCGA-HNSCC data, we established a nomogram based on the 17-NRGs signature and a variety of clinicopathological characteristics. The low-risk group exhibited a better effect when it came to immunotherapy. Conclusions: 17-NRGs signature and nomograms demonstrate excellent predictive performance and offer new perspectives for assessing pre-immune efficacy, which will facilitate future precision immuno-oncology research.

Cuprotosis Programmed-Cell-Death-Related lncRNA Signature Predicts Prognosis and Immune Landscape in PAAD Patients.

In terms of mortality and survival, pancreatic cancer is one of the worst malignancies. Known as a unique type of programmed cell death, cuprotosis contributes to tumor cell growth, angiogenesis, and metastasis. Cuprotosis programmed-cell-death-related lncRNAs (CRLs) have been linked to PAAD, although their functions in the tumor microenvironment and prognosis are not well understood. This study included data from the TCGA-PAAD cohort. Random sampling of PAAD data was conducted, splitting the data into two groups for use as a training set and test set (7:3). We searched for differentially expressed genes that were substantially linked to prognosis using univariate Cox and Lasso regression analysis. Through the use of multivariate Cox proportional risk regression, a risk-rating system for prognosis was developed. Correlations between the CRL signature and clinicopathological characteristics, tumor microenvironment, immunotherapy response, and chemotherapy sensitivity were further evaluated. Lastly, qRT-PCR was used to compare CRL expression in healthy tissues to that in tumors. Some CRLs are thought to have strong correlations with PAAD outcomes. These CRLs include AC005332.6, LINC02041, LINC00857, and AL117382.1. The CRL-based signature construction exhibited outstanding predictive performance and offers a fresh approach to evaluating pre-immune effectiveness, paving the way for future studies in precision immuno-oncology.

A novel anoikis-related gene signature predicts prognosis in patients with head and neck squamous cell carcinoma and reveals immune infiltration.

Background: Head and neck squamous cell carcinoma (HNSCC) is a highly aggressive disease with a poor prognosis for advanced tumors. Anoikis play a key role in cancer metastasis, facilitating the detachment and survival of cancer cells from the primary tumor site. However, few studies have focused on the role of anoikis in HNSC, especially on the prognosis. Methods: Anoikis-related genes (ANRGs) integrated from Genecards and Harmonizome portals were used to identify HNSCC subtypes and to construct a prognostic model for HNSCC patients. Also, we explored the immune microenvironment and enrichment pathways between different subtypes. Finally, we provide clinical experts with a novel nomogram based on ANRGs, with DCA curves indicating the potential clinical benefit of the model for clinical strategies. Results: We identified 69 survival-related HNSCC anoikis-related DEGs, from which 7 genes were selected to construct prognostic models. The prognostic risk score was identified as an independent prognostic factor. Functional analysis showed that these high and low risk groups had different immune status and drug sensitivity. Next risk scores were combined with HNSCC clinicopathological features together to construct a nomogram, and DCA analysis showed that the model could benefit patients from clinical treatment strategies. Conclusion: The predictive seven-gene signature and nomogram established in this study can assist clinicians in selecting personalized treatment for patients with HNSCC.

Machine learning to construct sphingolipid metabolism genes signature to characterize the immune landscape and prognosis of patients with uveal melanoma.

Background: Uveal melanoma (UVM) is the most common primary intraocular malignancy in adults and is highly metastatic, resulting in a poor patient prognosis. Sphingolipid metabolism plays an important role in tumor development, diagnosis, and prognosis. This study aimed to establish a reliable signature based on sphingolipid metabolism genes (SMGs), thus providing a new perspective for assessing immunotherapy response and prognosis in patients with UVM. Methods: In this study, SMGs were used to classify UVM from the TCGA-UVM and GEO cohorts. Genes significantly associated with prognosis in UVM patients were screened using univariate cox regression analysis. The most significantly characterized genes were obtained by machine learning, and 4-SMGs prognosis signature was constructed by stepwise multifactorial cox. External validation was performed in the GSE84976 cohort. The level of immune infiltration of 4-SMGs in high- and low-risk patients was analyzed by platforms such as CIBERSORT. The prediction of 4-SMGs on immunotherapy and immune checkpoint blockade (ICB) response in UVM patients was assessed by ImmuCellAI and TIP portals. Results: 4-SMGs were considered to be strongly associated with the prognosis of UVM and were good predictors of UVM prognosis. Multivariate analysis found that the model was an independent predictor of UVM, with patients in the low-risk group having higher overall survival than those in the high-risk group. The nomogram constructed from clinical characteristics and risk scores had good prognostic power. The high-risk group showed better results when receiving immunotherapy. Conclusions: 4-SMGs signature and nomogram showed excellent predictive performance and provided a new perspective for assessing pre-immune efficacy, which will facilitate future precision immuno-oncology studies.

Characterization of coagulation-related gene signature to predict prognosis and tumor immune microenvironment in skin cutaneous melanoma.

Backgroud: Skin cutaneous melanoma (SKCM) is an extremely metastatic form of skin cancer. However, there are few valuable molecular biomarkers, and accurate diagnosis is still a challenge. Hypercoagulable state encourages the infiltration and development of tumor cells and is significantly associated with poor prognosis in cancer patients. However, the use of a coagulation-related gene (CRG) signature for prognosis in SKCM, on the other hand, has yet to be determined. Method: We used data from The Cancer Genome Atlas (TCGA) and Genotype Tissue Expression (GTEx) databases to identify differentially expressed CRGs, then designed a prognostic model by using the LASSO algorithm, univariate and multivariate Cox regression analysis, and constructed a nomogram which was evaluated by calibration curves. Moreover, the Gene Expression Omnibus (GEO), GSE54467 was used as an independent validation. The correlation between risk score and clinicopathological characteristics, tumor microenvironment (TME), and immunotherapy was further analyzed. Results: To develop a prognostic model, seven CRGs in SKCM patients related to overall survival (OS) were selected: ANG, C1QA, CFB, DUSP6, KLKB1, MMP7, and RABIF. According to the Kaplan-Meier survival analysis, an increased OS was observed in the low-risk group than in the high-risk group (P<0.05). Immunotherapy was much more beneficial in the low-risk group, as per immune infiltration, functional enrichment, and immunotherapy analysis. Conclusions: The prognosis of SKCM patients may now be predicted with the use of a CRG prognostic model, thus guiding the development of treatment plans for SKCM patients and promoting OS rates.

Identification and validation of neurotrophic factor-related genes signature in HNSCC to predict survival and immune landscapes.

Background: Head and neck squamous cell carcinoma (HNSCC) is the seventh most common type of cancer worldwide. Its highly aggressive and heterogeneous nature and complex tumor microenvironment result in variable prognosis and immunotherapeutic outcomes for patients with HNSCC. Neurotrophic factor-related genes (NFRGs) play an essential role in the development of malignancies but have rarely been studied in HNSCC. The aim of this study was to develop a reliable prognostic model based on NFRGs for assessing the prognosis and immunotherapy of HNSCC patients and to provide guidance for clinical diagnosis and treatment. Methods: Based on the TCGA-HNSC cohort in the Cancer Genome Atlas (TCGA) database, expression profiles of NFRGs were obtained from 502 HNSCC samples and 44 normal samples, and the expression and prognosis of 2601 NFRGs were analyzed. TGCA-HNSC samples were randomly divided into training and test sets (7:3). GEO database of 97 tumor samples was used as the external validation set. One-way Cox regression analysis and Lasso Cox regression analysis were used to screen for differentially expressed genes significantly associated with prognosis. Based on 18 NFRGs, lasso and multivariate Cox proportional risk regression were used to construct a prognostic risk scoring system. ssGSEA was applied to analyze the immune status of patients in high- and low-risk groups. Results: The 18 NFRGs were considered to be closely associated with HNSCC prognosis and were good predictors of HNSCC. The multifactorial analysis found that the NFRGs signature was an independent prognostic factor for HNSCC, and patients in the low-risk group had higher overall survival (OS) than those in the high-risk group. The nomogram prediction map constructed from clinical characteristics and risk scores had good prognostic power. Patients in the low-risk group had higher levels of immune infiltration and expression of immune checkpoints and were more likely to benefit from immunotherapy. Conclusion: The NFRGs risk score model can well predict the prognosis of HNSCC patients. A nomogram based on this model can help clinicians classify HNSCC patients prognostically and identify specific subgroups of patients who may have better outcomes with immunotherapy and chemotherapy, and carry out personalized treatment for HNSCC patients.