A wide-range phylogenetic analysis of Zic proteins: implications for correlations between protein structure conservation and body plan complexity

We compared Zic homologues from a wide range of animals. Striking conservation was found in the zinc finger domains, in which an exon-intron boundary has been kept in all bilateralians but not cnidarians, suggesting that all of the bilateralian Zic genes are derived from a single gene in a bilateralian ancestor. There were additional conserved amino acid sequences, ZOC and ZF-NC. Combined analysis of the zinc finger, ZOC, and ZF-NC revealed the presence of two classes of Zic, based on the degree of protein structure conservation. The "conserved" class includes Zic proteins from the Arthropoda, Mollusca, Annelida, Echinodermata, and Chordata (vertebrates and cephalochordates), whereas the "diverged" class contains those from the Platyhelminthes, Cnidaria, Nematoda, and Chordata (urochordates). The result indicates that the ancestral bilateralian Zic protein had already acquired an entire set of conserved domains, but that this was lost and diverged in the platyhelminthes, nematodes, and urochordates.     

Phylogeny of the Metazoa Based on Morphological and 18S Ribosomal DNA Evidence

Cladistic analysis of traditional (i.e. morphological, developmental, ultrastructural) and molecular (18S rDNA) data sets (276+501 informative characters) provides a hypothesis about relationships of all meta-zoan higher taxa. Monophyly of Metazoa, Epith-eliozoa (= -03non-Porifera), Triploblastica, Mesozoa, Eutriploblastica (=Rhabditophora+Catenulida+“higher triploblasts”=Neotriploblastica, including Xeno- turbellida and Gnathostomulida), Rhabditophora, Syndermata (=“Rotifera”+Acanthocephala), Neotrichozoa (=Gastrotricha+Gnathostomulida), Nematozoa (=Nematoda+Nematomorpha), Panarthropoda (=Onychophora+Tardigrada+ Arthropoda), Cephalorhyncha, Deuterostomia, Ambulacralia (=Hemichordata+Echinodermata), Chordata, Phoronozoa (=Phoronida+“Brachiopoda”), Bryozoa, Trochozoa (=Eutrochozoa+Entoprocta+ Cycliophora), Eutrochozoa, and Chaetifera (=Annelida+ Pogonophora+Echiura) is strongly supported. Cnidaria (including Myxozoa), Ecdysozoa (=Cepha- lorhyncha + Nematozoa + Chaetognatha + Panarthropoda), Eucoelomata (=Bryozoa+Phoronozoa+Deuterostomia+Trochozoa, possibly including also Xenoturbellida), and Deuterostomia+Phoronozoa probably are monophyletic. Most traditional “phyla” are monophyletic, except for Porifera, Cnidaria (excluding Myxozoa), Platyhelminthes, Brachiopoda, and Rotifera. Three “hot” regions of the tree remain quite unresolved: basal Epitheliozoa, basal Triploblastica, and basal Neotriploblastica. A new phylogenetic classification of the Metazoa including 35 formally recognized phyla (Silicispongea, Calcispongea, Placozoa, Cnidaria, Ctenophora, Acoela, Nemertodermatida, Orthonecta, Rhombozoa, Rhabditophora, Catenulida, Syndermata, Gnathostomulida, Gastrotricha, Cephalorhyncha, Chaetognatha, Nematoda, Nematomorpha, Onychophora, Tardigrada, Arthropoda, Echinodermata, Hemichordata, Chordata, Phoronozoa, Bryozoa s. str., Xenoturbellida, Entoprocta, Cycliophora, Nemertea, Mollusca, Sipuncula, Echiura, Pogonophora, and Annelida) and few i ncertae sedis g roups (e.g. Myzostomida and Lobatocerebromorpha) is proposed.     

Chordate evolution and the three-phylum system

