GO enrichment analysis using non model organism

Question

[dago]

Hi everybody,

I have proteome data for a non model organism, and I would like to perform gene onthology (GO) enrichment analysis.

I would like to perform the enrichment for the up-regualted proteins using the all proteins encoded in the genome as background dataset. I downloaded the GO annotation for each protein directly from Uniprot and now I am kind of lost among the many tools available for doing the analysis.

I am looking for web tools or R-Bioconductor tools more appropriate form my purpose, because it seemed to my that most of the tools deal with model organism.

Any advice will be greatly appreciated.

[Renesh]

What is the organism you are working on?

Answer 1

[dago]

An update.

I used the Bioconductor package topGO and I think it is a great tool for performing enrichment for non model organisms.

[Ric]

Could please you share you code how you did it with topGO?

Answer 2

[Renesh]

If you are working on plant, use agriGO enrichment tool. It has most of the species listed. It supports gene IDs for performing GO enrichment. You can also try BinGO plugin of cytoscape.

If you organism is not listed in their default database, you can perform blast homology search of your proteins with closely related organims available in their database. Use these top hits resulted from blast for performing Go enrichment.

[dago]

I work with bacteria and I found the .goa file for the GO annotation on the Uniprot FTP server

[Renesh]

Use BinGO

You don't need to have GO IDs in BinGO. You should have gene IDs listed in their database. Check for your bacteria, if it is available there. If not then do as i suggested in above post.

[dago]

Indeed, BINGO seems to be quite interesting.

I read on its website that it is possible to perform analysis on Custom dataset (http://www.psb.ugent.be/cbd/papers/BiNGO/Customize.html), which would be what I need.

Do you have any idea how I can create my own ontology file or download an ontology file useful for bacteria?

Thanks again for your help