Problems with MACS ChIP-Seq peak calling
0
0
Entering edit mode
4.4 years ago
dcheng1 • 0

I downloaded raw data in GEO and analyzed it following general ChIP-Seq pipeline: trimming adapter, mapping and calling peaks. However, when compared with bigwig files in GEO, my results look like separate spikes instead of peaks.

Here is the command I used: macs14 --bw 200 -t -c -B -S --call-subpeaks

Does anyone know how did this happen? Thanks!!

enter image description here

ChIP-Seq MACS Peak calling • 1.7k views
ADD COMMENT
1
Entering edit mode

Need more information

We need to know the MACS2 output file you are working from. Do you load this file directly into IGV or create another file from loading it? Always include file extensions Please print the top 10 lines (including header if there is one) of the file you are loading into IGV

ADD REPLY
0
Entering edit mode
track type=bedGraph name="NA_treat_all" description="Extended tag pileup from MACS version 1.4.2 20120305"
chr1    10002   10003   1
chr1    10003   10004   2
chr1    10004   10008   3
chr1    10008   10014   4
chr1    10014   10015   5
chr1    10015   10016   7
chr1    10016   10020   8
chr1    10020   10021   9
chr1    10021   10026   10
ADD REPLY
1
Entering edit mode

you need to create a signal track. Calling peaks is one thing, visualising them is something else. MACS2 will have created 2 .bdg file called [prefix]_treat_pileup.bdg and [prefix]_control_lambda.bdg for your sample. These are used to create the type of file you want. The link describes how to create a FE or LogLR track but I think you need a ppois track.

$ macs2 bdgcmp -t [prefix]_treat_pileup.bdg -c [prefix]_control_lambda.bdg -o [prefix]_ppois.bdg -m ppois
ADD REPLY
1
Entering edit mode

I added markup to your post for increased readability. You can do this by selecting the text and clicking the 101010 button. When you compose or edit a post that button is in your toolbar, see image below:

101010 Button

ADD REPLY
0
Entering edit mode

Link for igv screenshot: https://ibb.co/dXZhvv

ADD REPLY
0
Entering edit mode

MACS call peaks (regions where there are peaks). BigWigs are just the coverage tracks (at least in this case). Those should not look similar.

ADD REPLY

Login before adding your answer.

Traffic: 2235 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6