An immediate-late gene expression module decodes ERK signal duration

Frequency and duration of sensory flicker controls astrocyte and neuron specific transcriptional profiles in 5xFAD mice

Background

Current clinical trials are investigating gamma frequency sensory stimulation as a potential therapeutic strategy for Alzheimer's disease, yet we lack a comprehensive picture of the effects of this stimulation on multiple aspects of brain function. While most prior research has focused on gamma frequency sensory stimulation, we previously showed that exposing mice to visual flickering stimulation increased MAPK and NFκB signaling in the visual cortex in a manner dependent on duration and frequency of sensory stimulation exposure. Because these pathways control multiple neuronal and glial functions and are differentially activated based on the duration and frequency of flicker stimulation, we aimed to define the transcriptional effects of different frequencies and durations of flicker stimulation on multiple brain functions.

Methods

We exposed 5xFAD mice to different frequencies of audio/visual flicker stimulation (constant light, 10Hz, 20Hz, 40Hz) for durations of 0.5hr, 1hr, or 4hr, then used bulk RNAseq to profile transcriptional changes within the visual cortex and hippocampus tissues. Using weighted gene co-expression network analysis, we identified modules of co-expressed genes controlled by frequency and/or duration of stimulation.

Results

Within the visual cortex, we found that all stimulation frequencies caused fast activation of a module of immune genes within 1hr and slower suppression of synaptic genes after 4hrs of stimulation. Interestingly, all frequencies of stimulation led to slow suppression of astrocyte specific gene sets, while activation of neuronal gene sets was frequency and duration specific. In contrast, in the hippocampus, immune and synaptic modules were suppressed based on the frequency of stimulation. Specifically,10Hz activated a module of genes associated with mitochondrial function, metabolism, and synaptic translation while 10Hz rapidly suppressed a module of genes linked to neurotransmitter activity.

Conclusion

Collectively, our data indicate that the frequency and duration of flicker stimulation controls immune, neuronal, and metabolic genes in multiple regions of the brain affected by Alzheimer's disease. Flicker stimulation may thus represent a potential therapeutic strategy that can be tuned based on the brain region and the specific cellular process to be modulated.

Deciphering the History of ERK Activity from Fixed-Cell Immunofluorescence Measurements

The Ras/ERK pathway drives cell proliferation and other oncogenic behaviors, and quantifying its activity in situ is of high interest in cancer diagnosis and therapy. Pathway activation is often assayed by measuring phosphorylated ERK. However, this form of measurement overlooks dynamic aspects of signaling that can only be observed over time. In this study, we combine a live, single-cell ERK biosensor approach with multiplexed immunofluorescence staining of downstream target proteins to ask how well immunostaining captures the dynamic history of ERK activity. Combining linear regression, machine learning, and differential equation models, we develop an interpretive framework for immunostains, in which Fra-1 and pRb levels imply long term activation of ERK signaling, while Egr-1 and c-Myc indicate recent activation. We show that this framework can distinguish different classes of ERK dynamics within a heterogeneous population, providing a tool for annotating ERK dynamics within fixed tissues.

Cis-regulatory control of transcriptional timing and noise in response to estrogen

Cis-regulatory elements control transcription levels, temporal dynamics, and cell-cell variation or transcriptional noise. However, the combination of regulatory features that control these different attributes is not fully understood. Here, we used single cell RNA-seq during an estrogen treatment time course and machine learning to identify predictors of expression timing and noise. We find that genes with multiple active enhancers exhibit faster temporal responses. We verified this finding by showing that manipulation of enhancer activity changes the temporal response of estrogen target genes. Analysis of transcriptional noise uncovered a relationship between promoter and enhancer activity, with active promoters associated with low noise and active enhancers linked to high noise. Finally, we observed that co-expression across single cells is an emergent property associated with chromatin looping, timing, and noise. Overall, our results indicate a fundamental tradeoff between a gene's ability to quickly respond to incoming signals and maintain low variation across cells.

A guide to ERK dynamics, part 2: downstream decoding

Signaling by the extracellular signal-regulated kinase (ERK) pathway controls many cellular processes, including cell division, death, and differentiation. In this second installment of a two-part review, we address the question of how the ERK pathway exerts distinct and context-specific effects on multiple processes. We discuss how the dynamics of ERK activity induce selective changes in gene expression programs, with insights from both experiments and computational models. With a focus on single-cell biosensor-based studies, we summarize four major functional modes for ERK signaling in tissues: adjusting the size of cell populations, gradient-based patterning, wave propagation of morphological changes, and diversification of cellular gene expression states. These modes of operation are disrupted in cancer and other related diseases and represent potential targets for therapeutic intervention. By understanding the dynamic mechanisms involved in ERK signaling, there is potential for pharmacological strategies that not only simply inhibit ERK, but also restore functional activity patterns and improve disease outcomes.

Regulation of ERK2 activity by dynamic S-acylation

