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

Aurora A Kinase Inhibition Is Synthetic Lethal with Loss of the RB1 Tumor Suppressor Gene

Loss-of-function mutations in the retinoblastoma gene RB1 are common in several treatment-refractory cancers such as small-cell lung cancer and triple-negative breast cancer. To identify drugs synthetic lethal with RB1 mutation (RB1 ^mut), we tested 36 cell-cycle inhibitors using a cancer cell panel profiling approach optimized to discern cytotoxic from cytostatic effects. Inhibitors of the Aurora kinases AURKA and AURKB showed the strongest RB1 association in this assay. LY3295668, an AURKA inhibitor with over 1,000-fold selectivity versus AURKB, is distinguished by minimal toxicity to bone marrow cells at concentrations active against RB1 ^mut cancer cells and leads to durable regression of RB1 ^mut tumor xenografts at exposures that are well tolerated in rodents. Genetic suppression screens identified enforcers of the spindle-assembly checkpoint (SAC) as essential for LY3295668 cytotoxicity in RB1-deficient cancers and suggest a model in which a primed SAC creates a unique dependency on AURKA for mitotic exit and survival. SIGNIFICANCE: The identification of a synthetic lethal interaction between RB1 and AURKA inhibition, and the discovery of a drug that can be dosed continuously to achieve uninterrupted inhibition of AURKA kinase activity without myelosuppression, suggest a new approach for the treatment of RB1-deficient malignancies, including patients progressing on CDK4/6 inhibitors.See related commentary by Dick and Li, p. 169.This article is highlighted in the In This Issue feature, p. 151.

Simultaneous zygotic inactivation of multiple genes in mouse through CRISPR/Cas9-mediated base editing

In vivo genetic mutation has become a powerful tool for dissecting gene function; however, multi-gene interaction and the compensatory mechanisms involved can make findings from single mutations, at best difficult to interpret, and, at worst, misleading. Hence, it is necessary to establish an efficient way to disrupt multiple genes simultaneously. CRISPR/Cas9-mediated base editing disrupts gene function by converting a protein-coding sequence into a stop codon; this is referred to as CRISPR-stop. Its application in generating zygotic mutations has not been well explored yet. Here, we first performed a proof-of-principle test by disrupting Atoh1, a gene crucial for auditory hair cell generation. Next, we individually mutated vGlut3 (Slc17a8), otoferlin (Otof) and prestin (Slc26a5), three genes needed for normal hearing function. Finally, we successfully disrupted vGlut3, Otof and prestin simultaneously. Our results show that CRISPR-stop can efficiently generate single or triple homozygous F0 mouse mutants, bypassing laborious mouse breeding. We believe that CRISPR-stop is a powerful method that will pave the way for high-throughput screening of mouse developmental and functional genes, matching the efficiency of methods available for model organisms such as Drosophila.

CRISPR knockout screening identifies combinatorial drug targets in pancreatic cancer and models cellular drug response

Predicting the response and identifying additional targets that will improve the efficacy of chemotherapy is a major goal in cancer research. Through large-scale in vivo and in vitro CRISPR knockout screens in pancreatic ductal adenocarcinoma cells, we identified genes whose genetic deletion or pharmacologic inhibition synergistically increase the cytotoxicity of MEK signaling inhibitors. Furthermore, we show that CRISPR viability scores combined with basal gene expression levels could model global cellular responses to the drug treatment. We develop drug response evaluation by in vivo CRISPR screening (DREBIC) method and validated its efficacy using large-scale experimental data from independent experiments. Comparative analyses demonstrate that DREBIC predicts drug response in cancer cells from a wide range of tissues with high accuracy and identifies therapeutic vulnerabilities of cancer-causing mutations to MEK inhibitors in various cancer types.

Predicting the response to chemotherapy is a major goal of cancer research. Here the authors use CRISPR knockout screens in pancreatic ductal adenocarcinoma cells to identify deletions synergistic with MEK inhibitors.

CRISPhieRmix: a hierarchical mixture model for CRISPR pooled screens

Pooled CRISPR screens allow researchers to interrogate genetic causes of complex phenotypes at the genome-wide scale and promise higher specificity and sensitivity compared to competing technologies. Unfortunately, two problems exist, particularly for CRISPRi/a screens: variability in guide efficiency and large rare off-target effects. We present a method, CRISPhieRmix, that resolves these issues by using a hierarchical mixture model with a broad-tailed null distribution. We show that CRISPhieRmix allows for more accurate and powerful inferences in large-scale pooled CRISPRi/a screens. We discuss key issues in the analysis and design of screens, particularly the number of guides needed for faithful full discovery.

Recent advances in combinatorial drug screening and synergy scoring

Treatment of complex diseases such as cancer, cardiovascular disease, diabetes or neurological disorders frequently warrants the utilization of drug combinations for therapeutic intervention. In fact, the most successful example is the current standard of care for HIV patients. However, identification of successful drug cocktails is not a simple task and is hampered by lack of standardization in terminology, experimental protocols and models as well as data analysis. Here we discuss the most recent developments in combinatorial drug screening by covering technological advancements in screening strategies, cellular model systems as well as novel drug classes. We believe the research progress being made provides promising basis to build on and identify, develop and optimize efficacious clinically relevant combinatorial drug treatments.

Fluorescence Activated Cell Sorting (FACS) in Genome-Wide Genetic Screening of Membrane Trafficking

About one-third of cellular proteins in eukaryotic cells are localized to membrane-enclosed organelles in the endomembrane system. Trafficking of these membrane proteins (including soluble lumenal proteins) among the organelles is mediated by small sac-like vesicles. Vesicle-mediated membrane trafficking regulates a broad range of biological processes, many of which are still poorly understood at the molecular level. A powerful approach to dissect a vesicle-mediated membrane trafficking pathway is unbiased genome-wide genetic screening, which only recently became possible in mammalian cells with the isolation of haploid human cell lines and the development of CRISPR-Cas9 genome editing. Here, we describe a FACS-based method to select populations of live mutant cells based on the surface levels of endogenous proteins or engineered reporters. Collection of these mutant populations enables subsequent deep sequencing and bioinformatics analysis to identify genes that regulate the trafficking pathway. This method can be readily adapted to genetically dissect a broad range of mammalian membrane trafficking processes using haploid genetics or CRISPR-Cas9 screens. © 2018 by John Wiley & Sons, Inc.

CRISPR/Cas9 System: A Bacterial Tailor for Genomic Engineering

