Search

The epidermal growth factor receptor, EGFR, is frequently activated in lung cancer and glioblastoma by genomic alterations including missense mutations. The different mutation spectra in these diseases are reflected in divergent responses to EGFR inhibition: significant patient benefit in lung cancer, but limited in glioblastoma. Here, we report a comprehensive mutational analysis of EGFR function. We perform saturation mutagenesis of EGFR and assess function of ~22,500 variants in a human EGFR-dependent lung cancer cell line. This approach reveals enrichment of erlotinib-insensitive variants of known and unknown significance in the dimerization, transmembrane, and kinase domains. Multiple EGFR extracellular domain variants, not associated with approved targeted therapies, are sensitive to afatinib and dacomitinib in vitro. Two glioblastoma patients with somatic EGFR G598V dimerization domain mutations show responses to dacomitinib treatment followed by within-pathway resistance mutation in one case. In summary, this comprehensive screen expands the landscape of functional EGFR variants and suggests broader clinical investigation of EGFR inhibition for cancers harboring extracellular domain mutations., (© 2024. The Author(s).)

The ability to understand and predict variable responses to therapeutic agents may improve outcomes in patients with cancer. We hypothesized that the basal gene-transcription state of cancer cell lines, coupled with cell viability profiles of small molecules, might be leveraged to nominate specific mechanisms of intrinsic resistance and to predict drug combinations that overcome resistance. We analyzed 564,424 sensitivity profiles to identify candidate gene-compound pairs, and validated nine such relationships. We determined the mechanism of a novel relationship, in which expression of the serine hydrolase enzymes monoacylglycerol lipase (MGLL) or carboxylesterase 1 (CES1) confers resistance to the histone lysine demethylase inhibitor GSK-J4 by direct enzymatic modification. Insensitive cell lines could be sensitized to GSK-J4 by inhibition or gene knockout. These analytical and mechanistic studies highlight the potential of integrating gene-expression features with small-molecule response to identify patient populations that are likely to benefit from treatment, to nominate rational candidates for combinations and to provide insights into mechanisms of action., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Although single-gene perturbation screens have revealed a number of new targets, vulnerabilities specific to frequently altered drivers have not been uncovered. An important question is whether the compensatory relationship between functionally redundant genes masks potential therapeutic targets in single-gene perturbation studies. To identify digenic dependencies, we developed a CRISPR paralog targeting library to investigate the viability effects of disrupting 3,284 genes, 5,065 paralog pairs and 815 paralog families. We identified that dual inactivation of DUSP4 and DUSP6 selectively impairs growth in NRAS and BRAF mutant cells through the hyperactivation of MAPK signaling. Furthermore, cells resistant to MAPK pathway therapeutics become cross-sensitized to DUSP4 and DUSP6 perturbations such that the mechanisms of resistance to the inhibitors reinforce this mechanism of vulnerability. Together, multigene perturbation technologies unveil previously unrecognized digenic vulnerabilities that may be leveraged as new therapeutic targets in cancer., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

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., (©2021 American Association for Cancer Research.)

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.

Purpose: Rhabdoid tumors are devastating pediatric cancers in need of improved therapies. We sought to identify small molecules that exhibit in vitro and in vivo efficacy against preclinical models of rhabdoid tumor., Experimental Design: We screened eight rhabdoid tumor cell lines with 481 small molecules and compared their sensitivity with that of 879 other cancer cell lines. Genome-scale CRISPR-Cas9 inactivation screens in rhabdoid tumors were analyzed to confirm target vulnerabilities. Gene expression and CRISPR-Cas9 data were queried across cell lines and primary rhabdoid tumors to discover biomarkers of small-molecule sensitivity. Molecular correlates were validated by manipulating gene expression. Subcutaneous rhabdoid tumor xenografts were treated with the most effective drug to confirm in vitro results., Results: Small-molecule screening identified the protein-translation inhibitor homoharringtonine (HHT), an FDA-approved treatment for chronic myelogenous leukemia (CML), as the sole drug to which all rhabdoid tumor cell lines were selectively sensitive. Validation studies confirmed the sensitivity of rhabdoid tumor to HHT was comparable with that of CML cell lines. Low expression of the antiapoptotic gene BCL2L1 , which encodes Bcl-XL, was the strongest predictor of HHT sensitivity, and HHT treatment consistently depleted Mcl-1, the synthetic-lethal antiapoptotic partner of Bcl-XL. Rhabdoid tumor cell lines and primary-tumor samples expressed low BCL2L1 , and overexpression of BCL2L1 induced resistance to HHT in rhabdoid tumor cells. Furthermore, HHT treatment inhibited rhabdoid tumor cell line and patient-derived xenograft growth in vivo ., Conclusions: Rhabdoid tumor cell lines and xenografts are highly sensitive to HHT, at least partially due to their low expression of BCL2L1 . HHT may have therapeutic potential against rhabdoid tumors., (©2020 American Association for Cancer Research.)

