A High-Resolution Genome-Wide CRISPR/Cas9 Viability Screen Reveals Structural Features and Contextual Diversity of the Human Cell-Essential Proteome

DDB1 engagement defines the selectivity of S656 analogs for cyclin K degradation over CDK inhibition

In efforts to identify additional therapeutic targets for Acute Myeloid Leukemia (AML), we performed a high-throughput screen that includes 56 primary specimens tested with 10,000 structurally diverse small molecules. One specific hit, called S656 acts as a molecular glue degrader (MGD), that mediates the CRL4-dependent proteolysis of cyclin K. Structurally, S656 features a moiety that binds to the ATP binding site of cyclin-dependent kinases (CDKs), allowing the recruitment of the CDK12-cyclin K complex, along with a binding site for DDB1 bridging the CRL4 complex. Structure activity relationship studies reveal that minimal modifications to the dimethylaniline moiety of S656 improve its cyclin K MGD function over CDK inhibition by promoting DDB1 engagement. This includes full occupation of the DDB1 pocket, preferably with hydrophobic terminal groups, and cation-π interaction with Arg928. Additionally, we demonstrate that despite structural diversity, cyclin K degraders exhibit similar functional activity in AML which is distinct from direct CDK12 inhibition.

Genome-aware annotation of CRISPR guides validates targets in variant cell lines and enhances discovery in screens

BACKGROUND

CRISPR-Cas9 technology has revolutionised genetic screens and can inform on gene essentiality and chemo-genetic interactions. It is easily deployed and widely supported with many pooled CRISPR libraries available commercially. However, discrepancies between the reference genomes used in the design of those CRISPR libraries and the cell line under investigation can lead to loss of signal or introduction of bias. The problem is particularly acute when dealing with variant cell lines such as cancer cell lines.

RESULTS

Here, we present an algorithm, EXOme-guided Re-annotation of nuCleotIde SEquences (Exorcise), which uses sequence search to detect and correct mis-annotations in CRISPR libraries. Exorcise verifies the presence of CRISPR targets in the target genome and applies corrections to CRISPR libraries using existing exome annotations. We applied Exorcise to re-annotate guides in pooled CRISPR libraries available on Addgene and found that libraries designed on a more permissive reference sequence had more mis-annotations. In simulated CRISPR screens, we modelled common mis-annotations and found that they adversely affect discovery of hits in the intermediate range. We then confirmed this by applying Exorcise on datasets from Dependency Map (DepMap) and the DNA Damage Response CRISPR Screen Viewer (DDRcs), where we found improved discovery power upon Exorcise while retaining the strongest hits.

CONCLUSIONS

Pooled CRISPR libraries map guide sequences to genes and these mappings might not be ready to use due to permissive library design or investigating a variant cell line. By re-annotating CRISPR guides, Exorcise focuses CRISPR experiments towards the genome of the cell line under investigation. Exorcise can be applied at the library design stage or the analysis stage and allows post hoc re-analysis of completed screens. It is available under a Creative Commons Zero v1.0 Universal licence at https://github.com/SimonLammmm/exorcise .

Targeting the dependence on PIK3C3-mTORC1 signaling in dormancy-prone breast cancer cells blunts metastasis initiation

Halting breast cancer metastatic relapses following primary tumor removal and the clinical dormant phase, remains challenging, due to a lack of specific vulnerabilities to target during dormancy. To address this, we conducted genome-wide CRISPR screens on two breast cancer cell lines with distinct dormancy properties: 4T1 (short-term dormancy) and 4T07 (prolonged dormancy). We discovered that loss of class-III PI3K, Pik3c3, revealed a unique vulnerability in 4T07 cells. Surprisingly, dormancy-prone 4T07 cells exhibited higher mTORC1 activity than 4T1 cells, due to lysosome-dependent signaling occurring at the cell periphery. Pharmacological inhibition of Pik3c3 counteracted this phenotype in 4T07 cells, and selectively reduced metastasis burden only in the 4T07 dormancy-prone model. This mechanism was also detected in human breast cancer cell lines in addition to a breast cancer patient-derived xenograft supporting that it may be relevant in humans. Our findings suggest dormant cancer cell-initiated metastasis may be prevented in patients carrying tumor cells that display PIK3C3-peripheral lysosomal signaling to mTORC1.

Statement of Significance

We reveal that dormancy-prone breast cancer cells depend on the class III PI3K to mediate a constant peripheral lysosomal positioning and mTORC1 hyperactivity. Targeting this pathway might blunt breast cancer metastasis.

Genome-scale exon perturbation screens uncover exons critical for cell fitness

CRISPR-Cas technology has transformed functional genomics, yet understanding of how individual exons differentially shape cellular phenotypes remains limited. Here, we optimized and conducted massively parallel exon deletion and splice-site mutation screens in human cell lines to identify exons that regulate cellular fitness. Fitness-promoting exons are prevalent in essential and highly expressed genes and commonly overlap with protein domains and interaction interfaces. Conversely, fitness-suppressing exons are enriched in nonessential genes, exhibiting lower inclusion levels, and overlap with intrinsically disordered regions and disease-associated mutations. In-depth mechanistic investigation of the screen-hit TAF5 alternative exon-8 revealed that its inclusion is required for assembly of the TFIID general transcription initiation complex, thereby regulating global gene expression output. Collectively, our orthogonal exon perturbation screens established a comprehensive repository of phenotypically important exons and uncovered regulatory mechanisms governing cellular fitness and gene expression.

The anticancer potential of the CLK kinases inhibitors 1C8 and GPS167 revealed by their impact on the epithelial-mesenchymal transition and the antiviral immune response

The diheteroarylamide-based compound 1C8 and the aminothiazole carboxamide-related compound GPS167 inhibit the CLK kinases, and affect the proliferation of a broad range of cancer cell lines. A chemogenomic screen previously performed with GPS167 revealed that the depletion of components associated with mitotic spindle assembly altered sensitivity to GPS167. Here, a similar screen performed with 1C8 also established the impact of components involved in mitotic spindle assembly. Accordingly, transcriptome analyses of cells treated with 1C8 and GPS167 indicated that the expression and RNA splicing of transcripts encoding mitotic spindle assembly components were affected. The functional relevance of the microtubule connection was confirmed by showing that subtoxic concentrations of drugs affecting mitotic spindle assembly increased sensitivity to GPS167. 1C8 and GPS167 impacted the expression and splicing of transcripts in pathways relevant to tumor progression, including MYC targets and the epithelial mesenchymal transition (EMT). Finally, 1C8 and GPS167 altered the expression and alternative splicing of transcripts involved in the antiviral immune response. Consistent with this observation, depleting the double-stranded RNA sensor DHX33 suppressed GPS167-mediated cytotoxicity on HCT116 cells. Our study uncovered molecular mechanisms through which 1C8 and GPS167 affect cancer cell proliferation as well as processes critical for metastasis.

A CRISPR base editing approach for the functional assessment of telomere biology disorder-related genes in human health and aging

Telomere Biology Disorders (TBDs) are a group of rare diseases characterized by the presence of short and/or dysfunctional telomeres. They comprise a group of bone marrow failure syndromes, idiopathic pulmonary fibrosis, and liver disease, among other diseases. Genetic alterations (variants) in the genes responsible for telomere homeostasis have been linked to TBDs. Despite the number of variants already identified as pathogenic, an even more significant number must be better understood. The study of TBDs is challenging since identifying these variants is difficult due to their rareness, it is hard to predict their impact on the disease onset, and there are not enough samples to study. Most of our knowledge about pathogenic variants comes from assessing telomerase activity from patients and their relatives affected by a TBD. However, we still lack a cell-based model to identify new variants and to study the long-term impact of such variants on the genes involved in TBDs. Herein, we present a cell-based model using CRISPR base editing to mutagenize the endogenous alleles of 21 genes involved in telomere biology. We identified key residues in the genes encoding 17 different proteins impacting cell growth. We provide functional evidence for variants of uncertain significance in patients with TBDs. We also identified variants resistant to telomerase inhibition that, similar to cells expressing wild-type telomerase, exhibited increased tumorigenic potential using an in vitro tumour growth assay. We believe that such cell-based approaches will significantly advance our understanding of the biology of TBDs and may contribute to the development of new therapies for this group of diseases.

