This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Clusters proteins at each taxonomic rank - OMA

Hello everybody,

Anyone knows how to clusters proteins ( for example 90% of sequence identity) at each taxonomic rank (for example phylum) from an orthoXML file downloaded from OMA or other files. This will create non-redundant proteins in all taxomomic rank, but will avoid to loose a phylum that contain sequences.

Thanks! Alejandro

oma clustering
• 1,388 views • link
updated by Adrian Altenhoff Adrian Altenhoff • written by anusblat • anusblat

1 answer

Hi Alejandro,

can you elaborate what you're trying to do? I'm one of the OMA developers, but I have a hard time to understand what you're doing.

In brief, OMA HOGs are groups of genes that all started diverging from a certain ancestral gene through speciation, so that gives you already gene clusters. If you want to subcluster them according to sequence identity, you would have to do that based on the available protein sequences. you could use pyham (https://github.com/DessimozLab/pyham) to extract the extent genes of any cluster you're interested.

• 0 views • link
written by Adrian Altenhoff Adrian Altenhoff

Log in to answer this question.

More posts like this