The generation of promoter-mediated transcriptional noise in bacteria

A genetic circuit on a single DNA molecule as an autonomous dissipative nanodevice

Realizing genetic circuits on single DNA molecules as self-encoded dissipative nanodevices is a major step toward miniaturization of autonomous biological systems. A circuit operating on a single DNA implies that genetically encoded proteins localize during coupled transcription-translation to DNA, but a single-molecule measurement demonstrating this has remained a challenge. Here, we use a genetically encoded fluorescent reporter system with improved temporal resolution and observe the synthesis of individual proteins tethered to a DNA molecule by transient complexes of RNA polymerase, messenger RNA, and ribosome. Against expectations in dilute cell-free conditions where equilibrium considerations favor dispersion, these nascent proteins linger long enough to regulate cascaded reactions on the same DNA. We rationally design a pulsatile genetic circuit by encoding an activator and repressor in feedback on the same DNA molecule. Driven by the local synthesis of only several proteins per hour and gene, the circuit dynamics exhibit enhanced variability between individual DNA molecules, and fluctuations with a broad power spectrum. Our results demonstrate that co-expressional localization, as a nonequilibrium process, facilitates single-DNA genetic circuits as dissipative nanodevices, with implications for nanobiotechnology applications and artificial cell design.

Transcription factor competition facilitates self-sustained oscillations in single gene genetic circuits

Genetic feedback loops can be used by cells to regulate internal processes or to keep track of time. It is often thought that, for a genetic circuit to display self-sustained oscillations, a degree of cooperativity is needed in the binding and unbinding of actor species. This cooperativity is usually modeled using a Hill function, regardless of the actual promoter architecture. Furthermore, genetic circuits do not operate in isolation and often transcription factors are shared between different promoters. In this work we show how mathematical modelling of genetic feedback loops can be facilitated with a mechanistic fold-change function that takes into account the titration effect caused by competing binding sites for transcription factors. The model shows how the titration effect facilitates self-sustained oscillations in a minimal genetic feedback loop: a gene that produces its own repressor directly without cooperative transcription factor binding. The use of delay-differential equations leads to a stability contour that predicts whether a genetic feedback loop will show self-sustained oscillations, even when taking the bursty nature of transcription into account.

Genetic, physiological, and cellular heterogeneities of bacterial pathogens in food matrices: Consequences for food safety

In complex food systems, bacteria live in heterogeneous microstructures, and the population displays phenotypic heterogeneities at the single-cell level. This review provides an overview of spatiotemporal drivers of phenotypic heterogeneity of bacterial pathogens in food matrices at three levels. The first level is the genotypic heterogeneity due to the possibility for various strains of a given species to contaminate food, each of them having specific genetic features. Then, physiological heterogeneities are induced within the same strain, due to specific microenvironments and heterogeneous adaptative responses to the food microstructure. The third level of phenotypic heterogeneity is related to cellular heterogeneity of the same strain in a specific microenvironment. Finally, we consider how these phenotypic heterogeneities at the single-cell level could be implemented in mathematical models to predict bacterial behavior and help ensure microbiological food safety.

Slow nucleosome dynamics set the transcriptional speed limit and induce RNA polymerase II traffic jams and bursts

Nucleosomes are recognized as key regulators of transcription. However, the relationship between slow nucleosome unwrapping dynamics and bulk transcriptional properties has not been thoroughly explored. Here, an agent-based model that we call the dynamic defect Totally Asymmetric Simple Exclusion Process (ddTASEP) was constructed to investigate the effects of nucleosome-induced pausing on transcriptional dynamics. Pausing due to slow nucleosome dynamics induced RNAPII convoy formation, which would cooperatively prevent nucleosome rebinding leading to bursts of transcription. The mean first passage time (MFPT) and the variance of first passage time (VFPT) were analytically expressed in terms of the nucleosome rate constants, allowing for the direct quantification of the effects of nucleosome-induced pausing on pioneering polymerase dynamics. The mean first passage elongation rate γ(h_c, h_o) is inversely proportional to the MFPT and can be considered to be a new axis of the ddTASEP phase diagram, orthogonal to the classical αβ-plane (where α and β are the initiation and termination rates). Subsequently, we showed that, for β = 1, there is a novel jamming transition in the αγ-plane that separates the ddTASEP dynamics into initiation-limited and nucleosome pausing-limited regions. We propose analytical estimates for the RNAPII density ρ, average elongation rate v, and transcription flux J and verified them numerically. We demonstrate that the intra-burst RNAPII waiting times t_in follow the time-headway distribution of a max flux TASEP and that the average inter-burst interval tIBI¯ correlates with the index of dispersion D_e. In the limit γ→0, the average burst size reaches a maximum set by the closing rate h_c. When α≪1, the burst sizes are geometrically distributed, allowing large bursts even while the average burst size NB¯ is small. Last, preliminary results on the relative effects of static and dynamic defects are presented to show that dynamic defects can induce equal or greater pausing than static bottle necks.

