ChIPpeakAnno getAnnotation request fail
0
1
Entering edit mode
igor ▴ 40
@igor
Last seen 8 months ago
United States

I am trying to use BioMart. Selecting a database works:

mart = useMart(host="useast.ensembl.org", biomart="ENSEMBL_MART_ENSEMBL", dataset="mmusculus_gene_ensembl", verbose=T)

If I try to get annotations using TSS featureType, it takes a few seconds, but works:

martAnnot = getAnnotation(mart=mart, featureType="TSS")

However, if I use Exon:

martAnnot = getAnnotation(mart=mart, featureType="Exon")

It fails after several minutes:

Error in value[[3L]](cond) : 

  Request to BioMart web service failed. Verify if you are still connected to the internet.  Alternatively the BioMart web service is temporarily down.

If the TSS works, then the connection and the service are both functional. I tried repeating both several times with same results. What is going on?

chippeakanno biomart ensembl • 955 views
ADD COMMENT
1
Entering edit mode

It looks like you lost the internet connection or your connection is not stable.

Maybe you want to have a try to use local annotations like this:

library(ChIPpeakAnno)

library(EnsDb.Mmusculus.v79)

exonAnno <- annoGR(EnsDb.Mmusculus.v79, feature="exon")

 

ADD REPLY
0
Entering edit mode

If it's the internet connection, then why do I not have problems with TSS? I repeated this multiple times over the course of several days.

Using EnsDb* packages limits me to just the three genomes that are available and just the build versions that are available.

ADD REPLY
1
Entering edit mode

I got it. It is a timeout error. It will be fix in next release version (2 or 3 days later). Thank you for reporting that.

ADD REPLY
1
Entering edit mode

Actually, it's pretty straightforward and easy to build EnsDb objects on any species/Ensembl release based on GTF files fetched from AnnotationHub or downloaded directly from the Ensembl ftp server:

library(ensembldb)
library(AnnotationHub)
## load AnnotationHub
ah <- AnnotationHub()

## Check what GTF files are available for Ensembl release 83, Mus musculus

> query(ah, c("Mus_musculus", "release-81"))
AnnotationHub with 7 records
# snapshotDate(): 2016-01-14
# $dataprovider: Ensembl
# $species: Mus musculus
# $rdataclass: FaFile, GRanges
# additional mcols(): taxonomyid, genome, description, tags, sourceurl,
#   sourcetype
# retrieve records with, e.g., 'object[["AH47973"]]'

            title                                 
  AH47973 | Mus_musculus.GRCm38.81.gtf            
  AH49243 | Mus_musculus.GRCm38.cdna.all.fa       
  AH49244 | Mus_musculus.GRCm38.dna_rm.toplevel.fa
  AH49245 | Mus_musculus.GRCm38.dna_sm.toplevel.fa
  AH49246 | Mus_musculus.GRCm38.dna.toplevel.fa   
  AH49247 | Mus_musculus.GRCm38.ncrna.fa          
  AH49248 | Mus_musculus.GRCm38.pep.all.fa     

## Fetch the GTF
> gtf <- ah[["AH47973"]]

## And build the EnsDb database.
> dbfile <- ensDbFromGRanges(gtf, organism="Mus_musculus",
 genomeVersion="GRCm38", version="81")
## Load the database.
> edb <- EnsDb(dbfile)

## you can use that straight away or build a package from it with makeEnsembldbPackage

## Note: in the devel version you can directly build the EnsDb from AnnotationHub; check the vignette for more information.
ADD REPLY
0
Entering edit mode

Great!. I will add this into development version of ChIPpeakAnno to build EnsDb or TxDb.

ADD REPLY

Login before adding your answer.

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