Goseq package in R:getting terms after using getgo

Question

[GabrielMontenegro]

I have a simple gene vector and trying to get the GO terms using the goseq package from Bioconductor in R.

In this example just one gene:

> unique(d$genes[1])
> [1] "SMARCC1"

> head(stack(getgo(unique(d$genes[1]),'hg19','geneSymbol',fetch.cats=c("GO:BP"))))

      values     ind
1 GO:0006139 SMARCC1
2 GO:0006325 SMARCC1
3 GO:0006333 SMARCC1
4 GO:0006337 SMARCC1
5 GO:0006338 SMARCC1
6 GO:0006351 SMARCC1

Now that I have the GO numbers, is there a way to get the GO terms? That is the exact biological process?

I saw this post which somehow had the GO terms after for each GO value, but I couldn't figure out where he got them.

Answer 1

[Parham]

I have an annotation file for the organism that I work with that holds GO numbers and their corresponding name. If you have something like that it should be easy to extract terms from that. Otherwise I guess you should look for such an annotation file. Probably it is available on Gene Ontology Consortium.

[GabrielMontenegro]

I was checking the GO consortium website, but couldn't find a table with that information. Would you mind showing me where I can get this info?

Answer 2

[Benn]

In the newest version of goseq, you get a table after using the goseq function which conatins the term column. You can use these to match with you head-stack table.

[GabrielMontenegro]

I just downloaded the latest version, but after running goget I get a list including the genes and the GO terms for every gene. Could you show me the exact function?

[Benn]

I was talking about the goseq function in the package, this is to do the enrichment analysis.

For example:

pwf=nullp(gene.vector,"hg19","geneSymbol")

GO.wall=goseq(pwf,"hg19","geneSymbol",use_genes_without_cat=TRUE)

The object GO.wall is a data frame which contains a column term. You can merge this further in R with your head-stack-getgo object.