Single cell analysis in head and neck cancer reveals potential immune evasion mechanisms during early metastasis.

Profiling tumors at single-cell resolution provides an opportunity to understand complexities underpinning lymph-node metastases in head and neck squamous-cell carcinoma. Single-cell RNAseq (scRNAseq) analysis of cancer-cell trajectories identifies a subpopulation of pre-metastatic cells, driven by actionable pathways including AXL and AURK. Blocking these two proteins blunts tumor invasion in patient-derived cultures. Furthermore, scRNAseq analyses of tumor-infiltrating CD8 + T-lymphocytes show two distinct trajectories to T-cell dysfunction, corroborated by their clonal architecture based on single-cell T-cell receptor sequencing. By determining key modulators of these trajectories, followed by validation using external datasets and functional experiments, we uncover a role for SOX4 in mediating T-cell exhaustion. Finally, interactome analyses between pre-metastatic tumor cells and CD8 + T-lymphocytes uncover a putative role for the Midkine pathway in immune-modulation and this is confirmed by scRNAseq of tumors from humanized mice. Aside from specific findings, this study demonstrates the importance of tumor heterogeneity analyses in identifying key vulnerabilities during early metastasis.

A chemical genetic screen identifies Aurora kinases as a therapeutic target in EGFR T790M negative, gefitinib-resistant head and neck squamous cell carcinoma (HNSCC).

BACKGROUND: Overexpression of epidermal growth factor receptor (EGFR), and downstream pathway activation appears to be a common oncogenic driver in the majority of head and neck squamous cell cancers (HNSCCs); yet targeting EGFR for the treatment of HNSCC has met with limited success. Apart from the anti-EGFR antibody cetuximab, no small molecule EGFR/tyrosine kinase inhibitors (TKIs) have progressed to routine clinical use. The aim of this study was to determine factors contributing to the lack of response to TKIs and identify alternative therapeutic vulnerabilities. METHODS: Genomic and transcriptomic sequencing, high-throughput compound screens, overexpression and siRNA knockdown, western blot, in vivo xenograft studies. FINDINGS: We derived three pairs of isogenic gefitinib (TKI)-sensitive and resistant patient-derived HNSCC cell lines. Genomic sequencing of gefitinib-resistant cell lines identified a lack of activating and resistance-associated EGFR mutations. Instead, transcriptomic sequencing showed upregulated EMT gene signature in the gefitinib-resistant cells with a corresponding increase in their migratory phenotype. Additionally, the resistant cell displayed reduced growth rate. Surprisingly, while gefitinib-resistant cells were independent of EGFR for survival, they nonetheless displayed activation of downstream ERK and AKT signalling. High-throughput screening (HTS) of druggable, small molecule libraries revealed that the gefitinib-resistant cells were particularly sensitive to inhibitors of genes involved in cell cycle and mitosis, such as Aurora kinase inhibitors (AKIs), cyclin-dependent kinase (CDK) inhibitors, and microtubule inhibitors. Notably our results showed that in the EGFR inhibited state, Aurora kinases are essential for cell survival. INTERPRETATION: Our study demonstrates that in the absence of activating EGFR mutations, HNSCCs may gain resistance to gefitinib through decreased cell proliferation, which makes them exceptionally vulnerable to cell-cycle inhibitors. FUNDING: Agency for Science, Technology, and Research (A*STAR), National Medical Research Council (NMRC), and the National Institutes of Health (NIH)/National Cancer Institute (NCI).

Pairing a prognostic target with potential therapeutic strategy for head and neck cancer.

OBJECTIVES: We have previously identified and validated a panel of molecular prognostic markers (ATP13A3, SSR3, and ANO1) for Head and Neck Squamous Cell Carcinoma (HNSCC). The aim of this study was to investigate the consequence of ATP13A3 dysregulation on signaling pathways, to aid in formulating a therapeutic strategy targeting ATP13A3-overexpressing HNSCC. MATERIALS AND METHODS: Gene Set Enrichment Analysis (GSEA) was performed on HNSCC microarray expression data (Internal local dataset [n = 92], TCGA [n = 232], EMBL [n = 81]) to identify pathways associated with high expression of ATP13A3. Validation was performed using immunohistochemistry (IHC) on tissue microarrays (TMAs) of head and neck cancers (n = 333), staining for ATP13A3 and phosphorylated Aurora kinase A (phospho-T288). Short interfering RNA was used to knockdown ATP13A3 expression in patient derived HNSCC cell lines. Protein expression of ATP13A3 and Aurora kinase A was then assessed by immunoblotting. RESULTS: GSEA identified Aurora kinase pathway to be associated with high expression of ATP13A3 (p = 0.026). The Aurora kinase pathway was also associated with a trend towards poor prognosis and tumor aggressiveness (p = 0.086, 0.094, respectively). Furthermore, the immunohistochemical staining results revealed a significant association between Aurora kinase activity and high ATP13A3 expression (p < 0.001). Knockdown of ATP13A3 in human head and neck cell lines showed decrease in Aurora kinase A levels. CONCLUSION: Tumors with high ATP13A3 are associated with high Aurora kinase activity. This suggests a potential therapeutic role of Aurora kinase inhibitors in a subset of poor prognosis HNSCC patients with overexpression of ATP13A3.

