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.

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

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

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

scRecover: Discriminating true and false zeros in single-cell RNA-seq data for imputation

Sccode: an R Package for Data-Specific Differentially Expressed Gene Detection on Single-Cell RNA-sequencing Data.

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

Observation weights unlock bulk RNA-seq tools for zero inflation and single-cell applications

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

eSVD-DE: cohort-wide differential expression in single-cell RNA-seq data using exponential-family embeddings

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

DEGseq: an R package for identifying differentially expressed genes from RNA-seq data

A comprehensive assessment of hurdle and zero-inflated models for single cell RNA-sequencing analysis

Imputing dropouts for single-cell RNA sequencing based on multi-objective optimization

Statistical inference for time course RNA-Seq data using a negative binomial mixed-effect model

scHD4E: Novel ensemble learning-based differential expression analysis method for single-cell RNA-sequencing data

A new bioinformatics tool to recover missing gene expression in single-cell RNA sequencing data

Robustness of single-cell RNA-seq for identifying differentially expressed genes

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

UMI-count modeling and differential expression analysis for single-cell RNA sequencing

A mathematical framework for raw counts of single-cell RNA-seq data analysis

Imputing single-cell RNA-seq data by considering cell heterogeneity and prior expression of dropouts