Next Article in Journal
Reflectance Confocal Microscopy and Dermoscopy of Facial Pigmented and Non-Pigmented Actinic Keratosis Features before and after Photodynamic Therapy Treatment
Next Article in Special Issue
Immune Escape Strategies in Head and Neck Cancer: Evade, Resist, Inhibit, Recruit
Previous Article in Journal
The Role of Gene Fusions in Thymic Epithelial Tumors
Previous Article in Special Issue
KEAP1/NRF2 Mutations in Stem Cells Define an Aggressive Subset of Head and Neck Cancer Patients Who Have a Poor Prognosis, Lung Metastasis, and Therapeutic Failure
 
 
Correction published on 21 October 2024, see Cancers 2024, 16(20), 3554.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identifying Predictive Biomarkers for Head and Neck Squamous Cell Carcinoma Response

1
Institute of Pathology, Rostock University Medical Center, 18057 Rostock, Germany
2
Department of Otorhinolaryngology, Head and Neck Surgery “Otto Koerner”, Rostock University Medical Center, 18057 Rostock, Germany
3
Department of Internal Medicine B, Cardiology, University Medicine Greifswald, 17475 Greifswald, Germany
4
Institute of Pathology, University of Luebeck, University Hospital Schleswig-Holstein, Campus Luebeck, 23538 Luebeck, Germany
5
Department of Oto-Rhino-Laryngology & Head and Neck Surgery, University of Lubeck, University Hospital Schleswig-Holstein, Campus Luebeck, 23538 Luebeck, Germany
6
Department of Internal Medicine, Medical Clinic III—Hematology, Oncology, Palliative Medicine, Rostock University Medical Center, 18057 Rostock, Germany
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(23), 5597; https://doi.org/10.3390/cancers15235597
Submission received: 19 October 2023 / Revised: 20 November 2023 / Accepted: 21 November 2023 / Published: 27 November 2023 / Corrected: 21 October 2024
(This article belongs to the Special Issue Head and Neck Cancers—Novel Approaches and Future Outlook)

Abstract

:

Simple Summary

Head and neck cancer patients have a poor prognosis, even with the best therapeutic options. This is partly due to the lack of reliable prognostic biomarkers. In this study, we aimed to unravel the prognostic and predictive value of immune-related biomarkers in the tumor microenvironment. We identified the PD-L1-stabilizing and membrane-associated protein CMTM6 (CKLF-like MARVEL transmembrane domain-containing protein) as beneficial for patients receiving radiotherapy. Co-expression with the regulatory T cell transcription factor FOXP3 was predictive of a response to radiotherapy. However, no association with survival was seen for other treatment regimens, including chemotherapy, radio-chemotherapy, and immunotherapy. In the latter, CTLA-4 was associated with worse outcomes, whereas PD-L1 was not predictive. Therefore, this study describes a strategy to identify predictive markers for treatment response in head and neck cancer patients and highlights the need for treatment-specific biomarker screening.

Abstract

The 5-year survival rate for head and neck squamous cell carcinoma (HNSCC) is approximately 65%. In addition to radio-chemotherapy, immunotherapy is an approach in the treatment of advanced HNSCC. A better understanding of the immune context would allow personalized treatment by identifying patients who are best suited for different treatment options. In our discovery cohort, we evaluated the expression profiles of CMTM6, PD-L1, CTLA-4, and FOXP3 in 177 HNSCCs from Caucasian patients of all tumor stages and different treatment regimens, correlating marker expression in tumor and immune cells with outcomes. Patients with CMTM6high-expressing tumors had a longer overall survival regardless of treatment. This prognostic benefit of CMTM6 in HNSCC was validated in an independent cohort. Focusing on the in the discovery cohort (n = 177), a good predictive effect of CMTM6high expression was seen in patients receiving radiotherapy (p = 0.07; log rank), but not in others. CMTM6 correlated with PD-L1, CTLA-4 and FOXP3 positivity, with patients possessing CMTM6high/FOXP3high tumors showing the longest survival regardless of treatment. In chemotherapy-treated patients, PD-L1 positivity was associated with longer progression-free survival (p < 0.05). In the 27 patients who received immunotherapy, gene expression analysis revealed lower levels of CTLA-4 and FOXP3 with either partial or complete response to this treatment, while no effect was observed for CMTM6 or PD-L1. The combination of these immunomodulatory markers seems to be an interesting prognostic and predictive signature for HNSCC patients with the ability to optimize individualized treatments.

1. Introduction

Head and neck squamous cell carcinoma (HNSCC) is a highly heterogeneous and common malignancy, with an increasing prevalence in the Western world. The majority of HNSCCs diagnosed worldwide are in men [1,2,3]. The major risk factors include chronic alcohol and tobacco use and high-risk human papillomavirus (HPV) infection [4,5]. The standard of care includes surgery, chemotherapy, and radiotherapy either alone or in combination with EGFR-targeting and/or immunotherapy. Although many patients initially respond well to platinum-based radio-chemotherapy, resistance and thus treatment refractoriness frequently occur. The heterogeneity of the disease and the associated individual patient response also make it difficult to standardize treatment protocols and estimate prognosis.
To date, HPV status is the most reliable prognostic factor [6,7], but other aspects, such as molecular and pathological characteristics, tumor location, and clinical stage, additionally influence patient outcomes. The introduction of the PD1-targeting antibodies pembrolizumab and nivolumab into the first- and second-line treatment of patients with recurrent and/or metastatic (r/m) HNSCC has raised hopes for improved outcomes, but objective response rates have been reported to be less than 20%. As a result, the prognosis has not improved significantly over the decades, and the five-year overall survival rate remains at 50–60% [8,9,10,11].
Positivity for programmed death-ligand 1 (PD-L1) is the only approved biomarker for immune checkpoint inhibition [12,13,14]. However, the limited predictive value is a major pitfall limiting immunotherapy success. We have recently identified CKLF-like MARVEL transmembrane domain 6 (CMTM6) as a novel prognostic biomarker with, so far, unknown predictive value. CMTM6 stabilizes PD-L1 on the tumor cells’ surface by inhibiting its ubiquitination and subsequent degradation via lysosomes, exhibiting opposing functions across tumors [15,16,17,18,19]. In colorectal cancer, high CMTM6 levels are associated with an inflamed tumor microenvironment, while in hepatocellular carcinoma, CMTM6 inhibits cell proliferation by preventing p21 ubiquitination [15,20]. In HNSCC, CMTM6 was previously reported to mediate cisplatin resistance and regulate AKT-mTORC1-dependent ribosome biogenesis [21,22]. Furthermore, tumor-secreted exosomal CMTM6 was reported to induce M2-like macrophage polarization, indicating an active crosstalk between tumor and immune cells [23]. To our knowledge, studies examining the impact of CMTM6 expression in HNSCC patients receiving immunotherapy are lacking.
Given this, and the fact that the mainstay of treatment for HNSCC patients includes platinum-based (radio-)therapy or immune checkpoint inhibition, the identification of reliable predictive biomarkers remains an unmet need. Increasing evidence suggests that immune status may differ between responders and non-responders. In a very recent report, a gene signature related to the inflammatory response was described, and OLR1 (a key receptor for oxidized low-density lipoprotein) and INHBA (a member of the transforming growth factor-β superfamily) were strongly associated with a poor prognosis of HNSCC patients [24]. In addition, a risk prognostic signature for predicting radiotherapy response based on the expression levels of therapeutic response-related genes has just been reported by Lin et al. [25]. In this study, we performed a detailed analysis using RNA and protein profiling of the tumor microenvironment obtained from primary and recurrent and/or metastatic HNSCC cases. We identified this marker as beneficial for patients receiving radiotherapy. Additionally, we confirmed the prognostic value of CMTM6 in an independent cohort.

2. Materials and Methods

2.1. Patient Cohort and Multi Tumor Tissue Blocks

Pretreatment samples from locally advanced and/or r/m HNSCC patients, partly treated with immune checkpoint inhibitors (ICI) (pembrolizumab, nivolumab), diagnosed between 2017 and 2021, were collected from the Rostock Pathology archives to establish our discovery cohort. If available, additional pre/post ICI treatment tumor biopsies were included. The patients’ clinicopathological characteristics, including response to treatment according to the Response Evaluation Criteria in Solid Tumors (RECIST) version 1.1 and progression-free survival (PFS), were obtained from clinical records. Using pseudonyms, follow-up data were obtained from the regional cancer registries. Tumor tissue was obtained from formalin-fixed, paraffin-embedded (FFPE) specimens. Tissue microarrays (TMAs) of three representative areas were constructed using 5 mm cores for samples from ICI-treated patients and 1 mm cores for all non-ICI-treated patients using a tissue chip microarrayer (Beecher Instruments, Silver Spring, MD, USA). The assembly of 5 mm cores preserved the tumor center and the invasion front, including infiltrating cells, for analysis and was performed with biopsy punches (Stiefel Laboratorium, Burgdorf, Germany). Tissues were collected with the patients’ consent. The institutional ethic committee at the University Hospital Rostock approved the study (number A2022-0120), which was conducted in accordance with the Declaration of Helsinki of 1975. An independent, well-described cohort from primary HNSCCs was used as a control group to validate our prognostic findings [26,27,28].

