Characterization of the immune cell function landscape in head and neck squamous carcinoma to assist in prognosis prediction and immunotherapy

Background: The malignant characteristics of cancer depend not only on intrinsic properties of cancer cells but also on the functions of infiltrating immune cells. In this study, we aimed to investigate the functional landscape of immune cells in head and neck squamous cell carcinoma (HNSCC). Methods: We employed single-sample gene set enrichment analysis to examine the immunophenotypes of HNSCC based on 29 immune cell functions (ICFs) in TCGA and GSE65858 datasets. We analyzed the clinical features, immune microenvironment, molecular profiles, and biological processes. Additionally, we developed and validated an ICF-based risk score for personalized prognosis prediction. We confirmed the value of the ICF score in our cohort using qRT-PCR and immunohistochemistry. Molecular docking was used to predict potential compounds for immunotherapy. Results: Three immunophenotypes (Immune-L, Immune-M, and Immune-H) were identified in 769 HNSCC samples. The characteristics of Immune-H were consistent with a “Hot” tumor, Immune-L was similar to a “Cold” tumor, and Immune-M exhibited intermediate features. The ICF risk score was associated with immune checkpoints, infiltrating immune cells, tumor mutation burden, and sensitivities to targeted/chemotherapeutic agents. Gene set variation analysis implicated the involvement of metabolic reprogramming pathways in the high-risk group. The combination of “Tumor Immune Dysfunction and Exclusion” and “Immunophenoscore” algorithms indicated that the low-risk group had a higher likelihood of benefiting from immunotherapy. Finally, we identified Eltrombopag and other compounds that may be beneficial for HNSCC immunotherapy. Conclusion: Our study provides a novel perspective on the tumor microenvironment of HNSCC, aiding in the understanding of HNSCC heterogeneity and the development of personalized/precision medicine.