Merging data from multiple cohorts and many donors with scVI

[adamgayoso]

Hi Phil,

@Phil_Bradley

is there a recommended minimum number of cells per donor ie batch (I am using the donor as the “batch” identifier)? I understand that there are suggestions to keep the number of cells greater than the number of genes, but I’m not sure if this applies within individual batches as well.

It’s hard to say. What you’re attempting to to is beyond what we have explicitly tested (and very cool!).
We are currently working on an extension related to this post:

Models for handling hundreds of batches

I think this is a great idea! Very curious to know what you find.

@Phil_Bradley

It looks like there are systematic technical differences in expression between the cohorts. Is there a way to include the cohort as an additional “batch” covariate for the model fitting? Does that even make sense, given that each donor has a cohort membership so the model already has freedom to fit those differences on a per-donor basis?

Yes you can put both the donor and the cohort as a key using categorical_covariate_keys here: https://docs.scvi-tools.org/en/stable/api/reference/scvi.data.setup_anndata.html#scvi.data.setup_anndata

@Phil_Bradley

Any other thoughts/suggestions you might have on integrating very many batches would be much appreciated! For example, would it be better to integrate a few big batches and then bring the rest into that latent space using something like scarches? Are there SCVI flags (use_layer_norm, use_batch_norm, etc?) that might be appropriate for handling many cells/batches?

I would be interested to learn a bit more about what you’re trying to do. What you stated are reasonable things to try though. Please feel free to email me (firstlast at berkeley dot edu) if you’d like to schedule a meeting!