SF3B1 mutations provide genetic vulnerability to copper ionophores in human acute myeloid leukemia

In a phenotypical screen of 56 acute myeloid leukemia (AML) patient samples and using a library of 10,000 compounds, we identified a hit with increased sensitivity toward SF3B1-mutated and adverse risk AMLs. Through structure-activity relationship studies, this hit was optimized into a potent, specific, and nongenotoxic molecule called UM4118. We demonstrated that UM4118 acts as a copper ionophore that initiates a mitochondrial-based noncanonical form of cell death known as cuproptosis. CRISPR-Cas9 loss-of-function screen further revealed that iron-sulfur cluster (ISC) deficiency enhances copper-mediated cell death. Specifically, we found that loss of the mitochondrial ISC transporter ABCB7 is synthetic lethal to UM4118. ABCB7 is misspliced and down-regulated in SF3B1-mutated leukemia, creating a vulnerability to copper ionophores. Accordingly, ABCB7 overexpression partially rescued SF3B1-mutated cells to copper overload. Together, our work provides mechanistic insights that link ISC deficiency to cuproptosis, as exemplified by the high sensitivity of SF3B1-mutated AMLs. We thus propose SF3B1 mutations as a biomarker for future copper ionophore-based therapies.

Predicting host-based, synthetic lethal antiviral targets from omics data

Traditional antiviral therapies often have limited effectiveness due to toxicity and the emergence of drug resistance. Host-based antivirals are an alternative, but can cause nonspecific effects. Recent evidence shows that virus-infected cells can be selectively eliminated by targeting synthetic lethal (SL) partners of proteins disrupted by viral infection. Thus, we hypothesized that genes depleted in CRISPR knockout (KO) screens of virus-infected cells may be enriched in SL partners of proteins altered by infection. To investigate this, we established a computational pipeline predicting antiviral SL drug targets. First, we identified SARS-CoV-2-induced changes in gene products via a large compendium of omics data. Second, we identified SL partners for each altered gene product. Last, we screened CRISPR KO data for SL partners required for cell viability in infected cells. Despite differences in virus-induced alterations detected by various omics data, they share many predicted SL targets, with significant enrichment in CRISPR KO-depleted datasets. Our comparison of SARS-CoV-2 and influenza infection data revealed potential broad-spectrum, host-based antiviral SL targets. This suggests that CRISPR KO data are replete with common antiviral targets due to their SL relationship with virus-altered states and that such targets can be revealed from analysis of omics datasets and SL predictions.

Active growth signaling promotes senescence and cancer cell sensitivity to CDK7 inhibition

Tumor growth is driven by continued cellular growth and proliferation. Cyclin-dependent kinase 7's (CDK7) role in activating mitotic CDKs and global gene expression makes it therefore an attractive target for cancer therapies. However, what makes cancer cells particularly sensitive to CDK7 inhibition (CDK7i) remains unclear. Here, we address this question. We show that CDK7i, by samuraciclib, induces a permanent cell-cycle exit, known as senescence, without promoting DNA damage signaling or cell death. A chemogenetic genome-wide CRISPR knockout screen identified that active mTOR (mammalian target of rapamycin) signaling promotes samuraciclib-induced senescence. mTOR inhibition decreases samuraciclib sensitivity, and increased mTOR-dependent growth signaling correlates with sensitivity in cancer cell lines. Reverting a growth-promoting mutation in PIK3CA to wild type decreases sensitivity to CDK7i. Our work establishes that enhanced growth alone promotes CDK7i sensitivity, providing an explanation for why some cancers are more sensitive to CDK inhibition than normally growing cells.

Visualizing the nucleoplasmic maturation of human pre-60S ribosomal particles

Eukaryotic ribosome assembly is a highly orchestrated process that involves over two hundred protein factors. After early assembly events on nascent rRNA in the nucleolus, pre-60S particles undergo continuous maturation steps in the nucleoplasm, and prepare for nuclear export. Here, we report eleven cryo-EM structures of the nuclear pre-60S particles isolated from human cells through epitope-tagged GNL2, at resolutions of 2.8-4.3 Å. These high-resolution snapshots provide fine details for several major structural remodeling events at a virtual temporal resolution. Two new human nuclear factors, L10K and C11orf98, were also identified. Comparative structural analyses reveal that many assembly factors act as successive place holders to control the timing of factor association/dissociation events. They display multi-phasic binding properties for different domains and generate complex binding inter-dependencies as a means to guide the rRNA maturation process towards its mature conformation. Overall, our data reveal that nuclear assembly of human pre-60S particles is generally hierarchical with short branch pathways, and a few factors display specific roles as rRNA chaperones by confining rRNA helices locally to facilitate their folding, such as the C-terminal domain of SDAD1.

Predicting host-based, synthetic lethal antiviral targets from omics data

Traditional antiviral therapies often have limited effectiveness due to toxicity and development of drug resistance. Host-based antivirals, while an alternative, may lead to non-specific effects. Recent evidence shows that virus-infected cells can be selectively eliminated by targeting synthetic lethal (SL) partners of proteins disrupted by viral infection. Thus, we hypothesized that genes depleted in CRISPR KO screens of virus-infected cells may be enriched in SL partners of proteins altered by infection. To investigate this, we established a computational pipeline predicting SL drug targets of viral infections. First, we identified SARS-CoV-2-induced changes in gene products via a large compendium of omics data. Second, we identified SL partners for each altered gene product. Last, we screened CRISPR KO data for SL partners required for cell viability in infected cells. Despite differences in virus-induced alterations detected by various omics data, they share many predicted SL targets, with significant enrichment in CRISPR KO-depleted datasets. Comparing data from SARS-CoV-2 and influenza infections, we found possible broad-spectrum, host-based antiviral SL targets. This suggests that CRISPR KO data are replete with common antiviral targets due to their SL relationship with virus-altered states and that such targets can be revealed from analysis of omics datasets and SL predictions.

Monoallelic variation in DHX9, the gene encoding the DExH-box helicase DHX9, underlies neurodevelopment disorders and Charcot-Marie-Tooth disease

DExD/H-box RNA helicases (DDX/DHX) are encoded by a large paralogous gene family; in a subset of these human helicase genes, pathogenic variation causes neurodevelopmental disorder (NDD) traits and cancer. DHX9 encodes a BRCA1-interacting nuclear helicase regulating transcription, R-loops, and homologous recombination and exhibits the highest mutational constraint of all DDX/DHX paralogs but remains unassociated with disease traits in OMIM. Using exome sequencing and family-based rare-variant analyses, we identified 20 individuals with de novo, ultra-rare, heterozygous missense or loss-of-function (LoF) DHX9 variant alleles. Phenotypes ranged from NDDs to the distal symmetric polyneuropathy axonal Charcot-Marie-Tooth disease (CMT2). Quantitative Human Phenotype Ontology (HPO) analysis demonstrated genotype-phenotype correlations with LoF variants causing mild NDD phenotypes and nuclear localization signal (NLS) missense variants causing severe NDD. We investigated DHX9 variant-associated cellular phenotypes in human cell lines. Whereas wild-type DHX9 was restricted to the nucleus, NLS missense variants abnormally accumulated in the cytoplasm. Fibroblasts from an individual with an NLS variant also showed abnormal cytoplasmic DHX9 accumulation. CMT2-associated missense variants caused aberrant nucleolar DHX9 accumulation, a phenomenon previously associated with cellular stress. Two NDD-associated variants, p.Gly411Glu and p.Arg761Gln, altered DHX9 ATPase activity. The severe NDD-associated variant p.Arg141Gln did not affect DHX9 localization but instead increased R-loop levels and double-stranded DNA breaks. Dhx9-/- mice exhibited hypoactivity in novel environments, tremor, and sensorineural hearing loss. All together, these results establish DHX9 as a critical regulator of mammalian neurodevelopment and neuronal homeostasis.

