Targeting TRIP13 in favorable histology Wilms tumor with nuclear export inhibitors synergizes with doxorubicin.

Wilms tumor (WT) is the most common renal malignancy of childhood. Despite improvements in the overall survival, relapse occurs in ~15% of patients with favorable histology WT (FHWT). Half of these patients will succumb to their disease. Identifying novel targeted therapies remains challenging in part due to the lack of faithful preclinical in vitro models. Here we establish twelve patient-derived WT cell lines and demonstrate that these models faithfully recapitulate WT biology using genomic and transcriptomic techniques. We then perform loss-of-function screens to identify the nuclear export gene, XPO1, as a vulnerability. We find that the FDA approved XPO1 inhibitor, KPT-330, suppresses TRIP13 expression, which is required for survival. We further identify synergy between KPT-330 and doxorubicin, a chemotherapy used in high-risk FHWT. Taken together, we identify XPO1 inhibition with KPT-330 as a potential therapeutic option to treat FHWTs and in combination with doxorubicin, leads to durable remissions in vivo.

RHOAL57V drives the development of diffuse gastric cancer through IGF1R-PAK1-YAP1 signaling.

Cancer-associated mutations in the guanosine triphosphatase (GTPase) RHOA are found at different locations from the mutational hotspots in the structurally and biochemically related RAS. Tyr42-to-Cys (Y42C) and Leu57-to-Val (L57V) substitutions are the two most prevalent RHOA mutations in diffuse gastric cancer (DGC). RHOAY42C exhibits a gain-of-function phenotype and is an oncogenic driver in DGC. Here, we determined how RHOAL57V promotes DGC growth. In mouse gastric organoids with deletion of Cdh1, which encodes the cell adhesion protein E-cadherin, the expression of RHOAL57V, but not of wild-type RHOA, induced an abnormal morphology similar to that of patient-derived DGC organoids. RHOAL57V also exhibited a gain-of-function phenotype and promoted F-actin stress fiber formation and cell migration. RHOAL57V retained interaction with effectors but exhibited impaired RHOA-intrinsic and GAP-catalyzed GTP hydrolysis, which favored formation of the active GTP-bound state. Introduction of missense mutations at KRAS residues analogous to Tyr42 and Leu57 in RHOA did not activate KRAS oncogenic potential, indicating distinct functional effects in otherwise highly related GTPases. Both RHOA mutants stimulated the transcriptional co-activator YAP1 through actin dynamics to promote DGC progression; however, RHOAL57V additionally did so by activating the kinases IGF1R and PAK1, distinct from the FAK-mediated mechanism induced by RHOAY42C. Our results reveal that RHOAL57V and RHOAY42C drive the development of DGC through distinct biochemical and signaling mechanisms.

Functional screening of amplification outlier oncogenes in organoid models of early tumorigenesis.

Somatic copy number gains are pervasive across cancer types, yet their roles in oncogenesis are insufficiently evaluated. This inadequacy is partly due to copy gains spanning large chromosomal regions, obscuring causal loci. Here, we employed organoid modeling to evaluate candidate oncogenic loci identified via integrative computational analysis of extreme copy gains overlapping with extreme expression dysregulation in The Cancer Genome Atlas. Subsets of "outlier" candidates were contextually screened as tissue-specific cDNA lentiviral libraries within cognate esophagus, oral cavity, colon, stomach, pancreas, and lung organoids bearing initial oncogenic mutations. Iterative analysis nominated the kinase DYRK2 at 12q15 as an amplified head and neck squamous carcinoma oncogene in p53-/- oral mucosal organoids. Similarly, FGF3, amplified at 11q13 in 41% of esophageal squamous carcinomas, promoted p53-/- esophageal organoid growth reversible by small molecule and soluble receptor antagonism of FGFRs. Our studies establish organoid-based contextual screening of candidate genomic drivers, enabling functional evaluation during early tumorigenesis.

Partial gene suppression improves identification of cancer vulnerabilities when CRISPR-Cas9 knockout is pan-lethal.

