Macroautophagy (hereafter referred to as autophagy) is an evolutionarily conserved catabolic process necessary for normal recycling of cellular constituents and for appropriate response to cellular stress. Although several genes belonging to the core molecular machinery involved in autophagosome formation have been discovered, relatively little is known about the nature of signaling networks controlling autophagy upon intracellular or extracellular stimuli. We discovered ATG8-like proteins (MAP1LC3B, GABARAP and GABARAPL1) as novel interactors of MAPK15/ERK8, a MAP kinase involved in cell proliferation and transformation. Based on the role of these proteins in the autophagic process, we demonstrated that MAPK15 is indeed localized to autophagic compartments and increased, in a kinase-dependent fashion, ATG8-like proteins lipidation, autophagosome formation and SQSTM1 degradation, while decreasing LC3B inhibitory phosphorylation. Interestingly, we also identified a conserved LC3-interacting region (LIR) in MAPK15 responsible for its interaction with ATG8-like proteins, for its localization to autophagic structures and, consequently, for stimulation of the formation of these compartments. Furthermore, we reveal that MAPK15 activity was induced in response to serum and amino-acid starvation and that this stimulus, in turn, required endogenous MAPK15 expression to induce the autophagic process. Altogether, these results suggested a new function for MAPK15 as a regulator of autophagy, acting through interaction with ATG8 family proteins. Also, based on the key role of this process in several human diseases, these results supported the use of this MAP kinase as a potential novel therapeutic target.

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.

Whole-exome sequencing of cell-free DNA (cfDNA) could enable comprehensive profiling of tumors from blood but the genome-wide concordance between cfDNA and tumor biopsies is uncertain. Here we report ichorCNA, software that quantifies tumor content in cfDNA from 0.1× coverage whole-genome sequencing data without prior knowledge of tumor mutations. We apply ichorCNA to 1439 blood samples from 520 patients with metastatic prostate or breast cancers. In the earliest tested sample for each patient, 34% of patients have ≥10% tumor-derived cfDNA, sufficient for standard coverage whole-exome sequencing. Using whole-exome sequencing, we validate the concordance of clonal somatic mutations (88%), copy number alterations (80%), mutational signatures, and neoantigens between cfDNA and matched tumor biopsies from 41 patients with ≥10% cfDNA tumor content. In summary, we provide methods to identify patients eligible for comprehensive cfDNA profiling, revealing its applicability to many patients, and demonstrate high concordance of cfDNA and metastatic tumor whole-exome sequencing.

We previously piloted the concept of a Connectivity Map (CMap), whereby genes, drugs, and disease states are connected by virtue of common gene-expression signatures. Here, we report more than a 1,000-fold scale-up of the CMap as part of the NIH LINCS Consortium, made possible by a new, low-cost, high-throughput reduced representation expression profiling method that we term L1000. We show that L1000 is highly reproducible, comparable to RNA sequencing, and suitable for computational inference of the expression levels of 81% of non-measured transcripts. We further show that the expanded CMap can be used to discover mechanism of action of small molecules, functionally annotate genetic variants of disease genes, and inform clinical trials. The 1.3 million L1000 profiles described here, as well as tools for their analysis, are available at https://clue.io.

Genetic and chemical perturbations impact diverse cellular phenotypes, including multiple indicators of cell health. These readouts reveal toxicity and antitumorigenic effects relevant to drug discovery and personalized medicine. We developed two customized microscopy assays, one using four targeted reagents and the other three targeted reagents, to collectively measure 70 specific cell health phenotypes including proliferation, apoptosis, reactive oxygen species, DNA damage, and cell cycle stage. We then tested an approach to predict multiple cell health phenotypes using Cell Painting, an inexpensive and scalable image-based morphology assay. In matched CRISPR perturbations of three cancer cell lines, we collected both Cell Painting and cell health data. We found that simple machine learning algorithms can predict many cell health readouts directly from Cell Painting images, at less than half the cost. We hypothesized that these models can be applied to accurately predict cell health assay outcomes for any future or existing Cell Painting dataset. For Cell Painting images from a set of 1500+ compound perturbations across multiple doses, we validated predictions by orthogonal assay readouts. We provide a web app to browse predictions: http://broad.io/cell-health-app. Our approach can be used to add cell health annotations to Cell Painting datasets.

