Epigenetic reprogramming shapes the cellular landscape of schwannoma.

Mechanisms specifying cancer cell states and response to therapy are incompletely understood. Here we show epigenetic reprogramming shapes the cellular landscape of schwannomas, the most common tumors of the peripheral nervous system. We find schwannomas are comprised of 2 molecular groups that are distinguished by activation of neural crest or nerve injury pathways that specify tumor cell states and the architecture of the tumor immune microenvironment. Moreover, we find radiotherapy is sufficient for interconversion of neural crest schwannomas to immune-enriched schwannomas through epigenetic and metabolic reprogramming. To define mechanisms underlying schwannoma groups, we develop a technique for simultaneous interrogation of chromatin accessibility and gene expression coupled with genetic and therapeutic perturbations in single-nuclei. Our results elucidate a framework for understanding epigenetic drivers of tumor evolution and establish a paradigm of epigenetic and metabolic reprograming of cancer cells that shapes the immune microenvironment in response to radiotherapy.

ENPP1 is an innate immune checkpoint of the anticancer cGAMP-STING pathway in breast cancer.

Ectonucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1) expression correlates with poor prognosis in many cancers, and we previously discovered that ENPP1 is the dominant hydrolase of extracellular cGAMP: a cancer-cell-produced immunotransmitter that activates the anticancer stimulator of interferon genes (STING) pathway. However, ENPP1 has other catalytic activities and the molecular and cellular mechanisms contributing to its tumorigenic effects remain unclear. Here, using single-cell RNA-seq, we show that ENPP1 in both cancer and normal tissues drives primary breast tumor growth and metastasis by dampening extracellular 2'3'-cyclic-GMP-AMP (cGAMP)-STING-mediated antitumoral immunity. ENPP1 loss-of-function in both cancer cells and normal tissues slowed primary tumor growth and abolished metastasis. Selectively abolishing the cGAMP hydrolysis activity of ENPP1 phenocopied ENPP1 knockout in a STING-dependent manner, demonstrating that restoration of paracrine cGAMP-STING signaling is the dominant anti-cancer mechanism of ENPP1 inhibition. Finally, ENPP1 expression in breast tumors deterministically predicated whether patients would remain free of distant metastasis after pembrolizumab (anti-PD-1) treatment followed by surgery. Altogether, ENPP1 blockade represents a strategy to exploit cancer-produced extracellular cGAMP for controlled local activation of STING and is therefore a promising therapeutic approach against breast cancer.

Integrative identification of non-coding regulatory regions driving metastatic prostate cancer.

Large-scale sequencing efforts of thousands of tumor samples have been undertaken to understand the mutational landscape of the coding genome. However, the vast majority of germline and somatic variants occur within non-coding portions of the genome. These genomic regions do not directly encode for specific proteins, but can play key roles in cancer progression, for example by driving aberrant gene expression control. Here, we designed an integrative computational and experimental framework to identify recurrently mutated non-coding regulatory regions that drive tumor progression. Application of this approach to whole-genome sequencing (WGS) data from a large cohort of metastatic castration-resistant prostate cancer (mCRPC) revealed a large set of recurrently mutated regions. We used (i) in silico prioritization of functional non-coding mutations, (ii) massively parallel reporter assays, and (iii) in vivo CRISPR-interference (CRISPRi) screens in xenografted mice to systematically identify and validate driver regulatory regions that drive mCRPC. We discovered that one of these enhancer regions, GH22I030351, acts on a bidirectional promoter to simultaneously modulate expression of U2-associated splicing factor SF3A1 and chromosomal protein CCDC157. We found that both SF3A1 and CCDC157 are promoters of tumor growth in xenograft models of prostate cancer. We nominated a number of transcription factors, including SOX6, to be responsible for higher expression of SF3A1 and CCDC157. Collectively, we have established and confirmed an integrative computational and experimental approach that enables the systematic detection of non-coding regulatory regions that drive the progression of human cancers.

ENPP1 is an innate immune checkpoint of the anticancer cGAMP-STING pathway.

