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.

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.

Identifying human pre-mRNA cleavage and polyadenylation factors by genome-wide CRISPR screens using a dual fluorescence readthrough reporter.

Messenger RNA precursors (pre-mRNA) generally undergo 3' end processing by cleavage and polyadenylation (CPA), which is specified by a polyadenylation site (PAS) and adjacent RNA sequences and regulated by a large variety of core and auxiliary CPA factors. To date, most of the human CPA factors have been discovered through biochemical and proteomic studies. However, genetic identification of the human CPA factors has been hampered by the lack of a reliable genome-wide screening method. We describe here a dual fluorescence readthrough reporter system with a PAS inserted between two fluorescent reporters. This system enables measurement of the efficiency of 3' end processing in living cells. Using this system in combination with a human genome-wide CRISPR/Cas9 library, we conducted a screen for CPA factors. The screens identified most components of the known core CPA complexes and other known CPA factors. The screens also identified CCNK/CDK12 as a potential core CPA factor, and RPRD1B as a CPA factor that binds RNA and regulates the release of RNA polymerase II at the 3' ends of genes. Thus, this dual fluorescence reporter coupled with CRISPR/Cas9 screens reliably identifies bona fide CPA factors and provides a platform for investigating the requirements for CPA in various contexts.

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.

N-terminal acetylation shields proteins from degradation and promotes age-dependent motility and longevity.

Most eukaryotic proteins are N-terminally acetylated, but the functional impact on a global scale has remained obscure. Using genome-wide CRISPR knockout screens in human cells, we reveal a strong genetic dependency between a major N-terminal acetyltransferase and specific ubiquitin ligases. Biochemical analyses uncover that both the ubiquitin ligase complex UBR4-KCMF1 and the acetyltransferase NatC recognize proteins bearing an unacetylated N-terminal methionine followed by a hydrophobic residue. NatC KO-induced protein degradation and phenotypes are reversed by UBR knockdown, demonstrating the central cellular role of this interplay. We reveal that loss of Drosophila NatC is associated with male sterility, reduced longevity, and age-dependent loss of motility due to developmental muscle defects. Remarkably, muscle-specific overexpression of UbcE2M, one of the proteins targeted for NatC KO-mediated degradation, suppresses defects of NatC deletion. In conclusion, NatC-mediated N-terminal acetylation acts as a protective mechanism against protein degradation, which is relevant for increased longevity and motility.

Golgi targeting of ARF-like GTPase Arl3p requires its Nalpha-acetylation and the integral membrane protein Sys1p.

Myristoylation of ARF family GTPases is required for their association with Golgi and endosomal membranes, where they regulate protein sorting and the lipid composition of these organelles. The Golgi-localized ARF-like GTPase Arl3p/ARP lacks a myristoylation signal, indicating that its targeting mechanism is distinct from myristoylated ARFs. We demonstrate that acetylation of the N-terminal methionine of Arl3p requires the NatC N(alpha)-acetyltransferase and that this modification is required for its Golgi localization. Chemical crosslinking and fluorescence microscopy experiments demonstrate that localization of Arl3p also requires Sys1p, a Golgi-localized integral membrane protein, which may serve as a receptor for acetylated Arl3p.

Ammonia-oxidizing archaea and ammonia-oxidizing bacteria in six full-scale wastewater treatment bioreactors.

In this study, dideoxy sequencing and 454 high-throughput sequencing were used to analyze diversities of the ammonia monooxygenase (amoA) genes and the 16S rRNA genes of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) in six municipal wastewater treatment plants. The results showed that AOB amoA genes were quite diverse in different wastewater treatment plants while the 16S rRNA genes were relatively conserved. Based on the observed complexity of amoA and 16S rRNA genes, most of the AOB can be assigned to the Nitrosomonas genus, with Nitrosomonas ureae, Nitrosomonas oligotropha, Nitrosomonas marina, and Nitrosomonas aestuarii being the four most dominant species. From the sequences of the AOA amoA genes, most AOA observed in this study belong to the CGI.1b group, i.e., the soil lineage. The AOB amoA and 16S rRNA genes were quantified by quantitative PCR and 454 high-throughput pyrosequencing, respectively. Although the results from the two approaches show some disconcordance, they both indicated that the abundance of AOB in activated sludge was very low.

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.

