Diverse Spatial Expression Patterns Emerge from Unified Kinetics of Transcriptional Bursting

Optogenetic dissection of transcriptional repression in a multicellular organism

Transcriptional control is fundamental to cellular function. However, despite knowing that transcription factors can repress or activate specific genes, how these functions are implemented at the molecular level has remained elusive, particularly in the endogenous context of developing animals. Here, we combine optogenetics, single-cell live-imaging, and mathematical modeling to study how a zinc-finger repressor, Knirps, induces switch-like transitions into long-lived quiescent states. Using optogenetics, we demonstrate that repression is rapidly reversible (~1 min) and memoryless. Furthermore, we show that the repressor acts by decreasing the frequency of transcriptional bursts in a manner consistent with an equilibrium binding model. Our results provide a quantitative framework for dissecting the in vivo biochemistry of eukaryotic transcriptional regulation.

A guide to studying 3D genome structure and dynamics in the kidney

The human genome is tightly packed into the 3D environment of the cell nucleus. Rapidly evolving and sophisticated methods of mapping 3D genome architecture have shed light on fundamental principles of genome organization and gene regulation. The genome is physically organized on different scales, from individual genes to entire chromosomes. Nuclear landmarks such as the nuclear envelope and nucleoli have important roles in compartmentalizing the genome within the nucleus. Genome activity (for example, gene transcription) is also functionally partitioned within this 3D organization. Rather than being static, the 3D organization of the genome is tightly regulated over various time scales. These dynamic changes in genome structure over time represent the fourth dimension of the genome. Innovative methods have been used to map the dynamic regulation of genome structure during important cellular processes including organism development, responses to stimuli, cell division and senescence. Furthermore, disruptions to the 4D genome have been linked to various diseases, including of the kidney. As tools and approaches to studying the 4D genome become more readily available, future studies that apply these methods to study kidney biology will provide insights into kidney function in health and disease.

Transcriptional bursting dynamics in gene expression

Gene transcription is a stochastic process that occurs in all organisms. Transcriptional bursting, a critical molecular dynamics mechanism, creates significant heterogeneity in mRNA and protein levels. This heterogeneity drives cellular phenotypic diversity. Currently, the lack of a comprehensive quantitative model limits the research on transcriptional bursting. This review examines various gene expression models and compares their strengths and weaknesses to guide researchers in selecting the most suitable model for their research context. We also provide a detailed summary of the key metrics related to transcriptional bursting. We compared the temporal dynamics of transcriptional bursting across species and the molecular mechanisms influencing these bursts, and highlighted the spatiotemporal patterns of gene expression differences by utilizing metrics such as burst size and burst frequency. We summarized the strategies for modeling gene expression from both biostatistical and biochemical reaction network perspectives. Single-cell sequencing data and integrated multiomics approaches drive our exploration of cutting-edge trends in transcriptional bursting mechanisms. Moreover, we examined classical methods for parameter estimation that help capture dynamic parameters in gene expression data, assessing their merits and limitations to facilitate optimal parameter estimation. Our comprehensive summary and review of the current transcriptional burst dynamics theories provide deeper insights for promoting research on the nature of cell processes, cell fate determination, and cancer diagnosis.

Nascent RNA kinetics with complex promoter architecture: Analytic results and parameter inference

Transcription is a stochastic process that involves several downstream operations which make it difficult to model and infer transcription kinetics from mature RNA numbers in individual cell. However, recent advances in single-cell technologies have enabled a more precise measurement of the fluctuations of nascent RNA that closely reflect transcription kinetics. In this paper we introduce a general stochastic model to mimic nascent RNA kinetics with complex promoter architecture. We derive the exact distribution and moments of nascent RNA using queuing theory techniques, which provide valuable insights into the effect of the molecular memory created by the multistep activation and deactivation on the stochastic kinetics of nascent RNA. Moreover, based on the analytical results, we develop a statistical method to infer the promoter memory from stationary nascent RNA distributions. Data analysis of synthetic data and a realistic example, the HIV-1 gene, verifies the validity of this inference method.

Transcription factor clusters as information transfer agents

Deciphering how genes interpret information from transcription factor (TFs) concentrations within the cell nucleus remains a fundamental question in gene regulation. Recent advancements have revealed the heterogeneous distribution of TF molecules, posing challenges to precisely decoding concentration signals. Using high-resolution single-cell imaging of the fluorescently tagged TF Bicoid in living Drosophila embryos, we show that Bicoid accumulation in submicron clusters preserves the spatial information of the maternal Bicoid gradient. These clusters provide precise spatial cues through intensity, size, and frequency. We further discover that gene targets of Bicoid, such as Hunchback and Eve, colocalize with these clusters in an enhancer binding affinity-dependent manner. Our modeling suggests that clustering offers a faster sensing mechanism for global nuclear concentrations than freely diffusing TF molecules detected by simple enhancers.

How Transcription Factor Clusters Shape the Transcriptional Landscape

In eukaryotic cells, gene transcription typically occurs in discrete periods of promoter activity, interspersed with intervals of inactivity. This pattern deviates from simple stochastic events and warrants a closer examination of the molecular interactions that activate the promoter. Recent studies have identified transcription factor (TF) clusters as key precursors to transcriptional bursting. Often, these TF clusters form at chromatin segments that are physically distant from the promoter, making changes in chromatin conformation crucial for promoter-TF cluster interactions. In this review, I explore the formation and constituents of TF clusters, examining how the dynamic interplay between chromatin architecture and TF clustering influences transcriptional bursting. Additionally, I discuss techniques for visualizing TF clusters and provide an outlook on understanding the remaining gaps in this field.

Gene activity fully predicts transcriptional bursting dynamics

Transcription commonly occurs in bursts, with alternating productive (ON) and quiescent (OFF) periods, governing mRNA production rates. Yet, how transcription is regulated through bursting dynamics remains unresolved. Here, we conduct real-time measurements of endogenous transcriptional bursting with single-mRNA sensitivity. Leveraging the diverse transcriptional activities in early fly embryos, we uncover stringent relationships between bursting parameters. Specifically, we find that the durations of ON and OFF periods are linked. Regardless of the developmental stage or body-axis position, gene activity levels predict individual alleles' average ON and OFF periods. Lowly transcribing alleles predominantly modulate OFF periods (burst frequency), while highly transcribing alleles primarily tune ON periods (burst size). These relationships persist even under perturbations of cis-regulatory elements or trans-factors and account for bursting dynamics measured in other species. Our results suggest a novel mechanistic constraint governing bursting dynamics rather than a modular control of distinct parameters by distinct regulatory processes.