2.2. Gene Expression Analysis

For RNA extraction, after reviewing all original tissue slides, FFPE tissue blocks were selected and recut for hematoxylin–eosin-stained sections to determine tumor surface area. After deparaffinization, RNA was extracted from unstained FFPE tumor sections using the RNeasy FFPE kit (Qiagen, Hilden, Germany). Microdissection was performed in cases showing <50% tumor content or inked borders. mRNA expression profiling was performed using the NanoString nCounter gene expression platform (NanoString Technologies, Seattle, WA, USA). To evaluate cancer-related genes, we applied the Nanostring™ PanCancer IO 360TM panel comprising 770 targets according to manufacturer’s instructions. In short, a hybridization reaction was performed using 100 ng of RNA per sample followed by further proceeding in the fully automated Prep Station with data acquisition conducted using the Digital Analyzer (performed at the Department of Internal Medicine B, Cardiology, University Medicine Greifswald, Germany). The data were exported as reporter code count (RCC) files and imported to the NanoString nSolver™ analysis software v4.0 for further analysis. Automatic quality control of mRNA was performed according to the software’s instructions.

2.3. Immunohistochemistry and Microscopic Evaluation

Sections of 4 μm were used. Heat-induced antigen retrieval was performed with a high pH buffer (20 min at 97 °C). The following steps were performed in an Autostainer link 48 instrument (Dako, Hamburg, Germany): 5 min of incubation in peroxidase-blocking buffer followed by 20 min of incubation with primary antibody (dilution and clone per antigen—AKT: 1:500; EPR1798, p-AKT (phosphor-S473): 1:50; EP2109Y, CMTM6: 1:1000; EPR23015-45, FOXP3: 1:200; 20034, all Abcam (Cambridge, UK); cytotoxic T-lymphocyte-associated protein 4 (CTLA-4): 1:200; UMAB249 (Rockville, MD, USA), D4B9C, PDL1: 1:100; 22C3, ß2-M: 1:800, A0072; all Dako) and 3,3′- diaminobenzidine (DAB) detection using the Dako-kit K8000 according to the manufacturer’s instructions. Slides were counterstained with hematoxylin. Appropriate positive and negative controls were used.
For assessment of immune cell distribution status (Figure S1), cytoplasmic and/or membranous reactivity of CTLA-4 and ß2 m as well as nuclear staining for FOXP3 were evaluated. Minimum and maximum expression per marker within all samples defined the ranges of the four applied categories (0: absent; 1: minimal; 2: moderate; 3: abundant), with 0/1 defined as low and 2/3 as high. Cytoplasmic and/or membranous reactivity for AKT and pAKT (S473) was scored semi-quantitatively by the percentage of positive cells and the staining intensity resulting in no/low- vs. high-expressing cases. Scores were built counting 3-5 high-power fields in the tumor center with an Axio-Cam 205 color microscope (field diameter 0.5 mm, Zeiss, Oberkochen, Germany).
PD-L1 and CMTM6 expression were scored analogous to the combined positive score (CPS) as described previously [29]. All scores were calculated within an area showing at least 100 viable tumor cells. Depending on PD-L1 and CMTM6 expression, samples were categorized as PD-L1-positive vs. -negative (cut off ≥ 1: positive) and CMTM6-low vs. -high (cut off ≥ 10: high).

2.4. Statistics

Based on D’Agostino, Pearson and Shapiro–Wilk tests for normality, Pearson’s or Spearman’s correlation coefficients (r) were used to analyze the linear association between two continuous variables and used the non-parametric t-test Mann–Whitney tests to compare groups. Overall survival (OS) and PFS curves were constructed using the Kaplan–Meier analysis, and statistical significance was determined using the log-rank test. Unadjusted univariable analysis was conducted. All statistical tests were two-sided, and p values below 0.05 were considered statistically significant. All statistical analyses were performed using GraphPadTM Prism® v7.0 software (supplier) and nSolver analysis software (version 4.0). After removing samples with the quality control (QC) flag, mRNA count normalization and log2 Fold Change (FC) calculation with Welch’s t-test p values were performed. Significant DEGs were defined by p < 0.05 and |log2FC| > 1 as the threshold.

3. Results

3.1. Patient and Tumor Characteristics

The discovery cohort included a total of 177 patients diagnosed with r/m HNSCC with a median age of 66 years (±9.8). Clinicopathologic characteristics are shown in Table 1, where stage refers to the time of treatment initiation. Approximately 45% of all patients were smokers (≥10 py) and >32% reported regular alcohol consumption. Thus, most of the cases included here were noxae-associated, either localized in the oral cavity or oropharynx. A total of 18 patients suffered from hypopharyngeal cancer and 15 from laryngeal cancer. One regional lymph node metastasis as well as four distant metastases were obtained, while the remaining samples in the cohort consisted of locally resected primary/recurrent tumors. A total of 20% of all cases were due to a previous high-risk HPV infection.
After diagnosis, patients received either cisplatin-based chemotherapy (n = 66), radiotherapy (n = 30), combined radio-chemotherapy (n = 13), ICI (n = 27), or ICI plus chemotherapy (n = 3). A small subset of patients were treated with other agents (n = 16) or received best supportive care only due to reduced general conditions (n = 22). Median PFS on treatment was 21 months and the OS was 29 months. PFS in the ICI-treated subgroup ranged from 1 to 18 months (median 4). Using the RECIST criteria, 1 patient showed a complete response, and 10 patients showed a partial response (=responders), while 5 patients were classified as mixed response/stable disease and the remaining 10 developed tumor progression.
A validation cohort of 286 treatment-naïve primary HNSCCs derived from Caucasian patients with comparable clinicopathological parameters was included in the initial survival analysis. Anonymyzed clinical information on this patient cohort was obtained from the cooperating institute (University of Lübeck) (Figure 1). Ref. [28] provides precise information about the validation cohort. Kaplan–Meier plots for the validation cohort stratified by CMTM6 are shown in Supplementary Figure S2.

3.2. CMTM6 Status Predicts Response to RCT but Not to Anti-PD1 Treatment

In a total of 463 patients from both cohorts combined, the mean CMTM6 CPS was 33.7 ± 32.3. PD-L1 staining was only available for the discovery cohort (n = 177), and mean PD-L1 CPS was 14.9 ± 26.2 (Table 2). Notably, 56 cases (31%) were either classified as “negative” (PD-L1 CPS < 1) or “low” (CMTM6 CPS < 10). In >40% of these cases, PD-L1 negativity correlated with low CMTM6 status.
Both markers had prognostic value, i.e., patients with PD-L1 positivity or high CMTM6 status showed a significantly better overall survival (Figure 2A,B left, p < 0.05). Hence, the prognostic value of PD-L1 and CMTM6 was formally confirmed in both patient cohorts. Stratification by HPV status additionally confirmed a significant prognostic impact for CMTM6 but not for PD-L1 (Figure 2A,B, middle and right).
We then examined the predictive value of these immunologic markers in n = 177 patients from the discovery cohort (Figure 2C,D). We focused on cisplatin, RT, RCT, and ICI to cover the most relevant treatment regimens. A predictive value for PD-L1 was only seen in the cisplatin subgroup. Here, patients with a CPS ≥ 1 showed a significantly better OS than their counterparts with CPS < 1 (Figure 2B). No statistically significant benefit was seen for the other regimens. For CMTM6, there was a trend towards a better outcome when the tumor was classified as high and the patient received RT (p = 0.07). Likewise, a slightly improved overall survival was seen in the other treatment groups, i.e., cisplatin or RCT. However, patients with high PD-L1 and/or CMTM6 tumors who received anti-PD1 ICI treatment had a comparable OS to their counterparts (Figure 2B,C). Hence, no predictive impact of either marker was seen in this subgroup.

3.3. Correlation of CMTM6 Status with Other Immune Biomarkers

