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

Understanding Biological Regulation Through Synthetic Biology

Engineering synthetic gene regulatory circuits proceeds through iterative cycles of design, building, and testing. Initial circuit designs must rely on often-incomplete models of regulation established by fields of reductive inquiry—biochemistry and molecular and systems biology. As differences in designed and experimentally observed circuit behavior are inevitably encountered, investigated, and resolved, each turn of the engineering cycle can force a resynthesis in understanding of natural network function. Here, we outline research that uses the process of gene circuit engineering to advance biological discovery. Synthetic gene circuit engineering research has not only refined our understanding of cellular regulation but furnished biologists with a toolkit that can be directed at natural systems to exact precision manipulation of network structure. As we discuss, using circuit engineering to predictively reorganize, rewire, and reconstruct cellular regulation serves as the ultimate means of testing and understanding how cellular phenotype emerges from systems-level network function.

Simulating biological processes: stochastic physics from whole cells to colonies

The last few decades have revealed the living cell to be a crowded spatially heterogeneous space teeming with biomolecules whose concentrations and activities are governed by intrinsically random forces. It is from this randomness, however, that a vast array of precisely timed and intricately coordinated biological functions emerge that give rise to the complex forms and behaviors we see in the biosphere around us. This seemingly paradoxical nature of life has drawn the interest of an increasing number of physicists, and recent years have seen stochastic modeling grow into a major subdiscipline within biological physics. Here we review some of the major advances that have shaped our understanding of stochasticity in biology. We begin with some historical context, outlining a string of important experimental results that motivated the development of stochastic modeling. We then embark upon a fairly rigorous treatment of the simulation methods that are currently available for the treatment of stochastic biological models, with an eye toward comparing and contrasting their realms of applicability, and the care that must be taken when parameterizing them. Following that, we describe how stochasticity impacts several key biological functions, including transcription, translation, ribosome biogenesis, chromosome replication, and metabolism, before considering how the functions may be coupled into a comprehensive model of a 'minimal cell'. Finally, we close with our expectation for the future of the field, focusing on how mesoscopic stochastic methods may be augmented with atomic-scale molecular modeling approaches in order to understand life across a range of length and time scales.

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.

Behavioral Variability and Phenotypic Diversity in Bacterial Chemotaxis

Living cells detect and process external signals using signaling pathways that are affected by random fluctuations. These variations cause the behavior of individual cells to fluctuate over time (behavioral variability) and generate phenotypic differences between genetically identical individuals (phenotypic diversity). These two noise sources reduce our ability to predict biological behavior because they diversify cellular responses to identical signals. Here, we review recent experimental and theoretical advances in understanding the mechanistic origin and functional consequences of such variation in Escherichia coli chemotaxis-a well-understood model of signal transduction and behavior. After briefly summarizing the architecture and logic of the chemotaxis system, we discuss determinants of behavior and chemotactic performance of individual cells. Then, we review how cell-to-cell differences in protein abundance map onto differences in individual chemotactic abilities and how phenotypic variability affects the performance of the population. We conclude with open questions to be addressed by future research.

Tuning Transcriptional Regulation through Signaling: A Predictive Theory of Allosteric Induction

Allosteric regulation is found across all domains of life, yet we still lack simple, predictive theories that directly link the experimentally tunable parameters of a system to its input-output response. To that end, we present a general theory of allosteric transcriptional regulation using the Monod-Wyman-Changeux model. We rigorously test this model using the ubiquitous simple repression motif in bacteria by first predicting the behavior of strains that span a large range of repressor copy numbers and DNA binding strengths and then constructing and measuring their response. Our model not only accurately captures the induction profiles of these strains, but also enables us to derive analytic expressions for key properties such as the dynamic range and [EC₅₀]. Finally, we derive an expression for the free energy of allosteric repressors that enables us to collapse our experimental data onto a single master curve that captures the diverse phenomenology of the induction profiles.

Antibiotic-induced population fluctuations and stochastic clearance of bacteria

Effective antibiotic use that minimizes treatment failures remains a challenge. A better understanding of how bacterial populations respond to antibiotics is necessary. Previous studies of large bacterial populations established the deterministic framework of pharmacodynamics. Here, characterizing the dynamics of population extinction, we demonstrated the stochastic nature of eradicating bacteria with antibiotics. Antibiotics known to kill bacteria (bactericidal) induced population fluctuations. Thus, at high antibiotic concentrations, the dynamics of bacterial clearance were heterogeneous. At low concentrations, clearance still occurred with a non-zero probability. These striking outcomes of population fluctuations were well captured by our probabilistic model. Our model further suggested a strategy to facilitate eradication by increasing extinction probability. We experimentally tested this prediction for antibiotic-susceptible and clinically-isolated resistant bacteria. This new knowledge exposes fundamental limits in our ability to predict bacterial eradication. Additionally, it demonstrates the potential of using antibiotic concentrations that were previously deemed inefficacious to eradicate bacteria.

Increased gene expression noise in human cancers is correlated with low p53 and immune activities as well as late stage cancer

