Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: OMA - without all vs all
1
gravatar for dtejadamartinez
3.1 years ago by
dtejadamartinez • 20
dtejadamartinez • 20 wrote:

I have a data set of 700 coding sequences (which I saved in a reference.fa file) and I need to obtain the orthology relationships (orthologous and paralogous) of these genes against 70 species. I have downloaded the genome cds from NCBI for each species. Now, is it possible to get only these comparisons without going through the all vs all phase? Could I start OMA without comparing all the genomes (cds) with each other? How do I launch these analyses ?

Thanks for your help,

orthology oma • 953 views
ADD COMMENT • link •
modified 3.1 years ago • written 3.1 years ago by dtejadamartinez • 20

Thanks for your help,

I think I would compute all vs all then.

ADD REPLY • link written 3.1 years ago by dtejadamartinez • 20
0
gravatar for Christophe Dessimoz
3.1 years ago by
Christophe Dessimoz • 700
University College London
Christophe Dessimoz • 700 wrote:

At the moment, it’s necessary to compute all-all among your 70 genomes as well to elucidate orthology relationships.

However, I see two fast alternatives:

1) If these 70 reference genomes are available in OMA, you could export them including all-against-all. Then, OMA would “just” need to compute the 700 CDS against the 70 genomes (and the CDS against themselves but that would be very fast.)

2) if you are mainly interested in functional annotation, you could use the new annotation feature of the OMA browser (https://omabrowser.org/oma/functions/). There is also a REST API function for that (https://omabrowser.org/api/docs#function)
ADD COMMENT • link modified 3.1 years ago • written 3.1 years ago by Christophe Dessimoz • 700

I am trying to run OMA (version OMA.2.4.2) on cluster using SLURM. I basically followed the procedure described in Oma Cheat Sheet, and in my DB directory I collected genomes for 20 species (coming from OMA using the export option) plus the genome of my species of interest. However, the array-job (splited into 1000 parallelized jobs) seems to take a lot of time (each job runs for over 8 hours!) Is there any way to speed up this step?

Thank you, Ilenia

ADD REPLY • link written 5 days ago by ilenia.urso • 0
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: 1013 users visited in the last hour