BACKGROUND: Hundreds of functional genomic screens have been performed across a diverse set of cancer contexts, as part of efforts such as the Cancer Dependency Map, to identify gene dependencies-genes whose loss of function reduces cell viability or fitness. Recently, large-scale screening efforts have shifted from RNAi to CRISPR-Cas9, due to superior efficacy and specificity. However, many effective oncology drugs only partially inhibit their protein targets, leading us to question whether partial suppression of genes using RNAi could reveal cancer vulnerabilities that are missed by complete knockout using CRISPR-Cas9. Here, we compare CRISPR-Cas9 and RNAi dependency profiles of genes across approximately 400 matched cancer cell lines. RESULTS: We find that CRISPR screens accurately identify more gene dependencies per cell line, but the majority of each cell line's dependencies are part of a set of 1867 genes that are shared dependencies across the entire collection (pan-lethals). While RNAi knockdown of about 30% of these genes is also pan-lethal, approximately 50% have selective dependency patterns across cell lines, suggesting they could still be cancer vulnerabilities. The accuracy of the unique RNAi selectivity is supported by associations to multi-omics profiles, drug sensitivity, and other expected co-dependencies. CONCLUSIONS: Incorporating RNAi data for genes that are pan-lethal knockouts facilitates the discovery of a wider range of gene targets than could be detected using the CRISPR dataset alone. This can aid in the interpretation of contrasting results obtained from CRISPR and RNAi screens and reinforce the importance of partial gene suppression methods in building a cancer dependency map.

An Atlas of Variant Effects to understand the genome at nucleotide resolution.

Sequencing has revealed hundreds of millions of human genetic variants, and continued efforts will only add to this variant avalanche. Insufficient information exists to interpret the effects of most variants, limiting opportunities for precision medicine and comprehension of genome function. A solution lies in experimental assessment of the functional effect of variants, which can reveal their biological and clinical impact. However, variant effect assays have generally been undertaken reactively for individual variants only after and, in most cases long after, their first observation. Now, multiplexed assays of variant effect can characterise massive numbers of variants simultaneously, yielding variant effect maps that reveal the function of every possible single nucleotide change in a gene or regulatory element. Generating maps for every protein encoding gene and regulatory element in the human genome would create an 'Atlas' of variant effect maps and transform our understanding of genetics and usher in a new era of nucleotide-resolution functional knowledge of the genome. An Atlas would reveal the fundamental biology of the human genome, inform human evolution, empower the development and use of therapeutics and maximize the utility of genomics for diagnosing and treating disease. The Atlas of Variant Effects Alliance is an international collaborative group comprising hundreds of researchers, technologists and clinicians dedicated to realising an Atlas of Variant Effects to help deliver on the promise of genomics.

Genome-scale functional genomics identify genes preferentially essential for multiple myeloma cells compared to other neoplasias.

Clinical progress in multiple myeloma (MM), an incurable plasma cell (PC) neoplasia, has been driven by therapies that have limited applications beyond MM/PC neoplasias and do not target specific oncogenic mutations in MM. Instead, these agents target pathways critical for PC biology yet largely dispensable for malignant or normal cells of most other lineages. Here we systematically characterized the lineage-preferential molecular dependencies of MM through genome-scale clustered regularly interspaced short palindromic repeats (CRISPR) studies in 19 MM versus hundreds of non-MM lines and identified 116 genes whose disruption more significantly affects MM cell fitness compared with other malignancies. These genes, some known, others not previously linked to MM, encode transcription factors, chromatin modifiers, endoplasmic reticulum components, metabolic regulators or signaling molecules. Most of these genes are not among the top amplified, overexpressed or mutated in MM. Functional genomics approaches thus define new therapeutic targets in MM not readily identifiable by standard genomic, transcriptional or epigenetic profiling analyses.

Mapping the landscape of genetic dependencies in chordoma.

Identifying the spectrum of genes required for cancer cell survival can reveal essential cancer circuitry and therapeutic targets, but such a map remains incomplete for many cancer types. We apply genome-scale CRISPR-Cas9 loss-of-function screens to map the landscape of selectively essential genes in chordoma, a bone cancer with few validated targets. This approach confirms a known chordoma dependency, TBXT (T; brachyury), and identifies a range of additional dependencies, including PTPN11, ADAR, PRKRA, LUC7L2, SRRM2, SLC2A1, SLC7A5, FANCM, and THAP1. CDK6, SOX9, and EGFR, genes previously implicated in chordoma biology, are also recovered. We find genomic and transcriptomic features that predict specific dependencies, including interferon-stimulated gene expression, which correlates with ADAR dependence and is elevated in chordoma. Validating the therapeutic relevance of dependencies, small-molecule inhibitors of SHP2, encoded by PTPN11, have potent preclinical efficacy against chordoma. Our results generate an emerging map of chordoma dependencies to enable biological and therapeutic hypotheses.