To perform specific functions, cells must express specific genes by copying the information in DNA into RNA via transcription. Structural proteins called nucleosomes are spaced every 200 base pairs along the length of a strand of DNA and play a crucial function in the regulation of gene activity by tightly binding DNA strands and condensing them into heterochromatin, preventing transcription by RNA polymerase II (RNAPII). Even on active genes where nucleosomes are loosely attached to DNA strands, the wrapping and unwrapping of nucleosomes pause transcription as RNAPII passes by. Previous mathematical models of transcription have compared this biological process to traffic on a one lane highway without obstructions. In contrast, our proposed model simulates transcription like traffic in a grid system where nucleosomes can be thought of as pedestrians or other vehicles crossing the road at regularly spaced intersections. Just as side street traffic and pedestrian crossings can cause cars to form convoys and cause jams limiting the max speed in an area, nucleosomes can cause RNAPII to form convoys that lead to bursts of mRNA production and limit the average polymerase flux through the gene.

Understanding the molecular mechanisms of transcriptional bursting

In recent years, it has been experimentally established that transcription, a fundamental biological process that involves the synthesis of messenger RNA molecules from DNA templates, does not proceed continuously as was expected. Rather, it exhibits a distinct dynamic behavior of alternating between productive phases when RNA molecules are actively synthesized and inactive phases when there is no RNA production at all. The bimodal transcriptional dynamics is now confirmed to be present in most living systems. This phenomenon is known as transcriptional bursting and it attracts significant amounts of attention from researchers in different fields. However, despite multiple experimental and theoretical investigations, the microscopic origin and biological functions of the transcriptional bursting remain unclear. Here we discuss the recent developments in uncovering the underlying molecular mechanisms of transcriptional bursting and our current understanding of them. Our analysis presents a physicochemical view of the processes that govern transcriptional bursting in living cells.

Probing Mechanisms of Transcription Elongation Through Cell-to-Cell Variability of RNA Polymerase

The process of transcription initiation and elongation are primary points of control in the regulation of gene expression. Although biochemical studies have uncovered the mechanisms involved in controlling transcription at each step, how these mechanisms manifest in vivo at the level of individual genes is still unclear. Recent experimental advances have enabled single-cell measurements of RNA polymerase (RNAP) molecules engaged in the process of transcribing a gene of interest. In this article, we use Gillespie simulations to show that measurements of cell-to-cell variability of RNAP numbers and interpolymerase distances can reveal the prevailing mode of regulation of a given gene. Mechanisms of regulation at each step, from initiation to elongation dynamics, produce qualitatively distinct signatures, which can further be used to discern between them. Most intriguingly, depending on the initiation kinetics, stochastic elongation can either enhance or suppress cell-to-cell variability at the RNAP level. To demonstrate the value of this framework, we analyze RNAP number distribution data for ribosomal genes in Saccharomyces cerevisiae from three previously published studies and show that this approach provides crucial mechanistic insights into the transcriptional regulation of these genes.

Dissecting the in vivo dynamics of transcription locking due to positive supercoiling buildup

Positive supercoiling buildup (PSB) is a pervasive phenomenon in the transcriptional programs of Escherichia coli. After finding a range of Gyrase concentrations where the inverse of the transcription rate of a chromosome-integrated gene changes linearly with the inverse of Gyrase concentration, we apply a LineWeaver-Burk plot to dissect the expected in vivo transcription rate in absence of PSB. We validate the estimation by time-lapse microscopy of single-RNA production kinetics of the same gene when single-copy plasmid-borne, shown to be impervious to Gyrase inhibition. Next, we estimate the fraction of time in locked states and number of transcription events prior to locking, which we validate by measurements under Gyrase inhibition. Replacing the gene of interest by one with slower transcription rate decreases the fraction of time in locked states due to PSB. Finally, we combine data from both constructs to infer a range of possible transcription initiation locking kinetics in a chromosomal location, obtainable by tuning the transcription rate. We validate with measurements of transcription activity at different induction levels. This strategy for dissecting transcription initiation locking kinetics due to PSB can contribute to resolve the transcriptional programs of E. coli and in the engineering of synthetic genetic circuits.

A Mechanochemical Model of Transcriptional Bursting

Populations of genetically identical cells generally show a large variability in cell phenotypes, which is typically associated with the stochastic nature of gene expression processes. It is widely believed that a significant source of such randomness is transcriptional bursting, which is when periods of active production of RNA molecules alternate with periods of RNA degradation. However, the molecular mechanisms of such strong fluctuations remain unclear. Recent studies suggest that DNA supercoiling, which happens during transcription, might be directly related to the bursting behavior. Stimulated by these observations, we developed a stochastic mechanochemical model of supercoiling-induced transcriptional bursting in which the RNA synthesis leads to the buildup of torsion in DNA. This slows down the RNA production until it is bound by the enzyme gyrase to DNA, which releases the stress and allows for the RNA synthesis to restart with the original rate. Using a thermodynamically consistent coupling between mechanical and chemical processes, the dynamic properties of transcription are explicitly evaluated. In addition, a first-passage method to evaluate the dynamics of transcription is developed. Theoretical analysis shows that transcriptional bursting is observed when both the supercoiling and the mechanical stress release due to gyrase are present in the system. It is also found that the overall RNA production rate is not constant and depends on the number of previously synthesized RNA molecules. A comparison with experimental data on bacteria allows us to evaluate the energetic cost of supercoiling during transcription. It is argued that the relatively weak mechanochemical coupling might allow transcription to be regulated most effectively.

