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.

Scywalker: scalable end-to-end data analysis workflow for nanopore single-cell transcriptome sequencing

Scywalker: Scalable end-to-end data analysis workflow for long-read single-cell transcriptome sequencing

Scanpy for analysis of large-scale single-cell gene expression data

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

SCALA: A complete solution for multimodal analysis of single-cell Next Generation Sequencing data

Interactive analysis of single-cell data using flexible workflows with SCTK2

Accelerating Single-Cell Sequencing Data Analysis with SciDAP: A User-Friendly Approach

Scater: pre-processing, quality control, normalization and visualization of single-cell RNA-seq data in R

scPower accelerates and optimizes the design of multi-sample single cell transcriptomic studies

scExplorer: A Comprehensive Web Server for Single-Cell RNA Sequencing Data Analysis

SingleCAnalyzer: Interactive Analysis of Single Cell RNA-Seq Data on the Cloud

scRNASequest: an ecosystem of scRNA-seq analysis, visualization, and publishing

Single-Cell Omics for Transcriptome CHaracterization (SCOTCH): isoform-level characterization of gene expression through long-read single-cell RNA sequencing

WASP: a versatile, web-accessible single cell RNA-Seq processing platform

Scdrake: a reproducible and scalable pipeline for scRNA-seq data analysis

A comprehensive workflow for allele-specific immune gene quantification and expression analysis in single-cell RNA-seq data

Accurate long-read transcript discovery and quantification at single-cell, pseudo-bulk and bulk resolution with Isosceles

sc-SPLASH provides ultra-efficient reference-free discovery in barcoded single-cell sequencing

AsaruSim: a single-cell and spatial RNA-Seq Nanopore long-reads simulation workflow

A flexible cross-platform single-cell data processing pipeline

Sainsc: a computational tool for segmentation-free analysis of in-situ capture