Cancer evolution poses a central obstacle to cure, as resistant clones expand under therapeutic selection pressures. Genome sequencing of relapsed disease can nominate genomic alterations conferring resistance but sample collection lags behind, limiting therapeutic innovation. Genome-wide screens offer a complementary approach to chart the compendium of escape genotypes, anticipating clinical resistance. We report genome-wide open reading frame (ORF) resistance screens for first- and third-generation epidermal growth factor receptor (EGFR) inhibitors and a MEK inhibitor. Using serial sampling, dose gradients, and mathematical modeling, we generate genotype-fitness maps across therapeutic contexts and identify alterations that escape therapy. Our data expose varying dose-fitness relationship across genotypes, ranging from complete dose invariance to paradoxical dose dependency where fitness increases in higher doses. We predict fitness with combination therapy and compare these estimates to genome-wide fitness maps of drug combinations, identifying genotypes where combination therapy results in unexpected inferior effectiveness. These data are applied to nominate combination optimization strategies to forestall resistant disease., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Kinases are involved in disease development and modulation of their activity can be therapeutically beneficial. Drug-resistant mutant kinases are valuable tools in drug discovery efforts, but the prediction of mutants across the kinome is challenging. Here, we generate deep mutational scanning data to identify mutant mammalian kinases that drive resistance to clinically relevant inhibitors. We aggregate these data with subsaturation mutagenesis data and use it to develop, test and validate a framework to prospectively identify residues that mediate kinase activity and drug resistance across the kinome. We validate predicted resistance mutations in CDK4, CDK6, ERK2, EGFR and HER2. Capitalizing on a highly predictable residue, we generate resistance mutations in TBK1, CSNK2A1 and BRAF. Unexpectedly, we uncover a potentially generalizable activation site that mediates drug resistance and confirm its impact in BRAF, EGFR, HER2 and MEK1. We anticipate that the identification of these residues will enable the broad interrogation of the kinome and its inhibitors.

High-grade serous ovarian cancer (HGSOC) is often sensitive to initial treatment with platinum and taxane combination chemotherapy, but most patients relapse with chemotherapy-resistant disease. To systematically identify genes modulating chemotherapy response, we performed pooled functional genomic screens in HGSOC cell lines treated with cisplatin, paclitaxel, or cisplatin plus paclitaxel. Genes in the intrinsic pathway of apoptosis were among the top candidate resistance genes in both gain-of-function and loss-of-function screens. In an open reading frame overexpression screen, followed by a mini-pool secondary screen, anti-apoptotic genes including BCL2L1 (BCL-XL) and BCL2L2 (BCL-W) were associated with chemotherapy resistance. In a CRISPR-Cas9 knockout screen, loss of BCL2L1 decreased cell survival whereas loss of proapoptotic genes promoted resistance. To dissect the role of individual anti-apoptotic proteins in HGSOC chemotherapy response, we evaluated overexpression or inhibition of BCL-2, BCL-XL, BCL-W, and MCL1 in HGSOC cell lines. Overexpression of anti-apoptotic proteins decreased apoptosis and modestly increased cell viability upon cisplatin or paclitaxel treatment. Conversely, specific inhibitors of BCL-XL, MCL1, or BCL-XL/BCL-2, but not BCL-2 alone, enhanced cell death when combined with cisplatin or paclitaxel. Anti-apoptotic protein inhibitors also sensitized HGSOC cells to the poly (ADP-ribose) polymerase inhibitor olaparib. These unbiased screens highlight anti-apoptotic proteins as mediators of chemotherapy resistance in HGSOC, and support inhibition of BCL-XL and MCL1, alone or combined with chemotherapy or targeted agents, in treatment of primary and recurrent HGSOC. IMPLICATIONS: Anti-apoptotic proteins modulate drug resistance in ovarian cancer, and inhibitors of BCL-XL or MCL1 promote cell death in combination with chemotherapy., (©2019 American Association for Cancer Research.)

