Project DRIVE: A Compendium of Cancer Dependencies and Synthetic Lethal Relationships Uncovered by Large-Scale, Deep RNAi Screening.

Elucidation of the mutational landscape of human cancer has progressed rapidly and been accompanied by the development of therapeutics targeting mutant oncogenes. However, a comprehensive mapping of cancer dependencies has lagged behind and the discovery of therapeutic targets for counteracting tumor suppressor gene loss is needed. To identify vulnerabilities relevant to specific cancer subtypes, we conducted a large-scale RNAi screen in which viability effects of mRNA knockdown were assessed for 7,837 genes using an average of 20 shRNAs per gene in 398 cancer cell lines. We describe findings of this screen, outlining the classes of cancer dependency genes and their relationships to genetic, expression, and lineage features. In addition, we describe robust gene-interaction networks recapitulating both protein complexes and functional cooperation among complexes and pathways. This dataset along with a web portal is provided to the community to assist in the discovery and translation of new therapeutic approaches for cancer.

DNA Repair Profiling Reveals Nonrandom Outcomes at Cas9-Mediated Breaks.

The repair outcomes at site-specific DNA double-strand breaks (DSBs) generated by the RNA-guided DNA endonuclease Cas9 determine how gene function is altered. Despite the widespread adoption of CRISPR-Cas9 technology to induce DSBs for genome engineering, the resulting repair products have not been examined in depth. Here, the DNA repair profiles of 223 sites in the human genome demonstrate that the pattern of DNA repair following Cas9 cutting at each site is nonrandom and consistent across experimental replicates, cell lines, and reagent delivery methods. Furthermore, the repair outcomes are determined by the protospacer sequence rather than genomic context, indicating that DNA repair profiling in cell lines can be used to anticipate repair outcomes in primary cells. Chemical inhibition of DNA-PK enabled dissection of the DNA repair profiles into contributions from c-NHEJ and MMEJ. Finally, this work elucidates a strategy for using "error-prone" DNA-repair machinery to generate precise edits.

Metabolic stress controls mTORC1 lysosomal localization and dimerization by regulating the TTT-RUVBL1/2 complex.

The metabolism of glucose and glutamine, primary carbon sources utilized by mitochondria to generate energy and macromolecules for cell growth, is directly regulated by mTORC1. We show that glucose and glutamine, by supplying carbons to the TCA cycle to produce ATP, positively feed back to mTORC1 through an AMPK-, TSC1/2-, and Rag-independent mechanism by regulating mTORC1 assembly and its lysosomal localization. We discovered that the ATP-dependent TTT-RUVBL1/2 complex was disassembled and repressed by energy depletion, resulting in its decreased interaction with mTOR. The TTT-RUVBL complex was necessary for the interaction between mTORC1 and Rag and formation of mTORC1 obligate dimers. In cancer tissues, TTT-RUVBL complex mRNAs were elevated and positively correlated with transcripts encoding proteins of anabolic metabolism and mitochondrial function-all mTORC1-regulated processes. Thus, the TTT-RUVBL1/2 complex responds to the cell's metabolic state, directly regulating the functional assembly of mTORC1 and indirectly controlling the nutrient signal from Rags to mTORC1.

Genome-wide CRISPR screen for PARKIN regulators reveals transcriptional repression as a determinant of mitophagy.

PARKIN, an E3 ligase mutated in familial Parkinson's disease, promotes mitophagy by ubiquitinating mitochondrial proteins for efficient engagement of the autophagy machinery. Specifically, PARKIN-synthesized ubiquitin chains represent targets for the PINK1 kinase generating phosphoS65-ubiquitin (pUb), which constitutes the mitophagy signal. Physiological regulation of PARKIN abundance, however, and the impact on pUb accumulation are poorly understood. Using cells designed to discover physiological regulators of PARKIN abundance, we performed a pooled genome-wide CRISPR/Cas9 knockout screen. Testing identified genes individually resulted in a list of 53 positive and negative regulators. A transcriptional repressor network including THAP11 was identified and negatively regulates endogenous PARKIN abundance. RNAseq analysis revealed the PARKIN-encoding locus as a prime THAP11 target, and THAP11 CRISPR knockout in multiple cell types enhanced pUb accumulation. Thus, our work demonstrates the critical role of PARKIN abundance, identifies regulating genes, and reveals a link between transcriptional repression and mitophagy, which is also apparent in human induced pluripotent stem cell-derived neurons, a disease-relevant cell type.

TMEM41B is a novel regulator of autophagy and lipid mobilization.

