Dissecting Tissue Compartment-Specific Protein Signatures in Primary and Metastatic Oropharyngeal Squamous Cell Carcinomas

Abstract

Head and neck squamous cell carcinoma (HNSCC) often presents with locoregional or distant disease, despite multimodal therapeutic approaches, which include surgical resection, chemoradiotherapy, and more recently, immunotherapy for metastatic or recurrent HNSCC. Therapies often target the primary and nodal regional HNSCC sites, and their efficacy at controlling occult distant sites remains poor. While our understanding of the tumor microenvironment conducive to effective therapies is increasing, the biology underpinning locoregional sites remains unclear. Here, we applied targeted spatial proteomic approaches to primary and lymph node metastasis from an oropharyngeal SCC (OPSCC) cohort to understand the expression of proteins within tumors, and stromal compartments of the respective sites in samples of both matched and unmatched patients. In unmatched analyses of n = 43 primary and 11 nodal metastases, our data indicated that tumor cells in nodal metastases had higher levels of Ki-67, PARP, BAD, and cleaved caspase 9, suggesting a role for increased proliferation, DNA repair, and apoptosis within these metastatic cells. Conversely, in matched analyses (n = 7), pro-apoptotic markers BIM and BAD were enriched in the stroma of primary tumors. Univariate, overall survival (OS) analysis indicated CD25 in tumor regions of primary tumors to be associated with reduced survival (HR = 3.3, p = 0.003), while progesterone receptor (PR) was associated with an improved OS (HR = 0.33, p = 0.015). This study highlights the utility of spatial proteomics for delineating the tumor and stromal compartment composition, and utility toward understanding these properties in locoregional metastasis. These findings indicate unique biological properties of lymph node metastases that may elucidate further understanding of distant metastatic in OPSCC.

Introduction

Head and neck squamous cell carcinoma (HNSCC) is the 7th leading cause of cancer worldwide, with approximately 890,000 new cases and 450,000 deaths (1). HNSCC is considered a heterogeneous malignancy arising from the upper aerodigestive tract, particularly from the squamous mucosal line. The lip, oral and nasal cavity, paranasal sinuses, larynx, nasopharynx, oropharynx, and hypopharynx are the areas involved in HNSCC (1). Several risk factors contribute to the development of HNSCC, including both tobacco and alcohol consumption. Viral-driven HNSCC is also found, with Epstein–Barr Virus (EBV) and human papillomavirus (HPV) responsible for the nasopharynx and oropharynx malignancies, respectively (1).

Oropharyngeal SCC (OPSCC) is responsible for a quarter of HNSCC (2), with a tendency to occur more commonly in non-smokers and frequent nodal involvement. Treatment of OPSCC is aimed at curing and organ preservation using a multimodality approach. In locoregional disease, patients are often treated with chemoradiation. Most notably, HPV-induced OPSCC has a better prognosis compared to HPV-negative OPSCC, namely, better radiation sensitivity and overall survival. The 5-year survival rate for patients with HPV-negative and HPV-positive OPSCC is 46 and 57.4%, respectively (3).

In approximately two-thirds of OPSCC patients, locoregional metastasis has been reported. Advanced nodal status, especially extranodal extension, is a poor prognostic predictor, and defining the molecular phenotypes of these multiple sites is important to understand not only the primary tumor, but also the influence of lymph node metastasis to develop effective therapies (4). To gain a deeper understanding of the primary and metastatic tissues, studies of the tumor microenvironment (TME) are needed, ideally by comparing primary and lymph node metastasis, and where possible, where these are from matched patients. The identification and characterization of potential TME biomarkers could have significant predictive and prognostic value in the treatment selection for OPSCC patients (5). However, this can be challenging due to inter- or intra-tumoral heterogeneity (6). The cellular composition and molecular interactions between the TME and the host immune cells could play a role in disease progression and treatment resistance. It is thought that locoregional metastasis are immune-cold and therefore treatment-resistant to any immuno-modulatory treatment strategies (1, 7).

