Abstract:Probabilistic models have provided the underpinnings for state-of-the-art performance in many single-cell omics data analysis tasks, including dimensionality reduction, clustering, differential expression, annotation, removal of unwanted variation, and integration across modalities. Many of the models being deployed are amenable to scalable stochastic inference techniques, and accordingly they are able to process single-cell datasets of realistic and growing sizes. However, the community-wide adoption of probabilistic approaches is hindered by a fractured software ecosystem resulting in an array of packages with distinct, and often complex interfaces. To address this issue, we developed scvi-tools (), a Python package that implements a variety of leading probabilistic methods. These methods, which cover many fundamental analysis tasks, are accessible through a standardized, easy-to-use interface with direct links to Scanpy, Seurat, and Bioconductor workflows. By standardizing the implementations, we were able to develop and reuse novel functionalities across different models, such as support for complex study designs through nonlinear removal of unwanted variation due to multiple covariates and reference-query integration via scArches. The extensible software building blocks that underlie scvi-tools also enable a developer environment in which new probabilistic models for single cell omics can be efficiently developed, benchmarked, and deployed. We demonstrate this through a code-efficient reimplementation of Stereoscope for deconvolution of spatial transcriptomics profiles. By catering to both the end user and developer audiences, we expect scvi-tools to become an essential software dependency and serve to formulate a community standard for probabilistic modeling of single cell omics. ### Competing Interest Statement O.C. is supported by the EPSRC Centre for Doctoral Training in Modern Statistics and Statistical Machine Learning (EP/S023151/1) and Novo Nordisk. V.S. is a full-time employee of Serqet Therapuetics and has ownership interest in Serqet Therapeutics. F.J.T. reports receiving consulting fees from Roche Diagnostics GmbH and Cellarity Inc., and ownership interest in Cellarity, Inc.

A unified framework of realistic in silico data generation and statistical model inference for single-cell and spatial omics

scDesign3 generates realistic in silico data for multimodal single-cell and spatial omics

scDesign2: a transparent simulator that generates high-fidelity single-cell gene expression count data with gene correlations captured

Scvi-Tools: a Library for Deep Probabilistic Analysis of Single-Cell Omics Data

Systematic evaluation with practical guidelines for single-cell and spatially resolved transcriptomics data simulation under multiple scenarios

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

Cell4D: a general purpose spatial stochastic simulator for cellular pathways

Simulating multiple variability in spatially resolved transcriptomics with scCube

Multi-task benchmarking of spatially resolved gene expression simulation models

Schi-Csim: a Flexible Simulator That Generates High-Fidelity Single-Cell Hi-C Data for Benchmarking

A Unified Probabilistic Framework for Modeling and Inferring Spatial Transcriptomic Data

High-throughput Single-Cell RNA-seq Data Imputation and Characterization with Surrogate-Assisted Automated Deep Learning

A Unified Statistical Framework for Single Cell and Bulk RNA Sequencing Data

Comprehensive single-cell RNA-seq analysis using deep interpretable generative modeling guided by biological hierarchy knowledge

scINSIGHT for interpreting single-cell gene expression from biologically heterogeneous data

Comparing AI versus Optimization Workflows for Simulation-Based Inference of Spatial-Stochastic Systems

Scminer: a Mutual Information-Based Framework for Identifying Hidden Drivers from Single-Cell Omics Data.

A novel framework for single-cell Hi-C clustering based on graph-convolution-based imputation and two-phase-based feature extraction

spaDesign: A Statistical Framework to Improve the Design of Sequencing-based Spatial Transcriptomics Experiments

Sequential design of single-cell experiments to identify discrete stochastic models for gene expression

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