bioRxivBioinformaticsPreprint
An assessment of normalization and differential expression methods for miRNA-seq analysis using a realistic benchmark dataset
Ernesto Aparicio‐PuertaiD, Andrea M BaraniD, John M. AshtoniD, Elizabeth M. PritchettiD, Anthony O. GacaiD, Jennifer Becker, Marc K. HalushkaiD, Seong-Hwan JuniD, Matthew N. McCalliD
Abstract
MicroRNAs are short noncoding RNAs that regulate gene expression and are commonly profiled by small RNA sequencing (miRNA-seq). Despite the widespread use of miRNA-seq, datasets are often analyzed with RNA-seq method such as DESeq2 or edgeR, which do not take into account the specific characteristics of miRNA-seq data. Here, we present a benchmark study of normalization and differential expression approaches using a realistic ground-truth dataset. By mixing mouse RNA of two organs, we generated expression trends while capturing biological and technical variability. Using monotonicity across the dataset and expected fold changes from the mixture design, we assessed normalization and differential expression methods. Normalization benchmarking showed that within-sample scaling, particularly Read Per Million (RPM), best preserved the expected monotonic trends, outperforming cross-sample methods such as TMM, rlog, and VST. These approaches sometimes recovered apparent monotonicity among abundant miRNAs, but inspection of individual profiles suggested likely over-correction. Regarding differential expression, edgeR consistently ranked among the best-performing methods across several metrics, including log2 fold-change estimation, with performance comparable to miRNA-seq-specific tools such as miRglmm and NBSR. DESeq2, edgeR-v4, and limma-based approaches tended to systematically underestimate log2 fold changes. Applying a common RPM-based normalization substantially improved the performance of cross-sample methods, highlighting the strong influence of normalization on differential expression analysis. Overall, our findings support within-sample scaling methods such as RPM for normalization, and edgeR, miRglmm, or NBSR for differential expression. The dataset has been made publicly available, providing a valuable resource for objective method comparison and future miRNA-seq software development.
- organism: mouse
- assay: bulk RNA-seq
- research focus: differential expression
- research focus: benchmarking
- approach: benchmark study
- Citations
- 0
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Improved differential expression analysis of miRNA-seq data by modeling competition to be counted
Seong-Hwan Jun, Marc K. Halushka, Matthew N. McCall·bioRxiv·May 7, 2024
- A Comparison of Methods: Normalizing High-Throughput RNA Sequencing Data
Rahul Reddy·bioRxiv·Sep 3, 2015
- Side-By-Side Analysis Of Alternative Approaches On Multi-Level RNA-Seq Data
Irina I Mohorianu·bioRxiv·Oct 6, 2018
- Cross-Platform Normalization Enables Machine Learning Model Training On Microarray And RNA-Seq Data Simultaneously
Steven M. Foltz, Casey S. Greene, Jaclyn Taroni·Communications BiologyPeer-reviewed·Feb 25, 2023
- An evaluation of RNA-seq differential analysis methods
Dongmei Li, Martin Stuart Zand, … Zidian Xie·PLoS ONEPeer-reviewed·Sep 16, 2022
- Evaluation of RNA-seq normalization methods using challenging datasets
Severin Uebbing·bioRxiv·Aug 27, 2018