p53 pulses lead to distinct patterns of gene expression albeit similar DNA-binding dynamics

Temporal regulation of gene expression through integration of p53 dynamics and modifications

The master regulator of the DNA damage response, the transcription factor p53, orchestrates multiple downstream responses and coordinates repair processes. In response to double-strand DNA breaks, p53 exhibits pulses of expression, but how it achieves temporal coordination of downstream responses remains unclear. Here, we show that p53's posttranslational modification state is altered between its first and second pulses of expression. We show that acetylations at two sites, K373 and K382, were reduced in the second pulse, and these acetylations differentially affected p53 target genes, resulting in changes in gene expression programs over time. This interplay between dynamics and modification may offer a strategy for cellular hubs like p53 to temporally organize multiple processes in individual cells.

Entrainment and multi-stability of the p53 oscillator in human cells

The tumor suppressor p53 responds to cellular stress and activates transcription programs critical for regulating cell fate. DNA damage triggers oscillations in p53 levels with a robust period. Guided by the theory of synchronization and entrainment, we developed a mathematical model and experimental system to test the ability of the p53 oscillator to entrain to external drug pulses of various periods and strengths. We found that the p53 oscillator can be locked and entrained to a wide range of entrainment modes. External periods far from p53's natural oscillations increased the heterogeneity between individual cells whereas stronger inputs reduced it. Single-cell measurements allowed deriving the phase response curves (PRCs) and multiple Arnold tongues of p53. In addition, multi-stability and non-linear behaviors were mathematically predicted and experimentally detected, including mode hopping, period doubling, and chaos. Our work revealed critical dynamical properties of the p53 oscillator and provided insights into understanding and controlling it. A record of this paper's transparent peer review process is included in the supplemental information.

Determinants of p53 DNA binding, gene regulation, and cell fate decisions

The extent to which transcription factors read and respond to specific information content within short DNA sequences remains an important question that the tumor suppressor p53 is helping us answer. We discuss recent insights into how local information content at p53 binding sites might control modes of p53 target gene activation and cell fate decisions. Significant prior work has yielded data supporting two potential models of how p53 determines cell fate through its target genes: a selective target gene binding and activation model and a p53 level threshold model. Both of these models largely revolve around an analogy of whether p53 is acting in a "smart" or "dumb" manner. Here, we synthesize recent and past studies on p53 decoding of DNA sequence, chromatin context, and cellular signaling cascades to elicit variable cell fates critical in human development, homeostasis, and disease.

Gene regulation by the tumor suppressor p53 - The omics era

The transcription factor p53 is activated in response to a variety of cellular stresses and serves as a prominent and potent tumor suppressor. Since its discovery, we have sought to understand how p53 functions as both a transcription factor and a tumor suppressor. Two decades ago, the field of gene regulation entered the omics era and began to study the regulation of entire genomes. The omics perspective has greatly expanded our understanding of p53 functions and has begun to reveal its gene regulatory network. In this mini-review, I discuss recent insights into the p53 transcriptional program from high-throughput analyses.

p53 rapidly restructures 3D chromatin organization to trigger a transcriptional response

Activation of the p53 tumor suppressor triggers a transcriptional program to control cellular response to stress. However, the molecular mechanisms by which p53 controls gene transcription are not completely understood. Here, we uncover the critical role of spatio-temporal genome architecture in this process. We demonstrate that p53 drives direct and indirect changes in genome compartments, topologically associating domains, and DNA loops prior to one hour of its activation, which escort the p53 transcriptional program. Focusing on p53-bound enhancers, we report 340 genes directly regulated by p53 over a median distance of 116 kb, with 74% of these genes not previously identified. Finally, we showcase that p53 controls transcription of distal genes through newly formed and pre-existing enhancer-promoter loops in a cohesin dependent manner. Collectively, our findings demonstrate a previously unappreciated architectural role of p53 as regulator at distinct topological layers and provide a reliable set of new p53 direct target genes that may help designs of cancer therapies.

Cell fate regulation governed by p53: Friends or reversible foes in cancer therapy

Cancer is a leading cause of death worldwide. Targeted therapies aimed at key oncogenic driver mutations in combination with chemotherapy and radiotherapy as well as immunotherapy have benefited cancer patients considerably. Tumor protein p53 (TP53), a crucial tumor suppressor gene encoding p53, regulates numerous downstream genes and cellular phenotypes in response to various stressors. The affected genes are involved in diverse processes, including cell cycle arrest, DNA repair, cellular senescence, metabolic homeostasis, apoptosis, and autophagy. However, accumulating recent studies have continued to reveal novel and unexpected functions of p53 in governing the fate of tumors, for example, functions in ferroptosis, immunity, the tumor microenvironment and microbiome metabolism. Among the possibilities, the evolutionary plasticity of p53 is the most controversial, partially due to the dizzying array of biological functions that have been attributed to different regulatory mechanisms of p53 signaling. Nearly 40 years after its discovery, this key tumor suppressor remains somewhat enigmatic. The intricate and diverse functions of p53 in regulating cell fate during cancer treatment are only the tip of the iceberg with respect to its equally complicated structural biology, which has been painstakingly revealed. Additionally, TP53 mutation is one of the most significant genetic alterations in cancer, contributing to rapid cancer cell growth and tumor progression. Here, we summarized recent advances that implicate altered p53 in modulating the response to various cancer therapies, including chemotherapy, radiotherapy, and immunotherapy. Furthermore, we also discussed potential strategies for targeting p53 as a therapeutic option for cancer.

Orthogonal inducible control of Cas13 circuits enables programmable RNA regulation in mammalian cells

