Question: how to filter ChIP-Seq peaks called by MACS based on the measurements?
gravatar for Nitin Narwade
20 months ago by
Nitin Narwade320
NCCS, Pune
Nitin Narwade320 wrote:

Hello everyone, I am a beginner in this field.

Actually, my question is about ChIP-Seq data analysis, I have tried ChIP-Seq data analysis for BRCA1 protein (Human) with two replicates. The data quality is good and I have successfully reached at the step of peak calling.

I have used bowtie for raw read alignment and MACS14 for peak calling. At the end of this I got 96483 peaks.

Is it possible that one successful ChIP-Seq experiment could have these many number of peaks? OR if I want to reduce the number of peaks which criteria is better? I have already tried FDR (%) [reported at]. but by using FDR (%) filter I hardly able to discard 10 peaks. In case of pValue criteria peaks have been already filtered by MACS14.

Any help would be really appreciated.

Thank you....!!

chip-seq • 864 views
ADD COMMENTlink modified 20 months ago by Petr Ponomarenko2.5k • written 20 months ago by Nitin Narwade320

In my opinion you should choose a p-value according with your data :)

ADD REPLYlink written 20 months ago by Lila M 390

If you have replicates then you should look into conducting some IDR analysis as done by ENCODE. Essentially you call peaks on each replicate individually and then you combine the 'true' peaks that overlap each other in both files.

ADD REPLYlink written 20 months ago by Sinji2.7k
gravatar for Petr Ponomarenko
20 months ago by
United States / Los Angeles /
Petr Ponomarenko2.5k wrote:

What is the protein you are analyzing, what is the binding area and how many there should be (compared to other known experiments)? Most likely you can find public data for the same species, protein, sequencing machine and chemistry at NCBI GEO and a corresponding paper in a good peer-reviewed journal. You can try your pipeline on that dataset, optimize parameters, compare with published data. This will give you a hint on how to analyze your own dataset and how to optimize parameters for it. 96483 peaks can be just right, or too little, or too few. Also, direct comparison with another replicate still can give you consistent artifacts that will be preserved between two different experiments just because of the systematic errors or (what is more likely) having certain areas of the genome called erroneously based on their sequence complexity.

ADD COMMENTlink written 20 months ago by Petr Ponomarenko2.5k
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: 1258 users visited in the last hour