Oncogenic KRAS Regulates Tumor Cell Signaling via Stromal Reciprocation.

Oncogenic mutations regulate signaling within both tumor cells and adjacent stromal cells. Here, we show that oncogenic KRAS (KRAS(G12D)) also regulates tumor cell signaling via stromal cells. By combining cell-specific proteome labeling with multivariate phosphoproteomics, we analyzed heterocellular KRAS(G12D) signaling in pancreatic ductal adenocarcinoma (PDA) cells. Tumor cell KRAS(G12D) engages heterotypic fibroblasts, which subsequently instigate reciprocal signaling in the tumor cells. Reciprocal signaling employs additional kinases and doubles the number of regulated signaling nodes from cell-autonomous KRAS(G12D). Consequently, reciprocal KRAS(G12D) produces a tumor cell phosphoproteome and total proteome that is distinct from cell-autonomous KRAS(G12D) alone. Reciprocal signaling regulates tumor cell proliferation and apoptosis and increases mitochondrial capacity via an IGF1R/AXL-AKT axis. These results demonstrate that oncogene signaling should be viewed as a heterocellular process and that our existing cell-autonomous perspective underrepresents the extent of oncogene signaling in cancer. VIDEO ABSTRACT.

Circulating extracellular vesicle-associated TGFß3 modulates response to cytotoxic therapy in head and neck squamous cell carcinoma.

Management of locally advanced head and neck squamous cell carcinoma (HNSCC) requires a multi-prong approach comprising surgery, radiation and/or chemotherapy, yet outcomes are limited. This is largely due to a paucity of biomarkers that can predict response to specific treatment modalities. Here, we evaluated TGFß3 protein levels in extracellular vesicles (EVs) released by HNSCC cells as a predictor for response to chemoradiation therapy (CRT). To this end, specific EV-fractions were isolated from cell lines or HNSCC patient plasma, and TGFß3 protein was quantified. In patients treated with CRT, TGFß3 levels were found to be significantly higher in plasma EV-fractions or non-responders compared with responders. High levels of TGFß3 levels in Annexin V-EVs were associated with the worst progression-free survival. In vitro experiments demonstrated that TGFß3 silencing sensitized HNSCC cells to cytotoxic therapies, and this phenotype could be rescued by treatment with exogenous. In addition, specific EV-fractions shed by cisplatin-resistant cells were sufficient to transfer the resistant phenotype to sensitive cells through activation of TGFß-signaling pathway. Therefore, our data show that TGFß3 transmitted through EV plays a significant role in response to cytotoxic therapy, which can be exploited as a potential biomarker for CRT response in HNSCC patients treated with curative intent.

Mouse RUNX1C regulates premegakaryocytic/erythroid output and maintains survival of megakaryocyte progenitors.

RUNX1 is crucial for the regulation of megakaryocyte specification, maturation, and thrombopoiesis. Runx1 possesses 2 promoters: the distal P1 and proximal P2 promoters. The major protein isoforms generated by P1 and P2 are RUNX1C and RUNX1B, respectively, which differ solely in their N-terminal amino acid sequences. RUNX1C is the most abundantly expressed isoform in adult hematopoiesis, present in all RUNX1-expressing populations, including the cKit+ hematopoietic stem and progenitor cells. RUNX1B expression is more restricted, being highly expressed in the megakaryocyte lineage but downregulated during erythropoiesis. We generated a Runx1 P1 knock-in of RUNX1B, termed P1-MRIPV This mouse line lacks RUNX1C expression but has normal total RUNX1 levels, solely comprising RUNX1B. Using this mouse line, we establish a specific requirement for the P1-RUNX1C isoform in megakaryopoiesis, which cannot be entirely compensated for by RUNX1B overexpression. P1 knock-in megakaryocyte progenitors have reduced proliferative capacity and undergo increased cell death, resulting in thrombocytopenia. P1 knock-in premegakaryocyte/erythroid progenitors demonstrate an erythroid-specification bias, evident from increased erythroid colony-forming ability and decreased megakaryocyte output. At a transcriptional level, multiple erythroid-specific genes are upregulated and megakaryocyte-specific transcripts are downregulated. In addition, proapoptotic pathways are activated in P1 knock-in premegakaryocyte/erythroid progenitors, presumably accounting for the increased cell death in the megakaryocyte progenitor compartment. Unlike in the conditional adult Runx1 null models, megakaryocytic maturation is not affected in the P1 knock-in mice, suggesting that RUNX1B can regulate endomitosis and thrombopoiesis. Therefore, despite the high degree of structural similarity, RUNX1B and RUNX1C isoforms have distinct and specific roles in adult megakaryopoiesis.

