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.

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

TWO-SIGMA-G: a new competitive gene set testing framework for scRNA-seq data accounting for inter-gene and cell–cell correlation

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

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

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

scTSSR: gene expression recovery for single-cell RNA sequencing using two-side sparse self-representation

Sclm: Automatic Detection of Consensus Gene Clusters Across Multiple Single-Cell Datasets

Bayesian segmented Gaussian copula factor model for single-cell sequencing data

MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data

Accounting for technical noise in differential expression analysis of single-cell RNA sequencing data

The Poisson distribution model fits UMI-based single-cell RNA-sequencing data

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

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

Efficient and accurate mixed model association tool for single-cell eQTL analysis

SSNMDI: a novel joint learning model of semi-supervised non-negative matrix factorization and data imputation for clustering of single-cell RNA-seq data

Zero-inflated Smoothing Spline (ZISS) Models for Individual-level Single-cell Temporal Data

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

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

CoMM-S2: a Collaborative Mixed Model Using Summary Statistics in Transcriptome-Wide Association Studies

JOINT for large-scale single-cell RNA-sequencing analysis via soft-clustering and parallel computing

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