Quantifying mRNA synthesis and decay rates using small RNAs

Tissue- and Species-Specific Patterns of RNA metabolism in Post-Mortem Mammalian Retina and Retinal Pigment Epithelium

Accurate analysis of gene expression in human tissues using RNA sequencing is dependent on the quality of source material. One major source of variation in mRNA quality is post-mortem time. While it is known that individual transcripts show differential post-mortem stability, few studies have directly and comprehensively analyzed mRNA stability following death, and in particular the extent to which tissue- and species-specific factors influence post-mortem mRNA stability are poorly understood. This knowledge is particularly important for ocular tissues studies, where tissues obtained post-mortem are frequently used for research or therapeutic applications. To directly investigate this question, we profiled mRNA levels in both neuroretina and retinal pigment epithelium (RPE) from mouse and baboon over a series of post-mortem intervals. We found substantial changes in gene expression as early as 15 minutes in the mouse and as early as three hours in the baboon eye tissues. Importantly, our findings demonstrate both tissue- and species- specific patterns of RNA metabolism, by identifying a set of genes that are either rapidly degraded or very stable in both species and/or tissues. Taken together, the data from this study lay the foundation for understanding RNA regulation post-mortem and provide novel insights into RNA metabolism in the tissues of the mammalian eye.

Quantitative Characteristic of ncRNA Regulation in Gene Regulatory Networks

RNA is mostly known for its role in protein synthesis, where it encodes information for protein sequence in its messenger RNA (mRNA) form (translation). Yet, RNA molecules regulate several cellular processes other than translation. Here, we present an overview of several mathematical models that help understanding and characterizing the role of noncoding RNA molecules (ncRNAs) in regulating gene expression and protein synthesis. First, we discuss relatively simple models where ncRNAs can modulate protein synthesis via targeting a mRNA. Then, we consider the case of feedback interactions between ncRNAs and their target proteins, and discuss several biological applications where these feedback architectures modulate a cellular phenotype and control the levels of intrinsic and extrinsic noise. Building from these simple circuit motifs, we examine feed-forward circuit motifs involving ncRNAs that generate precise spatial and temporal patterns of protein expression. Further, we investigate the competition between ncRNAs and other endogenous RNA molecules and show that the cross talk between coding and noncoding RNAs can form large genetic circuits that involve up to hundreds of chemical species. Finally, we discuss the role of ncRNAs in modulating cell-cell signaling pathways and therefore the dynamics of spatiotemporal pattern formation in a tissue.

Intrinsic noise of microRNA-regulated genes and the ceRNA hypothesis

MicroRNAs are small noncoding RNAs that regulate genes post-transciptionally by binding and degrading target eukaryotic mRNAs. We use a quantitative model to study gene regulation by inhibitory microRNAs and compare it to gene regulation by prokaryotic small non-coding RNAs (sRNAs). Our model uses a combination of analytic techniques as well as computational simulations to calculate the mean-expression and noise profiles of genes regulated by both microRNAs and sRNAs. We find that despite very different molecular machinery and modes of action (catalytic vs stoichiometric), the mean expression levels and noise profiles of microRNA-regulated genes are almost identical to genes regulated by prokaryotic sRNAs. This behavior is extremely robust and persists across a wide range of biologically relevant parameters. We extend our model to study crosstalk between multiple mRNAs that are regulated by a single microRNA and show that noise is a sensitive measure of microRNA-mediated interaction between mRNAs. We conclude by discussing possible experimental strategies for uncovering the microRNA-mRNA interactions and testing the competing endogenous RNA (ceRNA) hypothesis.

Complex degradation processes lead to non-exponential decay patterns and age-dependent decay rates of messenger RNA

Experimental studies on mRNA stability have established several, qualitatively distinct decay patterns for the amount of mRNA within the living cell. Furthermore, a variety of different and complex biochemical pathways for mRNA degradation have been identified. The central aim of this paper is to bring together both the experimental evidence about the decay patterns and the biochemical knowledge about the multi-step nature of mRNA degradation in a coherent mathematical theory. We first introduce a mathematical relationship between the mRNA decay pattern and the lifetime distribution of individual mRNA molecules. This relationship reveals that the mRNA decay patterns at steady state expression level must obey a general convexity condition, which applies to any degradation mechanism. Next, we develop a theory, formulated as a Markov chain model, that recapitulates some aspects of the multi-step nature of mRNA degradation. We apply our theory to experimental data for yeast and explicitly derive the lifetime distribution of the corresponding mRNAs. Thereby, we show how to extract single-molecule properties of an mRNA, such as the age-dependent decay rate and the residual lifetime. Finally, we analyze the decay patterns of the whole translatome of yeast cells and show that yeast mRNAs can be grouped into three broad classes that exhibit three distinct decay patterns. This paper provides both a method to accurately analyze non-exponential mRNA decay patterns and a tool to validate different models of degradation using decay data.

