Question: Harmony batch correction vs. regressing out donor effect in Seurat
0
gravatar for Lucy
8 months ago by
Lucy40
Lucy40 wrote:

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

ADD COMMENTlink modified 8 months ago by jared.andrews075.5k • written 8 months ago by Lucy40

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

ADD REPLYlink written 8 months ago by igor10k
1
gravatar for jared.andrews07
8 months ago by
jared.andrews075.5k
St. Louis, MO
jared.andrews075.5k wrote:

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 COMMENTlink written 8 months ago by jared.andrews075.5k

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

ADD REPLYlink written 8 months ago by Lucy40
1

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 REPLYlink written 8 months ago by jared.andrews075.5k

Great, thanks for the help!

ADD REPLYlink written 8 months ago by Lucy40

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 REPLYlink written 5 months ago by saeedabdurahman0

See the Seurat wrappers page.

ADD REPLYlink written 5 months ago by jared.andrews075.5k

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 REPLYlink written 3 months ago by 6232022150

Please post this as a new question.

ADD REPLYlink written 3 months ago by jared.andrews075.5k

Hi jared,

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

Best, Wei

ADD REPLYlink modified 3 months ago • written 3 months ago by 6232022150
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1019 users visited in the last hour