Multimodal transcriptional control of pattern formation in embryonic development

MONITTR allows real-time imaging of transcription and endogenous proteins in C. elegans

Maximizing cell survival under stress requires rapid and transient adjustments of RNA and protein synthesis. However, capturing these dynamic changes at both single-cell level and across an organism has been challenging. Here, we developed a system named MONITTR (MS2-embedded mCherry-based monitoring of transcription) for real-time simultaneous measurement of nascent transcripts and endogenous protein levels in C. elegans. Utilizing this system, we monitored the transcriptional bursting of fasting-induced genes and found that the epidermis responds to fasting by modulating the proportion of actively transcribing nuclei and transcriptional kinetics of individual alleles. Additionally, our findings revealed the essential roles of the transcription factors NHR-49 and HLH-30 in governing the transcriptional kinetics of fasting-induced genes under fasting. Furthermore, we tracked transcriptional dynamics during heat-shock response and ER unfolded protein response and observed rapid changes in the level of nascent transcripts under stress conditions. Collectively, our study provides a foundation for quantitatively investigating how animals spatiotemporally modulate transcription in various physiological and pathological conditions.

A simple MiMIC-based approach for tagging endogenous genes to visualise live transcription in Drosophila

Live imaging of transcription in the Drosophila embryo using the MS2 or PP7 systems is transforming our understanding of transcriptional regulation. However, insertion of MS2/PP7 stem-loops into endogenous genes requires laborious CRISPR genome editing. Here, we exploit the previously described Minos-mediated integration cassette (MiMIC) transposon system in Drosophila to establish a method for simply and rapidly inserting MS2/PP7 cassettes into any of the thousands of genes carrying a MiMIC insertion. In addition to generating a variety of stem-loop donor fly stocks, we have made new stocks expressing the complementary coat proteins fused to different fluorescent proteins. We show the utility of this MiMIC-based approach by MS2/PP7 tagging of endogenous genes and the long non-coding RNA roX1, then imaging their transcription in living embryos. We also present live transcription data from larval brains, the wing disc and ovary, thereby extending the tissues that can be studied using the MS2/PP7 system. Overall, this first high-throughput method for tagging mRNAs in Drosophila will facilitate the study of transcription dynamics of thousands of endogenous genes in a range of Drosophila tissues.

Bicoid-nucleosome competition sets a concentration threshold for transcription constrained by genome replication

Transcription factors (TFs) regulate gene expression despite constraints from chromatin structure and the cell cycle. Here we examine the concentration-dependent regulation of hunchback by the Bicoid morphogen through a combination of quantitative imaging, mathematical modeling and epigenomics in Drosophila embryos. By live imaging of MS2 reporters, we find that, following mitosis, the timing of transcriptional activation driven by the hunchback P2 (hb P2) enhancer directly reflects Bicoid concentration. We build a stochastic model that can explain in vivo onset time distributions by accounting for both the competition between Bicoid and nucleosomes at hb P2 and a negative influence of DNA replication on transcriptional elongation. Experimental modulation of nucleosome stability alters onset time distributions and the posterior boundary of hunchback expression. We conclude that TF-nucleosome competition is the molecular mechanism whereby the Bicoid morphogen gradient specifies the posterior boundary of hunchback expression.

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, in Berrocal et al., 2020, we showed that 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.

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.

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.

Dynamic microenvironments shape nuclear organization and gene expression

Live imaging has revealed that the regulation of gene expression is largely driven by transient interactions. For example, many regulatory proteins bind chromatin for just seconds, and loop-like genomic contacts are rare and last only minutes. These discoveries have been difficult to reconcile with our canonical models that are predicated on stable and hierarchical interactions. Proteomic microenvironments that concentrate nuclear factors may explain how brief interactions can still mediate gene regulation by creating conditions where reactions occur more frequently. Here, we summarize new imaging technologies and recent discoveries implicating microenvironments as a potential driver of nuclear function. Finally, we propose that key properties of proteomic microenvironments, such as their size, enrichment, and lifetimes, are directly linked to regulatory function.

Single-cell transcriptional dynamics in a living vertebrate

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

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.

