An AI-assisted integrated, scalable, single-cell phenomic-transcriptomic platform to elucidate intratumor heterogeneity against immune response.

We present a novel framework combining single-cell phenotypic data with single-cell transcriptomic analysis to identify factors underpinning heterogeneity in antitumor immune response. We developed a pairwise, tumor-immune discretized interaction assay between natural killer (NK-92MI) cells and patient-derived head and neck squamous cell carcinoma (HNSCC) cell lines on a microfluidic cell-trapping platform. Furthermore we generated a deep-learning computer vision algorithm that is capable of automating the acquisition and analysis of a large, live-cell imaging data set (>1 million) of paired tumor-immune interactions spanning a time course of 24 h across multiple HNSCC lines (n = 10). Finally, we combined the response data measured by Kaplan-Meier survival analysis against NK-mediated killing with downstream single-cell transcriptomic analysis to interrogate molecular signatures associated with NK-effector response. As proof-of-concept for the proposed framework, we efficiently identified MHC class I-driven cytotoxic resistance as a key mechanism for immune evasion in nonresponders, while enhanced expression of cell adhesion molecules was found to be correlated with sensitivity against NK-mediated cytotoxicity. We conclude that this integrated, data-driven phenotypic approach holds tremendous promise in advancing the rapid identification of new mechanisms and therapeutic targets related to immune evasion and response.

Single cell landscape of functionally cured chronic hepatitis B patients reveals activation of innate and altered CD4-CTL-driven adaptive immunity.

BACKGROUND & AIMS: Hepatitis B surface antigen (HBsAg) loss or functional cure (FC), is considered the desirable therapeutic outcome for chronic hepatitis B (CHB) patients. However, the immune-pathological biomarkers and underlying mechanisms remain unclear. In this study we comprehensively interrogate disease-associated cell states (DACS) identified within intra-hepatic tissue and matched PBMCs from either CHB or FC patients, at the resolution of single cells, to provide novel insights into putative mechanisms underlying FC. METHODS: We combined single cell transcriptomics (scRNA-seq) with multiparametric flow cytometry-based immune phenotyping, and multiplexed immunofluorescence to elucidate the immunopathological cell states associated with CHB vs FC. RESULTS: We find that the intra-hepatic environment of CHB and FC patients displays specific cell identities and molecular signatures that are distinct from those found in matched PBMCs. FC is associated with emergence of an altered adaptive immune response marked by CD4 cytotoxic T lymphocytes (CD4-CTLs), and an activated innate response represented by liver-resident natural killer (LR-NK) cells, specific Kupffer cell (KC) subtypes and marginated neutrophils. Surprisingly, we find that FC patients are also characterized by the presence of MHC class II-expressing hepatocytes and low but persistent levels of cccDNA and pgRNA, which may play an important role in achieving functional cure in HBV patients. CONCLUSIONS: Our study provides conceptually novel insights into the immuno-pathological control of HBV cure, and opens exciting new avenues for clinical management, biomarker discovery and therapeutic interventions. We believe that the discoveries from this study, as it relates to the activation of an innate and altered immune response that may facilitate sustained, low-grade inflammation, may have broader implications in the resolution of chronic viral hepatitis.

Rising global burden of cancer attributable to high BMI from 2010 to 2019.

BACKGROUND: High body mass index (BMI) is a major risk factor for cancer development, but its impact on the global burden of cancer remains unclear. METHODS: We estimated global and regional temporal trends in the burden of cancer attributable to high BMI, and the contributions of various cancer types using the framework of the Global Burden of Disease Study. RESULTS: From 2010 to 2019, there was a 35 % increase in deaths and a 34 % increase in disability-adjusted life-years from cancers attributable to high BMI. The age-standardized death rates for cancer attributable to high BMI increased over the study period (annual percentage change [APC] +0.48 %, 95 % CI 0.22 to 0.74 %). The greatest number of deaths from cancer attributable to high BMI occurred in Europe, but the fastest-growing age-standardized death rates and disability-adjusted life-years occurred in Southeast Asia. Liver cancer was the fastest-growing cause of cancer mortality (APC: 1.37 %, 95 % CI 1.25 to 1.49 %) attributable to high BMI. CONCLUSION: The global burden of cancer-related deaths attributable to high BMI has increased substantially from 2010 to 2019. The greatest increase in age-standardized death rates occurred in Southeast Asia, and liver cancer is the fastest-growing cause of cancer mortality attributable to high BMI. Urgent and sustained measures are required at a global and regional level to reverse these trends and slow the growing burden of cancer attributed to high BMI.

