Hello how to use ChIPpeakAnno for latest version of ensemble genome aligned data
1
0
Entering edit mode
Lucky • 0
@lucky-24920
Last seen 7 days ago
India

Hello

I did the alignment with the my Chip data with latest version of mouse genome from ensemble Mouse ensemble ref, how can use ChIPAnno package for this genome version, in the R work flow we have load library of Ensdb as mentioned. If google for mouse db I am getting latest version of library for mouse is EnsDb.Mmusculus.v79 is it fine to use this or is there any way to get the Ensdb for latest version of the genome

Thank you

ChIPpeakAnno • 1.6k views
ADD COMMENT
2
Entering edit mode
Kai Hu ▴ 70
@kai
Last seen 7 months ago
Worcester

The version of the annotation db must match with the genome used for mapping because the coordinates may differ for different genome releases. Therefore, if you used Mus_musculus.v103 for mapping, you'd best also use EnsDb.Mmusculus.v103 for annotation.

Below are 3 options to fetch the annotation files:

With AnnotationHub

library(AnnotationHub)

ah <- AnnotationHub()
EnsDb.Mmusculus <- query(ah, pattern = c("Mus musculus", "EnsDb"))
EnsDb.Mmusculus
EnsDb.Mmusculus.v101 <- EnsDb.Mmusculus[[length(EnsDb.Mmusculus)]]
class(EnsDb.Mmusculus.v101)

Note that AnnotationHub is not always up-to-date (it has v101 as the latest version as of 2021/03/06).

With ChIPpeakAnno::getAnnotation()

library(ChIPpeakAnno)
library(biomaRt)

listMarts()
mart <- useMart(biomart="ENSEMBL_MART_ENSEMBL", dataset="mmusculus_gene_ensembl")
Annotation <- getAnnotation(mart)

Note that getAnnotation() leverages biomart, which is always up-to-date.

Build your own EnsDb package (need extra effort)

There is a good tutorial on how to use makeEnsembldbPackage to generate EnsDb. It calls the Ensembl Perl API internally, and it may need some extra effort to set up the Perl environment.

ADD COMMENT
0
Entering edit mode

Okay Thank you for the reply I will Try above mentioned methods.

ADD REPLY

Login before adding your answer.

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