On the role of extrinsic noise in microRNA-mediated bimodal gene expression

Bioenergetic costs and the evolution of noise regulation by microRNAs

Noise control, together with other regulatory functions facilitated by microRNAs (miRNAs), is believed to have played important roles in the evolution of multicellular eukaryotic organisms. miRNAs can dampen protein fluctuations via enhanced degradation of messenger RNA (mRNA), but this requires compensation by increased mRNA transcription to maintain the same expression levels. The overall mechanism is metabolically expensive, leading to questions about how it might have evolved in the first place. We develop a stochastic model of miRNA noise regulation, coupled with a detailed analysis of the associated metabolic costs. Additionally, we calculate binding free energies for a range of miRNA seeds, the short sequences which govern target recognition. We argue that natural selection may have fine-tuned the Michaelis-Menten constant [Formula: see text] describing miRNA-mRNA affinity and show supporting evidence from analysis of experimental data. [Formula: see text] is constrained by seed length, and optimal noise control (minimum protein variance at a given energy cost) is achievable for seeds of 6 to 7 nucleotides in length, the most commonly observed types. Moreover, at optimality, the degree of noise reduction approaches the theoretical bound set by the Wiener-Kolmogorov linear filter. The results illustrate how selective pressure toward energy efficiency has potentially shaped a crucial regulatory pathway in eukaryotes.

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.

An integrated approach to identify bimodal genes associated with prognosis in câncer

Bimodal gene expression (where a gene expression distribution has two maxima) is associated with phenotypic diversity in different biological systems. A critical issue, thus, is the integration of expression and phenotype data to identify genuine associations. Here, we developed tools that allow both: i) the identification of genes with bimodal gene expression and ii) their association with prognosis in cancer patients from The Cancer Genome Atlas (TCGA). Bimodality was observed for 554 genes in expression data from 25 tumor types. Furthermore, 96 of these genes presented different prognosis when patients belonging to the two expression peaks were compared. The software to execute the method and the corresponding documentation are available at the Data access section.

Single-cell Transcriptomes Reveal Characteristics of MicroRNA 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.

One-to-one relationships between milk miRNA content and protein abundance in neonate duodenum support the potential for milk miRNAs regulating neonate development

Breast milk plays an essential role for offspring development; however, there lacks evidence of how specific milk components like nucleic acids mechanistically function to regulate neonate development. Previously, we found that maternal high-fat diet (HFD) not only significantly affected mRNA and miRNA content of the secreted milk transcriptome in mice but also affected the duodenal proteome of suckling pups. Here, we hypothesized that nucleic acids differentially expressed in milk of HFD fed dams are related to differentially abundant proteins in offspring duodenum nursed by HFD dams. We tested this hypothesis by analyzing one-to-one relationships in RNA-seq data of milk transcriptomes from control (10% kcal fat) and HFD (60% kcal fat) fed mice and liquid chromatography-tandem mass spectrometry (LC-MS/MS) duodenal proteome data from pups exposed to milk. Ten percent of differentially abundant duodenal proteins between controls and HFD-exposed pups had predicted upregulation or downregulation based on differential milk RNA content. Of these, 76% were targets of upregulated miRNA, and linear regression analysis indicated relationships (p < 0.05) between multiple milk miRNA counts and duodenal protein abundance. Duodenal proteins that were potential targets of milk miRNA enriched Gene Ontology (GO) terms and KEGG pathways related to cytoskeletal structure and neural development, suggesting potential regulation of pup enteric nervous system. One-to-one relationships between milk miRNA content and protein abundance in neonate duodenum support the potential for milk miRNAs regulating neonate development. Identification of milk miRNAs that changed in response to maternal diet will enable design of mechanistic studies that test effects on neonate.

microRNA-mediated noise processing in cells: A fight or a game?

In the past decades, microRNAs (miRNA) have much attracted the attention of researchers at the interface between life and theoretical sciences for their involvement in post-transcriptional regulation and related diseases. Thanks to the always more sophisticated experimental techniques, the role of miRNAs as "noise processing units" has been further elucidated and two main ways of miRNA noise-control have emerged by combinations of theoretical and experimental studies. While on one side miRNAs were thought to buffer gene expression noise, it has recently been suggested that miRNAs could also increase the cell-to-cell variability of their targets. In this Mini Review, we focus on the role of miRNAs in molecular noise processing and on the advantages as well as current limitations of theoretical modelling.

MicroRNAs organize intrinsic variation into stem cell states

Understanding how mammalian organisms achieve the full diversity of cell types in the adult organism is a central goal of developmental cell biology. Recent work has shown that some embryonic precursor cells can self-organize into developmental structures but the mechanisms of gene regulation that contribute to this process remain unknown. Here we show embryonic stem cells self-organize into distinct gene expression states that resemble developmental gene programs. We find that microRNAs, small noncoding regulators of gene expression, play a critical role in organizing fluctuations across gene networks to help achieve this organization into distinct expression states.

Pluripotent embryonic stem cells (ESCs) contain the potential to form a diverse array of cells with distinct gene expression states, namely the cells of the adult vertebrate. Classically, diversity has been attributed to cells sensing their position with respect to external morphogen gradients. However, an alternative is that diversity arises in part from cooption of fluctuations in the gene regulatory network. Here we find ESCs exhibit intrinsic heterogeneity in the absence of external gradients by forming interconverting cell states. States vary in developmental gene expression programs and display distinct activity of microRNAs (miRNAs). Notably, miRNAs act on neighborhoods of pluripotency genes to increase variation of target genes and cell states. Loss of miRNAs that vary across states reduces target variation and delays state transitions, suggesting variable miRNAs organize and propagate variation to promote state transitions. Together these findings provide insight into how a gene regulatory network can coopt variation intrinsic to cell systems to form robust gene expression states. Interactions between intrinsic heterogeneity and environmental signals may help achieve developmental outcomes.