CarD and RbpA modify the kinetics of initial transcription and slow promoter escape of the Mycobacterium tuberculosis RNA polymerase

The pathogen Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, enacts unique transcriptional regulatory mechanisms when subjected to host-derived stresses. Initiation of transcription by the Mycobacterial RNA polymerase (RNAP) has previously been shown to exhibit different open complex kinetics and stabilities relative to Escherichia coli (Eco) RNAP. However, transcription initiation rates also depend on the kinetics following open complex formation such as initial nucleotide incorporation and subsequent promoter escape. Here, using a real-time fluorescence assay, we present the first in-depth kinetic analysis of initial transcription and promoter escape for the Mtb RNAP. We show that in relation to Eco RNAP, Mtb displays slower initial nucleotide incorporation but faster overall promoter escape kinetics on the Mtb rrnAP3 promoter. Furthermore, in the context of the essential transcription factors CarD and RbpA, Mtb promoter escape is slowed via differential effects on initially transcribing complexes. Finally, based on their ability to increase the rate of open complex formation and decrease the rate of promoter escape, we suggest that CarD and RbpA are capable of activation or repression depending on the rate-limiting step of a given promoter's basal initiation kinetics.

Distribution of Initiation Times Reveals Mechanisms of Transcriptional Regulation in Single Cells

Transcription is the dominant point of control of gene expression. Biochemical studies have revealed key molecular components of transcription and their interactions, but the dynamics of transcription initiation in cells is still poorly understood. This state of affairs is being remedied with experiments that observe transcriptional dynamics in single cells using fluorescent reporters. Quantitative information about transcription initiation dynamics can also be extracted from experiments that use electron micrographs of RNA polymerases caught in the act of transcribing a gene (Miller spreads). Inspired by these data, we analyze a general stochastic model of transcription initiation and elongation and compute the distribution of transcription initiation times. We show that different mechanisms of initiation leave distinct signatures in the distribution of initiation times that can be compared to experiments. We analyze published data from micrographs of RNA polymerases transcribing ribosomal RNA genes in Escherichia coli and compare the observed distributions of interpolymerase distances with the predictions from previously hypothesized mechanisms for the regulation of these genes. Our analysis demonstrates the potential of measuring the distribution of time intervals between initiation events as a probe for dissecting mechanisms of transcription initiation in live cells.

Transcriptional Bursts in a Nonequilibrium Model for Gene Regulation by Supercoiling

We analyze transcriptional bursting within a stochastic nonequilibrium model, which accounts for the coupling between the dynamics of DNA supercoiling and gene transcription. We find a clear signature of bursty transcription when there is a separation between the timescales of transcription initiation and supercoiling dissipation (the latter may either be diffusive or mediated by topological enzymes, such as type I or type II topoisomerases). In multigenic DNA domains, we observe either bursty transcription or transcription waves; the type of behavior can be selected for by controlling gene activity and orientation. In the bursty phase, the statistics of supercoiling fluctuations at the promoter are markedly non-Gaussian.

Noise in bacterial gene expression

The expression level of a gene can fluctuate significantly between individuals within a population of genetically identical cells. The resultant phenotypic heterogeneity could be exploited by bacteria to adapt to changing environmental conditions. Noise is hence a genome-wide phenomenon that arises from the stochastic nature of the biochemical reactions that take place during gene expression and the relatively low abundance of the molecules involved. The production of mRNA and proteins therefore occurs in bursts, with alternating episodes of high and low activity during transcription and translation. Single-cell and single-molecule studies demonstrated that noise within gene expression is influenced by a combination of both intrinsic and extrinsic factors. However, our mechanistic understanding of this process at the molecular level is still rather limited. Further investigation is necessary that takes into account the detailed knowledge of gene regulation gained from biochemical studies.

Elucidating effects of reaction rates on dynamics of the lac circuit in Escherichia coli

Gene expression is regulated by a complex transcriptional network. It is of interest to quantify uncertainty of not knowing accurately reaction rates of underlying biochemical reactions, and to understand how they affect gene expression. Assuming a kinetic model of the lac circuit in Escherichia coli, regardless of how many reactions are involved in transcription regulation, transcription rate is shown to be the most important parameter affecting steady state production of mRNA and protein in the cell. In particular, doubling the transcription rate approximately doubles the number of mRNA synthesized at steady state for any rates of transcription inhibition and activation. On the other hand, increasing the rate of transcription inhibition by 10% reduces the average steady state count of mRNA by about 7%, whereas changes in the rate of transcription activation appear to have no such effect. Furthermore, for wide range of reaction rates in the kinetic model of the lac genetic switch considered, protein production was observed to always reach a maximum before the degradation reduces its count to zero, and this maximum was found to be always at least 27 protein molecules. Such value appears to be a fundamental structural property of genetic circuits making it very robust against changes in the internal and external conditions.

