Abstract:Large-scale comprehensive single-cell experiments are often resource-intensive and require the involvement of many laboratories and/or taking measurements at various times. This inevitably leads to batch effects, and systematic variations in the data that might occur due to different technology platforms, reagent lots, or handling personnel. Such technical differences confound biological variations of interest and need to be corrected during the data integration process. Data integration is a challenging task due to the overlapping of biological and technical factors, which makes it difficult to distinguish their individual contribution to the overall observed effect. Moreover, the choice of integration method may impact the downstream analyses, including searching for differentially expressed genes. From the existing data integration methods, we selected only those that return the full expression matrix. We evaluated six methods in terms of their influence on the performance of differential gene expression analysis in two single-cell datasets with the same biological study design that differ only in the way the measurement was done: one dataset manifests strong batch effects due to the measurements of each sample at a different time. Integrated data were visualized using the UMAP method. The evaluation was done both on individual gene level using parametric and non-parametric approaches for finding differentially expressed genes and on gene set level using gene set enrichment analysis. As an evaluation metric, we used two correlation coefficients, Pearson and Spearman, of the obtained test statistics between reference, test, and corrected studies. Visual comparison of UMAP plots highlighted ComBat-seq, limma, and MNN, which reduced batch effects and preserved differences between biological conditions. Most of the tested methods changed the data distribution after integration, which negatively impacts the use of parametric methods for the analysis. Two algorithms, MNN and Scanorama, gave very poor results in terms of differential analysis on gene and gene set levels. Finally, we highlight ComBat-seq as it led to the highest correlation of test statistics between reference and corrected dataset among others. Moreover, it does not distort the original distribution of gene expression data, so it can be used in all types of downstream analyses.

The effect of data transformation on low-dimensional integration of single-cell RNA-seq

An Introduction to the Analysis of Single-Cell RNA-Sequencing Data

Single Cells Make Big Data: New Challenges and Opportunities in Transcriptomics

Integration of Single-Cell RNA-Seq Datasets: A Review of Computational Methods

A comparison of deep learning-based pre-processing and clustering approaches for single-cell RNA sequencing data

A Comparison for Dimensionality Reduction Methods of Single-Cell RNA-seq Data

Influence of single-cell RNA sequencing data integration on the performance of differential gene expression analysis

Statistical Single Cell Multi-Omics Integration

Integration for single-cell RNA sequencing data based on the shared cell type assignment

How to Get Started with Single Cell RNA Sequencing Data Analysis

Practical bioinformatics pipelines for single-cell RNA-seq data analysis

The effect of background noise and its removal on the analysis of single-cell expression data

Optimization and Redevelopment of Single-Cell Data Analysis Workflow Based on Deep Generative Models

Processing pipelines and analytical methods for single-cell DNA methylation sequencing data

Critical downstream analysis steps for single-cell RNA sequencing data

The impact of package selection and versioning on single-cell RNA-seq analysis

A Review of Single-Cell RNA-Seq Annotation, Integration, and Cell-Cell Communication

Boosting the Power of Single-Cell Analysis

Single-Cell RNA Sequencing Analysis: A Step-by-Step Overview

Current best practices in single‐cell RNA‐seq analysis: a tutorial

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