Rxivist logo

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

By Christoph Hafemeister, Rahul Satija

Posted 14 Mar 2019
bioRxiv DOI: 10.1101/576827 (published DOI: 10.1186/s13059-019-1874-1)

Single-cell RNA-seq (scRNA-seq) data exhibits significant cell-to-cell variation due to technical factors, including the number of molecules detected in each cell, which can confound biological heterogeneity with technical effects. To address this, we present a modeling framework for the normalization and variance stabilization of molecular count data from scRNA-seq experiments. We propose that the Pearson residuals from 'regularized negative binomial regression', where cellular sequencing depth is utilized as a covariate in a generalized linear model, successfully remove the influence of technical characteristics from downstream analyses while preserving biological heterogeneity. Importantly, we show that an unconstrained negative binomial model may overfit scRNA-seq data, and overcome this by pooling information across genes with similar abundances to obtain stable parameter estimates. Our procedure omits the need for heuristic steps including pseudocount addition or log-transformation, and improves common downstream analytical tasks such as variable gene selection, dimensional reduction, and differential expression. Our approach can be applied to any UMI-based scRNA-seq dataset and is freely available as part of the R package sctransform (https://github.com/ChristophH/sctransform), with a direct interface to our single-cell toolkit Seurat.

Download data

  • Downloaded 21,782 times
  • Download rankings, all-time:
    • Site-wide: 308
    • In genomics: 12
  • Year to date:
    • Site-wide: 1,165
  • Since beginning of last month:
    • Site-wide: 1,857

Altmetric data

Downloads over time

Distribution of downloads per paper, site-wide


Sign up for the Rxivist weekly newsletter! (Click here for more details.)