Noise-mean relationship in mutated promoters

Epistatic impacts of cis- and trans-regulatory mutations on the distribution of mutational effects for gene expression in Saccharomyces cerevisiae

Epistasis can influence evolution by causing the distribution of phenotypic effects for new mutations to vary among genotypes. Here, we investigate how epistatic interactions between new mutations and an existing regulatory mutation might impact the evolution of gene expression using Saccharomyces cerevisiae. We do so by estimating the distribution of mutational effects for expression of a fluorescent reporter protein driven by the S. cerevisiae TDH3 promoter in a reference strain as well as in eight mutant strains. Each of the mutant strains differed from the reference strain by a single mutation affecting expression of the focal gene. We found that half of these regulatory mutations changed the variance and/or skewness of the distribution of mutational effects. A change in variance indicates a change in mutational robustness, and we found that one initial regulatory mutation increased mutational robustness while another decreased it. A change in skewness indicates a change in the relative frequency and/or effect size of mutations increasing or decreasing expression, and we found that the initial regulatory mutation in four strains had such an effect. Strikingly, in all four of these cases, the change in skewness increased the likelihood that new mutations would at least partially compensate for the effects of the initial regulatory mutation. If this form of epistatic impact on the distribution of mutational effects is common, it could provide a neutral mechanism reducing the divergence of gene expression and help explain the prevalence of alleles with compensatory effects in natural populations of S. cerevisiae.

The impact of phenotypic heterogeneity on fungal pathogenicity and drug resistance

Phenotypic heterogeneity in genetically clonal populations facilitates cellular adaptation to adverse environmental conditions while enabling a return to the basal physiological state. It also plays a crucial role in pathogenicity and the acquisition of drug resistance in unicellular organisms and cancer cells, yet the exact contributing factors remain elusive. In this review, we outline the current state of understanding concerning the contribution of phenotypic heterogeneity to fungal pathogenesis and antifungal drug resistance.

Exact burst-size distributions for gene-expression models with complex promoter structure

In prokaryotic and eukaryotic cells, most genes are transcribed in a bursty fashion on one hand and complex gene regulations may lead to complex promoter structure on the other hand. This raises an unsolved issue: how does promoter structure shape transcriptional bursting kinetics characterized by burst size and frequency? Here we analyze stochastic models of gene transcription, which consider complex regulatory mechanisms. Notably, we develop an efficient method to derive exact burst-size distributions. The analytical results show that if the promoter of a gene contains only one active state, the burst size indeed follows a geometric distribution, in agreement with the previous result derived under certain limiting conditions. However, if it contains a multitude of active states, the burst size in general obeys a non-geometric distribution, which is a linearly weighted sum of geometric distributions. This superposition principle reveals the essential feature of bursting kinetics in complex cases of transcriptional regulation although it seems that there has been no direct experimental confirmation. The derived burst-size distributions not only highlight the importance of promoter structure in regulating bursting kinetics, but can be also used in the exact inference of this kinetics based on experimental data.

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

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

Transcriptional noise, gene activation, and roles of SAGA and Mediator Tail measured using nucleotide recoding single-cell RNA-seq

We describe a time-resolved nascent single-cell RNA sequencing (RNA-seq) approach that measures gene-specific transcriptional noise and the fraction of active genes in S. cerevisiae. Most genes are expressed with near-constitutive behavior, while a subset of genes show high mRNA variance suggestive of transcription bursting. Transcriptional noise is highest in the cofactor/coactivator-redundant (CR) gene class (dependent on both SAGA and TFIID) and strongest in TATA-containing CR genes. Using this approach, we also find that histone gene transcription switches from a low-level, low-noise constitutive mode during M and M/G1 to an activated state in S phase that shows both an increase in the fraction of active promoters and a switch to a noisy and bursty transcription mode. Rapid depletion of cofactors SAGA and MED Tail indicates that both factors play an important role in stimulating the fraction of active promoters at CR genes, with a more modest role in transcriptional noise.

Living in a noisy world-origins of gene expression noise and its impact on cellular decision-making

The expression level of a gene can vary between genetically identical cells under the same environmental condition-a phenomenon referred to as gene expression noise. Several studies have now elucidated a central role of transcription factors in the generation of expression noise. Transcription factors, as the key components of gene regulatory networks, drive many important cellular decisions in response to cellular and environmental signals. Therefore, a very relevant question is how expression noise impacts gene regulation and influences cellular decision-making. In this Review, we summarize the current understanding of the molecular origins of expression noise, highlighting the role of transcription factors in this process, and discuss the ways in which noise can influence cellular decision-making. As advances in single-cell technologies open new avenues for studying expression noise as well as gene regulatory circuits, a better understanding of the influence of noise on cellular decisions will have important implications for many biological processes.

Analysis of developmental gene expression using smFISH and in silico staging of C. elegans embryos