Gene expression in metazoans is delicately organized. As genetic information transmits from DNA to RNA and protein, expression noise is inevitably generated. Recent studies begin to unveil the mechanisms of gene expression noise control, but the changes of gene expression precision in pathologic conditions like cancers are unknown. Here we analyzed the transcriptomic data of human breast, liver, lung and colon cancers, and found that the expression noise of more than 74.9% genes was increased in cancer tissues as compared to adjacent normal tissues. This suggested that gene expression precision controlling collapsed during cancer development. A set of 269 genes with noise increased more than 2-fold were identified across different cancer types. These genes were involved in cell adhesion, catalytic and metabolic functions, implying the vulnerability of deregulation of these processes in cancers. We also observed a tendency of increased expression noise in patients with low p53 and immune activity in breast, liver and lung caners but not in colon cancers, which indicated the contributions of p53 signaling and host immune surveillance to gene expression noise in cancers. Moreover, more than 53.7% genes had increased noise in patients with late stage than early stage cancers, suggesting that gene expression precision was associated with cancer outcome. Together, these results provided genomic scale explorations of gene expression noise control in human cancers.

The Chemical Fluctuation Theorem governing gene expression

Gene expression is a complex stochastic process composed of numerous enzymatic reactions with rates coupled to hidden cell-state variables. Despite advances in single-cell technologies, the lack of a theory accurately describing the gene expression process has restricted a robust, quantitative understanding of gene expression variability among cells. Here we present the Chemical Fluctuation Theorem (CFT), providing an accurate relationship between the environment-coupled chemical dynamics of gene expression and gene expression variability. Combined with a general, accurate model of environment-coupled transcription processes, the CFT provides a unified explanation of mRNA variability for various experimental systems. From this analysis, we construct a quantitative model of transcription dynamics enabling analytic predictions for the dependence of mRNA noise on the mRNA lifetime distribution, confirmed against stochastic simulation. This work suggests promising new directions for quantitative investigation into cellular control over biological functions by making complex dynamics of intracellular reactions accessible to rigorous mathematical deductions.

A unified framework to understand gene expression noise is still lacking. Here the authors derive a universal theorem relating the biological noise with dynamics of birth and death processes and present a model of transcription dynamics, allowing analytical prediction of the dependence of mRNA noise on mRNA lifetime variability.

Computational Investigation of Environment-Noise Interaction in Single-Cell Organisms: The Merit of Expression Stochasticity Depends on the Quality of Environmental Fluctuations

Organisms need to adapt to changing environments and they do so by using a broad spectrum of strategies. These strategies include finding the right balance between expressing genes before or when they are needed, and adjusting the degree of noise inherent in gene expression. We investigated the interplay between different nutritional environments and the inhabiting organisms’ metabolic and genetic adaptations by applying an evolutionary algorithm to an agent-based model of a concise bacterial metabolism. Our results show that constant environments and rapidly fluctuating environments produce similar adaptations in the organisms, making the predictability of the environment a major factor in determining optimal adaptation. We show that exploitation of expression noise occurs only in some types of fluctuating environment and is strongly dependent on the quality and availability of nutrients: stochasticity is generally detrimental in fluctuating environments and beneficial only at equal periods of nutrient availability and above a threshold environmental richness. Moreover, depending on the availability and nutritional value of nutrients, nutrient-dependent and stochastic expression are both strategies used to deal with environmental changes. Overall, we comprehensively characterize the interplay between the quality and periodicity of an environment and the resulting optimal deterministic and stochastic regulation strategies of nutrient-catabolizing pathways.

Characterization of a Lactococcus lactis promoter for heterologous protein production

•

Characterization of PTS-IIC, an endogenous constitutive promoter from L. lactis..

•

Cellobiose enhances activity from PTS-IIC promoter.

•

PTS-IIC promoter mediates protein expression in B. subtilis and E coli Nissle 1917.

Constitutively active promoter elements for heterologous protein production in Lactococcus lactis are scarce. Here, the promoter of the PTS-IIC gene cluster from L. lactis NZ3900 is described. This promoter was cloned upstream of an enhanced green fluorescent protein, GFPmut3a, and transformed into L. lactis. Transformants produced up to 13.5 μg of GFPmut3a per milliliter of log phase cells. Addition of cellobiose further increased the production of GFPmut3a by up to two-fold when compared to glucose. Analysis of mutations at two specific positions in the PTS-IIC promoter showed that a ‘T’ to ‘G’ mutation within the −35 element resulted in constitutive expression in glucose, while a ‘C’ at nucleotide 7 in the putative cre site enhanced promoter activity in cellobiose. Finally, this PTS-IIC promoter is capable of mediating protein expression in Bacillus subtilis and Escherichia coli Nissle 1917, suggesting the potential for future biotechnological applications of this element and its derivatives.

Engineered Control of Genetic Variability Reveals Interplay among Quorum Sensing, Feedback Regulation, and Biochemical Noise

Stochastic fluctuations in gene expression trigger both beneficial and harmful consequences for cell behavior. Therefore, achieving a desired mean protein expression level while minimizing noise is of interest in many applications, including robust protein production systems in industrial biotechnology. Here, we consider a synthetic gene circuit combining intracellular negative feedback and cell-to-cell communication based on quorum sensing. Accounting for both intrinsic and extrinsic noise, stochastic simulations allow us to analyze the capability of the circuit to reduce noise strength as a function of its parameters. We obtain mean expression levels and noise strengths for all species under different scenarios, showing good agreement with system-wide available experimental data of protein abundance and noise in Escherichia coli. Our in silico experiments, validated by preliminary in vivo results, reveal significant noise attenuation in gene expression through the interplay between quorum sensing and negative feedback and highlight the differential role that they play in regard to intrinsic and extrinsic noise.

