- Research article
- Open Access
The colonization of land by animals: molecular phylogeny and divergence times among arthropods
BMC Biology volume 2, Article number: 1 (2004)
The earliest fossil evidence of terrestrial animal activity is from the Ordovician, ~450 million years ago (Ma). However, there are earlier animal fossils, and most molecular clocks suggest a deep origin of animal phyla in the Precambrian, leaving open the possibility that animals colonized land much earlier than the Ordovician. To further investigate the time of colonization of land by animals, we sequenced two nuclear genes, glyceraldehyde-3-phosphate dehydrogenase and enolase, in representative arthropods and conducted phylogenetic and molecular clock analyses of those and other available DNA and protein sequence data. To assess the robustness of animal molecular clocks, we estimated the deuterostome-arthropod divergence using the arthropod fossil record for calibration and tunicate instead of vertebrate sequences to represent Deuterostomia. Nine nuclear and 15 mitochondrial genes were used in phylogenetic analyses and 61 genes were used in molecular clock analyses.
Significant support was found for the unconventional pairing of myriapods (millipedes and centipedes) with chelicerates (spiders, scorpions, horseshoe crabs, etc.) using nuclear and mitochondrial genes. Our estimated time for the divergence of millipedes (Diplopoda) and centipedes (Chilopoda) was 442 ± 50 Ma, and the divergence of insects and crustaceans was estimated as 666 ± 58 Ma. Our results also agree with previous studies suggesting a deep divergence (~1100 – 900 Ma) for arthropods and deuterostomes, considerably predating the Cambrian Explosion seen in the animal fossil record.
The consistent support for a close relationship between myriapods and chelicerates, using mitochondrial and nuclear genes and different methods of analysis, suggests that this unexpected result is not an artefact of analysis. We propose the name Myriochelata for this group of animals, which includes many that immobilize prey with venom. Our molecular clock analyses using arthropod fossil calibrations support earlier studies using vertebrate calibrations in finding that deuterostomes and arthropods diverged hundreds of millions of years before the Cambrian explosion. However, our molecular time estimate for the divergence of millipedes and centipedes is close to the divergence time inferred from fossils. This suggests that arthropods may have adapted to the terrestrial environment relatively late in their evolutionary history.
The terrestrial environment has been greatly altered by the actions of organisms over Earth's history. Prokaryotes were probably the first organisms to colonize land, and this occurred as early as 2.6 billion years ago [1–3]. The presence of organisms on exposed land will accelerate weathering through physical and chemical processes and may in turn affect the global atmosphere and climate . Therefore, it is of interest to know when different groups of organisms colonized land to better understand their effect on the biosphere.
The earliest undisputed fossils of terrestrial plants, animals, and fungi are all from the early Palaeozoic (Ordovician and Silurian; 495 – 417 million years ago, Ma) [5–9]. However, the earliest animal fossils are known from ~600 Ma [10–12] and the earliest representatives of the "plant lineage", such as green algae  and red algae , from even earlier (1000 – 750 Ma and 1200 Ma, respectively). This raises the possibility that land was colonized by multicellular eukaryotes prior to the Ordovician. A previous molecular clock analysis addressed the question of land colonization by plants and fungi, resulting in early time estimates of about 700 and 1000 Ma, respectively . However, molecular clock analyses have not addressed the colonization of land by animals.
Myriapods (centipedes, millipedes) and chelicerates (e.g., arachnids, horseshoe crabs) have figured prominently in the earliest evidence of terrestrial animals. For example, the first taxonomically identifiable body fossils of terrestrial animals are arachnids and chilopods from the late Silurian (~419 Ma) of England , and the oldest unambiguous evidence of sub-aerial animal activity is of arthropod (diplopod-like) trackways from ~450 Ma sediments . Older arthropod traces (possibly from the latest Cambrian) have been reported from terrestrial sediments, eolian dune deposits accumulated in a sandy beach environment of southern Canada , although the terrestrial nature of the trace makers is uncertain. Older, marine myriapod-like fossils are known , but no marine remains attributable to Chilopoda or Diplopoda have been found.
Here, we use molecular clock and phylogenetic methods to place temporal constraints on the early history of arthropods and the colonization of land by animals. Our focus is on the divergence of millipedes and centipedes because they represent the most ancient living lineages of terrestrial animals (based on the fossil record) whose common ancestor presumably was terrestrial. In turn, this provides a minimum estimate for the time when land was colonized by animals. With the use of a phylogenetic framework derived from analyses of all available nuclear genes, we also estimate times of other major divergences in the history of living arthropods, and the divergence of arthropods and deuterostomes. These additional data help to constrain a maximal time for the colonization of land by arthropods.
All phylogenetic analyses resulted in significant (>95%) clustering of myriapods and chelicerates (figure 1). Pancrustacea (insects and crustaceans) was found in most analyses, the only exceptions being minimum evolution with paralinear distances, and both weighted and unweighted parsimony analyses of nuclear + mitochondrial genes, and the unweighted parsimony analysis of only nuclear genes (Table 1). However, because these exceptional analyses used models that are less complex than recommended , it is possible that they were influenced by substitutional biases. Support for Pancrustacea was significant in most analyses of nuclear genes and the Bayesian analysis of nuclear + mitochondrial genes. The two representative crustaceans (branchiopods and malacostracans) formed a group in the nuclear gene analysis and the Bayesian all-gene analysis. However, insects joined malacostracans in the remaining all-gene analyses. This uncertainty of relationships within Pancrustacea has been encountered previously  and may be related to different rates of evolution among genes.
Divergence time estimates for the millipede-centipede split ranged from 357 – 537 Ma across different molecular clock methods, with an average of 442 ± 50 Ma (Table 2). Average time estimates for the other divergences among arthropods were 475 ± 53 Ma for Xiphosura-Arachnida, 642 ± 63 Ma for Myriapoda-Chelicerata, 614 ± 23 Ma for Branchiopoda-Malacostraca, 666 ± 58 Ma for Insecta-Crustacea, and 725 ± 46 Ma for Pancrustacea-(Chelicerata-Myriapoda). As expected, the use of different calibrations resulted in different time estimates, and although the penalized likelihood method (SGGPL) showed the greatest sensitivity to calibration differences, most methods yielded similar time estimates.
The use of exclusively fossil calibration points (Chilopoda-Diplopoda, Xiphosura-Arachnida) or molecular calibration points (Arthropoda-Deuterostomia) did not result in substantially different time estimates (table 2). Therefore, time estimates using the largest number of calibration points and proteins (averaged across methods) were chosen to summarize the time of divergence for each node (Table 2, Figure 2).
Using the assumption that the ancestral arthropod was aquatic (marine) and that terrestrialization is derived within arthropods, myriapods colonized land after the origin of the myriapod lineage (chelicerate-myriapod divergence; 642 ± 63 Ma) and before the millipede-centipede divergence (442 ± 50 Ma) or the earliest terrestrial fossils (~420 Ma). Time estimates for the deuterostome-arthropod divergence ranged from 913 – 1554 Ma across all methods and calibrations (average of 1130 ± 120 Ma), although SGGPL gave substantially higher estimates than the other methods (Table 3). The Bayesian (SGLDT and SGLMDT), multigene global (MGGLS), and supergene global (SGGLS) methods resulted in estimates of ~1100 – 900 Ma regardless of calibration point used or whether a vertebrate or non-vertebrate (tunicate) was used as the representative deuterostome.
A major limitation of this study, with respect to phylogenetic implications, is the sparse taxonomic sampling. As in most studies, there is a trade-off in terms of taxa and genes or proteins. In this case, we have emphasized a large number of proteins to increase the statistical resolution of relationships and time estimates at the expense of taxonomic representation. Nonetheless, our results agree with most previous molecular phylogenetic analyses in supporting a close relationship between insects and crustaceans (Pancrustacea) and between myriapods and chelicerates. Of the two groups, Pancrustacea has received the strongest support in the past [20–24]. Nonetheless, a myriapod-chelicerate grouping has been found previously with mitochondrial genes [23–25], nuclear ribosomal genes , and nuclear protein-coding genes . Additional evidence has come from hemocyanine structure .
We propose the name Myriochelata (in allusion to the joining of Myriapoda and Chelicerata) for the group containing myriapods and chelicerates, which otherwise is unnamed . Although we are unable to identify any morphological trait diagnostic of this clade, some trends are evident that might reflect the morphological or ecological nature of the ancestral myriochelatan. For example, many species of extant myriochelatans (e.g., spiders, scorpions, centipedes) inject and immobilize prey with a poison, albeit with structures that are not homologous. Envenomation of prey is also found among pancrustaceans, but it is less broadly distributed in that group. Certainly, envenomation has arisen multiple times in arthropods, associated mostly (but not exclusively) with terrestrial predation. The significance of this trait in arthropod evolution must await sequence evidence from a greater diversity of taxa (e.g., pycnogonids, remipedes) than is currently available, and a careful examination of the early fossil record of animals (especially from the Cambrian). In general, the difficulty in finding morphological characters diagnostic of these major clades of arthropods is probably the result of deep branching of the lineages and an early fossil record that shows great morphological diversification (the Cambrian Explosion) and some important gaps [26, 29, 30].
Timescale for animal evolution and colonization of land
Among animals, arthropods have been considered to be the earliest colonizers of land based on fossil evidence [5–7]. However, it is possible that other animal phyla colonized land even earlier. Among them, nematodes, tardigrades, and annelids are likely candidates given their current exploitation of terrestrial environments, yet these groups have relatively poor fossil records.
Determining the number of such colonization events requires a consideration of phylogeny, the fossil record, and morphological traits associated with terrestriality. For arthropods, at least four major colonization events are inferred, leading to the chilopods + diplopods, insects, arachnids, and isopod crustaceans. Additional events may have occurred in smaller lineages [24, 25, 31]. Moreover, the recent discovery of a basal marine hexapod fossil from the Devonian  suggests that some hexapod traits previously believed to have evolved as adaptations to land may have first appeared in a marine setting. Along the same lines, it is also possible that millipedes and centipedes colonized land independently. However, because the earliest fossils of those groups are presumably terrestrial  and our molecular time estimate is only 5% earlier than the age of those fossils, our assumption of a terrestrial common ancestor of millipedes and centipedes has little affect on the time of colonization.
Our relatively young time estimates for the millipede-centipede (442 ± 50 Ma) and xiphosuran-arachnid (475 ± 53 Ma) divergences contrast with the much older time estimate for the deuterostome-arthropod divergence (~1100-900 Ma) (Tables 1 and 2). The first two are close to the corresponding fossil record estimates whereas the third one is ~400 million years earlier than the earliest fossil evidence for animals . Most molecular clock analyses in the last three decades, including those using many genes, have resulted in similarly deep time estimates for the arthropod-deuterostome divergence . Currently a debate exists as to whether divergences among animals are best represented by molecular clock dates or the fossil record [33–37]. The molecular clock dates suggest that early animals fossilized poorly, possibly because they were small and soft-bodied. Alternatively, others have argued that molecular time estimates are older than the true times because of statistical biases, rate changes, and calibration biases [36, 38], although replies to those criticisms have been made .
The results of this study address several of these criticisms of molecular clocks in the following ways. (1) The time estimates reported here are not uniformly discordant with the fossil record suggesting that if statistical biases are present, they are not causing a directional and proportional bias in all time estimates. (2) Our time estimate for the deuterostome-arthropod divergence using the arthropod fossil record (exclusively) was similar to time estimates obtained in previous studies by using the vertebrate fossil record, suggesting that the vertebrate calibration is not obviously biasing time estimates. (3) Our use of tunicates instead of vertebrates for representing deuterostomes resulted in similar time estimates, further indicating that vertebrates (per se) are not causing a systematic bias in the time analysis. (4) Our use of a diversity of clock methodology, including Bayesian inference and likelihood-smoothing methods, did not alter the conclusions, indicating that the global clock methods used in previous studies were not responsible for biased time estimates. (5) The use of concatenated alignments (supergenes) yielded similar results to non-concatenated (multigene) analyses, indicating that the multigene analyses of previous studies were not responsible for the discordance between molecular clocks and the fossil record.
In summary, our molecular clock analysis resulted in some time estimates (e.g. millipede-centipede) in agreement with the fossil record and others (e.g., insect-crustacean) much earlier than fossil evidence (assuming that the crustaceans are, in fact, monophyletic). However, all studies, including this one, are limited by the small size of the available sequence data (5 – 10 proteins across major groups of arthropods). In the future, it should be possible to use several hundred proteins across a diversity of arthropod taxa, and such an analysis should greatly increase the precision of phylogenetic and molecular clock analyses.
We have found strong statistical support for the unconventional grouping of myriapods and chelicerates, a taxon that we herein name Myriochelata. We also note that many myriochelatans immobilize their prey with venom. Using only arthropod fossil calibrations, our molecular clock analyses support earlier studies that used vertebrate calibrations in finding a deep divergence of deuterostomes and arthropods, hundreds of millions of years before the Cambrian explosion. However, our much younger molecular time estimate for the millipede-centipede divergence is close to the divergence time inferred from fossils for that node. This suggests that the colonization of land by arthropods occurred relatively late in their evolutionary history.
Species and sequences
We sequenced the protein-coding region of two nuclear glycolytic enzymes, glyceraldehyde 3-phosphate dehydrogenase (G3PDH) and enolase, in an annelid (Nereis macrydi), a mollusk (Marisa sp.; not sequenced for enolase), and eight species of arthropods: a centipede (Lithobius sp.), millipede (Diplopoda sp.), blue crab (Callinectes sapidus), brine shrimp (Artemia sp.), water flea (Daphnia magna; not sequenced in G3PDH), muscle shrimp (Ostracoda sp.; possibly Cypridopsis sp.), Atlantic horseshoe crab (Limulus polyphemus), tarantula (Phormictopus sp.), and a scorpion (Centruroides sp.). Reverse-transcriptase PCR was used to amplify and sequence (both complementary strands) a total of approximately 1 100 base pairs (~363 amino acids) of the coding region of enolase and 890 base pairs (~294 amino acids) of the coding region of G3PDH. Additional sequence data were obtained from the public databases (GenBank). Genes were selected if there were nucleotide sequences available from an annelid (outgroup) and at least one representative of the following arthropod taxa: Branchiopoda (e.g., brine shrimp), Malacostraca (e.g., blue crabs), Diplopoda (e.g., millipedes), Chilopoda (e.g., centipedes), Insecta (e.g., insects), Xiphosura (e.g., horseshoe crabs), and Arachnida (e.g., spiders and scorpions). Sequences were aligned with Clustal X . Additional methods and sequence accession numbers are presented elsewhere (Additional file 1).
The alignments of all nine nuclear and 15 mitochondrial genes were concatenated (Additional file 1) because individual gene data sets were, for the most part, insufficient for statistical resolution of arthropod phylogeny (Additional file 2). This required the use of exemplars (representatives) from each of the major groups. Preliminary analyses using 16 – 54 taxa (Additional file 1) were used as a guide to choosing exemplars so that fast- or slow-evolving species were not selected.
Analyses of the concatenation of the nine nuclear genes, as well as a full data set of 24 genes, were carried out using minimum evolution, maximum parsimony, maximum likelihood, and Bayesian inference. Minimum evolution analyses were performed using Kimura two-parameter distances with transversions-only, gamma corrected Kimura two-parameters distances with transversions-only, and paralinear distances. For the paralinear distance analyses, the proportion of invariable sites in the considered alignment was estimated using likelihood and assuming a HKY85 + gamma + proportion of invariable sites model of DNA evolution. Identical sites were then removed proportionally to the base frequencies estimated from all sites (PAUP default settings). Maximum parsimony analyses were carried out with an equal weighting and a 2:1 weighting scheme favoring transversions over transitions. The model used for the maximum likelihood analyses (GTR + gamma+ proportion of invariable sites) was selected using Modeltest , whereas Bayesian inference was carried out under mixed models, and full parameter estimation was performed during tree search for each gene.
All minimum evolution analyses were performed using MEGA 2.1  and PAUP , while likelihood and parsimony analyses, as well as the estimation of the specific gamma parameters for the minimum evolution analyses were carried out using PAUP. Bayesian analyses were performed with MrBayes 3.0 . The branch and bound algorithm was used in the parsimony analyses, while in the likelihood analyses 100 heuristic searches were performed. In the latter, starting trees were obtained using random sequence addition and swapped using the tree bisection reconnection algorithm. In the Bayesian analysis, 300 000 generations were run and trees were used only after convergence was reached.
Robustness for the nodes in the minimum evolution, maximum parsimony, and maximum likelihood trees was evaluated using the bootstrap. PAUP settings for the bootstrap analyses were as follows: 2 000 replicates with full heuristic search (one random addition sequence) and the multiple tree option turned off. Trees were swapped using the TBR algorithm. Support for the groups recovered in the Bayesian analyses were expressed as their posterior probabilities. Gapped sites were removed prior to all analyses.
Molecular clock methods
Molecular clock analyses were performed using a diversity of methods : multigene global least squares (MGGLS) [43, 44], multigene local least squares (MGLLS) , supergene global least squares (SGGLS) , supergene local least squares (SGLLS), supergene local divtime (SGLDT; DivTime5b) , supergene local multidivtime (SGLMDT; Multidivtime) , and supergene local penalized likelihood (SGLPL; r8s version 1.5) [49, 50]. Least squares methods are distance based, SGLPL is a semi-parametric likelihood-smoothing method, and SGLDT and SGLMDT are Bayesian methods. Multigene methods use the mean or mode [37, 44, 51] of time estimates from individual proteins whereas supergene methods derive a single time estimate from the simultaneous analysis of all available proteins. SGLDT and SGLMDT are two different implementations of the Bayesian method of Thorne et al. . The difference between the two being that with SGLDT protein sequences are concatenated and considered as a single entity, whereas with SGLMDT each considered protein maintains its individuality during the analyses.
We used the divergence of centipedes (Chilopoda) and millipedes (Diplopoda) as the minimal (most recent) time of arthropod terrestrialization and the divergence of myriapods with their closest relatives (in this case, chelicerates) as the maximal (earliest) time for myriapod terrestrialization. The maximal time for arthropod terrestrialization was considered as the divergence between insects (Insecta) and crustaceans (Crustacea). This does not preclude the possibility that land was colonized even earlier by extinct species or groups of arthropods (or other animals).
Divergence times were estimated using two calibration points from the arthropod fossil record (Additional file 1): Xiphosura-Arachnida (~480 Ma)  and Chilopoda-Diplopoda (423 Ma; see Additional File 1 for explanation of date) [53, 54]. A third calibration point, the divergence of deuterostomes (vertebrates) and arthropods based on a molecular clock study using 50 proteins and calibrated with the vertebrate fossil record (993 Ma) , was used in some analyses and tested (with the arthropod fossil record and by substituting vertebrates with a tunicate, Ciona intestinalis) in others (see below). Analyses were performed using calibration points separately and simultaneously, for comparison. The millipede-centipede fossil divergence was not used to estimate the molecular clock time of that divergence.
GenBank was screened for nuclear encoded proteins (>100 amino acids) in which divergence times could be estimated between two or more arthropod groups, resulting in 84 proteins. For each, reciprocal BLAST analyses were used to identify and assemble sets of available sequences, and orthology was investigated through phylogenetic analyses. In those analyses, 61 proteins (27 291 amino acid positions) were selected for further analysis after determining that orthologous relationships (divergences corresponding to speciation events) could be distinguished from paralogous relationships (divergences corresponding to gene duplication events).
The least squares-based global and local clock methods (MGGLS, SGGLS, MGLLS, SGLLS) were performed on 49 proteins (Additional file 1) that passed the relative rate test  as implemented in PHYLTEST . Gamma-corrected Poisson distances were used for these analyses and specific gamma parameters were estimated with PAML . The Bayesian and the likelihood-smoothing local clock methods (SGLDT, SGLMDT, and SGLPL) were used with all 61 proteins that passed or did not pass the rate tests. Branch lengths for the penalized likelihood analyses were estimated with PAML (assuming a Gamma corrected JTT model of amino acid substitution). For the Bayesian analyses, branch length estimation was performed using the software AUTOestbranches  in the case of SGLDT, and Estbranches  in the case of SGLMDT. Standard errors for SGGLS and SGLLS were calculated from among-gene comparisons (i.e., MGGLS and MGLLS) and were not available for SGLPL.
Horodyski RJ, Knauth LP: Life on land in the Precambrian. Science. 1994, 263: 494-498.
Watanabe Y, Martini JEJ, Ohmoto H: Geochemical evidence for terrestrial ecosystems 2.6 billion years ago. Nature. 2000, 408: 574-578. 10.1038/35046052.
Prave AR: Life on land in the Proterozoic: Evidence from the Torridonian rocks of northwest Scotland. Geology. 2002, 30: 811-814. 10.1130/0091-7613(2002)030<0811:LOLITP>2.0.CO;2.
Schwartzman DW: Life, temperature, and the Earth. 1999, New York, Columbia University Press, 241-
Shear WA, Bonamo PM, Grierson JD, Rolfe WDI, Smith EL, Norton RA: Early land animals in North America. Science. 1984, 224: 492-494.
Jeram AJ, Selden PA, Edwards D: Land animals in the Silurian: arachnids and myriapods from Shropshire, England. Science. 1990, 250: 658-661.
Gray J, Shear W: Early life on land. American Scientist. 1992, 80: 444-456.
Redecker D, Kodner R, Graham LE: Glomalean fungi from the Ordovician. Science. 2000, 289: 1920-1921. 10.1126/science.289.5486.1920.
Wellman CH, Osterloff PL, Mohluddin U: Fragments of the earliest land plants. Nature. 2003, 425: 282-285. 10.1038/nature01884.
Xiao SH, Zhang Y, Knoll AH: Three-dimensional preservation of algae and animal embryos in a Neoproterozoic phosphorite. Nature. 1998, 391: 553-558. 10.1038/35318.
Li CW, Chen JY, Hua TE: Precambrian sponges with cellular structures. Science. 1998, 279: 879-882. 10.1126/science.279.5352.879.
Barfod GH, Albarede F, Knoll AH, Xiao SH, Telouk P, Frei R, Baker J: New Lu-Hf and Pb-Pb age constraints on the earliest animal fossils. Earth Planet Sc Lett. 2002, 201: 203-212. 10.1016/S0012-821X(02)00687-8.
Porter SM, Knoll AH: Testate amoebae in the Neoproterozoic Era: evidence from vase-shaped microfossils in the Chuar Group, Grand Canyon. Paleobiology. 2000, 26: 360-385.
Butterfield NJ: Bangiomorpha pubescens n. gen., n. sp.: implications for the evolution of sex, multicellularity, and the Mesoproterozoic/Neoproterozoic radiation of eukaryotes. Paleobiology. 2000, 26: 386-404.
Heckman DS, Geiser DM, Eidell BR, Stauffer RL, Kardos NL, Hedges SB: Molecular evidence for the early colonization of land by fungi and plants. Science. 2001, 293: 1129-1133. 10.1126/science.1061457.
Johnson EW, Briggs DEG, Suthren RJ, Wright JL, Tunnicliff SP: Non-marine arthropod traces from the subaereal Ordivician Borrowdale volcanic group, English Lake District. Geol Mag. 1994, 131: 395-406.
MacNaughton RB, Cole JM, Dalrymple RW, Braddy SJ, Briggs DEG, Lukie TD: First steps on land: Arthropod trackways in Cambrian-Ordovician eolian sandstone, southeastern Ontario, Canada. Geology. 2002, 30: 391-394. 10.1130/0091-7613(2002)030<0391:FSOLAT>2.0.CO;2.
Robison RA, Wiley EO: A new arthropod, Meristoma: more fallout from the Cambrian explosion. Journal of Paleontology. 1995, 69: 447-459.
Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
Shultz JW, Regier JC: Phylogenetic analysis of arthropods using two nuclear protein-encoding genes supports a crustacean + hexapod clade. Proc R Soc Lond B Biol Sci. 2000, 267: 1011-1019. 10.1098/rspb.2000.1104.
Friedrich M, Tautz D: Ribosomal DNA phylogeny of the major extant arthropod classes and the evolution of myriapods. Nature. 1995, 376: 165-167. 10.1038/376165a0.
Giribet G, Edgecombe GD, Wheeler WC: Arthropod phylogeny based on eight molecular loci and morphology. Nature. 2001, 413: 157-161. 10.1038/35093097.
Hwang UW, Friedrich M, Tautz D, Park CJ, Kim W: Mitochondrial protein phylogeny joins myriapods with chelicerates. Nature. 2001, 413: 154-157. 10.1038/35093090.
Nardi F, Spinsanti G, Boore JL, Carapelli A, Dallai R, Frati F: Hexapod origins: monophyletic or paraphyletic?. Science. 2003, 299: 1887-1889. 10.1126/science.1078607.
Delsuc F, Phillips MJ, Penny D: Comment on "Hexapod origins: monophyletic or paraphyletic?". Science. 2003, 301: 1482; author reply 1482-10.1126/science.1086558.
Cook CE, Smith ML, Telford MJ, Bastianello A, Akam M: Hox genes and the phylogeny of the arthropods. Curr Biol. 2001, 11: 759-763. 10.1016/S0960-9822(01)00222-6.
Kusche K, Burmester T: Diplopod hemocyanin sequence and the phylogenetic position of the Myriapoda. Mol Biol Evol. 2001, 18: 1566-1573.
Štys P, Zrzavý J: Phylogeny and classification of extant Arthropoda: a review of hypotheses and nomenclature. Eur J Entomol. 1994, 91: 257-275.
Fortey RA, Owens RM: Evolutionary radiations in the Trilobita. Major evolutionary radiations. Edited by: Taylor PD and Larwood GP. 1990, Oxford, Clarendon Press, 139-164.
Lieberman BS: Phylogenetic analysis of some basal early Cambrian trilobites, the biogeographic origins of the Eutrilobita, and the timing of the Cambrian radiation. J Paleontol. 2002, 76: 692-708.
Dunlop JA, Webster M: Fossil evidence, terrestrialization and arachnid phylogeny. J Arachnol. 1999, 27: 86-93.
Haas F, Waloszek D, Hartenberger R: Devonohexapodus bocksbergensis, a new marine hexapod from the Lower Devonian Hunsrück Slates, and the origin of Atelocerata and Hexapoda. Organisms, Diversity and Evolution. 2003, 3: 39-54.
Wang DY, Kumar S, Hedges SB: Divergence time estimates for the early history of animal phyla and the origin of plants, animals and fungi. Proc R Soc Lond B Biol Sci. 1999, 266: 163-171. 10.1098/rspb.1999.0617.
Wray GA, Levinton JS, Shapiro LH: Molecular evidence for deep Precambrian divergences among metazoan phyla. Science. 1996, 274: 568-573. 10.1126/science.274.5287.568.
Fortey RA, Briggs DEG, Wills MA: The Cambrian evolutionary 'explosion' recalibrated. Bioessays. 1997, 19: 429-434.
Benton MJ, Ayala FJ: Dating the tree of life. Science. 2003, 300: 1698-1700. 10.1126/science.1077795.
Hedges SB, Kumar S: Genomic clocks and evolutionary timescales. Trends in Genetics. 2003, 19: 200-206. 10.1016/S0168-9525(03)00053-2.
Benton MJ: Early origins of modern birds and mammals: molecules vs. morphology. Bioessays. 1999, 21: 1043-1051. 10.1002/(SICI)1521-1878(199912)22:1<1043::AID-BIES8>3.3.CO;2-2.
Thompson JD, Higgins DG, Gibson TJ: CLUSTALW: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Research. 1994, 22: 4673-4680.
Kumar S, Tamura S, Jakobsen IB, Nei M: MEGA2: Molecular Evolutionary Genetics Analysis software. 2001, Tempe, AZ, USA, Arizona State University
Swofford DL: PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4. 1998, Sunderland, Massachusetts, Sinauer Associates
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Hedges S. Blair, Parker Patrick H., Sibley Charles G., Kumar Sudhir: Continental breakup and the ordinal diversification of birds and mammals. Nature. 1996, 381: 226-229. 10.1038/381226a0.
Kumar S, Hedges SB: A molecular timescale for vertebrate evolution. Nature. 1998, 392: 917-920. 10.1038/31927.
Schubart Christoph D., Diesel Rudolf, Hedges S. Blair: Rapid evolution to terrestrial life in Jamaican crabs. Nature. 1998, 393: 363-365. 10.1038/30724.
Nei M, Xu P, Glazko G: Estimation of divergence times from multiprotein sequences for a few mammalian species and several distantly related organisms. Proceedings of the National Academy of Sciences (U.S.A.). 2001, 98: 2497-2502. 10.1073/pnas.051611498.
Thorne JL, Kishino H, Painter IS: Estimating the rate of evolution of the rate of molecular evolution. Molecular Biology and Evolution. 1998, 15: 1647-1657.
Thorne JL, Kishino H: Divergence time and evolutionary rate estimation with multilocus data. Systematic Biology. 2002, 51: 689-702. 10.1080/10635150290102456.
Sanderson MJ: Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Molecular Biology and Evolution. 2002, 19: 101-109.
Sanderson MJ: r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003, 19: 301-302. 10.1093/bioinformatics/19.2.301.
Hedges SB, Shah P: Comparison of mode estimation methods and application in molecular clock analysis. BMC Bioinformatics. 2003, 4: 31-10.1186/1471-2105-4-31.
Dunlop JA, Selden PA: The early history and phylogeny of the chelicerates. Arthropod relationships. Edited by: Fortey RA and Thomas RH. 1998, London, Chapman & Hall, 221-235.
Almond JE: The Silurian-Devonian fossil record of the Myriapoda. Phil Trans R Soc Lond B. 1985, 309: 227-237.
Shear WA, Selden PA: Rustling in the undergrowth: animals in early terrestrial ecosystems. Plants invade the land. Edited by: Gensel PG and Edwards D. 2001, New York, Columbia University Press, 29-51.
Takezaki Nauko, Rzhetsky Andrey, Nei Masatoshi: Phylogenetic test of the molecular clock and linearized trees. Molecular Biology and Evolution. 1995, 12: 823-833.
Kumar S: PHYLTEST. 1996, University Park, PA, USA., IMEG, Department of Biology, The Pennsylvania State University, 2.0.
Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. CABIOS. 1997, 13: 555-556.
We thank Frederick Schram for helpful advice concerning arthropod taxonomy and literature. This work was supported by grants to S.B.H. from the National Science Foundation (DBI-0112670) and the National Aeronautics and Space Administration (NCC2-1057).
DP carried out the bioinformatic research and drafted the manuscript. LLP and ML sequenced the new G3PDH and enolase sequences. LLP carried out preliminary analyses of the sequence data. SBH directed the research and assisted with drafting of the manuscript. All authors read and approved the final manuscript.