Cancer evolution poses a central obstacle to cure, as resistant clones expand under therapeutic selection pressures. Genome sequencing of relapsed disease can nominate genomic alterations conferring resistance but sample collection lags behind, limiting therapeutic innovation. Genome-wide screens offer a complementary approach to chart the compendium of escape genotypes, anticipating clinical resistance. We report genome-wide open reading frame (ORF) resistance screens for first- and third-generation epidermal growth factor receptor (EGFR) inhibitors and a MEK inhibitor. Using serial sampling, dose gradients, and mathematical modeling, we generate genotype-fitness maps across therapeutic contexts and identify alterations that escape therapy. Our data expose varying dose-fitness relationship across genotypes, ranging from complete dose invariance to paradoxical dose dependency where fitness increases in higher doses. We predict fitness with combination therapy and compare these estimates to genome-wide fitness maps of drug combinations, identifying genotypes where combination therapy results in unexpected inferior effectiveness. These data are applied to nominate combination optimization strategies to forestall resistant disease.

Despite advances in precision medicine, the clinical prospects for patients with ovarian and uterine cancers have not substantially improved. Here, we analyzed genome-scale CRISPR-Cas9 loss-of-function screens across 851 human cancer cell lines and found that frequent overexpression of SLC34A2-encoding a phosphate importer-is correlated with sensitivity to loss of the phosphate exporter XPR1, both in vitro and in vivo. In patient-derived tumor samples, we observed frequent PAX8-dependent overexpression of SLC34A2, XPR1 copy number amplifications and XPR1 messenger RNA overexpression. Mechanistically, in SLC34A2-high cancer cell lines, genetic or pharmacologic inhibition of XPR1-dependent phosphate efflux leads to the toxic accumulation of intracellular phosphate. Finally, we show that XPR1 requires the novel partner protein KIDINS220 for proper cellular localization and activity, and that disruption of this protein complex results in acidic "vacuolar" structures preceding cell death. These data point to the XPR1-KIDINS220 complex and phosphate dysregulation as a therapeutic vulnerability in ovarian cancer.

CRISPR-Cas9-based genetic screens have successfully identified cell type-dependent liabilities in cancer, including acute myeloid leukemia (AML), a devastating hematologic malignancy with poor overall survival. Because most of these screens have been performed in vitro using established cell lines, evaluating the physiologic relevance of these targets is critical. We have established a CRISPR screening approach using orthotopic xenograft models to validate and prioritize AML-enriched dependencies in vivo, including in CRISPR-competent AML patient-derived xenograft (PDX) models tractable for genome editing. Our integrated pipeline has revealed several targets with translational value, including SLC5A3 as a metabolic vulnerability for AML addicted to exogenous myo-inositol and MARCH5 as a critical guardian to prevent apoptosis in AML. MARCH5 repression enhanced the efficacy of BCL2 inhibitors such as venetoclax, further highlighting the clinical potential of targeting MARCH5 in AML. Our study provides a valuable strategy for discovery and prioritization of new candidate AML therapeutic targets. SIGNIFICANCE: There is an unmet need to improve the clinical outcome of AML. We developed an integrated in vivo screening approach to prioritize and validate AML dependencies with high translational potential. We identified SLC5A3 as a metabolic vulnerability and MARCH5 as a critical apoptosis regulator in AML, both of which represent novel therapeutic opportunities.This article is highlighted in the In This Issue feature, p. 275.

Human genome variation contributes to diversity in neurodevelopmental outcomes and vulnerabilities; recognizing the underlying molecular and cellular mechanisms will require scalable approaches. Here, we describe a "cell village" experimental platform we used to analyze genetic, molecular, and phenotypic heterogeneity across neural progenitor cells from 44 human donors cultured in a shared in vitro environment using algorithms (Dropulation and Census-seq) to assign cells and phenotypes to individual donors. Through rapid induction of human stem cell-derived neural progenitor cells, measurements of natural genetic variation, and CRISPR-Cas9 genetic perturbations, we identified a common variant that regulates antiviral IFITM3 expression and explains most inter-individual variation in susceptibility to the Zika virus. We also detected expression QTLs corresponding to GWAS loci for brain traits and discovered novel disease-relevant regulators of progenitor proliferation and differentiation such as CACHD1. This approach provides scalable ways to elucidate the effects of genes and genetic variation on cellular phenotypes.

Unlike most tumor suppressor genes, the most common genetic alterations in tumor protein p53 (TP53) are missense mutations1,2. Mutant p53 protein is often abundantly expressed in cancers and specific allelic variants exhibit dominant-negative or gain-of-function activities in experimental models3-8. To gain a systematic view of p53 function, we interrogated loss-of-function screens conducted in hundreds of human cancer cell lines and performed TP53 saturation mutagenesis screens in an isogenic pair of TP53 wild-type and null cell lines. We found that loss or dominant-negative inhibition of wild-type p53 function reliably enhanced cellular fitness. By integrating these data with the Catalog of Somatic Mutations in Cancer (COSMIC) mutational signatures database9,10, we developed a statistical model that describes the TP53 mutational spectrum as a function of the baseline probability of acquiring each mutation and the fitness advantage conferred by attenuation of p53 activity. Collectively, these observations show that widely-acting and tissue-specific mutational processes combine with phenotypic selection to dictate the frequencies of recurrent TP53 mutations.