Physical model of protein cluster positioning in growing bacteria

Chemotaxic receptors in bacteria form clusters at cell poles and also laterally, and this clustering plays an important role in signal transduction. These clusters were found to be periodically arranged on the surface of the bacterium Escherichia coli, independent of any known positioning mechanism. In this work we extend a model based on diffusion and aggregation to more realistic geometries and present a means based on "bursty" protein production to distinguish spontaneous positioning from an independently existing positioning mechanism. We also consider the case of isotropic cellular growth and characterize the degree of order arising spontaneously. Our model could also be relevant for other examples of periodically positioned protein clusters in bacteria.

Transcriptional Output Transiently Spikes Upon Mitotic Exit

The pulsatile nature of gene activity has recently emerged as a general property of the transcriptional process. It has been shown that the frequency and amplitude of transcriptional bursts can be subjected to extrinsic regulation. Here we have investigated if these parameters were constant throughout the cell cycle using the single molecule RNA FISH technique. We found evidence of transcriptional spikes upon mitotic exit in three different human cell lines. Recording of cell growth prior to hybridization and immuno-RNA FISH analysis revealed that these spikes were short-lived and subsided before completion of cytokinesis. The transient post-mitotic increase in transcriptional output was found to be the result of cells displaying a higher number of active alleles and/or an increased number of nascent transcripts per active allele, indicating that both the burst fraction and the amplitude of individual bursts can be increased upon mitotic exit. Our results further suggest that distinct regulatory mechanisms are at work shortly after mitotic exit and during the rest of interphase. We speculate that transcriptional spikes are associated with chromatin decondensation, a hallmark of post-mitotic cells that might alter the dynamics of transcriptional regulators and effectors.

Rate-limiting steps in transcription dictate sensitivity to variability in cellular components

Cell-to-cell variability in cellular components generates cell-to-cell diversity in RNA and protein production dynamics. As these components are inherited, this should also cause lineage-to-lineage variability in these dynamics. We conjectured that these effects on transcription are promoter initiation kinetics dependent. To test this, first we used stochastic models to predict that variability in the numbers of molecules involved in upstream processes, such as the intake of inducers from the environment, acts only as a transient source of variability in RNA production numbers, while variability in the numbers of a molecular species controlling transcription of an active promoter acts as a constant source. Next, from single-cell, single-RNA level time-lapse microscopy of independent lineages of Escherichia coli cells, we demonstrate the existence of lineage-to-lineage variability in gene activation times and mean RNA production rates, and that these variabilities differ between promoters and inducers used. Finally, we provide evidence that this can be explained by differences in the kinetics of the rate-limiting steps in transcription between promoters and induction schemes. We conclude that cell-to-cell and consequent lineage-to-lineage variability in RNA and protein numbers are both promoter sequence-dependent and subject to regulation.

What's Luck Got to Do with It: Single Cells, Multiple Fates, and Biological Nondeterminism

The field of single-cell biology has morphed from a philosophical digression at its inception, to a playground for quantitative biologists, to a major area of biomedical research. The last several years have witnessed an explosion of new technologies, allowing us to apply even more of the modern molecular biology toolkit to single cells. Conceptual progress, however, has been comparatively slow. Here, we provide a framework for classifying both the origins of the differences between individual cells and the consequences of those differences. We discuss how the concept of "random" differences is context dependent, and propose that rigorous definitions of inputs and outputs may bring clarity to the discussion. We also categorize ways in which probabilistic behavior may influence cellular function, highlighting studies that point to exciting future directions in the field.

Taking chances and making mistakes: non-genetic phenotypic heterogeneity and its consequences for surviving in dynamic environments

Natural selection has shaped the strategies for survival and growth of microorganisms. The success of microorganisms depends not only on slow evolutionary tuning but also on the ability to adapt to unpredictable changes in their environment. In principle, adaptive strategies range from purely deterministic mechanisms to those that exploit the randomness intrinsic to many cellular and molecular processes. Depending on the environment and selective pressures, particular strategies can lie somewhere along this continuum. In recent years, non-genetic cell-to-cell differences have received a lot of attention, not least because of their potential impact on the ability of microbial populations to survive in dynamic environments. Using several examples, we describe the origins of spontaneous and induced mechanisms of phenotypic adaptation. We identify some of the commonalities of these examples and consider the potential role of chance and constraints in microbial phenotypic adaptation.

Careful accounting of extrinsic noise in protein expression reveals correlations among its sources

In order to grow and replicate, living cells must express a diverse array of proteins, but the process by which proteins are made includes a great deal of inherent randomness. Understanding this randomness-whether it arises from the discrete stochastic nature of chemical reactivity ("intrinsic" noise), or from cell-to-cell variability in the concentrations of molecules involved in gene expression, or from the timings of important cell-cycle events like DNA replication and cell division ("extrinsic" noise)-remains a challenge. In this article we analyze a model of gene expression that accounts for several extrinsic sources of noise, including those associated with chromosomal replication, cell division, and variability in the numbers of RNA polymerase, ribonuclease E, and ribosomes. We then attempt to fit our model to a large proteomics and transcriptomics data set and find that only through the introduction of a few key correlations among the extrinsic noise sources can we accurately recapitulate the experimental data. These include significant correlations between the rate of mRNA degradation (mediated by ribonuclease E) and the rates of both transcription (RNA polymerase) and translation (ribosomes) and, strikingly, an anticorrelation between the transcription and the translation rates themselves.

