Our discovery cohort consisted of 94 patients from the HN TCGA cohort treated with surgery and postoperative radiation for whom detailed treatment and outcomes, including patterns of failure, are known, which is not available for the full TCGA cohort. The tumor characteristics of this subset are shown in Table 1. All patients had HPV/p16-negative tumors and were treated with ~ 60 Gy following surgical resection. Most patients had tumors of the oral cavity (63.8%) or larynx/hypopharynx (31.9%), with a small number of oropharyngeal tumors (4.3%), reflecting the higher proportion of HPV-positivity of these tumors as well as the usual nonoperative management of these tumors. The majority of tumors were staged as either T3 (24.5%) or T4 (59.6%), with close to half (48.9%) found to have nodal involvement following resection.
Table 1
Tumor characteristics of the discovery cohort.
Nodal stage | n | Percent |
0 | 44 | 46.8 |
1 | 20 | 21.3 |
2x | 1 | 1.1 |
2a | 2 | 2.1 |
2b | 13 | 13.8 |
2c | 8 | 8.5 |
3 | 2 | 2.1 |
x | 4 | 4.3 |
Tumor stage | | |
1 | 1 | 1.1 |
2 | 12 | 12.8 |
3 | 23 | 24.5 |
4 | 56 | 59.6 |
x | 2 | 2.1 |
Site | | |
Oral cavity | 60 | 63.8 |
OPX | 4 | 4.3 |
Larynx/hypopharynx | 30 | 31.9 |
Th2 cell type is significantly associated with clinical outcome in two similarly treated HNSCC cohorts.
To identify an association between outcome and pretreatment TILs and other cell types, we initially examined the discovery cohort and investigated the association of clinical variables (tumor stage, nodal stage and site) and cell types (as a continuous variable) with overall survival (OS). We found that nodal stage was significantly associated with OS (p = 0.03). The presence of monocytes (HR 2.1, 95% CI 1.008–4.34, p = 0.05), Th2 cells (HR 6.1, 95% CI 1.7–21.7, p = 0.005) and CD4 + memory T cells (HR 5.9, 95% CI 1.7–20.6, p = 0.006) was associated with worse OS in this cohort (Fig. 1a). On multivariable Cox regression analysis including both the immune cell types as well as nodal stage, nodal stage (p = 0.008), monocytes (p = 0.02), and Th2 cells (p = 0.008) remained significant.
We next evaluated the association between intratumoral cell populations and distant metastasis (DM) or locoregional recurrence (LRR). Time to DM was significantly associated with higher nodal stage (p = 6.4e-4). Additionally, we found that CD4 + memory T cells (HR 11.63, 95% CI 1.8–75.9, p = 0.01), sebocytes (HR 0.068, 95% CI 0.01–0.48, p = 0.007) and epithelial cells (HR 8.7e-3, 95% CI 3.56e-4–0.21, p = 0.01) were significantly associated with DM on univariable analysis (Fig. 1b). After adjusting for nodal status in a multivariable model, the presence of CD4 + memory cells remained significantly associated with DM (p = 0.012). We next examined LRR and found no significant association with any clinical variables evaluated. However, LRR was significantly associated with specific immune infiltrates, including common lymphoid progenitor cells (CLPs) (HR 9.5 95% CI 1.21–74.9, p = 0.03), Th2 cells (HR 9.0 1.4–57.0, p = 0.02), CD4 + memory T cells (HR 7.351843, 95% CI 1.3–41.3, p = 0.024) and immature dendritic cells (IDCs) (HR 0.34, 95% CI 0.17–0.67, p = 0.002) (Fig. 1c).
To validate the associations seen between intratumoral cell types and clinical outcome, we used gene expression data from an independent patient population treated in a similar fashion (e.g., surgery and postoperative radiation)11. We repeated the same analysis using xCell to predict cell type abundance in this validation cohort and focused only on the cell types identified to be significant in the discovery cohort. Regarding OS and DM, none of the cell types identified in the discovery cohort were found to be significant in the validation cohort. However, we found that the presence of Th2 cells was associated with significantly worse LRR both in the discovery and validation cohorts when the same cutoff value was used (lower tertile vs. others) (p = 0.019) (Fig. 1d).
Somatic mutations and gene expression associated with Th2 cell infiltration
To explore the potential relationship between tumor signaling and derived Th2 infiltrate, we initially examined differences in the presence of Th2 based upon tumor mutation in the entirety of the HPV-negative Head and Neck TCGA cohort for which both mutational and gene expression data were available (n = 406). Mutations in several genes, specifically CASP8, HRAS, and HLA-A, were associated with an increased presence of derived Th2 infiltrate (Fig. 2a). Although neither EP300 nor CREBBP mutation was individually associated with increased Th2 infiltrate, if taken together, mutations in these functionally similar histone (or lysine) aceytltransferases were also associated with increased Th2 infiltrate (p = 0.042, Fig. 2a).
We next examined the relationship between Th2 infiltrate and gene expression using RNA-seq data in the HPV-negative Head and Neck TCGA cohort, for which only gene expression data were required (n = 409). These data were analyzed using IPA, as described in the Methods, to identify pathways associated with derived Th2 infiltrate. This analysis identified several cancer- and immune-related pathways (hits related to cancer or immune response shown in Fig. 2b). For instance, the “BRCA1 in DRR” pathway positively and the “Tumor microenvironment (TME)” pathway negatively associated with Th2 infiltrate. Signaling nodes and expression patterns for both pathways are shown in Fig. 2c and Fig. 2d, respectively, with nodes significantly positively and negatively correlated with derived Th2 cells.
To confirm the relationship between these two pathways and derived Th2 infiltrate, pathway activation scores were generated as described in the Materials & Methods. As expected, a high degree of positive correlation was observed between derived Th2 infiltrate and BRCA1 in the DDR pathway (r = 0.73, p < 2.2e-16), while a significant negative correlation was observed with the TME pathway (r=-0.45, p < 2.2e-16) (Fig. 3a). We next generated pathway activation scores using gene expression data from a separate nonoverlapping HPV-negative HNSCC cohort and observed similar associations with derived Th2 infiltrate (Fig. 3b). Specifically, derived Th2 infiltrate was correlated with activation of BRCA1 in the DDR pathway (r = 0.77, p = 5.6e-9) and repression of the TME pathway (r=-0.44, p = 0.007). We also analyzed proteomic data from the separate HPV-negative Head and Neck CPTAC database and again found the same pathway associations, with derived Th2 infiltrate associated with activation of BRCA1 in DDR (r = 0.59, p < 2.2e-16) and repression of the TME pathway (r=-0.48, p = 1.6e-7) (Fig. 3c).