Hello,
I have just uploaded two bw files to IGV, they are from the same experiment, only the second one was sequenced more. When I load them more I notice that the file sequenced more has peaks that are lower, why is that in your opinion?
N.B. they were processed in the same way so the steps from fastq files to bw files are the same but the result is very different.
I attach the photo to make you understand.
What could be the explanation?
Thanks for helping me
•
written
by
michelafrancesconi9 •
1 answer
Independent IP trials from an experiment can behave differently. Signal to noise ratio in a data set is much more important than sequencing depth. If you look at figure 5 of Landt et al. (ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia Genome Res 2012), you'll see 4 replicate IPs of the same factor, 2 replicates each from two different experiments, and you'll see that the peaks are different heights, and each data set has unique characteristics in terms of IP enrichment. That paper also has a variety of recommendations for how to examine enrichment and compare replicates. Things like cross correlation plots, Fraction of Reads in Peaks (FrIP), and IDR analysis.
Log in to answer this question.
-
Troubleshooting poor DP values
written by manfi •
-
About batch effect in separate library preparation for multiple samples tested on same day
written by yujunfhf •
-
Problem with peak calling using hiddenDomains: entire chr12 lacking peaks.
written by carolin.mattausch •
-
No coverage in macs2 peak summit!
written by blur
-
deepTools computeMatrix need at least one array to concatenate error
written by leeb4 •
-
Are reads with the same length and position always "duplicates"
written by elcortegano
-
Problems with proper interpretation of visualized bam file from Illumina Miniseq platform
written by dominik.saburo •
-
Chip-Seq Peak height difference between BED and BW file formats.
written by morovatunc
-
Paired-end ChIP-Seq Peak Calling
written by simon.pearce •
-
Cisgenome Browser With Bar Files Vs Igv With Wig Files
written by Hamilton •
Hello.
What is the source of the BigWig files? Were they normalised? Do they contain duplicated reads/fragments? It could be the deeply sequenced sample contains more background (non-specific mapping).
I created bw files from bam files and filtered them so as to remove duplicates from both
I would use something like Deeptools bamCompare to generate bigWig files that are normalised by sequencing depth: https://deeptools.readthedocs.io/en/develop/content/tools/bamCompare.html
If the peaks are still shorter in the deeper sequenced sample then there is higher background (non-specific signal).
"from the same experiment"...are they two aliquots of sequencing from one library? Or are they replicates? (different IP samples but done within the same experiment). Replicates from ChIP seq (even within the same experiment) can vary quite a bit, such that sequencing depth would be beside the point.
They are two different replicates fro ChIP seq but in your opinion is normal?