4D nucleome equation predicts gene expression controlled by long-range enhancer-promoter interaction

Recent experimental evidence strongly supports that three-dimensional (3D) long-range enhancer-promoter (E-P) interactions have important influences on gene-expression dynamics, but it is unclear how the interaction information is translated into gene expression over time (4D). To address this question, we developed a general theoretical framework (named as a 4D nucleome equation), which integrates E-P interactions on chromatin and biochemical reactions of gene transcription. With this equation, we first present the distribution of mRNA counts as a function of the E-P genomic distance and then reveal a power-law scaling of the expression level in this distance. Interestingly, we find that long-range E-P interactions can induce bimodal and trimodal mRNA distributions. The 4D nucleome equation also allows for model selection and parameter inference. When this equation is applied to the mouse embryonic stem cell smRNA-FISH data and the E-P genomic-distance data, the predicted E-P contact probability and mRNA distribution are in good agreement with experimental results. Further statistical inference indicates that the E-P interactions prefer to modulate the mRNA level by controlling promoter activation and transcription initiation rates. Our model and results provide quantitative insights into both spatiotemporal gene-expression determinants (i.e., long-range E-P interactions) and cellular fates during development.

GRP78 acts as a cAMP/PKA signaling modulator through the MC4R pathway in porcine embryonic development

Glucose-regulated protein 78 (GRP78) binds to and stabilizes melanocortin 4 receptor (MC4R), which activates protein kinase A (PKA) by regulating G proteins. GRP78 is primarily used as a marker for endoplasmic reticulum stress; however, its other functions have not been well studied. Therefore, in this study, we aimed to investigate the function of GRP78 during porcine embryonic development. The developmental quality of porcine embryos, expression of cell cycle proteins, and function of mitochondria were evaluated by inhibiting the function of GRP78. Porcine oocytes were activated to undergo parthenogenesis, and blastocysts were obtained after 7 days of in vitro culture. GRP78 function was inhibited by adding 20 μM HA15 to the in vitro culture medium. The inhibition in GRP78 function led to a decrease in G proteins release, which subsequently downregulated the cyclic adenosine monophosphate (cAMP)/PKA pathway. Ultimately, inhibition of GRP78 function induced the inhibition of CDK1 and cyclin B expression and disruption of the cell cycle. In addition, inhibition of GRP78 function regulated DRP1 and SIRT1 expression, resulting in mitochondrial dysfunction. This study provides new insights into the role of GRP78 in porcine embryonic development, particularly its involvement in the regulation of the MC4R pathway and downstream cAMP/PKA signaling. The results suggest that the inhibition of GRP78 function in porcine embryos by HA15 treatment may have negative effects on embryo quality and development. This study also demonstrated that GRP78 plays a crucial role in the functioning of MC4R, which releases the G protein during porcine embryonic development.

A circadian-like gene network programs the timing and dosage of heterochronic miRNA transcription during C. elegans development

Development relies on the exquisite control of both the timing and the levels of gene expression to achieve robust developmental transitions. How cis- and trans-acting factors control both aspects simultaneously is unclear. We show that transcriptional pulses of the temporal patterning microRNA (miRNA) lin-4 are generated by two nuclear hormone receptors (NHRs) in C. elegans, NHR-85 and NHR-23, whose mammalian orthologs, Rev-Erb and ROR, function in the circadian clock. Although Rev-Erb and ROR antagonize each other to control once-daily transcription in mammals, NHR-85/NHR-23 heterodimers bind cooperatively to lin-4 regulatory elements to induce a single pulse of expression during each larval stage. Each pulse's timing, amplitude, and duration are dictated by the phased expression of these NHRs and the C. elegans Period ortholog, LIN-42, that binds to and represses NHR-85. Therefore, during nematode temporal patterning, an evolutionary rewiring of circadian clock components couples the timing of gene expression to the control of transcriptional dosage.

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.

Optogenetic control of YAP reveals a dynamic communication code for stem cell fate and proliferation

