RAD54-Like Protein 2 Is a Potential Diagnostic and Prognostic Biomarker in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a highly heterogeneous and aggressive malignancy with a poor prognosis. Pyroptosis triggered by gasdermins family proteins is reported vital for tumor microenvironment and cancer progression. However, pyroptosis-related gene expression and its relationship with immune infiltration and prognosis of HNSCC have not been fully defined. RNA-sequencing data of HNSCC patients were acquired from The Cancer Genome Atlas (TCGA) database. A pyroptosis-related gene expression signature and infiltrated immune cells were analyzed. Univariate, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression and nomogram analyses were used to construct a clinical-molecular risk model for survival prognosis. HNSCC was classified into three different molecular subtypes based on the expression information of pyroptosis-related genes. Immune cell infiltration was demonstrated to be distinct between the three subtypes. The segregation of patients into the high-risk group and low-risk group, were carried out using the signature of differential expression genes (DEGs) signature among the three molecular subtypes. The precision of this signature was corroborated by Receiver operating characteristic curve (ROC) analysis with the 3-year area under time-dependent ROC curve (AUC) reaching 0.711. The risk model was validated in another dataset from the Gene Expression Omnibus (GEO) database. Subsequently we established a clinical-molecular nomogram which combined the risk score with age and stage. The calibration plots for predicting the overall survival rate of 1-, 3-, and 5-years indicated that the nomogram performs well. The expression signature that encompasses pyroptosis-related genes could be used as molecular classification for HNSCC and pyroptosis might be a promising therapeutic target in HNSCC.

BackgroundThe homeobox gene Homeobox B7 (HOXB7) is overexpressed across a range of cancers and promotes tumorigenesis through varying effects on proliferation, survival, migration and invasion. However, its expression pattern and oncogenic role of HOXB7 in head and neck squamous cell carcinoma (HNSCC) remain largely unexplored. Here, we aimed to explore the expression pattern of HOXB7, its clinical significance as well as functional roles in HNSCC.MethodsHOXB7 mRNA expression in HNSCC was determined by data mining and analyses from TCGA (The Cancer Genome Atlas) and GEO (Gene Expression Omnibus) datasets. The protein abundance of HOXB7 was measured by immunohistochemistry in 119 primary HNSCC samples and associations between its expression and clinicopathological parameters and patient survival were evaluated. The pro-tumorigenic roles of HOXB7 in HNSCC were further delineated in vitro by loss-of-function assay. And a xenograft tumor model was established in nude mice to assess the role of HOXB7 in tumor growth. Connectivity Map (CMap) analysis was performed to identify bioactive small molecules which might be potential inhibitors for HOXB7.ResultsBioinformatics analyses showed that HOXB7 mRNA was significantly overexpressed in 8 independent HNSCC datasets from TCGA and GEO databases. HOXB7 protein was markedly upregulated in HNSCC samples as compared to normal counterparts and its overexpression significantly associated with high pathological grade, advanced clinical stage, cervical node metastasis (P = 0.0195, 0.0152, 0.0300) and reduced overall and disease-free survival (P = 0.0014, 0.0007). Univariate and multivariate Cox regression analyses further revealed HOXB7 as an independent prognostic factor for patients’ overall survival. Moreover, HOXB7 knockdown significantly inhibited cell proliferation, migration and invasion and induced cell apoptosis in HNSCC cells, and resulted in compromised tumour growth in vivo. Furthermore, CMap (Connectivity map) analysis has identified three potential bioactive small molecule inhibitors (NU-1025, thiamine, vinburnine) for HOXB7 targeted therapy in HNSCC.ConclusionsOur findings revealed that overexpression of HOXB7 was associates with tumour aggressiveness and unfavourable prognosis by serving a novel prognostic biomarker in HNSCC. Moreover, HOXB7 might be involved in the development and progression of HNSCC as an oncogene, and thereby might be a potential therapeutic target for HNSCC.

