Stochastic Expression Research Articles

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.

Stress-Induced Constraint on Expression Noise of Essential Genes in E. coli.

Gene expression is an inherently noisy process that is constrained by natural selection. Yet the condition dependence of constraint on expression noise remains unclear. Here, we address this problem by studying constraint on expression noise of E. coli genes in eight diverse growth conditions. In particular, we use variation in expression noise as an analog for constraint, examining its relationships to expression level and to the number of regulatory inputs from transcription factors across and within conditions. We show that variation in expression noise is negatively associated with expression level, implicating constraint to minimize expression noise of highly expressed genes. However, this relationship is condition dependent, with the strongest constraint observed when E. coli are grown in the presence of glycerol or ciprofloxacin, which result in carbon or antibiotic stress, respectively. In contrast, we do not observe evidence of constraint on expression noise of highly regulated genes, suggesting that highly expressed and highly regulated genes represent distinct classes of genes. Indeed, we find that essential genes are often highly expressed but not highly regulated, with elevated expression noise in glycerol and ciprofloxacin conditions. Thus, our findings support the hypothesis that selective constraint on expression noise is condition dependent in E. coli, illustrating how it may play a critical role in ensuring expression stability of essential genes in unstable environments.

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction.

Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1 , which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.

Exploiting phenotypic heterogeneity to improve production of glutathione by yeast

BackgroundGene expression noise (variation in gene expression among individual cells of a genetically uniform cell population) can result in heterogenous metabolite production by industrial microorganisms, with cultures containing both low- and high-producing cells. The presence of low-producing individuals may be a factor limiting the potential for high yields. This study tested the hypothesis that low-producing variants in yeast cell populations can be continuously counter-selected, to increase net production of glutathione (GSH) as an exemplar product.ResultsA counter-selection system was engineered in Saccharomyces cerevisiae based on the known feedback inhibition of gamma-glutamylcysteine synthetase (GSH1) gene expression, which is rate limiting for GSH synthesis: the GSH1 ORF and the counter-selectable marker GAP1 were expressed under control of the TEF1 and GSH-regulated GSH1 promoters, respectively. An 18% increase in the mean cellular GSH level was achieved in cultures of the engineered strain supplemented with D-histidine to counter-select cells with high GAP1 expression (i.e. low GSH-producing cells). The phenotype was non-heritable and did not arise from a generic response to D-histidine, unlike that with certain other test-constructs prepared with alternative markers.ConclusionsThe results corroborate that the system developed here improves GSH production by targeting low-producing cells. This supports the potential for exploiting end-product/promoter interactions to enrich high-producing cells in phenotypically heterogeneous populations, in order to improve metabolite production by yeast.

Mitochondrial protein heterogeneity stems from the stochastic nature of co-translational protein targeting in cell senescence

A decline in mitochondrial function is a hallmark of aging and neurodegenerative diseases. It has been proposed that changes in mitochondrial morphology, including fragmentation of the tubular mitochondrial network, can lead to mitochondrial dysfunction, yet the mechanism of this loss of function is unclear. Most proteins contained within mitochondria are nuclear-encoded and must be properly targeted to the mitochondria. Here, we report that sustained mRNA localization and co-translational protein delivery leads to a heterogeneous protein distribution across fragmented mitochondria. We find that age-induced mitochondrial fragmentation drives a substantial increase in protein expression noise across fragments. Using a translational kinetic and molecular diffusion model, we find that protein expression noise is explained by the nature of stochastic compartmentalization and that co-translational protein delivery is the main contributor to increased heterogeneity. We observed that cells primarily reduce the variability in protein distribution by utilizing mitochondrial fission-fusion processes rather than relying on the mitophagy pathway. Furthermore, we are able to reduce the heterogeneity of the protein distribution by inhibiting co-translational protein targeting. This research lays the framework for a better understanding of the detrimental impact of mitochondrial fragmentation on the physiology of cells in aging and disease.

Rational Design of High-Efficiency Synthetic Small Regulatory RNAs and Their Application in Robust Genetic Circuit Performance Through Tight Control of Leaky Gene Expression.

Synthetic sRNAs show promise as tools for targeted and programmable gene expression manipulation. However, the design of high-efficiency synthetic sRNAs is a challenging task that necessitates careful consideration of multiple factors. Therefore, this study aims to investigate rational design strategies that significantly and robustly enhance the efficiency of synthetic sRNAs. This is achieved by optimizing the following parameters: the sRNA scaffold, mRNA binding affinity, Hfq protein expression level, and mRNA secondary structure. By utilizing optimized synthetic sRNAs within a positive feedback circuit, we effectively addressed the issue of gene expression leakage─an enduring challenge in synthetic biology that undermines the reliability of genetic circuits in bacteria. Our designed synthetic sRNAs successfully prevented gene expression leakage, thus averting unintended circuit activation caused by initial expression noise, even in the absence of signal molecules. This result shows that high-efficiency synthetic sRNAs not only enable precise gene knockdown for metabolic engineering but also ensure the robust performance of synthetic circuits. The strategies developed here hold significant promise for broad applications across diverse biotechnological fields, establishing synthetic sRNAs as pivotal tools in advancing synthetic biology and gene regulation.

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.