RUNX1B Expression Is Highly Heterogeneous and Distinguishes Megakaryocytic and Erythroid Lineage Fate in Adult Mouse Hematopoiesis.

The Core Binding Factor (CBF) protein RUNX1 is a master regulator of definitive hematopoiesis, crucial for hematopoietic stem cell (HSC) emergence during ontogeny. RUNX1 also plays vital roles in adult mice, in regulating the correct specification of numerous blood lineages. Akin to the other mammalian Runx genes, Runx1 has two promoters P1 (distal) and P2 (proximal) which generate distinct protein isoforms. The activities and specific relevance of these two promoters in adult hematopoiesis remain to be fully elucidated. Utilizing a dual reporter mouse model we demonstrate that the distal P1 promoter is broadly active in adult hematopoietic stem and progenitor cell (HSPC) populations. By contrast the activity of the proximal P2 promoter is more restricted and its upregulation, in both the immature Lineage- Sca1high cKithigh (LSK) and bipotential Pre-Megakaryocytic/Erythroid Progenitor (PreMegE) populations, coincides with a loss of erythroid (Ery) specification. Accordingly the PreMegE population can be prospectively separated into "pro-erythroid" and "pro-megakaryocyte" populations based on Runx1 P2 activity. Comparative gene expression analyses between Runx1 P2+ and P2- populations indicated that levels of CD34 expression could substitute for P2 activity to distinguish these two cell populations in wild type (WT) bone marrow (BM). Prospective isolation of these two populations will enable the further investigation of molecular mechanisms involved in megakaryocytic/erythroid (Mk/Ery) cell fate decisions. Having characterized the extensive activity of P1, we utilized a P1-GFP homozygous mouse model to analyze the impact of the complete absence of Runx1 P1 expression in adult mice and observed strong defects in the T cell lineage. Finally, we investigated how the leukemic fusion protein AML1-ETO9a might influence Runx1 promoter usage. Short-term AML1-ETO9a induction in BM resulted in preferential P2 upregulation, suggesting its expression may be important to establish a pre-leukemic environment.

Correction: RUNX1B Expression Is Highly Heterogeneous and Distinguishes Megakaryocytic and Erythroid Lineage Fate in Adult Mouse Hematopoiesis.

[This corrects the article DOI: 10.1371/journal.pgen.1005814.].

Paracrine factors of human fetal MSCs inhibit liver cancer growth through reduced activation of IGF-1R/PI3K/Akt signaling.

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death in the world. The multikinase inhibitor sorafenib only demonstrated marginal improvement in overall survival for advanced disease prompted the search for alternative treatment options. Human mesenchymal stem cells (MSCs) have the ability to home to tumor cells. However, its functional roles on the tumor microenvironment remain controversial. Herein, we showed that conditioned media derived from human fetal MSC (CM-hfMSCs) expressed high level of the insulin growth factor binding proteins IGFBPs and can sequester free insulin-like growth factors (IGFs) to inhibit HCC cell proliferation. The inhibitory effect of IGFBPs on IGF signaling was further evident from the reduction of activated IGF-1R and PI3K/Akt, leading eventually to the induction of cell cycle arrest. We also demonstrated that CM-hfMSCs could enhance the therapeutic efficacy of sorafenib and sunitinib. To the best of our knowledge, this is the first report to show that CM-hfMSCs has a tumor-specific, antiproliferative effect that is not observed with normal human hepatocyte cells and patient-derived matched normal tissues. Our results thus suggest that CM-hfMSCs can provide a useful tool to design alternative/adjuvant treatment strategies for HCC, especially in related function to potentiate the effects of chemotherapeutic drugs.

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.

