Chip-seq peak annotation [Chipseeker]
1
1
Entering edit mode
4.1 years ago
anu014 ▴ 180

Hi all!

I was dealing with this package called 'Chipseeker' . I tried annotating a chip-seq peak calls fro Mycobacterium tuberculosis. As there is no 'org.Mtu.eg.db' package for annotation & 'TxDB' can be used in Chipseeker, so I created TxDB package & annotated as follows :

library(GenomicFeatures) txdb <- makeTxDbFromGFF("./m.tub1.gff3",circ_seqs=character(), organism="Mycobacterium tuberculosis") peakAnno <- annotatePeak(peak1, tssRegion=c(-3000, 3000), TxDb=txdb, genomicAnnotationPriority = c("Promoter", "5UTR", "3UTR", "Exon", "Intron", "Downstream", "Intergenic"))

Note: GFF was downloaded from ftp://ftp.ensemblgenomes.org/pub/bacteria/release-35/gff3/bacteria_0_collection/mycobacterium_tuberculosis_h37rv/Mycobacterium_tuberculosis_h37rv.ASM19595v2.35.gff3.gz

The chip-seq was done for 'Gyrase' but annotations from Chipseeker gave mostly 'promoters' . Is is correct or am I suppose to add some more steps? Please help me.

Thanking you in advance, Anupriya

ChIP-Seq chipseeker peak annotation • 2.3k views
ADD COMMENT
2
Entering edit mode

Check your ChIP-seq data on IGV and look at a couple of examples of peak calls. Do they mostly fall in the promoter region? If so, then the annotation is probably correct.

ADD REPLY
2
Entering edit mode
4.1 years ago
Guangchuang Yu ★ 2.5k

For virus/bacteria, tssRegion=c(-3000, 3000) is too large.

ADD COMMENT
1
Entering edit mode

Thank you Guangchuang for the reply :) So how can I optimise the parameter for 'bacteria'? And do parameter 'genomicAnnotationPriority' affects the annotation process in a way that Chipseeker will try to fit promoters first then rest of the types of annotation?

ADD REPLY
2
Entering edit mode

I don't know how to optimise the tssRegion parameter. c(-1000, 1000) or c(-3000, 3000) are commonly used for human. You may need to checkout some references for your species.

Yes, genomicAnnotationPriority parameter will affect the priority of genomic annotation and you can set lower priority of promoter annotation.

ADD REPLY
0
Entering edit mode

Can I use the same range for Mouse too? as it's almost similar to human..

ADD REPLY
1
Entering edit mode

yes. mouse is similar to huaman.

ADD REPLY
0
Entering edit mode

Okay. I'll be using same range for mouse now. Thank you Guangchuang.

ADD REPLY

Login before adding your answer.

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