Gene Expression Noise Research Articles

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.

Translation regulation by RNA stem-loops can reduce gene expression noise

BackgroundStochastic modelling plays a crucial role in comprehending the dynamics of intracellular events in various biochemical systems, including gene-expression models. Cell-to-cell variability arises from the stochasticity or noise in the levels of gene products such as messenger RNA (mRNA) and protein. The sources of noise can stem from different factors, including structural elements. Recent studies have revealed that the mRNA structure can be more intricate than previously assumed.ResultsHere, we focus on the formation of stem-loops and present a reinterpretation of previous data, offering new insights. Our analysis demonstrates that stem-loops that restrict translation have the potential to reduce noise.ConclusionsIn conclusion, we investigate a structured/generalised version of a stochastic gene-expression model, wherein mRNA molecules can be found in one of their finite number of different states and transition between them. By characterising and deriving non-trivial analytical expressions for the steady-state protein distribution, we provide two specific examples which can be readily obtained from the structured/generalised model, showcasing the model’s practical applicability.

Multimodal screen identifies noise-regulatory proteins.

Gene-expression noise can influence cell-fate choices across pathology and physiology. However, a crucial question persists: do regulatory proteins or pathways exist that control noise independently of mean expression levels? Our integrative approach, combining single-cell RNA sequencing with proteomics and regulator enrichment analysis, identifies 32 putative noise regulators. SON, a nuclear speckle-associated protein, alters transcriptional noise without changing mean expression levels. Furthermore, SON's noise control can propagate to the protein level. Long-read and total RNA sequencing shows that SON's noise control does not significantly change isoform usage or splicing efficiency. Moreover, SON depletion reduces state switching in pluripotent mouse embryonic stem cells and impacts their fate choice during differentiation. Collectively, we demonstrate a class of proteins that control noise orthogonally to mean expression levels. This work serves as a proof of concept that can identify other functional noise regulators throughout development and disease progression.

AFSC: A self-supervised augmentation-free spatial clustering method based on contrastive learning for identifying spatial domains

Recent research in spatial transcriptomics allows researchers to analyze gene expression without losing spatial information. Spatial information can assist in cell communication, identification of new cell subtypes, which provides important research methods for multiple fields such as microenvironment interactions and pathological processes of diseases. Identifying spatial domains is an important step in spatial transcriptomics analysis, and improving spatial clustering methods can benefit for identifying spatial domains. In addition to eliminating noise in original gene expression, how to use spatial information to assist clustering has also become a new problem. A variety of calculating methods have been applied to spatial clustering, including contrastive learning methods. However, existing spatial clustering methods based on contrastive learning use data augmentation to generate positive and negative pairs, which will inevitably destroy the biological meaning of the data. We propose a new self-supervised spatial clustering method based on contrastive learning, Augmentation-Free Spatial Clustering (AFSC), which integrates spatial information and gene expression to learn latent representations. We construct a contrastive learning module without negative pairs or data augmentation by designing Teacher and Student Encoder. We also design an unsupervised clustering module to make clustering and contrastive learning be trained together. Experiments on multiple spatial transcriptomics datasets at different resolutions demonstrate that our method performs well in self-supervised spatial clustering tasks. Furthermore, the learned representations can be used for various downstream tasks including visualization and trajectory inference.

Controlled Noise: Evidence of epigenetic regulation of Single-Cell expression variability.

Understanding single-cell expression variability (scEV) or gene expression noise among cells of the same type and state is crucial for delineating population-level cellular function. While epigenetic mechanisms are widely implicated in gene expression regulation, a definitive link between chromatin accessibility and scEV remains elusive. Recent advances in single-cell techniques enable the study of single-cell multiomics data that include the simultaneous measurement of scATAC-seq and scRNA-seq within individual cells, presenting an unprecedented opportunity to address this gap. This paper introduces an innovative testing pipeline to investigate the association between chromatin accessibility and scEV. With single-cell multiomics data of scATAC-seq and scRNA-seq, the pipeline hinges on comparing the prediction performance of scATAC-seq data on gene expression levels between highly variable genes (HVGs) and non-highly variable genes (non-HVGs). Applying this pipeline to paired scATAC-seq and scRNA-seq data from human hematopoietic stem and progenitor cells, we observed a significantly superior prediction performance of scATAC-seq data for HVGs compared to non-HVGs. Notably, there was substantial overlap between well-predicted genes and HVGs. The gene pathways enriched from well-predicted genes are highly pertinent to cell type-specific functions. Our findings support the notion that scEV largely stems from cell-to-cell variability in chromatin accessibility, providing compelling evidence for the epigenetic regulation of scEV and offering promising avenues for investigating gene regulation mechanisms at the single-cell level. The source code and data used in this paper can be found at https://github.com/SiweiCui/EpigeneticControlOfSingle-CellExpressionVariability. Supplementary data are available at Bioinformatics online.

