Abstract:Single-Cell RNA sequencing (scRNA-seq) measurements have facilitated genome-scale transcriptomic profiling of individual cells, with the hope of deconvolving cellular dynamic changes in corresponding cell sub-populations to better understand molecular mechanisms of different development processes. Several scRNA-seq analysis methods have been proposed to first identify cell sub-populations by clustering and then separately perform differential expression analysis to understand gene expression changes. Their corresponding statistical models and inference algorithms are often designed disjointly. We develop a new method -- SimCD -- that explicitly models cell heterogeneity and dynamic differential changes in one unified hierarchical gamma-negative binomial (hGNB) model, allowing simultaneous cell clustering and differential expression analysis for scRNA-seq data. Our method naturally defines cell heterogeneity by dynamic expression changes, which is expected to help achieve better performances on the two tasks compared to the existing methods that perform them separately. In addition, SimCD better models dropout (zero inflation) in scRNA-seq data by both cell- and gene-level factors and obviates the need for sophisticated pre-processing steps such as normalization, thanks to the direct modeling of scRNA-seq count data by the rigorous hGNB model with an efficient Gibbs sampling inference algorithm. Extensive comparisons with the state-of-the-art methods on both simulated and real-world scRNA-seq count data demonstrate the capability of SimCD to discover cell clusters and capture dynamic expression changes. Furthermore, SimCD helps identify several known genes affected by food deprivation in hypothalamic neuron cell subtypes as well as some new potential markers, suggesting the capability of SimCD for bio-marker discovery.

scDesign2: a transparent simulator that generates high-fidelity single-cell gene expression count data with gene correlations captured

scDesign3 generates realistic in silico data for multimodal single-cell and spatial omics

A unified framework of realistic in silico data generation and statistical model inference for single-cell and spatial omics

scPower accelerates and optimizes the design of multi-sample single cell transcriptomic studies

SimCH: simulation of single-cell RNA sequencing data by modeling cellular heterogeneity at gene expression level

Schi-Csim: a Flexible Simulator That Generates High-Fidelity Single-Cell Hi-C Data for Benchmarking

simCAS: an embedding-based method for simulating single-cell chromatin accessibility sequencing data

Systematic evaluation with practical guidelines for single-cell and spatially resolved transcriptomics data simulation under multiple scenarios

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

SCSsim: an Integrated Tool for Simulating Single-Cell Genome Sequencing Data.

scReadSim: a single-cell RNA-seq and ATAC-seq read simulator

CellSim: a novel software to calculate cell similarity and identify their co-regulation networks

Simulating multiple variability in spatially resolved transcriptomics with scCube

SCAN-ATAC-Sim: a scalable and efficient method for simulating single-cell ATAC-seq data from bulk-tissue experiments

Sequential design of single-cell experiments to identify discrete stochastic models for gene expression

Cytomulate: accurate and efficient simulation of CyTOF data

GLMsim: a GLM-based single cell RNA-seq simulator incorporating batch and biological effects

OneSC: A computational platform for recapitulating cell state transitions

SRTsim: spatial pattern preserving simulations for spatially resolved transcriptomics

Multi-task benchmarking of spatially resolved gene expression simulation models

TWO-SIGMA: A Novel Two-Component Single Cell Model-Based Association Method for Single-Cell RNA-seq Data