RNA plays an indispensable role in mammalian cell functions. Cas13, a class of RNA-guided ribonuclease, is a flexible tool for modifying and regulating coding and non-coding RNAs, with enormous potential for creating new cell functions. However, the lack of control over Cas13 activity has limited its cell engineering capability. Here, we present the CRISTAL (Control of RNA with Inducible SpliT CAs13 Orthologs and Exogenous Ligands) platform. CRISTAL is powered by a collection (10 total) of orthogonal split inducible Cas13 effectors that can be turned ON or OFF via small molecules in multiple cell types, providing precise temporal control. Also, we engineer Cas13 logic circuits that can respond to endogenous signaling and exogenous small molecule inputs. Furthermore, the orthogonality, low leakiness, and high dynamic range of our inducible Cas13d and Cas13b enable the design and construction of a robust incoherent feedforward loop, leading to near-perfect and tunable adaptation response. Finally, using our inducible Cas13 effectors, we achieve simultaneous multiplexed control of multiple genes in vitro and in mice. Together, our CRISTAL design represents a powerful platform for precisely regulating RNA dynamics to advance cell engineering and elucidate RNA biology.

Targeting class A GPCRs for hard tissue regeneration

G protein-coupled receptors (GPCRs) play important roles in various pathogeneses and physiological regulations. Owing to their functional diversity, GPCRs are considered one of the primary pharmaceutical targets. However, drugs targeting GPCRs have not been developed yet to regenerate hard tissues such as teeth and bones. Mesenchymal stromal cells (MSCs) have high proliferation and multi-lineage differentiation potential, which are essential for hard tissue regeneration. Here, we present a strategy for targeting class A GPCRs for hard tissue regeneration by promoting the differentiation of endogenous MSCs into osteogenic and odontogenic progenitor cells. Through in vitro screening targeted at class A GPCRs, we identified six target receptors (LPAR1, F2R, F2RL1, F2RL2, S1PR1, and ADORA2A) and candidate drugs with potent biomineralization effects. Through a combination of profiling whole transcriptome and accessible chromatin regions, we identified that p53 acts as a key transcriptional activator of genes that modulate the biomineralization process. Moreover, the therapeutic potential of class A GPCR-targeting drugs was demonstrated in tooth pulpotomy and calvarial defect models. The selected drugs revealed potent regenerative effects in both tooth and bone defects, represented by newly formed highly mineralized regions. Consequently, this study provides translational evidence for a new regenerative strategy for damaged hard tissue.

Energy metabolism modulates the regulatory impact of activators on gene expression

Gene expression is a regulated process fueled by ATP consumption. Therefore, regulation must be coupled to constraints imposed by the level of energy metabolism. Here, we explore this relationship both theoretically and experimentally. A stylized mathematical model predicts that activators of gene expression have variable impact depending on metabolic rate. Activators become less essential when metabolic rate is reduced and more essential when metabolic rate is enhanced. We find that, in the Drosophila eye, expression dynamics of the yan gene are less affected by loss of EGFR-mediated activation when metabolism is reduced, and the opposite effect is seen when metabolism is enhanced. The effects are also seen at the level of pattern regularity in the adult eye, where loss of EGFR-mediated activation is mitigated by lower metabolism. We propose that gene activation is tuned by energy metabolism to allow for faithful expression dynamics in the face of variable metabolic conditions.

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.

A molecular mechanism for the "digital" response of p53 to stress

The tumor suppressor protein p53 accumulates in response to cellular stress and consequently orchestrates the expression of multiple genes in a p53-level and time-dependent manner to overcome stress consequences, for which a molecular mechanism is currently unknown. Previously, we reported that DNA torsional flexibility distinguishes among p53 response elements (REs) and that transactivation at basal p53 levels is correlated with p53 REs flexibility. Here, we calculated the flexibility of ~200 p53 REs. By connecting functional outcomes of p53-target genes' activation to the calculated flexibility of their REs, we show that genes known to belong to pathways that are activated rapidly upon stress contain REs that are significantly more flexible relative to REs of genes known to be involved in pathways that are activated later in the response to stress. The global structural properties of several p53 REs belonging to different pathways were experimentally validated. Additionally, reporter-gene expression driven by flexible p53 REs occurred at lower p53 levels and with faster rates than expression from rigid REs. Furthermore, analysis of published endogenous mRNA levels of p53-target genes as a function of REs' flexibility showed that early versus late genes differ significantly in their flexibility properties of their REs and that highly flexible p53 REs enable high-activation level exclusively to early-response genes. Overall, we demonstrate that DNA flexibility of p53 REs contributes significantly to functional selectivity in the p53 system by facilitating the initial steps of p53-dependent target-genes expression, thereby contributing to survival versus death decisions in the p53 system.

Keratin 17 Impacts Global Gene Expression and Controls G2/M Cell Cycle Transition in Ionizing Radiation-Induced Skin Damage

Keratin 17 (K17) is a cytoskeletal protein that is part of the intermediate filaments in epidermal keratinocytes. In K17-/- mice, ionizing radiation induced more severe hair follicle damage, whereas the epidermal inflammatory response was attenuated compared with that in wild-type mice. Both p53 and K17 have a major impact on global gene expression because over 70% of the differentially expressed genes in the skin of wild-type mice showed no expression change in p53-/- or K17-/- skin after ionizing radiation. K17 does not interfere with the dynamics of p53 activation; rather, global p53 binding in the genome is altered in K17-/- mice. The absence of K17 leads to aberrant cell cycle progression and mitotic catastrophe in epidermal keratinocytes, which is due to nuclear retention, thus reducing the degradation of B-Myb, a key regulator of the G2/M cell cycle transition. These results expand our understanding of the role of K17 in regulating global gene expression and ionizing radiation-induced skin damage.

Proliferation history and transcription factor levels drive direct conversion

The sparse and stochastic nature of reprogramming has obscured our understanding of how transcription factors drive cells to new identities. To overcome this limit, we developed a compact, portable reprogramming system that increases direct conversion of fibroblasts to motor neurons by two orders of magnitude. We show that subpopulations with different reprogramming potentials are distinguishable by proliferation history. By controlling for proliferation history and titrating each transcription factor, we find that conversion correlates with levels of the pioneer transcription factor Ngn2, whereas conversion shows a biphasic response to Lhx3. Increasing the proliferation rate of adult human fibroblasts generates morphologically mature, induced motor neurons at high rates. Using compact, optimized, polycistronic cassettes, we generate motor neurons that graft with the murine central nervous system, demonstrating the potential for in vivo therapies.

