Abstract:Background Recently, rapid improvements in technology and decrease in sequencing costs have made RNA-Seq a widely used technique to quantify gene expression levels. Various normalization approaches have been proposed, owing to the importance of normalization in the analysis of RNA-Seq data. A comparison of recently proposed normalization methods is required to generate suitable guidelines for the selection of the most appropriate approach for future experiments. Results In this paper, we compared eight non-abundance (RC, UQ, Med, TMM, DESeq, Q, RPKM, and ERPKM) and two abundance estimation normalization methods (RSEM and Sailfish). The experiments were based on real Illumina high-throughput RNA-Seq of 35- and 76-nucleotide sequences produced in the MAQC project and simulation reads. Reads were mapped with human genome obtained from UCSC Genome Browser Database. For precise evaluation, we investigated Spearman correlation between the normalization results from RNA-Seq and MAQC qRT-PCR values for 996 genes. Based on this work, we showed that out of the eight non-abundance estimation normalization methods, RC, UQ, Med, TMM, DESeq, and Q gave similar normalization results for all data sets. For RNA-Seq of a 35-nucleotide sequence, RPKM showed the highest correlation results, but for RNA-Seq of a 76-nucleotide sequence, least correlation was observed than the other methods. ERPKM did not improve results than RPKM. Between two abundance estimation normalization methods, for RNA-Seq of a 35-nucleotide sequence, higher correlation was obtained with Sailfish than that with RSEM, which was better than without using abundance estimation methods. However, for RNA-Seq of a 76-nucleotide sequence, the results achieved by RSEM were similar to without applying abundance estimation methods, and were much better than with Sailfish. Furthermore, we found that adding a poly-A tail increased alignment numbers, but did not improve normalization results. Conclusion Spearman correlation analysis revealed that RC, UQ, Med, TMM, DESeq, and Q did not noticeably improve gene expression normalization, regardless of read length. Other normalization methods were more efficient when alignment accuracy was low; Sailfish with RPKM gave the best normalization results. When alignment accuracy was high, RC was sufficient for gene expression calculation. And we suggest ignoring poly-A tail during differential gene expression analysis.

Compositional Data Analysis is necessary for simulating and analyzing RNA-Seq data

Unifying the analysis of high-throughput sequencing datasets: characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis

Analysis and correction of compositional bias in sparse sequencing count data

Less is More: Relative Rank is More Informative Than Absolute Abundance for Compositional NGS Data.

Comparing the Normalization Methods for the Differential Analysis of Illumina High-Throughput RNA-Seq Data

Count-based differential expression analysis of RNA sequencing data using R and Bioconductor

A scaling normalization method for differential expression analysis of RNA-seq data

Degps is a Powerful Tool for Detecting Differential Expression in RNA-sequencing Studies

Normalization, testing, and false discovery rate estimation for RNA-sequencing data

Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences

Finding consistent patterns: a nonparametric approach for identifying differential expression in RNA-Seq data

Finding the centre: corrections for asymmetry in high-throughput sequencing datasets

A Powerful and Flexible Approach to the Analysis of RNA Sequence Count Data

Normalizing need not be the norm: count-based math for analyzing single-cell data

Evaluation of Normalization Methods for RNA-Seq Gene Expression Estimation

Comprehensive evaluation of differential expression analysis methods for RNA-seq data

Comparison of Scrna-Seq Data Analysis Method Combinations.

Statistical Modeling of RNA-Seq Data

Robust Differential Abundance Test in Compositional Data

RNA‐seq Data: Challenges in and Recommendations for Experimental Design and Analysis

A statistical normalization method and differential expression analysis for RNA-seq data between different species