Next Article in Journal
Multimodal Treatment of Pleural Mesothelioma with Cytoreductive Surgery and Hyperthermic Intrathoracic Chemotherapy: Impact of Additive Chemotherapy
Previous Article in Journal
Theranostic Risk Stratification for Thyroid Cancer in the Genomic Paradigm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Th2 Cells Are Associated with Tumor Recurrence Following Radiation

1
Department of Radiation Oncology, UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA 15232, USA
2
Department of Medicine, UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA 15232, USA
3
Department of Thoracic-Head and Neck Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
4
Department of Radiation Oncology, Stanford University, Stanford, CA 94305, USA
5
Department of Otolaryngology, UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA 15232, USA
6
Department of Surgery—Otolaryngology, Yale University, New Haven, CT 06520, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
These authors also contributed equally to this work.
Cancers 2024, 16(8), 1586; https://doi.org/10.3390/cancers16081586
Submission received: 28 March 2024 / Revised: 8 April 2024 / Accepted: 16 April 2024 / Published: 20 April 2024
(This article belongs to the Section Tumor Microenvironment)

Abstract

:

Simple Summary

Locoregional recurrence following the definitive treatment of head and neck cancers often leads to patient demise. Radiation is frequently used in the definitive treatment of these cancers. Resistance to radiation leads to increased locoregional recurrence, particularly in patients with human papillomavirus-negative head and neck squamous cell carcinomas. In this study, we show that tumor Th2 immune infiltrate is associated with increased locoregional recurrence in patients with human papillomavirus-negative head and neck squamous cell carcinoma. These findings can be used to develop novel immunotherapeutic approaches to improve the response to radiation and thus decrease locoregional recurrence and improve outcomes in patients with human papillomavirus-negative head and neck squamous cell carcinomas.

Abstract

The curative treatment of multiple solid tumors, including head and neck squamous cell carcinoma (HNSCC), utilizes radiation. The outcomes for HPV/p16-negative HNSCC are significantly worse than HPV/p16-positive tumors, with increased radiation resistance leading to worse locoregional recurrence (LRR) and ultimately death. This study analyzed the relationship between immune function and outcomes following radiation in HPV/p16-negative tumors to identify mechanisms of radiation resistance and prognostic immune biomarkers. A discovery cohort of 94 patients with HNSCC treated uniformly with surgery and adjuvant radiation and a validation cohort of 97 similarly treated patients were utilized. Tumor immune infiltrates were derived from RNAseq gene expression. The immune cell types significantly associated with outcomes in the discovery cohort were examined in the independent validation cohort. A positive association between high Th2 infiltration and LRR was identified in the discovery cohort and validated in the validation cohort. Tumor mutations in CREBBP/EP300 and CASP8 were significantly associated with Th2 infiltration. A pathway analysis of genes correlated with Th2 cells revealed the potential repression of the antitumor immune response and the activation of BRCA1-associated DNA damage repair in multiple cohorts. The Th2 infiltrates were enriched in the HPV/p16-negative HNSCC tumors and associated with LRR and mutations in CASP8, CREBBP/EP300, and pathways previously shown to impact the response to radiation.

1. Introduction

The curative treatment of nearly all locally advanced solid tumors involves radiation; however, there are few biologically derived drivers of response to radiation. Moreover, although the immune response impacts tumor radioresistance in preclinical models [1], the clinical impact of specific tumor-infiltrating lymphocytes (TILs) on the radiation response is unclear.
Head and neck squamous cell carcinoma (HNSCC) leads to the death of over 360,000 patients annually worldwide [2]. Outcomes are driven by tumor characteristics and human papillomavirus (HPV) status, with HPV-positive tumors exhibiting better outcomes [3]. In HNSCC, locoregional recurrence (LRR) following radiation, and not distant metastasis (DM), is the most common cause of cancer death [4]. Thus, radioresistance is a direct cause of death in HNSCC, particularly in HPV-negative disease. While there are no clinically utilized biomarkers of therapeutic resistance for HPV-negative tumors, genes such as TP53, CREBBP, EP300, and CASP8 may predict radioresistance [5,6,7,8].
Radiation relies on immune activation to exert maximal antitumor effects [9] via multiple pathways, including the release of tumor antigens and proinflammatory cytokines/chemokines and the activation of TILs and dendritic cells [10,11]. Furthermore, PD-L1 overexpression has been associated with worse LRR following radiation [12]. These observations, coupled with anti-PD1 activity in HNSCC [13,14], has led to clinical trials of concurrent immunotherapy and radiation [15,16,17]. Unfortunately, these studies did not meet their primary endpoints, although post hoc analyses have indicated some benefit in PD-L1-high patients [16,17,18,19]. Regardless, these results highlight the need for additional study of the intersection between curative radiation and the immune response.
In this study, we utilized two HPV-negative HNSCC populations treated with surgical resection and adjuvant radiation to examine TILs associated with LRR. Pretreatment tumors from both cohorts were digitally dissected for the abundance of TILs using an integration of gene set enrichment analysis and deconvolution approaches (xCell) [20]. Univariable and multivariable analyses for the association between immune cell type and clinical outcome were performed in the discovery cohort. Significant cell types from this analysis were then examined in the independent validation cohort. High-T helper 2 (Th2) infiltrate was identified as being significantly associated with LRR in both cohorts. Additional analysis was performed to identify mutations and pathways reliably associated with Th2 cell infiltration.

2. Materials and Methods

2.1. Discovery Cohort and Validation Cohort Clinical Outcome Datasets

