This is a test version of Biostars. For the public version, visit https://www.biostars.org.
multiple alignment of multiple genes to generate a phylogenetic tree

After we collected hundreds of one-copy orthologs among tens of specices, we want to generate a phylogenetic tree. But for the alignments of those sequences we are not sure what is the canonical steps?

(1) Should we concatenate protein or DNA sequences in each species then do the multiple alignment? Or should we do the multiple alignment for each ortholog and then concatenate the alignments of each orthologs? (2) Should DNA sequences be used or protein sequences? (3) What is the strategy to manually edit the alignments? Should the region showing a gap in any species be removed?

Thanks

phylogenetic tree
• 4,460 views • link
updated by Joe Joe • written by qwzhang0601 qwzhang0601

What kind of tree are you trying to build ? A species tree ?

• 0 views • link
written by Jean-Karim Heriche Jean-Karim Heriche

Yes. I want to generate a species tree

• 0 views • link
written by qwzhang0601 qwzhang0601

1 answer

You could concatenate the sequences but I'd be more inclined to calculate alignments and trees separately for each gene and then compute a consensus tree.

See ASTRAL-II for instance

You'd save yourself some work doing the alignments and concatenation them instead of aligning pre-concatenated sequences.

Whether you should use protein or NA depends on the question and genes involved. If they are potentially quite diverse, use the protein sequence or vice versa.

• 0 views • link
written by Joe Joe

Thank you for your suggestions. You mentioned "You could concatenate the sequences but I'd be more inclined to calculate alignments and trees separately for each gene and then compute a consensus tree." Would you please explain how to compute a consensus tree using the trees by separate genes?

Besides, I am using MEGA7, do you have suggestions on how to generate a phylogenetic tree of tens of species based on hundreds of single-copy orthologs among them?

Thanks

• 0 views • link
written by qwzhang0601 qwzhang0601

Computing the concensus tree is exactly what ASTRAL does so check it out. The input is simply the trees for each of the genes.

You should be able to use any program capable of producing standard format trees. I know Astral accepts newick trees. I tend to trust the trees out of RAxML personally but it depends on your sequences/alignments as to how good your trees will be

• 0 views • link
written by Joe Joe

Log in to answer this question.

More posts like this