PAF49: An RNA Polymerase I subunit essential for rDNA transcription and stabilization of PAF53

The application of genetic and biochemical techniques in yeast has informed our knowledge of transcription in mammalian cells. Such systems have allowed investigators to determine whether a gene was essential and to determine its function in rDNA transcription. However, there are significant differences in the nature of the transcription factors essential for transcription by Pol I in yeast and mammalian cells, and yeast RNA polymerase I contains 14 subunits while mammalian polymerase contains 13 subunits. We previously reported the adaptation of the auxin-dependent degron that enabled a combination of a "genetics-like" approach and biochemistry to study mammalian rDNA transcription. Using this system, we studied the mammalian orthologue of yeast RPA34.5, PAF49, and found that it is essential for rDNA transcription and cell division. The auxin-induced degradation of PAF49 induced nucleolar stress and the accumulation of P53. Interestingly, the auxin-induced degradation of AID-tagged PAF49 led to the degradation of its binding partner, PAF53, but not vice versa. A similar pattern of co-dependent expression was also found when we studied the non-essential, yeast orthologues. An analysis of the domains of PAF49 that are essential for rDNA transcription demonstrated a requirement for both the dimerization domain and an "arm" of PAF49 that interacts with PolR1B. Further, we demonstrate this interaction can be disrupted to inhibit Pol I transcription in normal and cancer cells which leads to the arrest of normal cells and cancer cell death. In summary, we have shown that both PAF53 and PAF49 are necessary for rDNA transcription and cell growth.

Genetic enhancers of partial PLK1 inhibition reveal hypersensitivity to kinetochore perturbations

Polo-like kinase 1 (PLK1) is a serine/threonine kinase required for mitosis and cytokinesis. As cancer cells are often hypersensitive to partial PLK1 inactivation, chemical inhibitors of PLK1 have been developed and tested in clinical trials. However, these small molecule inhibitors alone are not completely effective. PLK1 promotes numerous molecular and cellular events in the cell division cycle and it is unclear which of these events most crucially depend on PLK1 activity. We used a CRISPR-based genome-wide screening strategy to identify genes whose inactivation enhances cell proliferation defects upon partial chemical inhibition of PLK1. Genes identified encode proteins that are functionally linked to PLK1 in multiple ways, most notably factors that promote centromere and kinetochore function. Loss of the kinesin KIF18A or the outer kinetochore protein SKA1 in PLK1-compromised cells resulted in mitotic defects, activation of the spindle assembly checkpoint and nuclear reassembly defects. We also show that PLK1-dependent CENP-A loading at centromeres is extremely sensitive to partial PLK1 inhibition. Our results suggest that partial inhibition of PLK1 compromises the integrity and function of the centromere/kinetochore complex, rendering cells hypersensitive to different kinetochore perturbations. We propose that KIF18A is a promising target for combinatorial therapies with PLK1 inhibitors.

Using traditional machine learning and deep learning methods for on- and off-target prediction in CRISPR/Cas9: a review

CRISPR/Cas9 (Clustered Regularly Interspaced Short Palindromic Repeats and CRISPR-associated protein 9) is a popular and effective two-component technology used for targeted genetic manipulation. It is currently the most versatile and accurate method of gene and genome editing, which benefits from a large variety of practical applications. For example, in biomedicine, it has been used in research related to cancer, virus infections, pathogen detection, and genetic diseases. Current CRISPR/Cas9 research is based on data-driven models for on- and off-target prediction as a cleavage may occur at non-target sequence locations. Nowadays, conventional machine learning and deep learning methods are applied on a regular basis to accurately predict on-target knockout efficacy and off-target profile of given single-guide RNAs (sgRNAs). In this paper, we present an overview and a comparative analysis of traditional machine learning and deep learning models used in CRISPR/Cas9. We highlight the key research challenges and directions associated with target activity prediction. We discuss recent advances in the sgRNA-DNA sequence encoding used in state-of-the-art on- and off-target prediction models. Furthermore, we present the most popular deep learning neural network architectures used in CRISPR/Cas9 prediction models. Finally, we summarize the existing challenges and discuss possible future investigations in the field of on- and off-target prediction. Our paper provides valuable support for academic and industrial researchers interested in the application of machine learning methods in the field of CRISPR/Cas9 genome editing.

Viral protein engagement of GBF1 induces host cell vulnerability through synthetic lethality

Viruses co-opt host proteins to carry out their lifecycle. Repurposed host proteins may thus become functionally compromised; a situation analogous to a loss-of-function mutation. We term such host proteins as viral-induced hypomorphs. Cells bearing cancer driver loss-of-function mutations have successfully been targeted with drugs perturbing proteins encoded by the synthetic lethal (SL) partners of cancer-specific mutations. Similarly, SL interactions of viral-induced hypomorphs can potentially be targeted as host-based antiviral therapeutics. Here, we use GBF1, which supports the infection of many RNA viruses, as a proof-of-concept. GBF1 becomes a hypomorph upon interaction with the poliovirus protein 3A. Screening for SL partners of GBF1 revealed ARF1 as the top hit, disruption of which selectively killed cells that synthesize 3A alone or in the context of a poliovirus replicon. Thus, viral protein interactions can induce hypomorphs that render host cells selectively vulnerable to perturbations that leave uninfected cells otherwise unscathed. Exploiting viral-induced vulnerabilities could lead to broad-spectrum antivirals for many viruses, including SARS-CoV-2.

Bridging the splicing gap in human genetics with long-read RNA sequencing: finding the protein isoform drivers of disease

Aberrant splicing underlies many human diseases, including cancer, cardiovascular diseases and neurological disorders. Genome-wide mapping of splicing quantitative trait loci (sQTLs) has shown that genetic regulation of alternative splicing is widespread. However, identification of the corresponding isoform or protein products associated with disease-associated sQTLs is challenging with short-read RNA-seq, which cannot precisely characterize full-length transcript isoforms. Furthermore, contemporary sQTL interpretation often relies on reference transcript annotations, which are incomplete. Solutions to these issues may be found through integration of newly emerging long-read sequencing technologies. Long-read sequencing offers the capability to sequence full-length mRNA transcripts and, in some cases, to link sQTLs to transcript isoforms containing disease-relevant protein alterations. Here, we provide an overview of sQTL mapping approaches, the use of long-read sequencing to characterize sQTL effects on isoforms, the linkage of RNA isoforms to protein-level functions and comment on future directions in the field. Based on recent progress, long-read RNA sequencing promises to be part of the human disease genetics toolkit to discover and treat protein isoforms causing rare and complex diseases.

The histone methyltransferase SETD2 negatively regulates cell size

Cell size varies between cell types but is tightly regulated by cell intrinsic and extrinsic mechanisms. Cell size control is important for cell function, and changes in cell size are frequently observed in cancer. Here, we uncover a role for SETD2 in regulating cell size. SETD2 is a lysine methyltransferase and a tumor suppressor protein involved in transcription, RNA processing and DNA repair. At the molecular level, SETD2 is best known for associating with RNA polymerase II through its Set2-Rbp1 interacting (SRI) domain and methylating histone H3 on lysine 36 (H3K36) during transcription. Using multiple independent perturbation strategies, we identify SETD2 as a negative regulator of global protein synthesis rates and cell size. We provide evidence that overexpression of the H3K36 demethylase KDM4A or the oncohistone H3.3K36M also increase cell size. In addition, ectopic overexpression of a decoy SRI domain increased cell size, suggesting that the relevant substrate is engaged by SETD2 via its SRI domain. These data add a central role of SETD2 in regulating cellular physiology and warrant further studies on separating the different functions of SETD2 in cancer development.

CRISPR/Cas9 and genetic screens in malaria parasites: small genomes, big impact

