<p>You need bootstrap support using a model-based tree building algorhithm, via maximum likelihood (a few people use Bayes). The file format is relaxed phylip format (please submit a separate question if you have difficulties here - its a bit tricky).</p>

<p>I use standard RAxML <a href="https://github.com/stamatak/standard-RAxML" rel="nofollow noreferrer">here</a>, specifically raxmlHPC (easily downloadable and compiles on Linux and OSX). The codes are quite complicated and I've given them below. </p>

<p>A robust maximum likelihood tree is,</p>

<pre><code>raxmlHPC -m GTRGAMMA -n output.tre -o outgroup-p 10000 -s alignment_file.phy
</code></pre>

<p>To make a bootstrap tree,</p>

<pre><code>raxmlHPC -b 459 -m GTRGAMMA -n output.tre -o outgroup -p 10000 -s /path_to_alignment.phy -N 500
</code></pre>

<p>This tree will bootstrap for 500 replications, however to start I would use 100 replications.</p>

<p>Make a consensus tree of the bootstraps,</p>

<pre><code>raxmlHPC -f b -m GTRGAMMA -n output.tre -o outgroup -p 10000 -t /path_to_master_tree(best_tree) -z path_to_bootfile
</code></pre>

<p>You require bootstrap support >80% and please repeat this with and without 5-2/5-3 (it still looks long)</p>

<p>The cluster you have access to is fine for the calculation, it will take aound 24 hours to complete one bootstrap calculation for one data set and obviously you need to parallelise you calculation across the 22 contigs. </p>

<p>Viewing the tree, FigTree (for Mac OSX) is easy.</p>

<p><strong>Root</strong></p>

<p>Rooting can be complicated because I don't really know your bacteria.</p>

<hr>

<p>The recombination issue is more complicated, but I would construct 22 trees from your contigs and assess them for congruence. Panmixia is a concerned, which means too much recombination</p>

CC BY-SA 4.0
2020-03-14T20:18:52.947Z
2020-03-15T10:52:28.207Z
2020-03-15T10:52:28.207Z
M__
<p>I need to construct the phylogenetic tree of the whole genus of streptomyces. And <strong>my goal</strong> is to <strong>find out how my target strep species distributed in the whole strep genus and obtain the distribution pattern</strong>. Strep is a big genus, has 658 species(65% of them have the genome info), and my target species is about 25. <strong>But the traditional 16S method seems not helping when it comes to strep is a genetic recombination genus</strong>. And still lack of MLST scheme to distinguish my target species. I want to get the phylogeny for the whole genus, but I lack related genetics knowledge.</p>

<p>There is something I did before try to get the tree.</p>

<ol>
<li>I ran the orthomcl program for 400 strep genomes which I got from NCBI, and the result is: I got 22 orthologues</li>
<li>I ran the fastani program for 200 strep genomes(type material) and got the ani result.</li>
</ol>

<p>What can I do to achieve my goal? Thanks!</p>

CC BY-SA 4.0
2020-03-12T14:12:50.630Z
Hi @Tangli83, I will construct an answer in tomorrow evening (Saturday morning in China) based on our conversation
M__
CC BY-SA 4.0
2020-03-12T16:03:08.267Z
Oh, thanks,~ I thought my post didn't explain clearly.
TANGLi83
CC BY-SA 4.0
2020-03-13T04:37:32.887Z
that's fine I know your questions and research aims
M__
<fasta><phylogenetics><phylogeny><bacteria>
How to construct the phylogenetic tree of the "High Recombination" genus streptomyces?
# How to construct the phylogenetic tree of the "High Recombination" genus streptomyces?

(Asked by: username_0 on 2020-03-12)

I need to construct the phylogenetic tree of the whole genus of streptomyces. And **my goal** is to **find out how my target strep species distributed in the whole strep genus and obtain the distribution pattern**. Strep is a big genus, has 658 species(65% of them have the genome info), and my target species is about 25. **But the traditional 16S method seems not helping when it comes to strep is a genetic recombination genus**. And still lack of MLST scheme to distinguish my target species. I want to get the phylogeny for the whole genus, but I lack related genetics knowledge.

There is something I did before try to get the tree.

1. I ran the orthomcl program for 400 strep genomes which I got from NCBI, and the result is: I got 22 orthologues
2. I ran the fastani program for 200 strep genomes(type material) and got the ani result.

What can I do to achieve my goal? Thanks!

**Tags:** fasta, phylogenetics, phylogeny, bacteria
> **Comments:**
> **username_1**: Hi @Tangli83, I will construct an answer in tomorrow evening (Saturday morning in China) based on our conversation  
> **username_0**: Oh, thanks,~ I thought my post didn't explain clearly.  
> **username_1**: that's fine I know your questions and research aims  

---

## 1 answers

### Answer by username_1 on 2020-03-14. Score: 1

You need bootstrap support using a model-based tree building algorhithm, via maximum likelihood (a few people use Bayes). The file format is relaxed phylip format (please submit a separate question if you have difficulties here - its a bit tricky).
I use standard RAxML here, specifically raxmlHPC (easily downloadable and compiles on Linux and OSX). The codes are quite complicated and I've given them below. 
A robust maximum likelihood tree is,
raxmlHPC -m GTRGAMMA -n output.tre -o outgroup-p 10000 -s alignment_file.phy

To make a bootstrap tree,
raxmlHPC -b 459 -m GTRGAMMA -n output.tre -o outgroup -p 10000 -s /path_to_alignment.phy -N 500

This tree will bootstrap for 500 replications, however to start I would use 100 replications.
Make a consensus tree of the bootstraps,
raxmlHPC -f b -m GTRGAMMA -n output.tre -o outgroup -p 10000 -t /path_to_master_tree(best_tree) -z path_to_bootfile

You require bootstrap support >80% and please repeat this with and without 5-2/5-3 (it still looks long)
The cluster you have access to is fine for the calculation, it will take aound 24 hours to complete one bootstrap calculation for one data set and obviously you need to parallelise you calculation across the 22 contigs. 
Viewing the tree, FigTree (for Mac OSX) is easy.
Root
Rooting can be complicated because I don't really know your bacteria.

The recombination issue is more complicated, but I would construct 22 trees from your contigs and assess them for congruence. Panmixia is a concerned, which means too much recombination


---