Regulation of transcription during embryogenesis is key to development and differentiation. To study transcript expression throughout Caenorhabditis elegans embryogenesis at single-molecule resolution, we developed a high-throughput single-molecule fluorescence in situ hybridization (smFISH) method that relies on computational methods to developmentally stage embryos and quantify individual mRNA molecules in single embryos. We applied our system to sdc-2, a zygotically transcribed gene essential for hermaphrodite development and dosage compensation. We found that sdc-2 is rapidly activated during early embryogenesis by increasing both the number of mRNAs produced per transcription site and the frequency of sites engaged in transcription. Knockdown of sdc-2 and dpy-27, a subunit of the dosage compensation complex (DCC), increased the number of active transcription sites for the X chromosomal gene dpy-23 but not the autosomal gene mdh-1, suggesting that the DCC reduces the frequency of dpy-23 transcription. The temporal resolution from in silico staging of embryos showed that the deletion of a single DCC recruitment element near the dpy-23 gene causes higher dpy-23 mRNA expression after the start of dosage compensation, which could not be resolved using mRNAseq from mixed-stage embryos. In summary, we have established a computational approach to quantify temporal regulation of transcription throughout C. elegans embryogenesis and demonstrated its potential to provide new insights into developmental gene regulation.

A large sequenced mutant library - valuable reverse genetic resource that covers 98% of sorghum genes

Mutant populations are crucial for functional genomics and discovering novel traits for crop breeding. Sorghum, a drought and heat-tolerant C4 species, requires a vast, large-scale, annotated, and sequenced mutant resource to enhance crop improvement through functional genomics research. Here, we report a sorghum large-scale sequenced mutant population with 9.5 million ethyl methane sulfonate (EMS)-induced mutations that covered 98% of sorghum's annotated genes using inbred line BTx623. Remarkably, a total of 610 320 mutations within the promoter and enhancer regions of 18 000 and 11 790 genes, respectively, can be leveraged for novel research of cis-regulatory elements. A comparison of the distribution of mutations in the large-scale mutant library and sorghum association panel (SAP) provides insights into the influence of selection. EMS-induced mutations appeared to be random across different regions of the genome without significant enrichment in different sections of a gene, including the 5' UTR, gene body, and 3'-UTR. In contrast, there were low variation density in the coding and UTR regions in the SAP. Based on the Ka /Ks value, the mutant library (~1) experienced little selection, unlike the SAP (0.40), which has been strongly selected through breeding. All mutation data are publicly searchable through SorbMutDB (https://www.depts.ttu.edu/igcast/sorbmutdb.php) and SorghumBase (https://sorghumbase.org/). This current large-scale sequence-indexed sorghum mutant population is a crucial resource that enriched the sorghum gene pool with novel diversity and a highly valuable tool for the Poaceae family, that will advance plant biology research and crop breeding.

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.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian

Gene expression inherently gives rise to stochastic variation ("noise") in the production of gene products. Minimizing noise is crucial for ensuring reliable cellular functions. However, noise cannot be suppressed below a certain intrinsic limit. For constitutively expressed genes, this limit is typically assumed to be Poissonian noise, wherein the variance in mRNA numbers is equal to their mean. Here, we demonstrate that several cell division genes in fission yeast exhibit mRNA variances significantly below this limit. The reduced variance can be explained by a gene expression model incorporating multiple transcription and mRNA degradation steps. Notably, in this sub-Poissonian regime, distinct from Poissonian or super-Poissonian regimes, cytoplasmic noise is effectively suppressed through a higher mRNA export rate. Our findings redefine the lower limit of eukaryotic gene expression noise and uncover molecular requirements for achieving ultralow noise, which is expected to be important for vital cellular functions.

Contributions of mutation and selection to regulatory variation: lessons from the Saccharomyces cerevisiae TDH3 gene

Heritable variation in gene expression is common within and among species and contributes to phenotypic diversity. Mutations affecting either cis- or trans-regulatory sequences controlling gene expression give rise to variation in gene expression, and natural selection acting on this variation causes some regulatory variants to persist in a population for longer than others. To understand how mutation and selection interact to produce the patterns of regulatory variation we see within and among species, my colleagues and I have been systematically determining the effects of new mutations on expression of the TDH3 gene in Saccharomyces cerevisiae and comparing them to the effects of polymorphisms segregating within this species. We have also investigated the molecular mechanisms by which regulatory variants act. Over the past decade, this work has revealed properties of cis- and trans-regulatory mutations including their relative frequency, effects, dominance, pleiotropy and fitness consequences. Comparing these mutational effects to the effects of polymorphisms in natural populations, we have inferred selection acting on expression level, expression noise and phenotypic plasticity. Here, I summarize this body of work and synthesize its findings to make inferences not readily discernible from the individual studies alone. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Broad compatibility between yeast UAS elements and core promoters and identification of promoter elements that determine cofactor specificity

