Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 19 August 2021
Sec. Head and Neck Cancer

A 14-Marker Multiplexed Imaging Panel for Prognostic Biomarkers and Tumor Heterogeneity in Head and Neck Squamous Cell Carcinoma

Junichi MitsudaJunichi Mitsuda1Takahiro Tsujikawa,*Takahiro Tsujikawa1,2*Kanako YoshimuraKanako Yoshimura1Sumiyo SaburiSumiyo Saburi1Masaho SuetsuguMasaho Suetsugu1Kayo KitamotoKayo Kitamoto1Mari TakenakaMari Takenaka1Gaku OhmuraGaku Ohmura1Akihito AraiAkihito Arai1Hiroshi Ogi,Hiroshi Ogi3,4Kyoko ItohKyoko Itoh3Shigeru HiranoShigeru Hirano1
  • 1Department of Otolaryngology–Head and Neck Surgery, Kyoto Prefectural University of Medicine, Kyoto, Japan
  • 2Department of Cell, Developmental & Cancer Biology, Oregon Health and Science University, Portland, OR, United States
  • 3Department of Pathology and Applied Neurobiology, Kyoto Prefectural University of Medicine, Kyoto, Japan
  • 4SCREEN Holdings Co., Ltd., Kyoto, Japan

Recent advances made in treatment for head and neck squamous cell carcinoma (HNSCC) highlight the need for new prediction tools to guide therapeutic strategies. In this study, we aimed to develop a HNSCC-targeting multiplex immunohistochemical (IHC) panel that can evaluate prognostic factors and the intratumor heterogeneity of HNSCC. To identify IHC-based tissue biomarkers that constitute new multiplex IHC panel, a systematic review and meta-analysis were performed to analyze reported IHC biomarkers in laryngeal and pharyngeal SCC in the period of 2008–2018. The Cancer Genome Atlas (TCGA) and Reactome pathway databases were used to validate the prognostic and functional significance of the identified biomarkers. A 14-marker chromogenic multiplex IHC panel including identified biomarkers was used to analyze untreated HNSCC tissue. Forty-five high-quality studies and thirty-one candidate tissue biomarkers were identified (N = 7062). Prognostic validation in TCGA laryngeal and pharyngeal SCC cohort (N = 205) showed that β-catenin, DKK1, PINCH1, ADAM10, and TIMP1 were significantly associated with poor prognosis, which were related to functional categories such as immune system, cellular response, cell cycle, and developmental systems. Selected biomarkers were assembled to build a 14-marker panel, evaluating heterogeneity and polarized expression of tumor biomarkers in the tissue structures, which was particularly related to activation of Wnt/β-catenin pathway. Integrated IHC analysis based on a systemic review and meta-analysis provides an in situ proteomics tool to assess the aggressiveness and intratumor heterogeneity of HNSCC.

Introduction

The emergence of molecular targeting therapy and immune checkpoint blockade has revolutionized the clinical management of head and neck squamous cell carcinoma (HNSCC), but its mortality is still high at approximately 50%, leading to 330,000 deaths per year worldwide (1). Increased number of available therapeutic options urgently requires optimization of the sequence and combination of treatment modalities, thus requiring new tools for the prediction of therapeutic outcomes.

Immunohistochemistry (IHC) has served as the gold standard for histopathological evaluation of cancer. Since a large number of IHC-based biomarkers have been reported, integrative organization of reported biomarkers is needed to identify clinically meaningful biomarkers related to disease aggressiveness of HNSCC. However, systematic reviews and meta-analyses of IHC-based biomarkers have not been carried out for HNSCC, particularly for pharyngeal and laryngeal SCC, where decision making for treatment of these disease is directly associated with organ preservation, including voice and swallowing function.

Intratumor heterogeneity has been increasingly acknowledged since technological advances in single-cell sequencing methods were made, where the presence of subpopulations of cancer cells with distinct genotypic, phenotypic, and morphological profiles was observed within a single tumor, including HNSCC (24). Since intratumor heterogeneity is strongly associated with therapeutic response and resistance (5), better treatment strategies can be developed based on the characterization of intratumor heterogeneity. Given that transcriptomic profiles do not highly correlate with proteomic changes (6), intratumor heterogeneity needs to be evaluated at both the transcription and protein levels, highlighting the significance of revisiting IHC-based assessment in whole tumor tissue.

Recently, we established a 14-biomarker multiplex IHC platform to enable comprehensive profiling of tumor and immune infiltrates in HNSCC (7). In this study, we performed a systematic review and meta-analysis based on the preferred reporting items for systematic reviews and meta-analyses statement (PRISMA) (8) to identify IHC-based tissue biomarkers associated with the prognosis of pharyngeal and laryngeal SCC. Modified reporting recommendations for tumor marker prognostic studies (REMARK) criteria (9) were utilized for the selection of high-quality studies, followed by prognostic and functional validation using the Cancer Genome Atlas (TCGA) and the Reactome pathway database. This led to the identification of selective biomarkers including β-catenin, DKK1, PINCH1, ADAM10, TIMP1, HIF1α, and ZFX. Based on the above results, we developed a HNSCC-targeting multiplex IHC panel, enabling simultaneous assessment of 14 different tumor and stromal markers with preserved intratumor heterogeneity in a single formalin-fixed paraffin-embedded (FFPE) tissue. In combination with other single-cell and bulk tumor analyses, our chromogenic 14-plex tumor panel can provide spatial information, contributing to the further development of tissue-based biomarkers.

Materials and Methods

Systematic Review and Meta-Analysis (PRISMA)

A literature search was performed in the PubMed database (RRID : SCR_004846) on July 7, 2018, using the following search criteria: (head and neck) AND (cancer OR carcinoma) AND (prognos* OR surviv*) AND (immunohistoch*) AND (laryn* OR pharyn*). Two authors independently inspected the titles and abstracts of the retrieved articles between July 2008 and July 2018. Literature written in languages other than English and not analyzed using IHC data were excluded. To provide adequate measures of predictability and improve the quality of the tumor prognostic biomarker studies, modified REMARK criteria containing key items selected by REMARK criteria (9) were used (Supplementary Table 1).

The following information was obtained: author name, date of publication, number of patients, IHC protocols, associated multivariate hazard ratio (HR), 95% confidence intervals (CIs) and corresponding P values, overall survival (OS), and evaluated biomarkers. For proteins evaluated in a single study only, the summarized HR (including 95% CI) represented the value reported in that study. For proteins that were evaluated in multiple studies, meta-analysis was carried out using software “R version 3.3.3” and package “metafor”. HR and 95% CI were calculated using both the fixed-effect model and random-effect model. Heterogeneity was evaluated in multiple studies by calculating the Q statistic and I2 value. An HR greater than 1 indicated poor prognosis in the study group relative to the reference group and was considered statistically significant if the 95% CI did not overlap. The publication bias was assessed when there were ten or more studies through the inspection of funnel plots.

