Neuroanatomy of mud dragons: a comprehensive view of the nervous system in Echinoderes (Kinorhyncha) by confocal laser scanning microscopy

Background

The Scalidophora (Kinorhyncha, Loricifera and Priapulida) have an important phylogenetic position as early branching ecdysozoans, yet the architecture of their nervous organ systems is notably underinvestigated. Without such information, and in the absence of a stable phylogenetic context, we are inhibited from producing adequate hypotheses about the evolution and diversification of ecdysozoan nervous systems. Here, we utilize confocal laser scanning microscopy to characterize serotonergic, tubulinergic and FMRFamidergic immunoreactivity patterns in a comparative neuroanatomical study with three species of Echinoderes, the most speciose, abundant and diverse genus within Kinorhyncha.

Results

Neuroanatomy in Echinoderes as revealed by acetylated α-tubulin immunoreactivity includes a circumpharyngeal brain and ten neurite bundles in the head region that converge into five longitudinal nerves within the trunk. The ventral nerve cord is ganglionated, emerging from the brain with two connectives that converge in trunk segments 2–3, and diverge again within segment 8. The longitudinal nerves and ventral nerve cord are connected by two transverse neurites in segments 2–9. Differences among species correlate with the number, position and innervation of cuticular structures along the body. Patterns of serotoninergic and FMRFamidergic immunoreactivity correlate with the position of the brain neuropil and the ventral nerve cord. Distinct serotonergic and FMRFamidergic somata are associated with the brain neuropil and specific trunk segments along the ventral nerve cord.

Conclusions

Neural architecture is highly conserved across all three species, suggesting that our results reveal a pattern that is common to more than 40% of the species within Kinorhyncha. The nervous system of Echinoderes is segmented along most of the trunk; however, posterior trunk segments exhibit modifications that are likely associated with sensorial, motor or reproductive functions. Although all kinorhynchs show some evidence of an externally segmented trunk, it is unclear whether external segmentation matches internal segmentation of nervous and muscular organ systems across Kinorhyncha, as we observed in Echinoderes. The neuroanatomical data provided in this study not only expand the limited knowledge on kinorhynch nervous systems but also establish a comparative morphological framework within Scalidophora that will support broader inferences about the evolution of neural architecture among the deepest branching lineages of the Ecdysozoa.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1405-4) contains supplementary material, which is available to authorized users.

Miniaturization of tardigrades (water bears): Morphological and genomic perspectives

Tardigrades form a monophyletic group of microscopic ecdysozoans best known for surviving extreme environmental conditions. Due to their key phylogenetic position as a subgroup of the Panarthropoda, understanding tardigrade biology is important for comparative studies with related groups like Arthropoda. Panarthropods - and Ecdysozoa as a whole - likely evolved from macroscopic ancestors, with several taxa becoming secondarily miniaturized. Morphological and genomic evidence likewise points to a miniaturized tardigrade ancestor. The five-segmented tardigrade body typically measures less than 1 mm in length and consists of only about 1000 cells. Most organs comprise a relatively small number of cells, with the highest proportion belonging to the central nervous system, while muscles are reduced to a single cell each. Similarly, fully sequenced genomes of three tardigrade species - together with Hox gene expression data - point to extensive modifications, rearrangements, and major losses of genes and even a large body region. Parallels are evident with related ecdysozoans that may have also undergone genomic reductions, such as the nematode Caenorhabditis elegans. We interpret these data together as evidence of miniaturization in the tardigrade lineage, while cautioning that the effects of miniaturization may manifest in different ways depending on the organ or organ system under examination.

A transcriptome approach to ecdysozoan phylogeny