Characterization of colibactin-associated mutational signature in an Asian oral squamous cell carcinoma and in other mucosal tumor types.

Mutational signatures can reveal the history of mutagenic processes that cells were exposed to before and during tumorigenesis. We expect that as-yet-undiscovered mutational processes will shed further light on mutagenesis leading to carcinogenesis. With this in mind, we analyzed the mutational spectra of 36 Asian oral squamous cell carcinomas. The mutational spectra of two samples from patients who presented with oral bacterial infections showed novel mutational signatures. One of these novel signatures, SBS_AnT, is characterized by a preponderance of thymine mutations, strong transcriptional strand bias, and enrichment for adenines in the 4 bp 5' of mutation sites. The mutational signature described in this manuscript was shown to be caused by colibactin, a bacterial mutagen produced by E. coli carrying the pks-island. Examination of publicly available sequencing data revealed SBS_AnT in 25 tumors from several mucosal tissue types, expanding the list of tissues in which this mutational signature is observed.

Longitudinal single-cell RNA sequencing of patient-derived primary cells reveals drug-induced infidelity in stem cell hierarchy.

Chemo-resistance is one of the major causes of cancer-related deaths. Here we used single-cell transcriptomics to investigate divergent modes of chemo-resistance in tumor cells. We observed that higher degree of phenotypic intra-tumor heterogeneity (ITH) favors selection of pre-existing drug-resistant cells, whereas phenotypically homogeneous cells engage covert epigenetic mechanisms to trans-differentiate under drug-selection. This adaptation was driven by selection-induced gain of H3K27ac marks on bivalently poised resistance-associated chromatin, and therefore not expressed in the treatment-naïve setting. Mechanistic interrogation of this phenomenon revealed that drug-induced adaptation was acquired upon the loss of stem factor SOX2, and a concomitant gain of SOX9. Strikingly we observed an enrichment of SOX9 at drug-induced H3K27ac sites, suggesting that tumor evolution could be driven by stem cell-switch-mediated epigenetic plasticity. Importantly, JQ1 mediated inhibition of BRD4 could reverse drug-induced adaptation. These results provide mechanistic insights into the modes of therapy-induced cellular plasticity and underscore the use of epigenetic inhibitors in targeting tumor evolution.

GNA13 expression promotes drug resistance and tumor-initiating phenotypes in squamous cell cancers.

Treatment failure in solid tumors occurs due to the survival of specific subpopulations of cells that possess tumor-initiating (TIC) phenotypes. Studies have implicated G protein-coupled-receptors (GPCRs) in cancer progression and the acquisition of TIC phenotypes. Many of the implicated GPCRs signal through the G protein GNA13. In this study, we demonstrate that GNA13 is upregulated in many solid tumors and impacts survival and metastases in patients. GNA13 levels modulate drug resistance and TIC-like phenotypes in patient-derived head and neck squamous cell carcinoma (HNSCC) cells in vitro and in vivo. Blockade of GNA13 expression, or of select downstream pathways, using small-molecule inhibitors abrogates GNA13-induced TIC phenotypes, rendering cells vulnerable to standard-of-care cytotoxic therapies. Taken together, these data indicate that GNA13 expression is a potential prognostic biomarker for tumor progression, and that interfering with GNA13-induced signaling provides a novel strategy to block TICs and drug resistance in HNSCCs.

A three gene immunohistochemical panel serves as an adjunct to clinical staging of patients with head and neck cancer.

