Getting the highest taxonomy in the HOG database for a given OMA gene id
1
1
Entering edit mode
3.8 years ago
yogihooda88 ▴ 10

I am trying to find how conserved a given gene is using the OMA database. I have a list of 500 genes and for each gene, I want to identify the last ancestor that has the HOG family a given gene belongs to? Does anybody know how do I go about doing that?

  • Yogesh
OMA orthologs HOG • 865 views
ADD COMMENT
-1
Entering edit mode
3.8 years ago

There are several options. One way is to use the REST API. If you use e.g. P53_RAT as an input sequence:

https://omabrowser.org/api/hog/P53_RAT/?level=root

You can see that the level is Chordata (in the current Dec 2018 release).

If you use Python, you can use the omadb package

$ pip install omadb

import omadb
from omadb import Client 
c = Client() 
c.hogs[c.hogs['P53_RAT'].roothog_id].level

Alternatively, still in Python, you could use the pyHam library for this (https://github.com/DessimozLab/pyham), and load the HOGs dynamically from the OMA browser (or from an OrthoXML file):

import pyham

my_gene_query = 'P53_RAT' 
pyham_analysis = pyham.Ham(query_database=my_gene_query, use_data_from='oma')
pyham_analysis.get_list_top_level_hogs()[0].genome.name

If you use the R, you can use our package OmaDB:

install.packages('devtools')
library(devtools)
install_github('dessimozlab/omadb')

library(OmaDB)
getHOG(id = 'P53_RAT', level='root')$level
ADD COMMENT

Login before adding your answer.

Traffic: 1352 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6