Histone deacetylase Hos2 regulates protein expression noise by potentially modulating the protein translation machinery.

Non-genetic variations derived from expression noise at transcript or protein levels can result in cell-to-cell heterogeneity within an isogenic population. Although cells have developed strategies to reduce noise in some cellular functions, this heterogeneity can also facilitate varying levels of regulation and provide evolutionary benefits in specific environments. Despite several general characteristics of cellular noise having been revealed, the detailed molecular pathways underlying noise regulation remain elusive. Here, we established a dual-fluorescent reporter system in Saccharomyces cerevisiae and performed experimental evolution to search for mutations that increase expression noise. By analyzing evolved cells using bulk segregant analysis coupled with whole-genome sequencing, we identified the histone deacetylase Hos2 as a negative noise regulator. A hos2 mutant down-regulated multiple ribosomal protein genes and exhibited partially compromised protein translation, indicating that Hos2 may regulate protein expression noise by modulating the translation machinery. Treating cells with translation inhibitors or introducing mutations into several Hos2-regulated ribosomal protein genes-RPS9A, RPS28B and RPL42A-enhanced protein expression noise. Our study provides an effective strategy for identifying noise regulators and also sheds light on how cells regulate non-genetic variation through protein translation.

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.

Cis-regulatory control of transcriptional timing and noise in response to estrogen

Cis-regulatory elements control transcription levels, temporal dynamics, and cell-cell variation or transcriptional noise. However, the combination of regulatory features that control these different attributes is not fully understood. Here, we used single-cell RNA-seq during an estrogen treatment time course and machine learning to identify predictors of expression timing and noise. We found that genes with multiple active enhancers exhibit faster temporal responses. We verified this finding by showing that manipulation of enhancer activity changes the temporal response of estrogen target genes. Analysis of transcriptional noise uncovered a relationship between promoter and enhancer activity, with active promoters associated with low noise and active enhancers linked to high noise. Finally, we observed that co-expression across single cells is an emergent property associated with chromatin looping, timing, and noise. Overall, our results indicate a fundamental tradeoff between a gene's ability to quickly respond to incoming signals and maintain low variation across cells.

Molecular-Memory-Induced Counter-Intuitive Noise Attenuator in Protein Polymerization

Gene expression comprises many asymmetric and complex processes. Transcriptional details revealed by the whole genome indicate that genes resort to transcriptional bursting and accumulate molecular memory. However, it is still unclear how the interplay of transcriptional bursting and memory regulates robustness and expression noise. Here, we consider a model of multiple coupled processes of protein polymerization to focus on decoding the effect of molecular memory. Using non-Markovian transformation technology, we first define the memory index to measure the correlation window of expression to decipher the mechanism of regulation. The results indicate that memory from synthesis can amplify expression noise, while memory originating from polymerization can reduce the lower bound of the noise of gene products; that is, the memory from different sources plays distinct regulatory roles to induce non-symmetry. Moreover, it is counterintuitive that the dual regulation from memory and bursting expression can directly suppress system noise, violating the principle that transcriptional bursting enhances noise. Our results not only provide a theoretical framework for investigating the function of memory but also imply that expression noise is not part of a half-power relationship with, nor mediated by, memory.

A dual fluorescence protein expression system detects cell cycle dependent protein noise

Inherently identical cells exhibit significant phenotypic variation. It can be essential for many biological processes and is known to arise from stochastic, ‘noisy’, gene expression that is determined by intrinsic and extrinsic components. It is now obvious that the noise varies as a function of inducer concentration. However, its fluctuation over the cell cycle is limited. Applying dual colour fluorescence protein reporter system, Cyan Fluorescent Protein (CFP) and Yellow fluorescent protein (YFP) tagged multi-copy plasmids, we determine variation of the noise components over the phases in lac promoter induced by Isopropyl β-D-1-thiogalactopyranoside (IPTG) and in presence of additional Magnesium, Mg2+ ion. We, also, estimate the how such system deviates from observations of single-copy plasmid. Found 25 % difference between multi-copy system and single-copy system clarifies that observed noise is considerable and estimates population behaviour during the cell cycle. We show that total variation in cells induced with IPTG is determined by higher extrinsic than intrinsic noise. It increases from Lag to Exponential phase and decreases from Retardation to Stationary phase. By observing slow and fast dividing cells, we show that 5 mM Mg2+ increases population homogeneity compared to 2.5 mM Mg2+ in the environment. The experimental data obtained using dual colour fluorescence protein reporter system demonstrates that protein expression noise, depending on intra cellular ionic concentration, is tightly controlled by phase of the cell.