Stochastic Gene Expression in Proliferating Cells: Differing Noise Intensity in Single-Cell and Population Perspectives.

Random fluctuations (noise) in gene expression can be studied from two complementary perspectives: following expression in a single cell over time or comparing expression between cells in a proliferating population at a given time. Here, we systematically investigated scenarios where both perspectives lead to different levels of noise in a given gene product. We first consider a stable protein, whose concentration is diluted by cellular growth, and the protein inhibits growth at high concentrations, establishing a positive feedback loop. For a stochastic model with molecular bursting of gene products, we analytically predict and contrast the steady-state distributions of protein concentration in both frameworks. Although positive feedback amplifies the noise in expression, this amplification is much higher in the population framework compared to following a single cell over time. We also study other processes that lead to different noise levels even in the absence of such dilution-based feedback. When considering randomness in the partitioning of molecules between daughters during mitosis, we find that in the single-cell perspective, the noise in protein concentration is independent of noise in the cell cycle duration. In contrast, partitioning noise is amplified in the population perspective by increasing randomness in cell-cycle time. Overall, our results show that the commonly used single-cell framework that does not account for proliferating cells can, in some cases, underestimate the noise in gene product levels. These results have important implications for studying the inter-cellular variation of different stress-related expression programs across cell types that are known to inhibit cellular growth.

Fokker-Planck diffusion maps of multiple single cell microglial transcriptomes reveals radial differentiation into substates associated with Alzheimer's pathology.

The identification of microglia subtypes is important for understanding the role of innate immunity in neurodegenerative diseases. Current methods of unsupervised cell type identification assume a small noise-to-signal ratio of transcriptome measurements that would produce well-separated cell clusters. However, identification of subtypes is obscured by gene expression noise, diminishing the distances in transcriptome space between distinct cell types and blurring boundaries. Here we use Fokker-Planck (FP) diffusion maps to model cellular differentiation as a stochastic process whereby cells settle into local minima, corresponding to cell subtypes, in a potential landscape constructed from transcriptome data using a nearest neighbor graph approach. By applying critical transition fields, we identify individual cells on the verge of transitioning between subtypes, revealing microglial cells in inactivated, homeostatic state before radially transitioning into various specialized subtypes. Specifically, we show that cells from Alzheimer's disease patients are enriched in a microglia subtype associated to antigen presentation and T-cell recruitment.

Regulatory properties of transcription factors with diverse mechanistic function.

Transcription factors (TFs) regulate the process of transcription through the modulation of different kinetic steps. Although models can often describe the observed transcriptional output of a measured gene, predicting a TFs role on a given promoter requires an understanding of how the TF alters each step of the transcription process. In this work, we use a simple model of transcription to assess the role of promoter identity, and the degree to which TFs alter binding of RNAP (stabilization) and initiation of transcription (acceleration) on three primary characteristics: the range of steady-state regulation, cell-to-cell variability in expression, and the dynamic response time of a regulated gene. We find that steady state regulation and the response time of a gene behave uniquely for TFs that regulate incoherently, i.e that speed up one step but slow the other. We also find that incoherent TFs have dynamic implications, with one type of incoherent mode configuring the promoter to respond more slowly at intermediate TF concentrations. We also demonstrate that the noise of gene expression for these TFs is sensitive to promoter strength, with a distinct non-monotonic profile that is apparent under stronger promoters. Taken together, our work uncovers the coupling between promoters and TF regulatory modes with implications for understanding natural promoters and engineering synthetic gene circuits with desired expression properties.

Maintenance of steady-state mRNA levels by a microRNA-based feed forward loop in the presence of stochastic gene expression noise

