Question: [ChIPSeeker] Getting overlapping genomic feature annotation (e.g. peak locate in gene A & defined TSS region of gene B)
gravatar for chiefcat
4 months ago by
chiefcat70 wrote:

Hi guys,

I'm trying to use ChIPseeker but couldn't find a way to call the overlapping features during peak annotation.


TSS setting: -3kb +500bp

--> I would like to get the ChIP peak annotated with "Rrs1, 3'UTR" & "Adhfe1, TSS" I've tried a lot of different arguments but all of them just me ONLY "Rrs1, 3'UTR" enter image description here

Do anyone know how to do this?

And If a peak is overlapping with 2 different gene promoters, can ChIPseek annotate this?


ADD COMMENTlink modified 4 months ago by Guangchuang Yu1.7k • written 4 months ago by chiefcat70
gravatar for Guangchuang Yu
4 months ago by
Guangchuang Yu1.7k
China/Hong Kong/The University of Hong Kong
Guangchuang Yu1.7k wrote:

there is a genomicAnnotationPriority parameter documented in You can select to report "Rrs1, 3'UTR" OR "Adhfe1, TSS", but and "Rrs1, 3'UTR" & "Adhfe1, TSS" as annotatePeak only report the 1st hit according to the annotation priority.

Maybe setting the parameter addFlankGeneInfo = TRUE, flankDistance = 1000 will give you some useful information. The above parameters will report a new column of flanking genes with all the genes at the flanking region of the peak (e.g. [-1000, 1000] if flankDistance = 1000).

ADD COMMENTlink written 4 months ago by Guangchuang Yu1.7k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1671 users visited in the last hour