AP-1 Mediates Cellular Adaptation and Memory Formation During Therapy Resistance.

Cellular responses to environmental stimuli are typically thought to be governed by genetically encoded programs. We demonstrate that melanoma cells can form and maintain cellular memories during the acquisition of therapy resistance that exhibit characteristics of cellular learning and are dependent on the transcription factor AP-1. We show that cells exposed to a low dose of therapy adapt to become resistant to a high dose, demonstrating that resistance was not purely selective. The application of therapy itself results in the encoding of transient gene expression into cellular memory and that this encoding occurs for both transiently induced and probabilistically arising expression. Chromatin accessibility showed concomitant persistence. A two-color AP-1 reporter system showed that these memories are encoded in cis, constituting an example of activating cis epigenetics. Our findings establish the formation and maintenance of cellular memories as a critical aspect of gene regulation during the development of therapy resistance.

ERK hyperactivation serves as a unified mechanism of escape in intrinsic and acquired CDK4/6 inhibitor resistance in acral lentiginous melanoma.

Patients with metastatic acral lentiginous melanoma (ALM) suffer worse outcomes relative to patients with other forms of cutaneous melanoma (CM), and do not benefit as well to approved melanoma therapies. Identification of cyclin-dependent kinase 4 and 6 (CDK4/6) pathway gene alterations in >60% of ALMs has led to clinical trials of the CDK4/6 inhibitor (CDK4i/6i) palbociclib for ALM; however, median progression free survival with CDK4i/6i treatment was only 2.2 months, suggesting existence of resistance mechanisms. Therapy resistance in ALM remains poorly understood; here we report hyperactivation of MAPK signaling and elevated cyclin D1 expression serve as a mechanism of intrinsic early/adaptive CDK4i/6i resistance. ALM cells that have acquired CDK4i/6i resistance following chronic treatment exposure also exhibit hyperactivation of the MAPK pathway. MEK and/or ERK inhibition increases CDK4i/6i efficacy against therapy naïve and CDK4i/6i-resistant AM cells in xenograft and patient-derived xenograft (PDX) models and promotes a defective DNA repair, cell cycle arrested and apoptotic program. Notably, gene alterations poorly correlate with protein expression of cell cycle proteins in ALM or efficacy of CDK4i/6i, urging additional strategies when stratifying patients for CDK4i/6i trial inclusion. Concurrent targeting of the MAPK pathway and CDK4/6 represents a new approach for patients with metastatic ALM to improve outcomes.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity.

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

Diverse clonal fates emerge upon drug treatment of homogeneous cancer cells.

Even among genetically identical cancer cells, resistance to therapy frequently emerges from a small subset of those cells1-7. Molecular differences in rare individual cells in the initial population enable certain cells to become resistant to therapy7-9; however, comparatively little is known about the variability in the resistance outcomes. Here we develop and apply FateMap, a framework that combines DNA barcoding with single-cell RNA sequencing, to reveal the fates of hundreds of thousands of clones exposed to anti-cancer therapies. We show that resistant clones emerging from single-cell-derived cancer cells adopt molecularly, morphologically and functionally distinct resistant types. These resistant types are largely predetermined by molecular differences between cells before drug addition and not by extrinsic factors. Changes in the dose and type of drug can switch the resistant type of an initial cell, resulting in the generation and elimination of certain resistant types. Samples from patients show evidence for the existence of these resistant types in a clinical context. We observed diversity in resistant types across several single-cell-derived cancer cell lines and cell types treated with a variety of drugs. The diversity of resistant types as a result of the variability in intrinsic cell states may be a generic feature of responses to external cues.

Toward a systems-level probing of tumor clonality.

Cancer has been described as a genetic disease that clonally evolves in the face of selective pressures imposed by cell-intrinsic and extrinsic factors. Although classical models based on genetic data predominantly propose Darwinian mechanisms of cancer evolution, recent single-cell profiling of cancers has described unprecedented heterogeneity in tumors providing support for alternative models of branched and neutral evolution through both genetic and non-genetic mechanisms. Emerging evidence points to a complex interplay between genetic, non-genetic, and extrinsic environmental factors in shaping the evolution of tumors. In this perspective, we briefly discuss the role of cell-intrinsic and extrinsic factors that shape clonal behaviors during tumor progression, metastasis, and drug resistance. Taking examples of pre-malignant states associated with hematological malignancies and esophageal cancer, we discuss recent paradigms of tumor evolution and prospective approaches to further enhance our understanding of this spatiotemporally regulated process.

