This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Phylogenetic tree from denovo assemblies

Is there any neat way to do this? Ie to quickly get a core genome alignment from de novo assemblies and create a phylogenetic tree? I have hundreds of samples..

wgs phylogeny
• 3,856 views • link
updated by Sergey Naumenko Sergey Naumenko • written by samuel.lipworth • samuel.lipworth

Are you just looking for software that can create phylogenetic trees? If so, programs like SeaView work with alignment.

• 0 views • link
written by Hussain Ather Hussain Ather

Yeh but seaview cannot cope with hundreds of bacterial genomes, it crashes..

• 0 views • link
written by samuel.lipworth • samuel.lipworth

1 answer

Hello Samuel!

Not sure about neat and quickly. It depends on the biology of your species.

With very small genomes (viruses, 10k) you may try to

  1. make a multiple alignment of whole genomes with muscle (https://www.drive5.com/muscle/) or another aligner, and then
  2. build a phylogenetic tree with RAXML (https://sco.h-its.org/exelixis/software.html) or another program.
  3. visualize it with Dendroscope(http://dendroscope.org/)

If you are working with bacteria (5M) or eucariotic genomes (1M-3G), and you have different species, you may try to extract 16S/18S rRNA sequences and build trees upon them.

If you have individuals of the same species, or if you want more precise phylogeny, it is a longer story.

You have to annotate your genomes - mark up where genes are. It depends on species, for non-model organisms it may become a whole enterprise.

When you have you genes annotated, you have to

  1. cluster orthologous genes with Orthomcl(http://orthomcl.org/orthomcl/) (here paralogues are waiting to make your task even more challenging)
  2. make multiple alignments,
  3. concatenate them into one large alignment,
  4. infer the tree,
  5. visualize it.

You may use not all genes, but a smaller subset of 10-100 genes, the most conservative, housekeeping genes. The set is different for Bacteria/Eukarya.

Good luck!

SN

• 0 views • link
written by Sergey Naumenko Sergey Naumenko

This is a very useful comment! One addendum is that Orthofinder (a proteome clustering software similar to OrthoMCL but with more features and it accounts for some weirdnesses in the alignments) also gives you a file of one-to-one clusters as an easier start for the concatenated (MUSCLE?) aligment

Edit: Orthofinder also gives you a species tree but I have not really understood how it's being made except based on the distance matrix that is somehow generated from the clusters. That phylogeny does not contain bootstrap values or other confidence measurements, but it's a good start

• 0 views • link
written by Philipp Bayer Philipp Bayer

Log in to answer this question.

More posts like this