Unified bursting strategies in ectopic and endogenous even-skipped expression patterns

Transcription often occurs in bursts as gene promoters switch stochastically between active and inactive states. Enhancers can dictate transcriptional activity in animal development through the modulation of burst frequency, duration, or amplitude. Previous studies observed that different enhancers can achieve a wide range of transcriptional outputs through the same strategies of bursting control. For example, despite responding to different transcription factors, all even-skipped enhancers increase transcription by upregulating burst frequency and amplitude while burst duration remains largely constant. These shared bursting strategies suggest that a unified molecular mechanism constraints how enhancers modulate transcriptional output. Alternatively, different enhancers could have converged on the same bursting control strategy because of natural selection favoring one of these particular strategies. To distinguish between these two scenarios, we compared transcriptional bursting between endogenous and ectopic gene expression patterns. Because enhancers act under different regulatory inputs in ectopic patterns, dissimilar bursting control strategies between endogenous and ectopic patterns would suggest that enhancers adapted their bursting strategies to their trans-regulatory environment. Here, we generated ectopic even-skipped transcription patterns in fruit fly embryos and discovered that bursting strategies remain consistent in endogenous and ectopic even-skipped expression. These results provide evidence for a unified molecular mechanism shaping even-skipped bursting strategies and serve as a starting point to uncover the realm of strategies employed by other enhancers.

Comparative Single Cell Analysis of Transcriptional Bursting Reveals the Role of Genome Organization on de novo Transcript Origination

Spermatogenesis is a key developmental process underlying the origination of newly evolved genes. However, rapid cell type-specific transcriptomic divergence of the Drosophila germline has posed a significant technical barrier for comparative single-cell RNA-sequencing (scRNA-Seq) studies. By quantifying a surprisingly strong correlation between species-and cell type-specific divergence in three closely related Drosophila species, we apply a simple statistical procedure to identify a core set of 198 genes that are highly predictive of cell type identity while remaining robust to species-specific differences that span over 25-30 million years of evolution. We then utilize cell type classifications based on the 198-gene set to show how transcriptional divergence in cell type increases throughout spermatogenic developmental time, contrasting with traditional hourglass models of whole-organism development. With these cross-species cell type classifications, we then investigate the influence of genome organization on the molecular evolution of spermatogenesis vis-a-vis transcriptional bursting. We first demonstrate how mechanistic control of pre-meiotic transcription is achieved by altering transcriptional burst size while post-meiotic control is exerted via altered bursting frequency. We then report how global differences in autosomal vs. X chromosomal transcription likely arise in a developmental stage preceding full testis organogenesis by showing evolutionarily conserved decreases in X-linked transcription bursting kinetics in all examined somatic and germline cell types. Finally, we provide evidence supporting the cultivator model of de novo gene origination by demonstrating how the appearance of newly evolved testis-specific transcripts potentially provides short-range regulation of the transcriptional bursting properties of neighboring genes during key stages of spermatogenesis.

Real-time single-molecule imaging of transcriptional regulatory networks in living cells

Gene regulatory networks drive the specific transcriptional programmes responsible for the diversification of cell types during the development of multicellular organisms. Although our knowledge of the genes involved in these dynamic networks has expanded rapidly, our understanding of how transcription is spatiotemporally regulated at the molecular level over a wide range of timescales in the small volume of the nucleus remains limited. Over the past few decades, advances in the field of single-molecule fluorescence imaging have enabled real-time behaviours of individual transcriptional components to be measured in living cells and organisms. These efforts are now shedding light on the dynamic mechanisms of transcription, revealing not only the temporal rules but also the spatial coordination of underlying molecular interactions during various biological events.

T cell help induces Myc transcriptional bursts in germinal center B cells during positive selection

Antibody affinity maturation occurs in secondary lymphoid organs within germinal centers (GCs). At these sites, B cells mutate their antibody-encoding genes in the dark zone, followed by preferential selection of the high-affinity variants in the light zone by T cells. The strength of the T cell-derived selection signals is proportional to the B cell receptor affinity and to the magnitude of subsequent Myc expression. However, because the lifetime of Myc mRNA and its corresponding protein is very short, it remains unclear how T cells induce sustained Myc levels in positively selected B cells. Here, by direct visualization of mRNA and active transcription sites in situ, we found that an increase in transcriptional bursts promotes Myc expression during B cell positive selection in GCs. Elevated T cell help signals predominantly enhance the percentage of cells expressing Myc in GCs as opposed to augmenting the quantity of Myc transcripts per individual cell. Visualization of transcription start sites in situ revealed that T cell help promotes an increase in the frequency of transcriptional bursts at the Myc locus in GC B cells located primarily in the LZ apical rim. Thus, the rise in Myc, which governs positive selection of B cells in GCs, reflects an integration of transcriptional activity over time rather than an accumulation of transcripts at a specific time point.

Transcription factor exchange enables prolonged transcriptional bursts

Single-molecule imaging inside living cells has revealed that transcription factors (TFs) bind to DNA transiently, but a long-standing question is how this transient binding is related to transcription activation. Here, we devised a microscopy method to simultaneously measure transient TF binding at a single locus and the effect of these binding events on transcription. We show that DNA binding of the yeast TF Gal4 activates transcription of a target gene within a few seconds, with at least ∼20% efficiency and with a high initiation rate of ∼1 RNA/s. Gal4 DNA dissociation decreases transcription rapidly. Moreover, at a gene with multiple binding sites, individual Gal4 molecules only rarely stay bound throughout the entire burst but instead frequently exchange during a burst to increase the transcriptional burst duration. Our results suggest a mechanism for enhancer regulation in more complex eukaryotes, where TF cooperativity and exchange enable robust and responsive transcription regulation.

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.