The interleukin-3 receptor α subunit, CD123, is expressed in many hematologic malignancies including acute myeloid leukemia (AML) and blastic plasmacytoid dendritic cell neoplasm (BPDCN). Tagraxofusp (SL-401) is a CD123-targeted therapy consisting of interleukin-3 fused to a truncated diphtheria toxin payload. Factors influencing response to tagraxofusp other than CD123 expression are largely unknown. We interrogated tagraxofusp resistance in patients and experimental models and found that it was not associated with CD123 loss. Rather, resistant AML and BPDCN cells frequently acquired deficiencies in the diphthamide synthesis pathway, impairing tagraxofusp's ability to ADP-ribosylate cellular targets. Expression of DPH1, encoding a diphthamide pathway enzyme, was reduced by DNA CpG methylation in resistant cells. Treatment with the DNA methyltransferase inhibitor azacitidine restored DPH1 expression and tagraxofusp sensitivity. We also developed a drug-dependent ADP-ribosylation assay in primary cells that correlated with tagraxofusp activity and may represent an additional novel biomarker. As predicted by these results and our observation that resistance also increased mitochondrial apoptotic priming, we found that the combination of tagraxofusp and azacitidine was effective in patient-derived xenografts treated in vivo. These data have important implications for clinical use of tagraxofusp and led to a phase 1 study combining tagraxofusp and azacitidine in myeloid malignancies.

Mitochondrial apoptosis can be effectively targeted in lymphoid malignancies with the FDA-approved B cell lymphoma 2 (BCL-2) inhibitor venetoclax, but resistance to this agent is emerging. We show that venetoclax resistance in chronic lymphocytic leukemia is associated with complex clonal shifts. To identify determinants of resistance, we conducted parallel genome-scale screens of the BCL-2-driven OCI-Ly1 lymphoma cell line after venetoclax exposure along with integrated expression profiling and functional characterization of drug-resistant and engineered cell lines. We identified regulators of lymphoid transcription and cellular energy metabolism as drivers of venetoclax resistance in addition to the known involvement by BCL-2 family members, which were confirmed in patient samples. Our data support the implementation of combinatorial therapy with metabolic modulators to address venetoclax resistance., (Copyright © 2019 Elsevier Inc. All rights reserved.)

BET-bromodomain inhibition (BETi) has shown pre-clinical promise for MYC-amplified medulloblastoma. However, the mechanisms for its action, and ultimately for resistance, have not been fully defined. Here, using a combination of expression profiling, genome-scale CRISPR/Cas9-mediated loss of function and ORF/cDNA driven rescue screens, and cell-based models of spontaneous resistance, we identify bHLH/homeobox transcription factors and cell-cycle regulators as key genes mediating BETi's response and resistance. Cells that acquire drug tolerance exhibit a more neuronally differentiated cell-state and expression of lineage-specific bHLH/homeobox transcription factors. However, they do not terminally differentiate, maintain expression of CCND2, and continue to cycle through S-phase. Moreover, CDK4/CDK6 inhibition delays acquisition of resistance. Therefore, our data provide insights about the mechanisms underlying BETi effects and the appearance of resistance and support the therapeutic use of combined cell-cycle inhibitors with BETi in MYC-amplified medulloblastoma.