Microbes use diverse defence strategies that allow them to withstand exposure to a variety of genome invaders such as bacteriophages and plasmids. One such defence strategy is the use of RNA guided endonuclease called CRISPR-associated (Cas) 9 protein. The Cas9 protein, derived from type II CRISPR/Cas system, has been adapted as a versatile tool for genome targeting and engineering due to its simplicity and high efficiency over the earlier tools such as ZFNs and TALENs. With recent advancements, CRISPR/Cas9 technology has emerged as a revolutionary tool for modulating the genome in living cells and inspires innovative translational applications in different fields. In this paper we review the developments and its potential uses in the CRISPR/Cas9 technology as well as recent advancements in genome engineering using CRISPR/Cas9.

CRISPR/Cas9 cleavage efficiency regression through boosting algorithms and Markov sequence profiling

Motivation

CRISPR/Cas9 system is a widely used genome editing tool. A prediction problem of great interests for this system is: how to select optimal single-guide RNAs (sgRNAs), such that its cleavage efficiency is high meanwhile the off-target effect is low.

Results

This work proposed a two-step averaging method (TSAM) for the regression of cleavage efficiencies of a set of sgRNAs by averaging the predicted efficiency scores of a boosting algorithm and those by a support vector machine (SVM). We also proposed to use profiled Markov properties as novel features to capture the global characteristics of sgRNAs. These new features are combined with the outstanding features ranked by the boosting algorithm for the training of the SVM regressor. TSAM improved the mean Spearman correlation coefficiencies comparing with the state-of-the-art performance on benchmark datasets containing thousands of human, mouse and zebrafish sgRNAs. Our method can be also converted to make binary distinctions between efficient and inefficient sgRNAs with superior performance to the existing methods. The analysis reveals that highly efficient sgRNAs have lower melting temperature at the middle of the spacer, cut at 5'-end closer parts of the genome and contain more 'A' but less 'G' comparing with inefficient ones. Comprehensive further analysis also demonstrates that our tool can predict an sgRNA's cutting efficiency with consistently good performance no matter it is expressed from an U6 promoter in cells or from a T7 promoter in vitro.

Availability and implementation

Online tool is available at http://www.aai-bioinfo.com/CRISPR/. Python and Matlab source codes are freely available at https://github.com/penn-hui/TSAM.

Supplementary information

Supplementary data are available at Bioinformatics online.

Synthetic Lethal Networks for Precision Oncology: Promises and Pitfalls

Synthetic lethal interactions, in which the simultaneous loss of function of two genes produces a lethal phenotype, are being explored as a means to therapeutically exploit cancer-specific vulnerabilities and expand the scope of precision oncology. Currently, three Food and Drug Administration-approved drugs work by targeting the synthetic lethal interaction between BRCA1/2 and PARP. This review examines additional efforts to discover networks of synthetic lethal interactions and discusses both challenges and opportunities regarding the translation of new synthetic lethal interactions into the clinic.

Reduced intrinsic DNA curvature leads to increased mutation rate

BACKGROUND

Mutation rates vary across the genome. Many trans factors that influence mutation rates have been identified, as have specific sequence motifs at the 1-7-bp scale, but cis elements remain poorly characterized. The lack of understanding regarding why different sequences have different mutation rates hampers our ability to identify positive selection in evolution and to identify driver mutations in tumorigenesis.

RESULTS

Here, we use a combination of synthetic genes and sequences of thousands of isolated yeast colonies to show that intrinsic DNA curvature is a major cis determinant of mutation rate. Mutation rate negatively correlates with DNA curvature within genes, and a 10% decrease in curvature results in a 70% increase in mutation rate. Consistently, both yeast and humans accumulate mutations in regions with small curvature. We further show that this effect is due to differences in the intrinsic mutation rate, likely due to differences in mutagen sensitivity and not due to differences in the local activity of DNA repair.

CONCLUSIONS

Our study establishes a framework for understanding the cis properties of DNA sequence in modulating the local mutation rate and identifies a novel causal source of non-uniform mutation rates across the genome.

Interrogating the protein interactomes of RAS isoforms identifies PIP5K1A as a KRAS-specific vulnerability

In human cancers, oncogenic mutations commonly occur in the RAS genes KRAS, NRAS, or HRAS, but there are no clinical RAS inhibitors. Mutations are more prevalent in KRAS, possibly suggesting a unique oncogenic activity mediated by KRAS-specific interaction partners, which might be targeted. Here, we determine the specific protein interactomes of each RAS isoform by BirA proximity-dependent biotin identification. The combined interactomes are screened by CRISPR-Cas9 loss-of-function assays for proteins required for oncogenic KRAS-dependent, NRAS-dependent, or HRAS-dependent proliferation and censored for druggable proteins. Using this strategy, we identify phosphatidylinositol phosphate kinase PIP5K1A as a KRAS-specific interactor and show that PIP5K1A binds to a unique region in KRAS. Furthermore, PIP5K1A depletion specifically reduces oncogenic KRAS signaling and proliferation, and sensitizes pancreatic cancer cell lines to a MAPK inhibitor. These results suggest PIP5K1A as a potential target in KRAS signaling for the treatment of KRAS-mutant cancers.

RAS isoforms are frequently mutated in cancer but their inhibition remains challenging. By comparing the protein interactomes of the highly similar isoforms HRAS, NRAS and KRAS, the authors here identify PIP5K1A as a KRAS-specific interactor and a target to inhibit KRAS-driven cell growth.

Altered patterns of global protein synthesis and translational fidelity in RPS15-mutated chronic lymphocytic leukemia

Genomic studies have recently identified RPS15 as a new driver gene in aggressive and chemorefractory cases of chronic lymphocytic leukemia (CLL). RPS15 encodes a ribosomal protein whose conserved C-terminal domain extends into the decoding center of the ribosome. We demonstrate that mutations in highly conserved residues of this domain affect protein stability, by increasing its ubiquitin-mediated degradation, and cell-proliferation rates. On the other hand, we show that mutated RPS15 can be loaded into the ribosomes, directly impacting on global protein synthesis and/or translational fidelity in a mutation-specific manner. Quantitative mass spectrometry analyses suggest that RPS15 variants may induce additional alterations in the translational machinery, as well as a metabolic shift at the proteome level in HEK293T and MEC-1 cells. These results indicate that CLL-related RPS15 mutations might act following patterns known for other ribosomal diseases, likely switching from a hypo- to a hyperproliferative phenotype driven by mutated ribosomes. In this scenario, loss of translational fidelity causing altered cell proteostasis can be proposed as a new molecular mechanism involved in CLL pathobiology.

