Tumor- and circulating-free DNA methylation identifies clinically relevant small cell lung cancer subtypes.

Small cell lung cancer (SCLC) is an aggressive malignancy composed of distinct transcriptional subtypes, but implementing subtyping in the clinic has remained challenging, particularly due to limited tissue availability. Given the known epigenetic regulation of critical SCLC transcriptional programs, we hypothesized that subtype-specific patterns of DNA methylation could be detected in tumor or blood from SCLC patients. Using genomic-wide reduced-representation bisulfite sequencing (RRBS) in two cohorts totaling 179 SCLC patients and using machine learning approaches, we report a highly accurate DNA methylation-based classifier (SCLC-DMC) that can distinguish SCLC subtypes. We further adjust the classifier for circulating-free DNA (cfDNA) to subtype SCLC from plasma. Using the cfDNA classifier (cfDMC), we demonstrate that SCLC phenotypes can evolve during disease progression, highlighting the need for longitudinal tracking of SCLC during clinical treatment. These data establish that tumor and cfDNA methylation can be used to identify SCLC subtypes and might guide precision SCLC therapy.

PancanQTLv2.0: a comprehensive resource for expression quantitative trait loci across human cancers.

Expression quantitative trait locus (eQTL) analysis is a powerful tool used to investigate genetic variations in complex diseases, including cancer. We previously developed a comprehensive database, PancanQTL, to characterize cancer eQTLs using The Cancer Genome Atlas (TCGA) dataset, and linked eQTLs with patient survival and GWAS risk variants. Here, we present an updated version, PancanQTLv2.0 (https://hanlaboratory.com/PancanQTLv2/), with advancements in fine-mapping causal variants for eQTLs, updating eQTLs overlapping with GWAS linkage disequilibrium regions and identifying eQTLs associated with drug response and immune infiltration. Through fine-mapping analysis, we identified 58 747 fine-mapped eQTLs credible sets, providing mechanic insights of gene regulation in cancer. We further integrated the latest GWAS Catalog and identified a total of 84 592 135 linkage associations between eQTLs and the existing GWAS loci, which represents a remarkable ∼50-fold increase compared to the previous version. Additionally, PancanQTLv2.0 uncovered 659516 associations between eQTLs and drug response and identified 146948 associations between eQTLs and immune cell abundance, providing potentially clinical utility of eQTLs in cancer therapy. PancanQTLv2.0 expanded the resources available for investigating gene expression regulation in human cancers, leading to advancements in cancer research and precision oncology.

Autotaxin suppresses cytotoxic T cells via LPAR5 to promote anti-PD-1 resistance in non-small cell lung cancer.

Non-small cell lung cancers that harbor concurrent KRAS and TP53 (KP) mutations are immunologically warm tumors with partial responsiveness to anti-PD-(L)1 blockade; however, most patients observe little or no durable clinical benefit. To identify novel tumor-driven resistance mechanisms, we developed a panel of KP murine lung cancer models with intrinsic resistance to anti-PD-1 and queried differential gene expression between these tumors and anti-PD-1-sensitive tumors. We found that the enzyme autotaxin (ATX), and the metabolite it produces, lysophosphatidic acid (LPA), were significantly upregulated in resistant tumors and that ATX directly modulated antitumor immunity, with its expression negatively correlating with total and effector tumor-infiltrating CD8+ T cells. Pharmacological inhibition of ATX, or the downstream receptor LPAR5, in combination with anti-PD-1 was sufficient to restore the antitumor immune response and efficaciously control lung tumor growth in multiple KP tumor models. Additionally, ATX was significantly correlated with inflammatory gene signatures, including a CD8+ cytolytic score in multiple lung adenocarcinoma patient data sets, suggesting that an activated tumor-immune microenvironment upregulates ATX and thus provides an opportunity for cotargeting to prevent acquired resistance to anti-PD-1 treatment. These data reveal the ATX/LPA axis as an immunosuppressive pathway that diminishes the immune checkpoint blockade response in lung cancer.

Targeting immunosuppressive Ly6C+ classical monocytes reverses anti-PD-1/CTLA-4 immunotherapy resistance.

Introduction: Despite significant clinical advancement with the use of immune checkpoint blockade (ICB) in non-small cell lung cancer (NSCLC) there are still a major subset of patients that develop adaptive/acquired resistance. Understanding resistance mechanisms to ICB is critical to developing new therapeutic strategies and improving patient survival. The dynamic nature of the tumor microenvironment and the mutational load driving tumor immunogenicity limit the efficacy to ICB. Recent studies indicate that myeloid cells are drivers of ICB resistance. In this study we sought to understand which immune cells were contributing to resistance and if we could modify them in a way to improve response to ICB therapy. Results: Our results show that combination anti-PD-1/CTLA-4 produces an initial antitumor effect with evidence of an activated immune response. Upon extended treatment with anti-PD-1/CTLA-4 acquired resistance developed with an increase of the immunosuppressive populations, including T-regulatory cells, neutrophils and monocytes. Addition of anti-Ly6C blocking antibody to anti-PD-1/CTLA-4 was capable of completely reversing treatment resistance and restoring CD8 T cell activity in multiple KP lung cancer models and in the autochthonous lung cancer KrasLSL-G12D/p53fl/fl model. We found that there were higher classical Ly6C+ monocytes in anti-PD-1/CTLA-4 combination resistant tumors. B7 blockade illustrated the importance of dendritic cells for treatment efficacy of anti-Ly6C/PD-1/CTLA-4. We further determined that classical Ly6C+ monocytes in anti-PD-1/CTLA-4 resistant tumors are trafficked into the tumor via IFN-γ and the CCL2-CCR2 axis. Mechanistically we found that classical monocytes from ICB resistant tumors were unable to differentiate into antigen presenting cells and instead differentiated into immunosuppressive M2 macrophages or myeloid-derived suppressor cells (MDSC). Classical Ly6C+ monocytes from ICB resistant tumors had a decrease in both Flt3 and PU.1 expression that prevented differentiation into dendritic cells/macrophages. Conclusions: Therapeutically we found that addition of anti-Ly6C to the combination of anti-PD-1/CTLA-4 was capable of complete tumor eradication. Classical Ly6C+ monocytes differentiate into immunosuppressive cells, while blockade of classical monocytes drives dendritic cell differentiation/maturation to reinvigorate the anti-tumor T cell response. These findings support that immunotherapy resistance is associated with infiltrating monocytes and that controlling the differentiation process of monocytes can enhance the therapeutic potential of ICB.

CD70 is a therapeutic target upregulated in EMT-associated EGFR tyrosine kinase inhibitor resistance.

Effective therapeutic strategies are needed for non-small cell lung cancer (NSCLC) patients with epidermal growth factor receptor (EGFR) mutations that acquire resistance to EGFR tyrosine kinase inhibitors (TKIs) mediated by epithelial-to-mesenchymal transition (EMT). We investigate cell surface proteins that could be targeted by antibody-based or adoptive cell therapy approaches and identify CD70 as being highly upregulated in EMT-associated resistance. Moreover, CD70 upregulation is an early event in the evolution of resistance and occurs in drug-tolerant persister cells (DTPCs). CD70 promotes cell survival and invasiveness, and stimulation of CD70 triggers signal transduction pathways known to be re-activated with acquired TKI resistance. Anti-CD70 antibody drug conjugates (ADCs) and CD70-targeting chimeric antigen receptor (CAR) T cell and CAR NK cells show potent activity against EGFR TKI-resistant cells and DTPCs. These results identify CD70 as a therapeutic target for EGFR mutant tumors with acquired EGFR TKI resistance that merits clinical investigation.

Impad1 and Syt11 work in an epistatic pathway that regulates EMT-mediated vesicular trafficking to drive lung cancer invasion and metastasis.

Lung cancer is a highly aggressive and metastatic disease responsible for approximately 25% of all cancer-related deaths in the United States. Using high-throughput in vitro and in vivo screens, we have previously established Impad1 as a driver of lung cancer invasion and metastasis. Here we elucidate that Impad1 is a direct target of the epithelial microRNAs (miRNAs) miR-200 and miR∼96 and is de-repressed during epithelial-to-mesenchymal transition (EMT); thus, we establish a mode of regulation of the protein. Impad1 modulates Golgi apparatus morphology and vesicular trafficking through its interaction with a trafficking protein, Syt11. These changes in Golgi apparatus dynamics alter the extracellular matrix and the tumor microenvironment (TME) to promote invasion and metastasis. Inhibiting Impad1 or Syt11 disrupts the cancer cell secretome, regulates the TME, and reverses the invasive or metastatic phenotype. This work identifies Impad1 as a regulator of EMT and secretome-mediated changes during lung cancer progression.

Shared Nearest Neighbors Approach and Interactive Browser for Network Analysis of a Comprehensive Non-Small-Cell Lung Cancer Data Set.

PURPOSE: Advances in biological measurement technologies are enabling large-scale studies of patient cohorts across multiple omics platforms. Holistic analysis of these data can generate actionable insights for translational research and necessitate new approaches for data integration and mining. METHODS: We present a novel approach for integrating data across platforms on the basis of the shared nearest neighbors algorithm and use it to create a network of multiplatform data from the immunogenomic profiling of non-small-cell lung cancer project. RESULTS: Benchmarking demonstrates that the shared nearest neighbors-based network approach outperforms a traditional gene-gene network in capturing established interactions while providing new ones on the basis of the interplay between measurements from different platforms. When used to examine patient characteristics of interest, our approach provided signatures associated with and new leads related to recurrence and TP53 oncogenotype. CONCLUSION: The network developed offers an unprecedented, holistic view into immunogenomic profiling of non-small-cell lung cancer, which can be explored through the accompanying interactive browser that we built.

The microRNA-183/96/182 cluster inhibits lung cancer progression and metastasis by inducing an interleukin-2-mediated antitumor CD8+ cytotoxic T-cell response.

One of the mechanisms by which cancer cells acquire hyperinvasive and migratory properties with progressive loss of epithelial markers is the epithelial-to-mesenchymal transition (EMT). We have previously reported that in different cancer types, including nonsmall cell lung cancer (NSCLC), the microRNA-183/96/182 cluster (m96cl) is highly repressed in cells that have undergone EMT. In the present study, we used a novel conditional m96cl mouse to establish that loss of m96cl accelerated the growth of Kras mutant autochthonous lung adenocarcinomas. In contrast, ectopic expression of the m96cl in NSCLC cells results in a robust suppression of migration and invasion in vitro, and tumor growth and metastasis in vivo. Detailed immune profiling of the tumors revealed a significant enrichment of activated CD8+ cytotoxic T lymphocytes (CD8+ CTLs) in m96cl-expressing tumors, and m96cl-mediated suppression of tumor growth and metastasis was CD8+ CTL-dependent. Using coculture assays with naïve immune cells, we show that m96cl expression drives paracrine stimulation of CD8+ CTL proliferation and function. Using tumor microenvironment-associated gene expression profiling, we identified that m96cl elevates the interleukin-2 (IL2) signaling pathway and results in increased IL2-mediated paracrine stimulation of CD8+ CTLs. Furthermore, we identified that the m96cl modulates the expression of IL2 in cancer cells by regulating the expression of transcriptional repressors Foxf2 and Zeb1, and thereby alters the levels of secreted IL2 in the tumor microenvironment. Last, we show that in vivo depletion of IL2 abrogates m96cl-mediated activation of CD8+ CTLs and results in loss of metastatic suppression. Therefore, we have identified a novel mechanistic role of the m96cl in the suppression of lung cancer growth and metastasis by inducing an IL2-mediated systemic CD8+ CTL immune response.

Anti-tumor activity of cetuximab plus avelumab in non-small cell lung cancer patients involves innate immunity activation: findings from the CAVE-Lung trial.

BACKGROUND: We recently conducted Cetuximab-AVElumab-Lung (CAVE-Lung), a proof-of-concept, translational and clinical trial, to evaluate the combination of two IgG1 monoclonal antibodies (mAb): avelumab, an anti-PD-L1 drug, and cetuximab, an anti-epidermal growth factor receptor (EGFR) drug, as second- or third-line treatment in non-small cell lung cancer (NSCLC) patients. We have reported clinically relevant anti-tumor activity in 6/16 patients. Clinical benefit was accompanied by Natural Killer (NK) cell-mediated antibody-dependent cell cytotoxicity (ADCC). Among the 6 responding patients, 3 had progressed after initial response to a previous treatment with single agent anti-PD-1, nivolumab or pembrolizumab. METHODS: We report long-term clinical follow-up and additional findings on the anti-tumor activity and on the immune effects of cetuximab plus avelumab treatment for these 3 patients. RESULTS: As of November 30, 2021, 2/3 patients were alive. One patient was still on treatment from 34 months, while the other two patients had progression free survival (PFS) of 15 and 19 months, respectively. Analysis of serially collected peripheral blood mononuclear cells (PBMC) revealed long-term activation of NK cell-mediated ADCC. Comprehensive genomic profile analysis found somatic mutations and germline rare variants in DNA damage response (DDR) genes. Furthermore, by transcriptomic analysis of The Cancer Genome Atlas (TCGA) dataset we found that DDR mutant NSCLC displayed high STING pathway gene expression. In NSCLC patient-derived three-dimensional in vitro spheroid cultures, cetuximab plus avelumab treatment induced additive cancer cell growth inhibition as compared to single agent treatment. This effect was partially blocked by treatment with an anti-CD16 mAb, suggesting a direct involvement of NK cell activation. Furthermore, cetuximab plus avelumab treatment induced 10-, 20-, and 20-fold increase, respectively, in the gene expression of CCL5 and CXCL10, two STING downstream effector cytokines, and of interferon ß, as compared to untreated control samples. CONCLUSIONS: DDR mutations may contribute to DDR-induced STING pathway with sustained innate immunity activation following cetuximab plus avelumab combination in previously treated, PD-1 inhibitor responsive NSCLC patients.

Profiling of immune features to predict immunotherapy efficacy.

Immune checkpoint blockade (ICB) therapies exhibit substantial clinical benefit in different cancers, but relatively low response rates in the majority of patients highlight the need to understand mutual relationships among immune features. Here, we reveal overall positive correlations among immune checkpoints and immune cell populations. Clinically, patients benefiting from ICB exhibited increases for both immune stimulatory and inhibitory features after initiation of therapy, suggesting that the activation of the immune microenvironment might serve as the biomarker to predict immune response. As proof-of-concept, we demonstrated that the immune activation score (IS Δ) based on dynamic alteration of interleukins in patient plasma as early as two cycles (4-6 weeks) after starting immunotherapy can accurately predict immunotherapy efficacy. Our results reveal a systematic landscape of associations among immune features and provide a noninvasive, cost-effective, and time-efficient approach based on dynamic profiling of pre- and on-treatment plasma to predict immunotherapy efficacy.

Association of antibiotic treatment with immune-related adverse events in patients with cancer receiving immunotherapy.

BACKGROUND: To determine whether antibiotic treatment is a risk factor for immune-related adverse events (irAEs) across different patients with cancer receiving anti-PD-1/PD-L1 therapies. METHODS: The retrospective analysis includes clinical information from 767 patients with cancer treated at Hunan Cancer Hospital from 2017 to 2020. The pharmacovigilance data analysis includes individual cases of 38,705 safety reports from the US Food and Drug Administration Adverse Event Reporting System (FAERS) from 2014 to 2020, and 25,122 cases of safety reports from the World Health Organization database VigiBase from 2014 to 2019. All cases that received anti-PD-1/PD-L1 treatment were included. Multiomics data from patients across 25 cancer types were download from The Cancer Genome Atlas. Logistic regression and propensity score algorithm was employed to calculate OR of irAEs. RESULTS: Retrospective analysis of in-house patients showed that irAE potential risks are higher in all cancer (OR 2.12, 95% CI 1.38 to 3.22, false discovery rate (FDR) adjusted-p=1.93×10-3) and patients with lung cancer (OR 3.16, 95% CI 1.67 to 5.95, FDR adjusted-p=1.93×10-3) when using antibiotics. Potential risk of irAEs in patients with lung cancer with antibiotic treatment is significantly higher in FAERS (OR 1.39, 95% CI 1.21 to 1.59; FDR adjusted-p=1.62×10-5) and VigiBase (OR 1.32, 95% CI 1.09 to 1.59, FDR adjusted-p=0.05). Mechanistically, decreased microbial diversity caused by antibiotics use may increase the irAE risk through mediating the irAE-related factors. CONCLUSIONS: Our study is the first to comprehensively demonstrate the associations of irAEs and antibiotic during anti-PD-1/PD-L1 therapy across a wide spectrum of cancers by analyzing multisource data. Administration of antibiotics should be carefully evaluated in patients with cancer treated by anti-PD-1/PD-L1 to avoid potentially increasing irAE risk.

Enhanced Vulnerability of LKB1-Deficient NSCLC to Disruption of ATP Pools and Redox Homeostasis by 8-Cl-Ado.

Loss-of-function somatic mutations of STK11, a tumor suppressor gene encoding LKB1 that contributes to the altered metabolic phenotype of cancer cells, is the second most common event in lung adenocarcinomas and often co-occurs with activating KRAS mutations. Tumor cells lacking LKB1 display an aggressive phenotype, with uncontrolled cell growth and higher energetic and redox stress due to its failure to balance ATP and NADPH levels in response to cellular stimulus. The identification of effective therapeutic regimens for patients with LKB1-deficient non-small cell lung cancer (NSCLC) remains a major clinical need. Here, we report that LKB1-deficient NSCLC tumor cells displayed reduced basal levels of ATP and to a lesser extent other nucleotides, and markedly enhanced sensitivity to 8-Cl-adenosine (8-Cl-Ado), an energy-depleting nucleoside analog. Treatment with 8-Cl-Ado depleted intracellular ATP levels, raised redox stress, and induced cell death leading to a compensatory suppression of mTOR signaling in LKB1-intact, but not LKB1-deficient, cells. Proteomic analysis revealed that the MAPK/MEK/ERK and PI3K/AKT pathways were activated in response to 8-Cl-Ado treatment and targeting these pathways enhanced the antitumor efficacy of 8-Cl-Ado. IMPLICATIONS: Together, our findings demonstrate that LKB1-deficient tumor cells are selectively sensitive to 8-Cl-Ado and suggest that therapeutic approaches targeting vulnerable energy stores combined with signaling pathway inhibitors merit further investigation for this patient population.

GPEdit: the genetic and pharmacogenomic landscape of A-to-I RNA editing in cancers.

Altered A-to-I RNA editing has been widely observed in many human cancers and some editing sites are associated with drug sensitivity, implicating its therapeutic potential. Increasing evidence has demonstrated that a quantitative trait loci mapping approach is effective to understanding the genetic basis of RNA editing. We systematically performed RNA editing quantitative trait loci (edQTL) analysis in 33 human cancer types for >10 000 cancer samples and identified 320 029 edQTLs. We also identified 1688 ed-QTLs associated with patient overall survival and 4672 ed-QTLs associated with GWAS risk loci. Furthermore, we demonstrated the associations between RNA editing and >1000 anti-cancer drug response with ∼3.5 million significant associations. We developed GPEdit (https://hanlab.uth.edu/GPEdit/) to facilitate a global map of the genetic and pharmacogenomic landscape of RNA editing. GPEdit is a user-friendly and comprehensive database that provides an opportunity for a better understanding of the genetic impact and the effects on drug response of RNA editing in cancers.

Cold and heterogeneous T cell repertoire is associated with copy number aberrations and loss of immune genes in small-cell lung cancer.

Small-cell lung cancer (SCLC) is speculated to harbor complex genomic intratumor heterogeneity (ITH) associated with high recurrence rate and suboptimal response to immunotherapy. Here, using multi-region whole exome/T cell receptor (TCR) sequencing as well as immunohistochemistry, we reveal a rather homogeneous mutational landscape but extremely cold and heterogeneous TCR repertoire in limited-stage SCLC tumors (LS-SCLCs). Compared to localized non-small cell lung cancers, LS-SCLCs have similar predicted neoantigen burden and genomic ITH, but significantly colder and more heterogeneous TCR repertoire associated with higher chromosomal copy number aberration (CNA) burden. Furthermore, copy number loss of IFN-γ pathway genes is frequently observed and positively correlates with CNA burden. Higher mutational burden, higher T cell infiltration and positive PD-L1 expression are associated with longer overall survival (OS), while higher CNA burden is associated with shorter OS in patients with LS-SCLC.

The EMT activator ZEB1 accelerates endosomal trafficking to establish a polarity axis in lung adenocarcinoma cells.

Epithelial-to-mesenchymal transition (EMT) is a transcriptionally governed process by which cancer cells establish a front-rear polarity axis that facilitates motility and invasion. Dynamic assembly of focal adhesions and other actin-based cytoskeletal structures on the leading edge of motile cells requires precise spatial and temporal control of protein trafficking. Yet, the way in which EMT-activating transcriptional programs interface with vesicular trafficking networks that effect cell polarity change remains unclear. Here, by utilizing multiple approaches to assess vesicular transport dynamics through endocytic recycling and retrograde trafficking pathways in lung adenocarcinoma cells at distinct positions on the EMT spectrum, we find that the EMT-activating transcription factor ZEB1 accelerates endocytosis and intracellular trafficking of plasma membrane-bound proteins. ZEB1 drives turnover of the MET receptor tyrosine kinase by hastening receptor endocytosis and transport to the lysosomal compartment for degradation. ZEB1 relieves a plus-end-directed microtubule-dependent kinesin motor protein (KIF13A) and a clathrin-associated adaptor protein complex subunit (AP1S2) from microRNA-dependent silencing, thereby accelerating cargo transport through the endocytic recycling and retrograde vesicular pathways, respectively. Depletion of KIF13A or AP1S2 mitigates ZEB1-dependent focal adhesion dynamics, front-rear axis polarization, and cancer cell motility. Thus, ZEB1-dependent transcriptional networks govern vesicular trafficking dynamics to effect cell polarity change.

Targeting MYC-enhanced glycolysis for the treatment of small cell lung cancer.

INTRODUCTION: The transcription factor MYC is overexpressed in 30% of small cell lung cancer (SCLC) tumors and is known to modulate the balance between two major pathways of metabolism: glycolysis and mitochondrial respiration. This duality of MYC underscores the importance of further investigation into its role in SCLC metabolism and could lead to insights into metabolic targeting approaches. METHODS: We investigated differences in metabolic pathways in transcriptional and metabolomics datasets based on cMYC expression in patient and cell line samples. Metabolic pathway utilization was evaluated by flow cytometry and Seahorse extracellular flux methodology. Glycolysis inhibition was evaluated in vitro and in vivo using PFK158, a small molecular inhibitor of PFKFB3. RESULTS: MYC-overexpressing SCLC patient samples and cell lines exhibited increased glycolysis gene expression directly mediated by MYC. Further, MYC-overexpressing cell lines displayed enhanced glycolysis consistent with the Warburg effect, while cell lines with low MYC expression appeared more reliant on oxidative metabolism. Inhibition of glycolysis with PFK158 preferentially attenuated glucose uptake, ATP production, and lactate in MYC-overexpressing cell lines. Treatment with PFK158 in xenografts delayed tumor growth and decreased glycolysis gene expression. CONCLUSIONS: Our study highlights an in-depth characterization of SCLC metabolic programming and presents glycolysis as a targetable mechanism downstream of MYC that could offer therapeutic benefit in a subset of SCLC patients.

Structure-based classification predicts drug response in EGFR-mutant NSCLC.

Epidermal growth factor receptor (EGFR) mutations typically occur in exons 18-21 and are established driver mutations in non-small cell lung cancer (NSCLC)1-3. Targeted therapies are approved for patients with 'classical' mutations and a small number of other mutations4-6. However, effective therapies have not been identified for additional EGFR mutations. Furthermore, the frequency and effects of atypical EGFR mutations on drug sensitivity are unknown1,3,7-10. Here we characterize the mutational landscape in 16,715 patients with EGFR-mutant NSCLC, and establish the structure-function relationship of EGFR mutations on drug sensitivity. We found that EGFR mutations can be separated into four distinct subgroups on the basis of sensitivity and structural changes that retrospectively predict patient outcomes following treatment with EGFR inhibitors better than traditional exon-based groups. Together, these data delineate a structure-based approach for defining functional groups of EGFR mutations that can effectively guide treatment and clinical trial choices for patients with EGFR-mutant NSCLC and suggest that a structure-function-based approach may improve the prediction of drug sensitivity to targeted therapies in oncogenes with diverse mutations.

Targeting CDK4 overcomes EMT-mediated tumor heterogeneity and therapeutic resistance in KRAS-mutant lung cancer.

Lack of sustained response to therapeutic agents in patients with KRAS-mutant lung cancer poses a major challenge and arises partly due to intratumor heterogeneity that defines phenotypically distinct tumor subpopulations. To attain better therapeutic outcomes, it is important to understand the differential therapeutic sensitivities of tumor cell subsets. Epithelial-mesenchymal transition is a biological phenomenon that can alter the state of cells along a phenotypic spectrum and cause transcriptional rewiring to produce distinct tumor cell subpopulations. We utilized functional shRNA screens, in in vitro and in vivo models, to identify and validate an increased dependence of mesenchymal tumor cells on cyclin-dependent kinase 4 (CDK4) for survival, as well as a mechanism of resistance to MEK inhibitors. High zinc finger E-box binding homeobox 1 levels in mesenchymal tumor cells repressed p21, leading to perturbed CDK4 pathway activity. Increased dependence on CDK4 rendered mesenchymal cancer cells particularly vulnerable to selective CDK4 inhibitors. Coadministration of CDK4 and MEK inhibitors in heterogeneous tumors effectively targeted different tumor subpopulations, subverting the resistance to either single-agent treatment.