Autophagy maintains cellular homeostasis by targeting damaged organelles, pathogens, or misfolded protein aggregates for lysosomal degradation. The autophagic process is initiated by the formation of autophagosomes, which can selectively enclose cargo via autophagy cargo receptors. A machinery of well-characterized autophagy-related proteins orchestrates the biogenesis of autophagosomes; however, the origin of the required membranes is incompletely understood. Here, we have applied sensitized pooled CRISPR screens and identify the uncharacterized transmembrane protein TMEM41B as a novel regulator of autophagy. In the absence of TMEM41B, autophagosome biogenesis is stalled, LC3 accumulates at WIPI2- and DFCP1-positive isolation membranes, and lysosomal flux of autophagy cargo receptors and intracellular bacteria is impaired. In addition to defective autophagy, TMEM41B knockout cells display significantly enlarged lipid droplets and reduced mobilization and ß-oxidation of fatty acids. Immunostaining and interaction proteomics data suggest that TMEM41B localizes to the endoplasmic reticulum (ER). Taken together, we propose that TMEM41B is a novel ER-localized regulator of autophagosome biogenesis and lipid mobilization.

Macroautophagy is dispensable for growth of KRAS mutant tumors and chloroquine efficacy.

Macroautophagy is a key stress-response pathway that can suppress or promote tumorigenesis depending on the cellular context. Notably, Kirsten rat sarcoma (KRAS)-driven tumors have been reported to rely on macroautophagy for growth and survival, suggesting a potential therapeutic approach of using autophagy inhibitors based on genetic stratification. In this study, we evaluated whether KRAS mutation status can predict the efficacy to macroautophagy inhibition. By profiling 47 cell lines with pharmacological and genetic loss-of-function tools, we were unable to confirm that KRAS-driven tumor lines require macroautophagy for growth. Deletion of autophagy-related 7 (ATG7) by genome editing completely blocked macroautophagy in several tumor lines with oncogenic mutations in KRAS but did not inhibit cell proliferation in vitro or tumorigenesis in vivo. Furthermore, ATG7 knockout did not sensitize cells to irradiation or to several anticancer agents tested. Interestingly, ATG7-deficient and -proficient cells were equally sensitive to the antiproliferative effect of chloroquine, a lysosomotropic agent often used as a pharmacological tool to evaluate the response to macroautophagy inhibition. Moreover, both cell types manifested synergistic growth inhibition when treated with chloroquine plus the tyrosine kinase inhibitors erlotinib or sunitinib, suggesting that the antiproliferative effects of chloroquine are independent of its suppressive actions on autophagy.

Synthetic lethality of combined glutaminase and Hsp90 inhibition in mTORC1-driven tumor cells.

The mammalian target of rapamycin complex 1 (mTORC1) integrates multiple signals from growth factors, nutrients, and cellular energy status to control a wide range of metabolic processes, including mRNA biogenesis; protein, nucleotide, and lipid synthesis; and autophagy. Deregulation of the mTORC1 pathway is found in cancer as well as genetic disorders such as tuberous sclerosis complex (TSC) and sporadic lymphangioleiomyomatosis. Recent studies have shown that the mTORC1 inhibitor rapamycin and its analogs generally suppress proliferation rather than induce apoptosis. Therefore, it is critical to use alternative strategies to induce death of cells with activated mTORC1. In this study, a small-molecule screen has revealed that the combination of glutaminase (GLS) and heat shock protein 90 (Hsp90) inhibitors selectively triggers death of TSC2-deficient cells. At a mechanistic level, high mTORC1-driven translation rates in TSC1/2-deficient cells, unlike wild-type cells, sensitizes these cells to endoplasmic reticulum (ER) stress. Thus, Hsp90 inhibition drives accumulation of unfolded protein and ER stress. When combining proteotoxic stress with oxidative stress by depletion of the intracellular antioxidant glutathione by GLS inhibition, acute cell death is observed in cells with activated mTORC1 signaling. This study suggests that this combination strategy may have the potential to be developed into a therapeutic use for the treatment of mTORC1-driven tumors.

Functional epigenetics approach identifies BRM/SMARCA2 as a critical synthetic lethal target in BRG1-deficient cancers.

