Quantifying the Central Dogma in the p53 Pathway in Live Single Cells

Time will tell: comparing timescales to gain insight into transcriptional bursting

Recent imaging studies have captured the dynamics of regulatory events of transcription inside living cells. These events include transcription factor (TF) DNA binding, chromatin remodeling and modification, enhancer-promoter (E-P) proximity, cluster formation, and preinitiation complex (PIC) assembly. Together, these molecular events culminate in stochastic bursts of RNA synthesis, but their kinetic relationship remains largely unclear. In this review, we compare the timescales of upstream regulatory steps (input) with the kinetics of transcriptional bursting (output) to generate mechanistic models of transcription dynamics in single cells. We highlight open questions and potential technical advances to guide future endeavors toward a quantitative and kinetic understanding of transcription regulation.

Identification of genes with oscillatory expression in glioblastoma: the paradigm of SOX2

Quiescence, a reversible state of cell-cycle arrest, is an important state during both normal development and cancer progression. For example, in glioblastoma (GBM) quiescent glioblastoma stem cells (GSCs) play an important role in re-establishing the tumour, leading to relapse. While most studies have focused on identifying differentially expressed genes between proliferative and quiescent cells as potential drivers of this transition, recent studies have shown the importance of protein oscillations in controlling the exit from quiescence of neural stem cells. Here, we have undertaken a genome-wide bioinformatic inference approach to identify genes whose expression oscillates and which may be good candidates for controlling the transition to and from the quiescent cell state in GBM. Our analysis identified, among others, a list of important transcription regulators as potential oscillators, including the stemness gene SOX2, which we verified to oscillate in quiescent GSCs. These findings expand on the way we think about gene regulation and introduce new candidate genes as key regulators of quiescence.

Single-cell transcriptional dynamics in a living vertebrate

The ability to quantify transcriptional dynamics in individual cells via live imaging has revolutionized our understanding of gene regulation. However, such measurements are lacking in the context of vertebrate embryos. We addressed this deficit by applying MS2-MCP mRNA labeling to the quantification of transcription in zebrafish, a model vertebrate. We developed a platform of transgenic organisms, light sheet fluorescence microscopy, and optimized image analysis that enables visualization and quantification of MS2 reporters. We used these tools to obtain the first single-cell, real-time measurements of transcriptional dynamics of the segmentation clock. Our measurements challenge the traditional view of smooth clock oscillations and instead suggest a model of discrete transcriptional bursts that are organized in space and time. Together, these results highlight how measuring single-cell transcriptional activity can reveal unexpected features of gene regulation and how this data can fuel the dialogue between theory and experiment.

A dual role of RBM42 in modulating splicing and translation of CDKN1A/p21 during DNA damage response

p53-mediated cell cycle arrest during DNA damage is dependent on the induction of p21 protein, encoded by the CDKN1A gene. p21 inhibits cyclin-dependent kinases required for cell cycle progression to guarantee accurate repair of DNA lesions. Hence, fine-tuning of p21 levels is crucial to preserve genomic stability. Currently, the multilayered regulation of p21 levels during DNA damage is not fully understood. Herein, we identify the human RNA binding motif protein 42 (RBM42) as a regulator of p21 levels during DNA damage. Genome-wide transcriptome and interactome analysis reveals that RBM42 alters the expression of p53-regulated genes during DNA damage. Specifically, we demonstrate that RBM42 facilitates CDKN1A splicing by counteracting the splicing inhibitory effect of RBM4 protein. Unexpectedly, we also show that RBM42, underpins translation of various splicing targets, including CDKN1A. Concordantly, transcriptome-wide mapping of RBM42-RNA interactions using eCLIP further substantiates the dual function of RBM42 in regulating splicing and translation of its target genes, including CDKN1A. Collectively, our data show that RBM42 couples splicing and translation machineries to fine-tune gene expression during DNA damage response.

Differential cell-cycle control by oscillatory versus sustained Hes1 expression via p21

