A single-molecule view of transcription reveals convoys of RNA polymerases and multi-scale bursting

The Polycomb system sustains promoters in a deep OFF state by limiting pre-initiation complex formation to counteract transcription

The Polycomb system has fundamental roles in regulating gene expression during mammalian development. However, how it controls transcription to enable gene repression has remained enigmatic. Here, using rapid degron-based depletion coupled with live-cell transcription imaging and single-particle tracking, we show how the Polycomb system controls transcription in single cells. We discover that the Polycomb system is not a constitutive block to transcription but instead sustains a long-lived deep promoter OFF state, which limits the frequency with which the promoter can enter into a transcribing state. We demonstrate that Polycomb sustains this deep promoter OFF state by counteracting the binding of factors that enable early transcription pre-initiation complex formation and show that this is necessary for gene repression. Together, these important discoveries provide a rationale for how the Polycomb system controls transcription and suggests a universal mechanism that could enable the Polycomb system to constrain transcription across diverse cellular contexts.

Transcriptional bursting: from fundamentals to novel insights

Transcription occurs as irregular bursts in a very wide range of systems, including numerous different species and many genes within these. In this review, we examine the underlying theories, discuss how these relate to experimental measurements, and explore some of the discrepancies that have emerged among various studies. Finally, we consider more recent works that integrate novel concepts, such as the involvement of biomolecular condensates in enhancer-promoter interactions and their effects on the dynamics of transcriptional bursting.

The method in the madness: Transcriptional control from stochastic action at the single-molecule scale

Cell states result from the ordered activation of gene expression by transcription factors. Transcription factors face opposing design constraints: they need to be dynamic to trigger rapid cell state transitions, but also stable enough to maintain terminal cell identities indefinitely. Recent progress in live-cell single-molecule microscopy has helped define the biophysical principles underlying this paradox. Beyond transcription factor activity, single-molecule experiments have revealed that at nearly every level of transcription regulation, control emerges from multiple short-lived stochastic interactions, rather than deterministic, stable interactions typical of other biochemical pathways. This architecture generates consistent outcomes that can be rapidly choreographed. Here, we highlight recent results that demonstrate how order in transcription regulation emerges from the apparent molecular-scale chaos and discuss remaining conceptual challenges.

Single gene analysis in yeast suggests nonequilibrium regulatory dynamics for transcription

Fluctuations in the initiation rate of transcription, the first step in gene expression, ensue from the stochastic behavior of the molecular process that controls transcription. In steady state, the regulatory process is often assumed to operate reversibly, i.e., in equilibrium. However, reversibility imposes fundamental limits to information processing. For instance, the assumption of equilibrium is difficult to square with the precision with which the regulatory process executes its task in eukaryotes. Here we provide evidence - from microscopic analyses of the transcription dynamics at a single gene copy of yeast - that the regulatory process for transcription is cyclic and irreversible (out of equilibrium). The necessary coupling to reservoirs of free energy occurs via sequence-specific transcriptional activators and the recruitment, in part, of ATP-dependent chromatin remodelers. Our findings may help explain how eukaryotic cells reconcile the dual but opposing requirements for fast regulatory kinetics and high regulatory specificity.

Enhancer selectivity in space and time: from enhancer-promoter interactions to promoter activation

The primary regulators of metazoan gene expression are enhancers, originally functionally defined as DNA sequences that can activate transcription at promoters in an orientation-independent and distance-independent manner. Despite being crucial for gene regulation in animals, what mechanisms underlie enhancer selectivity for promoters, and more fundamentally, how enhancers interact with promoters and activate transcription, remain poorly understood. In this Review, we first discuss current models of enhancer-promoter interactions in space and time and how enhancers affect transcription activation. Next, we discuss different mechanisms that mediate enhancer selectivity, including repression, biochemical compatibility and regulation of 3D genome structure. Through 3D polymer simulations, we illustrate how the ability of 3D genome folding mechanisms to mediate enhancer selectivity strongly varies for different enhancer-promoter interaction mechanisms. Finally, we discuss how recent technical advances may provide new insights into mechanisms of enhancer-promoter interactions and how technical biases in methods such as Hi-C and Micro-C and imaging techniques may affect their interpretation.

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.

The HIV-1 Transcriptional Program: From Initiation to Elongation Control

A large body of work in the last four decades has revealed the key pillars of HIV-1 transcription control at the initiation and elongation steps. Here, I provide a recount of this collective knowledge starting with the genomic elements (DNA and nascent TAR RNA stem-loop) and transcription factors (cellular and the viral transactivator Tat), and later transitioning to the assembly and regulation of transcription initiation and elongation complexes, and the role of chromatin structure. Compelling evidence support a core HIV-1 transcriptional program regulated by the sequential and concerted action of cellular transcription factors and Tat to promote initiation and sustain elongation, highlighting the efficiency of a small virus to take over its host to produce the high levels of transcription required for viral replication. I summarize new advances including the use of CRISPR-Cas9, genetic tools for acute factor depletion, and imaging to study transcriptional dynamics, bursting and the progression through the multiple phases of the transcriptional cycle. Finally, I describe current challenges to future major advances and discuss areas that deserve more attention to both bolster our basic knowledge of the core HIV-1 transcriptional program and open up new therapeutic opportunities.