Large-scale investigation of the reasons why potentially important genes are ignored

Biomedical research has been previously reported to primarily focus on a minority of all known genes. Here, we demonstrate that these differences in attention can be explained, to a large extent, exclusively from a small set of identifiable chemical, physical, and biological properties of genes. Together with knowledge about homologous genes from model organisms, these features allow us to accurately predict the number of publications on individual human genes, the year of their first report, the levels of funding awarded by the National Institutes of Health (NIH), and the development of drugs against disease-associated genes. By explicitly identifying the reasons for gene-specific bias and performing a meta-analysis of existing computational and experimental knowledge bases, we describe gene-specific strategies for the identification of important but hitherto ignored genes that can open novel directions for future investigation.

PPP1R35 is a novel centrosomal protein that regulates centriole length in concert with the microcephaly protein RTTN

Centrosome structure, function, and number are finely regulated at the cellular level to ensure normal mammalian development. Here, we characterize PPP1R35 as a novel bona fide centrosomal protein and demonstrate that it is critical for centriole elongation. Using quantitative super-resolution microscopy mapping and live-cell imaging we show that PPP1R35 is a resident centrosomal protein located in the proximal lumen above the cartwheel, a region of the centriole that has eluded detailed characterization. Loss of PPP1R35 function results in decreased centrosome number and shortened centrioles that lack centriolar distal and microtubule wall associated proteins required for centriole elongation. We further demonstrate that PPP1R35 acts downstream of, and forms a complex with, RTTN, a microcephaly protein required for distal centriole elongation. Altogether, our study identifies a novel step in the centriole elongation pathway centered on PPP1R35 and elucidates downstream partners of the microcephaly protein RTTN.

RNAi/CRISPR Screens: from a Pool to a Valid Hit

High-throughput genetic screens interfering with gene expression are invaluable tools to identify gene function and phenotype-to-genotype interactions. Implementing such screens in the laboratory is challenging, and the choice between currently available technologies based on RNAi and CRISPR/Cas9 (CRISPR-associated protein 9) is not trivial. Identifying reliable candidate hits requires a streamlined experimental setup adjusted to the specific biological question. Here, we provide a critical assessment of the various RNAi/CRISPR approaches to pooled screens and discuss their advantages and pitfalls. We specify a set of best practices for key parameters enabling a reproducible screen and provide a detailed overview of analysis methods and repositories for identifying the best candidate gene hits.

SHLD2/FAM35A co-operates with REV7 to coordinate DNA double-strand break repair pathway choice

DNA double-strand breaks (DSBs) can be repaired by two major pathways: non-homologous end-joining (NHEJ) and homologous recombination (HR). DNA repair pathway choice is governed by the opposing activities of 53BP1, in complex with its effectors RIF1 and REV7, and BRCA1. However, it remains unknown how the 53BP1/RIF1/REV7 complex stimulates NHEJ and restricts HR to the S/G2 phases of the cell cycle. Using a mass spectrometry (MS)-based approach, we identify 11 high-confidence REV7 interactors and elucidate the role of SHLD2 (previously annotated as FAM35A and RINN2) as an effector of REV7 in the NHEJ pathway. FAM35A depletion impairs NHEJ-mediated DNA repair and compromises antibody diversification by class switch recombination (CSR) in B cells. FAM35A accumulates at DSBs in a 53BP1-, RIF1-, and REV7-dependent manner and antagonizes HR by limiting DNA end resection. In fact, FAM35A is part of a larger complex composed of REV7 and SHLD1 (previously annotated as C20orf196 and RINN3), which promotes NHEJ and limits HR Together, these results establish SHLD2 as a novel effector of REV7 in controlling the decision-making process during DSB repair.

[Searching for essential genes in cancer genomes]

The concept of essential genes, whose loss of functionality leads to cell death, is one of the fundamental concepts of genetics and is important for fundamental and applied research. This field is particularly promising in relation to oncology, since the search for genetic vulnerabilities of cancer cells allows us to identify new potential targets for antitumor therapy. The modern biotechnology capacities allow carrying out large-scale projects for sequencing somatic mutations in tumors, as well as directly interfering the genetic apparatus of cancer cells. They provided accumulation of a considerable body of knowledge about genetic variants and corresponding phenotypic manifestations in tumors. In the near future this knowledge will find application in clinical practice. This review describes the main experimental and computational approaches to the search for essential genes, concentrating on the application of these methods in the field of molecular oncology.

Rscreenorm: normalization of CRISPR and siRNA screen data for more reproducible hit selection

Background

Reproducibility of hits from independent CRISPR or siRNA screens is poor. This is partly due to data normalization primarily addressing technical variability within independent screens, and not the technical differences between them.

Results

We present “rscreenorm”, a method that standardizes the functional data ranges between screens using assay controls, and subsequently performs a piecewise-linear normalization to make data distributions across all screens comparable. In simulation studies, rscreenorm reduces false positives. Using two multiple-cell lines siRNA screens, rscreenorm increased reproducibility between 27 and 62% for hits, and up to 5-fold for non-hits. Using publicly available CRISPR-Cas screen data, application of commonly used median centering yields merely 34% of overlapping hits, in contrast with rscreenorm yielding 84% of overlapping hits. Furthermore, rscreenorm yielded at most 8% discordant results, whilst median-centering yielded as much as 55%.

Conclusions

Rscreenorm yields more consistent results and keeps false positive rates under control, improving reproducibility of genetic screens data analysis from multiple cell lines.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2306-z) contains supplementary material, which is available to authorized users.

Gene essentiality landscape and druggable oncogenic dependencies in herpesviral primary effusion lymphoma

Primary effusion lymphoma (PEL) is caused by Kaposi's sarcoma-associated herpesvirus. Our understanding of PEL is poor and therefore treatment strategies are lacking. To address this need, we conducted genome-wide CRISPR/Cas9 knockout screens in eight PEL cell lines. Integration with data from unrelated cancers identifies 210 genes as PEL-specific oncogenic dependencies. Genetic requirements of PEL cell lines are largely independent of Epstein-Barr virus co-infection. Genes of the NF-κB pathway are individually non-essential. Instead, we demonstrate requirements for IRF4 and MDM2. PEL cell lines depend on cellular cyclin D2 and c-FLIP despite expression of viral homologs. Moreover, PEL cell lines are addicted to high levels of MCL1 expression, which are also evident in PEL tumors. Strong dependencies on cyclin D2 and MCL1 render PEL cell lines highly sensitive to palbociclib and S63845. In summary, this work comprehensively identifies genetic dependencies in PEL cell lines and identifies novel strategies for therapeutic intervention.