The discovery cohort reflects 94 patients from the TCGA cohort with HPV-negative HNSCC treated with surgery and postoperative radiation, with detailed treatment and clinical outcomes [7]. The validation cohort consisted of 97 HPV-negative HNSCC patients treated with surgery and postoperative radiation for which expression data necessary for TIL prediction were available [12]. Tumors were staged using AJCC 7th edition.

2.2. Genomic and Transcriptomic Datasets

Tumor mutation and RNA-seq data available for the discovery cohort and the TCGA cohort have been previously described and downloaded from cBioPortal [21], while Illumina mRNA array data from the validation cohort have been published previously [22]. Unique (e.g., nonoverlapping with TCGA) RNA sequencing data from the International Cancer Genome Consortium (ICGC) data portal [23,24] (accessed on 1 October 2021) were used to validate the pathway analysis. A separate pathway validation was performed using additional RNA Seq data from the NCI’s Clinical Proteomic Tumor Analysis Consortium (CPTAC) HNSCC cohort (accessed on 1 November 2021) [25].

2.3. Digital Dissection of the Tumor Microenvironment for TIL Abundance

RNAseq data from both the discovery and validation cohorts were used to predict TIL abundance in tumors using xCell [20] (https://xcell.ucsf.edu/) (analyzed on 10 September 2021). xCell quantifies 64 immune, stromal, and tumor cell abundances in a heterogeneous tumor microenvironment from tissue RNA-seq data [20]. Specifically, -xCell was used to digitally deconvolute cell populations from RNAseq data using R library xcell:: xCellAnalysis with default settings. This method employs both ssGESA and deconvolution algorithms to compute the scores for each cell population in a sample, including Th2 cells, among others. All samples were organized based on their Th2 cell scores, which were generated using xCell. Subsequently, they were divided into three groups—T1, T2, and T3—with each group representing one-third of the total number of immune cell types present in at least 15 tumors and significantly associated with clinical outcome in the discovery cohort were selected for further validation using the validation cohort.

2.4. Identification of Somatic Mutations Associated with Th2 Cell Infiltrates

Somatic mutations deemed significant in the TCGA cohort according to MutSig [26] and identified as cancer genes in OncoKB [27] were determined. Average derived Th2 infiltrate was compared between mutant versus wild type for each mutation using one-way ANOVA and post hoc comparison (SPSS v27). Unadjusted p-values are shown for survival analysis. For pathway analysis, the multiple comparisons effect was adjusted using the Benjamini–Hochberg (BH) method [28].

2.5. Identification of Genes Positively or Negatively Correlated with Th2 Cell Infiltrates by Expression

Ingenuity pathway analysis (IPA) (Qiagen) was performed to determine enriched pathways and activated z-scores for genes with ≥1.5-fold difference between the upper and lower tertiles of derived Th2 infiltrate at an FDR of 0.05 in the HPV-negative TCGA cohort (n = 409). Significance of pathway enrichment was determined using the entire human genome as the reference set. The analysis of the IPA was performed with species selected as “Human, Mouse, Rat” and other parameters as default. The reference set in the IPA was set to “Ingenuity Knowledge Base (Genes only)” with analysis performed on 9 November 2021.

2.6. Survival Analysis

Overall survival (OS) was defined as the time from diagnosis until death or the last follow-up. Time to LRR or DM was defined as the time from diagnosis until either an event or the last follow-up. Clinical variables included site, nodal stage (0–1 versus ≥ 2a), and tumor stage (1–2 versus 3–4). xCell z-scores were transformed by log10 (score × 100 + 1) for the survival analysis. Univariable Cox proportional hazards (PH) model was used to test the associations between the abundance of each immune infiltrate and outcomes. Variables significant on univariable analysis were included in multivariable analysis where indicated. Kaplan–Meier analysis (SPSS (v27)) was performed using a log-rank test comparing the lower versus upper two tertiles of derived Th2 cells in both cohorts.

3. Results

The tumor characteristics of the discovery cohort are shown in Table 1. All the patients had HPV/p16-negative tumors and were treated with ~60 Gy following resection. Most of the patients had tumors of the oral cavity (63.8%) or larynx/hypopharynx (31.9%), with a small number of oropharyngeal tumors (4.3%). Most of the tumors were staged as either T3 (24.5%) or T4 (59.6%), with 48.9% of the patients found to exhibit nodal involvement following resection.

3.1. Th2 Cell Type Is Significantly Associated with Clinical Outcome in Two Similarly Treated HNSCC Cohorts

We examined the discovery cohort for an association between the clinical variables and cell types with OS. The 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) was not significantly associated with OS; however, both 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) were significantly associated with a worse OS (Figure 1A). On the multivariable Cox regression analysis, the nodal stage (p = 0.008), monocytes (p = 0.02), and Th2 cells (p = 0.008) remained significant.
We next evaluated the association between the intratumoral cell populations and DM or LRR. The time to DM was significantly associated with a higher nodal stage (p = 6.4 × 10−4). In addition, the 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.7 × 10−3, 95% CI 3.56 × 10−4–0.21, p = 0.01) were significantly associated with DM on the univariable analysis (Figure 1B). After adjusting for the 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 TILs, 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) (Figure 1C).
To validate the associations seen between the intratumoral cell types and clinical outcome, we used gene expression data from an independent patient population treated in a similar fashion [12]. We repeated the same analysis using xCell to predict the cell type abundance in this validation cohort and focused only on the cell types identified as being significant in the discovery cohort. None of the cell types identified in the discovery cohort were found to be significant for OS or DM 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) (Figure 1D).

3.2. Somatic Mutations and Gene Expression Associated with Th2 Cell Infiltration