Tissue-specific RNA Polymerase II promoter-proximal pause release and burst kinetics in a Drosophila embryonic patterning network

BACKGROUND

Formation of tissue-specific transcriptional programs underlies multicellular development, including dorsoventral (DV) patterning of the Drosophila embryo. This involves interactions between transcriptional enhancers and promoters in a chromatin context, but how the chromatin landscape influences transcription is not fully understood.

RESULTS

Here we comprehensively resolve differential transcriptional and chromatin states during Drosophila DV patterning. We find that RNA Polymerase II pausing is established at DV promoters prior to zygotic genome activation (ZGA), that pausing persists irrespective of cell fate, but that release into productive elongation is tightly regulated and accompanied by tissue-specific P-TEFb recruitment. DV enhancers acquire distinct tissue-specific chromatin states through CBP-mediated histone acetylation that predict the transcriptional output of target genes, whereas promoter states are more tissue-invariant. Transcriptome-wide inference of burst kinetics in different cell types revealed that while DV genes are generally characterized by a high burst size, either burst size or frequency can differ between tissues.

CONCLUSIONS

The data suggest that pausing is established by pioneer transcription factors prior to ZGA and that release from pausing is imparted by enhancer chromatin state to regulate bursting in a tissue-specific manner in the early embryo. Our results uncover how developmental patterning is orchestrated by tissue-specific bursts of transcription from Pol II primed promoters in response to enhancer regulatory cues.

Modulation of protein-DNA binding reveals mechanisms of spatiotemporal gene control in early Drosophila embryos

It is well known that enhancers regulate the spatiotemporal expression of their target genes by recruiting transcription factors (TFs) to the cognate binding sites in the region. However, the role of multiple binding sites for the same TFs and their specific spatial arrangement in determining the overall competency of the enhancer has yet to be fully understood. In this study, we utilized the MS2-MCP live imaging technique to quantitatively analyze the regulatory logic of the snail distal enhancer in early Drosophila embryos. Through systematic modulation of Dorsal and Twist binding motifs in this enhancer, we found that a mutation in any one of these binding sites causes a drastic reduction in transcriptional amplitude, resulting in a reduction in mRNA production of the target gene. We provide evidence of synergy, such that multiple binding sites with moderate affinities cooperatively recruit more TFs to drive stronger transcriptional activity than a single site. Moreover, a Hidden Markov-based stochastic model of transcription reveals that embryos with mutated binding sites have a higher probability of returning to the inactive promoter state. We propose that TF-DNA binding regulates spatial and temporal gene expression and drives robust pattern formation by modulating transcriptional kinetics and tuning bursting rates.

Regulatory Mechanisms for Transcriptional Bursting Revealed by an Event-Based Model

Gene transcription often occurs in discrete bursts, and it can be difficult to deduce the underlying regulatory mechanisms for transcriptional bursting with limited experimental data. Here, we categorize numerous states of single eukaryotic genes and identify 6 essential transcriptional events, each comprising a series of state transitions; transcriptional bursting is characterized as a sequence of 4 events, capable of being organized in various configurations, in addition to the beginning and ending events. By associating transcriptional kinetics with mean durations and recurrence probabilities of the events, we unravel how transcriptional bursting is modulated by various regulators including transcription factors. Through analytical derivation and numerical simulation, this study reveals key state transitions contributing to transcriptional sensitivity and specificity, typical characteristics of burst profiles, global constraints on intrinsic transcriptional noise, major regulatory modes in individual genes and across the genome, and requirements for fast gene induction upon stimulation. It is illustrated how biochemical reactions on different time scales are modulated to separately shape the durations and ordering of the events. Our results suggest that transcriptional patterns are essentially controlled by a shared set of transcriptional events occurring under specific promoter architectures and regulatory modes, the number of which is actually limited.

Lola-I is a promoter pioneer factor that establishes de novo Pol II pausing during development

While the accessibility of enhancers is dynamically regulated during development, promoters tend to be constitutively accessible and poised for activation by paused Pol II. By studying Lola-I, a Drosophila zinc finger transcription factor, we show here that the promoter state can also be subject to developmental regulation independently of gene activation. Lola-I is ubiquitously expressed at the end of embryogenesis and causes its target promoters to become accessible and acquire paused Pol II throughout the embryo. This promoter transition is required but not sufficient for tissue-specific target gene activation. Lola-I mediates this function by depleting promoter nucleosomes, similar to the action of pioneer factors at enhancers. These results uncover a level of regulation for promoters that is normally found at enhancers and reveal a mechanism for the de novo establishment of paused Pol II at promoters.

BurstDECONV: a signal deconvolution method to uncover mechanisms of transcriptional bursting in live cells

Monitoring transcription in living cells gives access to the dynamics of this complex fundamental process. It reveals that transcription is discontinuous, whereby active periods (bursts) are separated by one or several types of inactive periods of distinct lifetimes. However, decoding temporal fluctuations arising from live imaging and inferring the distinct transcriptional steps eliciting them is a challenge. We present BurstDECONV, a novel statistical inference method that deconvolves signal traces into individual transcription initiation events. We use the distribution of waiting times between successive polymerase initiation events to identify mechanistic features of transcription such as the number of rate-limiting steps and their kinetics. Comparison of our method to alternative methods emphasizes its advantages in terms of precision and flexibility. Unique features such as the direct determination of the number of promoter states and the simultaneous analysis of several potential transcription models make BurstDECONV an ideal analytic framework for live cell transcription imaging experiments. Using simulated realistic data, we found that our method is robust with regards to noise or suboptimal experimental designs. To show its generality, we applied it to different biological contexts such as Drosophila embryos or human cells.

Stochastic motion and transcriptional dynamics of pairs of distal DNA loci on a compacted chromosome

Chromosomes in the eukaryotic nucleus are highly compacted. However, for many functional processes, including transcription initiation, the pairwise motion of distal chromosomal elements such as enhancers and promoters is essential and necessitates dynamic fluidity. Here, we used a live-imaging assay to simultaneously measure the positions of pairs of enhancers and promoters and their transcriptional output while systematically varying the genomic separation between these two DNA loci. Our analysis reveals the coexistence of a compact globular organization and fast subdiffusive dynamics. These combined features cause an anomalous scaling of polymer relaxation times with genomic separation leading to long-ranged correlations. Thus, encounter times of DNA loci are much less dependent on genomic distance than predicted by existing polymer models, with potential consequences for eukaryotic gene expression.

