Genome-wide and targeted CRISPR screens identify RNF213 as a mediator of interferon gamma-dependent pathogen restriction in human cells

To define cellular immunity to the intracellular pathogen Toxoplasma gondii, we performed a genome-wide CRISPR loss-of-function screen to identify genes important for (interferon gamma) IFN-γ-dependent growth restriction. We revealed a role for the tumor suppressor NF2/Merlin for maximum induction of Interferon Stimulated Genes (ISG), which are positively regulated by the transcription factor IRF-1. We then performed an ISG-targeted CRISPR screen that identified the host E3 ubiquitin ligase RNF213 as necessary for IFN-γ-mediated control of T. gondii in multiple human cell types. RNF213 was also important for control of bacterial (Mycobacterium tuberculosis) and viral (Vesicular Stomatitis Virus) pathogens in human cells. RNF213-mediated ubiquitination of the parasitophorous vacuole membrane (PVM) led to growth restriction of T. gondii in response to IFN-γ. Moreover, overexpression of RNF213 in naive cells also impaired growth of T. gondii. Surprisingly, growth inhibition did not require the autophagy protein ATG5, indicating that RNF213 initiates restriction independent of a previously described noncanonical autophagy pathway. Mutational analysis revealed that the ATPase domain of RNF213 was required for its recruitment to the PVM, while loss of a critical histidine in the RZ finger domain resulted in partial reduction of recruitment to the PVM and complete loss of ubiquitination. Both RNF213 mutants lost the ability to restrict growth of T. gondii, indicating that both recruitment and ubiquitination are required. Collectively, our findings establish RNF213 as a critical component of cell-autonomous immunity that is both necessary and sufficient for control of intracellular pathogens in human cells.

A Genome-Wide CRISPR Interference Screen Reveals an StkP-Mediated Connection between Cell Wall Integrity and Competence in Streptococcus salivarius

Competence is one of the most efficient bacterial evolutionary and adaptative strategies by synchronizing production of antibacterial compounds and integration of DNA released by dead cells. In most streptococci, this tactic is orchestrated by the ComRS system, a pheromone communication device providing a short time window of activation in which only part of the population is responsive. Understanding how this developmental process integrates multiple inputs to fine-tune the adequate response is a long-standing question. However, essential genes involved in the regulation of ComRS have been challenging to study. In this work, we built a conditional mutant library using CRISPR interference and performed three complementary screens to investigate competence genetic regulation in the human commensal Streptococcus salivarius. We show that initiation of competence increases upon cell wall impairment, suggesting a connection between cell envelope stress and competence activation. Notably, we report a key role for StkP, a serine-threonine kinase known to regulate cell wall homeostasis. We show that StkP controls competence by a mechanism that reacts to peptidoglycan fragments. Together, our data suggest a key cell wall sensing mechanism coupling competence to cell envelope integrity. IMPORTANCE Survival of human commensal streptococci in the digestive tract requires efficient strategies which must be tightly and collectively controlled for responding to competitive pressure and drastic environmental changes. In this context, the autocrine signaling system ComRS controlling competence for natural transformation and predation in salivarius streptococci could be seen as a multi-input device integrating a variety of environmental stimuli. In this work, we revealed novel positive and negative competence modulators by using a genome-wide CRISPR interference strategy. Notably, we highlighted an unexpected connection between bacterial envelope integrity and competence activation that involves several cell wall sensors. Together, these results showcase how commensal streptococci can fine-tune the pheromone-based competence system by responding to multiple inputs affecting their physiological status in order to calibrate an appropriate collective behavior.

A Genome-Wide CRISPR Screen Identifies Factors Regulating Pluripotency Exit in Mouse Embryonic Stem Cells

Pluripotency maintenance and exit in embryonic stem cells is a focal topic in stem cell biology. However, the effects of screening under very stringent culture conditions (e.g., differentiation medium, no leukemia inhibitory factor, no chemical inhibitors such as PD0325901 and CHIR99021, and no feeder cells) and of prolonging culture for key factors that regulate pluripotency exit, have not yet been reported. Here, we used a genome-wide CRISPR library to perform such a screen in mouse embryonic stem cells. Naïve NANOG-GFP mESCs were first transfected with a mouse genome-wide CRISPR knockout library to obtain a mutant mESCs library, followed by screening for two months in a strict N2B27 differentiation medium. The clones that survived our stringent screening were analyzed to identify the inserted sgRNAs. In addition to identifying the enriched genes that were reported in previous studies (Socs3, Tsc1, Trp53, Nf2, Tcf7l1, Csnk1a1, and Dhx30), we found 17 unreported genes, among which Zfp771 and Olfr769 appeared to be involved in pluripotency exit. Furthermore, Zfp771 knockout ESCs showed a differentiation delay in embryonic chimera experiments, indicating Zfp771 played an important role in pluripotency exit. Our results show that stringent screening with the CRISPR library can reveal key regulators of pluripotency exit.

Retromer Complex and PI3K Complex II-Related Genes Mediate the Yeast (Saccharomyces cerevisiae) Sodium Metabisulfite Resistance Response