The interplay between different immune cell subtypes within the tumor microenvironment (TME) prompted us to check for additional biomarkers with potential prognostic and predictive value (Spearman correlation) in the discovery cohort. For this analysis, treatment-unrelated (=prognostic) and treatment-related (=predictive) correlations were conducted (Table 3). For the latter, only patients who received cisplatin, RT, or RCT were included. CMTM6 was analyzed for correlation with PD-L1, CTLA-4, FOXP3, AKT, and pAKT. Table 2 and Table 3 show the staining results and the correlations between the markers. Table 4 refers to staining pattern depending on the clinical parameters.
Independent of treatment, CMTM6 expression strongly correlated with FOXP3, a specific transcriptional regulator in Tregs. Moderate correlations were observed between CMTM6 and PD-L1 and CTLA-4. Notably, a strong correlation between CMTM6 and FOXP3 was seen across treatments. In all subgroups, CMTM6 correlated more strongly with CTLA-4 status than PD-L1. No association was found between CMTM6 and AKT or pAKT.
Then, the prognostic role of these immune markers was analyzed (Figure 3). High FOXP3 expression levels were associated with significantly improved overall survival in both HPVpos and HPVneg cases (Figure 3A, p < 0.05). To study if FOXP3 has comparable prognostic relevance to CMTM6, the overall survival of patients was analyzed in additional subgroups. Therefore, tumors were classified into FOXP3low/CMTM6low (n = 51), FOXP3high/CMTM6low (n = 6), FOXP3low/CMTM6high (n = 17), and FOXP3high/CMTM6high (n = 99). Patients in the latter group and those with FOXP3high/CMTM6low tumors had the best prognosis with survival rates exceeding 50% at the longest time point. This underlines the prognostic impact of both markers in HNSCC.
To see if this prognostic role could be transferred to the predictive model, we performed a subgroup analysis with respect to the individual treatments. A trend towards a better outcome was observed for all treatment regimens (Figure 3B), reaching statistical significance for patients receiving RCT (p < 0.05 vs. FOXP3low/CMTM6low). A closer look at the survival curve of ICI-treated patients revealed a delayed favorable response with a median OS of 43 months in the FOXP3high/CMTM6high subgroup (vs. 22 months FOXP3low/CMTM6low). Thus, although statistical significance was not achieved in this small patient cohort, we confirmed the predictive value by doubling the median OS of patients with r/m HNSCC.
CTLA-4 had no effect on median OS. Finally, the relevance of serine/threonine protein kinase AKT and its active phosphorylated form (pAKT) was analyzed, because AKT is a downstream target of CMTM6-mediated Wnt signaling and may be involved in cisplatin resistance. Hence, AKT is considered as a proto-oncogene. In our study, we identified a slightly improved median OS in the AKT high subgroup (vs. AKT low) and vice versa when pAKT was considered. Still, no significant survival benefit was seen in either group.
Further treatment-related analysis was only performed for the cisplatin subgroup but did not yield any prognostic value. Hence, the prognostic value of AKT/pAKT in our patient cohort is rather low.

3.4. Immunotherapy Response in r/m HNSCC Gene Signatures in Responders

Transcriptomic and immunological profiling was performed on pre- and post-treatment samples of r/m HNSCC patients receiving anti-PD1 immunotherapy. These samples were from the discovery cohort and only included one patient with high-risk oncogenic HPV infection; the remaining samples were HPV-negative. Gene expression analysis for anti-tumor immune activity showed two distinct subgroups (Figure 4A). Concerning the profiled gene signatures, responders had significantly lower mRNA expression levels of CTLA-4 than non-responders (Figure 4B). By plotting the receiver operator characteristics (ROCs), CTLA-4 had the highest capacity to predict response (area under the curve 0.753; PD-L1: 0.7) of 43 signatures. While the expression levels of genes associated with T cells, including exhausted CD8+ cytotoxic cells, as well as levels of the immune checkpoint receptor T cell immunoglobulin and ITIM domain (TIGIT) were lower in responders, this failed to be statistically significant in the multivariate analysis (Figure 4C). In line with our findings described above, the abundance of Tregs as measured by gene expression of FOXP3 was lower in these immunologically cold tumors, here showing a better treatment response to ICI, and vice versa. Patients benefiting longer from ICI had lower expression genes levels in tumors involved in signatures of CD8 T cells, TIGIT, and CTLA-4 compared to patients with shorter progression times (Figure 4D). Interestingly, we found no differences in PD-L1 expression between responders and non-responders. Genes involved in the pathway downstream of the IFN gamma axis were significantly higher in samples defined as “PD-L1-positive” by IHC (Figure 4E). In three patients with re-biopsy material at the time point of progression during ICI treatment, neither CMTM6 nor PD-L1 status by IHC changed according to the applied cut-offs, and Beta-2-microglobulin expression was preserved. Focusing on the profiled immune cell population in these patients (one classified as a responder, two as non-responders), the tumor inflammation signature (TIS), which predicts a response to pembrolizumab, decreased in the two non-responders (Figure 4F). On the other hand, CTLA-4 levels increased in the third patient, who progressed after 8 months of a partial response (Figure 4G). We also focused on the seven included patients with available tumor material (a) either before or (b) at disease progression during ICI treatment. Using IHC, PD-L1, CMTM6, FOXP3 and CTLA-4 status was stable in two patients with mixed changes in the remaining patients. Immune profiling generally revealed inhibition of immune cell abundance and antitumor immune activity; however, an increase in inhibitory immune signaling was observed when primary tumors were compared with the material from r/m disease, as previously described [28].

4. Discussion