Systematic efforts to sequence the cancer genome have identified large numbers of mutations and copy number alterations in human cancers. However, elucidating the functional consequences of these variants, and their interactions to drive or maintain oncogenic states, remains a challenge in cancer research. We developed REVEALER, a computational method that identifies combinations of mutually exclusive genomic alterations correlated with functional phenotypes, such as the activation or gene dependency of oncogenic pathways or sensitivity to a drug treatment. We used REVEALER to uncover complementary genomic alterations associated with the transcriptional activation of β-catenin and NRF2, MEK-inhibitor sensitivity, and KRAS dependency. REVEALER successfully identified both known and new associations, demonstrating the power of combining functional profiles with extensive characterization of genomic alterations in cancer genomes.

Oncogene-induced senescence functions to limit tumor development. However, a complete understanding of the signals that trigger this type of senescence is currently lacking. We found that mutations affecting NF1, Raf, and Ras induce a global negative feedback response that potently suppresses Ras and/or its effectors. Moreover, these signals promote senescence by inhibiting the Ras/PI3K pathway, which can impact the senescence machinery through HDM2 and FOXO. This negative feedback program is regulated in part by RasGEFs, Sprouty proteins, RasGAPs, and MKPs. Moreover, these signals function in vivo in benign human tumors. Thus, the ultimate response to the aberrant activation of the Ras pathway is a multifaceted negative feedback signaling network that terminates the oncogenic signal and participates in the senescence response.

Chronic HBV infection is a major cause of liver disease and cancer worldwide. Approaches for cure are lacking, and the knowledge of virus-host interactions is still limited. Here, we perform a genome-wide gain-of-function screen using a poorly permissive hepatoma cell line to uncover host factors enhancing HBV infection. Validation studies in primary human hepatocytes identified CDKN2C as an important host factor for HBV replication. CDKN2C is overexpressed in highly permissive cells and HBV-infected patients. Mechanistic studies show a role for CDKN2C in inducing cell cycle G1 arrest through inhibition of CDK4/6 associated with the upregulation of HBV transcription enhancers. A correlation between CDKN2C expression and disease progression in HBV-infected patients suggests a role in HBV-induced liver disease. Taken together, we identify a previously undiscovered clinically relevant HBV host factor, allowing the development of improved infectious model systems for drug discovery and the study of the HBV life cycle.

Many biological processes, such as cancer metastasis, organismal development, and acquisition of resistance to cytotoxic therapy, rely on the emergence of rare sub-clones from a larger population. Understanding how the genetic and epigenetic features of diverse clones affect clonal fitness provides insight into molecular mechanisms underlying selective processes. While large-scale barcoding with NGS readout has facilitated cellular fitness assessment at the population level, this approach does not support characterization of clones prior to selection. Single-cell genomics methods provide high biological resolution, but are challenging to scale across large populations to probe rare clones and are destructive, limiting further functional analysis of important clones.

The mitogen-activated protein kinase (MAPK) pathway is a critical effector of oncogenic RAS signaling, and MAPK pathway inhibition may be an effective combination treatment strategy. We performed genome-scale loss-of-function CRISPR-Cas9 screens in the presence of a MEK1/2 inhibitor (MEKi) in KRAS-mutant pancreatic and lung cancer cell lines and identified genes that cooperate with MEK inhibition. While we observed heterogeneity in genetic modifiers of MEKi sensitivity across cell lines, several recurrent classes of synthetic lethal vulnerabilities emerged at the pathway level. Multiple members of receptor tyrosine kinase (RTK)-RAS-MAPK pathways scored as sensitizers to MEKi. In particular, we demonstrate that knockout, suppression, or degradation of SHOC2, a positive regulator of MAPK signaling, specifically cooperated with MEK inhibition to impair proliferation in RAS-driven cancer cells. The depletion of SHOC2 disrupted survival pathways triggered by feedback RTK signaling in response to MEK inhibition. Thus, these findings nominate SHOC2 as a potential target for combination therapy.

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.