Variability of the innate immune response is globally constrained by transcriptional bursting

Transcription of almost all mammalian genes occurs in stochastic bursts, however the fundamental control mechanisms that allow appropriate single-cell responses remain unresolved. Here we utilise single cell genomics data and stochastic models of transcription to perform global analysis of the toll-like receptor (TLR)-induced gene expression variability. Based on analysis of more than 2000 TLR-response genes across multiple experimental conditions we demonstrate that the single-cell, gene-by-gene expression variability can be empirically described by a linear function of the population mean. We show that response heterogeneity of individual genes can be characterised by the slope of the mean-variance line, which captures how cells respond to stimulus and provides insight into evolutionary differences between species. We further demonstrate that linear relationships theoretically determine the underlying transcriptional bursting kinetics, revealing different regulatory modes of TLR response heterogeneity. Stochastic modelling of temporal scRNA-seq count distributions demonstrates that increased response variability is associated with larger and more frequent transcriptional bursts, which emerge via increased complexity of transcriptional regulatory networks between genes and different species. Overall, we provide a methodology relying on inference of empirical mean-variance relationships from single cell data and new insights into control of innate immune response variability.

Bayesian analysis dissects kinetic modulation during non-stationary gene expression

Labelling of nascent stem loops with fluorescent proteins has fostered the visualization of transcription in living cells. Quantitative analysis of recorded fluorescence traces can shed light on kinetic transcription parameters and regulatory mechanisms. However, existing methods typically focus on steady state dynamics. Here, we combine a stochastic process transcription model with a hierarchical Bayesian method to infer global as well locally shared parameters for groups of cells and recover unobserved quantities such as initiation times and polymerase loading of the gene. We apply our approach to the cyclic response of the yeast CUP1 locus to heavy metal stress. Within the previously described slow cycle of transcriptional activity on the scale of minutes, we discover fast time-modulated bursting on the scale of seconds. Model comparison suggests that slow oscillations of transcriptional output are regulated by the amplitude of the bursts. Several polymerases may initiate during a burst.

Transcriptional bursting: stochasticity in deterministic development

The transcription of DNA by RNA polymerase occurs as a discontinuous process described as transcriptional bursting. This bursting behavior is observed across species and has been quantified using various stochastic modeling approaches. There is a large body of evidence that suggests the bursts are actively modulated by transcriptional machinery and play a role in regulating developmental processes. Under a commonly used two-state model of transcription, various enhancer-, promoter- and chromatin microenvironment-associated features are found to differentially influence the size and frequency of bursting events - key parameters of the two-state model. Advancement of modeling and analysis tools has revealed that the simple two-state model and associated parameters may not sufficiently characterize the complex relationship between these features. The majority of experimental and modeling findings support the view of bursting as an evolutionarily conserved transcriptional control feature rather than an unintended byproduct of the transcription process. Stochastic transcriptional patterns contribute to enhanced cellular fitness and execution of proper development programs, which posit this mode of transcription as an important feature in developmental gene regulation. In this Review, we present compelling examples of the role of transcriptional bursting in development and explore the question of how stochastic transcription leads to deterministic organism development.

Dynamic epistasis analysis reveals how chromatin remodeling regulates transcriptional bursting

Transcriptional bursting has been linked to the stochastic positioning of nucleosomes. However, how bursting is regulated by the remodeling of promoter nucleosomes is unknown. Here, we use single-molecule live-cell imaging of GAL10 transcription in Saccharomyces cerevisiae to measure how bursting changes upon combined perturbations of chromatin remodelers, the transcription factor Gal4 and preinitiation complex components. Using dynamic epistasis analysis, we reveal how the remodeling of different nucleosomes regulates transcriptional bursting parameters. At the nucleosome covering the Gal4 binding sites, RSC and Gal4 binding synergistically facilitate each burst. Conversely, nucleosome remodeling at the TATA box controls only the first burst upon galactose induction. At canonical TATA boxes, the nucleosomes are displaced by TBP binding to allow for transcription activation even in the absence of remodelers. Overall, our results reveal how promoter nucleosome remodeling together with Gal4 and preinitiation complex binding regulates transcriptional bursting.

Modulation of transcription burst amplitude underpins dosage compensation in the Drosophila embryo

Dosage compensation, the balancing of X-linked gene expression between sexes and to the autosomes, is critical to an organism's fitness and survival. In Drosophila, dosage compensation involves hypertranscription of the male X chromosome. Here, we use quantitative live imaging and modeling at single-cell resolution to study X chromosome dosage compensation in Drosophila. We show that the four X chromosome genes studied undergo transcriptional bursting in male and female embryos. Mechanistically, our data reveal that transcriptional upregulation of male X chromosome genes is primarily mediated by a higher RNA polymerase II initiation rate and burst amplitude across the expression domain. In contrast, burst frequency is spatially modulated in nuclei within the expression domain in response to different transcription factor concentrations to tune the transcriptional response. Together, these data show how the local and global regulation of distinct burst parameters can establish the complex transcriptional outputs underpinning developmental patterning.

Transcriptional activators in the early Drosophila embryo perform different kinetic roles

Combinatorial regulation of gene expression by transcription factors (TFs) may in part arise from kinetic synergy-wherein TFs regulate different steps in the transcription cycle. Kinetic synergy requires that TFs play distinguishable kinetic roles. Here, we used live imaging to determine the kinetic roles of three TFs that activate transcription in the Drosophila embryo-Zelda, Bicoid, and Stat92E-by introducing their binding sites into the even-skipped stripe 2 enhancer. These TFs influence different sets of kinetic parameters, and their influence can change over time. All three TFs increased the fraction of transcriptionally active nuclei; Zelda also shortened the first-passage time into transcription and regulated the interval between transcription events. Stat92E also increased the lifetimes of active transcription. Different TFs can therefore play distinct kinetic roles in activating the transcription. This has consequences for understanding the composition and flexibility of regulatory DNA sequences and the biochemical function of TFs. A record of this paper's transparent peer review process is included in the supplemental information.

