A melanocyte lineage program confers resistance to MAP kinase pathway inhibition

BRAFV600E-mutant malignant melanomas depend on RAF/MEK/ERK (MAPK) signaling for tumor cell growth1. RAF and MEK inhibitors show remarkable clinical efficacy in BRAFV600E melanoma2, 3; however, resistance to these agents remains a formidable challenge2, 4. Global characterization of resistance mechanisms may inform the development of more effective therapeutic combinations. Here, we performed systematic gain-of-function resistance studies by expressing >15,500 genes individually in a BRAFV600E melanoma cell line treated with RAF, MEK, ERK, or combined RAF/MEK inhibitors. These studies revealed a cyclic AMP-dependent melanocytic signaling network not previously associated with drug resistance, including G-protein coupled receptors, adenyl cyclase, protein kinase A and cAMP response element binding protein (CREB). Preliminary analysis of biopsies from BRAFV600E melanoma patients revealed that phosphorylated (active) CREB was suppressed by RAF/MEK-inhibition but restored in relapsing tumors. Expression of transcription factors activated downstream of MAP kinase and cAMP pathways also conferred resistance, including c-FOS, NR4A1, NR4A2 and MITF. Combined treatment with MAP kinase pathway and histone deacetylase inhibitors suppressed MITF expression and cAMP-mediated resistance. Collectively, these data suggest that oncogenic dysregulation of a melanocyte lineage dependency can cause resistance to RAF/MEK/ERK inhibition, which may be overcome by combining signaling- and chromatin-directed therapeutics

Characterizing genomic alterations in cancer by complementary functional associations.

Systematic efforts to sequence the cancer genome have identified large numbers of mutations and copy number alterations in human cancers. However, elucidating the functional consequences of these variants, and their interactions to drive or maintain oncogenic states, remains a challenge in cancer research. We developed REVEALER, a computational method that identifies combinations of mutually exclusive genomic alterations correlated with functional phenotypes, such as the activation or gene dependency of oncogenic pathways or sensitivity to a drug treatment. We used REVEALER to uncover complementary genomic alterations associated with the transcriptional activation of β-catenin and NRF2, MEK-inhibitor sensitivity, and KRAS dependency. REVEALER successfully identified both known and new associations, demonstrating the power of combining functional profiles with extensive characterization of genomic alterations in cancer genomes

Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors

Whole-exome sequencing of cell-free DNA (cfDNA) could enable comprehensive profiling of tumors from blood but the genome-wide concordance between cfDNA and tumor biopsies is uncertain. Here we report ichorCNA, software that quantifies tumor content in cfDNA from 0.1× coverage whole-genome sequencing data without prior knowledge of tumor mutations. We apply ichorCNA to 1439 blood samples from 520 patients with metastatic prostate or breast cancers. In the earliest tested sample for each patient, 34% of patients have ≥10% tumor-derived cfDNA, sufficient for standard coverage whole-exome sequencing. Using whole-exome sequencing, we validate the concordance of clonal somatic mutations (88%), copy number alterations (80%), mutational signatures, and neoantigens between cfDNA and matched tumor biopsies from 41 patients with ≥10% cfDNA tumor content. In summary, we provide methods to identify patients eligible for comprehensive cfDNA profiling, revealing its applicability to many patients, and demonstrate high concordance of cfDNA and metastatic tumor whole-exome sequencing

Regulation of mTOR and Cell Growth in Response to Energy Stress by REDD1

The tuberous sclerosis tumor suppressors TSC1 and TSC2 regulate the mTOR pathway to control translation and cell growth in response to nutrient and growth factor stimuli. We have recently identified the stress response REDD1 gene as a mediator of tuberous sclerosis complex (TSC)-dependent mTOR regulation by hypoxia. Here, we demonstrate that REDD1 inhibits mTOR function to control cell growth in response to energy stress. Endogenous REDD1 is induced following energy stress, and REDD1(−)(/)(−) cells are highly defective in dephosphorylation of the key mTOR substrates S6K and 4E-BP1 following either ATP depletion or direct activation of the AMP-activated protein kinase (AMPK). REDD1 likely acts on the TSC1/2 complex, as regulation of mTOR substrate phosphorylation by REDD1 requires TSC2 and is blocked by overexpression of the TSC1/2 downstream target Rheb but is not blocked by inhibition of AMPK. Tetracycline-inducible expression of REDD1 triggers rapid dephosphorylation of S6K and 4E-BP1 and significantly decreases cellular size. Conversely, inhibition of endogenous REDD1 by short interfering RNA increases cell size in a rapamycin-sensitive manner, and REDD1(−)(/)(−) cells are defective in cell growth regulation following ATP depletion. These results define REDD1 as a critical transducer of the cellular response to energy depletion through the TSC-mTOR pathway

Comprehensive Mutational Analysis of the BRCA1-Associated DNA Helicase and Tumor-Suppressor FANCJ/BACH1/BRIP1

