Transcriptional Bursting Research Articles

Growth phase-dependent ribonucleic acid production dynamics

Transcriptional events play a crucial role in major cellular processes that specify the activity of an individual cells and influences cell population behavior in response to environment. Active (ON) and an inactive (OFF) states controls the transcriptional burst. Yet, the mechanism and kinetics of ON/OFF-state across the different growth phases of Escherichia coli remains elusive. Here, we have used a single mRNA detection method in live-cells to comprehend the ON/OFF mechanism of the first transcriptional (TF) and consecutive events (TC) controlled by lactose promoters, Plac and Plac/ara1. We determined that the duration of TF ON/OFF has different modes, exhibiting a close to inverse behavior to that of TC ON/OFF. Dynamics of ON/OFF states in fast and slow-dividing cells were affected by the promoter region during the initiation of transcription. Period of TF ON-state defines the behavior of TC by altering the number and the frequency of mRNAs formed. Furthermore, we have shown that delayed OFF-time in TF affects the dynamics of TC in both states, which is mainly determined by the upstream promoter region. Furthermore, using elongation arrest experiments, we independently validate that mRNA noise in TC is governed by the delayed OFF-period in TF. We have identified the position of the regulatory regions that plays a crucial role in noise (Fano) modulation. Taken together, our results suggest that the dynamics of the first transcriptional event, TF, pre-defines the diversity of the population.

Spectral neural approximations for models of transcriptional dynamics

The advent of high-throughput transcriptomics provides an opportunity to advance mechanistic understanding of transcriptional processes and their connections to cellular function at an unprecedented, genome-wide scale. These transcriptional systems, which involve discrete stochastic events, are naturally modeled using chemical master equations (CMEs), which can be solved for probability distributions to fit biophysical rates that govern system dynamics. While CME models have been used as standards in fluorescence transcriptomics for decades to analyze single-species RNA distributions, there are often no closed-form solutions to CMEs that model multiple species, such as nascent and mature RNA transcript counts. This has prevented the application of standard likelihood-based statistical methods for analyzing high-throughput, multi-species transcriptomic datasets using biophysical models. Inspired by recent work in machine learning to learn solutions to complex dynamical systems, we leverage neural networks and statistical understanding of system distributions to produce accurate approximations to a steady-state bivariate distribution for a model of the RNA life cycle that includes nascent and mature molecules. The steady-state distribution to this simple model has no closed-form solution and requires intensive numerical solving techniques: our approach reduces likelihood evaluation time by several orders of magnitude. We demonstrate two approaches, whereby solutions are approximated by 1) learning the weights of kernel distributions with constrained parameters or 2) learning both weights and scaling factors for parameters of kernel distributions. We show that our strategies, denoted by kernel weight regression and parameter-scaled kernel weight regression, respectively, enable broad exploration of parameter space and can be used in existing likelihood frameworks to infer transcriptional burst sizes, RNA splicing rates, and mRNA degradation rates from experimental transcriptomic data.

Dissection and integration of bursty transcriptional dynamics for complex systems

RNA velocity estimation is a potentially powerful tool to reveal the directionality of transcriptional changes in single-cell RNA-sequencing data, but it lacks accuracy, absent advanced metabolic labeling techniques. We developed an approach, TopicVelo, that disentangles simultaneous, yet distinct, dynamics by using a probabilistic topic model, a highly interpretable form of latent space factorization, to infer cells and genes associated with individual processes, thereby capturing cellular pluripotency or multifaceted functionality. Focusing on process-associated cells and genes enables accurate estimation of process-specific velocities via a master equation for a transcriptional burst model accounting for intrinsic stochasticity. The method obtains a global transition matrix by leveraging cell topic weights to integrate process-specific signals. In challenging systems, this method accurately recovers complex transitions and terminal states, while our use of first-passage time analysis provides insights into transient transitions. These results expand the limits of RNA velocity, empowering future studies of cell fate and functional responses.

A critical role for Pol II CTD phosphorylation in heterochromatic gene activation