Targeting the HTLV-I-Regulated BATF3/IRF4 Transcriptional Network in Adult T Cell Leukemia/Lymphoma

Adult T cell leukemia/lymphoma (ATLL) is a frequently incurable disease associated with the human lymphotropic virus type I (HTLV-I). RNAi screening of ATLL lines revealed that their proliferation depends on BATF3 and IRF4, which cooperatively drive ATLL-specific gene expression. HBZ, the only HTLV-I encoded transcription factor that is expressed in all ATLL cases, binds to an ATLL-specific BATF3 super-enhancer and thereby regulates the expression of BATF3 and its downstream targets, including MYC. Inhibitors of bromodomain-and-extra-terminal-domain (BET) chromatin proteins collapsed the transcriptional network directed by HBZ and BATF3, and were consequently toxic for ATLL cell lines, patient samples, and xenografts. Our study demonstrates that the HTLV-I oncogenic retrovirus exploits a regulatory module that can be attacked therapeutically with BET inhibitors.

Unsupervised correction of gene-independent cell responses to CRISPR-Cas9 targeting

BACKGROUND

Genome editing by CRISPR-Cas9 technology allows large-scale screening of gene essentiality in cancer. A confounding factor when interpreting CRISPR-Cas9 screens is the high false-positive rate in detecting essential genes within copy number amplified regions of the genome. We have developed the computational tool CRISPRcleanR which is capable of identifying and correcting gene-independent responses to CRISPR-Cas9 targeting. CRISPRcleanR uses an unsupervised approach based on the segmentation of single-guide RNA fold change values across the genome, without making any assumption about the copy number status of the targeted genes.

RESULTS

Applying our method to existing and newly generated genome-wide essentiality profiles from 15 cancer cell lines, we demonstrate that CRISPRcleanR reduces false positives when calling essential genes, correcting biases within and outside of amplified regions, while maintaining true positive rates. Established cancer dependencies and essentiality signals of amplified cancer driver genes are detectable post-correction. CRISPRcleanR reports sgRNA fold changes and normalised read counts, is therefore compatible with downstream analysis tools, and works with multiple sgRNA libraries.

CONCLUSIONS

CRISPRcleanR is a versatile open-source tool for the analysis of CRISPR-Cas9 knockout screens to identify essential genes.

Multiplexed assays of variant effects contribute to a growing genotype-phenotype atlas

Given the constantly improving cost and speed of genome sequencing, it is reasonable to expect that personal genomes will soon be known for many millions of humans. This stands in stark contrast with our limited ability to interpret the sequence variants which we find. Although it is, perhaps, easiest to interpret variants in coding regions, knowledge of functional impact is unknown for the vast majority of missense variants. While many computational approaches can predict the impact of coding variants, they are given a little weight in the current guidelines for interpreting clinical variants. Laboratory assays produce comparatively more trustworthy results, but until recently did not scale to the space of all possible mutations. The development of deep mutational scanning and other multiplexed assays of variant effect has now brought feasibility of this endeavour within view. Here, we review progress in this field over the last decade, break down the different approaches into their components, and compare methodological differences.

The shieldin complex mediates 53BP1-dependent DNA repair

53BP1 is a chromatin-binding protein that regulates the repair of DNA double-strand breaks by suppressing the nucleolytic resection of DNA termini1,2. This function of 53BP1 requires interactions with PTIP3 and RIF14-9, the latter of which recruits REV7 (also known as MAD2L2) to break sites10,11. How 53BP1-pathway proteins shield DNA ends is currently unknown, but there are two models that provide the best potential explanation of their action. In one model the 53BP1 complex strengthens the nucleosomal barrier to end-resection nucleases12,13, and in the other 53BP1 recruits effector proteins with end-protection activity. Here we identify a 53BP1 effector complex, shieldin, that includes C20orf196 (also known as SHLD1), FAM35A (SHLD2), CTC-534A2.2 (SHLD3) and REV7. Shieldin localizes to double-strand-break sites in a 53BP1- and RIF1-dependent manner, and its SHLD2 subunit binds to single-stranded DNA via OB-fold domains that are analogous to those of RPA1 and POT1. Loss of shieldin impairs non-homologous end-joining, leads to defective immunoglobulin class switching and causes hyper-resection. Mutations in genes that encode shieldin subunits also cause resistance to poly(ADP-ribose) polymerase inhibition in BRCA1-deficient cells and tumours, owing to restoration of homologous recombination. Finally, we show that binding of single-stranded DNA by SHLD2 is critical for shieldin function, consistent with a model in which shieldin protects DNA ends to mediate 53BP1-dependent DNA repair.

Synthetic Lethal Vulnerabilities in KRAS-Mutant Cancers

KRAS is the most commonly mutated oncogene in human cancer. Most KRAS-mutant cancers depend on sustained expression and signaling of KRAS, thus making it a high-priority therapeutic target. Unfortunately, development of direct small molecule inhibitors of KRAS function has been challenging. An alternative therapeutic strategy for KRAS-mutant malignancies involves targeting codependent vulnerabilities or synthetic lethal partners that are preferentially essential in the setting of oncogenic KRAS. KRAS activates numerous effector pathways that mediate proliferation and survival signals. Moreover, cancer cells must cope with substantial oncogenic stress conferred by mutant KRAS. These oncogenic signaling pathways and compensatory coping mechanisms of KRAS-mutant cancer cells form the basis for synthetic lethal interactions. Here, we review the compendium of previously identified codependencies in KRAS-mutant cancers, including the results of numerous functional genetic screens aimed at identifying KRAS synthetic lethal targets. Importantly, many of these vulnerabilities may represent tractable therapeutic opportunities.

Tackling resistance: emerging antimalarials and new parasite targets in the era of elimination

Malaria remains a significant contributor to global human mortality, and roughly half the world’s population is at risk for infection with Plasmodium spp. parasites. Aggressive control measures have reduced the global prevalence of malaria significantly over the past decade. However, resistance to available antimalarials continues to spread, including resistance to the widely used artemisinin-based combination therapies. Novel antimalarial compounds and therapeutic targets are greatly needed. This review will briefly discuss several promising current antimalarial development projects, including artefenomel, ferroquine, cipargamin, SJ733, KAF156, MMV048, and tafenoquine. In addition, we describe recent large-scale genetic and resistance screens that have been instrumental in target discovery. Finally, we highlight new antimalarial targets, which include essential transporters and proteases. These emerging antimalarial compounds and therapeutic targets have the potential to overcome multi-drug resistance in ongoing efforts toward malaria elimination.

