Leukemia stem cells (LSCs) represent a biologically distinct subpopulation of myeloid leukemias, with reduced cell cycle activity and increased resistance to therapeutic challenge. To better characterize key properties of LSCs, we employed a strategy based on identification of genes synergistically dysregulated by cooperating oncogenes. We hypothesized that such genes, termed "cooperation response genes" (CRGs), would represent regulators of LSC growth and survival. Using both a primary mouse model and human leukemia specimens, we show that CRGs comprise genes previously undescribed in leukemia pathogenesis in which multiple pathways modulate the biology of LSCs. In addition, our findings demonstrate that the CRG expression profile can be used as a drug discovery tool for identification of compounds that selectively target the LSC population. We conclude that CRG-based analyses provide a powerful means to characterize the basic biology of LSCs as well as to identify improved methods for therapeutic targeting.

Functional genomics attempts to understand the genome by perturbing the flow of information from DNA to RNA to protein, in order to learn how gene dysfunction leads to disease. CRISPR/Cas9 technology is the newest tool in the geneticist's toolbox, allowing researchers to edit DNA with unprecedented ease, speed and accuracy, and representing a novel means to perform genome-wide genetic screens to discover gene function. In this review, we first summarize the discovery and characterization of CRISPR/Cas9, and then compare it to other genome engineering technologies. We discuss its initial use in screening applications, with a focus on optimizing on-target activity and minimizing off-target effects. Finally, we comment on future challenges and opportunities afforded by this technology.

Huntington's disease, the most common inherited neurodegenerative disease, is characterized by a dramatic loss of deep-layer cortical and striatal neurons, as well as morbidity in midlife. Human genetic studies led to the identification of the causative gene, huntingtin. Recent genomic advances have also led to the identification of hundreds of potential interacting partners for huntingtin protein and many hypotheses as to the molecular mechanisms whereby mutant huntingtin leads to cellular dysfunction and death. However, the multitude of possible interacting partners and cellular pathways affected by mutant huntingtin has complicated efforts to understand the etiology of this disease, and to date no curative therapeutic exists. To address the general problem of identifying the disease-phenotype contributing genes from a large number of correlative studies, here we develop a synthetic lethal screening methodology for the mammalian central nervous system, called SLIC, for synthetic lethal in the central nervous system. Applying SLIC to the study of Huntington's disease, we identify the age-regulated glutathione peroxidase 6 (Gpx6) gene as a modulator of mutant huntingtin toxicity and show that overexpression of Gpx6 can dramatically alleviate both behavioral and molecular phenotypes associated with a mouse model of Huntington's disease. SLIC can, in principle, be used in the study of any neurodegenerative disease for which a mouse model exists, promising to reveal modulators of neurodegenerative disease in an unbiased fashion, akin to screens in simpler model organisms.

Clear-cell carcinomas (CCCs) are a histological group of highly aggressive malignancies commonly originating in the kidney and ovary. CCCs are distinguished by aberrant lipid and glycogen accumulation and are refractory to a broad range of anti-cancer therapies. Here we identify an intrinsic vulnerability to ferroptosis associated with the unique metabolic state in CCCs. This vulnerability transcends lineage and genetic landscape, and can be exploited by inhibiting glutathione peroxidase 4 (GPX4) with small-molecules. Using CRISPR screening and lipidomic profiling, we identify the hypoxia-inducible factor (HIF) pathway as a driver of this vulnerability. In renal CCCs, HIF-2α selectively enriches polyunsaturated lipids, the rate-limiting substrates for lipid peroxidation, by activating the expression of hypoxia-inducible, lipid droplet-associated protein (HILPDA). Our study suggests targeting GPX4 as a therapeutic opportunity in CCCs, and highlights that therapeutic approaches can be identified on the basis of cell states manifested by morphological and metabolic features in hard-to-treat cancers.

The formation of synapses, the specialized points of chemical communication between neurons, is a highly regulated developmental process fundamental to establishing normal brain circuitry. Perturbations of synapse formation and function causally contribute to human developmental and degenerative neuropsychiatric disorders, such as Alzheimer's disease, intellectual disability, and autism spectrum disorders. Many genes controlling synaptogenesis have been identified, but lack of facile experimental systems has made systematic discovery of regulators of synaptogenesis challenging. Thus, we created a high-throughput platform to study excitatory and inhibitory synapse development in primary neuronal cultures and used a lentiviral RNA interference library to identify novel regulators of synapse formation. This methodology is broadly applicable for high-throughput screening of genes and drugs that may rescue or improve synaptic dysfunction associated with cognitive function and neurological disorders.

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.

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.

