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.

ZINBMM: a general mixture model for simultaneous clustering and gene selection using single-cell transcriptomic data

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

A sparse negative binomial mixture model for clustering RNA-seq count data

BREM-SC: A Bayesian Random Effects Mixture Model for Joint Clustering Single Cell Multi-omics Data

An interpretable Bayesian clustering approach with feature selection for analyzing spatially resolved transcriptomics data

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

Integrating ZINB-Based Autoencoder with Adaptive Data Augmentation Graph Contrastive Learning for scRNA-seq Clustering

Scalable nonparametric clustering with unified marker gene selection for single-cell RNA-seq data

Model-based Clustering of Zero-Inflated Single-Cell RNA Sequencing Data via the EM Algorithm

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

Identifying Cell Subpopulations and Their Genetic Drivers from Single-Cell RNA-Seq Data Using a Biclustering Approach

A Cell Marker-Based Clustering Strategy (cmcluster) for Precise Cell Type Identification of Scrna-Seq Data

Joint CC and Bimax: A Biclustering Method for Single-Cell RNA-Seq Data Analysis

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

Single-cell RNA-seq clustering: datasets, models, and algorithms

Single-Cell Transcriptome Data Clustering via Multinomial Modeling and Adaptive Fuzzy K-Means Algorithm

Analysis of Single-Cell RNA-seq Data by Clustering Approaches

Binomial models uncover biological variation during feature selection of droplet-based single-cell RNA sequencing

ZINB-Based Graph Embedding Autoencoder for Single-Cell RNA-Seq Interpretations

scASGC: An adaptive simplified graph convolution model for clustering single-cell RNA-seq data

Gene selection and clustering of single-cell data based on Fisher score and genetic algorithm