The identification of reliable biomarkers for HNSCC is ongoing. In this study, we applied a classical pathological approach and examined the TME of HPV-related and -unrelated cases to elucidate the biological relevance of several known and novel biomarkers. We included cases from different sites and treatment regimens to cover the clinical heterogeneity and related challenges. The overall outcome was comparable to reports in the literature and confirmed the beneficial response to combined radio-chemotherapy compared to either monotherapy in the first-line treatment of HNSCC. However, in the r/m setting, the use of ICIs has little impact on PFS or OS with frequent relapse in all patients.
In our discovery cohort, we observed that CMTM6 positivity correlates with an improved OS, independently of adjuvant treatment modalities. We observed this beneficial prognostic impact in an independent validation cohort from treatment-naïve HNSCC patients with comparable clinical characteristics. CMTM6 is now well known to stabilize PD-L1 protein expression on the tumor cells’ surface [19,30]. It is widely expressed on both tumor and immune cells and plays an important role in regulating T cell activation and antitumor response. The prognostic impact of CMTM6 is comparable to PD-L1, underlying its overlapping biological functions. However, the predictive value of this novel biomarker is largely unknown [31]. By stratifying patients from the discovery cohort according to the most frequent treatment regimens, i.e., chemotherapy, radiotherapy (RT), radio-chemotherapy, or immune checkpoint inhibition, a predictive value was only seen in the RT group. Patients with high CMTM6 expression on tumor cells and/or in the TME showed a survival benefit compared to their CMTM6-low counterparts. While HPV is an approved biomarker, the prognostic impact was even higher here when tumors were also CMTM6-high. The prognostic impact of CMTM6 is in accordance with reports on colorectal (CRC) and hepatocellular carcinomas [15,20] but in contrast to the findings of Chen et al. [32,33]. In their cohort of 210 primary HNSCC samples collected from the School and Hospital of Stomatology, Wuhan University, China, the authors identified CMTM6 overexpression by IHC as a poor prognostic marker. Using digital image analysis, membranous and/or cytoplasmic CMTM6 positivity in tumor cells predicted a shorter OS independent of other clinicopathological parameters. Their scoring system for CMTM6, analogous to the tumor proportion score (TPS), excluded immune cells surrounding the tumor. In our study, we quantified CMTM6 positivity using the CPS. This score includes cells within the TME, which regulates tumor cell survival and confers anticancer effects [3,32]. Although the group describes similar staining patterns for PD-L1 and CMTM6 as seen here, the results are not comparable due to the different scoring methods. Furthermore, Chen et al. performed their study in an Asian cohort, where risk factors and outcomes differ from those in Caucasians [33]. With this in mind, the authors reviewed The Cancer Genome Atlas (TCGA) [34] to support their findings. While they do not provide information on the impact of CMTM6 RNA expression in the TCGA cohort, they show that both CMTM6 RNA expression and protein levels, when analyzed analogously to TPS, increased with higher pathological grade. Our cohort includes 21% HPV-related and therefore conventionally ungraded tumors. We did not see any correlation between tumor grade and CMTM6. However, another Chinese group investigated CMTM6 RNA levels in CRC samples and found a higher CMTM6 expression in early clinical stages (UICC I/II vs. III/IV) [15] (ref. Peng). Although CRC as an adenocarcinoma is not comparable to HNSCC, this underlies the importance of future studies to elucidate the impact of CMTM6 in solid cancers.
Still, no such correlation was seen in the other subgroups. The most unexpected finding was the missing predictive value for the ICI subgroup, which is described for lung carcinoma [35]. The same missing correlation was seen for PD-L1. Here, a beneficial effect was only obtained when patients received cisplatin but not with other therapies. Although this finding is somehow unexpected, it fits with recent reports in which ICI-based treatments yielded mixed responses and sometimes failed to improve the outcome of r/m HNSCC patients [36,37,38]. Hence, PD-L1 and/or CMTM6 alone may be insufficient to predict ICI responses in this patient cohort.
Here, we performed a more global approach by combining protein and gene expression-based analysis in patients from the discovery cohort. Firstly, we identified a strong correlation between CMTM6 and FOXP3. The latter is a specific transcription factor in regulatory T cells with known immunosuppressive effects and a correspondingly conflicting role in patient prognosis. In HNSCC, high FOXP3 levels have previously been reported to correlate with better outcomes [39,40] regardless of treatment [41]. Here, we further confirmed the prognostic impact independent of HPV status. Comparable results were shown in gastric cancer, where co-expression of FOXP3 and CMTM6 was related to a favorable prognosis [42]. FOXP3 was found to co-localize with tumor-infiltrating CD20+CD27+ B cells to promote effector and memory T cell differentiation and enhance B cell and NK cell activation and function. Therefore, the dual role of FOXP3 in cancer is advantageous for HNSCC patients, and this positive impact is probably enhanced by CMTM6. Adding to this, Jing et al. recently described an inflammatory response-related gene profile to predict the prognosis in HNSCC [24]. In a parallel study, a risk model was developed to predict therapy response and prognosis in response gene-defined subgroups of HNSCC [25]. With our study, we suggest additional predictive biomarkers that may be potential targets in the near future.
Another novel finding of our study was the strong correlation between CMTM6 and CTLA-4, which was even stronger than with PD-L1. This is intriguing and describes a new, previously unknown interaction of CMTM6 with immunoregulatory molecules [19]. Unlike CMTM6, CTLA-4 was not expressed in tumor cells, as reported recently by Hoffmann et al. [43]. In this study, lower CTLA4 promoter methylation correlated with a response to ICI and longer progression-free survival, but these findings were not confirmed at the protein level. This adds to Sholl’s conclusion that not only protein but also genomic and transcriptomic levels mediate the variable response to ICIs [44].
The interaction between the tumor and host immune-specific genetic and epigenetic factors is complex and it is now clear that PD-L1 protein expression by immunohistochemistry is neither sensitive nor specific enough to serve as a reliable predictive biomarker. Our study adds to this by showing that protein and gene expression levels of PD-L1 were insufficient predictors of a response to immunotherapy. Beyond PD-L1, the transcriptomic profile of this entity can identify possible correlations between an ICI response, TME, and CMTM6. By means of gene expression analysis, we found that the mRNA levels of CTLA-4 and FOXP3 were lower in the responder group of our patients treated with ICI. Since CTLA-4 is mainly expressed in tumor-draining lymph nodes, established tissue-based biomarkers for predicting its inhibitory capacity are lacking. In melanoma patients, the anti-CTLA-4 inhibitor ipilimumab is the standard of care regardless of PD-L1 status, and the combination with an anti-PD1 ICI improves its efficacy. Although some studies have failed to show a beneficial effect of the Nivo/Ipi combination, several phase II/III trials are currently investigating the efficacy of anti-CTLA-4 antibodies as a monotherapy or in combination with anti-PD1 antibodies in HNSCC (ClinicalTrials.gov Identifiers: NCT04080804, NCT04326257, NCT03624231, NCT03212469, NCT03799445). In addition, a lower CD8 profile indicated a trend towards an improved ICI response, which is in contrast to the data from Ayers et al. [45]. Using mRNA expression profiling, they showed that a T cell-inflamed microenvironment characterized by active IFN-γ signaling, cytotoxic effector molecules, and T cell-active cytokines is a common feature of the biology of tumors that respond to PD1 checkpoint blockade. With this knowledge, the authors established a focused set of genes to identify this PD1 checkpoint blockade-responsive biology. Although the genes analyzed were congruent with those studied here, possible explanations for the different findings may be that (1) our study cohort included less patients (n = 26 here vs. n = 43 by Ayers et al.), (2) only n = 8 patients received ICI as a first-line treatment, (3) the PFS in our cohort was shorter, and (4) PD-L1 CPS dis not predict patients’ responsiveness to treatment.
Combination therapies targeting PD1 and CTLA-4 may also fail due to primary or secondary resistance mechanisms. Beta-2-microglobulin (ß2M) is a major histocompatibility complex class I component and is essential for antigen presentation. Its loss, either by mutation or degradation, leads to ICI treatment failure [46,47]. Similarly, all tumors in the responder group expressed ß2M, whereas three non-responders lacked this protein on the tumor cell surface prior to treatment initiation. The three samples with re-biopsy material at the time of progression during ICI treatment had all retained ß2M expression, supporting alternative resistance mechanisms.
Finally, the relevance of the serine/threonine kinase AKT and its active form, pAKT, was studied. In oral cancer, pAKT is overexpressed because of reduced protein degradation and contributes to tumor growth, lymph node metastasis, and shorter survival time [48,49]. Other studies related pAKT expression to radio-resistance [50] or HPV positivity [51]. In our study, neither AKT nor its active form were associated with a significant survival benefit. This finding is rather expected and hence additional studies on larger patient cohorts will have to clarify the relevance of AKT/pAKT in HNSCC.
A limitation of our study was the small sample size of only n = 26 HNSCCs in the ICI cohort for gene expression profiling. However, as anti-PD1 treatment had only recently been approved for the treatment of HNSCC, the present cohort is of reasonable size, including eight patients without prior systemic therapy. Second, the TMA design weakens our results because separation into the tumor center, invasion front, and stromal compartment is more difficult than on whole mount sections. To address this, we used 1 mm cores in the construction of the TMA and punched tumor tissue multiple times to capture the heterogeneity of the TME. In addition, the biomarkers CMTM6 and PD-L1 were examined on whole slides. We also extracted mRNA from whole tumors, including immune cells. To increase the tumor cell content and decrease bystander cells, tumor microdissection was performed. In challenging cases, we could have used laser microdissection to separate tumor and immune cells, but the risk of losing mRNA-containing cells by focusing on the small amount of tumor in some of our cases seemed too high. Next, the measurement of CMTM6 was not well defined. To include CMTM6 expression on tumor and associated immune cells, we scored CMTM6 analogously to PD-L1 by calculating its CPS, which itself has its weaknesses. Receiver operating curve analysis revealed a cut-off of ≥10 for CMTM6. Future studies in independent cohorts will show whether another scoring system for CMTM6 protein expression is more valuable. Additionally, fluorescent labeling of the examined molecules, including immune cells, would prove the co-expression of CMTM6 with other markers in more detail. While membranous co-localization of CMTM6 and PD-L1 was shown via high-resolution florescence several times [35,52], the association with other co-stimulatory molecules, including FOXP3 and CTLA-4, especially in HNSCC, is less clear. Therefore, this question should be addressed in a future project, maybe even incorporating other co-stimulatory surface molecules like CD28.

5. Conclusions

In Caucasian HNSCC patients, CMTM6high tumors showed a longer OS in independent cohorts. CMTM6 correlates with PD-L1-, CTLA-4-, and FOXP3-positive Treg infiltration in HNSCCs with a better response to radiotherapy in CMTM6high tumors. While neither PD-L1 nor CMTM6 status predicted the treatment response to immune checkpoint inhibitors, and low CTLA-4 expression was associated with a good response. This argues for the inclusion of CTLA-4 inhibitors to target the antitumor immunity in HNSCC patients.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15235597/s1, Figure S1: Staining for immune checkpoint molecules and immune cells in high vs. low expression cases. Figure S2: Progression-free survival and overall survival in the validation cohort stratified by CMTM6 independent of treatment and for subgroups receiving either chemotherapy, radiation or chemoradiation.

Author Contributions

Conceptualization, A.-S.B. and C.M.; data curation, D.S.; formal analysis, A.H.Z. and B.S.; funding acquisition, A.-S.B. and C.M.; methodology, C.K., C.S., B.S. and C.R.; resources, J.R.-I. and C.I.; software, C.M.; writing—original draft, A.-S.B. and C.M.; writing—review and editing, A.-S.B., D.S., J.R.-I., C.I. and C.M. All authors have read and agreed to the published version of the manuscript.

Funding

The project was supported by a grant from the University of Medicine Rostock FORUN Program to ASB (Project no.: 889038).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board (A2022-0120).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We are grateful to Heike Zettl and Kloecking for providing clinical data from the cancer registry Mecklenburg-Western Pomerania.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

Cchemotherapy
CMTM6CKLF-like MARVEL transmembrane domain 6
CPScombined positive score
dday
HNSCChead and neck squamous cell carcinoma
HPVhuman papillomavirus
PD-L1programmed death-ligand 1
ICIimmune checkpoint inhibitors
OSoverall survival
PFS progression-free survival
pypack years
RCTradio-chemotherapy
RTradiotherapy
r/mrecurrent and/or metastatic
TIGITT cell immunoglobulin and ITIM domain