To explore the relationship between tumor signaling and Th2 infiltrate, we examined differences in the presence of Th2 based upon tumor mutation in the entirety of the HPV-negative 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 (Figure 2A). Although neither EP300 nor CREBBP mutations were individually associated with increased Th2 infiltrate, if taken together, mutations in these functionally similar histone acetyltransferases were associated with increased Th2 infiltrate (p = 0.042, Figure 2A). When correcting for multiple hypothesis testing, only CASP8 remained significant, with an FDR of 5%.
We next examined the relationship between the Th2 infiltrates and gene expression using the RNA-seq data of the HPV-negative TCGA cohort, for which only gene expression data were required (n = 409). These data were analyzed using IPA to identify pathways associated with the derived Th2 infiltrate. Several cancer- and immune-related pathways were identified (Figure 2B). For instance, the “BRCA1 in DRR” pathway and the “Tumor microenvironment (TME)” pathway were positively and negatively associated with the Th2 infiltrate, respectively. The signaling nodes and expression patterns for both pathways are shown in Figure A1 and Figure A2, respectively, with the nodes significantly positively or negatively correlated with the Th2 infiltrate.
To confirm the relationship between these two pathways and the derived Th2 infiltrate, pathway activation scores were generated as described in the Methods section. As expected, a high degree of positive correlation was observed between the derived Th2 infiltrate and BRCA1 in the DDR pathway (r = 0.73, p < 2.2 × 10−16), while a significant negative correlation was observed with the TME pathway (r = −0.45, p < 2.2 × 10−16) (Figure 3A). Pathway activation scores were generated using gene expression data from a separate nonoverlapping HPV-negative HNSCC cohort [22,23] and similar associations with the derived Th2 infiltrate were observed (Figure 3B). Specifically, the Th2 infiltrate was correlated with the activation of BRCA1 in the DDR pathway (r = 0.77, p = 5.6 × 10−9) and the repression of the TME pathway (r = −0.44, p = 0.007). We also analyzed RNA-Seq data from the separate HPV-negative CPTAC database and again found the same pathway associations, with the derived Th2 infiltrate associated with the activation of BRCA1 in DDR (r = 0.59, p < 2.2 × 10−16) and the repression of the TME pathway (r = −0.48, p = 1.6 × 10−7) (Figure 3C).
The correlation of the Th2 infiltrate with other pathways, specifically the Cell Cycle Control of Chromosomal Replication Pathway as well as the Kinetochore Metaphase Signaling Pathway was performed across the TCGA, ICGC, and CPTAC cohorts (Figure A3). In all the cohorts, both pathways are significantly positively correlated with the Th2 infiltrate.

4. Discussion

Patients with HNSCC are currently stratified by HPV/p16 status and clinical variables; however, there is a large variability in the response to therapy within these groups. Attempts to characterize HNSCC based on immune infiltrate and link these characterizations to outcomes have been hampered by small numbers, the heterogeneity of HNSCC itself, and a lack of uniformly treated patients [29,30,31].
Previously, our group demonstrated that high PD-L1 expression was associated with increased LRR after treatment with surgery and postoperative radiation in patients with HPV-negative HNSCC [12]. In contrast, patients with high CD8 TILs and low PD-L1 expression exhibited no LRR or death due to disease [12], suggesting that pretreatment TILs in the TME have a strong influence on outcomes. In this study, we build upon our previous work and investigate the TME in uniformly treated patients to elucidate additional biomarkers.
While the role of TILs in solid tumors is well known, the focus has primarily been on CD8+ or Th1 TILs as effector immune cells or regulatory T cells as suppressive immune cells. Conversely, Th2 cells and their signature cytokines (IL-4, IL-5, and IL-13) are largely characterized by their role in helminth infection or allergic responses and the promotion of IgE-mediated eosinophilic responses [32]. In the context of cancer, Th2 cells are considered to be protumorigenic [33], with potential protumorigenic actions of Th2 cells being reported, including the recruitment and repolarization of macrophages to suppressive (M2-like) lineages [34], the stimulation of the production of vascular endothelial growth factor (VEGF) [35], and the production of immunosuppressive cytokines [33]. The variability of the Th2 response in different preclinical tumor models may be in part due to plasticity between Th1 and Th2 polarization, with factors such as VEGF shown to repolarize T cells from a Th1 to a Th2 phenotype [36], while PD-1 blockade has been shown to revert the Th2 to Th1 phenotype [37]. The role of Th2 cells in the tumor microenvironment of clinical tumor samples is also unclear, with Th2 cells being associated with aggressive tumor histologies, including HNSCC [38], as well as histologies with better prognosis, such as Hodgkin’s lymphoma [39]. Analogous to our findings in this report, a recent study by Rui et al. demonstrated that high-Th2 infiltrate in prostate cancer was associated with increased recurrence after prostatectomy [40].
In this study, we also examined tumor genetic mutations associated with Th2 infiltrate. We found that mutations in HLA-A, CASP8, and HRAS correlated with the Th2 gene signature. Additionally, CREBBP/EP300 mutations were also associated with increased Th2 infiltrate. The CASP8 and CREBBP/EP300 mutations were of particular interest since we have shown that both mutations are associated with radioresistance in HNSCC [7,8]. However, caspase 8 has been shown to be critical in nuclear factor-κβ (NF-κβ) pathway activation, causing the IL-1 secretion required for the Th2 response [41]. This discordance between the dependence of Th2 function on caspase 8 and an association between CASP8 mutation and Th2 infiltrate may be at least partially explained by cancer-associated missense mutations in CASP8 associated with epithelial cancers that inhibit its apoptotic function but chronically activate NF-κβ signaling [42]. In patients whose tumors exhibit a Th2 gene signature, we found increased expression of genes related to the cell cycle and DNA repair, which can confer radioresistance. These effects may in part be generated through the Th2 cytokine IL-4, which has been shown to increase the expression of nonhomologous end joining DNA repair proteins and induce the proliferation of keratinocytes through c-myc-driven G0/G1-to-S phase cell cycle progression [43,44].
While our findings have revealed a novel association between outcomes and Th2 infiltrates in HPV-negative HNSCC, there are several limitations to our approach. Foremost, both our discovery and validation cohort analyses were performed in a retrospective fashion using prediction models for infiltrates based on mRNA expression data. Due to this limitation, further analyses of the TME in these samples, including direct TIL visualization, functional analyses of TIL infiltrates, or measurement of cytokine levels/activity, were not possible and will be performed in the future. Despite these limitations, our study did reveal a potential new prognostic biomarker for HPV-negative HNSCC that could be utilized to inform future treatment selection and clinical trial design.
The limitations of this study include the retrospective nature of the analyses, the limited sample size, and the inherent limitations of the transcriptomic, genomic, and proteomic analyses.