The ∼30 Mb genomes of the Plasmodium parasites that cause malaria each encode ∼5000 genes, but the functions of the majority remain unknown. This is due to a paucity of functional annotation from sequence homology, which is compounded by low genetic tractability compared with many model organisms. In recent years technical breakthroughs have made forward and reverse genome-scale screens in Plasmodium possible. Furthermore, the adaptation of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR-Associated protein 9 (CRISPR/Cas9) technology has dramatically improved gene editing efficiency at the single gene level. Here, we review the arrival of genetic screens in malaria parasites to analyse parasite gene function at a genome-scale and their impact on understanding parasite biology. CRISPR/Cas9 screens, which have revolutionised human and model organism research, have not yet been implemented in malaria parasites due to the need for more complex CRISPR/Cas9 gene targeting vector libraries. We therefore introduce the reader to CRISPR-based screens in the related apicomplexan Toxoplasma gondii and discuss how these approaches could be adapted to develop CRISPR/Cas9 based genome-scale genetic screens in malaria parasites. Moreover, since more than half of Plasmodium genes are required for normal asexual blood-stage reproduction, and cannot be targeted using knockout methods, we discuss how CRISPR/Cas9 could be used to scale up conditional gene knockdown approaches to systematically assign function to essential genes.

Identification of a glioma functional network from gene fitness data using machine learning

Glioblastoma multiforme (GBM) is an aggressive form of brain tumours that remains incurable despite recent advances in clinical treatments. Previous studies have focused on sub-categorizing patient samples based on clustering various transcriptomic data. While functional genomics data are rapidly accumulating, there exist opportunities to leverage these data to decipher glioma-associated biomarkers. We sought to implement a systematic approach to integrating data from high throughput CRISPR-Cas9 screening studies with machine learning algorithms to infer a glioma functional network. We demonstrated the network significantly enriched various biological pathways and may play roles in glioma tumorigenesis. From densely connected glioma functional modules, we further predicted 12 potential Wnt/β-catenin signalling pathway targeted genes, including AARSD1, HOXB5, ITGA6, LRRC71, MED19, MED24, METTL11B, SMARCB1, SMARCE1, TAF6L, TENT5A and ZNF281. Cox regression modelling with these targets was significantly associated with glioma overall survival prognosis. Additionally, TRIB2 was identified as a glioma neoplastic cell marker in single-cell RNA-seq of GBM samples. This work establishes novel strategies for constructing functional networks to identify glioma biomarkers for the development of diagnosis and treatment in clinical practice.

What sequencing technologies can teach us about innate immunity

For years, we have taken a reductionist approach to understanding gene regulation through the study of one gene in one cell at a time. While this approach has been fruitful it is laborious and fails to provide a global picture of what is occurring in complex situations involving tightly coordinated immune responses. The emergence of whole-genome techniques provides a system-level view of a response and can provide a plethora of information on events occurring in a cell from gene expression changes to splicing changes and chemical modifications. As with any technology, this often results in more questions than answers, but this wealth of knowledge is providing us with an unprecedented view of what occurs inside our cells during an immune response. In this review, we will discuss the current RNA-sequencing technologies and what they are helping us learn about the innate immune system.

Post-synaptic scaffold protein TANC2 in psychiatric and somatic disease risk

Understanding the shared genetic aetiology of psychiatric and medical comorbidity in neurodevelopmental disorders (NDDs) could improve patient diagnosis, stratification and treatment options. Rare tetratricopeptide repeat, ankyrin repeat and coiled-coil containing 2 (TANC2)-disrupting variants were disease causing in NDD patients. The post-synaptic scaffold protein TANC2 is essential for dendrite formation in synaptic plasticity and plays an unclarified but critical role in development. We here report a novel homozygous-viable Tanc2-disrupted function model in which mutant mice were hyperactive and had impaired sensorimotor gating consistent with NDD patient psychiatric endophenotypes. Yet, a multi-systemic analysis revealed the pleiotropic effects of Tanc2 outside the brain, such as growth failure and hepatocellular damage. This was associated with aberrant liver function including altered hepatocellular metabolism. Integrative analysis indicates that these disrupted Tanc2 systemic effects relate to interaction with Hippo developmental signalling pathway proteins and will increase the risk for comorbid somatic disease. This highlights how NDD gene pleiotropy can augment medical comorbidity susceptibility, underscoring the benefit of holistic NDD patient diagnosis and treatment for which large-scale preclinical functional genomics can provide complementary pleiotropic gene function information.

Summary: Disruption of mouse Tanc2 causes brain and liver abnormality, increasing psychiatric and somatic disease risk long term, highlighting the benefit of holistic diagnosis and treatment approaches for human neurodevelopmental disorder.

Loss of peptidase D binding restores the tumor suppressor functions of oncogenic p53 mutants

Tumor suppressor p53, a critical regulator of cell fate, is frequently mutated in cancer. Mutation of p53 abolishes its tumor-suppressing functions or endows oncogenic functions. We recently found that p53 binds via its proline-rich domain to peptidase D (PEPD) and is activated when the binding is disrupted. The proline-rich domain in p53 is rarely mutated. Here, we show that oncogenic p53 mutants closely resemble p53 in PEPD binding but are transformed into tumor suppressors, rather than activated as oncoproteins, when their binding to PEPD is disrupted by PEPD knockdown. Once freed from PEPD, p53 mutants undergo multiple posttranslational modifications, especially lysine 373 acetylation, which cause them to refold and regain tumor suppressor activities that are typically displayed by p53. The reactivated p53 mutants strongly inhibit cancer cell growth in vitro and in vivo. Our study identifies a cellular mechanism for reactivation of the tumor suppressor functions of oncogenic p53 mutants.

Chemical Genetics Screen Identifies COPB2 Tool Compounds That Alters ER Stress Response and Induces RTK Dysregulation in Lung Cancer Cells

Activating mutations in the epidermal growth factor receptor (EGFR) are common driver mutations in non-small cell lung cancer (NSCLC). First, second and third generation EGFR tyrosine kinase inhibitors (TKIs) are effective at inhibiting mutant EGFR NSCLC, however, acquired resistance is a major issue, leading to disease relapse. Here, we characterize a small molecule, EMI66, an analog of a small molecule which we previously identified to inhibit mutant EGFR signalling via a novel mechanism of action. We show that EMI66 attenuates receptor tyrosine kinase (RTK) expression and signalling and alters the electrophoretic mobility of Coatomer Protein Complex Beta 2 (COPB2) protein in mutant EGFR NSCLC cells. Moreover, we demonstrate that EMI66 can alter the subcellular localization of EGFR and COPB2 within the early secretory pathway. Furthermore, we find that COPB2 knockdown reduces the growth of mutant EGFR lung cancer cells, alters the post-translational processing of RTKs, and alters the endoplasmic reticulum (ER) stress response pathway. Lastly, we show that EMI66 treatment also alters the ER stress response pathway and inhibits the growth of mutant EGFR lung cancer cells and organoids. Our results demonstrate that targeting of COPB2 with EMI66 presents a viable approach to attenuate mutant EGFR signalling and growth in NSCLC.

A novel class of inhibitors that target SRSF10 and promote p53-mediated cytotoxicity on human colorectal cancer cells

The elevated expression of the splicing regulator SRSF10 in metastatic colorectal cancer (CRC) stimulates the production of the pro-tumorigenic BCLAF1-L splice variant. We discovered a group of small molecules with an aminothiazole carboxamide core (GPS167, GPS192 and others) that decrease production of BCLAF1-L. While additional alternative splicing events regulated by SRSF10 are affected by GPS167/192 in HCT116 cells (e.g. in MDM4, WTAP, SLK1 and CLK1), other events are shifted in a SRSF10-independent manner (e.g. in MDM2, NAB2 and TRA2A). GPS167/192 increased the interaction of SRSF10 with the CLK1 and CLK4 kinases, leading us to show that GPS167/192 can inhibit CLK kinases preferentially impacting the activity of SRSF10. Notably, GPS167 impairs the growth of CRC cell lines and organoids, inhibits anchorage-independent colony formation, cell migration, and promotes cytoxicity in a manner that requires SRSF10 and p53. In contrast, GPS167 only minimally affects normal colonocytes and normal colorectal organoids. Thus, GPS167 reprograms the tumorigenic activity of SRSF10 in CRC cells to elicit p53-dependent apoptosis.

