Question: how to filter ChIP-Seq peaks called by MACS based on the measurements?
gravatar for Nitin Narwade
3.1 years ago by
Nitin Narwade440
Nitin Narwade440 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 • 1.5k views
ADD COMMENTlink modified 3.1 years ago by Petr Ponomarenko2.6k • written 3.1 years ago by Nitin Narwade440

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

ADD REPLYlink written 3.1 years ago by Lila M 810

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 3.1 years ago by Sinji3.0k
gravatar for Petr Ponomarenko
3.1 years ago by
United States / Los Angeles /
Petr Ponomarenko2.6k 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 3.1 years ago by Petr Ponomarenko2.6k
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: 1340 users visited in the last hour