Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Towards a Supertree of Arthropoda: A Species-Level Supertree of the Spiny, Slipper and Coral Lobsters (Decapoda: Achelata)

Abstract

While supertrees have been built for many vertebrate groups (notably birds, mammals and dinosaurs), invertebrates have attracted relatively little attention. The paucity of supertrees of arthropods is particularly surprising given their economic and ecological importance, as well as their overwhelming contribution to biodiversity. The absence of comprehensive archives of machine-readable source trees, coupled with the need for software implementing repeatable protocols for managing them, has undoubtedly impeded progress. Here we present a supertree of Achelata (spiny, slipper and coral lobsters) as a proof of concept, constructed using new supertree specific software (the Supertree Toolkit; STK) and following a published protocol. We also introduce a new resource for archiving and managing published source trees. Our supertree of Achelata is synthesised from morphological and molecular source trees, and represents the most complete species-level tree of the group to date. Our findings are consistent with recent taxonomic treatments, confirming the validity of just two families: Palinuridae and Scyllaridae; Synaxidae were resolved within Palinuridae. Monophyletic Silentes and Stridentes lineages are recovered within Palinuridae, and all sub-families within Scyllaridae are found to be monophyletic with the exception of Ibacinae. We demonstrate the feasibility of building larger supertrees of arthropods, with the ultimate objective of building a complete species-level phylogeny for the entire phylum using a divide and conquer strategy.

Introduction

The Achelata are part of the highly diverse Decapoda, and are typically classified into five families; the extant Palinuridae (spiny lobsters), Scyllaridae (slipper lobsters) and Synaxidae (furry or coral lobsters), plus the extinct Cancrinidae and Tricarinidae. Despite their common names Achelata and Nephropidae (true lobsters) are not sister groups. Monophyly of Achelata is supported by all formal analyses, but the relationships and monophyly of the three constituent families remain equivocal. Morphological data support monophyly of all three extant families, while fossil data imply that Palinuridae are paraphyletic with respect to Synaxidae [1]. The most complete molecular phylogeny [1] concatenated data from five genes (18S, 28S, H3, 16S and COI) for 35 taxa, and found support for just two major clades; Scyllaridae and a group comprising Palinuridae and Synaxidae, with the latter placed within Palinuridae. Despite the generation of a number of smaller data sets and trees focussing on particular genera, the issue of deeper achelatan phylogeny has not been revisited; neither have more inclusive trees been inferred. Here, we synthesise the corpus of published trees using a supertree approach, specifically to investigate its efficacy and to test the suitability of our methods and new software [2] for application to further analyses.

Supertree methods offer a practicable means by which to synthesise large numbers of smaller trees with partially overlapping leaf sets. These “source trees” can themselves have been inferred from any type of data (e.g., morphology or molecules), each using the particular set of analytical approaches deemed most appropriate by their authors (e.g., distance methods, variously complex parsimony, likelihood or Bayesian frameworks). This inclusivity contrasts with supermatrix approaches containing diverse data types. In such cases, it is necessary either to apply a single analytical model across all characters [35], or to use the trees from the optimal analysis of one data type (e.g., the maximum likelihood analysis of molecular data) to construct a constraint tree or scaffold for the analysis of another (e.g., parsimony analysis of morphology) [6].

Large, inclusive and complete (i.e., all known species) phylogenies are vital for a variety of applications in evolutionary biology, ecology and conservation. Inter-nested patterns of shared evolutionary history mean that the traits and attributes of species are expected to share a complex, hierarchically correlated structure. Ecological or behavioural studies of trait correlations and interactions must therefore include phylogenetic information so that the conflated effects of shared ancestry can be factored out [712]. Similarly, the value of conserving a particular species may be at least partially determined by its “evolutionary distinctiveness”; the extent to which it is phylogenetically distant or distinct from its nearest relatives [13].

