HPV-associated head and neck cancer is characterized by distinct profiles of CD8+ T cells and myeloid-derived suppressor cells

Patients with HPV−-localized head and neck cancer (HNC) show inferior outcomes after surgery and radiochemotherapy compared to HPV-associated cancers. The underlying mechanisms remain elusive, but differences in immune status and immune activity may be implicated. In this study, we analyzed immune profiles of CD8+ T cells and myeloid-derived suppressor cells (MDSC) in HPV+ versus HPV− disease. The overall frequency of CD8+ T cells was reduced in HNC versus healthy donors but substantially increased after curative therapy (surgery and/or radiochemotherapy). In HPV+ patients, this increase was associated with significant induction of peripheral blood CD8+/CD45RA−/CD62L− effector memory cells. The frequency of HPV-antigen-specific CD8+ cells was low even in patients with virally associated tumors and dropped to background levels after curative therapy. Pre-therapeutic counts of circulating monocytic MDSC, but not PMN-MDSC, were increased in patients with HPV− disease. This increase was accompanied by reduced fractions of terminally differentiated CD8+ effector cells. HPV− tumors showed reduced infiltrates of CD8+ and CD45RO+ immune cells compared with HPV+ tumors. Importantly, frequencies of tumor tissue-infiltrating PMN-MDSC were increased, while percentages of Granzyme B+ and Ki-67+ CD8 T cells were reduced in patients with HPV− disease. We report differences in frequencies and relative ratios of MDSC and effector T cells in HPV− HNC compared with more immunogenic HPV-associated disease. Our data provide new insight into the immunological profiles of these two tumor entities and may be utilized for more tailored immunotherapeutic approaches in the future. Supplementary Information The online version contains supplementary material available at 10.1007/s00262-023-03571-8.


Introduction
Over the past decades, head and neck cancer (HNC), the worldwide 6th most common cancer, has experienced a significant change due to the identification of human papillomavirus (HPV)-associated disease, a molecular entity that is clearly distinct from HPV − HNC [1].In the oropharynx, the high-risk subtype HPV16 is-by far-the most frequent HPV subtype associated with the disease.Despite the fact that HPV + HNC is usually diagnosed at an advanced nodal infestation, therapeutic response and survival is significantly better in comparison with HPV − HNC patients [2][3][4].In some studies, the favorable prognosis was attributed to an increased sensitivity of HPV + HNC to radiation [5].However, recent evidence suggests that differences in immune profiles and response may contribute to the improved prognosis of HPV + HNC patients [4,6,7].In the peripheral blood, HNC is characterized by profound signs of immunosuppression associated with reduced T cells counts and induction of regulatory T cells as well as myeloid-derived suppressor cells (MDSC) [8][9][10][11].Although the prognostic relevance of circulating CD8 + T cells is less clear and demands further investigation, some studies already reported on the identification of CD8 + T cell subsets, which were associated with a more favorable outcome [12].Therefore, further investigations exploring potential differences in T cell subsets of HPV + and HPV − HNC patients are required, including functionality, antigen specificity, phenotypes, and-of utmost importance-possible causes for these differences.
Over the past decade, MDSC, which are induced in HNC and many other types of cancer [8], have emerged as important regulators of CD8 + T cell-mediated immunity.MDSC are a heterogeneous group of myeloid cells, and several mechanisms have been identified how these cells regulate and inhibit the activity of other immune cells, particularly T cells [13,14].Originally identified in mice [15], more recently, MDSC and subsets thereof have also been described in humans [16][17][18].
It was the aim of this study to investigate the frequency and subset distribution of CD8 + T cells and MDSC in patients with primary HNC.We focused our detailed longitudinal and comparative analysis of the CD8 + T cell compartment on a comparative analysis of HPV + versus HPV − HNC and hypothesized that differences in MDSC and CD8 + T cell subsets exist between these HNC subtypes.Our findings suggest different baseline levels and dynamics of MDSC and T cell phenotypes in HPV − versus HPV + HNC, which may contribute to differences in clinical outcome.

