NEBULA is a fast negative binomial mixed model for differential or co-expression analysis of large-scale multi-subject single-cell data

The increasing availability of single-cell data revolutionizes the understanding of biological mechanisms at cellular resolution. For differential expression analysis in multi-subject single-cell data, negative binomial mixed models account for both subject-level and cell-level overdispersions, but are computationally demanding. Here, we propose an efficient NEgative Binomial mixed model Using a Large-sample Approximation (NEBULA). The speed gain is achieved by analytically solving high-dimensional integrals instead of using the Laplace approximation. We demonstrate that NEBULA is orders of magnitude faster than existing tools and controls false-positive errors in marker gene identification and co-expression analysis. Using NEBULA in Alzheimer's disease cohort data sets, we found that the cell-level expression of APOE correlated with that of other genetic risk factors (including CLU, CST3, TREM2, C1q, and ITM2B) in a cell-type-specific pattern and an isoform-dependent manner in microglia. NEBULA opens up a new avenue for the broad application of mixed models to large-scale multi-subject single-cell data.

1. 1
   Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

   Genome Biology201497,097 citationsDataRank 1.7
2. 2
   Fitting Linear Mixed-Effects Models Using<b>lme4</b>

   Journal of Statistical Software201582,566 citationsDataRank 1.7
3. 3
   Integrating single-cell transcriptomic data across different conditions, technologies, and species

   Nature Biotechnology201814,465 citationsDataRank 1.4
4. 4
   glmmTMB Balances Speed and Flexibility Among Packages for Zero-inflated Generalized Linear Mixed Modeling

   The R Journal201712,014 citationsDataRank 1.4
5. 5
   Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets

   Cell20157,753 citationsDataRank 1.3