In this practical we will read in two replicate sets of peaks from the Myc Encode dataset and annotate these to gene using the ChIPseeker packge. We will visualise the results using ChIPseeker plotting functions and see how to combine samples into one plot.
In the second half of the practicals we will perform an enrichment analysis on the genes overlapping peak sets using the GOseq package.
First lets read in the two Myc Mel peakset replicates as we did for our previous exercise.
Now annotate the two peaksets to mm9 UCSC genes and draw the distribution of peaks in genomic features and distance of peaks to TSS regions.
ChIPseeker allows the user to plot multiple samples in one plot. Combine the csAnno objects for each sample into a list and name elements with their respective sample names. Use the list of csAnno objects with plotAnnoBar function.
Retrieve all unique names (as SYMBOLS) of genes with peaks in their TSS/promoter for replicate 1 and 2.
Create an object containing of Myc peaks in Mel replicate 1 overlapping Myc peaks in Mel replicate 2.