Protocol to quantify the rate of RNA polymerase II elongation with MS2/MCP- and PP7/PCP-based live-imaging systems in early Drosophila embryos

The MS2-PP7 two-color live-imaging system provides insights into the spatiotemporal dynamics of nascent transcripts at tagged loci. Here, we present a protocol to quantitatively measure the rate of RNA polymerase II elongation for each actively transcribing nucleus in living Drosophila embryos. The elongation rate is calculated by measuring the effective distance and the time elapsed between MS2 and PP7 trajectories. We describe steps for preparing embryo samples, performing live imaging, and measuring the elongation rate. For complete details on the use and execution of this protocol, please refer to Keller et al.1.

Identification of molecular determinants of gene-specific bursting patterns by high-throughput imaging screens

Stochastic transcriptional bursting is a universal property of active genes. While different genes exhibit distinct bursting patterns, the molecular mechanisms for gene-specific stochastic bursting are largely unknown. We have developed and applied a high-throughput-imaging based screening strategy to identify cellular factors and molecular mechanisms that determine the bursting behavior of human genes. Focusing on epigenetic regulators, we find that protein acetylation is a strong acute modulator of burst frequency, burst size and heterogeneity of bursting. Acetylation globally affects the Off-time of genes but has gene-specific effects on the On-time. Yet, these effects are not strongly linked to promoter acetylation, which do not correlate with bursting properties, and forced promoter acetylation has variable effects on bursting. Instead, we demonstrate acetylation of the Integrator complex as a key determinant of gene bursting. Specifically, we find that elevated Integrator acetylation decreases bursting frequency. Taken together our results suggest a prominent role of non-histone proteins in determining gene bursting properties, and they identify histone-independent acetylation of a transcription cofactor as an allosteric modulator of bursting via a far-downstream bursting checkpoint.

Solving stochastic gene-expression models using queueing theory: A tutorial review

Stochastic models of gene expression are typically formulated using the chemical master equation, which can be solved exactly or approximately using a repertoire of analytical methods. Here, we provide a tutorial review of an alternative approach based on queueing theory that has rarely been used in the literature of gene expression. We discuss the interpretation of six types of infinite-server queues from the angle of stochastic single-cell biology and provide analytical expressions for the stationary and nonstationary distributions and/or moments of mRNA/protein numbers and bounds on the Fano factor. This approach may enable the solution of complex models that have hitherto evaded analytical solution.

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

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

Mechanisms and Functions of the RNA Polymerase II General Transcription Machinery during the Transcription Cycle

Central to the development and survival of all organisms is the regulation of gene expression, which begins with the process of transcription catalyzed by RNA polymerases. During transcription of protein-coding genes, the general transcription factors (GTFs) work alongside RNA polymerase II (Pol II) to assemble the preinitiation complex at the transcription start site, open the promoter DNA, initiate synthesis of the nascent messenger RNA, transition to productive elongation, and ultimately terminate transcription. Through these different stages of transcription, Pol II is dynamically phosphorylated at the C-terminal tail of its largest subunit, serving as a control mechanism for Pol II elongation and a signaling/binding platform for co-transcriptional factors. The large number of core protein factors participating in the fundamental steps of transcription add dense layers of regulation that contribute to the complexity of temporal and spatial control of gene expression within any given cell type. The Pol II transcription system is highly conserved across different levels of eukaryotes; however, most of the information here will focus on the human Pol II system. This review walks through various stages of transcription, from preinitiation complex assembly to termination, highlighting the functions and mechanisms of the core machinery that participates in each stage.

Time will tell: comparing timescales to gain insight into transcriptional bursting

Recent imaging studies have captured the dynamics of regulatory events of transcription inside living cells. These events include transcription factor (TF) DNA binding, chromatin remodeling and modification, enhancer-promoter (E-P) proximity, cluster formation, and preinitiation complex (PIC) assembly. Together, these molecular events culminate in stochastic bursts of RNA synthesis, but their kinetic relationship remains largely unclear. In this review, we compare the timescales of upstream regulatory steps (input) with the kinetics of transcriptional bursting (output) to generate mechanistic models of transcription dynamics in single cells. We highlight open questions and potential technical advances to guide future endeavors toward a quantitative and kinetic understanding of transcription regulation.

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

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

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.

Exploring the role of transcriptional and post-transcriptional processes in mRNA co-expression

