Mutation and cell state compatibility is required and targetable in Ph+ acute lymphoblastic leukemia minimal residual disease.

Efforts to cure BCR::ABL1 B cell acute lymphoblastic leukemia (Ph+ ALL) solely through inhibition of ABL1 kinase activity have thus far been insufficient despite the availability of tyrosine kinase inhibitors (TKIs) with broad activity against resistance mutants. The mechanisms that drive persistence within minimal residual disease (MRD) remain poorly understood and therefore untargeted. Utilizing 13 patient-derived xenograft (PDX) models and clinical trial specimens of Ph+ ALL, we examined how genetic and transcriptional features co-evolve to drive progression during prolonged TKI response. Our work reveals a landscape of cooperative mutational and transcriptional escape mechanisms that differ from those causing resistance to first generation TKIs. By analyzing MRD during remission, we show that the same resistance mutation can either increase or decrease cellular fitness depending on transcriptional state. We further demonstrate that directly targeting transcriptional state-associated vulnerabilities at MRD can overcome BCR::ABL1 independence, suggesting a new paradigm for rationally eradicating MRD prior to relapse. Finally, we illustrate how cell mass measurements of leukemia cells can be used to rapidly monitor dominant transcriptional features of Ph+ ALL to help rationally guide therapeutic selection from low-input samples.

ABL kinases regulate the stabilization of HIF-1α and MYC through CPSF1.

The hypoxia-inducible factor 1-α (HIF-1α) enables cells to adapt and respond to hypoxia (Hx), and the activity of this transcription factor is regulated by several oncogenic signals and cellular stressors. While the pathways controlling normoxic degradation of HIF-1α are well understood, the mechanisms supporting the sustained stabilization and activity of HIF-1α under Hx are less clear. We report that ABL kinase activity protects HIF-1α from proteasomal degradation during Hx. Using a fluorescence-activated cell sorting (FACS)-based CRISPR/Cas9 screen, we identified HIF-1α as a substrate of the cleavage and polyadenylation specificity factor-1 (CPSF1), an E3-ligase which targets HIF-1α for degradation in the presence of an ABL kinase inhibitor in Hx. We show that ABL kinases phosphorylate and interact with CUL4A, a cullin ring ligase adaptor, and compete with CPSF1 for CUL4A binding, leading to increased HIF-1α protein levels. Further, we identified the MYC proto-oncogene protein as a second CPSF1 substrate and show that active ABL kinase protects MYC from CPSF1-mediated degradation. These studies uncover a role for CPSF1 in cancer pathobiology as an E3-ligase antagonizing the expression of the oncogenic transcription factors, HIF-1α and MYC.

Compressed phenotypic screens for complex multicellular models and high-content assays.

High-throughput phenotypic screens leveraging biochemical perturbations, high-content readouts, and complex multicellular models could advance therapeutic discovery yet remain constrained by limitations of scale. To address this, we establish a method for compressing screens by pooling perturbations followed by computational deconvolution. Conducting controlled benchmarks with a highly bioactive small molecule library and a high-content imaging readout, we demonstrate increased efficiency for compressed experimental designs compared to conventional approaches. To prove generalizability, we apply compressed screening to examine transcriptional responses of patient-derived pancreatic cancer organoids to a library of tumor-microenvironment (TME)-nominated recombinant protein ligands. Using single-cell RNA-seq as a readout, we uncover reproducible phenotypic shifts induced by ligands that correlate with clinical features in larger datasets and are distinct from reference signatures available in public databases. In sum, our approach enables phenotypic screens that interrogate complex multicellular models with rich phenotypic readouts to advance translatable drug discovery as well as basic biology.

MCB-613 exploits a collateral sensitivity in drug resistant EGFR-mutant non-small cell lung cancer through covalent inhibition of KEAP1.

