Abstract:Accounting for batch effects, especially latent batch effects, in differential expression (DE) analysis is critical for identifying true biological effects. Single-cell RNA sequencing (scRNA-seq) is a powerful tool for quantifying cell-to-cell variation in transcript abundance and characterizing cellular dynamics. Although many scRNA-seq DE analysis methods accommodate known batch variables, their performance has not been systematically evaluated. Moreover, the challenge of accounting for latent batch variables in scRNA-seq DE analysis is largely unmet. In contrast, many methods have been developed to account for batch variables (either known or latent) in other high-dimensional data, especially bulk RNA-seq. We extensively evaluate eleven methods for batch variables in different scRNA-seq DE analysis scenarios, with a primary focus on latent batch variables. We demonstrate that for known batch variables, incorporating them as covariates into a regression model outperformed approaches using batch-corrected matrix. For latent batches, fixed effects models have inflated FDRs, whereas aggregation-based methods and mixed effects models have significant power loss. Surrogate variable based methods generally control the FDR well while achieving good power with small group effects. However, their performance (except SVA) deteriorated substantially in scenarios involving large group effects and/or group label impurity. In these settings, SVA achieves relatively good performance despite occasionally inflated FDR (up to 0.2). Finally we make following recommendations for scRNA-seq DE analysis: 1) incorporating known batch variables instead of using batch-corrected data; 2) employing SVA for latent batch correction and 3) better methods are still needed to fully unleash the power of scRNA-seq.<ol class="links-for-figure"><li><a class="anchor download-link u-font-sans" href="https://ars.els-cdn.com/content/image/1-s2.0-S200103701930409X-ga1_lrg.jpg">Download : Download high-res image (88KB)</a></li><li><a class="anchor download-link u-font-sans" href="https://ars.els-cdn.com/content/image/1-s2.0-S200103701930409X-ga1.jpg">Download : Download full-size image</a></li></ol>

Trade-off between conservation of biological variation and batch effect removal in deep generative modeling for single-cell transcriptomics

Batch alignment of single-cell transcriptomics data using deep metric learning

Benchmarking deep learning methods for biologically conserved single-cell integration

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

Probabilistic harmonization and annotation of single‐cell transcriptomics data with deep generative models

HDMC: a novel deep learning-based framework for removing batch effects in single-cell RNA-seq data

A comparison of methods accounting for batch effects in differential expression analysis of UMI count based single cell RNA sequencing

A deep generative model for single-cell RNA sequencing with application to detecting differentially expressed genes

Mixed Effects Deep Learning for the interpretable analysis of single cell RNA sequencing data by quantifying and visualizing batch effects

Parameterizing Cell-to-cell Regulatory Heterogeneities Via Stochastic Transcriptional Profiles

A joint deep learning model enables simultaneous batch effect correction, denoising, and clustering in single-cell transcriptomics

Synthetic single cell RNA sequencing data from small pilot studies using deep generative models

Integrating single-cell RNA-seq datasets with substantial batch effects

Single-Cell RNA-Seq Debiased Clustering Via Batch Effect Disentanglement.

Deep feature extraction of single-cell transcriptomes by generative adversarial network

Dropout Imputation and Batch Effect Correction for Single-Cell RNA Sequencing Data

Deep Generative Autoencoder for Low-Dimensional Embeding Extraction from Single-Cell RNAseq Data

Simultaneous Deep Generative Modelling and Clustering of Single-Cell Genomic Data

A deep generative model for gene expression profiles from single-cell RNA sequencing

Model-based dimensionality reduction for single-cell RNA-seq using generalized bilinear models

Binomial models uncover biological variation during feature selection of droplet-based single-cell RNA sequencing