Reproducibility of Methods to Detect Differentially Expressed Genes from Single-Cell RNA Sequencing

Abstract

Detection of differentially expressed genes is a common task in single-cell RNA-seq (scRNA-seq) studies. Various methods based on both bulk-cell and single-cell approaches are in current use. Due to the unique distributional characteristics of singlecell data, it is important to compare these methods with rigorous statistical assessments. In this study, we assess the reproducibility of 9 tools for differential expression analysis in scRNA-seq data. These tools include four methods originally designed for scRNA-seq data, three popular methods originally developed for bulk-cell RNA-seq data but have been applied in scRNA-seq analysis, and two general statistical tests. Instead of comparing the performance across all genes, we compare the methods in terms of the rediscovery rates (RDRs) of top-ranked genes, separately for highly and lowly expressed genes. Three real and one simulated scRNA-seq data sets are used for the comparisons. The results indicate that some widely used methods, such as edgeR and monocle, have worse RDR performances compared to the other methods, especially for the top-ranked genes. For highly expressed genes, many bulk-cell-based methods can perform similarly to the methods designed for scRNA-seq data. But for the lowly expressed genes performance varies substantially; edgeR and monocle are too liberal and have poor control of false positives, while DESeq2 is too conservative and consequently loses sensitivity compared to the other methods. BPSC, Limma, DEsingle, MAST, t-test and Wilcoxon have similar performances in the real data sets. Overall, the scRNA-seq based method BPSC performs well against the other methods, particularly when there is a sufficient number of cells.

Dysregulation of brain and choroid plexus cell types in severe COVID-19

“…W model, and has demonstrated superior type-I error control without significantly sacrificing sensitivity [55][56][57][58] . First, we ensured that our data did not exhibit signs of confounding effects (Extended Data Fig.…”

Myonuclear transcriptional dynamics in response to exercise following satellite cell depletion

“…Differential gene expression analyses among datasets and among subsets within datasets were performed using the limma (v3.44.3) R package (Ritchie et al, 2015). We chose this package because previous studies report consistent performance in terms of sensitivity, specificity, and rediscovery rates for limma compared to numerous other available packages (Mou et al, 2019;Soneson and Robinson, 2018). Read counts were extracted from Seurat objects using the ''GetAssayData'' function and lowly expressed genes (total expression across all nuclei < 50 or having non-zero expression in fewer than 10 nuclei) were removed from analysis improve specificity because previous studies have found that too many lowly expressed genes contribute to increased false discovery (Soneson and Robinson, 2018).…”

Methods for nonparametric statistics in scientific research. Overview. Part 1.

“…Comparison of 441 daily smokers with 441 non-smokers was performed using the of Wilcoxon T-test and t-test. The study [47] evaluated the reproducibility of 9 instruments for the analysis of differential expression in scRNA-seq data. Statistical analysis was performed using the t-test and Wilcoxon T-test.…”