Targeted therapies have revolutionized cancer chemotherapy. Unfortunately, most patients develop multifocal resistance to these drugs within a matter of months. Here, we used a high-throughput phenotypic small molecule screen to identify MCB-613 as a compound that selectively targets EGFR-mutant, EGFR inhibitor-resistant non-small cell lung cancer (NSCLC) cells harboring diverse resistance mechanisms. Subsequent proteomic and functional genomic screens involving MCB-613 identified its target in this context to be KEAP1, revealing that this gene is selectively essential in the setting of EGFR inhibitor resistance. In-depth molecular characterization demonstrated that (1) MCB-613 binds KEAP1 covalently; (2) a single molecule of MCB-613 is capable of bridging two KEAP1 monomers together; and, (3) this modification interferes with the degradation of canonical KEAP1 substrates such as NRF2. Surprisingly, NRF2 knockout sensitizes cells to MCB-613, suggesting that the drug functions through modulation of an alternative KEAP1 substrate. Together, these findings advance MCB-613 as a new tool for exploiting the selective essentiality of KEAP1 in drug-resistant, EGFR-mutant NSCLC cells.

Spatially Resolved Single-Cell Assessment of Pancreatic Cancer Expression Subtypes Reveals Co-expressor Phenotypes and Extensive Intratumoral Heterogeneity.

Pancreatic ductal adenocarcinoma (PDAC) has been classified into classical and basal-like transcriptional subtypes by bulk RNA measurements. However, recent work has uncovered greater complexity to transcriptional subtypes than was initially appreciated using bulk RNA expression profiling. To provide a deeper understanding of PDAC subtypes, we developed a multiplex immunofluorescence (mIF) pipeline that quantifies protein expression of six PDAC subtype markers (CLDN18.2, TFF1, GATA6, KRT17, KRT5, and S100A2) and permits spatially resolved, single-cell interrogation of pancreatic tumors from resection specimens and core needle biopsies. Both primary and metastatic tumors displayed striking intratumoral subtype heterogeneity that was associated with patient outcomes, existed at the scale of individual glands, and was significantly reduced in patient-derived organoid cultures. Tumor cells co-expressing classical and basal markers were present in > 90% of tumors, existed on a basal-classical polarization continuum, and were enriched in tumors containing a greater admixture of basal and classical cell populations. Cell-cell neighbor analyses within tumor glands further suggested that co-expressor cells may represent an intermediate state between expression subtype poles. The extensive intratumoral heterogeneity identified through this clinically applicable mIF pipeline may inform prognosis and treatment selection for patients with PDAC. SIGNIFICANCE: A high-throughput pipeline using multiplex immunofluorescence in pancreatic cancer reveals striking expression subtype intratumoral heterogeneity with implications for therapy selection and identifies co-expressor cells that may serve as intermediates during subtype switching.

Deciphering the immunopeptidome in vivo reveals new tumour antigens.

Immunosurveillance of cancer requires the presentation of peptide antigens on major histocompatibility complex class I (MHC-I) molecules1-5. Current approaches to profiling of MHC-I-associated peptides, collectively known as the immunopeptidome, are limited to in vitro investigation or bulk tumour lysates, which limits our understanding of cancer-specific patterns of antigen presentation in vivo6. To overcome these limitations, we engineered an inducible affinity tag into the mouse MHC-I gene (H2-K1) and targeted this allele to the KrasLSL-G12D/+Trp53fl/fl mouse model (KP/KbStrep)7. This approach enabled us to precisely isolate MHC-I peptides from autochthonous pancreatic ductal adenocarcinoma and from lung adenocarcinoma (LUAD) in vivo. In addition, we profiled the LUAD immunopeptidome from the alveolar type 2 cell of origin up to late-stage disease. Differential peptide presentation in LUAD was not predictable by mRNA expression or translation efficiency and is probably driven by post-translational mechanisms. Vaccination with peptides presented by LUAD in vivo induced CD8+ T cell responses in naive mice and tumour-bearing mice. Many peptides specific to LUAD, including immunogenic peptides, exhibited minimal expression of the cognate mRNA, which prompts the reconsideration of antigen prediction pipelines that triage peptides according to transcript abundance8. Beyond cancer, the KbStrep allele is compatible with other Cre-driver lines to explore antigen presentation in vivo in the pursuit of understanding basic immunology, infectious disease and autoimmunity.