The monophyly of Ecdysozoa, which comprise molting phyla, has received strong support from several lines of evidence. However, the internal relationships of Ecdysozoa are still contended. We generated expressed sequence tags from a priapulid (penis worm), a kinorhynch (mud dragon), a tardigrade (water bear) and five chelicerate taxa by 454 transcriptome sequencing. A multigene alignment was assembled from 63 taxa, which comprised after matrix optimization 24,249 amino acid positions with high data density (2.6% gaps, 19.1% missing data). Phylogenetic analyses employing various models support the monophyly of Ecdysozoa. A clade combining Priapulida and Kinorhyncha (i.e. Scalidophora) was recovered as the earliest branch among Ecdysozoa. We conclude that Cycloneuralia, a taxon erected to combine Priapulida, Kinorhyncha and Nematoda (and others), are paraphyletic. Rather Arthropoda (including Onychophora) are allied with Nematoda and Tardigrada. Within Arthropoda, we found strong support for most clades, including monophyletic Mandibulata and Pancrustacea. The phylogeny within the Euchelicerata remained largely unresolved. There is conflicting evidence on the position of tardigrades: While Bayesian and maximum likelihood analyses of only slowly evolving genes recovered Tardigrada as a sister group to Arthropoda, analyses of the full data set, and of subsets containing genes evolving at fast and intermediate rates identified a clade of Tardigrada and Nematoda. Notably, the latter topology is also supported by the analyses of indel patterns.

The oldest known priapulid-like scalidophoran animal and its implications for the early evolution of cycloneuralians and ecdysozoans

Morphological phylogenetic analyses suggest that scalidophorans (priapulids, loriciferans, and kinorhynchs) and nematoids (nematodes and nematomorphs) form the ecdysozoan clade Cycloneuralia, which is a sister group to panarthropods. It has been proposed that extant priapulids and Cambrian priapulid-like scalidophorans, because of their conserved evolution, have the potential to illuminate the ancestral morphology, ecology, and developmental biology of highly derived ecdysozoans such as nematods and arthropods. As such, Cambrian fossils, particularly Markuelia and possibly olivooids, can inform the early evolution of scalidophorans, cycloneuralians, and ecdysozoans. However, the scalidophoran Markuelia is known exclusively as embryo fossils, and the olivooids have been alternatively interpreted as cnidarians or cycloneuralians. Here, we describe a post-embryonic scalidophoran fossil Eopriapulites sphinx new genus and species, which represents the oldest known scalidophoran, from the early Cambrian Period (∼535 Ma) in South China. E. sphinx is similar to modern scalidophorans in having an introvert armed with hollow scalids, a collar with coronal scalids, and a pharynx with pharyngeal teeth, but its scalids and pharyngeal teeth are arranged in a hexaradial pattern. Phylogenetically resolved as a stem-group scalidophoran, E. sphinx shares a hexaradial pattern with the hexaradial arrangement of certain anatomical structures in kinorhynchs, loriciferans, nematoids, and Cambrian fossils such as Eolympia pediculata, which could also be a scalidophoran. Thus, the bodyplan of ancestral cycloneuralians may have had a component of hexaradial symmetry (i.e., some but not necessarily all anatomical parts are hexaradially arranged). If panarthropods are nested within paraphyletic cycloneuralians, as several molecular phylogenetic analyses suggest, the ancestral ecdysozoans may have been a legless worm possibly with a component of hexaradial symmetry.

Comparative myoanatomy of Echinoderes (Kinorhyncha): a comprehensive investigation by CLSM and 3D reconstruction

INTRODUCTION

Kinorhyncha is a clade of marine invertebrate meiofauna. Their body plan includes a retractable introvert bearing rings of cuticular spines, and a limbless trunk with distinct segmentation of nervous, muscular and epidermal organ systems. As derived members within the basal branch of Ecdysozoa, kinorhynchs may provide an important example of convergence on the evolution of segmentation within one of three bilaterian superclades. We describe the myoanatomy of Echinoderes, the most specious kinorhynch genus, and build upon historical studies of kinorhynch ultrastructure and gross morphology. This is the first multi-species comparison of a complete organ system by confocal microscopy and three-dimensional reconstruction within Kinorhyncha.

RESULTS

Myoanatomy of adult Echinoderes is composed of the following: Head with two mouth cone circular muscles, nine pairs of oral style muscles, ten introvert retractors, one introvert circular muscle, and fourteen introvert circular muscle retractors; Neck with one circular muscle; Trunk showing distinct pairs of ventral and dorsal muscles within segments 1-10, dorsoventral muscles within segments 3-10, diagonal muscles within segments 1-8, longitudinal fibers spanning segments 1-9, three pairs of terminal spine muscles, and one pair of male penile spine muscles; Gut showing a pharynx with ten alternating rings of radial and circular muscle fibers enclosed in a complex sheath of protractors and retractors, an orthogonal grid of longitudinal and circular fibers surrounding the intestine, and paired hindgut dilators.

