Mapping of DDX11 genetic interactions defines sister chromatid cohesion as the major dependency.

DDX11/Chl1R is a conserved DNA helicase with roles in genome maintenance, DNA replication, and chromatid cohesion. Loss of DDX11 in humans leads to the rare cohesinopathy Warsaw breakage syndrome. DDX11 has also been implicated in human cancer where it has been proposed to have an oncogenic role and possibly to constitute a therapeutic target. Given the multiple roles of DDX11 in genome stability and its potential as an anticancer target, we set out to define a complete genetic interaction profile of DDX11 loss in human cell lines. Screening the human genome with clustered regularly interspaced short palindromic repeats (CRISPR) guide RNA drop out screens in DDX11-wildtype (WT) or DDX11-deficient cells revealed a strong enrichment of genes with functions related to sister chromatid cohesion. We confirm synthetic lethal relationships between DDX11 and the tumor suppressor cohesin subunit STAG2, which is frequently mutated in several cancer types and the kinase HASPIN. This screen highlights the importance of cohesion in cells lacking DDX11 and suggests DDX11 may be a therapeutic target for tumors with mutations in STAG2.

FAM72A antagonizes UNG2 to promote mutagenic repair during antibody maturation.

Activation-induced cytidine deaminase (AID) catalyses the deamination of deoxycytidines to deoxyuracils within immunoglobulin genes to induce somatic hypermutation and class-switch recombination1,2. AID-generated deoxyuracils are recognized and processed by subverted base-excision and mismatch repair pathways that ensure a mutagenic outcome in B cells3-6. However, why these DNA repair pathways do not accurately repair AID-induced lesions remains unknown. Here, using a genome-wide CRISPR screen, we show that FAM72A is a major determinant for the error-prone processing of deoxyuracils. Fam72a-deficient CH12F3-2 B cells and primary B cells from Fam72a-/- mice exhibit reduced class-switch recombination and somatic hypermutation frequencies at immunoglobulin and Bcl6 genes, and reduced genome-wide deoxyuracils. The somatic hypermutation spectrum in B cells from Fam72a-/- mice is opposite to that observed in mice deficient in uracil DNA glycosylase 2 (UNG2)7, which suggests that UNG2 is hyperactive in FAM72A-deficient cells. Indeed, FAM72A binds to UNG2, resulting in reduced levels of UNG2 protein in the G1 phase of the cell cycle, coinciding with peak AID activity. FAM72A therefore causes U·G mispairs to persist into S phase, leading to error-prone processing by mismatch repair. By disabling the DNA repair pathways that normally efficiently remove deoxyuracils from DNA, FAM72A enables AID to exert its full effects on antibody maturation. This work has implications in cancer, as the overexpression of FAM72A that is observed in many cancers8 could promote mutagenesis.

Paralogous synthetic lethality underlies genetic dependencies of the cancer-mutated gene STAG2.

STAG2, a component of the mitotically essential cohesin complex, is highly mutated in several different tumour types, including glioblastoma and bladder cancer. Whereas cohesin has roles in many cancer-related pathways, such as chromosome instability, DNA repair and gene expression, the complex nature of cohesin function has made it difficult to determine how STAG2 loss might either promote tumorigenesis or be leveraged therapeutically across divergent cancer types. Here, we have performed whole-genome CRISPR-Cas9 screens for STAG2-dependent genetic interactions in three distinct cellular backgrounds. Surprisingly, STAG1, the paralog of STAG2, was the only negative genetic interaction that was shared across all three backgrounds. We also uncovered a paralogous synthetic lethal mechanism behind a genetic interaction between STAG2 and the iron regulatory gene IREB2 Finally, investigation of an unusually strong context-dependent genetic interaction in HAP1 cells revealed factors that could be important for alleviating cohesin loading stress. Together, our results reveal new facets of STAG2 and cohesin function across a variety of genetic contexts.

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 understood1-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.

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.

Whole transcriptome analysis reveals differential gene expression profile reflecting macrophage polarization in response to influenza A H5N1 virus infection.

BACKGROUND: Avian influenza A H5N1 virus can cause lethal disease in humans. The virus can trigger severe pneumonia and lead to acute respiratory distress syndrome. Data from clinical, in vitro and in vivo suggest that virus-induced cytokine dysregulation could be a contributory factor to the pathogenesis of human H5N1 disease. However, the precise mechanism of H5N1 infection eliciting the unique host response are still not well understood. METHODS: To obtain a better understanding of the molecular events at the earliest time points, we used RNA-Seq to quantify and compare the host mRNA and miRNA transcriptomes induced by the highly pathogenic influenza A H5N1 (A/Vietnam/3212/04) or low virulent H1N1 (A/Hong Kong/54/98) viruses in human monocyte-derived macrophages at 1-, 3-, and 6-h post infection. RESULTS: Our data reveals that two macrophage populations corresponding to M1 (classically activated) and M2 (alternatively activated) macrophage subtypes respond distinctly to H5N1 virus infection when compared to H1N1 virus or mock infection, a distinction that could not be made from previous microarray studies. When this confounding variable is considered in our statistical model, a clear set of dysregulated genes and pathways emerges specifically in H5N1 virus-infected macrophages at 6-h post infection, whilst was not found with H1N1 virus infection. Furthermore, altered expression of genes in these pathways, which have been previously implicated in viral host response, occurs specifically in the M1 subtype. We observe a significant up-regulation of genes in the RIG-I-like receptor signaling pathway. In particular, interferons, and interferon-stimulated genes are broadly affected. The negative regulators of interferon signaling, the suppressors of cytokine signaling, SOCS-1 and SOCS-3, were found to be markedly up-regulated in the initial round of H5N1 virus replication. Elevated levels of these suppressors could lead to the eventual suppression of cellular antiviral genes, contributing to pathophysiology of H5N1 virus infection. CONCLUSIONS: Our study provides important mechanistic insights into the understanding of H5N1 viral pathogenesis and the multi-faceted host immune responses. The dysregulated genes could be potential candidates as therapeutic targets for treating H5N1 disease.

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.