References

  1. Lorenzoni, V.; Chaturvedi, A.K.; Vignat, J.; Laversanne, M.; Bray, F.; Vaccarella, S. The Current Burden of Oropharyngeal Cancer: A Global Assessment Based on GLOBOCAN 2020. Cancer Epidemiol. Biomarkers Prev. 2022, 31, 2054–2062. [Google Scholar] [CrossRef]
  2. Deneuve, S.; Pérol, O.; Dantony, E.; Guizard, A.V.; Bossard, N.; Virard, F.; Fervers, B. Diverging Incidence Trends of Oral Tongue Cancer Compared to Other Head and Neck Cancers in Young Adults in France. Int. J. Cancer 2022, 150, 1301–1309. [Google Scholar] [CrossRef]
  3. Deng, Y.; Wang, M.; Zhou, L.; Zheng, Y.; Li, N.; Tian, T.; Zhai, Z.; Yang, S.; Hao, Q.; Wu, Y.; et al. Global Burden of Larynx Cancer, 1990-2017: Estimates from the Global Burden of Disease 2017 Study. Aging 2020, 12, 2545–2583. [Google Scholar] [CrossRef]
  4. Leoncini, E.; Ricciardi, W.; Cadoni, G.; Arzani, D.; Petrelli, L.; Paludetti, G.; Brennan, P.; Luce, D.; Stucker, I.; Matsuo, K.; et al. Adult Height and Head and Neck Cancer: A Pooled Analysis within the INHANCE Consortium. Head Neck 2014, 36, 1391. [Google Scholar] [CrossRef]
  5. Facompre, N.D.; Rajagopalan, P.; Sahu, V.; Pearson, A.T.; Montone, K.T.; James, C.D.; Gleber-Netto, F.O.; Weinstein, G.S.; Jalaly, J.; Lin, A.; et al. Identifying Predictors of HPV-related Head and Neck Squamous Cell Carcinoma Progression and Survival through Patient-derived Models. Int. J. Cancer 2020, 147, 3236–3249. [Google Scholar] [CrossRef]
  6. Shinomiya, H.; Nibu, K. ichi Etiology, Diagnosis, Treatment, and Prevention of Human Papilloma Virus-Associated Oropharyngeal Squamous Cell Carcinoma. Int. J. Clin. Oncol. 2023, 28, 975–981. [Google Scholar] [CrossRef]
  7. Bates, J.E.; Steuer, C.E. HPV as a Carcinomic Driver in Head and Neck Cancer: A De-Escalated Future? Curr. Treat. Options Oncol. 2022, 23, 325–332. [Google Scholar] [CrossRef]
  8. Bragado, P.; Estrada, Y.; Sosa, M.S.; Avivar-Valderas, A.; Cannan, D.; Genden, E.; Teng, M.; Ranganathan, A.C.; Wen, H.-C.; Kapoor, A.; et al. Analysis of Marker-Defined HNSCC Subpopulations Reveals a Dynamic Regulation of Tumor Initiating Properties. PLoS ONE 2012, 7, e29974. [Google Scholar] [CrossRef]
  9. Eder, T.; Hess, A.K.; Konschak, R.; Stromberger, C.; Jöhrens, K.; Fleischer, V.; Hummel, M.; Balermpas, P.; von der Grün, J.; Linge, A.; et al. Interference of Tumour Mutational Burden with Outcome of Patients with Head and Neck Cancer Treated with Definitive Chemoradiation: A Multicentre Retrospective Study of the German Cancer Consortium Radiation Oncology Group. Eur. J. Cancer 2019, 116, 67–76. [Google Scholar] [CrossRef]
  10. Wienecke, A.; Kraywinkel, K. Epidemiology of Head and Neck Cancer in Germany. Onkologe 2019, 25, 190–200. [Google Scholar] [CrossRef]
  11. Grünwald, V.; Chirovsky, D.; Cheung, W.Y.; Bertolini, F.; Ahn, M.J.; Yang, M.H.; Castro, G.; Berrocal, A.; Sjoquist, K.; Kuyas, H.; et al. Global Treatment Patterns and Outcomes among Patients with Recurrent and/or Metastatic Head and Neck Squamous Cell Carcinoma: Results of the GLANCE H&N Study. Oral Oncol. 2020, 102, 104526. [Google Scholar] [CrossRef]
  12. Liu, Z.; Williams, M.; Stewart, J.; Glisson, B.S.; Fuller, C.; Roy-Chowdhuri, S. Evaluation of Programmed Death Ligand 1 Expression in Cytology to Determine Eligibility for Immune Checkpoint Inhibitor Therapy in Patients with Head and Neck Squamous Cell Carcinoma. Cancer Cytopathol. 2022, 130, 110–119. [Google Scholar] [CrossRef]
  13. Emancipator, K.; Huang, L.; Aurora-Garg, D.; Bal, T.; Cohen, E.E.W.; Harrington, K.; Soulières, D.; Le Tourneau, C.; Licitra, L.; Burtness, B.; et al. Comparing Programmed Death Ligand 1 Scores for Predicting Pembrolizumab Efficacy in Head and Neck Cancer. Mod. Pathol. 2020, 34, 532–541. [Google Scholar] [CrossRef]
  14. Sanchez-Canteli, M.; Granda-Díaz, R.; del Rio-Ibisate, N.; Allonca, E.; López-Alvarez, F.; Agorreta, J.; Garmendia, I.; Montuenga, L.M.; García-Pedrero, J.M.; Rodrigo, J.P. PD-L1 Expression Correlates with Tumor-Infiltrating Lymphocytes and Better Prognosis in Patients with HPV-Negative Head and Neck Squamous Cell Carcinomas. Cancer Immunol. Immunother. 2020, 69, 2089–2100. [Google Scholar] [CrossRef]
  15. Peng, Q.H.; Wang, C.H.; Chen, H.M.; Zhang, R.X.; Pan, Z.Z.; Lu, Z.H.; Wang, G.Y.; Yue, X.; Huang, W.; Liu, R.Y. CMTM6 and PD-L1 Coexpression Is Associated with an Active Immune Microenvironment and a Favorable Prognosis in Colorectal Cancer. J. Immunother. Cancer 2021, 9, e001638. [Google Scholar] [CrossRef]
  16. Yin, B.; Ding, J.; Hu, H.; Yang, M.; Huang, B.; Dong, W.; Li, F.; Han, L. Overexpressed CMTM6 Improves Prognosis and Associated With Immune Infiltrates of Ovarian Cancer. Front. Mol. Biosci. 2022, 9, 769032. [Google Scholar] [CrossRef]
  17. Zhang, C.; Zhao, S.; Wang, X. Co-Expression of CMTM6 and PD-L1: A Novel Prognostic Indicator of Gastric Cancer. Cancer Cell Int. 2021, 21, 78. [Google Scholar] [CrossRef]
  18. Dai, M.; Lan, T.; Li, X.; Xiao, B. High Expression of CMTM6 Is a Risk Factor for Poor Prognosis of Gastrointestinal Tumors: A Meta-Analysis. Asian J. Surg. 2022, 46, 66–72. [Google Scholar] [CrossRef]
  19. Zhao, Y.; Zhang, M.; Pu, H.; Guo, S.; Zhang, S.; Wang, Y. Prognostic Implications of Pan-Cancer CMTM6 Expression and Its Relationship with the Immune Microenvironment. Front. Oncol. 2021, 10, 585961. [Google Scholar] [CrossRef]
  20. Huang, Y.; Zhu, Y.; Yang, J.; Pan, Q.; Zhao, J.; Song, M.; Yang, C.; Han, Y.; Tang, Y.; Wang, Q.; et al. CMTM6 Inhibits Tumor Growth and Reverses Chemoresistance by Preventing Ubiquitination of P21 in Hepatocellular Carcinoma. Cell Death Dis. 2022, 13, 251. [Google Scholar] [CrossRef]
  21. Mohapatra, P.; Shriwas, O.; Mohanty, S.; Ghosh, A.; Smita, S.; Kaushik, S.R.; Arya, R.; Rath, R.; Das Majumdar, S.K.; Muduly, D.K.; et al. CMTM6 Drives Cisplatin Resistance by Regulating Wnt Signaling through the ENO-1/AKT/GSK3β Axis. JCI Insight 2021, 6, e143643. [Google Scholar] [CrossRef] [PubMed]
  22. Mohapatra, P.; Mohanty, S.; Ansari, S.A.; Shriwas, O.; Ghosh, A.; Rath, R.; Das Majumdar, S.K.; Swain, R.K.; Raghav, S.K.; Dash, R. CMTM6 Attenuates Cisplatin-Induced Cell Death in OSCC by Regulating AKT/c-Myc-Driven Ribosome Biogenesis. FASEB J. 2022, 36, e22566. [Google Scholar] [CrossRef] [PubMed]
  23. Pang, X.; Wang, S.S.; Zhang, M.; Jiang, J.; Fan, H.Y.; Wu, J.S.; Wang, H.F.; Liang, X.H.; Tang, Y. ling OSCC Cell-Secreted Exosomal CMTM6 Induced M2-like Macrophages Polarization via ERK1/2 Signaling Pathway. Cancer Immunol. Immunother. 2021, 70, 1015–1029. [Google Scholar] [CrossRef] [PubMed]
  24. Jing, S.-L.; Afshari, K.; Guo, Z.-C. Inflammatory Response-Related Genes Predict Prognosis in Patients with HNSCC. Immunol. Lett. 2023, 259, 46–60. [Google Scholar] [CrossRef] [PubMed]
  25. Lin, C.; Chen, Y.; Pan, J.; Lu, Q.; Ji, P.; Lin, S.; Liu, C.; Lin, S.; Li, M.; Zong, J. Identification of an Individualized Therapy Prognostic Signature for Head and Neck Squamous Cell Carcinoma. BMC Genomics 2023, 24, 221. [Google Scholar] [CrossRef] [PubMed]
  26. Klapper, L.; Idel, C.; Kuppler, P.; Jagomast, T.; von Bernuth, A.; Bruchhage, K.L.; Rades, D.; Offermann, A.; Kirfel, J.; Perner, S.; et al. TRIM24 Expression as an Independent Biomarker for Prognosis and Tumor Recurrence in HNSCC. J. Pers. Med. 2022, 12, 991. [Google Scholar] [CrossRef] [PubMed]
  27. von Bernuth, A.; Ribbat-Idel, J.; Klapper, L.; Jagomast, T.; Rades, D.; Leichtle, A.; Pries, R.; Bruchhage, K.L.; Perner, S.; Offermann, A.; et al. TRIM21 Expression as a Prognostic Biomarker for Progression-Free Survival in HNSCC. Int. J. Mol. Sci. 2023, 24, 5140. [Google Scholar] [CrossRef]
  28. Ribbat-Idel, J.; Perner, S.; Kuppler, P.; Klapper, L.; Krupar, R.; Watermann, C.; Paulsen, F.O.; Offermann, A.; Bruchhage, K.L.; Wollenberg, B.; et al. Immunologic “Cold” Squamous Cell Carcinomas of the Head and Neck Are Associated With an Unfavorable Prognosis. Front. Med. 2021, 8, 622330. [Google Scholar] [CrossRef]
  29. Schildhaus, H.U. Predictive Value of PD-L1 Diagnostics. Pathologe 2018, 39, 498–519. [Google Scholar] [CrossRef]
  30. Burr, M.L.; Sparbier, C.E.; Chan, Y.C.; Williamson, J.C.; Woods, K.; Beavis, P.A.; Lam, E.Y.N.; Henderson, M.A.; Bell, C.C.; Stolzenburg, S.; et al. CMTM6 Maintains the Expression of PD-L1 and Regulates Anti-Tumour Immunity. Nature 2017, 549, 101–105. [Google Scholar] [CrossRef]
  31. Liang, J.; Li, S.; Li, W.; Rao, W.; Xu, S.; Meng, H.; Zhu, F.; Zhai, D.; Cui, M.; Xu, D.; et al. CMTM6, a Potential Immunotherapy Target. J. Cancer Res. Clin. Oncol. 2022, 148, 47–56. [Google Scholar] [CrossRef] [PubMed]
  32. Wang, Q.; Shao, X.; Zhang, Y.; Zhu, M.; Wang, F.X.C.; Mu, J.; Li, J.; Yao, H.; Chen, K. Role of Tumor Microenvironment in Cancer Progression and Therapeutic Strategy. Cancer Med. 2023, 12, 11149–11165. [Google Scholar] [CrossRef]
  33. Chen, L.; Yang, Q.C.; Li, Y.C.; Yang, L.L.; Liu, J.F.; Li, H.; Xiao, Y.; Bu, L.L.; Zhang, W.F.; Sun, Z.J. Targeting CMTM6 Suppresses Stem Cell-Like Properties and Enhances Antitumor Immunity in Head and Neck Squamous Cell Carcinoma. Cancer Immunol. Res. 2020, 8, 179–191. [Google Scholar] [CrossRef] [PubMed]
  34. Chandrashekar, D.S.; Bashel, B.; Balasubramanya, S.A.H.; Creighton, C.J.; Ponce-Rodriguez, I.; Chakravarthi, B.V.S.K.; Varambally, S. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia 2017, 19, 649–658. [Google Scholar] [CrossRef] [PubMed]
  35. Zugazagoitia, J.; Liu, Y.; Toki, M.; McGuire, J.; Ahmed, F.S.; Henick, B.S.; Gupta, R.; Gettinger, S.N.; Herbst, R.S.; Schalper, K.A.; et al. Quantitative Assessment of CMTM6 in the Tumor Microenvironment and Association with Response to PD-1 Pathway Blockade in Advanced-Stage Non-Small Cell Lung Cancer. J. Thorac. Oncol. 2019, 14, 2084–2096. [Google Scholar] [CrossRef] [PubMed]
  36. Burtness, B.; Harrington, K.J.; Greil, R.; Soulières, D.; Tahara, M.; de Castro, G.; Psyrri, A.; Basté, N.; Neupane, P.; Bratland, Å.; 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]
  37. Saba, N.F.; Blumenschein, G.; Guigay, J.; Licitra, L.; Fayette, J.; Harrington, K.J.; Kiyota, N.; Gillison, M.L.; Ferris, R.L.; Jayaprakash, V.; et al. Nivolumab versus Investigator’s Choice in Patients with Recurrent or Metastatic Squamous Cell Carcinoma of the Head and Neck: Efficacy and Safety in CheckMate 141 by Age. Oral Oncol. 2019, 96, 7–14. [Google Scholar] [CrossRef]
  38. 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]
  39. de Ruiter, E.J.; Ooft, M.L.; Devriese, L.A.; Willems, S.M. The Prognostic Role of Tumor Infiltrating T-Lymphocytes in Squamous Cell Carcinoma of the Head and Neck: A Systematic Review and Meta-Analysis. Oncoimmunology 2017, 6, e1356148. [Google Scholar] [CrossRef]
  40. Kim, T.; Jung, S.H.; Kim, S.K.; Kwon, H.J. P16 Expression and Its Association with PD-L1 Expression and FOXP3-Positive Tumor Infiltrating Lymphocytes in Head and Neck Squamous Cell Carcinoma. Mol. Cell. Toxicol. 2019, 15, 137–143. [Google Scholar] [CrossRef]
  41. Seminerio, I.; Descamps, G.; Dupont, S.; de Marrez, L.; Laigle, J.A.; Lechien, J.R.; Kindt, N.; Journe, F.; Saussez, S. Infiltration of FoxP3+ Regulatory T Cells Is a Strong and Independent Prognostic Factor in Head and Neck Squamous Cell Carcinoma. Cancers 2019, 11, 227. [Google Scholar] [CrossRef]
  42. Yang, Y.; He, W.; Wang, Z.R.; Wang, Y.J.; Li, L.L.; Lu, J.Z.; Tao, Y.; Zhang, J.; Fu, S.J.; Wang, Z.P.; et al. Immune Cell Landscape in Gastric Cancer. Biomed Res. Int. 2021, 2021, 1930706. [Google Scholar] [CrossRef] [PubMed]
  43. Hoffmann, F.; Franzen, A.; de Vos, L.; Wuest, L.; Kulcsár, Z.; Fietz, S.; Maas, A.P.; Hollick, S.; Diop, M.Y.; Gabrielpillai, J.; et al. CTLA4 DNA Methylation Is Associated with CTLA-4 Expression and Predicts Response to Immunotherapy in Head and Neck Squamous Cell Carcinoma. Clin. Epigenetics 2023, 15, 112. [Google Scholar] [CrossRef] [PubMed]
  44. Sholl, L.M. Biomarkers of Response to Checkpoint Inhibitors beyond PD-L1 in Lung Cancer. Mod. Pathol. 2022, 35, 66–74. [Google Scholar] [CrossRef] [PubMed]
  45. Ayers, M.; Lunceford, J.; Nebozhyn, M.; Murphy, E.; Loboda, A.; Kaufman, D.R.; Albright, A.; Cheng, J.D.; Kang, S.P.; Shankaran, V.; et al. IFN-γ-Related MRNA Profile Predicts Clinical Response to PD-1 Blockade. J. Clin. Investig. 2017, 127, 2930–2940. [Google Scholar] [CrossRef]
  46. Wang, H.; Liu, B.; Wei, J. Beta2-Microglobulin(B2M) in Cancer Immunotherapies: Biological Function, Resistance and Remedy. Cancer Lett. 2021, 517, 96–104. [Google Scholar] [CrossRef]
  47. Amodio, V.; Mauri, G.; Reilly, N.M.; Sartore-bianchi, A.; Siena, S.; Bardelli, A.; Germano, G. Mechanisms of Immune Escape and Resistance to Checkpoint Inhibitor Therapies in Mismatch Repair Deficient Metastatic Colorectal Cancers. Cancers 2021, 13, 2638. [Google Scholar] [CrossRef]
  48. Pontes, H.A.R.; Pontes, F.S.C.; de Jesus, A.S.; Soares, M.C.P.; Gonçalves, F.L.N.; de Lucena Botelho, T.; do Carmo Ribeiro, J.; dos Santos Pinto, D. P-Akt and Its Relationship with Clinicopathological Features and Survival in Oral Squamous Cell Carcinoma: An Immunohistochemical Study. J. Oral Pathol. Med. 2015, 44, 532–537. [Google Scholar] [CrossRef]
  49. Archewa, P.; Pata, S.; Chotjumlong, P.; Supanchart, C.; Krisanaprakornkit, S.; Iamaroon, A. Akt2 and P-Akt Overexpression in Oral Cancer Cells Is Due to a Reduced Rate of Protein Degradation. J. Investig. Clin. Dent. 2017, 8, e12194. [Google Scholar] [CrossRef]
  50. Freudlsperger, C.; Horn, D.; Weißfuß, S.; Weichert, W.; Weber, K.J.; Saure, D.; Sharma, S.; Dyckhoff, G.; Grabe, N.; Plinkert, P.; et al. Phosphorylation of AKT(Ser473) Serves as an Independent Prognostic Marker for Radiosensitivity in Advanced Head and Neck Squamous Cell Carcinoma. Int. J. Cancer 2015, 136, 2775–2785. [Google Scholar] [CrossRef]
  51. Janecka-Widła, A.; Majchrzyk, K.; Mucha-Małecka, A.; Słonina, D.; Biesaga, B. Prognostic Potential of Akt, PAkt(Ser473) and PAkt(Thr308) Immunoreactivity in Relation to HPV Prevalence in Head and Neck Squamous Cell Carcinoma Patients. Pathol. Res. Pract. 2022, 229, 153684. [Google Scholar] [CrossRef]
  52. Martinez-Morilla, S.; Moutafi, M.; Rimm, D.L. Standardization of PD-L1 Immunohistochemistry. Mod. Pathol. 2021, 35, 294–295. [Google Scholar] [CrossRef]
