Question: Appropriate peak length of ChIP-seq. (Histone modification having different pattern at gene body & promoter & intergenic regions)
gravatar for chiefcat
21 months ago by
chiefcat100 wrote:

Hi everyone,

Nice to meet you all. I'm new to NGS data analysis. I have received H3K9me3 ChIP-seq. data, which peak calling was conducted by MACS2. I've observed that there are different enrichment pattern of the K9 signal (see below). I am concerning that the parameter for peak calling is not appropriate since the broad signal was called as separate small peaks. I know the --broad argument can be applied when handling broadly distributed histone modification signal, so that the nearby singal can be linked into a broad peak. However, if --broad is applied, will it affect the calling of those sharp peaks? Thanks.

: Relatively sharp peak at promoter

enter image description here Broad signal over gene body

Sharp peak at intergenic region: Sharp peak at intergenic region

ADD COMMENTlink modified 21 months ago by geek_y9.6k • written 21 months ago by chiefcat100

From my experience with ChIP data, it takes a lot of trial and error to find the correct parameters in order to identify your peaks. With the --broad parameter, I believe that there's an associate p-value cutoff parameter that's used when extending the peak from the initial peak region.

Looking at your ChIP signals, using the --broad parameter, you may end up extending peak regions too far as there appears to be a relatively large background signal. However, if your sharp peaks are well-resolved and have many reads mapping to them, it should work such that MACS2 identifies them perfectly even with the --broad parameter.

Also, consider running two types of analyses, if that works, i.e., one that identifies the narrow peaks and one that separately identifies the wider peaks.

Also, consider the other options such as HOMER, SICER, or bamCoverage.


ADD REPLYlink written 21 months ago by Kevin Blighe42k
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: 725 users visited in the last hour