Supertrees have been produced for a diverse array of major clades, although there has been considerable emphasis upon vertebrates; notably dinosaurs [14,15], crocodiles [16], mammals [17] and birds [18]. Surprisingly, there have been few published arthropod supertrees, despite the successful application of the approach to Adephaga (Coleoptera) at the level of genera [19], and at the family level to Odonata [20] and Hymenoptera [21]. Arthropoda probably contains between 2.5 and 10 million species, although estimates as high as 30 million have been obtained [2224]. Arthropods constitute more than 80% of described animal species, are abundant in virtually all habitats, and are of global economic [25] and ecological [26] importance. The ecosystem service value of arthropod pollinators alone is estimated at €153 billion per annum [25]. Given their pivotal ecological role in many habitats, a major loss of arthropod biodiversity would have catastrophic and far-reaching consequences [26]. At the same time, new species of arthropods (and especially insects) are being described at a rate exceeding that for all other animal, fungal and plant groups [27] with the conservation status of these species usually being unclear. An efficient means for synthesising the wealth of phylogenetic inferences for arthropods would therefore be valuable. We have chosen Achelata as an initial study clade since they are relatively small and well-documented, and contain representatives that have both consumptive (food source) and non-consumptive economic value (e.g., Palinuridae) [28].

Methods

Source Tree Collection

Potential source trees were identified from online resources. The Web of Knowledge Science Citation Index [29] was searched from 1980 to 2013 using the search terms: phylog*, taxonom*, systematic*, divers*, cryptic and clad* in conjunction with all scientific and common names for the Achelata from infra-order to sub-family level. All papers mentioning or implying the existence of a tree in their title or abstract were examined. In addition, the references cited by these papers were trawled for additional sources. All source trees and selected meta-data were digitised in their published form using TreeView [30] and the Supertree Toolkit (STK [2]). The latter is a fully integrated set of scripts designed to process trees and meta data, and to output matrices for MRP [31] supertree analysis or sets of trees for analysis using other supertree methods. The new version is either GUI or command line driven and offers much greater flexibility and functionality than its precursor that constituted a user defined processing pipeline [32]. Meta-data included bibliographic information, the types of characters used (e.g., molecular or morphological) and the methods used for tree inference. No corrections were made for synonyms or any other apparent errors or inconsistencies in the source trees prior to processing.

All the source tree data were deposited into our new resource at the Supertree Toolkit website [33]. This resource comprises a searchable, freely available database. All our source trees are archived here as they appear in their original published form, along with meta-data that allow further analyses re-purposing. This resource fills an important nîche, as few authors make their trees available in machine-readable form. All source trees curated for our arthropod supertrees will ultimately be archived here.

Data Processing

The tree presented here differs from all other previously published supertrees as it utilised the STK in order to standardise and partially automate the process of construction. It is vital to ensure that source trees are treated in a consistent and repeatable manner in assembling a supertree [34,35]. The STK was devised in order to increase the accuracy and uniformity of approach, as well as to speed data processing. We followed the protocol described by Davis and Page [18].

Once data collection and data entry were complete, we ensured that source trees met several criteria before inclusion in the analysis:

  1. Only trees presented by their authors explicitly as a reconstruction of evolutionary relationships were included. We therefore excluded taxonomies and informal phylogenies (i.e., we only included those derived from an explicit matrix of characters).
  2. Only phylogenies comprising clearly identified species, genera or higher taxa and clearly identifiable characters were included.
  3. Only trees derived from the analysis of a novel, independent dataset were included.

Non-independent studies were defined as those that utilised identical matrices (i.e., the same taxa and characters), or where one matrix was a subset of the other. In the former case, the “identical” source trees were weighted in inverse proportion to their number. In the latter case, the less inclusive tree was removed from the data set.