Abstract All vital functions of living cells rely on the production of various functional molecules through gene expression. The production periods are burst-like and stochastic due to the discrete nature of biochemical reactions. In certain contexts, the concentrations of RNA or protein require regulation to maintain a fine internal balance within the cell. Here we consider a motif of two types of RNA molecules – mRNA and an antagonistic microRNA – which are encoded by a shared coding sequence and form a feed forward loop (FFL). This control mechanism is shown to be perfectly adapting in the deterministic context. We demonstrate that the adaptation (of the mean value) becomes imperfect if production occurs in random bursts. The FFL nevertheless outperforms the benchmark feedback loop in terms of counterbalancing variations in the signal. Methodologically, we adapt a hybrid stochastic model, which has widely been used to model a single regulatory molecule, to the current case of a motif involving two species; the use of the Laplace transform thereby circumvents the problem of moment closure that arises owing to the mRNA–microRNA interaction. We expect that the approach can be applicable to other systems with nonlinear kinetics.

Effect of genomic and cellular environments on gene expression noise

BackgroundIndividual cells from isogenic populations often display large cell-to-cell differences in gene expression. This “noise” in expression derives from several sources, including the genomic and cellular environment in which a gene resides. Large-scale maps of genomic environments have revealed the effects of epigenetic modifications and transcription factor occupancy on mean expression levels, but leveraging such maps to explain expression noise will require new methods to assay how expression noise changes at locations across the genome.ResultsTo address this gap, we present Single-cell Analysis of Reporter Gene Expression Noise and Transcriptome (SARGENT), a method that simultaneously measures the noisiness of reporter genes integrated throughout the genome and the global mRNA profiles of individual reporter-gene-containing cells. Using SARGENT, we perform the first comprehensive genome-wide survey of how genomic locations impact gene expression noise. We find that the mean and noise of expression correlate with different histone modifications. We quantify the intrinsic and extrinsic components of reporter gene noise and, using the associated mRNA profiles, assign the extrinsic component to differences between the CD24+ “stem-like” substate and the more “differentiated” substate. SARGENT also reveals the effects of transgene integrations on endogenous gene expression, which will help guide the search for “safe-harbor” loci.ConclusionsTaken together, we show that SARGENT is a powerful tool to measure both the mean and noise of gene expression at locations across the genome and that the data generatd by SARGENT reveals important insights into the regulation of gene expression noise genome-wide.

Noise Reduction in Resource-Coupled Multi-Module Gene Circuits through Antithetic Feedback Control.

Gene circuits within the same host cell often experience coupling, stemming from the competition for limited resources during transcriptional and translational processes. This resource competition introduces an additional layer of noise to gene expression. Here we present three multi-module antithetic control strategies: negatively competitive regulation (NCR) controller, alongside local and global controllers, aimed at reducing the gene expression noise within the context of resource competition. Through stochastic simulations and fluctuation-dissipation theorem (FDT) analysis, our findings highlight the superior performance of the NCR antithetic controller in reducing noise levels. Our research provides an effective control strategy for attenuating resource-driven noise and offers insight into the development of robust gene circuits.

MicroRNA-mediated gene regulation and the resilience of multicellular animals.

MicroRNAs are small RNAs that enable parts of the genome to regulate the other parts of the genome by RNA::RNA complementarity. Genes that encode microRNAs function as trans-acting regulators of hundreds of other genes, primarily by inhibiting the production of protein from mRNAs to which the microRNAs can bind by base pairing. MicroRNAs and their Argonaute partner proteins constitute a regulatory complex (the miRISC) that exhibits astonishing regulatory versatility. microRNAs have been shown to perform diverse roles in genetic regulatory networks (GRNs) - to control developmental switches, to dampen gene expression noise, to coordinate multigene functional modules, and more broadly, to confer robustness and resilience to developmental and homeostatic processes. Genetic analysis reveals that the function of particular microRNAs can be conditional, such that the microRNA is required under particular environmental or physiological conditions, but relatively dispensable under other conditions. The diversity and versatility of microRNA function in animal systems reflects the many ways that miRISC can be regulated by cellular signaling pathways, and the structure-function interplay among microRNA, target, and Argonaute.

Living in a noisy world-origins of gene expression noise and its impact on cellular decision-making.