Extracellular signal-regulated kinases (ERK1/2) are key effector proteins of the mitogen-activated protein kinase pathway, choreographing essential processes of cellular physiology. Here, we discover that ERK1/2 are subject to S-acylation, a reversible lipid modification of cysteine residues, at C271/C254. The levels of ERK1/2 S-acylation are modulated by epidermal growth factor (EGF) signaling, mirroring its phosphorylation dynamics, and acylation-deficient ERK2 displays altered phosphorylation patterns. We show that ERK1/2 S-acylation is mediated by "writer" protein acyl transferases (PATs) and "eraser" acyl protein thioesterases (APTs) and that chemical inhibition of either lipid addition or removal alters ERK1/2's EGF-triggered transcriptional program. Finally, in a mouse model of metabolic syndrome, we find that ERK1/2 lipidation levels correlate with alterations in ERK1/2 lipidation writer/eraser expression, solidifying a link between ERK1/2 activity, ERK1/2 lipidation, and organismal health. This study describes how lipidation regulates ERK1/2 and offers insight into the role of dynamic S-acylation in cell signaling more broadly.

RAS-dependent RAF-MAPK hyperactivation by pathogenic RIT1 is a therapeutic target in Noonan syndrome-associated cardiac hypertrophy

RIT1 is a RAS guanosine triphosphatase (GTPase) that regulates different aspects of signal transduction and is mutated in lung cancer, leukemia, and in the germline of individuals with Noonan syndrome. Pathogenic RIT1 proteins promote mitogen-activated protein kinase (MAPK) hyperactivation; however, this mechanism remains poorly understood. Here, we show that RAF kinases are direct effectors of membrane-bound mutant RIT1 necessary for MAPK activation. We identify critical residues in RIT1 that facilitate interaction with membrane lipids and show that these are necessary for association with RAF kinases and MAPK activation. Although mutant RIT1 binds to RAF kinases directly, it fails to activate MAPK signaling in the absence of classical RAS proteins. Consistent with aberrant RAF/MAPK activation as a driver of disease, we show that pathway inhibition alleviates cardiac hypertrophy in a mouse model of RIT1 mutant Noonan syndrome. These data shed light on the function of pathogenic RIT1 and identify avenues for therapeutic intervention.

KRAS(G12D) drives lepidic adenocarcinoma through stem-cell reprogramming

Many cancers originate from stem or progenitor cells hijacked by somatic mutations that drive replication, exemplified by adenomatous transformation of pulmonary alveolar epithelial type II (AT2) cells1. Here we demonstrate a different scenario: expression of KRAS(G12D) in differentiated AT1 cells reprograms them slowly and asynchronously back into AT2 stem cells that go on to generate indolent tumours. Like human lepidic adenocarcinoma, the tumour cells slowly spread along alveolar walls in a non-destructive manner and have low ERK activity. We find that AT1 and AT2 cells act as distinct cells of origin and manifest divergent responses to concomitant WNT activation and KRAS(G12D) induction, which accelerates AT2-derived but inhibits AT1-derived adenoma proliferation. Augmentation of ERK activity in KRAS(G12D)-induced AT1 cells increases transformation efficiency, proliferation and progression from lepidic to mixed tumour histology. Overall, we have identified a new cell of origin for lung adenocarcinoma, the AT1 cell, which recapitulates features of human lepidic cancer. In so doing, we also uncover a capacity for oncogenic KRAS to reprogram a differentiated and quiescent cell back into its parent stem cell en route to adenomatous transformation. Our work further reveals that irrespective of a given cancer's current molecular profile and driver oncogene, the cell of origin exerts a pervasive and perduring influence on its subsequent behaviour.

ABL1/2 and DDR1 Drive MEKi Resistance in NRAS-Mutant Melanomas by Stabilizing RAF/MYC/ETS1 and Promoting RAF Homodimerization

Melanomas harboring NRAS mutations are a particularly aggressive and deadly subtype. If patients cannot tolerate or the melanomas are insensitive to immune checkpoint blockade, there are no effective 2nd-line treatment options. Drugs targeting the RAF/MEK/ERK pathway, which are used for BRAF-mutant melanomas, do little to increase progression-free survival (PFS). Here, using both loss-of-function and gain-of-function approaches, we show that ABL1/2 and DDR1 are critical nodes during NRAS-mutant melanoma intrinsic and acquired MEK inhibitor (MEKi) resistance. In some acquired resistance cells, ABL1/2 and DDR1 cooperate to stabilize RAF proteins, activate ERK cytoplasmic and nuclear signaling, repress p27/KIP1 expression, and drive RAF homodimerization. In contrast, other acquired resistance cells depend solely on ABL1/2 for their survival, and are sensitive to highly specific allosteric ABL1/2 inhibitors, which prevent β-catenin nuclear localization and destabilize MYC and ETS1 in an ERK-independent manner. Significantly, targeting ABL1/2 and DDR1 with an FDA-approved anti-leukemic drug, reverses intrinsic MEKi resistance, delays acquisition of acquired resistance, and doubles the survival time in a NRAS-mutant mouse model. These data indicate that repurposing FDA-approved drugs targeting ABL1/2 and DDR1 may be a novel and effective strategy for treating patients with treatment-refractory NRAS-driven melanomas.

Dusp1 regulates thermal tolerance limits in zebrafish by maintaining mitochondrial integrity