Estrogen-regulated feedback loop limits the efficacy of estrogen receptor-targeted breast cancer therapy

Endocrine therapy resistance invariably develops in advanced estrogen receptor-positive (ER+) breast cancer, but the underlying mechanisms are largely unknown. We have identified C-terminal SRC kinase (CSK) as a critical node in a previously unappreciated negative feedback loop that limits the efficacy of current ER-targeted therapies. Estrogen directly drives CSK expression in ER+ breast cancer. At low CSK levels, as is the case in patients with ER+ breast cancer resistant to endocrine therapy and with the poorest outcomes, the p21 protein-activated kinase 2 (PAK2) becomes activated and drives estrogen-independent growth. PAK2 overexpression is also associated with endocrine therapy resistance and worse clinical outcome, and the combination of a PAK2 inhibitor with an ER antagonist synergistically suppressed breast tumor growth. Clinical approaches to endocrine therapy-resistant breast cancer must overcome the loss of this estrogen-induced negative feedback loop that normally constrains the growth of ER+ tumors.

Explaining cancer type specific mutations with transcriptomic and epigenomic features in normal tissues

Most cancer driver genes are involved in generic cellular processes such as DNA repair, cell proliferation and cell adhesion, yet their mutations are often confined to specific cancer types. To resolve this paradox, we explained mutation frequencies of selected genes across tumor types with four features in the corresponding normal tissues from cancer-free subjects: mRNA expression and chromatin accessibility of mutated genes, mRNA expressions of their neighbors in curated pathways and the protein-protein interaction network. Encouragingly, these transcriptomic/epigenomic features in normal tissues were closely associated with mutational/functional characteristics in tumors. First, chromatin accessibility was a necessary but not sufficient condition for frequent mutations. Second, variations of mutation frequencies in selected genes across tissue types were significantly associated with all four features. Third, the genes possessing significant associations between mutation frequency variations and pathway gene expression were enriched with documented cancer genes. We further proposed a novel bivariate gene set enrichment analysis and confirmed that the pathway gene expression was the dominant factor in cancer gene enrichment. These findings shed lights on the functional roles of genes in normal tissues in shaping the mutational landscape during tumor genome evolution.

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.

Mitochondrial double-stranded RNA triggers antiviral signalling in humans

Mitochondria are descendants of endosymbiotic bacteria and retain essential prokaryotic features such as a compact circular genome. Consequently, in mammals, mitochondrial DNA is subjected to bidirectional transcription that generates overlapping transcripts, which are capable of forming long double-stranded RNA structures^1,2. However, to our knowledge, mitochondrial double-stranded RNA has not been previously characterized in vivo. Here we describe the presence of a highly unstable native mitochondrial double-stranded RNA species at single-cell level and identify key roles for the degradosome components mitochondrial RNA helicase SUV3 and polynucleotide phosphorylase PNPase in restricting the levels of mitochondrial double-stranded RNA. Loss of either enzyme results in massive accumulation of mitochondrial double-stranded RNA that escapes into the cytoplasm in a PNPase-dependent manner. This process engages an MDA5-driven antiviral signalling pathway that triggers a type I interferon response. Consistent with these data, patients carrying hypomorphic mutations in the gene PNPT1, which encodes PNPase, display mitochondrial double-stranded RNA accumulation coupled with upregulation of interferon-stimulated genes and other markers of immune activation. The localization of PNPase to the mitochondrial inter-membrane space and matrix suggests that it has a dual role in preventing the formation and release of mitochondrial double-stranded RNA into the cytoplasm. This in turn prevents the activation of potent innate immune defence mechanisms that have evolved to protect vertebrates against microbial and viral attack.

Uncovering the essential genes of the human malaria parasite Plasmodium falciparum by saturation mutagenesis

Severe malaria is caused by the apicomplexan parasite Plasmodium falciparum. Despite decades of research, the distinct biology of these parasites has made it challenging to establish high-throughput genetic approaches to identify and prioritize therapeutic targets. Using transposon mutagenesis of P. falciparum in an approach that exploited its AT-rich genome, we generated more than 38,000 mutants, saturating the genome and defining mutability and fitness costs for over 87% of genes. Of 5399 genes, our study defined 2680 genes as essential for optimal growth of asexual blood stages in vitro. These essential genes are associated with drug resistance, represent leading vaccine candidates, and include approximately 1000 Plasmodium-conserved genes of unknown function. We validated this approach by testing proteasome pathways for individual mutants associated with artemisinin sensitivity.

Mapping the Genetic Landscape of Human Cells

Seminal yeast studies have established the value of comprehensively mapping genetic interactions (GIs) for inferring gene function. Efforts in human cells using focused gene sets underscore the utility of this approach, but the feasibility of generating large-scale, diverse human GI maps remains unresolved. We developed a CRISPR interference platform for large-scale quantitative mapping of human GIs. We systematically perturbed 222,784 gene pairs in two cancer cell lines. The resultant maps cluster functionally related genes, assigning function to poorly characterized genes, including TMEM261, a new electron transport chain component. Individual GIs pinpoint unexpected relationships between pathways, exemplified by a specific cholesterol biosynthesis intermediate whose accumulation induces deoxynucleotide depletion, causing replicative DNA damage and a synthetic-lethal interaction with the ATR/9-1-1 DNA repair pathway. Our map provides a broad resource, establishes GI maps as a high-resolution tool for dissecting gene function, and serves as a blueprint for mapping the genetic landscape of human cells.

Conditional Degrons for Controlling Protein Expression at the Protein Level

The conditional depletion of a protein of interest (POI) is useful not only for loss-of-function studies, but also for the modulation of biological pathways. Technologies that work at the level of DNA, mRNA, and protein are available for temporal protein depletion. Compared with technologies targeting the pretranslation steps, direct protein depletion (or protein knockdown approaches) is advantageous in terms of specificity, reversibility, and time required for depletion, which can be achieved by fusing a POI with a protein domain called a degron that induces rapid proteolysis of the fusion protein. Conditional degrons can be activated or inhibited by temperature, small molecules, light, or the expression of another protein. The conditional degron-based technologies currently available are described and discussed.

