library(ChIPseeker)
library(TxDb.Hsapiens.UCSC.hg38.knownGene)
library(org.Hs.eg.db)
library(clusterProfiler)
options(ChIPseeker.ignore_1st_exon = T)
options(ChIPseeker.ignore_1st_intron = T)
options(ChIPseeker.ignore_downstream = T)
options(ChIPseeker.ignore_promoter_subcategory = T)
peak_regions <- annotatePeak(SeuratObj_ATAC@assays$ATAC@ranges, level = 'gene',tssRegion = c(-1000, 1000),
TxDb=TxDb.Hsapiens.UCSC.hg38.knownGene, annoDb = "org.Hs.eg.db")
head(peak_regions@anno)
plotAnnoPie(peak_regions,legend.position = 'top')
image.png