Defects in epigenetic regulation play a fundamental role in the development of cancer, and epigenetic regulators have recently emerged as promising therapeutic candidates. We therefore set out to systematically interrogate epigenetic cancer dependencies by screening an epigenome-focused deep-coverage design shRNA (DECODER) library across 58 cancer cell lines. This screen identified BRM/SMARCA2, a DNA-dependent ATPase of the mammalian SWI/SNF (mSWI/SNF) chromatin remodeling complex, as being essential for the growth of tumor cells that harbor loss of function mutations in BRG1/SMARCA4. Depletion of BRM in BRG1-deficient cancer cells leads to a cell cycle arrest, induction of senescence, and increased levels of global H3K9me3. We further demonstrate the selective dependency of BRG1-mutant tumors on BRM in vivo. Genetic alterations of the mSWI/SNF chromatin remodeling complexes are the most frequent among chromatin regulators in cancers, with BRG1/SMARCA4 mutations occurring in ∼10-15% of lung adenocarcinomas. Our findings position BRM as an attractive therapeutic target for BRG1 mutated cancers. Because BRG1 and BRM function as mutually exclusive catalytic subunits of the mSWI/SNF complex, we propose that such synthetic lethality may be explained by paralog insufficiency, in which loss of one family member unveils critical dependence on paralogous subunits. This concept of "cancer-selective paralog dependency" may provide a more general strategy for targeting other tumor suppressor lesions/complexes with paralogous subunits.

Genome-wide CRISPR screen identifies protein pathways modulating tau protein levels in neurons.

Aggregates of hyperphosphorylated tau protein are a pathological hallmark of more than 20 distinct neurodegenerative diseases, including Alzheimer's disease, progressive supranuclear palsy, and frontotemporal dementia. While the exact mechanism of tau aggregation is unknown, the accumulation of aggregates correlates with disease progression. Here we report a genome-wide CRISPR screen to identify modulators of endogenous tau protein for the first time. Primary screens performed in SH-SY5Y cells, identified positive and negative regulators of tau protein levels. Hit validation of the top 43 candidate genes was performed using Ngn2-induced human cortical excitatory neurons. Using this approach, genes and pathways involved in modulation of endogenous tau levels were identified, including chromatin modifying enzymes, neddylation and ubiquitin pathway members, and components of the mTOR pathway. TSC1, a critical component of the mTOR pathway, was further validated in vivo, demonstrating the relevance of this screening strategy. These findings may have implications for treating neurodegenerative diseases in the future.

CRISPR/Cas9-Based Chemogenomic Profiling in Mammalian Cells.

Chemogenomic profiling is a powerful and unbiased approach to elucidate pharmacological targets and the mechanism of bioactive compounds. It is based on identifying cellular hypersensitivity and resistance caused by individual gene modulations with genome-wide coverage. Due to the requirement of bar-coded, genome-wide deletion collections, high-resolution experiments of this nature have historically been limited to fungal systems. Pooled RNAi reagents have enabled similar attempts in mammalian cells but efforts have been hampered by significant off-target effects and experimental noise. The CRISPR/Cas9 system for the first time enables precise DNA editing at defined loci in a genome-wide fashion. Here we present the detailed protocol that leverages the CRISPR/Cas9 system for chemogenomic profiling and target identification of diverse chemical probes.

Genome-Scale CRISPR Screens Identify Human Pluripotency-Specific Genes.

Human pluripotent stem cells (hPSCs) generate a variety of disease-relevant cells that can be used to improve the translation of preclinical research. Despite the potential of hPSCs, their use for genetic screening has been limited by technical challenges. We developed a scalable and renewable Cas9 and sgRNA-hPSC library in which loss-of-function mutations can be induced at will. Our inducible mutant hPSC library can be used for multiple genome-wide CRISPR screens in a variety of hPSC-induced cell types. As proof of concept, we performed three screens for regulators of properties fundamental to hPSCs: their ability to self-renew and/or survive (fitness), their inability to survive as single-cell clones, and their capacity to differentiate. We identified the majority of known genes and pathways involved in these processes, as well as a plethora of genes with unidentified roles. This resource will increase the understanding of human development and genetics. This approach will be a powerful tool to identify disease-modifying genes and pathways.

YAP, but Not RSPO-LGR4/5, Signaling in Biliary Epithelial Cells Promotes a Ductular Reaction in Response to Liver Injury.

Biliary epithelial cells (BECs) form bile ducts in the liver and are facultative liver stem cells that establish a ductular reaction (DR) to support liver regeneration following injury. Liver damage induces periportal LGR5+ putative liver stem cells that can form BEC-like organoids, suggesting that RSPO-LGR4/5-mediated WNT/ß-catenin activity is important for a DR. We addressed the roles of this and other signaling pathways in a DR by performing a focused CRISPR-based loss-of-function screen in BEC-like organoids, followed by in vivo validation and single-cell RNA sequencing. We found that BECs lack and do not require LGR4/5-mediated WNT/ß-catenin signaling during a DR, whereas YAP and mTORC1 signaling are required for this process. Upregulation of AXIN2 and LGR5 is required in hepatocytes to enable their regenerative capacity in response to injury. Together, these data highlight heterogeneity within the BEC pool, delineate signaling pathways involved in a DR, and clarify the identity and roles of injury-induced periportal LGR5+ cells.