YAP is a transcriptional regulator that controls pluripotency, cell fate, and proliferation. How cells ensure the selective activation of YAP effector genes is unknown. This knowledge is essential to rationally control cellular decision-making. Here we leverage optogenetics, live-imaging of transcription, and cell fate analysis to understand and control gene activation and cell behavior. We reveal that cells decode the steady-state concentrations and timing of YAP activation to control proliferation, cell fate, and expression of the pluripotency regulators Oct4 and Nanog. While oscillatory YAP inputs induce Oct4 expression and proliferation optimally at frequencies that mimic native dynamics, cellular differentiation requires persistently low YAP levels. We identify the molecular logic of the Oct4 dynamic decoder, which acts through an adaptive change sensor. Our work reveals how YAP levels and dynamics enable multiplexing of information transmission for the regulation of developmental decision-making and establishes a platform for the rational control of these behaviors.

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.

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.

Dynamics of an incoherent feedforward loop drive ERK-dependent pattern formation in the early Drosophila embryo

Positional information in development often manifests as stripes of gene expression, but how stripes form remains incompletely understood. Here, we use optogenetics and live-cell biosensors to investigate the posterior brachyenteron (byn) stripe in early Drosophila embryos. This stripe depends on interpretation of an upstream ERK activity gradient and the expression of two target genes, tailless (tll) and huckebein (hkb), that exert antagonistic control over byn. We find that high or low doses of ERK signaling produce transient or sustained byn expression, respectively. Although tll transcription is always rapidly induced, hkb converts graded ERK inputs into a variable time delay. Nuclei thus interpret ERK amplitude through the relative timing of tll and hkb transcription. Antagonistic regulatory paths acting on different timescales are hallmarks of an incoherent feedforward loop, which is sufficient to explain byn dynamics and adds temporal complexity to the steady-state model of byn stripe formation. We further show that 'blurring' of an all-or-none stimulus through intracellular diffusion non-locally produces a byn stripe. Overall, we provide a blueprint for using optogenetics to dissect developmental signal interpretation in space and time.

A bistable autoregulatory module in the developing embryo commits cells to binary expression fates

Bistable autoactivation has been proposed as a mechanism for cells to adopt binary fates during embryonic development. However, it is unclear whether the autoactivating modules found within developmental gene regulatory networks are bistable, unless their parameters are quantitatively determined. Here, we combine in vivo live imaging with mathematical modeling to dissect the binary cell fate dynamics of the fruit fly pair-rule gene fushi tarazu (ftz), which is regulated by two known enhancers: the early (non-autoregulating) element and the autoregulatory element. Live imaging of transcription and protein concentration in the blastoderm revealed that binary Ftz fates are achieved as Ftz expression rapidly transitions from being dictated by the early element to the autoregulatory element. Moreover, we discovered that Ftz concentration alone is insufficient to activate the autoregulatory element, and that this element only becomes responsive to Ftz at a prescribed developmental time. Based on these observations, we developed a dynamical systems model and quantitated its kinetic parameters directly from experimental measurements. Our model demonstrated that the ftz autoregulatory module is indeed bistable and that the early element transiently establishes the content of the binary cell fate decision to which the autoregulatory module then commits. Further in silico analysis revealed that the autoregulatory element locks the Ftz fate quickly, within 35 min of exposure to the transient signal of the early element. Overall, our work confirms the widely held hypothesis that autoregulation can establish developmental fates through bistability and, most importantly, provides a framework for the quantitative dissection of cellular decision-making.

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.

Minimal synthetic enhancers reveal control of the probability of transcriptional engagement and its timing by a morphogen gradient

How enhancers interpret morphogen gradients to generate gene expression patterns is a central question in developmental biology. Recent studies have proposed that enhancers can dictate whether, when, and at what rate promoters engage in transcription, but the complexity of endogenous enhancers calls for theoretical models with too many free parameters to quantitatively dissect these regulatory strategies. To overcome this limitation, we established a minimal promoter-proximal synthetic enhancer in embryos of Drosophila melanogaster. Here, a gradient of the Dorsal activator is read by a single Dorsal DNA binding site. Using live imaging to quantify transcriptional activity, we found that a single binding site can regulate whether promoters engage in transcription in a concentration-dependent manner. By modulating the binding-site affinity, we determined that a gene's decision to transcribe and its transcriptional onset time can be explained by a simple model where the promoter traverses multiple kinetic barriers before transcription can ensue.

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.

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.