Added value of autoregulation and multi-step kinetics of transcription initiation

Bacterial gene expression regulation occurs mostly during transcription, which has two main rate-limiting steps: the close complex formation, when the RNA polymerase binds to an active promoter, and the subsequent open complex formation, after which it follows elongation. Tuning these steps' kinetics by the action of e.g. transcription factors, allows for a wide diversity of dynamics. For example, adding autoregulation generates single-gene circuits able to perform more complex tasks. Using stochastic models of transcription kinetics with empirically validated parameter values, we investigate how autoregulation and the multi-step transcription initiation kinetics of single-gene autoregulated circuits can be combined to fine-tune steady state mean and cell-to-cell variability in protein expression levels, as well as response times. Next, we investigate how they can be jointly tuned to control complex behaviours, namely, time counting, switching dynamics and memory storage. Overall, our finding suggests that, in bacteria, jointly regulating a single-gene circuit's topology and the transcription initiation multi-step dynamics allows enhancing complex task performance.

Bursting onto the scene? Exploring stochastic mRNA production in bacteria

Recent large-scale measurements of gene expression variability (or noise) in E. coli have led to the unexpected conclusion that the variability is in large part dictated by and increasing with the mean level of expression. Here we review the evidence for this apparent universal trend in variability, as well as for the related idea that transcription is fundamentally bursty. We examine recently proposed mechanisms for burstiness and universality and argue that they do not explain important features of observed data. Finally, we discuss potential limitations and pitfalls in the interpretation of experimental measurements of cell-to-cell variability.

Self-consistent theory of transcriptional control in complex regulatory architectures

Individual regulatory proteins are typically charged with the simultaneous regulation of a battery of different genes. As a result, when one of these proteins is limiting, competitive effects have a significant impact on the transcriptional response of the regulated genes. Here we present a general framework for the analysis of any generic regulatory architecture that accounts for the competitive effects of the regulatory environment by isolating these effects into an effective concentration parameter. These predictions are formulated using the grand-canonical ensemble of statistical mechanics and the fold-change in gene expression is predicted as a function of the number of transcription factors, the strength of interactions between the transcription factors and their DNA binding sites, and the effective concentration of the transcription factor. The effective concentration is set by the transcription factor interactions with competing binding sites within the cell and is determined self-consistently. Using this approach, we analyze regulatory architectures in the grand-canonical ensemble ranging from simple repression and simple activation to scenarios that include repression mediated by DNA looping of distal regulatory sites. It is demonstrated that all the canonical expressions previously derived in the case of an isolated, non-competing gene, can be generalised by a simple substitution to their grand canonical counterpart, which allows for simple intuitive incorporation of the influence of multiple competing transcription factor binding sites. As an example of the strength of this approach, we build on these results to present an analytical description of transcriptional regulation of the lac operon.

Models of life: epigenetics, diversity and cycles

This review emphasizes aspects of biology that can be understood through repeated applications of simple causal rules. The selected topics include perspectives on gene regulation, phage lambda development, epigenetics, microbial ecology, as well as model approaches to diversity and to punctuated equilibrium in evolution. Two outstanding features are repeatedly described. One is the minimal number of rules to sustain specific states of complex systems for a long time. The other is the collapse of such states and the subsequent dynamical cycle of situations that restitute the system to a potentially new metastable state.

Dissecting the stochastic transcription initiation process in live Escherichia coli

We investigate the hypothesis that, in Escherichia coli, while the concentration of RNA polymerases differs in different growth conditions, the fraction of RNA polymerases free for transcription remains approximately constant within a certain range of these conditions. After establishing this, we apply a standard model-fitting procedure to fully characterize the in vivo kinetics of the rate-limiting steps in transcription initiation of the Plac/ara-1 promoter from distributions of intervals between transcription events in cells with different RNA polymerase concentrations. We find that, under full induction, the closed complex lasts ∼788 s while subsequent steps last ∼193 s, on average. We then establish that the closed complex formation usually occurs multiple times prior to each successful initiation event. Furthermore, the promoter intermittently switches to an inactive state that, on average, lasts ∼87 s. This is shown to arise from the intermittent repression of the promoter by LacI. The methods employed here should be of use to resolve the rate-limiting steps governing the in vivo dynamics of initiation of prokaryotic promoters, similar to established steady-state assays to resolve the in vitro dynamics.

Deciphering Transcriptional Dynamics In Vivo by Counting Nascent RNA Molecules