Large panels of comprehensively characterized human cancer models, including the Cancer Cell Line Encyclopedia (CCLE), have provided a rigorous framework with which to study genetic variants, candidate targets, and small-molecule and biological therapeutics and to identify new marker-driven cancer dependencies. To improve our understanding of the molecular features that contribute to cancer phenotypes, including drug responses, here we have expanded the characterizations of cancer cell lines to include genetic, RNA splicing, DNA methylation, histone H3 modification, microRNA expression and reverse-phase protein array data for 1,072 cell lines from individuals of various lineages and ethnicities. Integration of these data with functional characterizations such as drug-sensitivity, short hairpin RNA knockdown and CRISPR-Cas9 knockout data reveals potential targets for cancer drugs and associated biomarkers. Together, this dataset and an accompanying public data portal provide a resource for the acceleration of cancer research using model cancer cell lines.

Combinatorial inhibition of MEK1/2 and CDK4/6 is currently undergoing clinical investigation in NRAS-mutant melanoma. To prospectively map the landscape of resistance to this investigational regimen, we utilized a series of gain- and loss-of-function forward genetic screens to identify modulators of resistance to clinical inhibitors of MEK1/2 and CDK4/6 alone and in combination. First, we identified NRAS-mutant melanoma cell lines that were dependent on NRAS for proliferation and sensitive to MEK1/2 and CDK4/6 combination treatment. We then used a genome-scale ORF overexpression screen and a CRISPR knockout screen to identify modulators of resistance to each inhibitor alone or in combination. These orthogonal screening approaches revealed concordant means of achieving resistance to this therapeutic modality, including tyrosine kinases, RAF, RAS, AKT, and PI3K signaling. Activated KRAS was sufficient to cause resistance to combined MEK/CDK inhibition and to replace genetic depletion of oncogenic NRAS. In summary, our comprehensive functional genetic screening approach revealed modulation of resistance to the inhibition of MEK1/2, CDK4/6, or their combination in NRAS-mutant melanoma. SIGNIFICANCE: These findings reveal that NRAS-mutant melanomas can acquire resistance to genetic ablation of NRAS or combination MEK1/2 and CDK4/6 inhibition by upregulating activity of the RTK-RAS-RAF and RTK-PI3K-AKT signaling cascade., (©2019 American Association for Cancer Research.)

Unlike most tumor suppressor genes, the most common genetic alterations in tumor protein p53 (TP53) are missense mutations 1,2 . Mutant p53 protein is often abundantly expressed in cancers and specific allelic variants exhibit dominant-negative or gain-of-function activities in experimental models 3-8 . To gain a systematic view of p53 function, we interrogated loss-of-function screens conducted in hundreds of human cancer cell lines and performed TP53 saturation mutagenesis screens in an isogenic pair of TP53 wild-type and null cell lines. We found that loss or dominant-negative inhibition of wild-type p53 function reliably enhanced cellular fitness. By integrating these data with the Catalog of Somatic Mutations in Cancer (COSMIC) mutational signatures database 9,10 , we developed a statistical model that describes the TP53 mutational spectrum as a function of the baseline probability of acquiring each mutation and the fitness advantage conferred by attenuation of p53 activity. Collectively, these observations show that widely-acting and tissue-specific mutational processes combine with phenotypic selection to dictate the frequencies of recurrent TP53 mutations.

Despite advances in cancer biology and therapeutics, drug resistance remains problematic. Resistance is often multifactorial, heterogeneous, and prone to undersampling. Nonetheless, many individual mechanisms of targeted therapy resistance may coalesce into a smaller number of convergences, including pathway reactivation (downstream re-engagement of original effectors), pathway bypass (recruitment of a parallel pathway converging on the same downstream output), and pathway indifference (development of a cellular state independent of the initial therapeutic target). Similar convergences may also underpin immunotherapy resistance. Such parsimonious, convergence-based frameworks may help explain resistance across tumor types and therapeutic categories and may also suggest strategies to overcome it., (Copyright © 2018 Elsevier Inc. All rights reserved.)

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.

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., Competing Interests: The authors have declared that no competing interests exist.

