Survival-based CRISPR genetic screens across a panel of permissive cell lines identify common and cell-specific SARS-CoV-2 host factors

SARS-CoV-2 depends on host cell components for infection and replication. Identification of virus-host dependencies offers an effective way to elucidate mechanisms involved in viral infection and replication. If druggable, host factor dependencies may present an attractive strategy for anti-viral therapy. In this study, we performed genome wide CRISPR knockout screens in Vero E6 cells and four human cell lines including Calu-3, UM-UC-4, HEK-293 and HuH-7 to identify genetic regulators of SARS-CoV-2 infection. Our findings identified only ACE2, the cognate SARS-CoV-2 entry receptor, as a common host dependency factor across all cell lines, while other host genes identified were largely cell line specific, including known factors TMPRSS2 and CTSL. Several of the discovered host-dependency factors converged on pathways involved in cell signalling, immune-related pathways, and chromatin modification. Notably, the chromatin modifier gene KMT2C in Calu-3 cells had the strongest impact in preventing SARS-CoV-2 infection when perturbed.

Cancer-selective metabolic vulnerabilities in MYC-amplified medulloblastoma

MYC-driven medulloblastoma (MB) is an aggressive pediatric brain tumor characterized by therapy resistance and disease recurrence. Here, we integrated data from unbiased genetic screening and metabolomic profiling to identify multiple cancer-selective metabolic vulnerabilities in MYC-driven MB tumor cells, which are amenable to therapeutic targeting. Among these targets, dihydroorotate dehydrogenase (DHODH), an enzyme that catalyzes de novo pyrimidine biosynthesis, emerged as a favorable candidate for therapeutic targeting. Mechanistically, DHODH inhibition acts on target, leading to uridine metabolite scarcity and hyperlipidemia, accompanied by reduced protein O-GlcNAcylation and c-Myc degradation. Pyrimidine starvation evokes a metabolic stress response that leads to cell-cycle arrest and apoptosis. We further show that an orally available small-molecule DHODH inhibitor demonstrates potent mono-therapeutic efficacy against patient-derived MB xenografts in vivo. The reprogramming of pyrimidine metabolism in MYC-driven medulloblastoma represents an unappreciated therapeutic strategy and a potential new class of treatments with stronger cancer selectivity and fewer neurotoxic sequelae.

Functional genomic landscape of cancer-intrinsic evasion of killing by T cells

The genetic circuits that allow cancer cells to evade destruction by the host immune system remain poorly understood^1-3. Here, to identify a phenotypically robust core set of genes and pathways that enable cancer cells to evade killing mediated by cytotoxic T lymphocytes (CTLs), we performed genome-wide CRISPR screens across a panel of genetically diverse mouse cancer cell lines that were cultured in the presence of CTLs. We identify a core set of 182 genes across these mouse cancer models, the individual perturbation of which increases either the sensitivity or the resistance of cancer cells to CTL-mediated toxicity. Systematic exploration of our dataset using genetic co-similarity reveals the hierarchical and coordinated manner in which genes and pathways act in cancer cells to orchestrate their evasion of CTLs, and shows that discrete functional modules that control the interferon response and tumour necrosis factor (TNF)-induced cytotoxicity are dominant sub-phenotypes. Our data establish a central role for genes that were previously identified as negative regulators of the type-II interferon response (for example, Ptpn2, Socs1 and Adar1) in mediating CTL evasion, and show that the lipid-droplet-related gene Fitm2 is required for maintaining cell fitness after exposure to interferon-γ (IFNγ). In addition, we identify the autophagy pathway as a conserved mediator of the evasion of CTLs by cancer cells, and show that this pathway is required to resist cytotoxicity induced by the cytokines IFNγ and TNF. Through the mapping of cytokine- and CTL-based genetic interactions, together with in vivo CRISPR screens, we show how the pleiotropic effects of autophagy control cancer-cell-intrinsic evasion of killing by CTLs and we highlight the importance of these effects within the tumour microenvironment. Collectively, these data expand our knowledge of the genetic circuits that are involved in the evasion of the immune system by cancer cells, and highlight genetic interactions that contribute to phenotypes associated with escape from killing by CTLs.

Systematic mapping of genetic interactions for de novo fatty acid synthesis identifies C12orf49 as a regulator of lipid metabolism

