scArches - single-cell architecture surgery¶
scArches is a package to integrate newly produced single-cell datasets into integrated reference atlases. Our method can facilitate large collaborative projects with decentralized training and integration of multiple datasets by different groups. scArches is compatible with scanpy. and hosts efficient implementations of all conditional generative models for single-cell data.
From version 0.3.0 we moved to PyTorch, however, you can still install the latest Tensoflow version by installing 0.2.4.
What can you do with scArches?¶
Construct single or multi-modal (CITE-seq) reference atlases and share the trained model and the data (if possible).
Download a pre-trained model for your atlas of interest, update it with new datasets and share with your collaborators.
Project and integrate query datasets on the top of a reference and use latent representation for downstream tasks, e.g.:diff testing, clustering, classification
What are the different models?¶
scArches is itself an algorithm to map to project query on the top of reference datasets and applies to different models. Here we provide a short explanation and hints on when to use which model. Our models are divided into three categories:
This class of algorithms require no cell type labels, meaning that you can create a reference and project a query without having access to cell type labels. We implemented two algorithms:
scVI (Lopez et al., 2018): Requires access to raw counts values for data integration and assumes
count distribution on the data (NB, ZINB, Poisson).
trVAE (Lotfollahi et al.,2020): It supports both normalized log transformed or count data as input and applies additional MMD loss to have better merging in the latent space.
- Supervised and Semi-supervised
This class of algorithms assumes the user has access to cell type labels when creating the reference data and usually perform better integration compared to. unsupervised methods. However, query data still can be unlabeled. In addition to integration, you can classify your query cells using these methods.
scANVI (Xu et al., 2019): It needs cell type labels for reference data. Your query data can be either unlabeled or labeled. In the case of unlabeled query data, you can use this method to also classify your query cells using reference labels.
scGen (Lotfollahi et al., 2019): This method requires cell-type labels for both reference building and query mapping. The query mapping for this method solely relies on the integrated reference and requre no fine-tuning.
- Bioligically informed
expiMap (Lotfollahi, Rybakov et al., 2022): This method takes prior knowledge from gene sets databases or users allowing to analyze your query data in the context of known gene programs.
These algorithms can be used to construct multi-modal references atlas and map query data from either modality on the top of the reference.
totalVI (Gayoso al., 2019): This model can be used to build multi-modal CITE-seq reference atalses. Query datasets can be either from sc-RNAseq or CITE-seq. In addition to integrating query with reference, one can use this model to impute the Proteins in the query datasets.
Which model to choose?¶
If your reference data is labeled (cell-type labels) and you have unlabeled or labeled query then use scArches scANVI.
If your reference data is labeled (cell-type labels) and you have labeled query then use scGen.
If your reference and query are both unlabeled our preferred model is scArches scVI and if it did not work for you try scArches trVAE.
If you have CITE-seq data and you want to integrate RNA-seq as query and denoise the proteins in the RNA-seq then use scArches totalVI.
Where to start?¶
To get a sense of how the model works please go through this tutorial. To find out how to construct and share or use pre-trained models example sections.
If scArches is useful in your research, please consider citing the preprint.