Co-expression of two or more genes at the single-cell level is usually associated with functional co-regulation. While mRNA co-expression-measured as the correlation in mRNA levels-can be influenced by both transcriptional and post-transcriptional events, transcriptional regulation is typically considered dominant. We review and connect the literature describing transcriptional and post-transcriptional regulation of co-expression. To enhance our understanding, we integrate four datasets spanning single-cell gene expression data, single-cell promoter activity data and individual transcript half-lives. Confirming expectations, we find that positive co-expression necessitates promoter coordination and similar mRNA half-lives. Surprisingly, negative co-expression is favored by differences in mRNA half-lives, contrary to initial predictions from stochastic simulations. Notably, this association manifests specifically within clusters of genes. We further observe a striking compensation between promoter coordination and mRNA half-lives, which additional stochastic simulations suggest might give rise to the observed co-expression patterns. These findings raise intriguing questions about the functional advantages conferred by this compensation between distal kinetic steps.

Regulation of the RNA polymerase II pre-initiation complex by its associated coactivators

The RNA polymerase II (Pol II) pre-initiation complex (PIC) is a critical node in eukaryotic transcription regulation, and its formation is the major rate-limiting step in transcriptional activation. Diverse cellular signals borne by transcriptional activators converge on this large, multiprotein assembly and are transduced via intermediary factors termed coactivators. Cryogenic electron microscopy, multi-omics and single-molecule approaches have recently offered unprecedented insights into both the structure and cellular functions of the PIC and two key PIC-associated coactivators, Mediator and TFIID. Here, we review advances in our understanding of how Mediator and TFIID interact with activators and affect PIC formation and function. We also discuss how their functions are influenced by their chromatin environment and selected cofactors. We consider how, through its multifarious interactions and functionalities, a Mediator-containing and TFIID-containing PIC can yield an integrated signal processing system with the flexibility to determine the unique temporal and spatial expression pattern of a given gene.

Mathematical Models of HIV-1 Dynamics, Transcription, and Latency

HIV-1 latency is a major barrier to curing infections with antiretroviral therapy and, consequently, to eliminating the disease globally. The establishment, maintenance, and potential clearance of latent infection are complex dynamic processes and can be best described with the help of mathematical models followed by experimental validation. Here, we review the use of viral dynamics models for HIV-1, with a focus on applications to the latent reservoir. Such models have been used to explain the multi-phasic decay of viral load during antiretroviral therapy, the early seeding of the latent reservoir during acute infection and the limited inflow during treatment, the dynamics of viral blips, and the phenomenon of post-treatment control. Finally, we discuss that mathematical models have been used to predict the efficacy of potential HIV-1 cure strategies, such as latency-reversing agents, early treatment initiation, or gene therapies, and to provide guidance for designing trials of these novel interventions.

Following the Birth, Life, and Death of mRNAs in Single Cells

Recent advances in single-molecule imaging of mRNAs in fixed and living cells have enabled the lives of mRNAs to be studied with unprecedented spatial and temporal detail. These approaches have moved beyond simply being able to observe specific events and have begun to allow an understanding of how regulation is coupled between steps in the mRNA life cycle. Additionally, these methodologies are now being applied in multicellular systems and animals to provide more nuanced insights into the physiological regulation of RNA metabolism.

Transcriptional Stochasticity as a Key Aspect of HIV-1 Latency

This review summarizes current advances in the role of transcriptional stochasticity in HIV-1 latency, which were possible in a large part due to the development of single-cell approaches. HIV-1 transcription proceeds in bursts of RNA production, which stem from the stochastic switching of the viral promoter between ON and OFF states. This switching is caused by random binding dynamics of transcription factors and nucleosomes to the viral promoter and occurs at several time scales from minutes to hours. Transcriptional bursts are mainly controlled by the core transcription factors TBP, SP1 and NF-κb, the chromatin status of the viral promoter and RNA polymerase II pausing. In particular, spontaneous variability in the promoter chromatin creates heterogeneity in the response to activators such as TNF-α, which is then amplified by the Tat feedback loop to generate high and low viral transcriptional states. This phenomenon is likely at the basis of the partial and stochastic response of latent T cells from HIV-1 patients to latency-reversing agents, which is a barrier for the development of shock-and-kill strategies of viral eradication. A detailed understanding of the transcriptional stochasticity of HIV-1 and the possibility to precisely model this phenomenon will be important assets to develop more effective therapeutic strategies.

Single-cell nascent RNA sequencing using click-chemistry unveils coordinated transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1-5. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations6-9. However, fundamental questions in the temporal regulation of transcription and enhancer-gene synchrony remain unanswered primarily due to the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq - a novel single-cell nascent RNA sequencing assay using click-chemistry - and unveil the coordinated transcription throughout the genome. scGRO-seq demonstrates the episodic nature of transcription, and estimates burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells. It reveals the co-transcription of functionally related genes and leverages the replication-dependent non-polyadenylated histone genes transcription to elucidate cell-cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq identifies networks of enhancers and genes and indicates that the bursting of transcription at super-enhancers precedes the burst from associated genes. By imparting insights into the dynamic nature of transcription and the origin and propagation of transcription signals, scGRO-seq demonstrates its unique ability to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

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.