Pulsed stimuli entrain p53 to synchronize single cells and modulate cell-fate determination

Entrainment to an external stimulus enables a synchronized oscillatory response across a population of cells, increasing coherent responses by reducing cell-to-cell heterogeneity. It is unclear whether the property of entrainability extends to systems where responses are intrinsic to the individual cell, rather than dependent on coherence across a population of cells. Using a combination of mathematical modeling, time-lapse fluorescence microscopy, and single-cell tracking, we demonstrated that p53 oscillations triggered by DNA double-strand breaks (DSBs) can be entrained with a periodic damage stimulus, despite such synchrony not known to function in effective DNA damage responses. Surprisingly, p53 oscillations were experimentally entrained over a wider range of DSB frequencies than predicted by an established computational model for the system. We determined that recapitulating the increased range of entrainment frequencies required, non-intuitively, a less robust oscillator and wider steady-state valley on the energy landscape. Further, we show that p53 entrainment can lead to altered expression dynamics of downstream targets responsible for cell fate in a manner dependent on target mRNA stability. Overall, this study demonstrates that entrainment can occur in a biological oscillator despite the apparent lack of an evolutionary advantage conferred through synchronized responses and highlights the potential of externally entraining p53 dynamics to reduce cellular variability and synchronize cell-fate responses for therapeutic outcomes.

Energy metabolism modulates the regulatory impact of activators on gene expression

Gene expression is a regulated process fueled by ATP consumption. Therefore, regulation must be coupled to constraints imposed by the level of energy metabolism. Here, we explore this relationship both theoretically and experimentally. A stylized mathematical model predicts that activators of gene expression have variable impact depending on metabolic rate. Activators become less essential when metabolic rate is reduced and more essential when metabolic rate is enhanced. We find that in the Drosophila eye, expression dynamics of the yan gene are less affected by loss of EGFR-mediated activation when metabolism is reduced, and the opposite effect is seen when metabolism is enhanced. The effects are also seen at the level of pattern regularity in the adult eye, where loss of EGFR-mediated activation is mitigated by lower metabolism. We propose that gene activation is tuned by energy metabolism to allow for faithful expression dynamics in the face of variable metabolic conditions.

Genome-wide CRISPR screening reveals ADCK3 as a key regulator in sensitizing endometrial carcinoma cells to MPA therapy

BACKGROUND

The effectiveness of conservative treatment of endometrial carcinoma (EC) with oral progesterone therapy, such as medroxyprogesterone acetate (MPA), can be blunted due to primary or acquired resistance, but the underlying mechanisms remain incompletely defined.

METHODS

Genome-wide CRISPR screening was performed to identify potential regulators in response to MPA in Ishikawa cells. Crystal violet staining, RT-qPCR, western blotting, ChIP-qPCR and luciferase assays were employed to elucidate the p53-AarF domain-containing kinase 3 (ADCK3) regulatory axis and its roles in sensitizing EC cells to MPA treatment.

RESULTS

ADCK3 is identified as a previously unrecognized regulator in response to MPA in EC cells. Loss of ADCK3 in EC cells markedly alleviated MPA-induced cell death. Mechanistically, loss of ADCK3 primarily suppresses MPA-mediated ferroptosis by abrogating arachidonate 15-lipoxygenase (ALOX15) transcriptional activation. Moreover, we validated ADCK3 as a direct downstream target of the tumor suppressor p53 in EC cells. By stimulating the p53-ADCK3 axis, the small-molecule compound Nutlin3A synergized with MPA to efficiently inhibit EC cell growth.

CONCLUSIONS

Our findings reveal ADCK3 as a key regulator of EC cells in response to MPA and shed light on a potential strategy for conservative EC treatment by activating the p53-ADCK3 axis to sensitize MPA-mediated cell death.

p53-dependent DNA repair during the DNA damage response requires actin nucleation by JMY

The tumour suppressor p53 is a nuclear transcription factor with key roles during DNA damage to enable a variety of cellular responses including cell cycle arrest, apoptosis and DNA repair. JMY is an actin nucleator and DNA damage-responsive protein whose sub-cellular localisation is responsive to stress and during DNA damage JMY undergoes nuclear accumulation. To gain an understanding of the wider role for nuclear JMY in transcriptional regulation, we performed transcriptomics to identify JMY-mediated changes in gene expression during the DNA damage response. We show that JMY is required for effective regulation of key p53 target genes involved in DNA repair, including XPC, XRCC5 (Ku80) and TP53I3 (PIG3). Moreover, JMY depletion or knockout leads to increased DNA damage and nuclear JMY requires its Arp2/3-dependent actin nucleation function to promote the clearance of DNA lesions. In human patient samples a lack of JMY is associated with increased tumour mutation count and in cells results in reduced cell survival and increased sensitivity to DNA damage response kinase inhibition. Collectively, we demonstrate that JMY enables p53-dependent DNA repair under genotoxic stress and suggest a role for actin in JMY nuclear activity during the DNA damage response.

Orthogonal inducible control of Cas13 circuits enables programmable RNA regulation in mammalian cells

RNA plays an indispensable role in mammalian cell functions. Cas13, a class of RNA-guided ribonuclease, is a flexible tool for modifying and regulating coding and non-coding RNAs, with enormous potential for creating new cell functions. However, the lack of control over Cas13 activity has limited its cell engineering capability. Here, we present the CRISTAL ( C ontrol of R NA with Inducible S pli T C A s13 Orthologs and Exogenous L igands) platform. CRISTAL is powered by a collection (10 total) of orthogonal split inducible Cas13s that can be turned ON or OFF via small molecules in multiple cell types, providing precise temporal control. Also, we engineered Cas13 logic circuits that can respond to endogenous signaling and exogenous small molecule inputs. Furthermore, the orthogonality, low leakiness, and high dynamic range of our inducible Cas13d and Cas13b enable the design and construction of a robust incoherent feedforward loop, leading to near-perfect and tunable adaptation response. Finally, using our inducible Cas13s, we achieve simultaneous multiplexed control of multiple genes in vitro and in mice. Together, our CRISTAL design represents a powerful platform for precisely regulating RNA dynamics to advance cell engineering and elucidate RNA biology.

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.

