Abstract
MicroRNAs (miRNAs) are short noncoding RNAs that are involved in post-transcriptional regulation of mRNAs. Microarrays have been employed to measure global miRNA expressions; however, because the number of miRNAs is much smaller than the number of mRNAs, it is not clear whether traditional normalization methods developed for mRNA arrays are suitable for miRNA. This is an important question, since normalization affects downstream analyses of the data. In this paper we develop a least-variant set (LVS) normalization method, which was previously shown to outperform other methods in mRNA analysis when standard assumptions are violated. The selection of the LVS miRNAs is based on a robust linear model fit of the probe-level data that takes into account the considerable differences in variances between probes. In a spike-in study, we show that the LVS has similar operating characteristics, in terms of sensitivity and specificity, compared with the ideal normalization, and it is better than no normalization, 75th percentile-shift, quantile, global median, VSN, and lowess normalization methods. We evaluate four expression-summary measures using a tissue data set; summarization from the robust model performs as well as the others. Finally, comparisons using expression data from two dissimilar tissues and two similar ones show that LVS normalization has better operating characteristics than other normalizations.