search

10XGenomics / cellranger

10x Genomics Single Cell Analysis
https://www.10xgenomics.com/support/software/cell-ranger
Other
340 stars 91 forks source link

cellranger aggr is unable to combine cellranger multi output #195

Closed dmalzl closed 1 year ago

dmalzl commented 1 year ago

Hi,

I am currently using cellranger to process some 10x datasets which were sequenced with Total-seq A. Unfortunately, the antibody to sample assignment was inconsistent with some even having different antibody combinations. I realized that for the cellranger multi runs the cmo reference is also added to the feature reference in the molecule_info.h5 which causes aggr to raise an error and abort. I for now resorted to concatenating the data manually using AnnData with the sample_feature_bc.h5 file but I keep getting very weird results when using scanpy's umap. I am yet to investigate why this is and I also do not completely understand how aggr works but I wonder if there is any magic happening when using aggr as compared to plain conatenation of the adatas.

dmalzl commented 1 year ago

In any case, I guess this might help future users who are in the same situation as me ;)

dmalzl commented 1 year ago

Nevermind, I had a bug in my config file generator which caused the same sample being analysed multiple times resulting in weird umap artifacts. Anyway I guess knowing that you also should use the same CMO reference might help someone (I hope)