OTUs (operational taxonomic units) were standardised to reduce the inclusion of higher taxa, and to remove synonyms and vernacular names (which were standardised using the freely available online WoRMS database [36]). Where authors used higher taxa as proxies for particular exemplars, we substituted those with the names of those genera or species. Where no exemplars were specified, higher taxa were removed from source trees by substituting those constituent taxa present in other source trees as a polytomy in the focal tree. This avoided artificial inflation of the data set. Definitions for higher taxa were derived from WoRMS [36].

Taxonomic overlap was checked once the nomenclature had been standardised. Each source tree required at least two taxa in common with at least one other source tree [37]. Overlap within our dataset was sufficient; therefore no source trees were removed and we were able to proceed to matrix creation without any further edits. See S1 File [38] for the source trees as they were included in the analysis, S2 File [38] for a reference list of all source trees and S3 File [38] for the STK data file. Source trees in their original form were deposited in the Supertree Toolkit website database [33].

Supertree construction

Achelata are a relatively small group (~150 species), and our data set contained 531 species (118 ingroup) from 55 source trees (S3 File [38]). Our supertree was inferred using Matrix Representation with Parsimony (MRP; [31]); the most commonly used and most tractable approach with medium to large data sets [39]. Source trees were encoded as a series of group inclusion characters using standard Baum and Ragan coding [31], and automated within the STK software. All taxa subtended by a given node in a source tree were scored as “1”, taxa not subtended from that node were scored as “0”, and taxa not present in that source tree were scored as “?”. Trees were rooted with a hypothetical, all-zero outgroup [40]. The resulting MRP matrix (S4 File) was analysed using standard parsimony algorithms in TNT [41]. We used the “xmult = 10” option, and ran 1000 replicates for the analysis, each using a different random starting point for the heuristic search. This improved exploratory coverage of the tree space, potentially avoiding local minima in the solutions.

Results

The analysis found 3000 MPTs of length 2889 steps. We then computed a Maximum Agreement Subtree (MAST) in PAUP* 4.0b10 [42]; the resulting tree comprised 82 ingroup taxa and was fully resolved. Fig 1 shows the complete tree, see S5 File [38] for the MAST in Newick format.

thumbnail
Fig 1. Maximum Agreement Subtree (MAST).

The MAST was calculated for 3000 MPTs of length 2889 steps with V scores superimposed (V+ not shown). Green circles denote nodes with V scores of 1, orange denote nodes with V support between 0.01–0.9, black denotes V support of 0 and red denotes negative V support. The tree was generated using FigTree [43].

https://doi.org/10.1371/journal.pone.0140110.g001

We calculated support for each bifurcating node in the supertree using the V index [44]. Each source tree was examined to determine if it contained a given supertree node. V indices range between +1 (where all source trees contain the node) and -1 (where no source trees contain it). Values over zero are consistent with support in the majority of source trees. A more relaxed index, V+, also takes permitted relationships (i.e., those consistent with polytomies in the source trees) into account, and therefore tends to yield higher values than V. All deep nodes in the supertree received positive V and V+ scores. Only nine nodes in the supertree received a negative V score, and just one had a negative V+. One additional node received a V score of zero. The negative values were all found towards the tips of the supertree and within the genera Panulirus and Thenus. This suggests that the fundamental splits in the tree are very well-supported and that only the more recent, species-level relationships show significant discordance in the source trees.

The supertree contains 56% of all described species of Achelata, and synthesises source trees published from 1992 to 2013 (no suitable trees were found from the period 1980–1991). Two clades were removed in the MAST consensus: Justitia/Linuparus and Palinurus. This suggests that these taxa are highly mobile within the source trees and could benefit from further study. Other taxa removed by the MAST were those that are poorly represented and/or poorly constrained within the source trees. For example, Ibacus pubescens is only present in one source tree.

Bininda-Emonds and Bryant [45] noted that the MRP method can lead to the creation of spurious clades and relationships that are not present in any of the source trees (“novel clades”). Although simulations have suggested that such anomalies are unlikely to be a significant problem [46], empirical studies have found an incidence of novel clades affecting up to 3% of taxa in the study [18]. However, no novel clades were found in this analysis.