Temperature tolerance restricts the distribution of a species. However, the molecular and cellular mechanisms that set the thermal tolerance limits of an organism are poorly understood. Here, we report on the function of dual-specificity phosphatase 1 (DUSP1) in thermal tolerance regulation. Notably, we found that dusp1 ^-/- zebrafish grew normally but survived within a narrowed temperature range. The higher susceptibility of these mutant fish to both cold and heat challenges was attributed to accelerated cell death caused by aggravated mitochondrial dysfunction and over-production of reactive oxygen species in the gills. The DUSP1-MAPK-DRP1 axis was identified as a key pathway regulating these processes in both fish and human cells. These observations suggest that DUSP1 may play a role in maintaining mitochondrial integrity and redox homeostasis. We therefore propose that maintenance of cellular redox homeostasis may be a key mechanism for coping with cellular thermal stress and that the interplay between signaling pathways regulating redox homeostasis in the most thermosensitive tissue (i.e., gills) may play an important role in setting the thermal tolerance limit of zebrafish.

Early response evaluation by single cell signaling profiling in acute myeloid leukemia

Aberrant pro-survival signaling is a hallmark of cancer cells, but the response to chemotherapy is poorly understood. In this study, we investigate the initial signaling response to standard induction chemotherapy in a cohort of 32 acute myeloid leukemia (AML) patients, using 36-dimensional mass cytometry. Through supervised and unsupervised machine learning approaches, we find that reduction of extracellular-signal-regulated kinase (ERK) 1/2 and p38 mitogen-activated protein kinase (MAPK) phosphorylation in the myeloid cell compartment 24 h post-chemotherapy is a significant predictor of patient 5-year overall survival in this cohort. Validation by RNA sequencing shows induction of MAPK target gene expression in patients with high phospho-ERK1/2 24 h post-chemotherapy, while proteomics confirm an increase of the p38 prime target MAPK activated protein kinase 2 (MAPKAPK2). In this study, we demonstrate that mass cytometry can be a valuable tool for early response evaluation in AML and elucidate the potential of functional signaling analyses in precision oncology diagnostics.

ERK-mediated NELF-A phosphorylation promotes transcription elongation of immediate-early genes by releasing promoter-proximal pausing of RNA polymerase II

Growth factor-induced, ERK-mediated induction of immediate-early genes (IEGs) is crucial for cell growth and tumorigenesis. Although IEG expression is mainly regulated at the level of transcription elongation by RNA polymerase-II (Pol-II) promoter-proximal pausing and its release, the role of ERK in this process remains unknown. Here, we identified negative elongation factor (NELF)-A as an ERK substrate. Upon growth factor stimulation, ERK phosphorylates NELF-A, which dissociates NELF from paused Pol-II at the promoter-proximal regions of IEGs, allowing Pol-II to resume elongation and produce full-length transcripts. Furthermore, we found that in cancer cells, PP2A efficiently dephosphorylates NELF-A, thereby preventing aberrant IEG expression induced by ERK-activating oncogenes. However, when PP2A inhibitor proteins are overexpressed, as is frequently observed in cancers, decreased PP2A activity combined with oncogene-mediated ERK activation conspire to induce NELF-A phosphorylation and IEG upregulation, resulting in tumor progression. Our data delineate previously unexplored roles of ERK and PP2A inhibitor proteins in carcinogenesis.

Relating individual cell division events to single-cell ERK and Akt activity time courses

Biochemical correlates of stochastic single-cell fates have been elusive, even for the well-studied mammalian cell cycle. We monitored single-cell dynamics of the ERK and Akt pathways, critical cell cycle progression hubs and anti-cancer drug targets, and paired them to division events in the same single cells using the non-transformed MCF10A epithelial line. Following growth factor treatment, in cells that divide both ERK and Akt activities are significantly higher within the S-G2 time window (~ 8.5-40 h). Such differences were much smaller in the pre-S-phase, restriction point window which is traditionally associated with ERK and Akt activity dependence, suggesting unappreciated roles for ERK and Akt in S through G2. Simple metrics of central tendency in this time window are associated with subsequent cell division fates. ERK activity was more strongly associated with division fates than Akt activity, suggesting Akt activity dynamics may contribute less to the decision driving cell division in this context. We also find that ERK and Akt activities are less correlated with each other in cells that divide. Network reconstruction experiments demonstrated that this correlation behavior was likely not due to crosstalk, as ERK and Akt do not interact in this context, in contrast to other transformed cell types. Overall, our findings support roles for ERK and Akt activity throughout the cell cycle as opposed to just before the restriction point, and suggest ERK activity dynamics may be more important than Akt activity dynamics for driving cell division in this non-transformed context.

P. aeruginosa type III and type VI secretion systems modulate early response gene expression in type II pneumocytes in vitro

BACKGROUND

Lung airway epithelial cells are part of innate immunity and the frontline of defense against bacterial infections. During infection, airway epithelial cells secrete proinflammatory mediators that participate in the recruitment of immune cells. Virulence factors expressed by bacterial pathogens can alter epithelial cell gene expression and modulate this response. Pseudomonas aeruginosa, a Gram-negative opportunistic pathogen, expresses numerous virulence factors to facilitate establishment of infection and evade the host immune response. This study focused on identifying the role of two major P. aeruginosa virulence factors, type III (T3SS) and type VI (T6SS) secretion systems, on the early transcriptome response of airway epithelial cells in vitro.

