[suggestion] what would be the appropriate pipeline to perform joint embedding of GEX and ATAC?

General Multiome
, ,
1

I have processed my snRNA of multiome using scanpy and snATAC using snapatac2. Both have been annotated as well.
I would like to perform joint embedding of rna and atac but the var attributes are different and so are the number of cells in both the modalities.
The overarching goal after embedding is obtain cell_data that will be used as an input for pycistopic->scenic+.

2

Hey @yojetsharma,

If that might help, this sounds like the setting described in the tutorials on mudata and muon, e.g. here.

Briefly, you could store two modalities in a MuData object, and use methods like MOFA+ and MultiVI to compute a join embedding, which you can store in the same MuData object.

1 Like