Three classes of yeast protein-coding genes are distinguished by their dependence on the transcription cofactors TFIID, SAGA, and Mediator (MED) Tail, but whether this dependence is determined by the core promoter, upstream activating sequences (UASs), or other gene features is unclear. Also unclear is whether UASs can broadly activate transcription from the different promoter classes. Here, we measure transcription and cofactor specificity for thousands of UAS-core promoter combinations and find that most UASs broadly activate promoters regardless of regulatory class, while few display strong promoter specificity. However, matching UASs and promoters from the same gene class is generally important for optimal expression. We find that sensitivity to rapid depletion of MED Tail or SAGA is dependent on the identity of both UAS and core promoter, while dependence on TFIID localizes to only the promoter. Finally, our results suggest the role of TATA and TATA-like promoter sequences in MED Tail function.

Complementary strategies for directing in vivo transcription factor binding through DNA binding domains and intrinsically disordered regions

DNA binding domains (DBDs) of transcription factors (TFs) recognize DNA sequence motifs that are highly abundant in genomes. Within cells, TFs bind a subset of motif-containing sites as directed by either their DBDs or DBD-external (nonDBD) sequences. To define the relative roles of DBDs and nonDBDs in directing binding preferences, we compared the genome-wide binding of 48 (∼30%) budding yeast TFs with their DBD-only, nonDBD-truncated, and nonDBD-only mutants. With a few exceptions, binding locations differed between DBDs and TFs, resulting from the cumulative action of multiple determinants mapped mostly to disordered nonDBD regions. Furthermore, TFs' preferences for promoters of the fuzzy nucleosome architecture were lost in DBD-only mutants, whose binding spread across promoters, implicating nonDBDs' preferences in this hallmark of budding yeast regulatory design. We conclude that DBDs and nonDBDs employ complementary DNA-targeting strategies, whose balance defines TF binding specificity along genomes.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian

Stochastic variation in gene products ("noise") is an inescapable by-product of gene expression. Noise must be minimized to allow for the reliable execution of cellular functions. However, noise cannot be suppressed beyond an intrinsic lower limit. For constitutively expressed genes, this limit is believed to be Poissonian, meaning that the variance in mRNA numbers cannot be lower than their mean. Here, we show that several cell division genes in fission yeast have mRNA variances significantly below this limit, which cannot be explained by the classical gene expression model for low-noise genes. Our analysis reveals that multiple steps in both transcription and mRNA degradation are essential to explain this sub-Poissonian variance. The sub-Poissonian regime differs qualitatively from previously characterized noise regimes, a hallmark being that cytoplasmic noise is reduced when the mRNA export rate increases. Our study re-defines the lower limit of eukaryotic gene expression noise and identifies molecular requirements for ultra-low noise which are expected to support essential cell functions.

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.

Transcription factor binding process is the primary driver of noise in gene expression

Noise in expression of individual genes gives rise to variations in activity of cellular pathways and generates heterogeneity in cellular phenotypes. Phenotypic heterogeneity has important implications for antibiotic persistence, mutation penetrance, cancer growth and therapy resistance. Specific molecular features such as the presence of the TATA box sequence and the promoter nucleosome occupancy have been associated with noise. However, the relative importance of these features in noise regulation is unclear and how well these features can predict noise has not yet been assessed. Here through an integrated statistical model of gene expression noise in yeast we found that the number of regulating transcription factors (TFs) of a gene was a key predictor of noise, whereas presence of the TATA box and the promoter nucleosome occupancy had poor predictive power. With an increase in the number of regulatory TFs, there was a rise in the number of cooperatively binding TFs. In addition, an increased number of regulatory TFs meant more overlaps in TF binding sites, resulting in competition between TFs for binding to the same region of the promoter. Through modeling of TF binding to promoter and application of stochastic simulations, we demonstrated that competition and cooperation among TFs could increase noise. Thus, our work uncovers a process of noise regulation that arises out of the dynamics of gene regulation and is not dependent on any specific transcription factor or specific promoter sequence.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

Gene regulation in Escherichia coli is commonly selected for both high plasticity and low noise

Bacteria often respond to dynamically changing environments by regulating gene expression. Despite this regulation being critically important for growth and survival, little is known about how selection shapes gene regulation in natural populations. To better understand the role natural selection plays in shaping bacterial gene regulation, here we compare differences in the regulatory behaviour of naturally segregating promoter variants from Escherichia coli (which have been subject to natural selection) to randomly mutated promoter variants (which have never been exposed to natural selection). We quantify gene expression phenotypes (expression level, plasticity and noise) for hundreds of promoter variants across multiple environments and show that segregating promoter variants are enriched for mutations with minimal effects on expression level. In many promoters, we infer that there is strong selection to maintain high levels of plasticity, and direct selection to decrease or increase cell-to-cell variability in expression. Taken together, these results expand our knowledge of how gene regulation is affected by natural selection and highlight the power of comparing naturally segregating polymorphisms to de novo random mutations to quantify the action of selection.

Kinetic characteristics of transcriptional bursting in a complex gene model with cyclic promoter structure