How gene activation works in heterochromatin, and how the mechanism might differ from the one used in euchromatin, has been largely unexplored. Previous work has shown that in SIR-regulated heterochromatin of Saccharomyces cerevisiae, gene activation occurs in the absence of covalent histone modifications and other alterations of chromatin commonly associated with transcription.Here we demonstrate that such activation occurs in a substantial fraction of cells, consistent with frequent transcriptional bursting, and this raises the possibility that an alternative activation pathway might be used. We address one such possibility, Pol II CTD phosphorylation, and explore this idea using a natural telomere-linked gene, YFR057w, as a model. Unlike covalent histone modifications, we find that Ser2, Ser5 and Ser7 CTD phosphorylated Pol II is prevalent at the drug-induced heterochromatic gene. Particularly enriched relative to the euchromatic state is Ser2 phosphorylation. Consistent with a functional role for Ser2P, YFR057w is negligibly activated in cells deficient in the Ser2 CTD kinases Ctk1 and Bur1 even though the gene is strongly stimulated when it is placed in a euchromatic context. Collectively, our results are consistent with a critical role for Ser2 CTD phosphorylation in driving Pol II recruitment and transcription of a natural heterochromatic gene – an activity that may supplant the need for histone epigenetic modifications.

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.

Disordered C-terminal domain drives spatiotemporal confinement of RNAPII to enhance search for chromatin targets.

Efficient gene expression requires RNA polymerase II (RNAPII) to find chromatin targets precisely in space and time. How RNAPII manages this complex diffusive search in three-dimensional nuclear space remains largely unknown. The disordered carboxy-terminal domain (CTD) of RNAPII, which is essential for recruiting transcription-associated proteins, forms phase-separated droplets in vitro, hinting at a potential role in modulating RNAPII dynamics. In the present study, we use single-molecule tracking and spatiotemporal mapping in living yeast to show that the CTD is required for confining RNAPII diffusion within a subnuclear region enriched for active genes, but without apparent phase separation into condensates. Both Mediator and global chromatin organization are required for sustaining RNAPII confinement. Remarkably, truncating the CTD disrupts RNAPII spatial confinement, prolongs target search, diminishes chromatin binding, impairs pre-initiation complex formation and reduces transcription bursting. The present study illuminates the pivotal role of the CTD in driving spatiotemporal confinement of RNAPII for efficient gene expression.

Yeast heterochromatin stably silences only weak regulatory elements by altering burst duration

Transcriptional silencing in Saccharomyces cerevisiae involves the generation of a chromatin state that stably represses transcription. Using multiple reporter assays, a diverse set of upstream activating sequence enhancers and core promoters were investigated for their susceptibility to silencing. We show that heterochromatin stably silences only weak and stress-induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements, and the partial repression of these elements did not result in bistable expression states. Permutation analysis of enhancers and promoters indicates that both elements are targets of repression. Chromatin remodelers help specific regulatory elements to resist repression, most probably by altering nucleosome mobility and changing transcription burst duration. The strong enhancers/promoters can be repressed if silencer-bound Sir1 is increased. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating-type gene regulatory elements but not strong housekeeping gene regulatory sequences.

Molecular-Memory-Induced Counter-Intuitive Noise Attenuator in Protein Polymerization

Gene expression comprises many asymmetric and complex processes. Transcriptional details revealed by the whole genome indicate that genes resort to transcriptional bursting and accumulate molecular memory. However, it is still unclear how the interplay of transcriptional bursting and memory regulates robustness and expression noise. Here, we consider a model of multiple coupled processes of protein polymerization to focus on decoding the effect of molecular memory. Using non-Markovian transformation technology, we first define the memory index to measure the correlation window of expression to decipher the mechanism of regulation. The results indicate that memory from synthesis can amplify expression noise, while memory originating from polymerization can reduce the lower bound of the noise of gene products; that is, the memory from different sources plays distinct regulatory roles to induce non-symmetry. Moreover, it is counterintuitive that the dual regulation from memory and bursting expression can directly suppress system noise, violating the principle that transcriptional bursting enhances noise. Our results not only provide a theoretical framework for investigating the function of memory but also imply that expression noise is not part of a half-power relationship with, nor mediated by, memory.

How transcription factor binding controls transcriptional bursting dynamics: A single-molecule view

In this issue, Pomp etal.1 simultaneously tracked transcriptional bursts of yeast gene GAL10 and transient binding of transcription factor Gal4 at the gene using novel methods. Dynamic exchange and infrequent long binding of Gal4 together enable prolonged transcriptional bursts of GAL10.

