Hi everyone,

I am a newbie at gene ontology analysis in R. I have a list of Uniprot ID from different organisms. These ids were found as a result of a blast search. I want to find GO term hierarchy of these proteins for Biological process category. Then I want to sum up all categories from leaves to root with number of proteins in corresponding category. In the end I want to have something like:

GO:0008150 (7 proteins)

-> Child 1 (7 proteins)

      -> Child 1.1 (5 proteins)
      -> Child 1.2 (2 proteins)

I found that GO.db library returns ancestors, children. However, they are not ordered by hierarchy. Uniprot IDs can be queried via UniProt.ws package. But I haven't tried it yet.

If you can help me, I would be appreciated.

Regards,

Canan