Recent genome sequencing efforts have identified millions of somatic mutations in cancer. However, the functional impact of most variants is poorly understood. Here we characterize 194 somatic mutations identified in primary lung adenocarcinomas. We present an expression-based variant-impact phenotyping (eVIP) method that uses gene expression changes to distinguish impactful from neutral somatic mutations. eVIP identified 69% of mutations analyzed as impactful and 31% as functionally neutral. A subset of the impactful mutations induces xenograft tumor formation in mice and/or confers resistance to cellular EGFR inhibition. Among these impactful variants are rare somatic, clinically actionable variants including EGFR S645C, ARAF S214C and S214F, ERBB2 S418T, and multiple BRAF variants, demonstrating that rare mutations can be functionally important in cancer.

Oncogenic mutations in the serine/threonine kinase B-RAF (also known as BRAF) are found in 50-70% of malignant melanomas. Pre-clinical studies have demonstrated that the B-RAF(V600E) mutation predicts a dependency on the mitogen-activated protein kinase (MAPK) signalling cascade in melanoma-an observation that has been validated by the success of RAF and MEK inhibitors in clinical trials. However, clinical responses to targeted anticancer therapeutics are frequently confounded by de novo or acquired resistance. Identification of resistance mechanisms in a manner that elucidates alternative 'druggable' targets may inform effective long-term treatment strategies. Here we expressed ∼600 kinase and kinase-related open reading frames (ORFs) in parallel to interrogate resistance to a selective RAF kinase inhibitor. We identified MAP3K8 (the gene encoding COT/Tpl2) as a MAPK pathway agonist that drives resistance to RAF inhibition in B-RAF(V600E) cell lines. COT activates ERK primarily through MEK-dependent mechanisms that do not require RAF signalling. Moreover, COT expression is associated with de novo resistance in B-RAF(V600E) cultured cell lines and acquired resistance in melanoma cells and tissue obtained from relapsing patients following treatment with MEK or RAF inhibitors. We further identify combinatorial MAPK pathway inhibition or targeting of COT kinase activity as possible therapeutic strategies for reducing MAPK pathway activation in this setting. Together, these results provide new insights into resistance mechanisms involving the MAPK pathway and articulate an integrative approach through which high-throughput functional screens may inform the development of novel therapeutic strategies.

Malignant melanomas harbouring point mutations (Val600Glu) in the serine/threonine-protein kinase BRAF (BRAF(V600E)) depend on RAF-MEK-ERK signalling for tumour cell growth. RAF and MEK inhibitors show remarkable clinical efficacy in BRAF(V600E) melanoma; however, resistance to these agents remains a formidable challenge. Global characterization of resistance mechanisms may inform the development of more effective therapeutic combinations. Here we carried out systematic gain-of-function resistance studies by expressing more than 15,500 genes individually in a BRAF(V600E) melanoma cell line treated with RAF, MEK, ERK or combined RAF-MEK inhibitors. These studies revealed a cyclic-AMP-dependent melanocytic signalling network not previously associated with drug resistance, including G-protein-coupled receptors, adenyl cyclase, protein kinase A and cAMP response element binding protein (CREB). Preliminary analysis of biopsies from BRAF(V600E) melanoma patients revealed that phosphorylated (active) CREB was suppressed by RAF-MEK inhibition but restored in relapsing tumours. Expression of transcription factors activated downstream of MAP kinase and cAMP pathways also conferred resistance, including c-FOS, NR4A1, NR4A2 and MITF. Combined treatment with MAPK-pathway and histone-deacetylase inhibitors suppressed MITF expression and cAMP-mediated resistance. Collectively, these data suggest that oncogenic dysregulation of a melanocyte lineage dependency can cause resistance to RAF-MEK-ERK inhibition, which may be overcome by combining signalling- and chromatin-directed therapeutics.

Fusion-transcription factors (fusion-TFs) represent a class of driver oncoproteins that are difficult to therapeutically target. Recently, protein degradation has emerged as a strategy to target these challenging oncoproteins. The mechanisms that regulate fusion-TF stability, however, are generally unknown. Using CRISPR-Cas9 screening, we discovered tripartite motif-containing 8 (TRIM8) as an E3 ubiquitin ligase that ubiquitinates and degrades EWS/FLI, a driver fusion-TF in Ewing sarcoma. Moreover, we identified TRIM8 as a selective dependency in Ewing sarcoma compared with >700 other cancer cell lines. Mechanistically, TRIM8 knockout led to an increase in EWS/FLI protein levels that was not tolerated. EWS/FLI acts as a neomorphic substrate for TRIM8, defining the selective nature of the dependency. Our results demonstrate that fusion-TF protein stability is tightly regulated and highlight fusion oncoprotein-specific regulators as selective therapeutic targets. This study provides a tractable strategy to therapeutically exploit oncogene overdose in Ewing sarcoma and potentially other fusion-TF-driven cancers.