Question: ATAC-Seq: Filter bams for differential analysis or not
gravatar for srhic
8 months ago by
srhic40 wrote:


I am doing some analysis to identify differential accessible regions between two conditions using atac-seq. The workflow I have followed uses peak calling with macs2 and then differential analysis of these peaks with deseq2. This gave me good results with expected genes changing accessibility .

However, after reading more about atac-seq I am a little confused if my analysis is correct. As atac-seq produces both nfr and nucleosomal fragments, it seems that some workflows recommend filtering the bams for nfr fragments <100bp and then using them for further analysis. However, some workflows also use all reads without filtering which is what I did.

My question is if my current analysis makes sense or if I should repeat it after filtering my bams to <100bp fragments? I am also concerned that in that case I am discarding more than half of all my data. What is the best practice in this case? Also if I do the filtering, should I use the filtered bams for both peak calling and and deseq2 or for only one of these steps?


atac-seq • 334 views
ADD COMMENTlink modified 8 months ago by Devon Ryan98k • written 8 months ago by srhic40
gravatar for Devon Ryan
8 months ago by
Devon Ryan98k
Freiburg, Germany
Devon Ryan98k wrote:

FYI, it's <120 bases that people are filtering at in most common organisms.

Whether to look at only nucleosome free regions or not depends on your goal, which only you can know. If you want to assess accessibility to transcription factors and such then you need to either look at shorter fragments or use an ATAC-specific peak caller (there are one or two out there). If instead you want to more generally measure compactness (e.g., changes between eu- and heterochromatin states) then not filtering makes sense. Of course in the latter case peak calling with MACS2 would be largely useless.

ADD COMMENTlink written 8 months ago by Devon Ryan98k

Thanks. I want to do a general comparison between two mouse cell lines one of which is a control and the other lacks a chromatin remodeling protein. Can you explain why you think MACS2 peak calling and differential analysis will not be useful for a general analysis? I have seen some papers using macs2 broad peaks for such analysis.

I am also looking at an atac-seq specific peak caller now (hmmratac) which takes into account nfrs and nucleosomes but again the confusion is once I have the peaks, should I use the unfiltered bams for deseq2 differential analysis?. I will try both and see how different the results are but I am not clear about this conceptually.

ADD REPLYlink written 8 months ago by srhic40

MACS2, even in broad mode, will miss larger-scale remodeling. Its broad mode works by merging relatively close called peaks together, which will typically fail for large changes since there will be few if any peaks to merge.

ADD REPLYlink written 8 months ago by Devon Ryan98k
Please log in to add an answer.


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