Implication of polymerase recycling for nascent transcript quantification by live cell imaging.

Transcription enables the production of RNA from a DNA template. Due to the highly dynamic nature of transcription, live-cell imaging methods play a crucial role in measuring the kinetics of this process. For instance, transcriptional bursts have been visualized using fluorescent phage-coat proteins that associate tightly with messenger RNA (mRNA) stem loops formed on nascent transcripts. To convert the signal emanating from a transcription site into meaningful estimates of transcription dynamics, the influence of various parameters on the measured signal must be evaluated. Here, the effect of gene length on the intensity of the transcription site focus was analyzed. Intuitively, a longer gene can support a larger number of transcribing polymerases, thus leading to an increase in the measured signal. However, measurements of transcription induced by hyper-osmotic stress responsive promoters display independence from gene length. A mathematical model of the stress-induced transcription process suggests that the formation of gene loops that favor the recycling of polymerase from the terminator to the promoter can explain the observed behavior. One experimentally validated prediction from this model is that the amount of mRNA produced from a short gene should be higher than for a long one as the density of active polymerase on the short gene will be increased by polymerase recycling. Our data suggest that this recycling contributes significantly to the expression output from a gene and that polymerase recycling is modulated by the promoter identity and the cellular state.

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.

Epigenomic states contribute to coordinated allelic transcriptional bursting in iPSC reprogramming.

Two alleles of a gene can be transcribed independently or coordinatedly, which can lead to temporal expression heterogeneity with potentially distinct impacts on cell fate. Here, we profiled genome-wide allelic transcriptional burst kinetics during the reprogramming of MEF to induced pluripotent stem cells. We show that the degree of coordination of allelic bursting differs among genes, and alleles of many reprogramming-related genes burst in a highly coordinated fashion. Notably, we show that the chromatin accessibility of the two alleles of highly coordinated genes is similar, unlike the semi-coordinated or independent genes, suggesting the degree of coordination of allelic bursting is linked to allelic chromatin accessibility. Consistently, we show that many transcription factors have differential binding affinity between alleles of semi-coordinated or independent genes. We show that highly coordinated genes are enriched with chromatin accessibility regulators such as H3K4me3, H3K4me1, H3K36me3, H3K27ac, histone variant H3.3, and BRD4. Finally, we demonstrate that enhancer elements are highly enriched in highly coordinated genes. Our study demonstrates that epigenomic states contribute to coordinated allelic bursting to fine-tune gene expression during induced pluripotent stem cell reprogramming.

The Mediator kinase module enhances polymerase activity to regulate transcriptional memory after heat stress in Arabidopsis

Plants are often exposed to recurring adverse environmental conditions in the wild. Acclimation to high temperatures entails transcriptional responses, which prime plants to better withstand subsequent stress events. Heat stress (HS)-induced transcriptional memory results in more efficient re-induction of transcription upon recurrence of heat stress. Here, we identified CDK8 and MED12, two subunits of the kinase module of the transcription co-regulator complex, Mediator, as promoters of heat stress memory and associated histone modifications in Arabidopsis. CDK8 is recruited to heat-stress memory genes by HEAT SHOCK TRANSCRIPTION FACTOR A2 (HSFA2). Like HSFA2, CDK8 is largely dispensable for the initial gene induction upon HS, and its function in transcriptional memory is thus independent of primary gene activation. In addition to the promoter and transcriptional start region of target genes, CDK8 also binds their 3ʹ-region, where it may promote elongation, termination, or rapid re-initiation of RNA polymerase II (Pol II) complexes during transcriptional memory bursts. Our work presents a complex role for the Mediator kinase module during transcriptional memory in multicellular eukaryotes, through interactions with transcription factors, chromatin modifications, and promotion of Pol II efficiency.

Brooklyn plots to identify co-expression dysregulation in single cell sequencing.