ERK Hyperactivation Serves as a Unified Mechanism of Escape in Intrinsic and Acquired CDK4/6 Inhibitor Resistance in Acral Lentiginous Melanoma.

Patients with metastatic acral lentiginous melanoma (ALM) suffer worse outcomes relative to patients with other forms of cutaneous melanoma (CM), and do not benefit as well to approved melanoma therapies. Identification of cyclin-dependent kinase 4 and 6 (CDK4/6) pathway gene alterations in > 60% of ALMs has led to clinical trials of the CDK4/6 inhibitor (CDK4i/6i) palbociclib for ALM; however, median progression free survival with CDK4i/6i treatment was only 2.2 months, suggesting existence of resistance mechanisms. Therapy resistance in ALM remains poorly understood; here we report hyperactivation of MAPK signaling and elevated cyclin D1 expression are a unified mechanism of both intrinsic and acquired CDK4i/6i resistance. MEK and/or ERK inhibition increases CDK4i/6i efficacy in a patient-derived xenograft (PDX) model of ALM and promotes a defective DNA repair, cell cycle arrested and apoptotic program. Notably, gene alterations poorly correlate with protein expression of cell cycle proteins in ALM or efficacy of CDK4i/6i, urging additional strategies when stratifying patients for CDK4i/6i trial inclusion. Concurrent targeting of the MAPK pathway and CDK4/6 represents a new approach to improve outcomes for patients with advanced ALM.

Unraveling non-genetic heterogeneity in cancer with dynamical models and computational tools.

Individual cells within an otherwise genetically homogenous population constantly undergo fluctuations in their molecular state, giving rise to non-genetic heterogeneity. Such diversity is being increasingly implicated in cancer therapy resistance and metastasis. Identifying the origins of non-genetic heterogeneity is therefore crucial for making clinical breakthroughs. We discuss with examples how dynamical models and computational tools have provided critical multiscale insights into the nature and consequences of non-genetic heterogeneity in cancer. We demonstrate how mechanistic modeling has been pivotal in establishing key concepts underlying non-genetic diversity at various biological scales, from population dynamics to gene regulatory networks. We discuss advances in single-cell longitudinal profiling techniques to reveal patterns of non-genetic heterogeneity, highlighting the ongoing efforts and challenges in statistical frameworks to robustly interpret such multimodal datasets. Moving forward, we stress the need for data-driven statistical and mechanistically motivated dynamical frameworks to come together to develop predictive cancer models and inform therapeutic strategies.

Changes in chromatin accessibility are not concordant with transcriptional changes for single-factor perturbations.

A major goal in the field of transcriptional regulation is the mapping of changes in the binding of transcription factors to the resultant changes in gene expression. Recently, methods for measuring chromatin accessibility have enabled us to measure changes in accessibility across the genome, which are thought to correspond to transcription factor-binding events. In concert with RNA-sequencing, these data in principle enable such mappings; however, few studies have looked at their concordance over short-duration treatments with specific perturbations. Here, we used tandem, bulk ATAC-seq, and RNA-seq measurements from MCF-7 breast carcinoma cells to systematically evaluate the concordance between changes in accessibility and changes in expression in response to retinoic acid and TGF-ß. We found two classes of genes whose expression showed a significant change: those that showed some changes in the accessibility of nearby chromatin, and those that showed virtually no change despite strong changes in expression. The peaks associated with genes in the former group had lower baseline accessibility prior to exposure to signal. Focusing the analysis specifically on peaks with motifs for transcription factors associated with retinoic acid and TGF-ß signaling did not reduce the lack of correspondence. Analysis of paired chromatin accessibility and gene expression data from distinct paths along the hematopoietic differentiation trajectory showed a much stronger correspondence, suggesting that the multifactorial biological processes associated with differentiation may lead to changes in chromatin accessibility that reflect rather than driving altered transcriptional status. Together, these results show many gene expression changes can happen independently of changes in the accessibility of local chromatin in the context of a single-factor perturbation.

Gene Networks with Transcriptional Bursting Recapitulate Rare Transient Coordinated High Expression States in Cancer.