CRISPR screens reveal genetic determinants of PARP inhibitor sensitivity and resistance in prostate cancer

Prostate cancer harboring BRCA1/2 mutations are often exceptionally sensitive to PARP inhibitors. However, genomic alterations in other DNA damage response genes have not been consistently predictive of clinical response to PARP inhibition. Here, we perform genome-wide CRISPR-Cas9 knockout screens in BRCA1/2-proficient prostate cancer cells and identify previously unknown genes whose loss has a profound impact on PARP inhibitor response. Specifically, MMS22L deletion, frequently observed (up to 14%) in prostate cancer, renders cells hypersensitive to PARP inhibitors by disrupting RAD51 loading required for homologous recombination repair, although this response is TP53-dependent. Unexpectedly, loss of CHEK2 confers resistance rather than sensitivity to PARP inhibition through increased expression of BRCA2, a target of CHEK2-TP53-E2F7-mediated transcriptional repression. Combined PARP and ATR inhibition overcomes PARP inhibitor resistance caused by CHEK2 loss. Our findings may inform the use of PARP inhibitors beyond BRCA1/2-deficient tumors and support reevaluation of current biomarkers for PARP inhibition in prostate cancer.

Coordination of MAPK and p53 dynamics in the cellular responses to DNA damage and oxidative stress

In response to different cellular stresses, the transcription factor p53 undergoes different dynamics. p53 dynamics, in turn, control cell fate. However, distinct stresses can generate the same p53 dynamics but different cell fate outcomes, suggesting integration of dynamic information from other pathways is important for cell fate regulation. To determine how MAPK activities affect p53-mediated responses to DNA breaks and oxidative stress, we simultaneously tracked p53 and either ERK, JNK, or p38 activities in single cells. While p53 dynamics were comparable between the stresses, cell fate outcomes were distinct. Combining MAPK dynamics with p53 dynamics was important for distinguishing between the stresses and for generating temporal ordering of cell fate pathways. Furthermore, cross-talk between MAPKs and p53 controlled the balance between proliferation and cell death. These findings provide insight into how cells integrate signaling pathways with distinct temporal patterns of activity to encode stress specificity and drive different cell fate decisions.

Divergent Molecular and Cellular Responses to Low and High-Dose Ionizing Radiation

Cancer risk after ionizing radiation (IR) is assumed to be linear with the dose; however, for low doses, definite evidence is lacking. Here, using temporal multi-omic systems analyses after a low (LD; 0.1 Gy) or a high (HD; 1 Gy) dose of X-rays, we show that, although the DNA damage response (DDR) displayed dose proportionality, many other molecular and cellular responses did not. Phosphoproteomics uncovered a novel mode of phospho-signaling via S12-PPP1R7, and large-scale dephosphorylation events that regulate mitotic exit control in undamaged cells and the G2/M checkpoint upon IR in a dose-dependent manner. The phosphoproteomics of irradiated DNA double-strand breaks (DSBs) repair-deficient cells unveiled extended phospho-signaling duration in either a dose-dependent (DDR signaling) or independent (mTOR-ERK-MAPK signaling) manner without affecting signal magnitude. Nascent transcriptomics revealed the transcriptional activation of genes involved in NRF2-regulated antioxidant defense, redox-sensitive ERK-MAPK signaling, glycolysis and mitochondrial function after LD, suggesting a prominent role for reactive oxygen species (ROS) in molecular and cellular responses to LD exposure, whereas DDR genes were prominently activated after HD. However, how and to what extent the observed dose-dependent differences in molecular and cellular responses may impact cancer development remain unclear, as the induction of chromosomal damage was found to be dose-proportional (10-200 mGy).

Characterization of chicken p53 transcriptional function via parallel genome-wide chromatin occupancy and gene expression analysis

The tumor suppressor p53, which acts primarily as a transcription factor, can regulate infections from various viruses in chickens. However, the underlying mechanisms of the antiviral functions of chicken p53 (chp53) remain unclear due to the lack of detailed information on its transcriptional regulation. Here, to gain comprehensive insights into chp53 transcriptional regulatory function in a global and unbiased manner, we determined the genome-wide chromatin occupancy of chp53 by chromatin immunoprecipitation, which was followed by sequencing and chp53-mediated gene expression profile by RNA sequencing using chemically immortalized leghorn male hepatoma (LMH) cells with ectopic expression of chp53 as the model. The integrated parallel genome-wide chromatin occupancy and gene expression analysis characterized chp53 chromatin occupancy and identified 754 direct target genes of chp53. Furthermore, functional annotation and cross-species comparative biological analyses revealed the conserved key biological functions and DNA binding motifs of p53 between chickens and humans, which may be due to the consensus amino acid sequence and structure of p53 DNA-binding domains. The present study, to our knowledge, provides the first comprehensive characterization of the chp53 transcriptional regulatory network, and can possibly help to improve our understanding of p53 transcriptional regulatory mechanisms and their antiviral functions in chickens.

Quantifying the phenotypic information in mRNA abundance

Quantifying the dependency between mRNA abundance and downstream cellular phenotypes is a fundamental open problem in biology. Advances in multimodal single-cell measurement technologies provide an opportunity to apply new computational frameworks to dissect the contribution of individual genes and gene combinations to a given phenotype. Using an information theory approach, we analyzed multimodal data of the expression of 83 genes in the Ca2+ signaling network and the dynamic Ca2+ response in the same cell. We found that the overall expression levels of these 83 genes explain approximately 60% of Ca2+ signal entropy. The average contribution of each single gene was 17%, revealing a large degree of redundancy between genes. Using different heuristics, we estimated the dependency between the size of a gene set and its information content, revealing that on average, a set of 53 genes contains 54% of the information about Ca2+ signaling. Our results provide the first direct quantification of information content about complex cellular phenotype that exists in mRNA abundance measurements.