Many chemotherapeutic drugs are differentially effective from one patient to the next. Understanding the causes of this variability is a critical step towards the development of personalized treatments and improvements to existing medications. Here, we investigate sensitivity to a group of anti-neoplastic drugs that target topoisomerase II using the model organism Caenorhabditis elegans. We show that wild strains of C. elegans vary in their sensitivity to these drugs, and we use an unbiased genetic approach to demonstrate that this natural variation is explained by a methionine-to-glutamine substitution in topoisomerase II (TOP-2). The presence of a non-polar methionine at this residue increases hydrophobic interactions between TOP-2 and its poison etoposide, as compared to a polar glutamine. We hypothesize that this stabilizing interaction results in increased genomic instability in strains that contain a methionine residue. The residue affected by this substitution is conserved from yeast to humans and is one of the few differences between the two human topoisomerase II isoforms (methionine in hTOPIIα and glutamine in hTOPIIβ). We go on to show that this amino acid difference between the two human topoisomerase isoforms influences cytotoxicity of topoisomerase II poisons in human cell lines. These results explain why hTOPIIα and hTOPIIβ are differentially affected by various poisons and demonstrate the utility of C. elegans in understanding the genetics of drug responses.

The CRISPR-Cas9 system has revolutionized gene editing both at single genes and in multiplexed loss-of-function screens, thus enabling precise genome-scale identification of genes essential for proliferation and survival of cancer cells. However, previous studies have reported that a gene-independent antiproliferative effect of Cas9-mediated DNA cleavage confounds such measurement of genetic dependency, thereby leading to false-positive results in copy number-amplified regions. We developed CERES, a computational method to estimate gene-dependency levels from CRISPR-Cas9 essentiality screens while accounting for the copy number-specific effect. In our efforts to define a cancer dependency map, we performed genome-scale CRISPR-Cas9 essentiality screens across 342 cancer cell lines and applied CERES to this data set. We found that CERES decreased false-positive results and estimated sgRNA activity for both this data set and previously published screens performed with different sgRNA libraries. We further demonstrate the utility of this collection of screens, after CERES correction, for identifying cancer-type-specific vulnerabilities.

CRISPR-Cas9 provides the means to perform genome editing and facilitates loss-of-function screens. However, we and others demonstrated that expression of the Cas9 endonuclease induces a gene-independent response that correlates with the number of target sequences in the genome. An alternative approach to suppressing gene expression is to block transcription using a catalytically inactive Cas9 (dCas9). Here we directly compare genome editing by CRISPR-Cas9 (cutting, CRISPRc) and gene suppression using KRAB-dCas9 (CRISPRi) in loss-of-function screens to identify cell essential genes. CRISPRc identified 98% of previously defined cell essential genes. After optimizing library construction by analysing transcriptional start sites (TSS), CRISRPi identified 92% of core cell essential genes and did not show a bias to regions involved in copy number alterations. However, bidirectional promoters scored as false positives in CRISRPi. We conclude that CRISPRc and CRISPRi have different off-target effects and combining these approaches provides complementary information in loss-of-function genetic screens.

DNA damage activates a signalling network that blocks cell-cycle progression, recruits DNA repair factors and/or triggers senescence or programmed cell death. Alterations in chromatin structure are implicated in the initiation and propagation of the DNA damage response. Here we further investigate the role of chromatin structure in the DNA damage response by monitoring ionizing-radiation-induced signalling and response events with a high-content multiplex RNA-mediated interference screen of chromatin-modifying and -interacting genes. We discover that an isoform of Brd4, a bromodomain and extra-terminal (BET) family member, functions as an endogenous inhibitor of DNA damage response signalling by recruiting the condensin II chromatin remodelling complex to acetylated histones through bromodomain interactions. Loss of this isoform results in relaxed chromatin structure, rapid cell-cycle checkpoint recovery and enhanced survival after irradiation, whereas functional gain of this isoform compacted chromatin, attenuated DNA damage response signalling and enhanced radiation-induced lethality. These data implicate Brd4, previously known for its role in transcriptional control, as an insulator of chromatin that can modulate the signalling response to DNA damage.

