Systematic identification of novel cancer genes through analysis of deep shRNA perturbation screens.

Systematic perturbation screens provide comprehensive resources for the elucidation of cancer driver genes. The perturbation of many genes in relatively few cell lines in such functional screens necessitates the development of specialized computational tools with sufficient statistical power. Here we developed APSiC (Analysis of Perturbation Screens for identifying novel Cancer genes) to identify genetic drivers and effectors in perturbation screens even with few samples. Applying APSiC to the shRNA screen Project DRIVE, APSiC identified well-known and novel putative mutational and amplified cancer genes across all cancer types and in specific cancer types. Additionally, APSiC discovered tumor-promoting and tumor-suppressive effectors, respectively, for individual cancer types, including genes involved in cell cycle control, Wnt/ß-catenin and hippo signalling pathways. We functionally demonstrated that LRRC4B, a putative novel tumor-suppressive effector, suppresses proliferation by delaying cell cycle and modulates apoptosis in breast cancer. We demonstrate APSiC is a robust statistical framework for discovery of novel cancer genes through analysis of large-scale perturbation screens. The analysis of DRIVE using APSiC is provided as a web portal and represents a valuable resource for the discovery of novel cancer genes.

Fully unsupervised deep mode of action learning for phenotyping high-content cellular images.

MOTIVATION: The identification and discovery of phenotypes from high content screening images is a challenging task. Earlier works use image analysis pipelines to extract biological features, supervised training methods or generate features with neural networks pretrained on non-cellular images. We introduce a novel unsupervised deep learning algorithm to cluster cellular images with similar Mode-of-Action (MOA) together using only the images' pixel intensity values as input. It corrects for batch effect during training. Importantly, our method does not require the extraction of cell candidates and works from the entire images directly. RESULTS: The method achieves competitive results on the labeled subset of the BBBC021 dataset with an accuracy of 97.09% for correctly classifying the MOA by nearest neighbors matching. Importantly, we can train our approach on unannotated datasets. Therefore, our method can discover novel MOAs and annotate unlabeled compounds. The ability to train end-to-end on the full resolution images makes our method easy to apply and allows it to further distinguish treatments by their effect on proliferation. AVAILABILITY AND IMPLEMENTATION: Our code is available at https://github.com/Novartis/UMM-Discovery. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ZNRF3 and RNF43 cooperate to safeguard metabolic liver zonation and hepatocyte proliferation.

AXIN2 and LGR5 mark intestinal stem cells (ISCs) that require WNT/ß-Catenin signaling for constant homeostatic proliferation. In contrast, AXIN2/LGR5+ pericentral hepatocytes show low proliferation rates despite a WNT/ß-Catenin activity gradient required for metabolic liver zonation. The mechanisms restricting proliferation in AXIN2+ hepatocytes and metabolic gene expression in AXIN2+ ISCs remained elusive. We now show that restricted chromatin accessibility in ISCs prevents the expression of ß-Catenin-regulated metabolic enzymes, whereas fine-tuning of WNT/ß-Catenin activity by ZNRF3 and RNF43 restricts proliferation in chromatin-permissive AXIN2+ hepatocytes, while preserving metabolic function. ZNRF3 deletion promotes hepatocyte proliferation, which in turn becomes limited by RNF43 upregulation. Concomitant deletion of RNF43 in ZNRF3 mutant mice results in metabolic reprogramming of periportal hepatocytes and induces clonal expansion in a subset of hepatocytes, ultimately promoting liver tumors. Together, ZNRF3 and RNF43 cooperate to safeguard liver homeostasis by spatially and temporally restricting WNT/ß-Catenin activity, balancing metabolic function and hepatocyte proliferation.

DOT1L inhibition is lethal for multiple myeloma due to perturbation of the endoplasmic reticulum stress pathway.

The histone 3 lysine 79 (H3K79) methyltransferase (HMT) DOT1L is known to play a critical role for growth and survival of MLL-rearranged leukemia. Serendipitous observations during high-throughput drug screens indicated that the use of DOT1L inhibitors might be expandable to multiple myeloma (MM). Through pharmacologic and genetic experiments, we could validate that DOT1L is essential for growth and viability of a subset of MM cell lines, in line with a recent report from another team. In vivo activity against established MM xenografts was observed with a novel DOT1L inhibitor. In order to understand the molecular mechanism of the dependency in MM, we examined gene expression changes upon DOT1L inhibition in sensitive and insensitive cell lines and discovered that genes belonging to the endoplasmic reticulum (ER) stress pathway and protein synthesis machinery were specifically suppressed in sensitive cells. Whole-genome CRISPR screens in the presence or absence of a DOT1L inhibitor revealed that concomitant targeting of the H3K4me3 methyltransferase SETD1B increases the effect of DOT1L inhibition. Our results provide a strong basis for further investigating DOT1L and SETD1B as targets in MM.