CONCLUSIONS

Myoanatomy is highly conserved between species of Echinoderes. Interspecific variation is observed in the arrangement and number of introvert fibers and the composition of pharyngeal muscles. Segmented trunk musculature facilitates the movements of articulated cuticular plates along the anterior-posterior axis. Intersegmental muscle fibers assist with dorsoventral and lateral trunk movements. Protractors, retractors and circular muscles coordinate eversion and retraction of the introvert and mouth cone, and relocation of the pharynx during locomotion and feeding behaviors. Pairs of posterior fibers suggest independent movements of terminal spines, and male penile spines. Within Scalidophora, myoanatomy is more similar between Kinorhyncha and Loricifera, than either group is to Priapulida. Kinorhynch myoanatomy may reflect a convergent transition from vermiform to segmented body plans during the early radiation of Ecdysozoa.

Comparative morphology of serotonergic-like immunoreactive elements in the central nervous system of kinorhynchs (Kinorhyncha, Cyclorhagida)

Cycloneuralian taxa exhibit similar organ system architectures, providing informative characters of metazoan evolution, yet very few modern comparative descriptions of cellular and molecular homologies within and among those taxa are available. We immunolabeled and characterized elements of the serotonergic nervous system in the kinorhynchs Echinoderes spinifurca, Antygomonas paulae, and Zelinkaderes brightae using confocal laser scanning microscopy. Fluorescent markers targeting DNA were combined with observations of auto-fluorescent structures to guide interpretations of the internal and external anatomy in each species. Results show a common pattern of the central nervous system with a circumenteric brain divided into ring-shaped anterior and posterior neuronal somata and a central neuropil connected to a multi-stringed, longitudinal ventral nerve cord. Structural similarities and differences in the nervous systems of these species were observed and described, stressing the incomplete ring nature of the anterior region of the kinorhynch brain, the functional relationship between the brain and the movable introvert, and the number and arrangement of nerve strings and somata of the ventral nerve cord. The ventral cord ends in two ventrolateral cell bodies in E. spinifurca, and forms a terminal loop associated with a midterminal spine in A. paulae and Z. brightae. The possible functional and phylogenetic significance of these features and arrangements are discussed.

Detection of a troponin I-like protein in non-striated muscle of the tardigrades (water bears)

Tardigrades, also known as water bears, have somatic muscle fibers that are responsible for movement of their body and legs. These muscle fibers contain thin and thick filaments in a non-striated pattern. However, the regulatory mechanism of muscle contraction in tardigrades is unknown. In the absence of extensive molecular and genomic information, we detected a protein of 31 kDa in whole lysates of tardigrades that cross-reacted with the antibody raised against nematode troponin I (TnI). TnI is a component of the troponin complex that regulates actin-myosin interaction in a Ca(2+)-dependent and actin-linked manner. This TnI-like protein was co-extracted with actin in a buffer containing ATP and EGTA, which is known to induce relaxation of a troponin-regulated contractile system. The TnI-like protein was specifically expressed in the somatic muscle fibers in adult animals and partially co-localized with actin filaments in a non-striated manner. Interestingly, the pharyngeal muscle did not express this protein. These observations suggest that the non-striated somatic muscle of tardigrades has an actin-linked and troponin-regulated system for muscle contraction.

A congruent solution to arthropod phylogeny: phylogenomics, microRNAs and morphology support monophyletic Mandibulata