The expression level of a gene can vary between genetically identical cells under the same environmental condition-a phenomenon referred to as gene expression noise. Several studies have now elucidated a central role of transcription factors in the generation of expression noise. Transcription factors, as the key components of gene regulatory networks, drive many important cellular decisions in response to cellular and environmental signals. Therefore, a very relevant question is how expression noise impacts gene regulation and influences cellular decision-making. In this Review, we summarize the current understanding of the molecular origins of expression noise, highlighting the role of transcription factors in this process, and discuss the ways in which noise can influence cellular decision-making. As advances in single-cell technologies open new avenues for studying expression noise as well as gene regulatory circuits, a better understanding of the influence of noise on cellular decisions will have important implications for many biological processes.

Feasibility of Gene Regulatory Network Inference Under Different Noise Conditions

One of the most difficult and pressing problems in computational cell biology is the inference of gene regulatory network structure from transcriptomic data. Although this topic has been widely studied for more than a decade, it is still unclear how well these methods work and how their performance differs under different conditions. We analyze the feasibility of network inference from mRNA abundance data using simulations of gene regulatory interactions, considering both intrinsic and extrinsic noise in the process of gene expression. We find that under conditions of only intrinsic noise in gene expression, the correlation between mRNA levels of genes in an activation relationship is quite low, suggesting that the task of network inference from transcriptomic data is very difficult under these conditions. By contrast, extrinsic noise affecting the process of gene expression, which could come from an upstream regulator, external stimulus, or change in the cell state, results in higher correlation between mRNA levels of these genes, potentially making the task of network inference from mRNA data more feasible. Lastly, we analyze the problem of false positives between genes that have no direct interaction but share a common upstream regulator, and explore a strategy for distinguishing between these false positives and true interactions based on noise profiles of mRNA expression levels.

Flagellar dynamics reveal fluctuations and kinetic limit in the Escherichia coli chemotaxis network

The Escherichia coli chemotaxis network, by which bacteria modulate their random run/tumble swimming pattern to navigate their environment, must cope with unavoidable number fluctuations (“noise”) in its molecular constituents like other signaling networks. The probability of clockwise (CW) flagellar rotation, or CW bias, is a measure of the chemotaxis network’s output, and its temporal fluctuations provide a proxy for network noise. Here we quantify fluctuations in the chemotaxis signaling network from the switching statistics of flagella, observed using time-resolved fluorescence microscopy of individual optically trapped E. coli cells. This approach allows noise to be quantified across the dynamic range of the network. Large CW bias fluctuations are revealed at steady state, which may play a critical role in driving flagellar switching and cell tumbling. When the network is stimulated chemically to higher activity, fluctuations dramatically decrease. A stochastic theoretical model, inspired by work on gene expression noise, points to CheY activation occurring in bursts, driving CW bias fluctuations. This model also shows that an intrinsic kinetic ceiling on network activity places an upper limit on activated CheY and CW bias, which when encountered suppresses network fluctuations. This limit may also prevent cells from tumbling unproductively in steep gradients.

Phase separation reduces cell-to-cell variability of transcriptional bursting

Gene expression is a stochastic and noisy process often occurring in "bursts". Experiments have shown that the compartmentalization of proteins by liquid-liquid phase separation is conducive to reducing the noise of gene expression. Therefore, an important goal is to explore the role of bursts in phase separation noise reduction processes. We propose a coupled model that includes phase separation and a two-state gene expression process. Using the timescale separation method, we obtain approximate solutions for the expectation, variance, and noise strength of the dilute phase. We find that a higher burst frequency weakens the ability of noise reduction by phase separation, but as the burst size increases, this ability first increases and then decreases. This study provides a deeper understanding of phase separation to reduce noise in the stochastic gene expression with burst kinetics.

Effects of microRNA-mediated negative feedback on gene expression noise

MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene expression post-transcriptionally in eukaryotes by binding with target mRNAs and preventing translation. miRNA-mediated feedback motifs are ubiquitous in various genetic networks that control cellular decision making. A key question is how such a feedback mechanism may affect gene expression noise. To answer this, we have developed a mathematical model to study the effects of a miRNA-dependent negative-feedback loop on mean expression and noise in target mRNAs. Combining analytics and simulations, we show the existence of an expression threshold demarcating repressed and expressed regimes in agreement with earlier studies. The steady-state mRNA distributions are bimodal near the threshold, where copy numbers of mRNAs and miRNAs exhibit enhanced anticorrelated fluctuations. Moreover, variation of negative-feedback strength shifts the threshold locations and modulates the noise profiles. Notably, the miRNA-mRNA binding affinity and feedback strength collectively shape the bimodality. We also compare our model with a direct auto-repression motif, where a gene produces its own repressor. Auto-repression fails to produce bimodal mRNA distributions as found in miRNA-based indirect repression, suggesting the crucial role of miRNAs in creating phenotypic diversity. Together, we demonstrate how miRNA-dependent negative feedback modifies the expression threshold and leads to a broader parameter regime of bimodality compared to the no-feedback case.

Uncovering the effect of RNA polymerase steric interactions on gene expression noise: Analytical distributions of nascent and mature RNA numbers.

The telegraph model is the standard model of stochastic gene expression, which can be solved exactly to obtain the distribution of mature RNA numbers per cell. A modification of this model also leads to an analytical distribution of nascent RNA numbers. These solutions are routinely used for the analysis of single-cell data, including the inference of transcriptional parameters. However, these models neglect important mechanistic features of transcription elongation, such as the stochastic movement of RNA polymerases and their steric (excluded-volume) interactions. Here we construct a model of gene expression describing promoter switching between inactive and active states, binding of RNA polymerases in the active state, their stochastic movement including steric interactions along the gene, and their unbinding leading to a mature transcript that subsequently decays. We derive the steady-state distributions of the nascent and mature RNA numbers in two important limiting cases: constitutive expression and slow promoter switching. We show that RNA fluctuations are suppressed by steric interactions between RNA polymerases, and that this suppression can in some instances even lead to sub-Poissonian fluctuations; these effects are most pronounced for nascent RNA and less prominent for mature RNA, since the latter is not a direct sensor of transcription. We find a relationship between the parameters of our microscopic mechanistic model and those of the standard models that ensures excellent consistency in their prediction of the first and second RNA number moments over vast regions of parameter space, encompassing slow, intermediate, and rapid promoter switching, provided the RNA number distributions are Poissonian or super-Poissonian. Furthermore, we identify the limitations of inference from mature RNA data, specifically showing that it cannot differentiate between highly distinct RNA polymerase traffic patterns on a gene.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian.

Gene expression inherently gives rise to stochastic variation ("noise") in the production of gene products. Minimizing noise is crucial for ensuring reliable cellular functions. However, noise cannot be suppressed below a certain intrinsic limit. For constitutively expressed genes, this limit is typically assumed to be Poissonian noise, wherein the variance in mRNA numbers is equal to their mean. Here, we demonstrate that several cell division genes in fission yeast exhibit mRNA variances significantly below this limit. The reduced variance can be explained by a gene expression model incorporating multiple transcription and mRNA degradation steps. Notably, in this sub-Poissonian regime, distinct from Poissonian or super-Poissonian regimes, cytoplasmic noise is effectively suppressed through a higher mRNA export rate. Our findings redefine the lower limit of eukaryotic gene expression noise and uncover molecular requirements for achieving ultralow noise, which is expected to be important for vital cellular functions.

The σB alternative sigma factor circuit modulates noise to generate different types of pulsing dynamics.

Single-cell approaches are revealing a high degree of heterogeneity, or noise, in gene expression in isogenic bacteria. How gene circuits modulate this noise in gene expression to generate robust output dynamics is unclear. Here we use the Bacillus subtilis alternative sigma factor σB as a model system for understanding the role of noise in generating circuit output dynamics. σB controls the general stress response in B. subtilis and is activated by a range of energy and environmental stresses. Recent single-cell studies have revealed that the circuit can generate two distinct outputs, stochastic pulsing and a single pulse response, but the conditions under which each response is generated are under debate. We implement a stochastic mathematical model of the σB circuit to investigate this and find that the system's core circuit can generate both response types. This is despite one response (stochastic pulsing) being stochastic in nature, and the other (single response pulse) being deterministic. We demonstrate that the main determinant for whichever response is generated is the degree with which the input pathway activates the core circuit, although the noise properties of the input pathway also biases the system towards one or the other type of output. Thus, our work shows how stochastic modelling can reveal the mechanisms behind non-intuitive gene circuit output dynamics.