Relating the Biogenesis and Function of P Bodies in Drosophila to Human Disease

Drosophila has been a premier model organism for over a century and many discoveries in flies have furthered our understanding of human disease. Flies have been successfully applied to many aspects of health-based research spanning from behavioural addiction, to dysplasia, to RNA dysregulation and protein misfolding. Recently, Drosophila tissues have been used to study biomolecular condensates and their role in multicellular systems. Identified in a wide range of plant and animal species, biomolecular condensates are dynamic, non-membrane-bound sub-compartments that have been observed and characterised in the cytoplasm and nuclei of many cell types. Condensate biology has exciting research prospects because of their diverse roles within cells, links to disease, and potential for therapeutics. In this review, we will discuss processing bodies (P bodies), a conserved biomolecular condensate, with a particular interest in how Drosophila can be applied to advance our understanding of condensate biogenesis and their role in disease.

Regulation of Pol II Pausing during Daily Gene Transcription in Mouse Liver

Cell autonomous circadian oscillation is present in central and various peripheral tissues. The intrinsic tissue clock and various extrinsic cues drive gene expression rhythms. Transcription regulation is thought to be the main driving force for gene rhythms. However, how transcription rhythms arise remains to be fully characterized due to the fact that transcription is regulated at multiple steps. In particular, Pol II recruitment, pause release, and premature transcription termination are critical regulatory steps that determine the status of Pol II pausing and transcription output near the transcription start site (TSS) of the promoter. Recently, we showed that Pol II pausing exhibits genome-wide changes during daily transcription in mouse liver. In this article, we review historical as well as recent findings on the regulation of transcription rhythms by the circadian clock and other transcription factors, and the potential limitations of those results in explaining rhythmic transcription at the TSS. We then discuss our results on the genome-wide characteristics of daily changes in Pol II pausing, the possible regulatory mechanisms involved, and their relevance to future research on circadian transcription regulation.

The Mediator complex regulates enhancer-promoter interactions

Enhancer-mediated gene activation generally requires physical proximity between enhancers and their target gene promoters. However, the molecular mechanisms by which interactions between enhancers and promoters are formed are not well understood. Here, we investigate the function of the Mediator complex in the regulation of enhancer-promoter interactions, by combining rapid protein depletion and high-resolution MNase-based chromosome conformation capture approaches. We show that depletion of Mediator leads to reduced enhancer-promoter interaction frequencies, which are associated with a strong decrease in gene expression. In addition, we find increased interactions between CTCF-binding sites upon Mediator depletion. These changes in chromatin architecture are associated with a redistribution of the Cohesin complex on chromatin and a reduction in Cohesin occupancy at enhancers. Together, our results indicate that the Mediator and Cohesin complexes contribute to enhancer-promoter interactions and provide insights into the molecular mechanisms by which communication between enhancers and promoters is regulated.

An expanded view of transcription

A new method expands chromatin to provide detailed images of transcription sites.

Bayesian analysis dissects kinetic modulation during non-stationary gene expression

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

Transcriptional bursting: stochasticity in deterministic development

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

Understanding spatiotemporal coupling of gene expression using single molecule RNA imaging technologies

Across all kingdoms of life, gene regulatory mechanisms underlie cellular adaptation to ever-changing environments. Regulation of gene expression adjusts protein synthesis and, in turn, cellular growth. Messenger RNAs are key molecules in the process of gene expression. Our ability to quantitatively measure mRNA expression in single cells has improved tremendously over the past decades. This revealed an unexpected coordination between the steps that control the life of an mRNA, from transcription to degradation. Here, we provide an overview of the state-of-the-art imaging approaches for measurement and quantitative understanding of gene expression, starting from the early visualizations of single genes by electron microscopy to current fluorescence-based approaches in single cells, including live-cell RNA-imaging approaches to FISH-based spatial transcriptomics across model organisms. We also highlight how these methods have shaped our current understanding of the spatiotemporal coupling between transcriptional and post-transcriptional events in prokaryotes. We conclude by discussing future challenges of this multidisciplinary field.Abbreviations: mRNA: messenger RNA; rRNA: ribosomal rDNA; tRNA: transfer RNA; sRNA: small RNA; FISH: fluorescence in situ hybridization; RNP: ribonucleoprotein; smFISH: single RNA molecule FISH; smiFISH: single molecule inexpensive FISH; HCR-FISH: Hybridization Chain-Reaction-FISH; RCA: Rolling Circle Amplification; seqFISH: Sequential FISH; MERFISH: Multiplexed error robust FISH; UTR: Untranslated region; RBP: RNA binding protein; FP: fluorescent protein; eGFP: enhanced GFP, MCP: MS2 coat protein; PCP: PP7 coat protein; MB: Molecular beacons; sgRNA: single guide RNA.

MEDIATOR SUBUNIT17 is required for transcriptional optimization of root system architecture in Arabidopsis