CRISPR screens identify genomic ribonucleotides as a source of PARP-trapping lesions

The observation that BRCA1- and BRCA2-deficient cells are sensitive to poly(ADP-ribose) polymerase (PARP) inhibitors spurred their development into cancer therapies that target homologous recombination (HR) deficiency1. The cytotoxicity of PARP inhibitors depends on PARP trapping, the formation of non-covalent protein-DNA adducts composed of inhibited PARP1 bound to DNA lesions of unclear origins1–4. To address the nature of such lesions and the cellular consequences of PARP trapping, we undertook three CRISPR screens to identify genes and pathways that mediate cellular resistance to olaparib, a clinically approved PARP inhibitor1. Here were present a high-confidence set of 73 genes whose mutation causes increased PARP inhibitor sensitivity. In addition to an expected enrichment for HR-related genes, we discovered that mutation in all three genes encoding RNase H2 sensitized cells to PARP inhibition. We establish that the underlying cause of the PARP inhibitor hypersensitivity of RNase H2-deficient cells is impaired ribonucleotide excision repair (RER)5. Embedded ribonucleotides, abundant in the genome of RER-deficient cells, are substrates for topoisomerase 1 cleavage, resulting in PARP-trapping lesions that impede DNA replication and endanger genome integrity. We conclude that genomic ribonucleotides are a hitherto unappreciated source of PARP-trapping DNA lesions, and that the frequent deletion of RNASEH2B in metastatic prostate cancer and chronic lymphocytic leukemia could provide an opportunity to exploit these findings therapeutically.

Correction of copy number induced false positives in CRISPR screens

Cell autonomous cancer dependencies are now routinely identified using CRISPR loss-of-function viability screens. However, a bias exists that makes it difficult to assess the true essentiality of genes located in amplicons, since the entire amplified region can exhibit lethal scores. These false-positive hits can either be discarded from further analysis, which in cancer models can represent a significant number of hits, or methods can be developed to rescue the true-positives within amplified regions. We propose two methods to rescue true positive hits in amplified regions by correcting for this copy number artefact. The Local Drop Out (LDO) method uses the relative lethality scores within genomic regions to assess true essentiality and does not require additional orthogonal data (e.g. copy number value). LDO is meant to be used in screens covering a dense region of the genome (e.g. a whole chromosome or the whole genome). The General Additive Model (GAM) method models the screening data as a function of the known copy number values and removes the systematic effect from the measured lethality. GAM does not require the same density as LDO, but does require prior knowledge of the copy number values. Both methods have been developed with single sample experiments in mind so that the correction can be applied even in smaller screens. Here we demonstrate the efficacy of both methods at removing the copy number effect and rescuing hits from some of the amplified regions. We estimate a 70-80% decrease of false positive hits with either method in regions of high copy number compared to no correction.

p53 inhibits CRISPR-Cas9 engineering in human pluripotent stem cells

CRISPR/Cas9 has revolutionized our ability to engineer genomes and conduct genome-wide screens in human cells1-3. Whereas some cell types are amenable to genome engineering, genomes of human pluripotent stem cells (hPSCs) have been difficult to engineer, with reduced efficiencies relative to tumour cell lines or mouse embryonic stem cells3-13. Here, using hPSC lines with stable integration of Cas9 or transient delivery of Cas9-ribonucleoproteins (RNPs), we achieved an average insertion or deletion (indel) efficiency greater than 80%. This high efficiency of indel generation revealed that double-strand breaks (DSBs) induced by Cas9 are toxic and kill most hPSCs. In previous studies, the toxicity of Cas9 in hPSCs was less apparent because of low transfection efficiency and subsequently low DSB induction3. The toxic response to DSBs was P53/TP53-dependent, such that the efficiency of precise genome engineering in hPSCs with a wild-type P53 gene was severely reduced. Our results indicate that Cas9 toxicity creates an obstacle to the high-throughput use of CRISPR/Cas9 for genome engineering and screening in hPSCs. Moreover, as hPSCs can acquire P53 mutations14, cell replacement therapies using CRISPR/Cas9-enginereed hPSCs should proceed with caution, and such engineered hPSCs should be monitored for P53 function.

MicroRNA-focused CRISPR-Cas9 library screen reveals fitness-associated miRNAs

MicroRNAs (miRNAs) are post-transcriptional gene regulators that play important roles in the control of cell fitness, differentiation, and development. The CRISPR-Cas9 gene-editing system is composed of the Cas9 nuclease in complex with a single guide RNA (sgRNA) and directs DNA cleavage at a predetermined site. Several CRISPR-Cas9 libraries have been constructed for genome-scale knockout screens of protein function; however, few libraries have included miRNA genes. Here we constructed a miRNA-focused CRISPR-Cas9 library that targets 1594 (85%) annotated human miRNA stem-loops. The sgRNAs in our LX-miR library are designed to have high on-target and low off-target activity, and each miRNA is targeted by four to five sgRNAs. We used this sgRNA library to screen for miRNAs that affect cell fitness of HeLa or NCI-N87 cells by monitoring the change in frequency of each sgRNA over time. By considering the expression in the tested cells and the dysregulation of the miRNAs in cancer specimens, we identified five HeLa pro-fitness and cervical cancer up-regulated miRNAs (miR-31-5p, miR-92b-3p, miR-146b-5p, miR-151a-3p, and miR-194-5p). Similarly, we identified six NCI-N87 pro-fitness and gastric cancer up-regulated miRNAs (miR-95-3p, miR-181a-5p, miR-188-5p, miR-196b-5p, miR-584-5p, and miR-1304-3p), as well as three anti-fitness and down-regulated miRNAs (let-7a-3p, miR-100-5p, and miR-149-5p). Some of those miRNAs are known to be oncogenic or tumor-suppressive, but others are novel. Taken together, the LX-miR library is useful for genome-wide unbiased screening to identify miRNAs important for cellular fitness and likely to be useful for other functional screens.

CK1α and IRF4 are essential and independent effectors of immunomodulatory drugs in primary effusion lymphoma

