- Open Access
Individual mRNA expression profiles reveal the effects of specific microRNAs
Genome Biology volume 9, Article number: R82 (2008)
MicroRNAs (miRNAs) are oligoribonucleotides with an important role in regulation of gene expression at the level of translation. Despite imperfect target complementarity, they can also significantly reduce mRNA levels. The validity of miRNA target gene predictions is difficult to assess at the protein level. We sought, therefore, to determine whether a general lowering of predicted target gene mRNA expression by endogenous miRNAs was detectable within microarray gene expression profiles.
The target gene sets predicted for each miRNA were mapped onto known gene expression data from a range of tissues. Whether considering mean absolute target gene expression, rank sum tests or 'ranked ratios', many miRNAs with significantly reduced target gene expression corresponded to those known to be expressed in the cognate tissue. Expression levels of miRNAs with reduced target mRNA levels were higher than those of miRNAs with no detectable effect on mRNA expression. Analysis of microarray data gathered after artificial perturbation of expression of a specific miRNA confirmed the predicted increase or decrease in influence of the altered miRNA upon mRNA levels. Strongest associations were observed with targets predicted by TargetScan.
We have demonstrated that the effect of a miRNA on its target mRNAs' levels can be measured within a single gene expression profile. This emphasizes the extent of this mode of regulation in vivo and confirms that many of the predicted miRNA-mRNA interactions are correct. The success of this approach has revealed the vast potential for extracting information about miRNA function from gene expression profiles.
MicroRNAs (miRNAs) are short oligonucleotides (approximately 22 bp) that regulate gene expression. Target genes are determined by sequence complementarity between the 3' untranslated region (UTR) and the mature miRNA, particularly in a 6 bp 'seed' region [1, 2]. A range of algorithms have been developed to predict the genes targeted by specific miRNAs . For example, 'TargetScan' [4, 5] searches for conserved 8-mer and 7-mer sites in 3' UTRs that match the seed region of a known miRNA. It is possible, therefore, to obtain lists of the potential target mRNAs for each miRNA. Plant miRNAs, which are often perfectly matched to their target sequences, act primarily by directing mRNA cleavage and degradation [6, 7]. In contrast, animal miRNAs have been shown to exert their effect largely via post-transcriptional inhibition of protein synthesis . However, it has been shown that miRNAs expressed in animal cells can affect mRNA levels, not only when they share almost complete complementarity with their target site , but more generally when base-pairing is partial [10–12]. When, for example, miR-124, which is known to be characteristic of neuronal tissue, was overexpressed, the genes that were down-regulated at the mRNA level included a preponderance of those expressed at lower levels in neuronal compared to other tissues . Conversely, silencing of miR-122 with a complementary, single-stranded RNA analogue, or 'antagomir', resulted in increased expression of mRNAs that were enriched in miR-122 recognition motifs  and miR-122 can direct cleavage of a reporter gene. Depletion of proteins required for miRNA processing has been shown to cause widespread alteration in mRNA levels [14, 15].
The suggestion that miRNAs can affect mRNA levels led to the prediction that a miRNA expressed at a high level in a specific tissue might leave a signature on the mRNA expression profile. Sood et al.  and Farh et al.  demonstrated that the predicted target genes of known tissue-specific miRNAs (for example, miR-122 in liver; miR-1 in heart/skeletal muscle and miR-7 in pituitary) were expressed at significantly lower levels, as determined by microarray analysis, in their cognate tissue relative to all other tissues.
The conclusive demonstrations that miRNAs can alter mRNA levels suggested to us that, within a specific tissue, the expression of genes predicted to be targeted by a specific mature miRNA might have a detectable inverse relationship with the expression level of that miRNA. This approach has been made feasible by advances in microarray technology and provision of comprehensive gene coverage, which have made global gene expression data increasingly reliable and reproducible [18, 19]. Concomitantly, public repositories such as Gene Expression Omnibus (GEO) [20, 21] and ArrayExpress  have made data from a huge range of tissues available to the scientific community. A method for extracting miRNA signatures from an mRNA expression dataset would be invaluable because it could immediately be applied to analyze miRNA activity in any situation for which microarray gene expression data are available.
Others have had limited success in detecting a significant effect of miRNAs within a single gene expression profile using a non-parametric approach based on gene expression ranking . However, by employing different predicted target and control datasets we were able to observe significant miRNA effects using a similar approach and by direct analysis of absolute target gene expression values (by the term 'target gene expression' we refer to the expression of predicted target genes). We were able to predict many of the previously characterized, highly expressed and/or tissue-specific miRNAs (for example, 14 of 25 in brain). This approach will facilitate investigation of the activity of miRNAs upon mRNA expression, without the need for ranking gene expression of each gene across a series of tissues [11, 23].
Results and discussion
Detection of miRNA signatures within endogenous gene expression profiles
miRNAs can down-regulate target mRNAs; therefore, one would expect that the target genes of a highly expressed miRNA might be expressed at a significantly lower level than those of a lowly expressed miRNA. In this case it might be possible to detect the presence of miRNAs from the relative expression of their predicted target genes. The profile of miRNAs expressed in one tissue differs from that in another and to test whether different 'signatures' were detectable, we first downloaded mRNA expression profiles for a range of tissues from GEO [20, 21]. A range of algorithms have been developed to predict miRNA target genes . However, the scarcity of experimentally confirmed interactions has made it difficult to develop reliable algorithms and validate existing methods. The relevance of existing rules is uncertain  and additional factors such as co-factor binding and relative positions of target sites [25, 26] undoubtedly play a role. Of the publicly available algorithms, we chose to initially use TargetScan  because its requirement for a perfect match to the seed region and cross-species conservation reduce the false-positive rate [3–5]. The resulting higher specificity of this algorithm maximizes the ability to detect effects on expression of real miRNA target genes. After detecting a signal we subsequently tested alternative miRNA target gene prediction algorithms (see below). For every mRNA expression dataset, the mRNA expression of predicted targets was mapped onto the respective miRNA families. The average number of predicted target genes for a single miRNA expressed in a given tissue was 134 (± 9 standard error; the number of predicted target genes for each miRNA expressed in all tissues is shown in Additional data file 1). We then tested the ability of three analytical approaches to detect the effects of variable endogenous miRNA expression on mRNA levels.
Wilcoxon rank sum test
Our first analysis followed the 'tissue-centric' approach described by Sood and colleagues . A vector of expression values for each set of specific miRNA target genes was compared to a vector of expression values for all predicted target genes. For all tissues, miRNAs with significantly low target gene expression were detected (Wilcoxon rank sum test), with lowest p-values ranging from 1.29 × 10-5 in brain to 7.23 × 10-3 in skeletal muscle. The results for all tissues are shown in Figure 1. It is notable that well characterized tissue-specific miRNAs, such as miR-122 in liver and miR-124 in brain, are all detected in the expected tissue and not elsewhere. This suggests that the 'signature' that a miRNA exerts upon mRNA expression can be detected within a single gene expression profile, without relation to levels in other tissues as previously reported .
In an alternative approach to analyzing the relative expression levels of all the predicted target genes of each miRNA within a particular tissue, we adapted the 'ranked ratio' (RR) described by Yu et al. . They first ranked the expression levels of each gene across a series of tissues. For each tissue the ranked genes were divided into two halves, one with high and one with low ranks. The RR values were then calculated by dividing the number of targeted genes in the 'low' ranked group by the 'high' ranked group. Instead of considering a range of tissues we ranked the targeted genes within a single expression dataset and for each miRNA calculated an RR value by dividing the number of predicted target genes with expression levels below the median absolute expression value by the number of predicted target genes above this value (comparison with other methods suggested that this was more effective than dividing genes into upper and lower halves - see below). This RR value is, therefore, an indicator of the distribution of a miRNA's target genes within a single mRNA population. A high RR indicates low expression in a greater proportion of target genes and is, therefore, indicative of miRNA expression in that tissue. The RR values for all miRNAs were calculated for all eight tissues and the ranked RR values for brain and liver are shown in Figure 2 (for all other tissues analyzed, see Additional data file 2). As expected, known tissue-specific miRNAs have high RR values in their cognate tissue.
Mean absolute expression
We next investigated whether an approach involving absolute target gene expression could be used to detect miRNA signatures. This could potentially identify miRNAs missed above, but runs the risk of being unduly influenced by single genes with a large change in expression. The technique is outlined in Figure 3a. The miRNAs were ordered by the mean expression value of their predicted target mRNAs, as shown for liver in Figure 3b. Of all the miRNAs in the liver, the lowest mean target gene expression value was that of miR-122a, a well characterized liver-specific miRNA . To determine the likelihood that this observed reduction in mRNA expression is due to the selection of mRNAs with specific miRNA targets, we calculated the probability (t-test) that these samples are drawn at random from amongst all those genes expressed in the tissue and that contain a predicted miRNA target sequence. The resulting probabilities for all tissues are plotted in Figure 3c and those miRNAs with low target gene expression include many known tissue-specific examples, such as miR-124a in brain (p = 6.2 × 10-4) and miR-1 in skeletal muscle (p = 1.9 × 10-2).
To test the reliability of this approach and the robustness of available microarray expression data, it was applied to an independent mouse expression dataset generated in several different laboratories (see Materials and methods). The sets of miRNAs predicted from the two datasets were very similar for all tissues and the extent of overlap is depicted in Figure 4. Mammalian miRNAs and their target sites are highly conserved; indeed, sequence conservation is a requirement of the TargetScan predictions . Accordingly, miRNA expression is conserved between species , at least for organisms with similar physiology  and miRNAs may have a role in reducing cross-species variation in mRNA expression . Analysis of mRNA expression profiles from human tissues (Additional data file 3) revealed that approximately one-third of the human miRNAs with low target gene expression corresponded to those predicted in equivalent murine tissues (30.6% and 35.5% for mouse datasets 1 and 2, respectively). For example, of 18 human miRNAs predicted in brain, 7 were common with mouse dataset 1 (Figure 4), rather than the lower number (approximately 2) expected if the groups of miRNAs were independent (the observed numbers were similarly high for all other tissues and the second mouse dataset). This provides further evidence for conserved miRNA expression and independent validation of the prediction method.
Comparison of miRNA signature detection methods
We next compared the results of the three methods, Wilcoxon rank sum test, RR and absolute expression t-test, using a 10% significance level and an equivalent number of miRNAs from the RR method. For all tissues there was significant overlap amongst predicted miRNAs (Figure 5), with the Wilcoxon rank sum test and absolute expression t-test in strongest agreement. To evaluate how well the miRNA signature detected in target gene expression predicts actual miRNA expression, we compared the tissue distribution of miRNAs predicted by at least two of the methods with that derived from experimental evidence (cloning and Northern blots). Table 1 illustrates the accordance between the tissues in which miRNA activity (upon target genes predicted by TargetScan) is computationally predicted and those for which there is experimental evidence of miRNA presence (particularly when more recently characterized miRNAs are excluded). This is supported by positive Matthews correlation coefficients (MCC)  for all tissues, ranging from 0.2-0.5 (average value 0.34; Additional data file 4).
Correlation of miRNA signatures with miRNA expression levels
miRNA microarrays are now available that provide a global indication of miRNA expression within a tissue. We therefore compared our predictions of miRNAs that alter mRNA expression with the actual expression of the miRNAs themselves, as determined by miRNA microarrays . For all tissues the expression levels of miRNAs with low target gene expression, determined by the absolute expression method (10% significance level), were significantly lower (t-test, p < 0.05) than those miRNAs having no detectable effect on their target genes (Figure 6). This provides further confirmation that the miRNA signatures we have detected are a consequence of miRNA expression in the cognate tissue. In addition to miRNAs with low target gene expression, we detected a set of miRNAs whose target genes were expressed at significantly higher levels than the background set (Figure 3b). The expression of these miRNAs, as determined by microarrays, was not significantly different from those with no effect on mRNA expression.
Recently, comprehensive miRNA expression data for human tissues determined by reverse transcription PCR (RT-PCR) have become available . This revealed an even clearer relationship between human miRNA copy number and level of predicted target gene expression (Figure 7). In an attempt to demonstrate the similarity between human and mouse miRNAs, the human orthologs of miRNAs predicted from analysis of murine data to have 'low' or 'mid' predicted target gene expression were selected. Surprisingly, significant differences were detected between the copy numbers in human tissues of these two groups of miRNAs, which had been selected based upon murine target gene expression (Additional data file 5). This is testament to the degree of conservation of miRNAs and their target genes between mice and humans and the accuracy of the RT-PCR measurements.
Some of those miRNAs with a significant effect on mRNA expression were highly expressed (for example, miR-122a, miR-124a, miR-125) and their observed lowering of mRNA levels could reasonably be attributed to a weak mRNA degradative activity secondary to their principal action directed at translation. However, other miRNAs that significantly affected target gene expression were not highly expressed, perhaps indicating a greater efficiency in mRNA degradation for these particular miRNAs. Therefore, the extent to which specific miRNAs cause mRNA degradation might be influencing our ability to detect their presence. We reasoned that the difference between miRNAs would be most marked between those highly expressed but having no detectable affect on mRNA expression and those expressed at a low level but with a significant impact on target mRNA expression. Other than extensive complementarity , the features of the miRNA-target interaction required for miRNAs to direct mRNA cleavage are unclear, although a number of features of site context, including position, local AU content and pairing with miRNA 3' residues have been shown to increase site efficacy . In a preliminary attempt to characterize the distinguishing properties of the potential classes of miRNAs described above, we analyzed the lengths of contiguous complementarity between miRNAs and predicted sites, but there were no significant differences.
Perturbation of miRNA expression affects cognate mRNA expression
In order to validate the miRNA signatures observed in mRNA profiles of normal tissues, we applied our approach to several publicly available gene expression datasets measured following an experimental perturbation of miRNA expression that resulted in either a decrease  or increase  in the activity of a specific miRNA. In each case the expected response was observed in target gene expression. Krutzfeldt et al. demonstrated that intravenous injection of chemically modified oligonucleotides, or 'antagomirs', complementary to miRNAs could specifically reduce the endogenous levels of the corresponding miRNA in mice. When we analyzed gene expression data from liver in which miR-122a had been silenced by use of an antagomir  the list of miRNAs with an influence on mRNA levels was very similar but miR-122a no longer had a detectable effect (Figure 8). This conclusively demonstrates that the signature of miR-122a on target gene expression that we observe is due to the physical presence of miR-122a rather than any evolutionary pressure on the expression of target genes co-expressed with their cognate miRNA.
Wang and Wang  transfected the HepG2 cell line with RNA duplexes that mimicked the miR-124 precursor. They demonstrated the effect of miR-124 over-expression by measuring mRNA expression profiles and showing that predicted targets of miR-124 were over-represented amongst the mRNAs down-regulated following this treatment, compared to a negative control RNA duplex. As expected, analysis of miRNA signatures in these expression datasets revealed miR-124 activity in cells in which it was overexpressed but not in controls (Figure 8).
Characteristics of predicted target gene groups
We could not find any increase in target complementarity for miRNAs with a significant effect on mRNA levels. However, we noticed that those miRNAs with low overall target gene expression in a specific tissue generally had a greater number of target genes expressed in that tissue. This was true across tissues and supports the observations of Farh et al. , who demonstrated that tissue-specific miRNA target genes are generally expressed in the cognate tissue, but at lower levels than other tissues.
In certain tissues some miRNA predicted target gene sets were expressed at significantly higher levels (10% probability) than those of other miRNAs. We observed that these often involved a highly expressed, tissue-specific miRNA in tissues in which the miRNA was not normally expressed, for example, miR-1 in liver (Figure 2). However, this phenomenon does not appear to be indicative of reduced suppression of miRNA levels due to low miRNA expression in these tissues (Figures 6 and 7).
The targeting of mRNAs by miRNAs has been associated with 3' UTR GC content and length [16, 35]. However, we found no consistent differences between the 3' UTR GC content of predicted target genes of miRNAs with significantly low target gene expression ('low' genes) and those of background genes. However, for all tissues, the predicted target genes of miRNAs with significantly high target gene expression ('high' genes) had significantly lower GC content (Additional data file 6). The 3' UTR length was significantly longer than background for 'low' genes, but 'high' genes were even longer. The biological significance of these observations is unclear. One might have expected genes with lower expression to be more AU rich because miRNAs have been implicated in degradation of mRNAs containing AU-rich elements . The longer 3' UTRs of 'low' genes would be predicted to have more miRNA target sites and, therefore, be subject to more degradation, but this is inconsistent with the even longer 3' UTRs observed for 'high' genes.
Use of alternative target prediction methods
The correlation demonstrated between expression of certain miRNAs and that of their predicted target genes indicates that at least some of the interactions suggested by TargetScan are valid. Random assignment of miRNA target genes rather than use of TargetScan predictions resulted in fewer miRNAs with low target gene expression and significantly less overlap between those miRNAs with low target gene expression in replicate expression datasets (p < 0.00001).
We next investigated whether we could detect variations in expression of the targets predicted by alternative algorithms. Many of the highly expressed, known tissue-specific miRNAs identified above (Table 1) were not observed when the equivalent analyses were performed with target gene sets predicted by miRanda [36, 37] or RNAhybrid [38, 39]. For example, in brain the TargetScan-based analyses predicted 14 of 24 known miRNAs whilst miRanda and RNAhybrid predicted 4 and 5, respectively. In Liver, TargetScan, miRanda and RNAhybrid predicted five, one and one, respectively, out of eight, and in skeletal muscle five, three and one out of five. The predicted targets of miR-122a were not lowly expressed in liver and, therefore, no effect could be observed following treatment with antagomir. Likewise, the expression levels of the predicted targets of miR-124 were not significantly lower following its overexpression. However, the miRNAs with low predicted target gene expression in a particular tissue were expressed at higher levels than others in that tissue (Additional data files 7 and 8), albeit less markedly than those miRNAs derived from equivalent analyses with TargetScan predictions (Additional data file 5). This suggests that these prediction methods do detect real miRNA targets, but that there may be more false predictions and/or that the target genes detected by these criteria are less susceptible to miRNA-mediated reduction in mRNA levels. The miRNAs detected in the analyses based upon miRanda or RNAhybrid target predictions are listed in Additional data file 9.
It has been suggested that target mRNAs can be repressed by miRNA binding sites in the 5' UTR . Analysis of all mouse 5' UTRs for target sites using RNAhybrid yielded fewer predictions than in 3' UTRs and there was no obvious link between miRNAs with low target gene expression and known, highly expressed, tissue-specific miRNAs such as miR-1 in heart or miR-122a in liver. However, there is evidence from the miRNA expression data of Liang et al.  that those miRNAs with low predicted target gene expression are more highly expressed (Additional data file 10). This suggests that functional interactions are occurring between miRNAs and 5' UTRs in vivo.
It is perhaps surprising that given the wide range of regulatory mechanisms acting on gene expression, we are able to detect the effect of miRNAs on their diverse collection of target genes. The effectiveness of three analytical approaches, including the use of absolute expression values that might be unevenly skewed by single extreme values, illustrates the extent of the effect and the correlation between independent datasets is remarkable. These findings provide further evidence that miRNAs play an important role in regulating mRNA expression and can affect the mRNA levels of multiple genes. Although the effects on mRNA expression profiles are highlighted following artificial manipulation of miRNA levels [12, 13], these are real in vivo effects because they are observed in data derived from endogenous situations. One limitation of the TargetScan predictions that we used, is that they cannot discriminate between miRNA family members with the same seed sequence (this has been addressed in release 4.0).
The mechanism by which miRNAs are causing degradation of their cognate mRNAs is unclear [8, 41] and may be a side-effect of the translation inhibition process and be related to high miRNA expression levels. In rare cases miRNAs can direct slicer-mediated cleavage if there is extensive complementarity with the target site, as reported for miR-196 and HOXB8 mRNA . This is not the case for the miRNA-target interactions reported here, but they may be inducing deadenylation and subsequent cleavage . The ability of specific miRNAs to reduce the level of certain target mRNAs very efficiently might explain why they have a detectable effect upon gene expression in tissues in which their expression has been shown to be very low, for example, miR-196 in brain. This raises the possibility that predicted miR-196 target sites included in this analysis may have distinctive features. An updated release of TargetScan (4.0) includes scores for features that contribute to predicted target site efficacy, which are summarized in a 'context score' . The context scores for the predicted targets of miR196 expressed in brain are not significantly higher than those of other miRNAs. However, there are several predicted target genes, including HOXB7 and HOXA5, with high context scores - efficient degradation of these targets could contribute to a detectable signal from low levels of miR-196. Analyses of miRNA signatures in gene expression profiles will help to pinpoint those miRNAs acting at the mRNA level and determine to what extent this results from high miRNA expression or reflects specific mechanisms related to certain miRNAs.
It is difficult to validate miRNA target predictions from effects on translation due to the technical difficulties of genome-wide quantitative protein profiling. Analysis of the effects of miRNAs on target mRNA expression provides a potential alternative. The target predictions we employed [27, 37, 39] were sufficiently accurate to enable detection of miRNA signatures within gene expression profiles. The ability to detect these miRNA signatures also attests to the sensitivity and reproducibility of current microarray technology. The strength of signal detected could be used to assess the effectiveness of different target prediction algorithms. The lack of significant effects in previous studies  might be because the target gene predictions employed were less effective at predicting miRNA-target interactions that result in mRNA degradation.
A unique advantage of this approach is that it is possible to infer miRNA activity from a single gene expression profile rather than having to first rank the expression of each gene across a series of experiments [11, 16, 17]. This will facilitate the analysis of individual datasets of interest from the vast repository of gene expression data that are steadily accumulating [20, 22]. It is applicable to data from any sufficiently sensitive and comprehensive platform and to gene expression profiles generated from biological sources too scarce to permit direct miRNA detection.
Our analysis of a limited number of gene expression profiles has demonstrated the major impact of miRNAs upon mRNA expression. The programs are amenable to provision of an online service that would enable individual researchers to gain an indication of miRNA activity in their tissue/cells of interest. This is a convenient approach to gain additional biological insight from gene expression profiles.
Materials and methods
Stand Alone Java and R programs , in conjugation with MS Access and MS Excel, were used for the following analyses.
Gene expression data
Mouse and human mRNA expression datasets that covered a range of tissues were downloaded from GEO [20, 21, 44]. In addition, a second mouse dataset covering the same tissues was compiled from data submitted by different laboratories. The tissues and GEO sample and platform accession numbers are listed in Additional data file 3. Evidence for tissue-specific miRNA expression supported by cloning, Northern hybridization or expressed sequence tag mapping was determined by literature review.
Calculation of miRNA target gene expression
Only single gene-specific Affymetrix probesets (suffix_at) were considered and only those expressed in each dataset analyzed (designated by a 'present' call if available, or alternatively a signal strength greater than the median value). Probeset IDs were converted to their cognate gene symbols, using the Biomart online suite . Lists of miRNA families and their predicted target genes published by Lewis et al.  were downloaded from the TargetScan website (version 3.1) . For mouse (Mouse Genome 430 2.0 genechip), 15,682 probesets were mapped to gene symbols and of these, between 3,120 and 3,792 were both predicted miRNA targets and expressed in the tissues analyzed (for human, 17,325 probesets and 4,135-5,906 expressed targets were analyzed). Additional miRNA target genes predicted by miRanda [36, 37] were downloaded from the miRBase website . Sites with scores and conservation values (P_org) greater than average were selected for analysis. The RNAhybrid algorithm  was run with the following parameters: p < 0.05, energy < -10 kcal/mol and with enforced binding to bases 2-7.
For each mRNA expression dataset, lists of the expressed genes predicted to be targeted by each miRNA and their respective expression levels were compiled (the average value was calculated for genes represented by > 1 probeset). Only miRNAs with > 50 predicted targets were considered further. These were sorted according to the average expression values (absolute or logarithmically transformed according to the primary microarray normalization algorithm used to analyze the data) of their predicted target genes and when necessary divided into upper and lower halves.
Three different analytical methods were used to compare the expression of the predicted targets of a specific miRNA in a particular tissue, with the set of expression values of all the genes predicted to be targeted by miRNAs in the same tissue. To analyze whether there was significant difference between the medians of the ranked gene expression values of the two sets, the nonparametric one-sided Wilcoxon rank sum test was employed. To test the hypothesis that the genes in each set are from the same population, we used a one-tailed Student's t-test (considering unequal variance and using log-normalized data). For every miRNA with p < 0.1, we considered that the expression of all the genes predicted to be targeted by that miRNA in that tissue was significantly more or less than the average expression of all the predicted target genes.
To calculate an RR value, predicted target genes were ordered by expression value and then divided into two groups by either the median gene rank or absolute expression value. For each miRNA, the number of genes targeted in the lower group was divided by the number in the higher to give an RR value.
MCCs  were calculated as a measure of performance.
Calculation of GC content and length of predicted target mRNA 3' UTRs
Mouse 3' UTR data were downloaded from Biomart . The GC content and the lengths of the UTRs were then calculated for the nonredundant sets of genes predicted to be targeted by groups of miRNAs of interest.
Analysis of miRNA-target site complementarity
For sets of miRNAs of interest, the number of predicted targeted genes having greater than or equal to x contiguous complementary nucleotides after the 6 bp seed was noted (where x ranged from 0-8, 0 being the seed itself and 8 the maximum complementarity observed after the seed region). The number of predicted target genes was calculated with and without consideration of multiple miRNA family members.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 is a figure showing the number of target genes predicted by TargetScan for each miRNA and expressed in each mouse tissue. Additional data file 2 is a figure displaying the RR values for all miRNAs in heart, kidney, lung, ovary, skeletal muscle and testes. Additional data file 3 is a table listing the sources of gene expression data from the NCBI GEO website. Additional data file 4 is a table listing the MCCs calculated for each tissue. Additional data file 5 is a series of box plots showing the correlation between mouse miRNAs with lowly expressed mouse target genes predicted by TargetScan and expression levels of orthologous human miRNAs detected in cognate tissues by RT-PCR. Additional data file 6 is a figure showing the GC content and lengths of the 3' UTRs of miRNA predicted target mRNAs in a range of tissues. Additional data file 7 is a series of box plots showing the correlation between miRNAs with lowly expressed target genes predicted by miRanda and miRNA expression levels detected by RT-PCR. Additional data file 8 is a series of box plots showing the correlation between miRNAs with lowly expressed target genes predicted by RNAhybrid and miRNA expression levels detected by RT-PCR. Additional data file 9 is a table listing miRNAs with significantly low target gene expression determined by miRanda and RNAhybrid.
Additional data file 10 is a series of box plots showing that miRNAs with lowly expressed predicted target genes, as defined by putative 5' UTR sites, are expressed at higher levels than those with no effect on target gene expression.
Gene Expression Omnibus
Matthews correlation coefficients
reverse transcription PCR
Bartel DP: MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.
Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: Determinants beyond seed pairing. Mol Cell. 2007, 27: 91-105. 10.1016/j.molcel.2007.06.017.
Maziere P, Enright AJ: Prediction of microRNA targets. Drug Discov Today. 2007, 12: 452-458. 10.1016/j.drudis.2007.04.002.
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.
Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115: 787-798. 10.1016/S0092-8674(03)01018-3.
Tang G, Reinhart BJ, Bartel DP, Zamore PD: A biochemical framework for RNA silencing in plants. Genes Dev. 2003, 17: 49-63. 10.1101/gad.1048103.
Llave C, Xie Z, Kasschau KD, Carrington JC: Cleavage of scarecrow-like mRNA targets directed by a class of arabidopsis miRNA. Science. 2002, 297: 2053-2056. 10.1126/science.1076311.
Pillai RS, Bhattacharyya SN, Filipowicz W: Repression of protein synthesis by miRNAs: How many mechanisms?. Trends Cell Biol. 2007, 17: 118-126. 10.1016/j.tcb.2006.12.007.
Yekta S, Shih IH, Bartel DP: MicroRNA-directed cleavage of HOXB8 mRNA. Science. 2004, 304: 594-596. 10.1126/science.1097434.
Bagga S, Bracht J, Hunter S, Massirer K, Holtz J, Eachus R, Pasquinelli AE: Regulation by let-7 and lin-4 miRNAs results in target mRNA degradation. Cell. 2005, 122: 553-563. 10.1016/j.cell.2005.07.031.
Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM: Microarray analysis shows that some microRNAs down-regulate large numbers of target mRNAs. Nature. 2005, 433: 769-773. 10.1038/nature03315.
Wang X, Wang X: Systematic identification of microRNA functions by combining target prediction and expression profiling. Nucleic Acids Res. 2006, 34: 1646-1652. 10.1093/nar/gkl068.
Krutzfeldt J, Rajewsky N, Braich R, Rajeev KG, Tuschl T, Manoharan M, Stoffel M: Silencing of microRNAs in vivo with 'antagomirs'. Nature. 2005, 438: 685-689. 10.1038/nature04303.
Schmitter D, Filkowski J, Sewer A, Pillai RS, Oakeley EJ, Zavolan M, Svoboda P, Filipowicz W: Effects of Dicer and Argonaute down-regulation on mRNA levels in human HEK293 cells. Nucleic Acids Res. 2006, 34: 4801-4815. 10.1093/nar/gkl646.
Rehwinkel J, Natalin P, Stark A, Brennecke J, Cohen SM, Izaurralde E: Genome-wide analysis of mRNAs regulated by Drosha and Argonaute proteins in Drosophila melanogaster . Mol Cell Biol. 2006, 26: 2965-2975. 10.1128/MCB.26.8.2965-2975.2006.
Sood P, Krek A, Zavolan M, Macino G, Rajewsky N: Cell-type-specific signatures of microRNAs on target mRNA expression. Proc Natl Acad Sci USA. 2006, 103: 2746-2751. 10.1073/pnas.0511045103.
Farh KK, Grimson A, Jan C, Lewis BP, Johnston WK, Lim LP, Burge CB, Bartel DP: The widespread impact of mammalian microRNAs on mRNA repression and evolution. Science. 2005, 310: 1817-1821. 10.1126/science.1121158.
MAQC Consortium, Shi L, Reid LH, Jones WD, Shippy R, Warrington JA, Baker SC, Collins PJ, de Longueville F, Kawasaki ES, Lee KY, Luo Y, Sun YA, Willey JC, Setterquist RA, Fischer GM, Tong W, Dragan YP, Dix DJ, Frueh FW, Goodsaid FM, Herman D, Jensen RV, Johnson CD, Lobenhofer EK, Puri RK, Schrf U, Thierry-Mieg J, Wang C, Wilson M, et al: The MicroArray quality control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol. 2006, 24: 1151-1161. 10.1038/nbt1239.
Kuo WP, Liu F, Trimarchi J, Punzo C, Lombardi M, Sarang J, Whipple ME, Maysuria M, Serikawa K, Lee SY, McCrann D, Kang J, Shearstone JR, Burke J, Park DJ, Wang X, Rector TL, Ricciardi-Castagnoli P, Perrin S, Choi S, Bumgarner R, Kim JH, Short GF, Freeman MW, Seed B, Jensen R, Church GM, Hovig E, Cepko CL, Park P, et al: A sequence-oriented comparison of gene expression measurements across different hybridization-based technologies. Nat Biotechnol. 2006, 24: 832-840. 10.1038/nbt1217.
Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: Mining tens of millions of expression profiles - database and tools update. Nucleic Acids Res. 2007, 35: D760-D765. 10.1093/nar/gkl887.
NCBI Gene Expression Omnibus. [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/projects/geo/]
Parkinson H, Kapushesky M, Shojatalab M, Abeygunawardena N, Coulson R, Farne A, Holloway E, Kolesnykov N, Lilja P, Lukk M, Mani R, Rayner T, Sharma A, William E, Sarkans U, Brazma A: ArrayExpress - a public database of microarray experiments and gene expression profiles. Nucleic Acids Res. 2007, 35: D747-D750. 10.1093/nar/gkl995.
Yu Z, Jian Z, Shen SH, Purisima E, Wang E: Global analysis of microRNA target gene expression reveals that miRNA targets are lower expressed in mature mouse and drosophila tissues than in the embryos. Nucleic Acids Res. 2007, 35: 152-164. 10.1093/nar/gkl1032.
Didiano D, Hobert O: Perfect seed pairing is not a generally reliable predictor for miRNA-target interactions. Nat Struct Mol Biol. 2006, 13: 849-851. 10.1038/nsmb1138.
Saetrom P, Heale BS, Snove O, Aagaard L, Alluin J, Rossi JJ: Distance constraints between microRNA target sites dictate efficacy and cooperativity. Nucleic Acids Res. 2007, 35: 2333-2342. 10.1093/nar/gkm133.
Hon LS, Zhang Z: The roles of binding site arrangement and combinatorial targeting in microRNA repression of gene expression. Genome Biol. 2007, 8: R166-10.1186/gb-2007-8-8-r166.
TargetScan: Prediction of MicroRNA Targets. [http://www.targetscan.org/]
Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, Tuschl T: Identification of tissue-specific microRNAs from mouse. Curr Biol. 2002, 12: 735-739. 10.1016/S0960-9822(02)00809-6.
Sempere LF, Freemantle S, Pitha-Rowe I, Moss E, Dmitrovsky E, Ambros V: Expression profiling of mammalian microRNAs uncovers a subset of brain-expressed microRNAs with possible roles in murine and human neuronal differentiation. Genome Biol. 2004, 5: R13-10.1186/gb-2004-5-3-r13.
Ason B, Darnell DK, Wittbrodt B, Berezikov E, Kloosterman WP, Wittbrodt J, Antin PB, Plasterk RH: Differences in vertebrate microRNA expression. Proc Natl Acad Sci USA. 2006, 103: 14385-14389. 10.1073/pnas.0603529103.
Cui Q, Yu Z, Purisima EO, Wang E: MicroRNA regulation and interspecific variation of gene expression. Trends Genet. 2007, 23: 372-375. 10.1016/j.tig.2007.04.003.
Matthews BW: Comparison of the predicted and observed secondary structure of T4 phage lysozyme. Biochim Biophys Acta. 1975, 405: 442-451.
Thomson JM, Parker J, Perou CM, Hammond SM: A custom microarray platform for analysis of microRNA gene expression. Nat Methods. 2004, 1: 47-53. 10.1038/nmeth704.
Liang Y, Ridzon D, Wong L, Chen C: Characterization of microRNA expression profiles in normal human tissues. BMC Genomics. 2007, 8: 166-10.1186/1471-2164-8-166.
Jing Q, Huang S, Guth S, Zarubin T, Motoyama A, Chen J, Di Padova F, Lin SC, Gram H, Han J: Involvement of microRNA in AU-rich element-mediated mRNA instability. Cell. 2005, 120: 623-634. 10.1016/j.cell.2004.12.038.
John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human microRNA targets. PLoS Biol. 2004, 2: e363-10.1371/journal.pbio.0020363.
Human MicroRNA Targets. [http://www.microrna.org/]
Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes. RNA. 2004, 10: 1507-1517. 10.1261/rna.5248604.
RNAhybrid Web Tool on the Bielefeld Bioinformatics Server. [http://bibiserv.techfak.uni-bielefeld.de/rnahybrid/]
Lytle JR, Yario TA, Steitz JA: Target mRNAs are repressed as efficiently by microRNA-binding sites in the 5' UTR as in the 3' UTR. Proc Natl Acad Sci USA. 2007, 104: 9667-9672. 10.1073/pnas.0703820104.
Jackson RJ, Standart N: How do microRNAs regulate gene expression?. Sci STKE. 2007, 2007: re1-10.1126/stke.3672007re1.
Giraldez AJ, Mishima Y, Rihel J, Grocock RJ, Van Dongen S, Inoue K, Enright AJ, Schier AF: Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science. 2006, 312: 75-79. 10.1126/science.1122689.
The Comprehensive R Archive Network. [http://cran.r-project.org/]
Edgar R, Domrachev M, Lash AE: Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.
The Biomart Project. [http://www.biomart.org/]
miRBase Targets. [http://microrna.sanger.ac.uk/targets/]
Gu J, He T, Pei Y, Li F, Wang X, Zhang J, Zhang X, Li Y: Primary transcripts and expressions of mammal intergenic microRNAs detected by mapping ESTs to their flanking sequences. Mamm Genome. 2006, 17: 1033-1041. 10.1007/s00335-006-0007-9.
Naraba H, Iwai N: Assessment of the microRNA system in salt-sensitive hypertension. Hypertens Res. 2005, 28: 819-826. 10.1291/hypres.28.819.
Zhao Y, Samal E, Srivastava D: Serum response factor regulates a muscle-specific microRNA that targets Hand2 during cardiogenesis. Nature. 2005, 436: 214-220. 10.1038/nature03817.
Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T: Identification of novel genes coding for small expressed RNAs. Science. 2001, 294: 853-858. 10.1126/science.1064921.
Lagos-Quintana M, Rauhut R, Meyer J, Borkhardt A, Tuschl T: New microRNAs from mouse and human. RNA. 2003, 9: 175-179. 10.1261/rna.2146903.
Hayashita Y, Osada H, Tatematsu Y, Yamada H, Yanagisawa K, Tomida S, Yatabe Y, Kawahara K, Sekido Y, Takahashi T: A polycistronic microRNA cluster, miR-17-92, is overexpressed in human lung cancers and enhances cell proliferation. Cancer Res. 2005, 65: 9628-9632. 10.1158/0008-5472.CAN-05-2352.
Yu Z, Raabe T, Hecht NB: MicroRNA Mirn122a reduces expression of the posttranscriptionally regulated germ cell transition protein 2 (Tnp2) messenger RNA (mRNA) by mRNA cleavage. Biol Reprod. 2005, 73: 427-433. 10.1095/biolreprod.105.040998.
We would like to acknowledge the contribution of all the researchers who have made their gene expression data publicly available and without whom this study would not have been possible. We are grateful to Drs Tom Gardiner, Gareth McKay, Tim Curtis and Colin Willoughby for helpful discussion. AA was supported by a SPUR studentship.
DACS conceived the study. AA collected the data and performed the data analyses. DACS and AA interpreted the results. DACS drafted the manuscript with assistance from AA. Both authors read and approved the final manuscript.
Electronic supplementary material
Additional data file 1: The number of predicted target genes (y-axis) for each miRNA (x-axis) is displayed with a different symbol for each tissue (grey circle, brain; green triangle, kidney; blue square, lung; green diamond, skeletal muscle; white circle, heart; red triangle, liver; purple square, ovary; yellow diamond, testes). miRNAs with less than 50 predicted target genes expressed in a specific tissue were not considered (see dashed line). (TIFF 204 KB)
Additional data file 2: Each panel displays the data from a single tissue. The miRNAs are arranged in ascending order according to RR value, shown on the left-hand y-axis and displayed as a red line. The higher values reflect lower target gene expression and are, therefore, indicative of miRNA activity. The number of predicted target genes is shown on the right-hand y-axis and the values for each miRNA indicated by a dashed line. (TIFF 1 MB)
Additional data file 5: Box plots depict the copy numbers (y-axis) of groups of microRNAs (x-axis) whose mouse orthologs have low, mid or high predicted target gene expression (as defined in the text). Where necessary to present the median and interquartile ranges effectively, up to two outliers were omitted. (TIFF 386 KB)
Additional data file 6: The graphs on each row show data from a single tissue, with the first column depicting 3' UTR length and the second GC content. The predicted target genes are divided into three groups, low (red), medium (orange) and high (yellow). Those in the low group are targeted by miRNAs with overall significantly low target gene expression, those in the medium group are targeted by miRNAs with overall target gene expression within the expected range and those in the high group are targeted by miRNAs with overall high target gene expression. The average lengths in bases or %GC content (y-axis) of the 3' UTRs of each group of genes (x-axis) are shown with standard error bars. (PNG 122 KB)
Additional data file 7: miRNAs were divided onto those with significantly lower than expected target mRNA expression (labeled 'low') and those with no detectable effect on their target expression (labeled 'mid'). The boxplots show the copy number (y-axis) of the miRNAs in each group (x-axis) and illustrate the significantly higher expression of miRNAs with low target gene expression. The scale of the y-axis (copy number) was chosen to facilitate visual comparison between groups and necessitated omission of up to three outliers from several graphs. (TIFF 336 KB)
Additional data file 8: miRNAs were divided into those with significantly lower than expected target mRNA expression (labeled 'low') and those with no detectable effect on the expression of their predicted targets (labeled 'mid'). The boxplots illustrate the significantly higher expression of miRNAs with low target gene expression. The scale of the y-axis (copy number) was chosen to facilitate visual comparison between groups and necessitated omission of up to three outliers from several graphs. (TIFF 342 KB)
Additional data file 10: miRNAs were divided into groups according to gene expression of targets predicted from potential 5' UTR miRNA binding sites. For each tissue (brain, heart, kidney, liver, lung, ovary, skeletal muscle and testes) a box plot shows the miRNA copy number, as determined by RT-PCR  (y-axis) of the miRNAs in each group ('low' and 'mid', x-axis). (TIFF 332 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Arora, A., Simpson, D.A. Individual mRNA expression profiles reveal the effects of specific microRNAs. Genome Biol 9, R82 (2008) doi:10.1186/gb-2008-9-5-r82
- miRNA Expression
- Additional Data File
- Target Gene Expression
- mRNA Expression Profile
- Predict Target Gene