Question: calling peaks without control
gravatar for jenyka_dm
8 months ago by
jenyka_dm10 wrote:

Hi again! I was wondering, when I call peaks with MACS using my treatment and control, I get 1600 peaks. However, when I call the peaks using treatment or control only, I get treatment (350 peaks) and control (600 peaks). I was under the impression that it should have given me a larger number of peaks for my treatment only sample, and that the treatment + control peaks are basically the peaks from the ;control sample only' subtracted from the 'treatment only' sample. Could you please explain me why is this so? Please forgive my ignorance, I am very new at this stuff. Thanks!

chip-seq • 369 views
ADD COMMENTlink modified 8 months ago by Jennifer Hillman Jackson25k • written 8 months ago by jenyka_dm10

I am also confused to why the peaks called separately for treatment do not reflect the peaks called for treatment against the control? Thanks!

ADD REPLYlink written 8 months ago by jenyka_dm10
gravatar for Jennifer Hillman Jackson
8 months ago by
United States
Jennifer Hillman Jackson25k wrote:


The control is reducing the number of "false peaks" (exist as a background) when used with the treated sample. Use a control whenever possible. You can play around with the parameters and compare the different results, along with a reference annotation - if you have one - to better interpret the differences. Try visualizing all results in a browser for a distinct genome region or two (zoom into a region) to make this comparison less complex.

The baseline peaks are called using mapped BAM(s) plus the parameters set for the tool. Calling peaks with just the input or treated sample alone can produce in spurious peak calling. Also, parameters work best when used in combinations that target a specific peak calling goal. Try the recommended combinations from the tool form help (including the tool manual) and compare the results to those obtained when using your own settings.

This is related to your other question and I encourage you to review the tool documentation and explore the MACS2 wiki and google group to better understand how it works. The full suite of MACS2 tools/functions are incorporated into Galaxy and line command usage can be replicated with the Galaxy wrappers around these tools/tool form options.

All of this assumes that mapping went well and that there are no technical errors (incorrect or different genome mapped against for one of the mapping runs, and other common usage issues).

Thanks! Jen, Galaxy team

ADD COMMENTlink modified 8 months ago • written 8 months ago by Jennifer Hillman Jackson25k
Please log in to add an answer.


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