In addition to the differences detected between HPV-positive and HPV-negative tumors, we observed substantial variability in the density of tumor-infiltrating lymphocytes and CD20+ B cell/CD8+ T cell interactions within the group of patients with HPV-associated tumors, splitting HPV-positive samples into “hot” and “cold” subgroups. Here, MS4A1 is linked to neoplasm.