Traditional metazoan phylogeny classifies the Vertebrata as a subphylum of the phylum Chordata, together with two other subphyla, the Urochordata (Tunicata) and the Cephalochordata. The Chordata, together with the phyla Echinodermata and Hemichordata, comprise a major group, the Deuterostomia. Chordates invariably possess a notochord and a dorsal neural tube. Although the origin and evolution of chordates has been studied for more than a century, few authors have intimately discussed taxonomic ranking of the three chordate groups themselves. Accumulating evidence shows that echinoderms and hemichordates form a clade (the Ambulacraria), and that within the Chordata, cephalochordates diverged first, with tunicates and vertebrates forming a sister group. Chordates share tadpole-type larvae containing a notochord and hollow nerve cord, whereas ambulacrarians have dipleurula-type larvae containing a hydrocoel. We propose that an evolutionary occurrence of tadpole-type larvae is fundamental to understanding mechanisms of chordate origin. Protostomes have now been reclassified into two major taxa, the Ecdysozoa and Lophotrochozoa, whose developmental pathways are characterized by ecdysis and trochophore larvae, respectively. Consistent with this classification, the profound dipleurula versus tadpole larval differences merit a category higher than the phylum. Thus, it is recommended that the Ecdysozoa, Lophotrochozoa, Ambulacraria and Chordata be classified at the superphylum level, with the Chordata further subdivided into three phyla, on the basis of their distinctive characteristics.     

An evolutionary comparison of leucine-rich repeat containing G protein-coupled receptors reveals a novel LGR subtype

Leucine-rich repeat containing G protein-coupled receptors or LGRs are receptors with important functions in development and reproduction. Belonging to this evolutionarily conserved group of receptors are the well-studied glycoprotein hormone receptors and relaxin receptors in mammals, as well as the bursicon receptor, which triggers cuticle hardening and tanning in freshly enclosed insects. In this study, the numerous LGR sequences in different animal phyla are analyzed and compared. Based on these data a phylogenetic tree was generated. This information sheds new light on structural and evolutionary aspects regarding this receptor group. Apart from vertebrates and insects, LGRs are also present in early chordates (Urochordata, Cephalochordata and Hyperoartia) and other arthropods (Arachnida and Branchiopoda) as well as in Mollusca, Echinodermata, Hemichordata, Nematoda, and even in ancient animal life forms, such as Cnidaria and Placozoa. Three distinct types of LGR exist, distinguishable by their number of leucine-rich repeats (LRRs), their type-specific hinge region and the presence or absence of an LDLa motif. Type C LGRs containing only one LDLa (C1 subtype) appear to be present in nearly all animal phyla. We here describe a second subtype, C2, containing multiple LDLa motifs, which was discovered in echinoderms, mollusks and in one insect species (Pediculus humanis corporis). In addition, eight putative LGRs can be predicted from the genome data of the placozoan species Trichoplax adhaerens. They may represent an ancient form of the LGRs, however, more genomic data will be required to confirm this hypothesis.     

The carboxylesterase/cholinesterase gene family in invertebrate deuterostomes

Carboxylesterase/cholinesterase family members are responsible for controlling the nerve impulse, detoxification and various developmental functions, and are a major target of pesticides and chemical warfare agents. Comparative structural analysis of these enzymes is thus important. The invertebrate deuterostomes (phyla Echinodermata and Hemichordata and subphyla Urochordata and Cephalochordata) lie in the transition zone between invertebrates and vertebrates, and are thus of interest to the study of evolution. Here we have investigated the carboxylesterase/cholinesterase gene family in the sequenced genomes of Strongylocentrotus purpuratus (Echinodermata), Saccoglossus kowalevskii (Hemichordata), Ciona intestinalis (Urochordata) and Branchiostoma floridae (Cephalochordata), using sequence analysis of the catalytic apparatus and oligomerisation domains, and phylogenetic analysis. All four genomes show blurring of structural boundaries between cholinesterases and carboxylesterases, with many intermediate enzymes. Non-enzymatic proteins are well represented. The Saccoglossus and Branchiostoma genomes show evidence of extensive gene duplication and retention. There is also evidence of domain shuffling, resulting in multidomain proteins consisting either of multiple carboxylesterase domains, or of carboxylesterase/cholinesterase domains linked to other domains, including RING finger, chitin-binding, immunoglobulin, fibronectin type 3, CUB, cysteine-rich-Frizzled, caspase activation and 7tm-1, amongst others. Such gene duplication and domain shuffling in the carboxylesterase/cholinesterase family appears to be unique to the invertebrate deuterostomes, and we hypothesise that these factors may have contributed to the evolution of the morphological complexity, particularly of the nervous system and neural crest, of the vertebrates.     

A novel amphioxus cadherin that localizes to epithelial adherens junctions has an unusual domain organization with implications for chordate phylogeny