CD38 marks the exhausted CD8+ tissue-resident memory T cells in hepatocellular carcinoma.

Introduction: Despite recent advances in immunotherapy for hepatocellular carcinoma (HCC), the overall modest response rate underscores the need for a better understanding of the tumor microenvironment (TME) of HCC. We have previously shown that CD38 is widely expressed on tumor-infiltrating leukocytes (TILs), predominantly on CD3+ T cells and monocytes. However, its specific role in the HCC TME remains unclear. Methods: In this current study, we used cytometry time-of-flight (CyTOF), bulk RNA sequencing on sorted T cells, and single-cell RNA (scRNA) sequencing to interrogate expression of CD38 and its correlation with T cell exhaustion in HCC samples. We also employed multiplex immunohistochemistry (mIHC) for validating our findings. Results: From CyTOF analysis, we compared the immune composition of CD38-expressing leukocytes in TILs, non-tumor tissue-infiltrating leukocytes (NIL), and peripheral blood mononuclear cells (PBMC). We identified CD8+ T cells as the dominant CD38-expressing TILs and found that CD38 expression was significantly higher in CD8+ TRM in TILs than in NILs. Furthermore, through transcriptomic analysis on sorted CD8+ TRM from HCC tumors, we observed a higher expression of CD38 along with T cell exhaustion genes, including PDCD1 and CTLA4, compared to the circulating memory CD8 T cells from PBMC. This was validated by scRNA sequencing that revealed co-expression of CD38 with PDCD1, CTLA4, and ITGAE (CD103) in T cells from HCC tumors. The protein co-expression of CD38 and PD-1 on CD8+ T cells was further demonstrated by mIHC on HCC FFPE tissues, marking CD38 as a T cell co-exhaustion marker in HCC. Lastly, the higher proportions of CD38+PD-1+ CD8+ T cells and CD38+PD-1+ TRM were significantly associated with the higher histopathological grades of HCC, indicating its role in the aggressiveness of the disease. Conclusion: Taken together, the concurrent expression of CD38 with exhaustion markers on CD8+ TRM underpins its role as a key marker of T cell exhaustion and a potential therapeutic target for restoring cytotoxic T cell function in HCC.

Identification of mechanism of cancer-cell-specific reactivation of hTERT offers therapeutic opportunities for blocking telomerase specifically in human colorectal cancer.

Transcriptional reactivation of hTERT is the limiting step in tumorigenesis. While mutations in hTERT promoter present in 19% of cancers are recognized as key drivers of hTERT reactivation, mechanisms by which wildtype hTERT (WT-hTERT) promoter is reactivated, in majority of human cancers, remain unknown. Using primary colorectal cancers (CRC) we identified Tert INTeracting region 2 (T-INT2), the critical chromatin region essential for reactivating WT-hTERT promoter in CRCs. Elevated ß-catenin and JunD level in CRC facilitates chromatin interaction between hTERT promoter and T-INT2 that is necessary to turn on hTERTexpression. Pharmacological screens uncovered salinomycin, which inhibits JunD mediated hTERT-T-INT2 interaction that is required for the formation of a stable transcription complex on the hTERT promoter. Our results showed for the first time how known CRC alterations, such as APC, lead to WT-hTERT promoter reactivation during stepwise-tumorigenesis and provide a new perspective for developing cancer-specific drugs.

Healthy and cancer cells harbor the same DNA sequence, but reactivation of the Human Telomerase Reverse Transcriptase (hTERT) gene is observed only in cancer cells. How does that happen was not known for over three decades of research? This study identifies a specific DNA structure that forms only in cancer cells and brings the necessary molecular machinery into the correct position to activate the hTERT gene. The detailed mechanism of hTERT activation provided in this study will be instrumental in designing cancer cell-specific hTERT inhibitors, especially since all the other ways of inhibiting telomerase failed in the clinic.