Hit screening with multivariate robust outlier detection.

Hit screening, which involves the identification of compounds or targets capable of modulating disease-relevant processes, is an important step in drug discovery. Some assays, such as image-based high-content screenings, produce complex multivariate readouts. To fully exploit the richness of such data, advanced analytical methods that go beyond the conventional univariate approaches should be employed. In this work, we tackle the problem of hit identification in multivariate assays. As with univariate assays, a hit from a multivariate assay can be defined as a candidate that yields an assay value sufficiently far away in distance from the mean or central value of inactives. Viewed another way, a hit is an outlier from the distribution of inactives. A method was developed for identifying multivariate hit in high-dimensional data sets based on principal components and robust Mahalanobis distance (the multivariate analogue to the Z- or T-statistic). The proposed method, termed mROUT (multivariate robust outlier detection), demonstrates superior performance over other techniques in the literature in terms of maintaining Type I error, false discovery rate and true discovery rate in simulation studies. The performance of mROUT is also illustrated on a CRISPR knockout data set from in-house phenotypic screening programme.

Therapeutic application of extracellular vesicular EGFR isoform D as a co-drug to target squamous cell cancers with tyrosine kinase inhibitors.

Targeting wild-type epidermal growth factor receptor (EGFR) using tyrosine kinase inhibitors (TKIs) never achieved its purported success in cancers such as head and neck squamous cell carcinoma, which are largely EGFR-dependent. We had previously shown that exceptional responders to TKIs have a genetic aberration that results in overexpression of an EGFR splice variant, isoform D (IsoD). IsoD lacks an integral transmembrane and kinase domain and is secreted in extracellular vesicles (EVs) in TKI-sensitive patient-derived cultures. Remarkably, the exquisite sensitivity to TKIs could be transferred to TKI-resistant tumor cells, and IsoD protein in the EV is necessary and sufficient to transfer the phenotype in vitro and in vivo across multiple models and drugs. This drug response requires an intact endocytic mechanism, binding to full-length EGFR, and signaling through Src-phosphorylation within the endosomal compartment. We propose a therapeutic strategy using EVs containing EGFR IsoD as a co-drug to expand the use of TKI therapy to EGFR-driven cancers.

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.

High-throughput immunoglobulin repertoire analysis distinguishes between human IgM memory and switched memory B-cell populations.

B-cell receptor (BCR) diversity is achieved centrally by rearrangement of Variable, Diversity, and Joining genes, and peripherally by somatic hypermutation and class-switching of the rearranged genes. Peripheral B-cell populations are subject to both negative and positive selection events in the course of their development that have the potential to shape the BCR repertoire. The origin of IgM(+)IgD(+)CD27(+) (IgM memory) cells is controversial. It has been suggested that they may be a prediversified, antigen-independent, population of cells or that they are a population of cells that develop in response to T-independent antigens. Most recently, it was suggested that the majority of IgM memory cells are directly related to switched memory cells and are early emigrants from the germinal center reaction. Advances in sequencing technology have enabled us to undertake large scale IGH repertoire analysis of transitional, naive, IgM memory and switched memory B-cell populations. We find that the memory B-cell repertoires differ from the transitional and naive repertoires, and that the IgM memory repertoire is distinct from that of class-switched memory. Thus we conclude that a large proportion of IgM memory cells develop in response to different stimuli than for class-switched memory cell development.

