I want to assign GO slim terms to some plant genes in non-model plants. So far as I know the only plant genomes relatively accurately annotated with GO slims are rice and Arabidopsis. So I am thinking of using best blast hit in rice or Arabidopsis (depends on whether the genome under study is monocot or dicot) to assign the GO slims. I'm working on genome scale dataset so cannot do manual curations, but I want be set up a reasonably accurate procedure. I'm new to this area and would appreciate if somebody can let me know if there are better ways to do this, or things other than blast similarity that I should take into account.
As you mention, you want to infer isofunctionality from homology, you you can use the homology hits to transfer rice and Arabidopsis annotations to your new species. An important caveat here is that if the gene has duplicated after speciacion, which is common in plants, your isofunctionality prediction may be less accurate, since plants tend to sub- and neofunctionalize their gene repertoires after duplication. So try to annotate unique one-to-one hits separately from one-top-many hits.
Hi Jingping, you could try the following two methods:
http://www.blast2go.com/b2ghome This method assigns function by blast search.
With this method, you first do a genome-wide InterPro domain scan, then assign GO terms based on InterPro2GO mapping.
One additional method is to first build family trees (e.g. http://www.phytozome.net/; http://www.pantherdb.org/genes/ ), then transfer annotations within the protein family. However this method involves manual curation so it might not work for you.
DHL On behalf of the Gene Ontology (GO) Help Desk
You might also be interested in the electronic GO annotations that are produced by the projection of manual annotations from Arabidopsis and rice to other plant species using Ensembl Compara ortholog data. The species that have been annotated using this method include poplar, grape, brachypodium, maize etc. and there are currently almost 270,000 annotations for these species.
The EBI's QuickGO browser allows you to view these annotations using this link; http://www.ebi.ac.uk/QuickGO/GAnnotation?a=&termUse=ancestor&relType=IPO%3D&customRelType=IPOR%2B-%3F%3D&ref=GO_REF%3A0000035
and also you will be able to use QuickGO (home page; www.ebi.ac.uk/QuickGO) to view these annotations in a GO slim, either a pre-defined one or one you make yourself.
The annotation method used to create these annotations is GO_REF:0000035, which is described on this page; http://www.geneontology.org/cgi-bin/references.cgi
If you need any help in using QuickGO to do your slim, please contact email@example.com
Rachael on behalf of the UniProt-GO annotation project.