Predictive modeling reveals that higher-order cooperativity drives transcriptional repression in a synthetic developmental enhancer

A challenge in quantitative biology is to predict output patterns of gene expression from knowledge of input transcription factor patterns and from the arrangement of binding sites for these transcription factors on regulatory DNA. We tested whether widespread thermodynamic models could be used to infer parameters describing simple regulatory architectures that inform parameter-free predictions of more complex enhancers in the context of transcriptional repression by Runt in the early fruit fly embryo. By modulating the number and placement of Runt binding sites within an enhancer, and quantifying the resulting transcriptional activity using live imaging, we discovered that thermodynamic models call for higher-order cooperativity between multiple molecular players. This higher-order cooperativity captures the combinatorial complexity underlying eukaryotic transcriptional regulation and cannot be determined from simpler regulatory architectures, highlighting the challenges in reaching a predictive understanding of transcriptional regulation in eukaryotes and calling for approaches that quantitatively dissect their molecular nature.

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.

Dynamics of Drosophila endoderm specification

During early Drosophila embryogenesis, a network of gene regulatory interactions orchestrates terminal patterning, playing a critical role in the subsequent formation of the gut. We utilized CRISPR gene editing at endogenous loci to create live reporters of transcription and light-sheet microscopy to monitor the individual components of the posterior gut patterning network across 90 min prior to gastrulation. We developed a computational approach for fusing imaging datasets of the individual components into a common multivariable trajectory. Data fusion revealed low intrinsic dimensionality of posterior patterning and cell fate specification in wild-type embryos. The simple structure that we uncovered allowed us to construct a model of interactions within the posterior patterning regulatory network and make testable predictions about its dynamics at the protein level. The presented data fusion strategy is a step toward establishing a unified framework that would explore how stochastic spatiotemporal signals give rise to highly reproducible morphogenetic outcomes.

Synthetic reconstruction of the hunchback promoter specifies the role of Bicoid, Zelda and Hunchback in the dynamics of its transcription

For over 40 years, the Bicoid-hunchback (Bcd-hb) system in the fruit fly embryo has been used as a model to study how positional information in morphogen concentration gradients is robustly translated into step-like responses. A body of quantitative comparisons between theory and experiment have since questioned the initial paradigm that the sharp hb transcription pattern emerges solely from diffusive biochemical interactions between the Bicoid transcription factor and the gene promoter region. Several alternative mechanisms have been proposed, such as additional sources of positional information, positive feedback from Hb proteins or out-of-equilibrium transcription activation. By using the MS2-MCP RNA-tagging system and analysing in real time, the transcription dynamics of synthetic reporters for Bicoid and/or its two partners Zelda and Hunchback, we show that all the early hb expression pattern features and temporal dynamics are compatible with an equilibrium model with a short decay length Bicoid activity gradient as a sole source of positional information. Meanwhile, Bicoid’s partners speed-up the process by different means: Zelda lowers the Bicoid concentration threshold required for transcriptional activation while Hunchback reduces burstiness and increases the polymerase firing rate.

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.

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.

Transcription Factor Dynamics

To predict transcription, one needs a mechanistic understanding of how the numerous required transcription factors (TFs) explore the nuclear space to find their target genes, assemble, cooperate, and compete with one another. Advances in fluorescence microscopy have made it possible to visualize real-time TF dynamics in living cells, leading to two intriguing observations: first, most TFs contact chromatin only transiently; and second, TFs can assemble into clusters through their intrinsically disordered regions. These findings suggest that highly dynamic events and spatially structured nuclear microenvironments might play key roles in transcription regulation that are not yet fully understood. The emerging model is that while some promoters directly convert TF-binding events into on/off cycles of transcription, many others apply complex regulatory layers that ultimately lead to diverse phenotypic outputs. Cracking this kinetic code is an ongoing and challenging task that is made possible by combining innovative imaging approaches with biophysical models.

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.

Single cell biology-a Keystone Symposia report

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

Molecular competition can shape enhancer activity in the Drosophila embryo

