This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Is it possible to visualize ChIP peaks of three different samples in one plot using ChIPseeker?

Hi,

I used MACS2 to call the peaks of three different ChIP samples. I have 3-bed files which I want to visualize.

I am trying to generate "Average Profile of ChIP peaks binding to TSS region" plot from ChIPseeker. Though I am able to generate three different "Average Profile of ChIP peaks binding to TSS region" plots for each bed file but, I want to generate one figure with will be having peaks of 3 bed files. In short I wanted to merge my three plots in one PNG.

I want to generate something like this: https://drive.google.com/file/d/1uGscEzuqmASboDLABGI2TcDdIWw3-b5B/view?usp=sharing

Any kind of help will be great. looking forward to hearing back from you. Thanks in advance Swadha

chip-seq chipseeker bdg
• 4,081 views • link
updated by ZZzzzzhong ZZzzzzhong • written by s.singh • s.singh

well HOMER has everything what you need. merge peaks from all the files and calculate enrichment over merged files.

• 1 views • link
written by Prakash Prakash

Thanks for getting back. I will definitely look into HOMER's kit

• 1 views • link
written by s.singh • s.singh

1 answer

This is easy to implement by CHIPseeker. Let's say 1_peak,2_peak and 3_peak are your peak_file.

files <- list(1_peak,2_peak,3_peak)
peakAnnoList <- lapply(files, annotatePeak, TxDb=txdb,tssRegion=c(-3000, 3000), verbose=FALSE)
plotDistToTSS(peakAnnoList) + theme(plot.title = element_text(hjust = 0.5))

Then you will get the picture you want.

• 1 views • link
written by ZZzzzzhong ZZzzzzhong

Thanks for getting back. Yes, I got the desired output. I am able to plot all the graphs in one plot. I have one more small question. How do I change the legends? By default, it labels the graphs as "peak1" peak2" peak3" and so forth. How do change the labels?

• 1 views • link
written by s.singh • s.singh

if you want to label 1_peak as name1(Any name you want) you can do thatfiles <- list(name1=1_peak,name2=2_peak,name3=3_peak)

• 1 views • link
written by ZZzzzzhong ZZzzzzhong

It is working now! Thank you!

• 1 views • link
written by s.singh • s.singh

Log in to answer this question.

More posts like this