Dear community,
I have gene expression data from several animals, which I'd like to compare. These data have UBERON identifiers for the tissues in which they're expressed, however, many tissues belong to different (yet, potentially) comparable anatomical entities (e.g. "kidney" and "adult mammalian kidney"). I know that UBERON, just like GO, has a hierarchical structure, so I thought that perhaps I could compare tissues belonging to the same group withinthe UBERON hierarchy. A way to access such data is via OLS, yet I haven't found an easy way to retrieve the UBERON relationships among multiple anatomical entities. I also investigated R package "rols", but I read that it might not be suitable for such task.
So my question is: Is it a easy-ish way to extract the hierarchical relationships between UBERON terms (from a list of UBERON ids, for example)?
Thanks a lot.
Best,
Hi Guille
It sounds like you are looking for a programmatic interface, and you might prefer R. The last time I looked I found a number of R generic ontology libraries were problematic. But as I recall there is a bgee library specifically for working with expression data annotated with uberon, you might find this useful.
We have a python library that can be used with Uberon
http://ontobio.readthedocs.org/
(you can use locally downloaded files, or it will transparently wrap a sparql service)
There are a number of REST endpoints that can be used as well. E.g. the OLS API, or the ontobee SPARQL endpoint.
How exactly do you intend to compare? Are you doing this for pairs of genes?
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot_______________________________________________
Obo-anatomy mailing list
Obo-a...@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/obo-anatomy
There are a number of ways to do this but I think I still need a bit more information.
For any pair of genes, you could get the bgee expression values for all classes, and then report values for shared superclasses. E.g. if G1 is expressed in the cerebellum, and G2 in the brain, you can roll up the the level of brain. This part should be straightforward using an ontology library or a pre-calculated table of ancestry relationships