BackgroundLysophosphatidic acid (LPA) and its receptors play a key role in regulating cancer progression. Upregulation of LPA receptor 2 (LPAR2) plays a role in carcinogenesis; however, the exact role of LPAR2 in tumors remains elusive. This study aims to explore the correlation between LPAR2 expression with tumor prognosis and immune infiltration in pan-cancers.Materials and methodsThe expression of LPAR2 in pan-cancers was analyzed using the Online Cancer Microarray Database (Oncomine), Tumor Immune Estimation Resource (TIMER), and UALCAN databases. The effects of LPAR2 on the clinical prognosis in pan-cancer were examined using the Kaplan–Meier plotter (KM plotter) as well as Gene Expression Profiling Interactive Analysis (GEPIA), UALCAN, and Human Protein Atlas (HPA) databases. Moreover, the R software program was applied for validation of expression and prognostic value of LPAR2 in tumor patients in the Cancer Genome Atlas (TCGA) dataset and the Gene Expression Omnibus (GEO) database. The relationship between the expression level of LPAR2 and the clinical and molecular criteria of head and neck squamous cell carcinoma (HNSC) and kidney renal clear cell carcinoma (KIRC) was analyzed using UALCAN, whereas the relationship between LPAR2 expression and prognosis in patients with HNSC and KIRC with different clinical characteristics was examined using the KM plotter. Furthermore, the correlation between LPAR2 expression and tumor immune infiltration was examined using TIMER. The correlation between LPAR2 expression and gene markers of tumor immune infiltrates was analyzed using TIMER and GEPIA. In addition, the cBioPortal for Cancer Genomics was used to calculate the mutations, methylations, and altered neighbor genes of LPAR2.ResultsThe expression of LPAR2 was significantly correlated with the outcome of multiple types of cancer, especially HNSC and KIRC. Furthermore, high expression of LPAR2 was significantly associated with various immune markers in the immune cell subsets of HNSC and KIRC.ConclusionsHigh expression of LPAR2 plays significantly different prognostic roles in HNSC and KIRC possibly owing to its association with different immune markers. LPAR2 is correlated with tumor immune cell infiltration and is a valuable prognostic biomarker for HNSC and KIRC. However, further experiments are required to validate these findings.

Background Globally, head and neck squamous cell carcinoma (HNSCC) is a common malignant tumor with high morbidity and mortality. Hence, it is important to find effective biomarkers for the diagnosis and prediction of the prognosis of patients with HNSCC. FAM3D had been proven to be vital in other cancers. However, its predictive and therapeutic value in HNSCC is unclear. Therefore, it is valuable to explore the association between the expression level of FAM3D and its impacts on the prognosis and tumor microenvironment in HNSCC. Methods The Cancer Genome Atlas (TCGA) dataset, Genotype-Tissue Expression (GTEx) dataset, the Clinical Proteomic Tumor Analysis Consortium (CPTAC) dataset, and The Human Protein Atlas (THPA) website were used to assess HNSCC expressions in tumor and nontumor tissues. Then, we further conducted immunohistochemistry experiment as internal cohort to validate the same results. The Cox regression analysis, Kaplan–Meier analysis, and nomograms were performed to find the predictive prognostic value of FAM3D in HNSCC patients and its relationship with the clinicopathological features in HNSCC. The Gene Expression Omnibus (GEO) dataset was utilized to externally verify the prognosis value of FAM3D in HNSCC. Gene Set Enrichment Analysis (GESA) was applied to search the molecular and biological functions of FAM3D. The association between FAM3D and immune cell infiltration was investigated with the Tumor Immune Estimating Resource, version 2 (TIMER2). The relationships between FAM3D expression and tumor microenvironment (TME) scores, immune checkpoints, and antitumor compound half-maximal inhibitory concentration predictions were also explored. Results In different datasets, FAM3D mRNA and protein levels were all significantly lower in HNSCC tissues than in normal tissues, and they were strongly inversely associated with tumor grade, stage, lymph node metastasis, and T stage. Patients with high-FAM3D-expression displayed better prognosis than those with low-FAM3D-expression. FAM3D was also determined to be a suitable biomarker for predicting the prognosis of patients with HNSCC. This was externally validated in the GEO dataset. As for gene and protein level, the functional and pathway research results of FAM3D indicated that it was enriched in alteration of immune-related pathways in HNSCC. The low-expression group had higher stromal and ESTIMATE scores by convention than the high-expression group. FAM3D expression were found to be positively correlated with immune infiltrating cells, such as cancer-associated fibroblasts, myeloid-derived suppressor cells, macrophage cells, T cell CD8+ cells, regulatory T cells, and T cell follicular helper cells. FAM3D's relationships with immune checkpoints and sensitivity to antitumor drugs were also investigated. Conclusion Our study explored the impact of FAM3D as a favorable prognostic marker for HNSCC on the tumor immune microenvironment from multiple perspectives. The results may provide new insights into HNSCC-targeted immunotherapy.