One-step multiplex toolkit for efficient generation of conditional gene silencing human cell lines

Loss-of-function analysis is one of the major arsenals we have for understanding gene functions in mammalian cells. For analysis of essential genes, the major challenge is to develop simple methodologies for tight and rapid inducible gene inactivation. One approach involves CRISPR-Cas9-mediated disruption of the endogenous locus in conjunction with the expression of a rescue construct, which can subsequently be turned off to produce a gene inactivation effect. Here we describe the development of a set of Sleeping Beauty transposon-based vectors for expressing auxin-inducible degron (AID)-tagged genes under the regulation of a tetracycline-controlled promoter. The dual transcriptional and degron-mediated post-translational regulation allows rapid and tight silencing of protein expression in mammalian cells. We demonstrated that both non-essential and essential genes could be targeted in human cell lines using a one-step transfection method. Moreover, multiple genes could be simultaneously or sequentially targeted, allowing inducible inactivation of multiple genes. These resources enable highly efficient generation of conditional gene silencing cell lines to facilitate functional studies of essential genes.

GuidePro: a multi-source ensemble predictor for prioritizing sgRNAs in CRISPR/Cas9 protein knockouts

MOTIVATION

The efficiency of CRISPR/Cas9-mediated protein knockout is determined by three factors: sequence-specific sgRNA activity, frameshift probability and the characteristics of targeted amino acids. A number of computational methods have been developed for predicting sgRNA efficiency from different perspectives. However, an integrative method that combines all three factors for rational sgRNA selection is still lacking.

RESULTS

We developed GuidePro, a two-layer ensemble predictor that enables the integration of multiple factors for the prioritization of sgRNAs in protein knockouts. Tested on independent datasets, GuidePro outperforms existing methods and demonstrates consistent superior performance in predicting phenotypes caused by protein loss-of-function, suggesting its robustness for prioritizing sgRNAs in various applications of CRISPR/Cas9 knockouts.

AVAILABILITY AND IMPLEMENTATION

GuidePro is available at https://github.com/MDhewei/GuidePro. A web application for prioritizing sgRNAs that target protein-coding genes in human, monkey and mouse genomes is available at https://bioinformatics.mdanderson.org/apps/GuidePro.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

SUMOylation- and GAR1-Dependent Regulation of Dyskerin Nuclear and Subnuclear Localization

The nuclear and subnuclear compartmentalization of the telomerase-associated protein and H/ACA ribonucleoprotein component dyskerin is an important although incompletely understood aspect of H/ACA ribonucleoprotein function. Four SUMOylation sites were previously identified in the C-terminal nuclear/nucleolar localization signal (N/NoLS) of dyskerin. We found that a cytoplasmic localized C-terminal truncation variant of dyskerin lacking most of the C-terminal N/NoLS represents an under-SUMOylated variant of dyskerin compared to wild-type dyskerin. We demonstrate that mimicking constitutive SUMOylation of dyskerin using a SUMO3 fusion construct can drive nuclear accumulation of this variant and that the SUMO site K467 in this N/NoLS is particularly important for the subnuclear localization of dyskerin to the nucleolus in a mature H/ACA complex assembly- and SUMO-dependent manner. We also characterize a novel SUMO-interacting motif in the mature H/ACA complex component GAR1 that mediates the interaction between dyskerin and GAR1. Mislocalization of dyskerin, either in the cytoplasm or excluded from the nucleolus, disrupts dyskerin function and leads to reduced interaction of dyskerin with the telomerase RNA. These data indicate a role for dyskerin C-terminal N/NoLS SUMOylation in regulating the nuclear and subnuclear localization of dyskerin, which is essential for dyskerin function as both a telomerase-associated protein and as an H/ACA ribonucleoprotein.

Structural basis for the final steps of human 40S ribosome maturation

Eukaryotic ribosomes consist of a small 40S and a large 60S subunit that are assembled in a highly coordinated manner. More than 200 factors ensure correct modification, processing and folding of ribosomal RNA and the timely incorporation of ribosomal proteins^1,2. Small subunit maturation ends in the cytosol, when the final rRNA precursor, 18S-E, is cleaved at site 3 by the endonuclease NOB1³. Previous structures of human 40S precursors have shown that NOB1 is kept in an inactive state by its partner PNO1⁴. The final maturation events, including the activation of NOB1 for the decisive rRNA-cleavage step and the mechanisms driving the dissociation of the last biogenesis factors have, however, remained unresolved. Here we report five cryo-electron microscopy structures of human 40S subunit precursors, which describe the compositional and conformational progression during the final steps of 40S assembly. Our structures explain the central role of RIOK1 in the displacement and dissociation of PNO1, which in turn allows conformational changes and activation of the endonuclease NOB1. In addition, we observe two factors, eukaryotic translation initiation factor 1A domain-containing protein (EIF1AD) and leucine-rich repeat-containing protein 47 (LRRC47), which bind to late pre-40S particles near RIOK1 and the central rRNA helix 44. Finally, functional data shows that EIF1AD is required for efficient assembly factor recycling and 18S-E processing. Our results thus enable a detailed understanding of the last steps in 40S formation in human cells and, in addition, provide evidence for principal differences in small ribosomal subunit formation between humans and the model organism Saccharomyces cerevisiae.

Viral protein engagement of GBF1 induces host cell vulnerability through synthetic lethality

Viruses co-opt host proteins to carry out their lifecycle. Repurposed host proteins may thus become functionally compromised; a situation analogous to a loss-of-function mutation. We term such host proteins viral-induced hypomorphs. Cells bearing cancer driver loss-of-function mutations have successfully been targeted with drugs perturbing proteins encoded by the synthetic lethal partners of cancer-specific mutations. Synthetic lethal interactions of viral-induced hypomorphs have the potential to be similarly targeted for the development of host-based antiviral therapeutics. Here, we use GBF1, which supports the infection of many RNA viruses, as a proof-of-concept. GBF1 becomes a hypomorph upon interaction with the poliovirus protein 3A. Screening for synthetic lethal partners of GBF1 revealed ARF1 as the top hit, disruption of which, selectively killed cells that synthesize poliovirus 3A. Thus, viral protein interactions can induce hypomorphs that render host cells vulnerable to perturbations that leave uninfected cells intact. Exploiting viral-induced vulnerabilities could lead to broad-spectrum antivirals for many viruses, including SARS-CoV-2.

SUMMARY

Using a viral-induced hypomorph of GBF1, Navare et al., demonstrate that the principle of synthetic lethality is a mechanism to selectively kill virus-infected cells.

SLC25A51 is a mammalian mitochondrial NAD+ transporter

Mitochondria require nicotinamide adenine dinucleotide (NAD⁺) in order to carry out the fundamental processes that fuel respiration and mediate cellular energy transduction. Mitochondrial NAD⁺ transporters have been identified in yeast and plants ^1,2 but their very existence is controversial in mammals ^3–5. Here we demonstrate that mammalian mitochondria are capable of taking up intact NAD⁺ and identify SLC25A51 (an essential ^6,7 mitochondrial protein of previously unknown function, also known as MCART1) as a mammalian mitochondrial NAD⁺ transporter. Loss of SLC25A51 decreases mitochondrial but not whole-cell NAD⁺ content, impairs mitochondrial respiration, and blocks the uptake of NAD⁺ into isolated mitochondria. Conversely, overexpression of SLC25A51 or a nearly identical paralog, SLC25A52, increases mitochondrial NAD⁺ levels and restores NAD⁺ uptake into yeast mitochondria lacking endogenous NAD⁺ transporters. Together, these findings identify SLC25A51 as the first transporter capable of importing NAD⁺ into mammalian mitochondria.