Competing constraints shape the nonequilibrium limits of cellular decision-making

Gene regulation is central to cellular function. Yet, despite decades of work, we lack quantitative models that can predict how transcriptional control emerges from molecular interactions at the gene locus. Thermodynamic models of transcription, which assume that gene circuits operate at equilibrium, have previously been employed with considerable success in the context of bacterial systems. However, the presence of ATP-dependent processes within the eukaryotic transcriptional cycle suggests that equilibrium models may be insufficient to capture how eukaryotic gene circuits sense and respond to input transcription factor concentrations. Here, we employ simple kinetic models of transcription to investigate how energy dissipation within the transcriptional cycle impacts the rate at which genes transmit information and drive cellular decisions. We find that biologically plausible levels of energy input can lead to significant gains in how rapidly gene loci transmit information but discover that the regulatory mechanisms underlying these gains change depending on the level of interference from noncognate activator binding. When interference is low, information is maximized by harnessing energy to push the sensitivity of the transcriptional response to input transcription factors beyond its equilibrium limits. Conversely, when interference is high, conditions favor genes that harness energy to increase transcriptional specificity by proofreading activator identity. Our analysis further reveals that equilibrium gene regulatory mechanisms break down as transcriptional interference increases, suggesting that energy dissipation may be indispensable in systems where noncognate factor interference is sufficiently large.

Stochastic motion and transcriptional dynamics of pairs of distal DNA loci on a compacted chromosome

Chromosomes in the eukaryotic nucleus are highly compacted. However, for many functional processes, including transcription initiation, the 3D pair-wise motion of distal chromosomal elements, such as enhancers and promoters, is essential and necessitates dynamic fluidity. Therefore, the interplay of chromosome organization and dynamics is crucial for gene regulation. Here, we use a live imaging assay to simultaneously measure the positions of pairs of enhancers and promoters and their transcriptional output in the developing fly embryo while systematically varying the genomic separation between these two DNA loci. Our analysis reveals a combination of a compact globular organization and fast subdiffusive dynamics. These combined features cause an anomalous scaling of polymer relaxation times with genomic separation and lead to long-ranged correlations compared to existing polymer models. This scaling implies that encounter times of DNA loci are much less dependent on genomic separation than predicted by existing polymer models, with potentially significant consequences for eukaryotic gene expression.

Modulation of protein-DNA binding reveals mechanisms of spatiotemporal gene control in early Drosophila embryos

It is well known that enhancers regulate the spatiotemporal expression of their target genes by recruiting transcription factors (TFs) to the cognate binding sites in the region. However, the role of multiple binding sites for the same TFs and their specific spatial arrangement in determining the overall competency of the enhancer has yet to be fully understood. In this study, we utilized the MS2-MCP live imaging technique to quantitatively analyze the regulatory logic of the snail distal enhancer in early Drosophila embryos. Through systematic modulation of Dorsal and Twist binding motifs in this enhancer, we found that a mutation in any one of these binding sites causes a drastic reduction in transcriptional amplitude, resulting in a reduction in total mRNA production of the target gene. We provide evidence of synergy, such that multiple binding sites with moderate affinities cooperatively recruit more TFs to drive stronger transcriptional activity than a single site. Moreover, a Hidden Markov-based stochastic model of transcription reveals that embryos with mutated binding sites have a higher probability of returning to the inactive promoter state. We propose that TF-DNA binding regulates spatial and temporal gene expression and drives robust pattern formation by modulating transcriptional kinetics and tuning bursting rates.

Genome-wide inference reveals that feedback regulations constrain promoter-dependent transcriptional burst kinetics

Gene expression in mammalian cells is highly variable and episodic, resulting in a series of discontinuous bursts of mRNAs. A challenge is to understand how static promoter architecture and dynamic feedback regulations dictate bursting on a genome-wide scale. Although single-cell RNA sequencing (scRNA-seq) provides an opportunity to address this challenge, effective analytical methods are scarce. We developed an interpretable and scalable inference framework, which combined experimental data with a mechanistic model to infer transcriptional burst kinetics (sizes and frequencies) and feedback regulations. Applying this framework to scRNA-seq data generated from embryonic mouse fibroblast cells, we found Simpson's paradoxes, i.e. genome-wide burst kinetics exhibit different characteristics in two cases without and with distinguishing feedback regulations. We also showed that feedbacks differently modulate burst frequencies and sizes and conceal the effects of transcription start site distributions on burst kinetics. Notably, only in the presence of positive feedback, TATA genes are expressed with high burst frequencies and enhancer-promoter interactions mainly modulate burst frequencies. The developed inference method provided a flexible and efficient way to investigate transcriptional burst kinetics and the obtained results would be helpful for understanding cell development and fate decision.

Quantifying how post-transcriptional noise and gene copy number variation bias transcriptional parameter inference from mRNA distributions

Transcriptional rates are often estimated by fitting the distribution of mature mRNA numbers measured using smFISH (single molecule fluorescence in situ hybridization) with the distribution predicted by the telegraph model of gene expression, which defines two promoter states of activity and inactivity. However, fluctuations in mature mRNA numbers are strongly affected by processes downstream of transcription. In addition, the telegraph model assumes one gene copy but in experiments, cells may have two gene copies as cells replicate their genome during the cell cycle. While it is often presumed that post-transcriptional noise and gene copy number variation affect transcriptional parameter estimation, the size of the error introduced remains unclear. To address this issue, here we measure both mature and nascent mRNA distributions of GAL10 in yeast cells using smFISH and classify each cell according to its cell cycle phase. We infer transcriptional parameters from mature and nascent mRNA distributions, with and without accounting for cell cycle phase and compare the results to live-cell transcription measurements of the same gene. We find that: (i) correcting for cell cycle dynamics decreases the promoter switching rates and the initiation rate, and increases the fraction of time spent in the active state, as well as the burst size; (ii) additional correction for post-transcriptional noise leads to further increases in the burst size and to a large reduction in the errors in parameter estimation. Furthermore, we outline how to correctly adjust for measurement noise in smFISH due to uncertainty in transcription site localisation when introns cannot be labelled. Simulations with parameters estimated from nascent smFISH data, which is corrected for cell cycle phases and measurement noise, leads to autocorrelation functions that agree with those obtained from live-cell imaging.