While transcription often occurs in a bursty manner, various possible regulations can lead to complex promoter patterns such as promoter cycles, giving rise to an important question: How do promoter kinetics shape transcriptional bursting kinetics? Here we introduce and analyze a general model of the promoter cycle consisting of multi-OFF states and multi-ON states, focusing on the effects of multi-ON mechanisms on transcriptional bursting kinetics. The derived analytical results indicate that burst size follows a mixed geometric distribution rather than a single geometric distribution assumed in previous studies, and ON and OFF times obey their own mixed exponential distributions. In addition, we find that the multi-ON mechanism can lead to bimodal burst-size distribution, antagonistic timing of ON and OFF, and diverse burst frequencies, each further contributing to cell-to-cell variability in the mRNA expression level. These results not only reveal essential features of transcriptional bursting kinetics patterns shaped by multi-state mechanisms but also can be used to the inferences of transcriptional bursting kinetics and promoter structure based on experimental data.

Silent transcription intervals and translational bursting lead to diverse phenotypic switching

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

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

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

Dynamic modulation of enhancer responsiveness by core promoter elements in living Drosophila embryos

Regulatory interactions between enhancers and core promoters are fundamental for the temporal and spatial specificity of gene expression in development. The central role of core promoters is to initiate productive transcription in response to enhancer's activation cues. However, it has not been systematically assessed how individual core promoter elements affect the induction of transcriptional bursting by enhancers. Here, we provide evidence that each core promoter element differentially modulates functional parameters of transcriptional bursting in developing Drosophila embryos. Quantitative live imaging analysis revealed that the timing and the continuity of burst induction are common regulatory steps on which core promoter elements impact. We further show that the upstream TATA also affects the burst amplitude. On the other hand, Inr, MTE and DPE mainly contribute to the regulation of the burst frequency. Genome editing analysis of the pair-rule gene fushi tarazu revealed that the endogenous TATA and DPE are both essential for its correct expression and function during the establishment of body segments in early embryos. We suggest that core promoter elements serve as a key regulatory module in converting enhancer activity into transcription dynamics during animal development.

Independent control of mean and noise by convolution of gene expression distributions

Gene expression noise can reduce cellular fitness or facilitate processes such as alternative metabolism, antibiotic resistance, and differentiation. Unfortunately, efforts to study the impacts of noise have been hampered by a scaling relationship between noise and expression level from individual promoters. Here, we use theory to demonstrate that mean and noise can be controlled independently by expressing two copies of a gene from separate inducible promoters in the same cell. We engineer low and high noise inducible promoters to validate this result in Escherichia coli, and develop a model that predicts the experimental distributions. Finally, we use our method to reveal that the response of a promoter to a repressor is less sensitive with higher repressor noise and explain this result using a law from probability theory. Our approach can be applied to investigate the effects of noise on diverse biological pathways or program cellular heterogeneity for synthetic biology applications.

The Stochastic Genome and Its Role in Gene Expression

Mammalian genomes have distinct levels of spatial organization and structure that have been hypothesized to play important roles in transcription regulation. Although much has been learned about these architectural features with ensemble techniques, single-cell studies are showing a new universal trend: Genomes are stochastic and dynamic at every level of organization. Stochastic gene expression, on the other hand, has been studied for years. In this review, we probe whether there is a causative link between the two phenomena. We specifically discuss the functionality of chromatin state, topologically associating domains (TADs), and enhancer biology in light of their stochastic nature and their specific roles in stochastic gene expression. We highlight persistent fundamental questions in this area of research.

Correlation between external regulators governs the mean-noise relationship in stochastic gene expression

Gene transcription in single cells is inherently a probabilistic process. The relationship between variance ($ \sigma^{2} $) and mean expression ($ \mu $) is of paramount importance for investigations into the evolutionary origins and consequences of noise in gene expression. It is often formulated as $ \log \left({{{\sigma}^{2}}}/{{{\mu}^{2}}}\; \right) = \beta\log\mu+\log\alpha $, where $ \beta $ is a key parameter since its sign determines the qualitative dependence of noise on mean. We reveal that the sign of $ \beta $ is controlled completely by external regulation, but independent of promoter structure. Specifically, it is negative if regulators as stochastic variables are independent but positive if they are correlated. The essential mechanism revealed here can well interpret diverse experimental phenomena underlying expression noise. Our results imply that external regulation rather than promoter sequence governs the mean-noise relationship.

Exploring dissipative sources of non-Markovian biochemical reaction systems

Many biological processes including important intracellular processes are governed by biochemical reaction networks. Usually, these reaction systems operate far from thermodynamic equilibrium, implying free-energy dissipation. On the other hand, single reaction events happen often in a memory manner, leading to non-Markovian kinetics. A question then arises: how do we calculate free-energy dissipation (defined as the entropy production rate) in this physically real case? We derive an analytical formula for calculating the energy consumption of a general reaction system with molecular memory characterized by nonexponential waiting-time distributions. It shows that this dissipation is composed of two parts: one from broken detailed balance of an equivalent Markovian system with the same topology and substrates, and the other from the direction-time dependence of waiting-time distributions. But, if the system is in a detailed balance and the waiting-time distribution is direction-time independent, there is no energy dissipation even in the non-Markovian case. These general results provide insights into the physical mechanisms underlying nonequilibrium processes. A continuous-time random-walk model and a generalized model of stochastic gene expression are chosen to clearly show dissipative sources and the relationship between energy dissipation and molecular memory.