Profiling transcription factor activity dynamics using intronic reads in time-series transcriptome data

Activities of transcription factors (TFs) are temporally modulated to regulate dynamic cellular processes, including development, homeostasis, and disease. Recent developments of bioinformatic tools have enabled the analysis of TF activities using transcriptome data. However, because these methods typically use exon-based target expression levels, the estimated TF activities have limited temporal accuracy. To address this, we proposed a TF activity measure based on intron-level information in time-series RNA-seq data, and implemented it to decode the temporal control of TF activities during dynamic processes. We showed that TF activities inferred from intronic reads can better recapitulate instantaneous TF activities compared to the exon-based measure. By analyzing public and our own time-series transcriptome data, we found that intron-based TF activities improve the characterization of temporal phasing of cycling TFs during circadian rhythm, and facilitate the discovery of two temporally opposing TF modules during T cell activation. Collectively, we anticipate that the proposed approach would be broadly applicable for decoding global transcriptional architecture during dynamic processes.

An oscillatory network controlling self-renewal of skeletal muscle stem cells

The balance between proliferation and differentiation of muscle stem cells is tightly controlled, ensuring the maintenance of a cellular pool needed for muscle growth and repair. Muscle stem cells can proliferate, they can generate differentiating cells, or they self-renew to produce new stem cells. Notch signaling plays a crucial role in this process. Recent studies revealed that expression of the Notch effector HES1 oscillates in activated muscle stem cells. The oscillatory expression of HES1 periodically represses transcription from the genes encoding the myogenic transcription factor MYOD and the Notch ligand DLL1, thereby driving MYOD and DLL1 oscillations. This oscillatory network allows muscle progenitor cells and activated muscle stem cells to remain in a proliferative and 'undecided' state, in which they can either differentiate or self-renew. When HES1 is downregulated, MYOD oscillations become unstable and are replaced by sustained expression, which drives the cells into terminal differentiation. During development and regeneration, proliferating stem cells contact each other and the stability of the oscillatory expression depends on regular DLL1 inputs provided by neighboring cells. In such communities of cells that receive and provide Notch signals, the appropriate timing of DLL1 inputs is important, as sustained DLL1 cannot replace oscillatory DLL1. Thus, in cell communities, DLL1 oscillations ensure the appropriate balance between self-renewal and differentiation. In summary, oscillations in myogenic cells are an important example of dynamic gene expression determining cell fate.

Progress and challenges in understanding the regulation and function of p53 dynamics

The dynamics of p53 expression provide a mechanism to increase differentiation between cellular stresses and specificity in appropriate responses. Here, we review recent advances in our understanding of the molecular mechanisms regulating p53 dynamics and the functions of the dynamics in the regulation of p53-dependent cell stress responses. We also compare dynamic encoding in the p53 system with that found in other important cell signaling systems, many of which can interact with the p53 network. Finally, we highlight some of the current challenges in understanding dynamic cell signaling within a larger cellular network context.

At a Crossroads to Cancer: How p53-Induced Cell Fate Decisions Secure Genome Integrity

P53 is known as the most critical tumor suppressor and is often referred to as the guardian of our genome. More than 40 years after its discovery, we are still struggling to understand all molecular details on how this transcription factor prevents oncogenesis or how to leverage current knowledge about its function to improve cancer treatment. Multiple cues, including DNA-damage or mitotic errors, can lead to the stabilization and nuclear translocation of p53, initiating the expression of multiple target genes. These transcriptional programs may be cell-type- and stimulus-specific, as is their outcome that ultimately imposes a barrier to cellular transformation. Cell cycle arrest and cell death are two well-studied consequences of p53 activation, but, while being considered critical, they do not fully explain the consequences of p53 loss-of-function phenotypes in cancer. Here, we discuss how mitotic errors alert the p53 network and give an overview of multiple ways that p53 can trigger cell death. We argue that a comparative analysis of different types of p53 responses, elicited by different triggers in a time-resolved manner in well-defined model systems, is critical to understand the cell-type-specific cell fate induced by p53 upon its activation in order to resolve the remaining mystery of its tumor-suppressive function.

Joint single-cell measurements of nuclear proteins and RNA in vivo

Identifying gene-regulatory targets of nuclear proteins in tissues is a challenge. Here we describe intranuclear cellular indexing of transcriptomes and epitopes (inCITE-seq), a scalable method that measures multiplexed intranuclear protein levels and the transcriptome in parallel across thousands of nuclei, enabling joint analysis of transcription factor (TF) levels and gene expression in vivo. We apply inCITE-seq to characterize cell state-related changes upon pharmacological induction of neuronal activity in the mouse brain. Modeling gene expression as a linear combination of quantitative protein levels revealed genome-wide associations of each TF and recovered known gene targets. TF-associated genes were coexpressed as distinct modules that each reflected positive or negative TF levels, showing that our approach can disentangle relative putative contributions of TFs to gene expression and add interpretability to inferred gene networks. inCITE-seq can illuminate how combinations of nuclear proteins shape gene expression in native tissue contexts, with direct applications to solid or frozen tissues and clinical specimens.

p53-Dependent Repression: DREAM or Reality?

Simple Summary

The tumor suppressor p53 is a complex cell signaling hub encompassing multiple transcription programs and governs a vast repertoire of biological responses. However, despite several decades of research, how p53 selects one program over another is still elusive. Recent attempts have used meta-analyses of p53 ChIP-seq data to determine the core p53 transcriptional program, conserved across different models and stimuli. This review highlights the complexity of the multiple layers of p53 regulation and the context specificity of p53 target genes. More specifically, we discuss the controversy over the mechanisms of p53-dependent transcriptional repression and its potential role in the flexibility of p53 response.