While a unique origin of the euarthropods is well established, relationships between the four euarthropod classes-chelicerates, myriapods, crustaceans and hexapods-are less clear. Unsolved questions include the position of myriapods, the monophyletic origin of chelicerates, and the validity of the close relationship of euarthropods to tardigrades and onychophorans. Morphology predicts that myriapods, insects and crustaceans form a monophyletic group, the Mandibulata, which has been contradicted by many molecular studies that support an alternative Myriochelata hypothesis (Myriapoda plus Chelicerata). Because of the conflicting insights from published molecular datasets, evidence from nuclear-coding genes needs corroboration from independent data to define the relationships among major nodes in the euarthropod tree. Here, we address this issue by analysing two independent molecular datasets: a phylogenomic dataset of 198 protein-coding genes including new sequences for myriapods, and novel microRNA complements sampled from all major arthropod lineages. Our phylogenomic analyses strongly support Mandibulata, and show that Myriochelata is a tree-reconstruction artefact caused by saturation and long-branch attraction. The analysis of the microRNA dataset corroborates the Mandibulata, showing that the microRNAs miR-965 and miR-282 are present and expressed in all mandibulate species sampled, but not in the chelicerates. Mandibulata is further supported by the phylogenetic analysis of a comprehensive morphological dataset covering living and fossil arthropods, and including recently proposed, putative apomorphies of Myriochelata. Our phylogenomic analyses also provide strong support for the inclusion of pycnogonids in a monophyletic Chelicerata, a paraphyletic Cycloneuralia, and a common origin of Arthropoda (tardigrades, onychophorans and arthropods), suggesting that previous phylogenies grouping tardigrades and nematodes may also have been subject to tree-reconstruction artefacts.

Myoanatomy of the marine tardigrade Halobiotus crispae (Eutardigrada: Hypsibiidae)

The muscular architecture of Halobiotus crispae (Eutardigrada: Hypsibiidae) was examined by means of fluorescent-coupled phalloidin in combination with confocal laser scanning microscopy and computer-aided three-dimensional reconstruction, in addition to light microscopy (Nomarski), scanning electron microscopy, and transmission electron microscopy (TEM). The somatic musculature of H. crispae is composed of structurally independent muscle fibers, which can be divided into a dorsal, ventral, dorsoventral, and a lateral musculature. Moreover, a distinct leg musculature is found. The number and arrangement of muscles differ in each leg. Noticeably, the fourth leg contains much fewer muscles when compared with the other legs. Buccopharyngeal musculature (myoepithelial muscles), intestinal musculature, and cloacal musculature comprise the animal's visceral musculature. TEM of stylet and leg musculature revealed ultrastructural similarities between these two muscle groups. Furthermore, microtubules are found in the epidermal cells of both leg and stylet muscle attachments. This would indicate that the stylet and stylet glands are homologues to the claw and claw glands, respectively. When comparing with previously published data on both heterotardigrade and eutardigrade species, it becomes obvious that eutardigrades possess very similar numbers and arrangement of muscles, yet differ in a number of significant details of their myoanatomy. This study establishes a morphological framework for the use of muscular architecture in elucidating tardigrade phylogeny.

Lophotrochozoa internal phylogeny: new insights from an up-to-date analysis of nuclear ribosomal genes

Resolving the relationships among animal phyla is a key biological problem that remains to be solved. Morphology is unable to determine the relationships among most phyla and although molecular data have unveiled a new evolutionary scenario, they have their own limitations. Nuclear ribosomal genes (18S and 28S rDNA) have been used effectively for many years. However, they are considered of limited use for resolving deep divergences such as the origin of the bilaterians, due to certain drawbacks such as the long-branch attraction (LBA) problem. Here, we attempt to overcome these pitfalls by combining several methods suggested in previous studies and routinely used in contemporary standard phylogenetic analyses but that have not yet been applied to any bilaterian phylogeny based on these genes. The methods used include maximum likelihood and Bayesian inference, the application of models with rate heterogeneity across sites, wide taxon sampling and compartmentalized analyses for each problematic clade. The results obtained show that the combination of the above-mentioned methodologies minimizes the LBA effect, and a new Lophotrochozoa phylogeny emerges. Also, the Acoela and Nemertodermatida are confirmed with maximum support as the first branching bilaterians. Ribosomal RNA genes are thus a reliable source for the study of deep divergences in the metazoan tree, provided that the data are treated carefully.

Comments on the eyes of tardigrades