Evolution of precision oncology-guided treatment paradigms.

Cancer treatment is gradually evolving from the classical use of nonspecific cytotoxic drugs targeting generic mechanisms of cell growth and proliferation. Instead, new "patient-specific treatment paradigms" that are based on an individual patient's tumor-specific molecular features are emerging, and these include "druggable" genomic alterations such as oncogenic driver mutations, downstream activities of cancer-signaling pathways, and the expression of specific genes involved in tumorigenesis and cancer progression. This evolving landscape of making evidence-based treatment decisions forms the foundation of precision oncology, which aims to deliver "the right drug, to the right patient and at the right time". The long-term vision for this approach is to maximize the treatment efficacy while minimizing exposure to ineffective therapy and reducing co-morbidity-related side effects. Successful clinical translation and implementation of this vision have the potential to revolutionize treatment paradigms from predominantly reactive, to more evidence-based, proactive and predictive care. In this article, we review the past and current approaches in precision oncology, and describe their remarkable power and limitations. We also speculate on the evolution of newly emerging methodologies of the future that can be used to address some of the key challenges associated with the existing paradigms. This article is categorized under: Cancer > Genetics/Genomics/Epigenetics Cancer > Molecular and Cellular Physiology Cancer > Computational Models.

Breakage fusion bridge cycles drive high oncogene copy number, but not intratumoral genetic heterogeneity or rapid cancer genome change.

Oncogene amplification is a major driver of cancer pathogenesis. Breakage fusion bridge (BFB) cycles, like extrachromosomal DNA (ecDNA), can lead to high copy numbers of oncogenes, but their impact on intratumoral heterogeneity, treatment response, and patient survival are not well understood due to difficulty in detecting them by DNA sequencing. We describe a novel algorithm that detects and reconstructs BFB amplifications using optical genome maps (OGMs), called OM2BFB. OM2BFB showed high precision (>93%) and recall (92%) in detecting BFB amplifications in cancer cell lines, PDX models and primary tumors. OM-based comparisons demonstrated that short-read BFB detection using our AmpliconSuite (AS) toolkit also achieved high precision, albeit with reduced sensitivity. We detected 371 BFB events using whole genome sequences from 2,557 primary tumors and cancer lines. BFB amplifications were preferentially found in cervical, head and neck, lung, and esophageal cancers, but rarely in brain cancers. BFB amplified genes show lower variance of gene expression, with fewer options for regulatory rewiring relative to ecDNA amplified genes. BFB positive (BFB (+)) tumors showed reduced heterogeneity of amplicon structures, and delayed onset of resistance, relative to ecDNA(+) tumors. EcDNA and BFB amplifications represent contrasting mechanisms to increase the copy numbers of oncogene with markedly different characteristics that suggest different routes for intervention.

Targeting the bicarbonate transporter SLC4A4 overcomes immunosuppression and immunotherapy resistance in pancreatic cancer.

Solid tumors are generally characterized by an acidic tumor microenvironment (TME) that favors cancer progression, therapy resistance and immune evasion. By single-cell RNA-sequencing analysis in individuals with pancreatic ductal adenocarcinoma (PDAC), we reveal solute carrier family 4 member 4 (SLC4A4) as the most abundant bicarbonate transporter, predominantly expressed by epithelial ductal cells. Functionally, SLC4A4 inhibition in PDAC cancer cells mitigates the acidosis of the TME due to bicarbonate accumulation in the extracellular space and a decrease in lactate production by cancer cells as the result of reduced glycolysis. In PDAC-bearing mice, genetic or pharmacological SLC4A4 targeting improves T cell-mediated immune response and breaches macrophage-mediated immunosuppression, thus inhibiting tumor growth and metastases. In addition, Slc4a4 targeting in combination with immune checkpoint blockade is able to overcome immunotherapy resistance and prolong survival. Overall, our data propose SLC4A4 as a therapeutic target to unleash an antitumor immune response in PDAC.

A comparative study of tumour-on-chip models with patient-derived xenografts for predicting chemotherapy efficacy in colorectal cancer patients.