Non-genetic transcriptional variability is a potential mechanism for therapy resistance in melanoma. Specifically, rare subpopulations of cells occupy a transient pre-resistant state characterized by coordinated high expression of several genes and survive therapy. How might these rare states arise and disappear within the population? It is unclear whether the canonical models of probabilistic transcriptional pulsing can explain this behavior, or if it requires special, hitherto unidentified mechanisms. We show that a minimal model of transcriptional bursting and gene interactions can give rise to rare coordinated high expression states. These states occur more frequently in networks with low connectivity and depend on three parameters. While entry into these states is initiated by a long transcriptional burst that also triggers entry of other genes, the exit occurs through independent inactivation of individual genes. Together, we demonstrate that established principles of gene regulation are sufficient to describe this behavior and argue for its more general existence. A record of this paper's transparent peer review process is included in the Supplemental Information.

A quantitative model of developmental RTK signaling.

Receptor tyrosine kinases (RTKs) control a wide range of developmental processes, from the first stages of embryogenesis to postnatal growth and neurocognitive development in the adult. A significant share of our knowledge about RTKs comes from genetic screens in model organisms, which provided numerous examples demonstrating how specific cell fates and morphologies are abolished when RTK activation is either abrogated or significantly reduced. Aberrant activation of such pathways has also been recognized in many forms of cancer. More recently, studies of human developmental syndromes established that excessive activation of RTKs and their downstream signaling effectors, most notably the Ras signaling pathway, can also lead to structural and functional defects. Given that both insufficient and excessive pathway activation can lead to abnormalities, mechanistic analysis of developmental RTK signaling must address quantitative questions about its regulation and function. Patterning events controlled by the RTK Torso in the early Drosophila embryo are well-suited for this purpose. This mini review summarizes current state of knowledge about Torso-dependent Ras activation and discusses its potential to serve as a quantitative model for studying the general principles of Ras signaling in development and disease.

How activating mutations affect MEK1 regulation and function.

The MEK1 kinase directly phosphorylates ERK2, after the activation loop of MEK1 is itself phosphorylated by Raf. Studies over the past decade have revealed a large number of disease-related mutations in the MEK1 gene that lead to tumorigenesis and abnormal development. Several of these mutations result in MEK1 constitutive activity, but how they affect MEK1 regulation and function remains largely unknown. Here, we address these questions focusing on two pathogenic variants of the Phe-53 residue, which maps to the well-characterized negative regulatory region of MEK1. We found that these variants are phosphorylated by Raf faster than the wild-type enzyme, and this phosphorylation further increases their enzymatic activity. However, the maximal activities of fully phosphorylated wild-type and mutant enzymes are indistinguishable. On the basis of available structural information, we propose that the activating substitutions destabilize the inactive conformation of MEK1, resulting in its constitutive activity and making it more prone to Raf-mediated phosphorylation. Experiments in zebrafish revealed that the effects of activating variants on embryonic development reflect the joint control of the negative regulatory region and activating phosphorylation. Our results underscore the complexity of the effects of activating mutations on signaling systems, even at the level of a single protein.

Silencing of retrotransposons by SETDB1 inhibits the interferon response in acute myeloid leukemia.

A propensity for rewiring genetic and epigenetic regulatory networks, thus enabling sustained cell proliferation, suppression of apoptosis, and the ability to evade the immune system, is vital to cancer cell propagation. An increased understanding of how this is achieved is critical for identifying or improving therapeutic interventions. In this study, using acute myeloid leukemia (AML) human cell lines and a custom CRISPR/Cas9 screening platform, we identify the H3K9 methyltransferase SETDB1 as a novel, negative regulator of innate immunity. SETDB1 is overexpressed in many cancers, and loss of this gene in AML cells triggers desilencing of retrotransposable elements that leads to the production of double-stranded RNAs (dsRNAs). This is coincident with induction of a type I interferon response and apoptosis through the dsRNA-sensing pathway. Collectively, our findings establish a unique gene regulatory axis that cancer cells can exploit to circumvent the immune system.

Parallel imaging of Drosophila embryos for quantitative analysis of genetic perturbations of the Ras pathway.

The Ras pathway patterns the poles of the Drosophila embryo by downregulating the levels and activity of a DNA-binding transcriptional repressor Capicua (Cic). We demonstrate that the spatiotemporal pattern of Cic during this signaling event can be harnessed for functional studies of mutations in the Ras pathway in human diseases. Our approach relies on a new microfluidic device that enables parallel imaging of Cic dynamics in dozens of live embryos. We found that although the pattern of Cic in early embryos is complex, it can be accurately approximated by a product of one spatial profile and one time-dependent amplitude. Analysis of these functions of space and time alone reveals the differential effects of mutations within the Ras pathway. Given the highly conserved nature of Ras-dependent control of Cic, our approach provides new opportunities for functional analysis of multiple sequence variants from developmental abnormalities and cancers.