Deciphering how the regulatory DNA sequence of a gene dictates its expression in response to intra and extracellular cues is one of the leading challenges in modern genomics. The development of novel single-cell sequencing and imaging techniques, as well as a better exploitation of currently available single-molecule imaging techniques, provides an avenue to interrogate the process of transcription and its dynamics in cells by quantifying the number of RNA polymerases engaged in the transcription of a gene (or equivalently the number of nascent RNAs) at a given moment in time. In this paper, we propose that measurements of the cell-to-cell variability in the number of nascent RNAs provide a mostly unexplored method for deciphering mechanisms of transcription initiation in cells. We propose a simple kinetic model of transcription initiation and elongation from which we calculate nascent RNA copy-number fluctuations. To demonstrate the usefulness of this approach, we test our theory against published nascent RNA data for twelve constitutively expressed yeast genes. Rather than transcription being initiated through a single rate limiting step, as it had been previously proposed, our single-cell analysis reveals the presence of at least two rate limiting steps. Surprisingly, half of the genes analyzed have nearly identical rates of transcription initiation, suggesting a common mechanism. Our analytical framework can be used to extract quantitative information about dynamics of transcription from single-cell sequencing data, as well as from single-molecule imaging and electron micrographs of fixed cells, and provides the mathematical means to exploit the quantitative power of these technologies.

Interplay between DNA supercoiling and transcription elongation

Transcription-coupled DNA supercoiling has been shown to be an important regulator of transcription that is broadly present in the cell. Here we review experimental work which shows that RNA polymerase is a powerful torsional motor that can alter DNA topology and structure, and DNA supercoiling in turn directly affects transcription elongation.

The Low Noise Limit in Gene Expression

Protein noise measurements are increasingly used to elucidate biophysical parameters. Unfortunately noise analyses are often at odds with directly measured parameters. Here we show that these inconsistencies arise from two problematic analytical choices: (i) the assumption that protein translation rate is invariant for different proteins of different abundances, which has inadvertently led to (ii) the assumption that a large constitutive extrinsic noise sets the low noise limit in gene expression. While growing evidence suggests that transcriptional bursting may set the low noise limit, variability in translational bursting has been largely ignored. We show that genome-wide systematic variation in translational efficiency can–and in the case of E. coli does–control the low noise limit in gene expression. Therefore constitutive extrinsic noise is small and only plays a role in the absence of a systematic variation in translational efficiency. These results show the existence of two distinct expression noise patterns: (1) a global noise floor uniformly imposed on all genes by expression bursting; and (2) high noise distributed to only a select group of genes.

Dynamic competition between transcription initiation and repression: Role of nonequilibrium steps in cell-to-cell heterogeneity

Transcriptional repression may cause transcriptional noise by a competition between repressor and RNA polymerase binding. Although promoter activity is often governed by a single limiting step, we argue here that the size of the noise strongly depends on whether this step is the initial equilibrium binding or one of the subsequent unidirectional steps. Overall, we show that nonequilibrium steps of transcription initiation systematically increase the cell-to-cell heterogeneity in bacterial populations. In particular, this allows also weak promoters to give substantial transcriptional noise.

Promoter architecture dictates cell-to-cell variability in gene expression

Variability in gene expression among genetically identical cells has emerged as a central preoccupation in the study of gene regulation; however, a divide exists between the predictions of molecular models of prokaryotic transcriptional regulation and genome-wide experimental studies suggesting that this variability is indifferent to the underlying regulatory architecture. We constructed a set of promoters in Escherichia coli in which promoter strength, transcription factor binding strength, and transcription factor copy numbers are systematically varied, and used messenger RNA (mRNA) fluorescence in situ hybridization to observe how these changes affected variability in gene expression. Our parameter-free models predicted the observed variability; hence, the molecular details of transcription dictate variability in mRNA expression, and transcriptional noise is specifically tunable and thus represents an evolutionarily accessible phenotypic parameter.

Mechanism of transcriptional bursting in bacteria

Transcription of highly expressed genes has been shown to occur in stochastic bursts. But the origin of such ubiquitous phenomenon has not been understood. Here, we present the mechanism in bacteria. We developed a high-throughput, in vitro, single-molecule assay to follow transcription on individual DNA templates in real time. We showed that positive supercoiling buildup on a DNA segment by transcription slows down transcription elongation and eventually stops transcription initiation. Transcription can be resumed upon gyrase binding to the DNA segment. Furthermore, using single-cell mRNA counting fluorescence in situ hybridization (FISH), we found that duty cycles of transcriptional bursting depend on the intracellular gyrase concentration. Together, these findings prove that transcriptional bursting of highly expressed genes in bacteria is primarily caused by reversible gyrase dissociation from and rebinding to a DNA segment, changing the supercoiling level of the segment.

FOUNDATIONS FOR MODELING THE DYNAMICS OF GENE REGULATORY NETWORKS: A MULTILEVEL-PERSPECTIVE REVIEW

