Abstract:Finding the right balance of quality and quantity can be important, and it is essential that project quality does not drop below the level where important main conclusions are missed or misstated. We use knock-out and over-expression studies as a simplification to test recovery of a known causal gene in RNA-Seq cell line experiments. When single-end RNA-Seq reads are aligned with STAR and quantified with htseq-count, we found potential value in testing the use of the Generalized Linear Model (GLM) implementation of edgeR with robust dispersion estimation more frequently for either single-variate or multi-variate 2-group comparisons (with the possibility of defining criteria less stringent than |fold-change| > 1.5 and FDR < 0.05). When considering a limited number of patient sample comparisons with larger sample size, there might be some decreased variability between methods (except for DESeq1). However, at the same time, the ranking of the gene identified using immunohistochemistry (for ER/PR/HER2 in breast cancer samples from The Cancer Genome Atlas) showed as possible shift in performance compared to the cell line comparisons, potentially highlighting utility for standard statistical tests and/or limma-based analysis with larger sample sizes. If this continues to be true in additional studies and comparisons, then that could be consistent with the possibility that it may be important to allocate time for potential methods troubleshooting for genomics projects. Analysis of public data presented in this study does not consider all experimental designs, and presentation of downstream analysis is limited. So, any estimate from this simplification would be an underestimation of the true need for some methods testing for every project. Additionally, this set of independent cell line experiments has a limitation in being able to determine the frequency of missing a highly important gene if the problem is rare (such as 10% or lower). For example, if there was an assumption that only one method can be tested for “initial” analysis, then it is not completely clear to the extent that using edgeR-robust might perform better than DESeq2 in the cell line experiments. Importantly, we do not wish to cause undue concern, and we believe that it should often be possible to define a gene expression differential expression workflow that is suitable for some purposes for many samples. Nevertheless, at the same time, we provide a variety of measures that we believe emphasize the need to critically assess every individual project and maximize confidence in published results.

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

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

scDEA: differential expression analysis in single-cell RNA-sequencing data via ensemble learning

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

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

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

A multitask clustering approach for single-cell RNA-seq analysis in Recessive Dystrophic Epidermolysis Bullosa

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

Identification of cell types, states and programs by learning gene set representations

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

Critical Differential Expression Assessment for Individual Bulk RNA-Seq Projects

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

Leveraging cross-source heterogeneity to improve the performance of bulk gene expression deconvolution

SCDC: Bulk Gene Expression Deconvolution by Multiple Single-Cell RNA Sequencing References

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

Comparative study on differential expression analysis methods for single-cell RNA sequencing data with small biological replicates: Based on single-cell transcriptional data of PBMCs from COVID-19 severe patients

Comprehensive evaluation of differential expression analysis methods for RNA-seq data

Demultiplexing of single-cell RNA-sequencing data using interindividual variation in gene expression

SDEvelo: a deep generative approach for transcriptional dynamics with cell-specific latent time and multivariate stochastic modeling

Detection of differentially abundant cell subpopulations in scRNA-seq data

SpatialDE2: Fast and localized variance component analysis of spatial transcriptomics