Oscillatory Hes1 expression activates cell proliferation, while high and sustained Hes1 expression induces quiescence, but the mechanism by which Hes1 differentially controls cell proliferation depending on its expression dynamics is unclear. Here, we show that oscillatory Hes1 expression down-regulates the expression of the cyclin-dependent kinase inhibitor p21 (Cdkn1a), which delays cell-cycle progression, and thereby activates the proliferation of mouse neural stem cells (NSCs). By contrast, sustained Hes1 overexpression up-regulates p21 expression and inhibits NSC proliferation, although it initially down-regulates p21 expression. Compared with Hes1 oscillation, sustained Hes1 overexpression represses Dusp7, a phosphatase for phosphorylated Erk (p-Erk), and increases the levels of p-Erk, which can up-regulate p21 expression. These results indicate that p21 expression is directly repressed by oscillatory Hes1 expression, but indirectly up-regulated by sustained Hes1 overexpression, suggesting that depending on its expression dynamics, Hes1 differentially controls NSC proliferation via p21.

Long-term p21 and p53 dynamics regulate the frequency of mitosis events and cell cycle arrest following radiation damage

Radiation exposure of healthy cells can halt cell cycle temporarily or permanently. In this work, we analyze the time evolution of p21 and p53 from two single cell datasets of retinal pigment epithelial cells exposed to several levels of radiation, and in particular, the effect of radiation on cell cycle arrest. Employing various quantification methods from signal processing, we show how p21 levels, and to a lesser extent p53 levels, dictate whether the cells are arrested in their cell cycle and how frequently these mitosis events are likely to occur. We observed that single cells exposed to the same dose of DNA damage exhibit heterogeneity in cellular outcomes and that the frequency of cell division is a more accurate monitor of cell damage rather than just radiation level. Finally, we show how heterogeneity in DNA damage signaling is manifested early in the response to radiation exposure level and has potential to predict long-term fate.

Encoding and Decoding of p53 Dynamics in Cellular Response to Stresses

In the cellular response to stresses, the tumor suppressor p53 is activated to maintain genomic integrity and fidelity. As a transcription factor, p53 exhibits rich dynamics to allow for discrimination of the type and intensity of stresses and to direct the selective activation of target genes involved in different processes including cell cycle arrest and apoptosis. In this review, we focused on how stresses are encoded into p53 dynamics and how the dynamics are decoded into cellular outcomes. Theoretical modeling may provide a global view of signaling in the p53 network by coupling the encoding and decoding processes. We discussed the significance of modeling in revealing the mechanisms of the transition between p53 dynamic modes. Moreover, we shed light on the crosstalk between the p53 network and other signaling networks. This review may advance the understanding of operating principles of the p53 signaling network comprehensively and provide insights into p53 dynamics-based cancer therapy.

A common pathway to cancer: Oncogenic mutations abolish p53 oscillations

The tumor suppressor p53 oscillates in response to DNA double-strand breaks, a behavior that has been suggested to be essential to its anti-cancer function. Nearly all human cancers have genetic alterations in the p53 pathway; a number of these alterations have been shown to be oncogenic by experiment. These alterations include somatic mutations and copy number variations as well as germline polymorphisms. Intriguingly, they exhibit a mixed pattern of interactions in tumors, such as co-occurrence, mutual exclusivity, and paradoxically, mutual antagonism. Using a differential equation model of p53-Mdm2 dynamics, we employ Hopf bifurcation analysis to show that these alterations have a common mode of action, to abolish the oscillatory competence of p53, thereby, we suggest, impairing its tumor suppressive function. In this analysis, diverse genetic alterations, widely associated with human cancers clinically, have a unified mechanistic explanation of their role in oncogenesis.

Inferring gene regulation from stochastic transcriptional variation across single cells at steady state

Regulatory relationships between transcription factors (TFs) and their target genes lie at the heart of cellular identity and function; however, uncovering these relationships is often labor-intensive and requires perturbations. Here, we propose a principled framework to systematically infer gene regulation for all TFs simultaneously in cells at steady state by leveraging the intrinsic variation in the transcriptional abundance across single cells. Through modeling and simulations, we characterize how transcriptional bursts of a TF gene are propagated to its target genes, including the expected ranges of time delay and magnitude of maximum covariation. We distinguish these temporal trends from the time-invariant covariation arising from cell states, and we delineate the experimental and technical requirements for leveraging these small but meaningful cofluctuations in the presence of measurement noise. While current technology does not yet allow adequate power for definitively detecting regulatory relationships for all TFs simultaneously in cells at steady state, we investigate a small-scale dataset to inform future experimental design. This study supports the potential value of mapping regulatory connections through stochastic variation, and it motivates further technological development to achieve its full potential.

Unraveling the effect of intra- and intercellular processes on acetaminophen-induced liver injury

