Question: Is it possible to visualize ChIP peaks of three different samples in one plot using ChIPseeker?
0
gravatar for swadha
10 weeks ago by
swadha20
University of California, United States
swadha20 wrote:

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

chipseeker chip-seq bdg • 243 views
ADD COMMENTlink modified 9 weeks ago by ZZzzzzhong150 • written 10 weeks ago by swadha20
1

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

ADD REPLYlink modified 10 weeks ago • written 10 weeks ago by prakash530

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

ADD REPLYlink written 10 weeks ago by swadha20
3
gravatar for ZZzzzzhong
9 weeks ago by
ZZzzzzhong150
ZZzzzzhong150 wrote:

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.

ADD COMMENTlink written 9 weeks ago by ZZzzzzhong150

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?

ADD REPLYlink written 9 weeks ago by swadha20
1

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)

ADD REPLYlink written 9 weeks ago by ZZzzzzhong150

It is working now! Thank you!

ADD REPLYlink written 8 weeks ago by swadha20
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: 568 users visited in the last hour