Tumor-specific genomic information has the potential to guide therapeutic strategies and revolutionize patient treatment. Currently, this approach is limited by an abundance of disease-associated mutants whose biological functions and impacts on therapeutic response are uncharacterized. To begin to address this limitation, we functionally characterized nearly all (99.84%) missense mutants of MAPK1/ERK2, an essential effector of oncogenic RAS and RAF. Using this approach, we discovered rare gain- and loss-of-function ERK2 mutants found in human tumors, revealing that, in the context of this assay, mutational frequency alone cannot identify all functionally impactful mutants. Gain-of-function ERK2 mutants induced variable responses to RAF-, MEK-, and ERK-directed therapies, providing a reference for future treatment decisions. Tumor-associated mutations spatially clustered in two ERK2 effector-recruitment domains yet produced mutants with opposite phenotypes. This approach articulates an allele-characterization framework that can be scaled to meet the goals of genome-guided oncology., (Copyright © 2016 The Author(s). Published by Elsevier Inc. All rights reserved.)

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.

To explore the distinct genotypic and phenotypic states of melanoma tumors, we applied single-cell RNA sequencing (RNA-seq) to 4645 single cells isolated from 19 patients, profiling malignant, immune, stromal, and endothelial cells. Malignant cells within the same tumor displayed transcriptional heterogeneity associated with the cell cycle, spatial context, and a drug-resistance program. In particular, all tumors harbored malignant cells from two distinct transcriptional cell states, such that tumors characterized by high levels of the MITF transcription factor also contained cells with low MITF and elevated levels of the AXL kinase. Single-cell analyses suggested distinct tumor microenvironmental patterns, including cell-to-cell interactions. Analysis of tumor-infiltrating T cells revealed exhaustion programs, their connection to T cell activation and clonal expansion, and their variability across patients. Overall, we begin to unravel the cellular ecosystem of tumors and how single-cell genomics offers insights with implications for both targeted and immune therapies., (Copyright © 2016, American Association for the Advancement of Science.)

We conducted a large-scale functional genetic study to characterize mechanisms of resistance to ALK inhibition in ALK-dependent lung cancer cells. We identify members of known resistance pathways and additional putative resistance drivers. Among the latter were members of the P2Y purinergic receptor family of G-protein-coupled receptors (P2Y1, P2Y2, and P2Y6). P2Y receptors mediated resistance in part through a protein-kinase-C (PKC)-dependent mechanism. Moreover, PKC activation alone was sufficient to confer resistance to ALK inhibitors, whereas combined ALK and PKC inhibition restored sensitivity. We observed enrichment of gene signatures associated with several resistance drivers (including P2Y receptors) in crizotinib-resistant ALK-rearranged lung tumors compared to treatment-naive controls, supporting a role for these identified mechanisms in clinical ALK inhibitor resistance., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Over the past decade, tremendous progress in high-throughput small molecule-screening methods has facilitated the rapid expansion of phenotype-based data. Parallel advances in genomic characterization methods have complemented these efforts by providing a growing list of annotated cell line features. Together, these developments have paved the way for feature-based identification of novel, exploitable cellular dependencies, subsequently expanding our therapeutic toolkit in cancer and other diseases. Here, we provide an overview of the evolution of phenotypic small-molecule profiling and discuss the most significant and recent profiling and analytical efforts, their impact on the field, and their clinical ramifications. We additionally provide a perspective for future developments in phenotypic profiling efforts guided by genomic science., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

Unlabelled: Most melanomas harbor oncogenic BRAF(V600) mutations, which constitutively activate the MAPK pathway. Although MAPK pathway inhibitors show clinical benefit in BRAF(V600)-mutant melanoma, it remains incompletely understood why 10% to 20% of patients fail to respond. Here, we show that RAF inhibitor-sensitive and inhibitor-resistant BRAF(V600)-mutant melanomas display distinct transcriptional profiles. Whereas most drug-sensitive cell lines and patient biopsies showed high expression and activity of the melanocytic lineage transcription factor MITF, intrinsically resistant cell lines and biopsies displayed low MITF expression but higher levels of NF-κB signaling and the receptor tyrosine kinase AXL. In vitro, these MITF-low/NF-κB-high melanomas were resistant to inhibition of RAF and MEK, singly or in combination, and ERK. Moreover, in cell lines, NF-κB activation antagonized MITF expression and induced both resistance marker genes and drug resistance. Thus, distinct cell states characterized by MITF or NF-κB activity may influence intrinsic resistance to MAPK pathway inhibitors in BRAF(V600)-mutant melanoma., Significance: Although most BRAF(V600)-mutant melanomas are sensitive to RAF and/or MEK inhibitors, a subset fails to respond to such treatment. This study characterizes a transcriptional cell state distinction linked to MITF and NF-κB that may modulate intrinsic sensitivity of melanomas to MAPK pathway inhibitors., (©2014 American Association for Cancer Research.)