5. Conclusions

An increased Th2 immune infiltrate was associated with locoregional recurrence in a discovery cohort of HPV-negative head and neck squamous cell carcinoma patients treated with surgery and postoperative radiation. This was confirmed in a separate validation cohort of similarly treated HPV-negative head and neck squamous cell carcinoma patients. These findings point to potential immunotherapeutic approaches to augment the response to radiation and decrease LRR in HPV-negative HNSCC patients.

Author Contributions

Conceptualization, H.D.S., M.K.A., R.B.P., C.R.P. and R.B.; methodology, H.D.S. and M.K.A.; validation, H.D.S. and R.B.; formal analysis, H.D.S. and R.B.; investigation, H.D.S. and M.K.A.; data curation, H.D.S., R.B., M.K.A., D.M. and J.M.; writing—original draft preparation, H.D.S. and M.K.A.; writing—review and editing, H.D.S., M.K.A., R.B., P.M.P., D.M., J.M., A.H., B.B., J.V.H., J.J.L., R.L.F., C.R.P., J.H.W. and R.B.P.; visualization, H.D.S. and M.K.A.; supervision, H.D.S. and C.R.P.; project administration, H.D.S.; funding acquisition, H.D.S., C.R.P. and J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by (i) the National Cancer Institute R01CA168485-08 (HS) and P50CA097190-15 (RF & HS), (ii) the National Institute for Dental and Craniofacial Research R01 DE028105 (HS) and R01 DE028061 (HS and CP), and (iii) T32 training grant supporting Pifer (5T32CA060397-27).

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki, and the data come from public data.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data used in this publication were generated by The Cancer Genome Atlas Program (TCGA) (https://portal.gdc.cancer.gov/) (accessed on 1 October 2021), International Cancer Genome Consortium (ICGC) (https://dcc.icgc.org/repositories) (accessed on 1 October 2021), and Clinical Proteomic Tumor Analysis Consortium (CPTAC) (https://proteomic.datacommons.cancer.gov/pdc/) (accessed on 1 November 2021). Validation cohort data have been previously published [7].

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Appendix A

Figure A1. “Role of BRCA1 in DDR” and “Tumor Microenvironment Pathway” from IPA are shown. Targets in red and green are significantly (FDR < 0.05, fold change > 1.5) positively and negatively correlated with derived Th2 infiltrate.
Figure A1. “Role of BRCA1 in DDR” and “Tumor Microenvironment Pathway” from IPA are shown. Targets in red and green are significantly (FDR < 0.05, fold change > 1.5) positively and negatively correlated with derived Th2 infiltrate.
Cancers 16 01586 g0a1
Figure A2. Signaling nodes and expression patterns for “BRCA1 in DDR” and “Tumor microenvironment pathway”.
Figure A2. Signaling nodes and expression patterns for “BRCA1 in DDR” and “Tumor microenvironment pathway”.
Cancers 16 01586 g0a2
Figure A3. Th2 infiltrate is significantly correlated with the Cell Cycle Control of Chromosomal Replication Pathway as well as the Kinetochore Metaphase Signaling Pathway across the TCGA, ICGC, and CPTAC cohorts.
Figure A3. Th2 infiltrate is significantly correlated with the Cell Cycle Control of Chromosomal Replication Pathway as well as the Kinetochore Metaphase Signaling Pathway across the TCGA, ICGC, and CPTAC cohorts.
Cancers 16 01586 g0a3aCancers 16 01586 g0a3b

References

  1. Sharabi, A.B.; Lim, M.; DeWeese, T.L.; Drake, C.G. Radiation and checkpoint blockade immunotherapy: Radiosensitisation and potential mechanisms of synergy. Lancet Oncol. 2015, 16, e498–e509. [Google Scholar] [CrossRef] [PubMed]
  2. Sung, H.; Ferlay, J.; Siegel, R.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA. Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed]
  3. Maxwell, J.H.; Grandis, J.R.; Ferris, R.L. HPV-Associated Head and Neck Cancer: Unique Features of Epidemiology and Clinical Management. Annu. Rev. Med. 2016, 67, 91–101. [Google Scholar] [CrossRef]
  4. Brockstein, B.; Haraf, D.J.; Rademaker, A.W.; Kies, M.S.; Stenson, K.M.; Rosen, F.; Mittal, B.B.; Pelzer, H.; Fung, B.B.; Witt, M.E.; et al. Patterns of failure, prognostic factors and survival in locoregionally advanced head and neck cancer treated with concomitant chemoradiotherapy: A 9-year, 337-patient, multi-institutional experience. Ann. Oncol. Off. J. Eur. Soc. Med. Oncol. 2004, 15, 1179–1186. [Google Scholar] [CrossRef]
  5. Poeta, M.L.; Manola, J.; Goldwasser, M.A.; Forastiere, A.; Benoit, N.; Califano, J.A.; Ridge, J.A.; Goddwin, J.; Kenady, D.; Saunders, J.; et al. TP53 Mutations and Survival in Squamous-Cell Carcinoma of the Head and Neck. N. Engl. J. Med. 2007, 357, 2552–2561. [Google Scholar] [CrossRef] [PubMed]
  6. Skinner, H.D.; Sandulache, V.C.; Ow, T.J.; Meyn, R.E.; Yordy, J.S.; Beadle, B.M.; Fitzgerald, A.L.; Giri, U.; Ang, K.K.; Meyers, J.N. TP53 disruptive mutations lead to head and neck cancer treatment failure through inhibition of radiation-induced senescence. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2012, 18, 290–300. [Google Scholar] [CrossRef] [PubMed]
  7. Kumar, M.; Molkentine, D.; Molkentine, J.; Bridges, K.; Xie, T.; Yang, L.; Hefner, A.; Gao, M.; Bahri, R.; Dhawan, A.; et al. Inhibition of histone acetyltransferase function radiosensitizes CREBBP/EP300 mutants via repression of homologous recombination, potentially targeting a gain of function. Nat. Commun. 2021, 12, 6340. [Google Scholar] [CrossRef] [PubMed]
  8. Uzunparmak, B.; Gao, M.; Lindemann, A.; Erikson, K.; Wang, L.; Eric, L.; Frank, S.J.; Gleber-Netto, F.O.; Zhao, M.; Skinner, H.D.; et al. Caspase 8 loss radiosensitizes head and neck squamous cell carcinoma to SMAC mimetic induced necroptosis. JCI Insight 2020, 5, e139837. [Google Scholar] [CrossRef]
  9. Lee, Y.; Auh, S.L.; Wang, Y.; Burnette, B.; Wang, Y.; Meng, Y.; Beckett, M.; Sharma, R.; Chin, R.; Tu, T.; et al. Therapeutic effects of ablative radiation on local tumor require CD8+ T cells: Changing strategies for cancer treatment. Blood 2009, 114, 589–595. [Google Scholar] [CrossRef]
  10. Golden, E.B.; Apetoh, L. Radiotherapy and immunogenic cell death. Semin. Radiat. Oncol. 2015, 25, 11–17. [Google Scholar] [CrossRef]
  11. Miyauchi, S.; Kim, S.S.; Pang, J.; Gold, K.A.; Gutkin, J.S.; Califano, J.A.; Mell, L.K.; Cohen, E.E.W.; Sharabi, A.B. Immune Modulation of Head and Neck Squamous Cell Carcinoma and the Tumor Microenvironment by Conventional Therapeutics. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2019, 25, 4211–4223. [Google Scholar] [CrossRef] [PubMed]
  12. Skinner, H.D.; Giri, U.; Yang, L.P.; Kumar, M.; Liu, Y.; Story, M.D.; Pickering, C.R.; Byers, L.A.; Williams, M.D.; Wang, J.; et al. Integrative Analysis Identifies a Novel AXL-PI3 Kinase-PD-L1 Signaling Axis Associated with Radiation Resistance in Head and Neck Cancer. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2017, 23, 2713–2722. [Google Scholar] [CrossRef] [PubMed]
  13. Ferris, R.L.; Blumenschein, G., Jr.; Fayette, J.; Guigay, J.; Colevas, A.D.; Licitra, L.; Harrington, K.; Kasper, S.; Vokes, E.E.; Even, C.; et al. Nivolumab for Recurrent Squamous-Cell Carcinoma of the Head and Neck. N. Engl. J. Med. 2016, 375, 1856–1867. [Google Scholar] [CrossRef] [PubMed]
  14. Burtness, B.; Harrington, K.J.; Greil, R.; Soulieres, D.; Tahara, M.; Castro, G., Jr.; Psyrri, A.; Baste, N.; Neupane, P.; Bratland, A.; et al. Pembrolizumab alone or with chemotherapy versus cetuximab with chemotherapy for recurrent or metastatic squamous cell carcinoma of the head and neck (KEYNOTE-048): A randomised, open-label, phase 3 study. Lancet 2019, 394, 1915–1928. [Google Scholar] [CrossRef] [PubMed]
  15. Cohen, E.E.W.; Bell, R.B.; Bifulco, C.B.; Burtness, B.; Gillison, M.L.; Harrington, K.J.; Le, Q.; Lee, N.Y.; Leidner, R.; Lewis, R.L.; et al. The Society for Immunotherapy of Cancer consensus statement on immunotherapy for the treatment of squamous cell carcinoma of the head and neck (HNSCC). J. Immunother. Cancer 2019, 7, 184. [Google Scholar] [CrossRef] [PubMed]
  16. Lee, N.Y.; Ferris, R.L.; Psyrri, A.; Haddad, R.I.; Tahara, M.; Bourhis, J.; Harrington, K.; Chang, P.M.; Lin, J.; Abdul Razaq, M.; et al. Avelumab plus standard-of-care chemoradiotherapy versus chemoradiotherapy alone in patients with locally advanced squamous cell carcinoma of the head and neck: A randomised, double-blind, placebo-controlled, multicentre, phase 3 trial. Lancet Oncol. 2021, 22, 450–462. [Google Scholar] [CrossRef] [PubMed]
  17. Tao, Y.; Biau, J.; Sun, X.S.; Sire, C.; Martin, L.; Alfonsi, M.; Prevost, J.B.; Modesto, A.; Lafond, C.; Tourani, J.M.; et al. Pembrolizumab versus cetuximab concurrent with radiotherapy in patients with locally advanced squamous cell carcinoma of head and neck unfit for cisplatin (GORTEC 2015-01 PembroRad): A multicenter, randomized, phase II trial. Ann. Oncol. Off. J. Eur. Soc. Med. Oncol. 2023, 34, 101–110. [Google Scholar] [CrossRef] [PubMed]
  18. Machiels, J.; Tao, Y.; Burtness, B.; Tahara, M.; Rischin, D.; Alves, G.V.; Lima, I.P.F.; Hughes, B.G.M.; Pointreau, Y.; Aksoy, S.; et al. Primary results of the phase 3 KEYNOTE-412 study: Pembrolizumab plus chemoradiation therapy (CRT) vs placebo plus CRT for locally advanced head and neck squamous cell carcinoma. Ann. Oncol. 2022, 33, S808–S869. [Google Scholar] [CrossRef]
  19. Mell, L.K.; Torres-Saavedra, P.A.; Wong, S.J.; Chang, S.S.; Kish, J.A.; Minn, A.; Jordan, R.C.; Liu, T.; Truong, M.T.; Winquist, E.W.; et al. Radiotherapy with Durvalumab vs. Cetuximab in Patients with Locoregionally Advanced Head and Neck Cancer and a Contraindication to Cisplatin: Phase II Results of NRGHN004. Int. J. Radiat. Oncol. Biol. Phys. 2022, 114, 1058. [Google Scholar] [CrossRef]
  20. Aran, D.; Hu, Z.; Butte, A.J. xCell: Digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017, 18, 220. [Google Scholar] [CrossRef]
  21. The Cancer Genome Atlas Network. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature 2015, 517, 576–582. [Google Scholar] [CrossRef] [PubMed]
  22. Skinner, H.D.; Giri, U.; Yang, L.; Woo, S.H.; Story, M.D.; Pickering, C.R.; Byers, L.A.; Williams, M.D.; El-Naggar, A.; Wang, J.; et al. Proteomic Profiling Identifies PTK2/FAK as a Driver of Radioresistance in HPV-negative Head and Neck Cancer. Clin. Cancer Res. 2016, 22, 4643–4650. [Google Scholar] [CrossRef] [PubMed]
  23. Zhang, J.; Bajari, R.; Andric, D.; Gerthoffert, F.; Lepsa, A.; Nahal-Bose, H.; Stein, L.D.; Ferretti, V. The International Cancer Genome Consortium Data Portal. Nat. Biotechnol. 2019, 37, 367–369. [Google Scholar] [CrossRef]
  24. Maitra, A.; Biswas, N.K.; Kumar, S.; Das, S.; Majumder, P.P.; Bagchi, I.; Bairagya, B.B.; Basu, A.; Das, D.; Dhar, R.; et al. Mutational landscape of gingivo-buccal oral squamous cell carcinoma reveals new recurrently-mutated genes and molecular subgroups. Nat. Commun. 2013, 4, 2873. [Google Scholar]
  25. Huang, C.; Chen, L.; Savage, S.R.; Eguez, R.V.; Dou, Y.; Li, Y.; Leprevost, F.V.; Jaehnig, E.J.; Lei, J.T.; Wen, B.; et al. Proteogenomic insights into the biology and treatment of HPV-negative head and neck squamous cell carcinoma. Cancer Cell. 2021, 39, 361–379.e16. [Google Scholar] [CrossRef] [PubMed]
  26. Lawrence, M.S.; Stojanov, P.; Polak, P.; Kryukov, G.V.; Cibulskis, K.; Sivachenko, A.; Carter, S.L.; Stewart, C.; Mermel, C.H.; Roberts, S.A.; et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature 2013, 499, 214–218. [Google Scholar] [CrossRef]
  27. Chakravarty, D.; Gao, J.; Phillips, S.M.; Kundra, R.; Zhang, H.; Wang, J.; Rudolph, J.E.; Yaeger, R.; Soumerai, T.; Nissan, M.H.; et al. OncoKB: A Precision Oncology Knowledge Base. JCO Precis. Oncol. 2017, 1, 1–16. [Google Scholar] [CrossRef] [PubMed]
  28. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  29. Wondergem, N.E.; Nauta, I.H.; Muijlwijk, T.; Leemans, C.R.; Van de Ven, R. The Immune Microenvironment in Head and Neck Squamous Cell Carcinoma: On Subsets and Subsites. Curr. Oncol. Rep. 2020, 22, 81. [Google Scholar] [CrossRef]
  30. Spector, M.E.; Bellile, E.; Amlani, L.; Zarins, K.; Smith, J.; Brenner, J.C.; Rozek, L.; Nguyen, A.; Thomas, D.; McHugh, J.B.; et al. Prognostic Value of Tumor-Infiltrating Lymphocytes in Head and Neck Squamous Cell Carcinoma. JAMA Otolaryngol. Neck Surg. 2019, 145, 1012–1019. [Google Scholar] [CrossRef]
  31. Mandal, R.; Senbabaoglu, Y.; Desrichard, A.; Havel, J.J.; Dalin, M.G.; Riaz, N.; Lee, K.; Ganly, I.; Hakimi, A.A.; Chan, T.A.; et al. The head and neck cancer immune landscape and its immunotherapeutic implications. JCI Insight 2016, 1, e89829. [Google Scholar] [CrossRef] [PubMed]
  32. Walker, J.A.; McKenzie, A.N.J. TH2 cell development and function. Nat. Rev. Immunol. 2018, 18, 121–133. [Google Scholar] [CrossRef] [PubMed]
  33. Hong, S.; Qian, J.; Yang, J.; Li, H.; Kwak, L.W.; Yi, Q. Roles of idiotype-specific t cells in myeloma cell growth and survival: Th1 and CTL cells are tumoricidal while Th2 cells promote tumor growth. Cancer Res. 2008, 68, 8456–8464. [Google Scholar] [CrossRef] [PubMed]
  34. Muraille, E.; Leo, O.; Moser, M. TH1/TH2 paradigm extended: Macrophage polarization as an unappreciated pathogen-driven escape mechanism? Front. Immunol. 2014, 5, 603. [Google Scholar] [CrossRef] [PubMed]
  35. Corrigan, C.J.; Wang, W.; Meng, Q.; Fang, C.; Wu, H.; Reay, V.; Lv, Z.; Fan, Y.; An, Y.; Wang, Y.; et al. T-helper cell type 2 (Th2) memory T cell-potentiating cytokine IL-25 has the potential to promote angiogenesis in asthma. Proc. Natl. Acad. Sci. USA 2011, 108, 1579–1584. [Google Scholar] [CrossRef] [PubMed]
  36. Grotz, T.E.; Jakub, J.W.; Mansfield, A.S.; Goldenstein, R.; Enninga, E.A.L.; Nevala, W.K.; Leontovich, A.A.; Markovic, S.N. Evidence of Th2 polarization of the sentinel lymph node (SLN) in melanoma. Oncoimmunology 2015, 4, e1026504. [Google Scholar] [CrossRef] [PubMed]
  37. Dulos, J.; Carven, G.J.; van Boxtel, S.J.; Evers, S.; Driessen-Engels, L.J.A.; Hobo, W.; Gorecka, M.A.; de Haan, A.F.J.; Mulders, P.; Punt, C.J.A.; et al. PD-1 blockade augments Th1 and Th17 and suppresses Th2 responses in peripheral blood from patients with prostate and advanced melanoma cancer. J. Immunother. 2012, 35, 169–178. [Google Scholar] [CrossRef] [PubMed]
  38. Sun, W.; Li, W.; Fu, Q.; Wu, C.; Lin, J.; Zhu, X.; Hou, W.; Wei, Y.; Wen, Y.; Wang, Y.; et al. Functionally distinct subsets of CD4+ regulatory T cells in patients with laryngeal squamous cell carcinoma are indicative of immune deregulation and disease progression. Oncol. Rep. 2015, 33, 354–362. [Google Scholar] [CrossRef] [PubMed]
  39. Schreck, S.; Friebel, D.; Buettner, M.; Distel, L.; Grabenbauer, G.; Young, L.S.; Niedobitek, G. Prognostic impact of tumour-infiltrating Th2 and regulatory T cells in classical Hodgkin lymphoma. Hematol. Oncol. 2009, 27, 31–39. [Google Scholar] [CrossRef]
  40. Rui, X.; Shao, S.; Wang, L.; Leng, J. Identification of recurrence marker associated with immune infiltration in prostate cancer with radical resection and build prognostic nomogram. BMC Cancer 2019, 19, 1179. [Google Scholar] [CrossRef]
  41. Qi, X.; Gurung, P.; Malireddi, R.K.S.; Karmaus, P.W.F.; Sharma, D.; Vogel, P.; Chi, H.; Green, D.R.; Kanneganti, T. Critical role of caspase-8-mediated IL-1 signaling in promoting Th2 responses during asthma pathogenesis. Mucosal Immunol. 2017, 10, 128–138. [Google Scholar] [CrossRef] [PubMed]
  42. Ando, M.; Kawazu, M.; Ueno, T.; Fukumura, K.; Yamato, A.; Soda, M.; Yamashita, Y.; Choi, Y.; Yamasoba, T.; Mano, H. Cancer-associated missense mutations of caspase-8 activate nuclear factor-κB signaling. Cancer Sci. 2013, 104, 1002–1008. [Google Scholar] [CrossRef] [PubMed]
  43. Yang, Y.; Yoo, H.M.; Choi, I.; Pyun, K.H.; Byun, S.M.; Ha, H. Interleukin 4-induced proliferation in normal human keratinocytes is associated with c-myc gene expression and inhibited by genistein. J. Investig. Dermatol. 1996, 107, 367–372. [Google Scholar] [CrossRef]
  44. Ciszewski, W.M.; Wagner, W.; Kania, K.D.; Dastych, J. Interleukin-4 enhances PARP-dependent DNA repair activity in vitro. J. Interferon Cytokine Res. 2014, 34, 734–740. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Association between immune subtype and outcome in HPV-negative HNSCC treated with radiation. (AC) Forest plot of univariable analysis for overall survival (A), distant metastasis (B), and locoregional recurrence (C) following surgery and adjuvant radiation. Hazard ratios and 95% confidence intervals are shown for each clinical outcome listed. (D) Kaplan–Meier curve showing locoregional recurrence as a function of high (upper two tertiles) versus low (lower tertile) Th2 infiltrate in the discovery (n = 94) and validation cohorts (n = 97). The Cox PH model was used in (AC). The log-rank test was used in (D).