Deciphering the signaling networks that underlie normal and disease processes remains a major challenge. Here, we report the discovery of signaling components involved in the Toll-like receptor (TLR) response of immune dendritic cells (DCs), including a previously unkown pathway shared across mammalian antiviral responses. By combining transcriptional profiling, genetic and small-molecule perturbations, and phosphoproteomics, we uncover 35 signaling regulators, including 16 known regulators, involved in TLR signaling. In particular, we find that Polo-like kinases (Plk) 2 and 4 are essential components of antiviral pathways in vitro and in vivo and activate a signaling branch involving a dozen proteins, among which is Tnfaip2, a gene associated with autoimmune diseases but whose role was unknown. Our study illustrates the power of combining systematic measurements and perturbations to elucidate complex signaling circuits and discover potential therapeutic targets.

Dysregulation of the phosphatidylinositol 3-kinase (PI3K) signaling pathway occurs frequently in human cancer. PTEN tumor suppressor or PIK3CA oncogene mutations both direct PI3K-dependent tumorigenesis largely through activation of the AKT/PKB kinase. However, here we show through phosphoprotein profiling and functional genomic studies that many PIK3CA mutant cancer cell lines and human breast tumors exhibit only minimal AKT activation and a diminished reliance on AKT for anchorage-independent growth. Instead, these cells retain robust PDK1 activation and membrane localization and exhibit dependency on the PDK1 substrate SGK3. SGK3 undergoes PI3K- and PDK1-dependent activation in PIK3CA mutant cancer cells. Thus, PI3K may promote cancer through both AKT-dependent and AKT-independent mechanisms. Knowledge of differential PI3K/PDK1 signaling could inform rational therapeutics in cancers harboring PIK3CA mutations.

Soft-tissue sarcomas, which result in approximately 10,700 diagnoses and 3,800 deaths per year in the United States, show remarkable histologic diversity, with more than 50 recognized subtypes. However, knowledge of their genomic alterations is limited. We describe an integrative analysis of DNA sequence, copy number and mRNA expression in 207 samples encompassing seven major subtypes. Frequently mutated genes included TP53 (17% of pleomorphic liposarcomas), NF1 (10.5% of myxofibrosarcomas and 8% of pleomorphic liposarcomas) and PIK3CA (18% of myxoid/round-cell liposarcomas, or MRCs). PIK3CA mutations in MRCs were associated with Akt activation and poor clinical outcomes. In myxofibrosarcomas and pleomorphic liposarcomas, we found both point mutations and genomic deletions affecting the tumor suppressor NF1. Finally, we found that short hairpin RNA (shRNA)-based knockdown of several genes amplified in dedifferentiated liposarcoma, including CDK4 and YEATS4, decreased cell proliferation. Our study yields a detailed map of molecular alterations across diverse sarcoma subtypes and suggests potential subtype-specific targets for therapy.

The application of RNA interference (RNAi) to mammalian cells has provided the means to perform phenotypic screens to determine the functions of genes. Although RNAi has revolutionized loss-of-function genetic experiments, it has been difficult to systematically assess the prevalence and consequences of off-target effects. The Connectivity Map (CMAP) represents an unprecedented resource to study the gene expression consequences of expressing short hairpin RNAs (shRNAs). Analysis of signatures for over 13,000 shRNAs applied in 9 cell lines revealed that microRNA (miRNA)-like off-target effects of RNAi are far stronger and more pervasive than generally appreciated. We show that mitigating off-target effects is feasible in these datasets via computational methodologies to produce a consensus gene signature (CGS). In addition, we compared RNAi technology to clustered regularly interspaced short palindromic repeat (CRISPR)-based knockout by analysis of 373 single guide RNAs (sgRNAs) in 6 cells lines and show that the on-target efficacies are comparable, but CRISPR technology is far less susceptible to systematic off-target effects. These results will help guide the proper use and analysis of loss-of-function reagents for the determination of gene function.

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.