USP9X mediates an acute adaptive response to MAPK suppression in pancreatic cancer but creates multiple actionable therapeutic vulnerabilities.

Pancreatic ductal adenocarcinomas (PDACs) frequently harbor KRAS mutations. Although MEK inhibitors represent a plausible therapeutic option, most PDACs are innately resistant to these agents. Here, we identify a critical adaptive response that mediates resistance. Specifically, we show that MEK inhibitors upregulate the anti-apoptotic protein Mcl-1 by triggering an association with its deubiquitinase, USP9X, resulting in acute Mcl-1 stabilization and protection from apoptosis. Notably, these findings contrast the canonical positive regulation of Mcl-1 by RAS/ERK. We further show that Mcl-1 inhibitors and cyclin-dependent kinase (CDK) inhibitors, which suppress Mcl-1 transcription, prevent this protective response and induce tumor regression when combined with MEK inhibitors. Finally, we identify USP9X as an additional potential therapeutic target. Together, these studies (1) demonstrate that USP9X regulates a critical mechanism of resistance in PDAC, (2) reveal an unexpected mechanism of Mcl-1 regulation in response to RAS pathway suppression, and (3) provide multiple distinct promising therapeutic strategies for this deadly malignancy.

Aberrant cell state plasticity mediated by developmental reprogramming precedes colorectal cancer initiation.

Cell state plasticity is carefully regulated in adult epithelia to prevent cancer. The aberrant expansion of the normally restricted capability for cell state plasticity in neoplasia is poorly defined. Using genetically engineered and carcinogen-induced mouse models of intestinal neoplasia, we observed that impaired differentiation is a conserved event preceding cancer development. Single-cell RNA sequencing (scRNA-seq) of premalignant lesions from mouse models and a patient with hereditary polyposis revealed that cancer initiates by adopting an aberrant transcriptional state characterized by regenerative activity, marked by Ly6a (Sca-1), and reactivation of fetal intestinal genes, including Tacstd2 (Trop2). Genetic inactivation of Sox9 prevented adenoma formation, obstructed the emergence of regenerative and fetal programs, and restored multilineage differentiation by scRNA-seq. Expanded chromatin accessibility at regeneration and fetal genes upon Apc inactivation was reduced by concomitant Sox9 suppression. These studies indicate that aberrant cell state plasticity mediated by unabated regenerative activity and developmental reprogramming precedes cancer development.

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.

Systematic Interrogation of Tumor Cell Resistance to Chimeric Antigen Receptor T-cell Therapy in Pancreatic Cancer.

Chimeric antigen receptor (CAR) T-cell therapy can lead to dramatic clinical responses in B-cell malignancies. However, early clinical trials with CAR T-cell therapy in non-B-cell malignancies have been disappointing to date, suggesting that tumor-intrinsic features contribute to resistance. To investigate tumor-intrinsic modes of resistance, we performed genome scale CRISPR-Cas9 screens in mesothelin (MSLN)-expressing pancreatic cancer cells. Co-culture with MSLN-targeting CAR T cells identified both antigen-dependent and antigen-independent modes of resistance. In particular, loss of the majority of the genes involved in the pathway responsible for GPI-anchor biosynthesis and attachment abrogated the ability of CAR T cells to target pancreatic cancer cells, suggesting that disruption of this pathway may permit MSLN CAR T-cell evasion in the clinic. Antigen-independent mediators of CAR T-cell response included members of the death receptor pathway as well as genes that regulate tumor transcriptional responses, including TFAP4 and INTS12. TFAP4-mediated CAR T resistance depended on the NFκB transcription factor p65, indicating that tumor resistance to CAR T-cell therapy likely involves alterations in tumor-intrinsic states. Overall, this study uncovers multiple antigen-dependent and -independent mechanisms of CAR T-cell evasion by pancreatic cancer, paving the way for overcoming resistance in this disease that is notoriously refractory to immunotherapy. SIGNIFICANCE: The identification and validation of key determinants of CAR T-cell response in pancreatic cancer provide insights into the landscape of tumor cell intrinsic resistance mechanisms and into approaches to improve therapeutic efficacy.

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.

