Gene Expression Noise Research Articles

Orthogonal control of mean and variability of endogenous genes in a human cell line

Stochastic fluctuations at the transcriptional level contribute to isogenic cell-to-cell heterogeneity in mammalian cell populations. However, we still have no clear understanding of the repercussions of this heterogeneity, given the lack of tools to independently control mean expression and variability of a gene. Here, we engineer a synthetic circuit to modulate mean expression and heterogeneity of transgenes and endogenous human genes. The circuit, a Tunable Noise Rheostat (TuNR), consists of a transcriptional cascade of two inducible transcriptional activators, where the output mean and variance can be modulated by two orthogonal small molecule inputs. In this fashion, different combinations of the inputs can achieve the same mean but with different population variability. With TuNR, we achieve low basal expression, over 1000-fold expression of a transgene product, and up to 7-fold induction of the endogenous gene NGFR. Importantly, for the same mean expression level, we are able to establish varying degrees of heterogeneity in expression within an isogenic population, thereby decoupling gene expression noise from its mean. TuNR is therefore a modular tool that can be used in mammalian cells to enable direct interrogation of the implications of cell-to-cell variability.

MicroRNA Based Feedforward Control of Intrinsic Gene Expression Noise.

Intrinsic noise, which arises in gene expression at low copy numbers, can be controlled by diverse regulatory motifs, including feedforward loops. Here, we study an example of a feedforward control system based on the interaction between an mRNA molecule and an antagonistic microRNA molecule encoded by the same gene, aiming to quantify the variability (or noise) in molecular copy numbers. Using linear noise approximation, we show that the mRNA noise is sub-Poissonian in case of non-bursty transcription, and exhibits a nonmonotonic response both to the species natural lifetime ratio and to the strength of the antagonistic interaction. Additionally, we use the Chemical Reaction Network Theory to prove that the mRNA copy number distribution is Poissonian in the absence of spontaneous mRNA decay channel. In case of transcriptional bursts, we show that feedforward control can attenuate the super-Poissonian gene-expression noise that is due to bursting, and that the effect is more considerable at the protein than at the mRNA level. Our results indicate that the strong coupling between mRNA and microRNA in the sense of burst stoichiometry and also of timing of production events renders the microRNA based feedforward motif an effective mechanism for the control of gene expression noise.

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

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

Gene Circuit Explorer (GeneEx): an interactive web-app for visualizing, simulating and analyzing gene regulatory circuits.

GeneEx is an interactive web-app that uses an ODE-based mathematical modeling approach to simulate, visualize and analyze gene regulatory circuits (GRCs) for an explicit kinetic parameter set or for a large ensemble of random parameter sets. GeneEx offers users the freedom to modify many aspects of the simulation such as the parameter ranges, the levels of gene expression noise and the GRC network topology itself. This degree of flexibility allows users to explore a variety of hypotheses by providing insight into the number and stability of attractors for a given GRC. Moreover, users have the option to upload, and subsequently compare, experimental gene expression data to simulated data generated from the analysis of a built or uploaded custom circuit. Finally, GeneEx offers a curated database that contains circuit motifs and known biological GRCs to facilitate further inquiry into these. Overall, GeneEx enables users to investigate the effects of parameter variation, stochasticity and/or topological changes on gene expression for GRCs using a systems-biology approach. GeneEx is available at https://geneex.jax.org. This web-app is released under the MIT license and is free and open to all users and there is no mandatory login requirement. Supplementary data are available at Bioinformatics online.

Copy Number Heterogeneity in the Virulence Plasmid of Salmonella enterica.