Abstract

p53 is a major tumor suppressor that integrates diverse types of signaling in mammalian cells. In response to a broad range of intra- or extra-cellular stimuli, p53 controls the expression of multiple target genes and elicits a vast repertoire of biological responses. The exact code by which p53 integrates the various stresses and translates them into an appropriate transcriptional response is still obscure. p53 is tightly regulated at multiple levels, leading to a wide diversity in p53 complexes on its target promoters and providing adaptability to its transcriptional program. As p53-targeted therapies are making their way into clinics, we need to understand how to direct p53 towards the desired outcome (i.e., cell death, senescence or other) selectively in cancer cells without affecting normal tissues or the immune system. While the core p53 transcriptional program has been proposed, the mechanisms conferring a cell type- and stimuli-dependent transcriptional outcome by p53 require further investigations. The mechanism by which p53 localizes to repressed promoters and manages its co-repressor interactions is controversial and remains an important gap in our understanding of the p53 cistrome. We hope that our review of the recent literature will help to stimulate the appreciation and investigation of largely unexplored p53-mediated repression.

Technologies for profiling the impact of genomic variants on transcription factor binding

Transcription factors (TFs) bind DNA in a sequence-specific manner and thereby regulate target gene expression. TF binding and its regulatory activity is highly context dependent, and is not only determined by specific cell types or differentiation stages but also relies on other regulatory mechanisms, such as DNA and chromatin modifications. Interactions between TFs and their DNA binding sites are critical mediators of phenotypic variation and play important roles in the onset of disease. A continuously growing number of studies therefore attempts to elucidate TF:DNA interactions to gain knowledge about regulatory mechanisms and disease-causing variants. Here we summarize how TF-binding characteristics and the impact of variants can be investigated, how bioinformatic tools can be used to analyze and predict TF:DNA binding, and what additional information can be obtained from the TF protein structure.

A Hill type equation can predict target gene expression driven by p53 pulsing

Many factors determine target gene expression dynamics under p53 pulsing. In this study, I sought to determine the mechanism by which duration, frequency, binding affinity and maximal transcription rate affect the expression dynamics of target genes. Using an analytical method to solve a simple model, I found that the fold change of target gene expression increases relative to the number of p53 pulses, and the optimal frequency, 0.18 h-1 , from two real p53 pulses drives the maximal fold change with a decay rate of 0.18 h-1 . Moreover, p53 pulses may also lead to a higher fold change than sustained p53. Finally, I discovered that a Hill-type equation, including these effect factors, can characterise target gene expression. The average error between the theoretical predictions and experiments was 23%. Collectively, this equation advances the understanding of transcription factor dynamics, where duration and frequency play a significant role in the fine regulation of target gene expression with higher binding affinity.

The Role of p53 Signaling in Colorectal Cancer

The transcription factor p53 functions as a critical tumor suppressor by orchestrating a plethora of cellular responses such as DNA repair, cell cycle arrest, cellular senescence, cell death, cell differentiation, and metabolism. In unstressed cells, p53 levels are kept low due to its polyubiquitination by the E3 ubiquitin ligase MDM2. In response to various stress signals, including DNA damage and aberrant growth signals, the interaction between p53 and MDM2 is blocked and p53 becomes stabilized, allowing p53 to regulate a diverse set of cellular responses mainly through the transactivation of its target genes. The outcome of p53 activation is controlled by its dynamics, its interactions with other proteins, and post-translational modifications. Due to its involvement in several tumor-suppressing pathways, p53 function is frequently impaired in human cancers. In colorectal cancer (CRC), the TP53 gene is mutated in 43% of tumors, and the remaining tumors often have compromised p53 functioning because of alterations in the genes encoding proteins involved in p53 regulation, such as ATM (13%) or DNA-PKcs (11%). TP53 mutations in CRC are usually missense mutations that impair wild-type p53 function (loss-of-function) and that even might provide neo-morphic (gain-of-function) activities such as promoting cancer cell stemness, cell proliferation, invasion, and metastasis, thereby promoting cancer progression. Although the first compounds targeting p53 are in clinical trials, a better understanding of wild-type and mutant p53 functions will likely pave the way for novel CRC therapies.

Connecting Timescales in Biology: Can Early Dynamical Measurements Predict Long-Term Outcomes?

Prediction of long-term outcomes from short-term measurements remains a fundamental challenge. Quantitative assessment of signaling dynamics, and the resulting transcriptomic and proteomic responses, has yielded fundamental insights into cellular outcomes. However, the utility of these measurements is limited by their short timescale (hours to days), while the consequences of these events frequently unfold over longer timescales. Here, we discuss the predictive power of static and dynamic measurements, drawing examples from fields that have harnessed the predictive capabilities of such measurements. We then explore potential approaches to close this timescale gap using complementary measurements and computational approaches, focusing on the example of dynamic measurements of signaling factors and their impacts on cellular outcomes.

Temporal gene regulation by p53 is associated with the rotational setting of its binding sites in nucleosomes

The tumor suppressor protein p53 is a DNA-binding transcription factor (TF) that, once activated, coordinates the expression of thousands of target genes. Increased p53 binding to gene promoters occurs shortly after p53 activation. Intriguingly, gene transcription exhibits differential kinetics with some genes being induced early (early genes) and others being induced late (late genes). To understand pre-binding factors contributing to the temporal gene regulation by p53, we performed time-course RNA sequencing experiments in human colon cancer cell line HCT116 treated with fluorouracil to identify early and late genes. Published p53 ChIP fragments co-localized with the early or late genes were used to uncover p53 binding sites (BS). We demonstrate that the BS associated with early genes are clustered around gene starts with decreased nucleosome occupancy. DNA analysis shows that these BS are likely exposed on nucleosomal surface if wrapped into nucleosomes, thereby facilitating stable interactions with and fast induction by p53. By contrast, p53 BS associated with late genes are distributed uniformly across the genes with increased nucleosome occupancy. Predicted rotational settings of these BS show limited accessibility. We therefore propose a hypothetical model in which the BS are fully, partially or not accessible to p53 in the nucleosomal context. The partial accessibility of the BS allows subunits of a p53 tetramer to bind, but the resulting p53-DNA complex may not be stable enough to recruit cofactors, which leads to delayed induction. Our work highlights the importance of DNA conformations of p53 BS in gene expression dynamics.

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.