Effects of Different Dietary Vegetable Lipid Sources on Health Status in Nile Tilapia (Oreochromis niloticus): Haematological Indices, Immune Response Parameters and Plasma Proteome

This study aimed to investigate the effects of DLs, including palm oil (PO; an SFAs), linseed oil (LO; n-3 PUFAs) and soybean oil (SBO; n-6 PUFAs) on the health status of Nile tilapia (Oreochromis niloticus) during adulthood. Three experimental diets incorporating PO, LO or SBO were fed to adult Nile tilapia for a period of 90 days, and haematological and innate immune parameters were evaluated. Proteome analysis was also conducted to evaluate the effects of DLs on plasma proteins. The tested DLs had no significant effects on red blood cell (RBC) count, haematocrit, haemoglobin, and total immunoglobulin and lysozyme activity. Dietary LO led to increased alternative complement 50 activity (ACH50), and proteome analysis revealed that PO and SBO enhanced A2ML, suggesting that different DLs promote immune system via different processes. Dietary LO or SBO increased the expression of several proteins involved in coagulation activity such as KNG1, HRG and FGG. Increased HPX in fish fed with PO suggests that SFAs are utilised in heme lipid-oxidation. Overall, DLs with distinct fatty acids (FAs) affect several parameters corresponding to health status in Nile tilapia, and dietary LO and SBO seemed to strengthen health in this species.

Cell Cycle Profiling Reveals Protein Oscillation, Phosphorylation, and Localization Dynamics

The cell cycle is a highly conserved process involving the coordinated separation of a single cell into two daughter cells. To relate transcriptional regulation across the cell cycle with oscillatory changes in protein abundance and activity, we carried out a proteome- and phospho-proteome-wide mass spectrometry profiling. We compared protein dynamics with gene transcription, revealing many transcriptionally regulated G2 mRNAs that only produce a protein shift after mitosis. Integration of CRISPR/Cas9 survivability studies further highlighted proteins essential for cell viability. Analyzing the dynamics of phosphorylation events and protein solubility dynamics over the cell cycle, we characterize predicted phospho-peptide motif distributions and predict cell cycle-dependent translocating proteins, as exemplified by the S-adenosylmethionine synthase MAT2A. Our study implicates this enzyme in translocating to the nucleus after the G1/S-checkpoint, which enables epigenetic histone methylation maintenance during DNA replication. Taken together, this data set provides a unique integrated resource with novel insights on cell cycle dynamics.

Genetic interaction mapping and exon-resolution functional genomics with a hybrid Cas9-Cas12a platform

Systematic mapping of genetic interactions (GIs) and interrogation of the functions of sizable genomic segments in mammalian cells represent important goals of biomedical research. To advance these goals, we present a CRISPR (clustered regularly interspaced short palindromic repeats)-based screening system for combinatorial genetic manipulation that employs coexpression of CRISPR-associated nucleases 9 and 12a (Cas9 and Cas12a) and machine-learning-optimized libraries of hybrid Cas9-Cas12a guide RNAs. This system, named Cas Hybrid for Multiplexed Editing and screening Applications (CHyMErA), outperforms genetic screens using Cas9 or Cas12a editing alone. Application of CHyMErA to the ablation of mammalian paralog gene pairs reveals extensive GIs and uncovers phenotypes normally masked by functional redundancy. Application of CHyMErA in a chemogenetic interaction screen identifies genes that impact cell growth in response to mTOR pathway inhibition. Moreover, by systematically targeting thousands of alternative splicing events, CHyMErA identifies exons underlying human cell line fitness. CHyMErA thus represents an effective screening approach for GI mapping and the functional analysis of sizable genomic regions, such as alternative exons.

Loss of UGP2 in brain leads to a severe epileptic encephalopathy, emphasizing that bi-allelic isoform-specific start-loss mutations of essential genes can cause genetic diseases

Developmental and/or epileptic encephalopathies (DEEs) are a group of devastating genetic disorders, resulting in early-onset, therapy-resistant seizures and developmental delay. Here we report on 22 individuals from 15 families presenting with a severe form of intractable epilepsy, severe developmental delay, progressive microcephaly, visual disturbance and similar minor dysmorphisms. Whole exome sequencing identified a recurrent, homozygous variant (chr2:64083454A > G) in the essential UDP-glucose pyrophosphorylase (UGP2) gene in all probands. This rare variant results in a tolerable Met12Val missense change of the longer UGP2 protein isoform but causes a disruption of the start codon of the shorter isoform, which is predominant in brain. We show that the absence of the shorter isoform leads to a reduction of functional UGP2 enzyme in neural stem cells, leading to altered glycogen metabolism, upregulated unfolded protein response and premature neuronal differentiation, as modeled during pluripotent stem cell differentiation in vitro. In contrast, the complete lack of all UGP2 isoforms leads to differentiation defects in multiple lineages in human cells. Reduced expression of Ugp2a/Ugp2b in vivo in zebrafish mimics visual disturbance and mutant animals show a behavioral phenotype. Our study identifies a recurrent start codon mutation in UGP2 as a cause of a novel autosomal recessive DEE syndrome. Importantly, it also shows that isoform-specific start-loss mutations causing expression loss of a tissue-relevant isoform of an essential protein can cause a genetic disease, even when an organism-wide protein absence is incompatible with life. We provide additional examples where a similar disease mechanism applies.

Imipridone Anticancer Compounds Ectopically Activate the ClpP Protease and Represent a New Scaffold for Antibiotic Development

The imipridones ONC201 and ONC212 selectively kill cancer cells and have been ascribed multiple mechanisms-of-action. Genome-wide CRISPR knockout screens revealed that loss of the mitochondrial proteases CLPP and MIPEP confer strong resistance to both compounds...

Systematic genetic interaction profiles can reveal the mechanisms-of-action of bioactive compounds. The imipridone ONC201, which is currently in cancer clinical trials, has been ascribed a variety of different targets. To investigate the genetic dependencies of imipridone action, we screened a genome-wide clustered regularly interspaced short palindromic repeats (CRISPR) knockout library in the presence of either ONC201 or its more potent analog ONC212. Loss of the mitochondrial matrix protease CLPP or the mitochondrial intermediate peptidase MIPEP conferred strong resistance to both compounds. Biochemical and surrogate genetic assays showed that impridones directly activate CLPP and that MIPEP is necessary for proteolytic maturation of CLPP into a catalytically competent form. Quantitative proteomic analysis of cells treated with ONC212 revealed degradation of many mitochondrial as well as nonmitochondrial proteins. Prompted by the conservation of ClpP from bacteria to humans, we found that the imipridones also activate ClpP from Escherichia coli, Bacillus subtilis, and Staphylococcus aureus in biochemical and genetic assays. ONC212 and acyldepsipeptide-4 (ADEP4), a known activator of bacterial ClpP, caused similar proteome-wide degradation profiles in S. aureus. ONC212 suppressed the proliferation of a number of Gram-positive (S. aureus, B. subtilis, and Enterococcus faecium) and Gram-negative species (E. coli and Neisseria gonorrhoeae). Moreover, ONC212 enhanced the ability of rifampin to eradicate antibiotic-tolerant S. aureus persister cells. These results reveal the genetic dependencies of imipridone action in human cells and identify the imipridone scaffold as a new entry point for antibiotic development.

Characterization of METTL16 as a cytoplasmic RNA binding protein