Although data are available from only vertebrates, urochordates, and three nonchordate animals, there are definite differences in the structures of classic cadherins between vertebrates plus urochordates and nonchordates. In this study we examined structural diversity of classic cadherins among bilaterian animals by obtaining new data from an amphioxus (Cephalochordata, Chordata), an acorn worm (Hemichordata), a sea star (Echinodermata), and an oyster (Mollusca). The structures of newly identified nonchordate cadherins are grouped together with those of the known sea urchin and Drosophila cadherins, whereas the structure of an amphioxus (Branchiostoma belcheri) cadherin, designated BbC, is differently categorized from those of other known chordate cadherins. BbC is identified as a cadherin by its cytoplasmic domain whose sequence is highly related to the cytoplasmic sequences of all known classic cadherins, but it lacks all of the five repeats constituting the extracellular homophilic-binding domain of other chordate cadherins. The ectodomains of BbC match the ectodomains found in nonchordate cadherins but not present in other chordate cadherins. We show that the BbC functions as a cell-cell adhesion molecule when expressed in Drosophila S2 cells and localizes to adherens junctions in the ectodermal epithelia in amphioxus embryos. We argue that BbC is the amphioxus homologue of the classic cadherins involved in the formation of epithelial adherens junctions. The structural relationships of the cadherin molecules allow us to propose a possibility that cephalochordates might be basal to the sister-groups vertebrates and urochordates.     

EvolMarkers: a database for mining exon and intron markers for evolution, ecology and conservation studies

Recent innovations in next-generation sequencing have lowered the cost of genome projects. Nevertheless, sequencing entire genomes for all representatives in a study remains expensive and unnecessary for most studies in ecology, evolution and conservation. It is still more cost-effective and efficient to target and sequence single-copy nuclear gene markers for such studies. Many tools have been developed for identifying nuclear markers, but most of these have focused on particular taxonomic groups. We have built a searchable database, EvolMarkers, for developing single-copy coding sequence (CDS) and exon-primed-intron-crossing (EPIC) markers that is designed to work across a broad range of phylogenetic divergences. The database is made up of single-copy CDS derived from BLAST searches of a variety of metazoan genomes. Users can search the database for different types of markers (CDS or EPIC) that are common to different sets of input species with different divergence characteristics. EvolMarkers can be applied to any taxonomic group for which genome data are available for two or more species. We included 82 genomes in the first version of EvolMarkers and have found the methods to be effective across Placozoa, Cnidaria, Arthropod, Nematoda, Annelida, Mollusca, Echinodermata, Hemichordata, Chordata and plants. We demonstrate the effectiveness of searching for CDS markers within annelids and show how to find potentially useful intronic markers within the lizard Anolis.     

Geologic events coupled with Pleistocene climatic oscillations drove genetic variation of Omei treefrog (Rhacophorus omeimontis) in southern China

Collagens provide structural support and guidance cues within the extracellular matrix of metazoans. Mammalian collagens XIII, XXIII and XXV form a unique subgroup of type II transmembrane proteins, each comprising a short N-terminal cytosolic domain, a transmembrane domain and a largely collagenous ectodomain. We name these collagens as MACITs (Membrane-Associated Collagens with Interrupted Triple-helices), and here investigate their evolution and conserved properties. To date, these collagens have been studied only in mammals. Knowledge of the representation of MACITs in other extant metazoans is lacking. This question is of interest for understanding structural/functional relationships in the MACIT family and also for insight into the evolution of MACITs in relation to the secreted, fibrillar collagens that are present throughout the metazoa. MACITs are restricted to bilaterians and are represented in the Ecdysozoa, Hemichordata, Urochordata and Vertebrata (Gnathostomata). They were not identified in available early-diverging metazoans, Lophotrochozoa, Echinodermata, Cephalochordata or Vertebrata (Cyclostomata). Whereas invertebrates encode a single MACIT, collagens XIII/XXIII/XXV of jawed vertebrates are paralogues that originated from the two rounds of en-bloc genome duplication occurring early in vertebrate evolution. MACITs have conserved domain architecture in which a juxta-membrane furin-cleavage site and the C-terminal 34 residues are especially highly conserved, whereas the cytoplasmic domains are weakly conserved. To study protein expression and function in a metazoan with a single MACIT gene, we focused on Caenorhabditis elegans and its col-99 gene. A col-99 cDNA was cloned and expressed as protein in mammalian CHO cells, two antibodies against COL-99 protein were generated, and a col-99-bearing fosmid gene construct col-99::egfp::flag was used to generate transgenic C. elegans lines. The encoded COL-99 polypeptide is 85 kDa in size and forms a trimeric protein. COL-99 is plasma membrane-associated and undergoes furin-dependent ectodomain cleavage and shedding. COL-99 is detected in mouth, pharynx, body wall and the tail, mostly in motor neurons and muscle systems and is enriched at neuromuscular junctions. Through identification of MACITs in multiple metazoan phyla we developed a model for the evolution of MACITs. The experimental data demonstrate conservation of MACIT molecular and cellular properties and tissue localisations in the invertebrate, C. elegans.     