Sodium metabisulfite (Na₂S₂O₅) is widely used as a preservative in the food and wine industry. However, it causes varying degrees of cellular damage to organisms. In order to improve our knowledge regarding its cyto-toxicity, a genome-wide screen using the yeast single deletion collection was performed. Additionally, a total of 162 Na₂S₂O₅-sensitive strains and 16 Na₂S₂O₅-tolerant strains were identified. Among the 162 Na₂S₂O₅ tolerance-related genes, the retromer complex was the top enriched cellular component. Further analysis demonstrated that retromer complex deletion leads to increased sensitivity to Na₂S₂O₅, and that Na₂S₂O₅ can induce mislocalization of retromer complex proteins. Notably, phosphatidylinositol 3-monophosphate kinase (PI3K) complex II, which is important for retromer recruitment to the endosome, might be a potential regulator mediating retromer localization and the yeast Na₂S₂O₅ tolerance response. Na₂S₂O₅ can decrease the protein expressions of Vps34, which is the component of PI3K complex. Therefore, Na₂S₂O₅-mediated retromer redistribution might be caused by the effects of decreased Vps34 expression levels. Moreover, both pharmaceutical inhibition of Vps34 functions and deletions of PI3K complex II-related genes affect cell tolerance to Na₂S₂O₅. The results of our study provide a global picture of cellular components required for Na₂S₂O₅ tolerance and advance our understanding concerning Na₂S₂O₅-induced cytotoxicity effects.

A synthetic lethal screen identifies HDAC4 as a potential target in MELK overexpressing cancers

Maternal embryonic leucine zipper kinase (MELK) is frequently overexpressed in cancer, but the role of MELK in cancer is still poorly understood. MELK was shown to have roles in many cancer-associated processes including tumor growth, chemotherapy resistance, and tumor recurrence. To determine whether the frequent overexpression of MELK can be exploited in therapy, we performed a high-throughput screen using a library of Saccharomyces cerevisiae mutants to identify genes whose functions become essential when MELK is overexpressed. We identified two such genes: LAG2 and HDA3. LAG2 encodes an inhibitor of the Skp, Cullin, F-box containing (SCF) ubiquitin-ligase complex, while HDA3 encodes a subunit of the HDA1 histone deacetylase complex. We find that one of these synthetic lethal interactions is conserved in mammalian cells, as inhibition of a human homolog of HDA3 (Histone Deacetylase 4, HDAC4) is synthetically toxic in MELK overexpression cells. Altogether, our work identified a novel potential drug target for tumors that overexpress MELK.

A genome-wide microRNA screen identifies the microRNA-183/96/182 cluster as a modulator of circadian rhythms

The regulatory mechanisms of circadian rhythms have been studied primarily at the level of the transcription-translation feedback loops of protein-coding genes. Regulatory modules involving noncoding RNAs are less thoroughly understood. In particular, emerging evidence has revealed the important role of microRNAs (miRNAs) in maintaining the robustness of the circadian system. To identify miRNAs that have the potential to modulate circadian rhythms, we conducted a genome-wide miRNA screen using U2OS luciferase reporter cells. Among 989 miRNAs in the library, 120 changed the period length in a dose-dependent manner. We further validated the circadian regulatory function of an miRNA cluster, miR-183/96/182, both in vitro and in vivo. We found that all three members of this miRNA cluster can modulate circadian rhythms. Particularly, miR-96 directly targeted a core circadian clock gene, PER2. The knockout of the miR-183/96/182 cluster in mice showed tissue-specific effects on circadian parameters and altered circadian rhythms at the behavioral level. This study identified a large number of miRNAs, including the miR-183/96/182 cluster, as circadian modulators. We provide a resource for further understanding the role of miRNAs in the circadian network and highlight the importance of miRNAs as a genome-wide layer of circadian clock regulation.

A Genome-Wide Screen in Saccharomyces cerevisiae Reveals a Critical Role for Oxidative Phosphorylation in Cellular Tolerance to Lithium Hexafluorophosphate

Lithium hexafluorophosphate (LiPF₆) is one of the leading electrolytes in lithium-ion batteries, and its usage has increased tremendously in the past few years. Little is known, however, about its potential environmental and biological impacts. In order to improve our understanding of the cytotoxicity of LiPF₆ and the specific cellular response mechanisms to it, we performed a genome-wide screen using a yeast (Saccharomyces cerevisiae) deletion mutant collection and identified 75 gene deletion mutants that showed LiPF₆ sensitivity. Among these, genes associated with mitochondria showed the most enrichment. We also found that LiPF₆ is more toxic to yeast than lithium chloride (LiCl) or sodium hexafluorophosphate (NaPF₆). Physiological analysis showed that a high concentration of LiPF₆ caused mitochondrial damage, reactive oxygen species (ROS) accumulation, and ATP content changes. Compared with the results of previous genome-wide screening for LiCl-sensitive mutants, we found that oxidative phosphorylation-related mutants were specifically hypersensitive to LiPF₆. In these deletion mutants, LiPF₆ treatment resulted in higher ROS production and reduced ATP levels, suggesting that oxidative phosphorylation-related genes were important for counteracting LiPF₆-induced toxicity. Taken together, our results identified genes specifically involved in LiPF₆-modulated toxicity, and demonstrated that oxidative stress and ATP imbalance maybe the driving factors in governing LiPF₆-induced toxicity.

A novel p53 regulator, C16ORF72/TAPR1, buffers against telomerase inhibition