High-throughput validation of ceRNA regulatory networks

Background

MicroRNAs (miRNAs) play multiple roles in tumor biology. Interestingly, reports from multiple groups suggest that miRNA targets may be coupled through competitive stoichiometric sequestration. Specifically, computational models predicted and experimental assays confirmed that miRNA activity is dependent on miRNA target abundance, and consequently, changes in the abundance of some miRNA targets lead to changes to the regulation and abundance of their other targets. The resulting indirect regulatory influence between miRNA targets resembles competition and has been dubbed competitive endogenous RNA (ceRNA). Recent studies have questioned the physiological relevance of ceRNA interactions, our ability to accurately predict these interactions, and the number of genes that are impacted by ceRNA interactions in specific cellular contexts.

Results

To address these concerns, we reverse engineered ceRNA networks (ceRNETs) in breast and prostate adenocarcinomas using context-specific TCGA profiles, and tested whether ceRNA interactions can predict the effects of RNAi-mediated gene silencing perturbations in PC3 and MCF7 cells._ENREF_22 Our results, based on tests of thousands of inferred ceRNA interactions that are predicted to alter hundreds of cancer genes in each of the two tumor contexts, confirmed statistically significant effects for half of the predicted targets.

Conclusions

Our results suggest that the expression of a significant fraction of cancer genes may be regulated by ceRNA interactions in each of the two tumor contexts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3790-7) contains supplementary material, which is available to authorized users.

Effect of transcription factor resource sharing on gene expression noise

Gene expression is intrinsically a stochastic (noisy) process with important implications for cellular functions. Deciphering the underlying mechanisms of gene expression noise remains one of the key challenges of regulatory biology. Theoretical models of transcription often incorporate the kinetics of how transcription factors (TFs) interact with a single promoter to impact gene expression noise. However, inside single cells multiple identical gene copies as well as additional binding sites can compete for a limiting pool of TFs. Here we develop a simple kinetic model of transcription, which explicitly incorporates this interplay between TF copy number and its binding sites. We show that TF sharing enhances noise in mRNA distribution across an isogenic population of cells. Moreover, when a single gene copy shares it’s TFs with multiple competitor sites, the mRNA variance as a function of the mean remains unaltered by their presence. Hence, all the data for variance as a function of mean expression collapse onto a single master curve independent of the strength and number of competitor sites. However, this result does not hold true when the competition stems from multiple copies of the same gene. Therefore, although previous studies showed that the mean expression follows a universal master curve, our findings suggest that different scenarios of competition bear distinct signatures at the level of variance. Intriguingly, the introduction of competitor sites can transform a unimodal mRNA distribution into a multimodal distribution. These results demonstrate the impact of limited availability of TF resource on the regulation of noise in gene expression.

Genetically identical cells, even when they are exposed to the same environmental conditions, display incredible diversity. Gene expression noise is attributed to be a key source of this phenotypic diversity. Transcriptional dynamics is a dominant source of expression noise. Although scores of theoretical and experimental studies have explored how noise is regulated at the level of transcription, most of them focus on the gene specific, cis regulatory elements, such as the number of transcription factor (TF) binding sites, their binding strength, etc. However, how the global properties of transcription, such as the limited availability of TFs impact noise in gene expression remains rather elusive. Here we build a theoretical model that incorporates the effect of limiting TF pool on gene expression noise. We find that competition between genes for TFs leads to enhanced variability in mRNA copy number across an isogenic population. Moreover, for gene copies sharing TFs with other competitor sites, mRNA variance as a function of the mean shows distinct imprints for one gene copy and multiple gene copies respectively. This stands in sharp contrast to the universal behavior found in mean expression irrespective of the different scenarios of competition. An interesting feature of competition is that introduction of competitor sites can transform a unimodal mRNA distribution into a multimodal distribution, which could lead to phenotypic variability.

Can We Predict Gene Expression by Understanding Proximal Promoter Architecture?

We review computational predictions of expression from the promoter architecture - the set of transcription factors that can bind the proximal promoter. We focus on spatial expression patterns in animals with complex body plans and many distinct tissue types. This field is ripe for change as functional genomics datasets accumulate for both expression and protein-DNA interactions. While there has been some success in predicting the breadth of expression (i.e., the fraction of tissue types a gene is expressed in), predicting tissue specificity remains challenging. We discuss how progress can be achieved through either machine learning or complementary combinatorial data mining. The likely impact of single-cell expression data is considered. Finally, we discuss the design of artificial promoters as a practical application.

Determining the Limitations and Benefits of Noise in Gene Regulation and Signal Transduction through Single Cell, Microscopy-Based Analysis

Stochastic fluctuations, termed "noise," in the level of biological molecules can greatly impact cellular functions. While biological noise can sometimes be detrimental, recent studies have provided an increasing number of examples in which biological noise can be functionally beneficial. Rather than provide an exhaustive review of the growing literature in this field, in this review, we focus on single-cell studies based on quantitative microscopy that have generated a deeper understanding of the sources, characteristics, limitations, and benefits of biological noise. Specifically, we highlight studies showing how noise can help coordinate the expression of multiple downstream target genes, impact the channel capacity of signaling networks, and interact synergistically with oscillatory dynamics to enhance the sensitivity of signal processing. We conclude with a discussion of current challenges and future opportunities.

