Biostar Beta. Not for public use.
Calling differential peaks on ATAC-seq data
Entering edit mode
13 months ago
ldyer2006 • 30

My lab has recently run an ATAC-seq analysis on 3 biological conditions (Day 0, day 1 and day 7) with two replicates assigned to each. Once the data is finally generated, I will need to use some differential peak calling software to identify regions of differential accessibility between the 3 conditions. So far, I've thought that MACS2 would be the best option for this (The bdgdiff tool), but after reading around, THOR seems like it may be easier to use and provide better results.

Can anyone give me some advice about the potential benefits or limitations of using these two programs with ATAC-seq data?

Entering edit mode


For ATAC-seq I have used MACS2 for peak calling, followed by edgeR to see if peaks are different between groups. edgeR works for counts data, which ATAC-seq data is when you count all the reads that map to the peak regions.

Entering edit mode

Have you completed this? I am trying to figure out how to do it now.

Check out this post: What are the best tools for differential detection between ATAC-seq samples?

I have used MACS2 and it took a while to get my head around it for ATAC but actually it turned out great. I have 4 conditions (each condition has 8 time points).

So now I am looking into getting differential expression.

Did you get your answer? - if so, what do you suggest, and what links were the biggest help?



Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1