Help generating multiple annotations under one H3K27me3 peak
Entering edit mode
9 weeks ago

Hi everyone, first post here!

I'm currently working on a ChIPseq data set enriching for H3K27me3 between four Arabidopsis genotypes.

My current problem is at the annotation step using annotatepeak from ChIPseeker. Because the peaks I have called (using MACS2) are broad, they cover multiple loci in the genome, but only become annotated with the closest TSS instead of say, three different loci. As a result, identifying differentially methylated regions between my WT and mutants is quite error prone when I look at the bed files in IGV. I've tried to show an example in the image attached.

In this example, the peak covering AT2G19000 is showing as not methylated in my global annotation for Mutant 1 because it has mapped only to the next gene across, AT2G19010

This is the line of code I use in R:

WT_consensus_anno <- annotatePeak(WT_covered_granges, tssRegion=c(-3000, 3000),
                      TxDb=txdb, annoDb="org.At.tair.db")

Where txdb is extracted from the TxDb.Athaliana.BioMart.plantsmart28 package, and WT_covered_granges is a GRanges object containing the called peaks from MACS2.

Any help would be much appreciated, I've been puzzling over this for a few days and searched everywhere for answers.

If I've left out any information please let me know

Thanks a lot!

ChIP-Seq R chipseeker • 188 views

Login before adding your answer.

Traffic: 1066 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6