The functional GRHL3-filaggrin axis maintains a tumor differentiation potential and influences drug sensitivity

IntroductionSalivary statherin, encoded by the STATH gene, maintains enamel integrity and oral microbial balance. Its downregulation has been linked to oral diseases, but its transcriptomic behavior in head and neck squamous cell carcinoma (HNSCC) remains largely unexplored. This study aimed to assess STATH as a diagnostic and prognostic biomarker in HNSCC and explore its potential biological role.Materials and methodsThe expression of STATH across 33 cancer types was analyzed via TIMER (Tumor Immune Estimation Resource), GEPIA (Gene Expression Profiling Interactive Analysis), and UALCAN (University of Alabama at Birmingham Cancer Data Analysis Portal), and promoter methylation and clinical correlations were explored in HNSCC. Prognostic significance was evaluated via Kaplan-Meier survival analysis (log-rank p<0.05). Immune cell infiltration was assessed by human papilloma virus (HPV) status via TIMER. Mutation data (n=674) were analyzed via cBioPortal. The findings were validated via the GEO GSE6631 dataset (n=22 paired samples). Co-expression and pathway enrichment analyses were conducted via UALCAN and Enrichr. Protein-protein interaction (PPI) networks were constructed via Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database.ResultsSTATH was significantly downregulated in HNSCC tissues compared with normal tissues according to TIMER (p<0.001) and GEPIA (p<0.01), with slightly reduced promoter methylation (median β=0.611; p=0.006) and rare genetic alterations (<1%). Higher STATH expression was associated with better overall survival (hazard ratio (HR)=0.914; log-rank p=0.033), especially in HPV-positive patients (HR=0.822; log-rank p=0.039). STATH expression was correlated with immune infiltration, particularly B cells, in HPV-positive tumors (r=0.355; p=0.001). STATH downregulation was validated (adjusted p=0.004). Co-expression analysis revealed enrichment in immune and oral health-associated pathways, including salivary protein's role in periodontitis (p=0.007; OR=118.8; LYZ, LTF) and dental caries (p=0.007; OR=59.36; LYZ, CCL28). PPIs interact with histatins and mucins.ConclusionSTATH is consistently downregulated in HNSCC, independent of promoter methylation or genetic mutations. Its suppression likely reflects tumor-driven remodeling of the salivary microenvironment, weakening host defenses and promoting disease progression. STATH has potential as a noninvasive diagnostic and prognostic biomarker.

BackgroundDNA methylation and miRNA-target genes play an important part in the early development of various tumors and have been studied as tumor biomarkers. Although previous studies have reported a cluster of molecular events (such as aberrant alterations of genomics and epigenetics), little is known of the potential biomarkers for early diagnosis and prognostic evaluation in head and neck squamous cell carcinoma (HNSCC).MethodsMultiple bioinformatics tools based on The Cancer Genome Atlas (TCGA) database and clinical samples were applied to evaluate the beneficial biomarkers in HNSCC. We focused on the role of plasminogen activator urokinase (PLAU), including diagnostic and prognostic significance, gene expression analysis, aberrant DNA methylation characteristics, interaction of miRNAs and associated signaling pathways.ResultsWe found that PLAU was aberrantly upregulated in HNSCC, regardless of the mRNA or protein level. The results of receiver operating characteristic (ROC) curve and Cox regression analysis revealed that PLAU was a diagnostic and independent prognostic factor for patients with HNSCC. Hypomethylation of PLAU was closely related to poor survival in HNSCC. Additionally, miR-23b-3p was predicted to target PLAU and was significantly downregulated in HNSCC tissues. Therefore, our findings suggested that PLAU functioned as a promoter in the pathological process of HNSCC. DNA hypomethylation and downregulation of miR-23b-3p were associated with PLAU overexpression. Finally, our findings provided evidence of a significant interaction between PLAU-target and miRNAs-target pathways, indicating that miR-23b-3p suppresses malignant properties of HNSCC by targeting PLAU via Ras/MAPK and Akt/mTOR signaling pathways.ConclusionsPLAU is overexpressed and may serve as an independent diagnostic and prognostic biomarker in HNSCC. Hypomethylation and downregulation of miR-23b-3p might account for the oncogenic role of PLAU in HNSCC.