Inter-patient and intra-tumour heterogeneity (ITH) have prompted the need for a more personalised approach to cancer therapy. Although patient-derived xenograft (PDX) models can generate drug response specific to patients, they are not sustainable in terms of cost and time and have limited scalability. Tumour Organ-on-Chip (OoC) models are in vitro alternatives that can recapitulate some aspects of the 3D tumour microenvironment and can be scaled up for drug screening. While many tumour OoC systems have been developed to date, there have been limited validation studies to ascertain whether drug responses obtained from tumour OoCs are comparable to those predicted from patient-derived xenograft (PDX) models. In this study, we established a multiplexed tumour OoC device, that consists of an 8 × 4 array (32-plex) of culture chamber coupled to a concentration gradient generator. The device enabled perfusion culture of primary PDX-derived tumour spheroids to obtain dose-dependent response of 5 distinct standard-of-care (SOC) chemotherapeutic drugs for 3 colorectal cancer (CRC) patients. The in vitro efficacies of the chemotherapeutic drugs were rank-ordered for individual patients and compared to the in vivo efficacy obtained from matched PDX models. We show that quantitative correlation analysis between the drug efficacies predicted via the microfluidic perfusion culture is predictive of response in animal PDX models. This is a first study showing a comparative framework to quantitatively correlate the drug response predictions made by a microfluidic tumour organ-on-chip (OoC) model with that of PDX animal models.

Single-cell and bulk transcriptome sequencing identifies two epithelial tumor cell states and refines the consensus molecular classification of colorectal cancer.

The consensus molecular subtype (CMS) classification of colorectal cancer is based on bulk transcriptomics. The underlying epithelial cell diversity remains unclear. We analyzed 373,058 single-cell transcriptomes from 63 patients, focusing on 49,155 epithelial cells. We identified a pervasive genetic and transcriptomic dichotomy of malignant cells, based on distinct gene expression, DNA copy number and gene regulatory network. We recapitulated these subtypes in bulk transcriptomes from 3,614 patients. The two intrinsic subtypes, iCMS2 and iCMS3, refine CMS. iCMS3 comprises microsatellite unstable (MSI-H) cancers and one-third of microsatellite-stable (MSS) tumors. iCMS3 MSS cancers are transcriptomically more similar to MSI-H cancers than to other MSS cancers. CMS4 cancers had either iCMS2 or iCMS3 epithelium; the latter had the worst prognosis. We defined the intrinsic epithelial axis of colorectal cancer and propose a refined 'IMF' classification with five subtypes, combining intrinsic epithelial subtype (I), microsatellite instability status (M) and fibrosis (F).

Trajectory of immune evasion and cancer progression in hepatocellular carcinoma.

Immune evasion is key to cancer initiation and later at metastasis, but its dynamics at intermediate stages, where potential therapeutic interventions could be applied, is undefined. Here we show, using multi-dimensional analyses of resected tumours, their adjacent non-tumour tissues and peripheral blood, that extensive immune remodelling takes place in patients with stage I to III hepatocellular carcinoma (HCC). We demonstrate the depletion of anti-tumoural immune subsets and accumulation of immunosuppressive or exhausted subsets along with reduced tumour infiltration of CD8 T cells peaking at stage II tumours. Corresponding transcriptomic modification occur in the genes related to antigen presentation, immune responses, and chemotaxis. The progressive immune evasion is validated in a murine model of HCC. Our results show evidence of ongoing tumour-immune co-evolution during HCC progression and offer insights into potential interventions to reverse, prevent or limit the progression of the disease.

Genome-wide screens identify specific drivers of mutant hTERT promoters.

Cancer-specific hTERT promoter mutations reported in 19% of cancers result in enhanced telomerase activity. Understanding the distinctions between transcriptional regulation of wild-type (WT) and mutant (Mut) hTERT promoters may open up avenues for development of inhibitors which specially block hTERT expression in cancer cells. To comprehensively identify physiological regulators of WT- or Mut-hTERT promoters, we generated several isogenic reporter cells driven by endogenous hTERT loci. Genome-wide CRISPR-Cas9 and small interfering RNA screens using these isogenic reporter lines identified specific regulators of Mut-hTERT promoters. We validate and characterize one of these hits, namely, MED12, a kinase subunit of mediator complex. We demonstrate that MED12 specifically drives expression of hTERT from the Mut-hTERT promoter by mediating long-range chromatin interaction between the proximal Mut-hTERT promoter and T-INT1 distal regulatory region 260 kb upstream. Several hits identified in our screens could serve as potential therapeutic targets, inhibition of which may specifically block Mut-hTERT promoter driven telomerase reactivation in cancers.