Quantitative PCR analysis shows that the virulence plasmid of Salmonella enterica serovar Typhimurium (pSLT) is a low-copy-number plasmid, with 1–2 copies per chromosome. However, fluorescence microscopy observation of pSLT labeled with a lacO fluorescent tag reveals cell-to-cell differences in the number of foci, which ranges from 1 to 8. As each focus must correspond to ≥1 plasmid copy, the number of foci can be expected to indicate the minimal number of pSLT copies per cell. A correlation is found between the number of foci and the bacterial cell volume. In contrast, heterogeneity in the number of foci appears to be independent of the cell volume and may have stochastic origin. As a consequence of copy number heterogeneity, expression of a pSLT-bone reporter gene shows high levels of cell-to-cell variation, especially in actively dividing cultures. These observations support the notion that low-copy-number plasmids can be a source of gene expression noise in bacterial populations.

Subpopulations of sensorless bacteria drive fitness in fluctuating environments.

Populations of bacteria often undergo a lag in growth when switching conditions. Because growth lags can be large compared to typical doubling times, variations in growth lag are an important but often overlooked component of bacterial fitness in fluctuating environments. We here explore how growth lag variation is determined for the archetypical switch from glucose to lactose as a carbon source in Escherichia coli. First, we show that single-cell lags are bimodally distributed and controlled by a single-molecule trigger. That is, gene expression noise causes the population before the switch to divide into subpopulations with zero and nonzero lac operon expression. While "sensorless" cells with zero preexisting lac expression at the switch have long lags because they are unable to sense the lactose signal, any nonzero lac operon expression suffices to ensure a short lag. Second, we show that the growth lag at the population level depends crucially on the fraction of sensorless cells and that this fraction in turn depends sensitively on the growth condition before the switch. Consequently, even small changes in basal expression can significantly affect the fraction of sensorless cells, thereby population lags and fitness under switching conditions, and may thus be subject to significant natural selection. Indeed, we show that condition-dependent population lags vary across wild E. coli isolates. Since many sensory genes are naturally low expressed in conditions where their inducer is not present, bimodal responses due to subpopulations of sensorless cells may be a general mechanism inducing phenotypic heterogeneity and controlling population lags in switching environments. This mechanism also illustrates how gene expression noise can turn even a simple sensory gene circuit into a bet hedging module and underlines the profound role of gene expression noise in regulatory responses.

Nonequilibrium models of optimal enhancer function

In prokaryotes, thermodynamic models of gene regulation provide a highly quantitative mapping from promoter sequences to gene-expression levels that is compatible with in vivo and in vitro biophysical measurements. Such concordance has not been achieved for models of enhancer function in eukaryotes. In equilibrium models, it is difficult to reconcile the reported short transcription factor (TF) residence times on the DNA with the high specificity of regulation. In nonequilibrium models, progress is difficult due to an explosion in the number of parameters. Here, we navigate this complexity by looking for minimal nonequilibrium enhancer models that yield desired regulatory phenotypes: low TF residence time, high specificity, and tunable cooperativity. We find that a single extra parameter, interpretable as the "linking rate," by which bound TFs interact with Mediator components, enables our models to escape equilibrium bounds and access optimal regulatory phenotypes, while remaining consistent with the reported phenomenology and simple enough to be inferred from upcoming experiments. We further find that high specificity in nonequilibrium models is in a trade-off with gene-expression noise, predicting bursty dynamics-an experimentally observed hallmark of eukaryotic transcription. By drastically reducing the vast parameter space of nonequilibrium enhancer models to a much smaller subspace that optimally realizes biological function, we deliver a rich class of models that could be tractably inferred from data in the near future.

Translational control of gene expression noise and its relationship to ageing in yeast

