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

Transcription templated assembly of the nucleolus in the C. elegans embryo

The nucleolus is a multicomponent structure made of RNA and proteins that serves as the site of ribosome biogenesis within the nucleus. It has been extensively studied as a prototype of a biomolecular condensate whose assembly is driven by phase separation. While the steady-state size of the nucleolus is quantitatively accounted for by the thermodynamics of phase separation, we show that experimental measurements of the assembly dynamics are inconsistent with a simple model of a phase-separating system relaxing to its equilibrium state. Instead, we show that the dynamics are well described by a model in which the transcription of ribosomal RNA actively drives nucleolar assembly. We find that our model of active transcription-templated assembly quantitatively accounts for the rapid kinetics observed in early embryos at different developmental stages, and for different RNAi perturbations of embryo size. Our model predicts a scaling of the time to assembly with the volume of the nucleus to the one-third power, which is confirmed by experimental data. Our study highlights the role of active processes such as transcription in controlling the placement and timing of assembly of membraneless organelles.

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.

Chromatin accessibility in the Drosophila embryo is determined by transcription factor pioneering and enhancer activation

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation.

Shared retinoic acid responsive enhancers coordinately regulate nascent transcription of Hoxb coding and non-coding RNAs in the developing mouse neural tube

Signaling pathways regulate the patterns of Hox gene expression that underlie their functions in the specification of axial identity. Little is known about the properties of cis-regulatory elements and underlying transcriptional mechanisms that integrate graded signaling inputs to coordinately control Hox expression. Here, we optimized a single molecule fluorescent in situ hybridization (smFISH) technique with probes spanning introns to evaluate how three shared retinoic acid response element (RARE)-dependent enhancers in the Hoxb cluster regulate patterns of nascent transcription in vivo at the level of single cells in wild-type and mutant embryos. We predominately detect nascent transcription of only a single Hoxb gene in each cell, with no evidence for simultaneous co-transcriptional coupling of all or specific subsets of genes. Single and/or compound RARE mutations indicate that each enhancer differentially impacts global and local patterns of nascent transcription, suggesting that selectivity and competitive interactions between these enhancers is important to robustly maintain the proper levels and patterns of nascent Hoxb transcription. This implies that rapid and dynamic regulatory interactions potentiate transcription of genes through combined inputs from these enhancers in coordinating the retinoic acid response.

Dynamic epistasis analysis reveals how chromatin remodeling regulates transcriptional bursting

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

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.

The transcription factor RUNX1 affects the maturation of porcine oocytes via the BMP15/TGF-β signaling pathway

Bone morphogenetic protein 15 (BMP15) is specifically expressed in oocytes in pigs at all stages from early stages to ovulation and has an important role in oocyte maturation. However, there are few reports on the molecular mechanisms by which BMP15 affects oocyte maturation. In this study, we identified the core promoter region of BMP15 using a dual luciferase activity assay and successfully predicted the DNA binding motif of the transcription factor RUNX1. The effect of BMP15 and RUNX1 on oocyte maturation was examined using the first polar body extrusion rate, a reactive oxygen species (ROS) assay and total glutathione (GSH) content at three time points of 12, 24 and 48 h of in vitro culture of porcine isolated oocytes. Subsequently, the effect of the transcription factor RUNX1 on the TGF-β signaling pathway (BMPR1B and ALK5) was further verified using RT-qPCR and Western blotting. We found that the overexpression of BMP15 significantly increased the first polar body extrusion rate (P < 0.01) and total glutathione content of oocytes cultured in vitro for 24 h and decreased reactive oxygen levels (P < 0.01), whereas interference with BMP15 decreased the first polar body extrusion rate (P < 0.01), increased reactive oxygen levels in oocytes cultured in vitro for 24 h (P < 0.01), and decreased glutathione content (P < 0.01). The dual luciferase activity assay and online software prediction showed that RUNX1 is a potential transcription factor binding to the core promoter region (-1203/-1423 bp) of BMP15. Overexpression of RUNX1 significantly increased the expression of BMP15 and oocyte maturation rate, while inhibition of RUNX1 decreased the expression of BMP15 and the oocyte maturation rate. Moreover, the expression of BMPR1B and ALK5 in the TGF-β signaling pathway increased significantly after overexpression of RUNX1, whereas their expression decreased after inhibition of RUNX1. Overall, our results suggest that the transcription factor RUNX1 positively regulates the expression of BMP15 and influences oocyte maturation through the TGF-β signaling pathway. This study provides a theoretical basis for further complementing the BMP15/TGF-β signaling pathway to regulate mammalian oocyte maturation.

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.

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.

Dynamic chromatin accessibility deploys heterotypic cis/trans-acting factors driving stomatal cell-fate commitment

Chromatin architecture and transcription factor (TF) binding underpin cell-fate specification during development, but their mutual regulatory relationships remain unclear. Here we report an atlas of dynamic chromatin landscapes during stomatal cell-lineage progression, in which sequential cell-state transitions are governed by lineage-specific bHLH TFs. Major reprogramming of chromatin accessibility occurs at the proliferation-to-differentiation transition. We discover novel co-cis regulatory elements (CREs) signifying the early precursor stage, BBR/BPC (GAGA) and bHLH (E-box) motifs, where master-regulatory bHLH TFs, SPEECHLESS and MUTE, consecutively bind to initiate and terminate the proliferative state, respectively. BPC TFs complex with MUTE to repress SPEECHLESS expression through a local deposition of repressive histone marks. We elucidate the mechanism by which cell-state-specific heterotypic TF complexes facilitate cell-fate commitment by recruiting chromatin modifiers via key co-CREs.

How does an organism extract relevant information from transcription factor concentrations?