Recoding RNA editing of AZIN1 predisposes to hepatocellular carcinoma.

A better understanding of human hepatocellular carcinoma (HCC) pathogenesis at the molecular level will facilitate the discovery of tumor-initiating events. Transcriptome sequencing revealed that adenosine-to-inosine (A→I) RNA editing of AZIN1 (encoding antizyme inhibitor 1) is increased in HCC specimens. A→I editing of AZIN1 transcripts, specifically regulated by ADAR1 (encoding adenosine deaminase acting on RNA-1), results in a serine-to-glycine substitution at residue 367 of AZIN1, located in ß-strand 15 (ß15) and predicted to cause a conformational change, induced a cytoplasmic-to-nuclear translocation and conferred gain-of-function phenotypes that were manifested by augmented tumor-initiating potential and more aggressive behavior. Compared with wild-type AZIN1 protein, the edited form has a stronger affinity to antizyme, and the resultant higher AZIN1 protein stability promotes cell proliferation through the neutralization of antizyme-mediated degradation of ornithine decarboxylase (ODC) and cyclin D1 (CCND1). Collectively, A→I RNA editing of AZIN1 may be a potential driver in the pathogenesis of human cancers, particularly HCC.

A network of multi-tasking proteins at the DNA replication fork preserves genome stability.

To elucidate the network that maintains high fidelity genome replication, we have introduced two conditional mutant alleles of DNA2, an essential DNA replication gene, into each of the approximately 4,700 viable yeast deletion mutants and determined the fitness of the double mutants. Fifty-six DNA2-interacting genes were identified. Clustering analysis of genomic synthetic lethality profiles of each of 43 of the DNA2-interacting genes defines a network (consisting of 322 genes and 876 interactions) whose topology provides clues as to how replication proteins coordinate regulation and repair to protect genome integrity. The results also shed new light on the functions of the query gene DNA2, which, despite many years of study, remain controversial, especially its proposed role in Okazaki fragment processing and the nature of its in vivo substrates. Because of the multifunctional nature of virtually all proteins at the replication fork, the meaning of any single genetic interaction is inherently ambiguous. The multiplexing nature of the current studies, however, combined with follow-up supporting experiments, reveals most if not all of the unique pathways requiring Dna2p. These include not only Okazaki fragment processing and DNA repair but also chromatin dynamics.

Global mapping of the yeast genetic interaction network.

A genetic interaction network containing approximately 1000 genes and approximately 4000 interactions was mapped by crossing mutations in 132 different query genes into a set of approximately 4700 viable gene yeast deletion mutants and scoring the double mutant progeny for fitness defects. Network connectivity was predictive of function because interactions often occurred among functionally related genes, and similar patterns of interactions tended to identify components of the same pathway. The genetic network exhibited dense local neighborhoods; therefore, the position of a gene on a partially mapped network is predictive of other genetic interactions. Because digenic interactions are common in yeast, similar networks may underlie the complex genetics associated with inherited phenotypes in other organisms.

Interaction between a Ras and a Rho GTPase couples selection of a growth site to the development of cell polarity in yeast.

Polarized cell growth requires the coupling of a defined spatial site on the cell cortex to the apparatus that directs the establishment of cell polarity. In the budding yeast Saccharomyces cerevisiae, the Ras-family GTPase Rsr1p/Bud1p and its regulators select the proper site for bud emergence on the cell cortex. The Rho-family GTPase Cdc42p and its associated proteins then establish an axis of polarized growth by triggering an asymmetric organization of the actin cytoskeleton and secretory apparatus at the selected bud site. We explored whether a direct linkage exists between the Rsr1p/Bud1p and Cdc42p GTPases. Here we show specific genetic interactions between RSR1/BUD1 and particular cdc42 mutants defective in polarity establishment. We also show that Cdc42p coimmunoprecipitated with Rsr1p/Bud1p from yeast extracts. In vitro studies indicated a direct interaction between Rsr1p/Bud1p and Cdc42p, which was enhanced by Cdc24p, a guanine nucleotide exchange factor for Cdc42p. Our findings suggest that Cdc42p interacts directly with Rsr1p/Bud1p in vivo, providing a novel mechanism by which direct contact between a Ras-family GTPase and a Rho-family GTPase links the selection of a growth site to polarity establishment.

A combined experimental and computational strategy to define protein interaction networks for peptide recognition modules.

Peptide recognition modules mediate many protein-protein interactions critical for the assembly of macromolecular complexes. Complete genome sequences have revealed thousands of these domains, requiring improved methods for identifying their physiologically relevant binding partners. We have developed a strategy combining computational prediction of interactions from phage-display ligand consensus sequences with large-scale two-hybrid physical interaction tests. Application to yeast SH3 domains generated a phage-display network containing 394 interactions among 206 proteins and a two-hybrid network containing 233 interactions among 145 proteins. Graph theoretic analysis identified 59 highly likely interactions common to both networks. Las17 (Bee1), a member of the Wiskott-Aldrich Syndrome protein (WASP) family of actin-assembly proteins, showed multiple SH3 interactions, many of which were confirmed in vivo by coimmunoprecipitation.