Gene expression noise influences organism evolution and fitness but is poorly understood. There is increasing evidence that the functional roles of components of the translation machinery influence noise intensity. In addition, modulation of the activities of at least some of these same components affects the replicative lifespan of a broad spectrum of organisms. In a novel comparative approach, we modulate the activities of the translation initiation factors eIFG1 and eIF4G2, both of which are involved in the process of recruiting ribosomal 43S preinitiation complexes to the 5' end of eukaryotic mRNAs. We show that tagging of the cell wall using a fluorescent dye allows us to follow gene expression noise as different yeast strains progress through successive cycles of replicative ageing. This procedure reveals a relationship between global protein synthesis rate and gene expression noise (cell-to-cell heterogeneity), which is accompanied by a parallel correlation between gene expression noise and the replicative age of mother cells. An alternative approach, based on microfluidics, confirms the interdependence between protein synthesis rate, gene expression noise and ageing. We additionally show that it is important to characterize the influence of the design of the microfluidic device on the nutritional state of the cells during such experiments. Analysis of the noise data derived from flow cytometry and fluorescence microscopy measurements indicates that both the intrinsic and the extrinsic noise components increase as a function of ageing.

Gene expression noise can promote the fixation of beneficial mutations in fluctuating environments.

Nongenetic phenotypic variation can either speed up or slow down adaptive evolution. We show that it can speed up evolution in environments where available carbon and energy sources change over time. To this end, we use an experimentally validated model of Escherichia coli growth on two alternative carbon sources, glucose and acetate. On the superior carbon source (glucose), all cells achieve high growth rates, while on the inferior carbon source (acetate) only a small fraction of the population manages to initiate growth. Consequently, populations experience a bottleneck when the environment changes from the superior to the inferior carbon source. Growth on the inferior carbon source depends on a circuit under the control of a transcription factor that is repressed in the presence of the superior carbon source. We show that noise in the expression of this transcription factor can increase the probability that cells start growing on the inferior carbon source. In doing so, it can decrease the severity of the bottleneck and increase mean population fitness whenever this fitness is low. A modest amount of noise can also enhance the fitness effects of a beneficial allele that increases the fraction of a population initiating growth on acetate. Additionally, noise can protect this allele from extinction, accelerate its spread, and increase its likelihood of going to fixation. Central to the adaptation-enhancing principle we identify is the ability of noise to mitigate population bottlenecks, particularly in environments that fluctuate periodically. Because such bottlenecks are frequent in fluctuating environments, and because periodically fluctuating environments themselves are common, this principle may apply to a broad range of environments and organisms.

Using fluorescence flow cytometry data for single-cell gene expression analysis in bacteria.

Fluorescence flow cytometry is increasingly being used to quantify single-cell expression distributions in bacteria in high-throughput. However, there has been no systematic investigation into the best practices for quantitative analysis of such data, what systematic biases exist, and what accuracy and sensitivity can be obtained. We investigate these issues by measuring the same E. coli strains carrying fluorescent reporters using both flow cytometry and microscopic setups and systematically comparing the resulting single-cell expression distributions. Using these results, we develop methods for rigorous quantitative inference of single-cell expression distributions from fluorescence flow cytometry data. First, we present a Bayesian mixture model to separate debris from viable cells using all scattering signals. Second, we show that cytometry measurements of fluorescence are substantially affected by autofluorescence and shot noise, which can be mistaken for intrinsic noise in gene expression, and present methods to correct for these using calibration measurements. Finally, we show that because forward- and side-scatter signals scale non-linearly with cell size, and are also affected by a substantial shot noise component that cannot be easily calibrated unless independent measurements of cell size are available, it is not possible to accurately estimate the variability in the sizes of individual cells using flow cytometry measurements alone. To aid other researchers with quantitative analysis of flow cytometry expression data in bacteria, we distribute E-Flow, an open-source R package that implements our methods for filtering debris and for estimating true biological expression means and variances from the fluorescence signal. The package is available at https://github.com/vanNimwegenLab/E-Flow.

Extrinsic noise prevents the independent tuning of gene expression noise and protein mean abundance in bacteria.