Sucrose and auxin are well-known determinants of root system architecture (RSA). However, the factors that connect the signaling pathways evoked by these two critical factors during root development are poorly understood. In this study, we report the role of MEDIATOR SUBUNIT17 (MED17) in RSA and its involvement in the transcriptional integration of sugar and auxin signaling pathways in Arabidopsis (Arabidopsis thaliana). Sucrose regulates root meristem activation through the TARGET OF RAPAMYCIN-E2 PROMOTER BINDING FACTOR A (TOR-E2FA) pathway, and auxin regulates lateral root (LR) development through AUXIN RESPONSE FACTOR-LATERAL ORGAN BOUNDARIES DOMAIN (ARF-LBDs). Both sucrose and auxin play a vital role during primary and LR development. However, there is no clarity on how sucrose is involved in the ARF-dependent regulation of auxin-responsive genes. This study establishes MED17 as a nodal point to connect sucrose and auxin signaling. Transcription of MED17 was induced by sucrose in an E2FA/B-dependent manner. Moreover, E2FA/B interacted with MED17, which can aid in the recruitment of the Mediator complex on the target promoters. Interestingly, E2FA/B and MED17 also occupied the promoter of ARF7, but not ARF19, leading to ARF7 expression, which then activates auxin signaling and thus initiates LR development. MED17 also activated cell division in the root meristem by occupying the promoters of cell-cycle genes, thus regulating their transcription. Thus, MED17 plays an important role in relaying the transcriptional signal from sucrose to auxin-responsive and cell-cycle genes to regulate primary and lateral root development, highlighting the role of the Mediator as the transcriptional processor for optimal root system architecture in Arabidopsis.

Analysis and design of single-cell experiments to harvest fluctuation information while rejecting measurement noise

Introduction: Despite continued technological improvements, measurement errors always reduce or distort the information that any real experiment can provide to quantify cellular dynamics. This problem is particularly serious for cell signaling studies to quantify heterogeneity in single-cell gene regulation, where important RNA and protein copy numbers are themselves subject to the inherently random fluctuations of biochemical reactions. Until now, it has not been clear how measurement noise should be managed in addition to other experiment design variables (e.g., sampling size, measurement times, or perturbation levels) to ensure that collected data will provide useful insights on signaling or gene expression mechanisms of interest. Methods: We propose a computational framework that takes explicit consideration of measurement errors to analyze single-cell observations, and we derive Fisher Information Matrix (FIM)-based criteria to quantify the information value of distorted experiments. Results and Discussion: We apply this framework to analyze multiple models in the context of simulated and experimental single-cell data for a reporter gene controlled by an HIV promoter. We show that the proposed approach quantitatively predicts how different types of measurement distortions affect the accuracy and precision of model identification, and we demonstrate that the effects of these distortions can be mitigated through explicit consideration during model inference. We conclude that this reformulation of the FIM could be used effectively to design single-cell experiments to optimally harvest fluctuation information while mitigating the effects of image distortion.

cis-Regulatory Elements in Plant Development, Adaptation, and Evolution

cis-Regulatory elements encode the genomic blueprints that ensure the proper spatiotemporal patterning of gene expression necessary for appropriate development and responses to the environment. Accumulating evidence implicates changes to gene expression as a major source of phenotypic novelty in eukaryotes, including acute phenotypes such as disease and cancer in mammals. Moreover, genetic and epigenetic variation affecting cis-regulatory sequences over longer evolutionary timescales has become a recurring theme in studies of morphological divergence and local adaptation. Here, we discuss the functions of and methods used to identify various classes of cis-regulatory elements, as well as their role in plant development and response to the environment. We highlight opportunities to exploit cis-regulatory variants underlying plant development and environmental responses for crop improvement efforts. Although a comprehensive understanding of cis-regulatory mechanisms in plants has lagged behind that in animals, we showcase several breakthrough findings that have profoundly influenced plant biology and shaped the overall understanding of transcriptional regulation in eukaryotes.

Women’s contribution in understanding how topoisomerases, supercoiling, and transcription control genome organization

One of the biggest paradoxes in biology is that human genome is roughly 2 m long, while the nucleus containing it is almost one million times smaller. To fit into the nucleus, DNA twists, bends and folds into several hierarchical levels of compaction. Still, DNA has to maintain a high degree of accessibility to be readily replicated and transcribed by proteins. How compaction and accessibility co-exist functionally in human cells is still a matter of debate. Here, we discuss how the torsional stress of the DNA helix acts as a buffer, regulating both chromatin compaction and accessibility. We will focus on chromatin supercoiling and on the emerging role of topoisomerases as pivotal regulators of genome organization. We will mainly highlight the major breakthrough studies led by women, with the intention of celebrating the work of this group that remains a minority within the scientific community.

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.

Functional partitioning of transcriptional regulators by patterned charge blocks

