Method of moments approach for generalized differential gene expression analysis (Perturb-seq)

Differential expression analysis of scRNA-seq data is central for characterizing how experimental factors affect the distribution of gene expression. However, distinguishing between biological and technical sources of cell-cell variability and assessing the statistical significance of quantitative comparisons between cell groups remains challenging. Here, we introduce memento, a tool designed to address these limitations by providing statistically robust and computationally efficient differential expression analysis of the mean, variability, and gene correlation from scRNA-seq. We applied memento to analyze 70,000 tracheal epithelial cells to identify interferon response genes with distinct variability and correlation patterns, 160,000 T cells perturbed with CRISPR-Cas9 to reconstruct gene-regulatory networks that control T cell activation, 1.2 million PMBCs to map cell-type-specific cis expression quantitative trait loci (eQTLs), and arbitrary cell groups within the entire 50 million cell CELLxGENE Discover data corpus. In all cases, memento identified more significant and reproducible differences in mean expression but also identified differences in variability and gene correlation that suggest distinct transcriptional regulation mechanisms imparted by cytokines, genetic perturbations, and natural genetic variation. These results demonstrate that memento is a first-in-class method for the quantitative analysis of scRNA-seq data, scalable to millions of cells and thousands of samples. Overall design: CD4+ T cells were isolated Primary human CD4+ T cells were isolated from peripheral blood mononuclear cells, edited using SLICE protocol and analyzed by scRNA-seq. Human tracheal epithelial cells were harvested from deceased organ donors according to established protocols (PMID: 1616056), stimulated by 4 different interferons followed by scRNA-seq.