Figure 1. Association between immune subtype and outcome in HPV-negative HNSCC treated with radiation. (AC) Forest plot of univariable analysis for overall survival (A), distant metastasis (B), and locoregional recurrence (C) following surgery and adjuvant radiation. Hazard ratios and 95% confidence intervals are shown for each clinical outcome listed. (D) Kaplan–Meier curve showing locoregional recurrence as a function of high (upper two tertiles) versus low (lower tertile) Th2 infiltrate in the discovery (n = 94) and validation cohorts (n = 97). The Cox PH model was used in (AC). The log-rank test was used in (D).
Cancers 16 01586 g001
Figure 2. Tumor alterations associated with derived Th2 infiltrate. (A) Volcano plot of mutations in HNSCC and Th2 infiltrate from the full HPV-negative TCGA HNSCC cohort with mutational and gene expression data (n = 406). Mutations exhibiting significantly altered Th2 infiltrate (unadjusted p-value and fold change (mean mutant/mean wild type) are shown, with orange indicating that the mutation is associated with increased infiltrate). (B) Ingenuity pathway analysis of genes associated with Th2 from the full TCGA Head and Neck cohort (n = 409), including z-scores and B-H adjusted p-values.
Figure 2. Tumor alterations associated with derived Th2 infiltrate. (A) Volcano plot of mutations in HNSCC and Th2 infiltrate from the full HPV-negative TCGA HNSCC cohort with mutational and gene expression data (n = 406). Mutations exhibiting significantly altered Th2 infiltrate (unadjusted p-value and fold change (mean mutant/mean wild type) are shown, with orange indicating that the mutation is associated with increased infiltrate). (B) Ingenuity pathway analysis of genes associated with Th2 from the full TCGA Head and Neck cohort (n = 409), including z-scores and B-H adjusted p-values.
Cancers 16 01586 g002
Figure 3. The DDR and TME scores are significantly associated with derived Th2 infiltrate in multiple clinical cohorts. BRCA1 in DDR and TME scores (as calculated in the Materials and Methods section) and derived Th2 infiltrate in the full HPV-negative TCGA cohort with available gene expression data (n = 409) (A), nonoverlapping HPV-negative HNSCC patients from the International Cancer Genome Consortium (ICGC) data (n = 40) (B), and the National Cancer Institute’s Clinical Proteomic Tumor Analysis Consortium (CPTAC) Head and Neck cohort (n = 110) (C). Correlation coefficient and p-value were computed by Spearman’s correlation in (AC). p-value < 0.05 indicates the correlation coefficient is significantly different from zero.
Figure 3. The DDR and TME scores are significantly associated with derived Th2 infiltrate in multiple clinical cohorts. BRCA1 in DDR and TME scores (as calculated in the Materials and Methods section) and derived Th2 infiltrate in the full HPV-negative TCGA cohort with available gene expression data (n = 409) (A), nonoverlapping HPV-negative HNSCC patients from the International Cancer Genome Consortium (ICGC) data (n = 40) (B), and the National Cancer Institute’s Clinical Proteomic Tumor Analysis Consortium (CPTAC) Head and Neck cohort (n = 110) (C). Correlation coefficient and p-value were computed by Spearman’s correlation in (AC). p-value < 0.05 indicates the correlation coefficient is significantly different from zero.
Cancers 16 01586 g003
Table 1. Tumor characteristics of the HPV-negative squamous cell carcinoma discovery cohort.
Table 1. Tumor characteristics of the HPV-negative squamous cell carcinoma discovery cohort.
Nodal StagenPercent
04446.8
12021.3
2x11.1
2a22.1
2b1313.8
2c88.5
322.1
x44.3
Tumor stage
111.1
21212.8
32324.5
45659.6
x22.1
Site
Oral cavity6063.8
OPX44.3
Larynx/hypopharynx3031.9
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Abdelhakiem, M.K.; Bao, R.; Pifer, P.M.; Molkentine, D.; Molkentine, J.; Hefner, A.; Beadle, B.; Heymach, J.V.; Luke, J.J.; Ferris, R.L.; et al. Th2 Cells Are Associated with Tumor Recurrence Following Radiation. Cancers 2024, 16, 1586. https://doi.org/10.3390/cancers16081586