ENPP1 expression correlates with poor prognosis in many cancers, and we previously discovered that ENPP1 is the dominant hydrolase of extracellular cGAMP: a cancer-cell-produced immunotransmitter that activates the anticancer STING pathway. However, ENPP1 has other catalytic activities and the molecular and cellular mechanisms contributing to its tumorigenic effects remain unclear. Here, using single cell RNA-seq (scRNA-seq), we show that ENPP1 overexpression drives primary breast tumor growth and metastasis by synergistically dampening extracellular cGAMP-STING mediated antitumoral immunity and activating immunosuppressive extracellular adenosine (eADO) signaling. In addition to cancer cells, stromal and immune cells in the tumor microenvironment (TME) also express ENPP1 that restrains their response to tumor-derived cGAMP. Enpp1 loss-of-function in both cancer cells and normal tissues slowed primary tumor initiation and growth and prevented metastasis in an extracellular cGAMP- and STING-dependent manner. Selectively abolishing the cGAMP hydrolysis activity of ENPP1 phenocopied total ENPP1 knockout, demonstrating that restoration of paracrine cGAMP-STING signaling is the dominant anti-cancer mechanism of ENPP1 inhibition. Strikingly, we find that breast cancer patients with low ENPP1 expression have significantly higher immune infiltration and improved response to therapeutics impacting cancer immunity upstream or downstream of the cGAMP-STING pathway, like PARP inhibitors and anti-PD1. Altogether, selective inhibition of ENPP1's cGAMP hydrolase activity alleviates an innate immune checkpoint to boost cancer immunity and is therefore a promising therapeutic approach against breast cancer that may synergize with other cancer immunotherapies.

IFITM proteins assist cellular uptake of diverse linked chemotypes.

The search for cell-permeable drugs has conventionally focused on low-molecular weight (MW), nonpolar, rigid chemical structures. However, emerging therapeutic strategies break traditional drug design rules by employing flexibly linked chemical entities composed of more than one ligand. Using complementary genome-scale chemical-genetic approaches we identified an endogenous chemical uptake pathway involving interferon-induced transmembrane proteins (IFITMs) that modulates the cell permeability of a prototypical biopic inhibitor of MTOR (RapaLink-1, MW: 1784 g/mol). We devised additional linked inhibitors targeting BCR-ABL1 (DasatiLink-1, MW: 1518 g/mol) and EIF4A1 (BisRoc-1, MW: 1466 g/mol), uptake of which was facilitated by IFITMs. We also found that IFITMs moderately assisted some proteolysis-targeting chimeras and examined the physicochemical requirements for involvement of this uptake pathway.

Highly efficient generation of isogenic pluripotent stem cell models using prime editing.

The recent development of prime editing (PE) genome engineering technologies has the potential to significantly simplify the generation of human pluripotent stem cell (hPSC)-based disease models. PE is a multicomponent editing system that uses a Cas9-nickase fused to a reverse transcriptase (nCas9-RT) and an extended PE guide RNA (pegRNA). Once reverse transcribed, the pegRNA extension functions as a repair template to introduce precise designer mutations at the target site. Here, we systematically compared the editing efficiencies of PE to conventional gene editing methods in hPSCs. This analysis revealed that PE is overall more efficient and precise than homology-directed repair of site-specific nuclease-induced double-strand breaks. Specifically, PE is more effective in generating heterozygous editing events to create autosomal dominant disease-associated mutations. By stably integrating the nCas9-RT into hPSCs we achieved editing efficiencies equal to those reported for cancer cells, suggesting that the expression of the PE components, rather than cell-intrinsic features, limit PE in hPSCs. To improve the efficiency of PE in hPSCs, we optimized the delivery modalities for the PE components. Delivery of the nCas9-RT as mRNA combined with synthetically generated, chemically-modified pegRNAs and nicking guide RNAs improved editing efficiencies up to 13-fold compared with transfecting the PE components as plasmids or ribonucleoprotein particles. Finally, we demonstrated that this mRNA-based delivery approach can be used repeatedly to yield editing efficiencies exceeding 60% and to correct or introduce familial mutations causing Parkinson's disease in hPSCs.