Structure and Function of Metazoan Ciliary Bands and Their Phylogenetic Significance

Ciliated epithelia, especially the ciliary bands used in swimming and filter feeding, of representatives of the following phyla have been investigated: Porifera, Cnidaria, Annelida, Mollusca, Sipuncula, Nemertini, Platyhelminthes, Entoprocta, Ectoprocta, Rotifera, Pterobranchia, Phoronida, Brachiopoda, Echinodermata and Enteropneusta. The trochaea theory predicts that Porifera and Cnidaria have only monociliate cells and lack ciliary bands used in filter-feeding, that the gastroneuralian larvae have downstream-collecting ciliary bands with prototroch and metatroch of compound cilia on multiciliate cells, and that notoneuralian larvae have an upstream-collecting neotroch on monociliate cells. The observations generally fit these predictions and the exceptions are discussed. In all the ciliated epithelia, except that of the sponge larva, each ciliated cell has an accessory centriole perpendicular to the basal body of the cilium and situated on its downstream side.     

Gluconeogenesis: An ancient biochemical pathway with a new twist

Synthesis of sugars from simple carbon sources is critical for survival of animals under limited nutrient availability. Thus, sugar-synthesizing enzymes should be present across the entire metazoan spectrum. Here, we explore the evolution of glucose and trehalose synthesis using a phylogenetic analysis of enzymes specific for the two pathways. Our analysis reveals that the production of trehalose is the more ancestral biochemical process, found in single cell organisms and primitive metazoans, but also in insects. The gluconeogenic-specific enzyme glucose-6-phosphatase (G6Pase) first appears in Cnidaria, but is also present in Echinodermata, Mollusca and Vertebrata. Intriguingly, some species of nematodes and arthropods possess the genes for both pathways. Moreover, expression data from Drosophila suggests that G6Pase and, hence, gluconeogenesis, initially had a neuronal function. We speculate that in insects—and possibly in some vertebrates—gluconeogenesis may be used as a means of neuronal signaling.     

A molecular analysis of the phylogenetic affinities of Saccoglossus cambrensis Brambell & Cole (Hemichordata).

Traditional approaches to phylogeny reconstruction have not allowed precise resolution of the evolutionary relationships between the major deuterostome phyla (chordates, hemichordates, echinoderms). Here we report the use of a molecular approach to investigate deuterostome phylogeny. We have used a polymerase chain reaction-based strategy to amplify, clone and sequence parts of the genes coding for 18S ribosomal RNA from Saccoglossus cambrensis (Hemichordata), Arbacia sp. (Echinodermata) and, for comparison, Mytilus edulis (Mollusca). We report the results of phylogenetic reconstructions using these, and homologous sequences from other eukaryotes. The results of our analyses are consistent with the hypothesis that S. cambrensis and vertebrates share a common ancestor not shared by echinoderms.     

The mitochondrial genome of <Emphasis Type="Italic">Phallusia mammillata</Emphasis> and <Emphasis Type="Italic">Phallusia fumigata</Emphasis> (Tunicata,Ascidiacea): high genome plasticity at intra-genus level

Background  

Within Chordata, the subphyla Vertebrata and Cephalochordata (lancelets) are characterized by a remarkable stability of the mitochondrial (mt) genome, with constancy of gene content and almost invariant gene order, whereas the limited mitochondrial data on the subphylum Tunicata suggest frequent and extensive gene rearrangements, observed also within ascidians of the same genus.     

Fertilization in inhibited in five diverse animal phyla by erythrosin B