Figure 1. Kaplan–Meier survival curves of the HNSCC patients. The discovery (n = 177) and validation cohort (n = 286) included HNSCC patients with comparable clinicopathological parameters receiving first- or second-line treatment or best supportive care. (A) Progression-free survival (PFS) and overall survival (OS) for the patients from both cohorts. Treatment regimens included radiotherapy (RT), radio-chemotherapy (RCT), chemotherapy, immune-oncology (IO), or ICI (immune checkpoint inhibition). (B) Treatment-related OS in the discovery cohort (n = 177) stratified according to the specific regimen. (left) Cisplatin-based chemotherapy (n = 64), RT (n = 29) or RCT (n = 13). (right) ICI mono (n = 14 and 4, respectively) or ICI-based combination regimens (n = 8). (C) Treatment-related OS in the validation cohort (n = 286). Log-rank analysis was performed to study differences between the individual regimens. n.s.—not significant.
Figure 1. Kaplan–Meier survival curves of the HNSCC patients. The discovery (n = 177) and validation cohort (n = 286) included HNSCC patients with comparable clinicopathological parameters receiving first- or second-line treatment or best supportive care. (A) Progression-free survival (PFS) and overall survival (OS) for the patients from both cohorts. Treatment regimens included radiotherapy (RT), radio-chemotherapy (RCT), chemotherapy, immune-oncology (IO), or ICI (immune checkpoint inhibition). (B) Treatment-related OS in the discovery cohort (n = 177) stratified according to the specific regimen. (left) Cisplatin-based chemotherapy (n = 64), RT (n = 29) or RCT (n = 13). (right) ICI mono (n = 14 and 4, respectively) or ICI-based combination regimens (n = 8). (C) Treatment-related OS in the validation cohort (n = 286). Log-rank analysis was performed to study differences between the individual regimens. n.s.—not significant.
Cancers 15 05597 g001
Figure 2. Kaplan–Meier survival curves of HNSCC patients. For the biomarker-driven and treatment-related prognosis, analysis was performed on n = 177 patients (discovery cohort) receiving first- or second-line treatment. Treatment regimens included chemotherapy, radiotherapy (RT), radio-chemotherapy (RCT), or ICI (immune checkpoint inhibition). (A,B) Overall survival (OS) according to PD-L1 or CMTM6 in all patients (left) and stratified according to HPV status (middle and right). (C) Treatment-related OS stratified according to the specific regimen and PD-L1 status. (D) Treatment-related OS stratified according to the specific regimen and CMTM6 status. Numbers of patients for each treatment group, i.e., cisplatin, RT, RCT, or ICI are given in their respective graphs. Log-rank analysis was performed to study any differences between the individual regimens. n.s.—not significant.
Figure 2. Kaplan–Meier survival curves of HNSCC patients. For the biomarker-driven and treatment-related prognosis, analysis was performed on n = 177 patients (discovery cohort) receiving first- or second-line treatment. Treatment regimens included chemotherapy, radiotherapy (RT), radio-chemotherapy (RCT), or ICI (immune checkpoint inhibition). (A,B) Overall survival (OS) according to PD-L1 or CMTM6 in all patients (left) and stratified according to HPV status (middle and right). (C) Treatment-related OS stratified according to the specific regimen and PD-L1 status. (D) Treatment-related OS stratified according to the specific regimen and CMTM6 status. Numbers of patients for each treatment group, i.e., cisplatin, RT, RCT, or ICI are given in their respective graphs. Log-rank analysis was performed to study any differences between the individual regimens. n.s.—not significant.
Cancers 15 05597 g002
Figure 3. Kaplan–Meier survival curves of HNSCC patients from the discovery cohort. For biomarker-driven and treatment-related prognosis, analysis was performed on n = 177 patients receiving first- or second-line treatment. Treatment regimens included chemotherapy, radiotherapy (RT), radio-chemotherapy (RCT), or ICI (immune checkpoint inhibition). (A) Overall survival (OS) of patients according to FOXP3 and combined FOXP3/CMTM6 status. (B) Treatment-related OS stratified according to the specific regimen and FOXP3/CMTM6 status. Cisplatin-based chemotherapy, RT, RCT, or ICI. Log-rank analysis was performed to study any differences between the individual regimens. n.s.—not significant.
Figure 3. Kaplan–Meier survival curves of HNSCC patients from the discovery cohort. For biomarker-driven and treatment-related prognosis, analysis was performed on n = 177 patients receiving first- or second-line treatment. Treatment regimens included chemotherapy, radiotherapy (RT), radio-chemotherapy (RCT), or ICI (immune checkpoint inhibition). (A) Overall survival (OS) of patients according to FOXP3 and combined FOXP3/CMTM6 status. (B) Treatment-related OS stratified according to the specific regimen and FOXP3/CMTM6 status. Cisplatin-based chemotherapy, RT, RCT, or ICI. Log-rank analysis was performed to study any differences between the individual regimens. n.s.—not significant.
Cancers 15 05597 g003
Figure 4. Gene signatures in HNSCC patients from the discovery group with their response to immune checkpoint inhibitors. (A) Heatmap representing the down- and upregulated genes involved in antitumor immune activity by unsupervised hierarchical clustering. Scores are scaled to have a mean of zero and a standard deviation of one. The standardized signature scores are truncated at ±3 standard deviations (99% of the data should fall within ±3 standard deviations of the mean). (B) Box plot indicating a significantly lower CTLA-4 signature in patients with a response to ICI as assessed by univariate analysis (p = 0.02). (C) Volcano plot of all gene signatures displays no significant difference between the ICI responders and non-responders in a multivariate analysis. Up indicates upregulated and down indicates downregulated genes. (D) Survival volcano plot of all gene signatures displays no significant difference between the signature’s hazard ratios in a multivariate analysis. Signatures further to the right are associated with a decreased risk of an event relative to the baseline. (E) Box plot indicating a significantly higher IFN downstream signature in patients whose tumors showed a CPS ≥ 1 for PD-L1 immunohistochemistry as assessed by univariate analysis (p = 0.04). (F) Spaghetti plot displaying a decrease in antitumor immune activity from initial diagnosis to r/m disease with start of ICI until the time point of progression under ICI in non-responders, and (G) an increase in CTLA-4 in recurrent tumor material from a patient who initially had a partial response to ICI. * p < 0.05. CPS—combined positive score; Pre-ICI—sample taken just before the start of ICI; Post-ICI—sample taken just after progression under ICI.
Figure 4. Gene signatures in HNSCC patients from the discovery group with their response to immune checkpoint inhibitors. (A) Heatmap representing the down- and upregulated genes involved in antitumor immune activity by unsupervised hierarchical clustering. Scores are scaled to have a mean of zero and a standard deviation of one. The standardized signature scores are truncated at ±3 standard deviations (99% of the data should fall within ±3 standard deviations of the mean). (B) Box plot indicating a significantly lower CTLA-4 signature in patients with a response to ICI as assessed by univariate analysis (p = 0.02). (C) Volcano plot of all gene signatures displays no significant difference between the ICI responders and non-responders in a multivariate analysis. Up indicates upregulated and down indicates downregulated genes. (D) Survival volcano plot of all gene signatures displays no significant difference between the signature’s hazard ratios in a multivariate analysis. Signatures further to the right are associated with a decreased risk of an event relative to the baseline. (E) Box plot indicating a significantly higher IFN downstream signature in patients whose tumors showed a CPS ≥ 1 for PD-L1 immunohistochemistry as assessed by univariate analysis (p = 0.04). (F) Spaghetti plot displaying a decrease in antitumor immune activity from initial diagnosis to r/m disease with start of ICI until the time point of progression under ICI in non-responders, and (G) an increase in CTLA-4 in recurrent tumor material from a patient who initially had a partial response to ICI. * p < 0.05. CPS—combined positive score; Pre-ICI—sample taken just before the start of ICI; Post-ICI—sample taken just after progression under ICI.
Cancers 15 05597 g004
Table 1. Clinico-pathological characteristics of the discovery cohort.
Table 1. Clinico-pathological characteristics of the discovery cohort.
Group CharacteristicsΣ n = 177
Female n (%)
Male n (%)
34 (19.2)
143 (80.8)
Median age (years ± SD)66 ± 9.8
Noxae
smoking (>10 py in %): yes/no
alcohol (>1 drink/d in %): yes/no
no information (smoking/alcohol)

