Visualising the number of overlapping peaks in ChIP-seq studies
I am usually asked if I can show the overlap between peaks in ChIP-seq studies using a venn diagram. However a venn diagram supposes the relationship is 1:1, where as a single peak in one sample could overlap multiple peaks in another sample. My question is how do you normally calculate the amount of peak overlap and how do you prefer to visualise/summarise this overlap?
• 9,234 views
•
link
updated
by
Ram
•
written
by
James Ashmore
•
written
by
James Ashmore
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Interpret Bedtools Overlap?
written by kstangline •
-
Finding common peaks between F-seq peak region files
written by a.rex
-
Alternatives to Venn Diagrams
written by jsneaththompson
-
Creating reproducible peak regions from multiple ChIP-seq studies
written by James Ashmore
-
determining peaks in a bedfile
written by Daniel James •
-
P-value for enrichment of ChIP-Seq peaks within a genomic annotation
written by User 7754
-
Venn Diagram Help
written by kevluv93
-
How do I view the height / distribution of ChIP-seq peaks?
written by James Ashmore
-
Chip-Seq Peak Calling Using Macs In Galaxy- Overlapping Peaks For Replicates
written by anon111 •
-
P Value Or Statistical Significance Of Real Peak Compared To Random Peak Overlaps
written by biorepine
Software To Find Overlaps Of Chip-Seq Peaks In Multiple Samples.
This is a common misconception that even the most seasoned analysts make.
A peak isn't a single entity. Its a region of the genome. You shouldn't use Venn Diagrams (or any other categorical visualization) for something with a beginning, end, and ideally intensity. Check out the deeptools package to see how to overlay positional information for 1 or more categories of sample/feature.
I haven't tried myself but venneuler for R seems like what you need: http://ieeexplore.ieee.org/xpls/icp.jsp?arnumber=5728808&tag=1
I have the same question as you have.
I have a post here about how to do this with mergePeaks in HOMER and the
VennDiagrampackage in R, and also a post here about how to do the same with UpSet plots.