A promising alternative for unraveling the principles under which the dynamic interactions among genes lead to cellular phenotypes relies on mathematical and computational models at different levels of abstraction, from the molecular level of protein-DNA interactions to the system level of functional relationships among genes. This review article presents, under a bottom–up perspective, a hierarchy of approaches to modeling gene regulatory network dynamics, from microscopic descriptions at the single-molecule level in the spatial context of an individual cell to macroscopic models providing phenomenological descriptions at the population-average level. The reviewed modeling approaches include Molecular Dynamics, Particle-Based Brownian Dynamics, the Master Equation approach, Ordinary Differential Equations, and the Boolean logic abstraction. Each of these frameworks is motivated by a particular biological context and the nature of the insight being pursued. The setting of gene network dynamic models from such frameworks involves assumptions and mathematical artifacts often ignored by the non-specialist. This article aims at providing an entry point for biologists new to the field and computer scientists not acquainted with some recent biophysically-inspired models of gene regulation. The connections promoting intuition between different abstraction levels and the role that approximations play in the modeling process are highlighted throughout the paper.

Phenotypic noise: effects of post-transcriptional regulatory processes affecting mRNA

The inherently stochastic nature of biomolecular processes is one of the main sources giving rise to cell-to-cell variations in protein and mRNA abundance, termed noise. Noise in isogenic populations can enhance survival under adverse conditions and stress, and has therefore played a fundamental role in evolution. On the other hand, noise may have detrimental effects and therefore cells must also display robustness to fluctuations and possess mechanisms of control in order to function properly. Noise can be introduced at every step in the cascade of intermediate events resulting in the production of functional proteins. While initial studies of noise focused on stochasticity introduced at the transcriptional level, recent years have witnessed a gradual shift of emphasis into the effects that post-transcriptional processes have on phenotypic noise. Here, we survey the insights that have been gained on the effects of processes that modify RNA transcript populations on phenotypic noise, including regulation by noncoding RNAs in prokaryotes and eukaryotes, alternative splicing and transcriptional interference.

Regulation of noise in gene expression

The biochemical processes leading to the synthesis of new proteins are random, as they typically involve a small number of diffusing molecules. They lead to fluctuations in the number of proteins in a single cell as a function of time and to cell-to-cell variability of protein abundances. These in turn can lead to phenotypic heterogeneity in a population of genetically identical cells. Phenotypic heterogeneity may have important consequences for the development of multicellular organisms and the fitness of bacterial colonies, raising the question of how it is regulated. Here we review the experimental evidence that transcriptional regulation affects noise in gene expression, and discuss how the noise strength is encoded in the architecture of the promoter region. We discuss how models based on specific molecular mechanisms of gene regulation can make experimentally testable predictions for how changes to the promoter architecture are reflected in gene expression noise.

Single-molecule methods for studying gene regulation in vivo

The recent emergence of new experimental tools employing sensitive fluorescence detection in vivo has made it possible to visualize various aspects of gene regulation at the single-molecule level in the native, intracellular context. In this review, we will first describe general considerations for in vivo, single-molecule fluorescence detection of DNA, mRNA, and protein molecules involved in gene regulation. We will then give an overview of the rapidly evolving suite of molecular tools available for observing gene regulation in vivo and discuss new insights they have brought into gene regulation.

Interplay of gene expression noise and ultrasensitive dynamics affects bacterial operon organization

Bacterial chromosomes are organized into polycistronic cotranscribed operons, but the evolutionary pressures maintaining them are unclear. We hypothesized that operons alter gene expression noise characteristics, resulting in selection for or against maintaining operons depending on network architecture. Mathematical models for 6 functional classes of network modules showed that three classes exhibited decreased noise and 3 exhibited increased noise with same-operon cotranscription of interacting proteins. Noise reduction was often associated with a decreased chance of reaching an ultrasensitive threshold. Stochastic simulations of the lac operon demonstrated that the predicted effects of transcriptional coupling hold for a complex network module. We employed bioinformatic analysis to find overrepresentation of noise-minimizing operon organization compared with randomized controls. Among constitutively expressed physically interacting protein pairs, higher coupling frequencies appeared at lower expression levels, where noise effects are expected to be dominant. Our results thereby suggest an important role for gene expression noise, in many cases interacting with an ultrasensitive switch, in maintaining or selecting for operons in bacterial chromosomes.

In some species, most notably bacteria, chromosomal genes are arranged into clusters called operons. In operons, the process of transcription is physically coupled: a single pass of the RNA polymerase enzyme reading that region of the chromosome simultaneously produces messenger RNA encoding multiple proteins. So far, we do not have a satisfying explanation for what evolutionary forces have maintained operons on bacterial chromosomes. We hypothesized that different types of interactions between operon-coded proteins affect how strongly operons are selected for between two genes. The proposed mechanism for this effect is that operons correlate gene expression noise, changing how it manifests in the post-translational network depending on the type of protein interaction. Mathematical models demonstrate that operons reduce noise for some types of interactions but not others. We found that operon-dependent noise reduction has an underlying dependence on surprisingly high sensitivity of the network to the ratio of proteins from each gene. Databases of genetic information show that E. coli has operons more frequently than random if operons reduce noise for the type of interaction various gene pairs have, but not otherwise. Our study thus provides an example of how the architecture of post-translational networks affects bacterial evolution.