A survey is given on the scarce information on the visual organs (eyes or ocelli) of Tardigrada. Many Eutardigrada and some Arthrotardigrada, namely the Echiniscidae, possess inverse pigment-cup ocelli, which are located in the outer lobe of the brain, and probably are of cerebral origin. Occurrence of such organs in tardigrades, suggested as being eyeless, has never been checked. Depending on the species, response to light (photokinesis) is negative, positive or indifferent, and may change during the ontogeny. The tardigrade eyes of the two eutardigrades examined up to now comprise a single pigment cup cell, one or two microvillous (rhabdomeric) sensory cells and ciliary sensory cell(s). In the eyes of the eutardigrade Milnesium tardigradum the cilia are differentiated in an outer branching segment and an inner (dendritic) segment. Because of the scarcity of information on the tardigrade eyes, their homology with the visual organs of other bilaterians is currently difficult to establish and further comparative studies are needed. Thus, the significance of these eyes for the evolution of arthropod visual systems is unclear yet.

The tardigrade Hypsibius dujardini, a new model for studying the evolution of development

Studying development in diverse taxa can address a central issue in evolutionary biology: how morphological diversity arises through the evolution of developmental mechanisms. Two of the best-studied developmental model organisms, the arthropod Drosophila and the nematode Caenorhabditis elegans, have been found to belong to a single protostome superclade, the Ecdysozoa. This finding suggests that a closely related ecdysozoan phylum could serve as a valuable model for studying how developmental mechanisms evolve in ways that can produce diverse body plans. Tardigrades, also called water bears, make up a phylum of microscopic ecdysozoan animals. Tardigrades share many characteristics with C. elegans and Drosophila that could make them useful laboratory models, but long-term culturing of tardigrades historically has been a challenge, and there have been few studies of tardigrade development. Here, we show that the tardigrade Hypsibius dujardini can be cultured continuously for decades and can be cryopreserved. We report that H. dujardini has a compact genome, a little smaller than that of C. elegans or Drosophila, and that sequence evolution has occurred at a typical rate. H. dujardini has a short generation time, 13-14 days at room temperature. We have found that the embryos of H. dujardini have a stereotyped cleavage pattern with asymmetric cell divisions, nuclear migrations, and cell migrations occurring in reproducible patterns. We present a cell lineage of the early embryo and an embryonic staging series. We expect that these data can serve as a platform for using H. dujardini as a model for studying the evolution of developmental mechanisms.

Segmental expression of Pax3/7 and engrailed homologs in tardigrade development

How morphological diversity arises through evolution of gene sequence is a major question in biology. In Drosophila, the genetic basis for body patterning and morphological segmentation has been studied intensively. It is clear that some of the genes in the Drosophila segmentation program are functioning similarly in certain other taxa, although many questions remain about when these gene functions arose and which taxa use these genes similarly to establish diverse body plans. Tardigrades are an outgroup to arthropods in the Ecdysozoa and, as such, can provide insight into how gene functions have evolved among the arthropods and their close relatives. We developed immunostaining methods for tardigrade embryos, and we used cross-reactive antibodies to investigate the expression of homologs of the pair-rule gene paired (Pax3/7) and the segment polarity gene engrailed in the tardigrade Hypsibius dujardini. We find that in H. dujardini embryos, Pax3/7 protein localizes not in a pair-rule pattern but in a segmentally iterated pattern, after the segments are established, in regions of the embryo where neurons later arise. Engrailed protein localizes in the posterior ectoderm of each segment before ectodermal segmentation is apparent. Together with previous results from others, our data support the conclusions that the pair-rule function of Pax3/7 is specific to the arthropods, that some of the ancient functions of Pax3/7 and Engrailed in ancestral bilaterians may have been in neurogenesis, and that Engrailed may have a function in establishing morphological boundaries between segments that is conserved at least among the Panarthropoda.

A comparative genomics analysis of codon reassignments reveals a link with mitochondrial proteome size and a mechanism of genetic code change via suppressor tRNAs

