[BioC] GOseq with ce6

Nadia Davidson nadia.davidson at mcri.edu.au
Sat Apr 12 03:31:42 CEST 2014

François Lefebvre <lefebvrf at ...> writes:

> Hi all,
> Our RNA-seq pipeline uses GOseq, and we are thankful to the authors for
> offering this great tools.
> GOseq crashed on us today on C. elegans data with the error":
>  "Couldn't grab GO categories automatically.  Please manually specify."
> The workaround is to use the non-native mechanism.
> After investigating, the error is thrown after the function getgo() cannot
> reliably identify the org package associated to the Ce6 genome.There is a
> collision with the yeast org package.
> So for fun I listed all other genomes for which that would not work:
> library(goseq,quietly = TRUE)
> genomes = sort(supportedGenomes()$db)
> supported = sapply(genomes,function(genome)
> {
>     orgstring = as.character(.ORG_PACKAGES[grep(gsub("[0-9]+", "", genome),
> names(.ORG_PACKAGES), ignore.case = TRUE)])
>        if (length(orgstring) != 1) {
>     #stop("Couldn't grab GO categories automatically.  Please manually
> specify.")
>     return(FALSE)
>        }else{
>     return(TRUE)
>     }
> })
> print(names(supported[!supported]))
> I just thought it was odd for the package not to work out of the box for
> the model organism C. elegans. It also wouldn't work for pig, and maybe
> more.
> Thank you!

Dear François,

Thank you for reporting this. You've found a bug in the getgo function. 
For c. elegans, two .ORG_PACKAGES were being returned because 
"sacCer" also has the string "ce" within it and grep was being used 
to match strings. I've fixed this in the development version of goseq. 

For pig, the issue was different, but related and has also been fixed. 
Many of the other genomes you listed above are still not supported, 
but this has to do with what organisms have "org.*" annotation 
packages in bioconductor. "supportedGenomes" is perhaps a 
bit misleading in this way, as it's just the gene lengths being 
automatically looked up that is supported.

Please let us know if you find any more issues.


More information about the Bioconductor mailing list