Text-based over-representation analysis of microarray gene lists with annotation bias.

A major challenge in microarray data analysis is the functional interpretation of gene lists. A common approach to address this is over-representation analysis (ORA), which uses the hypergeometric test (or its variants) to evaluate whether a particular functionally defined group of genes is represented more than expected by chance within a gene list. Existing applications of ORA have been largely limited to pre-defined terminologies such as GO and KEGG. We report our explorations of whether ORA can be applied to a wider mining of free-text. We found that a hitherto underappreciated feature of experimentally derived gene lists is that the constituents have substantially more annotation associated with them, as they have been researched upon for a longer period of time. This bias, a result of patterns of research activity within the biomedical community, is a major problem for classical hypergeometric test-based ORA approaches, which cannot account for such bias. We have therefore developed three approaches to overcome this bias, and demonstrate their usability in a wide range of published datasets covering different species. A comparison with existing tools that use GO terms suggests that mining PubMed abstracts can reveal additional biological insight that may not be possible by mining pre-defined ontologies alone.

CKAP2L Promotes Non-Small Cell Lung Cancer Progression through Regulation of Transcription Elongation.

Chromosomal instability (CIN) is a driver of clonal diversification and intratumor heterogeneity, providing genetic diversity that contributes to tumor progression. It is estimated that approximately 80% of solid cancers, including non-small cell lung cancer (NSCLC), exhibit features of CIN, which affects tumor growth and response to therapy. However, the molecular mechanisms connecting CIN to tumor progression are still poorly understood. Through an RNAi screen performed on genes involved in CIN and overexpressed in human lung adenocarcinoma samples, we identified the cytoskeleton-associated protein 2-like (CKAP2L) as a potential oncogene that promotes lung cancer proliferation and growth in vitro and in vivo. Mechanistically, CKAP2L directly interacted with RNA Pol II and regulated transcription elongation of key genes involved in spindle assembly checkpoint, chromosome segregation, cell cycle, and E2F signaling. Furthermore, depletion of CKAP2L increased the sensitivity of NSCLC cells to alvocidib, a pan-CDK inhibitor, leading to a significant reduction of cell proliferation and an increase in cell death. Altogether, these findings shed light on the molecular mechanisms through which CKAP2L, a protein involved in CIN, promotes cancer progression and suggest that its inhibition represents a novel therapeutic strategy in NSCLC. SIGNIFICANCE: These findings demonstrate the oncogenic function of CKAP2L through regulation of transcription elongation and suggest that targeting CKAP2L could enhance therapeutic response in patients with NSCLC.

A KRAS-responsive long non-coding RNA controls microRNA processing.

Wild-type KRAS (KRASWT) amplification has been shown to be a secondary means of KRAS activation in cancer and associated with poor survival. Nevertheless, the precise role of KRASWT overexpression in lung cancer progression is largely unexplored. Here, we identify and characterize a KRAS-responsive lncRNA, KIMAT1 (ENSG00000228709) and show that it correlates with KRAS levels both in cell lines and in lung cancer specimens. Mechanistically, KIMAT1 is a MYC target and drives lung tumorigenesis by promoting the processing of oncogenic microRNAs (miRNAs) through DHX9 and NPM1 stabilization while halting the biogenesis of miRNAs with tumor suppressor function via MYC-dependent silencing of p21, a component of the Microprocessor Complex. KIMAT1 knockdown suppresses not only KRAS expression but also KRAS downstream signaling, thereby arresting lung cancer growth in vitro and in vivo. Taken together, this study uncovers a role for KIMAT1 in maintaining a positive feedback loop that sustains KRAS signaling during lung cancer progression and provides a proof of principle that interfering with KIMAT1 could be a strategy to hamper KRAS-induced tumorigenesis.

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).

Pre-clinical activity of combined LSD1 and mTORC1 inhibition in MLL-translocated acute myeloid leukaemia.

