The Low Noise Limit in Gene Expression

Optimisation of gene expression noise for cellular persistence against lethal events

Bacterial cell persistence, crucial for survival under adverse conditions like antibiotic exposure, is intrinsically linked to stochastic fluctuations in gene expression. Certain genes, while inhibiting growth under normal circumstances, confer tolerance to antibiotics at elevated expression levels. The occurrence of antibiotic events lead to instantaneous cellular responses with varied survival probabilities correlated with gene expression levels. Notably, cells with lower protein concentrations face higher mortality rates. This study aims to elucidate an optimal strategy for protein expression conducive to cellular survival. Through comprehensive mathematical analysis, we determine the optimal burst size and frequency that maximise cell proliferation. Furthermore, we explore how the optimal expression distribution changes as the cost of protein expression to growth escalates. Our model reveals a hysteresis phenomenon, characterised by discontinuous transitions between deterministic and stochastic optima. Intriguingly, stochastic optima possess a noise floor, representing the minimal level of fluctuations essential for optimal cellular resilience.

N-terminal tagging of RNA Polymerase II shapes transcriptomes more than C-terminal alterations

RNA polymerase II (Pol II) has a C-terminal domain (CTD) that is unstructured, consisting of a large number of heptad repeats, and whose precise function remains unclear. Here, we investigate how altering the CTD's length and fusing it with protein tags affects transcriptional output on a genome-wide scale in mammalian cells at single-cell resolution. While transcription generally appears to occur in burst-like fashion, where RNA is predominantly made during short bursts of activity that are interspersed with periods of transcriptional silence, the CTD's role in shaping these dynamics seems gene-dependent; global patterns of bursting appear mostly robust to CTD alterations. Introducing protein tags with defined structures to the N terminus cause transcriptome-wide effects, however. We find the type of tag to dominate characteristics of the resulting transcriptomes. This is possibly due to Pol II-interacting factors, including non-coding RNAs, whose expression correlates with the tags. Proteins involved in liquid-liquid phase separation appear prominently.

A Kinetic Finite Volume Discretization of the Multidimensional PIDE Model for Gene Regulatory Networks

In this paper, a finite volume discretization scheme for partial integro-differential equations (PIDEs) describing the temporal evolution of protein distribution in gene regulatory networks is proposed. It is shown that the obtained set of ODEs can be formally represented as a compartmental kinetic system with a strongly connected reaction graph. This allows the application of the theory of nonnegative and compartmental systems for the qualitative analysis of the approximating dynamics. In this framework, it is straightforward to show the existence, uniqueness and stability of equilibria. Moreover, the computation of the stationary probability distribution can be traced back to the solution of linear equations. The discretization scheme is presented for one and multiple dimensional models separately. Illustrative computational examples show the precision of the approach, and good agreement with previous results in the literature.

Interrogating the Function of Bicistronic Translational Control Elements to Improve Consistency of Gene Expression

Context independent gene expression is required for genetic circuits to maintain consistent and predicable behavior. Previous efforts to develop context independent translation have leveraged the helicase activity of translating ribosomes via bicistronic design translational control elements (BCDs) located within an efficiently translated leader peptide. We have developed a series of bicistronic translational control elements with strengths that span several orders of magnitude, maintain consistent expression levels across diverse sequence contexts, and are agnostic to common ligation sequences used in modular cloning systems. We have used this series of BCDs to investigate several features of this design, including the spacing of the start and stop codons, the nucleotide identity upstream of the start codon, and factors affecting translation of the leader peptide. To demonstrate the flexibility of this architecture and their value as a generic modular expression control cassette for synthetic biology, we have developed a set of robust BCDs for use in several Rhodococcus species.

Epigenetic fluctuations underlie gene expression timescales and variability

Isogenic populations of mammalian cells exhibit significant gene expression variability. This variability can be separated into two sources, cis, or allele-specific sources, and trans and global processes. Furthermore, each source of variability has its own timescale. Fast timescales will result in rapid fluctuation of gene expression whereas slow timescales will result in longer persistence of gene expression levels over time. Here we investigated sources of gene expression that are intrinsic, i.e. coming from cis-regulatory factors and follow slow timescales. To do so, we developed a reporter system that isolates allele-specific variability and measures its persistence in imaging and long-term fluctuation analysis experiments. Our results identify a new source of gene expression variability that is allele-specific but that is fluctuating on timescales of days. We hypothesized that allele-specific fluctuations of epigenetic regulatory factors are responsible for the newly discovered allele-specific and slow source of gene expression variability. Using mathematical modeling we showed that the addition of this effect to the two-state model is sufficient to account for all empirical observation. Furthermore, using direct assays of chromatin markers we find fluctuation in H3K4me3 levels that match the observed changes in gene expression levels providing direct experimental support of our model. Collectively, our work shows that slow fluctuations of regulatory chromatin modifications contribute to the variability in gene expression.