Scaling of intrinsic noise in an autocratic reaction network

Biochemical reactions in living cells often produce stochastic trajectories due to the fluctuations of the finite number of the macromolecular species present inside the cell. A significant number of computational and theoretical studies have previously investigated stochasticity in small regulatory networks to understand its origin and regulation. At the systems level regulatory networks have been determined to be hierarchical resembling social networks. In order to determine the stochasticity in networks with hierarchical architecture, here we computationally investigated intrinsic noise in an autocratic reaction network in which only the upstream regulators regulate the downstream regulators. We studied the effects of the qualitative and quantitative nature of regulatory interactions on the stochasticity in the network. We established an unconventional scaling of noise with average abundance in which the noise passes through a minimum indicating that the network can be noisy both in the low and high abundance regimes. We determined that the bursty kinetics of the trajectories are responsible for such scaling. The scaling of noise remains intact for a mixed network that includes democratic subnetworks within the autocratic network.

Interplay between genetic, epigenetic, and gene expression variability: Considering complexity in evolvability

Genetic variability, epigenetic variability, and gene expression variability (noise) are generally considered independently in their relationship with phenotypic variation. However, they appear to be intrinsically interconnected and influence it in combination. The study of the interplay between genetic and epigenetic variability has the longest history. This article rather considers the introduction of gene expression variability in its relationships with the two others and reviews for the first time experimental evidences over the four relationships connected to gene expression noise. They show how introducing this third source of variability complicates the way of thinking evolvability and the emergence of biological novelty. Finally, cancer cells are proposed to be an ideal model to decipher the dynamic interplay between genetic, epigenetic, and gene expression variability when one of them is either experimentally increased or therapeutically targeted. This interplay is also discussed in an evolutionary perspective in the context of cancer cell drug resistance.

TATA and paused promoters active in differentiated tissues have distinct expression characteristics

Core promoter types differ in the extent to which RNA polymerase II (Pol II) pauses after initiation, but how this affects their tissue-specific gene expression characteristics is not well understood. While promoters with Pol II pausing elements are active throughout development, TATA promoters are highly active in differentiated tissues. We therefore used a genomics approach on late-stage Drosophila embryos to analyze the properties of promoter types. Using tissue-specific Pol II ChIP-seq, we found that paused promoters have high levels of paused Pol II throughout the embryo, even in tissues where the gene is not expressed, while TATA promoters only show Pol II occupancy when the gene is active. The promoter types are associated with different chromatin accessibility in ATAC-seq data and have different expression characteristics in single-cell RNA-seq data. The two promoter types may therefore be optimized for different properties: paused promoters show more consistent expression when active, while TATA promoters have lower background expression when inactive. We propose that tissue-specific genes have evolved to use two different strategies for their differential expression across tissues.

Explicit effect of stochastic reaction delay on gene expression

Apart from intrinsic stochastic variability, gene expression also involves stochastic reaction delay arising from heterogeneity and fluctuation processes, which can affect the efficiency of reactants (e.g., mRNA or protein) in exploring their environments. In contrast to the former that has been extensively investigated, the impact of the latter on gene expression remains not fully understood. Here, we analyze a non-Markovian model of bursty gene expression with general delay distribution. We analytically find that the effect of stochastic reaction delay is equivalent to the introduction of negative feedback, and stationary protein distribution only depends on the mean of the delay and is independent of its distribution. We numerically show that the stochastic reaction delay always slightly amplifies the mean protein level but remarkably reduces the protein noise (quantified by the ratio of the variance over the squared average). Our analysis indicates that stochastic reaction delay is an important factor affecting gene expression.

Identifying chromatin features that regulate gene expression distribution

Gene expression variability, differences in the number of mRNA per cell across a population of cells, is ubiquitous across diverse organisms with broad impacts on cellular phenotypes. The role of chromatin in regulating average gene expression has been extensively studied. However, what aspects of the chromatin contribute to gene expression variability is still underexplored. Here we addressed this problem by leveraging chromatin diversity and using a systematic investigation of randomly integrated expression reporters to identify what aspects of chromatin microenvironment contribute to gene expression variability. Using DNA barcoding and split-pool decoding, we created a large library of isogenic reporter clones and identified reporter integration sites in a massive and parallel manner. By mapping our measurements of reporter expression at different genomic loci with multiple epigenetic profiles including the enrichment of transcription factors and the distance to different chromatin states, we identified new factors that impact the regulation of gene expression distributions.

Recent advances in understanding RNA polymerase II structure and function

