This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to get output of annotated genes from ChIPseeker tool

Hi,

I am new at ChIP-Seq analysis. I used ChIPseeker package through R but I couldn't figure out how to get output of annoteted genes. Please help me!

Thank you.

chipseeker peakannotation
• 6,145 views • link
updated by kashiff007 kashiff007 • written by MS • MS

You can turn it into a data.frame using the generic function as.data.frame.

• 0 views • link
written by rpolicastro rpolicastro

Thank you! Is there any option in ChIPseeker to write overlapped genes?

• 0 views • link
written by MS • MS

Your question is unclear, but take a look at addFlankGeneInfo parameter

• 0 views • link
written by Aspire Aspire

Hi,

I already used that option as following:

peakAnnoList <- lapply(samplefiles, annotatePeak, TxDb=txdb, 
                       tssRegion=c(-1000, 1000), verbose=FALSE)

when I draw plots using;

genes= lapply(peakAnnoList, function(i) as.data.frame(i)$geneId)
vennplot(genes)

vennplot.peakfile(samplefiles)

vennplot(genes) finds 611 common genes between samplefiles and vennplot.peakfile(samplefiles) finds 674 peaks.

But when I check the output files I can not find same number common genes and peaks.

So that I want to write output of this common genes and peaks used at vennplot.

• 0 views • link
written by MS • MS

Hi EES,

You can use intersect() function to get overlapping genes/peaks.

e.g.

intersect(genes$A, genes$B)
• 0 views • link
written by gglim • gglim

2 answers

Annotating your peak and saving into your_peakAnno varable:

your_peakAnno <- annotatePeak(files[[7]], tssRegion=c(-2000, 0), TxDb=txdb)

Writing your_peakAnno table into your file:

write.table(as.data.frame(your_peakAnno@anno), file="your_peakAnno.txt", quote = F, row.names = F, sep = "\t")
• 0 views • link
written by kashiff007 kashiff007

Get annotation data frame

nanog_annot <- as.data.frame(peakAnnoList[["Nanog"]]@anno)
pou5f1_annot <- as.data.frame(peakAnnoList[["Pou5f1"]]@anno)

https://hbctraining.github.io/In-depth-NGS-Data-Analysis-Course/sessionV/lessons/12_annotation_functional_analysis.html

• 0 views • link
written by CB • CB

Log in to answer this question.

More posts like this