To delineate primary and nodal disease in OPSCC, we applied digital spatial profiling (DSP) of the TME using targeted-panel multiplex immunohistochemistry of tumor- and stromal-compartments. Our study found in unmatched analyses that the tumor compartments of the primary were enriched for ARG1, PD-L1, and nodal metastases were enriched for Ki-67, PARP, BAD, and cleaved caspase 9. Stromal compartments of the primary were enriched for VISTA and IDO1. In matched analyses, BIM and BAD were enriched in the stroma of primary tumors. Furthermore, protein signatures were identified to discriminate matched primary/nodal tissues, and survival associations were investigated compartmentally by tumor sample site.

Material and Methods

This study has the approval of the Queensland University of Technology Human Research Ethics Committee (UHREC #2000000494) and University of Queensland ratification. A tissue microarray (TMA) of OPSCC specimens was sourced from the Tristar Technologies Group (USA) and contained forty-three primary tumors, eleven nodal metastases, and seven matched pairs of primary tumor/nodal metastases with concordant clinicopathological annotations.

Nanostring GeoMx Digital Spatial Profiling (DSP)

The TMA slides from OPSCC samples were obtained and analyzed using the Nanostring GeoMX Digital Spatial Profiling (DSP) technology by the Systems Biology and Data Science Group at Griffith University (Gold Coast, Australia). Pan-cytokeratin and CD45 were the visualization markers used by the instrument to stain the tumor and lymphocytes, respectively. A protein panel of 68 antibodies was used, namely, immune activation, pan-tumor, immuno-oncology (IO) drug target, cell death, immune cell typing, human immune cell core panel, and PI3K/AKT panels. The slides were prepared according to the instructions of the manufacturer, and the tumor/stroma distinction was achieved by masking on PanCK⁺ or PanCK⁻ regions, respectively. Using the Nanostring nCounter^® platform, antibody barcodes were counted in accordance with the instructions of the manufacturer. In DSP analysis, external RNA Controls Consortium (ERCC) QC was employed to prepare the data for further bioinformatic analysis.

Bioinformatic Analysis

Data analysis was conducted in collaboration with the Queensland Cyber Infrastructure Foundation (QCIF, QLD, Australia). The quality of data was investigated using principal component analysis, and the suitability of the RUV-III normalization method was determined using coefficients of variation (8, 9). Differential analysis was carried out using Limma packages (10). Sparse partial least squares-discriminant analysis (sPLS-DA) within the mixOmics package was used to identify multivariate minimal protein signatures (11). The Kaplan–Meier survival analysis and Cox proportional hazards models were constructed within R studio (12) using Survival package (13) and plots generated by ggplot2 (14). Data shown are not adjusted for multiple testing. The false discovery rate (FDR) adjusted results were not significant within this cohort.

Results

OPSCC Patient Cohort

To investigate protein expression in primary tumors and nodal metastases, we evaluated tissue from a tissue microarray (TMA). The TMA included unmatched specimens from forty-three primary tumors and eleven nodal metastases, and seven matched primary and nodal metastases. All primary tumors were resected from oropharyngeal regions and had squamous cell carcinoma histology (Table 1). The TNM staging (8th edition) ranged from T2-4, N0-2, and M0 (Table 1). Only one patient within the cohort was HPV-16 positive.

Identification of Differentially Expressed Proteins by the Nanostring GeoMX DSP Assay

Nanostring Digital Spatial Profiler (DSP) was applied to investigate the protein expression of 68 TMA cores. By masking on PanCK⁺/PanCK⁻ regions, we compartmentalized protein expression within the tumor and stromal compartments (Figures 1A–D). Differential expression (DE) was performed on normalized protein counts to identify compartmental enrichment of proteins in unmatched primary tumors compared to nodal metastases and between patient-matched primary tumors and nodal metastases. We discovered significant differences in key deregulated proteins between the two comparison groups. Additionally, DE within the stroma between unmatched primary tumors and nodal metastasis revealed that V-domain IG suppressor of T cell activation (VISTA) and Indoleamine 2,3-dioxygenase 1 (IDO1) exhibited higher expression in primary tumors (Figures 2A, B). Parallel analysis of tumor compartments uncovered DE proteins. Programmed death-ligand 1 (PD-L1), Fibroblast activation protein-alpha (FAP-a), Poly (ADP-ribrose) polymerase (PARP), Ki-67, Progesterone receptor (PR), Arginase 1 (ARG1), CD56, BCL2-antagonist of cell death (BAD), and Cleaved Caspase 9. Accordingly, PD-L1, FAP-a, PR, and ARG1 had higher expression in primary tumors, while PARP, Ki-67, CD56, BAD, and Cleaved Caspase 9 were enriched in nodal metastases (Figures 2C, D). Furthermore, analysis of stromal compartments across matched primary and nodal metastasis specimens unveiled seven significant differentially expressed proteins, namely, smooth muscle actin (SMA), phosphatase and tensin homolog deleted on chromosome 10 (PTEN), CD163, Bcl-2-like protein 11 (BIM), BAD, PD-L1, and CD25. BIM, BAD, and CD25 were enriched in matched primary tumors, while SMA, PTEN, CD163, and PD-L1 were higher in matched nodal metastases (Figures 3A, B). Analysis of tumor compartments between matched primary and nodal metastasis specimens indicated no significant DE proteins.

Figure 1

Figure 2

Figure 3

Survival Associations of OPSCC Primary Tumor and Nodal Metastases

To assess the association between our protein expression and overall survival (OS), we performed a Cox proportional hazards model on all proteins. Analysis of the stromal compartment in primary tumor specimens revealed that the expression of NF1 (HR = 0.748, p = 0.025), CD27 (HR = 0.279, p = 0.035), and CD80 (HR = 0.703, p = 0.02) was associated with a better OS (Figure 4A). Moreover, tumoral compartment analysis across primary tumor samples indicated that PR (HR = 0.332, p = 0.015) was associated with better OS, while CD25 (HR = 3.311, p = 0.003) was associated with worse OS (Figure 4B). Survival analysis within the stromal compartment from nodal metastasis specimens showed that NY.ESO.1 (HR = 0.309, p = 0.016), and B7.H3 (HR = 0.249, p = 0.026) were associated with better OS (Figure 4C). Evaluation of tumoral compartment of nodal metastasis samples found several proteins associated with better OS, namely, SMA (HR = 0.685, p = 0.034), CD45 (HR = 0.463, p = 0.045), CD8 (HR = 0.251, p = 0.028), Fibronectin (HR = 0.472, p = 0.035), and STING (HR = 0.239, p = 0.024), however, BIM (HR = 1.979, p = 0.048), GZMA (HR = 4.332, p = 0.036), FOXP3 (HR = 3.258, p = 0.011), and PR (HR = 9.34, p = 0.017) were associated with worse OS (Figure 4D).

Figure 4

Multivariate Discrimination of OPSCC Primary Tumors From Nodal Metastases

Multivariate analysis by sparse partial least-squares discriminant analysis (sPLSDA) was employed to identify minimal protein signatures that collectively distinguish primary tumors from nodal metastases. Signatures within the stroma of matched primary vs nodal metastases stratified samples effectively (Figure 5A). The first signature (Figure 5B) included levels of SMA, PTEN, cleaved caspase 9, and CD25 (AUC = 0.979) (Figure 5C). A second signature comprised of CD95, CD80, and CD27 distinguished matched sample types as well (AUC = 0.918) (Figure 5D).

Figure 5

Discussion

HNSCC has a high risk of locoregional nodal metastasis, which affects patient prognosis and treatment outcomes (15). Patients with nodal metastasis are considered to have locoregionally advanced disease with a lower chance of remission (16). Currently, the predictors of nodal metastasis include tumor thickness and size, which have been shown to be unreliable predictors (17). It has been difficult to manage clinically negative neck nodes (N0) due to a lack of reliable predictors of occult metastasis (16, 18). In a study conducted by Shah et al., the authors found that there was a 40% chance of nodal metastasis in clinically node-negative neck dissections (19). For the purpose of distinguishing patients with a high risk of nodal metastasis, various pathological and clinical factors, namely, lymphovascular invasion, tumor differentiation, depth of invasion (DOI), and pattern of invasion (POI), have been reported (20, 21). Moreover, modern imaging modalities such as MRI, CT imaging, and PET/CT scanning, have been used to aid in the detection of locoregional nodal metastases. However, some radiographic features, namely, “subclinical”, or “microscopic”, or “occult” disease, remain difficult to diagnose using any of these approaches (22, 23). Therefore, companion diagnostics tools are needed to improve the prediction of the likelihood of the development of locoregional metastasis. Ideally, this would be possible by interrogating the primary tissue to determine its aggressiveness and propensity for metastasis.

Tumor tissue analysis by bulk expression or single cell RNA sequencing offers an overview of the molecular features of HNSCC tumors and their TME. These methods are incapable of revealing the spatial cellular properties required for the anti-tumor immune responses (24). Spatial proteomic approaches can provide compartment-specific tumor information to aid in delineating tumor composition. To garner insight into these properties that distinguish primary OPSCC tumors from their metastatic nodal counterparts, we have employed Digital Spatial Profiling to address a targeted profile of proteins and present this data as a first step in profiling HNSCC nodal involvement in an OPSCC cohort.

Unmatched analyses between primary tumors (n = 43) and nodal metastases (n = 11) provided insight into the potential dysregulation of several proteins, despite the inherent limitations associated with such a sampling strategy. VISTA, IDO1, and PD-L1 are key immune checkpoints, with VISTA and IDO1 appearing more abundant in the stromal compartment of primary tumors, whereas PD-L1 indicated higher expression in their respective tumor compartments. The V-domain Ig suppressor of T cell activation (VISTA) is an inhibitory immune checkpoint protein that is typically expressed on naïve CD4⁺ and Foxp3⁺ Tregs and functions by inhibiting T-cell proliferation and promoting naive to Treg conversion (25). VISTA was associated with several immune cell regions in the stroma but not in HNSCC tumors (26). Blockade of VISTA was found to boost anti-tumor immunity in the tumor microenvironment by increasing the number of activated dendritic cells (DCs) and decreasing the number of myeloid-derived suppressor cells (MDSCs). IDO1 induces T-cell apoptosis through activation of caspase 8 and releases mitochondrial cytochrome C, functioning in an immunosuppressive capacity (27). Programmed death-ligand 1 (PD-L1) is a canonical inhibitory immune checkpoint that binds PD-1 on the surface of tumors and immune cells (28). PD-L1 expression on the surface of HNSCC tumor cells is associated with a more robust anti-tumor immune response (29, 30).

Apoptotic pathways play an important role in tumorigenesis, and our results indicate that pro-apoptotic BAD and cleaved caspase 9 were enriched in tumor cells of unmatched nodal metastases. Conversely, in the matched analysis, BIM and BAD were enriched in the stroma of primary tumors. Under various physiological and patho-physiological conditions, Bcl-2 interacting mediator of cell death (BIM) promotes the intrinsic apoptotic pathway (31). Bcl2-associated agonist of cell death (BAD) is a member of the BCL2 family of proteins that act as pro-apoptotic regulators (32). The expression of BAD has been linked to chemoresistance in cancer patients (33, 34). Caspase 9 functions as a pro-apoptotic regulator, allowing the activation of effector caspases 3 and 7 (35). Caspase 9-induced apoptosis has been linked to chemotherapy response. Studies have shown that HNSCC tumors may be resistant to cisplatin if they have a reduced expression of caspase 9 (36).

It is interesting to note that despite an imbalance in samples in the unmatched analysis, several proteins appear enriched within nodal metastases relative to primary tumors. In addition to the pro-apoptotic markers above, Ki-67 and PARP appear enriched in nodal tumor cells. Ki-67 is an established proliferation marker (37), while PARP responds to DNA damage by recruiting effector proteins to repair single-strand breaks (38, 39). PARP inhibitors have been studied as a promising drug to overcome the limitations of conventional therapies that cause DNA damage, such as chemotherapy or radiotherapy (40). This pro-apoptotic, proliferative, and DNA damage phenotype of unmatched nodal metastatic cells is a novel finding in our data that requires further validation, perhaps indicating tumor evolution or response to changes in the cellular ecosystem of the lymph node.

Several other notable features of our data include increased expression of CD25 in the stroma of matched primary tumors. CD25, also known as the IL-2 receptor alpha, is a protein found on activated T cells, specifically Tregs (41). Interestingly, we found that it was CD25 expression within primary tumor regions, not stroma, that was associated with poorer OS. Additionally, PR expression appeared associated with better survival within primary tumors and was also enriched within their tumor regions relative to nodal metastases.

In addition to the differential expression of each individual protein, we applied a multivariate statistical model (sPLSDA) to further discern features that collectively discriminated between our matched patient samples. Of note, this model only performed effectively in stratifying these sample types by their stroma. Expression of CD25 and cleaved caspase 9 in the primary samples and PTEN and SMA in nodal samples could discern samples. Similarly, levels of CD80 and CD27 within primary samples and CD95 in nodal samples could separate these samples. This method offers an alternative to traditional differential expression that may provide insight into contributing differences in observed phenotypes using a multivariate approach.

A protein association with overall survival (OS) was investigated using the Univariate Cox proportional hazards model. Stromal expression of CD27 in primary tumor specimens, and NY-ESO.1 and B7-H3 expression in the nodal metastasis samples, was associated with improved OS. In nodal metastases tumor regions, CD8 and STING were associated with improved OS; however, GZMA was associated with poorer OS. Interestingly, the expression of PR in the tumoral compartment of primary and nodal metastasis tumors demonstrated a different survival pattern. Although the PR expression was associated with improved OS in primary tumors, it was associated with worse OS in nodal metastasis specimens. Immune response protein markers, CD27, and NY.ESO.1 within the stromal compartment of primary and nodal metastasis specimens, were associated with improved OS in our study. Cluster of differentiation 27 (CD27) belongs to the tumor necrosis factor (TNF) receptor superfamily and is involved in T and B cell co-stimulation (42). New York esophageal squamous cell carcinoma 1 (NY.ESO.1) is a member of the cancer testis antigen (CTA) family, which regulates both humoral and cellular immune responses. NY-ESO.1 expression has been linked to higher tumor differentiation grade and stage, and lymph node metastasis (43). In our study, we found that B7.H3 and STING protein expression, NF-κB pathway markers, were linked to improved OS. B7.H3, also known as CD276, promotes anti-tumor immune response by activating T and NK cells (44). STING, on the other hand, contributes to the immune response to tumor cells through the upregulation of interferon gamma 1 (IFN1) (45). Progesterone receptor (PR) is a type of androgen receptor and a member of the nuclear receptor family of transcription factors that regulates target gene expression networks in response to its ligand progesterone (46).

Our study has identified tumor and stromal compartment-specific proteins and signatures that may have predictive and prognostic implications for HNSCC and the development of nodal metastasis. Nevertheless, the study is impacted by the number of samples for each cohort, in particular the matched group. We propose further investigation to profile primary, locoregional, and distant metastasis from matched patient samples to understand the molecular features driving the development of metastasis in OPSCC.

Funding

This project is funded by the Garnett Passe and Rodney Williams Memorial Foundation (GPRWMF) conjoint grant for AK and BH. AK is supported by a fellowship from the NHMRC (1157741). KO’B and AK are supported by the PA Research Foundation.

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.

References

• 1

  JohnsonDEBurtnessBLeemansCRLuiVWYBaumanJEGrandisJR. Head and Neck Squamous Cell Carcinoma. Nat Rev Dis Primers (2020) 6(1):92. doi: 10.1038/s41572-020-00224-3

  • CrossRef
  • Google Scholar

• 2

  SabatiniMEChioccaS. Human Papillomavirus as a Driver of Head and Neck Cancers. Br J Cancer (2020) 122(3):306–14. doi: 10.1038/s41416-019-0602-7

  • CrossRef
  • Google Scholar

• 3

  ChenT-CWuC-TKoJ-YYangT-LLouP-JWangC-Pet al. Clinical Characteristics and Treatment Outcome of Oropharyngeal Squamous Cell Carcinoma in an Endemic Betel Quid Region. Sci Rep (2020) 10(1):526. doi: 10.1038/s41598-019-57177-1

  • CrossRef
  • Google Scholar

• 4

  ArgirisAKaramouzisMVRabenDFerrisRL. Head and Neck Cancer. Lancet (2008) 371(9625):1695–709. doi: 10.1016/S0140-6736(08)60728-X

  • CrossRef
  • Google Scholar

• 5

  Sadeghi RadHMonkmanJWarkianiMELadwaRO'ByrneKRezaeiNet al. Understanding the Tumor Microenvironment for Effective Immunotherapy. Medicinal Res Rev (2021) 41(3):1474–98. doi: 10.1002/med.21765

  • CrossRef
  • Google Scholar

• 6

  SternPLDalianisT. Oropharyngeal Squamous Cell Carcinoma Treatment in the Era of Immune Checkpoint Inhibitors. Viruses (2021) 13(7):1234. doi: 10.3390/v13071234

  • CrossRef
  • Google Scholar

• 7

  KulasingheATaheriTO’ByrneKHughesBGMKennyLPunyadeeraC. Highly Multiplexed Digital Spatial Profiling of the Tumor Microenvironment of Head and Neck Squamous Cell Carcinoma Patients. Front Oncol (2021) 10. doi: 10.3389/fonc.2020.607349

  • CrossRef
  • Google Scholar

• 8

  RissoDNgaiJSpeedTPDudoitS. Normalization of RNA-Seq Data Using Factor Analysis of Control Genes or Samples. Nat Biotechnol (2014) 32(9):896–902. doi: 10.1038/nbt.2931

  • CrossRef
  • Google Scholar

• 9

  MolaniaRGagnon-BartschJADobrovicASpeedTP. A New Normalization for Nanostring Ncounter Gene Expression Data. Nucleic Acids Res (2019) 47(12):6073–83. doi: 10.1093/nar/gkz433

  • CrossRef
  • Google Scholar

• 10

  RitchieMEPhipsonBWuDHuYLawCWShiWet al. Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res (2015) 43(7):e47. doi: 10.1093/nar/gkv007

  • CrossRef
  • Google Scholar

• 11

  RohartFGautierBSinghAKALC. Mixomics: An R Package for 'Omics Feature Selection and Multiple Data Integration. PloS Comput Biol (2017) 13(11):e1005752. doi: 10.1371/journal.pcbi.1005752

  • CrossRef
  • Google Scholar

• 12

  Team R. RStudio: Integrated Development for R. Boston, MA: RStudio (2020). Available at: http://www.rstudio.com/.

  • Google Scholar

• 13

  TherneauTM. A Package for Survival Analysis in R (2021). Available at: https://CRAN.R-project.org/package=survival.

  • Google Scholar

• 14

  HadleyW. Ggplot2: Elegrant Graphics for Data Analysis. Springer (2016).

  • Google Scholar

• 15

  ChowLQM. Head and Neck Cancer. N Engl J Med (2020) 382(1):60–72. doi: 10.1056/NEJMra1715715

  • CrossRef
  • Google Scholar

• 16

  AmpilFLNathanCASangsterGCalditoG. Head and Neck Cancer With Lower Neck Nodal Metastases: Management of 23 Cases and Review of the Literature. Oral Oncol (2012) 48(4):325–8. doi: 10.1016/j.oraloncology.2011.11.016

  • CrossRef
  • Google Scholar

• 17

  PuramSVTiroshIParikhASPatelAPYizhakKGillespieSet al. Single-Cell Transcriptomic Analysis of Primary and Metastatic Tumor Ecosystems in Head and Neck Cancer. Cell (2017) 171(7):1611–24.e24. doi: 10.1016/j.cell.2017.10.044

  • CrossRef
  • Google Scholar

• 18

  De SilvaRKSiriwardenaBSamaranayakaAAbeyasingheW. Tilakaratne WM. A Model to Predict Nodal Metastasis in Patients With Oral Squamous Cell Carcinoma. PloS One (2018) 13(8):e0201755. doi: 10.1371/journal.pone.0201755

  • CrossRef
  • Google Scholar

• 19

  ShahJP. Patterns of Cervical Lymph Node Metastasis From Squamous Carcinomas of the Upper Aerodigestive Tract. Am J Surg (1990) 160(4):405–9. doi: 10.1016/S0002-9610(05)80554-9

  • CrossRef
  • Google Scholar

• 20

  TaiSKLiWYYangMHChuPYWangYF. Perineural Invasion in T1 Oral Squamous Cell Carcinoma Indicates the Need for Aggressive Elective Neck Dissection. Am J Surg Pathol (2013) 37(8):1164–72. doi: 10.1097/PAS.0b013e318285f684

  • CrossRef
  • Google Scholar

• 21

  RajapaksheRMPallegamaRWJayasooriyaPRSiriwardenaBSAttygallaAMHewapathiranaSet al. A Retrospective Analysis to Determine Factors Contributing to the Survival of Patients With Oral Squamous Cell Carcinoma. Cancer Epidemiol (2015) 39(3):360–6. doi: 10.1016/j.canep.2015.02.011

  • CrossRef
  • Google Scholar

• 22

  SohnBKohYWKangWJLeeJHShinNYKimJ. Is There an Additive Value of 18 F-FDG PET-CT to CT/MRI for Detecting Nodal Metastasis in Oropharyngeal Squamous Cell Carcinoma Patients With Palpably Negative Neck? Acta Radiol (2016) 57(11):1352–9. doi: 10.1177/0284185115587544

  • CrossRef
  • Google Scholar

• 23

  KannBHAnejaSLoganadaneGVKellyJRSmithSMDeckerRHet al. Pretreatment Identification of Head and Neck Cancer Nodal Metastasis and Extranodal Extension Using Deep Learning Neural Networks. Sci Rep (2018) 8(1):14036. doi: 10.1038/s41598-018-32441-y

  • CrossRef
  • Google Scholar

• 24

  Sadeghi RadHBazazSRMonkmanJEbrahimi WarkianiMRezaeiNO'ByrneKet al. The Evolving Landscape of Predictive Biomarkers in Immuno-Oncology With a Focus on Spatial Technologies. Clin Trans Immunol (2020) 9(11):e1215. doi: 10.1002/cti2.1215

  • CrossRef
  • Google Scholar

• 25

  ElTanboulyMACroteauWNoelleRJLinesJL. VISTA: A Novel Immunotherapy Target for Normalizing Innate and Adaptive Immunity. Semin Immunol (2019) 42:101308. doi: 10.1016/j.smim.2019.101308

  • CrossRef
  • Google Scholar

• 26

  LinesJLPantaziEMakJSempereLFWangLO'ConnellSet al. VISTA is an Immune Checkpoint Molecule for Human T Cells. Cancer Res (2014) 74(7):1924–32. doi: 10.1158/0008-5472.CAN-13-1504

  • CrossRef
  • Google Scholar

• 27

  FallarinoFGrohmannUVaccaCBianchiROrabonaCSprecaAet al. T Cell Apoptosis by Tryptophan Catabolism. Cell Death Differ (2002) 9(10):1069–77. doi: 10.1038/sj.cdd.4401073

  • CrossRef
  • Google Scholar

• 28

  LieblMCHofmannTG. Identification of Responders to Immune Checkpoint Therapy: Which Biomarkers Have the Highest Value? J Eur Acad Dermatol Venereol (2019) 33:52–6. doi: 10.1111/jdv.15992

  • CrossRef
  • Google Scholar

• 29

  CohenEESoulièresDLe TourneauCDinisJLicitraLAhnM-Jet al. Pembrolizumab Versus Methotrexate, Docetaxel, or Cetuximab for Recurrent or Metastatic Head-and-Neck Squamous Cell Carcinoma (KEYNOTE-040): A Randomised, Open-Label, Phase 3 Study. Lancet (2019) 393(10167):156–67. doi: 10.1016/S0140-6736(18)31999-8

  • CrossRef
  • Google Scholar

• 30

  BrahmerJRTykodiSSChowLQHwuW-JTopalianSLHwuPet al. Safety and Activity of Anti–PD-L1 Antibody in Patients With Advanced Cancer. N Engl J Med (2012) 366(26):2455–65. doi: 10.1056/NEJMoa1200694

  • CrossRef
  • Google Scholar

• 31

  ShuklaSSaxenaSSinghBKKakkarP. BH3-Only Protein BIM: An Emerging Target in Chemotherapy. Eur J Cell Biol (2017) 96(8):728–38. doi: 10.1016/j.ejcb.2017.09.002

  • CrossRef
  • Google Scholar

• 32

  BoacBMAbbasiFIsmail-KhanRXiongYSiddiqueAParkHet al. Expression of the BAD Pathway is a Marker of Triple-Negative Status and Poor Outcome. Sci Rep (2019) 9(1):1–14. doi: 10.1038/s41598-019-53695-0

  • CrossRef
  • Google Scholar

• 33

  BansalNMarchionDCBicakuEXiongYChenNSticklesXBet al. BCL2 Antagonist of Cell Death Kinases, Phosphatases, and Ovarian Cancer Sensitivity to Cisplatin. J Gynecol Oncol (2012) 23(1):35–42. doi: 10.3802/jgo.2012.23.1.35

  • CrossRef
  • Google Scholar

• 34

  MarchionDCCottrillHMXiongYChenNBicakuEFulpWJet al. BAD Phosphorylation Determines Ovarian Cancer Chemosensitivity and Patient Survival. Clin Cancer Res (2011) 17(19):6356–66. doi: 10.1158/1078-0432.CCR-11-0735

  • CrossRef
  • Google Scholar

• 35

  ThornberryNALazebnikY. Caspases: Enemies Within. Science (1998) 281(5381):1312–6. doi: 10.1126/science.281.5381.1312

  • CrossRef
  • Google Scholar

• 36

  KuwaharaDTsutsumiKOyakeDOhtaTNishikawaHKoizukaI. Inhibition of Caspase-9 Activity and Apaf-1 Expression in Cisplatin-Resistant Head and Neck Squamous Cell Carcinoma Cells. Auris Nasus Larynx (2003) 30:85–8. doi: 10.1016/S0385-8146(02)00129-3

  • CrossRef
  • Google Scholar

• 37

  SunXKaufmanPD. Ki-67: More Than a Proliferation Marker. Chromosoma (2018) 127(2):175–86. doi: 10.1007/s00412-018-0659-8

  • CrossRef
  • Google Scholar

• 38

  LangelierM-FRiccioAAPascalJM. PARP-2 and PARP-3 are Selectively Activated by 5′ Phosphorylated DNA Breaks Through an Allosteric Regulatory Mechanism Shared With PARP-1. Nucleic Acids Res (2014) 42(12):7762–75. doi: 10.1093/nar/gku474

  • CrossRef
  • Google Scholar

• 39

  EustermannSWuW-FLangelierM-FYangJ-CEastonLERiccioAAet al. Structural Basis of Detection and Signaling of DNA Single-Strand Breaks by Human PARP-1. Mol Cell (2015) 60(5):742–54. doi: 10.1016/j.molcel.2015.10.032

  • CrossRef
  • Google Scholar

• 40

  LordCJAshworthA. PARP Inhibitors: Synthetic Lethality in the Clinic. Science (2017) 355(6330):1152–8. doi: 10.1126/science.aam7344

  • CrossRef
  • Google Scholar

• 41

  FlynnMJHartleyJA. The Emerging Role of Anti-CD 25 Directed Therapies as Both Immune Modulators and Targeted Agents in Cancer. Br J Haematol (2017) 179(1):20–35. doi: 10.1111/bjh.14770

  • CrossRef
  • Google Scholar

• 42

  StarzerAMBerghoffAS. New Emerging Targets in Cancer Immunotherapy: CD27 (Tnfrsf7). ESMO Open (2019) 4:e000629. doi: 10.1136/esmoopen-2019-000629

  • CrossRef
  • Google Scholar

• 43

  ThomasRAl-KhadairiGRoelandsJHendrickxWDermimeSBedognettiDet al. NY-ESO-1 Based Immunotherapy of Cancer: Current Perspectives. Front Immunol (2018) 9:947. doi: 10.3389/fimmu.2018.00947

  • CrossRef
  • Google Scholar

• 44

  LuoLChapovalAIFliesDBZhuGHiranoFWangSet al. B7-H3 Enhances Tumor Immunity In Vivo by Costimulating Rapid Clonal Expansion of Antigen-Specific CD8+ Cytolytic T Cells. J Immunol (2004) 173(9):5445–50. doi: 10.4049/jimmunol.173.9.5445

  • CrossRef
  • Google Scholar

• 45

  IshikawaHBarberGN. STING Is an Endoplasmic Reticulum Adaptor That Facilitates Innate Immune Signalling. Nature (2008) 455(7213):674–8. doi: 10.1038/nature07317

  • CrossRef
  • Google Scholar

• 46

  GrimmSLHartigSMEdwardsDP. Progesterone Receptor Signaling Mechanisms. J Mol Biol (2016) 428(19):3831–49. doi: 10.1016/j.jmb.2016.06.020

  • CrossRef
  • Google Scholar