ChIPseeker: Annotating Peaks with both upstream and downstream gene
0
0
Entering edit mode
ikumar2000 ▴ 10
@d32746b8
Last seen 2.1 years ago
United States

I am trying to annotate ChIPseq/ATAC peaks with ChIPseeker (AnnotatePeak) however it only reports the nearest gene (either upstream or downstream of the peak whichever is the closest). Is there any way to allow AnnotatePeak to get nearest genes in both direction i.e. two genes assigned to each peak ?

ChIPseeker • 1.7k views
ADD COMMENT
0
Entering edit mode

Take a look at the addFlankGeneInfo parameter. If TRUE, ChIPseeker returns all genes that overlap a flanking region around the peak. I am not 100% sure how overlapping is precisely defined (I think is an overlap of any part of the gene, not just of the TSS).

When you say that you are interested in the "two nearest genes", you do not specify if you are interested in the TSS, or just any part of a gene. For example, if a p[eak is close to the 3' UTR of one gene, and to the TSS of another - are you interested in both? If so, then addFlankGeneInfo might do the job.

ADD REPLY

Login before adding your answer.

Traffic: 718 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