Inhibitors that target the receptor tyrosine kinase (RTK)/Ras/mitogen-activated protein kinase (MAPK) pathway have led to clinical responses in lung and other cancers, but some patients fail to respond and in those that do resistance inevitably occurs (Balak et al., 2006; Kosaka et al., 2006; Rudin et al., 2013; Wagle et al., 2011). To understand intrinsic and acquired resistance to inhibition of MAPK signaling, we performed CRISPR-Cas9 gene deletion screens in the setting of BRAF, MEK, EGFR, and ALK inhibition. Loss of KEAP1, a negative regulator of NFE2L2/NRF2, modulated the response to BRAF, MEK, EGFR, and ALK inhibition in BRAF-, NRAS-, KRAS-, EGFR-, and ALK-mutant lung cancer cells. Treatment with inhibitors targeting the RTK/MAPK pathway increased reactive oxygen species (ROS) in cells with intact KEAP1, and loss of KEAP1 abrogated this increase. In addition, loss of KEAP1 altered cell metabolism to allow cells to proliferate in the absence of MAPK signaling. These observations suggest that alterations in the KEAP1/NRF2 pathway may promote survival in the presence of multiple inhibitors targeting the RTK/Ras/MAPK pathway.

The phagocyte oxidative burst, mediated by Nox2 NADPH oxidase-derived reactive oxygen species, confers host defense against a broad spectrum of bacterial and fungal pathogens. Loss-of-function mutations that impair function of the Nox2 complex result in a life-threatening immunodeficiency, and genetic variants of Nox2 subunits have been implicated in pathogenesis of inflammatory bowel disease (IBD). Thus, alterations in the oxidative burst can profoundly impact host defense, yet little is known about regulatory mechanisms that fine-tune this response. Here we report the discovery of regulatory nodes controlling oxidative burst by functional screening of genes within loci linked to human inflammatory disease. Implementing a multi-omics approach, we define transcriptional, metabolic and ubiquitin-cycling nodes controlled by Rbpj, Pfkl and Rnf145, respectively. Furthermore, we implicate Rnf145 in proteostasis of the Nox2 complex by endoplasmic reticulum-associated degradation. Consequently, ablation of Rnf145 in murine macrophages enhances bacterial clearance, and rescues the oxidative burst defects associated with Ncf4 haploinsufficiency.

To accomplish the remarkable task of lifelong infection, the Epstein-Barr virus (EBV) switches between four viral genome latency and lytic programmes to navigate the B-cell compartment and evade immune responses. The transforming programme, consisting of highly immunogenic EBV nuclear antigen (EBNA) and latent membrane proteins (LMPs), is expressed in newly infected B lymphocytes and in post-transplant lymphomas. On memory cell differentiation and in most EBV-associated Burkitt's lymphomas, all but one viral antigen are repressed for immunoevasion. To gain insights into the epigenetic mechanisms that restrict immunogenic oncoprotein expression, a genome-scale CRISPR-Cas9 screen was performed in EBV and Burkitt's lymphoma cells. Here, we show that the ubiquitin ligase ubiquitin-like PHD and RING finger domain-containing protein 1 (UHRF1) and its DNA methyltransferase partner DNA methyltransferase I (DNMT1) are critical for the restriction of EBNA and LMP expression. All UHRF1 reader and writer domains were necessary for silencing and DNMT3B was identified as an upstream viral genome CpG methylation initiator. Polycomb repressive complex I exerted a further layer of control over LMP expression, suggesting a second mechanism for latency programme switching. UHRF1, DNMT1 and DNMT3B are upregulated in germinal centre B cells, the Burkitt's lymphoma cell of origin, providing a molecular link between B-cell state and the EBV latency programme. These results suggest rational therapeutic targets to manipulate EBV oncoprotein expression.

Adoptive cell therapy represents a new paradigm in cancer immunotherapy, but it can be limited by the poor persistence and function of transferred T cells1. Here we use an in vivo pooled CRISPR-Cas9 mutagenesis screening approach to demonstrate that, by targeting REGNASE-1, CD8+ T cells are reprogrammed to long-lived effector cells with extensive accumulation, better persistence and robust effector function in tumours. REGNASE-1-deficient CD8+ T cells show markedly improved therapeutic efficacy against mouse models of melanoma and leukaemia. By using a secondary genome-scale CRISPR-Cas9 screening, we identify BATF as the key target of REGNASE-1 and as a rheostat that shapes antitumour responses. Loss of BATF suppresses the increased accumulation and mitochondrial fitness of REGNASE-1-deficient CD8+ T cells. By contrast, the targeting of additional signalling factors-including PTPN2 and SOCS1-improves the therapeutic efficacy of REGNASE-1-deficient CD8+ T cells. Our findings suggest that T cell persistence and effector function can be coordinated in tumour immunity and point to avenues for improving the efficacy of adoptive cell therapy for cancer.