A comprehensive and integrative reconstruction of evolutionary history for Anomura (Crustacea: Decapoda)
© Bracken-Grissom et al.; licensee BioMed Central Ltd. 2013
Received: 3 December 2012
Accepted: 21 May 2013
Published: 20 June 2013
The infraorder Anomura has long captivated the attention of evolutionary biologists due to its impressive morphological diversity and ecological adaptations. To date, 2500 extant species have been described but phylogenetic relationships at high taxonomic levels remain unresolved. Here, we reconstruct the evolutionary history—phylogeny, divergence times, character evolution and diversification—of this speciose clade. For this purpose, we sequenced two mitochondrial (16S and 12S) and three nuclear (H3, 18S and 28S) markers for 19 of the 20 extant families, using traditional Sanger and next-generation 454 sequencing methods. Molecular data were combined with 156 morphological characters in order to estimate the largest anomuran phylogeny to date. The anomuran fossil record allowed us to incorporate 31 fossils for divergence time analyses.
Our best phylogenetic hypothesis (morphological + molecular data) supports most anomuran superfamilies and families as monophyletic. However, three families and eleven genera are recovered as para- and polyphyletic. Divergence time analysis dates the origin of Anomura to the Late Permian ~259 (224–296) MYA with many of the present day families radiating during the Jurassic and Early Cretaceous. Ancestral state reconstruction suggests that carcinization occurred independently 3 times within the group. The invasion of freshwater and terrestrial environments both occurred between the Late Cretaceous and Tertiary. Diversification analyses found the speciation rate to be low across Anomura, and we identify 2 major changes in the tempo of diversification; the most significant at the base of a clade that includes the squat-lobster family Chirostylidae.
Our findings are compared against current classifications and previous hypotheses of anomuran relationships. Many families and genera appear to be poly- or paraphyletic suggesting a need for further taxonomic revisions at these levels. A divergence time analysis provides key insights into the origins of major lineages and events and the timing of morphological (body form) and ecological (habitat) transitions. Living anomuran biodiversity is the product of 2 major changes in the tempo of diversification; our initial insights suggest that the acquisition of a crab-like form did not act as a key innovation.
KeywordsAnomura Phylogeny Divergence times Diversification rates Molecular Morphology Character evolution Next-generation sequencing
The infraorder Anomura represents a highly diverse group of decapod crustaceans comprised of hermit crabs, mole crabs, king crabs, squat-lobsters and porcelain crabs. The fossil record contains representatives of nearly all extant families and spans the Norian/Rhaetian (Late Triassic)  to Holocene. Anomurans have colonized a wide variety of ecosystems including freshwater, anchialine cave, terrestrial and hydrothermal vent habitats, and are distributed from the ocean’s surface to depths more than 5000 m . Their morphological and ecological diversity are of doubtless scientific interest, but anomurans also represent an important economic commodity as evident in major commercial fisheries for some king crab and squat lobster genera [3–5] and the common use of hermit crabs as pets in the aquarium trade. Moreover, some species are threatened or endangered due to rarity in nature, e.g., Pylochelidae , overfishing, e.g., Lithodidae , or habitat loss, e.g., Aeglidae [8–10]. Thus, improved understanding of these groups bears not only on appreciation of their diversity and ecology, but also strategies for their conservation.
Anomuran classification has long been fraught with controversy [see reviews by [11–13]]. Early classifications from the 19th to the first half of the 20th centuries were based on adult morphological characters including mouthparts, antennae, gills, pleon type, and/or larval characteristics. These classifications often differed in higher-level composition and, in some cases, the infraordinal name (e.g. Anomura vs. Anomala). Since these studies, various researchers have proposed changes in the classification scheme [14–18], many of which remain actively debated. More recently, molecular and/or morphological data have been used to reevaluate anomuran relationships [19–21]. As currently defined, extant Anomura contains 7 superfamilies, 20 families, 335 genera, and more than 2500 species [17, 18, 22, 23]. Although the monophyly of Anomura is widely accepted [24–26], the elucidation of internal relationships among families, genera, and species using modern methods is dynamic and under continuous debate [11, 17, 18, 20, 23, 27].
One of the most debated evolutionary questions within Anomura is phylogenetic relationships between hermit and king crabs. Since the early 1800’s [e.g., [28, 29]], studies have suggested king crabs and hermit crabs are close relatives, despite first appearances to the contrary. King crabs are among the largest arthropods and have a crab-like body shape, whereas hermit crabs are relatively small and depend on a shell for protection. Despite glaring morphological differences as adults, an affinity between king crabs (lithodoids) and hermit crabs (paguroids) has been long suggested [30, 31]. Although most accept this claim, the evolutionary pathways and hypothesized ancestor of both groups has been debated for decades, with two major hypotheses being proposed. The first suggests that the lithodids (Lithodes or Paralithodes) evolved from a pagurid-like ancestor (Pagurus) (“hermit to king hypothesis”) while the second suggests the opposite evolutionary pathway (“king to hermit hypothesis”). Here we revisit these hypotheses in light of new phylogenetic data to test the “hermit to king”/“king to hermit” evolutionary pathway.
Additional controversy over anomuran relationships stems from apparently rampant examples of convergent and/or parallel evolution in body forms. Anomurans span an impressive array of body configurations that include: 1) crab-like forms 2) squat-lobster forms 3) hermit crab forms with pleonal (abdomen) symmetry (found in 1 hermit crab family) and 4) hermit crab forms with pleonal asymmetry (found in 4 hermit crab families). Recent studies suggest that the acquisition of a crab-like body form, known as carcinization [see,  for a review of this concept], has occurred multiple times during evolution of the group [20, 33]. Brachyura, all of which possess a “crab-like” body shape or slight modifications to this form, dominates decapod diversity with 6550+ species , and is considered the sister clade to Anomura [24–26, 35–38]. Given the success of brachyuran crabs, it has been hypothesized that acquisition of a crab-like form may have acted as a key innovation , possibly impacting diversification rates within these lineages. For the first time, we explore diversification patterns in Anomura and specifically test if carcinized lineages underwent unusually rapid diversification rates. If the emergence of the crab-like form promoted diversification we would expect the overall rate in carcinized lineages to be high compared to net of diversification across Anomura. Additionally, we test if the acquisition of different body forms (i.e., crab-like, squat-lobster-like, pleonal asymmetry and symmetry (hermit)) arose once or multiple times during the emergence of the anomurans and reconstruct the evolutionary pathways of these transitions.
Divergence dating is a powerful tool used to estimate the timing and origins of diversity, morphological traits, habitat shifts, and diversification. Although nearly all the family-level groups of Anomura are represented in the fossil record, the discovery has not been as frequent as that of other decapod groups (i.e., true crabs, lobsters). Two factors, variations in cuticular sclerotization and habitat preference, are likely responsible for the limited occurrence of anomuran fossils. Many taxa are weakly calcified, whereas others possess well-calcified claws and poorly calcified carapaces and pleons. In addition, habitats currently occupied by anomurans, including freshwater, terrestrial, intertidal marine, deep marine, and hydrothermal vent areas are strongly underrepresented in the fossil record. Despite these limitations, we incorporate 31 fossil calibrations to estimate the origin of lineages and major events during anomuran evolutionary history, including the transition of body forms and shift into freshwater and terrestrial environments.
Here, we present the taxonomically broadest and largest dataset yet assembled. We combine sequences generated by traditional Sanger and next-generation 454 sequencing methods with morphological characters, including 19/20 extant families and 137 species, to estimate phylogenetic relationships, character state evolution, divergence times, and diversification patterns among major lineages of this diverse clade of crustaceans. Our comprehensive sampling, in combination with modern integrative approaches, allows us to present the most complete evolutionary picture for the infraorder Anomura to date.
Taxonomy, voucher catalog numbers, and GenBank accession numbers for gene sequences used in this study
Aegla abtao Schmitt, 1942
Aegla alacalufi Jara & Lopez, 1981
Aegla camargoi Buckup & Rossi, 1977
Aegla cholchol Jara & Palacios, 1999
Aegla jarai Bond-Buckup & Buckup, 1994
Aegla papudo Schmitt, 1942
Aegla platensis Schmitt, 1942
Aegla uruguayana Schmitt, 1942
Aegla violacea Bond-Buckup & Buckup, 1994
Albunea catherinae Boyko, 2002
Albunea gibbesii Stimpson, 1859
Lepidopa californica Efford, 1971
Lepidopa dexterae Abele & Efford, 1972
Paraleucolepidopa myops (Stimpson, 1860)
Zygopa michaelis Holthuis, 1961
Blepharipoda occidentalis Randall, 1840
Chirostylus novaecaledoniae Baba, 1991
Gastroptychus novaezelandiae Baba, 1974
Gastroptychus rogeri Baba, 2000
Gastroptychus spinifer (A. Milne-Edwards, 1880)
Uroptychus nitidus (A. Milne-Edwards, 1880)
Uroptychus parvulus (Henderson, 1885)
Uroptychus scambus Benedict, 1902
Uroptychus spinirostris (Ahyong & Poore, 2004)
Birgus latro (Linnaeus, 1767)
Coenobita clypeatus (Fabricius, 1787)
Coenobita compressus (H. Milne Edwards)
Coenobita perlatus H. Milne Edwards, 1836
Areopaguristes hewatti (Wass, 1963)
Areopaguristes hewatti (Wass, 1963)
Areopaguristes hewatti (Wass, 1963)
Areopaguristes hummi (Wass, 1955)
Areopaguristes hummi (Wass, 1955)
Areopaguristes pilosus H. Milne Edwards, 1836
Calcinus laevimanus (Randall, 1840)
Calcinus obscurus Stimpson, 1859
Clibanarius albidigitus Nobili, 1901
Clibanarius antillensis Stimpson, 1859
Clibanarius corallinus (H. Milne-Edwards, 1848)
Clibanarius vittatus (Bosc, 1802)
Dardanus fucosus Biffar & Provenzano, 1972
Dardanus insignis (de Saussure, 1858)
Isocheles pilosus (Holmes, 1900)
Isocheles wurdemanni Stimpson, 1859
Paguristes cadenati Forest, 1954
Paguristes grayi (Benedict, 1901)
Paguristes nr. moorei
Paguristes punticeps Benedict, 1901
Paguristes sericeus A. Milne Edwards, 1880
Paguristes tortugae Schmitt, 1933
Paguristes triangulatus A. Milne-Edwards & Bouvier, 1893
Paguristes turgidus (Stimpson, 1857)
Petrochirus diogenes (Linnaeus, 1758)
Eumunida funambulus Gordon, 1930
Eumunida picta Smith, 1883
Eumunida picta Smith, 1883
Pseudomunida fragilis Haig, 1979
Alainius crosnieri Baba, 1991
MNHN:Norfolk I Stn DW 1703
Galathea rostrata A. Milne-Edwards, 1880
Hapalogaster mertensii Brandt, 1850
Oedignathus inermis (Stimpson, 1860)
Emerita brasiliensis Schmitt, 1935
Emerita emeritus (Linnaeus 1767)
Emerita talpoida (Say, 1817)
Kiwa hirsuta Macpherson, Jones & Segonzac, 2005
Glyptolithodes cristatipes (Faxon, 1893)
Lithodes santolla (Molina, 1782)
Lithodes santolla (Molina, 1782)
Lopholithodes mandtii Brandt, 1848
Paralithodes brevipes (H. Milne Edwards & Lucas, 1841)
Paralithodes camtschaticus Tilesius, 1815
Paralithodes platypus (Brandt, 1850)
Phyllolithodes papillosus Brandt, 1848
Lomis hirta (Lamarck, 1818)
Agononida procera Ahyong & Poore, 2004
Anoplonida inermis (Baba, 1994)
MNHN:SANTO Stn AT9
Babamunida kanaloa Schnabel, Martin & Moffitt, 2009
Bathymunida balssi Van Dam, 1838
MNHN:SANTO Stn AT5
Cervimunida johni Porter, 1903
Munida iris A. Milne-Edwards, 1880
Munida pusilla Benedict, 1902
Munida quadrispina Benedict, 1902
Munida gregaria (Fabricius, 1793)
Neonida grandis Baba & de Saint Laurent, 1996
MNHN:Lifou Stn CP2
Pleuroncodes monodon (H. Milne Edwards, 1837)
Galacantha rostrata A. Milne-Edwards, 1880
Galacantha valdiviae Balss, 1913
Leiogalathea laevirostris (Balss, 1913)
Munidopsis bairdii (Smith, 1884)
Munidopsis erinacea (A. Milne-Edwards, 1880)
Shinkaia crosnieri Baba & Williams, 1998
NTOU:Chan, Lee & Lee, 2000
Agaricochirus alexandri (A. Milne-Edwards & Bouvier, 1893)
Bythiopagurus macrocolus McLaughlin, 2003
Discorsopagurus schmitti (Stevens, 1925)
Goreopagurus piercei Wass, 1963
Iridopagurus caribbensis (A. Milne-Edwards & Bouvier, 1893)
Iridopagurus reticulatus García-Gómez, 1983
Labidochirus splendescens (Owen, 1839)
Manucomplanus ungulatus (Studer, 1883)
Pagurus bernhardus (Linnaeus, 1758)
Pagurus brevidactylus (Stimpson, 1859)
Pagurus bullisi Wass, 1963
Pagurus maclaughlinae García-Gómez, 1982
Pagurus nr. carolinensis
Pagurus pollicaris Say, 1817
Pagurus stimpsoni (A. Milne-Edwards & Bouvier, 1893)
Phimochirus holthuisi (Provenzano, 1961)
Phimochirus randalli (Provenzano, 1961)
Phimochirus randalli (Provenzano, 1961)
Porcellanopagurus filholi de Saint Laurent & McLaughlin, 2000
Pylopaguridium markhami McLaughlin & Lemaitre, 2001
Pylopagurus discoidalis (A. Milne-Edwards, 1880)
Tomopagurus merimaculosus (Glassell, 1937)
Xylopagurus cancellarius Walton, 1950
Parapagurus latimanus Henderson, 1888
Sympagurus acinops Lemaitre, 1989
Sympagurus dimorphus (Studer, 1883)
Sympagurus pictus Smith, 1883
KC7247/KC6835 (ULLZ 10849 16S/H3only)
Allopetrolisthes spinifrons (H. Milne Edwards, 1837)
Megalobrachium poeyi (Guérin-Méneville, 1855)
Neopisosoma angustifrons (Benedict, 1901)
Pachycheles ackleianus A. Milne-Edwards, 1880
Pachycheles haigae Rodrigues da Costa, 1960
Pachycheles pilosus (H. Milne Edwards, 1837)
Pachycheles rudis Stimpson, 1859
Pachycheles rugimanus A. Milne-Edwards, 1880
Parapetrolisthes tortugensis (Glassell, 1945)
Parapetrolisthes tortugensis (Glassell, 1945)
Petrolisthes armatus (Gibbes, 1850)
Petrolisthes armatus (Gibbes, 1850)
Petrolisthes galathinus (Bosc, 1802)
Petrolisthes laevigatus (Guérin, 1835)
Pisidia magdalenensis (Glassell, 1936)
Polyonyx gibbesii Haig, 1956
Porcellana sayana (Leach, 1820)
Pomatocheles jeffreysii Miers, 1879
Trizocheles spinosus (Henderson, 1888)
Xylocheles macrops Forest, 1987
Calappa gallus (Herbst, 1803)
Cyclograpsus cinereus (Dana, 1851)
Praebebalia longidactyla (Yokoya, 1933)
Chorilia longipes (Dana, 1852)
Cosmonotus grayi (White, 1848)
Lepidophthalmus louisianensis (Schmitt,1935)
Sergio mericeae (Manning & Felder, 1995)
Calaxius manningi Kensley et al., 2000
Calastacus crosnieri Kensley & Chan, 1998
Laomedia astacina de Haan, 1841
Thalassina anomala (Herbst, 1804)
Austinogebia narutensis (Sakai, 1986)
Hymenopenaeus debilis Smith, 1882
KC 4444/ULLZ 8531
KC 4454/ULLZ 6705
Latreutes fucorum (Fabricius, 1798)
KC 4498/ULLZ 9135
KC 4517/ULLZ 9174
Palaemonetes pugio Holthuis, 1949
KC 4542/ULLZ 7755
Alternative hypotheses regarding monophyly of the families Paguridae, Diogenidae, Hapalogastridae, Lithodidae, Munididae, Pylochelidae, and the superfamily Paguroidea [14–18, 23, 39] were tested using the Shimodaira-Hasegawa test (S-H). Three of the seven hypotheses were found to be significantly worse than our unconstrained topology (P < 0.05; MLbest = −68420.363272; MLDiogenidae = −68667.853268, MLPaguridae = −68497.123254; MLPaguroidea = −68825.722919). The remaining four hypotheses were not found to be significantly worse than our unconstrained topology (P > 0.05; MLHapalogastridae = −68432.438825; MLLithodidae = −68438.309801; MLMunididae = −68428.284597; MLPylochelidae = −68430.951016). Hypotheses that tested a “king to hermit” evolutionary pathway were all significantly worse than the alternative (i.e., “hermit to king”) as recovered in our best ML tree (P < 0.05; MLbest = −68420.363272; MLking-Paguroidea = −68777.179402; MLking-Paguridae = −68713.171227).
Divergence time analysis
Fossil calibrations used in BEAST divergence time analyses
Geological Age (MYA)
Aciculopoda mapesi Feldmann and Schweitzer, 2010
Late Devonian (Fammenian) 359-374
Pinnacaris dentata Garassino and Teruzzi, 1993
Late Triassic (Norian) 204-228
Palaeopalaemon newberryi Whitfield, 1880
Late Devonian 354-370
Callianassa s.l. bonjouri Étallon, 1861
Early Jurassic (Toarcian) 176-183
Upogebia s. l. obscura von Meyer, 1834
Early Triassic 245-251
Platykotta akaina Chablais, Feldmann and Schweitzer, 2011
Late Triassic (Norian/Rhaetian) 201.6-228
Protaegla miniscula Feldmann, Vega, Applegate, and Bishop, 1998
Early Cretaceous (Albian) 99.6-112
Pristinaspina gelasina Schweitzer and Feldmann, 2001
Late Cretaceous 65.5-99.6
Galatheites zitteli (Moericke, 1889)
Late Jurassic (Tithonian) 145.5-151
Shinkaia katapsyxis Schweitzer and Feldmann, 2008
Juracrista perculta Robins, Feldmann, and Schweitzer, 2012
Late Jurassic (Tithonian) 145.5-151
Munida konara Schweitzer and Feldmann, 2000
Sadayoshia pentacantha (Muller and Collins, 1991)
Late Eocene (Priabonian) 33.9-37.2
Based upon a drawing and description only, type material apparently destroyed: Palaeomunidopsis moutieri Van Straelen, 1925
Middle Jurassic (Bathonian) 168-165
Based upon actual fossil material: Gastrosacus wetzleri Von Meyer, 1851
Late Jurassic (Oxfordian-Tithonian) 161-145
Munidopsis foersteri Feldmann et al., 1993
Late Cretaceous (Campanian) 70.6-83.5
Jurellana tithonia Schweitzer and Feldmann, 2010
Late Jurassic (Tithonian) 145.5-151
Pachycheles dorsosulcatus Beschin, Busulini, De Angeli, and Tessier, 2007
Petrolisthes bittneri De Angeli and Garassino, 2002
Pisidia dorsosinuata De Angeli and Garassino, 2002
Praealbunea rickorum Fraaije, 2002
Late Cretaceous (Maastrichtian) 65.5-70.6
Albunea cuisiana Beschin and De Angeli, 1984
Zygopa galantensis De Angeli and Marangon, 2001
Lophomastix antiqua Schweitzer and Boyko, 2000
Paralomis debodeorum Feldmann, 1998
Based upon claws only: Palaeopagurus deslongchampsi Van Straelen, 1925
Early Jurassic (Pliensbachian) 190-183
Based upon carapace material: Diogenicheles theodorae Fraaije et al., 2012
Late Jurassic (Oxfordian) 161-156
Birgus latro Linnaeus, 1767
Annuntidiogenes ruizdegaonai Fraaije et al., 2008
Early Cretaceous (Albian) 99.6-112
Calcinus agnoensis Beschin et al., 2005
Pagurus malloryi Schweitzer and Feldmann 2001
Eotylaspis wehnerae van Bakel et al., 2008
Late Jurassic (Kimmeridgian) 151-156
Jurapylocheles malutka, Ammopylocheles mclaughlinae Van Bakel et al. 2008
Late Jurassic (Kimmeridgian) 151-156
The reworked version of MEDUSA  was used to detect whether any clade within the anomuran tree was best explained by independent diversification models, and to specifically address whether acquisition of the crab-like form resulted in an increase of diversification rates. The background tempo of diversification across the anomuran tree is characterized by a speciation rate lambda of 0.032572 lineages/Myr, and our results suggest that the diversification of anomurans is characterized by two periods where the tempo of diversification changes (Figure 4). A slow speciation rate is detected in the lineage leading to the monotypic and carcinized family Lomisidae, and an increase rate occurred in the squat-lobster family Chirostylidae. The ancient but species-depauperate branch leading to the monotypic family Lomisidae was optimally modelled separately with maximum likelihood estimate of lambda = 0 (rate reduction). The rate shift that occurred in the branch leading to the family Chirostylidae was characterized by a speciation rate lambda of 0.054182 (rate acceleration). All three resulting clade-specific diversification models were optimally fit as Yule models (AIC = 339.3032).
Recent studies on anomuran evolution have used molecular data [20, 21, 25, 38], morphological (including developmental) data [41–43], or a combination of the two  to resolve phylogenetic relationships. These studies have dramatically increased our understanding of anomuran relationships and resulted in several major changes within higher-level classification [17, 18, 27]. The instability of anomuran taxonomy in recent years highlights the need for continued phylogenetic study of this group at many levels, and we for now elect to follow the most recent and up-to-date classification scheme [17, 18, 22, 44–46].
Our total evidence approach combines 3669 molecular (nuclear and mitochondrial) and 156 morphological (adult, sperm and larval) characters from 137 species to recover the anomuran phylogeny (Figure 1). The addition of morphological data increased the support for many intra-familial and superfamily relationships that were poorly supported in the molecular-only phylogeny (Figure 2). As mentioned previously, anomurans have undergone dramatic changes in higher-level classification based on recent phylogenetic studies. Galatheoidea has been revised recently to exclude Aeglidae, Kiwaidae, and Chirostylidae [18, 23], and include only Galatheidae, Munididae, Munidopsidae, and Porcellanidae . With the recent revision of Galatheoidea, all superfamilies were recovered as monophyletic (i.e., Hippoidea, Aegloidea, Lomisoidea, Chirostyloidea, Galatheoidea, Lithodoidea), except for Paguroidea (Figures 1 and 2). We found Lithodoidea to be nested within Paguroidea, which is in accordance with all recent combined (molecular + morphology) and molecular-based phylogenetic studies [19–21, 33, 41]. An affinity between certain Lithodidae (Lithodes) and Paguridae (Pagurus) has been suggested since the early 1800’s [see  for review of literature], based on morphological characters including mouthparts, gills, and pleonal characters. However, the evolutionary pathways of the two groups continue to be debated (see also “Hermit to King, King to Hermit Evolutionary Hypotheses”) with all recent evidence pointing to a “hermit to king” hypothesis.
Family-level relationships were well resolved in the combined analysis (Figure 1) and in accordance with recent changes in classification [17, 19]. In 2010, Eumunida and Pseudomunida were removed from Chirostylidae and included in the newly erected Eumunididae, and the new family Munididae was erected on the basis of morphological and molecular evidence [17, 18]. Our results generally support these taxonomic revisions, recovering the Eumunididae as a monophyletic group, but finding Munididae to be possibly paraphyletic (Figure 1). Galatheidae was found nested inside Munididae, but alternative topologies that recovered Munididae as monophyletic were not significantly worse than our best estimate (see Results). Deeper sampling within both families is needed to resolve family and genus level relationships. The families Galatheidae, Munidopsidae, and Porcellanidae were all recovered as monophyletic with high support (Figure 1). The paraphyly and/or polyphyly of Diogenidae and Paguridae is consistent across the combined and molecular phylogenies and in accordance with recent phylogenies that have sampled sufficiently within these families [19–21, 25, 26, 41]. Alternative hypotheses proposing the monophyly of these families (i.e., Diogenidae, Paguridae) were rejected using S-H tests, confirming our findings (see Results). Coenobitidae (semi-terrestrial hermit crab) was deeply nested within Diogenidae (left-handed hermit crabs) while Paguristes and Areopaguristes are more closely related to Parapaguridae (deep-water hermit crabs) than to other members of Diogenidae (Figure 1). This relationship was first proposed by Boas , which he collectively called the Paguristinen. The families Pylochelidae and Hapalogastridae were found to be polyphyletic in the molecular analysis (Figure 2), but monophyletic in the combined analysis (Figure 1). Although we did not find Pylochelidae to be polyphyletic in our combined tree, alternative molecular-based and morphological phylogenies have recovered similar results that suggest a polyphyletic Pylochelidae [20, 48]. Additionally, there is morphological support for polyphyly among pylochelids separating Trizocheles and Mixtopagurus from the remaining pylochelid genera (based on form of ocular acicles, eye type and larval forms [20, 48]).
Generic relationships within Anomura seem to be much less resolved than superfamily and family level relationships. We found several genera to be poly- or paraphyletic (i.e., Munida, Munidopsis, Paguristes), in agreement with previous studies [19, 20, 41, 49]. Most instances of non-monophyly occur within highly speciose genera (i.e., Paguristes = ~115 spp., Pagurus = ~170 spp., Munida = ~240 spp.), suggesting deeper sampling and continued research needs to be undertaken on these groups.
King to hermit and hermit to king evolution: historical to recent hypotheses
Although past studies have shown an affinity between Paguridae (hermit crabs) and Lithodidae (king crabs), the evolutionary pathways and ancestry of these anomuran lineages have been debated for the past two centuries. The traditional and prevalent hypothesis posits that lithodids are free-living hermit crabs that abandoned shell use and underwent a series of morphological changes (carcinization) resulting in a crab-like form. It has been argued that the asymmetry of the lithodid female pleon, in particular, is evidence of asymmetrical hermit crab ancestry. Boas [31, 50] was the first to suggest the evolution of lithodids specifically from pagurid ancestors, and based on morphology proposed the ancestral pagurid to be closely related to Nematopagurus and Pylopagurus. Bouvier [51–53] similarly derived the lithodids from the pagurids, agreeing with Boas on the structural pleonal similarities between these two groups. However, Bouvier also proposed a series of gradual and linear progressive stages in the transformation of the pagurid pleon, starting from a pagurid precursor to various genera of hapalogastrids (Hapalogaster, Dermaturus) and lithodids (from Neolithodes, Paralithodes, Lithodes, Lopholithodes, Paralomis, Rhinolithodes, to Cryptolithodes). In modern times, this concept of pagurid and lithodid evolution was brought to attention when Cunningham et al.  coined the phrase “from hermit to king” in applying molecular analysis to study hermit crab and lithodid phylogeny, and was then widely popularized . A subsequent morphologically-based phylogenetic study by Richter and Scholtz  supported this same evolutionary view of pagurid and lithodid evolution. Recently, a study that examines the hemolymph vascular system in hermit and king crabs found close similarities in arterial systems of the dorsal cephalothorax .
An alternate, opposite view, often stated as the “king to hermit” evolutionary hypothesis, was proposed by McLaughlin and Lemaitre . Using morphological characters and an unusual application of cladistic methodology, McLaughlin and Lemaitre explored possible evolutionary pathways of carcinization across Anomura (rather than attempt to determine precise phylogenetic relationships among taxa). They acknowledged that the crab-like form might have arisen multiple times across the Anomura, but in the case of pagurid/lithodid evolution they concluded that the opposite evolutionary trajectory was more plausible, i.e., the transition was from a “crab-like” body form to a ‘hermit-crab” body form through a series of habitat change, calcium loss, and consequential adult morphological adaptations. Subsequent studies showed that the linear evolutionary scenario proposed by Bouvier did not correspond to the ontogenetic changes that take place in the megalopa to juvenile crab stages in at least 10 species of eight lithodid genera [43, 58]. Based on observations of the complex changes in pleonal tergites from megalopa to juvenile crab stages, these studies demonstrated that adult lithodid pleonal tergite structure in several species was the result of decalcification and sundering, not secondary calcification and fusion as had been proposed by Bouvier.
Our recent phylogenetic reconstruction of anomurans based on molecular and morphological data supports the traditional “hermit to king” hypothesis in congruence with all recent studies [19–21, 33, 59]. With the largest number of taxa and most robust molecular/morphological dataset ever used in a phylogenetic study of anomurans, our study once again shows Lithodoidea to be nested within Paguridae. Moreover, our conclusions are consistent with the fossil record, which suggest hermits are much older (Jurassic) than king crabs (Miocene, Table 2). Finally, topology testing rejects the “king to hermit” hypothesis and finds it as significantly worse than the alternative (P < 0.05) (i.e., “hermit to king”) (see Results).
While there is undeniable evidence of a close relationship between hermits and king crabs, it is less clear how morphological changes associated with carcinization may have proceeded within the Lithodoidea. A recent study comparing hermit and king crab circulatory systems identified several vascular changes that occurred as the result of carcinization, arguing for more comparative studies that look at morphology (both internal and external) and development . However, only with a clear phylogenetic hypothesis can many of these studies be correctly interpreted. Recent molecular or combined morphological-molecular phylogenies recover conflicting evolutionary relationships, although only three lithodoid genera (and not always the same, or excluding Hapalogastridae) have been used in previous analyses [19–21, 33]. Our phylogenetic reconstruction (Figure 3) shows the less carcinized and less calcified Hapalogastridae as sister to Lithodidae, in agreement with virtually every study since Bouvier’s in the 19th century. But within Lithodidae, and in contrast to Bouvier’s linear hypothesis, our study places Cryptolithodes, the most heavily calcified and carcinized lithodid, as an early branching lineage followed by more derived genera (see also McLaughlin and Lemaitre, 1997, Figure 2). It thus appears that the process of heavy calcification may have appeared at least twice in lithodid lineages. More lithodoid genera/species are needed to examine the process of carcinization within the Lithodoidea and to properly test Bouvier’s and Boas’ earlier hypotheses (explaining the transition of a shell-dwelling hermit crab to a fully calcified lithodid crab). In conclusion, while recent, modern studies, including ours, overwhelmingly and clearly support a “hermit to king” evolutionary scenario, it is also clear that the evolutionary process and concomitant morphological changes (particularly in pleonal tergal plates and pleopods) that occurred within the Lithodoidea to produce the various degrees of crab-like forms in that family, is at best poorly understood.
In our reaffirmation here of the “hermit to king” hypothesis, we revealed a close relationship between Lithodoidea and the pagurid, Discorsopagurus (Figures 1, 2, 3 and 4). Curiously, the same close relationship has surfaced in previous studies [21, 33]. This revelation is important to highlight because the “hermit [Paguridae] to king [Lithodidae]” hypothesis presupposes a distinctly asymmetrical shell-dwelling hermit crab-like ancestor from or close to the Paguridae, or more precisely Pagurus, as proposed by early  as well as modern studies . However, Pagurus is currently a taxonomic and paraphyletic conundrum of more than 160 species, and it remains unknown which of the different lineages within “Pagurus” could be the most likely candidate for lithodoid ancestry. The close relationship between Discorsopagurus and Lithodoidea may suggest a Discorsopagurus-like hermit crab as the precursor to the crab-like lithodoids. All species of Discorsopagurus are tube-dwellers, not shell-dwellers, and show pleonal asymmetry only in having unpaired pleopods. The genus is relatively small in size compared to the typically large-sized lithodoids with a distribution across both sides of the North Pacific, from the Sea of Japan to Puget Sound and the Straits of Juan de Fuca, Washington . The relationship between Discorsopagurus and lithodoids may not be coincidental in the North Pacific region where D. schmitti (this analysis) and all other Discorsopagurus species are found [61–63]. This region harbors the highest diversity of lithodoids, so it is plausible to expect closely related species (Discorsopagurus) in similar areas. Future studies with increased sampling within these groups will shed light into the evolutionary pathway of lithodoids from paguroid (possibly Discorsopagurus-like) ancestors.
Divergence times and character evolution
Our divergence dating analysis estimated the origin of Anomura to be in the Late Permian (~259 MYA) from a symmetrical crab-like ancestor (Figures 3 and 4). This is consistent with many higher-level decapod phylogenies finding Anomura and Brachyura as sister clades [26, 35, 37, 56], including the present study (Figure 1). Results estimate that the earliest diverging anomurans are the hippoids (~221 MYA), consistent with recent molecular estimates . Although this date is considerably older than the hippoid fossil record, closely related extinct forms extend into the Triassic and present day Hippoidea are found in substrates underrepresented in the fossil record. The superfamily Hippoidea containing Blepharipodidae, Hippidae, and Albuneidae, has been described as being similar to primitive brachyurans [20, 64], and ancestral reconstruction analysis confirms that the present day hippoids were derived from crab-like (brachyuran-like) predecessors (Figure 3A). The next radiation occurred in the Late Triassic, giving rise to the squat-lobsters and crab-like superfamilies Chirostyloidea and Galatheoidea, Aegloidea, Lomisoidea, and the hermit crab and crab-like superfamilies Paguroidea and Lithodoidea. Our results suggest these superfamily clades were derived from a squat-lobster-like ancestor approximately ~205 MYA (Figures 3A and 4). Interestingly, our divergence time and character reconstruction analyses (Figures 3A and 4) are consistent with fossil evidence, and more specifically, the discovery of Platykotta akaina, the oldest known anomuran fossil . Platykotta akaina, with a possibly squat-lobster-like body form, dates back to the Late Triassic (~201.6-228 MYA) and has strong morphological affinity with the superfamilies Chirostyloidea and Galatheoidea. This fossil was found as part of a biotic assemblage suggesting that Platykotta akaina thrived in tropical-subtropical waters and lived in the subtidal with connections to the open ocean [1, 65].
Around 137 MYA a squat-lobster like ancestor gave rise to a unique superfamily of anomurans, Aegloidea. Aegloid crabs represent the only freshwater anomuran family and can be found in caves, lakes, and streams throughout the Neotropical region of South America . Apart from a single species of freshwater hermit crab, Clibanarius fonticola, the transition into a completely freshwater environment only occurred in extant Aeglidae (Figure 3B). Fossil evidence suggests freshwater aeglids once inhabited marine waters, based on the fossil representative, Haumuriaegla glawssneri, found in New Zealand from Late Cretaceous rocks . In combination with our divergence time analyses, we hypothesize that the complete transition in freshwater occurred sometime between the Late Cretaceous and Miocene. This transition appears to have allowed for rapid diversification approximately 13 MYA (20–7.4 MYA).
From approximately 180 MYA to 147 MYA, the families of Galatheoidea radiated and diversified. These include the squat lobsters families Munidopsidae, Munididae and Galatheidae, and the porcelain crab family Porcellanidae. The porcellanids diverged in the Middle Jurassic (~172 MYA) from squat-lobster like ancestors, but a crab-like body form evolved by the Tithonian (~151-145.5 MYA) based on fossil evidence and ancestral reconstruction analyses. This was the first occurrence of carcinization from a squat-lobster or hermit-like ancestor within Anomura (Figures 3A and 4). Interestingly, Henderson  and Ortmann  suggested porcellanid crabs were derived galatheids despite the differences in body shape and form, and this is consistent with our current evolutionary hypothesis.
Lomisoidea and Chirostyloidea diverged around 122 MYA from a squat-lobster like ancestor. This body form was retained within the chirostyloids and underwent further carcinization, attaining a crab-like form in the monotypic Lomisidae, endemic to Australia.
Early hypotheses based on larval evidence proposed hermit crabs evolved as two independent lineages, Coenobitoidea and Paguroidea  and recent studies continue to explore superfamily and family level relationships [20, 23, 39]. In our combined analysis, the hermit crab families, Pylochelidae, Parapaguridae, Diogenidae, Coenobitidae, and Paguridae, formed a monophyletic group with the inclusion of Lithodidae or king crabs, and Hapalogastridae. We estimated these families arose early in the evolution of Anomura, approximately 205 MYA. The symmetrical hermit crabs, Pylochelidae, are unique with most having complete body symmetry and in utilizing broken gastropod shells, siboglinid tubes, and coral pieces for shelter and protection, in contrast to other hermit groups that commonly use coiled gastropod shells . Our analysis suggests pylochelids branched early in the evolution of hermit and king crabs, consistent with morphological and fossil evidence that place them as a basal primitive lineage . The oldest hermit crab fossils, Jurapylocheles malutka, Ammopylocheles mclaughlinae and Eotylaspis wehnerae of Kimmeridgian age (~151-156 MYA), belong to the families Pylochelidae and Parapaguridae (Table 2). This is consistent with our divergence time analysis, which recovers these families as early branching lineages. Diogenidae, Coenobitidae, and Paguridae typically possess an asymmetrical pleon accompanied by an enlarged right or left chela. According to our combined analysis, pleonal asymmetry in hermits appears to have been derived once in the evolution of the anomurans, most probably between 200–187 MYA. This contrasts with the results obtained by Tsang et al. , who proposed that the pleonal asymmetry evolved independently in two different hermit crab lineages, once in Parapaguridae, and a second time in Diogenidae, Coenobitidae, and Paguridae. These contrasting differences are the result of incongruent phylogenies based on total evidence (molecular + morphology, this paper) and molecular only approaches . Note, however, that our molecular-only analyses recover similar results to those of Tsang et al. . The semiterrestrial coenobitids colonized land from a marine ancestor sometime between 84 and 39 MYA (Figures 3B and 4). The emergence of Diogenidae (~173-167 MYA), Coenobitidae (~84 MYA), and Paguridae (~173 MYA) all predate their first appearance in the fossil record (Table 2, Figure 4). Carcinization occurred for the third time in the crab-like superfamily Lithodoidea between 29–18 MYA from an asymmetrical hermit-like ancestor. This estimation is consistent with other timing estimates of king crab carcinization .
The crab-like body form was recovered in our study as the ancestral state for all the anomurans. In our study, all alternative body forms were present (crab-like, squat lobster, symmetrical hermit, and asymmetrical hermit) early in the divergence of the anomurans. From these ancestral character states, carcinization occurred independently three times during the evolution of Anomura, once in the Lithodoidea through an asymmetrical hermit intermediate, and twice in Lomisidae and Porcellanidae through squat lobster intermediates (see AHI and SI, Figure 3A). These evolutionary pathways of the crab-like form, twice from squat lobster intermediaries and once through an asymmetrical ancestor, corroborates recent hypotheses . However, our tree differs significantly from Tsang et al.’s study  in the deep ancestral origins of carcinization. Tsang et al.’s hypothesis suggests a symmetrical hermit crab-like ancestor predated the squat lobster and asymmetrical intermediaries, whereas we recovered a crab-like ancestor to predate these intermediaries. We acknowledge that our analysis recovers two deep nodes that are unresolved, however symmetrical reconstruction at these nodes seems unlikely (Figure 3A). It must also be noted that the most recent common ancestor of Anomura is unresolved in the Tsang et al. analysis, although it appears to be a crab-like or symmetrical hermit ancestor. The major differences in the two analyses stems from the differences in phylogeny and more specifically the monophyly (our study) or polyphyly  of Paguroidea and families therein (i.e., Pylochelidae). There is agreement with Tsang et al. in the sister group relationship between Paguridae and Lithodoidea, although Tsang et al. used only four lithodid genera (vs. eight in our study) and did not include representatives of Hapalogastridae. In addition, both studies provide strong evidence for the intermediary ancestors directly predating carcinization across Anomura (twice through squat lobster (SI) and once through asymmetrical hermit (ASI), Figure 3A).
The multiple cases of carcinization among the anomurans have been noted since the early 1900s. Borradaile (1916) was the first to propose the term carcinization to explain the crab-like aspects of the hermit crab Porcellanopagurus and the tendency of anomurans to achieve this form, a phenomenon unique to Anomura not evident in other decapod lineages (e.g., lobsters, shrimp). The emergence of the crab-like form is not ‘evenly distributed’ across our phylogeny, first occurring in the older lineages Porcellanidae and Lomisidae and only more recently within the superfamily Lithodoidea. Some questions naturally arise. Why did carcinization occur independently three times during the evolution of the Anomura? Why did the presumably shell-dwelling asymmetrical hermit crab ancestors of lithodid king crabs forsake the use of shells for protection, which already provided them with survival advantage? Morrison et al.  suggest that the crab-like form might represent a key innovation that is associated with an evolutionary advantage, possibly due to the greater mobility and agility provided by this morphology. This seems to be evident within the true crabs, or Brachyura, which dominate decapods in terms of species richness [>6,559 species; 34] and have thrived in marine, freshwater, and terrestrial environments. Although diversification seems to be low in the crab-like anomurans when compared to the brachyurans, fossil evidence and divergence time analyses suggest crab-like anomurans are much younger when compared to the closely related true crabs (Table 2, Figure 4). Furthermore, the crab-like porcellanids are one of the oldest (~172 MYA, MRCA = 139 MYA) and most diverse families of anomurans [~247 species, 22]. Lithodids represent an even younger lineage, originating ~18 MYA, but comprising over 100 extant species. It is plausible that a crab-like form may hold some evolutionary advantage when considering age and diversification within Anomura, although this does not seem to hold true for all groups that underwent the crab-like transition (i.e., monotypic family Lomisidae). A second hypothesis explains the possible advantage of carcinization from a hermit-like ancestor. Previous studies have suggested a free-living body form may have a selective advantage in obtaining food resources when unconstrained by a gastropod shell [54, 73]. An example can be seen in the semi-terrestrial hermit crab, Birgus latro, a species that in the adult stages has lost dependence on shells as protection for the pleon, and instead has developed a calcified body .
The extraordinary morphological and ecological diversity of anomurans has long fascinated evolutionary biologists. Previous studies covering a wide range of faunas have shown how morphological or ecological factors may influence the course of subsequent evolutionary diversification [75–77], and in particular for anomurans it has been hypothesized that the acquisition of the crab-like form may have acted as a key-innovation .
Our analysis reports the pattern of diversification in Anomura to be characterized by a low net rate of diversification, with two major changes in the rates of speciation along its evolutionary history. The initial diversification of the group during the Late Permian was characterized by slow rates of diversification and it was not concomitant with major family radiations, which took place from the Jurassic onwards.
A significant change in the tempo of diversification was identified within the speciose squat-lobster family Chirostylidae, which has a higher speciation rate than the overall tempo of diversification across the anomuran tree (Figure 4). Recent studies based on the munidid squat-lobster genus Paramunida suggest that dramatic environmental change may provide great geological and habitat complexity, which in turn promotes isolation and rapid diversification . The fact that both families, chirostylids and munidids, diverged during the Late Triassic (see Figure 4) and currently occupy deep-sea habitats suggests that similar geological and environmental changes may also have driven major diversification within the Munididae, which shifted habitats at some point because the Jurassic forms are nearly all coral-reef associated. Currently, the family Chirostylidae accounts for 7% of all anomuran species, but the true diversity is underestimated and about 100 new species are in hand of taxonomists . Clearly, a more accurate phylogenetic framework is needed to interpret in detail the exceptionally high speciation rates reported here.
The monotypic family Lomisidae showed a strikingly lower rate than the overall tempo of diversification in Anomura. Lomis hirta is anomalous in its prolonged persistence despite an inferred speciation rate of zero (as recovered by the MEDUSA analysis, see Results). These taxa, old lineages with few extant species, have been reported in several invertebrates and vertebrates [40, 80, 81] and more recently in butterflies , suggesting that extremely low rates of diversification characterize these groups. High extinction rates could also account for this pattern; however, we report that a pure-birth Yule model best explains our data. Under a high-extinction scenario we would expect to see an overabundance of more recently arisen species that simply have not yet gone extinct; such a pattern is not observed in our phylogeny.
Our analysis failed to identify a correlation between the timing of branching events (speciation) and the evolutionary history of carcinized lineages, which suggests that the acquisition of a crab-like form did not play a major role in shaping extant anomuran biodiversity. However, a major limitation of the MEDUSA approach is that rate shifts cannot be assigned below the level of phylogenetic resolution , which prevents us from evaluating if the highly carcinized family Lithodidae underwent an unusual rapid diversification event. The lack of a rate shift in the branch leading to the three collapsed families (Paguridae, Hapalogastridae, and Lithodidae) does not necessarily imply that subclades within that group have not experienced changes in the tempo of diversification, which may be masked by the lack of taxonomic resolution among these taxa. Thus, further studies which focus on clarifying the systematics of the infraorder, with particular emphasis on the families Paguridae, Hapalogastridae, and Lithodidae, are necessary to examine the role of carcinization in anomuran diversification.
Anomuran relationships have been the topic of long debate, likely because of their extraordinary morphological and ecological diversity and their common targeting in fisheries. Here we estimate evolutionary relationships among 19 families, 7 superfamilies, and 137 species of anomurans based on morphological and molecular data to provide the most robust anomuran phylogeny to date. Many families and genera appear to be poly- or paraphyletic suggesting further taxonomic revisions at these levels. Carcinization evolved multiple times during the evolution of Anomura whereas transition into exclusively freshwater or semi-terrestrial environment occurred in the families Aeglidae and Coenobitidae, respectively. Divergence times date the origin of the group in the Late Permian, with subsequent radiations through the Jurassic and Cretaceous. Results suggest that anomurans diversified under low speciation rates with two major changes in the tempo of diversification. First insights suggest that the acquisition of the crab-like form did not play a major role in shaping the extant diversity of Anomura, but further examination is required in order to confirm this pattern.
Our study included extant representatives from 19 families, 77 genera, and 137 species of anomurans. The exceptionally rare family Pylojaquesidae is excluded for lack of molecular grade tissue samples. A total of 345 sequences from 76 of 144 anomuran specimens were new to this study, while sequences for all five genes from 68 taxa were obtained from GenBank. Newly included specimens were collected on cruise and field expeditions, from collaborators, or from the University of Louisiana at Lafayette Zoological (ULLZ) collection of molecular grade specimen and tissue samples (Table 1). Specimens were stored in 80% ethyl alcohol.
The sister group of Anomura is widely accepted to be Brachyura [24–26, 35–38], but because some molecular studies have recovered alternative arrangements [24, 25, 38] we included 18 outgroup taxa (see Table 1) spanning several decapod lineages. Different outgroups were included/excluded to explore sister relationships to Anomura and the impact of outgroup selection on anomuran relationships. They consist of representatives from infraorders Brachyura (5), Axiidea (4), Gebiidea (3), Caridea (4), and suborder Dendrobranchiata (2).
Our morphological data matrix consisted of 156 characters and 154 species (including outgroups) and was constructed in MacClade 4.0 (see Additional files 2 and 3). Citations of previously recognized characters and states are given following characters in Additional file 3. Codings for somatic morphological characters were scored based on examination of sequenced species (Table 1) supplemented by literature. For spermatozoal (130–143) and larval characters (144–156) that are highly conserved (but not available for every sequenced species), reasonable assumptions of monophyly were made in order to optimize the potential contribution of these data. Thus, for these spermatozoal characters, all members of a particular family for which data were available for some members were scored as uniform. For larval characters (primarily first zoeal stage), all members of a particular genus for which data were available for some members were scored as uniform. The larval characters that could be meaningfully scored across the breadth of taxa were included. Others are typically invariant within the family-level (and often superfamily-level) clades, as defined by recent revisionary classifications, and could have been effectively scored at family level. In deference to the possibility that some families might not be monophyletic, however, we took a more conservative, genus-level approach to larval character scoring. Monophyly (or not) of genera, however, with respect to the first zoeal characters employed does not affect results because of the level of generality of characters operating at low taxonomic levels.
Missing data were scored as unknown (?) and polymorphisms were scored as such rather than assuming a plesiomorphic state. Just as alignment gaps in molecular data have been variously treated as a fifth position or as missing in different studies, inapplicable character states in the morphological data may be scored as missing or as an additional character state, ‘inapplicable’ . We scored inapplicable character states as unknown (indicated by ‘-’), rather than an additional state, in order to avoid the possibility of nodes being supported by a non-existent character state .
DNA extraction, PCR, sequencing, and next-generation approaches
Total genomic DNA was extracted from the pleon or gills using the QiagenDNeasy® Blood and Tissue Kit Cat. No. 69582. Two partial mitochondrial genes, 16S and 12S, were amplified by PCR using the following primers, respectively: L2/L9 & 16S1472 or 16SF & 16S1472 [~580 bps, [85–87]] and 12S1F & 12SR or 12SF & 12S1R [~350 bps, ]. The nuclear large subunit 28S rRNA was amplified in sections by 1.3 F/4b, 3.25/4.4b, sA/5b, and 4.8/6b [~2200 bps, [89, 90]]. The nuclear small subunit 18S rRNA was amplified by A/L, C/Y, O/B [~1800 bps, [91, 92]] or by 1 F/2.9, 0.7/bi, 2.0/9R [89, 90], or by shorter internal primers (~1700 bps, B/D18s1R, D18s2F-D18s2R, D18s3F-D18s3R, D18s4F-D18s4R and D18s5F-A ). The histone H3 gene was amplified by H3AF/H3AR [~350 bps, ]. The majority of target gene regions were obtained through traditional Sanger sequencing and data for seven taxa were obtained through next-generation 454 sequencing (see below).
PCR amplifications were performed in 25 μl volumes containing 1 μl of Taq polymerase HotMaster or REDTaq, PCR buffer, 2.5 mM of deoxyribonucleotide triphosphate mix dNTPs, 0.5 μM forward and reverse primer, and extracted DNA. The thermal profile used an initial denaturation for 1 min at 94°C followed by 35–40 cycles of 30 sec at 94°C, 45 sec at 45-60°C depending on gene region, 1 min at 72º and a final extension of 10 min at 72°C. PCR products were purified using plate filters PrepEase™ PCR Purification 96-well Plate Kit, USB Corporation and sequenced with ABI BigDye® terminator mix (Applied Biosystems, Foster City, CA, USA). Cycle sequencing reactions were performed in an Applied Biosystems 9800 Fast Thermal Cycler (Applied Biosystems, Foster City, CA, USA), and sequencing products were run forward and reverse on an ABI 3730xl DNA Analyzer 96-capillary automated sequencer in the Brigham Young University (BYU) sequencing center.
Sequence data for seven taxa were obtained using a novel next-generation sequencing technique TAS: Targeted Amplification Sequencing on the 454 platform [95, 96]. The process required a two-step PCR to prepare selected DNA regions for targeted/directed sequencing. The first PCR used a locus specific primer (e.g., 16S, 12S, etc.) with a 22 bp adapter. These amplicons were cleaned using plate filters PrepEase™ PCR Purification 96-well Plate Kit, USB Corporation. One μl of cleaned PCR product was used as template for the second PCR. PCR II incorporated a 10 bp barcode multiplex identifier, MID, 4 bp key, and a 21 bp 454 Titanium primer. Samples were again cleaned using the Millipore system and subsequently combined in emulsion PCR and sequenced via 454 GS FLX Titanium pyrosequencing technology (Roche) at the BYU sequencing center. The bioinformatic pipeline, BarcodeCruncher, allowed us to exclude short reads, trim adapters, identify contamination, parse barcoded sequences, and assembly consensus sequences for phylogenetic reconstruction [for full description of methods see [95, 96]].
Sequences were cleaned and assembled using Sequencher 4.9 (GeneCodes, Ann Arbor, MI, USA). To check for pseudogenes, we followed suggestions by Song et al. (2008), which included extracting DNA from tissue with high amounts of mitochondrial gill tissue, translating protein-coding sequences H3 to check for indels and stop codons, comparing sequences among closely-related species, and building individual gene trees to ensure similar topologies . Comparing gene trees and BLAST searches helped identify contamination. Two datasets were assembled: 1) molecular dataset including all 5 gene regions 2) combined dataset including molecular + morphological data.
Individual gene alignments were performed using MAFFT, implementing the “E-INS-i” option. For non-protein coding genes 12S, 16S, 18S, 28S, GBlocks v0.91b were used to exclude regions of the alignment with questionable positional homology . The parameters used in GBlocks for 12S, 16S, 18S, 28S, were as follow: minimum number of sequences for a conserved position = 50/77/77/79; minimum number of sequences for a flanking position = 50/77/80/79; maximum number of contiguous non-conserved positions = 8/8/8/8; minimum length of a block = 5/5/5/5; allowed gap positions = half/half/half/half. Final alignments included 300, 474, 1632, and 931 base pairs for 12S, 16S, 18S, and 28S, respectively. After trimming for primer residue, the H3 alignment resulted in 332 base pairs. In MESQUITE , all genes were concatenated 3669 basepairs and partitioned for analysis. The final molecular dataset included 162 individuals as 3669 basepairs (5 genes) while the combined data set included the molecular dataset plus an additional 156 morphological characters.
The Maximum Likelihood (ML) analysis was conducted using RAxML Randomized Axelerated Maximum Likelihood [100–102]. Likelihood settings followed the General Time Reversible Model GTR with a gamma distribution and RAxML estimated all free parameters following a partitioned dataset. The first algorithm used in the analysis was the “-f a” option, for a rapid bootstrap analysis and search for the best tree in a single pass. The second algorithm implemented another search for the best tree implementing a “-f d” option of -#200 iterations of random starting trees. Likelihoods were compared to determine the best tree and bootstraps were mapped on the resulting topology. Confidence in the resulting topology was assessed using non-parametric bootstrap estimates  with 1000 pseudoreplicates and values > 50% are presented on the resulting phylogeny.
Bayesian analyses (BA) were performed in MrBayes v3.1.2b4  for the molecular and combined datasets morphology + molecular. We used the Markov k Mk,  model for the morphological characters equal state frequencies, combined with gamma distributes rates across sites. The model of evolution that best fit the individual datasets was determined by MODELTEST 3.7  and these parameters were applied to our molecular dataset. Three independent BA analyses were implemented each with 20 chains and a starting tree obtained from the ML analysis to help reach convergence. The molecular analysis ran for 30,000,000 generations, sampling one tree every 1000 generations. The combined analysis ran for 50,000,000 generations, sampling one tree every 5000 generations. To ensure that independent analyses converged on similar values, we graphically compared all likelihood parameters and scored means and variances using the program Tracer v1.4 . Burn-in and stationary distributions were determined by observing the likelihood -LnL scores and split frequencies for the data (~10 million generations). A 50% majority-rule consensus tree was obtained from the remaining saved trees, once the data reached convergence. Posterior probabilities Pp for clades were compared for congruence and post-burn-in-trees were combined between individual runs. Values > 0.5 are presented on the BA phylogram presented as percentages. All analyses were run on Marylou6 Dell PowerEdge M610 computing cluster at Brigham Young University. High support is defined as ≥ 95/70 Pp/bs, marginal support is ≥ 85/65 and low support is ≤ 84/64.
Alternative hypothesis testing
A partitioned S-H test  was used to test whether previous hypotheses of anomuran evolution implicit in modern, morphologically-based classifications [14–16] and morphological and/or molecular phylogenies [19, 20, 23] were significantly worse than our best ML tree. The test was implemented in RAxML using the same data partitions used to estimate our phylogeny. As in the ML analysis, the GTRGAMMA model was applied to each partition. Seven independent constrained tree topologies were constructed in Mesquite v.2.71 . Topological constraints were forced to the following monophyletic clades: Diogenidae, Hapalogastridae, Lithodidae, Munididae, Paguridae, Paguroidea, and Pylochelidae. These clades were tested to examine the validity of current generic assignments by testing the poly- and paraphyly of the families and superfamilies in the tree. Lastly, to test the king to hermit hypothesis, we forced the topologies: 1) king crabs (Lithodidae) ancestral to hermit crab superfamily Paguroidea and 2) king crabs (Lithodidae) ancestral to hermit crab family Paguridae to test if these hypotheses are significantly worse than a “hermit to king” evolutionary pathway.
We used ancestral state reconstruction (ASR) methods implemented in Mesquite v.2.71  to examine character evolution across the anomurans. We traced evolutionary pathways of two characters: body form and habitat. Body form was assigned as follows: 0: crab-like (carcinized), 1: squat lobster form, 2: asymmetrical pleon, 3: symmetrical pleon. Habitat was assigned as 0: semi-terrestrial, 1: freshwater, 2: marine. These characters were optimized across our best estimate of anomuran relationships (=combined (molecular + morphology) Bayesian phylogeny). Because the importance of employing different methods for ASR has been documented, we used both maximum parsimony and maximum likelihood methods . Likelihood methods are often preferred over parsimony reconstructions since they take into account branch lengths, all character state possibilities, and model evolution . The model of evolution used in the maximum likelihood analysis was the Markov k-state 1 (Mk1) parameter model, which allows equal probability for any character state change. All characters were scored and compiled based on specimen observation and/or literature searches. Reconstructions among all outgroup taxa are not shown.
Divergence time analyses
To estimate the relative timing of origins, diversification, body form evolution and transition in habitat, Bayesian molecular dating methods were implemented in BEAST v1.5.2 (Bayesian evolutionary analysis by sampling trees) . BEAST allows for missing data, multiple calibration points, relaxed clock models, and increased flexibility of model parameters when compared to other dating methods (e.g., Multidivtime). Substitution and clock models were unlinked and the dataset was partitioned by gene following models of evolution generated by MODELTEST v3.7. A relaxed uncorrelated lognormal clock model and Yule speciation tree prior were selected. We recognize that there are varied models to consider when using relaxed dating methods. Simulation studies that have compared accuracy of relaxed clock methods have recovered conflicting results, with some favoring uncorrelated models , others favoring autocorrelated models  and some favoring both (autocorrelated and uncorrelated) depending on the dataset . We chose an uncorrelated relaxed clock method due to the biological data under investigation and evidence from our divergence time analysis. It has been suggested that autocorrelation in life history traits (one biological assumption underlying autocorrelated relaxed clocks) would be less relevant in studies focused at high taxonomic levels, divergent taxa [112, 115], and sparely sampled datasets . We had little reason to believe we had autocorrelation in our anomuran dataset considering we were reconstructing a phylogeny across an infraorder that originated over 250 MYA (oldest fossil evidence = 201–228 MYA). Moreover, our sampling method focused at the superfamily and family level (and not within genera or species). Secondly, it has been suggested that rate autocorrelation can be measured by comparing the posterior and prior distributions of covariance in rates on neighboring branches [112, 116]. All covariance estimates in our divergence time analyses suggested we had no strong evidence for autocorrelation of rates in our phylogeny. The statistic measure between parent and child branches contained values that span zero suggesting branches with slow and fast rates are adjacent on the tree. For abovementioned reasons, we did not assume autocorrelated rates across our tree. Our best estimate of phylogeny (=combined (morphology + molecular) tree) was used as a starting tree and the tree searching parameters were removed from the BEAUTI xml file.
Using the non-parametric rate smoothing algorithm in r8s , we made branch lengths proportional to the timing chronogram rather than the substitutions per site phylogram. This allowed the tree to adhere to the topological and temporal constraints of using fossil calibrations. We implemented two runs in BEAST with MCMC chain length of 350 million generations logging parameters every 3,500 samples. To ensure that analyses converged on similar values with acceptable mixing, likelihood stationary and burn-in values, we used the program Tracer v1.4 . The runs were combined using LogCombiner . Estimates of the mean divergence times with 95% highest posterior density regions HPD and posterior probabilities represented as percentages are noted on the chronogram. The BEAST analyses were performed on the Marylou6 Dell PowerEdge M610 computing cluster at Brigham Young University.
Fossil and time calibrations
A total of 31 fossils were included in the analysis. We included fossils that represented the oldest known specimen for a particular family and/or genus (Table 2). Only fossils that could be confidently assigned to clades based on personal observation and/or previous literature were included in the analysis [, Table 2]. We followed recommendations by Parham et al.  when justifying fossil placement. Both deep and shallow fossil representatives were included. All fossils were placed at the crown (i.e., most recent common ancestor MRCA) or stem (i.e., node directly preceding the MRCA node of the clade). There was disagreement among authors about the familial assignment of Juracrista, as either a munidid or galatheid. The munidids and galatheids, however, are closely related so the affect on calibration is minimal. Therefore, we have followed the original taxonomy  and retained Juracrista in Munididae. Because divergence time should predate the fossil occurrence, all calibrations followed an exponential prior with the offset value set to the minimum calibration age. This distribution is suitable for modeling fossil calibrations, because it allows us to avoid a hard upper bound while providing an increased probability with the age range of fossil discovery .
We analysed patterns of diversification along the anomuran tree using a reworked version of MEDUSA http://www.webpages.uidaho.edu/~lukeh/software/index.html. This method infers clade-specific changes in the tempo of diversification (rate acceleration and rate reduction) across a tree from phylogenetic branch lengths and taxonomic extant species richness information, the latter to account for incomplete taxon sampling . MEDUSA utilized a stepwise AIC approach to adding clade-wise diversification models Yule or birth-death to a tree until the decrease in AIC failed to exceed a set threshold, which was dependent on tree size. Diversification analysis implemented in MEDUSA required that in the analysed tree terminal tips represented monophyletic taxonomic groups with known species richness. We attempted to resolve clades to the taxonomic level of family; however our best phylogenetic hypothesis did not support the monophyly of the families Munididae, Diogenidae, and Paguridae (Figure 2). This non-monophyly was also observed in molecular trees, (Additional file 4). Therefore, assuming the current family classification would lead us to infer shifts in diversification rates inconsistent with the phylogenetic relationships found in this study. Recent studies have suggested that further subdivisions are conceivable within Munididae , with the recognition of a well-differentiated clade including the genera Munida, Cervimunida, and Pleuroncodes and a second clade including the rest of the genera [19, 122, 123]. Our study highly supports the existence of the Munida, Cervimunida, and Pleuroncodes clade (Pp = 100), yielding also moderate support for the inclusion of the genus Sadayoshia (Pp = 91). Hence, species richness within the family Munididae was assigned to two different monophyletic subclades (see Additional file 4). In agreement with our data, Diogenidae has been recently recovered as paraphyletic . Previous studies have showed the genera Paguristes or Areopaguristes to be separated from other diogenids [19, 33, 95, 124], and closely related to Pseudopaguristes and Tetralobistes. Thus, in order to assign a known species richness value, we pruned this clade to include all the species belonging to those 4 genera, and recognized a second clade to accommodate the rest of species of Diogenidae plus Coenobitidae (See Additional file 4). Finally, the family Paguridae is recognized as one of the most challenging groups within Anomura and its monophyletic origin (in relation to Lithodidae and Hapalogastridae) has been debated [19–21]. Given the high species richness of this family (~ 542 species) and the lack of an appropriate taxonomic framework, it is not possible to account for the phylogenetic position of each unsampled extant species or all possible lineages. Thus, subdividing this clade into smaller subclades to assign species richness would be arbitrary, potentially leading to spurious results. Although suboptimal, the family Paguridae was collapsed together with the families Hapalaogastridae and Lithodidae in order to satisfy the monophyly assumption of MEDUSA.
We obtained information about species richness for major lineages of Anomura using the most recent published checklists [22, 44, 46, 126–130]. This information was assigned to 18 lineages of our phylogeny after pruning terminals belonging to the same monophyletic groups according to the aforementioned criteria (see Additional file 4).
Most recent common ancestor
Million years ago
We would like to thank a wide variety of friends and colleagues who have helped us collect anomurans from around the world. Staff and scientists at the Department of Invertebrate Zoology, National Museum of Natural History, Smithsonian Institution; Universidad Nacional Autónoma de México; Smithsonian Tropical Research Institute’ Marine Station, Bocas del Toro, and R/V Urraca; and Louisiana Universities Marine Consortium LUMCON, R/V Pelican, kindly assisted with research cruises, shipping of loans, or hosting our visits. We are deeply grateful to Joseph Brown and Matthew Pennell for their advice to perform diversification analyses and useful comments. We are thankful to T.Y Chan, Museum of Victoria, Melbourne (photo by M. Marmach), DLF, RL, S. DeJouannet, K. Schnabel (NIWA) whom provided anomuran photographs. Our work was partially supported by Brigham Young University and a Tree of Life grant from the US NSF EF-0531762 awarded to KAC, DLF, CES, RMF, Joel Martin, and Nikki Hannegan. This work was also partially funded by a grant from the U.S. National Science Foundation (NSF) PIRE for research on Patagonian Biodiversity to support collaboration among the following institutions (listed alphabetically): Brigham Young University, Centro Nacional Patagónico (AR), Dalhousie University, Instituto Bota´nico Darwinion (AR), Universidad Austral de Chile, Universidad Nacional del Comahue, Universidad de Concepción, and George Washington University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. PC was supported by a grant from the Fundación Caja Madrid. This is contribution no. 160 for the UL Laboratory for Crustacean Research and no. 915 for the Smithsonian Marine Station,-Ft. Pierce.
- Chablais J, Feldmann RM, Schweitzer CE: A new Triassic decapod, Platykotta akaina, from the Arabian shelf of the northern United Arab Emirates: earliest occurrence of the Anomura. Palaeontologische Zeitschrift. 2011, 85 (1): 93-102. 10.1007/s12542-010-0080-y.Google Scholar
- Macpherson E, Segonzac M: Species of the genus Munidopsis (Crustacea, Decapoda, Galatheidae) from the deep Atlantic Ocean, including cold-seep and hydrothermal vent areas. Zootaxa. 2005, 1095: 3-60.Google Scholar
- Dawson EW: King crabs of the world (Crustacea: Lithodidae) and their fisheries: a comprehensive bibliography. Misc Publ. 1989, 101: 338-Google Scholar
- Lovrich GA: La pesqueria mixta de las centollas Lithodes santolla y Paralomis gradulosa (Anomura: Lithodidae) en Tierra del Fuego, Argentina. Investig Mar. 1997, 25: 41-57.Google Scholar
- Wehrtmann IS, Acuna E: Squat Lobsters Fisheries. The Biology of Squat Lobsters. Edited by: Poore GCB, Ahyong S, Taylor J. 2011, Museum Victoria: CSIRO Publishing, 297-322.Google Scholar
- Macpherson E, Jones W, Segonzac M: A new squat lobster family of Galatheoidea (Crustacea, Decapoda: Anomura) from the hydrothermal vents of the Pacific-Antarctic Ridge. Zoosystema. 2005, 27: 709-723.Google Scholar
- Lovrich GA, Vinuesa JH: Reproductive potential of the lithodids Lithodes santolla and Paralomis granulosa (Anomura, Decapoda) in the Beagle Channel, Argentina. Sci Mar. 1999, 63: 355-360.Google Scholar
- Perez-Losada M, Jara CG, Bond-Buckup G, Crandall KA: Conservation phylogenetics of Chilean freshwater crabs Aegla (Anomura, Aeglidae): assigning priorities for aquatic habitat protection. Biol Conserv. 2002, 105 (3): 345-353. 10.1016/S0006-3207(01)00218-X.Google Scholar
- Jara C: Taxonomía, sistemática y zoogeográfica de las especies chilenas del género Aegla Leach (Crustacea: Decapoda: Anomura: Aeglidae). 1996, Concepción: University of Concepción, ChileGoogle Scholar
- Xu JW, Perez-Losada M, Jara CG, Crandall KA: Pleistocene glaciation leaves deep signature on the freshwater crab Aegla alacalufi in Chilean Patagonia. Mol Ecol. 2009, 18 (5): 904-918. 10.1111/j.1365-294X.2008.04070.x.PubMedGoogle Scholar
- Lemaitre R, McLaughlin PA: Recent advances and conflicts in concepts of anomuran phylogeny. Arthropod Systematics & Phylogeny. 2009, 67 (2): 119-135.Google Scholar
- Mclaughlin PA, Holthuis LB: Anomura Versus Anomala. Crustaceana. 1985, 49 (Sep): 204-209.Google Scholar
- Ahyong ST, Schnabel KE, Macpherson E: Phylogeny and fossil record of marine squat lobsters. The Biology of Squat Lobsters. Edited by: Poore GCB, Ahyong S, Taylor J. 2011, Museum Victoria: CSIRO Publishing, 73-104.Google Scholar
- Forest J: Ethology and distribution of pylochelidae (crustacea, decapoda coenobitoidea). Bulletin of Marine Science. 1987, 41 (2): 309-321.Google Scholar
- McLaughlin PA: A review of the phylogenetic position of the Lomidae (Crustacea, Decapoda, Anomala). J Crustacean Biol. 1983, 3 (3): 431-437.Google Scholar
- McLaughlin PA: Illustrated keys to families and genera of the superfamily Paguroidea (Crustacea: Decapoda: Anomura), with diagnoses of genera Paguridae. Memoirs of Museum Victoria. 2003, 60 (1): 111-114.Google Scholar
- Schnabel KE, Ahyong ST: A new classification of the Chirostyloidea (Crustacea: Decapoda: Anomura). Zootaxa. 2010, 2687: 56-64.Google Scholar
- Ahyong ST, Baba K, Macpherson E, Poore GCB: A new classification of the Galatheoidea (Crustacea: Decapoda: Anomura). Zootaxa. 2010, 2676: 57-68.Google Scholar
- Schnabel KE, Ahyong ST, Maas EW: Galatheoidea are not monophyletic - Molecular and morphological phylogeny of the squat lobsters (Decapoda: Anomura) with recognition of a new superfamily. Mol Phylogenet Evol. 2011, 58 (2): 157-168. 10.1016/j.ympev.2010.11.011.PubMedGoogle Scholar
- Tsang LM, Chan TY, Ahyong ST, Chu KH: Hermit to King, or Hermit to all: multiple transitions to Crab-like forms from Hermit crab ancestors. Syst Biol. 2011, 60 (5): 616-629. 10.1093/sysbio/syr063.PubMedGoogle Scholar
- Ahyong ST, Schnabel KE, Maas EW: Anomura phylogeny: New insights from molecular data. Decapod Crustacean Phylogenetics (Crustacean Issues 18). Edited by: Martin JW, Felder DL, Crandall KA. 2009, Bocan Raton, FL: CRC Press, 391-408.Google Scholar
- De Grave S, Pentcheff ND, Ahyong ST, Chan T-Y, Crandall KA, Dworschak PC, Felder DL, Feldman RM, Fransen C, Goulding L, et al: A classification of living and fossil genera of decapod crustaceans. Raffles B Zool. 2009, 21: 1-109.Google Scholar
- McLaughlin PA, Lemaitre R, Sorhannus U: Hermit crab phylogeny: a reappraisal and its “fall-out”. J Crustacean Biol. 2007, 27 (1): 97-115. 10.1651/S-2675.1.Google Scholar
- Porter ML, Perez-Losada M, Crandall KA: Model-based multi-locus estimation of decapod phylogeny and divergence times. Mol Phylogenet Evol. 2005, 37 (2): 355-369. 10.1016/j.ympev.2005.06.021.PubMedGoogle Scholar
- Bracken H, Toon A, Felder DL, Martin JW, Finley M, Rasmussen J, Palero F, Crandall KA: The Decapod Tree of Life: compiling the data and moving toward a consensus of decapod evolution. Arthropod Systematics and Phylogenetics. 2009, 67: 99-116.Google Scholar
- Tsang LM, Ma KY, Ahyong ST, Chan TY, Chu KH: Phylogeny of Decapoda using two nuclear protein-coding genes: Origin and evolution of the Reptantia. Mol Phylogenet Evol. 2008, 48 (1): 359-368. 10.1016/j.ympev.2008.04.009.PubMedGoogle Scholar
- McLaughlin PA, Lemaitre R: A New Classification for the Pylochelidae (Decapoda: Anomura: Paguroidea) and Descriptions of New Taxa. Raffles B Zool. 2009, Suppl. 20: 159-231.Google Scholar
- Milne Edwards H: Recherches sur l’organisation et la classification naturelle des Crustacéa Décapodes. Annales des Sciences Naturelles Zoologie, Paris. 1832, 25: 298-332.Google Scholar
- Milne Edwards H: Histoire naturelle des Crustacés, comprenant l’ana¬tomie, la physiologie et la classification de ces animaux. Librairie Encyclopédique de Roret, Paris. 1837, 2: 1-532. atlas, pp. 531–532, pls. 531–542Google Scholar
- Boas JEV: Studier over Decapodernes Slsegtskabsforhold. Slaegtskabsforhold Kongelige Danske Videnskabernes Selskabs Skrifter, Naturvidenskabelig of Mathema¬tisk Afdeling. 1880, 1 (2): 26-210.Google Scholar
- Boas JEV: Wissenschaftliche Mittheilungen. 1. Lithodes und Pagurus. Zool Anz. 1880, 3: 349-352.Google Scholar
- McLaughlin PA, Lemaitre R: Carcinization in the Anomura-fact or fiction? I. Evidence from adult morphology. Contrib Zool. 1997, 67 (2): 79-123.Google Scholar
- Morrison CL, Harvey AW, Lavery S, Tieu K, Huang Y, Cunningham CW: Mitochondrial gene rearrangements confirm the parallel evolution of the crab-like form. Proc R Soc London. 2002, 269: 345-350. 10.1098/rspb.2001.1886.Google Scholar
- Ng PKL, Guinot D, Davie PJ: Systema Brachyurorum: Part I. An annotated checklist of extant brachyuran crabs of the world. Raffles B Zool. 2008, 17: 1-286.Google Scholar
- Dixon CJ, Ahyong S, Schram FR: A new hypothesis of decapod phylogeny. Crustaceana. 2003, 76 (8): 935-975. 10.1163/156854003771997846.Google Scholar
- Schram FR, Dixon CJ: Decapod Phylogeny: addition of fossil evidence to a robust morphological cladistic data set. Bulletin of the Mitzunami Fossil Museum. 2004, 31: 1-19.Google Scholar
- Ahyong S, O’Meally D: Phylogeny of the Decapoda Reptantia: Resolution using three molecular loci and morphology. Raffles Bulletin Zoology. 2004, 52 (2): 673-693.Google Scholar
- Toon A, Finley M, Staples J, Crandall KA: Decapod phylogenetics and molecular evolution. Decapod Crustacean Phylogenetics (Crustacean Issues 18). Edited by: Martin JW, Crandall KA, Felder DL. 2009, Bocan Raton, FL: CRC Press, 14-28.Google Scholar
- McLaughlin PA: Hermit crabs - Are they really polyphyletic. J Crustacean Biol. 1983, 3 (4): 608-621.Google Scholar
- Alfaro ME, Santini F, Brock C, Alamillo H, Dornburg A, Rabosky DL, Carnevale G, Harmon LJ: Nine exceptional radiations plus high turnover explain species diversity in jawed vertebrates. Proc Natl Acad Sci. 2009, 106 (32): 13410-13414. 10.1073/pnas.0811087106.PubMed CentralPubMedGoogle Scholar
- Reimann A, Richter S, Scholtz G: Phylogeny of the Anomala (Crustacea, Decapoda, Reptantia) based on the ossicles of the foregut. Zool Anz. 2011, 250 (4): 316-342. 10.1016/j.jcz.2011.05.006.Google Scholar
- Lemaitre R, McLaughlin PA, Sorhannus U: Phylogenetic relationships within the Pylochelidae (Decapoda: Anomura: Paguroidea): A cladistic analysis based on morphological characters. Zootaxa. 2022, 2009: 1-14.Google Scholar
- McLaughlin PA, Lemaitre R, Tudge CC: Carcinization in the Anomura - fact or fiction? II. Evidence from larval, megalopal and early juvenile morphology. Contrib Zool. 2004, 73 (3): 165-205.Google Scholar
- Baba K, Macpherson E, Poore GCB, Ahyong ST, Bermudez A, Cabezas P, Lin CW, Nizinski M, Rodrigues C, Schnabel KE: Catalogue of squat lobsters of the world (Crustacea: Decapoda: Anomura - families Chirostylidae, Galatheidae and Kiwaidae). Zootaxa. 1905, 2008: 1-220.Google Scholar
- Martin JW, Davis GE: An updated classification of the Recent Crustacea. Natural History Museum of Los Angeles County, Science Series. 2001, 39: 1-124.Google Scholar
- McLaughlin PA, Boyko CB, Crandall KA, Komai T, Lemaitre R, Osawa M, Rahayu DL: Annotated checklist of anomuran decapod crustaceans of the world (exclusive of the kiwaoidea and families chirostylidae and galatheidae of the galatheoidea) - preamble and scope. Raffles B Zool. 2010, 1-4.Google Scholar
- Boas JEV: Die verwandtschaftliche Stellung der Gattung Lithodes. K Danske Vidensk Selsk Skr. 1924, 4: 1-34.Google Scholar
- Richter S, Scholtz G: Morphological evidence for a hermit-crab ancestry of lithodids (crustacea, decapoda, anomala, paguroidea). Zool Anz. 1994, 233 (5–6): 187-210.Google Scholar
- Ahyong ST, Andreakis N, Taylor J: Mitochondrial phylogeny of the deep-sea squat lobsters, Munidopsidae (Galatheoidea). Zool Anz. 2011, 250 (4): 367-377. 10.1016/j.jcz.2011.06.005.Google Scholar
- Boas FEV: Studier over decapodernes slaegtskabsforhold. Slaegtskabsforhold Kongelige Danske Videnskabernes Selskabs Skrifter, Naturvidenskabelig of Mathema¬tisk Afdeling. 1880, 5 (6): 25-210.Google Scholar
- Bouvier EL: Sur la transformation des Paguriens en crabes anomoures de la sousfamille des Lithodines. Comp Rendu heb Sean Acad Sciences, Paris. 1894, 119: 350-352.Google Scholar
- Bouvier EL: La transformation des Bernards l’Ermite en Lithodes. Naturaliste. 1897, 2 (19): 41-43.Google Scholar
- Bouvier EL: Sur les Lithodin^s hapalogastriques HapalogaslerelDermaturus. CR Soc Philom. 1985, 18: 56-Google Scholar
- Cunningham CW, Blackstone NW, Buss LW: Evolution of king crabs from hermit-crab ancestors. Nature. 1992, 355 (6360): 539-542. 10.1038/355539a0.PubMedGoogle Scholar
- Gould SJ: We are all monkeys’ uncles. Natural History. 1992, 6 (92): 14-21.Google Scholar
- Scholtz G, Richter S: Phylogenetic systematics of the reptantian decapoda (Crustacea, Malacostraca). Zoological J Linnean Soc. 1995, 113 (3): 289-328.Google Scholar
- Keiler J, Richter S, Wirkner CS: Evolutionary morphology of the hemolymph vascular system in hermit and king crabs (crustacea: decapoda: anomala). J Morphology. 2013, In pressGoogle Scholar
- McLaughlin PA, Lemaitre R: Aspects of Evolution in the anomuran superfamily Paguroidea: one larval prospective. Invertebrate Reproduction and Development. 2000, 38 (3): 159-169. 10.1080/07924259.2000.9652453.Google Scholar
- Keiler J, Richter S: Morphological diversity of setae on the grooming legs in Anomala (Decapoda: Reptantia) revealed by scanning electron microscopy. Zool Anz. 2011, 250 (4): 343-366. 10.1016/j.jcz.2011.04.004.Google Scholar
- Stevens BA: Hermit crabs of Friday Harbor. Publications of the Puget Sound Biological Station. 1925, 3: 273-308.Google Scholar
- Komai T: A new species of the genus Discorsopagurus (Crustacea: Deca¬poda: Paguridae) from Japan, previously known as D. schmitti (Stevens). Proc Biol Soc Washington. 1995, 108 (4): 617-628.Google Scholar
- Komai T: A new species of the hermit crab genus Discorsopagurus McLaughlin (Crustacea: Decapoda: Anomura: Paguridae) from Japan. Natural History Res. 2003, 7 (2): 181-192.Google Scholar
- Komai T, Takeda M: A new hermit crab of the genus Discorsopagurus (Crustacea: Anomura: Paguridae) from Japan. Species Diversity. 1996, 1: 75-85.Google Scholar
- Martin J, Ormsby B: A large brachyuran-like larva of the Hippidae (Crustacea: Decapoda: Anomura) from the Banda Sea, Indonesia: the largest known zoea. P Biol Soc Wash. 1991, 104 (3): 561-568.Google Scholar
- Maurer F, Rettori R, Martini R: Triassic stratigraphy, facies and evolution of the Arabian shelf in the northern United Arab Emirates. Int J Earth Sci. 2008, 97 (4): 765-784. 10.1007/s00531-007-0194-y.Google Scholar
- Perez-Losada M, Bond-Buckup G, Jara CG, Crandall KA: Molecular systematics and biogeography of the southern South American freshwater “crabs” Aegla (Decapoda : Anomura : Aeglidae) using multiple heuristic tree search approaches. Syst Biol. 2004, 53 (5): 767-780. 10.1080/10635150490522331.PubMedGoogle Scholar
- Mclaughlin PA, Murray T: Clibanarius-Fonticola, New Species (Anomura, Paguridae, Diogenidae), from a Fresh-Water Pool on Espiritu-Santo, Vanuatu. J Crustacean Biol. 1990, 10 (4): 695-702. 10.2307/1548413.Google Scholar
- Feldmann RM: Haumuriaegla glaessneri n. gen. and sp. (Decapoda, Anomura, Aeglidae) from Haumurian (Late Cretaceous) Rocks near Cheviot, New Zealand. New Zeal J Geol Geop. 1984, 27 (3): 379-385.Google Scholar
- Henderson JR: Report on the Anomura collected by H.M.S. Challenger during the years 1873–1876. In: Report of the Scientific Results of the Voyage of H.M.S. Challenger during the years 1873–76. Zoology. 1888, 27: 1-221. 221 plGoogle Scholar
- Ortman AE: Die Decapoden-Krebse des Strassburger Museums, mit besonderer Berücksichtigung der von Hern Dr. Doderlein bei Japan und den Liu-Kiu-Inseln gesammelten und z. Z. im Strassburger Museum aufbewahrten Formen IV. Theil. Die Abtheilungen Galatheidea and Paguridea. Zoologische Jahrbücher Abtheilung für Systematik. 1892, Jena: Geographie und Biologie der Thiere 6 Gustav Fisher, 241-326.Google Scholar
- MacDonald JD, Pike RB, Williamson DI: Larvae of the British species of Diogenes, Pagurus, Anapagurus and Lithodes (Crustaca, Decapoda). Proc R Soc London. 1957, 128: 209-257.Google Scholar
- van Bakel BWM, Fraaije RHB, Jagt JWM, Artal P: An unexpected diversity of Late Jurassic hermit crabs (Crustacea, Decapoda, Anomura) in Central Europe. Neues Jahrbuch fur Geologie und Palaontologie- Abhandlungen. 2008, 250 (2): 137-156. 10.1127/0077-7749/2008/0250-0137.Google Scholar
- Blackstone NW: Size, shell-living, and carcinization in geographic populations of a hermit crab, Pagurus hirsutiusculus. J Zool. 1989, 217: 477-490. 10.1111/j.1469-7998.1989.tb02503.x.Google Scholar
- Reese ES: Shell use: an adaptation for emigration from the sea by coconut crabs. Science. 1968, 161: 385-386. 10.1126/science.161.3839.385.PubMedGoogle Scholar
- Alfaro ME, Santini F, Brock CD: Do reefs drive diversification in marine teleosts? Evidence from the pufferfish and their allies (Order Tetraodontiformes). Evol. 2007, 61 (9): 2104-2126. 10.1111/j.1558-5646.2007.00182.x.Google Scholar
- Egan AN, Crandall KA: Divergence and diversification in North American Psoraleeae (Fabaceae) due to climate change. BMC Biol. 2008, 6 (1): 55-10.1186/1741-7007-6-55.PubMed CentralPubMedGoogle Scholar
- Rabosky DL, Donnellan SC, Talaba AL, Lovette IJ: Exceptional among-lineage variation in diversification rates during the radiation of Australia’s most diverse vertebrate clade. Proc R Soc B: Biol Sci. 2007, 274 (1628): 2915-2923. 10.1098/rspb.2007.0924.Google Scholar
- Cabezas P, Sanmartin I, Paulay G, Macpherson E, Machordom A: Deep under the sea: unraveling the evolutionary history of the deep-sea squat lobster paramunida (decapoda, munididae). Evol. 2012, 66 (6): 1878-1896. 10.1111/j.1558-5646.2011.01560.x.Google Scholar
- Macpherson E, Baba K: Taxonomy of squat lobsters. The Biology of Squat Lobsters. Edited by: Poore GCB, Ahyong S, Taylor J. 2011, Museum Victoria: CSIRO PublishingGoogle Scholar
- Avise JC, Nelson WS, Sugita H: A speciational history of “living fossils”: molecular evolutionary patterns in horseshoe crabs. Evol. 1994, 1986-2001.Google Scholar
- Simpson GG: Tempo and Mode in Evolution (Columbia Classics in Evolution Series). 1944, New York: Columbia University PressGoogle Scholar
- Heikkilä M, Kaila L, Mutanen M, Peña C, Wahlberg N: Cretaceous origin and repeated tertiary diversification of the redefined butterflies. Proc R Soc B: Biol Sci. 2012, 279 (1731): 1093-1099. 10.1098/rspb.2011.1430.Google Scholar
- Platnick NI, Griswold CE, Coddington JA: On missing entries in cladistic-analysis. Cladistics Int J Willi Hennig Soc. 1991, 7 (4): 337-343. 10.1111/j.1096-0031.1991.tb00042.x.Google Scholar
- Reid AL: Review of the Peripatopsidae (Onychophora) in Australia, with comments on peripatopsid relationships. Invertebr Taxon. 1996, 10 (4): 663-936. 10.1071/IT9960663.Google Scholar
- Palumbi S, Martin A, Romano S, McMillan WO, Stice L, Grabowski G: The Simple Fool’s Guide to PCR. 1991, Honolulu: Department of Zoology and Kewalo Marine Laboratory, University of HawaiiGoogle Scholar
- Crandall KA, Fitzpatrick JF: Crayfish molecular systematics: Using a combination of procedures to estimate phylogeny. Syst Biol. 1996, 45: 1-26. 10.1093/sysbio/45.1.1.Google Scholar
- Schubart CD, Neigel JE, Felder DL: Use of the mitochondrial16S rRNA gene for phylogenetic and population studies of Crustacea. Crustacean Issues. 2000, 12: 817-830.Google Scholar
- Buhay JE, Moni G, Mann N, Crandall KA: Molecular taxonomy in the dark: Evolutionary history, phylogeography, and diversity of cave crayfish in the subgenus Aviticambarus, genus Cambarus. Mol Phylogenet Evol. 2007, 42 (2): 435-448. 10.1016/j.ympev.2006.07.014.PubMedGoogle Scholar
- Whiting MF: Mecoptera is paraphyletic: multiple genes and phylogeny of Mecoptera and Siphonaptera. Zoologica Scripta. 2002, 31 (1): 93-104. 10.1046/j.0300-3256.2001.00095.x.Google Scholar
- Whiting MF, Carpenter JC, Wheeler QD, Wheeler WC: The strepsiptera problem: Phylogeny of the holometabolous insect orders inferred from 18S and 28S ribosomal DNA sequences and morphology. Syst Biol. 1997, 46 (1): 1-68.PubMedGoogle Scholar
- Apakupakul K, Siddall ME, Burreson EM: Higher level relationships of leeches (Annelida : Clitellata : Euhirudinea) based on morphology and gene sequences. Mol Phylogenet Evol. 1999, 12 (3): 350-359. 10.1006/mpev.1999.0639.PubMedGoogle Scholar
- Medlin LK, Elwood HJ, Stickel S, Sogin ML: The characterization of enzymatically amplified eukaryotic Ids-like rRNA coding regions. Gene. 1988, 71: 491-499. 10.1016/0378-1119(88)90066-2.PubMedGoogle Scholar
- Bracken HD, De Grave S, Felder DL: Phylogeny of the infraorder Caridea based on nuclear and mitochondrial genes (Crustacea: Decapoda). Decapod Crustacean Phylogenetics (Crustacean Issues 18). Edited by: Martin JW, Crandall KA, Felder DL. 2009, Bocan Raton: CRC Press, 274-300.Google Scholar
- Colgan DJ, McLauchlan A, Wilson GDF, Livingston SP, Edgecombe GD, Macaranas J, Cassis G, Gray MR: Histone 3 and U2 snRNA DNA sequences and arthropod molecular evolution. Aust J Zool. 1998, 46: 419-437. 10.1071/ZO98048.Google Scholar
- Bybee S, Bracken-Grissom HRH, Felder DL, Crandall KA: Directed next generation sequencing for phylogenetics: An example using Decapoda. Zool Anz. 2011, 250: 497-506. 10.1016/j.jcz.2011.05.010.Google Scholar
- Bybee SM, Bracken-Grissom HD, Haynes B, Hermansen R, Byers R, Clement M, Udall J, Crandall KA: A scalable next-gen approach to multilocus, multitaxa phylogenetics: An example of directed sequencing using Pancrustacea. Genome Biol Evol. 2011, 3: 1312-1323. 10.1093/gbe/evr106.PubMed CentralPubMedGoogle Scholar
- Song H, Buhay JE, Whiting MF, Crandall KA: Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. P Natl Acad Sci USA. 2008, 105 (36): 13486-13491. 10.1073/pnas.0803076105.Google Scholar
- Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17 (4): 540-552. 10.1093/oxfordjournals.molbev.a026334.PubMedGoogle Scholar
- Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. Version 273. 2010, http://mesquiteproject.org,Google Scholar
- Stamatakis A, Ludwig T, Meier H: RAxML-III: A fast program for maximum likelihood-based inference of large phylogenetic trees. Bioinformatics. 2005, 21 (4): 456-463. 10.1093/bioinformatics/bti191.PubMedGoogle Scholar
- Stamatakis A, Hoover P, Rougemont J: A Rapid Bootstrap Algorithm for the RAxML Web Servers. Syst Biol. 2008, 57 (5): 758-771. 10.1080/10635150802429642.PubMedGoogle Scholar
- Stamatakis A, Blagojevic F, Nikolopoulos DS, Antonopoulos CD: Exploring new search algorithms and hardware for phylogenetics: RAxML meets the IBM cell. J Vlsi Signal Proc Syst Signal Image Video Technol. 2007, 48 (3): 271-286. 10.1007/s11265-007-0067-4.Google Scholar
- Felsenstein J: Confidence-limits on phylogenies with a molecular clock. Syst Zool. 1985, 34 (2): 152-161. 10.2307/2413323.Google Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogeny. Biometrics. 2001, 17: 754-755.Google Scholar
- Lewis PO: A likelihod approach to estimating phylogeny from discrete morphological character data. Syst Biol. 2001, 50: 913-925. 10.1080/106351501753462876.PubMedGoogle Scholar
- Posada D, Crandall KA: MODELTEST: Testing the model of DNA substitution. Bioinformatics. 1998, 14 (9): 817-818. 10.1093/bioinformatics/14.9.817.PubMedGoogle Scholar
- Rambaut A, Drummond AJ: Tracer v1.4. Available from http://beast.bio.ed.ac.uk/Tracer 2007
- Shimodaira H, Hasegawa M: Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol Biol Evol. 1999, 16 (8): 1114-1116. 10.1093/oxfordjournals.molbev.a026201.Google Scholar
- Ekman S, Andersen HL, Wedin M: The limitations of ancestral state reconstruction and the evolution of the ascus in the Lecanorales (Lichenized ascomycota). Syst Biol. 2008, 57 (1): 141-156. 10.1080/10635150801910451.PubMedGoogle Scholar
- Cunningham CW, Omland KE, Oakley TH: Reconstructing ancestral character states: a critical reappraisal. Trends Ecol Evol. 1998, 13 (9): 361-366. 10.1016/S0169-5347(98)01382-2.PubMedGoogle Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMed CentralPubMedGoogle Scholar
- Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4 (5): 699-710.Google Scholar
- Lepage T, Bryant D, Philippe H, Lartillot N: A general comparison of relaxed molecular clock models. Mol Biol Evol. 2007, 24 (12): 2669-2680. 10.1093/molbev/msm193.PubMedGoogle Scholar
- Ronquist F, Klopfstein S, Vilhelmsen L, Schulmeister S, Murray DL, Rasnitsyn AP: A total-evidence approach to dating with fossils, applied to the early radiation of the hymenoptera. Syst Biol. 2012, 61 (6): 973-999. 10.1093/sysbio/sys058.PubMed CentralPubMedGoogle Scholar
- Gittleman JL, Kot M: Adaptation - Statistics and a Null Model for Estimating Phylogenetic Effects. Syst Zool. 1990, 39 (3): 227-241. 10.2307/2992183.Google Scholar
- Ho SYW: An examination of phylogenetic models of substitution rate variation among lineages. Biol Letters. 2009, 5 (3): 421-424. 10.1098/rsbl.2008.0729.Google Scholar
- Sanderson MJ: r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003, 19 (2): 301-302. 10.1093/bioinformatics/19.2.301.PubMedGoogle Scholar
- Schweitzer C, Feldmann R, Garassino A, Karasawa H, Schweigert G: Systematic List Of Fossil Decapod Crustacean Species, Crustaceana Monographs, Volume 10. 2010, Leiden, Boston: BrillGoogle Scholar
- Parham J, Donoghue P, Bell C, Calway T, Head J, Holroyd P, Inou J, Irmis R, et al: Best practices for justifying fossil calibrations. Syst Biol. 2012, 61 (2): 346-359. 10.1093/sysbio/syr107.PubMed CentralPubMedGoogle Scholar
- Robins CM, Feldmann RM, Schweitzer CE: The oldest Munididae (Decapoda: Anomura: Galatheoidea) from Ernstbrunn, Austria (Tithonian). Annalen des naturhistorischen museum in wien. 2012, 114: 289-300.Google Scholar
- Ho SYW: Calibrating molecular estimates of substitution rates and divergence times in birds. J Avian Biol. 2007, 38: 409-414.Google Scholar
- Machordom A, Macpherson E: Rapid radiation and cryptic speciation in squat lobsters of the genus Munida (Crustacea, Decapoda) and related genera in the South West Pacific: molecular and morphological evidence. Mol Phylogenet Evol. 2004, 33 (2): 259-279. 10.1016/j.ympev.2004.06.001.PubMedGoogle Scholar
- Cabezas P, Macpherson E, Machordom A: A new genus of squat lobster (Decapoda: Anomura: Galatheidae) from the South West Pacific and Indian Ocean inferred from morphological and molecular evidence. J Crustacean Biol. 2008, 28 (1): 68-75. 10.1651/07-2823R.1.Google Scholar
- Tsang LM, Lin FJ, Chu KH, Chan TY: Phylogeny of Thalassinidea (Crustacea, Decapoda) inferred from three rDNA sequences: implications for morphological evolution and superfamily classification. J Zoo Syst Evo Res. 2008, 46 (3): 216-223. 10.1111/j.1439-0469.2008.00459.x.Google Scholar
- Ayón-Parente M, Hendrickx ME: A new genus and new species of hermit crab (Crustacea: Anomura: Paguroidea: Diogenidae) from the eastern tropical Pacific. Zootaxa. 2010, 2677: 49-59.Google Scholar
- Osawa M, McLaughlin PA: Annotated checklist of anomuran decapod crustaceans of the world (exclusive of the Kiwaoidea and families Chirostylidae and Galatheidae of the Galatheoidea) Part II–Porcellanidae. Checklists of anomuran decapod crustaceans of the world (exclusive of the Kiwaoidea and families Chirostylidae and Galatheidae of the Galatheoidea) and marine lobsters of the world Raffles Bulletin of Zoology, Supplement. 2010, 23: 109-129.Google Scholar
- Macpherson E, Baba K: Taxonomy of squat lobsters. The Biology of Squat Lobsters. 2011, 39-Google Scholar
- Boyko CB, McLaughlin PA: Annotated Checklist of Anomuran Decapod Crustaceans of the World (Exclusive of the Kiwaoidea and Families Chirostylidae and Galatheidae of the Galatheoidea) Part Iv - Hippoidea. Raffles B Zool. 2010, Suppl. 23: 139-151.Google Scholar
- Mclaughlin PA, Komai T, Lemaitre R, Rahayu DL: Annotated checklist of anomuran decapod crustaceans of the world (exclusive of the Kiwaoidea and families Chirostylidae and Galatheidae of the Galatheoidea) Part I – Lithodoidea, Lomisoidea and Paguroidea. The Raffles Bulletin of Zoology. 2010, 23 (i): 5-107.Google Scholar
- Mclaughlin PA, Lemaitre R, Crandall KA: Annotated checklist of anomuran decapod crustaceans of the world (exclusive of the Kiwaoidea and families Chirostylidae and Galatheidae of the Galatheoidea) Part III – Aegloidea. The Raffles Bulletin of Zoology. 2010, 23 (i): 131-137.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.