An expanded universe of cancer targets.

The characterization of cancer genomes has provided insight into somatically altered genes across tumors, transformed our understanding of cancer biology, and enabled tailoring of therapeutic strategies. However, the function of most cancer alleles remains mysterious, and many cancer features transcend their genomes. Consequently, tumor genomic characterization does not influence therapy for most patients. Approaches to understand the function and circuitry of cancer genes provide complementary approaches to elucidate both oncogene and non-oncogene dependencies. Emerging work indicates that the diversity of therapeutic targets engendered by non-oncogene dependencies is much larger than the list of recurrently mutated genes. Here we describe a framework for this expanded list of cancer targets, providing novel opportunities for clinical translation.

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.

A Somatically Acquired Enhancer of the Androgen Receptor Is a Noncoding Driver in Advanced Prostate Cancer.

Increased androgen receptor (AR) activity drives therapeutic resistance in advanced prostate cancer. The most common resistance mechanism is amplification of this locus presumably targeting the AR gene. Here, we identify and characterize a somatically acquired AR enhancer located 650 kb centromeric to the AR. Systematic perturbation of this enhancer using genome editing decreased proliferation by suppressing AR levels. Insertion of an additional copy of this region sufficed to increase proliferation under low androgen conditions and to decrease sensitivity to enzalutamide. Epigenetic data generated in localized prostate tumors and benign specimens support the notion that this region is a developmental enhancer. Collectively, these observations underscore the importance of epigenomic profiling in primary specimens and the value of deploying genome editing to functionally characterize noncoding elements. More broadly, this work identifies a therapeutic vulnerability for targeting the AR and emphasizes the importance of regulatory elements as highly recurrent oncogenic drivers.

Organoid Modeling of the Tumor Immune Microenvironment.

In vitro cancer cultures, including three-dimensional organoids, typically contain exclusively neoplastic epithelium but require artificial reconstitution to recapitulate the tumor microenvironment (TME). The co-culture of primary tumor epithelia with endogenous, syngeneic tumor-infiltrating lymphocytes (TILs) as a cohesive unit has been particularly elusive. Here, an air-liquid interface (ALI) method propagated patient-derived organoids (PDOs) from >100 human biopsies or mouse tumors in syngeneic immunocompetent hosts as tumor epithelia with native embedded immune cells (T, B, NK, macrophages). Robust droplet-based, single-cell simultaneous determination of gene expression and immune repertoire indicated that PDO TILs accurately preserved the original tumor T cell receptor (TCR) spectrum. Crucially, human and murine PDOs successfully modeled immune checkpoint blockade (ICB) with anti-PD-1- and/or anti-PD-L1 expanding and activating tumor antigen-specific TILs and eliciting tumor cytotoxicity. Organoid-based propagation of primary tumor epithelium en bloc with endogenous immune stroma should enable immuno-oncology investigations within the TME and facilitate personalized immunotherapy testing.

Defining a Cancer Dependency Map.

Most human epithelial tumors harbor numerous alterations, making it difficult to predict which genes are required for tumor survival. To systematically identify cancer dependencies, we analyzed 501 genome-scale loss-of-function screens performed in diverse human cancer cell lines. We developed DEMETER, an analytical framework that segregates on- from off-target effects of RNAi. 769 genes were differentially required in subsets of these cell lines at a threshold of six SDs from the mean. We found predictive models for 426 dependencies (55%) by nonlinear regression modeling considering 66,646 molecular features. Many dependencies fall into a limited number of classes, and unexpectedly, in 82% of models, the top biomarkers were expression based. We demonstrated the basis behind one such predictive model linking hypermethylation of the UBB ubiquitin gene to a dependency on UBC. Together, these observations provide a foundation for a cancer dependency map that facilitates the prioritization of therapeutic targets.

KRAS and YAP1 converge to regulate EMT and tumor survival.

Cancer cells that express oncogenic alleles of RAS typically require sustained expression of the mutant allele for survival, but the molecular basis of this oncogene dependency remains incompletely understood. To identify genes that can functionally substitute for oncogenic RAS, we systematically expressed 15,294 open reading frames in a human KRAS-dependent colon cancer cell line engineered to express an inducible KRAS-specific shRNA. We found 147 genes that promoted survival upon KRAS suppression. In particular, the transcriptional coactivator YAP1 rescued cell viability in KRAS-dependent cells upon suppression of KRAS and was required for KRAS-induced cell transformation. Acquired resistance to Kras suppression in a Kras-driven murine lung cancer model also involved increased YAP1 signaling. KRAS and YAP1 converge on the transcription factor FOS and activate a transcriptional program involved in regulating the epithelial-mesenchymal transition (EMT). Together, these findings implicate transcriptional regulation of EMT by YAP1 as a significant component of oncogenic RAS signaling.