The majority of our source trees were derived from papers published post-2000 onwards (Fig 2), and the majority were derived from molecular characters, reflecting the increasing abundance of molecular studies in the literature (Fig 3). Although some data types have broad taxonomic coverage (e.g., adult morphology and ribosomal RNA genes), many are absent for large numbers of taxa (Fig 4A). Some of these (e.g., enzymes and mitochondrial DNA restriction sites) have been used only rarely in recent years, but even some otherwise commonly-sequenced genes (e.g., 12S RNA and COI) have not been sequenced for large numbers of taxa. Plotting taxon presence/absence in source trees and partitioning by family (Fig 4B) reveals that the source trees show a strong taxonomic bias; far better data coverage for Palinuridae than for Scyllaridae or Synaxidae. We also note that the size distribution of our source trees was strongly skewed towards small source trees with fewer than 20 taxa (Fig 5).

thumbnail
Fig 2. Number of phylogenies published by year.

The number of phylogenies published, and included in the supertree analysis, is heavily skewed towards recent years with relatively few trees from pre-2000.

https://doi.org/10.1371/journal.pone.0140110.g002

thumbnail
Fig 3. Source of character data by year.

There is a strong bias towards source trees derived from molecular data for all years and increasingly so from 2000 onwards.

https://doi.org/10.1371/journal.pone.0140110.g003

thumbnail
Fig 4. A) Characters sampled per taxon, partitioned by taxonomy. Palinuridae are better sampled than either Scyllaridae or Synaxidae for 50% of the characters used to build source trees. B) Number of taxa and presence/absence in source trees, partitioned by taxonomy. Palinuridae are better represented in the source trees than either Scyllaridae or Synaxidae.

https://doi.org/10.1371/journal.pone.0140110.g004

thumbnail
Fig 5. Source tree size shown by numbers of taxa present in source trees.

The data set is dominated by source trees with only a small (<20) number of taxa.

https://doi.org/10.1371/journal.pone.0140110.g005

Discussion

Achelatan Phylogeny

The Achelata supertree is fully-resolved and support is high throughout: 87% of nodes receive positive scores for both V and V+ indices (i.e., support or “permission” in the majority of the source trees for the V and V+ indices respectively). A well-supported, monophyletic Achelata is recovered, itself composed of large clades with high V indices; the Scyllaridae clade and the Palinuridae+Synaxidae clade. Synaxidae derive from within a paraphyletic Palinuridae, as reported by Palero et al. [1].

Within Palinuridae, Silentes and Stridentes are recovered as sister clades. The third proposed palinurid clade [1], comprising Linuparus and Justitia, does not have any taxa represented in the MAST phylogeny so cannot be assessed. Both Silentes and Stridentes receive high V and V+ scores (Silentes – 0.8; Stridentes – 0.55). All genera within Palinuridae are recovered as monophyletic with V and V+ scores of 1, with the exception of Panulirus which has a score of 0.8.

Within Scyllaridae all sub-families are recovered, with the exception of Ibacinae. The latter are split into two clades (Thenus and Evibacus/Parribacus) that are paraphyletic with respect to the Scyllarinae; a finding that is reflected in the source trees [1,47,48]. All sub-families and the two clades of Ibacinae have V indices of 1.00. All genera in the subfamilies Arctidinae, Ibacinae and Theninae are monophyletic, and also have V scores of 1.00. Within Scyllarinae all genera with the exception of Acantharctus, Scyllarus, Petrarctus and Galearctus are monophyletic, again all with postitive V and V+ scores. The polyphyly of Acantharctus is supported by the source data [48] clustering with Petrarctus and Scyllarus. The splitting of Scyllarus to cluster with Eduarctus and the inclusion of Acantharctus posteli are both relationships supported by the source trees [1,47,48]. Petrarctus is also split into two clades; again supported by the source trees [48,49]. The non-monophyly of Galearctus also reflects the source trees [49,50].