Transgenic reporters allow the measurement of regulatory DNA activity in vivo and consequently have long been useful tools for studying enhancers. Despite their utility, few studies have investigated the effects these reporters may have on the expression of other genes. Understanding these effects is required to accurately interpret reporter data and characterize gene regulatory mechanisms. By measuring the expression of Kruppel (Kr) enhancer reporters in live Drosophila embryos, we find reporters inhibit one another's expression and that of a nearby endogenous gene. Using synthetic transcription factor (TF) binding site arrays, we present evidence that competition for TFs is partially responsible for the observed transcriptional inhibition. We develop a simple thermodynamic model that predicts competition of the measured magnitude specifically when TF binding is restricted to distinct nuclear subregions. Our findings underline an unexpected role of the non-homogenous nature of the nucleus in regulating gene expression.

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

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

Stochastic pausing at latent HIV-1 promoters generates transcriptional bursting

Promoter-proximal pausing of RNA polymerase II is a key process regulating gene expression. In latent HIV-1 cells, it prevents viral transcription and is essential for latency maintenance, while in acutely infected cells the viral factor Tat releases paused polymerase to induce viral expression. Pausing is fundamental for HIV-1, but how it contributes to bursting and stochastic viral reactivation is unclear. Here, we performed single molecule imaging of HIV-1 transcription. We developed a quantitative analysis method that manages multiple time scales from seconds to days and that rapidly fits many models of promoter dynamics. We found that RNA polymerases enter a long-lived pause at latent HIV-1 promoters (>20 minutes), thereby effectively limiting viral transcription. Surprisingly and in contrast to current models, pausing appears stochastic and not obligatory, with only a small fraction of the polymerases undergoing long-lived pausing in absence of Tat. One consequence of stochastic pausing is that HIV-1 transcription occurs in bursts in latent cells, thereby facilitating latency exit and providing a rationale for the stochasticity of viral rebounds.

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.

Real-time single-cell characterization of the eukaryotic transcription cycle reveals correlations between RNA initiation, elongation, and cleavage

The eukaryotic transcription cycle consists of three main steps: initiation, elongation, and cleavage of the nascent RNA transcript. Although each of these steps can be regulated as well as coupled with each other, their in vivo dissection has remained challenging because available experimental readouts lack sufficient spatiotemporal resolution to separate the contributions from each of these steps. Here, we describe a novel application of Bayesian inference techniques to simultaneously infer the effective parameters of the transcription cycle in real time and at the single-cell level using a two-color MS2/PP7 reporter gene and the developing fruit fly embryo as a case study. Our method enables detailed investigations into cell-to-cell variability in transcription-cycle parameters as well as single-cell correlations between these parameters. These measurements, combined with theoretical modeling, suggest a substantial variability in the elongation rate of individual RNA polymerase molecules. We further illustrate the power of this technique by uncovering a novel mechanistic connection between RNA polymerase density and nascent RNA cleavage efficiency. Thus, our approach makes it possible to shed light on the regulatory mechanisms in play during each step of the transcription cycle in individual, living cells at high spatiotemporal resolution.

Live imaging and quantitation of nascent transcription using the MS2/MCP system in the Drosophila embryo

Visualizing transcription live in Drosophila is providing important new insights into the spatiotemporal regulation of transcription. Here, we describe a protocol to visualize and quantitate transcription from gene loci that are tagged with MS2 stem-loop sequences in the Drosophila embryo. MS2 stem-loop sequences are recognized by a coat protein fused to a fluorescent protein and visualized with microscopy. We also describe an analysis pipeline to extract and subsequently quantify transcription dynamics. For complete details on the use and execution of this protocol, please refer to Hoppe et al. (2020).

Kinetic sculpting of the seven stripes of the Drosophila even-skipped gene

We used live imaging to visualize the transcriptional dynamics of the Drosophila melanogaster even-skipped gene at single-cell and high-temporal resolution as its seven stripe expression pattern forms, and developed tools to characterize and visualize how transcriptional bursting varies over time and space. We find that despite being created by the independent activity of five enhancers, even-skipped stripes are sculpted by the same kinetic phenomena: a coupled increase of burst frequency and amplitude. By tracking the position and activity of individual nuclei, we show that stripe movement is driven by the exchange of bursting nuclei from the posterior to anterior stripe flanks. Our work provides a conceptual, theoretical and computational framework for dissecting pattern formation in space and time, and reveals how the coordinated transcriptional activity of individual nuclei shapes complex developmental patterns.

