Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Getting the highest taxonomy in the HOG database for a given OMA gene id
1
gravatar for yogihooda88
16 months ago by
yogihooda88 • 10
yogihooda88 • 10 wrote:

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
hog oma orthologs • 384 views
ADD COMMENT • link •
modified 16 months ago by Christophe Dessimoz • 660 • written 16 months ago by yogihooda88 • 10
-1
gravatar for Christophe Dessimoz
16 months ago by
Christophe Dessimoz • 660
University College London
Christophe Dessimoz • 660 wrote:

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 • link written 16 months ago by Christophe Dessimoz • 660
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 787 users visited in the last hour