The de novo synthesis of fatty acids has emerged as a therapeutic target for various diseases, including cancer. Because cancer cells are intrinsically buffered to combat metabolic stress, it is important to understand how cells may adapt to the loss of de novo fatty acid biosynthesis. Here, we use pooled genome-wide CRISPR screens to systematically map genetic interactions (GIs) in human HAP1 cells carrying a loss-of-function mutation in fatty acid synthase (FASN), whose product catalyses the formation of long-chain fatty acids. FASN-mutant cells show a strong dependence on lipid uptake that is reflected in negative GIs with genes involved in the LDL receptor pathway, vesicle trafficking and protein glycosylation. Further support for these functional relationships is derived from additional GI screens in query cell lines deficient in other genes involved in lipid metabolism, including LDLR, SREBF1, SREBF2 and ACACA. Our GI profiles also identify a potential role for the previously uncharacterized gene C12orf49 (which we call LUR1) in regulation of exogenous lipid uptake through modulation of SREBF2 signalling in response to lipid starvation. Overall, our data highlight the genetic determinants underlying the cellular adaptation associated with loss of de novo fatty acid synthesis and demonstrate the power of systematic GI mapping for uncovering metabolic buffering mechanisms in human cells.

Forward genetic screen in human podocytes identifies diphthamide biosynthesis genes as regulators of adhesion

Podocyte function is tightly linked to the complex organization of its cytoskeleton and adhesion to the underlying glomerular basement membrane. Adhesion of cultured podocytes to a variety of substrates is reported to correlate with podocyte health. To identify novel genes that are important for podocyte function, we designed an in vitro genetic screen based on podocyte adhesion to plates coated with either fibronectin or soluble Fms-like tyrosine kinase-1 (sFLT1)/Fc. A genome-scale pooled RNA interference screen on immortalized human podocytes identified 77 genes that increased adhesion to fibronectin, 101 genes that increased adhesion to sFLT1/Fc, and 44 genes that increased adhesion to both substrates when knocked down. Multiple shRNAs against diphthamide biosynthesis protein 1-4 (DPH1-DPH4) were top hits for increased adhesion. Immortalized human podocyte cells stably expressing these hairpins displayed increased adhesion to both substrates. We then used CRISPR-Cas9 to generate podocyte knockout cells for DPH1, DPH2, or DPH3, which also displayed increased adhesion to both fibronectin and sFLT1/Fc, as well as a spreading defect. Finally, we showed that Drosophila nephrocyte-specific knockdown of Dph1, Dph2, and Dph4 resulted in altered nephrocyte function. In summary, we report here a novel high-throughput method to identify genes important for podocyte function. Given the central role of podocyte adhesion as a marker of podocyte health, these data are a rich source of candidate regulators of glomerular disease.

Pooled CRISPR-Based Genetic Screens in Mammalian Cells

Genome editing using the CRISPR-Cas system has vastly advanced the ability to precisely edit the genomes of various organisms. In the context of mammalian cells, this technology represents a novel means to perform genome-wide genetic screens for functional genomics studies. Libraries of guide RNAs (sgRNA) targeting all open reading frames permit the facile generation of thousands of genetic perturbations in a single pool of cells that can be screened for specific phenotypes to implicate gene function and cellular processes in an unbiased and systematic way. CRISPR-Cas screens provide researchers with a simple, efficient, and inexpensive method to uncover the genetic blueprints for cellular phenotypes. Furthermore, differential analysis of screens performed in various cell lines and from different cancer types can identify genes that are contextually essential in tumor cells, revealing potential targets for specific anticancer therapies. Performing genome-wide screens in human cells can be daunting, as this involves the handling of tens of millions of cells and requires analysis of large sets of data. The details of these screens, such as cell line characterization, CRISPR library considerations, and understanding the limitations and capabilities of CRISPR technology during analysis, are often overlooked. Provided here is a detailed protocol for the successful performance of pooled genome-wide CRISPR-Cas9 based screens.

Evaluation and Design of Genome-Wide CRISPR/SpCas9 Knockout Screens