From muscles to nerves, our body is formed of many kinds of cells which can each respond slightly differently to the same harmful genetic changes. Understanding the exact relationship between mutations and cell-type specific function is essential to better grasp how conditions such as Parkinson's disease or amyotrophic lateral sclerosis progress and can be treated. Stem cells could be an important tool in that effort, as they can be directed to mature into many cell types in the laboratory. Yet it remains difficult to precisely introduce disease-relevant mutations in these cells. To remove this obstacle, Li et al. focused on prime editing, a cutting-edge 'search and replace' approach which can introduce new genetic information into a specific DNA sequence. However, it was unclear whether this technique could be used to efficiently create stem cell models of human diseases. A first set of experiments showed that prime editing is superior to conventional approaches when generating mutated genes in stem cells. Li et al. then further improved the efficiency and precision of the method by tweaking how prime editing components are delivered into the cells. The refined approach could be harnessed to quickly generate large numbers of stem cells carrying mutations associated with Parkinson's disease; crucially, prime editing could then also be used to revert a mutated gene back to its healthy form. The improved prime editing approach developed by Li et al. removes a major hurdle for scientists hoping to use stem cells to study genetic diseases. This could potentially help to unlock progress in how we understand and ultimately treat these conditions.

A campaign targeting a conserved Hsp70 binding site uncovers how subcellular localization is linked to distinct biological activities.

The potential of small molecules to localize within subcellular compartments is rarely explored. To probe this question, we measured the localization of Hsp70 inhibitors using fluorescence microscopy. We found that even closely related analogs had dramatically different distributions, with some residing predominantly in the mitochondria and others in the ER. CRISPRi screens supported this idea, showing that different compounds had distinct chemogenetic interactions with Hsp70s of the ER (HSPA5/BiP) and mitochondria (HSPA9/mortalin) and their co-chaperones. Moreover, localization seemed to determine function, even for molecules with conserved binding sites. Compounds with distinct partitioning have distinct anti-proliferative activity in breast cancer cells compared with anti-viral activity in cellular models of Dengue virus replication, likely because different sets of Hsp70s are required in these processes. These findings highlight the contributions of subcellular partitioning and chemogenetic interactions to small molecule activity, features that are rarely explored during medicinal chemistry campaigns.

A global cancer data integrator reveals principles of synthetic lethality, sex disparity and immunotherapy.

BACKGROUND: Advances in cancer biology are increasingly dependent on integration of heterogeneous datasets. Large-scale efforts have systematically mapped many aspects of cancer cell biology; however, it remains challenging for individual scientists to effectively integrate and understand this data. RESULTS: We have developed a new data retrieval and indexing framework that allows us to integrate publicly available data from different sources and to combine publicly available data with new or bespoke datasets. Our approach, which we have named the cancer data integrator (CanDI), is straightforward to implement, is well documented, and is continuously updated which should enable individual users to take full advantage of efforts to map cancer cell biology. We show that CanDI empowered testable hypotheses of new synthetic lethal gene pairs, genes associated with sex disparity, and immunotherapy targets in cancer. CONCLUSIONS: CanDI provides a flexible approach for large-scale data integration in cancer research enabling rapid generation of hypotheses. The CanDI data integrator is available at https://github.com/GilbertLabUCSF/CanDI .

An integrated functional and clinical genomics approach reveals genes driving aggressive metastatic prostate cancer.

Genomic sequencing of thousands of tumors has revealed many genes associated with specific types of cancer. Similarly, large scale CRISPR functional genomics efforts have mapped genes required for cancer cell proliferation or survival in hundreds of cell lines. Despite this, for specific disease subtypes, such as metastatic prostate cancer, there are likely a number of undiscovered tumor specific driver genes that may represent potential drug targets. To identify such genetic dependencies, we performed genome-scale CRISPRi screens in metastatic prostate cancer models. We then created a pipeline in which we integrated pan-cancer functional genomics data with our metastatic prostate cancer functional and clinical genomics data to identify genes that can drive aggressive prostate cancer phenotypes. Our integrative analysis of these data reveals known prostate cancer specific driver genes, such as AR and HOXB13, as well as a number of top hits that are poorly characterized. In this study we highlight the strength of an integrated clinical and functional genomics pipeline and focus on two top hit genes, KIF4A and WDR62. We demonstrate that both KIF4A and WDR62 drive aggressive prostate cancer phenotypes in vitro and in vivo in multiple models, irrespective of AR-status, and are also associated with poor patient outcome.

