How to use peak summits instead of TSS or transcript to plot heatmap for other factors
1
0
Entering edit mode
ret28 • 0
@user-24024
Last seen 4.0 years ago

Hi,

I am relatively new to R and analysing ChIPseq data sets and have a query.

I have used MACS2 callpeak for IP and input samples to find significantly enriched peaks for 2 proteins. I have the resulting summits.bed and .narrowPeak files. I can use these with e.g. ChIPseeker getPromoters and getTagMatrix to produce a heatmap or profile following the instructions to show TSS regions and read count frequency.

My question is how can i use the enriched peaks from one dataset as the windows for another set ? e.g. to see if a protein has a number of reads mapped within 3000 bases upstream or downstream of the peak centre from another protein.

Thanks

Rob

ChIPseqR ChIPs ChIPseeker ChIPSeq chipseq • 1.8k views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 1 day ago
United States

Generally speaking, you can read the bed files into R using rtracklayer, which will create a GRanges object. You could then use resize to make the GRanges object to increase the width of the peaks for the first protein, and then use findOverlaps to see what (expanded) peaks from protein 1 overlap the unexpanded peaks of protein 2.

You should read the vignettes for GenomicFeatures and rtracklayer for more background.

ADD COMMENT

Login before adding your answer.

Traffic: 845 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6