Background: Radioresistance in head and neck squamous cell carcinoma (HNSCC) patients means response failure to current treatment. In order to screen radioresistant biomarkers and mechanisms associated with HNSCC, differentially expressed genes (DEGs) associated with radioresistance in HNSCC were investigated.Methods: The HNSCC cell line with radioresistance, Hep2-R, was established and detected the radiosensitivity using MTT, colony formation assay and flow cytometry analysis. Clariom™ D chip was applied to compare DEGs between Hep2 and Hep2-R groups and build the differential gene expression profiles associated with radioresistance in HNSCC. Bioinformatic analysis were used to find biological functions and pathways that related to radioresistance in HNSCC, including cell adhesion, cytochrome P450 and drug metabolism. Gene Expression Omnibus (GEO) datasets were selected to verify DEGs between HNSCC radioresistant cells and tissues. The representation of DEGs were validated between HNSCC patients with complete response and post-operative radiation therapy failure. In addition, we evaluated the clinical prognosis of DEGs using The Cancer Genome Atlas (TCGA) database.Results: 2,360 DEGs (|Fold Change|>1.5, p < 0.05) were identified between Hep2 and Hep2-R, including 1,144 upregulated DEGs and 1,216 downregulated DEGs. They were further verified by HNSCC radioresistant cells and tissues in GEO. 13 radioresistant DEGs showed same difference in expression level between cells and tissues. By comparing 13 DEGs with HNSCC patients, upregulations of FN1, SOX4 and ETV5 were found identical with above results. Only FN1 was a prognostic indicator of HNSCC in TCGA.Conclusion: FN1 is the potential novel biomarker for predicting poor prognosis and radioresistance in HNSCC patients. Overexpression of FN1 plays an important role in the tumorigenesis, prognosis and radioresistance of HNSCC.

PurposeTo explore the potential of AVPR2 in the immunotherapy of head and neck squamous cell carcinoma (HNSCC), thus providing insights into a novel antitumour strategy.MethodsIn this study, we performed a comprehensive analysis of the AVPR2 gene in HNSCC using public datasets from The Cancer Genome Atlas and Gene Expression Omnibus. We explored the potential molecular mechanism of HNSCC in clinical prognosis and tumour immunity from the aspects of gene expression, prognosis, immune subtypes, and immune infiltration.ResultsAVPR2 expression was significantly downregulated in primary HNSCC tissue compared with normal tissue. HNSCC patients with high AVPR2 expression had a better prognosis. Moreover, the results of GSEA showed that immune subtype surface AVPR2 is involved in immune modulation. Furthermore, significant strong correlations between AVPR2 expression and infiltrating immune cells existed in HNSCC, and marker genes of infiltrating immune cells were also significantly related to AVPR2 expression in HNSCC. These results suggest that AVPR2 expression can influence the infiltration of tumour immune cells. Finally, we found that only high levels of B-cell infiltration, rather than those of other immune cells, can predict a longer overall survival in patients with HNSCC. Future studies are needed to explore the role of AVPR2 and tumour-infiltrating B cells in HNSCC.ConclusionsThe AVPR2 gene may be a prognostic biomarker of HNSCC. Moreover, AVPR2 may play a role in HNSCC immune modulation, and the regulation of tumour-infiltrating B cells by AVPR2 may be a key link.

Head and neck squamous cell carcinoma (HNSCC) is the sixth leading cancer worldwide, with a high recurrence rate and a low cure rate. Phosphoglycerate kinase 1 (PGK1), an essential enzyme in the aerobic glycolysis pathway, is a prognostic marker for a variety of cancers. However, it remains unclear whether a PGK1-based immune signature can be used as a prognostic biomarker in HNSCC patients. We explored the potential oncogenic mechanisms of PGK1 by multiple bioinformatics analyses combined with multiple databases, including the correlation between PGK1 and prognosis, and the infiltration of immune cells in HNSCC. Functional enrichment analyses were further performed to investigate the potential role of PGK1 in HNSCC. The expression of PGK1 was significantly higher in HNSCC tissues compared to normal tissues. High expression of PGK1 was associated with poor prognosis in HNSCC, and multivariate cox regression analysis showed that PGK1 could be an independent prognostic factor in HNSCC. Pathway analysis revealed that PGK1 may regulate the pathogenesis of HNSCC through the immune signaling pathway. Moreover, PGK1 expression significantly correlated with the infiltration level of 16 types of immune cells. The current study reports that PGK1 expression was increased in HNSCC and that high PGK1 expression was closely associated with poor prognosis and immune cell infiltration, which could serve as a promising independent prognostic biomarker and potential immunotherapeutic target for HNSCC.

