SUMMARY Most human epithelial tumors harbor numerous alterations, making it difficult to predict which genes are required for tumor survival. To systematically identify cancer dependencies, we analyzed 501 genome-scale loss-of-function screens performed in diverse human cancer cell lines. We developed DEMETER, an analytical framework that segregates on-from off-target effects of RNAi. 769 genes were differentially required in subsets of these cell lines at a threshold of six standard deviations from the mean. We found predictive models for 426 dependencies (55%) by nonlinear regression modeling considering 66,646 molecular features. Many dependencies fall into a limited number of classes, and unexpectedly, in 82% of models, the top biomarkers were expression-based. We demonstrated the basis behind one such predictive model linking hypermethylation of the UBB ubiquitin gene to a dependency on UBC. Together, these observations provide a foundation for a cancer dependency map that facilitates the prioritization of therapeutic targets.

The CRISPR-Cas9 system has revolutionized gene editing both on single genes and in multiplexed loss-of-function screens, enabling precise genome-scale identification of genes essential to proliferation and survival of cancer cells1,2. However, previous studies reported that a gene-independent anti-proliferative effect of Cas9-mediated DNA cleavage confounds such measurement of genetic dependency, leading to false positive results in copy number amplified regions3,4. We developed CERES, a computational method to estimate gene dependency levels from CRISPR-Cas9 essentiality screens while accounting for the copy-number-specific effect. As part of our efforts to define a cancer dependency map, we performed genome-scale CRISPR-Cas9 essentiality screens across 342 cancer cell lines and applied CERES to this dataset. We found that CERES reduced false positive results and estimated sgRNA activity for both this dataset and previously published screens performed with different sgRNA libraries. Here, we demonstrate the utility of this collection of screens, upon CERES correction, in revealing cancer-type-specific vulnerabilities.

Major efforts using loss-of-function genetic screens to systematically identify genes essential to the proliferation and survival of cancer cells have been reported [1][2][3][4][5][6][7][8][9] . Genes identified by these approaches may represent specific genetic vulnerabilities of cancer cells, suggesting treatment strategies and directing the development of novel therapeutics. The CRISPR-Cas9 genome editing system has proven to be a powerful tool to interrogate gene essentiality in cancer cell lines. Its relative ease of application, high rates of target validation, and increased specificity compared to RNA interference technology make it an ideal instrument for use in high-throughput functional genomic screening 10 .However, we and others have recently observed that measurements of genetic dependency in genome-scale CRISPR-Cas9 loss-of-function screens are influenced by the genomic copy number (CN) of the region targeted by the sgRNA-Cas9 complex [1][2][3][4] . Targeting Cas9 to DNA sequences within regions of high CN gain creates multiple DNA double-strand breaks (DSBs), inducing a gene-independent DNA damage response and a G2 cell-cycle arrest phenotype 2 .This systematic, sequence-independent effect due to DNA cleavage (copy-number effect)confounds the measurement of the consequences of gene deletion on cell viability (geneknockout effect) and is detectable even among low-level CN amplifications and deletions. In particular, this phenomenon hinders interpretation of CRISPR-Cas9 experiments in cancer cell

Although great progress has been made in the characterization of off-target effects of engineered nucleases, sensitive and unbiased genome-wide methods for the detection of off-target cleavage events and potential collateral damage are still lacking. Here we describe a linear amplification–mediated modification of a previously published high-throughput, genome-wide translocation sequencing (HTGTS) method that robustly detects DNA double-stranded breaks (DSBs) generated by engineered nucleases across the human genome based on their translocation to other endogenous or ectopic DSBs. HTGTS with different Cas9:sgRNA or TALEN-nucleases revealed off-target hotspots for given nucleases that ranged from a few or none to dozens or more, and extended the number of known off-targets for certain previously characterized nucleases by more than 10-fold. We also identified translocations between bona fide nuclease targets on homologous chromosomes, an undesired collateral effect that has not been described. Finally, HTGTS confirmed that the Cas9D10A paired nickase approach suppresses off-target cleavage genome-wide.