Telomere erosion in cells with insufficient levels of the telomerase reverse transcriptase (TERT), contributes to age-associated tissue dysfunction and senescence, and p53 plays a crucial role in this response. We undertook a genome-wide CRISPR screen to identify gene deletions that sensitized p53-positive human cells to telomerase inhibition. We uncovered a previously unannotated gene, C16ORF72, which we term Telomere Attrition and p53 Response 1 (TAPR1), that exhibited a synthetic-sick relationship with TERT loss. A subsequent genome-wide CRISPR screen in TAPR1-disrupted cells reciprocally identified TERT as a sensitizing gene deletion. Cells lacking TAPR1 or TERT possessed elevated p53 levels and transcriptional signatures consistent with p53 upregulation. The elevated p53 response in TERT- or TAPR1-deficient cells was exacerbated by treatment with the MDM2 inhibitor and p53 stabilizer nutlin-3a and coincided with a further reduction in cell fitness. Importantly, the sensitivity to treatment with nutlin-3a in TERT- or TAPR1-deficient cells was rescued by loss of p53. These data suggest that TAPR1 buffers against the deleterious consequences of telomere erosion or DNA damage by constraining p53. These findings identify C16ORF72/TAPR1 as new regulator at the nexus of telomere integrity and p53 regulation.

Systematic Identification of Regulators of Oxidative Stress Reveals Non-canonical Roles for Peroxisomal Import and the Pentose Phosphate Pathway

Reactive oxygen species (ROS) play critical roles in metabolism and disease, yet a comprehensive analysis of the cellular response to oxidative stress is lacking. To systematically identify regulators of oxidative stress, we conducted genome-wide Cas9/CRISPR and shRNA screens. This revealed a detailed picture of diverse pathways that control oxidative stress response, ranging from the TCA cycle and DNA repair machineries to iron transport, trafficking, and metabolism. Paradoxically, disrupting the pentose phosphate pathway (PPP) at the level of phosphogluconate dehydrogenase (PGD) protects cells against ROS. This dramatically alters metabolites in the PPP, consistent with rewiring of upper glycolysis to promote antioxidant production. In addition, disruption of peroxisomal import unexpectedly increases resistance to oxidative stress by altering the localization of catalase. Together, these studies provide insights into the roles of peroxisomal matrix import and the PPP in redox biology and represent a rich resource for understanding the cellular response to oxidative stress.

Identification of Required Host Factors for SARS-CoV-2 Infection in Human Cells

To better understand host-virus genetic dependencies and find potential therapeutic targets for COVID-19, we performed a genome-scale CRISPR loss-of-function screen to identify host factors required for SARS-CoV-2 viral infection of human alveolar epithelial cells. Top-ranked genes cluster into distinct pathways, including the vacuolar ATPase proton pump, Retromer, and Commander complexes. We validate these gene targets using several orthogonal methods such as CRISPR knockout, RNA interference knockdown, and small-molecule inhibitors. Using single-cell RNA-sequencing, we identify shared transcriptional changes in cholesterol biosynthesis upon loss of top-ranked genes. In addition, given the key role of the ACE2 receptor in the early stages of viral entry, we show that loss of RAB7A reduces viral entry by sequestering the ACE2 receptor inside cells. Overall, this work provides a genome-scale, quantitative resource of the impact of the loss of each host gene on fitness/response to viral infection.

A Genome-wide CRISPR Screen Identifies ZCCHC14 as a Host Factor Required for Hepatitis B Surface Antigen Production

A hallmark of chronic hepatitis B (CHB) virus infection is the presence of high circulating levels of non-infectious small lipid HBV surface antigen (HBsAg) vesicles. Although rare, sustained HBsAg loss is the idealized endpoint of any CHB therapy. A small molecule, RG7834, has been previously reported to inhibit HBsAg expression by targeting terminal nucleotidyltransferase proteins 4A and 4B (TENT4A and TENT4B). In this study, we describe a genome-wide CRISPR screen to identify other potential host factors required for HBsAg expression and to gain further insights into the mechanism of RG7834. We report more than 60 genes involved in regulating HBsAg and identify additional factors involved in RG7834 activity, including a zinc finger CCHC-type containing 14 (ZCCHC14) protein. We show that ZCCHC14, together with TENT4A/B, stabilizes HBsAg expression through HBV RNA tailing, providing a potential new therapeutic target to achieve functional cure in CHB patients.

Genome-Wide Screens Reveal that Resveratrol Induces Replicative Stress in Human Cells

Resveratrol is a natural product associated with wide-ranging effects in animal and cellular models, including lifespan extension. To identify the genetic target of resveratrol in human cells, we conducted genome-wide CRISPR-Cas9 screens to pinpoint genes that confer sensitivity or resistance to resveratrol. An extensive network of DNA damage response and replicative stress genes exhibited genetic interactions with resveratrol and its analog pterostilbene. These genetic profiles showed similarity to the response to hydroxyurea, an inhibitor of ribonucleotide reductase that causes replicative stress. Resveratrol, pterostilbene, and hydroxyurea caused similar depletion of nucleotide pools, inhibition of replication fork progression, and induction of replicative stress. The ability of resveratrol to inhibit cell proliferation and S phase transit was independent of the histone deacetylase sirtuin 1, which has been implicated in lifespan extension by resveratrol. These results establish that a primary impact of resveratrol on human cell proliferation is the induction of low-level replicative stress.

Genome-Wide Screen of the Hippocampus in Aged Rats Identifies Mitochondria, Metabolism and Aging Processes Implicated in Sevoflurane Anesthesia

