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.

The Parallel C++ Statistical Library for Bayesian Inference: QUESO

hIPPYlib-MUQ: A Bayesian Inference Software Framework for Integration of Data with Complex Predictive Models under Uncertainty

CUQIpy - Part I: computational uncertainty quantification for inverse problems in Python

Bayesian probabilistic propagation of hybrid uncertainties: Estimation of response expectation function, its variable importance and bounds

Scalable Uncertainty Quantification for Deep Operator Networks using Randomized Priors

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

CUQIpy - Part II: computational uncertainty quantification for PDE-based inverse problems in Python

NeuralUQ: A Comprehensive Library for Uncertainty Quantification in Neural Differential Equations and Operators

A Direct Coupled Solution Methodology for Efficient Robust Optimizations of Inverse Problems under Uncertainty

CUAOA: A Novel CUDA-Accelerated Simulation Framework for the QAOA

UQSA -- An R-Package for Uncertainty Quantification and Sensitivity Analysis for Biochemical Reaction Network Models

SOBER: Highly Parallel Bayesian Optimization and Bayesian Quadrature over Discrete and Mixed Spaces

Uncertainty Quantification of the 4th kind; optimal posterior accuracy-uncertainty tradeoff with the minimum enclosing ball

On Uncertainty Quantification for Near-Bayes Optimal Algorithms

Collaborative and Adaptive Bayesian Optimization for bounding variances and probabilities under hybrid uncertainties

Neural Message Passing for Objective-Based Uncertainty Quantification and Optimal Experimental Design.

A Quadrature Approach for General-Purpose Batch Bayesian Optimization via Probabilistic Lifting

Scalable Bayesian uncertainty quantification with data-driven priors for radio interferometric imaging

Uncertainty quantification for deep learning-based schemes for solving high-dimensional backward stochastic differential equations

Multilevel Sequential${}^2$ Monte Carlo for Bayesian Inverse Problems

A Nonparametric Bayesian Framework for Uncertainty Quantification in Stochastic Simulation