A matter of time: Using dynamics and theory to uncover mechanisms of transcriptional bursting

Eukaryotic transcription generally occurs in bursts of activity lasting minutes to hours; however, state-of-the-art measurements have revealed that many of the molecular processes that underlie bursting, such as transcription factor binding to DNA, unfold on timescales of seconds. This temporal disconnect lies at the heart of a broader challenge in physical biology of predicting transcriptional outcomes and cellular decision-making from the dynamics of underlying molecular processes. Here, we review how new dynamical information about the processes underlying transcriptional control can be combined with theoretical models that predict not only averaged transcriptional dynamics, but also their variability, to formulate testable hypotheses about the molecular mechanisms underlying transcriptional bursting and control.

Signaling Mechanism of Transcriptional Bursting: A Technical Resolution-Independent Study

Simple Summary

Following changing cellular signals, various genes adjust their activities and initiate transcripts with the right rates. The precision of such a transcriptional response has a fundamental role in the survival and development of lives. Quite unexpectedly, gene transcription has been uncovered to occur in sporadic bursts, rather than in a continuous manner. This has raised a provoking issue of how the bursting transmits regulatory signals, and it remains controversial whether the burst size, frequency, or both, take the role of signal transmission. Here, this study showed that only the burst frequency was subject to modulation by activators that carry the regulatory signals. A higher activator concentration led to a larger frequency, whereas the size remains unchanged. When very high, the burst cluster emerged, which may be mistaken as a large burst. This work thus supports the conclusion that transcription regulation is in a “digital” way.

Abstract

Gene transcription has been uncovered to occur in sporadic bursts. However, due to technical difficulties in differentiating individual transcription initiation events, it remains debated as to whether the burst size, frequency, or both are subject to modulation by transcriptional activators. Here, to bypass technical constraints, we addressed this issue by introducing two independent theoretical methods including analytical research based on the classic two-model and information entropy research based on the architecture of transcription apparatus. Both methods connect the signaling mechanism of transcriptional bursting to the characteristics of transcriptional uncertainty (i.e., the differences in transcriptional levels of the same genes that are equally activated). By comparing the theoretical predictions with abundant experimental data collected from published papers, the results exclusively support frequency modulation. To further validate this conclusion, we showed that the data that appeared to support size modulation essentially supported frequency modulation taking into account the existence of burst clusters. This work provides a unified scheme that reconciles the debate on burst signaling.

Quantitative dissection of transcription in development yields evidence for transcription-factor-driven chromatin accessibility

Thermodynamic models of gene regulation can predict transcriptional regulation in bacteria, but in eukaryotes, chromatin accessibility and energy expenditure may call for a different framework. Here, we systematically tested the predictive power of models of DNA accessibility based on the Monod-Wyman-Changeux (MWC) model of allostery, which posits that chromatin fluctuates between accessible and inaccessible states. We dissected the regulatory dynamics of hunchback by the activator Bicoid and the pioneer-like transcription factor Zelda in living Drosophila embryos and showed that no thermodynamic or non-equilibrium MWC model can recapitulate hunchback transcription. Therefore, we explored a model where DNA accessibility is not the result of thermal fluctuations but is catalyzed by Bicoid and Zelda, possibly through histone acetylation, and found that this model can predict hunchback dynamics. Thus, our theory-experiment dialogue uncovered potential molecular mechanisms of transcriptional regulatory dynamics, a key step toward reaching a predictive understanding of developmental decision-making.

Hematopoietic Stem Cell Transcription Factors in Cardiovascular Pathology