Components of transcriptional machinery are selectively partitioned into specific condensates, often mediated by protein disorder, yet we know little about how this specificity is achieved. Here, we show that condensates composed of the intrinsically disordered region (IDR) of MED1 selectively partition RNA polymerase II together with its positive allosteric regulators while excluding negative regulators. This selective compartmentalization is sufficient to activate transcription and is required for gene activation during a cell-state transition. The IDRs of partitioned proteins are necessary and sufficient for selective compartmentalization and require alternating blocks of charged amino acids. Disrupting this charge pattern prevents partitioning, whereas adding the pattern to proteins promotes partitioning with functional consequences for gene activation. IDRs with similar patterned charge blocks show similar partitioning and function. These findings demonstrate that disorder-mediated interactions can selectively compartmentalize specific functionally related proteins from a complex mixture of biomolecules, leading to regulation of a biochemical pathway.

From structure to molecular condensates: emerging mechanisms for Mediator function

Mediator is a large modular protein assembly whose function as a coactivator of transcription is conserved in all eukaryotes. The Mediator complex can integrate and relay signals from gene-specific activators bound at enhancers to activate the general transcription machinery located at promoters. It has thus been described as a bridge between these elements during initiation of transcription. Here, we review recent studies on Mediator relating to its structure, gene specificity and general requirement, roles in chromatin architecture as well as novel concepts involving phase separation and transcriptional bursting. We revisit the mechanism of action of Mediator and ultimately put forward models for its mode of action in gene activation.

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

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

Collective polymerase dynamics emerge from DNA supercoiling during transcription

All biological processes ultimately come from physical interactions. The mechanical properties of DNA play a critical role in transcription. RNA polymerase can over or under twist DNA (referred to as DNA supercoiling) when it moves along a gene, resulting in mechanical stresses in DNA that impact its own motion and that of other polymerases. For example, when enough supercoiling accumulates, an isolated polymerase halts, and transcription stops. DNA supercoiling can also mediate nonlocal interactions between polymerases that shape gene expression fluctuations. Here, we construct a comprehensive model of transcription that captures how RNA polymerase motion changes the degree of DNA supercoiling, which in turn feeds back into the rate at which polymerases are recruited and move along the DNA. Surprisingly, our model predicts that a group of three or more polymerases move together at a constant velocity and sustain their motion (forming what we call a polymeton), whereas one or two polymerases would have halted. We further show that accounting for the impact of DNA supercoiling on both RNA polymerase recruitment and velocity recapitulates empirical observations of gene expression fluctuations. Finally, we propose a mechanical toggle switch whereby interactions between genes are mediated by DNA twisting as opposed to proteins. Understanding the mechanical regulation of gene expression provides new insights into how endogenous genes can interact and informs the design of new forms of engineered interactions.

The Mediator complex as a master regulator of transcription by RNA polymerase II

The Mediator complex, which in humans is 1.4 MDa in size and includes 26 subunits, controls many aspects of RNA polymerase II (Pol II) function. Apart from its size, a defining feature of Mediator is its intrinsic disorder and conformational flexibility, which contributes to its ability to undergo phase separation and to interact with a myriad of regulatory factors. In this Review, we discuss Mediator structure and function, with emphasis on recent cryogenic electron microscopy data of the 4.0-MDa transcription preinitiation complex. We further discuss how Mediator and sequence-specific DNA-binding transcription factors enable enhancer-dependent regulation of Pol II function at distal gene promoters, through the formation of molecular condensates (or transcription hubs) and chromatin loops. Mediator regulation of Pol II reinitiation is also discussed, in the context of transcription bursting. We propose a working model for Mediator function that combines experimental results and theoretical considerations related to enhancer-promoter interactions, which reconciles contradictory data regarding whether enhancer-promoter communication is direct or indirect. We conclude with a discussion of Mediator's potential as a therapeutic target and of future research directions.

Technologies Enabling Single-Molecule Super-Resolution Imaging of mRNA

The transient nature of RNA has rendered it one of the more difficult biological targets for imaging. This difficulty stems both from the physical properties of RNA as well as the temporal constraints associated therewith. These concerns are further complicated by the difficulty in imaging endogenous RNA within a cell that has been transfected with a target sequence. These concerns, combined with traditional concerns associated with super-resolution light microscopy has made the imaging of this critical target difficult. Recent advances have provided researchers the tools to image endogenous RNA in live cells at both the cellular and single-molecule level. Here, we review techniques used for labeling and imaging RNA with special emphases on various labeling methods and a virtual 3D super-resolution imaging technique.

RNase HI Depletion Strongly Potentiates Cell Killing by Rifampicin in Mycobacteria