RESULTS

We performed RNA-seq analysis of the transcriptome response of type II pneumocytes during infection with P. aeruginosa in vitro. We observed that P. aeruginosa differentially upregulates immediate-early response genes and transcription factors that induce proinflammatory responses in type II pneumocytes. P. aeruginosa infection of type II pneumocytes was characterized by up-regulation of proinflammatory networks, including MAPK, TNF, and IL-17 signaling pathways. We also identified early response genes and proinflammatory signaling pathways whose expression change in response to infection with P. aeruginosa T3SS and T6SS mutants in type II pneumocytes. We determined that T3SS and T6SS modulate the expression of EGR1, FOS, and numerous genes that are involved in proinflammatory responses in epithelial cells during infection. T3SS and T6SS were associated with two distinct transcriptomic signatures related to the activation of transcription factors such as AP1, STAT1, and SP1, and the secretion of pro-inflammatory cytokines such as IL-6 and IL-8.

CONCLUSIONS

Taken together, transcriptomic analysis of epithelial cells indicates that the expression of immediate-early response genes quickly changes upon infection with P. aeruginosa and this response varies depending on bacterial viability and injectosomes. These data shed light on how P. aeruginosa modulates host epithelial transcriptome response during infection using T3SS and T6SS.

Signalling dynamics, cell decisions, and homeostatic control in health and disease

Cell signalling engenders cells with the capability to receive and process information from the intracellular and extracellular environments, trigger and execute biological responses, and communicate with each other. Ultimately, cell signalling is responsible for maintaining homeostasis at the cellular, tissue and systemic level. For this reason, cell signalling is a topic of intense research efforts aimed to elucidate how cells coordinate transitions between states in developing and adult organisms in physiological and pathological conditions. Here, we review current knowledge of how cell signalling operates at multiple spatial and temporal scales, focusing on how single-cell analytical techniques reveal mechanisms underpinning cell-to-cell variability, signalling plasticity, and collective cellular responses.

Restoration of DNA integrity and cell cycle by electric stimulation in planarian tissues damaged by ionizing radiation

Exposure to high levels of ionizing γ-radiation leads to irreversible DNA damage and cell death. Here, we establish that exogenous application of electric stimulation enables cellular plasticity to reestablish stem cell activity in tissues damaged by ionizing radiation. We show that sub-threshold direct current stimulation (DCS) rapidly restores pluripotent stem cell populations previously eliminated by lethally γ-irradiated tissues of the planarian flatworm Schmidtea mediterranea. Our findings reveal that DCS enhances DNA repair, transcriptional activity, and cell cycle entry in post-mitotic cells. These responses involve rapid increases in cytosolic [Ca2+] through the activation of L-type Cav channels and intracellular Ca2+ stores leading to the activation of immediate early genes and ectopic expression of stem cell markers in postmitotic cells. Overall, we show the potential of electric current stimulation to reverse the damaging effects of high dose γ-radiation in adult tissues. Furthermore, our results provide mechanistic insights describing how electric stimulation effectively translates into molecular responses capable of regulating fundamental cellular functions without the need for genetic or pharmacological intervention.

Analysis of the early response to spinal cord injury identified a key role for mTORC1 signaling in the activation of neural stem progenitor cells

Xenopus laevis are able to regenerate the spinal cord during larvae stages through the activation of neural stem progenitor cells (NSPCs). Here we use high-resolution expression profiling to characterize the early transcriptome changes induced after spinal cord injury, aiming to identify the signals that trigger NSPC proliferation. The analysis delineates a pathway that starts with a rapid and transitory activation of immediate early genes, followed by migration processes and immune response genes, the pervasive increase of NSPC-specific ribosome biogenesis factors, and genes involved in stem cell proliferation. Western blot and immunofluorescence analysis showed that mTORC1 is rapidly and transiently activated after SCI, and its pharmacological inhibition impairs spinal cord regeneration and proliferation of NSPC through the downregulation of genes involved in the G1/S transition of cell cycle, with a strong effect on PCNA. We propose that the mTOR signaling pathway is a key player in the activation of NPSCs during the early steps of spinal cord regeneration.

Short-term transcriptomic response to plasma membrane injury

Plasma membrane repair mechanisms are activated within seconds post-injury to promote rapid membrane resealing in eukaryotic cells and prevent cell death. However, less is known about the regeneration phase that follows and how cells respond to injury in the short-term. Here, we provide a genome-wide study into the mRNA expression profile of MCF-7 breast cancer cells exposed to injury by digitonin, a mild non-ionic detergent that permeabilizes the plasma membrane. We focused on the early transcriptional signature and found a time-dependent increase in the number of differentially expressed (> twofold, P < 0.05) genes (34, 114 and 236 genes at 20-, 40- and 60-min post-injury, respectively). Pathway analysis highlighted a robust and gradual three-part transcriptional response: (1) prompt activation of immediate-early response genes, (2) activation of specific MAPK cascades and (3) induction of inflammatory and immune pathways. Therefore, plasma membrane injury triggers a rapid and strong stress and immunogenic response. Our meta-analysis suggests that this is a conserved transcriptome response to plasma membrane injury across different cell and injury types. Taken together, our study shows that injury has profound effects on the transcriptome of wounded cells in the regeneration phase (subsequent to membrane resealing), which is likely to influence cellular status and has been previously overlooked.

