Stochastic Expression Research Articles

Single-cell Transcriptomes Reveal Characteristics of MicroRNAs in Gene Expression Noise Reduction

Isogenic cells growing in identical environments show cell-to-cell variations because of the stochasticity in gene expression. High levels of variation or noise can disrupt robust gene expression and result in tremendous consequences for cell behaviors. In this work, we showed evidence from single-cell RNA sequencing data analysis that microRNAs (miRNAs) can reduce gene expression noise at the mRNA level in mouse cells. We identified that the miRNA expression level, number of targets, target pool abundance, and miRNA–target interaction strength are the key features contributing to noise repression. miRNAs tend to work together in cooperative subnetworks to repress target noise synergistically in a cell type-specific manner. By building a physical model of post-transcriptional regulation and observing in synthetic gene circuits, we demonstrated that accelerated degradation with elevated transcriptional activation of the miRNA target provides resistance to extrinsic fluctuations. Together, through the integrated analysis of single-cell RNA and miRNA expression profiles, we demonstrated that miRNAs are important post-transcriptional regulators for reducing gene expression noise and conferring robustness to biological processes.

Pareto optimality between growth-rate and lag-time couples metabolic noise to phenotypic heterogeneity in Escherichia coli

Despite mounting evidence that in clonal bacterial populations, phenotypic variability originates from stochasticity in gene expression, little is known about noise-shaping evolutionary forces and how expression noise translates to phenotypic differences. Here we developed a high-throughput assay that uses a redox-sensitive dye to couple growth of thousands of bacterial colonies to their respiratory activity and show that in Escherichia coli, noisy regulation of lower glycolysis and citric acid cycle is responsible for large variations in respiratory metabolism. We found that these variations are Pareto optimal to maximization of growth rate and minimization of lag time, two objectives competing between fermentative and respiratory metabolism. Metabolome-based analysis revealed the role of respiratory metabolism in preventing the accumulation of toxic intermediates of branched chain amino acid biosynthesis, thereby supporting early onset of cell growth after carbon starvation. We propose that optimal metabolic tradeoffs play a key role in shaping and preserving phenotypic heterogeneity and adaptation to fluctuating environments.

Coordination of gene expression noise with cell size: analytical results for agent-based models of growing cell populations.

The chemical master equation and the Gillespie algorithm are widely used to model the reaction kinetics inside living cells. It is thereby assumed that cell growth and division can be modelled through effective dilution reactions and extrinsic noise sources. We here re-examine these paradigms through developing an analytical agent-based framework of growing and dividing cells accompanied by an exact simulation algorithm, which allows us to quantify the dynamics of virtually any intracellular reaction network affected by stochastic cell size control and division noise. We find that the solution of the chemical master equation-including static extrinsic noise-exactly agrees with the agent-based formulation when the network under study exhibits stochastic concentration homeostasis, a novel condition that generalizes concentration homeostasis in deterministic systems to higher order moments and distributions. We illustrate stochastic concentration homeostasis for a range of common gene expression networks. When this condition is not met, we demonstrate by extending the linear noise approximation to agent-based models that the dependence of gene expression noise on cell size can qualitatively deviate from the chemical master equation. Surprisingly, the total noise of the agent-based approach can still be well approximated by extrinsic noise models.

Dynamical analysis of the fission yeast cell cycle via Markov chain.

The cell cycle is a complex network involved in the regulation of cell growth and proliferation. Intrinsic molecular noise in gene expression in the cell cycle network can generate fluctuations in protein concentration. How the cell cycle network maintains its robust transitions between cell cycle phases in the presence of these fluctuations remains unclear. To understand the complex and robust behavior of the cell cycle system in the presence of intrinsic noise, we developed a Markov model for the fission yeast cell cycle system. We quantified the effect of noise on gene and protein activity and on the probability of transition between different phases of the cell cycle. Our analysis shows how network perturbations decide the fate of the cell. Our model predicts that the cell cycle pathway (subsequent transitions from [Formula: see text]) is the most robust and probable pathway among all possible trajectories in the cell cycle network. We performed a sensitivity analysis to find correlations between protein interaction weights and transition probabilities between cell cycle phases. The sensitivity analysis predicts how network perturbations affect the transition probability between different cell cycle phases and, consequently, affect different cell fates, thus, forming testable in vitro/in vivo hypotheses. Our simulation results agree with published experimental findings and reveal how noise in the cell cycle regulatory network can affect cell cycle progression.

Facial emotion detection in Vestibular Schwannoma patients with and without facial paresis

ABSTRACT This study investigates whether there exist differences in facial emotion detection accuracy in patients suffering from Vestibular Schwannoma (VS) due to their facial paresis. Forty-four VS patients, half of them with, and half of them without a facial paresis, had to classify pictures of facial expressions as being emotional or non-emotional. The visual information of images was systematically manipulated by adding different levels of visual noise. The study had a mixed design with emotional expression (happy vs. angry) and visual noise level (10% to 80%) as repeated measures and facial paresis (present vs. absent) and degree of facial dysfunction as between subjects’ factors. Emotion detection accuracy declined when visual information declined, an effect that was stronger for anger than for happy expressions. Overall, emotion detection accuracy for happy and angry faces did not differ between VS patients with or without a facial paresis, although exploratory analyses suggest that the ability to recognize emotions in angry facial expressions was slightly more impaired in patients with facial paresis. The findings are discussed in the context of the effects of facial paresis on emotion detection, and the role of facial mimicry, in particular, as an important mechanism for facial emotion processing and understanding.