The contribution of microRNA-mediated regulation to short- and long-term gene expression predictability

MicroRNAs are a class of short, noncoding RNAs which are essential for the coordination and timing of cell differentiation and embryonic development. However, despite their guiding role in development, microRNAs are dysregulated in many pathologies, including nearly all cases of cancer. While both development and oncogenesis can be thought of as extremes of phenotypic plasticity, they characteristically manifest on much different time scales: one taking place over a matter of weeks, the other typically requiring decades. Because microRNAs are believed to support this plasticity, a critically important question is how microRNAs affect phenotypic stability on different time scales, and what dynamical characteristics shift the balance between these two roles. To address this question, we extend a well-established mathematical model of transcriptional gene regulation to include translational regulation by microRNAs, and examine their effects on both short- and long-term gene expression predictability. Our findings show that microRNAs greatly improve short-term predictability for earlier, developmental phenotypes while causing a small decrease in long-term predictability, and that these effects are difficult to separate. In addition to providing a theoretical explanation for this seemingly duplicitous behavior, we describe some of the properties which determine the cost-benefit balance between short-term stabilization and long-term destabilization.

From Endogenous to Synthetic microRNA-Mediated Regulatory Circuits: An Overview

MicroRNAs are short non-coding RNAs that are evolutionarily conserved and are pivotal post-transcriptional mediators of gene regulation. Together with transcription factors and epigenetic regulators, they form a highly interconnected network whose building blocks can be classified depending on the number of molecular species involved and the type of interactions amongst them. Depending on their topology, these molecular circuits may carry out specific functions that years of studies have related to the processing of gene expression noise. In this review, we first present the different over-represented network motifs involving microRNAs and their specific role in implementing relevant biological functions, reviewing both theoretical and experimental studies. We then illustrate the recent advances in synthetic biology, such as the construction of artificially synthesised circuits, which provide a controlled tool to test experimentally the possible microRNA regulatory tasks and constitute a starting point for clinical applications.

Intrinsic and extrinsic noise of gene expression in lineage trees

Cell-to-cell heterogeneity is driven by stochasticity in intracellular reactions and the population dynamics. While these sources are usually studied separately, we develop an agent-based framework that accounts for both factors while tracking every single cell of a growing population. Apart from the common intrinsic variability, the framework also predicts extrinsic noise without the need to introduce fluctuating rate constants. Instead, extrinsic fluctuations are explained by cell cycle fluctuations and differences in cell age. We provide explicit formulas to quantify mean molecule numbers, intrinsic and extrinsic noise statistics in two-colour experiments. We find that these statistics differ significantly depending on the experimental setup used to observe the cells. We illustrate this fact using (i) averages over an isolated cell lineage tracked over many generations as observed in the mother machine, (ii) population snapshots with known cell ages as recorded in time-lapse microscopy, and (iii) snapshots with unknown cell ages as measured from static images or flow cytometry. Applying the method to models of stochastic gene expression and feedback regulation elucidates that isolated lineages, as compared to snapshot data, can significantly overestimate the mean number of molecules, overestimate extrinsic noise but underestimate intrinsic noise and have qualitatively different sensitivities to cell cycle fluctuations.

Kinetic Modelling of Competition and Depletion of Shared miRNAs by Competing Endogenous RNAs

Non-coding RNAs play a key role in the post-transcriptional regulation of mRNA translation and turnover in eukaryotes. miRNAs, in particular, interact with their target RNAs through protein-mediated, sequence-specific binding, giving rise to extended and highly heterogeneous miRNA-RNA interaction networks. Within such networks, competition to bind miRNAs can generate an effective positive coupling between their targets. Competing endogenous RNAs (ceRNAs) can in turn regulate each other through miRNA-mediated crosstalk. Albeit potentially weak, ceRNA interactions can occur both dynamically, affecting, e.g., the regulatory clock, and at stationarity, in which case ceRNA networks as a whole can be implicated in the composition of the cell's proteome. Many features of ceRNA interactions, including the conditions under which they become significant, can be unraveled by mathematical and in silico models. We review the understanding of the ceRNA effect obtained within such frameworks, focusing on the methods employed to quantify it, its role in the processing of gene expression noise, and how network topology can determine its reach.

Stochastic sequestration dynamics: a minimal model with extrinsic noise for bimodal distributions and competitors correlation

Many biological processes are known to be based on molecular sequestration. This kind of dynamics involves two types of molecular species, namely targets and sequestrants, that bind to form a complex. In the simple framework of mass-action law, key features of these systems appear to be threshold-like profiles of the amounts of free molecules as a function of the parameters determining their possible maximum abundance. However, biochemical processes are probabilistic and take place in stochastically fluctuating environments. How these different sources of noise affect the final outcome of the network is not completely characterised yet. In this paper we specifically investigate the effects induced by a source of extrinsic noise onto a minimal stochastic model of molecular sequestration. We analytically show how bimodal distributions of the targets can appear and characterise them as a result of noise filtering mediated by the threshold response. We then address the correlations between target species induced by the sequestrant and discuss how extrinsic noise can turn the negative correlation caused by competition into a positive one. Finally, we consider the more complex scenario of competitive inhibition for enzymatic kinetics and discuss the relevance of our findings with respect to applications.