Manipulation of the human tRNA pool reveals distinct tRNA sets that act in cellular proliferation or cell cycle arrest

Different subsets of the tRNA pool in human cells are expressed in different cellular conditions. The 'proliferation-tRNAs' are induced upon normal and cancerous cell division, while the 'differentiation-tRNAs' are active in non-dividing, differentiated cells. Here we examine the essentiality of the various tRNAs upon cellular growth and arrest. We established a CRISPR-based editing procedure with sgRNAs that each target a tRNA family. We measured tRNA essentiality for cellular growth and found that most proliferation-tRNAs are essential compared to differentiation- tRNAs in rapidly growing cell lines. Yet in more slowly dividing lines, the differentiation-tRNAs were more essential. In addition, we measured the essentiality of each tRNA family upon response to cell cycle arresting signals. Here we detected a more complex behavior with both proliferation-tRNAs and differentiation tRNAs showing various levels of essentiality. These results provide the so-far most comprehensive functional characterization of human tRNAs with intricate roles in various cellular states.

Life and Death Decision-Making by p53 and Implications for Cancer Immunotherapy

The tumor-suppressor protein p53 is mutated in approximately half of all cancers, whereas the p53 signaling network is perturbed in almost all cancers. In response to different stress stimuli, p53 selectively activates genes to elicit a cell survival or cell death response. How p53 makes the decision between life and death remains a fascinating question and an exciting field of research. Understanding how this decision is made has major implications for improving cancer treatments, particularly in recently evolved immune checkpoint inhibition therapy. We highlight progress and challenges in understanding the mechanisms governing the p53 life and death decision-making process, and discuss how this decision is relevant to immune system regulation. Finally, we discuss how knowledge of the p53 pro-survival and pro-death decision node can be applied to optimize immune checkpoint inhibitor therapy for cancer treatment.

Cell-Cycle-Dependent ERK Signaling Dynamics Direct Fate Specification in the Mammalian Preimplantation Embryo

Despite the noisy nature of single cells, multicellular organisms robustly generate different cell types from one zygote. This process involves dynamic cross regulation between signaling and gene expression that is difficult to capture with fixed-cell approaches. To study signaling dynamics and fate specification during preimplantation development, we generated a transgenic mouse expressing the ERK kinase translocation reporter and measured ERK activity in single cells of live embryos. Our results show primarily active ERK in both the inner cell mass and trophectoderm cells due to fibroblast growth factor (FGF) signaling. Strikingly, a subset of mitotic events results in a short pulse of ERK inactivity in both daughter cells that correlates with elevated endpoint NANOG levels. Moreover, endogenous tagging of Nanog in embryonic stem cells reveals that ERK inhibition promotes enhanced stabilization of NANOG protein after mitosis. Our data show that cell cycle, signaling, and differentiation are coordinated during preimplantation development.

IFI16 promotes human embryonic stem cell trilineage specification through interaction with p53

Transcriptional regulation plays an essential role in the self-renewal and differentiation of human embryonic stem cells (hESCs). However, how external signals disrupt the self-renewal regulatory network and further drive hESC differentiation remains largely unknown. Here, we found the immune regulative protein, gamma-interferon-inducible protein 16 (IFI16) was involved in the regulation of both self-renewal and differentiation gene expression during hESC trilineage specification through interaction with p53. IFI16 expression levels were upregulated through JNK activation. IFI16 knockdown delayed the downregulation of self-renewal gene expression and suppressed the upregulation of differentiation gene expression, while IFI16 overexpression accelerated trilineage specification. Furthermore, IFI16 stabilized p53-binding in the genome through IFI16-p53 interaction and differentially regulated self-renewal and differentiation gene expression. Together, our results suggest a particular role of IFI16 in differential gene expression regulation during trilineage specification of hESCs in a manner that is dependent on the genome-wide profile of p53-binding directed by IFI16-p53 interaction.

Regulation of circadian rhythm and sleep by miR-375-timeless interaction in Drosophila

MicroRNAs are important coordinators of circadian regulation that mediate the fine-tuning of gene expression. Although many studies have shown the effects of individual miRNAs on the circadian clock, the global functional miRNA-mRNA interaction network involved in the circadian system remains poorly understood. Here, we used CLEAR (Covalent Ligation of Endogenous Argonaute-bound RNAs)-CLIP (Cross-Linking and Immuno-Precipitation) to explore the regulatory functions of miRNAs in the circadian system by comparing the miRNA-mRNA interactions between Drosophila wild-type strain W¹¹¹⁸ and a mutant of the key circadian transcriptional regulator Clock (Clk^jrk ). This experimental approach unambiguously identified tens of thousands of miRNA-mRNA interactions in both the head and body. The miRNA-mRNA interactome showed dramatic changes in the Clk^jrk flies. Particularly, among ~300 miRNA-mRNA circadian relevant interactions, multiple interactions involving core clock genes pdp1, tim, and vri displayed distinct changes as a result of the Clk mutation. Based on the CLEAR-CLIP analysis, we found a novel regulation of the circadian rhythm and sleep by the miR-375-timeless interaction. The results indicated that Clk disruption abolished normal rhythmic expression of miR-375 and the functional regulation occurred in the l-LNv neurons, where miR-375 modulated the circadian rhythm and sleep via targeting timeless. This work provides the first global view of miRNA regulation in the circadian rhythm.

Rucaparib Treatment Alters p53 Oscillations in Single Cells to Enhance DNA-Double-Strand-Break-Induced Cell Cycle Arrest