BACKGROUND: Current management of head and neck squamous cell carcinoma (HNSCC) depends on tumor staging. Despite refinements in clinical staging algorithms, outcomes remain unchanged for the last two decades. In this study, we set out to identify a small, clinically applicable molecular panel to aid prognostication of patients with HNSCC. MATERIALS AND METHODS: Data from The Cancer Genome Atlas (TCGA) was used to derive copy number aberrations and expression changes to identify putative prognostic genes. To account for cross entity relevance of the biomarkers, HNSCC (n = 276), breast (n = 808) and lung cancer (n = 282) datasets were used to identify robust and reproducible markers with prognostic potential. Validation was performed using immunohistochemistry (IHC) on tissue microarrays of an independent cohort of HNSCC (n = 333). FINDINGS: Using GISTIC algorithm together with gene expression analysis, we identified six putative prognostic genes in at least two out of three cancers analyzed, of which four were successfully optimized for automated IHC. Of these, three were successfully validated; each molecular target being significantly prognostic on univariate analysis. Patients were differentially segregated into four prognostic groups based on the number of genes dysregulated (p < 0.001). The IHC panel remained an independent predictor of survival after adjusting for known survival covariates including clinical staging criteria in a multivariate Cox regression model (p < 0.001). . INTERPRETATION: We have identified and validated a clinically applicable IHC biomarker panel that is independently associated with overall survival. This panel is readily applicable, serving as a useful adjunct to current staging systems and provides novel targets for future therapeutic strategies.

Long noncoding RNA EGFR-AS1 mediates epidermal growth factor receptor addiction and modulates treatment response in squamous cell carcinoma.

Targeting EGFR is a validated approach in the treatment of squamous-cell cancers (SCCs), although there are no established biomarkers for predicting response. We have identified a synonymous mutation in EGFR, c.2361G>A (encoding p.Gln787Gln), in two patients with head and neck SCC (HNSCC) who were exceptional responders to gefitinib, and we showed in patient-derived cultures that the A/A genotype was associated with greater sensitivity to tyrosine kinase inhibitors (TKIs) as compared to the G/A and G/G genotypes. Remarkably, single-copy G>A nucleotide editing in isogenic models conferred a 70-fold increase in sensitivity due to decreased stability of the EGFR-AS1 long noncoding RNA (lncRNA). In the appropriate context, sensitivity could be recapitulated through EGFR-AS1 knockdown in vitro and in vivo, whereas overexpression was sufficient to induce resistance to TKIs. Reduced EGFR-AS1 levels shifted splicing toward EGFR isoform D, leading to ligand-mediated pathway activation. In co-clinical trials involving patients and patient-derived xenograft (PDX) models, tumor shrinkage was most pronounced in the context of the A/A genotype for EGFR-Q787Q, low expression of EGFR-AS1 and high expression of EGFR isoform D. Our study reveals how a 'silent' mutation influences the levels of a lncRNA, resulting in noncanonical EGFR addiction, and delineates a new predictive biomarker suite for response to EGFR TKIs.

Phenotype-driven precision oncology as a guide for clinical decisions one patient at a time.

Genomics-driven cancer therapeutics has gained prominence in personalized cancer treatment. However, its utility in indications lacking biomarker-driven treatment strategies remains limited. Here we present a "phenotype-driven precision-oncology" approach, based on the notion that biological response to perturbations, chemical or genetic, in ex vivo patient-individualized models can serve as predictive biomarkers for therapeutic response in the clinic. We generated a library of "screenable" patient-derived primary cultures (PDCs) for head and neck squamous cell carcinomas that reproducibly predicted treatment response in matched patient-derived-xenograft models. Importantly, PDCs could guide clinical practice and predict tumour progression in two n = 1 co-clinical trials. Comprehensive "-omics" interrogation of PDCs derived from one of these models revealed YAP1 as a putative biomarker for treatment response and survival in ~24% of oral squamous cell carcinoma. We envision that scaling of the proposed PDC approach could uncover biomarkers for therapeutic stratification and guide real-time therapeutic decisions in the future.Treatment response in patient-derived models may serve as a biomarker for response in the clinic. Here, the authors use paired patient-derived mouse xenografts and patient-derived primary culture models from head and neck squamous cell carcinomas, including metastasis, as models for high-throughput screening of anti-cancer drugs.

Serglycin expression: An independent marker of distant metastases in nasopharyngeal carcinoma.