mRNA modification by N6-methyladenosine (m6A) is involved in many post-transcriptional regulation processes including mRNA stability, splicing and promotion of translation. Accordingly, the recently identified mRNA methylation complex containing METTL3, METTL14, and WTAP has been the subject of intense study. However, METTL16 (METT10D) has also been identified as an RNA m6A methyltransferase that can methylate both coding and noncoding RNAs, but its biological role remains unclear. While global studies have identified many potential RNA targets of METTL16, only a handful, including the long noncoding RNA MALAT1, the snRNA U6, as well as the mRNA MAT2A have been verified and/or studied to any great extent. In this study we identified/verified METTL16 targets by immunoprecipitation of both endogenous as well as exogenous FLAG-tagged protein. Interestingly, exogenously overexpressed METTL16 differed from the endogenous protein in its relative affinity for RNA targets which prompted us to investigate METTL16's localization within the cell. Surprisingly, biochemical fractionation revealed that a majority of METTL16 protein resides in the cytoplasm of a number of cells. Furthermore, siRNA knockdown of METTL16 resulted in expression changes of a few mRNA targets suggesting that METTL16 may play a role in regulating gene expression. Thus, while METTL16 has been reported to be a nuclear protein, our findings suggest that METTL16 is also a cytoplasmic methyltransferase that may alter its RNA binding preferences depending on its cellular localization. Future studies will seek to confirm differences between cytoplasmic and nuclear RNA targets in addition to exploring the physiological role of METTL16 through long-term knockdown.

Conditional depletion of the RNA polymerase I subunit PAF53 reveals that it is essential for mitosis and enables identification of functional domains

Our knowledge of the mechanism of rDNA transcription has benefited from the combined application of genetic and biochemical techniques in yeast. Nomura's laboratory (Nogi, Y., Vu, L., and Nomura, M. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 7026-7030 and Nogi, Y., Yano, R., and Nomura, M. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 3962-3966) developed a system in yeast to identify genes essential for ribosome biogenesis. Such systems have allowed investigators to determine whether a gene was essential and to determine its function in rDNA transcription. However, there are significant differences in both the structures and components of the transcription apparatus and the patterns of regulation between mammals and yeast. Thus, there are significant deficits in our understanding of mammalian rDNA transcription. We have developed a system combining CRISPR/Cas9 and an auxin-inducible degron that enables combining a "genetics-like"approach with biochemistry to study mammalian rDNA transcription. We now show that the mammalian orthologue of yeast RPA49, PAF53, is required for rDNA transcription and mitotic growth. We have studied the domains of the protein required for activity. We have found that the C-terminal, DNA-binding domain (tandem-winged helix), the heterodimerization, and the linker domain were essential. Analysis of the linker identified a putative helix-turn-helix (HTH) DNA-binding domain. This HTH constitutes a second DNA-binding domain within PAF53. The HTH of the yeast and mammalian orthologues is essential for function. In summary, we show that an auxin-dependent degron system can be used to rapidly deplete nucleolar proteins in mammalian cells, that PAF53 is necessary for rDNA transcription and cell growth, and that all three PAF53 domains are necessary for its function.

Robust cullin-RING ligase function is established by a multiplicity of poly-ubiquitylation pathways

The cullin-RING ligases (CRLs) form the major family of E3 ubiquitin ligases. The prototypic CRLs in yeast, called SCF enzymes, employ a single E2 enzyme, Cdc34, to build poly-ubiquitin chains required for degradation. In contrast, six different human E2 and E3 enzyme activities, including Cdc34 orthologs UBE2R1 and UBE2R2, appear to mediate SCF-catalyzed substrate polyubiquitylation in vitro. The combinatorial interplay of these enzymes raises questions about genetic buffering of SCFs in human cells and challenges the dogma that E3s alone determine substrate specificity. To enable the quantitative comparisons of SCF-dependent ubiquitylation reactions with physiological enzyme concentrations, mass spectrometry was employed to estimate E2 and E3 levels in cells. In combination with UBE2R1/2, the E2 UBE2D3 and the E3 ARIH1 both promoted SCF-mediated polyubiquitylation in a substrate-specific fashion. Unexpectedly, UBE2R2 alone had negligible ubiquitylation activity at physiological concentrations and the ablation of UBE2R1/2 had no effect on the stability of SCF substrates in cells. A genome-wide CRISPR screen revealed that an additional E2 enzyme, UBE2G1, buffers against the loss of UBE2R1/2. UBE2G1 had robust in vitro chain extension activity with SCF, and UBE2G1 knockdown in cells lacking UBE2R1/2 resulted in stabilization of the SCF substrates p27 and CYCLIN E as well as the CUL2-RING ligase substrate HIF1α. The results demonstrate the human SCF enzyme system is diversified by association with multiple catalytic enzyme partners.

Application of CRISPR genetic screens to investigate neurological diseases

The adoption of CRISPR-Cas9 technology for functional genetic screens has been a transformative advance. Due to its modular nature, this technology can be customized to address a myriad of questions. To date, pooled, genome-scale studies have uncovered genes responsible for survival, proliferation, drug resistance, viral susceptibility, and many other functions. The technology has even been applied to the functional interrogation of the non-coding genome. However, applications of this technology to neurological diseases remain scarce. This shortfall motivated the assembly of a review that will hopefully help researchers moving in this direction find their footing. The emphasis here will be on design considerations and concepts underlying this methodology. We will highlight groundbreaking studies in the CRISPR-Cas9 functional genetics field and discuss strengths and limitations of this technology for neurological disease applications. Finally, we will provide practical guidance on navigating the many choices that need to be made when implementing a CRISPR-Cas9 functional genetic screen for the study of neurological diseases.

Genetic Interaction-Based Biomarkers Identification for Drug Resistance and Sensitivity in Cancer Cells