Using a comparative genomics approach we demonstrate a negative correlation between the number of codon reassignments undergone by 222 mitochondrial genomes and the mitochondrial genome size, the number of mitochondrial ORFs, and the sizes of the large and small subunit mitochondrial rRNAs. In addition, we show that the TGA-to-tryptophan codon reassignment, which has occurred 11 times in mitochondrial genomes, is found in mitochondrial genomes smaller than those which have not undergone the reassignment. We therefore propose that mitochondrial codon reassignments occur in a wide range of phyla, particularly in Metazoa, due to a reduced "proteomic constraint" on the mitochondrial genetic code, compared to the nuclear genetic code. The reduced proteomic constraint reflects the small size of the mitochondrial-encoded proteome and allows codon reassignments to occur with less likelihood of lethality. In addition, we demonstrate a striking link between nonsense codon reassignments and the decoding properties of naturally occurring nonsense suppressor tRNAs. This suggests that natural preexisting nonsense suppression facilitated nonsense codon reassignments and constitutes a novel mechanism of genetic code change. These findings explain for the first time the identity of the stop codons and amino acids reassigned in mitochondrial and nuclear genomes. Nonsense suppressor tRNAs provided the raw material for nonsense codon reassignments, implying that the properties of the tRNA anticodon have dictated the identity of nonsense codon reassignments.

Further use of nearly complete 28S and 18S rRNA genes to classify Ecdysozoa: 37 more arthropods and a kinorhynch

This work expands on a study from 2004 by Mallatt, Garey, and Shultz [Mallatt, J.M., Garey, J.R., Shultz, J.W., 2004. Ecdysozoan phylogeny and Bayesian inference: first use of nearly complete 28S and 18S rRNA gene sequences to classify the arthropods and their kin. Mol. Phylogenet. Evol. 31, 178-191] that evaluated the phylogenetic relationships in Ecdysozoa (molting animals), especially arthropods. Here, the number of rRNA gene-sequences was effectively doubled for each major group of arthropods, and sequences from the phylum Kinorhyncha (mud dragons) were also included, bringing the number of ecdysozoan taxa to over 80. The methods emphasized maximum likelihood, Bayesian inference and statistical testing with parametric bootstrapping, but also included parsimony and minimum evolution. Prominent findings from our combined analysis of both genes are as follows. The fundamental subdivisions of Hexapoda (insects and relatives) are Insecta and Entognatha, with the latter consisting of collembolans (springtails) and a clade of proturans plus diplurans. Our rRNA-gene data provide the strongest evidence to date that the sister group of Hexapoda is Branchiopoda (fairy shrimps, tadpole shrimps, etc.), not Malacostraca. The large, Pancrustacea clade (hexapods within a paraphyletic Crustacea) divided into a few basic subclades: hexapods plus branchiopods; cirripedes (barnacles) plus malacostracans (lobsters, crabs, true shrimps, isopods, etc.); and the basally located clades of (a) ostracods (seed shrimps) and (b) branchiurans (fish lice) plus the bizarre pentastomids (tongue worms). These findings about Pancrustacea agree with a recent study by Regier, Shultz, and Kambic that used entirely different genes [Regier, J.C., Shultz, J.W., Kambic, R.E., 2005a. Pancrustacean phylogeny: hexapods are terrestrial crustaceans and maxillopods are not monophyletic. Proc. R. Soc. B 272, 395-401]. In Malacostraca, the stomatopod (mantis shrimp) was not at the base of the eumalacostracans, as is widely claimed, but grouped instead with an euphausiacean (krill). Within centipedes, Craterostigmus was the sister to all other pleurostigmophorans, contrary to the consensus view. Our trees also united myriapods (millipedes and centipedes) with chelicerates (horseshoe crabs, spiders, scorpions, and relatives) and united pycnogonids (sea spiders) with chelicerates, but with much less support than in the previous rRNA-gene study. Finally, kinorhynchs joined priapulans (penis worms) at the base of Ecdysozoa.

The anatomy, affinity, and phylogenetic significance of Markuelia

