代码库12-ChIPseeker注释基因组区域

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')