Previous studies have shown multiple mechanisms and pathophysiological changes after anesthesia, and genome-wide studies have been implemented in the studies of brain aging and neurodegenerative diseases. However, the genome-wide gene expression patterns and modulation networks after general anesthesia remains to be elucidated. Therefore, whole transcriptome microarray analysis was used to explore the coding gene expression patterns in the hippocampus of aged rats after sevoflurane anesthesia. Six hundred and thirty one upregulated and 183 downregulated genes were screened out, then 44 enriched terms of biological process, 16 of molecular function and 18 of the cellular components were identified by Gene Ontology (GO) and KEGG analysis. Among them, oxidative stress, metabolism, aging, and neurodegeneration were the most enriched biological processes and changed functions. Thus, involved genes of these processes were selected for qPCR verification and a good consistency was confirmed. The potential signaling pathways were further constructed including mitochondrion and oxidative stress-related Hifs-Prkcd-Akt-Nfe2l2-Sod1 signaling, multiple metabolism signaling (Scd2, Scap-Hmgcs2, Aldh18a1-Glul and Igf1r), as well as aging and neurodegeneration related signaling (Spidr-Ercc4-Cdkn1a-Pmaip1 and Map1lc3b). These results provide potential therapeutic gene targets for brain function modulation and memory formation process after inhaled anesthesia in the elderly, which could be valuable for preventing postoperative brain disorders and diseases, such as perioperative neurocognitive disorders (PND), from the genetic level in the future.

A Genome-wide ER-phagy Screen Highlights Key Roles of Mitochondrial Metabolism and ER-Resident UFMylation

Selective autophagy of organelles is critical for cellular differentiation, homeostasis, and organismal health. Autophagy of the ER (ER-phagy) is implicated in human neuropathy but is poorly understood beyond a few autophagosomal receptors and remodelers. By using an ER-phagy reporter and genome-wide CRISPRi screening, we identified 200 high-confidence human ER-phagy factors. Two pathways were unexpectedly required for ER-phagy. First, reduced mitochondrial metabolism represses ER-phagy, which is opposite of general autophagy and is independent of AMPK. Second, ER-localized UFMylation is required for ER-phagy to repress the unfolded protein response via IRE1α. The UFL1 ligase is brought to the ER surface by DDRGK1 to UFMylate RPN1 and RPL26 and preferentially targets ER sheets for degradation, analogous to PINK1-Parkin regulation during mitophagy. Our data provide insight into the cellular logic of ER-phagy, reveal parallels between organelle autophagies, and provide an entry point to the relatively unexplored process of degrading the ER network.

Doxycycline-Dependent Self-Inactivation of CRISPR-Cas9 to Temporally Regulate On- and Off-Target Editing

Exome and deep sequencing of cells treated with a panel of lentiviral guide RNA demonstrate that both on- and off-target editing proceed in a time-dependent manner. Thus, methods to temporally control Cas9 activity would be beneficial. To address this need, we describe a "self-inactivating CRISPR (SiC)" system consisting of a single guide RNA that deactivates the Streptococcus pyogenes Cas9 nuclease in a doxycycline-dependent manner. This enables defined, temporal control of Cas9 activity in any cell type and also in vivo. Results show that SiC may enable a reduction in off-target editing, with less effect on on-target editing rates. This tool facilitates diverse applications including (1) the timed regulation of genetic knockouts in hard-to-transfect cells using lentivirus, including human leukocytes for the identification of glycogenes regulating leukocyte-endothelial cell adhesion; (2) genome-wide lentiviral sgRNA (single guide RNA) library applications where Cas9 activity is ablated after allowing pre-determined editing times. Thus, stable knockout cell pools are created for functional screens; and (3) temporal control of Cas9-mediated editing of myeloid and lymphoid cells in vivo, both in mouse peripheral blood and bone marrow. Overall, SiC enables temporal control of gene editing and may be applied in diverse application including studies that aim to reduce off-target genome editing.

Mitochondrial Reprogramming Underlies Resistance to BCL-2 Inhibition in Lymphoid Malignancies

Mitochondrial apoptosis can be effectively targeted in lymphoid malignancies with the FDA-approved B cell lymphoma 2 (BCL-2) inhibitor venetoclax, but resistance to this agent is emerging. We show that venetoclax resistance in chronic lymphocytic leukemia is associated with complex clonal shifts. To identify determinants of resistance, we conducted parallel genome-scale screens of the BCL-2-driven OCI-Ly1 lymphoma cell line after venetoclax exposure along with integrated expression profiling and functional characterization of drug-resistant and engineered cell lines. We identified regulators of lymphoid transcription and cellular energy metabolism as drivers of venetoclax resistance in addition to the known involvement by BCL-2 family members, which were confirmed in patient samples. Our data support the implementation of combinatorial therapy with metabolic modulators to address venetoclax resistance.

A Genome-Wide Knockout Screen in Human Macrophages Identified Host Factors Modulating Salmonella Infection

A genome-scale CRISPR knockout library screen of THP-1 human macrophages was performed to identify loss-of-function mutations conferring resistance to Salmonella uptake. The screen identified 183 candidate genes, from which 14 representative genes involved in actin dynamics (ACTR3, ARPC4, CAPZB, TOR3A, CYFIP2, CTTN, and NHLRC2), glycosaminoglycan metabolism (B3GNT1), receptor signaling (PDGFB and CD27), lipid raft formation (CLTCL1), calcium transport (ATP2A2 and ITPR3), and cholesterol metabolism (HMGCR) were analyzed further. For some of these pathways, known chemical inhibitors could replicate the Salmonella resistance phenotype, indicating their potential as targets for host-directed therapy. The screen indicated a role for the relatively uncharacterized gene NHLRC2 in both Salmonella invasion and macrophage differentiation. Upon differentiation, NHLRC2 mutant macrophages were hyperinflammatory and did not exhibit characteristics typical of macrophages, including atypical morphology and inability to interact and phagocytose bacteria/particles. Immunoprecipitation confirmed an interaction of NHLRC2 with FRYL, EIF2AK2, and KLHL13.IMPORTANCE Salmonella exploits macrophages to gain access to the lymphatic system and bloodstream to lead to local and potentially systemic infections. With an increasing number of antibiotic-resistant isolates identified in humans, Salmonella infections have become major threats to public health. Therefore, there is an urgent need to identify alternative approaches to anti-infective therapy, including host-directed therapies. In this study, we used a simple genome-wide screen to identify 183 candidate host factors in macrophages that can confer resistance to Salmonella infection. These factors may be potential therapeutic targets against Salmonella infections.

