Hello,
I am working with a single-cell RNA-seq dataset and would like to use Harmony for batch correction (sc.external.pp.harmony_integrate in Scanpy). I am not sure which key to use.
My experimental design includes multiple biological conditions (control, injury, and control or injury + treatment at different timepoints), totaling 12 demultiplexed samples. Samples were multiplexed using hashtag oligos (HTO) and demultiplexed after sequencing. The libraries were generated in 6 pools (each pool containing 2 biological samples), and each pool was prepared separately. Some pools were sequenced in a separate runs.
Should the Harmony batch key correspond to the 12 individual demultiplexed samples, or to the 6 original multiplexed pools? If I overcorrect the data would I realize?
Thank you
1 answer
My best practice would be to do the 12 individual samples. After all, you have two major sources of variation: a) the library preps of each of the 12, and the conditions of each of these. Both you do typically not want reflected in clustering and visualization. You want celltypes cluster by celltypes. DE analysis, using the raw data will then take care of the condition differences. That's at least how I have been doing this over the years. I don't think overcorrection is really.a thing here. After all, you probably will have a good idea which celltypes to expect, so check whether these can be found and separated. You could start by assigning celltypes using reference labels, e.g. using SingleR first, so this whole clustering-UMAP typical workflow is deemphasized and you immediately see which celltypes are present.
Log in to answer this question.