The fossil record provides a paucity of data on the development of extinct organisms, particularly for their embryology. The recovery of fossilized embryos heralds new insight into the evolution of development but advances are limited by an almost complete absence of phylogenetic constraint. Markuelia is an exception to this, known from cleavage and pre-hatchling stages as a vermiform and profusely annulated direct-developing bilaterian with terminal circumoral and posterior radial arrays of spines. Phylogenetic analyses have hitherto suggested assignment to stem-Scalidophora (phyla Kinorhyncha, Loricifera, Priapulida). We test this assumption with additional data and through the inclusion of additional taxa. The available evidence supports stem-Scalidophora affinity, leading to the conclusion that scalidophorans, cyclonerualians, and ecdysozoans are primitive direct developers, and the likelihood that scalidophorans are primitively metameric.

The eutardigrade Thulinia stephaniae has an indeterminate development and the potential to regulate early blastomere ablations

We present a detailed analysis of the cell lineage of the tardigrade Thulinia stephaniae with a 4D-microscopy system (3D time-lapse recording). The recording, of the entire development from embryogenesis until hatching, allowed us to analyze the fate of single descendants from early blastomeres up to germ layer formation and tissue development. The embryo undergoes an irregular indeterminate cleavage pattern without early fate restriction. During gastrulation, mesodermal and endodermal precursors, and a pair of primordial germ cells migrate through a blastopore at the prospective position of the mouth. Our results are not consistent with earlier descriptions of mesoderm formation by enterocoely in tardigrades. The mesoderm in Thulinia stephaniae originates from a variable number of blastomeres, which form mesodermal bands that later produce the serial somites. The nervous system is formed by neural progenitor cells, which delaminate from the neurogenic ectoderm. Early embryogenesis of Thulinia stephaniae is highly regulative, even after laser ablations of blastomeres at the two- and four-cell stages 'normal' juveniles are formed. This has never been observed before for a protostome. Germ cell specification occurs late during development between the sixth and seventh cell generation. Comparing the development of other protostomes with that of the Tardigrada, which occupy a basal position within the Arthropoda, suggests that an indeterminate cleavage and regulatory development is not only part of the ground pattern of the Arthropoda, but probably of the entire Ecdysozoa.

Ecdysozoan phylogeny and Bayesian inference: first use of nearly complete 28S and 18S rRNA gene sequences to classify the arthropods and their kin

Relationships among the ecdysozoans, or molting animals, have been difficult to resolve. Here, we use nearly complete 28S+18S ribosomal RNA gene sequences to estimate the relations of 35 ecdysozoan taxa, including newly obtained 28S sequences from 25 of these. The tree-building algorithms were likelihood-based Bayesian inference and minimum-evolution analysis of LogDet-transformed distances, and hypotheses were tested wth parametric bootstrapping. Better taxonomic resolution and recovery of established taxa were obtained here, especially with Bayesian inference, than in previous parsimony-based studies that used 18S rRNA sequences (or 18S plus small parts of 28S). In our gene trees, priapulan worms represent the basal ecdysozoans, followed by nematomorphs, or nematomorphs plus nematodes, followed by Panarthropoda. Panarthropoda was monophyletic with high support, although the relationships among its three phyla (arthropods, onychophorans, tardigrades) remain uncertain. The four groups of arthropods-hexapods (insects and related forms), crustaceans, chelicerates (spiders, scorpions, horseshoe crabs), and myriapods (centipedes, millipedes, and relatives)-formed two well-supported clades: Hexapoda in a paraphyletic crustacea (Pancrustacea), and 'Chelicerata+Myriapoda' (a clade that we name 'Paradoxopoda'). Pycnogonids (sea spiders) were either chelicerates or part of the 'chelicerate+myriapod' clade, but not basal arthropods. Certain clades derived from morphological taxonomy, such as Mandibulata, Atelocerata, Schizoramia, Maxillopoda and Cycloneuralia, are inconsistent with these rRNA data. The 28S gene contained more signal than the 18S gene, and contributed to the improved phylogenetic resolution. Our findings are similar to those obtained from mitochondrial and nuclear (e.g., elongation factor, RNA polymerase, Hox) protein-encoding genes, and should revive interest in using rRNA genes to study arthropod and ecdysozoan relationships.

Bayesian Inference of the Metazoan Phylogeny

