4
$\begingroup$

I made a phylogenetic tree for about 300 species using a region of the COX1 gene. The tree looks visually meaningful: species of the same genus are all clustered together.

Then I performed a 350-replicate bootstrapping and I analyzed the support values. Surprisingly I found very low values for nodes that are more close to the root and high support values for nodes that are more close to the leaves.

Almost all nodes that split into 2 leaves have Bootstrap >95% but the ones close to the root have Bootstrap<20%.

Does that mean I chose a highly variable DNA region?
How would you re-build the tree to make it more consistent on more ancient nodes? I was thinking of using a concatenation of sequences with both high and low variability regions, but i am not sure it would work.

$\endgroup$
1
  • $\begingroup$ What are the families and the Order under investigation? Is this a nucleotide or protein tree? If you are considering the 'control region of mtDNA it suggests this is a nucleotide tree. $\endgroup$
    – M__
    Commented Apr 2 at 10:05

1 Answer 1

2
$\begingroup$

This is the expected result. We can be much more certain about the terminal nodes of the tree than we can be about the tree's deep structure, trivially (but not exclusively) because of uncertainty propagation from the tips down into the internal nodes (a deep node is only as good as all of the nodes that lead to it). You can see a similar phenomenon in the various phylogenies in this classic paper.

It is also possible that your gene is not variable enough, not just that it is too variable. If there is a subset of highly informative positions in your gene, bootstrapping samples that miss those positions would be expected to yield noisy results.

In terms of how to resolve these, the traditional answer is "more data". This data can come in either the form of more species/tips or more sequences per species (ideally both). It is probably unsurprising that you aren't able to get very confident with a single gene. Even a rather older approach uses 31 universally conserved genes. As is traditional in statistics, we get more certain about our estimates as we use more data to arrive at them. In these methods, the genes (proteins) are indeed concatenated.

Note also that single gene trees can be positively misleading about the true phylogenetic relationship, due to incomplete lineage sorting or other such phenomena. In this case, you could have excellent bootstrap support for a wrong phylogeny.

In another paper, you will find the following statement:

In molecular phylogenetic analysis, it is often the case that some relationships are robust, and relatively "easy" to recover while others are difficult to resolve, leading to phylogenetic hypotheses that consist of a patchwork of well and poorly supported nodes. When difficult nodes are encountered, the next logical step is to add taxa and/or data under the reasonable assumption that additional taxa or characters might enable resolution and/or provide support for poorly supported nodes...

The amount of data required for resolution of difficult phylogenetic problems associated with short internodes, especially those deep in a tree can represent a particularly difficult challenge [5–7] that often requires massive amounts of sequence data to resolve.

Deep nodes are just harder than shallow nodes! The possible exception here is the very oldest splits, where enormous amounts of evolutionary time might separate outgroups quite clearly from a relatively messy subtree (see the Fishbein paper for more examples of this).

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

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