Question: Can be merged biological replicates?
gravatar for star
12 months ago by
star190 wrote:

Can we merge/ combine biological replicates from same sample together in ChIP-seq / ATAC-seq data when we want to do peak calling (or like RNA seq we should process each replicate separately)? because I am seeing in some papers, the data were merged for peak calling and some where else processed separately!

Thanks in advance!

ADD COMMENTlink modified 11 months ago by ZZzzzzhong200 • written 12 months ago by star190

I'd perform peak calling on individual samples and then use bedtools intersect or mergePeaks from Homer package to assess the overlap among replicates.

ADD REPLYlink modified 12 months ago • written 12 months ago by Sej Modha4.3k

I have used that approach before for ATACseq. So first call peaks for all bam files at once, and then quantify how many reads were mapped to the peak regions for each sample (bam file) for statistical analysis in edgeR.

ADD REPLYlink written 12 months ago by Benn7.5k
gravatar for ZZzzzzhong
11 months ago by
ZZzzzzhong200 wrote:

For the experiment with several replicates, it is recommended to concatenate several ChIP-seq treatment files into a single file. To do this, under Unix/Mac or Cygwin (for windows OS), type:

$ cat replicate1.bed replicate2.bed replicate3.bed > all_replicates.bed

For BAM or SAM files, samtools can be used to combine replicates.

at MACS2 Manual( recommend to merge biological replicates.But you can also call peak separately then use overlapping peaks

ADD COMMENTlink modified 11 months ago by ATpoint21k • written 11 months ago by ZZzzzzhong200

If you merge replicates prior to sequencing, it might be useful to downsample the data to have equal readcounts. Otherwise, the peak calling is dominated by the sample with the highest counts. See this post on BioC for details.

Also @ZZzzzzhong, if you copy statements from a website (the first sentence of your answer), please indicate it by using the Blockquote option in the formatting bar and putting in a link to the reference directly behind it.

ADD REPLYlink modified 11 months ago • written 11 months ago by ATpoint21k
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: 765 users visited in the last hour