Defining a Cancer Dependency Map - PubMed (original) (raw)

. 2017 Jul 27;170(3):564-576.e16.

doi: 10.1016/j.cell.2017.06.010.

Francisca Vazquez 2, Phil G Montgomery 1, Barbara A Weir 2, Gregory Kryukov 2, Glenn S Cowley 1, Stanley Gill 2, William F Harrington 1, Sasha Pantel 1, John M Krill-Burger 1, Robin M Meyers 1, Levi Ali 1, Amy Goodale 1, Yenarae Lee 1, Guozhi Jiang 1, Jessica Hsiao 1, William F J Gerath 1, Sara Howell 1, Erin Merkel 1, Mahmoud Ghandi 1, Levi A Garraway 3, David E Root 1, Todd R Golub 4, Jesse S Boehm 1, William C Hahn 5

Affiliations

Defining a Cancer Dependency Map

Aviad Tsherniak et al. Cell. 2017.

Abstract

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 SDs 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.

Keywords: RNAi screens; cancer dependencies; cancer targets; genetic vulnerabilities; genomic biomarkers; precision medicine; predictive modeling; seed effects; shRNA.

Copyright © 2017 Elsevier Inc. All rights reserved.

PubMed Disclaimer

Figures

Figure 1

Figure 1. Computational segregation of on- and off-target effects of RNAi

(A) Tumor types by their contribution to cancer mortality (left) and their cancer cell line representation in the reported dataset (right). (B) Distributions of Pearson correlation coefficients for pairs of shRNA viability profiles before (left) and after (right) removal of inferred seed effects and selection of effective shRNAs (n > 12,000) by DEMETER. Pairs of shRNAs selected randomly (blue lines), targeting the same gene (orange) and sharing a seed sequence (green). (C) Schematic representation of DEMETER and its computational model. Gene- and seed-related effects are estimated from shRNA depletion data. The color of inner circles represents the shRNA target gene and the color of outer circles represents the shRNAs seed sequence. (D) For the top 0.1% most depleted shRNA readouts and the top 0.1% DEMETER gene dependency scores across the whole dataset, the fraction of data points corresponding to a cell line not expressing the target gene. (E) A heatmap depicts the dependency scores (rows) across 501 screened lines (columns) for PIK3CA and the 7 genes that have significantly correlated dependency profiles (z-score > 3). These data were used to plot a gene network, with each edge representing a significant correlation between a pair of dependency profiles. Genes are colored by functional classes. The same analysis was used to generate gene networks for PTK2 (F) and MED12 (G).

Figure 2

Figure 2. The landscape of genetic dependencies in 501 cancer cell lines

(A) Histograms of gene dependency scores for the indicated genes for all cell lines (x-axis). (B) For each differential dependency strength (line color), and for each number of cell lines (xaxis), the number of genes that are differential dependencies is shown (y-axis). (C) Distribution of the number of 6σ dependencies per cell line. (D) Distribution of 6σ dependencies by protein classes. (E) The number of 6σ dependencies annotated as druggable by either being included in DGIdb or IUPHAR/BPS Guide to Pharmacology. (F) The fraction of cell lines (y-axis) that have a 6σ differential dependency on at least one gene in a set of a given size (x-axis). Blue line – considering all 6σ dependencies; orange line – considering only druggable ones.

Figure 3

Figure 3. Prediction of differential dependencies using molecular markers

(A) The number of 6σ dependencies with predictive models built using all features (Unbiased, blue), features of genes related to the dependency gene (Related, red) and those falling into one of the four identified dependency classes (green). (B) Cumulative fraction of 6σ and non-6σ dependencies with predictive models (y-axis) using all features (red bars), plus related features (green), plus those in the four dependency classes (blue). (C) The proportion of the top predictive feature type (copy number, orange; expression, green; mutation, blue) in all unbiased models of 6σ dependencies (D) Top five features of predictive models for 3 gene dependencies in white circles. Circle size is proportional to the relative importance of each feature to the model’s predictive power. (E) Predictive accuracy of ATLANTIS models using only single features (black and colored bars) and using all features (gray bars). (F) Four classes of MDPs, each with a representative example and the top 10 predictable dependencies. Red dotted circles highlight the most dependent cell lines. (Top left) A histogram of GNAI2 dependency scores (x-axis). The two cell lines most dependent on GNAI2 harbor the same indel mutation. (Top right) POU2F2 dependency scores (x-axis) and expression levels (y-axis). Cell lines over-expressing POU2F2 are the most dependent lines. (Bottom left) RPL17 dependency (x-axis) and copy number (y-axis) illustrating a CYCLOPS dependency. (Bottom right) FERMT1 dependency (x-axis) and FERMT2 expression levels (y-axis) for cell lines with low expression of FERMT3 (log2RPKM < 3). Cell lines most dependent on FERMT1 do not express FERMT2.