In Silico Analysis

To evaluate the prognostic significance of the extracted markers, TCGA data (RRID : SCR_003193) was used. The data on 205 laryngeal and pharyngeal SCC cases containing information on the primary site, survival information, OS, and values of biomarker expression were downloaded from the UCSC Xena project (http://xena.ucsc.edu) on September 5, 2019. Based on the median biomarker expression, all cases were divided into high-expression (N = 102) and low-expression groups (N = 103). The 5-year OS was analyzed using the Kaplan-Meier and log-rank tests. No statistical power analyses were used to determine the sample size. The sample size was based on the number of available and qualified tumor samples for this study. A p value of less than 0.05 was considered statistically significant. The functional significance of the extracted biomarkers was evaluated using the Reactome pathway database (https://reactome.org) (RRID : SCR_003485), where molecular reactions are systematically described in molecular detail to generate a network of molecular transformations (10).

Multiplex IHC

Sections (5 µm) were obtained from FFPE samples of untreated tongue, oropharyngeal, hypopharyngeal and laryngeal SCC (N = 4). Multiplex IHC was performed as described previously (7). Briefly, FFPE tumor sections were subjected to sequential detection using 14-different antibodies for tumor and stromal markers. Stained slides were digitally scanned using Nanozoomer S60 (Hamamatsu Photonics) at 20 × magnification. Tissue sections were then stripped of antibody and chromogen by cleaning with alcohol and heat-mediated antigen retrieval, followed by subsequent rounds of antibody staining and imaging. A complete list of antibodies and conditions used for staining are provided in Supplementary Table 2.

Digital image processing was performed as previously described (7). Digital images reflecting the antibody panel were co-registered and aligned. A sequential gating strategy of image cytometry was applied to identify tumor cell phenotypes based on negative cell staining by using FCS Express 6 Image Cytometry Version 6.03.0011 (De Novo Software). All studies involving human tissues were approved by the Institutional Review Board at Kyoto Prefectural University of Medicine (ERB-C-43-4), and written informed consent was obtained from all patients.

Results

A Systematic Review and Meta-Analysis of IHC-Based Studies Identifies 31 Biomarkers Related to Poor Prognosis in Laryngeal and Pharyngeal Cancer

To identify tumor biomarkers that were functionally associated with the prognosis of HNSCC, we performed a systematic review and meta-analysis of previously reported IHC markers in laryngeal and pharyngeal cancers. Our search of the laryngeal or pharyngeal SCC IHC prognostic literature yielded 431 manuscripts for consideration, as shown in Figure 1. 56 studies were not related to pharyngeal or laryngeal SCC. 45 studies were written in other languages. 12 manuscripts were excluded, as protein expression was not assessed using IHC. 12 case reports, 9 reviews, 2 meta-analyses, and a preliminary study were excluded. The remaining 297 studies, which used IHC to measure protein expression in laryngeal and pharyngeal SCC, were assessed in terms of study design. Modified REMARK criteria were applied to assess the completeness of reporting (Supplementary Table 1). Of the 297 cohort studies, 71 failed to publish the HRs and 95% CIs. 67 studies did not report OS. 33 studies did not perform a multivariable log rank analysis. 26 studies were unrelated to the biomarkers. 32 studies were excluded due to a lack of information such as age, gender. Finally, 45 high-quality cohort studies reporting multivariable-adjusted estimates for 45 biomarkers in laryngeal or pharyngeal SCC were extracted (Table 1. N = 7062). Since 5 out of 45 biomarkers (metadherin, ZEB2, p53, Ki-67, and p16) were reported by two or more studies, a meta-analysis was performed based on the fixed-effect model (Figures 2A–E) and random effect model (Supplementary Figures 1A–E), which identified ZEB2 and metadherin as poor prognostic indicators (Figures 2A, B). p16 was associated with favorable OS (HR = 0.58, P = 0.017) (Figure 2E), in agreement with current clinical management, where p16 was utilized as a highly sensitive surrogate marker for prognostically favorable human papilloma virus-related oropharyngeal cancer (16, 17). Publication bias was not evaluable due to the limited number of studies available for funnel plot inspections (56). Following the exclusion of 5 biomarkers with favorable prognosis and 13 biomarkers that did not meet statistical significance, our systematic review and meta-analysis identified 31 biomarkers that were significantly associated with short OS in laryngeal and pharyngeal HNSCC.

FIGURE 1
www.frontiersin.org

Figure 1 A PRISMA flow diagram of the systematic review of IHC-based prognostic biomarkers and multiplex IHC panel selection. (A) A prisma flow diagram presents the PubMed searches, the number of manuscripts excluded, and the identified high-quality cohort studies. (B) Selection strategies based on prognostic and functional significance are applied for establishment of the 14-biomaker multiplex IHC panel.

TABLE 1
www.frontiersin.org

Table 1 IHC biomarkers related to overall survival of laryngeal or pharyngeal SCC.

FIGURE 2
www.frontiersin.org

Figure 2 Forest plots of the effects of five biomarkers evaluated by two or more independent studies on laryngeal and pharyngeal HNSCC survival. (A–E) Forest plots of hazard ratio (HR) for overall survival (OS) of metadherin (A), ZEB2 (B), p53 (C), Ki-67 (D), and p16 (E) are shown based on a fixed effect model. Bars present 95% confidence intervals (CI) of HR, and the center of the lozenge gives the combined HR. Combined fixed effect HRs and tests for heterogeneity (I2) are based on the generic inverse variance (I-V) method.

Functional and Prognostic Significance of Identified IHC Biomarkers Reveals a Broad Spectrum of Tumor Characteristics of HNSCC

The biomarkers identified by the systematic review and meta-analysis were functionally classified into the hallmarks of cancer, such as angiogenesis, evasion of growth suppressors, genome instability and mutation, immune evasion, invasion and metastasis, metabolism, replicable immortality, resisting cell death, sustaining proliferation, and tumor promoting inflammation (57) (Table 1). Details of representative markers classified by hallmarks are described in the Supplementary Information, showing that prognostic IHC biomarkers were associated with a wide range of hallmarks, particularly with increased invasion and metastasis as well as sustained proliferation. Notably, multiple markers, such as DKK1, ADAM10, and ZFX were related to the Wnt/β-catenin pathway (5860).

Furthermore, the functional significance of the identified IHC biomarkers was evaluated using the Reactome pathway database (61), visualizing signaling and metabolic pathways associated with reference gene datasets. A total of 31 genes corresponding to the 31 IHC biomarkers identified by the systematic review and meta-analysis were utilized as a reference dataset, and the results were visualized as a Voronoi diagram based on the ReacFoam function (10), which exhibited an overview of pathways highly related to the identified biomarkers (Figure 3 and Supplementary Figure 2). The immune system, cell cycle, cellular response to external stimuli, and developmental systems were related to the selected 31 biomarkers. β-catenin and related markers were observed in a wide range of functional categories, such as signal transduction, gene expression, and disease, suggesting the possible involvement of Wnt/β-catenin pathways in the progression of HNSCC (Figure 3 and Supplementary Figure 2).

FIGURE 3
www.frontiersin.org

Figure 3 Functional significance of the selected IHC biomarkers identified by the systematic review and meta-analysis. Functional significance of the identified IHC biomarkers is visualized by using the Reactome pathway database (https://reactome.org), reporting on biological pathways associated with dataset of genes. The results are visualized as a Voronoi diagram, showing an overview of the biological pathways related to the 31 biomarkers identified by the systematic review and meta-analysis. Significantly enriched pathways are shown with a color scale from dark to light yellow. IHC biomarkers associated with each pathway are shown.

Since the selected IHC-based biomarkers were related to a broad range of functional properties, we next validated the prognostic significance of these biomarkers using a TCGA cohort in laryngeal and pharyngeal cancer (N = 205). OS was evaluated in relation to representative 31 biomarkers (Figures 4A–G, Supplementary Figures 3A–X and 4), showing that β-catenin, DKK1, TIMP1, PINCH1, and ADAM10 were significantly associated with poor prognosis in patients with laryngeal and pharyngeal SCC (Figures 4A–E). HIF1α and ZFX tended to be associated with a poor prognosis. Interestingly, high expression of MCM7 or COX2 was associated with favorable prognosis in transcriptome analysis, potentially due to the discrepancy between transcriptional and protein levels. Overall, five biomarkers showed prognostic significance in IHC-based analysis as well as during validation by TCGA analysis, which served as candidates for robust predictive biomarkers for the prognosis of laryngeal and pharyngeal cancer.

FIGURE 4
www.frontiersin.org

Figure 4 Prognostic significance of the selected IHC biomarkers identified by the systematic review and meta-analysis. (A–G) Kaplan-Meier analyses of overall survival (OS) in the laryngeal and pharyngeal HNSCC cohort of the Cancer Genome Atlas (TCGA) (N = 205) stratified by gene expression of the identified biomarkers are shown. Median is used for the cutoff values. Statistical significance is determined using log-rank test.

14-Biomarker Multiplex IHC Panel Highlights Intratumor Heterogeneity in HNSCC

To expand our multiplex IHC platform, we established an HNSCC-targeting multiplex IHC panel assembling candidate biomarkers derived from our systematic review and meta-analysis. Based on a chromogenic sequential IHC method enabling analysis of more than 12 proteins in one tissue section (7), a 14-biomarker multiplex IHC panel was assembled from β-catenin, DKK1, PINCH1, ADAM10, and TIMP1 used as prognostic markers. In addition to hematoxylin staining, HIF1α as a hypoxia marker, ZFX as a transcriptional regulator of self-renewal and stemness, ZEB2 as an epithelial to mesenchymal transformation (EMT) marker, and Ki67 as a proliferation marker were included in the panel, considering the functional significance of these markers. To simultaneously analyze tumor microenvironment-related factors, the panel also included pan-cytokeratin (pCK) as an epithelial marker, CD3 and CD68 as immune markers, and alpha smooth muscle actin (aSMA) as a fibroblastic and vascular endothelial marker (Figure 5A). Negative controls were used for confirmation of complete antibody and signal stripping by omitting the primary antibody in the next cycle (data not shown). This panel enabled simultaneous evaluation of tumor, immune, and stroma biomarkers in a single slide, visualizing the opposite distribution of HIF1α and TIMP1 in the intratumor and invasive margins, respectively, (Figures 5A, B) and accumulation of immune markers (CD3 and CD68) and stroma markers (aSMA) in adjacent non-malignant regions (Figure 5C and Supplementary Figure 5).

FIGURE 5
www.frontiersin.org

Figure 5 Fourteen-marker multiplex IHC panel to visualize prognostic tumor cell biomarkers and immune/stromal cells in a single FFPE tissue. A FFPE tissue of tongue SCC is visualized by the established 14-marker multiplex IHC. (A) Representative 12-color composite images are shown. Biomarkers and colors are shown at the bottom. Scale bar = 100 µm. (B) Top and bottom panels present expression of HIF1α and TIMP1 in the intratumor and invasive margin regions, respectively. Blue and white arrows depict HIF1α+ cells and TIMP1+ cells, respectively. (C) Immune markers (CD3 and CD68) and aSMA+ cells in stroma region are shown. White, red, and green arrows present CD3+, CD68+, and aSMA+ cells, respectively. Biomarkers and colors are shown at the right. Scale bars = 10 µm. Corresponding single-marker images are shown in Supplementary Figure 3.

Next, a single quantitative assessment was performed using an image cytometry (Figure 6A). Since image cytometry analysis enables single-cell analysis with preserved spatial information (62), cell densities and phenotypes of cancer, immune and stromal cells were spatially evaluated at different proximities to the invasive margin of stroma versus intratumor regions (Figures 6B–E). Spatial analysis revealed that some tumor cell biomarkers exhibited polarized localization in the tissue structure. For example, increased expression of HIF1α was observed in a distant area from the invasive margin (Figure 6B), suggesting the presence of hypoxic conditions at the center of the tumor. In contrast, strong expression of ZFX was observed at the invasive margin (Figure 6C), suggesting the presence of heterogeneity in increased survival and self-renewing capabilities of tumor cells via ZFX (63). The expression of β-catenin and ZEB2 exhibited characteristic polarized localization in the center and invasive margin, respectively (Figures 6B, C). Simultaneously, spatial analysis of T cells (CD3+), tumor-associated macrophages (CD68+), and aSMA+ cells revealed heterogeneity of immune and stromal cell densities (Figure 6D).

FIGURE 6
www.frontiersin.org

Figure 6 The 14-biomarker multiplex IHC panel quantitatively assesses intratumoral heterogeneity of HNSCC. (A) The gating strategy used for single cell-based chromogenic signal intensity is shown where cell size/area, and location are utilized for density plots similar to flow cytometry. X and y axes were shown on a logarithmic scale. The gating strategies identified cell populations of pCK+ tumor cells with selected tumor cell biomarkers such as β-catenin, DKK1, ADAM10, ZEB2, ZFX, HIF1α, TIMP1, and PINCH1, together with aSMA+ cells (pCKaSMA+), T cells (pCKaSMACD3+) and macrophages (pCKaSMACD68+) in a single tissue. (B–D) Left panels depict contour mapping based on cell density of each marker. Dash lines present the invasive margin between tumor cell nests and stroma. Right panels depict cell densities of identified lineages in every 200 µm-width from stroma to the tumor side. Tumor cell biomarkers are classified into tumor core-polarized (B), and invasive margin-polarized markers, in comparison with immune and stroma markers (D). (E) IHC image and pixel density contour graph indicate distribution of colocalized markers at the invasive margin. Line graph show two colocalized markers. Black line shows β-catenin+ DKK1+ ADAM10+ ZFX+ cells and green line showed β-catenin+ DKK1+ ADAM10+ ZFX+ ZEB2+ cells.

Furthermore, this panel revealed the potential colocalization of multiple tumor cell markers such as β-catenin, ADAM10, ZFX, and ZEB2 (Figure 6E). Interestingly, ZEB2, an EMT marker, was strongly colocalized with Wnt/β-catenin related markers at the invasive margin (Figure 6E). Given that the Wnt/β-catenin pathway is associated with cancer cell aggressiveness, and ZEB2 is presumably related to metastatic capability, multiple biomarker-positive cancer cells might be related to the presence of aggressive subpopulations of cancer cells, which are potentially localized at the invasive margin of HNSCC tumors. Taken together, these findings indicate that our 14-biomarker multiplex IHC panel provides a platform for the quantitative assessment of intratumor heterogeneity based on multiple tissue biomarkers identified by a systematic review and meta-analysis.

Discussion

In this study, we developed a 14-marker multiplexed imaging panel for prognostic biomarkers identified by a systematic review and meta-analysis of prognostically and functionally-relevant IHC-based biomarkers in laryngeal and pharyngeal SCC. IHC analysis using this 14-marker multiplexed panel revealed intratumor heterogeneity that was especially related to an aggressive subpopulation associated with Wnt/β-catenin pathway in HNSCC.

While some of the identified biomarkers such as p53 and Ki-67 are shared among many other cancer types (6468), our systematic review of laryngeal and pharyngeal HNSCC characteristically identified the prognostic significance of Wnt/β-catenin pathway-related markers such as β-catenin, DKK1, ADAM10, and ZFX (Table 1 and Figure 4) (27, 47, 58, 67, 68). This trend was more pronounced in laryngeal SCC as observed in Supplementary Figure 4. These Wnt/β-catenin-related markers are involved in a wide range of functional categories in the Reactome pathway analysis (Figure 3 and Supplementary Figure 2), suggesting that activated β-catenin might be related to disease aggressiveness of HNSCC, particularly in laryngeal SCC.

Although the importance of intratumor heterogeneity has been increasingly acknowledged, the vast majority of IHC-based studies have not considered the significance of heterogeneous expression of tumor biomarkers. Intratumor heterogeneity can be associated with poor clinical outcomes due to the presence of subclones of cancer cells with aggressive characteristics, such as rapid proliferation, metastatic potential, and therapeutic resistance (2). Recent studies based on single-cell RNA sequencing have revealed that HNSCC has remarkable intratumor heterogeneity in terms of clonal diversity (69) and locally activated partial EMT program at the leading edge of primary tumors (70). In agreement with these findings, our multiplex IHC analysis demonstrated the localization of aggressive cancer cell populations expressing Wnt/β-catenin related markers and ZEB2 in the invasive margin (Figure 6E). Given that TIMP1 and macrophages were reported to induce the expression of EMT transcriptional factors, including ZEB2 (26, 71), increased expression of TIMP1 and CD68 in the invasive margin might be related to the regulation of EMT and metastatic potential. Multiplex IHC-based simultaneous detection in a single tissue enables quantitative evaluation of multiple overlapping biomarkers together with spatial distribution in the tumor architecture, contributing to the understanding of intratumor heterogeneity of HNSCC.

Immune checkpoint blockade has gained increasing importance in the treatment of HNSCC, and many studies have investigated the association between specific immune cell populations and clinical outcomes of HNSCC (72, 73). Despite the single-case analysis in our study, increased expression of β-catenin in the center of the tumor might be associated with the low density of CD3+ T cells (Figure 6D), which is potentially associated with T-cell exclusion mechanisms driven by Wnt/β-catenin pathway activation (74). These findings suggest that our 14-biomarker multiplex IHC panel can serve as a platform to characterize the intratumor heterogeneity of cancer cell phenotypes in association with immune cell densities.

In this study, we observed discrepancies between IHC-based results and transcriptomic analyses in terms of prognostic significance. Notably, 15 out of 45 biomarkers exhibited poor prognostic significance in IHC studies, but not in transcriptomic analysis in TCGA (Table 1 and Figure 4), leading to the notion that transcriptomic profiles did not necessarily predict protein expression (75). These observations provide a rationale for revisiting IHC-based proteomics, which can provide a clinically applicable approach for analyzing the intensity and functional properties of protein expression without significant cost (76).

Although extensive research was performed, a potential limitation of this study is publication bias. Due to a great diversity of IHC biomarkers, identified biomarkers in our study were based on a small number of the original studies (Figure 2), in analogous to previous IHC meta-analyses (64, 66). To avoid bias due to any potential heterogeneity, we utilized TCGA database to validate the prognostic significance of the identified biomarkers (Figure 4). Another limitation in this study is the lack of subgroup analyses stratified by the HNSCC subsites such as larynx, oropharynx, and hypopharynx. Since we observed some interesting differences between laryngeal and pharyngeal SCC (Supplementary Figure 4), further studies are required for exploring subsite-specific biomarkers associated with clinical outcomes.

In summary, the present study, based on a systematic review and meta-analysis of laryngeal and pharyngeal SCC, identified IHC-based prognostic biomarkers. Established multiplex IHC panels enable assessment of intratumor heterogeneity of HNSCC, which is particularly related to activation of the Wnt/β-catenin pathway. In combination with other single-cell and bulk tumor analyses, our tissue-based biomarker analysis can provide new insights into diagnostic and treatment strategies for HNSCC.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by The Institutional Review Board at Kyoto Prefectural University of Medicine (ERB-C-43-4). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

JM performed manuscript preparation. JM and TT performed data analysis. JM, MS, and KK collected data and performed IHC. TT conceived of the research and designed research studies. TT, MT, GO, and AA obtained research funding. TT, GO, and AA participated in sample preparation. HO and KI provided digital analytic support for multiplex IHC. TT and SH supervised this study. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the Japanese ministry of education, culture, sports, science and technology (17H07016, 19K18814, 19K18739, and 19K23893), and the public promoting association Asano foundation for studies on medicine, and by the research promotion award from the Oto-Rhino-Laryngological Society of Japan, Inc. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.

Conflict of Interest

HO is an employee for SCREEN Holdings Co., Ltd. KI has received research funding from SCREEN Holdings Co., Ltd.

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

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.713561/full#supplementary-material

References

1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2018) 68:394–424. doi: 10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Mroz EA, Tward AD, Pickering CR, Myers JN, Ferris RL, Rocco JW. High Intratumor Genetic Heterogeneity Is Related to Worse Outcome in Patients With Head and Neck Squamous Cell Carcinoma. Cancer (2013) 119:3034–42. doi: 10.1002/cncr.28150

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Puram SV, Rocco JW. Molecular Aspects of Head and Neck Cancer Therapy. Hematol Oncol Clin North Am (2015) 29:971–92. doi: 10.1016/j.hoc.2015.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Mroz EA, Rocco JW. Intra-Tumor Heterogeneity in Head and Neck Cancer and Its Clinical Implications. World J Otorhinolaryngol Neck Surg (2016) 2:60–7. doi: 10.1016/j.wjorl.2016.05.007

CrossRef Full Text | Google Scholar

5. Easwaran H, Tsai HC, Baylin SB. Cancer Epigenetics: Tumor Heterogeneity, Plasticity of Stem-Like States, and Drug Resistance. Mol Cell (2014) 54:716–27. doi: 10.1016/j.molcel.2014.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Vogel C, Marcotte EM. Insights Into the Regulation of Protein Abundance From Proteomic and Transcriptomic Analyses. Nat Rev Genet (2012) 13:227–32. doi: 10.1038/nrg3185

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Tsujikawa T, Kumar S, Borkar RN, Azimi V, Thibault G, Chang YH, et al. Quantitative Multiplex Immunohistochemistry Reveals Myeloid-Inflamed Tumor-Immune Complexity Associated With Poor Prognosis. Cell Rep (2017) 19:203–17. doi: 10.1016/j.celrep.2017.03.037

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Moher D, Liberati A, Tetzlaff J, Altman DG, Altman D, Antes G, et al. Preferred Reporting Items for Systematic Reviews and Meta-Analyses: The PRISMA Statement. PloS Med (2009) 6:e1000097. doi: 10.1371/journal.pmed.1000097

PubMed Abstract | CrossRef Full Text | Google Scholar

9. McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M, Clark GM. Reporting Recommendations for Tumour Marker Prognostic Studies (REMARK). Br J Cancer (2005) 93:387–91. doi: 10.1038/sj.bjc.6602678

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Jassal B, Matthews L, Viteri G, Gong C, Lorente P, Fabregat A, et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res (2020) 48:D498–503. doi: 10.1093/nar/gkz1031

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Li J-J, Yang X-M, Wang S-H, Tang Q-L. Prognostic Role of Epidermal Growth Factor-Like Domain 7 Protein Expression in Laryngeal Squamous Cell Carcinoma. J Laryngol Otol (2011) 125:1152–7. doi: 10.1017/S0022215111002441

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Almadori G, Lauriola L, Coli A, Bussu F, Gallus R, Scannone D, et al. Minichromosome Maintenance Protein 7 and Geminin Expression: Prognostic Value in Laryngeal Squamous Cell Carcinoma in Patients Treated With Radiotherapy and Cetuximab. Head Neck (2017) 39:684–93. doi: 10.1002/hed.24670

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Schrijvers ML, van der laan BFAM, de bock GH, Pattje WJ, Mastik MF, Menkema L, et al. Overexpression of Intrinsic Hypoxia Markers HIF1α and CA-IX Predict for Local Recurrence in Stage T1-T2 Glottic Laryngeal Carcinoma Treated With Radiotherapy. IntJRadiation Oncol BiolPhys (2008) 72(1):161–9. doi: 10.1016/j.ijrobp.2008.05.025

CrossRef Full Text | Google Scholar

14. Zhu M, Yin F, Yang L, Chen S, Chen R, Zhou W, et al. Contribution of TIP30 to Chemoresistance in Laryngeal Carcinoma. Cell Death Dis (2014) 5:e1468–10. doi: 10.1038/cddis.2014.424

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Baek SK, Jung KY, Lee SH, Woo JS, Kwon SY, Chung EJ, et al. Prognostic Significance of Vascular Endothelial Growth Factor-C Expression and Lymphatic Vessel Density in Supraglottic Squamous Cell Carcinoma. Laryngoscope (2009) 119:1325–30. doi: 10.1002/lary.20483

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Salazar CR, Anayannis N, Smith RV, Wang Y, Jr Haigentz M, Garg M, et al. Combined P16 and Human Papillomavirus Testing Predicts Head and Neck Cancer Survival. Int J Cancer (2014) 135:2404–12. doi: 10.1002/ijc.28876

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Golusiński P, Pazdrowski J, Szewczyk M, Misiołek M, Pietruszewska W, Klatka J, et al. Is Immunohistochemical Evaluation of P16 in Oropharyngeal Cancer Enough to Predict the HPV Positivity? Rep Pract Oncol Radiother (2017) 22:237–42. doi: 10.1016/j.rpor.2017.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Wang H, Zhang Z, Sun R, Lin H, Gong L, Fang M, et al. HPV Infection and Anemia Status Stratify the Survival of Early T2 Laryngeal Squamous Cell Carcinoma. J Voice (2015) 29:356–62. doi: 10.1016/j.jvoice.2014.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Lehnerdt GF, Bachmann HS, Adamzik M, Panic A, Köksal E, Weller P, et al. AQP1, AQP5, Bcl-2 and P16 in Pharyngeal Squamous Cell Carcinoma. J Laryngol Otol (2015) 129:580–6. doi: 10.1017/S002221511500119X

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Young RJ, Urban D, Angel C, Corry J, Lyons B, Vallance N, et al. Frequency and Prognostic Significance of P16 INK4A Protein Overexpression and Transcriptionally Active Human Papillomavirus Infection in Laryngeal Squamous Cell Carcinoma. Br J Cancer (2015) 112:1098–104. doi: 10.1038/bjc.2015.59

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Descamps G, Karaca Y, Lechien JR, Kindt N, Decaestecker C, Remmelink M, et al. Classical Risk Factors, But Not HPV Status, Predict Survival After Chemoradiotherapy in Advanced Head and Neck Cancer Patients. J Cancer Res Clin Oncol (2016) 142:2185–96. doi: 10.1007/s00432-016-2203-7

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Qiu X, You Y, Huang J, Wang X, Zhu H, Wang Z. LAMP3 and TP53 Overexpression Predicts Poor Outcome in Laryngeal Squamous Cell Carcinoma. Int J Clin Exp Pathol (2015) 8:5519–27.

PubMed Abstract | Google Scholar

23. de Miguel-Luken MJ, Chaves-Conde M, de Miguel-Luken V, Muñoz-Galván S, López-Guerra JL, Mateos JC, et al. MAP17 (PDZKIP1) as a Novel Prognostic Biomarker for Laryngeal Cancer. Oncotarget (2015) 6:12625–36. doi: 10.18632/oncotarget.3470

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Re M, Zizzi A, Ferrante L, Stramazzotti D, Goteri G, Gioacchini FM, et al. P63 and Ki-67 Immunostainings in Laryngeal Squamous Cell Carcinoma Are Related to Survival. Eur Arch Oto-Rhino-Laryngol (2014) 271:1641–51. doi: 10.1007/s00405-013-2833-1

CrossRef Full Text | Google Scholar

25. Han L, Jiang B, Wu H, Zhang S, Lu X. Expression and Prognostic Value of MAGE-A9 in Laryngeal Squamous Cell Carcinoma. Int J Clin Exp Pathol (2014) 7:6734–42.

PubMed Abstract | Google Scholar

26. Zhu GJ, Song PP, Zhou H, Shen XH, Wang JG, Ma XF, et al. Role of Epithelial-Mesenchymal Transition Markers E-Cadherin, N-Cadherin, β-Catenin and ZEB2 in Laryngeal Squamous Cell Carcinoma. Oncol Lett (2018) 15:3472–81. doi: 10.3892/ol.2018.7751

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Shi Y, Gong H-L, Zhou L, Tian J, Wang Y. Dickkopf-1 Is a Novel Prognostic Biomarker for Laryngeal Squamous Cell Carcinoma. Acta Otolaryngol (2014) 134:753–9. doi: 10.3109/00016489.2014.894251

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Zhang Q, Zhang W, Zhang J, Xu H, You Y. Aberrant Kif2a and Ki67 Expression Predicts Poor Survival in Laryngeal Squamous Cell Carcinoma. Auris Nasus Larynx (2016) 43:433–9. doi: 10.1016/j.anl.2015.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Lu M, Zhu H, Wang X, Zhang D, Xiong L, Zhu J, et al. LAMP1 Expression Is Associated With Malignant Behaviours and Predicts Unfavourable Prognosis in Laryngeal Squamous Cell Carcinoma. Pathology (2016) 48:684–90. doi: 10.1016/j.pathol.2016.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Ma J, Wang J, Fan W, Pu X, Zhang D, Fan C, et al. Upregulated TIMP-1 Correlates With Poor Prognosis of Laryngeal Squamous Cell Carcinoma. Int J Clin Exp Pathol (2014) 7:246–54.

PubMed Abstract | Google Scholar

31. Wu H, Xu H, Zhang S, Wang X, Zhu H, Zhang H, et al. Potential Therapeutic Target and Independent Prognostic Marker of TROP2 in Laryngeal Squamous Cell Carcinoma. Head Neck (2013) 35:1373–8. doi: 10.1002/hed.23138

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Jouppila-Mättö A, Mannermaa A, Sironen R, Kosma V-M, Soini Y, Pukkila M. SIP1 Predicts Progression and Poor Prognosis in Pharyngeal Squamous Cell Carcinoma. Histol Histopathol (2015) 30:569–79. doi: 10.14670/HH-30.569

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Yu C, Liu Y, Tan H, Li G, Su Z, Ren S, et al. Metadherin Regulates Metastasis of Aquamous Cell Carcinoma of the Head and Neck via AKT Signaling Pathway-Mediated Epithelial-Mesenchymal Transition. Cancer Lett (2014) 343:258–67. doi: 10.1016/j.canlet.2013.09.033

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Liu Y, Su Z, Li G, Yu C, Ren S, Huang D, et al. Increased Expression of Metadherin Protein Predicts Worse Disease-Free and Overall Survival in Laryngeal Squamous Cell Carcinoma. Int J Cancer (2013) 133:671—679. doi: 10.1002/ijc.28071

CrossRef Full Text | Google Scholar

35. Zhang Q, Xu H, You Y, Zhang J, Chen R. High Gpx1 Expression Predicts Poor Survival in Laryngeal Squamous Cell Carcinoma. Auris Nasus Larynx (2018) 45:13–9. doi: 10.1016/j.anl.2017.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kaira K, Toyoda M, Shimizu A, Imai H, Sakakura K, Nikkuni O, et al. Decreasing Expression of Glucose-Regulated Protein GRP78/Bip as a Significant Prognostic Predictor in Patients With Advanced Laryngeal Squamous Cell Carcinoma. Head Neck (2016) 38:1539–44. doi: 10.1002/hed.24471

PubMed Abstract | CrossRef Full Text | Google Scholar

37. García-Carracedo D, Ángeles Villaronga M, Álvarez-Teijeiro S, Hermida-Prado F, Santamaría I, Allonca E, et al. Impact of PI3K/AKT/Mtor Pathway Activation on the Prognosis of Patients With Head and Neck Squamous Cell Carcinomas. Oncotarget (2016) 7:29780–93. doi: 10.18632/oncotarget.8957

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Schrijvers ML, Pattje WJ, Slagter-Menkema L, Mastik MF, Gibcus JH, Langendijk JA, et al. FADD Expression as a Prognosticator in Early-Stage Glottic Squamous Cell Carcinoma of the Larynx Treated Primarily With Radiotherapy. Int J Radiat Oncol Biol Phys (2012) 83:1220–6. doi: 10.1016/j.ijrobp.2011.09.060

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Ma J, Liu S, Zhang W, Zhang F, Wang S, Wu L, et al. High Expression of NDRG3 Associates With Positive Lymph Node Metastasis and Unfavourable Overall Survival in Laryngeal Squamous Cell Carcinoma. Pathology (2016) 48:691–6. doi: 10.1016/j.pathol.2016.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Shen B, Li D, Dong P, Gao S. Expression of ABC Transporters Is an Unfavorable Prognostic Factor in Laryngeal Squamous Cell Carcinoma. Ann Otol Rhinol Laryngol (2011) 120:820–7. doi: 10.1177/000348941112001208

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Schrader CH, Kolb M, Zaoui K, Flechtenmacher C, Grabe N, Weber KJ, et al. Kallikrein-Related Peptidase 6 Regulates Epithelial-to-Mesenchymal Transition and Serves as Prognostic Biomarker for Head and Neck Squamous Cell Carcinoma Patients. Mol Cancer (2015) 14:1–14. doi: 10.1186/s12943-015-0381-6

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Tsinias G, Nikou S, Papadas T, Pitsos P, Papadaki H, Bravou V. High PINCH1 Expression in Human Laryngeal Carcinoma Associates With Poor Prognosis. Anal Cell Pathol (2018) 2018:1–12. doi: 10.1155/2018/2989635

CrossRef Full Text | Google Scholar

43. Liu XK, Li Q, Xu LH, Hu LJ, Liao WG, Zhang XR, et al. Expression and Clinical Significance of SIAH in Laryngeal Squamous Cell Carcinoma. Med Oncol (2013) 30:1–8. doi: 10.1007/s12032-013-0485-z

CrossRef Full Text | Google Scholar

44. Tang X, Shen X, Li L, Zhang Y, Chen G. SOX2 Overexpression Correlates With Poor Prognosis in Laryngeal Squamous Cell Carcinoma. Auris Nasus Larynx (2013) 40:481–6. doi: 10.1016/j.anl.2013.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Yang F, Ma H, Feng L, Lian M, Wang R, Fan E, et al. Zinc Finger Protein X-Linked (ZFX) Contributes to Patient Prognosis, Cell Proliferation and Apoptosis in Human Laryngeal Squamous Cell Carcinoma. Int J Clin Exp Pathol (2015) 8:13886–99.

PubMed Abstract | Google Scholar

46. Ma H, Du X, Zhang S, Wang Q, Yin Y, Qiu X, et al. Achaete-Scute Complex Homologue-1 Promotes Development of Laryngocarcinoma via Facilitating the Epithelial-Mesenchymal Transformation. Tumor Biol (2017) 39:1010428317705752. doi: 10.1177/1010428317705752

CrossRef Full Text | Google Scholar

47. You B, Gu M, Cao X, Li X, Shi S, Shan Y, et al. Clinical Significance of ADAM10 Expression in Laryngeal Carcinoma. Oncol Lett (2017) 13:1353–9. doi: 10.3892/ol.2016.5546

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Wang R, Guo Y, Ma H, Feng L, Wang Q, Chen X, et al. Tumor Necrosis Factor Superfamily Member 13 Is a Novel Biomarker for Diagnosis and Prognosis and Promotes Cancer Cell Proliferation in Laryngeal Squamous Cell Carcinoma. Tumor Biol (2016) 37:2635–45. doi: 10.1007/s13277-015-4016-8

CrossRef Full Text | Google Scholar

49. Zhu G, Liu L, She L, Tan H, Wei M, Chen C, et al. Elevated Expression of Histone Demethylase PHF8 Associates With Adverse Prognosis in Patients of Laryngeal and Hypopharyngeal Squamous Cell Carcinoma. Epigenomics (2015) 7:143–53. doi: 10.2217/epi.14.82

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Yang L, Wang H, Wang Y, He Z, Chen H, Liang S, et al. Prostate Tumor Overexpressed-1, in Conjunction With Human Papillomavirus Status, Predicts Outcome in Early-Stage Human Laryngeal Squamous Cell Carcinoma. Oncotarget (2016) 7:31878–91. doi: 10.18632/oncotarget.8103

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Zhou X, Liu Y, Tan G. Prognostic Value of Elevated SHIP2 Expression in Laryngeal Squamous Cell Carcinoma. Arch Med Res (2011) 42:589–95. doi: 10.1016/j.arcmed.2011.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Yang H, Jiang WQ, Cao Y, Sun YA, Wei J, An X, et al. Elevated ZNF703 Protein Expression Is an Independent Unfavorable Prognostic Factor for Survival of the Patients With Head and Neck Squamous Cell Carcinoma. Dis Markers (2015) 2015:1–178. doi: 10.1155/2015/640263

CrossRef Full Text | Google Scholar

53. Zhang J, Cheng L, Zhou L. Prognostic Significance of Γ-H2AX in Laryngeal Squamous Cell Carcinoma After Surgery. Chin (Engl) (2014) 127:2664–7.

Google Scholar

54. Chen YF, Luo RZ, Li Y, Cui BK, Song M, Yang AK, et al. High Expression Levels of COX-2 and P300 Are Associated With Unfavorable Survival in Laryngeal Squamous Cell Carcinoma. Eur Arch Oto-Rhino-Laryngol (2013) 270:1009–17. doi: 10.1007/s00405-012-2275-1

CrossRef Full Text | Google Scholar

55. Yan M, Xu Q, Zhang P, Zhou XJ, Zhang ZY, Chen WT. Correlation of NF-κb Signal Pathway With Tumor Metastasis of Human Head and Neck Squamous Cell Carcinoma. BMC Cancer (2010) 10:1–13. doi: 10.1186/1471-2407-10-437

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Sutton AJ, Duval SJ, Tweedie RL, Abrams KR, Jones DR. Empirical Assessment of Effect of Publication Bias on Meta-Analyses. Br Med J (2000) 320:1574–7. doi: 10.1136/bmj.320.7249.1574

CrossRef Full Text | Google Scholar

57. Hanahan D, Weinberg RA. Hallmarks of Cancer: The Next Generation. Cell (2011) 144:646–74. doi: 10.1016/j.cell.2011.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Chen ZY, Du Y, Wang L, Liu XH, Guo J, Weng XD. MiR-543 Promotes Cell Proliferation and Metastasis of Renal Cell Carcinoma by Targeting Dickkopf 1 Through the Wnt/β-Catenin Signaling Pathway. J Cancer (2018) 9:3660–8. doi: 10.7150/jca.27124

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Anders L, Mertins P, Lammich S, Murgia M, Hartmann D, Saftig P, et al. Furin-, ADAM 10-, and Γ-Secretase-Mediated Cleavage of a Receptor Tyrosine Phosphatase and Regulation of β-Catenin's Transcriptional Activity. Mol Cell Biol (3917) 2006) 26:3917–34. doi: 10.1128/MCB.26.10.3917-3934.2006

CrossRef Full Text | Google Scholar

60. Wang C, Fu SY, da Wang M, Yu WB, Cui QS, Wang HR, et al. Zinc Finger Protein X-Linked Promotes Expansion of Epcam+ Cancer Stem-Like Cells in Hepatocellular Carcinoma. Mol Oncol (2017) 11:455–69. doi: 10.1002/1878-0261.12036

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Joshi-Tope G, Gillespie M, Vastrik I, D’Eustachio P, Schmidt E, de Bono B, et al. Reactome: A Knowledgebase of Biological Pathways. Nucleic Acids Res (2005) 33:428–32. doi: 10.1093/nar/gki072

CrossRef Full Text | Google Scholar

62. Tsujikawa T, Thibault G, Azimi V, Sivagnanam S, Banik G, Means C, et al. Robust Cell Detection and Segmentation for Image Cytometry Reveal Th17 Cell Heterogeneity. Cytom Part A (2019) 95:389–98. doi: 10.1002/cyto.a.23726

CrossRef Full Text | Google Scholar

63. Yin J, Jiang Y, Wu H, Wang J, Zhang S, Liu H. Overexpression of ZFX and Its Involvement in Squamous Cell Carcinoma of the Tongue. Oncol Rep (2015) 33:141–8. doi: 10.3892/or.2014.3572

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Oliveira LR, Ribeiro-Silva A. Prognostic Significance of Immunohistochemical Biomarkers in Oral Squamous Cell Carcinoma. Int J Oral Maxillofac Surg (2011) 40:298–307. doi: 10.1016/j.ijom.2010.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Rainsbury JW, Ahmed W, Williams HK, Roberts S, Paleri V, Mehanna H. Prognostic Biomarkers of Survival in Oropharyngeal Squamous Cell Carcinoma: Systematic Review and Meta-Analysis. Head Neck (2013) 35:1048–55. doi: 10.1002/hed.22950

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Jamieson NB, Carter CR, McKay CJ, Oien KA. Tissue Biomarkers for Prognosis in Pancreatic Ductal Adenocarcinoma: A Systematic Review and Meta-Analysis. Clin Cancer Res (2011) 17:3316–31. doi: 10.1158/1078-0432.CCR-10-3284

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Bouvet M, Ellis LM, Nishizaki M, Fujiwara T, Liu W, Bucana CD, et al. Adenovirus-Mediated Wild-Type P53 Gene Transfer Down-Regulates Vascular Endothelial Growth Factor Expression and Inhibits Angiogenesis in Human Colon Cancer. Cancer Res (1998) 58(11):2288–92.

PubMed Abstract | Google Scholar

68. Yerushalmi R, Woods R, Ravdin PM, Hayes MM, Gelmon KA. Ki67 in Breast Cancer: Prognostic and Predictive Potential. Lancet Oncol (2010) 11:174–83. doi: 10.1016/S1470-2045(09)70262-1

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Zhang XC, Xu C, Mitchell RM, Zhang B, Zhao D, Li Y, et al. Tumor Evolution and Intratumor Heterogeneity of an Oropharyngeal Squamous Cell Carcinoma Revealed by Whole-Genome Sequencing. Neoplasia (United States) (2013) 15:1371–8. doi: 10.1593/neo.131400

CrossRef Full Text | Google Scholar

70. Puram SV, Tirosh I, Parikh AS, Patel AP, Yizhak K, Gillespie S, et al. Single-Cell Transcriptomic Analysis of Primary and Metastatic Tumor Ecosystems in Head and Neck Cancer. Cell (2017) 171:1611–24.e24. doi: 10.1016/j.cell.2017.10.044

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Lee SH, Koo BS, Kim JM, Huang S, Rho YS, Bae WJ, et al. Wnt/β-Catenin Signalling Maintains Self-Renewal and Tumourigenicity of Head and Neck Frontsquamous Cell Carcinoma Stem-Like Cells by Activating Oct4. J Pathol (2014) 234:99–107. doi: 10.1002/path.4383

PubMed Abstract | CrossRef Full Text | Google Scholar

72. de Ruiter EJ, Ooft ML, Devriese LA, Willems SM. 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:1–10. doi: 10.1080/2162402X.2017.1356148

CrossRef Full Text | Google Scholar

73. De Meulenaere A, Vermassen T, Aspeslagh S, Vandecasteele K, Rottey S, Ferdinande L. Tils in Head and Neck Cancer: Ready for Clinical Implementation and Why (Not)? Head Neck Pathol (2017) 11:354–63. doi: 10.1007/s12105-016-0776-8

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Luke JJ, Bao R, Sweis RF, Spranger S, Gajewski TF. WNT/B-Catenin Pathway Activation Correlates With Immune Exclusion Across Human Cancers. Clin Cancer Res (2019) 25:3074–83. doi: 10.1158/1078-0432.CCR-18-1942

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Akbani R, Ng PKS, Werner HMJ, Shahmoradgoli M, Zhang F, Ju Z, et al. A Pan-Cancer Proteomic Perspective on the Cancer Genome Atlas. Nat Commun (2014) 5:1–15. doi: 10.1038/ncomms4887

CrossRef Full Text | Google Scholar

76. Uhlén M, Björling E, Agaton C, Szigyarto CAK, Amini B, Andersen E, et al. A Human Protein Atlas for Normal and Cancer Tissues Based on Antibody Proteomics. Mol Cell Proteomics (2005) 4:1920–32. doi: 10.1074/mcp.M500279-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: meta-analysis, biomarker, immunohistochemistry, tumor heterogeneity, head and neck squamous cell carcinoma (HNSCC)

Citation: Mitsuda J, Tsujikawa T, Yoshimura K, Saburi S, Suetsugu M, Kitamoto K, Takenaka M, Ohmura G, Arai A, Ogi H, Itoh K and Hirano S (2021) A 14-Marker Multiplexed Imaging Panel for Prognostic Biomarkers and Tumor Heterogeneity in Head and Neck Squamous Cell Carcinoma. Front. Oncol. 11:713561. doi: 10.3389/fonc.2021.713561

Received: 08 June 2021; Accepted: 15 July 2021;
Published: 19 August 2021.

Edited by:

Makoto Tahara, National Cancer Center Hospital East, Japan

Reviewed by:

Ryuji Yasumatsu, Kyushu University, Japan
Leandro Luongo Matos, Universidade de São Paulo, Brazil

Copyright © 2021 Mitsuda, Tsujikawa, Yoshimura, Saburi, Suetsugu, Kitamoto, Takenaka, Ohmura, Arai, Ogi, Itoh and Hirano. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Takahiro Tsujikawa, tu-ji@koto.kpu-m.ac.jp

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.