In high dosages, acetaminophen (APAP) can cause severe liver damage, but susceptibility to liver failure varies across individuals and is influenced by factors such as health status. Because APAP-induced liver injury and recovery is regulated by an intricate system of intra- and extracellular molecular signaling, we here aim to quantify the importance of specific modules in determining the outcome after an APAP insult and of potential targets for therapies that mitigate adversity. For this purpose, we integrated hepatocellular acetaminophen metabolism, DNA damage response induction and cell fate into a multiscale mechanistic liver lobule model which involves various cell types, such as hepatocytes, residential Kupffer cells and macrophages. Our model simulations show that zonal differences in metabolism and detoxification efficiency are essential determinants of necrotic damage. Moreover, the extent of senescence, which is regulated by intracellular processes and triggered by extracellular signaling, influences the potential to recover. In silico therapies at early and late time points after APAP insult indicated that prevention of necrotic damage is most beneficial for recovery, whereas interference with regulation of senescence promotes regeneration in a less pronounced way.

Model-based translation of DNA damage signaling dynamics across cell types

Interindividual variability in DNA damage response (DDR) dynamics may evoke differences in susceptibility to cancer. However, pathway dynamics are often studied in cell lines as alternative to primary cells, disregarding variability. To compare DDR dynamics in the cell line HepG2 with primary human hepatocytes (PHHs), we developed a HepG2-based computational model that describes the dynamics of DDR regulator p53 and targets MDM2, p21 and BTG2. We used this model to generate simulations of virtual PHHs and compared the results to those for PHH donor samples. Correlations between baseline p53 and p21 or BTG2 mRNA expression in the absence and presence of DNA damage for HepG2-derived virtual samples matched the moderately positive correlations observed for 50 PHH donor samples, but not the negative correlations between p53 and its inhibitor MDM2. Model parameter manipulation that affected p53 or MDM2 dynamics was not sufficient to accurately explain the negative correlation between these genes. Thus, extrapolation from HepG2 to PHH can be done for some DDR elements, yet our analysis also reveals a knowledge gap within p53 pathway regulation, which makes such extrapolation inaccurate for the regulator MDM2. This illustrates the relevance of studying pathway dynamics in addition to gene expression comparisons to allow reliable translation of cellular responses from cell lines to primary cells. Overall, with our approach we show that dynamical modeling can be used to improve our understanding of the sources of interindividual variability of pathway dynamics.

Susceptibility to develop cancer varies among people, partially due to differences in genetic background. Ideally, healthy human-derived cells are used to investigate intracellular signaling pathways and their interindividual variability contributing to cancer susceptibility. Because cells from healthy human tissue are difficult to obtain and culture for periods longer than a few days, cell lines are often used as substitute. However, it is unclear to what extent signaling dynamics in cell lines represent dynamics in healthy human tissue. We asked whether we could reproduce interindividual variability in DNA damage response gene expression in a set of 50 human liver cell donors. Therefore, we built a mathematical model that simulates temporal expression dynamics of the DNA damage response in the HepG2 liver cell line upon chemical activation and used the simulations to create virtual donors. Our virtual donors displayed similar relations between genes as the samples from human donors, provided that we adjusted the strength of specific molecular interactions. Thus, our approach can be used to examine the applicability of widely used cell systems to healthy human tissue in terms of their dynamic responses.

Louder for longer: Myc amplifies gene expression by extended transcriptional bursting

Transcription is a complex, dynamic process. Using live single-cell measurements, Patange et al. show, in a recent issue of Cell Reports, that elevated levels of the transcription factor MYC enhance target gene RNA production by increasing the duration but not frequency of transcriptional bursts.

The Origins of Phenotypic Heterogeneity in Cancer

Heterogeneity is a pervasive feature of cancer, and understanding the sources and regulatory mechanisms underlying heterogeneity could provide key insights to help improve the diagnosis and treatment of cancer. In this review, we discuss the origin of heterogeneity in the phenotype of individual cancer cells. Genotype-phenotype (G-P) maps are widely used in evolutionary biology to represent the complex interactions of genes and the environment that lead to phenotypes that impact fitness. Here, we present the rationale of an extended G-P (eG-P) map with a cone structure in cancer. The eG-P cone is formed by cells that are similar at the genome layer but gradually increase variability in the epigenome, transcriptome, proteome, metabolome and signalome layers to produce large variability at the phenome layer. Experimental evidence from single-cell -omics analyses supporting the cancer eG-P cone concept is presented, and the impact of epimutations and the interaction of cancer and tumor microenvironmental eG-P cones are integrated with the current understanding of cancer biology. The eG-P cone concept uncovers potential therapeutic strategies to reduce cancer evolution and improve cancer treatment. More methods to study phenotypes in single cells will be key to better understand cancer cell fitness in tumor biology and therapeutics.