It is generally accepted that prokaryotes can tune gene expression noise independently of protein mean abundance by varying the relative levels of transcription and translation. Here, we address this question quantitatively, using a custom-made library of 40 Bacillus subtilis strains expressing a fluorescent protein under the control of different transcription and translation control elements. We quantify noise and mean protein abundance by fluorescence microscopy and show that for most of the natural transcription range of B. subtilis, expression noise is equally sensitive to variations in the transcription or translation rate because of the prevalence of extrinsic noise. In agreement, analysis of whole-genome transcriptomic and proteomic datasets suggests that noise optimization through transcription and translation tuning during evolution may only occur in a regime of weak transcription. Therefore, independent control of mean abundance and noise can rarely be achieved, which has strong implications for both genome evolution and biological engineering.

Profiling the Non-genetic Origins of Cancer Drug Resistance with a Single-Cell Functional Genomics Approach Using Predictive Cell Dynamics.

Non-genetic heterogeneity observed in clonal cell populations is an immediate cause of drug resistance that remains challenging to profile because of its transient nature. Here, we coupled three single-cell technologies to link the predicted drug response of a cell to its own genome-wide transcriptomic profile. As a proof of principle, we analyzed the response to tumor-necrosis-factor-related apoptosis-inducing ligand (TRAIL) in HeLa cells to demonstrate that cell dynamics can discriminate the transient transcriptional states at the origin of cell decisions such as sensitivity and resistance. Our same-cell approach, named fate-seq, can reveal the molecular factors regulating the efficacy of a drug in clonal cells, providing therapeutic targets of non-genetic drug resistance otherwise confounded in gene expression noise. A record of this paper's transparent peer review process is included in the Supplemental Information.

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

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

Uncoupling gene expression noise along the central dogma using genome engineered human cell lines

Eukaryotic protein synthesis is an inherently stochastic process. This stochasticity stems not only from variations in cell content between cells but also from thermodynamic fluctuations in a single cell. Ultimately, these inherently stochastic processes manifest as noise in gene expression, where even genetically identical cells in the same environment exhibit variation in their protein abundances. In order to elucidate the underlying sources that contribute to gene expression noise, we quantify the contribution of each step within the process of protein synthesis along the central dogma. We uncouple gene expression at the transcriptional, translational, and post-translational level using custom engineered circuits stably integrated in human cells using CRISPR. We provide a generalized framework to approximate intrinsic and extrinsic noise in a population of cells expressing an unbalanced two-reporter system. Our decomposition shows that the majority of intrinsic fluctuations stem from transcription and that coupling the two genes along the central dogma forces the fluctuations to propagate and accumulate along the same path, resulting in increased observed global correlation between the products.

Contributions to the \u2018noise floor\u2019 in gene expression in a population of dividing cells

Experiments with cells reveal the existence of a lower bound for protein noise, the noise floor, in highly expressed genes. Its origins are still debated. We propose a minimal model of gene expression in a proliferating bacterial cell population. The model predicts the existence of a noise floor and it semi-quantitatively reproduces the curved shape of the experimental noise vs. mean protein concentration plots. When the cell volume increases in a different manner than does the mean protein copy number, the noise floor level is determined by the cell population’s age structure and by the dependence of the mean protein concentration on cell age. Additionally, the noise floor level may depend on a biological limit for the mean number of bursts in the cell cycle. In that case, the noise floor level depends on the burst size distribution width but it is insensitive to the mean burst size. Our model quantifies the contributions of each of these mechanisms to gene expression noise.

A Stochastic Model of Gene Expression with Polymerase Recruitment and Pause Release

Transcriptional bursting is a major source of noise in gene expression. The telegraph model of gene expression, whereby transcription switches between on and off states, is the dominant model for bursting. Recently, it was shown that the telegraph model cannot explain a number of experimental observations from perturbation data. Here, we study an alternative model that is consistent with the data and which explicitly describes RNA polymerase recruitment and polymerase pause release, two steps necessary for messenger RNA (mRNA) production. We derive the exact steady-state distribution of mRNA numbers and an approximate steady-state distribution of protein numbers, which are given by generalized hypergeometric functions. The theory is used to calculate the relative sensitivity of the coefficient of variation of mRNA fluctuations for thousands of genes in mouse fibroblasts. This indicates that the size of fluctuations is mostly sensitive to the rate of burst initiation and the mRNA degradation rate. Furthermore, we show that 1) the time-dependent distribution of mRNA numbers is accurately approximated by a modified telegraph model with a Michaelis-Menten like dependence of the effective transcription rate on RNA polymerase abundance, and 2) the model predicts that if the polymerase recruitment rate is comparable or less than the pause release rate, then upon gene replication, the mean number of RNA per cell remains approximately constant. This gene dosage compensation property has been experimentally observed and cannot be explained by the telegraph model with constant rates.