51.9/34.5
36.7/63.3
11.3/2.3
p16/HPV status (%)
positive
negative

21
79
Tumor stage (UICC stage 8th edition) (n)
1/2/3/4

37/56/29/55
Localization (n)
oral cavity
oropharynx
hypopharynx
larynx

79
65
18
15
Treatment (n)
C/RT/RCT/ICI/ICI + C
Other/no adjuvant therapy

66/33/13/27/3
16/19
Median survival (month)
PFS
OS

21.0 ± 20.6
29.0 ± 21.4
Values are given as absolute/relative numbers and mean ± SD unless otherwise stated. Abbreviations: py—pack years; d—day; HPV—human papilloma virus; RT—radiotherapy; RCT—radio-chemotherapy; C—chemotherapy (not specified); ICI—immune checkpoint inhibition (nivolumab or pembrolizumab); PFS—progression-free survival; OS—overall survival.
Table 2. Immunohistochemistry staining results for the markers in the discovery group.
Table 2. Immunohistochemistry staining results for the markers in the discovery group.
Defined Criteria for Evaluation
Cut-Off *Low/NegativeHigh/PositiveMedian
n = 177n (%)n (%)n (%)n (%)
CMTM6(CPS) 1057 (32.2)120 (67.8)33.7
PD-L1(CPS) 158 (32.8)119 (67.2)14.9
CTLA-4(density of positive cells) low vs. high93 (52.5)84 (47.5)-
FOXP3(density of positive cells) low vs. high77 (43.5)100 (56.5)-
AKT(% positive cells × staining intensity) no/low vs. high119 (67.2)58 (32.8)-
pAKT(% positive cells × staining intensity) no/low vs. high152 (85.9)25 (14.1) 67.2-
* in (brackets): scheme for evaluation. Abbreviations: CPS, combined positive score.
Table 3. Spearman correlation for markers in the discovery cohort including samples from patients receiving no adjuvant or either chemotherapy/radiation/chemoradiation.
Table 3. Spearman correlation for markers in the discovery cohort including samples from patients receiving no adjuvant or either chemotherapy/radiation/chemoradiation.
CMTM6 vs.Biomarker
PD-L1CTLA-4FOXP3AKTpAKT
Treatment
Unrelated
n = 151 #
r0.260.630.850.310.22
95% CI0.098 to 0.410.52 to 0.720.80 to 0.890.15 to 0.460.047 to 0.38
p (two-tailed)0.0014<0.0001<0.00010.00020.0110
p value summary**************
RT
n = 33
r0.320.760.890.190.16
95% CI−0.061 to 0.620.54 to 0.880.78 to 0.95−0.20 to 0.54−0.26 to 0.53
p (two-tailed)0.0877<0.0001<0.00010.32210.4350
p value summaryns********nsns
Cisplatin
n = 66
r0.360.550.850.030.13
95% CI−0.12 to 0.560.34 to 0.700.76 to 0.91−0.25 to 0.31−0.16 to 0.40
p (two-tailed)0.0037<0.0001<0.00010.83110.3636
p value summary**********nsns
RCT
n = 13
r0.510.840.710.290.31
95% CI−0.10 to 0.850.53 to 0.950.25 to 0.91−0.39 to 0.77−0.38 to 0.77
p (two-tailed)0.0890.00060.00780.41820.3549
p value summaryns*****nsns
Interpretation: ≥0.70: very strong relationship; 0.40–0.69: strong relationship; 0.30–0.39: moderate relationship; 0.20–0.29: weak relationship; 0.01–0.19: no relationship. ns—not significant. # PD-L1 n = 148; AKT n = 134; pAKT n = 130. * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001.
Table 4. Contingency table for the markers in the discovery group.
Table 4. Contingency table for the markers in the discovery group.
Biomarker
n = 177CMTM6 (n) PD-L1 (n) FOXP3 (n) CTLA-4 (n)
LowHighpNegativePositivepLowHighpLowHighp
Gender ns ns ns ns
   Male38105 4895 6083 7568
   Female925 1024 1717 1816