Predicting heterogeneity in clone-specific therapeutic vulnerabilities using single-cell transcriptomic signatures.

While understanding molecular heterogeneity across patients underpins precision oncology, there is increasing appreciation for taking intra-tumor heterogeneity into account. Based on large-scale analysis of cancer omics datasets, we highlight the importance of intra-tumor transcriptomic heterogeneity (ITTH) for predicting clinical outcomes. Leveraging single-cell RNA-seq (scRNA-seq) with a recommender system (CaDRReS-Sc), we show that heterogeneous gene-expression signatures can predict drug response with high accuracy (80%). Using patient-proximal cell lines, we established the validity of CaDRReS-Sc's monotherapy (Pearson r>0.6) and combinatorial predictions targeting clone-specific vulnerabilities (>10% improvement). Applying CaDRReS-Sc to rapidly expanding scRNA-seq compendiums can serve as in silico screen to accelerate drug-repurposing studies. Availability: https://github.com/CSB5/CaDRReS-Sc .

LncRNA SFTA1P mediates positive feedback regulation of the Hippo-YAP/TAZ signaling pathway in non-small cell lung cancer.

Long non-coding RNAs (lncRNAs) regulate numerous biological processes involved in both development and carcinogenesis. Hippo-YAP/TAZ signaling, a critical pathway responsible for organ size control, is often dysregulated in a variety of cancers. However, the nature and function of YAP/TAZ-regulated lncRNAs during tumorigenesis remain largely unexplored. By profiling YAP/TAZ-regulated lncRNAs, we identified SFTA1P as a novel transcriptional target and a positive feedback regulator of YAP/TAZ signaling. Using non-small cell lung cancer (NSCLC) cell lines, we show that SFTA1P is transcriptionally activated by YAP/TAZ in a TEAD-dependent manner. Functionally, knockdown of SFTA1P in NSCLC cell lines inhibited proliferation, induced programmed cell death, and compromised their tumorigenic potential. Mechanistically, SFTA1P knockdown decreased TAZ protein abundance and consequently, the expression of YAP/TAZ transcriptional targets. We provide evidence that this phenomenon could potentially be mediated via its interaction with TAZ mRNA to regulate TAZ translation. Our results reveal SFTA1P as a positive feedback regulator of Hippo-YAP/TAZ signaling, which may serve as the molecular basis for lncRNA-based therapies against YAP/TAZ-driven cancers.

Molecular docking-aided identification of small molecule inhibitors targeting ß-catenin-TCF4 interaction.

Here we report a molecular docking-based approach to identify small molecules that can target the ß-catenin (ß-cat)-TCF4 protein-protein interaction (PPI), a key effector complex for nuclear Wnt signaling activity. Specifically, we developed and optimized a computational model of ß-cat using publicly available ß-cat protein crystal structures, and existing ß-cat-TCF4 interaction inhibitors as the training set. Using our computational model to an in silico screen predicted 27 compounds as good binders to ß-cat, of which 3 were identified to be effective against a Wnt-responsive luciferase reporter. In vitro functional validation experiments revealed GB1874 as an inhibitor of the Wnt pathway that targets the ß-cat-TCF4 PPI. GB1874 also affected the proliferation and stemness of Wnt-addicted colorectal cancer (CRC) cells in vitro. Encouragingly, GB1874 inhibited the growth of CRC tumor xenografts in vivo, thus demonstrating its potential for further development into therapeutics against Wnt-associated cancer indications.

Single-Cell RNA-seq Reveals Angiotensin-Converting Enzyme 2 and Transmembrane Serine Protease 2 Expression in TROP2+ Liver Progenitor Cells: Implications in Coronavirus Disease 2019-Associated Liver Dysfunction.

