Abstract:We provide additional plots and tables for the results section in Supplementary Information. We have developed an R package “RZiMM”, to implement our proposed method, and scripts to perform simulational and real applications. Applications of single-cell RNA sequencing in various biomedical research areas have been blooming. This new technology provides unprecedented opportunities to study disease heterogeneity at the cellular level. However, unique characteristics of scRNA-seq data, including large dimensionality, high dropout rates, and possibly batch effects, bring great difficulty into the analysis of such data. Not appropriately addressing these issues obstructs true scientific discovery. Herein we propose a unified Regularized Zero-inflated Mixture Model framework, designed for scRNA-seq data (RZiMM-scRNA), to simultaneously detect cell subgroups and identify gene differential expression based on a developed importance score, accounting for both dropouts and batch effects. We conduct extensive simulation studies in which we evaluate the performance of RZiMM-scRNA and compare it with several popular methods, including Seurat, SC3, K-means, and hierarchical clustering. Simulation results show that RZiMM-scRNA demonstrates superior clustering performance and enhanced biomarker detection accuracy, compared to alternative methods, especially when cell subgroups are less distinct, verifying the robustness of our method. Our empirical investigations focus on two brain tumor studies dealing with astrocytoma of various grades, including the most malignant of all brain tumors, glioblastoma multiforme (GBM). Our goal is to delineate cell heterogeneity and identify driving biomarkers associated with these tumors. Notably, RZiMM-scNRA successfully identifies a small group of oligodendrocyte cells, which has drawn much attention in biomedical literature on brain cancers. In addition, our method discovers several new biomarkers which are not discussed in the original studies, including PLP1, BCAN, and PTPRZ1—all associated with the development and malignant growth of glioma—as well as CAMK2B, which is downregulated in glioma and GBM and implicated in neurodevelopment, brain function, learning and memory processes.

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

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

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

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

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

Cellular Similarity based Imputation for Single cell RNA Sequencing Data

Scimc: a Platform for Benchmarking Comparison and Visualization Analysis of Scrna-Seq Data Imputation Methods.

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

The shaky foundations of simulating single-cell RNA sequencing data

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

SimuSCoP: Reliably Simulate Illumina Sequencing Data Based on Position and Context Dependent Profiles

Single-Cell Transcriptome Profiling Simulation Reveals the Impact of Sequencing Parameters and Algorithms on Clustering

DIMM-SC: a Dirichlet mixture model for clustering droplet-based single cell transcriptomic data

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

scCGImpute: An Imputation Method for Single-Cell RNA Sequencing Data Based on Similarities between Cells and Relationships among Genes

A feasible roadmap to identifying significant intercellular genomic heterogeneity in deep sequencing data

Machine learning and statistical methods for clustering single-cell RNA-sequencing data

RZiMM-scRNA: A regularized zero-inflated mixture model framework for single-cell RNA-seq data

scAMACE: Model-based approach to the joint analysis of single-cell data on chromatin accessibility, gene expression and methylation

scRNA-seq mixology: towards better benchmarking of single cell RNA-seq analysis methods

A Network Propagation Based Approach for Measuring Cell-Cell Similarity.