Data integration single cell using Harmony

Question

[_deb]

Dear all,

I have six samples from a 10x experiment: three under control conditions and three subjected to a treatment. I am performing QC preprocessing on each sample separately. Now, I want to integrate the samples using Harmony before proceeding with downstream analyses. My question is: Is it better to run the RunHarmony() function on all preprocessed and merged control samples, then do the same for the preprocessed treated samples, and finally merge the two resulting integrated objects? Or should I first merge all six preprocessed samples into a single Seurat object and then run RunHarmony once on the combined dataset? Any suggestions on which way is better? Those are very large files...

[Arup Ghosh]

Harmony can integrate over multiple covariates, allowing you to provide both control/treatment and batch information, and run it once. Running Harmony separately will result in incompatible embeddings.

https://portals.broadinstitute.org/harmony/articles/quickstart.html#harmony-with-two-or-more-covariates

Answer 1

[yura.grabovska]

Harmony requires PCA reduction as an input so you should merge your data and carry out PCA reduction as a minimum. If you're using Seurat v5 it will now normalise each individual layer separately until you run JoinLayers() after getting a good integrated reduction from Harmony or another method like RPCA, sciVI etc