SFPQ Depletion Is Synthetically Lethal with BRAFV600E in Colorectal Cancer Cells

Oncoproteins such as the BRAF^V600E kinase endow cancer cells with malignant properties, but they also create unique vulnerabilities. Targeting of BRAF^V600E-driven cytoplasmic signaling networks has proved ineffective, as patients regularly relapse with reactivation of the targeted pathways. We identify the nuclear protein SFPQ to be synthetically lethal with BRAF^V600E in a loss-of-function shRNA screen. SFPQ depletion decreases proliferation and specifically induces S-phase arrest and apoptosis in BRAF^V600E-driven colorectal and melanoma cells. Mechanistically, SFPQ loss in BRAF-mutant cancer cells triggers the Chk1-dependent replication checkpoint, results in decreased numbers and reduced activities of replication factories, and increases collision between replication and transcription. We find that BRAF^V600E-mutant cancer cells and organoids are sensitive to combinations of Chk1 inhibitors and chemically induced replication stress, pointing toward future therapeutic approaches exploiting nuclear vulnerabilities induced by BRAF^V600E.

A rationally engineered decoder of transient intracellular signals

Cells can encode information about their environment by modulating signaling dynamics and responding accordingly. Yet, the mechanisms cells use to decode these dynamics remain unknown when cells respond exclusively to transient signals. Here, we approach design principles underlying such decoding by rationally engineering a synthetic short-pulse decoder in budding yeast. A computational method for rapid prototyping, TopoDesign, allowed us to explore 4122 possible circuit architectures, design targeted experiments, and then rationally select a single circuit for implementation. This circuit demonstrates short-pulse decoding through incoherent feedforward and positive feedback. We predict incoherent feedforward to be essential for decoding transient signals, thereby complementing proposed design principles of temporal filtering, the ability to respond to sustained signals, but not to transient signals. More generally, we anticipate TopoDesign to help designing other synthetic circuits with non-intuitive dynamics, simply by assembling available biological components.

Deciphering key regulators involved in epilepsy-induced cardiac damage through whole transcriptome and proteome analysis in a rat model

OBJECTIVE

Sudden unexpected death in epilepsy (SUDEP) is a major outcome of cardiac dysfunction in patients with epilepsy. In continuation of our previous work, the present study was envisaged to explore the key regulators responsible for cardiac damage associated with chronic seizures using whole transcriptome and proteome analysis in a rat model of temporal lobe epilepsy.

METHODS

A standard lithium-pilocarpine protocol was used to induce recurrent seizures in rats. The isolated rat heart tissue was subjected to transcriptomic and proteomic analysis. An integrated approach of RNA-Seq, proteomics, and system biology analysis was used to identify key regulators involved in seizure-linked cardiac changes. The analyzed differential expression patterns and network interactions were supported by gene and protein expression studies.

RESULTS

Altogether, 1157 differentially expressed genes and 1264 proteins were identified in the cardiac tissue of epileptic animals through RNA-Seq and liquid chromatography with tandem mass spectrometry-based proteomic analysis, respectively. The network analysis revealed seven critical genes-STAT3, Myc, Fos, Erbb2, Erbb3, Notch1, and Mapk8-that could play a role in seizure-mediated cardiac changes. The LC-MS/MS analysis supported the activation of the transforming growth factor β (TGF-β) pathway in the heart of epileptic animals. Furthermore, our gene and protein expression studies established a key role of STAT3, Erbb, and Mapk8 to develop cardiac changes linked with recurrent seizures.

SIGNIFICANCE

The present multi-omics study identified STAT3, Mapk8, and Erbb as key regulators involved in seizure-associated cardiac changes. It provided a deeper understanding of molecular, cellular, and network-level operations of the identified regulators that lead to cardiac changes in epilepsy.

Kinetic modeling of stem cell transcriptome dynamics to identify regulatory modules of normal and disturbed neuroectodermal differentiation

Thousands of transcriptome data sets are available, but approaches for their use in dynamic cell response modelling are few, especially for processes affected simultaneously by two orthogonal influencing variables. We approached this problem for neuroepithelial development of human pluripotent stem cells (differentiation variable), in the presence or absence of valproic acid (signaling variable). Using few basic assumptions (sequential differentiation states of cells; discrete on/off states for individual genes in these states), and time-resolved transcriptome data, a comprehensive model of spontaneous and perturbed gene expression dynamics was developed. The model made reliable predictions (average correlation of 0.85 between predicted and subsequently tested expression values). Even regulations predicted to be non-monotonic were successfully validated by PCR in new sets of experiments. Transient patterns of gene regulation were identified from model predictions. They pointed towards activation of Wnt signaling as a candidate pathway leading to a redirection of differentiation away from neuroepithelial cells towards neural crest. Intervention experiments, using a Wnt/beta-catenin antagonist, led to a phenotypic rescue of this disturbed differentiation. Thus, our broadly applicable model allows the analysis of transcriptome changes in complex time/perturbation matrices.