Silent transcription intervals and translational bursting lead to diverse phenotypic switching

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

Disentangling Intrinsic and Extrinsic Gene Expression Noise in Growing Cells

Gene expression is a stochastic process. Despite the increase of protein numbers in growing cells, the protein concentrations are often found to be confined within small ranges throughout the cell cycle. Generally, the noise in protein concentration can be decomposed into an intrinsic and an extrinsic component, where the former vanishes for high expression levels. Considering the time trajectory of protein concentration as a random walker in the concentration space, an effective restoring force (with a corresponding "spring constant") must exist to prevent the divergence of concentration due to random fluctuations. In this work, we prove that the magnitude of the effective spring constant is directly related to the fraction of intrinsic noise in the total protein concentration noise. We show that one can infer the magnitude of intrinsic, extrinsic, and measurement noises of gene expression solely based on time-resolved data of protein concentration, without any a priori knowledge of the underlying gene expression dynamics. We apply this method to experimental data of single-cell bacterial gene expression. The results allow us to estimate the average copy numbers and the translation burst parameters of the studied proteins.

3 '-5 ' crosstalk contributes to transcriptional bursting

BACKGROUND

Transcription in mammalian cells is a complex stochastic process involving shuttling of polymerase between genes and phase-separated liquid condensates. It occurs in bursts, which results in vastly different numbers of an mRNA species in isogenic cell populations. Several factors contributing to transcriptional bursting have been identified, usually classified as intrinsic, in other words local to single genes, or extrinsic, relating to the macroscopic state of the cell. However, some possible contributors have not been explored yet. Here, we focus on processes at the 3 ' and 5 ' ends of a gene that enable reinitiation of transcription upon termination.

RESULTS

Using Bayesian methodology, we measure the transcriptional bursting in inducible transgenes, showing that perturbation of polymerase shuttling typically reduces burst size, increases burst frequency, and thus limits transcriptional noise. Analysis based on paired-end tag sequencing (PolII ChIA-PET) suggests that this effect is genome wide. The observed noise patterns are also reproduced by a generative model that captures major characteristics of the polymerase flux between the ends of a gene and a phase-separated compartment.

CONCLUSIONS

Interactions between the 3 ' and 5 ' ends of a gene, which facilitate polymerase recycling, are major contributors to transcriptional noise.

Contributions to the 'noise floor' 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.

Ab initio scaling laws between noise and mean of gene expression

Gene expression is inherently noisy due to fluctuations occurring at the molecular level. From a top-down perspective, noise has been traditionally decomposed into an intrinsic component that scales inversely with the mean expression level and an extrinsic component that is constant in absence of regulatory changes. Here, we adopt a bottom-up approach to reveal that extrinsic noise, by itself, can follow the aforementioned decomposition, which entails that one component of it can be confounded with intrinsic noise. Analytical expressions of the noise-mean relationship were derived for different scenarios, which were in part supported by numerical simulations.

Mammalian gene expression variability is explained by underlying cell state

Gene expression variability in mammalian systems plays an important role in physiological and pathophysiological conditions. This variability can come from differential regulation related to cell state (extrinsic) and allele-specific transcriptional bursting (intrinsic). Yet, the relative contribution of these two distinct sources is unknown. Here, we exploit the qualitative difference in the patterns of covariance between these two sources to quantify their relative contributions to expression variance in mammalian cells. Using multiplexed error robust RNA fluorescent in situ hybridization (MERFISH), we measured the multivariate gene expression distribution of 150 genes related to Ca2+ signaling coupled with the dynamic Ca2+ response of live cells to ATP. We show that after controlling for cellular phenotypic states such as size, cell cycle stage, and Ca2+ response to ATP, the remaining variability is effectively at the Poisson limit for most genes. These findings demonstrate that the majority of expression variability results from cell state differences and that the contribution of transcriptional bursting is relatively minimal.