Towards a Supertree of Arthropoda?

Our species-level supertree of Achelata was constructed using a protocol [18] implemented in newly updated and freely available software; the Supertree Toolkit (STK) [2]. The STK allows users to input trees and meta-data via a full GUI, and implements numerous functions including the standardisation of taxonomy, the substitution of higher taxa, checks for adequate overlap between source trees and the down-weighting of non-independent source trees (many of these outputting graphical summaries). The protocol and software allowed us to build this tree relatively quickly, and also archived data in a form (XML) that can be readily re-analysed and re-purposed by other workers. We will utilise this approach to generate further supertrees of Arthropoda: a phylum that has received surprisingly little attention from supertree workers.

Our supertree represents the source data well. Despite containing only 56% of described species, it is nevertheless the most complete phylogeny of Achelata produced to date. We highlight taxa that appear to be in need of further study as well as areas of the tree that are well-supported.

Supertrees are now widely accepted as a valid means of obtaining large, complete phylogenies relatively quickly, and without the need to collate and analyse primary data. As such, they maximise the value of source trees already in the literature. Each of these are constructed by experts on their focal groups, and usually represent the investment of considerable analytical time and computational resources. Although the size of supermatrices that can be analysed within tractable time frames is increasing, there is still disagreement and incomplete overlap between the results of phylogenomic and other large studies. Supertree methods are therefore likely to remain important as a means to synthesise trees resulting from these largest of analyses.

Supertrees have many potential applications in comparative biology and macroevolutionary studies; published examples include diversification rates through deep time, origins of modern taxa, and origins of species richness [14,20,21]. However they are particularly useful in the fields of comparative trait analysis and conservation, where large and inclusive cladograms are needed in order to remove the effects of phylogenetic correlation and to quantify evolutionary distinctiveness respectively [51,52]. Yet, their utility could be further enhanced in a number of ways.

  1. Although supertree methods readily allow the inclusion of trees derived using any method from any type of data, there are still likely to be many described species that are not included in any primary phylogenies. The inclusivity of supertrees is therefore limited by the progress of the wider systematic community, although the largest trees are likely (almost by definition) to be supertrees at any given time. We will explore methods for utilising taxonomic information, thereby allowing supertrees to include all described species prior to their inclusion in published phylogenies. This will yield complete but unresolved supertrees that can be refined as new phylogenies are published.
  2. At present, supertrees comprising hundreds or thousands of taxa can only be built within tractable search times using MRP methods. These parsimony supertrees contain no valid branch length information, and require time-consuming post hoc calibration (e.g., using fossils) in order to set them against any absolute or even relative time frame.
  3. Although we found no rogue taxa in this study, the phenomenon of spurious clades in MRP supertrees is well-documented [45]. Bayesian and Maximum Likelihood supertree methods that may obviate such problems are being developed [53,54], but currently such approaches cannot handle data sets with more than a few tens of terminals.

Addressing these issues will enable the building of larger and more complete species trees without the problems associated with current supertree methods. Eventually supertree methods are likely to be rendered obsolete as computing power increases and more sequence data becomes available but they are still necessary for the foreseeable future; one way forward may be the combination of supertree approaches with supermatrix methods [55,56] as a means to “divide and conquer” large data sets.

Supporting Information

S1 File. Achelata source trees in their original published form.

https://doi.org/10.1371/journal.pone.0140110.s001

(TRE)

S2 File. Bibliography for Achelata source trees.

https://doi.org/10.1371/journal.pone.0140110.s002

(BIB)

S3 File. Supertree Toolkit file with Achelata trees and meta-data.

https://doi.org/10.1371/journal.pone.0140110.s003

(PHYML)

Acknowledgments

We are grateful to all those who supplied source trees and to two reviewers whose comments helped us to improve and refine this manuscript.

Author Contributions

Conceived and designed the experiments: KED MAW. Performed the experiments: KED TWH. Analyzed the data: KED TWH CD. Wrote the paper: KED MAW.

