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

Summary Due to genome instability, most cancers exhibit loss of regions containing tumor suppressor genes and collateral loss of other genes. To identify cancer-specific vulnerabilities that are the result of copy-number losses, we performed integrated analyses of genome-wide copy-number and RNAi profiles and identified 56 genes for which gene suppression specifically inhibited the proliferation of cells harboring partial copy-number loss of that gene. These CYCLOPS (Copy-number alterations Yielding Cancer Liabilities Owing to Partial losS) genes are enriched for spliceosome, proteasome and ribosome components. One CYCLOPS gene, PSMC2, encodes an essential member of the 19S proteasome. Normal cells express excess PSMC2, which resides in a complex with PSMC1, PSMD2, and PSMD5 and acts as a reservoir protecting cells from PSMC2 suppression. Cells harboring partial PSMC2 copy-number loss lack this complex and die after PSMC2 suppression. These observations define a distinct class of cancer-specific liabilities resulting from genome instability.

Brain metastases from lung adenocarcinoma (BM-LUAD) cause significant patient mortality. To identify genomic alterations that promote brain metastases, we performed whole-exome sequencing of 73 BM-LUAD cases. Using case-control analyses, we discovered candidate drivers of brain metastasis by identifying genes with more frequent copy-number aberrations in BM-LUAD compared to 503 primary lung adenocarcinomas. We identified three regions with significantly higher amplification frequencies in BM-LUAD, including MYC (12% vs 6%), YAP1 (7% vs 0.8%), and MMP13 (10% vs 0.6%) and significantly more frequent deletions in CDKN2A/B (27% vs 13%). We confirmed that amplification frequencies of MYC and YAP1 / MMP13 were elevated in an independent cohort of 105 patients. Functional assessment in patient-derived xenograft mouse models validated that MYC , YAP1 or MMP13 overexpression increased the brain metastasis incidence. These results demonstrate that somatic alterations contribute to brain metastases and that genomic sequencing of a large number of metastatic tumors can reveal novel metastatic drivers.