The phylogenetic analysis comprises an astonishing 73,060 terminal eukaryotic taxa, 9535 molecular characters and, for good measure, they threw in 604 morphological characters. It is therefore the largest phylogenetic analysis published to date and almost six times larger than the former world record. Such feat presented many technical challenges. The logistics required the automatizing of every step in the analysis, via computer scripts, to retrieve and sort thousands of GenBank entries, to align the sequences to construct the data matrix, to perform the actual searches for the optimal solutions, and to interpretation of the mammoth-size phylogenetic trees. The crux of the analysis, the search for the optimal phylogenetic trees, was done with the powerful parsimony phylogenetic program TNT running in parallel in three multi-processor computers for 2.5 months.

Fig. 1. Pruned strict consensus tree for the combined data set (seven trees, 1879 taxa excluded). The bar shows the span of 5000 species.

The resulting phylogeny recovers most traditional taxonomic groups. This is interesting for various reasons. First, as noted about, our understanding of the tree of life is the results of many taxonomically localized efforts that have been informally pasted together¹. This is the first time a phylogeny has been reconstructed from scratch, letting the data speak unconstrained for itself without assuming that certain evolutionary relationships most be true a priori. Second, it shows that there is enough historical information contained in the data so that the optimal solution is not a complete mess or largely unresolved answer– consider that there are 9 X 10^345,593 possible tree combinations for the number of terminals included. Third, that we do have the current capacity, both in terms of software and hardware, to carry out such a large analysis. And last, but related to the previous two points, that parsimony methods for phylogenetic reconstruction are up for the task. The latter point is worth noting because early simulations, based on just a few taxa (a grand total of four actually) scared systematists into thinking that parsimony methods may result in erroneous reconstructions. Later studies using real data and a much larger collection of species has shown that this is not the case, and this 73,060 taxa analysis serves as the largest of these test cases.

The authors are no strangers when it comes to computer implementation of phylogenetic methods. James S. Farris is a pioneer in the field who developed the algorithmic foundations and produced the some of the first phylogenetic programs in the late 1960s, when the character information for each taxon to be analyzed was contained in a punch card and random addition sequence for the phylogenetic tree construction meant that the set of cards was shuffled by hand before feeding them into the terminal connected to the mainframe. Likewise, Pablo A. Goloboff has been responsible for many of the rapid search techniques developed during the 1990s up to the present, that seek to cover the searchable tree-space in a fast and efficient way.

It seems that, for phylogenetics, the only limit that remains is the availability of data.

References and notes

Goloboff, P., Catalano, S., Marcos Mirande, J., Szumik, C., Salvador Arias, J., Källersjö, M., & Farris, J. (2009). Phylogenetic analysis of 73 060 taxa corroborates major eukaryotic groups Cladistics DOI: 10.1111/j.1096-0031.2009.00255.x