Microenvironment drives cell state, plasticity, and drug response in pancreatic cancer.

Prognostically relevant RNA expression states exist in pancreatic ductal adenocarcinoma (PDAC), but our understanding of their drivers, stability, and relationship to therapeutic response is limited. To examine these attributes systematically, we profiled metastatic biopsies and matched organoid models at single-cell resolution. In vivo, we identify a new intermediate PDAC transcriptional cell state and uncover distinct site- and state-specific tumor microenvironments (TMEs). Benchmarking models against this reference map, we reveal strong culture-specific biases in cancer cell transcriptional state representation driven by altered TME signals. We restore expression state heterogeneity by adding back in vivo-relevant factors and show plasticity in culture models. Further, we prove that non-genetic modulation of cell state can strongly influence drug responses, uncovering state-specific vulnerabilities. This work provides a broadly applicable framework for aligning cell states across in vivo and ex vivo settings, identifying drivers of transcriptional plasticity and manipulating cell state to target associated vulnerabilities.

Flow induces barrier and glycocalyx-related genes and negative surface charge in a lab-on-a-chip human blood-brain barrier model.

Microfluidic lab-on-a-chip (LOC) devices allow the study of blood-brain barrier (BBB) properties in dynamic conditions. We studied a BBB model, consisting of human endothelial cells derived from hematopoietic stem cells in co-culture with brain pericytes, in an LOC device to study fluid flow in the regulation of endothelial, BBB and glycocalyx-related genes and surface charge. The highly negatively charged endothelial surface glycocalyx functions as mechano-sensor detecting shear forces generated by blood flow on the luminal side of brain endothelial cells and contributes to the physical barrier of the BBB. Despite the importance of glycocalyx in the regulation of BBB permeability in physiological conditions and in diseases, the underlying mechanisms remained unclear. The MACE-seq gene expression profiling analysis showed differentially expressed endothelial, BBB and glycocalyx core protein genes after fluid flow, as well as enriched pathways for the extracellular matrix molecules. We observed increased barrier properties, a higher intensity glycocalyx staining and a more negative surface charge of human brain-like endothelial cells (BLECs) in dynamic conditions. Our work is the first study to provide data on BBB properties and glycocalyx of BLECs in an LOC device under dynamic conditions and confirms the importance of fluid flow for BBB culture models.

Polar Infection of Echovirus-30 Causes Differential Barrier Affection and Gene Regulation at the Blood-Cerebrospinal Fluid Barrier.

Echovirus-30 (E-30) is responsible for the extensive global outbreaks of meningitis in children. To gain access to the central nervous system, E-30 first has to cross the epithelial blood-cerebrospinal fluid barrier. Several meningitis causing bacteria preferentially infect human choroid plexus papilloma (HIBCPP) cells in a polar fashion from the basolateral cell side. Here, we investigated the polar infection of HIBCPP cells with E-30. Both apical and basolateral infections caused a significant decrease in the transepithelial electrical resistance of HIBCPP cells. However, to reach the same impact on the barrier properties, the multiplicity of infection of the apical side had to be higher than that of the basolateral infection. Furthermore, the number of infected cells at respective time-points after basolateral infection was significantly higher compared to apical infection. Cytotoxic effects of E-30 on HIBCPP cells during basolateral infection were observed following prolonged infection and appeared more drastically compared to the apical infection. Gene expression profiles determined by massive analysis of cDNA ends revealed distinct regulation of specific genes depending on the side of HIBCPP cells' infection. Altogether, our data highlights the polar effects of E-30 infection in a human in vitro model of the blood-cerebrospinal fluid barrier leading to central nervous system inflammation.

Contribution of brain pericytes in blood-brain barrier formation and maintenance: a transcriptomic study of cocultured human endothelial cells derived from hematopoietic stem cells.