Erythrosin B inhibited fertilization in species from the following phyla: Echiuroida, Annelida, Mollusca, Echinodermata, and Chordata. This effect was not due to diminished gamete viability but may reflect a general interference of sperm-egg membrane fusion by Erythrosin B.     

Advanced formulation of base pair changes in the stem regions of ribosomal RNAs; its application to mitochondrial rRNAs for resolving the phylogeny of animals

The ribosomal RNAs (rRNAs) of animal mitochondria, especially those of arthropod mitochondria, have a higher content of G:U and U:G base pairs in their stem regions than the nuclear rRNAs. Thus, the theoretical formulation of base pair changes is extended to incorporate the faster base pair changes A:U<-->G:U<-->G:C and U:A<-->U:G<-->C:G into the previous formulation of the slower base pair changes between A:U, G:C, C:G and U:A. The relative base pair change probability containing the faster and slower base pair changes is theoretically derived to estimate the divergence time of rRNAs under the influence of selection for these base pairs. Using the cartilaginous fish-teleost fish divergence and the crustacean-insect divergence as calibration points, the present method successfully predicts the divergence times of the main branches of animals: Deuterostomia and Protostomia diverged 9.2 x 10(8) years ago, the divergence of Echinodermata, Hemichordata and Cephalochordata succeedingly occurred during the period from 8 x 10(8) to 6 x 10(8) years ago, while Arthropoda, Annelida and Mollusca diverged almost concomitantly about 7 x 10(8) years ago. The dating for the divergence of Platyhelminthes and Cnidaria is traced back to 1.2 x 10(9) years ago. This result is consistent with the fossil records in the Stirling Range Formation of southwestern Australia, the Ediacara and Avalon faunas and the Cambrian Burgess Shale. Thus, the present method may be useful for estimating the divergence times of animals ranging from 10(8) to 10(9) years ago, resolving the difficult problems, e.g. deviation from rate constancy and large sampling variances, in the usual methods of treating apparent change rates between individual bases and/or base pairs.     

Phylogenèse des Phoronida

Phylogenese of Phoronida. Lophophorata and the Archimerata concept The main phylogenetic characteristics of Phoronida and other Lophophorates are discussed: 1. Archimeric segmentation of the body; 2. Egg cleavage of radial (or biradial) type, coeloblastula and gastrulation by invagination (emboly); 3. Mesoderm formation by a derived enteroccelous method (primitive stage of enteroc? ly); 4. Bringing of the anus anteriorly to lie rocoelous method (primitive stage of enterocoely); 4. Bringing of the anus anteriorly to lie sence of a true lophophore; 6. Larva not of trochophoral type, but actinotroch related to Tornaria-Dipleurula; 7. Nervous system basi-epithelial with primitive neurulation in Phoronopsis, without any orthogons; 8. Adult nervous ganglion neo-formed, not issuing from the apicale plate; in Phoronida this ganglion is located in the prosome and not in the mesosome; 9. Lack of cephalization. All these characteristics are closely related to that of the primitive phyla of the Chordata assemblage. The only exception is the presence of protonephridia with solenocytes in actinotroch, but such organs are also described in Cephalochordata. The Lophophorata (Phoronida, Brachiopoda, Bryozoa) are undoubtedly a primitive group on the Chordata trend and the Phoronida appear to be the most phylogenetically evolved phylum of this group with predominating position. The validity of placing the Lophophorata within the Echinoderm-Hemichordata assemblage is demonstrated. The term Archic?lomata appears not suitable and its substitution by Archimerata, assemblage at the base of the Chordata, is here proposed. The Archimerata concept brings together the Lophophorata, Echinodermata and Hemidiordata and is considered as a phylogenetic stage and a natural systematic unit.     

Molecular Phylogeny of Metazoan Intermediate Filament Proteins