Treatment of BRAF-mutant melanoma with combined dabrafenib and trametinib, which target RAF and the downstream MAP-ERK kinase (MEK)1 and MEK2 kinases, respectively, improves progression-free survival and response rates compared with dabrafenib monotherapy. Mechanisms of clinical resistance to combined RAF/MEK inhibition are unknown. We performed whole-exome sequencing (WES) and whole-transcriptome sequencing (RNA-seq) on pretreatment and drug-resistant tumors from five patients with acquired resistance to dabrafenib/trametinib. In three of these patients, we identified additional mitogen-activated protein kinase (MAPK) pathway alterations in the resistant tumor that were not detected in the pretreatment tumor, including a novel activating mutation in MEK2 (MEK2(Q60P)). MEK2(Q60P) conferred resistance to combined RAF/MEK inhibition in vitro, but remained sensitive to inhibition of the downstream kinase extracellular signal-regulated kinase (ERK). The continued MAPK signaling-based resistance identified in these patients suggests that alternative dosing of current agents, more potent RAF/MEK inhibitors, and/or inhibition of the downstream kinase ERK may be needed for durable control of BRAF-mutant melanoma.

Most patients with BRAF(V600)-mutant metastatic melanoma develop resistance to selective RAF kinase inhibitors. The spectrum of clinical genetic resistance mechanisms to RAF inhibitors and options for salvage therapy are incompletely understood. We performed whole-exome sequencing on formalin-fixed, paraffin-embedded tumors from 45 patients with BRAF(V600)-mutant metastatic melanoma who received vemurafenib or dabrafenib monotherapy. Genetic alterations in known or putative RAF inhibitor resistance genes were observed in 23 of 45 patients (51%). Besides previously characterized alterations, we discovered a "long tail" of new mitogen-activated protein kinase (MAPK) pathway alterations (MAP2K2, MITF) that confer RAF inhibitor resistance. In three cases, multiple resistance gene alterations were observed within the same tumor biopsy. Overall, RAF inhibitor therapy leads to diverse clinical genetic resistance mechanisms, mostly involving MAPK pathway reactivation. Novel therapeutic combinations may be needed to achieve durable clinical control of BRAF(V600)-mutant melanoma. Integrating clinical genomics with preclinical screens may model subsequent resistance studies.

Malignant melanomas harbouring point mutations (Val600Glu) in the serine/threonine-protein kinase BRAF (BRAF(V600E)) depend on RAF-MEK-ERK signalling for tumour cell growth. RAF and MEK inhibitors show remarkable clinical efficacy in BRAF(V600E) melanoma; however, resistance to these agents remains a formidable challenge. Global characterization of resistance mechanisms may inform the development of more effective therapeutic combinations. Here we carried out systematic gain-of-function resistance studies by expressing more than 15,500 genes individually in a BRAF(V600E) melanoma cell line treated with RAF, MEK, ERK or combined RAF-MEK inhibitors. These studies revealed a cyclic-AMP-dependent melanocytic signalling 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 BRAF(V600E) melanoma patients revealed that phosphorylated (active) CREB was suppressed by RAF-MEK inhibition but restored in relapsing tumours. 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 MAPK-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 signalling- and chromatin-directed therapeutics.

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