More than 50 years after the identification of RNA polymerase II, the enzyme responsible for the transcription of most eukaryotic genes, studies have continued to reveal fresh aspects of its structure and regulation. New technologies, coupled with years of development of a vast catalog of RNA polymerase II accessory proteins and activities, have led to new revelations about the transcription process. The maturation of cryo-electron microscopy as a tool for unraveling the detailed structure of large molecular machines has provided numerous structures of the enzyme and its accessory factors. Advances in biophysical methods have enabled the observation of a single polymerase’s behavior, distinct from work on aggregate population averages. Other recent work has revealed new properties and activities of the general initiation factors that RNA polymerase II employs to accurately initiate transcription, as well as chromatin proteins that control RNA polymerase II’s firing frequency, and elongation factors that facilitate the enzyme’s departure from the promoter and which control sequential steps and obstacles that must be navigated by elongating RNA polymerase II. There has also been a growing appreciation of the physical properties conferred upon many of these proteins by regions of each polypeptide that are of low primary sequence complexity and that are often intrinsically disordered. This peculiar feature of a surprisingly large number of proteins enables a disordered region of the protein to morph into a stable structure and creates an opportunity for pathway participants to dynamically partition into subcompartments of the nucleus. These subcompartments host designated portions of the chemical reactions that lead to mRNA synthesis. This article highlights a selection of recent findings that reveal some of the resolved workings of RNA polymerase II and its ensemble of supporting factors.

Capturing and Understanding the Dynamics and Heterogeneity of Gene Expression in the Living Cell

The regulation of gene expression is a fundamental process enabling cells to respond to internal and external stimuli or to execute developmental programs. Changes in gene expression are highly dynamic and depend on many intrinsic and extrinsic factors. In this review, we highlight the dynamic nature of transient gene expression changes to better understand cell physiology and development in general. We will start by comparing recent in vivo procedures to capture gene expression in real time. Intrinsic factors modulating gene expression dynamics will then be discussed, focusing on chromatin modifications. Furthermore, we will dissect how cell physiology or age impacts on dynamic gene regulation and especially discuss molecular insights into acquired transcriptional memory. Finally, this review will give an update on the mechanisms of heterogeneous gene expression among genetically identical individual cells. We will mainly focus on state-of-the-art developments in the yeast model but also cover higher eukaryotic systems.

Regulation Between HSF1 Isoforms and HSPs Contributes to the Variation in Thermal Tolerance Between Two Oyster Congeners

Heat shock transcription factor 1 (HSF1) plays an important role in regulating heat shock, which can activate heat shock proteins (HSPs). HSPs can protect organisms from thermal stress. Oysters in the intertidal zone can tolerate thermal stress. The Pacific oyster (Crassostrea gigas gigas) and Fujian oyster (C. gigas angulata)—allopatric subspecies with distinct thermal tolerances—make good study specimens for analyzing and comparing thermal stress regulation. We cloned and compared HSF1 isoforms, which is highly expressed under heat shock conditions in the two subspecies. The results revealed that two isoforms (HSF1a and HSF1d) respond to heat shock in both Pacific and Fujian oysters, and different heat shock conditions led to various combinations of isoforms. Subcellular localization showed that isoforms gathered in the nucleus when exposed to heat shock. The co-immunoprecipitation revealed that HSF1d can be a dimer. In addition, we selected HSPs that are expressed under the heat shock response, according to the RNA-seq and proteomic analyses. For the HSPs, we analyzed the coding part and the promoter sequences. The result showed that the domains of HSPs are conserved in two subspecies, but the promoters are significantly different. The Dual-Luciferase assay showed that the induced expression isoform HSF1d had the highest activity in C. gigas gigas, while the constitutively-expressed HSF1a was most active in C. gigas angulata. In addition, variation in the level of HSP promoters appeared to be correlated with gene expression. We argue that this gene is regulated based on the different expression levels between the two subspecies’ responses to heat shock. In summary, various stress conditions can yield different HSF1 isoforms and respond to heat shock in both oyster subspecies. Differences in how the isoforms and promoter are activated may contribute to their differential expressions. Overall, the results comparing C. gigas gigas and C. gigas angulata suggest that these isoforms have a regulatory relationship under heat shock, providing valuable information on the thermal tolerance mechanism in these commercially important oyster species.

Buffering noise: KAT2A modular contributions to stabilization of transcription and cell identity in cancer and development

KAT2A is a histone acetyltransferase recently identified as a vulnerability in at least some forms of Acute Myeloid Leukemia (AML). Its loss or inhibition prompts leukemia stem cells out of self-renewal and into differentiation with ultimate exhaustion of the leukemia pool. We have recently linked the Kat2a requirement in AML to control of transcriptional noise, reflecting an evolutionary-conserved role of Kat2a in promoting burst-like promoter activity and stabilizing gene expression. We suggest that through this role, Kat2a contributes to preservation of cell identity. KAT2A exerts its acetyltransferase activity in the context of two macromolecular complexes, Spt-Ada-Gcn5-Acetyltransferase (SAGA) and Ada-Two-A-Containing (ATAC), but the specific contribution of each complex to stabilization of gene expression is currently unknown. By reviewing specific gene targets and requirements of the two complexes in cancer and development, we suggest that SAGA regulates lineage-specific programs, and ATAC maintains biosynthetic activity through control of ribosomal protein and translation-associated genes, on which cells may be differentially dependent. While our data suggest that KAT2A-mediated regulation of transcriptional noise in AML may be exerted through ATAC, we discuss potential caveats and probe general vs. complex-specific contributions of KAT2A to transcriptional stability, with implications for control and perturbation of cell identity.