A Ubiquitination Cascade Regulating the Integrated Stress Response and Survival in Carcinomas.

Systematic identification of signaling pathways required for the fitness of cancer cells will facilitate the development of new cancer therapies. We used gene essentiality measurements in 1,086 cancer cell lines to identify selective coessentiality modules and found that a ubiquitin ligase complex composed of UBA6, BIRC6, KCMF1, and UBR4 is required for the survival of a subset of epithelial tumors that exhibit a high degree of aneuploidy. Suppressing BIRC6 in cell lines that are dependent on this complex led to a substantial reduction in cell fitness in vitro and potent tumor regression in vivo. Mechanistically, BIRC6 suppression resulted in selective activation of the integrated stress response (ISR) by stabilization of the heme-regulated inhibitor, a direct ubiquitination target of the UBA6/BIRC6/KCMF1/UBR4 complex. These observations uncover a novel ubiquitination cascade that regulates ISR and highlight the potential of ISR activation as a new therapeutic strategy. SIGNIFICANCE: We describe the identification of a heretofore unrecognized ubiquitin ligase complex that prevents the aberrant activation of the ISR in a subset of cancer cells. This provides a novel insight on the regulation of ISR and exposes a therapeutic opportunity to selectively eliminate these cancer cells. See related commentary Leli and Koumenis, p. 535. This article is highlighted in the In This Issue feature, p. 517.

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.

VRK1 as a synthetic lethal target in VRK2 promoter-methylated cancers of the nervous system.

Collateral lethality occurs when loss of a gene/protein renders cancer cells dependent on its remaining paralog. Combining genome-scale CRISPR/Cas9 loss-of-function screens with RNA sequencing in over 900 cancer cell lines, we found that cancers of nervous system lineage, including adult and pediatric gliomas and neuroblastomas, required the nuclear kinase vaccinia-related kinase 1 (VRK1) for their survival in vivo. VRK1 dependency was inversely correlated with expression of its paralog VRK2. VRK2 knockout sensitized cells to VRK1 loss, and conversely, VRK2 overexpression increased cell fitness in the setting of VRK1 loss. DNA methylation of the VRK2 promoter was associated with low VRK2 expression in human neuroblastomas and adult and pediatric gliomas. Mechanistically, depletion of VRK1 reduced barrier-to-autointegration factor phosphorylation during mitosis, resulting in DNA damage and apoptosis. Together, these studies identify VRK1 as a synthetic lethal target in VRK2 promoter-methylated adult and pediatric gliomas and neuroblastomas.

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 .

CREB5 reprograms FOXA1 nuclear interactions to promote resistance to androgen receptor-targeting therapies.

Metastatic castration-resistant prostate cancers (mCRPCs) are treated with therapies that antagonize the androgen receptor (AR). Nearly all patients develop resistance to AR-targeted therapies (ARTs). Our previous work identified CREB5 as an upregulated target gene in human mCRPC that promoted resistance to all clinically approved ART. The mechanisms by which CREB5 promotes progression of mCRPC or other cancers remains elusive. Integrating ChIP-seq and rapid immunoprecipitation and mass spectroscopy of endogenous proteins, we report that cells overexpressing CREB5 demonstrate extensive reprogramming of nuclear protein-protein interactions in response to the ART agent enzalutamide. Specifically, CREB5 physically interacts with AR, the pioneering actor FOXA1, and other known co-factors of AR and FOXA1 at transcription regulatory elements recently found to be active in mCRPC patients. We identified a subset of CREB5/FOXA1 co-interacting nuclear factors that have critical functions for AR transcription (GRHL2, HOXB13) while others (TBX3, NFIC) regulated cell viability and ART resistance and were amplified or overexpressed in mCRPC. Upon examining the nuclear protein interactions and the impact of CREB5 expression on the mCRPC patient transcriptome, we found that CREB5 was associated with Wnt signaling and epithelial to mesenchymal transitions, implicating these pathways in CREB5/FOXA1-mediated ART resistance. Overall, these observations define the molecular interactions among CREB5, FOXA1, and pathways that promote ART resistance.