Efficient genome-wide first-generation phenotypic screening system in mice using the piggyBac transposon

Genome-wide phenotypic screens provide an unbiased way to identify genes involved in particular biological traits, and have been widely used in lower model organisms. However, cost and time have limited the utility of such screens to address biological and disease questions in mammals. Here we report a highly efficient piggyBac (PB) transposon-based first-generation (F1) dominant screening system in mice that enables an individual investigator to conduct a genome-wide phenotypic screen within a year with fewer than 300 cages. The PB screening system uses visually trackable transposons to induce both gain- and loss-of-function mutations and generates genome-wide distributed new insertions in more than 55% of F1 progeny. Using this system, we successfully conducted a pilot F1 screen and identified 5 growth retardation mutations. One of these mutants, a Six1/4 ^PB/+ mutant, revealed a role in milk intake behavior. The mutant animals exhibit abnormalities in nipple recognition and milk ingestion, as well as developmental defects in cranial nerves V, IX, and X. This PB F1 screening system offers individual laboratories unprecedented opportunities to conduct affordable genome-wide phenotypic screens for deciphering the genetic basis of mammalian biology and disease pathogenesis.

A CRISPR Activation Screen Identifies Genes That Protect against Zika Virus Infection

Zika virus (ZIKV) is an arthropod-borne emerging pathogen causing febrile illness. ZIKV is associated Guillain-Barré syndrome and other neurological complications. Infection during pregnancy is associated with pregnancy complications and developmental and neurological abnormalities collectively defined as congenital Zika syndrome. There is still no vaccine or specific treatment for ZIKV infection. To identify host factors that can rescue cells from ZIKV infection, we used a genome-scale CRISPR activation screen. Our highly ranking hits included a short list of interferon-stimulated genes (ISGs) previously reported to have antiviral activity. Validation of the screen results highlighted interferon lambda 2 (IFN-λ2) and interferon alpha-inducible protein 6 (IFI6) as genes providing high levels of protection from ZIKV. Activation of these genes had an effect on an early stage in viral infection. In addition, infected cells expressing single guide RNAs (sgRNAs) for both of these genes displayed lower levels of cell death than did the controls. Furthermore, the identified genes were significantly induced in ZIKV-infected placenta explants. Thus, these results highlight a set of ISGs directly relevant for rescuing cells from ZIKV infection or its associated cell death and substantiate CRISPR activation screens as a tool to identify host factors impeding pathogen infection.IMPORTANCE Zika virus (ZIKV) is an emerging vector-borne pathogen causing a febrile disease. ZIKV infection might also trigger Guillain-Barré syndrome, neuropathy, and myelitis. Vertical transmission of ZIKV can cause fetus demise, stillbirth, or severe congenital abnormalities and neurological complications. There is no vaccine or specific antiviral treatment against ZIKV. We used a genome-wide CRISPR activation screen, where genes are activated from their native promoters to identify host cell factors that protect cells from ZIKV infection or associated cell death. The results provide a better understanding of key host factors that protect cells from ZIKV infection and might assist in identifying novel antiviral targets.

A Genome-Wide Haploid Genetic Screen Identifies Heparan Sulfate-Associated Genes and the Macropinocytosis Modulator TMED10 as Factors Supporting Vaccinia Virus Infection

Vaccinia virus is a promising viral vaccine and gene delivery candidate and has historically been used as a model to study poxvirus-host cell interactions. We employed a genome-wide insertional mutagenesis approach in human haploid cells to identify host factors crucial for vaccinia virus infection. A library of mutagenized HAP1 cells was exposed to modified vaccinia virus Ankara (MVA). Deep-sequencing analysis of virus-resistant cells identified host factors involved in heparan sulfate synthesis, Golgi organization, and vesicular protein trafficking. We validated EXT1, TM9SF2, and TMED10 (TMP21/p23/p24δ) as important host factors for vaccinia virus infection. The critical roles of EXT1 in heparan sulfate synthesis and vaccinia virus infection were confirmed. TM9SF2 was validated as a player mediating heparan sulfate expression, explaining its contribution to vaccinia virus infection. In addition, TMED10 was found to be crucial for virus-induced plasma membrane blebbing and phosphatidylserine-induced macropinocytosis, presumably by regulating the cell surface expression of the TAM receptor Axl.IMPORTANCE Poxviruses are large DNA viruses that can infect a wide range of host species. A number of these viruses are clinically important to humans, including variola virus (smallpox) and vaccinia virus. Since the eradication of smallpox, zoonotic infections with monkeypox virus and cowpox virus are emerging. Additionally, poxviruses can be engineered to specifically target cancer cells and are used as a vaccine vector against tuberculosis, influenza, and coronaviruses. Poxviruses rely on host factors for most stages of their life cycle, including attachment to the cell and entry. These host factors are crucial for virus infectivity and host cell tropism. We used a genome-wide knockout library of host cells to identify host factors necessary for vaccinia virus infection. We confirm a dominant role for heparin sulfate in mediating virus attachment. Additionally, we show that TMED10, previously not implicated in virus infections, facilitates virus uptake by modulating the cellular response to phosphatidylserine.