Non-genetic diversity modulates population performance

Biological functions are typically performed by groups of cells that express predominantly the same genes, yet display a continuum of phenotypes. While it is known how one genotype can generate such non-genetic diversity, it remains unclear how different phenotypes contribute to the performance of biological function at the population level. We developed a microfluidic device to simultaneously measure the phenotype and chemotactic performance of tens of thousands of individual, freely swimming Escherichia coli as they climbed a gradient of attractant. We discovered that spatial structure spontaneously emerged from initially well-mixed wild-type populations due to non-genetic diversity. By manipulating the expression of key chemotaxis proteins, we established a causal relationship between protein expression, non-genetic diversity, and performance that was theoretically predicted. This approach generated a complete phenotype-to-performance map, in which we found a nonlinear regime. We used this map to demonstrate how changing the shape of a phenotypic distribution can have as large of an effect on collective performance as changing the mean phenotype, suggesting that selection could act on both during the process of adaptation.

Transcriptional bursting is intrinsically caused by interplay between RNA polymerases on DNA

Cell-to-cell variability plays a critical role in cellular responses and decision-making in a population, and transcriptional bursting has been broadly studied by experimental and theoretical approaches as the potential source of cell-to-cell variability. Although molecular mechanisms of transcriptional bursting have been proposed, there is little consensus. An unsolved key question is whether transcriptional bursting is intertwined with many transcriptional regulatory factors or is an intrinsic characteristic of RNA polymerase on DNA. Here we design an in vitro single-molecule measurement system to analyse the kinetics of transcriptional bursting. The results indicate that transcriptional bursting is caused by interplay between RNA polymerases on DNA. The kinetics of in vitro transcriptional bursting is quantitatively consistent with the gene-nonspecific kinetics previously observed in noisy gene expression in vivo. Our kinetic analysis based on a cellular automaton model confirms that arrest and rescue by trailing RNA polymerase intrinsically causes transcriptional bursting.

Transcriptional bursting is a potential source of cell-to-cell variability but the molecular mechanisms are unclear. Here the authors use single molecule imaging to analyse the kinetics of bursting on DNA and observe that bursting is an intrinsic property of RNA polymerases on DNA.

Mechanical bounds to transcriptional noise

Over the past several decades it has been increasingly recognized that stochastic processes play a central role in transcription. Although many stochastic effects have been explained, the source of transcriptional bursting (one of the most well-known sources of stochasticity) has continued to evade understanding. Recent results have pointed to mechanical feedback as the source of transcriptional bursting, but a reconciliation of this perspective with preexisting views of transcriptional regulation is lacking. In this article, we present a simple phenomenological model that is able to incorporate the traditional view of gene expression within a framework with mechanical limits to transcription. By introducing a simple competition between mechanical arrest and relaxation copy number probability distributions collapse onto a shared universal curve under shifting and rescaling and a lower limit of intrinsic noise for any mean expression level is found.

Temperature-Dependent Model of Multi-step Transcription Initiation in Escherichia coli Based on Live Single-Cell Measurements

Transcription kinetics is limited by its initiation steps, which differ between promoters and with intra- and extracellular conditions. Regulation of these steps allows tuning both the rate and stochasticity of RNA production. We used time-lapse, single-RNA microscopy measurements in live Escherichia coli to study how the rate-limiting steps in initiation of the P_lac/ara-1 promoter change with temperature and induction scheme. For this, we compared detailed stochastic models fit to the empirical data in maximum likelihood sense using statistical methods. Using this analysis, we found that temperature affects the rate limiting steps unequally, as nonlinear changes in the closed complex formation suffice to explain the differences in transcription dynamics between conditions. Meanwhile, a similar analysis of the P_tetA promoter revealed that it has a different rate limiting step configuration, with temperature regulating different steps. Finally, we used the derived models to explore a possible cause for why the identified steps are preferred as the main cause for behavior modifications with temperature: we find that transcription dynamics is either insensitive or responds reciprocally to changes in the other steps. Our results suggests that different promoters employ different rate limiting step patterns that control not only their rate and variability, but also their sensitivity to environmental changes.

Temperature affects the behavior of cells, such as their growth rate. However, it is not well understood how these changes result from the changes at the single molecule level. We observed the production of individual RNA molecules in live cells under a wide range of temperatures. This allowed us to determine not only how fast they are produced, but also how much variability there is in this process. Next, we fit a stochastic model to the data to identify which rate-limiting steps during RNA production are responsible for the observed differences between conditions. We found that genes differ in how their RNA production is limited by different steps and in how these are affected by the temperature, which explains why different genes respond differently to temperature fluctuations.

Ribosome biogenesis in replicating cells: Integration of experiment and theory