Cell Size-Based Decision-Making of a Viral Gene Circuit

Latently infected T cells able to reinitiate viral propagation throughout the body remain a major barrier to curing HIV. Distinguishing between latently infected cells and uninfected cells will advance efforts for viral eradication. HIV decision-making between latency and active replication is stochastic, and drug cocktails that increase bursts of viral gene expression enhance reactivation from latency. Here, we show that a larger host-cell size provides a natural cellular mechanism for enhancing burst size of viral expression and is necessary to destabilize the latent state and bias viral decision-making. Latently infected Jurkat and primary CD4+ T cells reactivate exclusively in larger activated cells, while smaller cells remain silent. In addition, reactivation is cell-cycle dependent and can be modulated with cell-cycle-arresting compounds. Cell size and cell-cycle dependent decision-making of viral circuits may guide stochastic design strategies and applications in synthetic biology and may provide important determinants to advance diagnostics and therapies.

Bohn-Wippert et al. investigate reactivation of T cells latently infected with HIV. They discover that only larger cells exit latency, while smaller cells remain silent. Viral expression bursts are cell size and cell-cycle dependent, presenting dynamic cell states, capable of active control, as sources of viral fate determination.

Homeostasis of protein and mRNA concentrations in growing cells

Many experiments show that the numbers of mRNA and protein are proportional to the cell volume in growing cells. However, models of stochastic gene expression often assume constant transcription rate per gene and constant translation rate per mRNA, which are incompatible with these experiments. Here, we construct a minimal gene expression model to fill this gap. Assuming ribosomes and RNA polymerases are limiting in gene expression, we show that the numbers of proteins and mRNAs both grow exponentially during the cell cycle and that the concentrations of all mRNAs and proteins achieve cellular homeostasis; the competition between genes for the RNA polymerases makes the transcription rate independent of the genome number. Furthermore, by extending the model to situations in which DNA (mRNA) can be saturated by RNA polymerases (ribosomes) and becomes limiting, we predict a transition from exponential to linear growth of cell volume as the protein-to-DNA ratio increases.

Probing Intercell Variability Using Bulk Measurements

The measurement of noise is critical when assessing the design and function of synthetic biological systems. Cell-to-cell variability can be quantified experimentally using single-cell measurement techniques such as flow cytometry and fluorescent microscopy. However, these approaches are costly and impractical for high-throughput parallelized experiments, which are frequently conducted using plate-reader devices. In this paper we describe reporter systems that allow estimation of the cell-to-cell variability in a biological system's output using only measurements of a cell culture's bulk properties. We analyze one potential implementation of such a system that is based upon a fluorescent protein FRET reporter pair, finding that with typical parameters from the literature it is able to reliably estimate variability. We also briefly describe an alternate implementation based upon an activating sRNA circuit. The feasible region of parameter values for which the reporter system can function is assessed, and the dependence of its performance on both extrinsic and intrinsic noise is investigated. Experimental realization of these constructs can yield novel reporter systems that allow measurement of a synthetic gene circuit's output, as well as the intrapopulation variability of this output, at little added cost.

The Pathway To Establishing HIV Latency Is Critical to How Latency Is Maintained and Reversed

HIV infection requires lifelong antiretroviral therapy because of the persistence of latently infected CD4⁺ T cells. The induction of virus expression from latently infected cells occurs following T cell receptor (TCR) activation, but not all latently infected cells respond to TCR stimulation. We compared two models of latently infected cells using an enhanced green fluorescent protein (EGFP) reporter virus to infect CCL19-treated resting CD4⁺ (rCD4⁺) T cells (preactivation latency) or activated CD4⁺ T cells that returned to a resting state (postactivation latency). We isolated latently infected cells by sorting for EGFP-negative (EGFP^-) cells after infection. These cells were cultured with antivirals and stimulated with anti-CD3/anti-CD28, mitogens, and latency-reversing agents (LRAs) and cocultured with monocytes and anti-CD3. Spontaneous EGFP expression was more frequent in postactivation than in preactivation latency. Stimulation of latently infected cells with monocytes/anti-CD3 resulted in an increase in EGFP expression compared to that for unstimulated controls using the preactivation latency model but led to a reduction in EGFP expression in the postactivation latency model. The reduced EGFP expression was not associated with reductions in the levels of viral DNA or T cell proliferation but depended on direct contact between monocytes and T cells. Monocytes added to the postactivation latency model during the establishment of latency reduced spontaneous virus expression, suggesting that monocyte-T cell interactions at an early time point postinfection can maintain HIV latency. This direct comparison of pre- and postactivation latency suggests that effective strategies needed to reverse latency will depend on how latency is established.IMPORTANCE One strategy being evaluated to eliminate latently infected cells that persist in HIV-infected individuals on antiretroviral therapy (ART) is to activate HIV expression or production with the goal of inducing virus-mediated cytolysis or immune-mediated clearance of infected cells. The gold standard for the activation of latent virus is T cell receptor stimulation with anti-CD3/anti-CD28. However, this stimulus activates only a small proportion of latently infected cells. We show clear differences in the responses of latently infected cells to activating stimuli based on how latent infection is established, an observation that may potentially explain the persistence of noninduced intact proviruses in HIV-infected individuals on ART.