Formation, maintenance, and repair of the blood-brain barrier (BBB) are critical for central nervous system homeostasis. The interaction of endothelial cells (ECs) with brain pericytes is known to induce BBB characteristics in brain ECs during embryogenesis and can be used to differentiate human ECs from stem cell source in in vitro BBB models. However, the molecular events involved in BBB maturation are not fully understood. To this end, human ECs derived from hematopoietic stem cells were cultivated with either primary bovine or cell line-derived human brain pericytes to induce BBB formation. Subsequently, the transcriptomic profiles of solocultured vs. cocultured ECs were analysed over time by Massive Analysis of cDNA Ends (MACE) technology. This RNA sequencing method is a 3'-end targeted, tag-based, reduced representation transcriptome profiling technique, that can reliably quantify all polyadenylated transcripts including those with low expression. By analysing the generated transcriptomic profiles, we can explore the molecular processes responsible for the functional changes observed in ECs in coculture with brain pericytes (e.g. barrier tightening, changes in the expression of transporters and receptors). Our results identified several up- and downregulated genes and signaling pathways that provide a valuable data source to further delineate complex molecular processes that are involved in BBB formation and BBB maintenance. In addition, this data provides a source to identify novel targets for central nervous system drug delivery strategies.

Using antagonistic pleiotropy to design a chemotherapy-induced evolutionary trap to target drug resistance in cancer.

Local adaptation directs populations towards environment-specific fitness maxima through acquisition of positively selected traits. However, rapid environmental changes can identify hidden fitness trade-offs that turn adaptation into maladaptation, resulting in evolutionary traps. Cancer, a disease that is prone to drug resistance, is in principle susceptible to such traps. We therefore performed pooled CRISPR-Cas9 knockout screens in acute myeloid leukemia (AML) cells treated with various chemotherapies to map the drug-dependent genetic basis of fitness trade-offs, a concept known as antagonistic pleiotropy (AP). We identified a PRC2-NSD2/3-mediated MYC regulatory axis as a drug-induced AP pathway whose ability to confer resistance to bromodomain inhibition and sensitivity to BCL-2 inhibition templates an evolutionary trap. Across diverse AML cell-line and patient-derived xenograft models, we find that acquisition of resistance to bromodomain inhibition through this pathway exposes coincident hypersensitivity to BCL-2 inhibition. Thus, drug-induced AP can be leveraged to design evolutionary traps that selectively target drug resistance in cancer.

Optofluidic real-time cell sorter for longitudinal CTC studies in mouse models of cancer.

Circulating tumor cells (CTCs) play a fundamental role in cancer progression. However, in mice, limited blood volume and the rarity of CTCs in the bloodstream preclude longitudinal, in-depth studies of these cells using existing liquid biopsy techniques. Here, we present an optofluidic system that continuously collects fluorescently labeled CTCs from a genetically engineered mouse model (GEMM) for several hours per day over multiple days or weeks. The system is based on a microfluidic cell sorting chip connected serially to an unanesthetized mouse via an implanted arteriovenous shunt. Pneumatically controlled microfluidic valves capture CTCs as they flow through the device, and CTC-depleted blood is returned back to the mouse via the shunt. To demonstrate the utility of our system, we profile CTCs isolated longitudinally from animals over 4 days of treatment with the BET inhibitor JQ1 using single-cell RNA sequencing (scRNA-Seq) and show that our approach eliminates potential biases driven by intermouse heterogeneity that can occur when CTCs are collected across different mice. The CTC isolation and sorting technology presented here provides a research tool to help reveal details of how CTCs evolve over time, allowing studies to credential changes in CTCs as biomarkers of drug response and facilitating future studies to understand the role of CTCs in metastasis.

Systematic mapping of BCL-2 gene dependencies in cancer reveals molecular determinants of BH3 mimetic sensitivity.