Primary effusion lymphoma (PEL) is an aggressive cancer with few treatment options. The immunomodulatory drugs (IMiDs) lenalidomide and pomalidomide have recently been shown to kill PEL cell lines, and lenalidomide is in clinical trials against PEL. IMiDs bind to the CRL4^CRBN E3 ubiquitin ligase complex, leading to the acquisition of the Ikaros family zinc finger proteins 1 and 3 (IKZF1 and IKZF3), casein kinase 1 α (CK1α), and zinc finger protein 91 (ZFP91) as neosubstrates. IMiDs are effective against multiple myeloma because of degradation of IKZF1 and IKZF3 and the consequent loss of interferon regulatory factor 4 (IRF4) and MYC expression. Lenalidomide is also effective in chromosome 5q deletion-associated myelodysplastic syndrome as a result of degradation of CK1α. An essential IKZF1-IRF4-MYC axis has recently been proposed to underlie the toxicity of IMiDs in PEL. Here, we further investigate IMiD effectors in PEL cell lines, based on genome-wide CRISPR/Cas9 screens for essential human genes. These screens and extensive validation experiments show that, of the 4 neosubstrates, only CK1α is essential for the survival of PEL cell lines. In contrast, IKZF1 and IKZF3 are dispensable, individually or in combination. IRF4 was critical in all 8 PEL cell lines tested, and surprisingly, IMiDs triggered downregulation of IRF4 expression independently of both IKZF1 and IKZF3. Reexpression of CK1α and/or IRF4 partially rescued PEL cell lines from IMiD-mediated toxicity. In conclusion, IMiD toxicity in PEL cell lines is independent of IKZF1 and IKZF3 but proceeds through degradation of the neosubstrate CK1α and downregulation of IRF4.

DeepCRISPR: optimized CRISPR guide RNA design by deep learning

A major challenge for effective application of CRISPR systems is to accurately predict the single guide RNA (sgRNA) on-target knockout efficacy and off-target profile, which would facilitate the optimized design of sgRNAs with high sensitivity and specificity. Here we present DeepCRISPR, a comprehensive computational platform to unify sgRNA on-target and off-target site prediction into one framework with deep learning, surpassing available state-of-the-art in silico tools. In addition, DeepCRISPR fully automates the identification of sequence and epigenetic features that may affect sgRNA knockout efficacy in a data-driven manner. DeepCRISPR is available at http://www.deepcrispr.net/ .

A multiprotein supercomplex controlling oncogenic signalling in lymphoma

B cell receptor (BCR) signalling has emerged as a therapeutic target in B cell lymphomas, but inhibiting this pathway in diffuse large B cell lymphoma (DLBCL) has benefited only a subset of patients¹. Gene expression profiling identified two major subtypes of DLBCL, known as germinal centre B cell-like and activated B cell-like (ABC)^2,3, that show poor outcomes after immunochemotherapy in ABC. Autoantigens drive BCR-dependent activation of NF-κB in ABC DLBCL through a kinase signalling cascade of SYK, BTK and PKCβ to promote the assembly of the CARD11-BCL10-MALT1 adaptor complex, which recruits and activates IκB kinase^4-6. Genome sequencing revealed gain-of-function mutations that target the CD79A and CD79B BCR subunits and the Toll-like receptor signalling adaptor MYD88^5,7, with MYD88(L265P) being the most prevalent isoform. In a clinical trial, the BTK inhibitor ibrutinib produced responses in 37% of cases of ABC¹. The most striking response rate (80%) was observed in tumours with both CD79B and MYD88(L265P) mutations, but how these mutations cooperate to promote dependence on BCR signalling remains unclear. Here we used genome-wide CRISPR-Cas9 screening and functional proteomics to determine the molecular basis of exceptional clinical responses to ibrutinib. We discovered a new mode of oncogenic BCR signalling in ibrutinib-responsive cell lines and biopsies, coordinated by a multiprotein supercomplex formed by MYD88, TLR9 and the BCR (hereafter termed the My-T-BCR supercomplex). The My-T-BCR supercomplex co-localizes with mTOR on endolysosomes, where it drives pro-survival NF-κB and mTOR signalling. Inhibitors of BCR and mTOR signalling cooperatively decreased the formation and function of the My-T-BCR supercomplex, providing mechanistic insight into their synergistic toxicity for My-T-BCR⁺ DLBCL cells. My-T-BCR supercomplexes characterized ibrutinib-responsive malignancies and distinguished ibrutinib responders from non-responders. Our data provide a framework for the rational design of oncogenic signalling inhibitors in molecularly defined subsets of DLBCL.

RIOK1 kinase activity is required for cell survival irrespective of MTAP status

Genotype specific vulnerabilities of cancer cells constitute a promising strategy for the development of new therapeutics. Deletions of non-essential genes in tumors can generate unique vulnerabilities which could be exploited therapeutically. The MTAP gene is recurrently deleted in human cancers because of its chromosomal proximity to the tumor suppressor gene CDKN2A. Recent studies have uncovered an increased dependency of MTAP-deleted cancer cells on the function of a PRMT5 containing complex, including WDR77, PRMT5 and the kinase RIOK1. As RIOK1 kinase activity constitutes a potential therapeutic target, we wanted to test if MTAP deletion confers increased sensitivity to RIOK1 inhibition. Using CRISPR/Cas9-mediated genome engineering we generated analog sensitive alleles of RIOK1 in isogenic cell lines differing only by MTAP status. While we were able to independently confirm an increased dependency of MTAP-deleted cells on PRMT5, we did not detect a differential requirement for RIOK1 kinase activity between MTAP-proficient and deficient cells. Our results reveal that the kinase activity of RIOK1 is required for the survival of cancer cell lines irrespective of their MTAP status and cast doubt on the therapeutic exploitability of RIOK1 in the context of MTAP-deleted cancers.

The multi-subunit GID/CTLH E3 ubiquitin ligase promotes cell proliferation and targets the transcription factor Hbp1 for degradation

In yeast, the glucose-induced degradation-deficient (GID) E3 ligase selectively degrades superfluous gluconeogenic enzymes. Here, we identified all subunits of the mammalian GID/CTLH complex and provide a comprehensive map of its hierarchical organization and step-wise assembly. Biochemical reconstitution demonstrates that the mammalian complex possesses inherent E3 ubiquitin ligase activity, using Ube2H as its cognate E2. Deletions of multiple GID subunits compromise cell proliferation, and this defect is accompanied by deregulation of critical cell cycle markers such as the retinoblastoma (Rb) tumor suppressor, phospho-Histone H3 and Cyclin A. We identify the negative regulator of pro-proliferative genes Hbp1 as a bonafide GID/CTLH proteolytic substrate. Indeed, Hbp1 accumulates in cells lacking GID/CTLH activity, and Hbp1 physically interacts and is ubiquitinated in vitro by reconstituted GID/CTLH complexes. Our biochemical and cellular analysis thus demonstrates that the GID/CTLH complex prevents cell cycle exit in G1, at least in part by degrading Hbp1.