Eukaryotic gene regulation at equilibrium, or non?

Models of transcriptional regulation that assume equilibrium binding of transcription factors have been less successful at predicting gene expression from sequence in eukaryotes than in bacteria. This could be due to the non-equilibrium nature of eukaryotic regulation. Unfortunately, the space of possible non-equilibrium mechanisms is vast and predominantly uninteresting. The key question is therefore how this space can be navigated efficiently, to focus on mechanisms and models that are biologically relevant. In this review, we advocate for the normative role of theory-theory that prescribes rather than just describes-in providing such a focus. Theory should expand its remit beyond inferring mechanistic models from data, towards identifying non-equilibrium gene regulatory schemes that may have been evolutionarily selected, despite their energy consumption, because they are precise, reliable, fast, or otherwise outperform regulation at equilibrium. We illustrate our reasoning by toy examples for which we provide simulation code.

Analysis of a non-Markov transcription model with nuclear RNA export and RNA nuclear retention

Transcription involves gene activation, nuclear RNA export (NRE) and RNA nuclear retention (RNR). All these processes are multistep and biochemical. A multistep reaction process can create memories between reaction events, leading to non-Markovian kinetics. This raises an unsolved issue: how does molecular memory affect stochastic transcription in the case that NRE and RNR are simultaneously considered? To address this issue, we analyze a non-Markov model, which considers multistep activation, multistep NRE and multistep RNR can interpret many experimental phenomena. In order to solve this model, we introduce an effective transition rate for each reaction. These effective transition rates, which explicitly decode the effect of molecular memory, can transform the original non-Markov issue into an equivalent Markov one. Based on this technique, we derive analytical results, showing that molecular memory can significantly affect the nuclear and cytoplasmic mRNA mean and noise. In addition to the results providing insights into the role of molecular memory in gene expression, our modeling and analysis provide a paradigm for studying more complex stochastic transcription processes.

Differential regulation of alternative promoters emerges from unified kinetics of enhancer-promoter interaction

Many eukaryotic genes contain alternative promoters with distinct expression patterns. How these promoters are differentially regulated remains elusive. Here, we apply single-molecule imaging to quantify the transcriptional regulation of two alternative promoters (P1 and P2) of the Bicoid (Bcd) target gene hunchback in syncytial blastoderm Drosophila embryos. Contrary to the previous notion that Bcd only activates P2, we find that Bcd activates both promoters via the same two enhancers. P1 activation is less frequent and requires binding of more Bcd molecules than P2 activation. Using a theoretical model to relate promoter activity to enhancer states, we show that the two promoters follow common transcription kinetics driven by sequential Bcd binding at the two enhancers. Bcd binding at either enhancer primarily activates P2, while P1 activation relies more on Bcd binding at both enhancers. These results provide a quantitative framework for understanding the kinetic mechanisms of complex eukaryotic gene regulation.

Alternative promoters differ in their expression patterns, whose mechanisms are not well understood. Here the authors show that alternative promoters of a Drosophila embryonic gene hunchback are regulated by different action modes of two enhancers.

The Impact of Space and Time on the Functional Output of the Genome

Over the past two decades, it has become clear that the multiscale spatial and temporal organization of the genome has important implications for nuclear function. This review centers on insights gained from recent advances in light microscopy on our understanding of transcription. We discuss spatial and temporal aspects that shape nuclear order and their consequences on regulatory components, focusing on genomic scales most relevant to function. The emerging picture is that spatiotemporal constraints increase the complexity in transcriptional regulation, highlighting new challenges, such as uncertainty about how information travels from molecular factors through the genome and space to generate a functional output.

Nup98-dependent transcriptional memory is established independently of transcription

Cellular ability to mount an enhanced transcriptional response upon repeated exposure to external cues is termed transcriptional memory, which can be maintained epigenetically through cell divisions and can depend on a nuclear pore component Nup98. The majority of mechanistic knowledge on transcriptional memory has been derived from bulk molecular assays. To gain additional perspective on the mechanism and contribution of Nup98 to memory, we used single-molecule RNA FISH (smFISH) to examine the dynamics of transcription in Drosophila cells upon repeated exposure to the steroid hormone ecdysone. We combined smFISH with mathematical modeling and found that upon hormone exposure, cells rapidly activate a low-level transcriptional response, but simultaneously begin a slow transition into a specialized memory state characterized by a high rate of expression. Strikingly, our modeling predicted that this transition between non-memory and memory states is independent of the transcription stemming from initial activation. We confirmed this prediction experimentally by showing that inhibiting transcription during initial ecdysone exposure did not interfere with memory establishment. Together, our findings reveal that Nup98's role in transcriptional memory is to stabilize the forward rate of conversion from low to high expressing state, and that induced genes engage in two separate behaviors - transcription itself and the establishment of epigenetically propagated transcriptional memory.

Scalable inference of transcriptional kinetic parameters from MS2 time series data

MOTIVATION

The MS2-MCP (MS2 coat protein) live imaging system allows for visualization of transcription dynamics through the introduction of hairpin stem-loop sequences into a gene. A fluorescent signal at the site of nascent transcription in the nucleus quantifies mRNA production. Computational modelling can be used to infer the promoter states along with the kinetic parameters governing transcription, such as promoter switching frequency and polymerase loading rate. However, modelling of the fluorescent trace presents a challenge due its persistence; the observed fluorescence at a given time point depends on both current and previous promoter states. A compound state Hidden Markov Model (cpHMM) was recently introduced to allow inference of promoter activity from MS2-MCP data. However, the computational time for inference scales exponentially with gene length and the cpHMM is therefore not currently practical for application to many eukaryotic genes.

RESULTS

We present a scalable implementation of the cpHMM for fast inference of promoter activity and transcriptional kinetic parameters. This new method can model genes of arbitrary length through the use of a time-adaptive truncated compound state space. The truncated state space provides a good approximation to the full state space by retaining the most likely set of states at each time during the forward pass of the algorithm. Testing on MS2-MCP fluorescent data collected from early Drosophila melanogaster embryos indicates that the method provides accurate inference of kinetic parameters within a computationally feasible timeframe. The inferred promoter traces generated by the model can also be used to infer single-cell transcriptional parameters.

