Harmony batch correction vs. regressing out donor effect in Seurat
1
5
Entering edit mode
5.2 years ago
Lucy ▴ 160

Hi,

I am performing combined analysis of three scRNA-seq samples from different donors generated using 10x Genomics technology. I am trying to decide whether to use Harmony (https://www.biorxiv.org/content/10.1101/461954v2) to remove donor effects or whether to regress out donor during the sctransform step (https://satijalab.org/seurat/v3.0/sctransform_vignette.html) in the Seurat pipeline.

What are the relative advantages and disadvantages of the two approaches?

Best wishes,

Lucy

Seurat scRNA-seq Harmony RNA-seq batch-effect • 7.8k views
ADD COMMENT
0
Entering edit mode

You might also find this recent discussion helpful: about batch correction in scRNA-seq

ADD REPLY
4
Entering edit mode
5.2 years ago

You should not use SCTransform to regress out batch effects. Rather, you should use one of Seurat's integration methods. Seurat also has a number of wrappers around different integration methods, including Harmony.

Personally, fastMNN has worked well for me, but it's usually worth trying a few methods, as they don't all perform similarly across all datasets.

ADD COMMENT
0
Entering edit mode

Thank you. What is the reason to not use SCTransform to regress out batch effects?

ADD REPLY
1
Entering edit mode

The authors don't recommend it, as it's only meant to remove minor variation/biases due to differences in cell cycle, reads in mitochondrial genes, etc. rather than batch effects that affect every cell in a much more systemic manner. Their docs are a bit of a mess, so I can't find a specific reference, but if you dig around their github issues page, you can find threads that talk about it. Additionally, you can read their preprint about sctransform and/or their paper on integration specifically that shed more light on this.

ADD REPLY
0
Entering edit mode

Great, thanks for the help!

ADD REPLY
0
Entering edit mode

Hi @Jared.andrews07,

It is an interesting response for me. I have now tried both Harmony and Seurat Integration. But how did you try FastMNN? For eg; is there a package which I can use? Or It will be great if you can share a public github link that was useful for you or any notebook that you have prepared yourself?

Thanks in advance!

ADD REPLY
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Hi,

I have 3 patients with normal and tumor tissue sample(10× technology), there are six samples in total. I want combine them and find whether there is some difference between normal and tumor sample in specific cell type (such as immune cell). I am struggle to choose integrate method. For the seurat integrate method, SCtransform method, harmony method, which one is more suitable for my case, further, should I regress out donor when perform scaleData in seurat? Thanks a lot, hope to get your suggestions!

Best, Wei

ADD REPLY
0
Entering edit mode

Please post this as a new question.

ADD REPLY
0
Entering edit mode

Hi jared,

I have post a new question, hope to get your suggestions!Thanks a lot! Here is the new post

Best, Wei

ADD REPLY

Login before adding your answer.

Traffic: 1885 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6