Regulation by small RNAs via coupled degradation: mean-field and variational approaches

Regulatory genes called small RNAs (sRNAs) are known to play critical roles in cellular responses to changing environments. For several sRNAs, regulation is effected by coupled stoichiometric degradation with messenger RNAs (mRNAs). The nonlinearity inherent in this regulatory scheme indicates that exact analytical solutions for the corresponding stochastic models are intractable. Here, we present a variational approach to analyze a well-studied stochastic model for regulation by sRNAs via coupled degradation. The proposed approach is efficient and provides accurate estimates of mean mRNA levels as well as higher-order terms. Results from the variational ansatz are in excellent agreement with data from stochastic simulations for a wide range of parameters, including regions of parameter space where mean-field approaches break down. The proposed approach can be applied for quantitative modeling of stochastic gene expression in complex regulatory networks.

Connecting protein and mRNA burst distributions for stochastic models of gene expression

The intrinsic stochasticity of gene expression can lead to large variability in protein levels for genetically identical cells. Such variability in protein levels can arise from infrequent synthesis of mRNAs which in turn give rise to bursts of protein expression. Protein expression occurring in bursts has indeed been observed experimentally and recent studies have also found evidence for transcriptional bursting, i.e. production of mRNAs in bursts. Given that there are distinct experimental techniques for quantifying the noise at different stages of gene expression, it is of interest to derive analytical results connecting experimental observations at different levels. In this work, we consider stochastic models of gene expression for which mRNA and protein production occurs in independent bursts. For such models, we derive analytical expressions connecting protein and mRNA burst distributions which show how the functional form of the mRNA burst distribution can be inferred from the protein burst distribution. Additionally, if gene expression is repressed such that observed protein bursts arise only from single mRNAs, we show how observations of protein burst distributions (repressed and unrepressed) can be used to completely determine the mRNA burst distribution. Assuming independent contributions from individual bursts, we derive analytical expressions connecting means and variances for burst and steady-state protein distributions. Finally, we validate our general analytical results by considering a specific reaction scheme involving regulation of protein bursts by small RNAs. For a range of parameters, we derive analytical expressions for regulated protein distributions that are validated using stochastic simulations. The analytical results obtained in this work can thus serve as useful inputs for a broad range of studies focusing on stochasticity in gene expression.

Hierarchical genetic networks and noncoding RNAs

In eukaryotic cells, many genes are transcribed into noncoding RNAs. Such RNAs may associate with mRNAs and inhibit their translation and facilitate degradation. To clarify what may happen in this case, we propose a kinetic model describing the effect of noncoding RNAs on a mRNA-protein network with the hierarchical three-layer architecture. For positive regulation of the layers, our model predicts either bistability with a fairly narrow hysteresis loop or a unique steady state. For negative or mixed regulation, the steady state is found to be unique.

Evidence for an autonomous 5' target recognition domain in an Hfq-associated small RNA

The abundant class of bacterial Hfq-associated small regulatory RNAs (sRNAs) parallels animal microRNAs in their ability to control multiple genes at the posttranscriptional level by short and imperfect base pairing. In contrast to the universal length and seed pairing mechanism of microRNAs, the sRNAs are heterogeneous in size and structure, and how they regulate multiple targets is not well understood. This paper provides evidence that a 5' located sRNA domain is a critical element for the control of a large posttranscriptional regulon. We show that the conserved 5' end of RybB sRNA recognizes multiple mRNAs of Salmonella outer membrane proteins by ≥7-bp Watson-Crick pairing. When fused to an unrelated sRNA, the 5' domain is sufficient to guide target mRNA degradation and maintain σ(E)-dependent envelope homeostasis. RybB sites in mRNAs are often conserved and flanked by 3' adenosine. They are found in a wide sequence window ranging from the upstream untranslated region to the deep coding sequence, indicating that some targets might be repressed at the level of translation, whereas others are repressed primarily by mRNA destabilization. Autonomous 5' domains seem more common in sRNAs than appreciated and might improve the design of synthetic RNA regulators.