Altered open chromatin regions, impacting gene expression, is a feature of some human disorders. We discovered it is possible to detect global changes in genomically-related adjacent gene co-expression within single cell RNA sequencing (scRNA-seq) data. We built a software package to generate and test non-randomness using 'Brooklyn plots' to identify the percent of genes significantly co-expressed from the same chromosome in ∼10MB intervals across the genome. These plots establish an expected low baseline of co-expression in scRNA-seq from most cell types, but, as seen in dilated cardiomyopathy cardiomyocytes, altered patterns of open chromatin appear. These may relate to larger regions of transcriptional bursting, observable in single cell, but not bulk datasets.

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

BackgroundFormation 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.ResultsHere 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.ConclusionsThe 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.

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.

Power-law behavior of transcriptional bursting regulated by enhancer-promoter communication.

Revealing how transcriptional bursting kinetics are genomically encoded is challenging because genome structures are stochastic at the organization level and are suggestively linked to gene transcription. To address this challenge, we develop a generic theoretical framework that integrates chromatin dynamics, enhancer-promoter (E-P) communication, and gene-state switching to study transcriptional bursting. The theory predicts that power law can be a general rule to quantitatively describe bursting modulations by E-P spatial communication. Specifically, burst frequency and burst size are up-regulated by E-P communication strength, following power laws with positive exponents. Analysis of the scaling exponents further reveals that burst frequency is preferentially regulated. Bursting kinetics are down-regulated by E-P genomic distance with negative power-law exponents, and this negative modulation desensitizes at large distances. The mutual information between burst frequency (or burst size) and E-P spatial distance further reveals essential characteristics of the information transfer from E-P communication to transcriptional bursting kinetics. These findings, which are in agreement with experimental observations, not only reveal fundamental principles of E-P communication in transcriptional bursting but also are essential for understanding cellular decision-making.

Enhancer-promoter communication in Drosophila developmental gene transcription.

Enhancers play an essential role in gene regulation by receiving cues from transcription factors and relaying these signals to modulate transcription from target promoters. Enhancer-promoter communications occur across large linear distances of the genome and with high specificity. The molecular mechanisms that underlie enhancer-mediated control of transcription remain unresolved. In this review, we focus on research in Drosophila uncovering the molecular mechanisms governing enhancer-promoter communication and discuss the current understanding of developmental gene regulation. The functions of protein acetylation, pausing of RNA polymerase II, transcriptional bursting, and the formation of nuclear hubs in the induction of tissue-specific programs of transcription during zygotic genome activation are considered.

Modeling Control of Supercoiling Dynamics and Transcription Using DNA-Binding Proteins.

Nearly all natural and synthetic gene networks rely on the fundamental process of transcription to enact biological feedback, genetic programs, and living circuitry. In this work, we investigate the efficacy of controlling transcription using a new biophysical mechanism, control of localized supercoiling near a gene of interest. We postulate a basic reaction network model for describing the general phenomenon of transcription and introduce a separate set of equations to describe the dynamics of supercoiling. We show that supercoiling and transcription introduce a shared reaction flux term in the model dynamics and illustrate how the modulation of supercoiling can be used to control transcription rates. We show the supercoiling-transcription model can be written as a nonlinear state-space model, with a radial basis function nonlinearity to capture the empirical relationship between supercoiling and transcription rates. We show the system admits a single, globally exponentially stable equilibrium point. Notably, we show that mRNA steady-state levels can be controlled directly by increasing a length-scale parameter for genetic spacing. Finally, we build a mathematical model to explore the use of a DNA binding protein, to define programmable boundary conditions on supercoiling propagation, which we show can be used to control transcriptional bursting or pulsatile transcriptional response. We show there exists a stabilizing control law for mRNA tracking, using the method of control Lyapunov functions and illustrate these results with numerical simulations.

Phase separation reduces cell-to-cell variability of transcriptional bursting

Gene expression is a stochastic and noisy process often occurring in "bursts". Experiments have shown that the compartmentalization of proteins by liquid-liquid phase separation is conducive to reducing the noise of gene expression. Therefore, an important goal is to explore the role of bursts in phase separation noise reduction processes. We propose a coupled model that includes phase separation and a two-state gene expression process. Using the timescale separation method, we obtain approximate solutions for the expectation, variance, and noise strength of the dilute phase. We find that a higher burst frequency weakens the ability of noise reduction by phase separation, but as the burst size increases, this ability first increases and then decreases. This study provides a deeper understanding of phase separation to reduce noise in the stochastic gene expression with burst kinetics.