Genome-scale identification of cellular pathways required for cell surface recognition

Interactions mediated by cell surface receptors initiate important instructive signaling cues but can be difficult to detect in biochemical assays because they are often highly transient and membrane-embedded receptors are difficult to solubilize in their native conformation. Here, we address these biochemical challenges by using a genome-scale, cell-based genetic screening approach using CRISPR gene knockout technology to identify cellular pathways required for specific cell surface recognition events. By using high-affinity monoclonal antibodies and low-affinity ligands, we determined the necessary screening parameters, including the importance of establishing binding contributions from the glycocalyx, that permitted the unequivocal identification of genes encoding directly interacting membrane-embedded receptors with high statistical confidence. Importantly, we show that this genome-wide screening approach additionally identified receptor-specific pathways that are required for functional display of receptors on the cell surface that included chaperones, enzymes that add post-translational modifications, trafficking proteins, and transcription factors. Finally, we demonstrate the utility of the approach by identifying IGF2R (insulin like growth factor 2 receptor) as a binding partner for the R2 subunit of GABA_B receptors. We show that this interaction is direct and is critically dependent on mannose-6-phosphate, providing a mechanism for the internalization and regulation of GABA_B receptor signaling. We conclude that this single approach can reveal both the molecular nature and the genetic pathways required for functional cell surface display of receptors recognized by antibodies, secreted proteins, and membrane-embedded ligands without the need to make any prior assumptions regarding their biochemical properties.

Genetic Network Complexity Shapes Background-Dependent Phenotypic Expression

The phenotypic consequences of a given mutation can vary across individuals. This so-called background effect is widely observed, from mutant fitness of loss-of-function variants in model organisms to variable disease penetrance and expressivity in humans; however, the underlying genetic basis often remains unclear. Taking insights gained from recent large-scale surveys of genetic interaction and suppression analyses in yeast, we propose that the genetic network context for a given mutation may shape its propensity of exhibiting background-dependent phenotypes. We argue that further efforts in systematically mapping the genetic interaction networks beyond yeast will provide not only key insights into the functional properties of genes, but also a better understanding of the background effects and the (un)predictability of traits in a broader context.

The International Mouse Phenotyping Consortium (IMPC): a functional catalogue of the mammalian genome that informs conservation

The International Mouse Phenotyping Consortium (IMPC) is building a catalogue of mammalian gene function by producing and phenotyping a knockout mouse line for every protein-coding gene. To date, the IMPC has generated and characterised 5186 mutant lines. One-third of the lines have been found to be non-viable and over 300 new mouse models of human disease have been identified thus far. While current bioinformatics efforts are focused on translating results to better understand human disease processes, IMPC data also aids understanding genetic function and processes in other species. Here we show, using gorilla genomic data, how genes essential to development in mice can be used to help assess the potentially deleterious impact of gene variants in other species. This type of analyses could be used to select optimal breeders in endangered species to maintain or increase fitness and avoid variants associated to impaired-health phenotypes or loss-of-function mutations in genes of critical importance. We also show, using selected examples from various mammal species, how IMPC data can aid in the identification of candidate genes for studying a condition of interest, deliver information about the mechanisms involved, or support predictions for the function of genes that may play a role in adaptation. With genotyping costs decreasing and the continued improvements of bioinformatics tools, the analyses we demonstrate can be routinely applied.

MLL-fusion-driven leukemia requires SETD2 to safeguard genomic integrity

MLL-fusions represent a large group of leukemia drivers, whose diversity originates from the vast molecular heterogeneity of C-terminal fusion partners of MLL. While studies of selected MLL-fusions have revealed critical molecular pathways, unifying mechanisms across all MLL-fusions remain poorly understood. We present the first comprehensive survey of protein-protein interactions of seven distantly related MLL-fusion proteins. Functional investigation of 128 conserved MLL-fusion-interactors identifies a specific role for the lysine methyltransferase SETD2 in MLL-leukemia. SETD2 loss causes growth arrest and differentiation of AML cells, and leads to increased DNA damage. In addition to its role in H3K36 tri-methylation, SETD2 is required to maintain high H3K79 di-methylation and MLL-AF9-binding to critical target genes, such as Hoxa9. SETD2 loss synergizes with pharmacologic inhibition of the H3K79 methyltransferase DOT1L to induce DNA damage, growth arrest, differentiation, and apoptosis. These results uncover a dependency for SETD2 during MLL-leukemogenesis, revealing a novel actionable vulnerability in this disease.

Methionine metabolism influences genomic architecture and gene expression through H3K4me3 peak width

Nutrition and metabolism are known to influence chromatin biology and epigenetics through post-translational modifications, yet how this interaction influences genomic architecture and connects to gene expression is unknown. Here we consider, as a model, the metabolically-driven dynamics of H3K4me3, a histone methylation mark that is known to encode information about active transcription, cell identity, and tumor suppression. We analyze the genome-wide changes in H3K4me3 and gene expression in response to alterations in methionine availability in both normal mouse physiology and human cancer cells. Surprisingly, we find that the location of H3K4me3 peaks is largely preserved under methionine restriction, while the response of H3K4me3 peak width encodes almost all aspects of H3K4me3 biology including changes in expression levels, and the presence of cell identity and cancer-associated genes. These findings may reveal general principles for how nutrient availability modulates specific aspects of chromatin dynamics to mediate biological function.

Interrogation of Mammalian Protein Complex Structure, Function, and Membership Using Genome-Scale Fitness Screens

Protein complexes are assemblies of subunits that have co-evolved to execute one or many coordinated functions in the cellular environment. Functional annotation of mammalian protein complexes is critical to understanding biological processes, as well as disease mechanisms. Here, we used genetic co-essentiality derived from genome-scale RNAi- and CRISPR-Cas9-based fitness screens performed across hundreds of human cancer cell lines to assign measures of functional similarity. From these measures, we systematically built and characterized functional similarity networks that recapitulate known structural and functional features of well-studied protein complexes and resolve novel functional modules within complexes lacking structural resolution, such as the mammalian SWI/SNF complex. Finally, by integrating functional networks with large protein-protein interaction networks, we discovered novel protein complexes involving recently evolved genes of unknown function. Taken together, these findings demonstrate the utility of genetic perturbation screens alone, and in combination with large-scale biophysical data, to enhance our understanding of mammalian protein complexes in normal and disease states.