HIF1A employs CDK8-mediator to stimulate RNAPII elongation in response to hypoxia.

The transcription factor HIF1A is a key mediator of the cellular response to hypoxia. Despite the importance of HIF1A in homeostasis and various pathologies, little is known about how it regulates RNA polymerase II (RNAPII). We report here that HIF1A employs a specific variant of the Mediator complex to stimulate RNAPII elongation. The Mediator-associated kinase CDK8, but not the paralog CDK19, is required for induction of many HIF1A target genes. HIF1A induces binding of CDK8-Mediator and the super elongation complex (SEC), containing AFF4 and CDK9, to alleviate RNAPII pausing. CDK8 is dispensable for HIF1A chromatin binding and histone acetylation, but it is essential for binding of SEC and RNAPII elongation. Global analysis of active RNAPII reveals that hypoxia-inducible genes are paused and active prior to their induction. Our results provide a mechanistic link between HIF1A and CDK8, two potent oncogenes, in the cellular response to hypoxia.

Structure-function analysis of the SHOC2-MRAS-PP1C holophosphatase complex.

Receptor tyrosine kinase (RTK)-RAS signalling through the downstream mitogen-activated protein kinase (MAPK) cascade regulates cell proliferation and survival. The SHOC2-MRAS-PP1C holophosphatase complex functions as a key regulator of RTK-RAS signalling by removing an inhibitory phosphorylation event on the RAF family of proteins to potentiate MAPK signalling1. SHOC2 forms a ternary complex with MRAS and PP1C, and human germline gain-of-function mutations in this complex result in congenital RASopathy syndromes2-5. However, the structure and assembly of this complex are poorly understood. Here we use cryo-electron microscopy to resolve the structure of the SHOC2-MRAS-PP1C complex. We define the biophysical principles of holoenzyme interactions, elucidate the assembly order of the complex, and systematically interrogate the functional consequence of nearly all of the possible missense variants of SHOC2 through deep mutational scanning. We show that SHOC2 binds PP1C and MRAS through the concave surface of the leucine-rich repeat region and further engages PP1C through the N-terminal disordered region that contains a cryptic RVXF motif. Complex formation is initially mediated by interactions between SHOC2 and PP1C and is stabilized by the binding of GTP-loaded MRAS. These observations explain how mutant versions of SHOC2 in RASopathies and cancer stabilize the interactions of complex members to enhance holophosphatase activity. Together, this integrative structure-function model comprehensively defines key binding interactions within the SHOC2-MRAS-PP1C holophosphatase complex and will inform therapeutic development .

EGFR in limbo.

The epidermal growth factor receptor (EGFR) is normally activated by ligand-induced dimerization. Oncogenic mutations in EGFR promote activation in a largely ligand-independent manner. Shan et al. uncover a partially disordered state of EGFR kinase, providing evidence that oncogenic mutations counteract this intrinsic structural instability to promote dimerization and aberrant activation.

ß-Catenin-driven cancers require a YAP1 transcriptional complex for survival and tumorigenesis.

Wnt/ß-catenin signaling plays a key role in the pathogenesis of colon and other cancers; emerging evidence indicates that oncogenic ß-catenin regulates several biological processes essential for cancer initiation and progression. To decipher the role of ß-catenin in transformation, we classified ß-catenin activity in 85 cancer cell lines in which we performed genome-scale loss-of-function screens and found that ß-catenin active cancers are dependent on a signaling pathway involving the transcriptional regulator YAP1. Specifically, we found that YAP1 and the transcription factor TBX5 form a complex with ß-catenin. Phosphorylation of YAP1 by the tyrosine kinase YES1 leads to localization of this complex to the promoters of antiapoptotic genes, including BCL2L1 and BIRC5. A small-molecule inhibitor of YES1 impeded the proliferation of ß-catenin-dependent cancers in both cell lines and animal models. These observations define a ß-catenin-YAP1-TBX5 complex essential to the transformation and survival of ß-catenin-driven cancers.

Cancer vulnerabilities unveiled by genomic loss.