Ribosomes-the primary macromolecular machines responsible for translating the genetic code into proteins-are complexes of precisely folded RNA and proteins. The ways in which their production and assembly are managed by the living cell is of deep biological importance. Here we extend a recent spatially resolved whole-cell model of ribosome biogenesis in a fixed volume [Earnest et al., Biophys J 2015, 109, 1117-1135] to include the effects of growth, DNA replication, and cell division. All biological processes are described in terms of reaction-diffusion master equations and solved stochastically using the Lattice Microbes simulation software. In order to determine the replication parameters, we construct and analyze a series of Escherichia coli strains with fluorescently labeled genes distributed evenly throughout their chromosomes. By measuring these cells' lengths and number of gene copies at the single-cell level, we could fit a statistical model of the initiation and duration of chromosome replication. We found that for our slow-growing (120 min doubling time) E. coli cells, replication was initiated 42 min into the cell cycle and completed after an additional 42 min. While simulations of the biogenesis model produce the correct ribosome and mRNA counts over the cell cycle, the kinetic parameters for transcription and degradation are lower than anticipated from a recent analytical time dependent model of in vivo mRNA production. Describing expression in terms of a simple chemical master equation, we show that the discrepancies are due to the lack of nonribosomal genes in the extended biogenesis model which effects the competition of mRNA for ribosome binding, and suggest corrections to parameters to be used in the whole-cell model when modeling expression of the entire transcriptome. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 735-751, 2016.

Stochastic Kinetics of Nascent RNA

The stochastic kinetics of transcription is typically inferred from the distribution of RNA numbers in individual cells. However, cellular RNA reflects additional processes downstream of transcription, hampering this analysis. In contrast, nascent (actively transcribed) RNA closely reflects the kinetics of transcription. We present a theoretical model for the stochastic kinetics of nascent RNA, which we solve to obtain the probability distribution of nascent RNA per gene. The model allows us to evaluate the kinetic parameters of transcription from single-cell measurements of nascent RNA. The model also predicts surprising discontinuities in the distribution of nascent RNA, a feature which we verify experimentally.

Fluctuations in reactive networks subject to extrinsic noise studied in the framework of the chemical Langevin equation

Theoretical and experimental studies have shown that the fluctuations of in vivo systems break the fluctuation-dissipation theorem. One can thus ask what information is contained in the correlation functions of protein concentrations and how they relate to the response of the reactive network to a perturbation. Answers to these questions are of prime importance to extract meaningful parameters from the in vivo fluorescence correlation spectroscopy data. In this paper we study the fluctuations of the concentration of a reactive species involved in a cyclic network that is in a nonequilibrium steady state perturbed by a noisy force, taking into account both the breaking of detailed balance and extrinsic noises. Using a generic model for the network and the extrinsic noise, we derive a chemical Langevin equation that describes the dynamics of the system, we determine the expressions of the correlation functions of the concentrations, and we estimate the deviation of the fluctuation-dissipation theorem and the range of parameters in which an effective temperature can be defined.

Scaling-Up Signaling Pathway Analysis

A new technique for simultaneously measuring the activities of many signaling pathways unravels interconnected signaling networks.

Transcript degradation and noise of small RNA-controlled genes in a switch activated network in Escherichia coli

Post-transcriptional regulatory processes may change transcript levels and affect cell-to-cell variability or noise. We study small-RNA downregulation to elucidate its effects on noise in the iron homeostasis network of Escherichia coli. In this network, the small-RNA RyhB undergoes stoichiometric degradation with the transcripts of target genes in response to iron stress. Using single-molecule fluorescence in situ hybridization, we measured transcript numbers of the RyhB-regulated genes sodB and fumA in individual cells as a function of iron deprivation. We observed a monotonic increase of noise with iron stress but no evidence of theoretically predicted, enhanced stoichiometric fluctuations in transcript numbers, nor of bistable behavior in transcript distributions. Direct detection of RyhB in individual cells shows that its noise is much smaller than that of these two targets, when RyhB production is significant. A generalized two-state model of bursty transcription that neglects RyhB fluctuations describes quantitatively the dependence of noise and transcript distributions on iron deprivation, enabling extraction of in vivo RyhB-mediated transcript degradation rates. The transcripts’ threshold-linear behavior indicates that the effective in vivo interaction strength between RyhB and its two target transcripts is comparable. Strikingly, the bacterial cell response exhibits Fur-dependent, switch-like activation instead of a graded response to iron deprivation.

Regulation of cell-to-cell variability in divergent gene expression

Cell-to-cell variability (noise) is an important feature of gene expression that impacts cell fitness and development. The regulatory mechanism of this variability is not fully understood. Here we investigate the effect on gene expression noise in divergent gene pairs (DGPs). We generated reporters driven by divergent promoters, rearranged their gene order, and probed their expressions using time-lapse fluorescence microscopy and single-molecule fluorescence in situ hybridization (smFISH). We show that two genes in a co-regulated DGP have higher expression covariance compared with the separate, tandem and convergent configurations, and this higher covariance is caused by more synchronized firing of the divergent transcriptions. For differentially regulated DGPs, the regulatory signal of one gene can stochastically ‘leak' to the other, causing increased gene expression noise. We propose that the DGPs' function in limiting or promoting gene expression noise may enhance or compromise cell fitness, providing an explanation for the conservation pattern of DGPs.

Gene expression noise affects cell fitness and development. Here, Yan et al. show that co-regulated divergent gene pairs (DGPs) suppress uncorrelated gene expression noise due to more synchronized transcription firing, and differentially regulated DGPs enhance gene expression noise due to transcription leakage.

Measurement of gene regulation in individual cells reveals rapid switching between promoter states

In vivo mapping of transcription-factor binding to the transcriptional output of the regulated gene is hindered by probabilistic promoter occupancy, the presence of multiple gene copies, and cell-to-cell variability. We demonstrate how to overcome these obstacles in the lysogeny maintenance promoter of bacteriophage lambda, P(RM). We simultaneously measured the concentration of the lambda repressor CI and the number of messenger RNAs (mRNAs) from P(RM) in individual Escherichia coli cells, and used a theoretical model to identify the stochastic activity corresponding to different CI binding configurations. We found that switching between promoter configurations is faster than mRNA lifetime and that individual gene copies within the same cell act independently. The simultaneous quantification of transcription factor and promoter activity, followed by stochastic theoretical analysis, provides a tool that can be applied to other genetic circuits.