PWD/Ph-Encoded Genetic Variants Modulate the Cellular Wnt/β-Catenin Response to Suppress Apc Min-Triggered Intestinal Tumor Formation

Genetic predisposition affects the penetrance of tumor-initiating mutations, such as APC mutations that stabilize β-catenin and cause intestinal tumors in mice and humans. However, the mechanisms involved in genetically predisposed penetrance are not well understood. Here, we analyzed tumor multiplicity and gene expression in tumor-prone Apc ^Min/+ mice on highly variant C57BL/6J (B6) and PWD/Ph (PWD) genetic backgrounds. (B6 × PWD) F1 APC ^Min offspring mice were largely free of intestinal adenoma, and several chromosome substitution (consomic) strains carrying single PWD chromosomes on the B6 genetic background displayed reduced adenoma numbers. Multiple dosage-dependent modifier loci on PWD chromosome 5 each contributed to tumor suppression. Activation of β-catenin-driven and stem cell-specific gene expression in the presence of Apc ^Min or following APC loss remained moderate in intestines carrying PWD chromosome 5, suggesting that PWD variants restrict adenoma initiation by controlling stem cell homeostasis. Gene expression of modifier candidates and DNA methylation on chromosome 5 were predominantly cis controlled and largely reflected parental patterns, providing a genetic basis for inheritance of tumor susceptibility. Human SNP variants of several modifier candidates were depleted in colorectal cancer genomes, suggesting that similar mechanisms may also affect the penetrance of cancer driver mutations in humans. Overall, our analysis highlights the strong impact that multiple genetic variants acting in networks can exert on tumor development. SIGNIFICANCE: These findings in mice show that, in addition to accidental mutations, cancer risk is determined by networks of individual gene variants.

SPEED2: inferring upstream pathway activity from differential gene expression

Extracting signalling pathway activities from transcriptome data is important to infer mechanistic origins of transcriptomic dysregulation, for example in disease. A popular method to do so is by enrichment analysis of signature genes in e.g. differentially regulated genes. Previously, we derived signatures for signalling pathways by integrating public perturbation transcriptome data and generated a signature database called SPEED (Signalling Pathway Enrichment using Experimental Datasets), for which we here present a substantial upgrade as SPEED2. This web server hosts consensus signatures for 16 signalling pathways that are derived from a large number of transcriptomic signalling perturbation experiments. When providing a gene list of e.g. differentially expressed genes, the web server allows to infer signalling pathways that likely caused these genes to be deregulated. In addition to signature lists, we derive ‘continuous’ gene signatures, in a transparent and automated fashion without any fine-tuning, and describe a new algorithm to score these signatures.

Genome-wide dynamics of RNA synthesis, processing, and degradation without RNA metabolic labeling

The quantification of the kinetic rates of RNA synthesis, processing, and degradation are largely based on the integrative analysis of total and nascent transcription, the latter being quantified through RNA metabolic labeling. We developed INSPEcT−, a computational method based on the mathematical modeling of premature and mature RNA expression that is able to quantify kinetic rates from steady-state or time course total RNA-seq data without requiring any information on nascent transcripts. Our approach outperforms available solutions, closely recapitulates the kinetic rates obtained through RNA metabolic labeling, improves the ability to detect changes in transcript half-lives, reduces the cost and complexity of the experiments, and can be adopted to study experimental conditions in which nascent transcription cannot be readily profiled. Finally, we applied INSPEcT− to the characterization of post-transcriptional regulation landscapes in dozens of physiological and disease conditions. This approach was included in the INSPEcT Bioconductor package, which can now unveil RNA dynamics from steady-state or time course data, with or without the profiling of nascent RNA.

Systems-Level Properties of EGFR-RAS-ERK Signaling Amplify Local Signals to Generate Dynamic Gene Expression Heterogeneity

Intratumoral heterogeneity is associated with aggressive tumor behavior, therapy resistance, and poor patient outcomes. Such heterogeneity is thought to be dynamic, shifting over periods of minutes to hours in response to signaling inputs from the tumor microenvironment. However, models of this process have been inferred from indirect or post-hoc measurements of cell state, leaving the temporal details of signaling-driven heterogeneity undefined. Here, we developed a live-cell model system in which microenvironment-driven signaling dynamics can be directly observed and linked to variation in gene expression. Our analysis reveals that paracrine signaling between two cell types is sufficient to drive continual diversification of gene expression programs. This diversification emerges from systems-level properties of the EGFR-RAS-ERK signaling cascade, including intracellular amplification of amphiregulin-mediated paracrine signals and differential kinetic filtering by target genes including Fra-1, c-Myc, and Egr1. Our data enable more precise modeling of paracrine-driven transcriptional variation as a generator of gene expression heterogeneity. A record of this paper's transparent peer review process is included in the Supplemental Information.

Temporal perturbation of ERK dynamics reveals network architecture of FGF2/MAPK signaling