Due to genome instability, most cancers exhibit loss of regions containing tumor suppressor genes and collateral loss of other genes. To identify cancer-specific vulnerabilities that are the result of copy number losses, we performed integrated analyses of genome-wide copy number and RNAi profiles and identified 56 genes for which gene suppression specifically inhibited the proliferation of cells harboring partial copy number loss of that gene. These CYCLOPS (copy number alterations yielding cancer liabilities owing to partial loss) genes are enriched for spliceosome, proteasome, and ribosome components. One CYCLOPS gene, PSMC2, encodes an essential member of the 19S proteasome. Normal cells express excess PSMC2, which resides in a complex with PSMC1, PSMD2, and PSMD5 and acts as a reservoir protecting cells from PSMC2 suppression. Cells harboring partial PSMC2 copy number loss lack this complex and die after PSMC2 suppression. These observations define a distinct class of cancer-specific liabilities resulting from genome instability.

Biologically informed deep neural network for prostate cancer discovery.

The determination of molecular features that mediate clinically aggressive phenotypes in prostate cancer remains a major biological and clinical challenge1,2. Recent advances in interpretability of machine learning models as applied to biomedical problems may enable discovery and prediction in clinical cancer genomics3-5. Here we developed P-NET-a biologically informed deep learning model-to stratify patients with prostate cancer by treatment-resistance state and evaluate molecular drivers of treatment resistance for therapeutic targeting through complete model interpretability. We demonstrate that P-NET can predict cancer state using molecular data with a performance that is superior to other modelling approaches. Moreover, the biological interpretability within P-NET revealed established and novel molecularly altered candidates, such as MDM4 and FGFR1, which were implicated in predicting advanced disease and validated in vitro. Broadly, biologically informed fully interpretable neural networks enable preclinical discovery and clinical prediction in prostate cancer and may have general applicability across cancer types.

BET bromodomain inhibition as a therapeutic strategy to target c-Myc.

MYC contributes to the pathogenesis of a majority of human cancers, yet strategies to modulate the function of the c-Myc oncoprotein do not exist. Toward this objective, we have targeted MYC transcription by interfering with chromatin-dependent signal transduction to RNA polymerase, specifically by inhibiting the acetyl-lysine recognition domains (bromodomains) of putative coactivator proteins implicated in transcriptional initiation and elongation. Using a selective small-molecule bromodomain inhibitor, JQ1, we identify BET bromodomain proteins as regulatory factors for c-Myc. BET inhibition by JQ1 downregulates MYC transcription, followed by genome-wide downregulation of Myc-dependent target genes. In experimental models of multiple myeloma, a Myc-dependent hematologic malignancy, JQ1 produces a potent antiproliferative effect associated with cell-cycle arrest and cellular senescence. Efficacy of JQ1 in three murine models of multiple myeloma establishes the therapeutic rationale for BET bromodomain inhibition in this disease and other malignancies characterized by pathologic activation of c-Myc.

Binding of TMPRSS2-ERG to BAF Chromatin Remodeling Complexes Mediates Prostate Oncogenesis.

Chromosomal rearrangements resulting in the fusion of TMPRSS2, an androgen-regulated gene, and the ETS family transcription factor ERG occur in over half of prostate cancers. However, the mechanism by which ERG promotes oncogenic gene expression and proliferation remains incompletely understood. Here, we identify a binding interaction between ERG and the mammalian SWI/SNF (BAF) ATP-dependent chromatin remodeling complex, which is conserved among other oncogenic ETS factors, including ETV1, ETV4, and ETV5. We find that ERG drives genome-wide retargeting of BAF complexes in a manner dependent on binding of ERG to the ETS DNA motif. Moreover, ERG requires intact BAF complexes for chromatin occupancy and BAF complex ATPase activity for target gene regulation. In a prostate organoid model, BAF complexes are required for ERG-mediated basal-to-luminal transition, a hallmark of ERG activity in prostate cancer. These observations suggest a fundamental interdependence between ETS transcription factors and BAF chromatin remodeling complexes in cancer.

Discovery of antibodies and cognate surface targets for ovarian cancer by surface profiling.