AMA Style

Abdelhakiem MK, Bao R, Pifer PM, Molkentine D, Molkentine J, Hefner A, Beadle B, Heymach JV, Luke JJ, Ferris RL, et al. Th2 Cells Are Associated with Tumor Recurrence Following Radiation. Cancers. 2024; 16(8):1586. https://doi.org/10.3390/cancers16081586

Chicago/Turabian Style

Abdelhakiem, Mohamed K., Riyue Bao, Phillip M. Pifer, David Molkentine, Jessica Molkentine, Andrew Hefner, Beth Beadle, John V. Heymach, Jason J. Luke, Robert L. Ferris, and et al. 2024. "Th2 Cells Are Associated with Tumor Recurrence Following Radiation" Cancers 16, no. 8: 1586. https://doi.org/10.3390/cancers16081586

APA Style

Abdelhakiem, M. K., Bao, R., Pifer, P. M., Molkentine, D., Molkentine, J., Hefner, A., Beadle, B., Heymach, J. V., Luke, J. J., Ferris, R. L., Pickering, C. R., Wang, J. H., Patel, R. B., & Skinner, H. D. (2024). Th2 Cells Are Associated with Tumor Recurrence Following Radiation. Cancers, 16(8), 1586. https://doi.org/10.3390/cancers16081586

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop