1
$\begingroup$

I ran a calculation in RAxML to determine the majority consensus phylogeny of a maximum likelihood bootstrap (How to show bootstrap values on a phylogenetic tree constructed with RAxML), and I got three output files:

  1. RAxML_bipartitions.output_bootstrap.tre
  2. RAxML_bipartitionsBranchLabels.output_bootstrap.tre
  3. RAxML_info.output_bootstrap.tre

What is the difference between the RAxML_bipartitions.output_bootstrap.tre file and RAxML_bipartitionsBranchLabels.output_bootstrap.tre?

Thank you for your help!

$\endgroup$
  • $\begingroup$ Done it, I managed to meet my own deadline :-) $\endgroup$ – Michael Dec 8 '19 at 20:04
  • 1
    $\begingroup$ Great! Thank you very much! ; ) $\endgroup$ – Leah Dec 10 '19 at 1:11
1
$\begingroup$

In summary,

RAxML_bipartitions.output_bootstrap.tre

Is the only file of interest. The reason this is true in this context is really complicated and you have to understand the statistics of likelihood and how they are interpreted within phylogeny to understand why. This file is simply the final output of a non-parametric bootstrap analysis performed by maximum likelihood.

What on earth is a non-parametric boostrap?

A non-parametric bootstrap is resampling each alignment position with replacement. Thus if we have alignment positions 1,2,3,4,5 A bootstrap resample for 2 replicates might be,

Replicate 1

1,1,3,5,2

Replicate 2

4,2,5,2,1

The ML algorithm will make trees of replicates 1 and 2 and find the consensus between them. If you think about it in any other context a bootstrap replicate is pretty meaningless because it no longer reflects the true biological sequence. Thus information on how the consensus was derived, are not really of interest to us providing we are confident this has been done correctly, viz. RAxML_bipartitionsBranchLabels.output_bootstrap.tre and RAxML_bipartitionsBranchLabels.output_bootstrap.tre

So why is this output of limited use?

There are situations to some investigators this information is useful, but assess the robustness of a tree topology its not needed. The only thing we want is a phylogram (bestTree) with the bootstrap values superimposed on them. We really don't need complicated stuff such as the tree to be represented for example as a polytomy (non-bifurcating tree) because we can just read the bootstraps to make that deduction (values >> 75%). In addition, there is not perfect consensus what boostrap value constitutes robustness, but generally most agree >80% is robust.


What output files have useful information in them?

The information that is important are the files associated with "bestTree", that was the single maximum likelihood tree performed on the intact native sequence. The "info" file for this contains 3 really important parameters:

  • -lnL ... very important!!
  • Gamma distribution parameter "alpha",
  • PINVAR, proportion of invariant sites,

-lnL is the highest log-likelihood (probability) of the phylogeny. It is usually a very small number for which where is an enormous amount of theory over it.

Alpha parameter of the gamma distribution this is the shape parameter of the mutation rate, if it is very low (<1) the distribution of mutations across the alignment is very tight clustered and approximates to a negative binomial distribution. This means some sites don't mutate at all and a small number of sites mutate alot. If it is very large >200 (which is never observed) it approximates to the Poisson distribution, meaning the mutation distribution is randomised across the alignment.

PINVAR this is a straight percentage/frequency and simply means the sites that don't mutate.

How are they calculated?

PINVAR and alpha are not emperically calculated, i.e. if you look at an alignment and say 'no mutations at that position', PINVAR would of course agree but may consider other invariant depending on the phylogeny. These parameters are calculated by maximum likelihood and you can begin to see why the calculation takes so long ... alpha and PINVAR affect the tree topology (which affects -lnL), but the topology affects alpha and PINVAR. Thus, is a multidimensional search of tree and parameter space.

So what stuff do I report in my Results?

Anyway reporting -lnL is good technique and shows the reader you've done maximum likelihood, citing PINVAR and alpha from gamma distribution helps ('Methods' parameters were calculated reiteratively under maximum likelihood). This is only useful for bestTree. The -lnL, PINVAR and gamma's alpha are also calculated for every single bootstrap replicate, but these values are of limited to use, because we have resampled the data, only the consensus tree counts... Obviously presenting the bootstrapped phylogram is extremely important.

Welcome to the technical world of phylogeny!

The amino acid matrix you used BTW .. LG is in vogue right now.

How do I do it?

When I do this stuff its via Biopython and ETE3, I capture the values within the pipeline and don't examine the output files of RAxML because I generate my own.

|improve this answer|||||
$\endgroup$
  • 1
    $\begingroup$ Thanks a lot for such a complete answer @Michael G; I greatly appreciate your help! $\endgroup$ – Leah Dec 10 '19 at 1:10
  • $\begingroup$ Hi Leah, no problem at all. If anything else arises please pop a question on site (a note on a question helps). The only trick now is making sure all the methods are correctly referenced (LG matrix is an MBE paper around 2010, but there has been some impressive applications). You might want to acknowledge Bioinfo SE. $\endgroup$ – Michael Dec 10 '19 at 13:36
  • 1
    $\begingroup$ Thanks @Michael G, will do! $\endgroup$ – Leah Dec 13 '19 at 1:41

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.