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.

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.

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.

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.

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.

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.

Genome-wide association studies (GWAS) have identified thousands of variants associated with human diseases and traits. However, the majority of GWAS-implicated variants are in non-coding regions of the genome and require in depth follow-up to identify target genes and decipher biological mechanisms. Here, rather than focusing on causal variants, we have undertaken a pooled loss-of-function screen in primary hematopoietic cells to interrogate 389 candidate genes contained in 75 loci associated with red blood cell traits. Using this approach, we identify 77 genes at 38 GWAS loci, with most loci harboring 1-2 candidate genes. Importantly, the hit set was strongly enriched for genes validated through orthogonal genetic approaches. Genes identified by this approach are enriched in specific and relevant biological pathways, allowing regulators of human erythropoiesis and modifiers of blood diseases to be defined. More generally, this functional screen provides a paradigm for gene-centric follow up of GWAS for a variety of human diseases and traits.

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.

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.

Identification of host genes essential for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection may reveal novel therapeutic targets and inform our understanding of coronavirus disease 2019 (COVID-19) pathogenesis. Here we performed genome-wide CRISPR screens in Vero-E6 cells with SARS-CoV-2, Middle East respiratory syndrome CoV (MERS-CoV), bat CoV HKU5 expressing the SARS-CoV-1 spike, and vesicular stomatitis virus (VSV) expressing the SARS-CoV-2 spike. We identified known SARS-CoV-2 host factors, including the receptor ACE2 and protease Cathepsin L. We additionally discovered pro-viral genes and pathways, including HMGB1 and the SWI/SNF chromatin remodeling complex, that are SARS lineage and pan-coronavirus specific, respectively. We show that HMGB1 regulates ACE2 expression and is critical for entry of SARS-CoV-2, SARS-CoV-1, and NL63. We also show that small-molecule antagonists of identified gene products inhibited SARS-CoV-2 infection in monkey and human cells, demonstrating the conserved role of these genetic hits across species. This identifies potential therapeutic targets for SARS-CoV-2 and reveals SARS lineage-specific and pan-CoV host factors that regulate susceptibility to highly pathogenic CoVs.

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.

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.

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.

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.

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.

Genetic interactions mediate the emergence of phenotype from genotype, but initial technologies for multiplex genetic perturbation in mammalian cells suffer from inefficiency and are challenging to scale. Recent focus on paralog synthetic lethality in cancer cells offers an opportunity to evaluate different CRISPR/Cas multiplexing technologies and improve on the state of the art. Here we report a meta-analysis of CRISPR genetic interactions screens, identifying a candidate set of background-independent paralog synthetic lethals, and find that the CRISPR/enCas12a platform provides superior sensitivity and assay replicability. We demonstrate that enCas12a can independently target up to four genes from a single guide array, and build on this knowledge by constructing a one-component library that expresses arrays of four guides per clone, a platform we call 'in4mer'. Our genome-scale human library, with only 44k clones, is substantially smaller than a typical CRISPR/Cas9 monogenic library while also targeting more than two thousand paralog pairs, triples, and quads. Proof of concept screens in two cell lines demonstrate discrimination of core and context-dependent essential genes similar to that of state of the art CRISPR/Cas9 libraries, as well as detection of synthetic lethal and masking (also known as buffering) genetic interactions between paralogs of various family sizes, a capability not offered by any extant library. Importantly, the in4mer platform offers a fivefold reduction in the number of clones required to assay genetic interactions, dramatically improving the cost and effort required for these studies.

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.

Stimulator of interferon genes (STING) is an intracellular sensor of cyclic di-nucleotides involved in the innate immune response against pathogen- or self-derived DNA. STING trafficking is tightly linked to its function, and its dysregulation can lead to disease. Here, we systematically characterize genes regulating STING trafficking and examine their impact on STING-mediated responses. Using proximity-ligation proteomics and genetic screens, we demonstrate that an endosomal sorting complex required for transport (ESCRT) complex containing HGS, VPS37A and UBAP1 promotes STING degradation, thereby terminating STING-mediated signaling. Mechanistically, STING oligomerization increases its ubiquitination by UBE2N, forming a platform for ESCRT recruitment at the endosome that terminates STING signaling via sorting in the lysosome. Finally, we show that expression of a UBAP1 mutant identified in patients with hereditary spastic paraplegia and associated with disrupted ESCRT function, increases steady-state STING-dependent type I IFN responses in healthy primary monocyte-derived dendritic cells and fibroblasts. Based on these findings, we propose that STING is subject to a tonic degradative flux and that the ESCRT complex acts as a homeostatic regulator of STING signaling.

Recognition of Cas9 and other proteins encoded in delivery vectors has limited CRISPR technology in vivo. Here, we present a protocol for genome engineering using selective CRISPR antigen removal (SCAR) lentiviral vectors in Renca mouse model. This protocol describes how to conduct an in vivo genetic screen with a sgRNA library and SCAR vectors that can be applied to different cell lines and contexts. For complete details on the use and execution of this protocol, please refer to Dubrot et al. (2021).1.