While inhibitors of BCL-2 family proteins (BH3 mimetics) have shown promise as anti-cancer agents, the various dependencies or co-dependencies of diverse cancers on BCL-2 genes remain poorly understood. Here we develop a drug screening approach to define the sensitivity of cancer cells from ten tissue types to all possible combinations of selective BCL-2, BCL-XL, and MCL-1 inhibitors and discover that most cell lines depend on at least one combination for survival. We demonstrate that expression levels of BCL-2 genes predict single mimetic sensitivity, whereas EMT status predicts synergistic dependence on BCL-XL+MCL-1. Lastly, we use a CRISPR/Cas9 screen to discover that BFL-1 and BCL-w promote resistance to all tested combinations of BCL-2, BCL-XL, and MCL-1 inhibitors. Together, these results provide a roadmap for rationally targeting BCL-2 family dependencies in diverse human cancers and motivate the development of selective BFL-1 and BCL-w inhibitors to overcome intrinsic resistance to BH3 mimetics.

Mapping Effector-Phenotype Landscapes in KRAS-Driven Cancers.

Oncogenic KRAS can activate numerous effector pathways to drive malignant progression. However, the relationships between specific effectors and oncogenic phenotypes, and the extent to which these relationships vary across heterogeneous tumors, are incompletely understood. Recently in Cell Reports, a team of scientists described an innovative, combinatorial siRNA-based approach to functionally link KRAS effectors and phenotypes in a large panel of cancer cell lines. Central to this work was the identification of two major subtypes of KRAS-mutant cancers with distinct effector landscapes and tractable therapeutic vulnerabilities.

Evaluation of epidermal growth factor receptor signaling effects in gastric cancer cell lines by detailed motility-focused phenotypic characterization linked with molecular analysis.

BACKGROUND: Gastric cancers frequently overexpress the epidermal growth factor receptor (EGFR), which has been implicated in pathological processes including tumor cell motility, invasion and metastasis. Targeting EGFR with the inhibitory antibody cetuximab may affect the motile and invasive behavior of tumor cells. Here, we evaluated the effects of EGFR signaling in gastric cancer cell lines to link the phenotypic behavior of the cells with their molecular characteristics. METHODS: Phenotypic effects were analyzed in four gastric cancer cell lines (AGS, Hs746T, LMSU and MKN1) by time-lapse microscopy and transwell invasion assay. Effects on EGFR signaling were detected using Western blot and proteome profiler analyses. A network was constructed linking EGFR signaling to the regulation of cellular motility. RESULTS: The analysis of the effects of treatment with epidermal growth factor (EGF) and cetuximab revealed that only one cell line (MKN1) was sensitive to cetuximab treatment in all phenotypic assays, whereas the other cell lines were either not responsive (Hs746T, LMSU) or sensitive only in certain tests (AGS). Cetuximab inhibited EGFR, MAPK and AKT activity and associated components of the EGFR signaling pathway to different degrees in cetuximab-sensitive MKN1 cells. In contrast, no such changes were observed in Hs746T cells. Thus, the different phenotypic behaviors of the cells were linked to their molecular response to treatment. Genetic alterations had different associations with response to treatment: while PIK3CA mutations and KRAS mutation or amplification were not obstructive, the MET mutation was associated with non-response. CONCLUSION: These results identify components of the EGFR signaling network as important regulators of the phenotypic and molecular response to cetuximab treatment.

Transforming Growth Factor-ß Drives the Transendothelial Migration of Hepatocellular Carcinoma Cells.

The entry of malignant hepatocytes into blood vessels is a key step in the dissemination and metastasis of hepatocellular carcinoma (HCC). The identification of molecular mechanisms involved in the transmigration of malignant hepatocytes through the endothelial barrier is of high relevance for therapeutic intervention and metastasis prevention. In this study, we employed a model of hepatocellular transmigration that mimics vascular invasion using hepatic sinusoidal endothelial cells and malignant hepatocytes evincing a mesenchymal-like, invasive phenotype by transforming growth factor (TGF)-ß. Labelling of respective cell populations with various stable isotopes and subsequent mass spectrometry analyses allowed the "real-time" detection of molecular changes in both transmigrating hepatocytes and endothelial cells. Interestingly, the proteome profiling revealed 36 and 559 regulated proteins in hepatocytes and endothelial cells, respectively, indicating significant changes during active transmigration that mostly depends on cell-cell interaction rather than on TGF-ß alone. Importantly, matching these in vitro findings with HCC patient data revealed a panel of common molecular alterations including peroxiredoxin-3, epoxide hydrolase, transgelin-2 and collectin 12 that are clinically relevant for the patient's survival. We conclude that hepatocellular plasticity induced by TGF-ß is crucially involved in blood vessel invasion of HCC cells.