Changes in Sef levels influence auditory brainstem development and function.

During development of the CNS, secreted morphogens of the fibroblast growth factor (FGF) family have multiple effects on cell division, migration, and survival depending on where, when, and how much FGF signal is received. The consequences of misregulating the FGF pathway were studied in a mouse with decreased levels of the FGF antagonist Sef. To uncover effects in the nervous system, we focused on the auditory system, which is accessible to physiological analysis. We found that the mitogen-activated protein kinase pathway is active in the rhombic lip, a germinal zone that generates diverse types of neurons, including the cochlear nucleus complex of the auditory system. Sef is expressed immediately adjacent to the rhombic lip, overlapping with FGF15 and FGFR1, which is also present in the lip itself. This pattern suggests that Sef may normally function in non-rhombic lip cells and prevent them from responding to FGF ligand in the vicinity. Consistent with this idea, overexpression of Sef in chicks decreased the size of the auditory nuclei. Cochlear nucleus defects were also apparent in mice with reduced levels of Sef, with 13% exhibiting grossly dysmorphic cochlear nuclei and 26% showing decreased amounts of GFAP in the cochlear nucleus. Additional evidence for cochlear nucleus defects was obtained by electrophysiological analysis of Sef mutant mice, which have normal auditory thresholds but abnormal auditory brainstem responses. These results show both increases and decreases in Sef levels affect the assembly and function of the auditory brainstem.

Mutational analysis of betaCOP (Sec26p) identifies an appendage domain critical for function.

BACKGROUND: The appendage domain of the gammaCOP subunit of the COPI vesicle coat bears a striking structural resemblance to adaptin-family appendages despite limited primary sequence homology. Both the gammaCOP appendage domain and an equivalent region on betaCOP contain the FxxxW motif; the conservation of this motif suggested the existence of a functional appendage domain in betaCOP. RESULTS: Sequence comparisons in combination with structural prediction tools show that the fold of the COOH-terminus of Sec26p is strongly predicted to closely mimic that of adaptin-family appendages. Deletion of the appendage domain of Sec26p results in inviability in yeast, over-expression of the deletion construct is dominant negative and mutagenesis of this region identifies residues critical for function. The ArfGAP Glo3p was identified via suppression screening as a potential downstream modulator of Sec26p in a manner that is independent of the GAP activity of Glo3p but requires the presence of the COOH-terminal ISS motifs. CONCLUSION: Together, these results indicate an essential function for the predicted betaCOP appendage and suggest that both COPI appendages perform a biologically active regulatory role with a structure related to adaptin-family appendage domains.

TRRAP is a central regulator of human multiciliated cell formation.

The multiciliated cell (MCC) is an evolutionarily conserved cell type, which in vertebrates functions to promote directional fluid flow across epithelial tissues. In the conducting airway, MCCs are generated by basal stem/progenitor cells and act in concert with secretory cells to perform mucociliary clearance to expel pathogens from the lung. Studies in multiple systems, including Xenopus laevis epidermis, murine trachea, and zebrafish kidney, have uncovered a transcriptional network that regulates multiple steps of multiciliogenesis, ultimately leading to an MCC with hundreds of motile cilia extended from their apical surface, which beat in a coordinated fashion. Here, we used a pool-based short hairpin RNA screening approach and identified TRRAP, an essential component of multiple histone acetyltransferase complexes, as a central regulator of MCC formation. Using a combination of immunofluorescence, signaling pathway modulation, and genomic approaches, we show that (a) TRRAP acts downstream of the Notch2-mediated basal progenitor cell fate decision and upstream of Multicilin to control MCC differentiation; and (b) TRRAP binds to the promoters and regulates the expression of a network of genes involved in MCC differentiation and function, including several genes associated with human ciliopathies.

p53 inhibits CRISPR-Cas9 engineering in human pluripotent stem cells.

