Agora-generic pipeline seems to be “stuck” for several days after adding 1 more genome about agora

[diekei]

Hi @alouis72,

Many thanks for the super swift response! and the two possible suggestions.
We're mostly interested on reconstructing ALG for the oldest ancestors, so I don't think I can use the first solution approach.
I've tried to run it using the modified agora-generic.py scripts as suggested, and it worked!

Thank you!!
Best,
Arif

[alouis72]

Hi maxnest,
sorry for the late answer.
It seems that agora-generic is "stuck" when trying to process agora on the constrained families "1.0-1.0".

It's probably due to the fact that there are very few ancestral genes in this category for one ancestor after adding the new species.
Will it be possible for you to give me the number of conserved pairs for each ancestors in "pairwise/pairs-size-1.0-1.0/log" (without the name of the ancestors if it's private of course)?
by doing this command:

grep "conserved pairs for" pairwise/pairs-size-1.0-1.0/log

On the data before the insertion of the new genome and after the insertion?

you can contact me directly at [email protected] if you want.
regards,
Alex

[diekei]

Hi @alouis72 and colleagues,

Thanks for the great tool! I also got the same problem with my dataset, but instead that the problem appeared after adding one more genome, I couldn't run the generic pipeline in the first place (but works for the basic). I wonder what is the problem and whether you could help solving it. Thank you!

[alouis72]

Hi diekei,
I think that the problem is that when adding new genomes, the parameter of constrained genes 1.0-1.0 is to stringent, leading to graphs too complex to be parsed when the algorithm tries to fillin-in the first skeleton of ancestral blocks…

Maybe the thing to do, is to try generic pipeline only on specific ancestors (with option -target= ) and loop.
and stay on basic for the oldests ancestors.
Or
edit the agora-generic.py scripts, and remove the parameter 1.0-1.0 on line 29:

for sizeParams in [(1.0,1.0), (0.9,1.1), (0.77,1.33)]:
Becoming
for sizeParams in [(0.9,1.1), (0.77,1.33)]:

Hope this will work,
regards,
Alex