Metazoan phylogeny remains one of evolutionary biology's major unsolved problems. Molecular and morphological data, as well as different analytical approaches, have produced highly conflicting results due to homoplasy resulting from more than 570 million years of evolution. To date, parsimony has been the only feasible combined approach but is highly sensitive to long-branch attraction. Recent development of stochastic models for discrete morphological characters and computationally efficient methods for Bayesian inference has enabled combined molecular and morphological data analysis with rigorous statistical approaches less prone to such inconsistencies. We present the first statistically founded analysis of a metazoan data set based on a combination of morphological and molecular data and compare the results with a traditional parsimony analysis. Interestingly, the Bayesian analyses demonstrate a high degree of congruence between morphological and molecular data, and both data sets contribute to the result of the combined analysis. Additionally, they resolve several irregularities obtained in previous studies and show high credibility values for controversial groups such as the ecdysozoans and lophotrochozoans. Parsimony, on the contrary, shows conflicting results, with morphology being congruent to the Bayesian results and the molecular data set producing peculiarities that are largely reflected in the combined analysis.

Molecules, development and fossils in the study of metazoan evolution; Articulata versus Ecdysozoa revisited

Two conflicting hypotheses of protostome relationships, Articulata and Ecdysozoa, are reviewed by evaluating the evidence in favor and against each one of them. Understanding early embryonic development and segmentation in non-arthropod non-annelid protostomes seems crucial to the debate. New ways of coding metazoan matrices, avoiding ground-patterns and higher taxa, and incorporating fossil evidence seems the best way to avoid circular debates. Molecular data served as the catalyzer for the Ecdysozoa hypothesis, although morphological support had been implicitly suggested. Most molecular analyses published so far have shown some support for Ecdysozoa, whereas none has ever supported Articulata. Here, new analyses of up to four nuclear loci, including 18S rRNA, myosin heavy chain II, histone H3 and elongation factor 1-alpha are conducted to test the molecular support for Ecdysozoa, and, at least under some parameter sets, most data sets show a clade formed by the molting animals. In contrast, support for Articulata is not found under any analytical conditions.

Testing the new animal phylogeny: first use of combined large-subunit and small-subunit rRNA gene sequences to classify the protostomes

Although the small-subunit ribosomal RNA (SSU rRNA) gene is widely used in the molecular systematics, few large-subunit (LSU) rRNA gene sequences are known from protostome animals, and the value of the LSU gene for invertebrate systematics has not been explored. The goal of this study is to test whether combined LSU and SSU rRNA gene sequences support the division of protostomes into Ecdysozoa (molting forms) and Lophotrochozoa, as was proposed by Aguinaldo et al. (1997) (Nature 387:489) based on SSU rRNA sequences alone. Nearly complete LSU gene sequences were obtained, and combined LSU + SSU sequences were assembled, for 15 distantly related protostome taxa plus five deuterostome outgroups. When the aligned LSU + SSU sequences were analyzed by tree-building methods (minimum evolution analysis of LogDet-transformed distances, maximum likelihood, and maximum parsimony) and by spectral analysis of LogDet distances, both Ecdysozoa and Lophotrochozoa were indeed strongly supported (e.g., bootstrap values >90%), with higher support than from the SSU sequences alone. Furthermore, with the LogDet-based methods, the LSU + SSU sequences resolved some accepted subgroups within Ecdysozoa and Lophotrochozoa (e.g., the polychaete sequence grouped with the echiuran, and the annelid sequences grouped with the mollusc and lophophorates)-subgroups that SSU-based studies do not reveal. Also, the mollusc sequence grouped with the sequences from lophophorates (brachiopod and phoronid). Like SSU sequences, our LSU + SSU sequences contradict older hypotheses that grouped annelids with arthropods as Articulata, that said flatworms and nematodes were basal bilateralians, and considered lophophorates, nemerteans, and chaetognaths to be deuterostomes. The position of chaetognaths within protostomes remains uncertain: our chaetognath sequence associated with that of an onychophoran, but this was unstable and probably artifactual. Finally, the benefits of combining LSU with SSU sequences for phylogenetic analyses are discussed: LSU adds signal, it can be used at lower taxonomic levels, and its core region is easy to align across distant taxa-but its base frequencies tend to be nonstationary across such taxa. We conclude that molecular systematists should use combined LSU + SSU rRNA genes rather than SSU alone.