Next-generation characterization of the Cancer Cell Line Encyclopedia.

Large panels of comprehensively characterized human cancer models, including the Cancer Cell Line Encyclopedia (CCLE), have provided a rigorous framework with which to study genetic variants, candidate targets, and small-molecule and biological therapeutics and to identify new marker-driven cancer dependencies. To improve our understanding of the molecular features that contribute to cancer phenotypes, including drug responses, here we have expanded the characterizations of cancer cell lines to include genetic, RNA splicing, DNA methylation, histone H3 modification, microRNA expression and reverse-phase protein array data for 1,072 cell lines from individuals of various lineages and ethnicities. Integration of these data with functional characterizations such as drug-sensitivity, short hairpin RNA knockdown and CRISPR-Cas9 knockout data reveals potential targets for cancer drugs and associated biomarkers. Together, this dataset and an accompanying public data portal provide a resource for the acceleration of cancer research using model cancer cell lines.

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.

Project DRIVE: A Compendium of Cancer Dependencies and Synthetic Lethal Relationships Uncovered by Large-Scale, Deep RNAi Screening.

Elucidation of the mutational landscape of human cancer has progressed rapidly and been accompanied by the development of therapeutics targeting mutant oncogenes. However, a comprehensive mapping of cancer dependencies has lagged behind and the discovery of therapeutic targets for counteracting tumor suppressor gene loss is needed. To identify vulnerabilities relevant to specific cancer subtypes, we conducted a large-scale RNAi screen in which viability effects of mRNA knockdown were assessed for 7,837 genes using an average of 20 shRNAs per gene in 398 cancer cell lines. We describe findings of this screen, outlining the classes of cancer dependency genes and their relationships to genetic, expression, and lineage features. In addition, we describe robust gene-interaction networks recapitulating both protein complexes and functional cooperation among complexes and pathways. This dataset along with a web portal is provided to the community to assist in the discovery and translation of new therapeutic approaches for cancer.

Resistance mechanisms to TP53-MDM2 inhibition identified by in vivo piggyBac transposon mutagenesis screen in an Arf-/- mouse model.

Inhibitors of double minute 2 protein (MDM2)-tumor protein 53 (TP53) interaction are predicted to be effective in tumors in which the TP53 gene is wild type, by preventing TP53 protein degradation. One such setting is represented by the frequent CDKN2A deletion in human cancer that, through inactivation of p14ARF, activates MDM2 protein, which in turn degrades TP53 tumor suppressor. Here we used piggyBac (PB) transposon insertional mutagenesis to anticipate resistance mechanisms occurring during treatment with the MDM2-TP53 inhibitor HDM201. Constitutive PB mutagenesis in Arf-/- mice provided a collection of spontaneous tumors with characterized insertional genetic landscapes. Tumors were allografted in large cohorts of mice to assess the pharmacologic effects of HDM201. Sixteen out of 21 allograft models were sensitive to HDM201 but ultimately relapsed under treatment. A comparison of tumors with acquired resistance to HDM201 and untreated tumors identified 87 genes that were differentially and significantly targeted by the PB transposon. Resistant tumors displayed a complex clonality pattern suggesting the emergence of several resistant subclones. Among the most frequent alterations conferring resistance, we observed somatic and insertional loss-of-function mutations in transformation-related protein 53 (Trp53) in 54% of tumors and transposon-mediated gain-of-function alterations in B-cell lymphoma-extra large (Bcl-xL), Mdm4, and two TP53 family members, resulting in expression of the TP53 dominant negative truncations ΔNTrp63 and ΔNTrp73. Enhanced BCL-xL and MDM4 protein expression was confirmed in resistant tumors, as well as in HDM201-resistant patient-derived tumor xenografts. Interestingly, concomitant inhibition of MDM2 and BCL-xL demonstrated significant synergy in p53 wild-type cell lines in vitro. Collectively, our findings identify several potential mechanisms by which TP53 wild-type tumors may escape MDM2-targeted therapy.