Question: Comparison of replicates from NGS : merge before or after peak calling ?
0
gravatar for anais1396
17 months ago by
anais139620
Brussels
anais139620 wrote:

Hello everyone !!

I have a question for you : I have biological replicates from Atac-seq experience and I want to perform a peak calling and then a comparison of the overlaps between my files. I have 5 samples for healthy people and 9 samples for patients and I have 2 strategies :

  • first, I have bam files from alignment and I merge the 5 bam for healthy together and make a second merge for the 9 samples of patients (with samtools merge) and thus I obtain 2 files (bam) and I perform peak calling on each and then search for overlap between the two sets of peaks with intersectBed

  • second, I have bam files from alignment (so I have a total of 14 bam files), I perform peak calling on each of the 14 samples and then I make a merge (with bedtools merge) of the 5 peak files of healthy and a second merge of the 9 peaks files of patients and then I search for overlap between the 2 sets with intersectBed.

Which one is the best strategy and why ?? Should I merge before or after the peak calling ?

Thank you in advance !!

Anaïs

merge ngs peak calling • 733 views
ADD COMMENTlink modified 17 months ago by h.mon28k • written 17 months ago by anais139620
0
gravatar for h.mon
17 months ago by
h.mon28k
Brazil
h.mon28k wrote:

It depends on the sequencing coverage of each sample, but ideally each sample have been sequenced to an appropriate coverage, map and call peaks separately, and then merge peaks. If you have low sequencing coverage per sample (I hope you don't), merge before calling.

For more extensive answers, this issue has been covered before, more than once. Searching for merge peak calling would have led you to several posts, such as:

ChIP-Seq: Calling peaks with replicates

how to pool together biological replicates?

ChipSeq: merge bam file before peak calling

Chip-Seq merging peak files

ADD COMMENTlink written 17 months ago by h.mon28k

How can I simply check the coverage of my samples ?

ADD REPLYlink written 17 months ago by anais139620
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2283 users visited in the last hour