Stimulation of PC-12 cells with epidermal (EGF) versus nerve (NGF) growth factors (GFs) biases the distribution between transient and sustained single-cell ERK activity states, and between proliferation and differentiation fates within a cell population. We report that fibroblast GF (FGF2) evokes a distinct behavior that consists of a gradually changing population distribution of transient/sustained ERK signaling states in response to increasing inputs in a dose response. Temporally controlled GF perturbations of MAPK signaling dynamics applied using microfluidics reveal that this wider mix of ERK states emerges through the combination of an intracellular feedback, and competition of FGF2 binding to FGF receptors (FGFRs) and heparan sulfate proteoglycan (HSPG) co-receptors. We show that the latter experimental modality is instructive for model selection using a Bayesian parameter inference. Our results provide novel insights into how different receptor tyrosine kinase (RTK) systems differentially wire the MAPK network to fine-tune fate decisions at the cell population level.

Gene Regulatory Strategies that Decode the Duration of NFκB Dynamics Contribute to LPS- versus TNF-Specific Gene Expression

Pathogen-derived lipopolysaccharide (LPS) and cytokine tumor necrosis factor (TNF) activate NFκB with distinct duration dynamics, but how immune response genes decode NFκB duration to produce stimulus-specific expression remains unclear. Here, detailed transcriptomic profiling of combinatorial and temporal control mutants identified 81 genes that depend on stimulus-specific NFκB duration for their stimulus-specificity. Combining quantitative experimentation with mathematical modeling, we found that for some genes a long mRNA half-life allowed effective decoding, but for many genes this was insufficient to account for the data; instead, we found that chromatin mechanisms, such as a slow transition rate between inactive and RelA-bound enhancer states, could also decode NFκB dynamics. Chromatin-mediated decoding is favored by genes acting as immune effectors (e.g., tissue remodelers and T cell recruiters) rather than immune regulators (e.g., signaling proteins and monocyte recruiters). Overall, our results delineate two gene regulatory strategies that decode stimulus-specific NFκB dynamics and determine distinct biological functions.

Temporal perturbation of ERK dynamics reveals network architecture of FGF2/MAPK signaling

Stimulation of PC-12 cells with epidermal (EGF) versus nerve (NGF) growth factors (GFs) biases the distribution between transient and sustained single-cell ERK activity states, and between proliferation and differentiation fates within a cell population. We report that fibroblast GF (FGF2) evokes a distinct behavior that consists of a gradually changing population distribution of transient/sustained ERK signaling states in response to increasing inputs in a dose response. Temporally controlled GF perturbations of MAPK signaling dynamics applied using microfluidics reveal that this wider mix of ERK states emerges through the combination of an intracellular feedback, and competition of FGF2 binding to FGF receptors (FGFRs) and heparan sulfate proteoglycan (HSPG) co-receptors. We show that the latter experimental modality is instructive for model selection using a Bayesian parameter inference. Our results provide novel insights into how different receptor tyrosine kinase (RTK) systems differentially wire the MAPK network to fine-tune fate decisions at the cell population level.

Cell type-dependent differential activation of ERK by oncogenic KRAS in colon cancer and intestinal epithelium

Oncogenic mutations in KRAS or BRAF are frequent in colorectal cancer and activate the ERK kinase. Here, we find graded ERK phosphorylation correlating with cell differentiation in patient-derived colorectal cancer organoids with and without KRAS mutations. Using reporters, single cell transcriptomics and mass cytometry, we observe cell type-specific phosphorylation of ERK in response to transgenic KRAS^G12V in mouse intestinal organoids, while transgenic BRAF^V600E activates ERK in all cells. Quantitative network modelling from perturbation data reveals that activation of ERK is shaped by cell type-specific MEK to ERK feed forward and negative feedback signalling. We identify dual-specificity phosphatases as candidate modulators of ERK in the intestine. Furthermore, we find that oncogenic KRAS, together with β-Catenin, favours expansion of crypt cells with high ERK activity. Our experiments highlight key differences between oncogenic BRAF and KRAS in colorectal cancer and find unexpected heterogeneity in a signalling pathway with fundamental relevance for cancer therapy.

KRAS^G12V and BRAF^V600E are oncogenic mutations that activate ERK signalling. Here, the authors use single cell analysis in intestinal organoids and show that BRAF^V600E activates ERK in all intestinal cell types, while KRAS^G12V induces ERK activation in only a subset of cells, depending on cell differentiation state.

Communication codes in developmental signaling pathways

A handful of core intercellular signaling pathways play pivotal roles in a broad variety of developmental processes. It has remained puzzling how so few pathways can provide the precision and specificity of cell-cell communication required for multicellular development. Solving this requires us to quantitatively understand how developmentally relevant signaling information is actively sensed, transformed and spatially distributed by signaling pathways. Recently, single cell analysis and cell-based reconstitution, among other approaches, have begun to reveal the 'communication codes' through which information is represented in the identities, concentrations, combinations and dynamics of extracellular ligands. They have also revealed how signaling pathways decipher these features and control the spatial distribution of signaling in multicellular contexts. Here, we review recent work reporting the discovery and analysis of communication codes and discuss their implications for diverse developmental processes.

Cancer mutations and targeted drugs can disrupt dynamic signal encoding by the Ras-Erk pathway

