Question: Chip-seq peak annotation [Chipseeker]
1
gravatar for anu014
2.4 years ago by
anu014170
India
anu014170 wrote:

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

ADD COMMENTlink modified 2.3 years ago by Guangchuang Yu2.2k • written 2.4 years ago by anu014170
2

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 REPLYlink written 2.4 years ago by Sinji2.9k
2
gravatar for Guangchuang Yu
2.3 years ago by
Guangchuang Yu2.2k
China/Guangzhou/Southern Medical University
Guangchuang Yu2.2k wrote:

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

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by Guangchuang Yu2.2k
1

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 REPLYlink modified 2.3 years ago • written 2.3 years ago by anu014170
2

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 REPLYlink modified 2.3 years ago • written 2.3 years ago by Guangchuang Yu2.2k

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

ADD REPLYlink written 19 months ago by anu014170
1

yes. mouse is similar to huaman.

ADD REPLYlink written 19 months ago by Guangchuang Yu2.2k

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

ADD REPLYlink written 19 months ago by anu014170
Please log in to add an answer.

Help
Access

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