Role of Episamarcandin in Promoting the Apoptosis of Human Colon Cancer HCT116 Cells through the PI3K-Akt Signaling Pathway

This study identifies the active ingredients of Ferula sinkiangensis and investigates the role and mechanism of episamarcandin in colon cancer cells. The silica gel column chromatography was utilized to separate the chemical components of Ferula sinkiangensis. Sephadex LH-20 and semipreparative HPLC were adopted for further separation and purification. The compound episamarcandin showed good anticolon cancer activity among the 13 monomeric compounds obtained. Its effects on the apoptosis, cell cycle, and invasion and migration of colon cancer HCT 116 cells and PI3K-Akt signaling pathway were further investigated. The results showed that, similar to positive control cisplatin, episamarcandin inhibited the proliferation, promoted the apoptosis, arrested cells at G0/G1 phase, and suppressed migration and invasion of HCT 116 cells. A large number of apoptotic HCT 116 cells were observed under a transmission electron microscope. Fluorescence real-time quantitative PCR and western blot analysis showed that episamarcandin increased the expression of PTEN, p53, and Bax and decreased the expression of P-Akt, Akt, mTOR, Bcl-xl, and Bcl-2. Conclusively, episamarcandin may inhibit cell proliferation, migration, and invasion and promote the apoptosis of human colon cancer HCT 116 cells possibly through the PI3K-Akt signaling pathway.

Quantitative imaging of RNA polymerase II activity in plants reveals the single-cell basis of tissue-wide transcriptional dynamics

The responses of plants to their environment are often dependent on the spatiotemporal dynamics of transcriptional regulation. While live-imaging tools have been used extensively to quantitatively capture rapid transcriptional dynamics in living animal cells, the lack of implementation of these technologies in plants has limited concomitant quantitative studies in this kingdom. Here, we applied the PP7 and MS2 RNA-labelling technologies for the quantitative imaging of RNA polymerase II activity dynamics in single cells of living plants as they respond to experimental treatments. Using this technology, we counted nascent RNA transcripts in real time in Nicotiana benthamiana (tobacco) and Arabidopsis thaliana. Examination of heat shock reporters revealed that plant tissues respond to external signals by modulating the proportion of cells that switch from an undetectable basal state to a high-transcription state, instead of modulating the rate of transcription across all cells in a graded fashion. This switch-like behaviour, combined with cell-to-cell variability in transcription rate, results in mRNA production variability spanning three orders of magnitude. We determined that cellular heterogeneity stems mainly from stochasticity intrinsic to individual alleles instead of variability in cellular composition. Together, our results demonstrate that it is now possible to quantitatively study the dynamics of transcriptional programs in single cells of living plants.

The Genomic Response to TGF-β1 Dictates Failed Repair and Progression of Fibrotic Disease in the Obstructed Kidney

Tubulointerstitial fibrosis is a common and diagnostic hallmark of a spectrum of chronic renal disorders. While the etiology varies as to the causative nature of the underlying pathology, persistent TGF-β1 signaling drives the relentless progression of renal fibrotic disease. TGF-β1 orchestrates the multifaceted program of kidney fibrogenesis involving proximal tubular dysfunction, failed epithelial recovery or re-differentiation, capillary collapse and subsequent interstitial fibrosis eventually leading to chronic and ultimately end-stage disease. An increasing complement of non-canonical elements function as co-factors in TGF-β1 signaling. p53 is a particularly prominent transcriptional co-regulator of several TGF-β1 fibrotic-response genes by complexing with TGF-β1 receptor-activated SMADs. This cooperative p53/TGF-β1 genomic cluster includes genes involved in cellular proliferative control, survival, apoptosis, senescence, and ECM remodeling. While the molecular basis for this co-dependency remains to be determined, a subset of TGF-β1-regulated genes possess both p53- and SMAD-binding motifs. Increases in p53 expression and phosphorylation, moreover, are evident in various forms of renal injury as well as kidney allograft rejection. Targeted reduction of p53 levels by pharmacologic and genetic approaches attenuates expression of the involved genes and mitigates the fibrotic response confirming a key role for p53 in renal disorders. This review focuses on mechanisms underlying TGF-β1-induced renal fibrosis largely in the context of ureteral obstruction, which mimics the pathophysiology of pediatric unilateral ureteropelvic junction obstruction, and the role of p53 as a transcriptional regulator within the TGF-β1 repertoire of fibrosis-promoting genes.