The adaptation of CRISPR/SpCas9 technology to mammalian cell lines is transforming the study of human functional genomics. Pooled libraries of CRISPR guide RNAs (gRNAs) targeting human protein-coding genes and encoded in viral vectors have been used to systematically create gene knockouts in a variety of human cancer and immortalized cell lines, in an effort to identify whether these knockouts cause cellular fitness defects. Previous work has shown that CRISPR screens are more sensitive and specific than pooled-library shRNA screens in similar assays, but currently there exists significant variability across CRISPR library designs and experimental protocols. In this study, we reanalyze 17 genome-scale knockout screens in human cell lines from three research groups, using three different genome-scale gRNA libraries. Using the Bayesian Analysis of Gene Essentiality algorithm to identify essential genes, we refine and expand our previously defined set of human core essential genes from 360 to 684 genes. We use this expanded set of reference core essential genes, CEG2, plus empirical data from six CRISPR knockout screens to guide the design of a sequence-optimized gRNA library, the Toronto KnockOut version 3.0 (TKOv3) library. We then demonstrate the high effectiveness of the library relative to reference sets of essential and nonessential genes, as well as other screens using similar approaches. The optimized TKOv3 library, combined with the CEG2 reference set, provide an efficient, highly optimized platform for performing and assessing gene knockout screens in human cell lines.

High-Resolution CRISPR Screens Reveal Fitness Genes and Genotype-Specific Cancer Liabilities

The ability to perturb genes in human cells is crucial for elucidating gene function and holds great potential for finding therapeutic targets for diseases such as cancer. To extend the catalog of human core and context-dependent fitness genes, we have developed a high-complexity second-generation genome-scale CRISPR-Cas9 gRNA library and applied it to fitness screens in five human cell lines. Using an improved Bayesian analytical approach, we consistently discover 5-fold more fitness genes than were previously observed. We present a list of 1,580 human core fitness genes and describe their general properties. Moreover, we demonstrate that context-dependent fitness genes accurately recapitulate pathway-specific genetic vulnerabilities induced by known oncogenes and reveal cell-type-specific dependencies for specific receptor tyrosine kinases, even in oncogenic KRAS backgrounds. Thus, rigorous identification of human cell line fitness genes using a high-complexity CRISPR-Cas9 library affords a high-resolution view of the genetic vulnerabilities of a cell.

High-Resolution CRISPR Screens Reveal Fitness Genes and Genotype-Specific Cancer Liabilities

The ability to perturb genes in human cells is crucial for elucidating gene function and holds great potential for finding therapeutic targets for diseases such as cancer. To extend the catalog of human core and context-dependent fitness genes, we have developed a high-complexity second-generation genome-scale CRISPR-Cas9 gRNA library and applied it to fitness screens in five human cell lines. Using an improved Bayesian analytical approach, we consistently discover 5-fold more fitness genes than were previously observed. We present a list of 1,580 human core fitness genes and describe their general properties. Moreover, we demonstrate that context-dependent fitness genes accurately recapitulate pathway-specific genetic vulnerabilities induced by known oncogenes and reveal cell-type-specific dependencies for specific receptor tyrosine kinases, even in oncogenic KRAS backgrounds. Thus, rigorous identification of human cell line fitness genes using a high-complexity CRISPR-Cas9 library affords a high-resolution view of the genetic vulnerabilities of a cell.

Essential gene profiles in breast, pancreatic, and ovarian cancer cells

Genomic analyses are yielding a host of new information on the multiple genetic abnormalities associated with specific types of cancer. A comprehensive description of cancer-associated genetic abnormalities can improve our ability to classify tumors into clinically relevant subgroups and, on occasion, identify mutant genes that drive the cancer phenotype ("drivers"). More often, though, the functional significance of cancer-associated mutations is difficult to discern. Genome-wide pooled short hairpin RNA (shRNA) screens enable global identification of the genes essential for cancer cell survival and proliferation, providing a "functional genomic" map of human cancer to complement genomic studies. Using a lentiviral shRNA library targeting ~16,000 genes and a newly developed, dynamic scoring approach, we identified essential gene profiles in 72 breast, pancreatic, and ovarian cancer cell lines. Integrating our results with current and future genomic data should facilitate the systematic identification of drivers, unanticipated synthetic lethal relationships, and functional vulnerabilities of these tumor types.

SIGNIFICANCE

This study presents a resource of genome-scale, pooled shRNA screens for 72 breast, pancreatic, and ovarian cancer cell lines that will serve as a functional complement to genomics data, facilitate construction of essential gene profiles, help uncover synthetic lethal relationships, and identify uncharacterized genetic vulnerabilities in these tumor types.

SIGNIFICANCE

This study presents a resource of genome-scale, pooled shRNA screens for 72 breast, pancreatic, and ovarian cancer cell lines that will serve as a functional complement to genomics data, facilitate construction of essential gene profiles, help uncover synthetic lethal relationships, and identify uncharacterized genetic vulnerabilities in these tumor types.

