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.
Display omitted
► Partial copy number loss renders cancer cells dependent on CYCLOPS genes ► CYCLOPS genes are enriched for components of the proteasome, including PSMC2 ► Copy number loss of PSMC2 leads to lower expression and loss of a PSMC2 reservoir ► Targeting PSMC2 disrupts tumor maintenance in a genotype-specific manner
Reducing the mRNA levels of genes that have suffered copy number loss due to genomic instability in cancer cells leads to tumor-cell-specific growth inhibition. Thus, probing cancer genomes for genes that have copy number losses can reveal new therapeutic targets.
Abstract
Background
Hundreds of functional genomic screens have been performed across a diverse set of cancer contexts, as part of efforts such as the Cancer Dependency Map, to identify gene ...dependencies—genes whose loss of function reduces cell viability or fitness. Recently, large-scale screening efforts have shifted from RNAi to CRISPR-Cas9, due to superior efficacy and specificity. However, many effective oncology drugs only partially inhibit their protein targets, leading us to question whether partial suppression of genes using RNAi could reveal cancer vulnerabilities that are missed by complete knockout using CRISPR-Cas9. Here, we compare CRISPR-Cas9 and RNAi dependency profiles of genes across approximately 400 matched cancer cell lines.
Results
We find that CRISPR screens accurately identify more gene dependencies per cell line, but the majority of each cell line’s dependencies are part of a set of 1867 genes that are shared dependencies across the entire collection (pan-lethals). While RNAi knockdown of about 30% of these genes is also pan-lethal, approximately 50% have selective dependency patterns across cell lines, suggesting they could still be cancer vulnerabilities. The accuracy of the unique RNAi selectivity is supported by associations to multi-omics profiles, drug sensitivity, and other expected co-dependencies.
Conclusions
Incorporating RNAi data for genes that are pan-lethal knockouts facilitates the discovery of a wider range of gene targets than could be detected using the CRISPR dataset alone. This can aid in the interpretation of contrasting results obtained from CRISPR and RNAi screens and reinforce the importance of partial gene suppression methods in building a cancer dependency map.
Clear cell renal cell carcinoma (ccRCC) is the most common form of kidney cancer and is often linked to loss of chromosome 3p, which harbors the VHL tumor suppressor gene, loss of chromosome 14q, ...which includes HIF1A, and gain of chromosome 5q. The relevant target(s) on chromosome 5q is not known. Here, we show that 5q amplification leads to overexpression of the SQSTM1 oncogene in ccRCC lines and tumors. Overexpression of SQSTM1 in ccRCC lines promoted resistance to redox stress and increased soft agar growth, while downregulation of SQSTM1 decreased resistance to redox stress, impaired cellular fitness, and decreased tumor formation. Therefore, the selection pressure to amplify 5q in ccRCC is driven, at least partly, by SQSTM1.
•Copy number gains of chromosome 5q drive SQSTM1 overexpression•SQSTM1 gene product p62 activates NRF2 and promotes resistance to redox stress•mTOR acts downstream or parallel of p62 in renal cancer cells•p62 promotes renal cancer growth in vitro and in vivo
The mitogen-activated protein kinase (MAPK) pathway is a critical effector of oncogenic RAS signaling, and MAPK pathway inhibition may be an effective combination treatment strategy. We performed ...genome-scale loss-of-function CRISPR-Cas9 screens in the presence of a MEK1/2 inhibitor (MEKi) in KRAS-mutant pancreatic and lung cancer cell lines and identified genes that cooperate with MEK inhibition. While we observed heterogeneity in genetic modifiers of MEKi sensitivity across cell lines, several recurrent classes of synthetic lethal vulnerabilities emerged at the pathway level. Multiple members of receptor tyrosine kinase (RTK)-RAS-MAPK pathways scored as sensitizers to MEKi. In particular, we demonstrate that knockout, suppression, or degradation of SHOC2, a positive regulator of MAPK signaling, specifically cooperated with MEK inhibition to impair proliferation in RAS-driven cancer cells. The depletion of SHOC2 disrupted survival pathways triggered by feedback RTK signaling in response to MEK inhibition. Thus, these findings nominate SHOC2 as a potential target for combination therapy.
Display omitted
•RTK-RAS-MAPK pathway members score strongly in genome-scale MEKi modifier screens•Depletion of SHOC2 potently sensitizes RAS-driven cells to MEK inhibition•SHOC2 loss impairs RTK-mediated adaptive reactivation of MAPK signaling induced by MEKi•A model of SHOC2 degradation suggests a combination therapeutic strategy with MEKi
Sulahian, Kwon, and Walsh et al. performed several loss-of-function CRISPR-Cas9 screens in KRAS-mutant cancer cells treated with a MEK inhibitor and define the landscape of modifiers of MEK inhibitor sensitivity while highlighting that SHOC2 is a potent synthetic lethal target that serves as a critical signaling node to mediate MAP kinase pathway reactivation upon MEK inhibition.
Genome-scale RNAi libraries enable the systematic interrogation of gene function. However, the interpretation of RNAi screens is complicated by the observation that RNAi reagents designed to suppress ...the mRNA transcripts of the same gene often produce a spectrum of phenotypic outcomes due to differential on-target gene suppression or perturbation of off-target transcripts. Here we present a computational method, Analytic Technique for Assessment of RNAi by Similarity (ATARiS), that takes advantage of patterns in RNAi data across multiple samples in order to enrich for RNAi reagents whose phenotypic effects relate to suppression of their intended targets. By summarizing only such reagent effects for each gene, ATARiS produces quantitative, gene-level phenotype values, which provide an intuitive measure of the effect of gene suppression in each sample. This method is robust for data sets that contain as few as 10 samples and can be used to analyze screens of any number of targeted genes. We used this analytic approach to interrogate RNAi data derived from screening more than 100 human cancer cell lines and identified HNF1B as a transforming oncogene required for the survival of cancer cells that harbor HNF1B amplifications. ATARiS is publicly available at http://broadinstitute.org/ataris.
Ewing sarcoma is an aggressive bone and soft tissue tumor in children and adolescents, with treatment remaining a clinical challenge. This disease is mediated by somatic chromosomal translocations of ...the EWS gene and a gene encoding an ETS transcription factor, most commonly, FLI1. While direct targeting of aberrant transcription factors remains a pharmacological challenge, identification of dependencies incurred by EWS/FLI1 expression would offer a new therapeutic avenue. We used a combination of super-enhancer profiling, near-whole genome shRNA-based and small-molecule screening to identify cyclin D1 and CDK4 as Ewing sarcoma-selective dependencies. We revealed that super-enhancers mark Ewing sarcoma specific expression signatures and EWS/FLI1 target genes in human Ewing sarcoma cell lines. Particularly, a super-enhancer regulates cyclin D1 and promotes its expression in Ewing sarcoma. We demonstrated that Ewing sarcoma cells require CDK4 and cyclin D1 for survival and anchorage-independent growth. Additionally, pharmacologic inhibition of CDK4 with selective CDK4/6 inhibitors led to cytostasis and cell death of Ewing sarcoma cell lines in vitro and growth delay in an in vivo Ewing sarcoma xenograft model. These results demonstrated a dependency in Ewing sarcoma on CDK4 and cyclin D1 and support exploration of CDK4/6 inhibitors as a therapeutic approach for patients with this disease.
Through an shRNA screen, we identified the protein arginine methyltransferase Prmt1 as a vulnerable intervention point in murine p53/Rb-null osteosarcomas, the human counterpart of which lacks ...effective therapeutic options. Depletion of Prmt1 in p53-deficient cells impaired tumor initiation and maintenance
and
Mechanistic studies reveal that translation-associated pathways were enriched for Prmt1 downstream targets, implicating Prmt1 in translation control. In particular, loss of Prmt1 led to a decrease in arginine methylation of the translation initiation complex, thereby disrupting its assembly and inhibiting translation. p53/Rb-null cells were sensitive to p53-induced translation stress, and analysis of human cancer cell line data from Project Achilles further revealed that Prmt1 and translation-associated pathways converged on the same functional networks. We propose that targeted therapy against Prmt1 and its associated translation-related pathways offer a mechanistic rationale for treatment of osteosarcomas and other cancers that exhibit dependencies on translation stress response.
.
Genomic instability is a hallmark of human cancer, and results in widespread somatic copy number alterations. We used a genome-scale shRNA viability screen in human cancer cell lines to ...systematically identify genes that are essential in the context of particular copy-number alterations (copy-number associated gene dependencies). The most enriched class of copy-number associated gene dependencies was CYCLOPS (Copy-number alterations Yielding Cancer Liabilities Owing to Partial losS) genes, and spliceosome components were the most prevalent. One of these, the pre-mRNA splicing factor
, is also frequently mutated in cancer. We validated
as a CYCLOPS gene and found that human cancer cells harboring partial
copy-loss lack a reservoir of SF3b complex that protects cells with normal
copy number from cell death upon partial
suppression. These data provide a catalog of copy-number associated gene dependencies and identify partial copy-loss of wild-type
as a novel, non-driver cancer gene dependency.
Despite significant progress in the molecular understanding of medulloblastoma, stratification of risk in patients remains a challenge. Focus has shifted from clinical parameters to molecular ...markers, such as expression of specific genes and selected genomic abnormalities, to improve accuracy of treatment outcome prediction. Here, we show how integration of high-level clinical and genomic features or risk factors, including disease subtype, can yield more comprehensive, accurate, and biologically interpretable prediction models for relapse versus no-relapse classification. We also introduce a novel Bayesian nomogram indicating the amount of evidence that each feature contributes on a patient-by-patient basis.
A Bayesian cumulative log-odds model of outcome was developed from a training cohort of 96 children treated for medulloblastoma, starting with the evidence provided by clinical features of metastasis and histology (model A) and incrementally adding the evidence from gene-expression-derived features representing disease subtype-independent (model B) and disease subtype-dependent (model C) pathways, and finally high-level copy-number genomic abnormalities (model D). The models were validated on an independent test cohort (n = 78).
On an independent multi-institutional test data set, models A to D attain an area under receiver operating characteristic (au-ROC) curve of 0.73 (95% CI, 0.60 to 0.84), 0.75 (95% CI, 0.64 to 0.86), 0.80 (95% CI, 0.70 to 0.90), and 0.78 (95% CI, 0.68 to 0.88), respectively, for predicting relapse versus no relapse.
The proposed models C and D outperform the current clinical classification schema (au-ROC, 0.68), our previously published eight-gene outcome signature (au-ROC, 0.71), and several new schemas recently proposed in the literature for medulloblastoma risk stratification.
The identification of somatic genetic alterations that confer sensitivity to pharmacologic inhibitors has led to new cancer therapies. To identify mutations that confer an exceptional dependency, ...shRNA-based loss-of-function data were analyzed from a dataset of numerous cell lines to reveal genes that are essential in a small subset of cancer cell lines. Once these cell lines were determined, detailed genomic characterization from these cell lines was utilized to ascertain the genomic aberrations that led to this extreme dependency. This method, in a large subset of lung cancer cell lines, yielded a single lung adenocarcinoma cell line, NCI-H1437, which is sensitive to RNA interference of MAP2K1 expression. Notably, NCI-H1437 is the only lung cancer cell line included in the dataset with a known activating mutation in MAP2K1 (Q56P). Subsequent validation using shRNA and CRISPR-Cas9 confirmed MAP2K1 dependency. In vitro and in vivo inhibitor studies established that NCI-H1437 cells are sensitive to MEK1 inhibitors, including the FDA-approved drug trametinib. Like NCI-H1437 cells, the MAP2K1-mutant cell lines SNU-C1 (colon) and OCUM-1 (gastric) showed decreased viability after MAP2K1 depletion via Cas9-mediated gene editing. Similarly, these cell lines were particularly sensitive to trametinib treatment compared with control cell lines. On the basis of these data, cancers that harbor driver mutations in MAP2K1 could benefit from treatment with MEK1 inhibitors. Furthermore, this functional data mining approach provides a general method to experimentally test genomic features that confer dependence in tumors.
Cancers with an activated RAS/MAPK pathway driven by oncogenic MAP2K1 mutations may be particularly sensitive to MEK1 inhibitor treatments.