BACKGROUND: Nasopharyngeal carcinoma (NPC) has a high propensity for metastasis. The purpose of this study was for us to determine whether serglycin expression can be used to predict distant metastases. METHODS: Serglycin expression of tumor tissue of 112 patients with NPC was assessed based on percentage of tumor cells expressing serglycin, staining intensity, percentage of tumor-infiltrated lymphocyte (TIL) expressing serglycin and TIL-staining intensity. RESULTS: Risk factors for distant metastases include sex, smoking status, tumor intensity, and TIL percentage for serglycin. The odds of distant metastases was 4.13 and 0.18 in patients with strong tumor intensity and >50% TIL percentage, respectively. Based on a nomogram incorporating predictors, patients were stratified into 2 prognostic groups. The proportion of distant metastases in the high-risk group (strong tumor intensity and ≤50% TIL percentage) was 78% versus 19% in the low risk group (p < .001). CONCLUSION: Patients with NPC with tumors showing strong tumor intensity and low TIL percentage with serglycin may be at high risk for distant metastases.

Mutational landscapes of tongue carcinoma reveal recurrent mutations in genes of therapeutic and prognostic relevance.

BACKGROUND: Carcinoma of the oral tongue (OTSCC) is the most common malignancy of the oral cavity, characterized by frequent recurrence and poor survival. The last three decades has witnessed a change in the OTSCC epidemiological profile, with increasing incidence in younger patients, females and never-smokers. Here, we sought to characterize the OTSCC genomic landscape and to determine factors that may delineate the genetic basis of this disease, inform prognosis and identify targets for therapeutic intervention. METHODS: Seventy-eight cases were subjected to whole-exome (n = 18) and targeted deep sequencing (n = 60). RESULTS: While the most common mutation was in TP53, the OTSCC genetic landscape differed from previously described cohorts of patients with head and neck tumors: OTSCCs demonstrated frequent mutations in DST and RNF213, while alterations in CDKN2A and NOTCH1 were significantly less frequent. Despite a lack of previously reported NOTCH1 mutations, integrated analysis showed enrichments of alterations affecting Notch signaling in OTSCC. Importantly, these Notch pathway alterations were prognostic on multivariate analyses. A high proportion of OTSCCs also presented with alterations in drug targetable and chromatin remodeling genes. Patients harboring mutations in actionable pathways were more likely to succumb from recurrent disease compared with those who did not, suggesting that the former should be considered for treatment with targeted compounds in future trials. CONCLUSIONS: Our study defines the Asian OTSCC mutational landscape, highlighting the key role of Notch signaling in oral tongue tumorigenesis. We also observed somatic mutations in multiple therapeutically relevant genes, which may represent candidate drug targets in this highly lethal tumor type.

An eleven gene molecular signature for extra-capsular spread in oral squamous cell carcinoma serves as a prognosticator of outcome in patients without nodal metastases.

OBJECTIVES: Extracapsular spread (ECS) is an important prognostic factor for oral squamous cell carcinoma (OSCC) and is used to guide management. In this study, we aimed to identify an expression profile signature for ECS in node-positive OSCC using data derived from two different sources: a cohort of OSCC patients from our institution (National Cancer Centre Singapore) and The Cancer Genome Atlas (TCGA) head and neck squamous cell carcinoma (HNSCC) cohort. We also sought to determine if this signature could serve as a prognostic factor in node negative cancers. MATERIALS AND METHODS: Patients with a histological diagnosis of OSCC were identified from an institutional database and fresh tumor samples were retrieved. RNA was extracted and gene expression profiling was performed using the Affymetrix GeneChip Human Genome U133 Plus 2.0 microarray platform. RNA sequence data and corresponding clinical data for the TCGA HNSCC cohort were downloaded from the TCGA Data Portal. All data analyses were conducted using R package and SPSS. RESULTS: We identified an 11 gene signature (GGH, MTFR1, CDKN3, PSRC1, SMIM3, CA9, IRX4, CPA3, ZSCAN16, CBX7 and ZFP3) which was robust in segregating tumors by ECS status. In node negative patients, patients harboring this ECS signature had a significantly worse overall survival (p=0.04). CONCLUSIONS: An eleven gene signature for ECS was derived. Our results also suggest that this signature is prognostic in a separate subset of patients with no nodal metastasis Further validation of this signature on other datasets and immunohistochemical studies are required to establish utility of this signature in stratifying early stage OSCC patients.

Tongue carcinoma infrequently harbor common actionable genetic alterations.