The histone demethylase lysine-specific demethylase 1 (LSD1 or KDM1A) has emerged as a candidate therapeutic target in acute myeloid leukaemia (AML); tranylcypromine-derivative inhibitors induce loss of clonogenic activity and promote differentiation, in particular in the MLL-translocated molecular subtype of AML. In AML, the use of drugs in combination often delivers superior clinical activity. To identify genes and cellular pathways that collaborate with LSD1 to maintain the leukaemic phenotype, and which could be targeted by combination therapies, we performed a genome-wide CRISPR-Cas9 dropout screen. We identified multiple components of the amino acid sensing arm of mTORC1 signalling-RRAGA, MLST8, WDR24 and LAMTOR2-as cellular sensitizers to LSD1 inhibition. Knockdown of mTORC1 components, or mTORC1 pharmacologic inhibition, in combination with LSD1 inhibition enhanced differentiation in both cell line and primary cell settings, in vitro and in vivo, and substantially reduced the frequency of clonogenic primary human AML cells in a modelled minimal residual disease setting. Synergistic upregulation of a set of transcription factor genes associated with terminal monocytic lineage differentiation was observed. Thus, dual mTORC1 and LSD1 inhibition represents a candidate combination approach for enhanced differentiation in MLL-translocated AML which could be evaluated in early phase clinical trials.

Profiling of Circulating Free DNA Using Targeted and Genome-wide Sequencing in Patients with SCLC.

INTRODUCTION: SCLC accounts for approximately 250,000 deaths worldwide each year. Acquisition of adequate tumor biopsy samples is challenging, and liquid biopsies present an alternative option for patient stratification and response monitoring. METHODS: We applied whole genome next-generation sequencing to circulating free DNA (cfDNA) from 39 patients with limited-stage (LS) SCLC and 30 patients with extensive-stage SCLC to establish genome-wide copy number aberrations and also performed targeted mutation analysis of 110 SCLC associated genes. Quantitative metrics were calculated for copy number aberrations, including percent genome amplified (PGA [the percentage of genomic regions amplified]), Z-score (a measure of standard deviation), and Moran's I (a measure of spatial autocorrelation). In addition CellSearch, an epitope-dependent enrichment platform, was used to enumerate circulating tumor cells (CTCs) from a parallel blood sample. RESULTS: Genome-wide and targeted cfDNA sequencing data identified tumor-related changes in 94% of patients with LS SCLC and 100% of patients with extensive-stage SCLC. Parallel analysis of CTCs based on at least 1 CTC/7.5 mL of blood increased tumor detection frequencies to 95% for LS SCLC. Both CTC counts and cfDNA readouts correlated with disease stage and overall survival. CONCLUSIONS: We demonstrate that a simple cfDNA genome-wide copy number approach provides an effective means of monitoring patients through treatment and show that targeted cfDNA sequencing identifies potential therapeutic targets in more than 50% of patients. We are now incorporating this approach into additional studies and trials of targeted therapies.

Vulnerability of drug-resistant EML4-ALK rearranged lung cancer to transcriptional inhibition.

A subset of lung adenocarcinomas is driven by the EML4-ALK translocation. Even though ALK inhibitors in the clinic lead to excellent initial responses, acquired resistance to these inhibitors due to on-target mutations or parallel pathway alterations is a major clinical challenge. Exploring these mechanisms of resistance, we found that EML4-ALK cells parental or resistant to crizotinib, ceritinib or alectinib are remarkably sensitive to inhibition of CDK7/12 with THZ1 and CDK9 with alvocidib or dinaciclib. These compounds robustly induce apoptosis through transcriptional inhibition and downregulation of anti-apoptotic genes. Importantly, alvocidib reduced tumour progression in xenograft mouse models. In summary, our study takes advantage of the transcriptional addiction hypothesis to propose a new treatment strategy for a subset of patients with acquired resistance to first-, second- and third-generation ALK inhibitors.

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.