Question: ChIPseeker peak annotation tssRegion
gravatar for lifeng.liu.biology
7 months ago by
lifeng.liu.biology10 wrote:

I was confused about the argument of tssRegion in peak annotation function of ChIPseeker. It seems that there is no change of the annotation output when I changed the tssRegion setting, from tssRegion=c(-3000, 3000) to tssRegion=c(-2000, 0). I am working on a plant species, Brachypodium distachyon (~272 Mb). I thought this is to define the promoter region, which supposed to be upstream of TSS or 5UTR and setting upstream=3000, downstream=3000 is just for plotting around TSS. So, what is the point of this argument here?

And, does anyone know how the promoter was defined in genomic annotation categories (promoter, 5/3 UTR, exon, intron, downstream, intergenic)? For downstream, does it start after 3' UTR? Can I manually set the downstream range, e.g. 500bp?

peakAnno <- annotatePeak(files[[4]], tssRegion=c(-3000, 3000), TxDb=txdb)


ADD COMMENTlink modified 7 months ago • written 7 months ago by lifeng.liu.biology10
gravatar for Guangchuang Yu
7 months ago by
Hong Kong
Guangchuang Yu990 wrote:

do you have any reproducible example of change the `tssRegion` setting has no effect of the result?

ADD COMMENTlink written 7 months ago by Guangchuang Yu990

R script used here:




txdb <- loadDb("Bdistachyon_v314.sqlite")
peak_H3K27me3 <- readPeakFile(file_H3K27me3)

peakAnno <- annotatePeak(peak_H3K27me3, tssRegion=c(-3000, 3000), TxDb=txdb, genomicAnnotationPriority = c( "5UTR",  "Exon", "Intron", "3UTR", "Promoter", "Downstream","Intergenic"))


upsetplot(peakAnno,"grey23", matrix.color = "black", = "grey23")

write.csv(peakAnno, "H3K27me3_all_peaks_annotation_U3000_D3000.csv")

peakAnno <- annotatePeak(peak_H3K27me3, tssRegion=c(-2000, 0), TxDb=txdb, genomicAnnotationPriority = c( "5UTR",  "Exon", "Intron", "3UTR", "Promoter", "Downstream","Intergenic"))


upsetplot(peakAnno,"grey23", matrix.color = "black", = "grey23")

write.csv(peakAnno, "H3K27me3_all_peaks_annotation_U2000_D0.csv")

lifengl-dm:scripts lifengliu$ diff H3K27me3_all_peaks_annotation_U2000_D0.csv H3K27me3_all_peaks_annotation_U3000_D3000.csv

lifengl-dm:scripts lifengliu$

And the numbers in AnnoPie or AnnoUpset are exactly the same. 

so the tssRegion argument here is to define the promoter region or something else? If you need the bed file and txdb file, I can send them to you.

ADD REPLYlink modified 7 months ago • written 7 months ago by lifeng.liu.biology10
gravatar for lifeng.liu.biology
7 months ago by
lifeng.liu.biology10 wrote:

Figured it out. The problem is the priority of the genomic annotation. Using the priority of ​genomicAnnotationPriority = c( "5UTR",  "Exon", "Intron", "3UTR", "Promoter", "Downstream","Intergenic") gives me same result with tssRegion=c(-3000, 3000)​ or tssRegion=c(-2000, 0)​. But if I change it to genomicAnnotationPriority = c( "Promoter", "5UTR",  "Exon", "Intron", "3UTR",  "Downstream","Intergenic"), it gives me quite different annotation and the annotation (genomicAnnotationPriority = c( "Promoter", "5UTR",  "Exon", "Intron", "3UTR",  "Downstream","Intergenic"), tssRegion=c(-2000, 0)​​) makes more sense in my case. 

ADD COMMENTlink written 7 months ago by lifeng.liu.biology10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 421 users visited in the last hour