Sequential and distinct roles of the cadherin domain-containing protein Axl2p in cell polarization in yeast cell cycle.

Polarization of cell growth along a defined axis is essential for the generation of cell and tissue polarity. In the budding yeast Saccharomyces cerevisiae, Axl2p plays an essential role in polarity-axis determination, or more specifically, axial budding in MATa or alpha cells. Axl2p is a type I membrane glycoprotein containing four cadherin-like motifs in its extracellular domain. However, it is not known when and how Axl2p functions together with other components of the axial landmark, such as Bud3p and Bud4p, to direct axial budding. Here, we show that the recruitment of Axl2p to the bud neck after S/G2 phase of the cell cycle depends on Bud3p and Bud4p. This recruitment is mediated via an interaction between Bud4p and the central region of the Axl2p cytoplasmic tail. This region of Axl2p, together with its N-terminal region and its transmembrane domain, is sufficient for axial budding. In addition, our work demonstrates a previously unappreciated role for Axl2p. Axl2p interacts with Cdc42p and other polarity-establishment proteins, and it regulates septin organization in late G1 independently of its role in polarity-axis determination. Together, these results suggest that Axl2p plays sequential and distinct roles in the regulation of cellular morphogenesis in yeast cell cycle.

The GTPase Arf1p and the ER to Golgi cargo receptor Erv14p cooperate to recruit the golgin Rud3p to the cis-Golgi.

Rud3p is a coiled-coil protein of the yeast cis-Golgi. We find that Rud3p is localized to the Golgi via a COOH-terminal domain that is distantly related to the GRIP domain that recruits several coiled-coil proteins to the trans-Golgi by binding the small Arf-like GTPase Arl1p. In contrast, Rud3p binds to the GTPase Arf1p via this COOH-terminal "GRIP-related Arf-binding" (GRAB) domain. Deletion of RUD3 is lethal in the absence of the Golgi GTPase Ypt6p, and a screen of other mutants showing a similar genetic interaction revealed that Golgi targeting of Rud3p also requires Erv14p, a cargo receptor that cycles between the endoplasmic reticulum and Golgi. The one human protein with a GRAB domain, GMAP-210 (CEV14/Trip11/Trip230), is known to be on the cis-Golgi, but the COOH-terminal region that contains the GRAB domain has been reported to bind to centrosomes and gamma-tubulin (Rios, R.M, A. Sanchis, A.M. Tassin, C. Fedriani, and M. Bornens. 2004. Cell. 118:323-335). In contrast, we find that this region binds to the Golgi in a GRAB domain-dependent manner, suggesting that GMAP-210 may not link the Golgi to gamma-tubulin and centrosomes.

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.

Pooled Lentiviral CRISPR-Cas9 Screens for Functional Genomics in Mammalian Cells.

CRISPR-Cas9 technology provides a simple way to introduce targeted mutations into mammalian cells to induce loss-of-function phenotypes. The CRISPR-Cas9 system has now successfully been applied for genetic screens in many cell types, providing a powerful tool for functional genomics with manifold applications. Genome-wide guide-RNA (gRNA) libraries allow facile generation of a pool of cells, each harboring a gene knockout mutation that can be used for the study of gene function, pathway analysis or the identification of genes required for cellular fitness. Furthermore, CRISPR genetic screens can be applied for the discovery of genes whose knockout sensitizes cells to drug treatments or mediates drug resistance. Here, we provide a detailed protocol discussing the necessary steps for the successful performance of pooled CRISPR-Cas9 screens.

Essential Gene Profiles for Human Pluripotent Stem Cells Identify Uncharacterized Genes and Substrate Dependencies.

Human pluripotent stem cells (hPSCs) provide an invaluable tool for modeling diseases and hold promise for regenerative medicine. For understanding pluripotency and lineage differentiation mechanisms, a critical first step involves systematically cataloging essential genes (EGs) that are indispensable for hPSC fitness, defined as cell reproduction in this study. To map essential genetic determinants of hPSC fitness, we performed genome-scale loss-of-function screens in an inducible Cas9 H1 hPSC line cultured on feeder cells and laminin to identify EGs. Among these, we found FOXH1 and VENTX, genes that encode transcription factors previously implicated in stem cell biology, as well as an uncharacterized gene, C22orf43/DRICH1. hPSC EGs are substantially different from other human model cell lines, and EGs in hPSCs are highly context dependent with respect to different growth substrates. Our CRISPR screens establish parameters for genome-wide screens in hPSCs, which will facilitate the characterization of unappreciated genetic regulators of hPSC biology.

