Abstract:BackgroundThe analysis of single-cell RNA sequencing (scRNAseq) data plays an important role in understanding the intrinsic and extrinsic cellular processes in biological and biomedical research. One significant effort in this area is the detection of differentially expressed (DE) genes. scRNAseq data, however, are highly heterogeneous and have a large number of zero counts, which introduces challenges in detecting DE genes. Addressing these challenges requires employing new approaches beyond the conventional ones, which are based on a nonzero difference in average expression. Several methods have been developed for differential gene expression analysis of scRNAseq data. To provide guidance on choosing an appropriate tool or developing a new one, it is necessary to evaluate and compare the performance of differential gene expression analysis methods for scRNAseq data.ResultsIn this study, we conducted a comprehensive evaluation of the performance of eleven differential gene expression analysis software tools, which are designed for scRNAseq data or can be applied to them. We used simulated and real data to evaluate the accuracy and precision of detection. Using simulated data, we investigated the effect of sample size on the detection accuracy of the tools. Using real data, we examined the agreement among the tools in identifying DE genes, the run time of the tools, and the biological relevance of the detected DE genes.ConclusionsIn general, agreement among the tools in calling DE genes is not high. There is a trade-off between true-positive rates and the precision of calling DE genes. Methods with higher true positive rates tend to show low precision due to their introducing false positives, whereas methods with high precision show low true positive rates due to identifying few DE genes. We observed that current methods designed for scRNAseq data do not tend to show better performance compared to methods designed for bulk RNAseq data. Data multimodality and abundance of zero read counts are the main characteristics of scRNAseq data, which play important roles in the performance of differential gene expression analysis methods and need to be considered in terms of the development of new methods.

A Two-Part Mixed Model for Differential Expression Analysis in Single-Cell High-Throughput Gene Expression Data.

A two-step strategy for detecting differential gene expression in cDNA microarray data

Differential expression analysis for paired RNA-seq data

Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation

A robust model for cell type-specific interindividual variation in single-cell RNA sequencing data

Assessing transcriptomic heterogeneity of single-cell RNASeq data by bulk-level gene expression data

Differential Expression Analyses for Single-Cell RNA-Seq: Old Questions on New Data

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

Accounting for technical noise in differential expression analysis of single-cell RNA sequencing data

Comparative analysis of differential gene expression analysis tools for single-cell RNA sequencing data

Unit-Free and Robust Detection of Differential Expression from RNA-Seq Data

Modeling expression ranks for noise-tolerant differential expression analysis of scRNA-seq data

Data Exploration, Quality Control and Testing in Single-Cell qPCR-Based Gene Expression Experiments

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

Bias, robustness and scalability in differential expression analysis of single-cell RNA-seq data

DEsingle for Detecting Three Types of Differential Expression in Single-Cell RNA-seq Data

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

A flexible model for correlated count data, with application to multi-condition differential expression analyses of single-cell RNA sequencing data

Network Module-Based Model in the Differential Expression Analysis for RNA-seq.

SimCD: Simultaneous Clustering and Differential expression analysis for single-cell transcriptomic data

A Unified Model for Differential Expression Analysis of RNA-seq Data Via L1-Penalized Linear Regression