Although antibodies targeting specific tumor-expressed antigens are the standard of care for some cancers, the identification of cancer-specific targets amenable to antibody binding has remained a bottleneck in development of new therapeutics. To overcome this challenge, we developed a high-throughput platform that allows for the unbiased, simultaneous discovery of antibodies and targets based on phenotypic binding profiles. Applying this platform to ovarian cancer, we identified a wide diversity of cancer targets including receptor tyrosine kinases, adhesion and migration proteins, proteases and proteins regulating angiogenesis in a single round of screening using genomics, flow cytometry, and mass spectrometry. In particular, we identified BCAM as a promising candidate for targeted therapy in high-grade serous ovarian cancers. More generally, this approach provides a rapid and flexible framework to identify cancer targets and antibodies.

Synthetic lethal interaction between oncogenic KRAS dependency and STK33 suppression in human cancer cells.

An alternative to therapeutic targeting of oncogenes is to perform "synthetic lethality" screens for genes that are essential only in the context of specific cancer-causing mutations. We used high-throughput RNA interference (RNAi) to identify synthetic lethal interactions in cancer cells harboring mutant KRAS, the most commonly mutated human oncogene. We find that cells that are dependent on mutant KRAS exhibit sensitivity to suppression of the serine/threonine kinase STK33 irrespective of tissue origin, whereas STK33 is not required by KRAS-independent cells. STK33 promotes cancer cell viability in a kinase activity-dependent manner by regulating the suppression of mitochondrial apoptosis mediated through S6K1-induced inactivation of the death agonist BAD selectively in mutant KRAS-dependent cells. These observations identify STK33 as a target for treatment of mutant KRAS-driven cancers and demonstrate the potential of RNAi screens for discovering functional dependencies created by oncogenic mutations that may enable therapeutic intervention for cancers with "undruggable" genetic alterations.

Unexpected pieces to the senescence puzzle.

Although initially described as the end state of cells after extended rounds of division in culture, it is now clear that cellular senescence induced by different stimuli plays an important role in tumor suppression in vivo. Three recent studies in Cell report that secreted proteins play an important role in enforcing the senescence response (Acosta et al., 2008; Kuilman et al., 2008; Wajapeyee et al., 2008). These new studies identify unanticipated contributors to this tumor-suppressing cell state.

Dependency of a therapy-resistant state of cancer cells on a lipid peroxidase pathway.

Plasticity of the cell state has been proposed to drive resistance to multiple classes of cancer therapies, thereby limiting their effectiveness. A high-mesenchymal cell state observed in human tumours and cancer cell lines has been associated with resistance to multiple treatment modalities across diverse cancer lineages, but the mechanistic underpinning for this state has remained incompletely understood. Here we molecularly characterize this therapy-resistant high-mesenchymal cell state in human cancer cell lines and organoids and show that it depends on a druggable lipid-peroxidase pathway that protects against ferroptosis, a non-apoptotic form of cell death induced by the build-up of toxic lipid peroxides. We show that this cell state is characterized by activity of enzymes that promote the synthesis of polyunsaturated lipids. These lipids are the substrates for lipid peroxidation by lipoxygenase enzymes. This lipid metabolism creates a dependency on pathways converging on the phospholipid glutathione peroxidase (GPX4), a selenocysteine-containing enzyme that dissipates lipid peroxides and thereby prevents the iron-mediated reactions of peroxides that induce ferroptotic cell death. Dependency on GPX4 was found to exist across diverse therapy-resistant states characterized by high expression of ZEB1, including epithelial-mesenchymal transition in epithelial-derived carcinomas, TGFß-mediated therapy-resistance in melanoma, treatment-induced neuroendocrine transdifferentiation in prostate cancer, and sarcomas, which are fixed in a mesenchymal state owing to their cells of origin. We identify vulnerability to ferroptic cell death induced by inhibition of a lipid peroxidase pathway as a feature of therapy-resistant cancer cells across diverse mesenchymal cell-state contexts.

Discovery of a selective inhibitor of doublecortin like kinase 1.

Doublecortin like kinase 1 (DCLK1) is an understudied kinase that is upregulated in a wide range of cancers, including pancreatic ductal adenocarcinoma (PDAC). However, little is known about its potential as a therapeutic target. We used chemoproteomic profiling and structure-based design to develop a selective, in vivo-compatible chemical probe of the DCLK1 kinase domain, DCLK1-IN-1. We demonstrate activity of DCLK1-IN-1 against clinically relevant patient-derived PDAC organoid models and use a combination of RNA-sequencing, proteomics and phosphoproteomics analysis to reveal that DCLK1 inhibition modulates proteins and pathways associated with cell motility in this context. DCLK1-IN-1 will serve as a versatile tool to investigate DCLK1 biology and establish its role in cancer.