FANCJ (BRIP1/BACH1) is a hereditary breast and ovarian cancer (HBOC) gene encoding a DNA helicase. Similar to HBOC genes, BRCA1 and BRCA2, FANCJ is critical for processing DNA inter-strand crosslinks (ICL) induced by chemotherapeutics, such as cisplatin. Consequently, cells deficient in FANCJ or its catalytic activity are sensitive to ICL-inducing agents. Unfortunately, the majority of FANCJ clinical mutations remain uncharacterized, limiting therapeutic opportunities to effectively use cisplatin to treat tumors with mutated FANCJ. Here, we sought to perform a comprehensive screen to identify FANCJ loss-of-function (LOF) mutations. We developed a FANCJ lentivirus mutation library representing approximately 450 patient-derived FANCJ nonsense and missense mutations to introduce FANCJ mutants into FANCJ knockout (K/O) HeLa cells. We performed a high-throughput screen to identify FANCJ LOF mutants that, as compared with wild-type FANCJ, fail to robustly restore resistance to ICL-inducing agents, cisplatin or mitomycin C (MMC). On the basis of the failure to confer resistance to either cisplatin or MMC, we identified 26 missense and 25 nonsense LOF mutations. Nonsense mutations elucidated a relationship between location of truncation and ICL sensitivity, as the majority of nonsense mutations before amino acid 860 confer ICL sensitivity. Further validation of a subset of LOF mutations confirmed the ability of the screen to identify FANCJ mutations unable to confer ICL resistance. Finally, mapping the location of LOF mutations to a new homology model provides additional functional information. IMPLICATIONS: We identify 51 FANCJ LOF mutations, providing important classification of FANCJ mutations that will afford additional therapeutic strategies for affected patients

Creation of Novel Protein Variants with CRISPR/Cas9-Mediated Mutagenesis: Turning a Screening By-Product into a Discovery Tool

<div><p>CRISPR/Cas9 screening has proven to be a versatile tool for genomics research. Based on unexpected results from a genome-wide screen, we developed a CRISPR/Cas9-mediated approach to mutagenesis, exploiting the allelic diversity generated by error-prone non-homologous end-joining (NHEJ) to identify novel gain-of-function and drug resistant alleles of the MAPK signaling pathway genes MEK1 and BRAF. We define the parameters of a scalable technique to easily generate cell populations containing thousands of endogenous allelic variants to map gene functions. Further, these results highlight an unexpected but important phenomenon, that Cas9-induced gain-of-function alleles are an inherent by-product of normal Cas9 loss-of-function screens and should be investigated during analysis of data from large-scale positive selection screens.</p></div

CloneSifter: enrichment of rare clones from heterogeneous cell populations

Abstract Background Many biological processes, such as cancer metastasis, organismal development, and acquisition of resistance to cytotoxic therapy, rely on the emergence of rare sub-clones from a larger population. Understanding how the genetic and epigenetic features of diverse clones affect clonal fitness provides insight into molecular mechanisms underlying selective processes. While large-scale barcoding with NGS readout has facilitated cellular fitness assessment at the population level, this approach does not support characterization of clones prior to selection. Single-cell genomics methods provide high biological resolution, but are challenging to scale across large populations to probe rare clones and are destructive, limiting further functional analysis of important clones. Results Here, we develop CloneSifter, a methodology for tracking and enriching rare clones throughout their response to selection. CloneSifter utilizes a CRISPR sgRNA-barcode library that facilitates the isolation of viable cells from specific clones within the barcoded population using a sequence-specific retrieval reporter. We demonstrate that CloneSifter can measure clonal fitness of cancer cell models in vitro and retrieve targeted clones at abundance as low as 1 in 1883 in a heterogeneous cell population. Conclusions CloneSifter provides a means to track and access specific and rare clones of interest across dynamic changes in population structure to comprehensively explore the basis of these changes

MicroSCALE Screening Reveals Genetic Modifiers of Therapeutic Response in Melanoma

Cell microarrays are a promising tool for performing large-scale functional genomic screening in mammalian cells at reasonable cost, but owing to technical limitations they have been restricted for use with a narrow range of cell lines and short-term assays. Here, we describe MicroSCALE (Microarrays of Spatially Confined Adhesive Lentiviral Features), a cell microarray–based platform that enables application of this technology to a wide range of cell types and longer-term assays. We used MicroSCALE to uncover kinases that when overexpressed partially desensitized B-RAF[superscript V600E]–mutant melanoma cells to inhibitors of the mitogen-activated protein kinase kinase kinase (MAPKKK) RAF, the MAPKKs MEK1 and 2 (MEK1/2, mitogen-activated protein kinase kinase 1 and 2), mTOR (mammalian target of rapamycin), or PI3K (phosphatidylinositol 3-kinase). These screens indicated that cells treated with inhibitors acting through common mechanisms were affected by a similar profile of overexpressed proteins. In contrast, screens involving inhibitors acting through distinct mechanisms yielded unique profiles, a finding that has potential relevance for small-molecule target identification and combination drugging studies. Further, by integrating large-scale functional screening results with cancer cell line gene expression and pharmacological sensitivity data, we validated the nuclear factor κB pathway as a potential mediator of resistance to MAPK pathway inhibitors. The MicroSCALE platform described here may enable new classes of large-scale, resource-efficient screens that were not previously feasible, including those involving combinations of cell lines, perturbations, and assay outputs or those involving limited numbers of cells and limited or expensive reagents.Broad Institute of MIT and Harvard (Scientific Planning and Allocation of Resources Committee Grant

Large variation generated by an individual sgRNA can be used to select for gain-of-function alleles.

<p>(A) Schematic of the experimental approach to examine variants generated by sgRNA targeting and identify gain-of-function alleles. (B) For each of three conditions, variants detected by deep sequencing of the MAP2K1 locus were ranked by their abundance and the reads per million is plotted. Left panel shows all the variants, right panel enlarges one region. (C) Alignment and abundances of selected MAP2K1 alleles generated in untreated or vemurafenib-selected A375 cells by an individual sgRNA targeting MAP2K1 at the site encoding K59. Mismatched and inserted nucleotides are shown in red. Variants that are studied in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0170445#pone.0170445.g002" target="_blank">Fig 2</a> are labeled in blue. Variants are ranked by within-sample abundance.</p