Using synthetic biology to make cells tomorrow's test tubes

The main tenet of physical biology is that biological phenomena can be subject to the same quantitative and predictive understanding that physics has afforded in the context of inanimate matter. However, the inherent complexity of many of these biological processes often leads to the derivation of complex theoretical descriptions containing a plethora of unknown parameters. Such complex descriptions pose a conceptual challenge to the establishment of a solid basis for predictive biology. In this article, we present various exciting examples of how synthetic biology can be used to simplify biological systems and distill these phenomena down to their essential features as a means to enable their theoretical description. Here, synthetic biology goes beyond previous efforts to engineer nature and becomes a tool to bend nature to understand it. We discuss various recent and classic experiments featuring applications of this synthetic approach to the elucidation of problems ranging from bacteriophage infection, to transcriptional regulation in bacteria and in developing embryos, to evolution. In all of these examples, synthetic biology provides the opportunity to turn cells into the equivalent of a test tube, where biological phenomena can be reconstituted and our theoretical understanding put to test with the same ease that these same phenomena can be studied in the in vitro setting.

A biophysical model of supercoiling dependent transcription predicts a structural aspect to gene regulation

Background

Transcription in Escherichia coli generates positive supercoiling in the DNA, which is relieved by the enzymatic activity of gyrase. Recently published experimental evidence suggests that transcription initiation and elongation are inhibited by the buildup of positive supercoiling. It has therefore been proposed that intermittent binding of gyrase plays a role in transcriptional bursting. Considering that transcription is one of the most fundamental cellular processes, it is desirable to be able to account for the buildup and release of positive supercoiling in models of transcription.

Results

Here we present a detailed biophysical model of gene expression that incorporates the effects of supercoiling due to transcription. By directly linking the amount of positive supercoiling to the rate of transcription, the model predicts that highly transcribed genes’ mRNA distributions should substantially deviate from Poisson distributions, with enhanced density at low mRNA copy numbers. Additionally, the model predicts a high degree of correlation between expression levels of genes inside the same supercoiling domain.

Conclusions

Our model, incorporating the supercoiling state of the gene, makes specific predictions that differ from previous models of gene expression. Genes in the same supercoiling domain influence the expression level of neighboring genes. Such structurally dependent regulation predicts correlations between genes in the same supercoiling domain. The topology of the chromosome therefore creates a higher level of gene regulation, which has broad implications for understanding the evolution and organization of bacterial genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13628-016-0027-0) contains supplementary material, which is available to authorized users.

Single-cell analysis of transcription kinetics across the cell cycle

Transcription is a highly stochastic process. To infer transcription kinetics for a gene-of-interest, researchers commonly compare the distribution of mRNA copy-number to the prediction of a theoretical model. However, the reliability of this procedure is limited because the measured mRNA numbers represent integration over the mRNA lifetime, contribution from multiple gene copies, and mixing of cells from different cell-cycle phases. We address these limitations by simultaneously quantifying nascent and mature mRNA in individual cells, and incorporating cell-cycle effects in the analysis of mRNA statistics. We demonstrate our approach on Oct4 and Nanog in mouse embryonic stem cells. Both genes follow similar two-state kinetics. However, Nanog exhibits slower ON/OFF switching, resulting in increased cell-to-cell variability in mRNA levels. Early in the cell cycle, the two copies of each gene exhibit independent activity. After gene replication, the probability of each gene copy to be active diminishes, resulting in dosage compensation.

Effects of DNA replication on mRNA noise

There are several sources of fluctuations in gene expression. Here we study the effects of time-dependent DNA replication, itself a tightly controlled process, on noise in mRNA levels. Stochastic simulations of constitutive and regulated gene expression are used to analyze the time-averaged mean and variation in each case. The simulations demonstrate that to capture mRNA distributions correctly, chromosome replication must be realistically modeled. Slow relaxation of mRNA from the low copy number steady state before gene replication to the high steady state after replication is set by the transcript's half-life and contributes significantly to the shape of the mRNA distribution. Consequently both the intrinsic kinetics and the gene location play an important role in accounting for the mRNA average and variance. Exact analytic expressions for moments of the mRNA distributions that depend on the DNA copy number, gene location, cell doubling time, and the rates of transcription and degradation are derived for the case of constitutive expression and subsequently extended to provide approximate corrections for regulated expression and RNA polymerase variability. Comparisons of the simulated models and analytical expressions to experimentally measured mRNA distributions show that they better capture the physics of the system than previous theories.

Mitochondria and the non-genetic origins of cell-to-cell variability: More is different

Gene expression activity is heterogeneous in a population of isogenic cells. Identifying the molecular basis of this variability will improve our understanding of phenomena like tumor resistance to drugs, virus infection, or cell fate choice. The complexity of the molecular steps and machines involved in transcription and translation could introduce sources of randomness at many levels, but a common constraint to most of these processes is its energy dependence. In eukaryotic cells, most of this energy is provided by mitochondria. A clonal population of cells may show a large variability in the number and functionality of mitochondria. Here, we discuss how differences in the mitochondrial content of each cell contribute to heterogeneity in gene products. Changes in the amount of mitochondria can also entail drastic alterations of a cell's gene expression program, which ultimately leads to phenotypic diversity. Also watch the Video Abstract.