References

  1. 1. Palero F, Crandall KA, Abelló P, Macpherson E, Pascual M. Phylogenetic relationships between spiny, slipper and coral lobsters (Crustacea, Decapoda, Achelata). Molecular Phylogenetics & Evolution 2009;50:152–162.
  2. 2. Hill J, Davis KE. The Supertree Toolkit 2. A new and improved software package with a Graphical User Interface for supertree construction. Biodiversity Data Journal 2014;2:e1053. pmid:24891820
  3. 3. Bernhard D, Ribera I, Komarek A, Beutel RG. Phylogenetic analysis of Hydrophiloidea (Coleoptera: Polyphaga) based on molecular data and morphological characters of adults and immature stages. Insect Systematics & Evolution 2009;40:3–31.
  4. 4. Dabert J, Dabert M, Mironov S. Phylogeny of feather mite subfamily Avenzoariinae (Acari: Analgoidea: Avenzoariidae) inferred from combined analyses of molecular and morphological data. Molecular Phylogenetics & Evolution 2001;201:124–135.
  5. 5. Lee MSY, Hugall AF, Lawson R, Scanlon JD. Phylogeny of snakes (Serpentes): combining morphological and molecular data in likelihood, Bayesian and parsimony analyses. Systematics & Biodiversity 2007;5:371–389.
  6. 6. Geisler JH, McGowen MR, Yang G, Gatesy J. A supermatrix analysis of genomic, morphological, and paleontological data from crown Cetacea. BMC Evolutionary Biology 2011;11:112. pmid:21518443
  7. 7. Altermatt F, Pearse IS. Similarity and Specialization of the Larval versus Adult Diet of European Butterflies and Moths. American Naturalist 2011;178:372–82. pmid:21828993
  8. 8. Higgins L, Coddington J, Goodnight C, Kuntner M. Testing ecological and developmental hypotheses of mean and variation in adult size in nephilid orb weaving spiders. Evolutionary Ecology 2011;25:1289–306.
  9. 9. Higginson DM, Miller KB, Segraves KA, Pitnick S. Female reproductive tract form drives the evolution of complex sperm morphology. Proceedings of the National Academy of Sciences of the United States of America 2012;109:4538–43. pmid:22323584
  10. 10. Lind EM, Barbosa P. Life history traits predict relative abundance in an assemblage of forest caterpillars. Ecology 2010;91:3274–83. pmid:21141188
  11. 11. Poeyry J, Luoto M, Heikkinen RK, Kuussaari M, Saarinen K. Species traits explain recent range shifts of Finnish butterflies. Global Change Biology 2009;15:732–43.
  12. 12. Sekar S. A meta-analysis of the traits affecting dispersal ability in butterflies: can wingspan be used as a proxy? Journal of Animal Ecology 2012;81:174–84. pmid:21988561
  13. 13. Isaac NJB, Turvey ST, Collen B, Waterman C, Baillie JEM. Mammals on the EDGE: Conservation Priorities Based on Threat and Phylogeny, PLOS ONE 2007;2:e296. pmid:17375184
  14. 14. Lloyd GT, Davis KE, Pisani D, Tarver JE, Ruta M, Sakamoto M, et al. Dinosaurs and the Cretaceous Terrestrial Revolution. Proceedings of the Royal Society B—Biological Sciences 2008;275:2483–2490.
  15. 15. Pisani D, Yates AM, Langer MC, Benton MJ. A genus-level supertree of the Dinosauria. Proceedings of the Royal Society B—Biological Sciences 2002;269:915–921.
  16. 16. Gatesy J, Baker RH, Hayashi C. Inconsistencies in arguments for the supertree approach: supermatrices versus supertrees of Crocodylia. Systematic Biology 2004;53:342–55. pmid:15205058
  17. 17. Bininda-Emonds ORP, Cardillo M, Jones KE, MacPhee RDE, Beck RMD, Grenyer R, et al. The delayed rise of present-day mammals. Nature 2007;446:507–51. pmid:17392779
  18. 18. Davis KE, Page RDM. Reweaving the Tapestry: a Supertree of Birds. PLOS Currents Tree of Life 2014; Jun 9. Edition 1.
  19. 19. Rolf G. Beutel RG, Ribera I, Bininda-Emonds ORP. A genus-level supertree of Adephaga (Coleoptera). Organisms Diversity & Evolution 2008;7:255–269.
  20. 20. Davis RB, Nicholson DB, Saunders ELR, Mayhew PJ. Fossil gaps inferred from phylogenies alter the apparent nature of diversification in dragonflies and their relatives. BMC Evolutionary Biology 2011;11:252. pmid:21917167
  21. 21. Davis RB, Baldauf SL, Peter J Mayhew PJ. The origins of species richness in the Hymenoptera: insights from a family-level supertree BMC Evolutionary Biology 2010;10:109. pmid:20423463
  22. 22. Adis J. Thirty million arthropod species—too many or too few? Journal of Tropical Ecology 1990;6:115–118.
  23. 23. May RM. The dimensions of life on Earth. In: Raven PH, editor. Nature and Human Society: The Quest for a Sustainable World. Washington, D.C.: National Academy Press; 2000. p. 30–45.
  24. 24. Wilson EO. 2000. The Creation of Biodiversity. In: Raven PH, editor. Nature and Human Society: The Quest for a Sustainable World. Washington, D.C.: National Academy Press; 2000. p. 22–29.
  25. 25. Gallai N, Salles J-M, Settele J, Vaissière BE. Economic valuation of the vulnerability of world agriculture confronted with pollinator decline. Ecological Economics 2009;68:810–82.
  26. 26. Redak RA. Arthropods and Multispecies Habitat Conservation Plans: Are We Missing Something? Environmental Management 2000;26:S97–S107. pmid:10801993
  27. 27. Costello MJ, May RM, Stork NE. Can We Name Earth's Species Before They Go Extinct? Science 2013;339:413–416. pmid:23349283
  28. 28. Rudd MA. The non-consumptive economic value of spiny lobster, Panulirus argus, in the Turks and Caicos Islands. Environmental Conservation 2001;28:226–234.
  29. 29. Web of Knowledge Science Citation Index [Internet]–[cited 3rd March 2015]. Available: http://wok.mimas.ac.uk.
  30. 30. Page RDM. Tree View: An application to display phylogenetic trees on personal computers. Computer Applications in the Biosciences 1996;12:357–358. pmid:8902363
  31. 31. Baum BR, Ragan MA. The MRP method. In: Bininda-Emonds ORP, editor. Phylogenetic Supertrees: Combining information to reveal the tree of life. Dordrecht, The Netherlands: Kluwer Academic; 2004. p. 17–34.
  32. 32. Davis KE, Hill J. The supertree tool kit. BMC research notes 2010; 3:95. pmid:20377857
  33. 33. Hill J, Tovar J, Davis KE The Supertree Toolkit database [Internet]–[cited 3rd March 2015]. Available: http://supertreetoolkit.org/
  34. 34. Bininda-Emonds ORP, Jones KE, Price SA, Cardillo M, Grenyer R, Purvis A. Garbage in, garbage out: data issues in supertree construction. In Bininda-Emonds ORP, editor. Phylogenetic Supertrees: Combining information to reveal the tree of life. Dordrecht, The Netherlands: Kluwer Academic; 2004. p. 267–280.
  35. 35. Gatesy J, Baker CRH, Hayashi C. Inconsistencies in arguments for the supertree approach: supermatrices versus supertrees of Crocodylia. Systematic Biology 2004;53:342–355. pmid:15205058
  36. 36. WoRMS Editorial Board [Internet]. World Register of Marine Species–[cited 2015-03-03]. Available: http://www.marinespecies.org.
  37. 37. Sanderson MJ, Purvis A, Henze C. Phylogenetic supertrees: assembling the tree of life. Trends in Ecology & Evolution 1998;13:105–109.
  38. 38. Davis KE, Hesketh TW, Delmer C, Wills MA. Figshare [Internet]–[cited 29th July 2015]. Available: share.com/articles/Achelata_Supertree_Data/1321250.
  39. 39. Cotton JA, Slater CSC, Wilkinson M. Discriminating supported and unsupported relationships in supertrees using triplets. Systematic Biology 2006;55:345–350. pmid:16611604
  40. 40. Ragan MA. Phylogenetic inference based on matrix representation of trees. Molecular Phylogenetics & Evolution 1992;1:53–58.
  41. 41. Goloboff PA, Farris JS, Nixon KC. TNT, a free program for phylogenetic analysis. Cladistics 2008;24:774–786.
  42. 42. Swofford DL. PAUP*. Phylogenetic Analysis Using Parsimony (*and other methods). Sunderland, Massachusetts: Sinauer Associates; 2002.
  43. 43. Rambaut, A. FigTree [Internet]–[cited 3rd March 2015]. Available: http://tree.bio.ed.ac.uk/software/figtree/
  44. 44. Wilkinson M, Pisani D, Cotton J, Corfe I. Measuring support and finding unsupported relationships in supertrees. Systematic Biology 2005;54:823–831. pmid:16243766
  45. 45. Bininda-Emonds ORP, Bryant HN. Properties of matrix representation with parsimony analysis. Systematic Biology 1998;47:497–508. pmid:12066690
  46. 46. Bininda-Emonds ORP. Novel versus unsupported clades: assessing the qualitative support for clades in MRP supertrees. Systematic Biology 2003;52:839–848. pmid:14668120
  47. 47. Bracken HDB, Toon A, Felder DL, Martin JW, Finley M, Rasmussen J, et al. The decapod tree of life: Compiling the data and moving toward a consensus of decapod evolution. Arthropod Systematics & Phylogeny 1008;67:99–116.
  48. 48. Yang C-H, Bracken-Grissom H, Kim , Crandall KA, Chan T-Y. Phylogenetic relationships, character evolution, and taxonomic implications within the slipper lobsters (Crustacea: Decapoda: Scyllaridae). Molecular Phylogenetics & Evolution 2012;62:237–250.
  49. 49. Yang C-H, Chen I-S, Chan T-Y. A new slipper lobster of the genus Petrarctus (Crustacea: Decapoda: Scyllaridae) from the West Pacific. The Raffles Bulletin of Zoology Supplement 2008;19:71–81.
  50. 50. Yang C-H, Chen I-S, Chan T-Y. A new slipper lobster of the genus Galearctus Holthuis, 2002 (Crustacea, Decapoda, Scyllaridae) from New Caledonia. Zoosystema 2011;33:207–217.
  51. 51. Huang D. Threatened Reef Corals of the World. PLOS ONE 2012 7(3): e34459. pmid:22479633
  52. 52. Nyakatura K, Bininda-Emonds Olaf RP. Updating the evolutionary history of Carnivora (Mammalia): a new species-level supertree complete with divergence time estimates. BMC Biology 2012;10:12. pmid:22369503
  53. 53. Oliveira Martins L, Mallo D, Posada D. A Bayesian Supertree Model for Genome-Wide Species Tree Reconstruction. Systematic Biology 2015;
  54. 54. Akanni WA, Creevey CJ, Wilkinson M, Pisani D. L.U.St: a tool for approximated maximum likelihood supertree reconstruction. BMC Bioinformatics 2014;15:183. pmid:24925766
  55. 55. Bininda-Emonds ORP. The future of supertrees: bridging the gap with supermatrices. Palaeodiversity 2010;3 Supplement:99–106.
  56. 56. Von Haeseler A. Do we still need supertrees? BMC Biology 2012;10:13. pmid:22369571