Multidrug-resistant (MDR) tuberculosis (TB) is defined by the resistance of Mycobacterium tuberculosis, the causative organism, to the first-line antibiotics rifampicin and isoniazid. Mitigating or reversing resistance to these drugs offers a means of preserving and extending their use in TB treatment. R-loops are RNA/DNA hybrids that are formed in the genome during transcription, and they can be lethal to the cell if not resolved. RNase HI is an enzyme that removes R-loops, and this activity is essential in M. tuberculosis: knockouts of rnhC, the gene encoding RNase HI, are nonviable. This essentiality makes it a candidate target for the development of new antibiotics. In the model organism Mycolicibacterium smegmatis, RNase HI activity is provided by two enzymes, RnhA and RnhC. We show that the partial depletion of RNase HI activity in M. smegmatis, by knocking out either of the genes encoding RnhA or RnhC, led to the accumulation of R-loops. The sensitivity of the knockout strains to the antibiotics moxifloxacin, streptomycin, and rifampicin was increased, the latter by a striking near 100-fold. We also show that R-loop accumulation accompanies partial transcriptional inhibition, suggesting a mechanistic basis for the synergy between RNase HI depletion and rifampicin. A model of how transcriptional inhibition can potentiate R-loop accumulation is presented. Finally, we identified four small molecules that inhibit recombinant RnhC activity and that also potentiated rifampicin activity in whole-cell assays against M. tuberculosis, supporting an on-target mode of action and providing the first step in developing a new class of antimycobacterial drug.

A transcriptional cycling model recapitulates chromatin-dependent features of noisy inducible transcription

Activation of gene expression in response to environmental cues results in substantial phenotypic heterogeneity between cells that can impact a wide range of outcomes including differentiation, viral activation, and drug resistance. An important source of gene expression noise is transcriptional bursting, or the process by which transcripts are produced during infrequent bursts of promoter activity. Chromatin accessibility impacts transcriptional bursting by regulating the assembly of transcription factor and polymerase complexes on promoters, suggesting that the effect of an activating signal on transcriptional noise will depend on the initial chromatin state at the promoter. To explore this possibility, we simulated transcriptional activation using a transcriptional cycling model with three promoter states that represent chromatin remodeling, polymerase binding and pause release. We initiated this model over a large parameter range representing target genes with different chromatin environments, and found that, upon increasing the polymerase pause release rate to activate transcription, changes in gene expression noise varied significantly across initial promoter states. This model captured phenotypic differences in activation of latent HIV viruses integrated at different chromatin locations and mediated by the transcription factor NF-κB. Activating transcription in the model via increasing one or more of the transcript production rates, as occurs following NF-κB activation, reproduced experimentally measured transcript distributions for four different latent HIV viruses, as well as the bimodal pattern of HIV protein expression that leads to a subset of reactivated virus. Importantly, the parameter ‘activation path’ differentially affected gene expression noise, and ultimately viral activation, in line with experimental observations. This work demonstrates how upstream signaling pathways can be connected to biological processes that underlie transcriptional bursting, resulting in target gene-specific noise profiles following stimulation of a single upstream pathway.

Many genes are transcribed in infrequent bursts of mRNA production through a process called transcriptional bursting, which contributes to variability in responses between cells. Heterogeneity in cell responses can have important biological impacts, such as whether a cell supports viral replication or responds to a drug, and thus there is an effort to describe this process with mathematical models to predict biological outcomes. Previous models described bursting as a transition between an “OFF” state or an “ON” state, an elegant and simple mathematical representation of complex molecular mechanisms, but one which failed to capture how upstream activation signals affected bursting. To address this, we added an additional promoter state to better reflect biological mechanisms underlying bursting. By fitting this model to variable activation of quiescent HIV infections in T cells, we showed that our model more accurately described viral expression variability across cells in response to an upstream stimulus. Our work highlights how mathematical models can be further developed to understand complex biological mechanisms and suggests ways to connect transcriptional bursting to upstream activation pathways.

Label and quantify mRNA molecules in live cell experiments using SunRISER and dNEMO

Visualization of mRNA molecules in single cells has revealed their core mechanisms as well as sources of cell-to-cell and spatiotemporal heterogeneity. Here, we describe a protocol to label, visualize, and quantify mRNA molecules by time-lapse imaging with the capability of resolving mRNA molecules over durations of hours to days. We provide links to mRNA-labeling plasmids as well as free software for a semi-automated image analysis pipeline.

For complete details on the use and execution of this protocol, please refer to Guo and Lee (2022) and Kowalczyk et al. (2021).

•

SunRISER is a toolbox optimized for single mRNA live-cell imaging experiments

•

Rapid detection of labeled single mRNA molecules in fluorescence microscopy images

•

Automated single-cell segmentation and mRNA quantification

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Visualization, Quantification, and Modeling of Endogenous RNA Polymerase II Phosphorylation at a Single-copy Gene in Living Cells