The recent coronavirus disease 2019 (COVID-19) pandemic is caused by severe acute respiratory syndrome coronavirus 2. COVID-19 was first reported in China (December 2019) and is now prevalent across the globe. Entry of severe acute respiratory syndrome coronavirus 2 into mammalian cells requires the binding of viral Spike (S) proteins to the angiotensin-converting enzyme 2 receptor. Once entered, the S protein is primed by a specialized serine protease, transmembrane serine protease 2 in the host cell. Importantly, besides the respiratory symptoms that are consistent with other common respiratory virus infections when patients become viremic, a significant number of COVID-19 patients also develop liver comorbidities. We explored whether a specific target cell-type in the mammalian liver could be implicated in disease pathophysiology other than the general deleterious response to cytokine storms. Here, we used single-cell RNA-seq to survey the human liver and identified potentially implicated liver cell-type for viral ingress. We analyzed ~300,000 single cells across five different (i.e., human fetal, healthy, cirrhotic, tumor, and adjacent normal) liver tissue types. This study reports on the co-expression of angiotensin-converting enzyme 2 and transmembrane serine protease 2 in a TROP2+ liver progenitor population. Importantly, we detected enrichment of this cell population in the cirrhotic liver when compared with tumor tissue. These results indicated that in COVID-19-associated liver dysfunction and cell death, a viral infection of TROP2+ progenitors in the liver might significantly impair liver regeneration in patients with liver cirrhosis.

DUSP16 promotes cancer chemoresistance through regulation of mitochondria-mediated cell death.

Drug resistance is a major obstacle to the treatment of most human tumors. In this study, we find that dual-specificity phosphatase 16 (DUSP16) regulates resistance to chemotherapy in nasopharyngeal carcinoma, colorectal cancer, gastric and breast cancer. Cancer cells expressing higher DUSP16 are intrinsically more resistant to chemotherapy-induced cell death than cells with lower DUSP16 expression. Overexpression of DUSP16 in cancer cells leads to increased resistance to cell death upon chemotherapy treatment. In contrast, knockdown of DUSP16 in cancer cells increases their sensitivity to treatment. Mechanistically, DUSP16 inhibits JNK and p38 activation, thereby reducing BAX accumulation in mitochondria to reduce apoptosis. Analysis of patient survival in head & neck cancer and breast cancer patient cohorts supports DUSP16 as a marker for sensitivity to chemotherapy and therapeutic outcome. This study therefore identifies DUSP16 as a prognostic marker for the efficacy of chemotherapy, and as a therapeutic target for overcoming chemoresistance in cancer.

miR-582-5p Is a Tumor Suppressor microRNA Targeting the Hippo-YAP/TAZ Signaling Pathway in Non-Small Cell Lung Cancer.

The Hippo-YAP/TAZ signaling pathway is an evolutionarily conserved signaling pathway involved in a broad spectrum of biological processes, including tumorigenesis. Whilst aberrant Hippo-YAP/TAZ signaling is frequently reported in various cancers, the genetic alterations of this pathway are relatively rare, suggesting regulation at the post-transcriptional level. MicroRNAs play key role in tumorigenesis by regulating gene expression post-transcriptionally. Amongst the cancer-relevant microRNAs, miR-582-5p suppresses cell growth and tumorigenesis by inhibiting the expression of oncogenes, including AKT3, MAP3K2 and NOTCH1. Given the oncogenic role of YAP/TAZ in solid tumors, we scrutinized the possible interplay between miR-582-5p and Hippo-YAP/TAZ signaling. Correlation analysis in NSCLC cells revealed a positive relationship between the expression of mature miR-582-5p and the proportion of phosphorylated YAP/TAZ. Intriguingly, YAP/TAZ knockdown reduced the expression of mature miR-582-5p but increased that of primary miR-582. Overexpression of miR-582-5p resulted in increased phosphorylation of YAP/TAZ with a concomitant reduction in cell proliferation and enhanced apoptosis. Mechanistically, we find that miR-582-5p targets actin regulators NCKAP1 and PIP5K1C, which may be responsible for the observed alteration in F-actin, known to modulate YAP/TAZ. We postulate that regulation of the actin cytoskeleton by miR-582-5p may attenuate YAP/TAZ activity. Altogether, this study reveals a novel mechanism of YAP/TAZ regulation by miR-582-5p in a cytoskeleton-dependent manner and suggests a negative feedback loop, highlighting the therapeutic potential of restoring miR-582-5p expression in treating NSCLC.

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