Quantitative single-cell ion-channel gene expression profiling through an improved qRT-PCR technique combined with whole cell patch clamp

Cellular excitability originates from a concerted action of different ion channels. The genomic diversity of ion channels (over 100 different genes) underlies the functional diversity of neurons in the central nervous system (CNS) and even within a specific type of neurons large differences in channel expression have been observed. Patch-clamp is a powerful technique to study the electrophysiology of excitability at the single cell level, allowing exploration of cell-to-cell variability. Only a few attempts have been made to link electrophysiological profiling to mRNA transcript levels and most suffered from experimental noise precluding conclusive quantitative correlations. Here we describe a refinement to the technique that combines patch-clamp analysis with quantitative real-time (qRT) PCR at the single cell level. Hereto the expression of a housekeeping gene was used to normalize for cell-to-cell variability in mRNA isolation and the subsequent processing steps for performing qRT-PCR. However, the mRNA yield from a single cell was insufficient for performing a valid qRT-PCR assay; this was resolved by including a RNA amplification step. The technique was validated on a stable Ltk(-) cell line expressing the Kv2.1 channel and on embryonic dorsal root ganglion (DRG) cells probing for the expression of Kv2.1. Current density and transcript quantity displayed a clear correlation when the qRT-PCR assay was done in twofold and the data normalized to the transcript level of the housekeeping gene GAPD. Without this normalization no significant correlation was obtained. This improved technique should prove very valuable for studying the molecular background of diversity in cellular excitability.

General properties of transcriptional time series in Escherichia coli

Gene activity is described by the time-series of discrete, stochastic mRNA production events. This transcriptional time-series exhibits intermittent, bursty behavior. One consequence of this temporal intricacy is that gene expression can be tuned by varying different features of the time-series. What schemes for varying the transcriptional time-series are observed in the cell? Are the observed properties of these time-series optimized for cellular function? To address these questions, we characterize mRNA copy-number statistics at single-molecule resolution from multiple Escherichia coli promoters. We find that the degree of burstiness depends only on the gene expression level, while being independent of the details of gene regulation. The observed behavior is explained by the underlying variation in the duration of bursting events. Using information theory, we find that the properties of the transcriptional time series allow the cell to efficiently map the extracellular concentration of inducer molecules to intracellular levels of mRNA and proteins.

Identification of Vibrio harveyi as a causative bacterium for a tail rot disease of sea bream Sparus aurata from research hatchery in Malta

A bacterial disease was reported from gilthead sea bream (Sparus aurata) within a hatchery environment in Malta. Symptoms included complete erosion of tail, infection in the eye, mucous secretion and frequent mortality. A total of 540 strains were initially isolated in marine agar from different infected body parts and culture water sources. Subsequently 100 isolates were randomly selected, identified biochemically and all were found to be Vibrio harveyi-related organisms; finally from 100 isolates a total of 13 numbers were randomly selected and accurately identified as V. harveyi by 16S rRNA gene sequencing and species-specific PCR. Ribotyping of these strains with HindIII revealed total of six clusters. In vivo challenge study with representative isolates from each cluster proved two clusters each were highly pathogenic, moderately pathogenic and non-pathogenic. All 13 isolates were positive for hemolysin gene, a potential virulence factor. Further analysis revealed probably a single copy of this gene was encoded in all isolates, although not in the same locus in the genome. Although V. harveyi was reported to be an important pathogen for many aquatic organisms, to our knowledge this might be the first report of disease caused by V. harveyi and their systematic study in the sea bream hatchery from Malta.

On the spontaneous stochastic dynamics of a single gene: complexity of the molecular interplay at the promoter

BACKGROUND

Gene promoters can be in various epigenetic states and undergo interactions with many molecules in a highly transient, probabilistic and combinatorial way, resulting in a complex global dynamics as observed experimentally. However, models of stochastic gene expression commonly consider promoter activity as a two-state on/off system. We consider here a model of single-gene stochastic expression that can represent arbitrary prokaryotic or eukaryotic promoters, based on the combinatorial interplay between molecules and epigenetic factors, including energy-dependent remodeling and enzymatic activities.

RESULTS

We show that, considering the mere molecular interplay at the promoter, a single-gene can demonstrate an elaborate spontaneous stochastic activity (eg. multi-periodic multi-relaxation dynamics), similar to what is known to occur at the gene-network level. Characterizing this generic model with indicators of dynamic and steady-state properties (including power spectra and distributions), we reveal the potential activity of any promoter and its influence on gene expression. In particular, we can reproduce, based on biologically relevant mechanisms, the strongly periodic patterns of promoter occupancy by transcription factors (TF) and chromatin remodeling as observed experimentally on eukaryotic promoters. Moreover, we link several of its characteristics to properties of the underlying biochemical system. The model can also be used to identify behaviors of interest (eg. stochasticity induced by high TF concentration) on minimal systems and to test their relevance in larger and more realistic systems. We finally show that TF concentrations can regulate many aspects of the stochastic activity with a considerable flexibility and complexity.

