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.

Gene count normalization in single-cell imaging-based spatially resolved transcriptomics

Gene count normalization in single-cell imaging-based spatially resolved transcriptomics

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

Cgcorrect: a Method to Correct for Confounding Cell-Cell Variation Due to Cell Growth in Single-Cell Transcriptomics

Use of normalization methods for analysis of microarrays containing a high degree of gene effects

Transcriptional output, cell-type densities, and normalization in spatial transcriptomics

Library size normalization affects spatial domain identification

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

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

Evaluation of Normalization Methods for RNA-Seq Gene Expression Estimation

A Systematic Comparison of Normalization Methods for Eqtl Analysis

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

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

Normalization of Real-Time Quantitative Reverse Transcription-PCR Data: A Model-Based Variance Estimation Approach to Identify Genes Suited for Normalization, Applied to Bladder and Colon Cancer Data Sets

Variation-preserving normalization unveils blind spots in gene expression profiling

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

Library size confounds biology in spatial transcriptomics data

Bayesian inference of gene expression states from single-cell RNA-seq data

A multi-platform normalization method for meta-analysis of gene expression data

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