Genome-Scale CRISPRa Screen Identifies Novel Factors for Cellular Reprogramming

Primed epiblast stem cells (EpiSCs) can be reverted to a pluripotent embryonic stem cell (ESC)-like state by expression of single reprogramming factor. We used CRISPR activation to perform a genome-scale, reprogramming screen in EpiSCs and identified 142 candidate genes. Our screen validated a total of 50 genes, previously not known to contribute to reprogramming, of which we chose Sall1 for further investigation. We show that Sall1 augments reprogramming of mouse EpiSCs and embryonic fibroblasts and that these induced pluripotent stem cells are indeed fully pluripotent including formation of chimeric mice. We also demonstrate that Sall1 synergizes with Nanog in reprogramming and that overexpression in ESCs delays their conversion back to EpiSCs. Lastly, using RNA sequencing, we identify and validate Klf5 and Fam189a2 as new downstream targets of Sall1 and Nanog. In summary, our work demonstrates the power of using CRISPR technology in understanding molecular mechanisms that mediate complex cellular processes such as reprogramming.

•

Genome-scale CRISPRa screen in mouse EpiSCs identifies novel reprogramming factors

•

50 novel genes, including Sall1 and Fam189a2, identified to mediate reprogramming

•

Sall1 synergizes with Nanog to increase reprogramming efficiency in EpiSCs and MEFs

•

RNA-seq provides insight into downstream pathways of Sall1 and Nanog-mediated reprogramming

CRISPR-Cas9-based genome-wide screening of Toxoplasma gondii

Apicomplexan parasites, such as Toxoplasma gondii, cause extensive morbidity and mortality in humans and livestock, highlighting the need for a deeper understanding of their molecular biology. Although techniques for the generation of targeted gene disruptions have long been available for apicomplexans, such methods are not readily scalable to the entire genome. We recently used CRISPR-Cas9 to disrupt all nuclear protein-coding genes in T. gondii using a pooled format. The method relies on transfection of a guide RNA library into parasites constitutively expressing Cas9. Here, we present the complete workflow of such a screen, including preparation of the guide RNA library, growth and testing of the recipient strain, generation of the mutant population, culture conditions for the screen, preparation of genomic DNA libraries, next-generation sequencing of the guide RNA loci, and analysis to detect fitness-conferring genes. This method can be deployed to study how culture conditions affect the repertoire of genes needed by parasites, which will enable studies of their metabolic needs, host specificity, and drug-resistance mechanisms. In addition, by manipulating the background in which the screen is performed, researchers will be able to investigate genetic interactions, which may help uncover redundancy or epistasis in the parasite genome. Using this method, a genome-wide screen and its analysis can be completed in 3 weeks, after ∼1 month of preparation to generate the library and grow the cells needed, making it a powerful tool for uncovering functionally important genes in apicomplexan parasites.

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

In mitophagy, damaged mitochondria are targeted for disposal by the autophagy machinery. PARKIN promotes signaling of mitochondrial damage to the autophagy machinery for engagement, and PARKIN mutations cause Parkinson’s disease, possibly because damaged mitochondria accumulate in neurons. Because regulation of PARKIN abundance and the impact on signaling are poorly understood, we performed a genetic screen to identify PARKIN abundance regulators. Both positive and negative regulators were identified and will help us to further understand mitophagy and Parkinson’s disease. We show that some of the identified genes negatively regulate PARKIN gene expression, which impacts signaling of mitochondrial damage in mitophagy. This link between transcriptional repression and mitophagy is also apparent in neurons in culture, bearing implications for disease.

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.

Explicit Modeling of siRNA-Dependent On- and Off-Target Repression Improves the Interpretation of Screening Results

