Question: how to pool together biological replicates?
2
gravatar for ahsan.raza.rana
3.6 years ago by
ahsan.raza.rana40 wrote:

I have 2 biological replicates (Chip-Seq) and i want to know how to pool them together so that after pooling them together i may apply Homer to find peaks. I am doing it for the very 1st time thats y asking such a q.

Now my q is : I have found the peaks of replicates separately by Homer.

1.How to find the peaks higher than 1 read per million (rpm) from each replcate and then merge the peaks of replicates, So, how can I merge the peaks of replicates???????????????

2.Can I use Normalized tag count column to check for peaks higher than 1rpm.If yes,then what would be the interpretation of 145.7 normalized tag count. also what would be the value to check for >than 1rpm.

  1. how can i merge the peaks ?by mergeBed with -n option will be ok or not

Waiting for reply

chip-seq • 3.6k views
ADD COMMENTlink modified 3.6 years ago • written 3.6 years ago by ahsan.raza.rana40

I would also consider calling peaks on the separate replicates and then intersect the two peak sets to get a consensus.

ADD REPLYlink written 3.6 years ago by dariober10k
5
gravatar for robm9119
3.6 years ago by
robm9119150
United States
robm9119150 wrote:

You can pool the reads and map them to get a single bam file. Alternatively you can merge bam files from the two replicates using samtools merge.

By the way, why would you want to merge the replicates? If the replicates are adequately sequenced, a better type of analysis would be to call peaks for the replicates separately and only keep the peaks that agree between the two replicates.

ADD COMMENTlink written 3.6 years ago by robm9119150
1

IDR and other measures of reproducibility are discussed in the ENCODE guidelines for ChIP-seq

ADD REPLYlink written 3.6 years ago by fanli.gcb690

How can I keep the common peaks between replicates? how to do this ? I have no idea?

ADD REPLYlink written 3.6 years ago by ahsan.raza.rana40
2

Bedtools intersect:

http://bedtools.readthedocs.io/en/latest/content/tools/intersect.html

ADD REPLYlink written 3.6 years ago by robm9119150

Or mergePeaks from HOMER.

ADD REPLYlink written 2.9 years ago by boczniak767680
1
gravatar for EVR
3.6 years ago by
EVR570
Earth
EVR570 wrote:

you can either use cat function to pool all the fastq files or samtools merge to merge the bam files.

ADD COMMENTlink written 3.6 years ago by EVR570

thanks alot.I have sam files. Can i use samtools merge function to merge the sam files?

ADD REPLYlink written 3.6 years ago by ahsan.raza.rana40

I am not sure whether to use samtools merge as it help says only bam files. My safe bet would pool the fastq files and generate sam file with appropriate program

ADD REPLYlink written 3.6 years ago by EVR570
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: 1848 users visited in the last hour