Question: [OMA] When input sequence data are DNA, OMA conversion in AA ?
0
gravatar for pourbiostars
13 days ago by
pourbiostars0 wrote:

Hello,

I'm actually using OMA with this parameter:

  • " # Type of input sequence data, has to be either 'DNA' or 'AA' InputDataType := 'DNA'; "

So OMA knows that I'm using DNA sequences. But I find this parameter inside the parameter file:

  • " # alignments which have a score lower than MinScore will not be considered. # The scores are in Gonnet PAM matrices units. MinScore := 181; "

Gonnet PAM matrices being amino acid based, I would like to know if OMA convert DNA into AA before the allVSall step ?

Thank you.

Best regards.

dna aa orthologous oma orthologs • 121 views
ADD COMMENTlink modified 11 days ago by Christophe Dessimoz310 • written 13 days ago by pourbiostars0
0
gravatar for Christophe Dessimoz
11 days ago by
University College London
Christophe Dessimoz310 wrote:

No, OMA does not translate DNA into AA.

While it's true that Gonnet PAM matrices were defined on amino-acid, it's also possible to derive empirical scoring matrices from DNA in the same way. In OMA standalone, if the input data type is DNA, we use such matrices to perform alignment and evolutionary distance inference.

ADD COMMENTlink written 11 days ago by Christophe Dessimoz310

Thank you for your answer,

When comparing species (from the same plant taxonomy family) does using CDS sequence instead of AA matter when the clustering is done ?

ADD REPLYlink written 10 days ago by pourbiostars0

We recommend using protein sequences, because we have done much more benchmarking on this type of data. There are practical advantages too: public genomes exported from OMA, including precomputed all-against-all, are only available at for protein sequences. By contrast, the DNA-level comparison support is still experimental in OMA. We have started using it for non-coding regions.

So for coding regions, I strongly recommend that you use the translated sequences.

ADD REPLYlink written 7 days ago by Christophe Dessimoz310
Please log in to add an answer.

Help
Access

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