Figure 4

Figure 4. Oncogene and expression addiction MDPs are enriched in lineage-specific transcription factors

(A) Percentage of transcription factors (TF) among all genes and the four dependency classes. * - P-value < 0.05, *** - P-value < 10−15, Fisher’s exact test. (B) Lineage enrichment (odds ratio; y-axis) of mutation- and expression-driven TF dependencies (N=50) for lineages (x-axis) with significant enrichment (Fisher’s exact FDR < 0.05) in a single (blue) or multiple (orange) lineages. (C) Distributions of 6σ TF dependencies overrepresented in non-essential lineages (ovary, breast, prostate, multiple myeloma and melanoma) compared to known mutation-driven dependencies (BRAF, PIK3CA); dots depict dependency scores greater than 4σ.

Figure 5

Figure 5. UBB/UBC as a paralog deficiency MDP in ovarian cancer cell lines

(A) UBC dependency scores (x-axis) versus UBB expression levels (y-axis). UBB expression (y-axis) versus promoter methylation (x-axis; Fraction) in (B) ovarian cell lines (CCLE data) and (C) tumors (TCGA data). (D) GFP viability competition assay in UBBlow and UBBhigh ovarian cell lines using 4 shRNAs targeting UBC. Log2 fold change of shUBC expressing cells relative to negative controls is shown. (E) Time course of relative viability upon UBC suppression with or without ectopic expression of monoubiquitin (UBB) in a UBBlow cell line (SNU8). Data represent fold change relative to day 1 normalized to pLKO_TRC005-nullT. Error bars represent SD. (F) Levels of conjugated ubiquitin upon UBC suppression in UBBlow (SNU8) and UBBhigh (TOV112D) cell lines.

Figure 6

Figure 6. Effects of scale on the completeness of a Cancer Dependency Map

(A) For each differential dependency with a significant predictive model, the predictive power of the best model (y-axis) and its MDP class (color) along with the strength of the dependency in the most dependent cell line (x-axis). (B) Discovery status of a curated set of 39 mutation- and expression-related dependencies in the dataset. We computed the correlations of each marker with all the differential dependencies and categorized them as (1) Discovered, (2) Not concordant (3) Insufficient context or (4) No differential dependency (see Methods). (C) Fraction of predictable 6σ dependencies, summarized by the number of 6σ dependent cell lines. (D) Results of a down-sampling analysis showing the number of 6σ differential dependencies identified (y-axis) in randomly-sampled subsets of the screened cell lines (x-axis). The blue, orange, green and red lines correspond to dependencies observed in at least 1, 5, 10 or 20 cell lines, respectively.

Similar articles

Cited by

References

    1. Aguirre AJ, Meyers RM, Weir BA, Vazquez F, Zhang CZ, Ben-David U, Cook A, Ha G, Harrington WF, Doshi MB, et al. Genomic Copy Number Dictates a Gene-Independent Cell Response to CRISPR/Cas9 Targeting. Cancer Discov 2016 - PMC - PubMed
    1. Aksoy BA, Demir E, Babur O, Wang W, Jing X, Schultz N, Sander C. Prediction of individualized therapeutic vulnerabilities in cancer from genomic profiles. Bioinformatics. 2014;30:2051–2059. - PMC - PubMed
    1. Ashton JM, Balys M, Neering SJ, Hassane DC, Cowley G, Root DE, Miller PG, Ebert BL, McMurray HR, Land H, et al. Gene sets identified with oncogene cooperativity analysis regulate in vivo growth and survival of leukemia stem cells. Cell Stem Cell. 2012;11:359–372. - PMC - PubMed
    1. Barbie DA, Tamayo P, Boehm JS, Kim SY, Moody SE, Dunn IF, Schinzel AC, Sandy P, Meylan E, Scholl C, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009;462:108–112. - PMC - PubMed
    1. Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, Wilson CJ, Lehar J, Kryukov GV, Sonkin D, et al. The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483:603–607. - PMC - PubMed

MeSH terms

Substances

Grants and funding

LinkOut - more resources