Transcription factors as multifaceted modulators of gene expression that play a central role in cell proliferation, differentiation, lineage commitment, and disease progression. They interact among themselves and create complex spatiotemporal gene regulatory networks that modulate hematopoiesis, cardiogenesis, and conditional differentiation of hematopoietic stem cells into cells of cardiovascular lineage. Additionally, bone marrow-derived stem cells potentially contribute to the cardiovascular cell population and have shown potential as a therapeutic approach to treat cardiovascular diseases. However, the underlying regulatory mechanisms are currently debatable. This review focuses on some key transcription factors and associated epigenetic modifications that modulate the maintenance and differentiation of hematopoietic stem cells and cardiac progenitor cells. In addition to this, we aim to summarize different potential clinical therapeutic approaches in cardiac regeneration therapy and recent discoveries in stem cell-based transplantation.

Robust and efficient gene regulation through localized nuclear microenvironments

Developmental enhancers drive gene expression in specific cell types during animal development. They integrate signals from many different sources mediated through the binding of transcription factors, producing specific responses in gene expression. Transcription factors often bind low-affinity sequences for only short durations. How brief, low-affinity interactions drive efficient transcription and robust gene expression is a central question in developmental biology. Localized high concentrations of transcription factors have been suggested as a possible mechanism by which to use these enhancer sites effectively. Here, we discuss the evidence for such transcriptional microenvironments, mechanisms for their formation and the biological consequences of such sub-nuclear compartmentalization for developmental decisions and evolution.

Modulation of the Promoter Activation Rate Dictates the Transcriptional Response to Graded BMP Signaling Levels in the Drosophila Embryo

Morphogen gradients specify cell fates during development, with a classic example being the bone morphogenetic protein (BMP) gradient's conserved role in embryonic dorsal-ventral axis patterning. Here, we elucidate how the BMP gradient is interpreted in the Drosophila embryo by combining live imaging with computational modeling to infer transcriptional burst parameters at single-cell resolution. By comparing burst kinetics in cells receiving different levels of BMP signaling, we show that BMP signaling controls burst frequency by regulating the promoter activation rate. We provide evidence that the promoter activation rate is influenced by both enhancer and promoter sequences, whereas Pol II loading rate is primarily modulated by the enhancer. Consistent with BMP-dependent regulation of burst frequency, the numbers of BMP target gene transcripts per cell are graded across their expression domains. We suggest that graded mRNA output is a general feature of morphogen gradient interpretation and discuss how this can impact on cell-fate decisions.

Shadow enhancers can suppress input transcription factor noise through distinct regulatory logic

Shadow enhancers, groups of seemingly redundant enhancers, are found in a wide range of organisms and are critical for robust developmental patterning. However, their mechanism of action is unknown. We hypothesized that shadow enhancers drive consistent expression levels by buffering upstream noise through a separation of transcription factor (TF) inputs at the individual enhancers. By measuring the transcriptional dynamics of several Kruppel shadow enhancer configurations in live Drosophila embryos, we showed that individual member enhancers act largely independently. We found that TF fluctuations are an appreciable source of noise that the shadow enhancer pair can better buffer than duplicated enhancers. The shadow enhancer pair is also uniquely able to maintain low levels of expression noise across a wide range of temperatures. A stochastic model demonstrated the separation of TF inputs is sufficient to explain these findings. Our results suggest the widespread use of shadow enhancers is partially due to their noise suppressing ability.

Effects of cell cycle variability on lineage and population measurements of messenger RNA abundance

Many models of gene expression do not explicitly incorporate a cell cycle description. Here, we derive a theory describing how messenger RNA (mRNA) fluctuations for constitutive and bursty gene expression are influenced by stochasticity in the duration of the cell cycle and the timing of DNA replication. Analytical expressions for the moments show that omitting cell cycle duration introduces an error in the predicted mean number of mRNAs that is a monotonically decreasing function of η, which is proportional to the ratio of the mean cell cycle duration and the mRNA lifetime. By contrast, the error in the variance of the mRNA distribution is highest for intermediate values of η consistent with genome-wide measurements in many organisms. Using eukaryotic cell data, we estimate the errors in the mean and variance to be at most 3% and 25%, respectively. Furthermore, we derive an accurate negative binomial mixture approximation to the mRNA distribution. This indicates that stochasticity in the cell cycle can introduce fluctuations in mRNA numbers that are similar to the effect of bursty transcription. Finally, we show that for real experimental data, disregarding cell cycle stochasticity can introduce errors in the inference of transcription rates larger than 10%.