In eukaryotic cells, RNA Polymerase II (RNAP2) is the enzyme in charge of transcribing mRNA from DNA. RNAP2 possesses an extended carboxy-terminal domain (CTD) that gets dynamically phosphorylated as RNAP2 progresses through the transcription cycle, therefore regulating each step of transcription from recruitment to termination. Although RNAP2 residue-specific phosphorylation has been characterized in fixed cells by immunoprecipitation-based assays, or in live cells by using tandem gene arrays, these assays can mask heterogeneity and limit temporal and spatial resolution. Our protocol employs multi-colored complementary fluorescent antibody-based (Fab) probes to specifically detect the CTD of the RNAP2 (CTD-RNAP2), and its phosphorylated form at the serine 5 residue (Ser5ph-RNAP2) at a single-copy HIV-1 reporter gene. Together with high-resolution fluorescence microscopy, single-molecule tracking analysis, and rigorous computational modeling, our system allows us to visualize, quantify, and predict endogenous RNAP2 phosphorylation dynamics and mRNA synthesis at a single-copy gene, in living cells, and throughout the transcription cycle. Graphical abstract: Schematic of the steps for visualizing, quantifying, and predicting RNAP2 phosphorylation at a single-copy gene.

Activation domains can decouple the mean and noise of gene expression

Regulatory mechanisms set a gene's average level of expression, but a gene's expression constantly fluctuates around that average. These stochastic fluctuations, or expression noise, play a role in cell-fate transitions, bet hedging in microbes, and the development of chemotherapeutic resistance in cancer. An outstanding question is what regulatory mechanisms contribute to noise. Here, we demonstrate that, for a fixed mean level of expression, strong activation domains (ADs) at low abundance produce high expression noise, while weak ADs at high abundance generate lower expression noise. We conclude that differences in noise can be explained by the interplay between a TF's nuclear concentration and the strength of its AD's effect on mean expression, without invoking differences between classes of ADs. These results raise the possibility of engineering gene expression noise independently of mean levels in synthetic biology contexts and provide a potential mechanism for natural selection to tune the noisiness of gene expression.

Long-term imaging of individual mRNA molecules in living cells

Single-cell imaging of individual mRNAs has revealed core mechanisms of the central dogma. However, most approaches require cell fixation or have limited sensitivity for live-cell applications. Here, we describe SunRISER (SunTag-based reporter for imaging signal-enriched mRNA), a computationally and experimentally optimized approach for unambiguous detection of single mRNA molecules in living cells. When viewed by epifluorescence microscopy, SunRISER-labeled mRNAs show strong signal to background and resistance to photobleaching, which together enable long-term mRNA imaging studies. SunRISER variants, using 8× and 10× stem-loop arrays, demonstrate effective mRNA detection while significantly reducing alterations to target mRNA sequences. We characterize SunRISER to observe mRNA inheritance during mitosis and find that stressors enhance diversity among post-mitotic sister cells. Taken together, SunRISER enables a glimpse into living cells to observe aspects of the central dogma and the role of mRNAs in rare and dynamical trafficking events.

Illuminating RNA biology through imaging

RNA processing plays a central role in accurately transmitting genetic information into functional RNA and protein regulators. To fully appreciate the RNA life-cycle, tools to observe RNA with high spatial and temporal resolution are critical. Here we review recent advances in RNA imaging and highlight how they will propel the field of RNA biology. We discuss current trends in RNA imaging and their potential to elucidate unanswered questions in RNA biology.

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.

Bromodomains regulate dynamic targeting of the PBAF chromatin-remodeling complex to chromatin hubs

Chromatin remodelers actively target arrays of acetylated nucleosomes at select enhancers and promoters to facilitate or shut down the repeated recruitment of RNA polymerase II during transcriptional bursting. It is poorly understood how chromatin remodelers such as PBAF dynamically target different chromatin states inside a live cell. Our live-cell single-molecule fluorescence microscopy study reveals chromatin hubs throughout the nucleus where PBAF rapidly cycles on and off the genome. Deletion of PBAF's bromodomains impairs targeting and stable engagement of chromatin in hubs. Dual color imaging reveals that PBAF targets both euchromatic and heterochromatic hubs with distinct genome-binding kinetic profiles that mimic chromatin stability. Removal of PBAF's bromodomains stabilizes H3.3 binding within chromatin, indicating that bromodomains may play a direct role in remodeling of the nucleosome. Our data suggests that PBAF's dynamic bromodomain-mediated engagement of a nucleosome may reflect the chromatin-remodeling potential of differentially bound chromatin states.

Walking a tightrope: The complex balancing act of R-loops in genome stability

Although transcription is an essential cellular process, it is paradoxically also a well-recognized cause of genomic instability. R-loops, non-B DNA structures formed when nascent RNA hybridizes to DNA to displace the non-template strand as single-stranded DNA (ssDNA), are partially responsible for this instability. Yet, recent work has begun to elucidate regulatory roles for R-loops in maintaining the genome. In this review, we discuss the cellular contexts in which R-loops contribute to genomic instability, particularly during DNA replication and double-strand break (DSB) repair. We also summarize the evidence that R-loops participate as an intermediate during repair and may influence pathway choice to preserve genomic integrity. Finally, we discuss the immunogenic potential of R-loops and highlight their links to disease should they become pathogenic.