Cancer cells generally harbor hundreds of alterations in the cancer genomes and act as crucial factors in the development and progression of cancer. Gene alterations in the cancer genome form genetic interactions, which affect the response of patients to drugs. We developed an algorithm that mines copy number alteration and whole-exome mutation profiles from The Cancer Genome Atlas (TCGA), as well as functional screen data generated to identify potential genetic interactions for specific cancer types. As a result, 4,529 synthetic viability (SV) interactions and 10,637 synthetic lethality (SL) interactions were detected. The pharmacogenomic datasets revealed that SV interactions induced drug resistance in cancer cells and that SL interactions mediated drug sensitivity in cancer cells. Deletions of HDAC1 and DVL1, both of which participate in the Notch signaling pathway, had an SV effect in cancer cells, and deletion of DVL1 induced resistance to HDAC1 inhibitors in cancer cells. In addition, patients with low expression of both HDAC1 and DVL1 had poor prognosis. Finally, by integrating current reported genetic interactions from other studies, the Cancer Genetic Interaction database (CGIdb) (http://www.medsysbio.org/CGIdb) was constructed, providing a convenient retrieval for genetic interactions in cancer.

Link clustering explains non-central and contextually essential genes in protein interaction networks

Recent studies have shown that many essential genes (EGs) change their essentiality across various contexts. Finding contextual EGs in pathogenic conditions may facilitate the identification of therapeutic targets. We propose link clustering as an indicator of contextual EGs that are non-central in protein-protein interaction (PPI) networks. In various human and yeast PPI networks, we found that 29–47% of EGs were better characterized by link clustering than by centrality. Importantly, non-central EGs were prone to change their essentiality across different human cell lines and between species. Compared with central EGs and non-EGs, non-central EGs had intermediate levels of expression and evolutionary conservation. In addition, non-central EGs exhibited a significant impact on communities at lower hierarchical levels, suggesting that link clustering is associated with contextual essentiality, as it depicts locally important nodes in network structures.

Unbiased Boolean analysis of public gene expression data for cell cycle gene identification

Cell proliferation is essential for the development and maintenance of all organisms and is dysregulated in cancer. Using synchronized cells progressing through the cell cycle, pioneering microarray studies defined cell cycle genes based on cyclic variation in their expression. However, the concordance of the small number of synchronized cell studies has been limited, leading to discrepancies in definition of the transcriptionally regulated set of cell cycle genes within and between species. Here we present an informatics approach based on Boolean logic to identify cell cycle genes. This approach used the vast array of publicly available gene expression data sets to query similarity to CCNB1, which encodes the cyclin subunit of the Cdk1-cyclin B complex that triggers the G2-to-M transition. In addition to highlighting conservation of cell cycle genes across large evolutionary distances, this approach identified contexts where well-studied genes known to act during the cell cycle are expressed and potentially acting in nondivision contexts. An accessible web platform enables a detailed exploration of the cell cycle gene lists generated using the Boolean logic approach. The methods employed are straightforward to extend to processes other than the cell cycle.

CRISPR-Cas in Streptococcus pyogenes

The discovery and characterization of the prokaryotic CRISPR-Cas immune system has led to a revolution in genome editing and engineering technologies. Despite the fact that most applications emerged after the discovery of the type II-A CRISPR-Cas9 system of Streptococcus pyogenes, its biological importance in this organism has received little attention. Here, we provide a comprehensive overview of the current knowledge about CRISPR-Cas systems from S. pyogenes. We discuss how the interplay between CRISPR-mediated immunity and horizontal gene transfer might have modeled the evolution of this pathogen. We review the current literature about the CRISPR-Cas systems present in S. pyogenes (types I-C and II-A), and describe their distinctive biochemical and functional features. Finally, we summarize the main biotechnological applications that have arisen from the discovery of the CRISPR-Cas9 system in S. pyogenes.

The BioGRID interaction database: 2019 update

The Biological General Repository for Interaction Datasets (BioGRID: https://thebiogrid.org) is an open access database dedicated to the curation and archival storage of protein, genetic and chemical interactions for all major model organism species and humans. As of September 2018 (build 3.4.164), BioGRID contains records for 1 598 688 biological interactions manually annotated from 55 809 publications for 71 species, as classified by an updated set of controlled vocabularies for experimental detection methods. BioGRID also houses records for >700 000 post-translational modification sites. BioGRID now captures chemical interaction data, including chemical-protein interactions for human drug targets drawn from the DrugBank database and manually curated bioactive compounds reported in the literature. A new dedicated aspect of BioGRID annotates genome-wide CRISPR/Cas9-based screens that report gene-phenotype and gene-gene relationships. An extension of the BioGRID resource called the Open Repository for CRISPR Screens (ORCS) database (https://orcs.thebiogrid.org) currently contains over 500 genome-wide screens carried out in human or mouse cell lines. All data in BioGRID is made freely available without restriction, is directly downloadable in standard formats and can be readily incorporated into existing applications via our web service platforms. BioGRID data are also freely distributed through partner model organism databases and meta-databases.

Genome-wide CRISPR-dCas9 screens in E. coli identify essential genes and phage host factors

High-throughput genetic screens are powerful methods to identify genes linked to a given phenotype. The catalytic null mutant of the Cas9 RNA-guided nuclease (dCas9) can be conveniently used to silence genes of interest in a method also known as CRISPRi. Here, we report a genome-wide CRISPR-dCas9 screen using a starting pool of ~ 92,000 sgRNAs which target random positions in the chromosome of E. coli. To benchmark our method, we first investigate its utility to predict gene essentiality in the genome of E. coli during growth in rich medium. We could identify 79% of the genes previously reported as essential and demonstrate the non-essentiality of some genes annotated as essential. In addition, we took advantage of the intermediate repression levels obtained when targeting the template strand of genes to show that cells are very sensitive to the expression level of a limited set of essential genes. Our data can be visualized on CRISPRbrowser, a custom web interface available at crispr.pasteur.fr. We then apply the screen to discover E. coli genes required by phages λ, T4 and 186 to kill their host, highlighting the involvement of diverse host pathways in the infection process of the three tested phages. We also identify colanic acid capsule synthesis as a shared resistance mechanism to all three phages. Finally, using a plasmid packaging system and a transduction assay, we identify genes required for the formation of functional λ capsids, thus covering the entire phage cycle. This study demonstrates the usefulness and convenience of pooled genome-wide CRISPR-dCas9 screens in bacteria and paves the way for their broader use as a powerful tool in bacterial genomics.

Over the past few years, CRISPR-Cas technologies have emerged as powerful tools to edit genomes and modulate gene expression. They have been applied to perform high-throughput genetic screens with the purpose to understand the function of genes in a systematic manner, but the application of these screens to bacteria have so far remained limited. Here, we present the use of a library of ~92,000 guide RNAs directing the dCas9 protein to silence one by one all the genes in the chromosome of E. coli. To benchmark our method, we first investigate the performance of the technique to identify essential genes, highlighting several non-essential genes also found to be essential by other methods. We then apply our method to detect bacterial genes required by three different bacteriophages to kill E. coli and for the production of functional progeny by phage λ. Our screens highlight previously known and new genetic interactions between phages and their host’s pathways and emphasize the importance of bacterial capsule in the resistance to multiple phages. Altogether, our results demonstrate the usefulness of genome-wide CRISPR-dCas9 screens in bacteria to uncover genes involved in various phenotypes.

CRISPRO: identification of functional protein coding sequences based on genome editing dense mutagenesis

CRISPR/Cas9 pooled screening permits parallel evaluation of comprehensive guide RNA libraries to systematically perturb protein coding sequences in situ and correlate with functional readouts. For the analysis and visualization of the resulting datasets, we develop CRISPRO, a computational pipeline that maps functional scores associated with guide RNAs to genomes, transcripts, and protein coordinates and structures. No currently available tool has similar functionality. The ensuing genotype-phenotype linear and three-dimensional maps raise hypotheses about structure-function relationships at discrete protein regions. Machine learning based on CRISPRO features improves prediction of guide RNA efficacy. The CRISPRO tool is freely available at gitlab.com/bauerlab/crispro .

Integrating genetic and protein-protein interaction networks maps a functional wiring diagram of a cell

Systematic experimental approaches have led to construction of comprehensive genetic and protein-protein interaction networks for the budding yeast, Saccharomyces cerevisiae. Genetic interactions capture functional relationships between genes using phenotypic readouts, while protein-protein interactions identify physical connections between gene products. These complementary, and largely non-overlapping, networks provide a global view of the functional architecture of a cell, revealing general organizing principles, many of which appear to be evolutionarily conserved. Here, we focus on insights derived from the integration of large-scale genetic and protein-protein interaction networks, highlighting principles that apply to both unicellular and more complex systems, including human cells. Network integration reveals fundamental connections involving key functional modules of eukaryotic cells, defining a core network of cellular function, which could be elaborated to explore cell-type specificity in metazoans.

Metabolic Engineering of Escherichia coli for Enhanced Production of Naringenin 7-Sulfate and Its Biological Activities

Flavonoids are one of the predominant groups of plant polyphenols, and these compounds have significant effects on human health and nutrition. Sulfated flavonoids have more favorable attributes compared to their parent compounds such as increased solubility, stability, and bioavailability. In this research, we developed a microbial system to produce sulfated naringenin using Escherichia coli expressing a sulfotransferase (ST) from Arabidopsis thaliana (At2g03770). This wild-type strain was used as a model system for testing clustered regularly interspaced short palindromic repeats (CRISPR) interference (CRISPRi) metabolic engineering strategies. Using synthetic sgRNA to mediate transcriptional repression of cysH, a gene encoding 3'-phosphoadenosine-5'-phosphosulfate (PAPS) ST, which is involved in sulfur metabolism, resulted in an increase in intracellular PAPS accumulation by over 3.28-fold without impairing cell growth. Moreover, naringenin 7-sulfate production by engineering E. coli with its cysH gene repressed in the open reading frame through CRISPRi was enhanced by 2.83-fold in compared with the wild-type control. To improve the efficiency of biotransformation, the concentration of SO42- , glucose, and substrate were optimized. The bioproductivity of naringenin 7-sulfate was 135.49 μM [∼143.1 mg (47.7 mg L^-1)] in a 3-L fermenter at 36 h. These results demonstrated that the CRISPRi system was successfully applied for the first time in E. coli to develop an efficient microbial strain for production of a sulfated flavonoid. In addition, antibacterial and anticancer activities of naringenin 7-sulfate were investigated and found to be higher than the parent compound.