DNA double strand breaks induce oscillatory expression of the transcription factor p53 that is dependent on ataxia telangiectasia mutated (ATM) activity and the rate of double strand break resolution. Although p53 dynamics are known to play a role in the regulation of cell fate determination, the consequences of the variability in dynamics associated with differences in repair rates and utilized repair pathways are unknown. Using single-cell time-lapse microscopy, we found that disruption of specific repair pathways has distinct impacts on p53 dynamics. The small-molecule rucaparib, an inhibitor of the alternative end-joining-associated protein poly (ADP-ribose) polymerase (PARP), increased p53 pulse duration, altering the temporal expression of multiple p53 target genes. As a result, combination treatments of the radiomimetic drug neocarzinostatin with rucaparib drove prolonged growth arrest beyond that of DNA damage alone. This study highlights how pharmacological manipulation of DNA repair pathways may be used to alter p53 dynamics to enhance therapeutic regimens.

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

Oscillatory Dynamics of p53-Mdm2 Circuit in Response to DNA Damage Caused by Ionizing Radiation

Although the dynamical behavior of the p53-Mdm2 loop has been extensively studied, the understanding of the mechanism underlying the regulation of this pathway still remains limited. Herein, we developed an integrated model with five basic components and three ubiquitous time delays for the p53-Mdm2 interaction in response to DNA damage following ionizing radiation (IR). We showed that a sufficient amount of activated ATM level can initiate the p53 oscillations with nearly the same amplitude over a wide range of the ATM level; a proper range of p53 level is also required for generating the oscillations, for too high or too low levels it would fail to generate the oscillations; and increased Mdm2 level leads to decreased amplitude of the p53 oscillation and reduced expression of the p53 activity. Moreover, we found that the negative feedback loop formed between p53 and nuclear Mdm2 plays a dominant role in determining the p53 dynamics, whereas when interaction strength of the negative feedback loop becomes weaker, the positive feedback loop formed between p53 and cytoplasmatic Mdm2 can induce different types of dynamics. Furthermore, we demonstrated that the total time delay required for protein production and nuclear translocation of Mdm2 can induce p53 oscillations even when the p53 level is at a certain stable high steady state or at a certain stable low steady state. In addition, the two important features of the oscillatory dynamics-amplitude and period-can be controlled by such time delay. These results are in agreement with multiple experimental observations and may enrich our understanding of the dynamics of the p53 network.

Pinning Control for the p53-Mdm2 Network Dynamics Regulated by p14ARF

p53 regulates the cellular response to genotoxic damage and prevents carcinogenic events. Theoretical and experimental studies state that the p53-Mdm2 network constitutes the core module of regulatory interactions activated by cellular stress induced by a variety of signaling pathways. In this paper, a strategy to control the p53-Mdm2 network regulated by p14ARF is developed, based on the pinning control technique, which consists into applying local feedback controllers to a small number of nodes (pinned ones) in the network. Pinned nodes are selected on the basis of their importance level in a topological hierarchy, their degree of connectivity within the network, and the biological role they perform. In this paper, two cases are considered. For the first case, the oscillatory pattern under gamma-radiation is recovered; afterward, as the second case, increased expression of p53 level is taken into account. For both cases, the control law is applied to p14ARF (pinned node based on a virtual leader methodology), and overexpressed Mdm2-mediated p53 degradation condition is considered as carcinogenic initial behavior. The approach in this paper uses a computational algorithm, which opens an alternative path to understand the cellular responses to stress, doing it possible to model and control the gene regulatory network dynamics in two different biological contexts. As the main result of the proposed control technique, the two mentioned desired behaviors are obtained.

A Switch in p53 Dynamics Marks Cells That Escape from DSB-Induced Cell Cycle Arrest

Cellular responses to stimuli can evolve over time, resulting in distinct early and late phases in response to a single signal. DNA damage induces a complex response that is largely orchestrated by the transcription factor p53, whose dynamics influence whether a damaged cell will arrest and repair the damage or will initiate cell death. How p53 responses and cellular outcomes evolve in the presence of continuous DNA damage remains unknown. Here, we have found that a subset of cells switches from oscillating to sustained p53 dynamics several days after undergoing damage. The switch results from cell cycle progression in the presence of damaged DNA, which activates the caspase-2-PIDDosome, a complex that stabilizes p53 by inactivating its negative regulator MDM2. This work defines a molecular pathway that is activated if the canonical checkpoints fail to halt mitosis in the presence of damaged DNA.

Hybrid analysis of gene dynamics predicts context-specific expression and offers regulatory insights

Motivation

To understand the regulatory pathways underlying diseases, studies often investigate the differential gene expression between genetically or chemically differing cell populations. Differential expression analysis identifies global changes in transcription and enables the inference of functional roles of applied perturbations. This approach has transformed the discovery of genetic drivers of disease and possible therapies. However, differential expression analysis does not provide quantitative predictions of gene expression in untested conditions. We present a hybrid approach, termed Differential Expression in Python (DiffExPy), that uniquely combines discrete, differential expression analysis with in silico differential equation simulations to yield accurate, quantitative predictions of gene expression from time-series data.

Results

To demonstrate the distinct insight provided by DiffExpy, we applied it to published, in vitro, time-series RNA-seq data from several genetic PI3K/PTEN variants of MCF10a cells stimulated with epidermal growth factor. DiffExPy proposed ensembles of several minimal differential equation systems for each differentially expressed gene. These systems provide quantitative models of expression for several previously uncharacterized genes and uncover new regulation by the PI3K/PTEN pathways. We validated model predictions on expression data from conditions that were not used for model training. Our discrete, differential expression analysis also identified SUZ12 and FOXA1 as possible regulators of specific groups of genes that exhibit late changes in expression. Our work reveals how DiffExPy generates quantitatively predictive models with testable, biological hypotheses from time-series expression data.

Availability and implementation

DiffExPy is available on GitHub (https://github.com/bagherilab/diffexpy).

Supplementary information

Supplementary data are available at Bioinformatics online.