ScANVI is a variational inference model for single-cell RNA-seq data that can learn an underlying
latent space, integrate technical batches and impute dropouts.
In addition, to scVI, ScANVI is a semi-supervised model that can leverage labeled data to learn a
cell-type classifier in the latent space and afterward predict cell types of new data.
The learned low-dimensional latent representation of the data can be used for visualization and
clustering.
scANVI takes as input a scRNA-seq gene expression matrix with cells and genes as well as a
cell-type annotation for a subset of cells.
We provide an extensive
user guide.
- See our original manuscript for further details of the model:
scANVI manuscript.
- See our manuscript on scvi-hub
how to leverage pre-trained models.
This model can be used for fine tuning on new data using our Arches framework:
Arches tutorial.
Model Description
Tabula Sapiens is a benchmark, first-draft human cell atlas of nearly 500,000 cells from 24 organs of 15 normal human subjects.
Metrics
We provide here key performance metrics for the uploaded model, if provided by the data uploader.
Coefficient of variation
The cell-wise coefficient of variation summarizes how well variation between different cells is
preserved by the generated model expression. Below a squared Pearson correlation coefficient of 0.4
, we would recommend not to use generated data for downstream analysis, while the generated latent
space might still be useful for analysis.
Cell-wise Coefficient of Variation:
| Metric | Training Value | Validation Value |
|---|
| Mean Absolute Error | 0.71 | 0.76 |
| Pearson Correlation | 0.92 | 0.92 |
| Spearman Correlation | 0.94 | 0.94 |
| R² (R-Squared) | 0.83 | 0.82 |
The gene-wise coefficient of variation summarizes how well variation between different genes is
preserved by the generated model expression. This value is usually quite high.
Gene-wise Coefficient of Variation:
| Metric | Training Value |
|---|
| Mean Absolute Error | 12.95 |
| Pearson Correlation | 0.74 |
| Spearman Correlation | 0.78 |
| R² (R-Squared) | -0.99 |
Differential expression metric
The differential expression metric provides a summary of the differential expression analysis
between cell types or input clusters. We provide here the F1-score, Pearson Correlation
Coefficient of Log-Foldchanges, Spearman Correlation Coefficient, and Area Under the Precision
Recall Curve (AUPRC) for the differential expression analysis using Wilcoxon Rank Sum test for each
cell-type.
Differential expression:
| Index | gene_f1 | lfc_mae | lfc_pearson | lfc_spearman | roc_auc | pr_auc | n_cells |
|---|
| kidney epithelial cell | 0.96 | 0.46 | 0.91 | 0.99 | 0.61 | 0.30 | 9278.00 |
| CD8-positive, alpha-beta T cell | 0.92 | 1.83 | 0.65 | 0.93 | 0.03 | 0.04 | 537.00 |
| B cell | 0.93 | 2.04 | 0.67 | 0.93 | 0.03 | 0.04 | 494.00 |
| macrophage | 0.78 | 1.17 | 0.62 | 0.91 | 0.86 | 0.71 | 339.00 |
| CD4-positive, alpha-beta T cell | 0.93 | 2.55 | 0.64 | 0.89 | 0.03 | 0.03 | 299.00 |
| natural killer cell | 0.86 | 2.99 | 0.59 | 0.86 | 0.08 | 0.03 | 201.00 |
| endothelial cell | 0.68 | 1.45 | 0.69 | 0.92 | 0.51 | 0.35 | 101.00 |
| monocyte | 0.71 | 4.62 | 0.64 | 0.77 | 0.17 | 0.02 | 54.00 |
| non-classical monocyte | 0.62 | 5.20 | 0.63 | 0.72 | 0.15 | 0.02 | 27.00 |
| neutrophil | 0.74 | 6.67 | 0.49 | 0.46 | 0.40 | 0.03 | 22.00 |
| T cell | 0.67 | 7.36 | 0.43 | 0.47 | 0.33 | 0.02 | 17.00 |
| plasma cell | 0.32 | 7.64 | 0.40 | 0.40 | 0.22 | 0.02 | 5.00 |
| intermediate monocyte | 0.21 | 5.75 | 0.55 | 0.55 | 0.27 | 0.02 | 2.00 |
Model Properties
We provide here key parameters used to setup and train the model.
Model Parameters
These provide the settings to setup the original model:
1{
2 "n_hidden": 128,
3 "n_latent": 20,
4 "n_layers": 3,
5 "dropout_rate": 0.05,
6 "dispersion": "gene",
7 "gene_likelihood": "nb",
8 "use_observed_lib_size": true,
9 "linear_classifier": false,
10 "datamodule": null,
11 "latent_distribution": "normal",
12 "use_batch_norm": "none",
13 "use_layer_norm": "both",
14 "encode_covariates": true
15}
Setup Data Arguments
Arguments passed to setup_anndata of the original model:
1{
2 "labels_key": "cell_type",
3 "unlabeled_category": "unknown",
4 "layer": "counts",
5 "batch_key": "donor_assay",
6 "size_factor_key": null,
7 "categorical_covariate_keys": null,
8 "continuous_covariate_keys": null,
9 "use_minified": false
10}
Data Registry
Registry elements for AnnData manager:
| Registry Key | scvi-tools Location |
|---|
| X | adata.layers['counts'] |
| batch | adata.obs['_scvi_batch'] |
| labels | adata.obs['_scvi_labels'] |
| latent_qzm | adata.obsm['scanvi_latent_qzm'] |
| latent_qzv | adata.obsm['scanvi_latent_qzv'] |
| minify_type | adata.uns['_scvi_adata_minify_type'] |
| observed_lib_size | adata.obs['observed_lib_size'] |
Summary Statistics
| Summary Stat Key | Value |
|---|
| n_batch | 2 |
| n_cells | 11376 |
| n_extra_categorical_covs | 0 |
| n_extra_continuous_covs | 0 |
| n_labels | 14 |
| n_latent_qzm | 20 |
| n_latent_qzv | 20 |
| n_vars | 3000 |
Training
Training data url: Not provided by uploader
If provided by the original uploader, for those interested in understanding or replicating the
training process, the code is available at the link below.
References
The Tabula Sapiens Consortium. The Tabula Sapiens: A multiple-organ, single-cell transcriptomic atlas of humans. Science, May 2022. doi:10.1126/science.abl4896