A detailed understanding of the mechanisms by which tumors acquire resistance to targeted anticancer agents should speed the development of treatment strategies with lasting clinical efficacy. RAF inhibition in BRAF-mutant melanoma exemplifies the promise and challenge of many targeted drugs; although response rates are high, resistance invariably develops. Here, we articulate overarching principles of resistance to kinase inhibitors, as well as a translational approach to characterize resistance in the clinical setting through tumor mutation profiling. As a proof of principle, we performed targeted, massively parallel sequencing of 138 cancer genes in a tumor obtained from a patient with melanoma who developed resistance to PLX4032 after an initial dramatic response. The resulting profile identified an activating mutation at codon 121 in the downstream kinase MEK1 that was absent in the corresponding pretreatment tumor. The MEK1(C121S) mutation was shown to increase kinase activity and confer robust resistance to both RAF and MEK inhibition in vitro. Thus, MEK1(C121S) or functionally similar mutations are predicted to confer resistance to combined MEK/RAF inhibition. These results provide an instructive framework for assessing mechanisms of acquired resistance to kinase inhibition and illustrate the use of emerging technologies in a manner that may accelerate personalized cancer medicine.

Functional characterization of the human genome requires tools for systematically modulating gene expression in both loss-of-function and gain-of-function experiments. We describe the production of a sequence-confirmed, clonal collection of over 16,100 human open-reading frames (ORFs) encoded in a versatile Gateway vector system. Using this ORFeome resource, we created a genome-scale expression collection in a lentiviral vector, thereby enabling both targeted experiments and high-throughput screens in diverse cell types.

Oncogenic mutations in the serine/threonine kinase B-RAF (also known as BRAF) are found in 50-70% of malignant melanomas. Pre-clinical studies have demonstrated that the B-RAF(V600E) mutation predicts a dependency on the mitogen-activated protein kinase (MAPK) signalling cascade in melanoma-an observation that has been validated by the success of RAF and MEK inhibitors in clinical trials. However, clinical responses to targeted anticancer therapeutics are frequently confounded by de novo or acquired resistance. Identification of resistance mechanisms in a manner that elucidates alternative 'druggable' targets may inform effective long-term treatment strategies. Here we expressed ∼600 kinase and kinase-related open reading frames (ORFs) in parallel to interrogate resistance to a selective RAF kinase inhibitor. We identified MAP3K8 (the gene encoding COT/Tpl2) as a MAPK pathway agonist that drives resistance to RAF inhibition in B-RAF(V600E) cell lines. COT activates ERK primarily through MEK-dependent mechanisms that do not require RAF signalling. Moreover, COT expression is associated with de novo resistance in B-RAF(V600E) cultured cell lines and acquired resistance in melanoma cells and tissue obtained from relapsing patients following treatment with MEK or RAF inhibitors. We further identify combinatorial MAPK pathway inhibition or targeting of COT kinase activity as possible therapeutic strategies for reducing MAPK pathway activation in this setting. Together, these results provide new insights into resistance mechanisms involving the MAPK pathway and articulate an integrative approach through which high-throughput functional screens may inform the development of novel therapeutic strategies.

Inactivating mutations in NF1 underlie the prevalent familial cancer syndrome neurofibromatosis type 1 [1]. The NF1-encoded protein is a Ras GTPase-activating protein (RasGAP) [2]. Accordingly, Ras is aberrantly activated in NF1-deficient tumors; however, it is unknown which effector pathways critically function in tumor development. Here we provide in vivo evidence that TORC1/mTOR activity is essential for tumorigenesis. Specifically, we show that the mTOR inhibitor rapamycin potently suppresses the growth of aggressive NF1-associated malignancies in a genetically engineered murine model. However, in these tumors rapamycin does not function via mechanisms generally assumed to mediate tumor suppression, including inhibition of HIF-1alpha and indirect suppression of AKT, but does suppress the mTOR target Cyclin D1 [3]. These results demonstrate that mTOR inhibitors may be an effective targeted therapy for this commonly untreatable malignancy. Moreover, they indicate that mTOR inhibitors do not suppress all tumor types via the same mechanism, suggesting that current biomarkers that rely on HIF-1alpha suppression may not be informative for all cancers. Finally, our results reveal important differences between the effects of mTOR inhibition on the microvasculature in genetically engineered versus xenograft models and indicate that the former may be required for effective preclinical screening with this class of inhibitors.