Mixture distributions in a stochastic gene expression model with delayed feedback: a WKB approximation approach

Noise in gene expression can be substantively affected by the presence of production delay. Here we consider a mathematical model with bursty production of protein, a one-step production delay (the passage of which activates the protein), and feedback in the frequency of bursts. We specifically focus on examining the steady-state behaviour of the model in the slow-activation (i.e. large-delay) regime. Using a formal asymptotic approach, we derive an autonomous ordinary differential equation for the inactive protein that applies in the slow-activation regime. If the differential equation is monostable, the steady-state distribution of the inactive (active) protein is approximated by a single Gaussian (Poisson) mode located at the globally stable fixed point of the differential equation. If the differential equation is bistable (due to cooperative positive feedback), the steady-state distribution of the inactive (active) protein is approximated by a mixture of Gaussian (Poisson) modes located at the stable fixed points; the weights of the modes are determined from a WKB approximation to the stationary distribution. The asymptotic results are compared to numerical solutions of the chemical master equation.

Orthogonal tuning of gene expression noise using CRISPR-Cas.

The control of gene expression noise is important for improving drug treatment and the performance of synthetic biological systems. Previous work has tuned gene expression noise by changing the rate of transcription initiation, mRNA degradation, and mRNA translation. However, these methods are invasive: they require changes to the target genetic components. Here, we create an orthogonal system based on CRISPR-dCas9 to tune gene expression noise. Specifically, we modulate the gene expression noise of a reporter gene in Escherichia coli by incorporating CRISPR activation and repression (CRISPRar) simultaneously in a single cell. The CRISPRar uses a single dCas9 that recognizes two different single guide RNAs (sgRNA). We build a library of sgRNA variants with different expression activation and repression strengths. We find that expression noise and mean of a reporter gene can be tuned independently by CRISPRar. Our results suggest that the expression noise is tuned by the competition between two sgRNAs that modulate the binding of RNA polymerase to promoters. The CRISPRar may change how we tune expression noise at the genomic level. Our work has broad impacts on the study of gene functions, phenotypical heterogeneity, and genetic circuit control.

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

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

Enhancement of gene expression noise from transcription\xa0factor binding to genomic decoy sites

The genome contains several high-affinity non-functional binding sites for transcription factors (TFs) creating a hidden and unexplored layer of gene regulation. We investigate the role of such “decoy sites” in controlling noise (random fluctuations) in the level of a TF that is synthesized in stochastic bursts. Prior studies have assumed that decoy-bound TFs are protected from degradation, and in this case decoys function to buffer noise. Relaxing this assumption to consider arbitrary degradation rates for both bound/unbound TF states, we find rich noise behaviors. For low-affinity decoys, noise in the level of unbound TF always monotonically decreases to the Poisson limit with increasing decoy numbers. In contrast, for high-affinity decoys, noise levels first increase with increasing decoy numbers, before decreasing back to the Poisson limit. Interestingly, while protection of bound TFs from degradation slows the time-scale of fluctuations in the unbound TF levels, the decay of bound TFs leads to faster fluctuations and smaller noise propagation to downstream target proteins. In summary, our analysis reveals stochastic dynamics emerging from nonspecific binding of TFs and highlights the dual role of decoys as attenuators or amplifiers of gene expression noise depending on their binding affinity and stability of the bound TF.