chip-seq analysis (bowtie2, MACS2)
2
0
Entering edit mode
7.5 years ago

Could you please explain where I made a mistake analyzing ChiP-seq data.

I took fastq file from modencode (set 2639), then mapped it to dm3 genome (bowtie2). Then I took .sam file from its output and called peaks (MACS2).

Final .bed file (with narrow peaks) is completely different from gff3 file from modencode:

Number of peaks differs about 2 times.

bowtie2 macs2 alignment fastq • 6.1k views
ADD COMMENT
2
Entering edit mode
7.5 years ago

It really depends on , what tools for mapping and peak calling were used along with the parameters. Slight change can cause differences, make sure the peakset you downloaded, has been analysed the same way.

ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia

ADD COMMENT
0
Entering edit mode

So, if I have only raw data from modencode, should I use their guidelines to analyze it? I know my question seems stupid, sorry for that, I'm a newbie in bioinformatics.

ADD REPLY
0
Entering edit mode

It's not stupid, you can use any method for analysing, if you know what you are doing. But using different methods to get the exactly same result is difficult. Mapping might be still fine (using different mapping algorithms) but peak callers can generate huge variability, depending on parameters and use of controls.

ADD REPLY
0
Entering edit mode

Ok, I understood you. To be clear: my goal is to define which genes are regulated by the TF. I'm trying to determine binding sites for TF. But if results of peak calling are completely different I suppose that the final results of analysis will differ. Am I right or not?

And could you please advise how to determine which parameters I should use in MACS2 (mfold, q-value, p-values etc.) and in bowtie2 (alignment options, input options etc.)? Maybe any manual exists or learning videos. I tried using only default parameters earlier.

ADD REPLY
2
Entering edit mode
7.5 years ago
Ming Tang ★ 2.7k

The results can be very different with different analyzing methods. No surprise at all.

ADD COMMENT

Login before adding your answer.

Traffic: 2460 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