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.

Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression

Feature selection followed by a novel residuals-based normalization that includes variance stabilization simplifies and improves single-cell gene expression analysis

Feature selection followed by a novel residuals-based normalization simplifies and improves single-cell gene expression analysis

Analytic Pearson residuals for normalization of single-cell RNA-seq UMI data

Non-linear Normalization for Non-UMI Single Cell RNA-Seq

Feature selection and dimension reduction for single-cell RNA-Seq based on a multinomial model

scKWARN: Kernel-weighted-average robust normalization for single-cell RNA-seq data

Normalization of Single-cell RNA-seq Data Using Partial Least Squares with Adaptive Fuzzy Weight

bayNorm: Bayesian gene expression recovery, imputation and normalization for single-cell RNA-sequencing data

Comparison and evaluation of statistical error models for scRNA-seq

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

Root Causal Inference from Single Cell RNA Sequencing with the Negative Binomial

Model-based dimensionality reduction for single-cell RNA-seq using generalized bilinear models

Negative binomial count splitting for single-cell RNA sequencing data

Normalizing need not be the norm: count-based math for analyzing single-cell data

BASiCS workflow: a step-by-step analysis of expression variability using single cell RNA sequencing data

Data normalization for addressing the challenges in the analysis of single-cell transcriptomic datasets

Compound models and Pearson residuals for single-cell RNA-seq data without UMIs

A statistical normalization method and differential expression analysis for RNA-seq data between different species

Scater: pre-processing, quality control, normalization and visualization of single-cell RNA-seq data in R

RCRnorm: An integrated system of random-coefficient hierarchical regression models for normalizing NanoString nCounter data