Abstract:Integrated analysis of human gene expression data from multiple studies has become essential in genomics research for complex traits. However, integrating data generated from different cohorts with different platforms, such as microarray and RNA-seq, often requires data preprocessing, including normalization. In this study, we empirically evaluate 9 commonly used cross-platform normalization methods. We classify these methods into two main types: joint and separate normalization. Joint methods normalize multiple datasets together, while separate methods normalize each dataset independently. We further divide these methods into unsupervised and supervised approaches depending on whether they use outcomes during their normalization process. Examples of joint unsupervised methods include Quantile Normalization (QN), while Rank-in serves as an example of a joint supervised method. Training Distribution Matching (TDM) is an example of a separate unsupervised method. We assess each method's ability to cluster samples, predict outcomes, and detect differentially expressed (DE) genes using three real datasets and simulated data. First, our real data analysis suggests that while joint supervised methods can cluster sample groups better than the other two method groups, they double use the outcome data with artificially inflated clustering performance. Their biases are further demonstrated by their inflated type I error in DE analysis and clustering results from simulated data with no DE genes. Second, for outcome prediction, supervised normalization is no longer applicable. Among the unsupervised methods, QN significantly outperforms the other approaches, regardless of whether RNA-seq data is used to predict microarray outcomes or vice versa. Finally, we compare normalization methods on downstream DE analysis using simulation. In addition to direct DE analysis on the combined normalized data using the non-parametric Wilcoxon rank-sum test, we also perform a meta-analysis that combines p-values of the DE analysis from each individual data. For DE analysis, the meta-analysis consistently achieves the best balance between controlling type I error and maximizing power for DE gene detection. Our research suggests that while normalization is critical for the integrated analysis of transcriptomics data, simple QN is the most efficient and unbiased normalization approach for outcome prediction, and meta-analysis is the most appropriate for DE analysis.

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

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

RNAflow: An Effective and Simple RNA-Seq Differential Gene Expression Pipeline Using Nextflow

Improving the Diversity of Captured Full-Length Isoforms Using a Normalized Single-Molecule RNA-sequencing Method

Explicit Scale Simulation for analysis of RNA-sequencing with ALDEx2

Depth Normalization of Small RNA Sequencing: Using Data and Biology to Select a Suitable Method

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

A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis

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

A Unified Model for Joint Normalization and Differential Gene Expression Detection in RNA-Seq Data.

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

Normalization of RNA-Seq data using adaptive trimmed mean with multi-reference

Evaluating Cross-Platform Normalization Methods for Integrated Microarray and RNA-seq Data Analysis

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

Variation-preserving normalization unveils blind spots in gene expression profiling

Data normalization for addressing the challenges in the analysis of single-cell transcriptomic datasets

Replacing Normalizations with Interval Assumptions Improves the Rigor and Robustness of Differential Expression and Differential Abundance Analyses

Feature selection followed by a novel residuals-based normalization simplifies and improves single-cell gene expression analysis

Feature selection followed by a novel residuals-based normalization that includes variance stabilization simplifies and improves single-cell gene expression analysis

Non-linear Normalization for Non-UMI Single Cell RNA-Seq

Beyond Normalization: Incorporating Scale Uncertainty in Microbiome and Gene Expression Analysis