Head and neck squamous cell carcinoma (HNSCC) is a prevalent cancer that originates from the squamous cells. The role of the replication factor C subunit 3 (RFC3) in HNSCC progression remains elusive. The aim of this study was to uncover RFC3 significance in HNSCC. The Cancer Genome Atlas (TCGA-HNSCC) dataset was initially used to assess RFC3 expression and its association with HNSCC clinical features. Subsequently, quantitative reverse transcription PCR (RT-qPCR) confirmed RFC3 mRNA expression in oral squamous cell carcinoma (OSCC), a primary HNSCC type. Survival rates were evaluated using the Kaplan-Meier plot, while the Tumor Immune Estimation Resource (TIMER) database probed RFC3-immune cell interaction. Additionally, in silico tools were used to examine the RFC3 protein network and engagement in HNSCC pathways. RFC3 expression is significantly upregulated in HNSCC, including OSCC. Upregulated RFC3 expression was significantly correlated with the clinicopathological features of HNSCC, including tumor stage, grade, metastasis, and patient survival. RFC3 is also associated with immune cell infiltration. Functional analysis has highlighted its involvement in DNA replication, mismatch repair, and cell cycle pathways. Interestingly, RFC3 high expression is linked to well-known oncogenic signaling pathways, such as MYC/MYCN, HIPPO, and mTOR. In conclusion, RFC3 can be considered a novel prognostic biomarker for HNSCC, and further studies on its functional mechanisms may help to use RFC3 as a therapeutic target for HNSCC. The clinical relevance of this study lies in identifying RFC3 as a novel biomarker and prognostic indicator for HNSCC, offering insights that could impact future clinical approaches.

Telomere dysfunction has been identified as a biological marker of cancer progression in several types of cancer, including Head and Neck Squamous Cell Carcinoma (HNSCC). This study aimed to characterize the telomere maintenance genes (TMG)-related signature in prognosis and treatment response in HNSCC. The transcriptome and clinical data of HNSCC were obtained from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus databases, respectively. Non-negative matrix factorization (NMF) was used to identify molecular subtypes derived from TMG. Gene set enrichment analysis (GSEA) was performed to analyze the differentially expressed pathways between subtypes, and a risk score model derived from TMG was established. Kaplan-Meier survival analysis was used to evaluate inter-group prognostic features, and the correlation between TMG-derived molecular subtypes and risk score model with immune infiltration, immunotherapy, and chemosensitivity was assessed. Two HNSCC subtypes were identified based on 59 TMG-related genes, which exhibit significant heterogeneity in prognosis, immune cell infiltration, and treatment response. Additionally, a TMG-derived risk signature containing 9 genes was developed to assess the prognosis of HNSCC patients. The signature had significant predictive ability for HNSCC prognosis and was significantly correlated with immune cell infiltration and immunotherapy response. A nomogram integrating the risk signature, N stage and radiotherapy was constructed to predict 1-, 3-, and 5-year overall survival (OS) of HNSCC patients, which had better performance than other prognostic models and included TMG-derived risk score, radiotherapy, and N stage. This study identified TMG-derived molecular subtypes in HNSCC and developed a novel prognostic score model, highlighting the potential value of TMG in HNSCC prognosis and immunotherapy.