How does an organism regulate its genes? The involved regulation typically occurs in terms of a signal processing chain: an externally applied stimulus or a maternally supplied transcription factor leads to the expression of some downstream genes, which, in turn, are transcription factors for further genes. Especially during development, these transcription factors are frequently expressed in amounts where noise is still important; yet, the signals that they provide must not be lost in the noise. Thus, the organism needs to extract exactly relevant information in the signal. New experimental approaches involving single-molecule measurements at high temporal precision as well as increased precision in manipulations directly on the genome are allowing us to tackle this question anew. These new experimental advances mean that also from the theoretical side, theoretical advances should be possible. In this review, I will describe, specifically on the example of fly embryo gene regulation, how theoretical approaches, especially from inference and information theory, can help in understanding gene regulation. To do so, I will first review some more traditional theoretical models for gene regulation, followed by a brief discussion of information-theoretical approaches and when they can be applied. I will then introduce early fly development as an exemplary system where such information-theoretical approaches have traditionally been applied and can be applied; I will specifically focus on how one such method, namely the information bottleneck approach, has recently been used to infer structural features of enhancer architecture.

Eukaryotic gene regulation at equilibrium, or non?

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

Structures and consequences of pioneer factor binding to nucleosomes

Pioneer transcription factors are able to bind a partially exposed motif on the surface of a nucleosome, enabling the proteins to target sites in silent regions of chromatin that have been compacted by linker histone. The targeting of nucleosomal DNA by pioneer factors has been observed in vitro and in vivo, where binding can promote local nucleosome exposure that allows other transcription factors, nucleosome remodelers, and histone modifiers to engage the chromatin and elicit gene activation or further repression. Pioneer factors thereby establish new gene expression programs during cell fate changes that occur during embryonic development, regeneration, and cancer. Here, we review recent biophysical studies that reveal the structural features and strategies used by pioneer factors to accomplish nucleosome binding and the consequential changes to nucleosomes that can lead to DNA accessibility.

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.

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.

Notch-dependent and -independent transcription are modulated by tissue movements at gastrulation

Cells sense and integrate external information from diverse sources that include mechanical cues. Shaping of tissues during development may thus require coordination between mechanical forces from morphogenesis and cell-cell signalling to confer appropriate changes in gene expression. By live-imaging Notch-induced transcription in real time, we have discovered that morphogenetic movements during Drosophila gastrulation bring about an increase in activity-levels of a Notch-responsive enhancer. Mutations that disrupt the timing of gastrulation resulted in concomitant delays in transcription up-regulation that correlated with the start of mesoderm invagination. As a similar gastrulation-induced effect was detected when transcription was elicited by the intracellular domain NICD, it cannot be attributed to forces exerted on Notch receptor activation. A Notch-independent vnd enhancer also exhibited a modest gastrulation-induced activity increase in the same stripe of cells. Together, these observations argue that gastrulation-associated forces act on the nucleus to modulate transcription levels. This regulation was uncoupled when the complex linking the nucleoskeleton and cytoskeleton (LINC) was disrupted, indicating a likely conduit. We propose that the coupling between tissue-level mechanics, arising from gastrulation, and enhancer activity represents a general mechanism for ensuring correct tissue specification during development and that Notch-dependent enhancers are highly sensitive to this regulation.

From complex datasets to predictive models of embryonic development

Modern studies of embryogenesis are increasingly quantitative, powered by rapid advances in imaging, sequencing and genome manipulation technologies. Deriving mechanistic insights from the complex datasets generated by these new tools requires systematic approaches for data-driven analysis of the underlying developmental processes. Here, we use data from our work on signal-dependent gene repression in the Drosophila embryo to illustrate how computational models can compactly summarize quantitative results of live imaging, chromatin immunoprecipitation and optogenetic perturbation experiments. The presented computational approach is ideally suited for integrating rapidly accumulating quantitative data and for guiding future studies of embryogenesis.

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.

Logic and lineage impacts on functional transcription factor deployment for T-cell fate commitment

Transcription factors are the major agents that read the regulatory sequence information in the genome to initiate changes in expression of specific genes, both in development and in physiological activation responses. Their actions depend on site-specific DNA binding and are largely guided by their individual DNA target sequence specificities. However, their action is far more conditional in a real developmental context than would be expected for simple reading of local genomic DNA sequence, which is common to all cells in the organism. They are constrained by slow-changing chromatin states and by interactions with other transcription factors, which affect their occupancy patterns of potential sites across the genome. These mechanisms lead to emergent discontinuities in function even for transcription factors with minimally changing expression. This is well revealed by diverse lineages of blood cells developing throughout life from hematopoietic stem cells, which use overlapping combinations of transcription factors to drive strongly divergent gene regulation programs. Here, using development of T lymphocytes from hematopoietic multipotent progenitor cells as a focus, recent evidence is reviewed on how binding specificity and dynamics, transcription factor cooperativity, and chromatin state changes impact the effective regulatory functions of key transcription factors including PU.1, Runx1, Notch-RBPJ, and Bcl11b.

Gene Regulation and Cellular Metabolism: An Essential Partnership

It is recognized that cell metabolism is tightly connected to other cellular processes such as regulation of gene expression. Metabolic pathways not only provide the precursor molecules necessary for gene expression, but they also provide ATP, the primary fuel driving gene expression. However, metabolic conditions are highly variable since nutrient uptake is not a uniform process. Thus, cells must continually calibrate gene expression to their changing metabolite and energy budgets. This review discusses recent advances in understanding the molecular and biophysical mechanisms that connect metabolism and gene regulation as cells navigate their growth, proliferation, and differentiation. Particular focus is given to these mechanisms in the context of organismal development.

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.