Gene autoregulation via intronic microRNAs and its functions
- Carla Bosia†1, 2Email author,
- Matteo Osella†3, 4Email author,
- Mariama El Baroudi5,
- Davide Corà2, 6 and
- Michele Caselle2, 7
© Bosia et al.; licensee BioMed Central Ltd. 2012
Received: 11 June 2012
Accepted: 1 October 2012
Published: 10 October 2012
MicroRNAs, post-transcriptional repressors of gene expression, play a pivotal role in gene regulatory networks. They are involved in core cellular processes and their dysregulation is associated to a broad range of human diseases. This paper focus on a minimal microRNA-mediated regulatory circuit, in which a protein-coding gene (host gene) is targeted by a microRNA located inside one of its introns.
Autoregulation via intronic microRNAs is widespread in the human regulatory network, as confirmed by our bioinformatic analysis, and can perform several regulatory tasks despite its simple topology. Our analysis, based on analytical calculations and simulations, indicates that this circuitry alters the dynamics of the host gene expression, can induce complex responses implementing adaptation and Weber’s law, and efficiently filters fluctuations propagating from the upstream network to the host gene. A fine-tuning of the circuit parameters can optimize each of these functions. Interestingly, they are all related to gene expression homeostasis, in agreement with the increasing evidence suggesting a role of microRNA regulation in conferring robustness to biological processes. In addition to model analysis, we present a list of bioinformatically predicted candidate circuits in human for future experimental tests.
The results presented here suggest a potentially relevant functional role for negative self-regulation via intronic microRNAs, in particular as a homeostatic control mechanism of gene expression. Moreover, the map of circuit functions in terms of experimentally measurable parameters, resulting from our analysis, can be a useful guideline for possible applications in synthetic biology.
microRNAs (miRNAs) are small (about 22 nucleotides) single-strand RNAs able to interfere post-transcriptionally with the protein production of their targets. Targeting a vast proportion of protein-coding genes[1–3], miRNA-mediated regulation composes an important layer in gene regulatory networks. The implication of miRNAs in several core cellular processes[4–7] as well as in many human diseases[8, 9] further confirms their biological importance.
Approximately half of the miRNA genes can be found in intergenic regions (between genes), whereas the intragenic miRNAs (inside genes) are predominantly located inside introns and usually oriented on the same DNA strand of the host gene (a trend further confirmed by our bioinformatic analysis shown in a following section). Intergenic miRNA genes present their own promoter region[11, 12] and their expression is expected to be regulated by the same molecular mechanisms that control the expression of protein-coding genes. On the other hand, experimental and computational results are consistent with the idea that same-strand intronic miRNAs are co-transcribed with their host gene[13–17], and then processed to finally become mature functional miRNAs[18, 19] (although exceptions to this common scheme of co-transcripton have been reported[20–22]).
The host-miRNA co-expression can have a specific functional role. In fact, an intronic miRNA can support the function of its host gene by silencing genes that are functionally antagonistic to the host, or more generally act synergistically with the host by coordinating the expression of genes with related functions.
In addition to this “cooperative” miRNA-host relation, different studies showed that intronic miRNAs can directly regulate the expression of their host gene, establishing a negative feedback regulation[10, 25, 26]. In particular, instances of negative autoregulatory feedbacks via intronic miRNAs were firstly found by expression analysis in human. More recently, two independent large-scale bioinformatic analysis, based on different algorithms of target prediction, claimed that the occurency of intronic miRNA-mediated self-loops (iMSLs) in the human regulatory network is significantly higher than expected by chance[10, 26]. The over-representation of such regulatory module can be interpreted as a sign of evolutionary positive selection that has led to an accumulation of a specific topology able to perform useful elementary regulatory tasks. In addition, two iMSL circuits have been confirmed experimentally: regulation of EGFL7 by its intronic miRNA miR-126[28, 29] and regulation of ARPP-21 by miR-128b. Both regulations were associated to relevant biological functions, the former playing a role in cancer proliferation, while the latter in vertebrate brain physiology.
The combination of all these pieces of evidence suggests that iMSLs are an often exploited and presumably functionally relevant regulatory circuitry. The open question concerns the peculiar functions that an iMSL can accomplish and that could have thus driven their pervasive spreading in the human regulatory network. Moreover, it would be interesting to understand what specificities of post-transcriptional autoregulation by miRNAs can make them better suited to fullfil certain tasks with respect to the trascriptional self-regulation, so widely used in bacteria. In this paper we address these questions by modeling the dynamical and stochastic behaviour of the iMSL circuit and comparing its properties to those of alternative regulatory strategies such as constitutive expression and transcriptional self repression.
Our results show that, despite of its minimal topology, the iMSL circuitry can implement different biological functions. It can speed-up the host gene protein production in response to an activating signal, while delaying its switching-off kinetics when the activation drops; it can buffer fluctuations coming from the upstream network, and generate complex behaviours like a host gene expression response obeying “Weber’s law” (i.e. the magnitude of the response depends only on the fold change of the input signal). While these different functions can be optimized individually, by tuning parameters like molecular production/degradation rates, it will be shown that they all represent different ways of making the host gene expression robust to external fluctuations. Therefore, autoregulation via intronic miRNAs can generally represent an efficient homeostatic control of the host gene expression, in agreement with the observation that miRNAs are often involved in signaling networks to ensure homeostasis and gene expression robustness[31–34]. In addition to model analysis, we present our own bioinformatical search for iMSLs in human to further assess their statistical over-representation and to propose the best predicted candidates for possible future experimental tests.
Besides the understanding of the role of endogenous iMSLs, our results can be useful for the growing field of synthetic biology[35, 36], which has succesfully started to make use of RNA-based post-transcriptional regulations[37, 38]. The function-topology map presented in this paper can contribute to draw up the manual of biological circuits that carry out specific functions for synthetic engineering, adding a simple and efficient wiring strategy that can increase systems’ robustness in different conditions. A synthetic realization of an iMSL has been indeed recently produced and proven to be effective in reducing the expression dependency on gene dosage. Therefore, the potential additional functions we will show associated to iMSLs could be tested in the near future.
Results and discussion
Outline of the model
On the other hand, miRNAs can exert their action repressing translation or inducing degradation of their target mRNAs. We construct our model supposing an action on target translation. While most of the results shown in this paper are independent of this choice, some dynamical properties of the circuit can actually change if miRNA action is mainly due to induction of mRNA degradation. This issue is discussed in more detail in Additional file1.
where r and p are the levels of host gene mRNA and protein products, s is the level of miRNAs, and g i represents the degradation rate of the molecular species i. As discussed in the introduction, intronic miRNAs (same-strand with the host) are expected to be co-transcribed with their host gene, therefore their production rate k r (q) has the same dependence on the input TF level.
A different representation was introduced in the context of bacterial sRNA regulation[45–47] and subsequently applied with slight modifications to eukaryotic miRNA regulation. In this representation, the degree of catalicity, i.e. the ability of a miRNA to affect multiple mRNAs without being degraded, was parametrized explicitly. The use of an effective phenomenological function (like the one in Equation 2) implicitly assumes a catalytic action, as commonly believed for miRNAs.
Moreover, miRNA biogenesis is a complex and highly regulated multi-step process (see[5, 18, 49] for a review) that finally leads to a mature miRNA loaded into the RNA-induced silencing complex (RISC) which is actually the active form that can downregulate the target mRNAs. Recent models of regulation by small RNAs have tried to build a more detailed modeling framework taking into account these processing steps[39, 50], for example including miRNA incorporation into the RISC complex as a reversible binding reaction. Such detailed modeling approaches have succesfully explained results from synthetic biology experiments. However, a more comprehensive mathematical representation of the biological processes involved comes at the expense of an increased number of free parameters, making more difficult an extensive exploration of the parameter space. While in our modeling strategy the dynamics of miRNA biogenesis is not included, in the supporting information (Additional file1) we report the analysis of the possible effects of the delay that the biogenesis process can introduce between miRNA transcription and miRNA-mediated repression. More generally, the relations between different possible models of miRNA regulation are discussed in detail in the supporting information (Additional file1), where it is shown that most of the results that will be presented in the following are essentially independent on the modeling strategy, provided that certain generic conditions on the parameters are satisfied.
In an analogous manner, it is possible to model the two circuits that we will use for comparison: a gene simply activated by the TF (sTF) without any feedback regulation (scheme in Figure1B) and a transcriptional self-loop (tSL), in which the negative feedback is realized through transcriptional repression (scheme in Figure1C). The properties of each circuit will be compared using a so called mathematically controlled comparison: all the common parameters will be kept to equivalent values, constraining the remainders so as to achieve the same steady state of protein concentration.
A deterministic description based on ordinary differential equations can effectively describe the mean kinetic behaviour of genetic circuits, thus its predictions can be tested with experiments based on averages over cell populations. In fact, equivalent mathematical treatments have correctly predicted the dynamical features of several endogenous and synthetics circuits[27, 30]. However, since gene expression is inherently a stochastic process[51–53], we will also make use of a stochastic description based on a master equation approach, that has Equations 3 as a “mean-field” limit (complete model in Additional file1). To compare the stochastic properties and the noise susceptibility of the three regulatory strategies in Figure1, we calculated analytically the relative fluctuations in protein level p at steady state and confirmed our results with Gillespie simulations (see the Methods section for details on simulations).
Autoregulation via intronic miRNAs has many of the structural properties of miRNA-mediated incoherent feedforward loops, that represent a diffused and functionally relevant motif in regulatory networks[25, 44, 54–56]. In fact, iMSLs can be thought as a mimimal feedforwad topology with perfect co-expression of the target gene and the miRNA buffering node, and thus, as previously observed, can be considered a special case of post-transcriptional incoherent feedforward loops. As a matter of fact, many of the functions that we will show the iMSL can perform are consistent with the feedforward nature of this circuit, and the analysis presented in this paper could be easily generalized to miRNA-mediated incoherent feedforward loops, adding new pieces to our understanding of microRNA regulation in simple circuits.
Response times to external signals are altered by autoregulation via intronic microRNAs
The response of a transcriptional unit to a stimulus, such as a change in a TF concentration, is steered by the lifetime of its mRNA and protein products. A fast protein turnover speeds up the kinetics, but with a consequent high metabolic cost, while in the case of long-living proteins the timescale of changes in concentration can be comparable to the cell cycle time[30, 57], which can be of several hours. However, the dynamics of a gene expression also depends strongly on the regulatory circuitry in which the gene is embedded. For example, it has been proven that negative transcriptional self regulation (like the one in Figure1C) and incoherent feedforward loops speed up the expression rise-time after induction[57, 58], while coherent feedforwad loops introduce delays.
We address in this section the question of how the host gene kinetics is changed by being a target of its intronic miRNA. To this aim, we consider two opposite simplified situations: (i) a sudden activating signal that fully saturates the promoter, and (ii) the opposite case of an istantaneous drop of the activating signal that completely switches off transcription. Case (i) can be studied assuming that at t=0 the transcription rate k r (q) switches from zero to its maximum value k r , and measuring the response time T ON defined as the time needed to reach half of the final protein steady state. In other words, we integrate numerically Equations 3 to calculate the time T ON such that p(T ON )/p ss = 0.5 (where p ss is the final steady-state protein level), starting from the condition r(0) = s(0) = p(0) = 0. In case (ii), in which we assume a drop of the activating signal at t = 0, we can similarly define a response time T OFF looking at the decrease of p(t) after a switch of the transcription rate from k r to zero at time t = 0. The same analysis is performed on a sTF (scheme in Figure1B) and a tSL (scheme in Figure1C) for comparison. The response time T0 of the simple transcription unit sTF is used as a normalization, since TON(OFF)/T0 is a measure of how much a circuit can alter the response time with respect to an unregulated gene.
As the repression increases, the response acceleration to an activating signal relies more and more on an overshoot of protein concentration, well above the final steady state, both for iMSLs and tSLs. If the input signal has to drive the host gene to its functional steady state, a large overshoot can be unwanted since it represents an unnecessary metabolic cost and a possible source of toxic effects. Thus, there is probably a limitation in the repression strength that can be applied to minimize the time separation between two functional steady states. On the other hand, a regime of strong repression makes the iMSL a pulse generator, a dynamical feature that can eventually lead to adaptation as will be discussed in a following section. The observation that the iMSL can function as a pulse generator is fully consistent with the fact that this circuit is a particular case of incoherent feedforwad loops that were indeed previously associated to pulse generation.
While the speeding up of activation is a property shared by incoherent feedforward loops (and specifically by iMSLs) and tSLs, an interesting peculiarity of the miRNA-mediated regulation in iMSLs emerges looking at the time required for p concentration to reach zero, starting from a constitutive level (Figure2D reports an example of this dynamics). The iMSL can delay the switch-off kinetics of the host in the same repression regime where it can accelerate the activation and the extent of the introduced delay is again dependent on the mRNA to miRNA lifetime ratio (Figure2B). This apparently counterintuitive behaviour can be easily qualitatively explained. When a constitutively expressed gene senses a transcription stop signal, the velocity of protein concentration decrease is established only by protein and mRNA degradation rates. For example, long living mRNAs are more persistent and can be translated for a longer time after the stop of transcription, and long living proteins are obviously more resilient. The same is true for tSLs or transcriptional feedforward loops: as the transcription is switched off, transcriptional repressors cannot exert any regulation and the protein level simply undergoes the exponential decrease dictated by mRNA and protein degradation. On the other hand, thanks to the post-transcriptional regulation in iMSLs, for each single miRNA that is degraded the still present mRNAs sense an increase in their translation rate. This increase clearly depends on the repression strength that miRNAs can exert (thus on the repression fold p ss /p0) and on the relative stability of mRNAs and miRNAs (τ r /τ s ), as a fast miRNA turnover leads to a higher translation rate of the remaining mRNAs. Eventually, the general increase in mRNA translation rate for each miRNA degradation event can lead to a temporary boost in protein concentration above the original steady state (see Figure2D).
It is important to notice that the dynamics just described can be altered if the miRNA acts mostly on mRNA degradation and depends on the timescale of miRNA-mRNA binding-unbinding. While the iMSL can always speed up the host gene expression in activation, the delay in the switch-off dynamics can vanish in case of fast miRNA-mediated induction of mRNA degradation. This issue is discussed in more details in the Additional file1.
The circuit response dynamics can robustly keep the host gene in a high-expression state
This property can be measured more quantitatively by the distance d from the ON-steady-state that is reached by the target protein level during a temporary absence of the input activator lasting a time T∗. As shown in Figure3C, the iMSL regulation keeps the host gene protein product close to its steady state in presence of input fluctuations that would almost switch-off a gene transcriptionally self-regulated or constitutively expressed.
Intronic microRNAs, targeting their host gene, can implement adaptation and Weber’s law
Adaptation is defined as the ability of a system to respond to a change in the input but subsequently return to the original level, even if the stimulus persists. Adaptation is ubiquitous in signaling systems. Examples of nearly perfect adaptation range from chemotaxis in bacteria to sensor cells in higher organisms. In all these systems, the benefit of adaptation can be summarized as the possibility of signal detection irrespective of the background level, thus widening the range of accessible signals and keeping the system robust in presence of perturbations.
Simple network topologies, as negative feedback loops with a buffering node or incoherent feedforward loops, can be at the basis of the cellular implementation of adaptation. As a special case of incoherent feedforward loops, also iMSLs are expected to be suitable to implement adaptation. In this section, we investigate whether and in what conditions the minimal topology of a post-transcriptional self-regulation through intronic miRNAs can perform adaptation.
It is easy to show analytically (see Additional file1) that in the regime of strong repression (s/h≫1 in the Michaelis-Menten function in Equation 2) the steady state of p concentration is independent of the input level q, which is clearly a hallmark of perfect adaptation: after an eventual dynamical response to a change in q, the system always returns to its original equilibrium level. On the other hand, it is impossible to achieve such an independence on the input level at equilibrium using a tSL, as confirmed by the fact that in general circuits with just two molecular species are not adaptive.
The general requirement of strong repression for an effective implementation of an adaptive response is in agreement with the results of an analogous theoretical analysis previously performed with detailed models of different feedforward loop circuits, including iMSLs. In particular, the prediction that stronger repression should lead to better adaptivity was tested experimentally with a synthetic transcriptional feedforward loop, using the number of DNA molecules coding for the circuit to modulate the input signal.
Certain adaptive systems, besides the ability to return to their original value after a signal response, present also a degree of response that is proportional to the relative change in the input signal and not to its absolute value. This feature is known as Weber’s law, originally introduced in the context of human sensory response. Recently, this dependence on input fold-change was demonstrated experimentally in eukaryotic signaling systems[74, 75], and theoretically the feedforward loop topology was proposed as a candidate to Weber’s law implementation in gene regulatory networks. Once again, it is natural to examine in what conditions also the minimal iMSL circuit can satisfy Weber’s law.
It is possible to show analytically (see Additional file1) that iMSLs are responsive to input fold-change if three conditions are satisfied:
Strong repression: s/h ≫ 1 ⇒ k p (s) ≈ k p h/s(condition for perfect adaptation),
Almost linear promoter activation k r (q) ≈ q k r /h r ,
Fast mRNA dynamics (short mRNA half-life with respect to miRNA and protein ones): r(t)→r ss .
Parameter space of adaptation and Weber’s law
Using the observables defined in Equations 4 and 5, it is possible to explore the conditions in which adaptation and Weber’s law are successfully performed by iMSLs. An illustrative example is depicted in Figure4C, where two effective parameters are varied: the effective promoter activation q/h r , and 1/h which measures the repression strength since h is the number of miRNAs necessary to reduce to one half the target translation rate. The grey region depicts the parameter space where precise adaptation is performed (P > 10), while in the excluded red region the dynamical response of the circuit is not able to go beyond the noise level (S < 2C V p ). The E value is reported with the color code in the legend when the minimal condition E < 0.1 is satisfied, i.e the two steps of the input produce the same response within 10%. Adaptation and Weber’s law can be encoded by iMSLs in a parameter region that span several orders of magnitude of the effective parameters. Therefore, the only constraint is that the effective parameters have to approach the appropriate limits, without the need of fine-tuning.
It should be noticed that the general condition of strong repression required for both functions is limited by the circuit sensitivity. This is partially due to the fact that a too strong repression can rise the noise level of the circuit (see next section) making the achievement of a signal significantly above fluctuations harder.
It is interesting to consider what are the functional advantages that these two functions can provide to the host gene. Both adaptation and Weber’s law can bestow robustness to the expression program of the host gene. An expression state that is not influenced by constant inputs thanks to adaptation is robust with respect to the ubiquitous cell-to-cell variability in TF concentrations, but it is still responsive to signals that induce dynamical variations of TF levels. When additionally Weber’s law is implemented, also the dynamical response can be kept homogeneous in a cell population. In fact, in this case the response profile is only due to the input fold-change and not on its absolute value that is affected by the potentially variable background level. Moreover, Weber’s law naturally encodes a noise filter. In fact, since the noise level is expected to scale with the background TF concentration, a dependence on fold-change can rescale appropriately the threshold at which the response is triggered, thus allowing a better signal/noise discrimination in different background conditions.
Autoregulation via intronic microRNAs reduces the host gene expression fluctuations
It is interesting to notice that, even if iMSLs and tSLs show similar noise reduction properties, the miRNA-mediated self-regulation actually performs better than the transcriptional self-regulation. As it is possible to see in Figure5A (where histograms and continuous lines are respectively the result of Gillespie simulations with full nonlinear dynamics and gamma distributions with analytically calculated moments), the probability distributions of the target protein level for the three circuits are different. Both autoregulatory circuits lead to a target distribution less sparse than a sTF, showing that they effectively reduce fluctuations, but the iMSL distribution is clearly more peaked than the tSL one. Similarly, both self-regulation strategies show an optimal noise buffering for an intermediate repression strength, but again the attenuation is larger in the miRNA-mediated case (see Figure5B). This is more clearly shown in Figure5C,D, where the noise reduction (withrepresenting the target noise in the case of a simple transcription unit producing the same mean amount of proteins) is reported for the two autoregulatory circuits. Noise reduction is explored for different levels of transcriptional activation (〈q〉/h r ) and target repression (〈p〉/〈p0〉, where 〈p0〉 is the target mean value in absence of repression) to shed light into noise control and target suppression interdependence. Averages are here intended at steady state, thus the repression level measure 〈p〉/〈p0〉 for the stochastic model is perfectly equivalent to the one used in the response time analysis (see Figure2). In the regime where the target is more sensitive to TF fluctuations, i.e. q is far from saturating the promoter, the iMSL can reduce the fluctuations up to a factor 0.5 (Figure5C), while the tSL (Figure5D) is much less effective. Moreover, the heat maps in Figures5C,D indicate that the iMSL can buffer fluctuations over a wider range of conditions as well as to a greater extent.
As pointed out in, an optimal miRNA-mediated noise buffering does not necessarily require a strong repression. Indeed, Figure5C shows that a reduction of the mean protein expression to 50% of its constitutive level is sufficient to reduce the noise by approximately 40%. This means that the intronic miRNA can keep the expression of its host gene in its homeostatic regime, while filtering out fluctuations, without exerting a strong reduction of its concentration. This result agrees well with the observation that miRNAs act often to fine-tune their targets rather than to switch them off completely.
Sketch of the one-to-many topology-function map
This section summarizes the functions found to be associated to intronic miRNA-mediated self-loops into a qualitative “map of functions”, showing the different, although overlapping, ranges of biochemical parameters in which each specific function is optimized. The emerging map between parameter values and functions can be useful to understand the presence of the iMSL architecture in different biological contexts and gives general guidelines for the design of synthetic circuits with a desired behaviour, well beyond the simple suggestion of a network topology.
Looking at a finer scale the strong repression regime (〈p〉/〈p0〉 < 1/2), a smooth transition in the dynamical behaviour of the circuit can be observed (see Figure6B). At first, the host gene is able to fastly transit between two well distinct steady states after induction. When the repression is further increased, this fast ON-activation relies increasingly on an overshoot well above the final equilibrium at which the dynamics asympotically relaxes. Therefore, the concentration profile resembles a pulse. Finally, for high enough repression the system returns to the initial steady state after the pulse, a necessary condition for the implementation of adaptation and Weber’s law.
The relative half-life of the molecules involved, in particular of miRNAs and mRNAs, is another ingredient that can strongly influence the dynamical behaviour (see Figure6C). For example, a miRNA half-life comparable to the mRNA one allows a trade-off between acceleration of the ON-dynamics and delay of the OFF-dynamics, making the state of high expression of the host gene robust to fluctuations. On the other hand, mRNA lifetime must be short with respect to the other molecules lifetimes for a dynamical response following Weber’s law.
The present analysis of the iMSL functions considers the circuit as isolated, while realistically a single microRNA can target hundreds of genes. As recently pointed out, the degree of repression of a target depends on the level of expression of all possible target genes[45, 81], since their mRNAs can dilute the pool of available miRNAs. Therefore, the expression profile of alternative miRNA targets is a variable that can potentially alter the dynamics of iMSLs (as shown for incoherent feedforward loops), and thus have to be carefully taken into account in experimental tests on endogenous iMSLs.
Identification of intronic miRNA-mediated self-loops in human
This study presents a fairly comprehensive survey of the possibile functions associated to a miRNA-mediated circuit composed of one protein-coding gene (host gene) negatively regulated by a miRNA located in one of its introns. In particular, we have shown that, thanks to the miRNA-mediated self-regulation, the host gene expression responds to input changes with an altered timing, its response can be adaptative and follow Weber’s law, and fluctuations propagating from the upstream network are buffered at steady state. Each of these functions can confer robustness to the expression program of the host gene, suggesting that miRNA-mediated self-loops represents a simple homeostatic control. For example, adaptation makes the host gene expression level at equilibrium independent of the cell-to-cell variability of transcription factor expression without compromising its sensitivity to input changes. Similarly, the host expression dynamics, as modified by miRNA autorepression, can mantain the host gene in a high-expression state in the face of downwards fluctuations in activators’ concentration. The association of miRNA-mediated self-loops with functions with different specificities, but apparently the same final aim, suggests that, depending on the desired level of the host gene expression and on the type of fluctuations that have to be more frequently filtered out, the details of the regulatory interactions and the characteristics of the molecules involved could have been fine-tuned over evolutionary timescales accordingly. Such a fine-tuning of expression parameters has been shown to be possible even over short timescales in in vitro evolutionary experiments.
The comparison with an unregulated transcriptional unit and with a transcriptional negative feedback indicates that the specificities of miRNA regulation makes the post-transcriptional circuit better suited to implement a homeostatic control. This result is in line with the accumulating clues that miRNAs can help the cell to function reliably in presence of perturbations[31–34].
Finally, our systematic analysis of the constraints on biochemical parameters necessary to optimize each function can guide the realization of synthetic versions of miRNA-mediated self-loops, as well as contribute to the understanding of the role of their many occurrences in endogenous networks. In this perspective, we also provide a list of bioinformatically predicted miRNA-mediated self-loops in human for future experimental tests.
Simulations were implemented by using Gillespie’s first reaction algorithm. The reactions simulated are those presented in Figure1 with additional transcription, translation and degradations for the input transcription factor q. Reactions that depend on a regulator were allowed to have as rates the corresponding full nonlinear functions. Results in Figure5 are at steady state, which is assumed to be reached when the deterministic evolution of the system in analysis is at a distance from the steady state (its asymptotic value) smaller than its 0.05% (more than 10 times the protein half-life). Each data point or histogram is the result of 100000 trials.
In order to identify human intragenic miRNAs, and associate them to their host genes, we used Ensembl-release 57 database. We collected the data for all human known protein coding genes (consisting in a total of 22.257 entries with a stable Ensembl Gene Identifier (ENSG)). For each gene we then retained for further analysis only the longest Transcript Identifier (ENST). The data on human miRNAs were extracted from Ensembl v.57, that includes miRBase v.13 (Table S1). To identify the iMSLs, we used eight tools for miRNA/target gene interaction predictions: TargetScan human v. 5.0[3, 82], miRanda - release 2008[83, 84], RNA22, PITA-4way, MirTarget2, PicTar, Diana microT v.3 and TargetMiner v.1. To test the over-representation of the putative iMSLs, we performed two different types of randomization strategies. Specifically, we randomly permuted 1000 times the intronic miRNA/host-gene link and the union of miRNA/target gene interactions datasets predicted by the different algorithms. In both cases we created, according to the two reshuffling strategies, 1000 independent reshuffled copies of the original network. Then for each of them we evaluated the number of iMSLs confirmed by at least one algorithm and obtained in this way the two histograms plotted in Figure7C. Then for both reshuffling strategies a Z score can be defined as :, where x is the actual number of iMSLs in the network (i.e. 77), while <x> r and σ r are the mean value and the standard deviation of the distribution of iMSLs in the reshuffled samples. These Z scores turned out to be rather large: Z = 4.63 for the first test and Z = 5.52 for the second one (see Figure7C). These values show that the number of iMSLs in the network is definitely larger than random and thus suggest, in agreement with what already observed in[20, 26], that this particular class of network motifs is most probably under positive selection.
CB, MO and MC designed research. CB, MO, ME, DC performed research. CB, MO, ME and MC wrote the paper. All authors read and approved the final manuscript.
We would like to thank Marco Cosentino Lagomarsino for useful comments on the manuscript. MO acknowledges support from the Human Frontier Science Program Organization (Grant RGY0069/2009-C). DC acknowledges support by Fondazione Piemontese per la Ricerca sul Cancro-ONLUS and Associazione Italiana per la Ricerca sul Cancro (Intramural Grant 5x1000 2008).
- Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120: 15-20. 10.1016/j.cell.2004.12.035. [http://dx.doi.org/10.1016/j.cell.2004.12.035], 10.1016/j.cell.2004.12.035View ArticleGoogle Scholar
- Flynt AS, Lai EC: Biological principles of microRNA-mediated regulation: shared themes amid diversity. Nat Rev Genet. 2008, 9 (11): 831-842. 10.1038/nrg2455. [http://dx.doi.org/10.1038/nrg2455], 10.1038/nrg2455View ArticleGoogle Scholar
- Friedman RC, Farh KKH, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19: 92-105. [http://dx.doi.org/10.1101/gr.082701.108],View ArticleGoogle Scholar
- Ambros V: The functions of animal microRNAs. Nature. 2004, 431 (7006): 350-355. 10.1038/nature02871. [http://dx.doi.org/10.1038/nature02871], 10.1038/nature02871View ArticleGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.View ArticleGoogle Scholar
- Bushati N, Cohen SM: microRNA functions. Annu Rev Cell Dev Biol. 2007, 23: 175-205. 10.1146/annurev.cellbio.23.090506.123406. [http://dx.doi.org/10.1146/annurev.cellbio.23.090506.123406], 10.1146/annurev.cellbio.23.090506.123406View ArticleGoogle Scholar
- Stefani G, Slack FJ: Small non-coding RNAs in animal development. Nat Rev Mol Cell Biol. 2008, 9 (3): 219-230. 10.1038/nrm2347. [http://dx.doi.org/10.1038/nrm2347], 10.1038/nrm2347View ArticleGoogle Scholar
- Alvarez-Garcia I, Miska EA: MicroRNA functions in animal development and human disease. Development. 2005, 132 (21): 4653-4662. 10.1242/dev.02073. [http://dx.doi.org/10.1242/dev.02073], 10.1242/dev.02073View ArticleGoogle Scholar
- Esquela-Kerscher A, Slack FJ: Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006, 6 (4): 259-269. 10.1038/nrc1840. [http://dx.doi.org/10.1038/nrc1840], 10.1038/nrc1840View ArticleGoogle Scholar
- Hinske LCG, Galante PAF, Kuo WP, Ohno-Machado L: A potential role for intragenic miRNAs on their hosts’ interactome. BMC Genomics. 2010, 11: 533-10.1186/1471-2164-11-533. [http://dx.doi.org/10.1186/1471-2164-11-533], 10.1186/1471-2164-11-533View ArticleGoogle Scholar
- Li SC, Tang P, Lin WC: Intronic microRNA: discovery and biological implications. DNA Cell Biol. 2007, 26 (4): 195-207. 10.1089/dna.2006.0558. [http://dx.doi.org/10.1089/dna.2006.0558], 10.1089/dna.2006.0558View ArticleGoogle Scholar
- Saini HK, Griffiths-Jones S, Enright AJ: Genomic analysis of human microRNA transcripts. Proc Natl Acad Sci U S A. 2007, 104 (45): 17719-17724. 10.1073/pnas.0703890104. [http://dx.doi.org/10.1073/pnas.0703890104], 10.1073/pnas.0703890104View ArticleGoogle Scholar
- Rodriguez A, Griffiths-Jones S, Ashurst JL, Bradley A: Identification of mammalian microRNA host genes and transcription units. Genome Res. 2004, 14 (10A): 1902-1910. 10.1101/gr.2722704. [http://dx.doi.org/10.1101/gr.2722704], 10.1101/gr.2722704View ArticleGoogle Scholar
- Baskerville S, Bartel DP: Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. RNA. 2005, 11 (3): 241-247. 10.1261/rna.7240905. [http://dx.doi.org/10.1261/rna.7240905], 10.1261/rna.7240905View ArticleGoogle Scholar
- Lin SL, Miller JD, Ying SY: Intronic microRNA (miRNA). J Biomed Biotechnol. 2006, 2006 (4): 26818-[http://dx.doi.org/10.1155/JBB/2006/26818],Google Scholar
- Kim YK, Kim VN: Processing of intronic microRNAs. EMBO J. 2007, 26 (3): 775-783. 10.1038/sj.emboj.7601512. [http://dx.doi.org/10.1038/sj.emboj.7601512], 10.1038/sj.emboj.7601512View ArticleGoogle Scholar
- Ma N, Wang X, Qiao Y, Li F, Hui Y, Zou C, Jin J, Lv G, Peng Y, Wang L, Huang H, Zhou L, Zheng X, Gao X: Coexpression of an intronic microRNA and its host gene reveals a potential role for miR-483-5p as an IGF2 partner. Mol Cell Endocrinol. 2011, 333: 96-101. 10.1016/j.mce.2010.11.027. [http://dx.doi.org/10.1016/j.mce.2010.11.027], 10.1016/j.mce.2010.11.027View ArticleGoogle Scholar
- Kim VN, Han J, Siomi MC: Biogenesis of small RNAs in animals. Nat Rev Mol Cell Biol. 2009, 10 (2): 126-139. 10.1038/nrm2632. [http://dx.doi.org/10.1038/nrm2632], 10.1038/nrm2632View ArticleGoogle Scholar
- Morlando M, Ballarino M, Gromak N, Pagano F, Bozzoni I, Proudfoot NJ: Primary microRNA transcripts are processed co-transcriptionally. Nat Struct Mol Biol. 2008, 15 (9): 902-909. 10.1038/nsmb.1475.View ArticleGoogle Scholar
- Ozsolak F, Poling LL, Wang Z, Liu H, Liu XS, Roeder RG, Zhang X, Song JS, Fisher DE: Chromatin structure analyses identify miRNA promoters. Genes Dev. 2008, 22 (22): 3172-3183. 10.1101/gad.1706508. [http://dx.doi.org/10.1101/gad.1706508], 10.1101/gad.1706508View ArticleGoogle Scholar
- Monteys AM, Spengler RM, Wan J, Tecedor L, Lennox KA, Xing Y, Davidson BL: Structure and activity of putative intronic miRNA promoters. RNA. 2010, 16 (3): 495-505. 10.1261/rna.1731910. [http://dx.doi.org/10.1261/rna.1731910], 10.1261/rna.1731910View ArticleGoogle Scholar
- He C, Li Z, Chen P, Huang H, Hurst LD, Chen J: Young intragenic miRNAs are less coexpressed with host genes than old ones: implications of miRNA-host gene coevolution. Nucleic Acids Res. 2012, 40 (9): 4002-4012. 10.1093/nar/gkr1312. [http://dx.doi.org/10.1093/nar/gkr1312], 10.1093/nar/gkr1312View ArticleGoogle Scholar
- Barik S: An intronic microRNA silences genes that are functionally antagonistic to its host gene. Nucleic Acids Res. 2008, 36 (16): 5232-5241. 10.1093/nar/gkn513. [http://dx.doi.org/10.1093/nar/gkn513], 10.1093/nar/gkn513View ArticleGoogle Scholar
- Lutter D, Marr C, Krumsiek J, Lang EW: Intronic microRNAs support their host genes by mediating synergistic and antagonistic regulatory effects. BMC Genomics. 2010, 11: 224-10.1186/1471-2164-11-224. [http://dx.doi.org/10.1186/1471-2164-11-224], 10.1186/1471-2164-11-224View ArticleGoogle Scholar
- Tsang J, Zhu J, van Oudenaarden A: MicroRNA-mediated feedback and feedforward loops are recurrent network motifs in mammals. Mol Cell. 2007, 26 (5): 753-767. 10.1016/j.molcel.2007.05.018. [http://dx.doi.org/10.1016/j.molcel.2007.05.018], 10.1016/j.molcel.2007.05.018View ArticleGoogle Scholar
- Megraw M, Sethupathy P, Gumireddy K, Jensen S, Huang Q: Isoform specific gene auto-regulation via miRNAs: a case study on miR-128b and ARPP-21. Theor Chem Acc. 2009, 125: 593-598.View ArticleGoogle Scholar
- Alon U: Network motifs: theory and experimental approaches. Nat Rev Genet. 2007, 8 (6): 450-461. 10.1038/nrg2102. [http://dx.doi.org/10.1038/nrg2102], 10.1038/nrg2102View ArticleGoogle Scholar
- Sun Y, Bai Y, Zhang F, Wang Y, Guo Y: miR-126 inhibits non-small cell lung cancer cells proliferation by targeting EGFL7. Biochem Biophys Res Commun. 2010, 391 (3): 1483-1489. 10.1016/j.bbrc.2009.12.098. [http://dx.doi.org/10.1016/j.bbrc.2009.12.098], 10.1016/j.bbrc.2009.12.098View ArticleGoogle Scholar
- Nikolic I, Plate KH, Schmidt MH: EGFL7 meets miRNA-126: an angiogenesis alliance. J Angiogenes Res. 2010, 2: 9-10.1186/2040-2384-2-9. [http://dx.doi.org/10.1186/2040-2384-2-9], 10.1186/2040-2384-2-9View ArticleGoogle Scholar
- Alon U: An Introduction to Systems Biology: Design Principles of Biological Circuits. 2007, Boca Raton; FL: Chapman & Hall/CRC Mathematical/Taylor & FrancisGoogle Scholar
- Inui M, Martello G, Piccolo S: MicroRNA control of signal transduction. Nat Rev Mol Cell Biol. 2010, 11 (4): 252-263. [http://dx.doi.org/10.1038/nrm2868],View ArticleGoogle Scholar
- Li X, Cassidy JJ, Reinke CA, Fischboeck S, Carthew RW: A microRNA imparts robustness against environmental fluctuation during development. Cell. 2009, 137 (2): 273-282. 10.1016/j.cell.2009.01.058. [http://dx.doi.org/10.1016/j.cell.2009.01.058], 10.1016/j.cell.2009.01.058View ArticleGoogle Scholar
- Li Y, Wang F, Lee JA: MicroRNA-9a ensures the precise specification of sensory organ precursors in Drosophila. Genes Dev. 2006, 20 (20): 2793-2805. 10.1101/gad.1466306. [http://dx.doi.org/10.1101/gad.1466306], 10.1101/gad.1466306View ArticleGoogle Scholar
- Ebert MS, Sharp PA: Roles for MicroRNAs in Conferring Robustness to Biological Processes. Cell. 2012, 149 (3): 515-524. 10.1016/j.cell.2012.04.005. [http://dx.doi.org/10.1016/j.cell.2012.04.005], 10.1016/j.cell.2012.04.005View ArticleGoogle Scholar
- Haynes KA, Silver PA: Eukaryotic systems broaden the scope of synthetic biology. J Cell Biol. 2009, 187 (5): 589-596. 10.1083/jcb.200908138. [http://dx.doi.org/10.1083/jcb.200908138], 10.1083/jcb.200908138View ArticleGoogle Scholar
- Khalil AS, Collins JJ: Synthetic biology: applications come of age. Nat Rev Genet. 2010, 11 (5): 367-379. [http://dx.doi.org/10.1038/nrg2775],View ArticleGoogle Scholar
- Deans TL, Cantor CR, Collins JJ: A tunable genetic switch based on RNAi and repressor proteins for regulating gene expression in mammalian cells. Cell. 2007, 130 (2): 363-372. 10.1016/j.cell.2007.05.045. [http://dx.doi.org/10.1016/j.cell.2007.05.045], 10.1016/j.cell.2007.05.045View ArticleGoogle Scholar
- Rinaudo K, Bleris L, Maddamsetti R, Subramanian S, Weiss R, Benenson Y: A universal RNAi-based logic evaluator that operates in mammalian cells. Nat Biotechnol. 2007, 25 (7): 795-801. 10.1038/nbt1307. [http://dx.doi.org/10.1038/nbt1307], 10.1038/nbt1307View ArticleGoogle Scholar
- Bleris L, Xie Z, Glass D, Adadey A, Sontag E, Benenson Y: Synthetic incoherent feedforward circuits show adaptation to the amount of their genetic template. Mol Syst Biol. 2011, 7: 519-[http://dx.doi.org/10.1038/msb.2011.49],View ArticleGoogle Scholar
- Bintu L, Buchler NE, Garcia HG, Gerland U, Hwa T, Kondev J, Phillips R: Transcriptional regulation by the numbers: models. Curr Opin Genet Dev. 2005, 15 (2): 116-124. 10.1016/j.gde.2005.02.007. [http://dx.doi.org/10.1016/j.gde.2005.02.007], 10.1016/j.gde.2005.02.007View ArticleGoogle Scholar
- Valencia-Sanchez MA, Liu J, Hannon GJ, Parker R: Control of translation and mRNA degradation by miRNAs and siRNAs. Genes Dev. 2006, 20 (5): 515-524. 10.1101/gad.1399806. [http://dx.doi.org/10.1101/gad.1399806], 10.1101/gad.1399806View ArticleGoogle Scholar
- Cuccato G, Polynikis A, Siciliano V, Graziano M, di Bernardo M, di Bernardo D: Modeling RNA interference in mammalian cells. BMC Syst Biol. 2011, 5: 19-10.1186/1752-0509-5-19. [http://dx.doi.org/10.1186/1752-0509-5-19], 10.1186/1752-0509-5-19View ArticleGoogle Scholar
- Komorowski M, Miekisz J, Kierzek AM: Translational repression contributes greater noise to gene expression than transcriptional repression. Biophys J. 2009, 96 (2): 372-384. 10.1016/j.bpj.2008.09.052. [http://dx.doi.org/10.1016/j.bpj.2008.09.052], 10.1016/j.bpj.2008.09.052View ArticleGoogle Scholar
- Osella M, Bosia C, Corà D, Caselle M: The role of incoherent microRNA-mediated feedforward loops in noise buffering. PLoS Comput Biol. 2011, 7 (3): e1001101-10.1371/journal.pcbi.1001101. [http://dx.doi.org/10.1371/journal.pcbi.1001101], 10.1371/journal.pcbi.1001101View ArticleGoogle Scholar
- Levine E, Zhang Z, Kuhlman T, Hwa T: Quantitative characteristics of gene regulation by small RNA. PLoS Biol. 2007, 5 (9): e229-10.1371/journal.pbio.0050229. [http://dx.doi.org/10.1371/journal.pbio.0050229], 10.1371/journal.pbio.0050229View ArticleGoogle Scholar
- Shimoni Y, Friedlander G, Hetzroni G, Niv G, Altuvia S, Biham O, Margalit H: Regulation of gene expression by small non-coding RNAs: a quantitative view. Mol Syst Biol. 2007, 3: 138-[http://dx.doi.org/10.1038/msb4100181],View ArticleGoogle Scholar
- Mehta P, Goyal S, Wingreen NS: A quantitative comparison of sRNA-based and protein-based gene regulation. Mol Syst Biol. 2008, 4: 221-[http://dx.doi.org/10.1038/msb.2008.58],View ArticleGoogle Scholar
- Mukherji S, Ebert MS, Zheng GXY, Tsang JS, Sharp PA, van Oudenaarden A: MicroRNAs can generate thresholds in target gene expression. Nat Genet. 2011, 43 (9): 854-859. 10.1038/ng.905. [http://dx.doi.org/10.1038/ng.905], 10.1038/ng.905View ArticleGoogle Scholar
- Krol J, Loedige I, Filipowicz W: The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet. 2010, 11 (9): 597-610. [http://dx.doi.org/10.1038/nrg2843],Google Scholar
- Tigges M, Dénervaud N, Greber D, Stelling J, Fussenegger M: A synthetic low-frequency mammalian oscillator. Nucleic Acids Res. 2010, 38 (8): 2702-2711. 10.1093/nar/gkq121. [http://dx.doi.org/10.1093/nar/gkq121], 10.1093/nar/gkq121View ArticleGoogle Scholar
- Kaern M, Elston TC, Blake WJ, Collins JJ: Stochasticity in gene expression: from theories to phenotypes. Nat Rev Genet. 2005, 6 (6): 451-464. 10.1038/nrg1615. [http://dx.doi.org/10.1038/nrg1615], 10.1038/nrg1615View ArticleGoogle Scholar
- Maheshri N, O’Shea EK: Living with noisy genes: how cells function reliably with inherent variability in gene expression. Annu Rev Biophys Biomol Struct. 2007, 36: 413-434. 10.1146/annurev.biophys.36.040306.132705. [http://dx.doi.org/10.1146/annurev.biophys.36.040306.132705], 10.1146/annurev.biophys.36.040306.132705View ArticleGoogle Scholar
- Raj A, van Oudenaarden A: Cell. 2008, 135 (2): 216-226. 10.1016/j.cell.2008.09.050. [http://dx.doi.org/10.1016/j.cell.2008.09.050], 10.1016/j.cell.2008.09.050View ArticleGoogle Scholar
- Re A, Corà D, Taverna D, Caselle M: Genome-wide survey of microRNA-transcription factor feed-forward regulatory circuits in human. Mol Biosyst. 2009, 5 (8): 854-867. 10.1039/b900177h. [http://dx.doi.org/10.1039/b900177h], 10.1039/b900177hView ArticleGoogle Scholar
- El Baroudi M, Corà D, Bosia C, Osella M, Caselle M: A curated database of miRNA mediated feed-forward loops involving MYC as master regulator. PLoS One. 2011, 6 (3): e14742-10.1371/journal.pone.0014742. [http://dx.doi.org/10.1371/journal.pone.0014742], 10.1371/journal.pone.0014742View ArticleGoogle Scholar
- Eduati F, Camillo BD, Karbiener M, Scheideler M, Corà D, Caselle M, Toffolo G: Dynamic modeling of miRNA-mediated feed-forward loops. J Comput Biol. 2012, 19 (2): 188-199. 10.1089/cmb.2011.0274. [http://dx.doi.org/10.1089/cmb.2011.0274], 10.1089/cmb.2011.0274View ArticleGoogle Scholar
- Rosenfeld N, Elowitz MB, Alon U: Negative autoregulation speeds the response times of transcription networks. J Mol Biol. 2002, 323 (5): 785-793. 10.1016/S0022-2836(02)00994-4.View ArticleGoogle Scholar
- Mangan S, Itzkovitz S, Zaslaver A, Alon U: The incoherent feed-forward loop accelerates the response-time of the gal system of Escherichia coli. J Mol Biol. 2006, 356 (5): 1073-1081. 10.1016/j.jmb.2005.12.003. [http://dx.doi.org/10.1016/j.jmb.2005.12.003], 10.1016/j.jmb.2005.12.003View ArticleGoogle Scholar
- Mangan S, Zaslaver A, Alon U: The coherent feedforward loop serves as a sign-sensitive delay element in transcription networks. J Mol Biol. 2003, 334 (2): 197-204. 10.1016/j.jmb.2003.09.049.View ArticleGoogle Scholar
- Fan J, Yang X, Wang W, Wood WH, Becker KG, Gorospe M: Global analysis of stress-regulated mRNA turnover by using cDNA arrays. Proc Natl Acad Sci U S A. 2002, 99 (16): 10611-10616. 10.1073/pnas.162212399. [http://dx.doi.org/10.1073/pnas.162212399], 10.1073/pnas.162212399View ArticleGoogle Scholar
- Chen CYA, Ezzeddine N, Shyu AB: Messenger RNA half-life measurements in mammalian cells. Methods Enzymol. 2008, 448: 335-357. [http://dx.doi.org/10.1016/S0076-6879(08)02617-7],View ArticleGoogle Scholar
- Yang E, van Nimwegen E, Zavolan M, Rajewsky N, Schroeder M, Magnasco M, Darnell JE: Decay rates of human mRNAs: correlation with functional characteristics and sequence attributes. Genome Res. 2003, 13 (8): 1863-1872. [http://dx.doi.org/10.1101/gr.1272403],Google Scholar
- Sharova LV, Sharov AA, Nedorezov T, Piao Y, Shaik N, Ko MSH: Database for mRNA half-life of 19 977 genes obtained by DNA microarray analysis of pluripotent and differentiating mouse embryonic stem cells. DNA Res. 2009, 16: 45-58. 10.1093/dnares/dsn030. [http://dx.doi.org/10.1093/dnares/dsn030], 10.1093/dnares/dsn030View ArticleGoogle Scholar
- Yen HCS, Xu Q, Chou DM, Zhao Z, Elledge SJ: Global protein stability profiling in mammalian cells. Science. 2008, 322 (5903): 918-923. 10.1126/science.1160489. [http://dx.doi.org/10.1126/science.1160489], 10.1126/science.1160489View ArticleGoogle Scholar
- van Rooij E, Sutherland LB, Qi X, Richardson JA, Hill J, Olson EN: Control of stress-dependent cardiac growth and gene expression by a microRNA. Science. 2007, 316 (5824): 575-579. 10.1126/science.1139089. [http://dx.doi.org/10.1126/science.1139089], 10.1126/science.1139089View ArticleGoogle Scholar
- Kai ZS, Pasquinelli AE: MicroRNA assassins: factors that regulate the disappearance of miRNAs. Nat Struct Mol Biol. 2010, 17: 5-10. 10.1038/nsmb.1762. [http://dx.doi.org/10.1038/nsmb.1762], 10.1038/nsmb.1762View ArticleGoogle Scholar
- Sethi P, Lukiw WJ: Micro-RNA abundance and stability in human brain: specific alterations in Alzheimer’s disease temporal lobe neocortex. Neurosci Lett. 2009, 459 (2): 100-104. 10.1016/j.neulet.2009.04.052. [http://dx.doi.org/10.1016/j.neulet.2009.04.052], 10.1016/j.neulet.2009.04.052View ArticleGoogle Scholar
- Chatterjee S, Grosshans H: Active turnover modulates mature microRNA activity in Caenorhabditis elegans. Nature. 2009, 461 (7263): 546-549. 10.1038/nature08349. [http://dx.doi.org/10.1038/nature08349], 10.1038/nature08349View ArticleGoogle Scholar
- Kollmann M, Lovdok L, Bartholomé K, Timmer J, Sourjik V: Design principles of a bacterial signalling network. Nature. 2005, 438 (7067): 504-507. 10.1038/nature04228. [http://dx.doi.org/10.1038/nature04228], 10.1038/nature04228View ArticleGoogle Scholar
- Matthews HR, Reisert J: Calcium, the two-faced messenger of olfactory transduction and adaptation. Curr Opin Neurobiol. 2003, 13 (4): 469-475. 10.1016/S0959-4388(03)00097-7.View ArticleGoogle Scholar
- Ma W, Trusina A, El-Samad H, Lim WA, Tang C: Defining network topologies that can achieve biochemical adaptation. Cell. 2009, 138 (4): 760-773. 10.1016/j.cell.2009.06.013. [http://dx.doi.org/10.1016/j.cell.2009.06.013], 10.1016/j.cell.2009.06.013View ArticleGoogle Scholar
- Sontag ED: Remarks on feedforward circuits, adaptation, and pulse memory. IET Syst Biol. 2010, 4: 39-51. 10.1049/iet-syb.2008.0171. [http://dx.doi.org/10.1049/iet-syb.2008.0171], 10.1049/iet-syb.2008.0171View ArticleGoogle Scholar
- François P, Siggia ED: A case study of evolutionary computation of biochemical adaptation. Phys Biol. 2008, 5 (2): 026009-10.1088/1478-3975/5/2/026009. [http://dx.doi.org/10.1088/1478-3975/5/2/026009], 10.1088/1478-3975/5/2/026009View ArticleGoogle Scholar
- Goentoro L, Kirschner MW: Evidence that fold-change, and not absolute level, of beta-catenin dictates Wnt signaling. Mol Cell. 2009, 36 (5): 872-884. 10.1016/j.molcel.2009.11.017. [http://dx.doi.org/10.1016/j.molcel.2009.11.017], 10.1016/j.molcel.2009.11.017View ArticleGoogle Scholar
- Cohen-Saidon C, Cohen AA, Sigal A, Liron Y, Alon U: Dynamics and variability of ERK2 response to EGF in individual living cells. Mol Cell. 2009, 36 (5): 885-893. 10.1016/j.molcel.2009.11.025. [http://dx.doi.org/10.1016/j.molcel.2009.11.025], 10.1016/j.molcel.2009.11.025View ArticleGoogle Scholar
- Goentoro L, Shoval O, Kirschner MW, Alon U: The incoherent feedforward loop can provide fold-change detection in gene regulation. Mol Cell. 2009, 36 (5): 894-899. 10.1016/j.molcel.2009.11.018. [http://dx.doi.org/10.1016/j.molcel.2009.11.018], 10.1016/j.molcel.2009.11.018View ArticleGoogle Scholar
- Dublanche Y, Michalodimitralis K, Kummerer N, Foglierini M, Serrano L: Noise in transcription negative feedback loops: simulation and experimental analysis. Molecular Systems Biology. 2006, 2: 41-View ArticleGoogle Scholar
- Shahrezaei V, Ollivier J, PS S: Colored extrinsic fluctuations and stochastic gene expression. Mol Syst Biol. 2008, 4: 196-View ArticleGoogle Scholar
- Singh A, Hespanha JP: Optimal feedback strength for noise suppression in autoregulatory gene networks. Biophys J. 2009, 96 (10): 4013-4023. 10.1016/j.bpj.2009.02.064. [http://dx.doi.org/10.1016/j.bpj.2009.02.064], 10.1016/j.bpj.2009.02.064View ArticleGoogle Scholar
- Baek D, Villén J, Shin C, Camargo FD, Gygi SP, Bartel DP: The impact of microRNAs on protein output. Nature. 2008, 455 (7209): 64-71. 10.1038/nature07242. [http://dx.doi.org/10.1038/nature07242], 10.1038/nature07242View ArticleGoogle Scholar
- Arvey A, Larsson E, Sander C, Leslie CS, Marks DS: Target mRNA abundance dilutes microRNA and siRNA activity. Mol Syst Biol. 2010, 6: 363-[http://dx.doi.org/10.1038/msb.2010.24],View ArticleGoogle Scholar
- Lewis BP, hung Shih I, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115 (7): 787-798. 10.1016/S0092-8674(03)01018-3.View ArticleGoogle Scholar
- John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets. PLoS Biol. 2004, 2 (11): e363-10.1371/journal.pbio.0020363. [http://dx.doi.org/10.1371/journal.pbio.0020363], 10.1371/journal.pbio.0020363View ArticleGoogle Scholar
- Betel D, Wilson M, Gabow A, Marks DS, Sander C: The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008, 36 (Database issue): D149-D153. [http://dx.doi.org/10.1093/nar/gkm995],Google Scholar
- Miranda KC, Huynh T, Tay Y, Ang YS, Tam WL, Thomson AM, Lim B, Rigoutsos I: A pattern-based method for the identification of MicroRNA binding sites and their corresponding heteroduplexes. Cell. 2006, 126 (6): 1203-1217. 10.1016/j.cell.2006.07.031. [http://dx.doi.org/10.1016/j.cell.2006.07.031], 10.1016/j.cell.2006.07.031View ArticleGoogle Scholar
- Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition. Nat Genet. 2007, 39 (10): 1278-1284. 10.1038/ng2135. [http://dx.doi.org/10.1038/ng2135], 10.1038/ng2135View ArticleGoogle Scholar
- Wang X, Naqa IME: Prediction of both conserved and nonconserved microRNA targets in animals. Bioinformatics. 2008, 24 (3): 325-332. 10.1093/bioinformatics/btm595. [http://dx.doi.org/10.1093/bioinformatics/btm595], 10.1093/bioinformatics/btm595View ArticleGoogle Scholar
- Krek A, Grün D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N: Combinatorial microRNA target predictions. Nat Genet. 2005, 37 (5): 495-500. 10.1038/ng1536. [http://dx.doi.org/10.1038/ng1536], 10.1038/ng1536View ArticleGoogle Scholar
- Maragkakis M, Reczko M, Simossis VA, Alexiou P, Papadopoulos GL, Dalamagas T, Giannopoulos G, Goumas G, Koukis E, Kourtis K, Vergoulis T, Koziris N, Sellis T, Tsanakas P, Hatzigeorgiou AG: DIANA-microT web server: elucidating microRNA functions through target prediction. Nucleic Acids Res. 2009, 37 (Web Server issue): W273-W276. [http://dx.doi.org/10.1093/nar/gkp292],View ArticleGoogle Scholar
- Bandyopadhyay S, Mitra R: TargetMiner: microRNA target prediction with systematic identification of tissue-specific negative examples. Bioinformatics. 2009, 25 (20): 2625-2631. 10.1093/bioinformatics/btp503. [http://dx.doi.org/10.1093/bioinformatics/btp503], 10.1093/bioinformatics/btp503View ArticleGoogle Scholar
- Dekel E, Alon U: Optimality and evolutionary tuning of the expression level of a protein. Nature. 2005, 436 (7050): 588-592. 10.1038/nature03842. [http://dx.doi.org/10.1038/nature03842], 10.1038/nature03842View ArticleGoogle Scholar
- Gillespie D: A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. J Comput Phys. 1976, 22: 403-434. 10.1016/0021-9991(76)90041-3.View ArticleGoogle Scholar