Oncogene-induced senescence functions to limit tumor development. However, a complete understanding of the signals that trigger this type of senescence is currently lacking. We found that mutations affecting NF1, Raf, and Ras induce a global negative feedback response that potently suppresses Ras and/or its effectors. Moreover, these signals promote senescence by inhibiting the Ras/PI3K pathway, which can impact the senescence machinery through HDM2 and FOXO. This negative feedback program is regulated in part by RasGEFs, Sprouty proteins, RasGAPs, and MKPs. Moreover, these signals function in vivo in benign human tumors. Thus, the ultimate response to the aberrant activation of the Ras pathway is a multifaceted negative feedback signaling network that terminates the oncogenic signal and participates in the senescence response.

p63 is essential for normal epithelial development and is overexpressed in the vast majority of squamous cell carcinomas (SCC). Recent work had shown that DeltaNp63alpha is essential for survival of SCC cells, raising the possibility that the p63 pathway may be an attractive therapeutic target in these tumors. Nevertheless, it is unknown whether a therapeutic window exists for inhibiting p63 in tumor cells versus normal epithelia. Here, we show that SCC cells are uniquely dependent on DeltaNp63alpha for survival, unlike normal p63-expressing epithelial cells, and that dependence is mediated through tumor-specific up-regulation of the related protein p73. In normal primary human keratinocytes, we find that inhibition of endogenous p63 by RNA interference (RNAi) induces p21(CIP1) expression, inhibits cell cycle progression, and ultimately promotes cellular senescence. In contrast, p63 inhibition in SCC cells induces proapoptotic bcl-2 family members and rapidly triggers apoptosis. Expression of p73 is low in uncultured basal keratinocytes but is markedly up-regulated in both SCC cell lines and primary tumors in vivo. Whereas p21(CIP1) induction following loss of p63 in normal cells is independent of p53 and p73, both proapoptotic gene induction and cell death following p63 RNAi in tumor cells are p73 dependent. Finally, ectopic p73 expression in primary keratinocytes does not affect baseline cell proliferation but is sufficient to trigger cell death following loss of p63. Together, these findings define a specific molecular mechanism of p63 dependence through p73 up-regulation, and they provide a rationale for targeting the p63 pathway as a therapeutic strategy in SCCs.

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.

Loss-of-function mutations in the NF1 tumor suppressor gene underlie the familial cancer syndrome neurofibromatosis type I (NF1). The NF1-encoded protein, neurofibromin, functions as a Ras-GTPase activating protein (RasGAP). Accordingly, deregulation of Ras is thought to contribute to NF1 development. However, the critical effector pathways involved in disease pathogenesis are still unknown. We show here that the mTOR pathway is tightly regulated by neurofibromin. mTOR is constitutively activated in both NF1-deficient primary cells and human tumors in the absence of growth factors. This aberrant activation depends on Ras and PI3 kinase, and is mediated by the phosphorylation and inactivation of the TSC2-encoded protein tuberin by AKT. Importantly, tumor cell lines derived from NF1 patients, and a genetically engineered cell system that requires Nf1-deficiency for transformation, are highly sensitive to the mTOR inhibitor rapamycin. Furthermore, while we show that the activation of endogenous Ras leads to constitutive mTOR signaling in this disease state, we also demonstrate that in normal cells Ras is differentially required for mTOR signaling in response to various growth factors. Thus, these findings identify the NF1 tumor suppressor as an indispensable regulator of TSC2 and mTOR. Furthermore, our results also demonstrate that Ras plays a critical role in the activation of mTOR in both normal and tumorigenic settings. Finally, these data suggest that rapamycin, or its derivatives, may represent a viable therapy for NF1.

We identified REDD1 as a novel transcriptional target of p53 induced following DNA damage. During embryogenesis, REDD1 expression mirrors the tissue-specific pattern of the p53 family member p63, and TP63 null embryos show virtually no expression of REDD1, which is restored in mouse embryo fibroblasts following p63 expression. In differentiating primary keratinocytes, TP63 and REDD1 expression are coordinately downregulated, and ectopic expression of either gene inhibits in vitro differentiation. REDD1 appears to function in the regulation of reactive oxygen species (ROS); we show that TP63 null fibroblasts have decreased ROS levels and reduced sensitivity to oxidative stress, which are both increased following ectopic expression of either TP63 or REDD1. Thus, REDD1 encodes a shared transcriptional target that implicates ROS in the p53-dependent DNA damage response and in p63-mediated regulation of epithelial differentiation.