Question: Bacterial whole genome tree
gravatar for jeccy.J
5.6 years ago by
jeccy.J50 wrote:

Hello all,

I was trying to draw whole genome phylogeny about 200 similar strains of bacterial genome. I align all the genome using Mauve aligner and got tree file and i can see it with any tree viewer. But when i gone through few recent paper i saw after alignment they usually extract the conserved block and draw the tree separately.  Can anyone suggest me which will be optimal solution.




whole genome phylogeny • 2.6k views
ADD COMMENTlink modified 5.0 years ago by Biostar ♦♦ 20 • written 5.6 years ago by jeccy.J50
gravatar for Asaf
5.6 years ago by
Asaf8.3k wrote:

You can draw trees using several biological data. Usually people use highly conserved proteins like ribosomal proteins, align them, select blocks of alignments and draw a tree using this data. Another option is rRNA or other markers like insertion sequences and more. You should take into account that evolution of proteins might be different from evolution of rRNA or non-coding DNA.

That being said, every tree you build is valid, it just represents the data you used to build it. Since you are dealing with closely related bacteria it might not be such a bad idea to use DNA rather than highly conserved proteins.

Another issue is the algorithm used by Mauve, I'm not familiar with it so I can't say anything about it. You should have a way to validate your tree so you could trust it.

Hope I helped.

ADD COMMENTlink written 5.6 years ago by Asaf8.3k
gravatar for 5heikki
5.6 years ago by
5heikki8.9k wrote:

I believe Gblocks and trimal are the most widely used alignment trimmers. What software did they use in the few recent papers you saw? Since you're aligning complete genomes, the alignment is probably extremely long and thus not suitable for Bayesian tree-building approaches. RAxML would probably suit your needs well but if you want something much faster, then start with FastTree. Also, if you want a rooted tree, remember that you'll need an outgroup..

ADD COMMENTlink modified 5.6 years ago • written 5.6 years ago by 5heikki8.9k

Thanks for reply. I saw the paper used Gblock to find conserved block and tree build by Figtree. If i got properly :acceding to your opinion my approach should be 1. align all genome using mauve 2. and draw tree using FastTree or RAxML. No need to find the conserved block after alignment.


ADD REPLYlink written 5.6 years ago by jeccy.J50

You read the article wrong because FigTree is for illustration of trees that were made with whatever algorithm. You also read my post wrong. RAxML and FastTree can make tree files (e.g. nexus format), but you can't illustrate said files with these programs.

ADD REPLYlink modified 5.6 years ago • written 5.6 years ago by 5heikki8.9k

sorry little modified tree was build by Mauve alignment and tree view by FigTree. 

for fast tree i saw in manual they mention i have to use like 
$FastTree -gtr -nt < alignment.file > tree_file 

Could you please let me know in alignment file can i use a .aln file generated by mauve or else i have to build with other tool such as Mugsy ??

ADD REPLYlink written 5.6 years ago by jeccy.J50
Please log in to add an answer.


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