Imaging mRNA In Vivo, from Birth to Death

RNA is the fundamental information transfer system in the cell. The ability to follow single messenger RNAs (mRNAs) from transcription to degradation with fluorescent probes gives quantitative information about how the information is transferred from DNA to proteins. This review focuses on the latest technological developments in the field of single-mRNA detection and their usage to study gene expression in both fixed and live cells. By describing the application of these imaging tools, we follow the journey of mRNA from transcription to decay in single cells, with single-molecule resolution. We review current theoretical models for describing transcription and translation that were generated by single-molecule and single-cell studies. These methods provide a basis to study how single-molecule interactions generate phenotypes, fundamentally changing our understating of gene expression regulation.

Macromolecular Crowding Induces Spatial Correlations That Control Gene Expression Bursting Patterns

Recent superresolution microscopy studies in E. coli demonstrate that the cytoplasm has highly variable local concentrations where macromolecular crowding plays a central role in establishing membrane-less compartmentalization. This spatial inhomogeneity significantly influences molecular transport and association processes central to gene expression. Yet, little is known about how macromolecular crowding influences gene expression bursting-the episodic process where mRNA and proteins are produced in bursts. Here, we simultaneously measured mRNA and protein reporters in cell-free systems, showing that macromolecular crowding decoupled the well-known relationship between fluctuations in the protein population (noise) and mRNA population statistics. Crowded environments led to a 10-fold increase in protein noise even though there were only modest changes in the mRNA population and fluctuations. Instead, cell-like macromolecular crowding created an inhomogeneous spatial distribution of mRNA ("spatial noise") that led to large variability in the protein production burst size. As a result, the mRNA spatial noise created large temporal fluctuations in the protein population. These results highlight the interplay between macromolecular crowding, spatial inhomogeneities, and the resulting dynamics of gene expression, and provide insights into using these organizational principles in both cell-based and cell-free synthetic biology.

A Post-Transcriptional Feedback Mechanism for Noise Suppression and Fate Stabilization

Diverse biological systems utilize fluctuations ("noise") in gene expression to drive lineage-commitment decisions. However, once a commitment is made, noise becomes detrimental to reliable function, and the mechanisms enabling post-commitment noise suppression are unclear. Here, we find that architectural constraints on noise suppression are overcome to stabilize fate commitment. Using single-molecule and time-lapse imaging, we find that-after a noise-driven event-human immunodeficiency virus (HIV) strongly attenuates expression noise through a non-transcriptional negative-feedback circuit. Feedback is established through a serial cascade of post-transcriptional splicing, whereby proteins generated from spliced mRNAs auto-deplete their own precursor unspliced mRNAs. Strikingly, this auto-depletion circuitry minimizes noise to stabilize HIV's commitment decision, and a noise-suppression molecule promotes stabilization. This feedback mechanism for noise suppression suggests a functional role for delayed splicing in other systems and may represent a generalizable architecture of diverse homeostatic signaling circuits.

Fine-tuning of noise in gene expression with nucleosome remodeling