CRISPR/Cas9 has revolutionized our ability to engineer genomes and conduct genome-wide screens in human cells1-3. Whereas some cell types are amenable to genome engineering, genomes of human pluripotent stem cells (hPSCs) have been difficult to engineer, with reduced efficiencies relative to tumour cell lines or mouse embryonic stem cells3-13. Here, using hPSC lines with stable integration of Cas9 or transient delivery of Cas9-ribonucleoproteins (RNPs), we achieved an average insertion or deletion (indel) efficiency greater than 80%. This high efficiency of indel generation revealed that double-strand breaks (DSBs) induced by Cas9 are toxic and kill most hPSCs. In previous studies, the toxicity of Cas9 in hPSCs was less apparent because of low transfection efficiency and subsequently low DSB induction3. The toxic response to DSBs was P53/TP53-dependent, such that the efficiency of precise genome engineering in hPSCs with a wild-type P53 gene was severely reduced. Our results indicate that Cas9 toxicity creates an obstacle to the high-throughput use of CRISPR/Cas9 for genome engineering and screening in hPSCs. Moreover, as hPSCs can acquire P53 mutations14, cell replacement therapies using CRISPR/Cas9-enginereed hPSCs should proceed with caution, and such engineered hPSCs should be monitored for P53 function.

Identification of a novel NAMPT inhibitor by CRISPR/Cas9 chemogenomic profiling in mammalian cells.

Chemogenomic profiling is a powerful and unbiased approach to elucidate pharmacological targets and the mechanism of bioactive compounds. Until recently, genome-wide, high-resolution experiments of this nature have been limited to fungal systems due to lack of mammalian genome-wide deletion collections. With the example of a novel nicotinamide phosphoribosyltransferase (NAMPT) inhibitor, we demonstrate that the CRISPR/Cas9 system enables the generation of transient homo- and heterozygous deletion libraries and allows for the identification of efficacy targets and pathways mediating hypersensitivity and resistance relevant to the compound mechanism of action.

Functional CRISPR screening identifies the ufmylation pathway as a regulator of SQSTM1/p62.

SQSTM1 is an adaptor protein that integrates multiple cellular signaling pathways and whose expression is tightly regulated at the transcriptional and post-translational level. Here, we describe a forward genetic screening paradigm exploiting CRISPR-mediated genome editing coupled to a cell selection step by FACS to identify regulators of SQSTM1. Through systematic comparison of pooled libraries, we show that CRISPR is superior to RNAi in identifying known SQSTM1 modulators. A genome-wide CRISPR screen exposed MTOR signalling and the entire macroautophagy machinery as key regulators of SQSTM1 and identified several novel modulators including HNRNPM, SLC39A14, SRRD, PGK1 and the ufmylation cascade. We show that ufmylation regulates SQSTM1 by eliciting a cell type-specific ER stress response which induces SQSTM1 expression and results in its accumulation in the cytosol. This study validates pooled CRISPR screening as a powerful method to map the repertoire of cellular pathways that regulate the fate of an individual target protein.

Disordered methionine metabolism in MTAP/CDKN2A-deleted cancers leads to dependence on PRMT5.

5-Methylthioadenosine phosphorylase (MTAP) is a key enzyme in the methionine salvage pathway. The MTAP gene is frequently deleted in human cancers because of its chromosomal proximity to the tumor suppressor gene CDKN2A. By interrogating data from a large-scale short hairpin RNA-mediated screen across 390 cancer cell line models, we found that the viability of MTAP-deficient cancer cells is impaired by depletion of the protein arginine methyltransferase PRMT5. MTAP-deleted cells accumulate the metabolite methylthioadenosine (MTA), which we found to inhibit PRMT5 methyltransferase activity. Deletion of MTAP in MTAP-proficient cells rendered them sensitive to PRMT5 depletion. Conversely, reconstitution of MTAP in an MTAP-deficient cell line rescued PRMT5 dependence. Thus, MTA accumulation in MTAP-deleted cancers creates a hypomorphic PRMT5 state that is selectively sensitized toward further PRMT5 inhibition. Inhibitors of PRMT5 that leverage this dysregulated metabolic state merit further investigation as a potential therapy for MTAP/CDKN2A-deleted tumors.

Signaling to the Rho GTPases: networking with the DH domain.

The Dbl homology (DH) domain was first identified in the Dbl oncogene product as the limit region required for mediating guanine nucleotide exchange on the Rho family GTPase Cdc42. Since the initial biochemical characterization of the DH domain, this conserved motif has been identified in a large family of proteins. In each case, a pleckstrin homology (PH) domain immediately follows the DH domain and this tandem DH-PH module is the signature motif of the Dbl family of guanine nucleotide exchange factors (GEFs). Recent structural studies have provided significant insight into the molecular basis of guanine nucleotide exchange by Dbl family GEFs, opening the door for understanding the specificity of the DH/GTPase interaction as well as providing a starting point for understanding how the exchange activity of these proteins is modulated to achieve specific biological outcomes in the cell.