RNAi is broadly used to map gene regulatory networks, but the identification of genes that are responsible for the observed phenotypes is challenging, as small interfering RNAs (siRNAs) simultaneously downregulate the intended on targets and many partially complementary off targets. Additionally, the scarcity of publicly available control datasets hinders the development and comparative evaluation of computational methods for analyzing the data. Here, we introduce PheLiM (https://github.com/andreariba/PheLiM), a method that uses predictions of siRNA on- and off-target downregulation to infer gene-specific contributions to phenotypes. To assess the performance of PheLiM, we carried out siRNA- and CRISPR/Cas9-based genome-wide screening of two well-characterized pathways, bone morphogenetic protein (BMP) and nuclear factor κB (NF-κB), and we reanalyzed publicly available siRNA screens. We demonstrate that PheLiM has the overall highest accuracy and most reproducible results compared to other available methods. PheLiM can accommodate various methods for predicting siRNA off targets and is broadly applicable to the identification of genes underlying complex phenotypes.

The New Role for an Old Kinase: Protein Kinase CK2 Regulates Metal Ion Transport

The pleiotropic serine/threonine protein kinase CK2 was the first kinase discovered. It is renowned for its role in cell proliferation and anti-apoptosis. The complexity of this kinase is well reflected by the findings of past decades in terms of its heterotetrameric structure, subcellular location, constitutive activity and the extensive catalogue of substrates. With the advent of non-biased high-throughput functional genomics such as genome-wide deletion mutant screening, novel aspects of CK2 functionality have been revealed. Our recent discoveries using the model organism Saccharomyces cerevisiae and mammalian cells demonstrate that CK2 regulates metal toxicity. Extensive literature search reveals that there are few but elegant works on the role of CK2 in regulating the sodium and zinc channels. As both CK2 and metal ions are key players in cell biology and oncogenesis, understanding the details of CK2’s regulation of metal ion homeostasis has a direct bearing on cancer research. In this review, we aim to garner the recent data and gain insights into the role of CK2 in metal ion transport.

Genome-wide screen uncovers novel pathways for tRNA processing and nuclear-cytoplasmic dynamics

In this resource, Wu et al. present the first comprehensive unbiased analysis of the role of nearly an entire proteome in tRNA biology and describe 162 novel and 12 previously known Saccharomyces cerevisiae gene products that function in tRNA processing, turnover, and subcellular movement. The findings from this genome-wide screen describe putative novel pathways for tRNA nuclear export and extensive links between tRNA biology and other aspects of cell physiology.

Transfer ribonucleic acids (tRNAs) are essential for protein synthesis. However, key gene products involved in tRNA biogenesis and subcellular movement remain to be discovered. We conducted the first comprehensive unbiased analysis of the role of nearly an entire proteome in tRNA biology and describe 162 novel and 12 previously known Saccharomyces cerevisiae gene products that function in tRNA processing, turnover, and subcellular movement. tRNA nuclear export is of particular interest because it is essential, but the known tRNA exporters (Los1 [exportin-t] and Msn5 [exportin-5]) are unessential. We report that mutations of CRM1 (Exportin-1), MEX67/MTR2 (TAP/p15), and five nucleoporins cause accumulation of unspliced tRNA, a hallmark of defective tRNA nuclear export. CRM1 mutation genetically interacts with los1Δ and causes altered tRNA nuclear–cytoplasmic distribution. The data implicate roles for the protein and mRNA nuclear export machineries in tRNA nuclear export. Mutations of genes encoding actin cytoskeleton components and mitochondrial outer membrane proteins also cause accumulation of unspliced tRNA, likely due to defective splicing on mitochondria. Additional gene products, such as chromatin modification enzymes, have unanticipated effects on pre-tRNA end processing. Thus, this genome-wide screen uncovered putative novel pathways for tRNA nuclear export and extensive links between tRNA biology and other aspects of cell physiology.

Spatial sequestration and detoxification of Huntingtin by the ribosome quality control complex

Huntington disease (HD) is a neurological disorder caused by polyglutamine expansions in mutated Huntingtin (mHtt) proteins, rendering them prone to form inclusion bodies (IB). We report that in yeast, such IB formation is a factor-dependent process subjected to age-related decline. A genome-wide, high-content imaging approach, identified the E3 ubiquitin ligase, Ltn1 of the ribosome quality control complex (RQC) as a key factor required for IB formation, ubiquitination, and detoxification of model mHtt. The failure of ltn1∆ cells to manage mHtt was traced to another RQC component, Tae2, and inappropriate control of heat shock transcription factor, Hsf1, activity. Moreover, super-resolution microscopy revealed that mHtt toxicity in RQC-deficient cells was accompanied by multiple mHtt aggregates altering actin cytoskeletal structures and retarding endocytosis. The data demonstrates that spatial sequestration of mHtt into IBs is policed by the RQC-Hsf1 regulatory system and that such compartmentalization, rather than ubiquitination, is key to mHtt detoxification.

DOI:http://dx.doi.org/10.7554/eLife.11792.001

Huntington’s disease is a neurological disease that is caused by mutations in the gene that encodes a protein called Htt. Individuals with this mutation gradually lose neurons as they age, resulting in declines in muscle coordination and mental abilities. The mutant Htt proteins tend to form clumps inside cells, but it is not clear if these clumps are the cause of the disease symptoms or whether they have a protective effect.

Yang et al. used yeast as a model to investigate whether the mutant Htt proteins need other molecules to allow them to form clumps. The experiments identified several new molecules that are required for mutated Htt to form clumps. Some of these are components of a system called the Ribosome Quality Control (RQC) complex, which monitors newly made proteins and labels abnormal ones for destruction. However, Yang et al.’s findings suggest that the RQC complex regulates the formation of Htt clumps through a different pathway involving a protein called heat shock factor 1. In this case, cells would need to fine-tune heat shock factor 1 activity to make mutant Htt proteins clump together to protect cells from damage.

Future experiments should expand Yang et al.’s findings to animal models of Huntington’s disease and identify which other molecules contribute to the formation of Htt clumps. One challenge will be to find out why older neurons fail to form clumps of Htt proteins, and whether this can be overcome by drugs that boost the activity of the molecules that Yang et al. identified.

DOI:http://dx.doi.org/10.7554/eLife.11792.002

Genome-wide CRISPR screen in a mouse model of tumor growth and metastasis

Genetic screens are powerful tools for identifying genes responsible for diverse phenotypes. Here we describe a genome-wide CRISPR/Cas9-mediated loss-of-function screen in tumor growth and metastasis. We mutagenized a non-metastatic mouse cancer cell line using a genome-scale library with 67,405 single-guide RNAs (sgRNAs). The mutant cell pool rapidly generates metastases when transplanted into immunocompromised mice. Enriched sgRNAs in lung metastases and late-stage primary tumors were found to target a small set of genes, suggesting that specific loss-of-function mutations drive tumor growth and metastasis. Individual sgRNAs and a small pool of 624 sgRNAs targeting the top-scoring genes from the primary screen dramatically accelerate metastasis. In all of these experiments, the effect of mutations on primary tumor growth positively correlates with the development of metastases. Our study demonstrates Cas9-based screening as a robust method to systematically assay gene phenotypes in cancer evolution in vivo.

Genetic basis for Saccharomyces cerevisiae biofilm in liquid medium

Biofilm-forming microorganisms switch between two forms: free-living planktonic and sessile multicellular. Sessile communities of yeast biofilms in liquid medium provide a primitive example of multicellularity and are clinically important because biofilms tend to have other growth characteristics than free-living cells. We investigated the genetic basis for yeast, Saccharomyces cerevisiae, biofilm on solid surfaces in liquid medium by screening a comprehensive deletion mutant collection in the Σ1278b background and found 71 genes that were essential for biofilm development. Quantitative northern blots further revealed that AIM1, ASG1, AVT1, DRN1, ELP4, FLO8, FMP10, HMT1, KAR5, MIT1, MRPL32, MSS11, NCP1, NPR1, PEP5, PEX25, RIM8, RIM101, RGT1, SNF8, SPC2, STB6, STP22, TEC1, VID24, VPS20, VTC3, YBL029W, YBL029C-A, YFL054C, YGR161W-C, YIL014C-A, YIR024C, YKL151C, YNL200C, YOR034C-A, and YOR223W controlled biofilm through FLO11 induction. Almost all deletion mutants that were unable to form biofilms in liquid medium also lost the ability to form surface-spreading biofilm colonies (mats) on agar and 69% also lost the ability to grow invasively. The protein kinase A isoform Tpk3p functioned specifically in biofilm and mat formation. In a tpk3 mutant, transcription of FLO11 was induced three-fold compared with wild-type, but biofilm development and cell-cell adhesion was absent, suggesting that Tpk3p regulates FLO11 positive posttranscriptionally and negative transcriptionally.The study provides a resource of biofilm-influencing genes for additional research on biofilm development and suggests that the regulation of FLO11 is more complex than previously anticipated.

Enhancing the power to detect low-frequency variants in genome-wide screens

In genetic association studies a conventional test statistic is proportional to the correlation coefficient between the trait and the variant, with the result that it lacks power to detect association for low-frequency variants. Considering the link between the conventional association test statistics and the linkage disequilibrium measure r(2), we propose a test statistic analogous to the standardized linkage disequilibrium D' to increase the power of detecting association for low-frequency variants. By both simulation and real data analysis we show that the proposed D' test is more powerful than the conventional methods for detecting association for low-frequency variants in a genome-wide setting. The optimal coding strategy for the D' test and its asymptotic properties are also investigated. In summary, we advocate using the D' test in a dominant model as a complementary approach to enhancing the power of detecting association for low-frequency variants with moderate to large effect sizes in case-control genome-wide association studies.

Saccharomyces cerevisiae genes involved in survival of heat shock

The heat-shock response in cells, involving increased transcription of a specific set of genes in response to a sudden increase in temperature, is a highly conserved biological response occurring in all organisms. Despite considerable attention to the processes activated during heat shock, less is known about the role of genes in survival of a sudden temperature increase. Saccharomyces cerevisiae genes involved in the maintenance of heat-shock resistance in exponential and stationary phase were identified by screening the homozygous diploid deletants in nonessential genes and the heterozygous diploid mutants in essential genes for survival after a sudden shift in temperature from 30 to 50°. More than a thousand genes were identified that led to altered sensitivity to heat shock, with little overlap between them and those previously identified to affect thermotolerance. There was also little overlap with genes that are activated or repressed during heat-shock, with only 5% of them regulated by the heat-shock transcription factor. The target of rapamycin and protein kinase A pathways, lipid metabolism, vacuolar H⁺-ATPase, vacuolar protein sorting, and mitochondrial genome maintenance/translation were critical to maintenance of resistance. Mutants affected in l-tryptophan metabolism were heat-shock resistant in both growth phases; those affected in cytoplasmic ribosome biogenesis and DNA double-strand break repair were resistant in stationary phase, and in mRNA catabolic processes in exponential phase. Mutations affecting mitochondrial genome maintenance were highly represented in sensitive mutants. The cell division transcription factor Swi6p and Hac1p involved in the unfolded protein response also play roles in maintenance of heat-shock resistance.

Genome-wide screen of human bromodomain-containing proteins identifies Cecr2 as a novel DNA damage response protein

The formation of γ-H2AX foci after DNA double strand breaks (DSBs) is crucial for the cellular response to this lethal DNA damage. We previously have shown that BRG1, a chromatin remodeling enzyme, facilitates DSB repair by stimulating γ-H2AX formation, and this function of BRG1 requires the binding of BRGI to acetylated histone H3 on γ-H2AX-containing nucleosomes using its bromodomain (BRD), a protein module that specifically recognizes acetyl-Lys moieties. We also have shown that the BRD of BRG1, when ectopically expressed in cells, functions as a dominant negative inhibitor of the BRG1 activity to stimulate γ-H2AX and DSB repair. Here, we found that BRDs from a select group of proteins have no such activity, suggesting that the γ-H2AX inhibition activity of BRG1 BRD is specific. This finding led us to search for more BRDs that exhibit γ-H2AX inhibition activity in the hope of finding additional BRD-containing proteins involved in DNA damage responses. We screened a total of 52 individual BRDs present in 38 human BRD-containing proteins, comprising 93% of all human BRDs. We identified the BRD of cat eye syndrome chromosome region candidate 2 (Cecr2), which recently was shown to form a novel chromatin remodeling complex with unknown cellular functions, as having a strong γ-H2AX inhibition activity. This activity of Cecr2 BRD is specific because it depends on the chromatin binding affinity of Cecr2 BRD. Small interfering RNA knockdown experiments showed that Cecr2 is important for γ-H2AX formation and DSB repair. Therefore, our genomewide screen identifies Cecr2 as a novel DNA damage response protein.