CONCLUSIONS

This tight promoter-mediated control of stochasticity may constitute a powerful asset for the cell. Remarkably, a strongly periodic activity that demonstrates a complex TF concentration-dependent control is obtained when molecular interactions have typical characteristics observed on eukaryotic promoters (high mobility, functional redundancy, many alternate states/pathways). We also show that this regime results in a direct and indirect energetic cost. Finally, this model can constitute a framework for unifying various experimental approaches. Collectively, our results show that a gene - the basic building block of complex regulatory networks - can itself demonstrate a significantly complex behavior.

Noise management by molecular networks

Fluctuations in the copy number of key regulatory macromolecules ("noise") may cause physiological heterogeneity in populations of (isogenic) cells. The kinetics of processes and their wiring in molecular networks can modulate this molecular noise. Here we present a theoretical framework to study the principles of noise management by the molecular networks in living cells. The theory makes use of the natural, hierarchical organization of those networks and makes their noise management more understandable in terms of network structure. Principles governing noise management by ultrasensitive systems, signaling cascades, gene networks and feedback circuitry are discovered using this approach. For a few frequently occurring network motifs we show how they manage noise. We derive simple and intuitive equations for noise in molecule copy numbers as a determinant of physiological heterogeneity. We show how noise levels and signal sensitivity can be set independently in molecular networks, but often changes in signal sensitivity affect noise propagation. Using theory and simulations, we show that negative feedback can both enhance and reduce noise. We identify a trade-off; noise reduction in one molecular intermediate by negative feedback is at the expense of increased noise in the levels of other molecules along the feedback loop. The reactants of the processes that are strongly (cooperatively) regulated, so as to allow for negative feedback with a high strength, will display enhanced noise.

Analysis of stochasticity in promoter activation by using a dual-fluorescence reporter system

Stochastic dynamics of promoter activity in bacterial cells were studied by using a dual-fluorescence reporter system of protein expression. The dual-fluorescence reporter system enabled us to calculate the amplitude of intrinsic noise generated during transcription and translation. By fitting the experimental data to a simple stochastic model of protein expression, we could estimate parameters representing the stochastic transition between the active and inactive states of a promoter. Using the system, we analyzed the stochastic dynamics of promoter activation of genes in the lysine biosynthesis pathway in Escherichia coli. We found that the promoter of lysA has a significantly slower transition rate between active and inactive states than other promoters in the lysine biosynthesis pathway. The infrequent switching between active and inactive states can be a dominant source of noise in lysA expression. Analysis using the dual-fluorescence reporter system provided a better understanding of stochastic dynamics in promoter activation.

Elongation dynamics shape bursty transcription and translation

Cells in isogenic populations may differ substantially in their molecular make up because of the stochastic nature of molecular processes. Stochastic bursts in process activity have a great potential for generating molecular noise. They are characterized by (short) periods of high process activity followed by (long) periods of process silence causing different cells to experience activity periods varying in size, duration, and timing. We present an analytically solvable model of bursts in molecular networks, originally developed for the analysis of telecommunication networks. We define general measures for model-independent characterization of bursts (burst size, significance, and duration) from stochastic time series. Inspired by the discovery of bursts in mRNA and protein production by others, we use those indices to investigate the role of stochastic motion of motor proteins along biopolymer chains in determining burst properties. Collisions between neighboring motor proteins can attenuate bursts introduced at the initiation site on the chain. Pausing of motor proteins can give rise to bursts. We investigate how these effects are modulated by the length of the biopolymer chain and the kinetic properties of motion. We discuss the consequences of those results for transcription and translation.

Dynamical analysis on gene activity in the presence of repressors and an interfering promoter

Transcription is regulated through interplay among transcription factors, an RNA polymerase (RNAP), and a promoter. Even for a simple repressive transcription factor that disturbs promoter activity at initial binding of RNAP, its repression level is not determined solely by the dissociation constant of transcription factor but is sensitive to timescales of processes in RNAP. We first analyze the promoter activity under strong repression by a slow binding repressor, in which case transcription events occur in bursts, followed by long quiescent periods while a repressor binds to the operator; the number of transcription events, bursting, and quiescent times are estimated by reaction rates. We then examine interference effect from an opposing promoter, using the correlation function of initiation events for a single promoter. The interference is shown to de-repress the promoter because RNAPs from the opposing promoter most likely encounter the repressor and remove it in case of strong repression. This de-repression mechanism should be especially prominent for the promoters that facilitate fast formation of open complex with the repressor whose binding rate is slower than approximately 1/s. Finally, we discuss possibility of this mechanism for high activity of promoter PR in the hyp-mutant of lambda-phage.