Patients with human papillomavirus-negative head and neck squamous cell carcinoma (HPV-negative HNSCC) have worse outcomes than HPV-positive HNSCC. In our study, we used a published dataset and investigated the microbes enriched in molecularly classified tumor groups. We showed that microbial signatures could distinguish Hypoxia/Immune phenotypes similar to the gene expression signatures. Furthermore, we identified three highly-correlated microbes with immune processes that are crucial for immunotherapy response. The survival of patients in a molecularly heterogenous group shows significant differences based on the co-abundance of the three microbes. Overall, we present evidence that tumor-associated microbiota are critical components of the tumor ecosystem that may impact tumor microenvironment and immunotherapy response. The results of our study warrant future investigation to experimentally validate the conclusions, which have significant impacts on clinical decision-making, such as treatment selection.

Background: Solid tumors vary by the immunogenic potential of the tumor microenvironment (TME) and the likelihood of response to immunotherapy. The emerging literature has identified key immune cell populations that significantly impact immune activation or suppression within the TME. This study investigated candidate T-cell populations and their differential infiltration within different tumor types as estimated from mRNA co-expression levels of the corresponding cellular markers. Methods: We analyzed the mRNA co-expression levels of cellular biomarkers that define stem-like tumor-infiltrating lymphocytes (TILs), tissue-resident memory T-cells (TRM), early dysfunctional T-cells, late dysfunctional T-cells, activated-potentially anti-tumor (APA) T-cells and Butyrophilin 3A (BTN3A) isoforms, utilizing clinical and transcriptomic data from 1892 patients diagnosed with melanoma, bladder, ovarian, or pancreatic carcinomas. Real-world data were collected under the Total Cancer Care Protocol and the Avatar^® project (NCT03977402) across 18 cancer centers. Furthermore, we compared the survival outcomes following immune checkpoint inhibitors (ICIs) based on immune cell gene expression. Results: In melanoma and bladder cancer, the estimated infiltration of APA T-cells differed significantly (p = 4.67 × 10⁻¹² and p = 5.80 × 10⁻¹², respectively) compared to ovarian and pancreatic cancers. Ovarian cancer had lower TRM T-cell infiltration than melanoma, bladder, and pancreatic (p = 2.23 × 10⁻⁸, 3.86 × 10⁻²⁸, and 7.85 × 10⁻⁹, respectively). Similar trends were noted with stem-like, early, and late dysfunctional T-cells. Melanoma and ovarian expressed BTN3A isoforms more than other malignancies. Higher densities of stem-like TILs; TRM, early and late dysfunctional T-cells; APA T-cells; and BTN3A isoforms were associated with increased survival in melanoma (p = 0.0075, 0.00059, 0.013, 0.005, 0.0016, and 0.041, respectively). The TRM gene signature was a moderate predictor of survival in the melanoma cohort (AUROC = 0.65), with similar findings in testing independent public datasets of ICI-treated patients with melanoma (AUROC 0.61–0.64). Conclusions: Key cellular elements related to immune activation are more heavily infiltrated within ICI-responsive versus non-responsive malignancies, supporting a central role in anti-tumor immunity. In melanoma patients treated with ICIs, higher densities of stem-like TILs, TRM T-cells, early dysfunctional T-cells, late dysfunctional T-cells, APA T-cells, and BTN3A isoforms were associated with improved survival.