We have cloned cytoplasmic intermediate filament (IF) proteins from a large number of invertebrate phyla using cDNA probes, the monoclonal antibody IFA, peptide sequence information, and various RT-PCR procedures. Novel IF protein sequences reported here include the urochordata and nine protostomic phyla, i.e., Annelida, Brachiopoda, Chaetognatha, Echiura, Nematomorpha, Nemertea, Platyhelminthes, Phoronida, and Sipuncula. Taken together with the wealth of data on IF proteins of vertebrates and the results on IF proteins of Cephalochordata, Mollusca, Annelida, and Nematoda, two IF prototypes emerge. The L-type, which includes 35 sequences from 11 protostomic phyla, shares with the nuclear lamins the long version of the coil 1b subdomain and, in most cases, a homology segment of some 120 residues in the carboxyterminal tail domain. The S-type, which includes all four subfamilies (types I to IV) of vertebrate IF proteins, lacks 42 residues in the coil 1b subdomain and the carboxyterminal lamin homology segment. Since IF proteins from all three phyla of the chordates have the 42-residue deletion, this deletion arose in a progenitor prior to the divergence of the chordates into the urochordate, cephalochordate, and vertebrate lineages, possibly already at the origin of the deuterostomic branch. Four phyla recently placed into the protostomia on grounds of their 18S rDNA sequences (Brachiopoda, Nemertea, Phoronida, and Platyhelminthes) show IF proteins of the L-type and fit by sequence identity criteria into the lophotrochozoic branch of the protostomia. Received: 2 April 1998 / Accepted: 19 June 1998     

Isolation of Hox and Parahox genes in the hemichordate Ptychodera flava and the evolution of deuterostome Hox genes

Because of their importance for proper development of the bilaterian embryo, Hox genes have taken center stage for investigations into the evolution of bilaterian metazoans. Taxonomic surveys of major protostome taxa have shown that Hox genes are also excellent phylogenetic markers, as specific Hox genes are restricted to one of the two great protostome clades, the Lophotrochozoa or the Ecdysozoa, and thus support the phylogenetic relationships as originally deduced by 18S rDNA studies. Deuterostomes are the third major group of bilaterians and consist of three major phyla, the echinoderms, the hemichordates, and the chordates. Most morphological studies have supported Hemichordata+Chordata, whereas molecular studies support Echinodermata+Hemichordata, a clade known as Ambulacraria. To test these competing hypotheses, complete or near complete cDNAs of eight Hox genes and four Parahox genes were isolated from the enteropneust hemichordate Ptychodera flava. Only one copy of each Hox gene was isolated suggesting that the Hox genes of P. flava are arranged in a single cluster. Of particular importance is the isolation of three posterior or Abd-B Hox genes; these genes are only shared with echinoderms, and thus support the monophyly of Ambulacraria.     

Triploblastic relationships with emphasis on the acoelomates and the position of Gnathostomulida,Cycliophora, Plathelminthes,and Chaetognatha: a combined approach of 18S rDNA sequences and morphology

Triploblastic relationships were examined in the light of molecular and morphological evidence. Representatives for all triploblastic "phyla" (except Loricifera) were represented by both sources of phylogenetic data. The 18S ribosomal (rDNA) sequence data for 145 terminal taxa and 276 morphological characters coded for 36 supraspecific taxa were combined in a total evidence regime to determine the most consistent picture of triploblastic relationships for these data. Only triploblastic taxa are used to avoid rooting with distant outgroups, which seems to happen because of the extreme distance that separates diploblastic from triploblastic taxa according to the 18S rDNA data. Multiple phylogenetic analyses performed with variable analysis parameters yield largely inconsistent results for certain groups such as Chaetognatha, Acoela, and Nemertodermatida. A normalized incongruence length metric is used to assay the relative merit of the multiple analyses. The combined analysis having the least character incongruence yields the following scheme of relationships of four main clades: (1) Deuterostomia [((Echinodermata + Enteropneusta) (Cephalochordata (Urochordata + Vertebrata)))]; (2) Ecdysozoa [(((Priapulida + Kinorhyncha) (Nematoda + Nematomorpha)) ((Onychophora + Tardigrada) Arthropoda))]; (3) Trochozoa [((Phoronida + Brachiopoda) (Entoprocta (Nemertea (Sipuncula (Mollusca (Pogonophora (Echiura + Annelida)))))))]; and (4) Platyzoa [((Gnathostomulida (Cycliophora + Syndermata)) (Gastrotricha + Plathelminthes))]. Chaetognatha, Nemertodermatida, and Bryozoa cannot be assigned to any one of these four groups. For the first time, a data analysis recognizes a clade of acoelomates, the Platyzoa (sensu Cavalier-Smith, Biol. Rev. 73:203-266, 1998). Other relationships that corroborate some morphological analyses are the existence of a clade that groups Gnathostomulida + Syndermata (= Gnathifera), which is expanded to include the enigmatic phylum Cycliophora, as sister group to Syndermata.