Patients and control subjects and HLA/HPV testing
Blood samples were taken from 89 head and neck cancer patients before the onset of therapy (pre-therapy sample) and again beginning from the first tumor-free check-up (post-therapy samples).Pre-therapy blood samples were used for HLA typing (performed at the Institute for Transfusion Medicine, University Hospital Essen).HPV status was determined in primary tumor material by an HPV family gene probe (INFORM HPV III Family 16, Roche Diagnostics, Mannheim, Germany) followed by a (PGMY09/11 & GP5 + /6 + ) nested PCR (performed by the Institute of Pathology).In addition, surplus tumor tissue samples were stained for pRb and p16 (see below).This staining pattern was used as a surrogate marker for HPV RNA-positivity according to a method described by Holzinger et al. [19].Blood samples from 18 healthy control donors were taken and analyzed accordingly.Experiments were approved by the local ethics committee, and informed written consent was obtained from each individual.

Preparations of tumor microarrays (TMAs)
Formalin-fixed and paraffin-embedded tumor tissue from HNC patients was cut into 1.5-µm-thick whole tissue sections.Hematoxylin and eosin (H&E)-stained sections were used to identify and mark representative tumor regions by an experienced pathologist (AB).Marked tumor regions were 1 3 punched out (3 mm diameter) and transferred into empty recipient blocks (spaced 1 mm apart).Tissue punches from kidney, liver, and healthy mucosa were added as controls.Blocks were sealed, 3-µm-thick sections were cut, transferred on slides (SuperFrost, R. Langenbrinck Labor-und Medizintechnik, Emmendingen, Germany), and dried over night at 60 °C.
For pRb and p16 INK4a the following score was used: Staining intensity was scored by using a 1 to 4 scale: 1 = none, 2 = low, 3 = medium, and 4 = high staining intensity.A sample was considered HPV16 RNA + when p16 INK4a was  "high" (scores 3 or 4) and pRB was "low" (scores 1 or 2).Any other combination was considered HPV16 RNA − [19].For CD8β and CD45RO the following score was used: The amount of stained tissue was scored on a scale from 1 to 4: 1 ≤ 10% stained, 2 = 10-50%, 3 = 50-80%, and 4 ≥ 80% stained.The score was applied to the epithelial tumor cell nests (tumor) and to the intratumoral stroma regions adjacent to the tumor cell areas.For Fig. 7E, we combined CD8β and CD45RO in tumor and stroma into an "immunoscorelike grading" analogous to a recently proposed method [22].A score of 2 or more (10% or more of tissue stained) was considered "high" for the epithelial tumor region.A score of 3 or more (50% or more of tissue stained) was considered "high" for the stromal regions.An IG0 (Immunograde 0) means that no "high" scores were assigned to the tissue; an IG4 means that four "high" scores were assigned.Immunogrades 0 and 1 were defined as "low infiltrate," and the immunogrades 2-4 were defined as "high infiltrate" group.IHC analysis was performed with a Zeiss Axioscope 2 microscope at 100 × and 200 × magnification (Zeiss, Jena, Germany).

Statistical analysis
Results are depicted as box plots with medians (and 5th, 25th, 75th and 95th percentiles shown) or scatter plots (with lines connecting samples from the same patient).Statistical analysis was performed by Mann−Whitney rank sum test, Kruskal-Wallis one-way analysis of variance followed by Dunn's method post hoc, calculation of linear regression or Chi-squared test.A p value ≤ 0.05 was taken as the level of significance.The statistical analysis was performed with Systat Sigmaplot 11.0 software (Systat Software, Erkrath, Germany).