Engineering stochastic fluctuations of gene expression (or “noise”) is integral to precisely bias cellular-fate decisions and statistical phenotypes in both single-cell and multi-cellular systems. Epigenetic regulation has been shown to constitute a large source of noise, and thus, engineering stochasticity is deeply intertwined with epigenetics. Here, utilizing chromatin remodeling, we report that Caffeic acid phenethyl ester (CA) and Pyrimethamine (PYR), two inhibitors of BAF250a, a subunit of the Brahma-associated factor (BAF) nucleosome remodeling complex, enable differential and tunable control of noise in transcription and translation from the human immunodeficiency virus long terminal repeat promoter in a dose and time-dependent manner. CA conserves noise levels while increasing mean abundance, resulting in direct tuning of the transcriptional burst size, while PYR strictly increases transcriptional initiation frequency while conserving a constant transcriptional burst size. Time-dependent treatment with CA reveals non-continuous tuning with noise oscillating at a constant mean abundance at early time points and the burst size increasing for treatments after 5 h. Treatments combining CA and Protein Kinase C agonists result in an even larger increase of abundance while conserving noise levels with a highly non-linear increase in variance of up to 63× untreated controls. Finally, drug combinations provide non-antagonistic combinatorial tuning of gene expression noise and map a noise phase space for future applications with viral and synthetic gene vectors. Active remodeling of nucleosomes and BAF-mediated control of gene expression noise expand a toolbox for the future design and engineering of stochasticity in living systems.

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.

Translation initiation events on structured eukaryotic mRNAs generate gene expression noise

Gene expression stochasticity plays a major role in biology, creating non-genetic cellular individuality and influencing multiple processes, including differentiation and stress responses. We have addressed the lack of knowledge about posttranscriptional contributions to noise by determining cell-to-cell variations in the abundance of mRNA and reporter protein in yeast. Two types of structural element, a stem–loop and a poly(G) motif, not only inhibit translation initiation when inserted into an mRNA 5΄ untranslated region, but also generate noise. The noise-enhancing effect of the stem–loop structure also remains operational when combined with an upstream open reading frame. This has broad significance, since these elements are known to modulate the expression of a diversity of eukaryotic genes. Our findings suggest a mechanism for posttranscriptional noise generation that will contribute to understanding of the generally poor correlation between protein-level stochasticity and transcriptional bursting. We propose that posttranscriptional stochasticity can be linked to cycles of folding/unfolding of a stem–loop structure, or to interconversion between higher-order structural conformations of a G-rich motif, and have created a correspondingly configured computational model that generates fits to the experimental data. Stochastic events occurring during the ribosomal scanning process can therefore feature alongside transcriptional bursting as a source of noise.

Effects of growth rate and promoter activity on single-cell protein expression

Protein expression in a single cell depends on its global physiological state. Moreover, genetically-identical cells exhibit variability (noise) in protein expression, arising from the stochastic nature of biochemical processes, cell growth and division. While it is well understood how cellular growth rate influences mean protein expression, little is known about the relationship between growth rate and noise in protein expression. Here we quantify this relationship in Bacillus subtilis by a novel combination of experiments and theory. We measure the effects of promoter activity and growth rate on the expression of a fluorescent protein in single cells. We disentangle the observed protein expression noise into protein-specific and systemic contributions, using theory and variance decomposition. We find that noise in protein expression depends solely on mean expression levels, regardless of whether expression is set by promoter activity or growth rate, and that noise increases linearly with growth rate. Our results can aid studies of (synthetic) gene circuits of single cells and their condition dependence.

Resource Sharing Controls Gene Expression Bursting

Episodic gene expression, with periods of high expression separated by periods of no expression, is a pervasive biological phenomenon. This bursty pattern of expression draws from a finite reservoir of expression machinery in a highly time variant way, i.e., requiring no resources most of the time but drawing heavily on them during short intense bursts, that intimately links expression bursting and resource sharing. Yet, most recent investigations have focused on specific molecular mechanisms intrinsic to the bursty behavior of individual genes, while little is known about the interplay between resource sharing and global expression bursting behavior. Here, we confine Escherichia coli cell extract in both cell-sized microfluidic chambers and lipid-based vesicles to explore how resource sharing influences expression bursting. Interestingly, expression burst size, but not burst frequency, is highly sensitive to the size of the shared transcription and translation resource pools. The intriguing implication of these results is that expression bursts are more readily amplified than initiated, suggesting that burst formation occurs through positive feedback or cooperativity. When extrapolated to prokaryotic cells, these results suggest that large translational bursts may be correlated with large transcriptional bursts. This correlation is supported by recently reported transcription and translation bursting studies in E. coli. The results reported here demonstrate a strong intimate link between global expression burst patterns and resource sharing, and they suggest that bursting plays an important role in optimizing the use of limited, shared expression resources.