Does Gene Expression Noise Play a Functional Role in Plants?

Gene expression in individual cells can be surprisingly noisy. In unicellular organisms this noise can be functional; for example, by allowing a subfraction of the population to prepare for environmental stress. The role of gene expression noise in multicellular organisms has, however, remained unclear. In this review, we discuss how new techniques are revealing an unexpected level of variability in gene expression between and within genetically identical plants. We describe recent progress as well as speculate on the function of transcriptional noise as a mechanism for generating functional phenotypic diversity in plants.

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

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

SIR2 Expression Noise Can Generate Heterogeneity in Viability but Does Not Affect Cell-to-Cell Epigenetic Silencing of Subtelomeric URA3 in Yeast

Chromatin structure clearly modulates gene expression noise, but the reverse influence has never been investigated, namely how the cell-to-cell expression heterogeneity of chromatin modifiers may generate variable rates of epigenetic modification. Sir2 is a well-characterized histone deacetylase of the Sirtuin family. It strongly influences chromatin silencing, especially at telomeres, subtelomeres and rDNA. This ability to influence epigenetic landscapes makes it a good model to study the largely unexplored interplay between gene expression noise and other epigenetic processes leading to phenotypic diversification. Here, we addressed this question by investigating whether noise in the expression of SIR2 was associated with cell-to-cell heterogeneity in the frequency of epigenetic silencing at subtelomeres in Saccharomyces cerevisiae Using cell sorting to isolate subpopulations with various expression levels, we found that heterogeneity in the cellular concentration of Sir2 does not lead to heterogeneity in the epigenetic silencing of subtelomeric URA3 between these subpopulations. We also noticed that SIR2 expression noise can generate cell-to-cell variability in viability, with lower levels being associated with better viability. This work shows that SIR2 expression fluctuations are not sufficient to generate cell-to-cell heterogeneity in the epigenetic silencing of URA3 at subtelomeres in Saccharomyces cerevisiae but can strongly affect cellular viability.

Stochastic allelic expression as trigger for contractile imbalance in hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM), the most common inherited cardiac disease, is caused by several mostly heterozygous mutations in sarcomeric genes. Hallmarks of HCM are cardiomyocyte and myofibrillar disarray and hypertrophy and fibrosis of the septum and the left ventricle. To date, a pathomechanism common to all mutations remains elusive. We have proposed that contractile imbalance, an unequal force generation of neighboring cardiomyocytes, may contribute to development of HCM hallmarks. At the same calcium concentration, we found substantial differences in force generation between individual cardiomyocytes from HCM patients with mutations in β-MyHC (β-myosin heavy chain). Variability among cardiomyocytes was significantly larger in HCM patients as compared with donor controls. We assume that this heterogeneity in force generation among cardiomyocytes may lead to myocardial disarray and trigger hypertrophy and fibrosis. We provided evidence that burst-like transcription of the MYH7-gene, encoding for β-MyHC, is associated with unequal fractions of mutant per wild-type mRNA from cell to cell (cell-to-cell allelic imbalance). This will presumably lead to unequal fractions of mutant per wild-type protein from cell to cell which may underlie contractile imbalance. In this review, we discuss molecular mechanisms of burst-like transcription with regard to contractile imbalance and disease development in HCM.

Cell-to-cell expression dispersion of B-cell surface proteins is linked to genetic variants in humans

Variability in gene expression across a population of homogeneous cells is known to influence various biological processes. In model organisms, natural genetic variants were found that modify expression dispersion (variability at a fixed mean) but very few studies have detected such effects in humans. Here, we analyzed single-cell expression of four proteins (CD23, CD55, CD63 and CD86) across cell lines derived from individuals of the Yoruba population. Using data from over 30 million cells, we found substantial inter-individual variation of dispersion. We demonstrate, via de novo cell line generation and subcloning experiments, that this variation exceeds the variation associated with cellular immortalization. We detected a genetic association between the expression dispersion of CD63 and the rs971 SNP. Our results show that human DNA variants can have inherently-probabilistic effects on gene expression. Such subtle genetic effects may participate to phenotypic variation and disease outcome.

Triqueneaux, Burny, Symmons et al. show association between gene expression noise and genotypes, using single-cell expression of four proteins across human-derived lymphoblastoid cell lines. This study suggests that very subtle regulatory effects of human DNA variants may contribute to phenotypic variation and disease outcome.

Transcription in Living Cells: Molecular Mechanisms of Bursting

Transcription in several organisms from certain bacteria to humans has been observed to be stochastic in nature: toggling between active and inactive states. Periods of active nascent RNA synthesis known as bursts represent individual gene activation events in which multiple polymerases are initiated. Therefore, bursting is the single locus illustration of both gene activation and repression. Although transcriptional bursting was originally observed decades ago, only recently have technological advances enabled the field to begin elucidating gene regulation at the single-locus level. In this review, we focus on how biochemical, genomic, and single-cell data describe the regulatory steps of transcriptional bursts.

