- Research article
- Open Access
High numbers of activated helper T cells are associated with better clinical outcome in early stage vulvar cancer, irrespective of HPV or p53 status
Journal for ImmunoTherapy of Cancervolume 7, Article number: 236 (2019)
Vulvar squamous cell carcinoma (VSCC) has been suggested to consist of three subtypes; HPV-positive, HPV-negative mutated TP53 or HPV-negative TP53 wildtype, with different clinical courses. To analyze the immune infiltrate in these molecular subtypes and its impact on clinical outcome, an in-depth study of the tumor immune microenvironment was performed.
Sixty-five patients with invasive VSCC matched for age, FIGO stage and treatment modality, were grouped according to the presence of HPV and p53 protein expression status. Archived tissues were analyzed for intraepithelial and stromal expression of CD3, CD8, Foxp3, PD-1, and pan-keratin in randomly selected areas using immunofluorescence. Additional phenotyping of T cells was performed ex-vivo on VSCC (n = 14) and blood samples by flow cytometry. Healthy vulvar samples and blood served as controls.
Based on T-cell infiltration patterns about half of the VSCC were classified as inflamed or altered-excluded while one-third was immune-deserted. High intraepithelial helper T cell infiltration was observed in 78% of the HPV-induced VSCC, 60% of the HPVnegVSCC/p53wildtype and 40% of the HPVnegVSCC with abnormal p53 expression. A high intraepithelial infiltration with activated (CD3+PD-1+), specifically helper T cells (CD3+CD8−Foxp3−), was associated with a longer recurrence-free period and overall survival, irrespective of HPV and p53 status. Flow cytometry confirmed the tumor-specific presence of activated (CD4+PD-1++CD161−CD38+HLA-DR+ and CD8+CD103+CD161−NKG2A+/−PD1++CD38++HLA-DR+) effector memory T cells.
This is the first study demonstrating an association between intraepithelial T cells and clinical outcome in VSCC. Our data suggest that abnormal p53 expressing VSCCs mostly are cold tumors whereas HPV-driven VSCCs are strongly T-cell infiltrated.
Immunotherapy of cancer has established itself as a new breakthrough approach that offers long-term durable clinical responses in patients with advanced cancers. As the initiation and regulation of the immune response to tumors is complex and multistep in nature, inspection of the different processes involved is required to provide the optimal (combination) of immunotherapeutic modalities available . This is highly relevant for vulvar squamous cell carcinoma (VSCC), the predominant histologic subtype of vulvar cancers, for which new treatment options are urgently needed. Because despite current treatment, consisting of radical surgery and/or (chemo) radiotherapy which causes impressive morbidity, lymphedema, sexual and psychological dysfunction and wound healing disorders [2, 3], 46% of VSCC patients still develop recurrences after 10-years .
At this point however, little is known about the role and impact of cellular immunity on the clinical outcome of VSCC. Both CD4 and CD8 T cells as well as B cells infiltrate VSCC [5,6,7]. The CD4 cells comprise CD4+ helper T cells and regulatory T cells (Tregs). Often a strong infiltration with one type of T cells is paralleled by the others [5, 6, 8]. In three studies focussing on the prognostic role of CD4+ and/or CD8+ T cells or Tregs no impact on clinical outcome was found [6, 9, 10]. On the one hand, these analyses may have been influenced by the heterogeneity of the study group with respect to tumor etiology, stage and treatment. Furthermore, enumeration of all T cells, irrespective of their location in the tumor , as well as preselection of highly infiltrated areas only [6, 10], may also have influenced study outcomes. On the other hand, the impact of T cells may be nullified by the presence of immune regulatory mechanisms, as VSCC can be massively infiltrated with M2 macrophages and Tregs . Moreover, VSCC can express the immunoregulatory enzyme, indoleamine 2,3-dioxygenase (IDO) or PD-L1, both of which were shown to negatively influence clinical outcome [10, 11]. Notably, PD-L1 was mainly expressed in lymphocyte rich areas , suggesting that it functioned as an adaptive escape mechanism , and implying that in some VSCC a functionally active antitumor response is present. This notion is sustained by the observation that the intraepithelial presence of Granzyme B-positive cells is related to better overall survival (OS) in patients with localized VSCC .
At present, three distinct etiologic pathways in the development of VSCC have been proposed. One type is driven by the overexpression of high-risk human papilloma virus oncogenes E6 and E7 (HPVposVSCC). The second type is not related to HPV and can be categorized based on the mutational status of the tumor suppressor gene TP53 associated with high protein levels of p53 (HPVnegVSCC/p53abn). We recently identified a third type as a substantial group of patients with a HPV-negative VSCC displaying normal expression levels of p53 protein (HPVnegVSCC/p53wt) but frequently bearing other mutations . Importantly, HPV-driven VSCC display better OS and a longer recurrence-free period (RFP) than HPVnegVSCC [14,15,16,17]. Interestingly among the latter group, local recurrences more often occurred after treatment in patients with HPVnegVSCC/p53abn than in HPVnegVSCC/p53wt . With the first reports showing an influence of different oncogenic pathways on local immunity [18, 19], we asked the question if the differences in RFP and OS observed between the three groups of VSCC driven by different oncogenic pathways could be explained by the local immune response. Bearing in mind the limitations of previous studies, we selected three cohorts of VSCCs based on their HPV and p53 protein (abn/wt) status which were highly matched for clinicopathological parameters and enumerated different types of intraepithelial and stromal T cells in randomly selected fields of VSCC, using multiplex immunofluorescence. In-depth analysis of T cells was performed on freshly dispersed tissue by flow cytometry. Our study revealed a strong impact of intraepithelial activated T cells on clinical outcome, in particular a dense infiltration with intraepithelial CD4+ T cells was highly associated with RFP and OS, irrespective of HPV or p53 status. Moreover, the percentage of tumors highly infiltrated with these T cells varied between the three different subtypes, with HPV-induced VSCC most often strongly infiltrated (78%) followed by the HPVnegVSCC/p53wt (60%) and the lowest infiltration in the HPVnegVSCC/p53abn group (40%).
Material and methods
Archived formalin-fixed paraffin-embedded (FFPE) tumor tissue from VSCC patients was selected from a larger cohort with known HPV and p53 status. HPV presence was tested by HPV-PCR and p16 IHC . Tumors that were positive in both tests were assigned as HPVposVSCC. When both tests were negative, tumors were scored as HPVnegVSCC. The HPVnegVSCC were further sub-classified based on the wildtype or abnormal expression of p53 (HPVnegVSCC/p53wt and HPVnegVSCC/p53abn) as previously described . In addition, archived FFPE healthy HPV-negative vulvar tissue from 10 women who underwent labial reduction surgery served as controls. Fresh tumor tissue (n = 14) and blood samples (n = 34) were obtained from 38 patients participating in the large observational CIRCLE study. Women with histologically proven p16ink4a-negative VSCC were included in this study investigating cellular immunity against anogenital lesions [21, 22]. Tumor staging was done according to FIGO staging 2009. Patients were included after signing informed consent. The study was conducted in accordance with the Declaration of Helsinki and approved by the local medical ethical committee of the Leiden University Medical Center (P08.197 and B16.024) and in agreement with the Dutch law. The materials were used according to the Dutch Federation of Medical Research Association guidelines. The patients received standard-of-care treatment consisting of primary surgery.
Multiplex six color staining, image acquisition and analysis by VECTRA
For the direct and indirect immunofluorescent six color staining and detection, 4 μm FFPE tissue sections were deparaffinized and prepared with heat-induced antigen citrate (10 mM, pH 6.0) retrieval as described previously . Antibody specificity and optimal conditions for antigen retrieval were assessed by single-plex IHC using tonsils as a positive control . After incubation with superblock buffer (Thermo Fisher Scientific, Waltham, MA, USA), the primary antibodies and isotype/species-specific secondary fluorescent antibodies were applied (Additional file 1). Nuclear counterstain was obtained with DAPI. Tissue slides were imaged at 20x magnification with the Vectra 3.0 Automated Quantitative Pathology Imaging System (Perkin Elmer). Imaging analysis and spectral separation of dyes was performed with the InForm Cell Analysis software (Perkin Elmer) by using spectral libraries defined with single-marker immunofluorescence detection. Five random multispectral imaging fields of interest were selected for acquisition from each tumor, depending on its size. Tissue and cellular segmentation was done as described before . The following phenotypes were identified for the T cell panel: total T cells (CD3+), CD8+ T cells (CD3+CD8+Foxp3−), helper T cells (CD3+CD8−Foxp3−), Tregs (CD3+CD8−Foxp3+), PD-1 expressing T cells (CD3+PD1+). All images were visually inspected to confirm the correct attribution and quantification of phenotypes, and segmentation of tissue. Because PD-1 could be expressed by CD3+CD8− and CD3+CD8+ cells, the CD3+PD1+ phenotypes were separately analyzed. All phenotypes in both areas were normalized by tissue area (number of cells/mm2). In addition, ten HPVposVSCC, six HPVnegVSCC/p53wt and five HPVnegVSCC/p53abn VSCC samples were used to study Tbet (Santa Cruz, clone H-210, dilution 1:100) expressing CD3+ cells with immunofluorescence.
Blood and tumor cell isolation and culturing
Venous blood samples were drawn prior to surgery, and peripheral blood mononuclear cells (PBMC) were isolated using Ficoll density gradient centrifugation as described previously [24, 25]. VSCC tumor material was obtained and handled as described [24, 25]. First, tumor material was cut into small pieces. One-third of the tumor pieces was incubated for 60 min at 37 °C in Iscove’s Modified Dulbecco’s Medium (IMDM, Gibco by life technologies, ThermoFisher Scientific, Lonza, Verviers, Belgium) with 10% human AB serum (Capricorn Scientific, Esdorfergrund, Germany) and supplemented with high dose of antibiotics (50 μg/ml Gentamycin (Gibco/ Thermo Fisher Scientific (TFS), Bleiswijk, the Netherlands), 25 μg/ml Fungizone (Gibco/Thermo Fisher Scientific), after which the tumor pieces were put in culture in IMDM supplemented with 10% human AB serum (IMDM complete) and 1000 IU/ml human recombinant IL-2 (Aldesleukin, Novartis, Arnhem, the Netherlands). Cultures (n = 14) were replenished every 2–3 days with fresh IMDM complete and IL-2 to a final concentration of 1000 IU/ml. After 2–4 weeks, when sufficient T cells were obtained, the cells were cryopreserved and stored in liquid nitrogen until use. Approximately two-third of the tumor pieces was incubated for 15 min at 37 °C in IMDM dissociation mixture containing 10% human AB serum, high dose of antibiotics (as above) and 0.38 mg/ml of the commercially available Liberase enzymes (Liberase TL, research grade, Roche). Following incubation, cell suspension was put on a 70 μm cell strainer (Falcon, Durham, NC, USA) to obtain a single cell suspension, counted using trypan blue exclusion (Sigma, St Louis, MO, USA), and cryopreserved at approximately 2 million cells/vial. All cells were stored in the vapor phase of liquid nitrogen until further use.
Flow cytometry and data analysis
Cryopreserved PBMC (n = 34) and/or cryopreserved freshly isolated single cell tumor samples (n = 12) were thawed and assessed by flow cytometry as described before [26, 27]. In brief, samples were thawed according to standard operation procedures and stained with the LIVE-DEAD® Fixable yellow dead cell stain kit (ThermoFisher Scientific) for 20 min at room temperature to identify dead cells. Following incubation, the cells were washed, incubated with PBS/0.5%BSA/10%FCS for 10 min on ice to block Fc receptors. After the cells were washed, the cells were stained for 30 min on ice and in the dark with fluorochrome-conjugated antibodies. Intracytoplasmic/intranuclear staining was conducted with the BD Pharmingen Transcription Factor Buffer set (BD Biosciences) according to manufacturers’ protocol. Details on antibodies used are listed in Additional file 1. Acquisition of cells was done on a BD LSR Fortessa. Data was analyzed by either manual gating using DIVA software (version 8.02; BD Biosciences) or by high-dimensional single cell data analysis using hierarchical Stochastical Neighbor Embedding (HSNE)  in Cytosplore. To automatically discover stratifying biological signatures at the single cell level, we used the fully automated hierarchical clustering (unsupervised) tool CITRUS in the cloud-based cytobank software (Fluidigm Sciences) with an FDR of 1%.
Cytokine production of phytohemagglutinin (PHA)-stimulated TIL
To determine capacity of tumor infiltrating lymphocyte (TIL) batches from HPVnegVSCC tumors to produce cytokines in response to mitogenic stimulation, cultured TIL batches (n = 14) were stimulated with 0.5 μg/ml PHA (HA16 Remel; ThermoFischer Scientific) for 4 days, after which supernatants were harvested and analyzed by Cytometric Bead Array (CBA, Th1/Th2 kit, BD Bioscience, Breda, the Netherlands) according to the manufacturer’s instructions. The cut-off value for cytokine production was 20 pg/ml, except for IFN-γ for which it was 100 pg/ml. Positive cytokine production was defined as at least twice above that of the unstimulated cells [25, 29].
For data analysis the statistical software package SPSS 23.0 (SPSS Inc., Chicago, IL) was used. Group comparisons of categorical data were performed by chi-square test. The non-parametric Mann-Whitney U test was used for continuous variables when comparing two groups. For the survival analysis, patients were categorized into two groups based on numerical immune cell count. First, the median cell count was used as cut-off value. To optimize the chance to detect a relationship between T-cell subsets and clinical outcome in a relatively small group of patients, the best cut-off value for the different T cell subsets was determined using receiver operating characteristics (ROC) curve analysis. The T-cell subset values with the best accuracy (i.e. with greatest sensitivity and specificity) were selected as the most optimal cut-off value for (OS or RFP). Based on these cut-off values, the immune cell counts were categorized in two groups and a log-rank test was performed to calculate the difference in OS or RFP. The RFP was censored for lost-to-follow up and death. Two sided p-values < 0.05 were considered significant. GraphPad Prism 7 (GraphPad Software Inc., LA Jolla, CA, USA) was used to illustrate the data by graphs and figures.
A cohort of 65 primary VSCC samples, divided in HPVposVSCC (n = 23), HPVnegVSCC/p53wt (n = 20), and HPVnegVSCC/p53abn (n = 22) was analyzed. All cases were matched for age (40–85 years), FIGO stage, absence of lymph node and distant metastasis, ≥8 mm tumor-free margins, no use of immunosuppression, and no documented medical history. However, HPVposVSCC were younger than the other groups despite matching because younger women are more likely to have HPVposVSCC than HPVnegVSCC . An overview of patient characteristics and treatment is given in Additional file 2. In line with current literature [15, 17, 31], the group of patients with HPVposVSCC displayed a better OS and RFP than those with HPVnegVSCC (Additional file 3). Furthermore, the recurrence rate increased from HPVposVSCC (13%), HPVnegVSCC/p53wt (40%) to 59% in HPVnegVSCC/p53abn (Additional file 2). Together this confirms our selection of a representative cohort of patients for our study.
The HPVposVSCCs are most often strongly infiltrated with T cells
The archived tissues sections were simultaneously analyzed for the expression of CD3, CD8, Foxp3, PD-1, and pan-keratin by multispectral immunofluorescence VECTRA analysis, both in the epithelial and stromal compartments (Additional file 4). Quantification of the T cells per square mm of tumor revealed that the stroma of VSCC was more densely infiltrated with CD3+ T cells, CD3+CD8−Foxp3− T cells, CD3+CD8−Foxp3+ Tregs, and CD3+CD8+Foxp3− T cells than healthy controls. The number of intraepithelial Tregs was also higher in VSCC (Fig. 1; Additional file 4; Additional file 5). Comparison of the three subgroups revealed a strong difference in T cell infiltration between HPVposVSCC and HPVnegVSCC/p53abn. The majority of HPVposVSCC was well infiltrated whereas the HPVnegVSCC/p53abn most often displayed a low T-cell infiltration. The group of HPVnegVSCC/p53wt showed a more variable pattern, with low and high T-cell infiltrated tumors (Fig. 1; Additional file 4; Additional file 5). The number of tumor-infiltrating intraepithelial cells was highly correlated to the other intraepithelial T-cell subsets and to their numbers in the stroma (Additional file 6). This suggests a coordinated response of CD3+CD8−Foxp3− and CD3+CD8+Foxp3− T cells in VSCC. Quantification of intraepithelial CD3+Tbet+ T cells, representing type 1 immunity , revealed higher numbers in HPVposVSCC compared to both HPVnegVSCC subtypes (Fig. 1; Additional file 5).
Immune inflamed, altered-excluded, altered-immunosuppressed and deserted VSCC
Based on the previously published categories of T-cell infiltration patterns , the VSCC were characterized (Fig. 2) as immune-deserted (n = 19), −altered (n = 41) or -inflamed (n = 5). The immune-altered group was the largest and could be subdivided based on two distinct patterns of T cells in the stroma . The altered-excluded tumors (n = 24) showed more stromal T cells at the invasive border whereas in the altered-immunosuppressed VSCC (n = 17) T cells were dispersed throughout the whole stroma (Fig. 2). Notably, the number of CD3+ T cells at the invasive border (Fig. 2c) was highly correlated with the intraepithelial CD3+ T-cell count (p = 0.000; Additional file 6) in the altered-excluded VSCC. Moreover, the average number of intraepithelial CD3+ T cells in the altered-excluded was higher than in the altered-immune suppressed (mean 612 ± SD 539 vs mean 157 ± SD 92, p < 0.001, respectively). To evaluate the impact of these four VSCC categories on survival, a Kaplan-Meier analysis was performed. The immune-inflamed group displayed a superior RFP and OS (Additional file 7). Interestingly, the group of altered-excluded VSCC showed a similar OS whereas the RFP was less good when compared to that of the immune-inflamed group. Therefore, the immune inflamed and altered-excluded VSCC were classified as hot tumors.
The intratumoral CD3+CD8−Foxp3− T cell count is an independent prognostic marker for RFP and OS irrespective of VSCC type
The better OS of the two categories displaying stronger intraepithelial infiltration than altered-immune suppressed and immune-desert VSCC suggested an important role for intraepithelial T cells on clinical outcome. For each T-cell subset the median cell count (Additional file 8) and the optimized cut-off point, as determined by ROC curve analysis, was used to categorize the patients tumor into low or high infiltrated and subsequently its impact on clinical outcome was determined. High intraepithelial infiltration with CD3+ T cells, in particular of CD3+CD8−Foxp3− T cells was strongly associated with longer RFP (p = 0.001) and OS (p = 0.004). A strong infiltration with CD3+PD1+ T cells was also associated with a longer RFP (p = 0.032). The intraepithelial infiltration with CD3+CD8+ T cells or CD3+CD8−Foxp3+ Tregs was not prognostic for clinical outcome (Fig. 3; Additional file 9). Importantly, the prognostic power of CD3+CD8−Foxp3− T cells for RFP was retained when only the HPVnegVSCC were analyzed (Additional file 9). To sustain this notion, the impact of tumor infiltrating CD3+CD8−Foxp3− T cells in clinical outcome was corrected for age, and p53 and HPV status (Additional file 10). In the univariate analysis, only high CD3+CD8−Foxp3− counts and age correlated with RFP. In the multivariate analysis, high infiltration with CD3+CD8−Foxp3− T cells but not age was associated with longer RFP (HR 3.30 (1.22–8.94), p = 0.018). Thus, the CD3+CD8−Foxp3− T-cell infiltration is expected to be an important prognostic marker for clinical outcome, irrespective of whether these VSCC are caused by the HPV-derived oncogenes or other oncogenic pathways (e.g. p53 mutation).
HPVnegVSCC are infiltrated by activated CD8+ and CD4+ effector memory T cells
The vast majority (~ 80%) of vulvar cancers are not induced by HPV.  While 78% (18/23) of HPVposVSCC were strongly infiltrated with CD3+CD8−Foxp3− T cells, there was also a considerable fraction of HPVnegVSCC displaying evidence of their immunogenicity with 60% (12/20) of the HPVnegVSCC/p53wt and 40% (9/22) of the HPVnegVSCC/p53abn showing high intraepithelial CD3+CD8−Foxp3− T-cell infiltration and longer RFP and OS. In order to gain a better understanding of these tumor-infiltrating T cells in HPVnegVSCC, a serie of fresh HPVnegVSCC tumor biopsies was used to culture tumor-infiltrating lymphocytes (TIL; n = 14) and for ex-vivo phenotypic analysis of freshly dissociated and directly liquid nitrogen stored tumor-infiltrating T cells by flow cytometry (n = 12). Upon mitogenic stimulation, the growing TILs predominantly produced the type I cytokine IFN-γ and the type 2 cytokine IL-5 at on average similar concentrations, suggesting the presence of both Th1 and Th2 cells in these tumors (Fig. 4). Only low concentrations of TNF-α, IL-4 and IL-10 were found.
To analyze the tumor-infiltrating T cells an antibody mix against CD45, CD3, CD4, CD8, CCR7, CD45RA, CD103, CD161, PD-1, CD38, HLA-DR and NKG2A was used to stain the fresh VSCC digests. In addition, PBMC of healthy female controls (n = 11) and PBMC of HPVnegVSCC (n = 29) were stained. A combined hierarchical Stochastical Neighbor Embedding (HSNE) analysis of the high-dimensional single cell data obtained from staining the blood and tumor samples resulted in the identification of several distinct immune populations (clusters), which were more prominently present or absent in the tumors or PBMC of the VSCC patients (Fig. 5a). In order to automatically discover stratifying biological signatures within VSCC blood and tumor samples, we made use of the automated and data-driven CITRUS platform, as an unbiased and thorough correlation-based tool for mining and inspection of cell subsets at the single cell level nested within high-dimensional datasets . The CITRUS analysis resulted in ten distinctive (groups of) lymphocyte populations significantly higher present in one or more of the three different types of samples (Fig. 5b). The CD4 and CD8 distribution within the total CD3+ T-cell population did not differ between PBMC and tumors of VSCC patients (Fig. 5c). Of the five identified CD8+ T-cell subsets, populations 2, 4, and 5 were significantly underrepresented in VSCC tumors (Fig. 5d). Population 2 comprised CD8+CD161+PD1+CD38+ HLA-DR− effector memory RA+ T cells (Temra), whereas population 4 (CD8+ Tcm/em) and population 5 (CD8+ Temra) did not display these markers (Additional file 11). The CD8+ T-cell population (#1) which was clearly overrepresented in VSCC tumors, consisted of CD8+CD103+CD161−NKG2A+/−PD1++CD38++HLA-DR+ Tem cells. Of the five different CD4+ T-cell subsets identified, population 6 (CD4+PD-1−CD161−CD38++HLA-DR− naïve T cells) and population 10 (CD4+PD-1+CD38−CD161−HLA-DR− Tem/cm) were lower in VSCC than PBMC. In contrast, two populations of CD4+ effector T cells were found at significantly higher percentages and comprised CD4+PD-1++CD161−CD38+HLA-DR+Tem (#7) as well as CD4+PD-1−CD161−CD38−HLA-DR−Tcm/em (#9). The co-expression of PD-1, CD38 and HLA-DR is indicative for T-cell activation. As such the strong tumor-specific infiltration of HPVnegVSCC with activated CD8+ (population 1) and CD4+ (population 7) effector T cells sustains the notion these tumors are immunogenic and explain why their presence is associated with better clinical outcome.
For half (6/12) of the freshly digested VSCC samples enough material was available to characterize the T cell infiltrate with a second antibody mix against CD3, CD4, CD8, CD25, CD127, Foxp3, Tim-3, Lag-3 and Tbet. These samples were analyzed for the presence of different types of Tregs, Tbet+ cells and the two checkpoint molecules according to the strategy shown in Additional file 12. Similar to what was found in the FFPE tissue samples, a tumor-specific increase in activated and proliferating (Ki67+) Tregs was observed (Fig. 6a). Furthermore, tumor-specific increases in the percentages of Tim-3 and Lag-3 Tregs, CD8+ and non-Treg CD4+ T cells were observed (Fig. 6b), confirming that part of the tumor-infiltrating T cells has been activated. Last but not least, on average 30% of the CD8+ and 20% of the non-Treg CD4+ T cells expressed the transcription factor Tbet, which is in line with the IFN-γ production of the cultured TILs. Finally, only a small percentage of the Tregs expressed Tbet (Fig. 6b).
Based on the cytokine production and the expression of several checkpoints, transcription factors and activation markers, we conclude that HPVnegVSCC are infiltrated with variable numbers of activated type 1 and 2 CD8+ and CD4+ effector T cells as well as Tregs.
We asked the question if the observed differences in recurrence rate and survival between the three subtypes of VSCC [14, 17], classified by the presence of HPV, the overexpression of p53 or the absence of both, may have an immunological background. Our study is the first to show that a strong infiltration of the tumor cell nests with helper (CD3+CD8−Foxp3−) T cells is important for clinical outcome after primary surgery, irrespective of whether the VSCC are caused by HPV or other oncogenic pathways, including TP53 mutations. Probable reasons for the failure to detect this association in previous studies [6, 9,10,11] are related to the importance of the location of the T cells in the tumor and the homogeneity in stage and treatment of the VSCC patients analyzed. In line with the RFP and the percentage of recurrences found in each of the three subtypes in VSCC, the percentage of tumors with high intraepithelial helper T-cell infiltration was the highest in the HPV-driven VSCC (78%), followed by VSCC not associated with HPV or p53 overexpression (60%), and the lowest in VSCC with abnormal p53 expression (40%). Importantly, these data suggest that T-cell infiltration of VSCC may be influenced by the oncogenic pathway underlying the development of a particular VSCC.
Although the helper T cells display the strongest relationship with clinical outcome, this does not mean that CD8+ T cells are not important in VSCC. Whenever there is a strong intraepithelial infiltration with T cells, this is because both subsets of T cells are present in large numbers. Furthermore, we noted a positive association between the intraepithelial presence of CD3+PD-1+ T cells and clinical outcome. In-depth flow cytometric analysis revealed that this PD1+ T-cell population comprised activated CD4+PD-1++CD161−CD38+HLA-DR+ and CD8+CD103+CD161−NKG2A+/−PD1++CD38++HLA-DR+ effector memory T cells. Potentially, helper T cells play an important role in VSCC because a substantial fraction of VSCC can partially downregulate HLA class I expression whilst the levels of tumor-expressed HLA class II may go up . Based on the percentages of T cells co-expressing the transcription factor Tbet, as counted in the tumor sections and measured in fresh VSCC by flow cytometry, and by the detection of IFN-γ and IL-5 in the supernatants of stimulated TIL, the VSCC-infiltrating T cells are both of a type 1 and 2 phenotype.
Current treatment of VSCC does not take into account the differences in etiology and clinical outcome . The relationship between T-cell infiltration and clinical outcome suggests that immunotherapy may form a new treatment option for VSCC, as in other tumor types this was associated with a better response to immunotherapy [32, 33]. These other cancer types were categorized into immune-inflamed (hot), − altered (excluded or suppressed) and -deserted (cold) tumors in order to define which immunotherapeutic (combination) approach may work best. For instance, hot tumors show the best response to checkpoint blockade (e.g. anti-PD1 and anti-CTLA-4) [32, 33]. Also VSCC could be divided according to these four immune phenotypes. Only a few (5 of 65) VSCC were categorized as truly immune-inflamed while a substantial portion (37%, n = 24) displayed the immune altered-excluded phenotype. However, the latter group displayed a significant stronger intraepithelial T cell infiltration when compared to the immune altered-suppressed phenotype and showed a better OS. Patients with inflamed and immune altered-excluded VSCC may be selected for treatment with immunotherapy. In our study high percentages of intratumoral T cells expressed PD-1. Others found variable percentages of cases in which the VSCC (> 30%) or VSCC-infiltrating immune cells (> 90%) expressed PD-L1 [11, 36, 37], congruent with our observation that there are varying numbers of tumor-infiltrating T cells which can produce IFN-γ, as indicated by expression of Tbet, and may lead to adaptive PD-L1 expression . Altogether, this makes a strong case for the treatment of inflamed and immune altered-excluded VSCC with PD-1/PD-L1 checkpoint therapy. Indeed, a first successfully treated case with advanced stage recurrent vulvar cancer has been reported with PD-L1 blockade . A tumor-specific increase of the CD4+ T-cell response seems more likely to be achieved by CTLA4-blockade than by targeting PD-1 , arguing for a combination of PD-L1 and CTLA-4 blockade to reinvigorate the tumor-specific CD4+ T-cell response. Another option would be the use of an agonistic antibody to OX-40 , which in combination with PD-L1 blockade displayed synergistic effects on CD4+ T-cell reactivity . Moreover, treatment of the generally T-cell infiltrated HPVposVSCC may include therapeutic HPV16 vaccination as HPV-driven oropharyngeal cancers responded well to the combination of checkpoint therapy and therapeutic vaccination .
One-third of our VSCC were phenotyped as immune-deserted or cold tumors which may exist because of a lack of antigens or their presentation (immune ignorance), or because of several deficits leading to a lack of priming or to tolerance [32, 42]. Apart from the HPV-induced VSCC, the viral proteins of which may drive a strong T-cell response similar to what is seen in HPV-driven oropharyngeal cancer , 40–60% of the HPVnegVSCC display strong intraepithelial T cell infiltration. This suggests that also in these tumors immunogenic tumor antigens are expressed and presented. Currently, the antigens recognized by T cells in HPVnegVSCC are unknown but the majority of primary VSCC express for instance the well-known tumor antigens MAGEA1 and MAGEA4 , but it is unknown if these antigens function as target for the VSCC-infiltrating T-cells as this still needs to be studied. We have no data on a direct relation between the level of T-cell infiltration and the availability of strong tumor antigens expressed by VSCC, but our work on HPV-specific immunity in different cancers , as well as that of others in for instance melanoma [44,45,46], suggest that this is not the case. More likely, a lack of inflammation or danger signals has played a role in HPVnegVSCC. Intratumoral activation of the Stimulator of Interferon Genes (STING) pathway , the use of oncolytic viruses , but also intratumoral injections of toll-like receptor (TLR)-agonists  have been shown to sensitize cold tumors to checkpoint blockade. For cold VSCC tumors, the TLR7/8-agonist imiquimod may be a promising topically applied therapeutic agent. Imiquimod elevates numerous genes involved in the regulation of innate immunity, resulting in the migration of DC to the application site, and subsequently the activation of a type 1 T-cell response . Patients with a precancerous lesion of HPVposVSCC responded very well to imiquimod therapy . Notably, imiquimod treatment of breast metastases in the skin not only converted them from cold to hot, as demonstrated by a profound infiltration with CD4+ and CD8+ T cells, but also led to tumor regression .
There are several limitations to our study. The correlation between better clinical outcome and intraepithelial T-cell infiltration was found in a highly homogeneous patient group with early stage cancer and treated with surgery. If this relationship also exists in locally advanced cancer patients treated with (chemo) radiotherapy has to be determined. Furthermore, next to the median cell count, we optimized the chance to detect a statistically significant relation between T-cell subsets and clinical outcome in a relatively small group of patients. Hence our results need to be validated in a larger cohort. In addition, our data suggests that the etiology of the VSCC may have an impact on its immunogenicity. While this would fit with the concept that different oncogenic pathways may influence local immunity [18, 19], the numbers of VSCC analyzed are such that the outcome can only be used for hypothesis generation. Moreover, we have not analyzed the myeloid cell component, which on itself may impact prognosis and T-cell function. Finally, less than 20% of VSCC are induced by HPV. As such, we merely collected small pieces of fresh tumor tissue and PBMCs from HPVnegVSCC patients. A more extensive comparison of the tumor-infiltrating immune phenotypes between HPV-induced and HPVnegVSCC, therefore, was not possible.
In conclusion, our observation that a strong coordinated intraepithelial infiltration with T cells is highly associated with a better clinical course of early stage VSCC after surgery and suggests that this group of patients may benefit from immunotherapy as an alternative to potential mutilating surgery in this delicate anatomical area . In parallel to the use of the different categories of immune cell infiltrated tumors in other types of cancer [32, 33], these tumor classifications can also be used to tailor immunotherapy approaches in VSCC. Near future studies should focus on the effects of checkpoint blockade in patients classified with inflamed and altered-excluded tumors while patients diagnosed with an immune-deserted (cold) VSCC may benefit more from therapies that induce acute inflammation. Furthermore, future studies should assess the mutational landscape of VSCC as this will reveal if a tumor-specific (neo)-antigen T-cell repertoire could be harnessed to treat the HPVnegVSCC as well as assess if the HPV-specific T-cell repertoire can be exploited to treat HPVposVSCC.
Availability of data and materials
All data generated or analyzed during this study are included in this published article (and its additional files) and are available from the corresponding author on reasonable request.
Anti-cytotoxic T-lymphocyte-associated protein 4
Anti-programmed cell death protein 1
Bovine serum albumin
Cytotoxic T-lymphocyte-associated protein 4
Fetal calf serum
Formalin-fixed paraffin embedded
Human leukocyte antigen
Human papilloma virus
Human papilloma virus negative vulvar squamous cell carcinoma
Human papilloma virus negative vulvar squamous cell carcinoma with p53 abnormal expression
Human papilloma virus negative vulvar squamous cell carcinoma with p53 wildtype expression
Human papilloma virus positive vulvar squamous cell carcinoma
Melanoma-associated antigen 1
Melanoma-associated antigen 4
Peripheral blood mononuclear cells
Phosphate buffered saline
Programmed cell death protein 1
Programmed cell death ligand 1
Receiver operating characteristics
Stimulator of interferon genes
Central memory T cell
Effector memory T cell
Effector memory RA+ T cells
Tumor infiltrating lymphocyte
Tumor necrosis factor alpha
Regulatory T cell
Vulvar squamous cell carcinoma
Mellman I, Hubbard-Lucey VM, Tontonoz MJ, Kalos MD, Chen DS, Allison JP, et al. De-risking immunotherapy: report of a consensus workshop of the Cancer immunotherapy consortium of the Cancer Research Institute. Cancer Immunol Res. 2016;4(4):279–88. https://doi.org/10.1158/2326-6066.Cir-16-0045.
van der Velden J, Fons G, Lawrie TA. Primary groin irradiation versus primary groin surgery for early vulvar cancer. Cochrane Database Syst Rev. 2011;(5):Cd002224. https://doi.org/10.1002/14651858.CD002224.pub2.
Gaarenstroom K, Kenter G, Trimbos J, Agous I, Amant F, Peters AA, et al. Postoperative complications after vulvectomy and inguinofemoral lymphadenectomy using separate groin incisions. Int J Gynecol Cancer. 2003;13(4):522–7.
Te Grootenhuis NC, Pouwer AW, de Bock GH, Hollema H, Bulten J, van der Zee AGJ, et al. Prognostic factors for local recurrence of squamous cell carcinoma of the vulva: A systematic review. Gynecol Oncol. 2017;148(3):622–31. https://doi.org/10.1016/j.ygyno.2017.11.006.
Raspollini MR, Asirelli G, Taddei GL. Analysis of lymphocytic infiltrate does not help in prediction of vulvar squamous cell carcinoma arising in a background of lichen sclerosus. Int J Gynaecol Obstet. 2008;100(2):190–1. https://doi.org/10.1016/j.ijgo.2007.07.019.
Sznurkowski JJ, Zawrocki A, Emerich J, Biernat W. Prognostic significance of CD4+ and CD8+ T cell infiltration within cancer cell nests in vulvar squamous cell carcinoma. Int J Gynecol Cancer. 2011;21(4):717–21. https://doi.org/10.1097/IGC.0b013e3182131f36.
van Esch EM, van Poelgeest MI, Kouwenberg S, Osse EM, Trimbos JB, Fleuren GJ, et al. Expression of coinhibitory receptors on T cells in the microenvironment of usual vulvar intraepithelial neoplasia is related to proinflammatory effector T cells and an increased recurrence-free survival. Int J Cancer. 2015;136(4):E95–106. https://doi.org/10.1002/ijc.29174.
van Esch EM, van Poelgeest MI, Trimbos JB, Fleuren GJ, Jordanova ES, van der Burg SH. Intraepithelial macrophage infiltration is related to a high number of regulatory T cells and promotes a progressive course of HPV-induced vulvar neoplasia. Int J Cancer. 2015;136(4):E85–94. https://doi.org/10.1002/ijc.29173.
de Jong RA, Toppen NL, Ten Hoor KA, Boezen HM, Kema IP, Hollema H, et al. Status of cellular immunity lacks prognostic significance in vulvar squamous carcinoma. Gynecol Oncol. 2012;125(1):186–93. https://doi.org/10.1016/j.ygyno.2011.12.416.
Sznurkowski JJ, Zawrocki A, Emerich J, Sznurkowska K, Biernat W. Expression of indoleamine 2,3-dioxygenase predicts shorter survival in patients with vulvar squamous cell carcinoma (vSCC) not influencing on the recruitment of FOXP3-expressing regulatory T cells in cancer nests. Gynecol Oncol. 2011;122(2):307–12. https://doi.org/10.1016/j.ygyno.2011.04.050.
Hecking T, Thiesler T, Schiller C, Lunkenheimer JM, Ayub TH, Rohr A, et al. Tumoral PD-L1 expression defines a subgroup of poor-prognosis vulvar carcinomas with non-viral etiology. Oncotarget. 2017;8(54):92890–903. https://doi.org/10.18632/oncotarget.21641.
Taube JM, Anders RA, Young GD, Xu H, Sharma R, McMiller TL, et al. Colocalization of inflammatory response with B7-h1 expression in human melanocytic lesions supports an adaptive resistance mechanism of immune escape. Sci Transl Med. 2012;4(127):127ra37. https://doi.org/10.1126/scitranslmed.3003689.
Sznurkowski JJ, Zawrocki A, Biernat W. Subtypes of cytotoxic lymphocytes and natural killer cells infiltrating cancer nests correlate with prognosis in patients with vulvar squamous cell carcinoma. Cancer Immunol Immunother. 2014;63(3):297–303. https://doi.org/10.1007/s00262-013-1511-x.
Nooij LS, Ter Haar NT, Ruano D, Rakislova N, van Wezel T, Smit V, et al. Genomic characterization of vulvar (pre) cancers identifies distinct molecular subtypes with prognostic significance. Clin Cancer Res. 2017;23(22):6781–9. https://doi.org/10.1158/1078-0432.ccr-17-1302.
van de Nieuwenhof HP, van Kempen LC, de Hullu JA, Bekkers RL, Bulten J, Melchers WJ, et al. The etiologic role of HPV in vulvar squamous cell carcinoma fine tuned. Cancer Epidemiol Biomark Prev. 2009;18(7):2061–7. https://doi.org/10.1158/1055-9965.epi-09-0209.
Dong F, Kojiro S, Borger DR, Growdon WB, Oliva E. Squamous cell carcinoma of the vulva: a subclassification of 97 cases by Clinicopathologic, Immunohistochemical, and molecular features (p16, p53, and EGFR). Am J Surg Pathol. 2015;39(8):1045–53. https://doi.org/10.1097/pas.0000000000000454.
Hinten F, Molijn A, Eckhardt L, Massuger L, Quint W, Bult P, et al. Vulvar cancer: two pathways with different localization and prognosis. Gynecol Oncol. 2018;149(2):310–7. https://doi.org/10.1016/j.ygyno.2018.03.003.
Spranger S, Bao R, Gajewski TF. Melanoma-intrinsic beta-catenin signalling prevents anti-tumour immunity. Nature. 2015;523(7559):231–5. https://doi.org/10.1038/nature14404.
Peng W, McKenzie JA, Hwu P. Complementing T-cell function: an inhibitory role of the complement system in T-cell-mediated antitumor immunity. Cancer Discov. 2016;6(9):953–5. https://doi.org/10.1158/2159-8290.Cd-16-0698.
de Sanjose S, Alemany L, Ordi J, Tous S, Alejo M, Bigby SM, et al. Worldwide human papillomavirus genotype attribution in over 2000 cases of intraepithelial and invasive lesions of the vulva. Eur J Cancer. 2013;49(16):3450–61. https://doi.org/10.1016/j.ejca.2013.06.033.
Piersma SJ, Welters MJ, van der Hulst JM, Kloth JN, Kwappenberg KM, Trimbos BJ, et al. Human papilloma virus specific T cells infiltrating cervical cancer and draining lymph nodes show remarkably frequent use of HLA-DQ and -DP as a restriction element. Int J Cancer. 2008;122(3):486–94. https://doi.org/10.1002/ijc.23162.
de Vos van Steenwijk PJ, Heusinkveld M, Ramwadhdoebe TH, Lowik MJ, van der Hulst JM, Goedemans R, et al. An unexpectedly large polyclonal repertoire of HPV-specific T cells is poised for action in patients with cervical cancer. Cancer Res. 2010;70(7):2707–17. https://doi.org/10.1158/0008-5472.can-09-4299.
Ijsselsteijn ME, Brouwer TP, Abdulrahman Z, Reidy E, Ramalheiro A, Heeren AM, et al. Cancer immunophenotyping by 7 colour multispectral imaging without tyramide signal amplification. J Pathol Clin Res. 2018;5(1):3–11. https://doi.org/10.1002/cjp2.113.
Welters MJP, Ma W, Santegoets S, Goedemans R, Ehsan I, Jordanova ES, et al. Intratumoral HPV16-specific T cells constitute a type I-oriented tumor microenvironment to improve survival in HPV16-driven oropharyngeal Cancer. Clin Cancer Res. 2018;24(3):634–47. https://doi.org/10.1158/1078-0432.ccr-17-2140.
Santegoets SJ, van Ham VJ, Ehsan I, Charoentong P, Duurland CL, van Unen V, et al. The Anatomical Location Shapes the Immune Infiltrate in Tumors of Same Etiology and Affects Survival. Clin Cancer Res. 2018;25(1):240–52. https://doi.org/10.1158/1078-0432.Ccr-18-1749.
Santegoets SJ, Dijkgraaf EM, Battaglia A, Beckhove P, Britten CM, Gallimore A, et al. Monitoring regulatory T cells in clinical samples: consensus on an essential marker set and gating strategy for regulatory T cell analysis by flow cytometry. Cancer Immunol Immunother. 2015;64(10):1271–86. https://doi.org/10.1007/s00262-015-1729-x.
Santegoets SJ, Duurland CL, Jordanova ES, van Ham JJ, Ehsan I, van Egmond SL, et al. Tbet-positive regulatory T cells accumulate in oropharyngeal cancers with ongoing tumor-specific type 1 T cell responses. J Immunother Cancer. 2019;7(1):14. https://doi.org/10.1186/s40425-019-0497-0.
Pezzotti N, Höllt T, Lelieveldt B, Eisemann E, Vilanova A. Hierarchical Stochastic Neighbor Embedding. Compt Graph Forum. 2016;35:21–30.
Welters MJP, Ma W, Santegoets SJ, Goedemans R, Ehsan I, Jordanova ES, et al. Intratumoral HPV16-specific T-cells constitute a type 1 oriented tumor microenvironment to improve survival in HPV16-driven oropharyngeal cancer. Clin Cancer Res. 2017;24(3):634–47. https://doi.org/10.1158/1078-0432.ccr-17-2140.
Rumbold AR, Tan SE, Condon JR, Taylor-Thomson D, Nickels M, Tabrizi SN, et al. Investigating a cluster of vulvar cancer in young women: a cross-sectional study of genital human papillomavirus prevalence. BMC Infect Dis. 2012;12:243. https://doi.org/10.1186/1471-2334-12-243.
Monk BJ, Burger RA, Lin F, Parham G, Vasilev SA, Wilczynski SP. Prognostic significance of human papillomavirus DNA in vulvar carcinoma. Obstet Gynecol. 1995;85(5 Pt 1):709–15.
Chen DS, Mellman I. Elements of cancer immunity and the cancer-immune set point. Nature. 2017;541(7637):321–30. https://doi.org/10.1038/nature21349.
Galon J, Bruni D. Approaches to treat immune hot, altered and cold tumours with combination immunotherapies. Nat Rev Drug Discov. 2019;18(3):197–218. https://doi.org/10.1038/s41573-018-0007-y.
Bruggner RV, Bodenmiller B, Dill DL, Tibshirani RJ, Nolan GP. Automated identification of stratifying signatures in cellular subpopulations. Proc Natl Acad Sci U S A. 2014;111(26):E2770–7. https://doi.org/10.1073/pnas.1408792111.
van Esch EM, Tummers B, Baartmans V, Osse EM, Ter Haar N, Trietsch MD, et al. Alterations in classical and nonclassical HLA expression in recurrent and progressive HPV-induced usual vulvar intraepithelial neoplasia and implications for immunotherapy. Int J Cancer. 2014;135(4):830–42. https://doi.org/10.1002/ijc.28713.
Howitt BE, Sun HH, Roemer MG, Kelley A, Chapuy B, Aviki E, et al. Genetic basis for PD-L1 expression in squamous cell carcinomas of the cervix and vulva. JAMA Oncol. 2016;2(4):518–22. https://doi.org/10.1001/jamaoncol.2015.6326.
Thangarajah F, Morgenstern B, Pahmeyer C, Schiffmann LM, Puppe J, Mallmann P, et al. Clinical impact of PD-L1 and PD-1 expression in squamous cell cancer of the vulva. J Cancer Res Clin Oncol. 2019;145(6):1651–60. https://doi.org/10.1007/s00432-019-02915-1.
Garcia-Diaz A, Shin DS, Moreno BH, Saco J, Escuin-Ordinas H, Rodriguez GA, et al. Interferon receptor signaling pathways regulating PD-L1 and PD-L2 expression. Cell Rep. 2017;19(6):1189–201. https://doi.org/10.1016/j.celrep.2017.04.031.
Shields LBE, Gordinier ME. Pembrolizumab in Recurrent Squamous Cell Carcinoma of the Vulva: Case Report and Review of the Literature. Gynecol Obstet Investig. 2018;84:1–5. https://doi.org/10.1159/000491090.
Wei SC, Levine JH, Cogdill AP, Zhao Y, Anang NAS, Andrews MC, et al. Distinct Cellular Mechanisms Underlie Anti-CTLA-4 and Anti-PD-1 Checkpoint Blockade. Cell. 2017;170(6):1120–1133.e17. https://doi.org/10.1016/j.cell.2017.07.024.
Hirschhorn-Cymerman D, Budhu S, Kitano S, Liu C, Zhao F, Zhong H, et al. Induction of tumoricidal function in CD4+ T cells is associated with concomitant memory and terminally differentiated phenotype. J Exp Med. 2012;209(11):2113–26. https://doi.org/10.1084/jem.20120532.
Massarelli E, William W, Johnson F, Kies M, Ferrarotto R, Guo M, et al. Combining immune checkpoint blockade and tumor-specific vaccine for patients with incurable human papillomavirus 16-related Cancer: a phase 2 clinical trial. JAMA Oncol. 2019;5(1):67–73. https://doi.org/10.1001/jamaoncol.2018.4051.
Sznurkowski JJ, Zawrocki A, Karczewska J, Emerich J, Biernat W. Should we consider cancer/testis antigens NY-ESO-1, MAGE-A4 and MAGE-A1 as potential targets for immunotherapy in vulvar squamous cell carcinoma? Histopathology. 2011;58(3):481–3. https://doi.org/10.1111/j.1365-2559.2011.03772.x.
Corrales L, McWhirter SM, Dubensky TW Jr, Gajewski TF. The host STING pathway at the interface of cancer and immunity. J Clin Invest. 2016;126(7):2404–11. https://doi.org/10.1172/jci86892.
Bourgeois-Daigneault MC, Roy DG, Aitken AS, El Sayes N, Martin NT, Varette O, et al. Neoadjuvant oncolytic virotherapy before surgery sensitizes triple-negative breast cancer to immune checkpoint therapy. Sci Transl Med. 2018;10(422):eaao1641. https://doi.org/10.1126/scitranslmed.aao1641.
Sato-Kaneko F, Yao S, Ahmadi A, Zhang SS, Hosoya T, Kaneda MM, et al. Combination immunotherapy with TLR agonists and checkpoint inhibitors suppresses head and neck cancer. JCI Insight. 2017;2(18):93397. https://doi.org/10.1172/jci.insight.93397.
Schon MP, Schon M. Imiquimod: mode of action. Br J Dermatol. 2007;157(Suppl 2):8–13. https://doi.org/10.1111/j.1365-2133.2007.08265.x.
van Seters M, van Beurden M, ten Kate FJ, Beckmann I, Ewing PC, Eijkemans MJ, et al. Treatment of vulvar intraepithelial neoplasia with topical imiquimod. N Engl J Med. 2008;358(14):1465–73. https://doi.org/10.1056/NEJMoa072685.
Adams S, Kozhaya L, Martiniuk F, Meng TC, Chiriboga L, Liebes L, et al. Topical TLR7 agonist imiquimod can induce immune-mediated rejection of skin metastases in patients with breast cancer. Clin Cancer Res. 2012;18(24):6748–57. https://doi.org/10.1158/1078-0432.ccr-12-1149.
We gratefully thank all the patients and healthy individuals who participated in this study. Furthermore Sandra van den Broek-Veldstra and Margriet Löwik for including the patients in the Circle study.
This study was financially supported by grants from the Dutch Cancer Society 2016–10168 to MIEvP, TB, and SHvdB.
Ethics approval and consent to participate
The study was conducted in accordance with the Declaration of Helsinki and approved by the local medical ethical committee of the Leiden University Medical Center (P08.197 and B16.024) and in agreement with the Dutch law. The materials were used according to the Dutch Federation of Medical Research Association guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Antibody panels. (DOCX 19 kb)
Patient characteristics of FFPE cohort. (DOCX 23 kb)
HPVposVSCChave better overall survival compared to HPVnegVSCC. (PDF 260 kb)
Tissue segmentation and image analysis by VECTRA and total T cell infiltrate in VSCC subtypes and healthy controls. (PDF 419 kb)
Statistical differences in T-cell infiltration between VSCC subtypes and healthy controls. (DOCX 16 kb)
Pearson correlations between the numbers of intraepithelial and stromal T cells. (DOCX 14 kb)
Differences in survival for the four immune categories of VSCC. (PDF 300 kb)
Median and optimized cut-off point by ROC curve analysis per phenotype and outcome. (DOCX 17 kb)
Clinical impact of several subsets of intraepithelial T cells in the total group of VSCC and in HPVnegVSCCpatients only. (PDF 399 kb)
Uni- and multivariate analysis for recurrence-free period. (DOCX 18 kb)
Clustering analysis using CITRUS revealed 10 distinctive populations of CD4+ and CD8+ T cells. (PDF 643 kb)
Gating strategy for CD8+, non-Treg CD4+ and regulatory (Treg) T-cell populations. (PDF 263 kb)