Gene Expression at a Single Molecule Level: Implications for MDS and AML

Non-genetic heterogeneity, or gene expression stochasticity, is an important source of variability in biological systems. With the advent and improvement of single molecule resolution technologies, it has been shown that transcription dynamics and resultant transcript number fluctuations generate significant cell-to-cell variability that has important biological effects and may contribute substantially to both tissue homeostasis and disease. In this respect, the pathophysiology of stem cell-derived malignancies such as AML and MDS, which has historically been studied at the ensemble level, may require re-evaluation. To that end, it is our aim in this review to highlight the results of recent single-molecule, biophysical, and systems studies of gene expression dynamics, with the explicit purpose of demonstrating how the insights from these basic science studies may help inform and progress the field of leukemia biology and, ultimately, research into novel therapies.

The guardian's choice: how p53 enables context-specific decision-making in individual cells

p53 plays a central role in defending the genomic integrity of our cells. In response to genotoxic stress, this tumour suppressor orchestrates the expression of hundreds of target genes, which induce a variety of cellular outcomes ranging from damage repair to induction of apoptosis. In this review, we examine how the p53 response is regulated on several levels in individual cells to allow precise and context-specific fate decisions. We discuss that the p53 response is not only controlled by its canonical regulators but also controlled by interconnected signalling pathways that influence the dynamics of p53 accumulation upon damage and modulate its transcriptional activity at target gene promoters. Additionally, we consider how the p53 response is diversified through a variety of mechanisms at the promoter level and beyond to induce context-specific outcomes in individual cells. These layers of regulation allow p53 to react in a stimulus-specific manner and fine-tune its signalling according to the individual needs of a given cell, enabling it to take the right decision on survival or death.

Tumor suppressor p53: from engaging DNA to target gene regulation

The p53 transcription factor confers its potent tumor suppressor functions primarily through the regulation of a large network of target genes. The recent explosion of next generation sequencing protocols has enabled the study of the p53 gene regulatory network (GRN) and underlying mechanisms at an unprecedented depth and scale, helping us to understand precisely how p53 controls gene regulation. Here, we discuss our current understanding of where and how p53 binds to DNA and chromatin, its pioneer-like role, and how this affects gene regulation. We provide an overview of the p53 GRN and the direct and indirect mechanisms through which p53 affects gene regulation. In particular, we focus on delineating the ubiquitous and cell type-specific network of regulatory elements that p53 engages; reviewing our understanding of how, where, and when p53 binds to DNA and the mechanisms through which these events regulate transcription. Finally, we discuss the evolution of the p53 GRN and how recent work has revealed remarkable differences between vertebrates, which are of particular importance to cancer researchers using mouse models.

First Responders Shape a Prompt and Sharp NF-κB-Mediated Transcriptional Response to TNF-α

Nuclear factor (NF)-κB controls the transcriptional response to inflammatory signals by translocating into the nucleus, but we lack a single-cell characterization of the resulting transcription dynamics. Here we show that upon tumor necrosis factor (TNF)-α transcription of NF-κB target genes is heterogeneous in individual cells but results in an average nascent transcription profile that is prompt (i.e., occurs almost immediately) and sharp (i.e., increases and decreases rapidly) compared with NF-κB nuclear localization. Using an NF-κB-controlled MS2 reporter we show that the single-cell nascent transcription is more heterogeneous than NF-κB translocation dynamics, with a fraction of synchronized “first responders” that shape the average transcriptional profile and are more prone to respond to multiple TNF-α stimulations. A mathematical model combining NF-κB-mediated gene activation and a gene refractory state is able to reproduce these features. Our work shows how the expression of target genes induced by transcriptional activators can be heterogeneous across single cells and yet time resolved on average.

•

Nascent transcription upon TNF-α is heterogeneous, with a subset of “first responders”

•

The average nascent transcription is prompt and sharper than NF-κB response

•

First responders do not depend on NF-κB dynamics and respond more to pulsed stimuli

•

A model including NF-κB and a gene refractory state reproduces these observations