Differential recruitment of activating and inhibitory Fc gamma RII during phagocytosis

Human myeloid cells express both activating and inhibitory receptors of the FcgammaRII family. FcgammaRIIA mediates processes associated with cell activation, including phagocytosis of IgG-opsonized particles, whereas coengagement of the inhibitory FcgammaRIIB downregulates such signaling. We analyzed the relative recruitment of these two receptors during phagocytosis of IgG-coated particles by ts20 Chinese hamster fibroblast cells cotransfected with both receptors carrying distinguishable fluorescent protein tags. We found that FcgammaRIIA is substantially enriched at sites of particle binding relative to its inhibitory counterpart, with a greater than 2-fold increase in the local ratio of activating to inhibitory receptor compared with that for the plasma membrane as a whole. Experiments with chimeric receptors revealed that the preferential enrichment of FcgammaRIIA results from differences between the extracellular domains of the receptors, and indicated that the lesser recruitment of FcgammaRIIB limits its ability to effectively inhibit FcgammaRIIA-mediated phagocytosis. Mutagenesis studies indicated that FcgammaRIIA residues leucine 132 and phenylalanine 160, which lie in IgG-binding regions of FcgammaRIIA and which differ in FcgammaRIIB, both contribute to the local relative enrichment of FcgammaRIIA by increasing its affinity for IgG1 relative to that of FcgammaRIIB. In human monocytes, engagement of approximately equal amounts of FcgammaRIIB was required to substantially inhibit FcgammaRIIA-mediated phagocytosis. These results demonstrate that differences in affinity for IgG between activating and inhibitory FcgammaR can result in substantial local changes in their relative concentrations during phagocytosis, with important functional consequences.

Prolonging microtubule dysruption enhances the immunogenicity of chronic lymphocytic leukaemia cells

Cytotoxic chemotherapies do not usually mediate the expression of an immunogenic gene programme in tumours, despite activating many of the signalling pathways employed by highly immunogenic cells. Concomitant use of agents that modulate and complement stress-signalling pathways activated by chemotherapeutic agents may then enhance the immunogenicity of cancer cells, increase their susceptibility to T cell-mediated controls and lead to higher clinical remission rates. Consistent with this hypothesis, the microtubule inhibitor, vincristine, caused chronic lymphocytic leukaemia (CLL) cells to die rapidly, without increasing their immunogenicity. Protein kinase C (PKC) agonists (such as bryostatin) delayed the death of vincristine-treated CLL cells and made them highly immunogenic, with increased stimulatory abilities in mixed lymphocyte responses, production of proinflammatory cytokines, expression of co-stimulatory molecules and activation of c-Jun N-terminal kinase (JNK), p38 and nuclear factor kappa B (NF-kappaB) signalling pathways. This phenotype was similar to the result of activating CLL cells through Toll-like receptors (TLRs), which communicate 'danger' signals from infectious pathogens. Use of PKC agonists and microtubule inhibitors to mimic TLR-signalling, and increase the immunogenicity of CLL cells, has implications for the design of chemo-immunotherapeutic strategies.

Phosphorylation-independent ubiquitylation and endocytosis of Fc gammaRIIA

Endocytosis of the Fc receptor Fc gammaRIIA depends on a functional ubiquitin conjugation system, and the receptor becomes ubiquitylated upon ligand binding. Phosphorylation of tyrosines in Fc gammaRIIA by Src family kinases is thought to be the initiating event in its signaling. However, although the Src family kinase inhibitor PP1 inhibited both ligand-induced phosphorylation of Fc gammaRIIA and phagocytosis in ts20 cells expressing Fc gammaRIIA, it did not inhibit receptor ubiquitylation or endocytosis of soluble ligands. Conversely, genistein and the proteasomal inhibitor MG132 did not inhibit receptor phosphorylation but strongly inhibited both receptor ubiquitylation and endocytosis. A region of the receptor lying within the immunoreceptor tyrosine-based activation motif was found to be necessary for both ubiquitylation and endocytosis. Ubiquitylation occurs at the plasma membrane before internalization. Endocytosis of Fc gammaRIIA is dependent on clathrin but independent of the adaptor protein AP-2. These findings point to a novel mechanism for ubiquitylation and endocytosis of this immunoreceptor.