The Ras-Erk (extracellular signal-regulated kinase) pathway encodes information in its dynamics; the duration and frequency of Erk activity can specify distinct cell fates. To enable dynamic encoding, temporal information must be accurately transmitted from the plasma membrane to the nucleus. We used optogenetic profiling to show that both oncogenic B-Raf mutations and B-Raf inhibitors can cause corruption of this transmission, so that short pulses of input Ras activity are distorted into abnormally long Erk outputs. These changes can reshape downstream transcription and cell fates, resulting in improper decisions to proliferate. These findings illustrate how altered dynamic signal transmission properties, and not just constitutively increased signaling, can contribute to cell proliferation and perhaps cancer, and how optogenetic profiling can dissect mechanisms of signaling dysfunction in disease.

Microenvironmental Signals and Biochemical Information Processing: Cooperative Determinants of Intratumoral Plasticity and Heterogeneity

Intra-tumor cellular heterogeneity is a major challenge in cancer therapy. Tumors are composed of multiple phenotypic subpopulations that vary in their ability to initiate metastatic tumors and in their sensitivity to chemotherapy. In many cases, cells can transition between these subpopulations, not by genetic mutation, but instead through reversible changes in signal transduction or gene expression programs. This plasticity begins at the level of the microenvironment where local autocrine and paracrine signals, exosomes, tumor–stroma interactions, and extracellular matrix (ECM) composition create a signaling landscape that varies over space and time. The integration of this complex array of signals engages signaling pathways that control gene expression. The resulting modulation of gene expression programs causes individual cells to sample a wide array of phenotypic states that support tumor growth, dissemination, and therapeutic resistance. In this review, we discuss how information flows dynamically within the microenvironmental landscape to inform cell state decisions and to create intra-tumoral heterogeneity. We address the role of plasticity in the acquisition of transient and prolonged drug resistant states and discuss how targeted pharmacological modification of the signaling landscape may be able to constrain phenotypic plasticity, leading to improved treatment responses.

Linear Integration of ERK Activity Predominates over Persistence Detection in Fra-1 Regulation

ERK signaling regulates the expression of target genes, but it is unclear how ERK activity dynamics are interpreted. Here, we investigate this question using simultaneous, live, single-cell imaging of two ERK activity reporters and expression of Fra-1, a target gene controlling epithelial cell identity. We find that Fra-1 is expressed in proportion to the amplitude and duration of ERK activity. In contrast to previous "persistence detector" and "selective filter" models in which Fra-1 expression only occurs when ERK activity persists beyond a threshold duration, our observations demonstrate that the network regulating Fra-1 expression integrates total ERK activity and responds to it linearly. However, exploration of a generalized mathematical model of the Fra-1 coherent feedforward loop demonstrates that it can perform either linear integration or persistence detection, depending on the basal mRNA production rate and protein production delays. Our data indicate that significant basal expression and short delays cause Fra-1 to respond linearly to integrated ERK activity.

Widespread activation of antisense transcription of the host genome during herpes simplex virus 1 infection

Background

Herpesviruses can infect a wide range of animal species. Herpes simplex virus 1 (HSV-1) is one of the eight herpesviruses that can infect humans and is prevalent worldwide. Herpesviruses have evolved multiple ways to adapt the infected cells to their needs, but knowledge about these transcriptional and post-transcriptional modifications is sparse.

Results

Here, we show that HSV-1 induces the expression of about 1000 antisense transcripts from the human host cell genome. A subset of these is also activated by the closely related varicella zoster virus. Antisense transcripts originate either at gene promoters or within the gene body, and they show different susceptibility to the inhibition of early and immediate early viral gene expression. Overexpression of the major viral transcription factor ICP4 is sufficient to turn on a subset of antisense transcripts. Histone marks around transcription start sites of HSV-1-induced and constitutively transcribed antisense transcripts are highly similar, indicating that the genetic loci are already poised to transcribe these novel RNAs. Furthermore, an antisense transcript overlapping with the BBC3 gene (also known as PUMA) transcriptionally silences this potent inducer of apoptosis in cis.

Conclusions

We show for the first time that a virus induces widespread antisense transcription of the host cell genome. We provide evidence that HSV-1 uses this to downregulate a strong inducer of apoptosis. Our findings open new perspectives on global and specific alterations of host cell transcription by viruses.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1329-5) contains supplementary material, which is available to authorized users.

An immediate-late gene expression module decodes ERK signal duration

The RAF‐MEK‐ERK signalling pathway controls fundamental, often opposing cellular processes such as proliferation and apoptosis. Signal duration has been identified to play a decisive role in these cell fate decisions. However, it remains unclear how the different early and late responding gene expression modules can discriminate short and long signals. We obtained both protein phosphorylation and gene expression time course data from HEK293 cells carrying an inducible construct of the proto‐oncogene RAF. By mathematical modelling, we identified a new gene expression module of immediate–late genes (ILGs) distinct in gene expression dynamics and function. We find that mRNA longevity enables these ILGs to respond late and thus translate ERK signal duration into response amplitude. Despite their late response, their GC‐rich promoter structure suggested and metabolic labelling with 4SU confirmed that transcription of ILGs is induced immediately. A comparative analysis shows that the principle of duration decoding is conserved in PC12 cells and MCF7 cells, two paradigm cell systems for ERK signal duration. Altogether, our findings suggest that ILGs function as a gene expression module to decode ERK signal duration.