Dynamics of simple gene-network motifs subject to extrinsic fluctuations

Cellular processes do not follow deterministic rules; even in identical environments genetically identical cells can make random choices leading to different phenotypes. This randomness originates from fluctuations present in the biomolecular interaction networks. Most previous work has been focused on the intrinsic noise (IN) of these networks. Yet, especially for high-copy-number biomolecules, extrinsic or environmental noise (EN) has been experimentally shown to dominate the variation. Here, we develop an analytical formalism that allows for calculation of the effect of EN on gene-expression motifs. We introduce a method for modeling bounded EN as an auxiliary species in the master equation. The method is fully generic and is not limited to systems with small EN magnitudes. We focus our study on motifs that can be viewed as the building blocks of genetic switches: a nonregulated gene, a self-inhibiting gene, and a self-promoting gene. The role of the EN properties (magnitude, correlation time, and distribution) on the statistics of interest are systematically investigated, and the effect of fluctuations in different reaction rates is compared. Due to its analytical nature, our formalism can be used to quantify the effect of EN on the dynamics of biochemical networks and can also be used to improve the interpretation of data from single-cell gene-expression experiments.

Universal protein distributions in a model of cell growth and division

Protein distributions measured under a broad set of conditions in bacteria and yeast were shown to exhibit a common skewed shape, with variances depending quadratically on means. For bacteria these properties were reproduced by temporal measurements of protein content, showing accumulation and division across generations. Here we present a stochastic growth-and-division model with feedback which captures these observed properties. The limiting copy number distribution is calculated exactly, and a single parameter is found to determine the distribution shape and the variance-to-mean relation. Estimating this parameter from bacterial temporal data reproduces the measured distribution shape with high accuracy and leads to predictions for future experiments.

Phenotypic heterogeneity and temporal expression of the capsular polysaccharide in Staphylococcus aureus

Bacteria respond to ever-changing environments through several adaptive strategies. This includes mechanisms leading to a high degree of phenotypic variability within a genetically homogeneous population. In Staphylococcus aureus, the capsular polysaccharide (CP) protects against phagocytosis, but also impedes adherence to endothelial cells and/or matrix proteins. We analysed the regulation of core biosynthesis genes (capA-P) necessary for CP synthesis using single-cell assays (immunofluorescence and promoter-activity). In persistent human carriers, we found a distinct subpopulation of nasal S. aureus to be CP positive. In vitro, cap expression is also heterogeneous and strongly growth-phase dependent. We asked whether this peculiar expression pattern (earlyOff/lateHeterogen) is orchestrated by the quorum system Agr. We show that the Agr-driven effector molecule RNAIII promotes cap expression largely via inactivation of the repressor Rot. High NaCl, deletion of CodY or Sae also resulted in higher cap expression but did not change the earlyOFF/lateHeterogen expression pattern. Activity of the quorum system itself is largely homogenous and does not account for the observed heterogeneity of cap expression or the strictly growth phase dependent expression. Our findings are in contrast to the prevailing view that quorum sensing is the main driving force for virulence gene expression when bacterial cell densities increase.

Noise in gene expression is coupled to growth rate

Genetically identical cells exposed to the same environment display variability in gene expression (noise), with important consequences for the fidelity of cellular regulation and biological function. Although population average gene expression is tightly coupled to growth rate, the effects of changes in environmental conditions on expression variability are not known. Here, we measure the single-cell expression distributions of approximately 900 Saccharomyces cerevisiae promoters across four environmental conditions using flow cytometry, and find that gene expression noise is tightly coupled to the environment and is generally higher at lower growth rates. Nutrient-poor conditions, which support lower growth rates, display elevated levels of noise for most promoters, regardless of their specific expression values. We present a simple model of noise in expression that results from having an asynchronous population, with cells at different cell-cycle stages, and with different partitioning of the cells between the stages at different growth rates. This model predicts non-monotonic global changes in noise at different growth rates as well as overall higher variability in expression for cell-cycle–regulated genes in all conditions. The consistency between this model and our data, as well as with noise measurements of cells growing in a chemostat at well-defined growth rates, suggests that cell-cycle heterogeneity is a major contributor to gene expression noise. Finally, we identify gene and promoter features that play a role in gene expression noise across conditions. Our results show the existence of growth-related global changes in gene expression noise and suggest their potential phenotypic implications.

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.

Expression noise facilitates the evolution of gene regulation

Although it is often tacitly assumed that gene regulatory interactions are finely tuned, how accurate gene regulation could evolve from a state without regulation is unclear. Moreover, gene expression noise would seem to impede the evolution of accurate gene regulation, and previous investigations have provided circumstantial evidence that natural selection has acted to lower noise levels. By evolving synthetic Escherichia coli promoters de novo, we here show that, contrary to expectations, promoters exhibit low noise by default. Instead, selection must have acted to increase the noise levels of highly regulated E. coli promoters. We present a general theory of the interplay between gene expression noise and gene regulation that explains these observations. The theory shows that propagation of expression noise from regulators to their targets is not an unwanted side-effect of regulation, but rather acts as a rudimentary form of regulation that facilitates the evolution of more accurate regulation.