BACKGROUND: Oral tongue squamous cell carcinomas (TSCC) are a unique subset of head and neck cancers with a distinct demographic profile, where up to half of the cases are never smokers. A small proportion of patients with OSCC are known to respond to EGFR TKI. We used a high-sensitivity mass spectrometry-based mutation profiling platform to determine the EGFR mutation status, as well as other actionable alterations in a series of Asian TSCC. METHODS: 66 TSCC patients treated between 1998-2009 with complete clinico-pathologic data were included in this study. Somatic mutation profiling was performed using Sequenom LungCarta v1.0, and correlated with clinical parameters. RESULTS: Mutations were identified in 20/66(30.3%) of samples and involved TP53, STK11, MET, PIK3CA, BRAF and NRF2. No activating EGFR mutations or KRAS mutations were discovered in our series, where just over a third were never smokers. The most common mutations were in p53 (10.6%; n = 7) and MET (10.6%, n = 11) followed by STK11 (9.1%, n = 6) and PIK3CA (4.5%, n = 3). BRAF and NRF2 mutations, which are novel in TSCC, were demonstrated in one sample each. There was no significant correlation between overall mutation status and smoking history (p = 0.967) or age (p = 0.360). Positive MET alteration was associated with poorer loco-regional recurrence free survival (LRFS) of 11 months [vs 90 months in MET-negative group (p = 0.008)]. None of the other mutations were significantly correlated with LRFS or overall survival. Four of these tumors were propagated as immortalized cell lines and demonstrated the same mutations as the original tumor. CONCLUSIONS: Using the Sequenom multiplexed LungCarta panel, we identified mutations in 6 genes, TP53, STK11, MET, PIK3CA, BRAF and NRF2, with the notable absence of EGFR and HER2 mutations in our series of Asian OSCC. Primary cell line models recapitulated the mutation profiles of the original primary tumours and provide an invaluable resource for experimental cancer therapeutics.

Targeting cancer stem cell plasticity through modulation of epidermal growth factor and insulin-like growth factor receptor signaling in head and neck squamous cell cancer.

Emerging data suggest that cancer stem cells (CSCs) exist in equilibrium with differentiated cells and that stochastic transitions between these states can account for tumor heterogeneity and drug resistance. The aim of this study was to establish an in vitro system that recapitulates stem cell plasticity in head and neck squamous cell cancers (HNSCCs) and identify the factors that play a role in the maintenance and repopulation of CSCs. Tumor spheres were established using patient-derived cell lines via anchorage-independent cell culture techniques. These tumor spheres were found to have higher aldehyde dehydrogenase (ALD) cell fractions and increased expression of Kruppel-like factor 4, SRY (sex determining region Y)-box 2, and Nanog and were resistant to γ-radiation, 5-fluorouracil, cisplatin, and etoposide treatment compared with monolayer culture cells. Monolayer cultures were subject to single cell cloning to generate clones with high and low ALD fractions. ALDHigh clones showed higher expression of stem cell and epithelial-mesenchymal transition markers compared with ALDLow clones. ALD fractions, representing stem cell fractions, fluctuated with serial passaging, equilibrating at a level specific to each cell line, and could be augmented by the addition of epidermal growth factor (EGF) and/or insulin. ALDHigh clones showed increased EGF receptor (EGFR) and insulin-like growth factor-1 receptor (IGF-1R) phosphorylation, with increased activation of downstream pathways compared with ALDLow clones. Importantly, blocking these pathways using specific inhibitors against EGFR and IGF-1R reduced stem cell fractions drastically. Taken together, these results show that HNSCC CSCs exhibit plasticity, with the maintenance of the stem cell fraction dependent on the EGFR and IGF-1R pathways and potentially amenable to targeted therapeutics.

Dragon TF Association Miner: a system for exploring transcription factor associations through text-mining.

We present Dragon TF Association Miner (DTFAM), a system for text-mining of PubMed documents for potential functional association of transcription factors (TFs) with terms from Gene Ontology (GO) and with diseases. DTFAM has been trained and tested in the selection of relevant documents on a manually curated dataset containing >3000 PubMed abstracts relevant to transcription control. On our test data the system achieves sensitivity of 80% with specificity of 82%. DTFAM provides comprehensive tabular and graphical reports linking terms to relevant sets of documents. These documents are color-coded for easier inspection. DTFAM complements the existing biological resources by collecting, assessing, extracting and presenting associations that can reveal some of the not so easily observable connections among the entities found which could explain the functions of TFs and help decipher parts of gene transcriptional regulatory networks. DTFAM summarizes information from a large volume of documents saving time and making analysis simpler for individual users. DTFAM is freely available for academic and non-profit users at http://research.i2r.a-star.edu.sg/DRAGON/TFAM/.