Question: macs2 normalization on broadPeak
2
gravatar for firatuyulur
12 weeks ago by
firatuyulur170
firatuyulur170 wrote:

Hi,

I have a number of atacseq bams and their total read numbers differ. I did peak calling with and without --SPMR for normalization and there is no difference between the two ......broad_peaks.broadPeak files. yet there is a difference in the resulting two ....broad_treat_pileup.bdg files. While using callpeak, does --SPMR have no effect on the signalValues? If that is the case, should I do another peak calling from the bdg files?

My end goal is to see what happened to some specific regions, whether they got condensed or decondensed and I cannot make sure whether this read depth difference effect the signalValues.

macs2 callpeak atacseq • 278 views
ADD COMMENTlink modified 12 weeks ago by James Ashmore2.1k • written 12 weeks ago by firatuyulur170
2
gravatar for James Ashmore
12 weeks ago by
James Ashmore2.1k
UK/Edinburgh/MRC Centre for Regenerative Medicine
James Ashmore2.1k wrote:

Whether you use the SPMR flag or not, MACS2 will normalise the data internally to call peaks. The SPMR flag only effects the signal track produced. With the flag present the signal will be normalised to reads per million, this is for comparison with other samples which have been sequenced to different read depths.

ADD COMMENTlink modified 12 weeks ago • written 12 weeks ago by James Ashmore2.1k

So what you are saying is the signalValues that I see in broadPeaks file are independent of the read depth as those are the values obtained from a normalised peak calling and I can compare two samples according to their signalValues although they have different read depths?

ADD REPLYlink written 12 weeks ago by firatuyulur170

The signalValue in the broadPeak file produced by MACS2 is not a measure of how many reads are aligned within the peak. It is a measure of fold-change above background.

ADD REPLYlink written 12 weeks ago by James Ashmore2.1k

Isnt fold change related to the number of reads aligned to the peak?

ADD REPLYlink written 12 weeks ago by firatuyulur170

James thanks for the help. After this quick conversation I started looking at the regions of peaks on igv. I put the spmr called bedgraphs to igv and checked broadpeaks beds by eye. Still cannot strongly argue about it, I suppose you are right about peak calling normalizes while doing its job.

ADD REPLYlink written 12 weeks ago by firatuyulur170
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: 1550 users visited in the last hour