Uncoupling neurogenic gene networks in the Drosophila embryo.

The EGF signaling pathway specifies neuronal identities in the Drosophila embryo by regulating developmental patterning genes such as intermediate neuroblasts defective (ind). EGFR is activated in the ventral midline and neurogenic ectoderm by the Spitz ligand, which is processed by the Rhomboid protease. CRISPR/Cas9 was used to delete defined rhomboid enhancers mediating expression at each site of Spitz processing. Surprisingly, the neurogenic ectoderm, not the ventral midline, was found to be the dominant source of EGF patterning activity. We suggest that Drosophila is undergoing an evolutionary transition in central nervous system (CNS)-organizing activity from the ventral midline to the neurogenic ectoderm.

Divergent effects of intrinsically active MEK variants on developmental Ras signaling.

Germline mutations in Ras pathway components are associated with a large class of human developmental abnormalities, known as RASopathies, that are characterized by a range of structural and functional phenotypes, including cardiac defects and neurocognitive delays. Although it is generally believed that RASopathies are caused by altered levels of pathway activation, the signaling changes in developing tissues remain largely unknown. We used assays with spatiotemporal resolution in Drosophila melanogaster (fruit fly) and Danio rerio (zebrafish) to quantify signaling changes caused by mutations in MAP2K1 (encoding MEK), a core component of the Ras pathway that is mutated in both RASopathies and cancers in humans. Surprisingly, we discovered that intrinsically active MEK variants can both increase and reduce the levels of pathway activation in vivo. The sign of the effect depends on cellular context, implying that some of the emerging phenotypes in RASopathies may be caused by increased, as well as attenuated, levels of Ras signaling.

In vivo severity ranking of Ras pathway mutations associated with developmental disorders.

Germ-line mutations in components of the Ras/MAPK pathway result in developmental disorders called RASopathies, affecting about 1/1,000 human births. Rapid advances in genome sequencing make it possible to identify multiple disease-related mutations, but there is currently no systematic framework for translating this information into patient-specific predictions of disease progression. As a first step toward addressing this issue, we developed a quantitative, inexpensive, and rapid framework that relies on the early zebrafish embryo to assess mutational effects on a common scale. Using this assay, we assessed 16 mutations reported in MEK1, a MAPK kinase, and provide a robust ranking of these mutations. We find that mutations found in cancer are more severe than those found in both RASopathies and cancer, which, in turn, are generally more severe than those found only in RASopathies. Moreover, this rank is conserved in other zebrafish embryonic assays and Drosophila-specific embryonic and adult assays, suggesting that our ranking reflects the intrinsic property of the mutant molecule. Furthermore, this rank is predictive of the drug dose needed to correct the defects. This assay can be readily used to test the strengths of existing and newly found mutations in MEK1 and other pathway components, providing the first step in the development of rational guidelines for patient-specific diagnostics and treatment of RASopathies.

Role of dopant concentration, crystal phase and particle size on microbial inactivation of Cu-doped TiO2 nanoparticles.

The properties of Cu-doped TiO(2) nanoparticles (NPs) were independently controlled in a flame aerosol reactor by varying the molar feed ratios of the precursors, and by optimizing temperature and time history in the flame. The effect of the physico-chemical properties (dopant concentration, crystal phase and particle size) of Cu-doped TiO(2) nanoparticles on inactivation of Mycobacterium smegmatis (a model pathogenic bacterium) was investigated under three light conditions (complete dark, fluorescent light and UV light). The survival rate of M. smegmatis (in a minimal salt medium for 2 h) exposed to the NPs varied depending on the light irradiation conditions as well as the dopant concentrations. In dark conditions, pristine TiO(2) showed insignificant microbial inactivation, but inactivation increased with increasing dopant concentration. Under fluorescent light illumination, no significant effect was observed for TiO(2). However, when TiO(2) was doped with copper, inactivation increased with dopant concentration, reaching more than 90% (>3 wt% dopant). Enhanced microbial inactivation by TiO(2) NPs was observed only under UV light. When TiO(2) NPs were doped with copper, their inactivation potential was promoted and the UV-resistant cells were reduced by over 99%. In addition, the microbial inactivation potential of NPs was also crystal-phase-and size-dependent under all three light conditions. A lower ratio of anatase phase and smaller sizes of Cu-doped TiO(2) NPs resulted in decreased bacterial survival. The increased inactivation potential of doped TiO(2) NPs is possibly due to both enhanced photocatalytic reactions and leached copper ions.