Revealing molecular pathways for cancer cell fitness through a genetic screen of the cancer translatome.

The major cap-binding protein eukaryotic translation initiation factor 4E (eIF4E), an ancient protein required for translation of all eukaryotic genomes, is a surprising yet potent oncogenic driver. The genetic interactions that maintain the oncogenic activity of this key translation factor remain unknown. In this study, we carry out a genome-wide CRISPRi screen wherein we identify more than 600 genetic interactions that sustain eIF4E oncogenic activity. Our data show that eIF4E controls the translation of Tfeb, a key executer of the autophagy response. This autophagy survival response is triggered by mitochondrial proteotoxic stress, which allows cancer cell survival. Our screen also reveals a functional interaction between eIF4E and a single anti-apoptotic factor, Bcl-xL, in tumor growth. Furthermore, we show that eIF4E and the exon-junction complex (EJC), which is involved in many steps of RNA metabolism, interact to control the migratory properties of cancer cells. Overall, we uncover several cancer-specific vulnerabilities that provide further resolution of the cancer translatome.

Mapping cancer genetics at single-cell resolution.

CRISPR genetic interaction mapping reveals networks of genes that drive cancer phenotypes.

Pharmaceutical-Grade Rigosertib Is a Microtubule-Destabilizing Agent.

We recently used CRISPRi/a-based chemical-genetic screens and cell biological, biochemical, and structural assays to determine that rigosertib, an anti-cancer agent in phase III clinical trials, kills cancer cells by destabilizing microtubules. Reddy and co-workers (Baker et al., 2020, this issue of Molecular Cell) suggest that a contaminating degradation product in commercial formulations of rigosertib is responsible for the microtubule-destabilizing activity. Here, we demonstrate that cells treated with pharmaceutical-grade rigosertib (>99.9% purity) or commercially obtained rigosertib have qualitatively indistinguishable phenotypes across multiple assays. The two formulations have indistinguishable chemical-genetic interactions with genes that modulate microtubule stability, both destabilize microtubules in cells and in vitro, and expression of a rationally designed tubulin mutant with a mutation in the rigosertib binding site (L240F TUBB) allows cells to proliferate in the presence of either formulation. Importantly, the specificity of the L240F TUBB mutant for microtubule-destabilizing agents has been confirmed independently. Thus, rigosertib kills cancer cells by destabilizing microtubules, in agreement with our original findings.

Clonal ZEB1-Driven Mesenchymal Transition Promotes Targetable Oncologic Antiangiogenic Therapy Resistance.

Glioblastoma (GBM) responses to bevacizumab are invariably transient with acquired resistance. We profiled paired patient specimens and bevacizumab-resistant xenograft models pre- and post-resistance toward the primary goal of identifying regulators whose targeting could prolong the therapeutic window, and the secondary goal of identifying biomarkers of therapeutic window closure. Bevacizumab-resistant patient specimens and xenografts exhibited decreased vessel density and increased hypoxia versus pre-resistance, suggesting that resistance occurs despite effective therapeutic devascularization. Microarray analysis revealed upregulated mesenchymal genes in resistant tumors correlating with bevacizumab treatment duration and causing three changes enabling resistant tumor growth in hypoxia. First, perivascular invasiveness along remaining blood vessels, which co-opts vessels in a VEGF-independent and neoangiogenesis-independent manner, was upregulated in novel biomimetic 3D bioengineered platforms modeling the bevacizumab-resistant microenvironment. Second, tumor-initiating stem cells housed in the perivascular niche close to remaining blood vessels were enriched. Third, metabolic reprogramming assessed through real-time bioenergetic measurement and metabolomics upregulated glycolysis and suppressed oxidative phosphorylation. Single-cell sequencing of bevacizumab-resistant patient GBMs confirmed upregulated mesenchymal genes, particularly glycoprotein YKL-40 and transcription factor ZEB1, in later clones, implicating these changes as treatment-induced. Serum YKL-40 was elevated in bevacizumab-resistant versus bevacizumab-naïve patients. CRISPR and pharmacologic targeting of ZEB1 with honokiol reversed the mesenchymal gene expression and associated stem cell, invasion, and metabolic changes defining resistance. Honokiol caused greater cell death in bevacizumab-resistant than bevacizumab-responsive tumor cells, with surviving cells losing mesenchymal morphology. Employing YKL-40 as a resistance biomarker and ZEB1 as a target to prevent resistance could fulfill the promise of antiangiogenic therapy. SIGNIFICANCE: Bevacizumab resistance in GBM is associated with mesenchymal/glycolytic shifts involving YKL-40 and ZEB1. Targeting ZEB1 reduces bevacizumab-resistant GBM phenotypes. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/7/1498/F1.large.jpg.