Inorganic phosphate deprivation causes tRNA nuclear accumulation via retrograde transport in Saccharomyces cerevisiae.

Nuclear export of tRNA is an essential eukaryotic function, yet the one known yeast tRNA nuclear exporter, Los1, is nonessential. Moreover recent studies have shown that tRNAs can move retrograde from the cytosol to the nucleus by an undefined process. Therefore, additional gene products involved in tRNA nucleus-cytosol dynamics have yet to be identified. Synthetic genetic array (SGA) analysis was employed to identify proteins involved in Los1-independent tRNA transport and in regulating tRNA nucleus-cytosol distribution. These studies uncovered synthetic interactions between los1Delta and pho88Delta involved in inorganic phopsphate uptake. Further analysis revealed that inorganic phosphate deprivation causes transient, temperature-dependent nuclear accumulation of mature cytoplasmic tRNA within nuclei via a Mtr10- and retrograde-dependent pathway, providing a novel connection between tRNA subcellular dynamics and phosphate availability.

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.

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.

Synthetic lethal analysis implicates Ste20p, a p21-activated potein kinase, in polarisome activation.

The p21-activated kinases Ste20p and Cla4p carry out undefined functions that are essential for viability during budding in Saccharomyces cerevisiae. To gain insight into the roles of Ste20p, we have used a synthetic lethal mutant screen to identify additional genes that are required in the absence of Cla4p. Altogether, we identified 65 genes, including genes with roles in cell polarity, mitosis, and cell wall maintenance. Herein, we focus on a set that defines a function carried out by Bni1p and several of its interacting proteins. We found that Bni1p and a group of proteins that complex with Bni1p (Bud6p, Spa2p, and Pea2p) are essential in a cla4delta mutant background. Bni1p, Bud6p, Spa2, and Pea2p are members of a group of polarity determining proteins referred to as the polarisome. Loss of polarisome proteins from a cla4delta strain causes cells to form elongated buds that have mislocalized septin rings. In contrast, other proteins that interact with or functionally associate with Bni1p and have roles in nuclear migration and cytokinesis, including Num1p and Hof1p, are not essential in the absence of Cla4p. Finally, we have found that Bni1p is phosphorylated in vivo, and a substantial portion of this phosphorylation is dependent on STE20. Together, these results suggest that one function of Ste20p may be to activate the polarisome complex by phosphorylation of Bni1p.

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.

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.

De Novo Genome and Transcriptome Assembly of the Canadian Beaver (Castor canadensis).

The Canadian beaver (Castor canadensis) is the largest indigenous rodent in North America. We report a draft annotated assembly of the beaver genome, the first for a large rodent and the first mammalian genome assembled directly from uncorrected and moderate coverage (< 30 ×) long reads generated by single-molecule sequencing. The genome size is 2.7 Gb estimated by k-mer analysis. We assembled the beaver genome using the new Canu assembler optimized for noisy reads. The resulting assembly was refined using Pilon supported by short reads (80 ×) and checked for accuracy by congruency against an independent short read assembly. We scaffolded the assembly using the exon-gene models derived from 9805 full-length open reading frames (FL-ORFs) constructed from the beaver leukocyte and muscle transcriptomes. The final assembly comprised 22,515 contigs with an N50 of 278,680 bp and an N50-scaffold of 317,558 bp. Maximum contig and scaffold lengths were 3.3 and 4.2 Mb, respectively, with a combined scaffold length representing 92% of the estimated genome size. The completeness and accuracy of the scaffold assembly was demonstrated by the precise exon placement for 91.1% of the 9805 assembled FL-ORFs and 83.1% of the BUSCO (Benchmarking Universal Single-Copy Orthologs) gene set used to assess the quality of genome assemblies. Well-represented were genes involved in dentition and enamel deposition, defining characteristics of rodents with which the beaver is well-endowed. The study provides insights for genome assembly and an important genomics resource for Castoridae and rodent evolutionary biology.