AVAILABILITY AND IMPLEMENTATION

Python implementation is available at https://github.com/ManchesterBioinference/burstInfer, along with code to reproduce the examples presented here.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mean-field theory accurately captures the variation of copy number distributions across the mRNA life cycle

We consider a stochastic model where a gene switches between two states, an mRNA transcript is released in the active state, and subsequently it undergoes an arbitrary number of sequential unimolecular steps before being degraded. The reactions effectively describe various stages of the mRNA life cycle such as initiation, elongation, termination, splicing, export, and degradation. We construct a mean-field approach that leads to closed-form steady-state distributions for the number of transcript molecules at each stage of the mRNA life cycle. By comparison with stochastic simulations, we show that the approximation is highly accurate over all the parameter space, independent of the type of expression (constitutive or bursty) and of the shape of the distribution (unimodal, bimodal, and nearly bimodal). The theory predicts that in a population of identical cells, any bimodality is gradually washed away as the mRNA progresses through its life cycle.

Silent transcription intervals and translational bursting lead to diverse phenotypic switching

For complex process of gene expression, we use theoretical analysis and stochastic simulations to study the phenotypic diversity induced by silent transcription intervals and translational bursting.

Exact distributions for stochastic gene expression models with arbitrary promoter architecture and translational bursting

Gene expression in individual cells is inherently variable and sporadic, leading to cell-to-cell variability in mRNA and protein levels. Recent single-cell and single-molecule experiments indicate that promoter architecture and translational bursting play significant roles in controlling gene expression noise and generating the phenotypic diversity that life exhibits. To quantitatively understand the impact of these factors, it is essential to construct an accurate mathematical description of stochastic gene expression and find the exact analytical results, which is a formidable task. Here, we develop a stochastic model of bursty gene expression, which considers the complex promoter architecture governing the variability in mRNA expression and a general distribution characterizing translational burst. We derive the analytical expression for the corresponding protein steady-state distribution and all moment statistics of protein counts. We show that the total protein noise can be decomposed into three parts: the low-copy noise of protein due to probabilistic individual birth and death events, the noise due to stochastic switching between promoter states, and the noise resulting from translational busting. The theoretical results derived provide quantitative insights into the biochemical mechanisms of stochastic gene expression.

Two promoters integrate multiple enhancer inputs to drive wild-type knirps expression in the Drosophila melanogaster embryo

Proper development depends on precise spatiotemporal gene expression patterns. Most developmental genes are regulated by multiple enhancers and often by multiple core promoters that generate similar transcripts. We hypothesize that multiple promoters may be required either because enhancers prefer a specific promoter or because multiple promoters serve as a redundancy mechanism. To test these hypotheses, we studied the expression of the knirps locus in the early Drosophila melanogaster embryo, which is mediated by multiple enhancers and core promoters. We found that one of these promoters resembles a typical "sharp" developmental promoter, while the other resembles a "broad" promoter usually associated with housekeeping genes. Using synthetic reporter constructs, we found that some, but not all, enhancers in the locus show a preference for one promoter, indicating that promoters provide both redundancy and specificity. By analyzing the reporter dynamics, we identified specific burst properties during the transcription process, namely burst size and frequency, that are most strongly tuned by the combination of promoter and enhancer. Using locus-sized reporters, we discovered that enhancers with no promoter preference in a synthetic setting have a preference in the locus context. Our results suggest that the presence of multiple promoters in a locus is due both to enhancer preference and a need for redundancy and that "broad" promoters with dispersed transcription start sites are common among developmental genes. They also imply that it can be difficult to extrapolate expression measurements from synthetic reporters to the locus context, where other variables shape a gene's overall expression pattern.

Single-molecule conformational dynamics of a transcription factor reveals a continuum of binding modes controlling association and dissociation

Binding and unbinding of transcription factors to DNA are kinetically controlled to regulate the transcriptional outcome. Control of the release of the transcription factor NF-κB from DNA is achieved through accelerated dissociation by the inhibitor protein IκBα. Using single-molecule FRET, we observed a continuum of conformations of NF-κB in free and DNA-bound states interconverting on the subseconds to minutes timescale, comparable to in vivo binding on the seconds timescale, suggesting that structural dynamics directly control binding kinetics. Much of the DNA-bound NF-κB is partially bound, allowing IκBα invasion to facilitate DNA dissociation. IκBα induces a locked conformation where the DNA-binding domains of NF-κB are too far apart to bind DNA, whereas a loss-of-function IκBα mutant retains the NF-κB conformational ensemble. Overall, our results suggest a novel mechanism with a continuum of binding modes for controlling association and dissociation of transcription factors.

Cis-Regulatory Logic Produces Gene-Expression Noise Describing Phenotypic Heterogeneity in Bacteria

Gene transcriptional process is random. It occurs in bursts and follows single-molecular kinetics. Intermittent bursts are measured based on their frequency and size. They influence temporal fluctuations in the abundance of total mRNA and proteins by generating distinct transcriptional variations referred to as “noise”. Noisy expression induces uncertainty because the association between transcriptional variation and the extent of gene expression fluctuation is ambiguous. The promoter architecture and remote interference of different cis-regulatory elements are the crucial determinants of noise, which is reflected in phenotypic heterogeneity. An alternative perspective considers that cellular parameters dictating genome-wide transcriptional kinetics follow a universal pattern. Research on noise and systematic perturbations of promoter sequences reinforces that both gene-specific and genome-wide regulation occur across species ranging from bacteria and yeast to animal cells. Thus, deciphering gene-expression noise is essential across different genomics applications. Amidst the mounting conflict, it is imperative to reconsider the scope, progression, and rational construction of diversified viewpoints underlying the origin of the noise. Here, we have established an indication connecting noise, gene expression variations, and bacterial phenotypic variability. This review will enhance the understanding of gene-expression noise in various scientific contexts and applications.

Understanding the molecular mechanisms of transcriptional bursting