Permutational analysis of Saccharomyces cerevisiae regulatory elements

Gene expression in Saccharomyces cerevisiae is regulated at multiple levels. Genomic and epigenomic mapping of transcription factors and chromatin factors has led to the delineation of various modular regulatory elements—enhancers (upstream activating sequences), core promoters, 5′ untranslated regions (5′ UTRs) and transcription terminators/3′ untranslated regions (3′ UTRs). However, only a few of these elements have been tested in combinations with other elements and the functional interactions between the different modular regulatory elements remain under explored. We describe a simple and rapid approach to build a combinatorial library of regulatory elements and have used this library to study 26 different enhancers, core promoters, 5′ UTRs and transcription terminators/3′ UTRs to estimate the contribution of individual regulatory parts in gene expression. Our combinatorial analysis shows that while enhancers initiate gene expression, core promoters modulate the levels of enhancer-mediated expression and can positively or negatively affect expression from even the strongest enhancers. Principal component analysis (PCA) indicates that enhancer and promoter function can be explained by a single principal component while UTR function involves multiple functional components. The PCA also highlights outliers and suggest differences in mechanisms of regulation by individual elements. Our data also identify numerous regulatory cassettes composed of different individual regulatory elements that exhibit equivalent gene expression levels. These data thus provide a catalog of elements that could in future be used in the design of synthetic regulatory circuits.

Genome-wide kinetic properties of transcriptional bursting in mouse embryonic stem cells

Transcriptional bursting is the stochastic activation and inactivation of promoters, contributing to cell-to-cell heterogeneity in gene expression. However, the mechanism underlying the regulation of transcriptional bursting kinetics (burst size and frequency) in mammalian cells remains elusive. In this study, we performed single-cell RNA sequencing to analyze the intrinsic noise and mRNA levels for elucidating the transcriptional bursting kinetics in mouse embryonic stem cells. Informatics analyses and functional assays revealed that transcriptional bursting kinetics was regulated by a combination of promoter- and gene body-binding proteins, including the polycomb repressive complex 2 and transcription elongation factors. Furthermore, large-scale CRISPR-Cas9-based screening identified that the Akt/MAPK signaling pathway regulated bursting kinetics by modulating transcription elongation efficiency. These results uncovered the key molecular mechanisms underlying transcriptional bursting and cell-to-cell gene expression noise in mammalian cells.

The route to transcription initiation determines the mode of transcriptional bursting in E. coli

Transcription is fundamentally noisy, leading to significant heterogeneity across bacterial populations. Noise is often attributed to burstiness, but the underlying mechanisms and their dependence on the mode of promotor regulation remain unclear. Here, we measure E. coli single cell mRNA levels for two stress responses that depend on bacterial sigma factors with different mode of transcription initiation (σ⁷⁰ and σ⁵⁴). By fitting a stochastic model to the observed mRNA distributions, we show that the transition from low to high expression of the σ⁷⁰-controlled stress response is regulated via the burst size, while that of the σ⁵⁴-controlled stress response is regulated via the burst frequency. Therefore, transcription initiation involving σ⁵⁴ differs from other bacterial systems, and yields bursting kinetics characteristic of eukaryotic systems.

Transcription noise in bacteria is often attributed to burstiness, but the mechanisms are unclear. Here, the authors show that the transition from low to high expression can be regulated via burst size or burst frequency, depending on the mode of transcription initiation determined by different sigma factors.

Allele-specific single-cell RNA sequencing reveals different architectures of intrinsic and extrinsic gene expression noises

Gene expression noise refers to the variation of the expression level of a gene among isogenic cells in the same environment, and has two sources: extrinsic noise arising from the disparity of the cell state and intrinsic noise arising from the stochastic process of gene expression in the same cell state. Due to the low throughput of the existing method for measuring the two noise components, the architectures of intrinsic and extrinsic expression noises remain elusive. Using allele-specific single-cell RNA sequencing, we here estimate the two noise components of 3975 genes in mouse fibroblast cells. Our analyses verify predicted influences of several factors such as the TATA-box and microRNA targeting on intrinsic or extrinsic noises and reveal gene function-associated noise trends implicating the action of natural selection. These findings unravel differential regulations, optimizations, and biological consequences of intrinsic and extrinsic noises and can aid the construction of desired synthetic circuits.

Extrinsic Noise and Heavy-Tailed Laws in Gene Expression

Noise in gene expression is one of the hallmarks of life at the molecular scale. Here we derive analytical solutions to a set of models describing the molecular mechanisms underlying transcription of DNA into RNA. Our ansatz allows us to incorporate the effects of extrinsic noise-encompassing factors external to the transcription of the individual gene-and discuss the ramifications for heterogeneity in gene product abundance that has been widely observed in single cell data. Crucially, we are able to show that heavy-tailed distributions of RNA copy numbers cannot result from the intrinsic stochasticity in gene expression alone, but must instead reflect extrinsic sources of variability.