Results
Various types of squamous cell carcinoma (SCC) induce systemic immune suppression, and alterations in the T cell compartment have been reported in patients with HNC [23,24].Therefore, we studied a series of 89 patients with localized or locally advanced HNC treated at the West German Cancer Center.Treatment strategies (surgery alone, surgery followed by adjuvant therapy and primary radiochemotherapy) were proposed by the multidisciplinary tumor board for head and neck cancers, and treatments were conducted using state-of-the art technique and facilities (Table 1).First, we analyzed the frequency, phenotype, and subset distribution of peripheral blood CD8 + T cells in our cohort of HNC patients compared to healthy donors (HD).We found that the relative frequency of CD8 + T cells within the T cell compartment was reduced in cancer patients (Fig. 1A).CCR7, CD62L, PD-1, and TIM-3 are important markers of T cell functionality.Analyzing expression of these markers by CD8 + T cells, we observed a significant relative reduction of CCR7 + and CD62L + CD8 cells in cancer patients (Fig. 1B,  C).Considering the expression of activation and exhaustion markers PD-1 and TIM-3, our analysis revealed no significant differences in CD8 + T cells between HNSCC patients and HD (Fig. 1D), with double negative cells representing the highest fractions in both groups.We next analyzed the subset composition of CD8 cells (naïve, central memory, effector memory, terminally differentiated effector cells) in HNC patients versus HD.Following the protocol delineated in Fig. 1E we applied CD62L and CD45RA as subset markers.While the frequency of T CM , T EM , and T TD was similar in HNC patients and HD, there was a striking and statistically significant reduction of naïve CD8 + T cells in HNC patients (Fig. 1E-F).These data support previous reports on low frequencies of CD8 + T cell populations in HNC patients and reveal further differences between the phenotypic T cell properties of HNC patients and HD.
Next, we conducted a longitudinal follow-up analysis of patients and explored therapy effects on the peripheral CD8 + T cell compartment.As illustrated in Fig. 2A, the relative frequency of CD8 + T cells strongly increased in patients after therapy.Interestingly, the percentage of CD8 + T cells with CCR7 expression further decreased after therapy (Fig. 2B).The latter finding prompted us to investigate treatment effects on CD8 + T cell subsets.In agreement with the reduced expression of CCR7, we observed a significant treatment-associated relative increase in CD8 + T EM and T TD cells (left versus right columns in Fig. 2C).Surprisingly, no change in immune Fig. 2 Pre-therapy to post-therapy changes in CD8 + T cells of HNC patients.Pre-therapy samples were obtained from tumor-bearing patients before the onset of therapy.Post-therapy samples were obtained at the first monitoring visit during routine tumor aftercare.Relative frequency of pre-therapy and post-therapy samples of CD8 + cells (A) and CCR7 + /CD8 + cells (B) was determined in HNC patients and healthy donors.CD8 + T cell subsets (C) and expression of PD-1 and TIM-3 (D) were quantified by flow cytometry (DN = double negative, DP = double positive).In paired columns of the same color the left column represents the pre-therapy sample and the right column represents the post-therapy sample.Note the post-therapy increase in T EM and T TD in (C).(E-G) shows a representative example of an extended longitudinal monitoring in a HNC patient.Blood was taken over a period of 15 month.Data are presented as box plots showing median with percentiles (25th, 75th) and 5th/95th whiskers.The Mann-Whitney Rank Sum test (A-B) and the Kruskal-Wallis ANOVA on ranks [post hoc: Dunn's method] (C-D) were used for statistical analysis.* indicates a p value < 0.05 ◂ checkpoint expression (PD-1, TIM-3) was observed in pre-versus post-therapeutic samples (Fig. 2D).For the majority of patients, multiple post-therapeutic time points were available for analysis and representative patient trajectories suggesting durable subset shifts are shown in Fig. 2E-G.
HPV-positivity in HNC is associated with improved treatment response and better outcomes [25].It has been suggested that improved anti-tumor immunity may contribute to this effect [26].We, therefore, compared the treatment effect on CD8 + T cells and HPV-specific immune responses in HPV + versus HPV − patients.Therapy-associated relative increases in CD8 + T cells and reduction of CCR7 + /CD8 + T cells could be confirmed for both HPV + and HPV − disease (Fig. 3A, B).However, a detailed analysis of the CD8 + subset composition revealed HPV-associated differences.A clear and statistically significant increase in peripheral blood CD8 + T EM post-therapy could only be observed in patients with HPV + tumors, whereas the frequency of CD8 + T EM remained largely unchanged in patients with HPV − disease (Fig. 3C).Interestingly, the frequency of PD-1 + and TIM3 + CD8 + T cells remained low after therapy for both patient groups.In sum, these data suggest a more pronounced treatment effect on CD8 + T cell phenotypes in HPV + HNC.
In previous studies, T cells with specificity for HPV antigens have been reported in HPV-associated HNC.These HPV-specific T cells were mainly isolated from tissue samples [27,28].However, some studies reported the existence of HPV-specific T cells also in the peripheral blood, albeit at very low frequencies [29,30].We devised the highly sensitive dextramer technology as well as a multi-peptide ELIS-POT to enumerate HPV-specific T cells pre-and post-therapy [31].While CD8 + T cells with HPV16 E7-specificity could hardly be detected in HPV − patients, the majority of patients with HPV16 + tumors displayed HPV-reactive cells at frequencies of more than 1 in 10,000 peripheral blood CD8 + T cells (Fig. 3E).Interestingly, associated with the removal of the primary tumor by curative therapy (primary or adjuvant radiotherapy), the frequency of these HPV16specific CD8 + T cells decreased in all patients tested in our study (Fig. 3E, lower panel).Functional reactivity to HPV16 E6/E7 multi-peptides was tested by ELISPOT technology.The frequency of HPV-responsive IFN-γ producing cells was very low in pretherapeutic samples.Low-level specific responses could be detected in some patients after therapy, but no clear correlation with the HPV status was evident (Fig. 3F).
Recent molecular analyses suggested the importance of genetically HPV-driven oncogenesis (HPV RNA positivity) compared to the presence of HPV DNA only [32].We used a method developed by Holzinger et al. [19] as a surrogate for RNA positivity.To this end, HPV DNA + samples from patients with oropharyngeal SCC were stained for p16 and pRb, and tissues with high expression of p16 and simultaneous low expression of pRb were classified as HPV16 RNA + (Fig. 4A).Using this classification, no significant differences between patient stratification according to DNA positivity (Fig. 3) and HPV RNA positivity (Fig. 4) were found.Thus, our data suggest that therapy-associated changes in CD8 + T Fig. 3 Analysis of CD8 + T cells in HPV − versus HPV + patients.HPV16 status was determined by a HPV family gene probe followed by a nested PCR to determine viral subtype if HPV + .The relative frequency of CD8 + cells (A) and CCR7 + /CD8 + cells (B) was determined in pre-therapy and post-therapy samples by flow cytometry.CD8 + T cell subsets (C) and expression of PD-1 and TIM-3 (D) were quantified by flow cytometry (DN = double negative, DP = double positive).In paired columns of the same color, the left column represents the pre-therapy sample and the right column represents the post-therapy sample.Note the post-therapy increase of T EM in HPV + patients (C).(E) The frequency of HPV16 E7-specific CD8 + T cells was determined by flow cytometric dextramer analysis.Note the decrease of HPV-specific T cells in post-therapy samples of HPV + patients.(F) IFN-γ ELISpot was used to determine the HPV16-specific cytokine response of CD8 + T cells, positive control: Phytohemagglutinin-L (PHA).Data are presented as box plots showing median with percentiles (25th, 75th) and 5th/95th whiskers.The Mann-Whitney rank sum test (A-B, E) and the Kruskal-Wallis ANOVA on ranks [post hoc: Dunn's method] (C-D) were used for statistical analysis.* indicates a p value < 0.05  cells occur in both HPV DNA + and RNA + patients and are independent of the mode of therapy (Fig. 4B, C).MDSCs regulate the activation and function of T cells in multiple murine cancer models.While the identity and function of MDSC are less clear in cancer patients, a similar function for MDSC has also been proposed in humans.Employing a previously established protocol [16,20], we determined the frequency of the three major human MDSC subsets in the peripheral blood of our patient cohort (Fig. 5A, B).In agreement with our previous studies [8], we found an increase in both monocytic MDSC (M-MDSC) and polymorphonuclear MDSC (PMN-MDSC) in HNC patients as compared to HD.While the frequency of PMN-MDSC was similar between patients with HPV + and HPV − tumors, high levels of M-MDSC (above 1%) were almost exclusively detected in patients with HPV − tumors.This finding was further supported by a correlation of high M-MDSC frequencies and higher expression of CD62L on CD8 + T cells (Fig. 5C).Also, a high number of M-MDSC moderately correlated with reduced fractions of T TD effector cells (Fig. 5D).These findings establish a correlation between MDSC frequencies and CD8 + T cells phenotypes in HNC In order to support and enhance our findings from the peripheral blood, we performed an analysis of the tumor microenvironment in a second cohort (Table 2) and investigated MDSC and T cells in tumor tissue islands as well as surrounding stroma.Because no robust consensus markers for tissue M-MDSC exist, only PMN-MDSC were investigated by co-staining of CD66b and LOX1 (lectin-type oxidized LDL receptor 1) expression [14].To further investigate CD8 + T cell phenotypes, Granzyme B as major cytotoxic CD8 + T cell molecule was chosen for analysis of T cell functionality and Ki67 as surrogate for T cell proliferation activity.Of note, the frequency and density of LOX1 + PMN-MDSC was increased in tumor islands of HPV − patients in comparison with the surrounding stroma, a phenomenon that was not observed in HPV + patients (Fig. 6A).Conversely, the prevalence of CD8 + T cells that express Granzyme B + and Ki67 + was notably greater within the tumor stroma of HPV + patients compared to HPV − patients (Fig. 6B).Consequently, the ratio of PMN-MDSC (CD66b + /LOX1 + ) and proliferating (Ki67 + ) CD8 + T cells as well as cytotoxic Granzyme B + /CD8 + T cells, was higher in HPV − patients (Fig. 6C) suggesting a higher level of intratumoral immunosuppression in those patients.Within the HPV + group, patients exhibiting low levels of intratumoral CD66b + / LOX1 + PMN-MDSC (falling below the median) demonstrated significantly improved survival rates (Fig. 6D).
Finally, we analyzed potential differences in the CD8 + T cell infiltrate of HPV + versus HPV − tumors in patients.We focused on CD8 + and CD45RO + cells as these markers are currently evaluated in the context of their immunoscore [33] and have been shown to be of prognostic relevance in HNC and other types of cancer [34][35][36].We found that CD8 + T cells were evenly distributed between the epithelial tumor islands (iTu, intratumoral) and the stromal regions (St, stromal) with the majority of patients showing less than 10% CD8 + T cells in these regions (score 1) (Fig. 7A).CD45RO + T cells were more frequent than CD8 + T cells and preferentially located in the stroma.Interestingly, the intratumoral region was more densely infiltrated by T cells in HPV + tumors compared to HPV − (Fig. 7B-D).When we combined the CD8 + T cell score and the CD45RO + T cell score in the intratumoral and stromal region into a comprehensive "immunoscore-like grade" (see material and methods for details), an increased infiltration of HPV + tumors with CD8 + and CD45RO + cells was confirmed (Fig. 7E).Collectively, these data show a denser T cell infiltrate in HPV + tumors, which is mainly caused by a stronger infiltration of the epithelial tumor islands (iT region).

Conclusion
It has been suggested that immune-mediated mechanisms contribute to the differential prognosis of patients with HPV + versus HPV − HNC.In accordance with this concept, several differences between the immune status of HPVrelated and non-related HNC have been described [6,27,37].Herein, we carried out the first detailed longitudinal analysis of T cells and MDSC in HPV-related HNC and compared pre-and post-therapy samples from 89/61 patients from whom material was available for this type of analysis.We found that pre-therapy levels of circulating monocytic MDSC were higher in HPV − disease.Importantly high numbers of MDSC correlated with low fractions of terminally differentiated CD8 effector T cells in the cohort of HPV − patients.HPV + tumors showed an increased infiltrate with CD8 + /CD45RO + T cells, Granzyme B + /CD8 + , and Ki-67 + /CD8 + .This increase in CD8 effector cells was accompanied by reduced infiltration of tumor tissue by PMN-MDSC.Therapy of HPV + patients was associated with an increase in CD8 + T EM in the peripheral blood, which was not evident in patients with HPV − disease.
In contrast to peripheral blood T cells, the clinical relevance of tumor-infiltrating CD8 + T cells is better established.A number of studies have associated high TIL counts with favorable prognosis in HNC [38][39][40].Controversial data exist, however, on a potential correlation between TIL frequency and HPV status.While in some studies such a correlation was found [27], it was missing in other studies [40].In our study we found a stronger infiltration of HPV + tumors with CD8 + and CD45RO + T cells.Considering the different biology of TIL in the tumor stroma versus tumor cell islands, we separately enumerated CD8 + T cells in stromal versus tumor regions and included markers originally used in the immunoscore [22].Combining both markers and intratumoral localizations, the HPV + tumors displayed a higher immune cell infiltrate.This was mainly caused by a stronger Fig. 6 PMN MDSC fractions and CD8 + T cell ratios differ in HPV + and HPV − patients and are associated with survival in the HPV + cohort.Analysis of 61 tumor tissue samples from patients with p16 + or p16 − oropharyngeal cancer (second cohort, see Table 2) by multi-color immunofluorescence.The frequency and density of CD66b +/ LOX1 + PMN-MDSC (A) and frequency of Granzyme B + cytotoxic or Ki67 + proliferating (B) CD8 + T cells was determined by digital pathology using tissue studio software modules.(C) shows PMN-MDSC to effector T cell ratios.The prognostic value of the frequency of PMN-MDSC was analyzed in HPV + and HPV − patients using Kaplan-Meier curves and overall survival (D).Data are presented as median dot plots.The Mann-Whitney rank sum test and the Kruskal-Wallis ANOVA on ranks [post hoc: Dunn's method] were used for statistical analysis.Survival analysis is depicted in Kaplan-Meier curves.* indicates a p value < 0.05 ◂ infiltration of the tumor cell islands in HPV + tumors.One could speculate that this deeper infiltration into the tumor center is associated with improved anti-tumor immune activity.However, in a study by Oguejiofor only CD8 + stromal and not CD8 + T cell infiltration of the tumor area was associated with improved clinical outcome [39].
Another novel aspect of our study is the longitudinal analysis which allowed us to examine therapy-associated changes.In a previous study of 20 patients with HPV + HNC, Parikh et al. found a decrease in HPV-peptide-induced IFN-γ in bulk cultures of PBMC when comparing pre-and posttherapy samples [41].While the authors did not determine the frequency of HPV-specific T cells in their system, they suggested that reduced IFN-γ was caused by immune inhibitory mechanisms at post-therapy time points and induced via chemoradiotherapy.Our data show a distinct reduction of HPV-specific T cells in post-therapy samples.This reduction of HPV-specific immunity may be caused by a reduction of tumor burden and viral load.Similar observations have been made for anti-HPV serum antibodies [42].Interestingly, in our study, this reduction in HPV-specific T cells was paralleled by an increase in the fraction of CD8 + T effector memory cells (T EM ) at post-therapy time points, which is indicative of a shift in the functional, anti-tumor effective T cell compartment of these patients.
Heusinkveld et al. reported that the frequency and functionality of HPV-specific T cells are higher in the tumor as compared to the peripheral blood [28].The exploitation of HPV-specific immunity may require therapeutic costimulation in order to achieve clinically meaningful responses [43].Alternatively, recent advances in ex vivo expansion of HPV-specific T cells may even make adoptive T cell transfer possible [44].We found a relatively low frequency (0.5-2 in 10.000 CD8 + T cells) of HPV-specific T cells in the peripheral blood of HPV + cancer patients.This frequency further declined after therapy.
Interaction of PD-1 with its ligand PD-L1 has been implicated in the regulation of the immune response in HNC and other types of cancer.Our data show a low expression of PD-1 on CD8 + peripheral blood T cells in HNC patients.Low expression of PD-1 on peripheral blood has also been reported by Lyford-Pike et al. and by Partlova et al. [27,45].This is in contrast to the higher expression of PD-1 on tumor-infiltrating lymphocytes noted in those studies.Although not statistically significant in the entire cohort, we noted a small increase in the percentage of PD-1 + CD8 + T cells in the course of therapy in some patients (Fig. 2F).Originally considered an exhaustion marker, under certain conditions, PD-1 expression may also represent antigenexperienced T cells with the unique capacity to recognize tumor-specific neoantigens [46,47].Thus, an increase in the course of therapy, as observed in our study, could be indicative of post-therapeutic immune activation.
Recent evidence suggests that MDSC significantly contribute to the regulation of anti-tumor immunity and in particular regulation of T cell activity in murine and human malignant disease.For the first time, we comprehensively compared levels of MDSC subsets in HPV + versus HPV − human HNC.Interestingly the frequency of M-MDSC was significantly higher in the peripheral blood of HPV − patients.In addition, in HPV − disease, high levels of MDSC were associated with lower fractions of terminally differentiated peripheral blood effector T cells.In murine models of vaccination against HPV + tumors, MDSC partially controlled the activity of anti-tumor CD8 + T cells [48].Some earlier studies suggested that PMN-MDSC induced antigen-specific tolerance, while monocytic MDSC is involved in non-specific suppression [49].Interestingly, in a recent multicenter analysis of circulating MDSC subsets across multiple diseases within the Mye-EUNITER (www.mye-eunit er.eu) consortium, programmed death-ligand 1 (PD-L1) was found to be primarily expressed in M-MDSC [50].As the PD-1 and PD-L1 interaction results in T cell anergy and inhibits T cell proliferation, the increased PD-L1 expressing M-MDSC fractions in HPV − HNC patients potentially contribute to our observed differences in T cell phenotypes with lower T EM cells and lower immunoscores for these patients.Accordingly, Spehner and colleagues identified M-MDSC levels in another HPV-dependent squamous cell carcinoma entity as the major influence on TH1-adaptive immune responses and patients' survival [51].
Additionally, our investigations show that differences between HPV + and HPV − tumor patients are also present in the tumor tissue.In previous work we provided evidence for in situ intratumoral activity by LOX1 + PMN-MDSC [14].These findings are in line with reports by Chai and colleagues who identified ARG1/iNOS-dependent inhibition of T cell proliferation by LOX1 + PMN-MDSC in glioblastoma [52].In this study we report reduced Ki67 and Granzyme B expression in CD8 + T cells of HPV − patients that also have higher PMN-MDSC / CD8 + Granzyme B + ratios.In HPV + patients, we observed improved clinical outcome in patients with lower PMN-MDSC fractions.An association of MDSC and T cell frequencies was also observed by Liang and colleagues who reported-in cervical cancer patients-a negative correlation of PMN-MDSC and densities of CD8 + T cells, influencing clinicopathologic parameters and prognosis [53].Although the immunoregulatory mechanisms of MDSC in cancer patients are still far from being understood, these data support a link between MDSC expansion and T cell differentiation and function, which could further explain the somewhat poor prognosis of patients with HPV − disease.
In sum, our analyses revealed important differences in the MDSC and T cell compartment in HPV + versus HPV − disease.Patients with HPV + HNC showed a post-therapeutic expansion of CD8 + T EM .These patients also demonstrated higher immune cell infiltration of their tumors as defined by the immunoscore markers CD8 and CD45RO.In contrast, patients with HPV − tumors had higher pre-treatment levels of M-MDSC and this correlated with a reduced frequency of CD8 + T effector cells in the blood.HPV − patients also showed higher frequencies of intratumoral PMN-MDSC and lower frequencies of proliferating and cytotoxic T cells.In HPV + patients, lower intratumoral PMN-MDSC levels were associated with an improved survival.These differences in MDSC and T cell functionality may further contribute to and explain differences in clinical outcome of patients with HPV + versus HPV − HNC.

Fig. 1
Fig. 1 Differences in the CD8 + T cell compartment between HNC patients and healthy donors.PBMC from HNC patients and healthy donors were isolated by density gradient centrifugation, stored frozen until use, and analyzed by flow cytometry.Data are depicted as the relative frequency of CD8 + in CD3 + (A), CCR7 + in CD8 + (B), and CD62L + in CD8 (C).In (D) expression of activation and exhaustion markers PD-1 and TIM-3 on CD8 cells was quantified by flow cytometry (DP = double positive, DN = double negative).(E) Depicts the gating strategy for CD8 + T cell phenotype subpopulations and representative examples of a HNC patient and a healthy donor.Gating was

Fig. 4
Fig. 4 Analysis of CD8 + T cells in HPV16-driven versus non-HPVdriven HNC.HPV16 RNA positivity was determined by immunohistochemistry for the surrogate markers pRb and p16.Whole tissue slides (FFPE) from each patient were stained for pRb and p16 with representative examples shown in (A).High p16 and low pRb indicate RNA positivity and HPV16-driven HNC.Low p16 and high pRb indicate RNA negativity.The relative frequency of CD8 + cells (B) and CCR7 + /CD8 + cells (C) was determined in pre-therapy and post-therapy (surgery, surgery + R(C)Tx, primary R(C)Tx) samples by flow cytometry.CD8 + T cell subsets (D) and expression of PD-1

Fig. 5
Fig. 5 Peripheral blood MDSC frequency is linked to HPV status and correlates with T cell phenotype.MDSC analysis was performed in pre-therapy patient and healthy donor samples by flow cytometry.Polymorphonuclear MDSC (PMN-MDSC), monocytic MDSC (M-MDSC), and immature or early-stage MDSC (e-MDSC) were determined in HPV − , HPV + HNC patients, and healthy donors.Note the increased M-MDSC frequency in HPV − patients (A).The gat-

Fig. 7 T
Fig.7T cell infiltration in the tumor microenvironment differs between HPV − and HPV + patients.Tumor microarrays (including several samples for each patient) were immunohistochemically stained for CD8 and CD45RO, respectively.Epithelial tumor cell nests (iT) and intratumoral stromal (St) regions were scored independently.A score based on the approximate percentage of stained cells (1: 0-10%, 2:10-50%, 3:50-80%, and 4:80-100%) was assigned to each region by three blinded, independent, and trained investigators.Score distribution was compared between all patients, HPV − and

Table 1
Clinical information for patients in cohort 1

Table 2
Clinical information for patients in cohort 2