A Landscape of Therapeutic Cooperativity in KRAS Mutant Cancers Reveals Principles for Controlling Tumor Evolution.

Combinatorial inhibition of effector and feedback pathways is a promising treatment strategy for KRAS mutant cancers. However, the particular pathways that should be targeted to optimize therapeutic responses are unclear. Using CRISPR/Cas9, we systematically mapped the pathways whose inhibition cooperates with drugs targeting the KRAS effectors MEK, ERK, and PI3K. By performing 70 screens in models of KRAS mutant colorectal, lung, ovarian, and pancreas cancers, we uncovered universal and tissue-specific sensitizing combinations involving inhibitors of cell cycle, metabolism, growth signaling, chromatin regulation, and transcription. Furthermore, these screens revealed secondary genetic modifiers of sensitivity, yielding a SRC inhibitor-based combination therapy for KRAS/PIK3CA double-mutant colorectal cancers (CRCs) with clinical potential. Surprisingly, acquired resistance to combinations of growth signaling pathway inhibitors develops rapidly following treatment, but by targeting signaling feedback or apoptotic priming, it is possible to construct three-drug combinations that greatly delay its emergence.

Precision medicine for hepatocelluar carcinoma using molecular pattern diagnostics: results from a preclinical pilot study.

The aim of this study was to design a road map for personalizing cancer therapy in hepatocellular carcinoma (HCC) by using molecular pattern diagnostics. As an exploratory study, we investigated molecular patterns of tissues of two tumors from individual HCC patients, which in previous experiments had shown contrasting reactions to the phase 2 transforming growth factor beta receptor 1 inhibitor galunisertib. Cancer-driving molecular patterns encompass - inter alias - altered transcription profiles and somatic mutations in coding regions differentiating tumors from their respective peritumoral tissues and from each other. Massive analysis of cDNA ends and all-exome sequencing demonstrate a highly divergent transcriptional and mutational landscape, respectively, for the two tumors, that offers potential explanations for the tumors contrasting responses to galunisertib. Molecular pattern diagnostics (MPDs) suggest alternative, individual-tumor-specific therapies, which in both cases deviate from the standard sorafenib treatment and from each other. Suggested personalized therapies use kinase inhibitors and immune-focused drugs as well as low-toxicity natural compounds identified using an advanced bioinformatics routine included in the MPD protocol. The MPD pipeline we describe here for the prediction of suitable drugs for treatment of two contrasting HCCs may serve as a blueprint for the design of therapies for various types of cancer.

Epstein-Barr virus ensures B cell survival by uniquely modulating apoptosis at early and late times after infection.

Latent Epstein-Barr virus (EBV) infection is causally linked to several human cancers. EBV expresses viral oncogenes that promote cell growth and inhibit the apoptotic response to uncontrolled proliferation. The EBV oncoprotein LMP1 constitutively activates NFκB and is critical for survival of EBV-immortalized B cells. However, during early infection EBV induces rapid B cell proliferation with low levels of LMP1 and little apoptosis. Therefore, we sought to define the mechanism of survival in the absence of LMP1/NFκB early after infection. We used BH3 profiling to query mitochondrial regulation of apoptosis and defined a transition from uninfected B cells (BCL-2) to early-infected (MCL-1/BCL-2) and immortalized cells (BFL-1). This dynamic change in B cell survival mechanisms is unique to virus-infected cells and relies on regulation of MCL-1 mitochondrial localization and BFL-1 transcription by the viral EBNA3A protein. This study defines a new role for EBNA3A in the suppression of apoptosis with implications for EBV lymphomagenesis.