Exploring genetic interaction manifolds constructed from rich single-cell phenotypes.

How cellular and organismal complexity emerges from combinatorial expression of genes is a central question in biology. High-content phenotyping approaches such as Perturb-seq (single-cell RNA-sequencing pooled CRISPR screens) present an opportunity for exploring such genetic interactions (GIs) at scale. Here, we present an analytical framework for interpreting high-dimensional landscapes of cell states (manifolds) constructed from transcriptional phenotypes. We applied this approach to Perturb-seq profiling of strong GIs mined from a growth-based, gain-of-function GI map. Exploration of this manifold enabled ordering of regulatory pathways, principled classification of GIs (e.g., identifying suppressors), and mechanistic elucidation of synergistic interactions, including an unexpected synergy between CBL and CNN1 driving erythroid differentiation. Finally, we applied recommender system machine learning to predict interactions, facilitating exploration of vastly larger GI manifolds.

DNA-Dependent Protein Kinase Drives Prostate Cancer Progression through Transcriptional Regulation of the Wnt Signaling Pathway.

PURPOSE: Protein kinases are known to play a prominent role in oncogenic progression across multiple cancer subtypes, yet their role in prostate cancer progression remains underexplored. The purpose of this study was to identify kinases that drive prostate cancer progression.Experimental Design: To discover kinases that drive prostate cancer progression, we investigated the association between gene expression of all known kinases and long-term clinical outcomes in tumor samples from 545 patients with high-risk disease. We evaluated the impact of genetic and pharmacologic inhibition of the most significant kinase associated with metastatic progression in vitro and in vivo. RESULTS: DNA-dependent protein kinase (DNAPK) was identified as the most significant kinase associated with metastatic progression in high-risk prostate cancer. Inhibition of DNAPK suppressed the growth of both AR-dependent and AR-independent prostate cancer cells. Gene set enrichment analysis nominated Wnt as the top pathway associated with DNAPK. We found that DNAPK interacts with the Wnt transcription factor LEF1 and is critical for LEF1-mediated transcription. CONCLUSIONS: Our data show that DNAPK drives prostate cancer progression through transcriptional regulation of Wnt signaling and is an attractive therapeutic target in aggressive prostate cancer.

KRASG12C inhibition produces a driver-limited state revealing collateral dependencies.

Inhibitors targeting KRASG12C, a mutant form of the guanosine triphosphatase (GTPase) KRAS, are a promising new class of oncogene-specific therapeutics for the treatment of tumors driven by the mutant protein. These inhibitors react with the mutant cysteine residue by binding covalently to the switch-II pocket (S-IIP) that is present only in the inactive guanosine diphosphate (GDP)-bound form of KRASG12C, sparing the wild-type protein. We used a genome-scale CRISPR interference (CRISPRi) functional genomics platform to systematically identify genetic interactions with a KRASG12C inhibitor in cellular models of KRASG12C mutant lung and pancreatic cancer. Our data revealed genes that were selectively essential in this oncogenic driver-limited cell state, meaning that their loss enhanced cellular susceptibility to direct KRASG12C inhibition. We termed such genes "collateral dependencies" (CDs) and identified two classes of combination therapies targeting these CDs that increased KRASG12C target engagement or blocked residual survival pathways in cells and in vivo. From our findings, we propose a framework for assessing genetic dependencies induced by oncogene inhibition.