How to get the list of genes from DiffBind dba.report or dba.peakset
3
0
Entering edit mode
Bordiya ▴ 10
@bordiya-16743
Last seen 6.2 years ago

Hi all,

I am analyzing H2AZ ChIP seq data. I have the dba.report and dba.peakset which give chromosome number and peak position of differential peaks. I would like to get list of genes corresponding to those peaks. How can I do this?

Thank you!

diffbind • 2.9k views
ADD COMMENT
1
Entering edit mode
Rory Stark ★ 5.2k
@rory-stark-5741
Last seen 5 weeks ago
Cambridge, UK

There are a number of annotation packages available that can work on the GRanges object returned by dba.report() or a bed file. Within Bioconductor, ChIPpeakAnno is popular -- many DiffBind users feed the results into ChIPpeakAnno.

ADD COMMENT
0
Entering edit mode

Thanks a lot Rory!

ADD REPLY
1
Entering edit mode
Bordiya ▴ 10
@bordiya-16743
Last seen 6.2 years ago

I used ChIPpeakAnno in combination with other packages, please find the script in the link below. It works great!!

 

 

 

ADD COMMENT
1
Entering edit mode
mdidish ▴ 10
@cbee3b30
Last seen 3 months ago
France

Three years after, if someone has the same question, here is a proposition with ChIPseeker:

library(ChIPseeker)
library(TxDb.Mmusculus.UCSC.mm10.knownGene)
library(org.Mm.eg.db)

peak.counted <- dba.count(peak.dba, summits=T, score = DBA_SCORE_SUMMIT_POS)
# retrieve summit list (GRanges object)
summits <- dba.peakset(kdm6b.counted, bRetrieve=TRUE)
# annotate with ChIPseeker
peak.anno = annotatePeak(
  peak = summits,
  tssRegion = c(-1000, 1000),
  TxDb = TxDb.Mmusculus.UCSC.mm10.knownGene,
  annoDb = "org.Mm.eg.db"
)
plotAnnoPie(peak.anno)
dfPA = as.data.frame(peak.anno)
ADD COMMENT

Login before adding your answer.

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