Head and neck squamous cell carcinoma (HNSCC) contributes significantly to global health challenges, presenting primarily in the oral cavity, pharynx, nasopharynx, and larynx. HNSCC has a high propensity for lymphatic metastasis. Diffuse large B-cell lymphoma (DLBCL), the most common subtype of non-Hodgkin lymphoma, exhibits significant heterogeneity and aggressive behavior, leading to high mortality rates. Epstein-Barr virus (EBV) is notably associated with DLBCL and certain types of HNSCC. The purpose of this study is to elucidate the molecular and immune interplay between HNSCC and DLBCL using bioinformatics and machine learning (ML) to identify shared biomarkers and potential therapeutic targets. Differentially expressed genes (DEGs) were identified using the "limma" package in R from the HNSCC dataset in The Cancer Genome Atlas (TCGA) database, and relevant modules were selected through weighted gene co-expression network analysis (WGCNA) from a DLBCL dataset in the Gene Expression Omnibus (GEO) database. Based on their intersection genes, functional enrichment analyses were conducted using Gene Ontology (GO), Disease Ontology, and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Protein-protein interaction (PPI) networks and ML algorithms were employed to screen for biomarkers. The prognostic value of these biomarkers was evaluated using Kaplan-Meier (K-M) survival analysis and receiver operating characteristic (ROC) curve analyses. The Human Protein Atlas (HPA) database facilitated the examination of messenger RNA (mRNA) and protein expressions. Further analyses of mutations, immune infiltration, drug predictions, and pan-cancer impacts were performed. Additionally, single-cell RNA sequencing (scRNA-seq) data analysis at the cell type level was conducted to provide deeper insights into the tumor microenvironment. From 2,040 DEGs and 1,983 module-related genes, 85 shared genes were identified. PPI analysis with six algorithms proposed 21 prospective genes, followed ML examination yielded 16 candidates. Survival and ROC analyses pinpointed four hub genes-ACACB, MMP8, PAX5, and TNFAIP6-as significantly associated with patient outcomes, demonstrating high predictive capabilities. Evaluations of mutations and immune infiltration, coupled with drug prediction and a comprehensive cancer analysis, highlighted these biomarkers' roles in tumor immune response and treatment efficacy. The scRNA-seq data analysis revealed an increased abundance of fibroblasts, epithelial cells and mononuclear phagocyte system (MPs) in HNSCC tissues compared to lymphoid tissues. MMP8 showed higher expression in five cell types in HNSCC tissues, while TNFAIP6 and PAX5 exhibited higher expression in specific cell types. Leveraging bioinformatics and ML, this study identified four pivotal genes with significant diagnostic capabilities for DLBCL and HNSCC. The survival analysis corroborates their diagnostic accuracy, supporting the development of a diagnostic nomogram to assist in clinical decision-making.

Bioinformatics analyses proposed xenotropic and polytropic retrovirus receptor 1 as a potential diagnostic and prognostic biomarker and immunotherapeutic target in head and neck squamous cell carcinoma

Background: Glutathione S-transferase kappa 1 (GSTK1) is critical in sarcoma and breast cancer (BRCA) development. However, the clinical significance of GSTK1 in head and neck squamous cell carcinoma (HNSC) remains unclear. This study is the first investigation into the role of GSTK1 in HNSC. Methods: All original data were downloaded from the Cancer Genome Atlas (TCGA) dataset and verified by R Base Package 4.2.0. The expression of GSTK1 in various cancers was explored with TIMER and TCGA databases. Prognostic value of GSTK1 was analyzed via survival module of Kaplan-Meier plotter and Human Protein Atlas database and Cox regression analysis. The association between GSTK1 and clinical features was evaluated by Wilcoxon signed-rank test and logistic regression analysis. The relationship between GSTK1 and immune infiltration and methylation level was further explored. The expression of GSTK1 and its correlation with immune cell infiltration was verified by Immunohistochemical staining (IHC). Results: GSTK1 was lower in HNSC, BRCA, Lung squamous cell carcinoma, and Thyroid carcinoma than in para-carcinoma. Low GSTK1 expression was associated with worse overall survival in Bladder urothelial carcinoma, Kidney renal papillary cell carcinoma, BRCA, and HNSC. However, only in BRCA and HNSC, GSTK1 expression in tumors was lower than that in normal tissues. Cox regression analyses confirmed that GSKT1 was an independent prognostic factor of overall survival in HNSC patients. The decrease in GSTK1 expression in HNSC was significantly correlated with high T stage and smoker history. IHC showed that the expression level of GSTK1 in HNSC was lower than that in para-carcinoma. In addition, GSEA showed that three pathways related to immune infiltration were positively correlated, while two pathways related to DNA methylation were negatively correlated with expression of GSTK1. Further analysis showed that GSTK1 was moderately positively correlated with the infiltration level of T cells and Cytotoxic cells, which was further confirmed by IHC. The methylation level of GSTK1 was associated with prognosis in patients with HNSC. Conclusion: Low GSTK1 expression may be a potential molecular marker for poor prognosis in HNSC and provide new insight for the development of diagnostic marker or therapeutic target.