Smoking # ns ns ns ns
   Yes3359 3854 3755 5636
   No2936 2045 3530 3332
Alcohol ## ns ns ns ns
   Yes3031 2833 2437 2833
   No4369 4171 4963 6547
HPV
   Pos.1028 929 1326 1524
   Neg.4792 4891 6474 6969
TNM ns ns ns *
   T1928 1225 1819 1027
   T21838 1541 2630 2531
   T31118 1316 920 1514
   T41936 1837 2431 3421
   N−/N+23/3466/54ns23/3566/53ns37/4048/52ns40/4449/44ns
   M0/M143/14102/17ns43/15103/16ns63/1483/17ns65/1981/12ns
Recurrence
   No3492*2491****4891****25100**
   Yes2328 3428 279 2131
# for n = 157; ## for n = 173; * p < 0.05; ** p < 0.01; **** p < 0.0001. Abbreviations: pp value; HPV—human papilloma virus; ns—not significant.
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

Becker, A.-S.; Kluge, C.; Schofeld, C.; Zimpfer, A.H.; Schneider, B.; Strüder, D.; Redwanz, C.; Ribbat-Idel, J.; Idel, C.; Maletzki, C. Identifying Predictive Biomarkers for Head and Neck Squamous Cell Carcinoma Response. Cancers 2023, 15, 5597. https://doi.org/10.3390/cancers15235597

AMA Style

Becker A-S, Kluge C, Schofeld C, Zimpfer AH, Schneider B, Strüder D, Redwanz C, Ribbat-Idel J, Idel C, Maletzki C. Identifying Predictive Biomarkers for Head and Neck Squamous Cell Carcinoma Response. Cancers. 2023; 15(23):5597. https://doi.org/10.3390/cancers15235597

Chicago/Turabian Style

Becker, Anne-Sophie, Cornelius Kluge, Carsten Schofeld, Annette Helene Zimpfer, Björn Schneider, Daniel Strüder, Caterina Redwanz, Julika Ribbat-Idel, Christian Idel, and Claudia Maletzki. 2023. "Identifying Predictive Biomarkers for Head and Neck Squamous Cell Carcinoma Response" Cancers 15, no. 23: 5597. https://doi.org/10.3390/cancers15235597

APA Style

Becker, A. -S., Kluge, C., Schofeld, C., Zimpfer, A. H., Schneider, B., Strüder, D., Redwanz, C., Ribbat-Idel, J., Idel, C., & Maletzki, C. (2023). Identifying Predictive Biomarkers for Head and Neck Squamous Cell Carcinoma Response. Cancers, 15(23), 5597. https://doi.org/10.3390/cancers15235597

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