Question: Batch and combining reads from two flow cells for the same ChIP/ATAC-seq library
0
gravatar for reskejak
2.2 years ago by
reskejak20
Michigan State University
reskejak20 wrote:

I'm sure this has been answered previously, but I have been unable to find a directly-similar question here for ChIP-seq (or similar) data.

We have a library which has been sequenced twice, on two flow cells, corresponding to one biological replicate sample from an experiment (each flow cell is an independent experiment). This library has the lowest complexity output from the group of samples assayed, so we would like to combine reads from the two flow cells.

Qualitatively, we have observed that reads from both flow cells appear highly similar (in terms of mapping stats, duplicates, complexity, signal-to-noise, FRiP, peaks called, UCSC genome browser, etc.) and plotting PCA with DiffBind shows these samples to be nearly identical with or without modeling batch.

Can we directly combine these reads by way of cat or similar? What about combining bam files as in samtools merge? How else can we assess the appropriateness of combining these libraries, and furthermore which would be the most practical method to doing so?

ADD COMMENTlink modified 2.1 years ago by Devon Ryan96k • written 2.2 years ago by reskejak20
3
gravatar for Devon Ryan
2.1 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

As a general principal, if you just sequenced the same library twice then go ahead and merge the results (it doesn't much matter if you cat the fastq files together or use samtools merge on the BAM files, you'll get essentially the same results regardless).

In your particular case you've also nicely assessed the comparability of the technical replicates. That's great, well done! There's no hard and fast rule for exactly how similar technical replicates need to be to merge them. But if you look at the QC metrics you created and the runs look as similar as yours sound then you can safely merge the runs.

ADD COMMENTlink written 2.1 years ago by Devon Ryan96k
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: 1068 users visited in the last hour