In recent years, it has been experimentally established that transcription, a fundamental biological process that involves the synthesis of messenger RNA molecules from DNA templates, does not proceed continuously as was expected. Rather, it exhibits a distinct dynamic behavior of alternating between productive phases when RNA molecules are actively synthesized and inactive phases when there is no RNA production at all. The bimodal transcriptional dynamics is now confirmed to be present in most living systems. This phenomenon is known as transcriptional bursting and it attracts significant amounts of attention from researchers in different fields. However, despite multiple experimental and theoretical investigations, the microscopic origin and biological functions of the transcriptional bursting remain unclear. Here we discuss the recent developments in uncovering the underlying molecular mechanisms of transcriptional bursting and our current understanding of them. Our analysis presents a physicochemical view of the processes that govern transcriptional bursting in living cells.

PRC1 drives Polycomb-mediated gene repression by controlling transcription initiation and burst frequency

The Polycomb repressive system plays a fundamental role in controlling gene expression during mammalian development. To achieve this, Polycomb repressive complexes 1 and 2 (PRC1 and PRC2) bind target genes and use histone modification-dependent feedback mechanisms to form Polycomb chromatin domains and repress transcription. The inter-relatedness of PRC1 and PRC2 activity at these sites has made it difficult to discover the specific components of Polycomb chromatin domains that drive gene repression and to understand mechanistically how this is achieved. Here, by exploiting rapid degron-based approaches and time-resolved genomics, we kinetically dissect Polycomb-mediated repression and discover that PRC1 functions independently of PRC2 to counteract RNA polymerase II binding and transcription initiation. Using single-cell gene expression analysis, we reveal that PRC1 acts uniformly within the cell population and that repression is achieved by controlling transcriptional burst frequency. These important new discoveries provide a mechanistic and conceptual framework for Polycomb-dependent transcriptional control.

Gene expression at a single-molecule level: implications for myelodysplastic syndromes and acute myeloid leukemia

Nongenetic 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 acute myeloid leukemia and myelodysplastic syndromes, which has historically been studied at the ensemble level, may require reevaluation. 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.

Approximating gene transcription dynamics using steady-state formulas

Understanding how genes in a single cell respond to dynamically changing signals has been a central question in stochastic gene transcription research. Recent studies have generated massive steady-state or snapshot mRNA distribution data of individual cells, and inferred a large spectrum of kinetic transcription parameters under varying conditions. However, there have been few algorithms to convert these static data into the temporal variation of kinetic rates. Real-time imaging has been developed to monitor stochastic transcription processes at the single-cell level, but the immense technicality has prevented its application to most endogenous loci in mammalian cells. In this article, we introduced a stochastic gene transcription model with variable kinetic rates induced by unstable cellular conditions. We approximated the transcription dynamics using easily obtained steady-state formulas in the model. We tested the approximation against experimental data in both prokaryotic and eukaryotic cells and further solidified the conditions that guarantee the robustness of the method. The method can be easily implemented to provide convenient tools for quantifying dynamic kinetics and mechanisms underlying the widespread static transcription data, and may shed a light on circumventing the limitation of current bursting data on transcriptional real-time imaging.

Quantitative imaging of transcription in living Drosophila embryos reveals the impact of core promoter motifs on promoter state dynamics

Genes are expressed in stochastic transcriptional bursts linked to alternating active and inactive promoter states. A major challenge in transcription is understanding how promoter composition dictates bursting, particularly in multicellular organisms. We investigate two key Drosophila developmental promoter motifs, the TATA box (TATA) and the Initiator (INR). Using live imaging in Drosophila embryos and new computational methods, we demonstrate that bursting occurs on multiple timescales ranging from seconds to minutes. TATA-containing promoters and INR-containing promoters exhibit distinct dynamics, with one or two separate rate-limiting steps respectively. A TATA box is associated with long active states, high rates of polymerase initiation, and short-lived, infrequent inactive states. In contrast, the INR motif leads to two inactive states, one of which relates to promoter-proximal polymerase pausing. Surprisingly, the model suggests pausing is not obligatory, but occurs stochastically for a subset of polymerases. Overall, our results provide a rationale for promoter switching during zygotic genome activation.

The early Drosophila embryo as a model system for quantitative biology

With the rise of new tools, from controlled genetic manipulations and optogenetics to improved microscopy, it is now possible to make clear, quantitative and reproducible measurements of biological processes. The humble fruit fly Drosophila melanogaster, with its ease of genetic manipulation combined with excellent imaging accessibility, has become a major model system for performing quantitative in vivo measurements. Such measurements are driving a new wave of interest from physicists and engineers, who are developing a range of testable dynamic models of active systems to understand fundamental biological processes. The reproducibility of the early Drosophila embryo has been crucial for understanding how biological systems are robust to unavoidable noise during development. Insights from quantitative in vivo experiments in the Drosophila embryo are having an impact on our understanding of critical biological processes, such as how cells make decisions and how complex tissue shape emerges. Here, to highlight the power of using Drosophila embryogenesis for quantitative biology, I focus on three main areas: (1) formation and robustness of morphogen gradients; (2) how gene regulatory networks ensure precise boundary formation; and (3) how mechanical interactions drive packing and tissue folding. I further discuss how such data has driven advances in modelling.

Exact results for queuing models of stochastic transcription with memory and crosstalk

Gene transcription is a complex multistep biochemical process, which can create memory between individual reaction events. On the other hand, many inducible genes, when activated by external cues, are often coregulated by several competitive pathways with crosstalk. This raises an unexplored question: how do molecular memory and crosstalk together affect gene expressions? To address this question, we introduce a queuing model of stochastic transcription, where two crossing signaling pathways are used to direct gene activation in response to external signals and memory functions to model multistep reaction processes involved in transcription. We first establish, based on the total probability principle, the chemical master equation for this queuing model, and then we derive, based on the binomial moment approach, exact expressions for statistical quantities (including distributions) of mRNA, which provide insights into the roles of crosstalk and memory in controlling the mRNA level and noise. We find that molecular memory of gene activation decreases the mRNA level but increases the mRNA noise, and double activation pathways always reduce the mRNA noise in contrast to a single pathway. In addition, we find that molecular memory can make the mRNA bimodality disappear.