Expression levels of microRNAs are not associated with their regulatory activities
© Liang et al; licensee BioMed Central Ltd. 2011
Received: 24 May 2011
Accepted: 19 September 2011
Published: 19 September 2011
MicroRNAs (miRNAs) regulate their targets by triggering mRNA degradation or translational repression. The negative relationship between miRNAs and their targets suggests that the regulatory effect of a miRNA could be determined from the expression levels of its targets. Here, we investigated the relationship between miRNA activities determined by computational programs and miRNA expression levels by using data in which both mRNA and miRNA expression from the same samples were measured. We found that different from the intuitive expectation one might have, miRNA activity shows very weak correlation with miRNA expression, which indicates complex regulating mechanisms between miRNAs and their target genes.
This manuscript was reviewed by an anonymous reviewer and Dr Yuriy Gusev.
KeywordsMicroRNA microRNA activity microRNA expression
MicroRNAs (miRNAs) as small non-coding RNAs have been widely recognized as crucial regulators in a broad range of biological processes . miRNAs regulate gene expressions by base-paring with the 3' untranslated region of target mRNAs, mainly leading to mRNA degradation or translational repression. Since one miRNA could regulate more than one target gene and one gene might be regulated by more than one miRNAs, the relationship between miRNAs and their target genes is quite complicated. Facing this challenge, several computational programs, including miReduce , MIR , Sylamer  and mirAct  have been developed to evaluate miRNA activity changes between different biological states according to gene-expression data. All these programs are based on a general consideration that an up-regulation of a miRNA's target expression should result from the decrease of its inhibitory activity, or vice versa. One could intuitively argue that the miRNA activities might be directly reflected by up- or down-regulation of miRNA expression levels. In order to test this argument, we investigated the relationship between miRNA activities determined by these four programs and miRNA expression levels by using data in which both mRNA and miRNA expression from the same samples were measured. Accordingly, two microarray data sets denoted as PNAS05  and Nature05  were analyzed. The data set PNAS05 contains matched expression data from 9 papillary thyroid carcinoma (PTC) samples and 9 unaffected thyroid tissue samples. The data set Nature05 contains matched mRNA and miRNA expression data from 8 different cancerous tissues and corresponding normal tissues.
Correlation coefficients between the raw p-values and the absolute of miRNA log-ratios
The present results suggest that the regulatory activity of miRNAs should depend on many different factors besides their expression levels. It has been shown that RNA accessory proteins such as RNA binding proteins (RBPs) can regulate miRNA activity without impacting miRNA expression levels and blur the relationship between miRNA activity and miRNA expression . For instance, a human-RBP Dnd1 regulated miRNA activity by binding mRNAs and determining access for miRNA regulators . Similarly, in E. coli, a protein Hfq was found to bind to mRNAs and determines sRNA activity . Furthermore, it was demonstrated that miRNA activity is also influenced by target abundance . Therefore, these discovered various factors impacting miRNA activity make the situation more complex.
In summary, different from the intuitive expectation one might have, miRNA activity shows very weak correlation with miRNA expression, which indicates complex regulating mechanisms between miRNAs and their target genes.
Reviewer's report 1
Using a computer program generated p-value as an indicator of miRNA regulatory activity and analysis of correlation, the authors conclude that the association between expression levels of miRNAs and their regulatory activities is very weak. One could argue that it is possible to reach the same conclusion by considering available evidences alone without having to rely on the statistical (correlation) analysis of two variables as done in this article. It should be noted that although p-values are generated by computer programs that analyze activity change of miRNAs, it is not clear what these values actually represent: either biologically or statistically. Biologically, it would be more meaningful to focus analyses only on those miRNAs that have similar regulatory activities, rather than lumping together all miRNAs in one analysis. It is also statistically problematic because the analyses of simple correlation failed to consider potential confounding factors such as age and sex for tissues that provide data for computer programs and the resultant p-values. The importance of considering confounding factors in this type of data analysis has been articulated in Potter (2003)* Therefore, the conclusion from the analysis should be considered superficial, and unnecessary.
*Potter J. Epidemiology, cancer genetics and microarrays: making correct inferences, using appropriate designs. Trends Genet 2003, 19(12):690-695.
Dear reviewer, thank you for your careful reviewing on our manuscript. The following are our responses to your comments.
First, the statistical analysis performed in the manuscript aims to identify the general relationship between miRNA activity and miRNA expression. Although some available evidences support the absence of association between the two, however some support the opposite. Also, current evidences are limited and most of them are case-by-case examples. So we think it might not be appropriate to make a general conclusion by considering available evidences alone, at least at the present time. Especially, it is difficult to measure the biological activity of a miRNA via direct experimental techniques. One effective and generally accepted alternative is to evaluate activity of a miRNA by assessing the expression of its targets. The computational determined activities make it possible to systematically evaluate the association between miRNA expression levels and biological activities.
Second, the p-values generated by the programs reflect the significance of miRNA activity changes. In fact, all the scoring systems of the four computational programs used in the manuscript have explicit definition and meaning. For instance, the mirAct web server infers the regulatory effect of a miRNA via a two-step procedure. First, a sample score measuring the activity of a miRNA in a sample is obtained by comparing the expression levels of its non-targets with those of targets. In the case of rank transformation of expression values, the difference of the average ranks between a miRNA's non-targets and targets is used. In the case of z-score representation of expression values, the two-sample t-statistic is applied. Then, the miRNA activity changes across different classes of samples are investigated by examining the sample scores via Kruskal-Wallis test, which tests the null hypothesis that all classes have identical miRNA activity.
Third, we highly appreciate your preciseness on the confounding factors that might influence data analysis. We completely agree your opinion that care must be taken when dealing with high-throughput data. That's why we performed the analysis using expression data from nine different tissues and four different computational programs. We think our results should be reliable. On one hand, as stated by the data providers that they had taken strategies to control the confounding factors that might influence their experiments. On the other hand, we repeated the analysis using different tissues and programs and obtained stable results.
Reviewer, Dr Yuriy Gusev, Georgetown University Medical Center, United States of America
This reviewer provided no comments for publication.
List of abbreviations
Acknowledgements and Funding
This work is supported by Chinese National Fundamental Research Project (Grants 2011CB910200), National Natural Science Foundation of China (30900270, 30821065), Knowledge Innovation Program of the Chinese Academy of Sciences (KSCX1-YW-02)]
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.PubMedView ArticleGoogle Scholar
- Sood P, Krek A, Zavolan M, Macino G, Rajewsky N: Cell-type-specific signatures of microRNA on target mRNA expression. Proc Natl Acad Sci USA. 2006, 103: 2746-2751. 10.1073/pnas.0511045103.PubMedPubMed CentralView ArticleGoogle Scholar
- Cheng C, Li LM: Inferring microRNA activities by combining gene expression with microRNA target prediction. PLoS One. 2008, 3: e1989-10.1371/journal.pone.0001989.PubMedPubMed CentralView ArticleGoogle Scholar
- van Dongen S, Abreu-Goodger C, Enright AJ: Detecting microRNA binding and siRNA off-target effects from expression data. Nature Methods. 2008, 5: 1023-1025. 10.1038/nmeth.1267.PubMedPubMed CentralView ArticleGoogle Scholar
- Liang Z, Zhou H, He Z, Zheng HR, Wu JR: mirAct: a web tool for evaluating microRNA activity based on gene expression data. Nucleic Acids Res. 2011.Google Scholar
- He H, Jazdzewski K, Li W, Liyanarachchi S, Nagy R, Volinia S, Calin G, Liu C, Franssila K, Suster S, Kloos RT, Croce CM, de la Chapelle A: The role of microRNA genes in papillary thyroid carcinoma. Proc Natl Acad Sci USA. 2005, 102: 19075-19080. 10.1073/pnas.0509603102.PubMedPubMed CentralView ArticleGoogle Scholar
- Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, Downing JR, Jacks T, Horvitz HR, Golub TR: MicroRNA expression profiles classify human cancers. Nature. 2005, 435: 834-838. 10.1038/nature03702.PubMedView ArticleGoogle Scholar
- Israel A, Sharan R, Ruppin E, Galun E: Increased microRNA activity in human cancers. PLoS One. 2009, 4 (6): e6045-10.1371/journal.pone.0006045.PubMedPubMed CentralView ArticleGoogle Scholar
- Cheng C, Fu X, Alves P, Gerstein M: mRNA expression profiles show differential regulatory effects of microRNAs between estrogen receptor-positive and estrogen receptor-negative breast cancer. Genome Biol. 2009, 10: R90-10.1186/gb-2009-10-9-r90.PubMedPubMed CentralView ArticleGoogle Scholar
- Fu X, Xue C, Huang Y, Xie Y, Li Y: The activity and expression of microRNAs in prostate cancers. Mol Biosyst. 2010, 6: 2561-2572. 10.1039/c0mb00100g.PubMedView ArticleGoogle Scholar
- Filipowicz W, Bhattacharyya SN, Sonenberg N: Mechanisms of post-transcirptional regulation by microRNAs: are the answers in sight. Nat Rev Genet. 2009, 9: 102-114.Google Scholar
- Kedde M, Strasser MJ, Boldajipour B, Oude Vrielink JA, Slanchev K, le Sage C, Nagel R, Voorhoeve PM, van Duijse J, Ørom UA, Lund AH, Perrakis A, Raz E, Agami R: RNA-binding protein Dnd1 inhibits microRNA access to target mRNA. Cell. 2007, 131: 1273-1286. 10.1016/j.cell.2007.11.034.PubMedView ArticleGoogle Scholar
- Geissmann TA, Touati D: Hfq, a new chaperoning role: binding to messenger RNA determines access for small RNA regulator. EMBO J. 2004, 23: 396-405. 10.1038/sj.emboj.7600058.PubMedPubMed CentralView ArticleGoogle Scholar
- Arvey A, Larsson E, Sander C, Leslie CS, Marks DS: Target mRNA abundance dilutes microRNA and siRNA activity. Mol Sys Biol. 2010, 6: 1-7.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.