|Home | About | Journals | Submit | Contact Us | Français|
Melanoma associated antigens (MAGE) are potential targets for immunotherapy and have been associated with poor overall survival (OS) in head and neck squamous cell carcinoma (HNSCC). However, little is known about MAGE in lymph node metastases (LNM) and recurrent disease (RD) of HNSCC.
To assess whether MAGE expression increases with metastasis or recurrence, a tissue microarray (TMA) of 552 primary tumors (PT), 219 LNM and 75 RD was evaluated by immunohistochemistry for MAGE antigens using three monoclonal antibodies to multiple MAGE family members. Mean expression intensity (MEI) was obtained from triplicates of each tumor specimen.
The median MEI compared between PT, LNM and RD was significantly higher in LNM and RD. In paired samples, MEI was comparable in PT to respective LNM, but significantly different from RD. Up to 25% of patients were negative for pan-MAGE or MAGE-A3/A4 in PT, but positive in RD. The prognostic impact of MAGE expression was validated in the TMA cohort and also in TCGA data (mRNA). OS was significantly lower for patients expressing pan-MAGE or MAGE-A3/A4 in both independent cohorts.
MAGE expression was confirmed as a prognostic marker in HNSCC and may be important for immunotherapeutic strategies as a shared antigen.
Head and neck squamous cell carcinoma (HNSCC) of the oral cavity, pharynx and larynx combined account for 60,000 new cases of malignant tumors and 12,000 tumor related deaths annually in the USA, with comparable numbers in Europe [1–3]. The prognosis of HNSCC patients remains unsatisfactory with 5-year overall survival rates of approximately 50–65% without significant improvement in the last decades . Recent advances in the field of cancer immunology and cancer immunotherapy imply that immunotherapy may be key to improve this outcome. Notably, a subset of HNSCC patients treated with checkpoint blockade therapy experience clinical benefit [5–7]. Characterizing tumor antigens recognized by the immune system and potentially targeted during immunotherapy is therefore likely to yield important biomarkers of prognostic or predictive value. Cancer-testis antigens (CTA) are a group of immunogenic tumor antigens, which can be specifically detected in different cancers, but not in healthy normal tissue except for human testis [8, 9]. Their restricted expression and immunogenicity make them good candidates to assess as shared targets of immunotherapy.
Melanoma-associated antigens (MAGE) are among the most frequently expressed CTA in HNSCC [10–13]. We previously demonstrated in a large cohort of surgically treated HNSCC patients that expression of pan-MAGE and MAGE-A3/A4 is associated with reduced overall survival (OS). The prognostic value was strongest if the staining pattern was considered. Simultaneous cytoplasmic (cyt) and nuclear (nuc) antigen expression (cyt+nuc) was identified as an independent prognostic marker .
Because MAGE expression is often heterogeneous within primary tumors (PT), dynamic changes and clonal evolution of MAGE expression should be also assessed upon tumor progression. Yet, very little is known about MAGE expression in lymph node metastases (LNM) or recurrent disease (RD) compared to PT. Since the most frequent types of treatment failure are local and regional recurrences , it is important to correlate expression of MAGE in primary tumors, lymph node metastases, and recurrent tumors within the same patients.
Mechanistic data from bladder cancer imply that MAGE-A3 is predominantly expressed in a cancer stem cell like subpopulation . Further mechanistic data from HNSCC cell lines showed reduced treatment efficacy in MAGE expressing cell lines treated with different chemotherapeutics . These factors would suggest an enrichment of MAGE expressing cancer cells during cancer evolution, especially after unsuccessful cancer treatment. At the same time, immune responses to MAGE have been described in HNSCC and may result in reduced MAGE expression during tumor evolution . Based on these previous mechanistic studies, we hypothesized that MAGE expression may underlie a selection pressure during tumor evolution resulting in differential expression in primary tumors, respective lymph node metastases and recurrences. We took advantage of a unique large set of patients with paired tissues collected from primaries, metastases, and recurrences, and used them to assess MAGE evolution and to validate the prognostic impact of MAGE expression in independent HNSCC cohorts.
The cohort comprised 129 oral cavity, 175 oropharyngeal, 66 hypopharyngeal and 182 laryngeal primaries. Mean OS was 105.4 months, mean recurrence-free survival (RFS) was 82.2 months. The majority of patients were treated surgically (n = 458), 74 patients were treated by primary chemoradiation and 20 patients were treated with palliative intention. Clinico-pathological criteria of the patient cohort separately for each primary site can be found in Table Table11.
As previously reported, we observed three different expression patterns of MAGE for pan-MAGE (M3H67) and MAGE-A3/A4 (57B) expression: cyt, nuc and cyt+nuc. MAGE-A1 (MA454) is always solely expressed cytoplasmically . Representative staining examples of negative samples and the different expression patterns (cyt, nuc, cyt+nuc) are given in Figure 1A–1D. Additionally, examples for different mean expression intensities (MEI) are presented in Figure 1E–1H.
Visual scoring of antigen expression and expression pattern was performed microscopically by two independent analysts. The expression was considered positive if any of the three triplicates for each sample stained positive for the respective antigen. Pan-MAGE was found in 173/ 423 (41%) of primary tumors, 97/212 (46%) of lymph node metastases and 18/49 (37%) of recurrences respectively. MAGE-A3/A4 was positive in 166/423 (39%) of primary tumors, 90/212 (42%) of lymph node metastases and 18/49 (37%) of recurrences. MAGE-A1 was less frequently found with 78/424 (18%) in primary tumors, 38/212 (18%) in lymph node metastases and 12/49 (25%) in recurrent tumors. Notably, the frequency of combined cyt+nuc expression seems to be higher in RD. Detailed numbers and the distribution of the expression pattern for each antigen are presented in Table Table22.
MAGE expression frequency in primary tumors was significantly different based on the primary site for pan-MAGE (p = 0.04), MAGE-A3/A4 (p = 0.04) and MAGE-A1 (p = 0.002) whereas the lowest expression rate was found in oral cavity primaries and the highest expression rate in hypopharyngeal primaries. MAGE expression was less frequently found in HPV positive cases based on HPV-DNA for pan-MAGE (HPV-DNAnegative/pan-MAGEpositive: 165/387 (42.6%), HPV-DNApositive/pan-MAGEpositive: 8/36 (22.2%) p = 0.018) as well as for MAGE-A3/A4 (HPV-DNAnegative/MAGE-A3/A4positive: 159/387 (41.0%), HPV-DNApositive/MAGE-A3/A4positive: 6/36 (16.7%), p = 0.013). Based on HPV DNA and p16 status similar results were found for pan-MAGE (HPV(DNA+p16)negative/pan-MAGEpositive 168/393 (42.7%), HPV(DNA+p16)positive/pan-MAGEpositive 2/26 (7.7%), p = 0.002) and MAGE-A3/A4 (HPV(DNA+p16)negative/ MAGE-A3/A4positive 161/393 (40.9%), HPV(DNA+p16)positive/pan-MAGEpositive 3/26 (11.5%), p = 0.011). Detailed numbers including the distribution of the expression pattern for MAGE positive cases are provided in Table Table22.
We first analyzed data in each type of samples (PT, LNM, RD), focusing only on cases with MEI ≥ 0.1, and found that the median MEI of the groups differed. Pan-MAGE: PT median MEI = 0.24 (n = 120), LNM median MEI = 0.36 (n = 68), RD median MEI = 0.37 (n = 14). MAGE-A3/A4: primary PT median MEI = 0.27 (n = 101), LNM median MEI = 0.31 (n = 67), RD median MEI = 0.40 (n = 12). MAGE-A1: PT median MEI=0.21 (n = 26), LNM median MEI = 0.26 (n = 20), RD median MEI = 0.33 (n = 7).
The MEI medians of the respective groups (primary tumors, lymph node metastases, recurrences) were compared by Kruskal-Wallis test. For pan-MAGE (Kruskal-Wallis: 8.06, p = 0.018) and MAGE-A1 (Kruskal-Wallis: 6.40, p = 0.041), a significant difference between the median MEI was found. The differences between the median MEI for MAGE-A3/A4 expression did not reach significance (Kruskal-Wallis: 4.39, p = 0.111). MEI distribution and the median MEI with interquartile range for the three groups (PT, LNM, RD) are depicted for each antigen in Figure Figure22.
In order to compare the MEI between primary tumors (PT), lymph node metastases (LNM) and recurrent disease (RD) in individual patients, paired samples were analyzed. Paired samples were compared using Wilcoxon signed rank test for related samples.
The median of differences between the MEI of PT and respective LNM was not different for pan-MAGE (n = 166, p = 0.900), MAGE-A3/A4 (n = 153, p = 0.594) or MAGE-A1 (n = 150, 0.601). Individual MEI values for each paired sample are depicted in Figure 3A–3C. The graph is divided into four quadrants by a cut-off of MEI = 0.1. The lower left quadrants (Q1: PT and LNM negative (MEI < 0.1)) and upper right quadrants (Q3: PT positive and LNM positive (MEI ≥ 0.1)) show concordant samples. The upper left quadrants (Q2: PT negative (MEI < 0.1) and LNM positive (MEI ≥ 0.1)) and lower right quadrants (Q4: PT positive (MEI ≥ 0.1) and LNM negative (MEI < 0.1)) show discordant samples. Most related sample MEI values are concordant: pan-MAGE: 129/166 (77.7%); MAGE-A3/A4: 119/153 (77.8%); MAGE-A1: 136/150 (90.7%)). For pan-MAGE, 37/166 (22.3%) cases showed discordant expression between PT and LNM (Q2: 15/166 (9%); Q4: 22/166 (13.3%)). For MAGE-A3/A4, 34/153 (22.2%) (Q2: 19/153 (12.4%); Q4: 15/153 (9.8%)) and for MAGE-A1 14/150 (9.3%) (Q2: 9/150 (6%); Q4: 5/150 (3.3%) cases respectively were discordant between PT and LNM.
Comparing the median of differences between the MEI of PT and respective RD significant differences were observed for pan-MAGE (n = 23, p = 0.034), MAGE-A3/A4 (n = 20, p = 0.041), but not for MAGE-A1 (n = 18, p = 0.144). Individual MEI values for each paired sample are depicted in Figure 3D–3F. The graph is also divided into four quadrants (by a cut-off of MEI = 0.1 for PT and RD respectively. This graph shows a higher number of RD with MEI ≥ 0.1 that showed MEI < 0.1 in the respective PT (upper left quadrant: pan-MAGE: 4/23 (17.4%); MAGE-A3/A4: 5/20 (25%); MAGE-A1: 3/18 (16.7%)) than the other way around (lower right quadrant: pan-MAGE: 1/23 (4.3%); MAGE-A3/A4: 1/20 (5%); MAGE-A1: 0/18 (0%)).
Known prognostic markers were evaluated for their impact on OS by Cox regression analysis. The analysis was performed for the whole cohort with the following known prognostic parameters: T-, N- and M-classification, grading, primary treatment, human papilloma virus (HPV) status, sex, primary site and smoking habits. Results are presented in Table Table3.3. Cox regression analysis revealed the highest HR for the primary treatment (surgical vs. non-surgical; HR = 3.012 (95% CI 2.189–4.146, p < 0.001). This can be attributed to the treatment algorithms in Germany. Most centers and the majority of patients prefer primary surgical treatment in the curative setting for eligible patients. Only a minority of resectable patients receive primary (chemo-)radiotherapy. Thus, elderly patients, patients with high comorbidities or unresectable cancers are overrepresented in the non-surgical cohort. To eliminate this confounder, the following known prognostic parameters were separately evaluated in the cohort of patients who received primary surgical treatment with curative intent: T- and N-classification, resection margin status, lymphangiosis, vascular invasion, grading, extracapsular extension (ECE), HPV status, sex, primary site and smoking history. Results are presented in Table Table44.
A Pearson correlation revealed significant correlations between pan-MAGE expression and T-classification, N-status, HPV-status and smoking, but all correlation coefficients were R < 0.2. For MAGE-A3/A4 a significant correlation was found with T-classification, HPV-status and smoking. All correlation coefficients were R < 0.2. MAGE-A1 was significantly correlated with T-classification, N-status and HPV-status. All correlation coefficients were R < 0.2. All correlation coefficients and significance values are provided in Table Table55.
Due to the very low correlation coefficients, the significant correlations with known prognostic markers can be considered irrelevant.
In order to validate our previously reported prognostic impact of pan-MAGE and MAGE-A3/A4 expression on OS in surgically treated patients , we performed survival analyses for OS using the Kaplan-Meier method. Patients with available survival information were filtered for patients who received primary surgical treatment with curative intent. OS was significantly lower for patients showing expression of pan-MAGE (mean OS pan-MAGEneg (n = 202) = 122.8 months vs. mean OS pan-MAGEpos (n = 143) = 88.7 months, p = 0.025) and MAGE-A3/A4 (mean OS MAGE-A3/A4neg (n = 208) = 125.5 months vs. mean OS MAGE-A3/A4pos (n = 138) = 66.9 months, p = 0.012) (Figure 4A, 4C).
Next, we evaluated the impact of the respective expression pattern on OS. Compared to negative cases, OS for pan-MAGEcytORnuc (n = 73, mean OS=71.3 months) was lower, but did not reach significance (p = 0.173). In contrast, OS in pan-MAGEcyt+nuc (n = 70, mean OS=82.9 months) was significantly reduced (p = 0.034). Compared to negative cases, OS for MAGE-A3/A4cytORnuc (n = 67, mean OS = 68.6 months) was lower, but did not reach significance (p = 0.111). In contrast, OS in MAGE-A3/A4cyt+nuc (n = 71, mean OS=65.9 months) was significantly lower (p = 0.017) (Figure 4B, 4D).
In the subgroup of oral cavity primaries, mean OS was significantly lower for pan-MAGEcyt+nuc (31.5 months, p = 0.004) as well as pan-MAGEcytORnuc (29.7 months, p = 0.049) compared to negative cases (72.7 months). Mean OS was also significantly lower for MAGE-A3/A4cyt+nuc (28.6 months, p = 0.009) compared to negative cases (70.8 months). Mean OS of MAGE-A3/A4cytORnuc was also lower (48.2 months), but did not reach significance compared to negative cases (p = 0.246).
In the subgroup of oropharyngeal primaries, mean OS was not significantly different for pan-MAGEcyt+nuc (93.9 months, p = 0.125) or pan-MAGEcytORnuc (64.2 months, p = 0.281) compared to negative cases (87.0 months). However, mean OS was lower for MAGE-A3/A4cyt+nuc (55.7 months, p = 0.08) compared to negative cases (125.9 months). Mean OS of MAGE-A3/A4cytORnuc was significantly lower (58.7 months, p = 0.043) compared to negative cases.
In the subgroups of hypopharyngeal and laryngeal primaries, the survival was not significantly lower for pan-MAGE or MAGE-A3/A4 expression (p > 0.05).
None of the evaluated antigens was significantly associated with disease-free survival (DFS; data not shown).
A multivariate Cox regression analysis for OS taking into account T-classification, HPV-status (DNA), primary tumor site, pan-MAGE expression and MAGE-A3/A4 expression was performed. For this analysis 344 patients were available. Multivariate cox regression revealed the following markers as independently prognostic: T-classification (T2 HR = 2.2; T3 HR = 3.02; T4 HR = 3.3), HPV-status (HPV negative: HR = 3.5), primary site compared to laryngeal primaries (oral cavity: HR = 3.2; oropharynx: HR = 1.7; hypopharynx: HR = 1.1) and MAGE-A3/A4 expression (HR = 1.4).
As control cohort, TCGA data were analyzed for OS. For the gene products recognized by the M3H67 antibody (pan-MAGE: MAGE-A1, -A3, -A4, -A8, -A10, -B2, -C2), 117/522 (22 %) of cases had upregulated mRNA expression for any of the genes. OS was significantly reduced for cases with increased expression of pan-MAGE gene mRNA (negative cases: 64.78 months vs. positive cases: 30.91 months; p = 0.0057). For the gene products recognized by the 57B antibody (MAGE-A3/A4: MAGE-A3, -A4, -A6, -A12), 93/522 (18 %) of the cases showed upregulated mRNA expression for any of the genes. OS was also significantly reduced for cases expressing MAGE-A3/A4 genes (negative cases: 57.88 months vs. positive cases: 30.91 months; p = 0.0094). Survival curves for the queried genes from the TCGA dataset are given in Supplementary Figure 1. DFS was not significantly different for any of the queried gene sets in the TCGA cohort (pan-MAGE: p = 0.679; MAGE-A3/A4: p = 0.361; data not shown).
This is to the best of our knowledge, the first study examining the expression of MAGE antigens in a large cohort of HNSCC patients with focus on the comparison of expression frequency and expression intensity in paired samples of primary tumors, lymph node metastases, and recurrences.
The MAGE expression frequency of pan-MAGE and MAGE-A1 in PT in this analysis was comparable to our previously published results in a similarly large independent cohort of HNSCC patients from Hamburg, Germany, but MAGE-A3/A4 expression was found more frequently (TMA Hamburg: 27.7% vs. TMA Bonn: 39%) .
In this analysis, a significant correlation with known prognostic markers was found. However, the correlation coefficients were weak (R < 0.2). Therefore we consider this correlation irrelevant. Previously, two other studies have reported no significant correlation with clinicopathological criteria [10, 13]. In contrast to our results, one group reported a higher prevalence of MAGE-A3/A4 expression (57B) in T3/T4 laryngeal cancer and of MAGE-A1 (MA454) and pan-MAGE (M3H67) in cases which developed regional recurrence, but the corresponding correlation coefficients were not reported .
The main goal of this study was the comparison of MAGE expression in PT, LNM and RD. The use of digital pathology enabled us to obtain continuous values for the expression intensity from 0 to 1. This way we were able to obtain a mean expression intensity (MEI) for the available triplicates of each tumor sample and to compare the median of the MEI between the three groups. Significant differences between the MEI of PT, LNM and RD were observed for pan-MAGE and MAGE-A1. The graphed data (Figure (Figure2)2) indicate an increase of the MEI when comparing PT with LNM or RD, respectively. This may indicate an increasing expression intensity of MAGE during tumor evolution although the expression frequencies were similar between PT, LNM and RD (Table (Table2).2). Mechanistic data in HNSCC cell lines implied reduced responses to different chemotherapeutic drugs . MAGE-A has been shown to impair p53 function  which may lead to reduced chemosensitivity as well as reduced radiosensitivity . Another group showed increased MAGE-A3 RNA expression in a cancer-stem cell population in bladder cancer . Zamuner et al. reported reduced incidence of recurrence in MAGE-A3/A6 positive tumors . However, Zamuner et al. did not assess MAGE expression in samples of recurrent tumors. Thus, these data do not contradict our results that in recurrent tumors a higher median MEI of pan-MAGE and MAGE-A1 was found.
The comparison of the MEI in paired samples from primary tumors and recurrences revealed a significantly different MEI distribution. A higher number of recurrences with a MEI ≥ 0.1 and a MEI < 0.1 in the respective PT was observed. Published data from gastrointestinal stromal tumors identified CTA expression as a marker for early recurrence , which was validated prospectively in a subsequent study . These facts further support our hypothesis that MAGE expression intensity increases in recurrent cancer. Most published data indicate a selection advantage of MAGE expression for cancer cells and support our hypothesis of differential MAGE expression during tumor evolution. In clinical trials with MAGE-specific vaccination for RD, the expression should therefore be assessed in a sample from the recurrent tumor and not from archived tissue.
The analysis of paired samples of PT and LNM enables us to draw further conclusions. Although the MEI in paired samples of PT and LNM was well concordant, mismatches were also observed in ≈20% for pan-MAGE and MAGE-A3/A4.
Regarding our previously reported prognostic implications of pan-MAGE and MAGE-A3/A4, we were able to validate the lower OS of pan-MAGE and MAGE-A3/A4 expressing patients as well as the impaired survival in cases with cyt+nuc expression in this independent cohort . However, concerning the prognostic value of pan-MAGE and MAGE-A3/A4, the impact on prognosis in this cohort was smaller than in our previously analyzed cohort. This may be a result of the different distribution of primary sites and patient characteristics in the two cohorts. Compared to the TMA Hamburg, this cohort consisted of a smaller proportion of oral cavity cancers and a larger number of oropharyngeal tumors . One of the main prognostic factors in HNSCC is HPV. However, in the TMA Bonn, the rate of HPV DNA positive cases was low with 39/552 (7%) cases. There was a lower frequency of MAGE expression in HPV positive cases compared to negative cases. But in this cohort, the overall number of HPV DNA positive cases (n = 39) is too low to generalize these findings regarding MAGE expression in HPV positive HNSCC. The prognostic impact of MAGE expression on OS was pronounced in the oral cavity cohort in both independent patient cohorts and oral cavity squamous cell carcinoma is rarely associated with HPV infection (TMA Bonn: 4.7%). Due to the low rate of HPV positive cases and results of the multivariate Cox regression analysis, we can rule out that the prognostic disadvantage of pan-MAGE and MAGE-A3/A4 expression was actually based on the survival advantage of HPV positive patients. On the contrary, a subgroup analysis of HPV positive patients revealed that there was a non-significant trend to improved survival of MAGE-positive cases in the HPV positive subgroup (no events in the MAGE-positive patients compared to a death rate of 20% in MAGE negative cases after 5 years, Supplementary Figure 2). This may be due to increased immune-surveillance in HPV positive HNSCC [23–25].
Furthermore the laryngeal cohort from Bonn had a higher rate of early glottic laryngeal primaries known to have a good prognosis and few events to analyze. Thus, in the laryngeal subgroup of patients, we did not find reduced survival for neither pan-MAGE nor MAGE-A3/A4 expression. These facts may explain why we were not able to confirm pan-MAGE expression as an independent prognostic marker by multivariate cox regression. However, MAGE-A3/A4 expression was confirmed as an independent prognostic marker in this cohort.
The analysis of the independent cohort from the TCGA databank further confirmed the prognostic significance of pan-MAGE and MAGE-A3/A4 genes using RNA sequencing data. Thus, the prognostic implications are based on >1000 individual patients and two different methods (RNA sequencing, IHC).
Another group reported better DFS for MAGE-A3/A6 positive compared to negative patients based on polymerase chain reaction data for MAGE-A3/A6 RNA . MAGE-A6 is a homologue protein to MAGE-A3 and the 57B antibody binds to MAGE-A6 as well, as we have previously shown by ELISA . In our data no significant differences for DFS according to MAGE expression were evident. However, there was a tendency to reduced DFS for MAGE-A3/A4 positive cases, which did not reach significance. The TCGA data also support our finding that DFS is not significantly different based on the expression of MAGE. Thus, the impact of MAGE expression on the incidence of recurrence and DFS will have to be finally validated prospectively. One reason for the discrepancies between these different cohorts may lie in the extent of CTA specific immunity and whether/how these immune reactions were counterbalanced by immune modulatory checkpoints. However, a prospective validation of the prognostic impact is needed to establish MAGE expression as a prognostic marker for clinical routine.
Antigens of the MAGE family are interesting shared targets for vaccination due to their tumor specificity and immunogenicity . For HNSCC, no large trials with a relevant patient number have been performed and published to date. However, in a large phase III MAGE-A3 vaccination trial for NSCLC, no improvement of DFS was found . Central tolerance to these self antigens is discussed as an explanation for the observation of low-avidity T cells to CTA [27, 28] .
At the same time, new opportunities for immunotherapy have emerged with the successful development of different immune checkpoint modulators targeting inhibitory molecules such as cytotoxic T lymphocyte-associated antigen 4 (CTLA-4), programmed death 1 (PD-1) and its ligand PD-L1 or immune stimulatory molecules such as OX40 or CD137. Data from melanoma treated with immune checkpoint modulators suggest that a high mutational load is associated with a good clinical response to immune checkpoint blockade . This is explained by a higher probability of mutations leading to the expression of so-called mutational neoantigens which can be recognized by the immune system [29–31]. In the studies demonstrating immune responses to neoantigens in patients treated with immune checkpoint blockade, other non-mutational antigens such as cancer testis antigens (CTA) were, due to the methodology of the antigen screening process, neglected thus far. However, in melanoma patients treated with ipilimumab an increased rate of NY-ESO-1 specific immunity has been associated with improved clinical benefit of treatment, especially in patients developing both NY-ESO-1 specific antibodies and specific CD8 T cells . Thus, the addition of an antigen-specific treatment in the form of a vaccine may increase the rate of responders to immune checkpoint targeted therapy. The combination of immune-checkpoint modulators and a vaccine, composed of several different CTA, may be an opportunity to improve the benefit of CTA vaccination. To evaluate whether this strategy may be a valid approach to improve OS of CTA-positive HNSCC patients, a systematic evaluation of spontaneous immunity to CTA in HNSCC patients before and during conventional curative treatment should be performed.
Due to tumor heterogeneity, the tissue sampling process for TMA generation may represent a bias. We reduced sampling bias by sampling three separate cores of each tumor sample in the TMA construction progress. Furthermore, the addition of digital pathology to microscopic evaluation enabled us to address tumor heterogeneity by obtaining MEI as a continuous semiquantitative variable for protein expression derived from the analysis of the three triplicates.
This analysis was retrospectively performed. Thus, we cannot exclude a selection bias. Selection bias was addressed by the selection of every patient with available tissue specimen who was treated for HNSCC between 1997 and 2011 at the University Hospital of Bonn.
In conclusion, the prognostic disadvantage of MAGE expression and the increasing MEI in recurrent disease imply that the intrinsic functions of MAGE genes support an aggressive malignant phenotype. It seems that during conventional therapy, the immune system fails to outbalance this aggressiveness in MAGE positive tumors, resulting in reduced OS. A systematic analysis of MAGE-specific immunity and MAGE gene function in HNSCC is needed to shed light on the causal relationship of the prognostic disadvantage. Also, the impact of MAGE expression in relation to HPV status needs to be validated in a larger cohort of HPV positive patients. Immunotherapy, incorporating immune checkpoint modulation in combination with a vaccine may improve MAGE-specific immunity and improve OS. Still, prospective trials are needed to determine the value of MAGE expression as a prognostic marker for clinical routine.
Patients were treated according to local treatment guidelines between 1997 and 2011 at the University Hospital of Bonn. A tissue microarray (TMA) was constructed from archived formalin-fixed paraffin-embedded tissue samples after clinical diagnostics were finished. Representative 0.6mm cores were assembled into TMA blocks. The TMA included 552 primary tumors, 219 lymph node metastases and 75 recurrences. Each tumor specimen was sampled in triplicates. The study was approved by the internal review board of the University Hospital of Bonn (#174/13).
For an independent control cohort, TCGA database was accessed on May 1st, 2016. At this time 522 HNSCC cases with RNA Seq V2 data were available for analysis. The cBioPortal for Cancer Genomics (cbioportal.org) was used to access and analyze TCGA database [33, 34]. Antigens for the query were based on the previously established antibody specificities of the respective antibodies. M3H67 (pan-MAGE): MAGE-A1, -A3, -A4, -A8, -A10, -B2, -C2 and 57B (MAGE-A3/A4): MAGE-A3, -A4, -A6, -A12. The z-score was set at 2.0.
Immunohistochemistry (IHC) for MAGE-antigens was performed using three antibodies (pan-MAGE: clone M3H67, MAGE-A3/A4: clone 57B, MAGE-A1: clone MA454) on 4μm TMA sections as described previously [10, 35]. Antibody specificities have also been previously analyzed by enzyme linked immunosorbent assay . Based on these previously reported results, we will refer to M3H67 reactivity as pan-MAGE expression, to 57B as MAGE-A3/A4 and to MA454 as MAGE-A1. The following primary antibody concentrations were used: M3H67 = dilution of 1:3324 (0.74 μg/ml), 57B = dilution of 1:150, MA454 = dilution of 1:589 (5.7 μg/ml).
Each tumor core was microscopically evaluated by two independent investigators and was scored as positive for the respective antigen if specific staining was present in ≥ 10% of tumor cells. The tumor sample was classified positive if at least one core of the triplicate was scored as positive in agreement between the two investigators. Additionally, the expression pattern was noted (cytoplasmic OR nuclear (cytORnuc), cytoplasmic AND nuclear (cyt+nuc)).
TMAs were then digitized using the Zeiss MIRAX DESK scanner. Staining intensity per tumor area was determined using a semiautomated, quantitative image analysis software (Definiens Architect XD 1.2, Definiens, Munich, Germany) in order to obtain a continuous spectrum of average staining intensity per tumor area in arbitrary units (maximum range of readout 0.0–1.0) for the complete tumor areas of each core. Thus, expression intensity per area contains both, the stained area of the tumor in relation to the total tumor area and the expression intensity averaged over the total tumor area. The mean expression intensity (MEI) was derived from the mean of up to three available cores for each tumor specimen. This way intratumoral heterogeneity was addressed by taking into account both, the stained tumor area and the staining intensity in relation to the total tumor area for each tumor in three separate tumor cores. The MEI was determined separately for the triplicates of primary tumors, lymph node metastases and recurrences.
IBM SPSS statistics version 21.0 was used for statistical analysis and graphing of results unless otherwise indicated. Cox regression analysis for OS was performed to determine hazard ratios (HR) with a 95% confidence interval (CI). For multivariate cox regression backward likelihood ratio was applied. A Pearson correlation was used to quantify correlations between MAGE expression and known prognostic markers. OS was defined as the time interval from diagnosis until death. OS was analyzed using the Kaplan-Meier method. Groups were compared by log-rank-test. TCGA data were analyzed for OS using the analysis tool cbioportal.org [33, 34].
MEI of primary tumors, lymph node metastases and recurrences was compared and graphed using Kruskal-Wallis test with Graph Pad Prism version 6. MEI in paired samples was compared using the Wilcoxon signed rank test for related samples. Significance was assumed for p < 0.05.
The results published or shown here are in part based upon data generated by the TCGA Research Network:http://cancergenome.nih.gov/.
The authors thank Dr. Gerd Ritter of the Ludwig Institute for Cancer Research at Memorial Sloan Kettering Cancer Center for providing the antibodies for pan-MAGE (M3H67) and MAGE-A1 (MA454).
CONFLICTS OF INTEREST
Simon Laban is a member of an advisory board for Astra Zeneca. Thomas K. Hoffmann is a member of advisory boards for Merck Serono and Merck Sharp & Dohme. Patrick J. Schuler is a member of the advisory board for Bristol Myers Squibb.
Authors’ contributionsS. Laban, Thomas K. Hoffmann and S. Perner planned and supervised the experiments, analyzed and graphed data, performed statistics and wrote the manuscript.
G. Giebel performed the experiments and analyzed data.
G. Giebel and N. Klümper performed the digital image analysis and analyzed data.
S. Perner and A. Schröck established the TMA. G. Kristiansen provided the tissue for TMA generation. A. Schröck and F. Bootz provided clinical data.
J. Thierauf , M.N.Theodoraki, J. Doescher, R. Remark, S. Gnjatic, A. G. Sikora, G. Litjens, N. Grabe, G. Spagnoli, R. Krupar, P. J. Schuler, C. Brunner and J. Brägelmann contributed to the writing of the manuscript, analyzed data and performed statistics.
All authors proof-read the manuscript and contributed to the discussion of results.