Modeling expression ranks for noise-tolerant differential expression analysis of scRNA-seq data

Systematic delineation of complex biological systems is an ever-challenging and resource-intensive process. Single-cell transcriptomics allows us to study cell-to-cell variability in complex tissues at an unprecedented resolution. Accurate modeling of gene expression plays a critical role in the statistical determination of tissue-specific gene expression patterns. In the past few years, considerable efforts have been made to identify appropriate parametric models for single-cell expression data. The zero-inflated version of Poisson/negative binomial and log-normal distributions have emerged as the most popular alternatives owing to their ability to accommodate high dropout rates, as commonly observed in single-cell data. Although the majority of the parametric approaches directly model expression estimates, we explore the potential of modeling expression ranks, as robust surrogates for transcript abundance. Here we examined the performance of the discrete generalized beta distribution (DGBD) on real data and devised a Wald-type test for comparing gene expression across two phenotypically divergent groups of single cells. We performed a comprehensive assessment of the proposed method to understand its advantages compared with some of the existing best-practice approaches. We concluded that besides striking a reasonable balance between Type I and Type II errors, ROSeq, the proposed differential expression test, is exceptionally robust to expression noise and scales rapidly with increasing sample size. For wider dissemination and adoption of the method, we created an R package called ROSeq and made it available on the Bioconductor platform.

Massively parallel gene expression variation measurement of a synonymous codon library

BackgroundCell-to-cell variation in gene expression strongly affects population behavior and is key to multiple biological processes. While codon usage is known to affect ensemble gene expression, how codon usage influences variation in gene expression between single cells is not well understood.ResultsHere, we used a Sort-seq based massively parallel strategy to quantify gene expression variation from a green fluorescent protein (GFP) library containing synonymous codons in Escherichia coli. We found that sequences containing codons with higher tRNA Adaptation Index (TAI) scores, and higher codon adaptation index (CAI) scores, have higher GFP variance. This trend is not observed for codons with high Normalized Translation Efficiency Index (nTE) scores nor from the free energy of folding of the mRNA secondary structure. GFP noise, or squared coefficient of variance (CV2), scales with mean protein abundance for low-abundant proteins but does not change at high mean protein abundance.ConclusionsOur results suggest that the main source of noise for high-abundance proteins is likely not originating at translation elongation. Additionally, the drastic change in mean protein abundance with small changes in protein noise seen from our library implies that codon optimization can be performed without concerning gene expression noise for biotechnology applications.

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

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

Generation of Genetic Tools for Gauging Multiple-Gene Expression at the Single-Cell Level.

Key microbial processes in many bacterial species are heterogeneously expressed in single cells of bacterial populations. However, the paucity of adequate molecular tools for live, real-time monitoring of multiple-gene expression at the single-cell level has limited the understanding of phenotypic heterogeneity. To investigate phenotypic heterogeneity in the ubiquitous opportunistic pathogen Pseudomonas aeruginosa, a genetic tool that allows gauging multiple-gene expression at the single-cell level has been generated. This tool, named pRGC, consists of a promoter-probe vector for transcriptional fusions that carries three reporter genes coding for the fluorescent proteins mCherry, green fluorescent protein (GFP), and cyan fluorescent protein (CFP). The pRGC vector has been characterized and validated via single-cell gene expression analysis of both constitutive and iron-regulated promoters, showing clear discrimination of the three fluorescence signals in single cells of a P. aeruginosa population without the need for image processing for spectral cross talk correction. In addition, two pRGC variants have been generated for either (i) integration of the reporter gene cassette into a single neutral site of P. aeruginosa chromosome that is suitable for long-term experiments in the absence of antibiotic selection or (ii) replication in bacterial genera other than Pseudomonas The easy-to-use genetic tools generated in this study will allow rapid and cost-effective investigation of multiple-gene expression in populations of environmental and pathogenic bacteria, hopefully advancing the understanding of microbial phenotypic heterogeneity.IMPORTANCE Within a bacterial population, single cells can differently express some genes, even though they are genetically identical and experience the same chemical and physical stimuli. This phenomenon, known as phenotypic heterogeneity, is mainly driven by gene expression noise and results in the emergence of bacterial subpopulations with distinct phenotypes. The analysis of gene expression at the single-cell level has shown that phenotypic heterogeneity is associated with key bacterial processes, including competence, sporulation, and persistence. In this study, new genetic tools have been generated that allow easy cloning of up to three promoters upstream of distinct fluorescent genes, making it possible to gauge multiple-gene expression at the single-cell level by fluorescence microscopy without the need for advanced image-processing procedures. A proof of concept has been provided by investigating iron uptake and iron storage gene expression in response to iron availability in P. aeruginosa.

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.

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.

Analysis of microRNA Regulation in Single Cells.

MicroRNAs (miRNAs) regulate gene expression by binding to mRNAs. Consequently, they reduce target gene expression levels and expression variability, also known as "noise." Single-cell RNA sequencing (scRNA-seq) technology has been used to study miRNA and mRNA expression in single cells, and has demonstrated its strength in quantifying cell-to-cell variation. Here we describe how to investigate miRNA regulation using data with both mRNA and miRNA expression in single cell format. We show that miRNAs reduce the expression levels and also expression noise of target genes in single cells. Finally, we also discuss potential improvements in experimental design and computational analysis of scRNA-seq in order to reduce or partition the technical noise.

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.

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.

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.