Origins of binary gene expression in post-transcriptional regulation by microRNAs

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.

Emergent correlations in gene expression dynamics as footprints of resource competition

Genetic circuits need a cellular environment to operate in, which naturally couples the circuit function with the overall functionality of gene regulatory network. To execute their functions, all gene circuits draw resources in the form of RNA polymerases, ribosomes, and tRNAs. Recent experiments pointed out that the role of resource competition on synthetic circuit outputs could be immense. However, the effect of complexity of the circuit architecture on resource sharing dynamics is yet unexplored. In this paper, we employ mathematical modelling and in-silico experiments to identify the sources of resource trade-off and to quantify its impact on the function of a genetic circuit, keeping our focus on regulation of immediate downstream proteins, which are often used as protein read-outs. We show that estimating gene expression dynamics from readings of downstream protein data might be unreliable when the resource is limited and ribosome affinities are asymmetric. We focus on the impact of mRNA copy number and ribosome binding site (RBS) strength on the nonlinear isocline that emerges with two regimes, prominently separated by a tipping point, and study how correlation and competition dominate each other depending on various circuit parameters. Focusing further on genetic toggle circuit, we have identified major effects of resource competition in this model motif and quantified the observations. The observations are testable in wet-lab experiments, as all the parameters chosen are experimentally relevant.

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.

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

Several studies highlighted the relevance of extrinsic noise in shaping cell decision making and differentiation in molecular networks. Bimodal distributions of gene expression levels provide experimental evidence of phenotypic differentiation, where the modes of the distribution often correspond to different physiological states of the system. We theoretically address the presence of bimodal phenotypes in the context of microRNA (miRNA)-mediated regulation. MiRNAs are small noncoding RNA molecules that downregulate the expression of their target mRNAs. The nature of this interaction is titrative and induces a threshold effect: below a given target transcription rate almost no mRNAs are free and available for translation. We investigate the effect of extrinsic noise on the system by introducing a fluctuating miRNA-transcription rate. We find that the presence of extrinsic noise favours the presence of bimodal target distributions which can be observed for a wider range of parameters compared to the case with intrinsic noise only and for lower miRNA-target interaction strength. Our results suggest that combining threshold-inducing interactions with extrinsic noise provides a simple and robust mechanism for obtaining bimodal populations without requiring fine tuning. Furthermore, we characterise the protein distribution's dependence on protein half-life.