Animal phylogeny and the ancestry of bilaterians: inferences from morphology and 18S rDNA gene sequences

SUMMARY Insight into the origin and early evolution of the animal phyla requires an understanding of how animal groups are related to one another. Thus, we set out to explore animal phylogeny by analyzing with maximum parsimony 138 morphological characters from 40 metazoan groups, and 304 18S rDNA sequences, both separately and together. Both types of data agree that arthropods are not closely related to annelids: the former group with nematodes and other molting animals (Ecdysozoa), and the latter group with molluscs and other taxa with spiral cleavage. Furthermore, neither brachiopods nor chaetognaths group with deuterostomes; brachiopods are allied with the molluscs and annelids (Lophotrochozoa), whereas chaetognaths are allied with the ecdysozoans. The major discordance between the two types of data concerns the rooting of the bilaterians, and the bilaterian sister-taxon. Morphology suggests that the root is between deuterostomes and protostomes, with ctenophores the bilaterian sister-group, whereas 18S rDNA suggests that the root is within the Lophotrochozoa with acoel flatworms and gnathostomulids as basal bilaterians, and with cnidarians the bilaterian sister-group. We suggest that this basal position of acoels and gnathostomulids is artifactal because for 1000 replicate phylogenetic analyses with one random sequence as outgroup, the majority root with an acoel flatworm or gnathostomulid as the basal ingroup lineage. When these problematic taxa are eliminated from the matrix, the combined analysis suggests that the root lies between the deuterostomes and protostomes, and Ctenophora is the bilaterian sister-group. We suggest that because chaetognaths and lophophorates, taxa traditionally allied with deuterostomes, occupy basal positions within their respective protostomian clades, deuterostomy most likely represents a suite of characters plesiomorphic for bilaterians.     

Phylogeny of protostome moulting animals (Ecdysozoa) inferred from 18 and 28S rRNA gene sequences

Reliability of reconstruction of phylogenetic relationships within a group of protostome moulting animals was evaluated by means of comparison of 18 and 28S rRNA gene sequences sets both taken separately and combined. Reliability of reconstructions was evaluated by values of the bootstrap support of major phylogenetic tree nodes and by degree of congruence of phylogenetic trees inferred by various methods. By both criteria, phylogenetic trees reconstructed from the combined 18 and 28S rRNA gene sequences were better than those inferred from 18 and 28S sequences taken separately. Results obtained are consistent with phylogenetic hypothesis separating protostome animals into two major clades, moulting Ecdysozoa (Priapulida + Kinorhyncha, Nematoda + Nematomorpha, Onychophora + Tardigrada, Myriapoda + Chelicerata, Crustacea + Hexapoda) and unmoulting Lophotrochozoa (Plathelminthes, Nemertini, Annelida, Mollusca, Echiura, Sipuncula). Clade Cephalorhyncha does not include nematomorphs (Nematomorpha). Conclusion was taken that it is necessary to use combined 18 and 28S data in phylogenetic studies.     

Phylogeny of Molting Protostomes (Ecdysozoa) as Inferred from 18S and 28S rRNA Gene Sequences

Phylogenetic relationships within the group of molting protostomes were reconstructed by comparing the sets of 18S and 28S rRNA gene sequences considered either separately or in combination. The reliability of reconstructions was estimated from the bootstrap indices for major phylogenetic tree nodes and from the degree of congruence of phylogenetic trees obtained by different methods. By either criterion, the phylogenetic trees reconstructed on the basis of both 18 and 28S rRNA gene sequences were better than those based on the 18S or 28S sequences alone. The results of reconstruction are consistent with the phylogenetic hypothesis classifying protostomes into two major clades: molting Ecdysozoa (Priapulida + Kinorhyncha, Nematoda + Nematomorpha, Onychophora + Tardigrada, Myriapoda + Chelicerata, and Crustacea + Hexapoda) and nonmolting Lophotrochozoa (Plathelminthes, Nemertini, Annelida, Mollusca, Echiura, and Sipuncula). Nematomorphs (Nematomorpha) do not belong to the clade Cephalorhyncha (Priapulida + Kinorhyncha). It is concluded that combined data on the 18S and 28S rRNA gene sequences provide a more reliable basis for phylogenetic inferences.__________Translated from Molekulyarnaya Biologiya, Vol. 39, No. 4, 2005, pp. 590–601.Original Russian Text Copyright © 2005 by Petrov, Vladychenskaya.     

Evaluating hypotheses of deuterostome phylogeny and chordate evolution with new LSU and SSU ribosomal DNA data

We investigated evolutionary relationships among deuterostome subgroups by obtaining nearly complete large-subunit ribosomal RNA (LSU rRNA)-gene sequences for 14 deuterostomes and 3 protostomes and complete small-subunit (SSU) rRNA-gene sequences for five of these animals. With the addition of previously published sequences, we compared 28 taxa using three different data sets (LSU only, SSU only, and combined LSU + SSU) under minimum evolution (with LogDet distances), maximum likelihood, and maximum parsimony optimality criteria. Additionally, we analyzed the combined LSU + SSU sequences with spectral analysis of LogDet distances, a technique that measures the amount of support and conflict within the data for every possible grouping of taxa. Overall, we found that (1) the LSU genes produced a tree very similar to the SSU gene tree, (2) adding LSU to SSU sequences strengthened the bootstrap support for many groups above the SSU-only values (e.g., hemichordates plus echinoderms as Ambulacraria; lancelets as the sister group to vertebrates), (3) LSU sequences did not support SSU-based hypotheses of pterobranchs evolving from enteropneusts and thaliaceans evolving from ascidians, and (4) the combined LSU + SSU data are ambiguous about the monophyly of chordates. No tree-building algorithm united urochordates conclusively with other chordates, although spectral analysis did so, providing our only evidence for chordate monophyly. With spectral analysis, we also evaluated several major hypotheses of deuterostome phylogeny that were constructed from morphological, embryological, and paleontological evidence. Our rRNA-gene analysis refutes most of these hypotheses and thus advocates a rethinking of chordate and vertebrate origins.     

Phylogenomic analyses of lophophorates (brachiopods, phoronids and bryozoans) confirm the Lophotrochozoa concept

Based on embryological and morphological evidence, Lophophorata was long considered to be the sister or paraphyletic stem group of Deuterostomia. By contrast, molecular data have consistently indicated that the three lophophorate lineages, Ectoprocta, Brachiopoda and Phoronida, are more closely related to trochozoans (annelids, molluscs and related groups) than to deuterostomes. For this reason, the lophophorate groups and Trochozoa were united to Lophotrochozoa. However, the relationships of the lophophorate lineages within Lophotrochozoa are still largely unresolved. Maximum-likelihood and Bayesian analyses were performed based on a dataset comprising 11,445 amino acid positions derived from 79 ribosomal proteins of 39 metazoan taxa including new sequences obtained from a brachiopod and a phoronid. These analyses show that the three lophophorate lineages are affiliated with trochozoan rather than deuterostome phyla. All hypotheses claiming that they are more closely related to Deuterostomia than to Protostomia can be rejected by topology testing. Monophyly of lophophorates was not recovered but that of Bryozoa including Ectoprocta and Entoprocta and monophyly of Brachiozoa including Brachiopoda and Phoronida were strongly supported. Alternative hypotheses that are refuted include (i) Brachiozoa as the sister group of Mollusca, (ii) ectoprocts as sister to all other Lophotrochozoa including Platyzoa, and (iii) ectoprocts as sister or to all other protostomes except chaetognaths.     

Investigation of molluscan phylogeny using large-subunit and small-subunit nuclear rRNA sequences

The Mollusca represent one of the most morphologically diverse animal phyla, prompting a variety of hypotheses on relationships between the major lineages within the phylum based upon morphological, developmental, and paleontological data. Analyses of small-ribosomal RNA (SSU rRNA) gene sequence have provided limited resolution of higher-level relationships within the Mollusca. Recent analyses suggest large-subunit (LSU) rRNA gene sequences are useful in resolving deep-level metazoan relationships, particularly when combined with SSU sequence. To this end, LSU (approximately 3.5 kb in length) and SSU (approximately 2 kb) sequences were collected for 33 taxa representing the major lineages within the Mollusca to improve resolution of intraphyletic relationships. Although the LSU and combined LSU+SSU datasets appear to hold potential for resolving branching order within the recognized molluscan classes, low bootstrap support was found for relationships between the major lineages within the Mollusca. LSU+SSU sequences also showed significant levels of rate heterogeneity between molluscan lineages. The Polyplacophora, Gastropoda, and Cephalopoda were each recovered as monophyletic clades with the LSU+SSU dataset. While the Bivalvia were not recovered as monophyletic clade in analyses of the SSU, LSU, or LSU+SSU, the Shimodaira-Hasegawa test showed that likelihood scores for these results did not differ significantly from topologies where the Bivalvia were monophyletic. Analyses of LSU sequences strongly contradict the widely accepted Diasoma hypotheses that bivalves and scaphopods are closely related to one another. The data are consistent with recent morphological and SSU analyses suggesting scaphopods are more closely related to gastropods and cephalopods than to bivalves. The dataset also presents the first published DNA sequences from a neomeniomorph aplacophoran, a group considered critical to our understanding of the origin and early radiation of the Mollusca.     

Phylogenetic analysis of the Metastrongyloidea (Nematoda: Strongylida) inferred from ribosomal RNA gene sequences

Phylogenetic relationships among nematodes of the strongylid superfamily Metastrongyloidea were analyzed using partial sequences from the large-subunit ribosomal RNA (LSU rRNA) and small-subunit ribosomal RNA (SSU rRNA) genes. Regions of nuclear ribosomal DNA (rDNA) were amplified by polymerase chain reaction, directly sequenced, aligned, and phylogenies inferred using maximum parsimony. Phylogenetic hypotheses inferred from the SSU rRNA gene supported the monophyly of representative taxa from each of the 7 currently accepted metastrongyloid families. Metastrongyloid taxa formed the sister group to representative trichostrongyloid sequences based on SSU data. Sequences from either the SSU or LSU RNA regions alone provided poor resolution for relationships within the Metastrongyloidea. However, a combined analysis using sequences from all rDNA regions yielded 3 equally parsimonious trees that represented the abursate Filaroididae as polyphyletic, Parafilaroides decorus as the sister species to the monophyletic Pseudaliidae, and a sister group relationship between Oslerus osleri and Metastrongylus salmi. Relationships among 3 members of the Crenosomatidae, and 1 representative of the Skrjabingylidae (Skrjabingylus chitwoodorum) were not resolved by these combined data. However, members of both these groups were consistently resolved as the sister group to the other metastrongyloid families. These relationships are inconsistent with traditional classifications of the Metastrongyloidea and existing hypotheses for their evolution.     

Lophotrochozoan phylogeny assessed with LSU and SSU data: evidence of lophophorate polyphyly

Of the three major bilaterian clades, Lophotrochozoa has the greatest diversity and disparity of body forms and is the least understood in terms of phylogenetic history. Within this clade, small nuclear ribosomal subunit (SSU or 18S) studies have failed to provide resolution and other molecular markers have insufficient taxon sampling. To examine relationships within Lophotrochozoa, we collected and complied complete SSU data and nearly complete (>90%) large nuclear ribosomal subunit (LSU or 28S) data totaling approximately 5kb per taxon, for 36 lophotrochozoans. Results of LSU and combined SSU+LSU likelihood analyses provide topologies more consistent with morphological data than analyses of SSU data alone. Namely, most phyla recognized on morphological grounds are recovered as monophyletic entities when the LSU data is considered (contra SSU data alone). These new data show with significant support that "Lophophorata" (traditionally recognized to include Brachiopoda, Phoronida, and Bryozoa) is not a monophyletic entity. Further, the data suggest that Platyzoa is real and may be derived within lophotrochozans rather than a basal or sister taxon. The recently discovered Cycliophora are allied to entoprocts, consistent with their initial placement based on morphology. Additional evidence for Syndermata (i.e., Rotifera+Acanthocephala) is also found. Although relationships among groups with trochophore-like larvae could not be resolved and nodal support values are generally low, the addition of LSU data is a considerable advance in our understanding of lophotrochozoan phylogeny from the molecular perspective.     

Complete sequence and structure of ribosomal RNA gene of Heterosporis anguillarum

The ribosomal RNA (rRNA) gene region of the microsporidium Heterosporis anguillarum has been examined. Complete DNA sequence data (4060 bp, GenBank Accession No. AF402839) of the rRNA gene of H. anguillarum are presented for the small subunit gene (SSU rRNA: 1359 bp), the internal transcribed spacer (ITS: 37 bp), and the large subunit gene (LSU rRNA: 2664 bp). The secondary structures of the H. anguillarum SSU and LSU rRNA genes are constructed and described. This is the first complete sequence of an rRNA gene published for a fish-infecting microsporidian species. In the phylogenetic analysis, the sequences, including partial SSU rRNA, ITS, and partial LSU rRNA sequences of the fish-infecting microsporidia, were aligned and analysed. The taxonomic position of H. anguillarum as suggested by Lom et al. (2000; Dis Aquat Org 43:225-231) is confirmed in this paper.     

Phylogenetic analysis of ribosomal RNA operons from uncultivated coastal marine bacterioplankton

Analyses of small subunit ribosomal RNA genes (SSU rDNAs) have significantly influenced our understanding of the composition of aquatic microbial assemblages. Unfortunately, SSU rDNA sequences often do not have sufficient resolving power to differentiate closely related species. To address this general problem for uncultivated bacterioplankton taxa, we analysed and compared sequences of polymerase chain reaction (PCR)-generated and bacterial artificial chromosome (BAC)-derived clones that contained most of the SSU rDNAs, the internal transcribed spacer (ITS) and the large subunit ribosomal RNA gene (LSU rDNA). The phylogenetic representation in the rRNA operon PCR library was similar to that reported previously in coastal bacterioplankton SSU rDNA libraries. We observed good concordance between the phylogenetic relationships among coastal bacterioplankton inferred from SSU or LSU rDNA sequences. ITS sequences confirmed the close intragroup relationships among members of the SAR11, SAR116 and SAR86 clades that were predicted by SSU and LSU rDNA sequence analyses. We also found strong support for homologous recombination between the ITS regions of operons from the SAR11 clade.     

Analysis of the complete mitochondrial DNA sequence of the brachiopod terebratulina retusa places Brachiopoda within the protostomes

Brachiopod phylogeny is still a controversial subject. Analyses using nuclear 18SrRNA and mitochondrial 12SrDNA sequences place them within the protostomes but some recent interpretations of morphological data support a relationship with deuterostomes. In order to investigate brachiopod affinities within the metazoa further, we compared the gene arrangement on the brachiopod mitochondrial genome with several metazoan taxa. The complete (15 451 bp) mitochondrial DNA (mtDNA) sequence of the articulate brachiopod Terebratulina retusa was determined from two overlapping long polymerase chain reaction products. All the genes are encoded on the same strand and gene order comparisons showed that.only one major rearrangement is required to interconvert the T. retusa and Katharina tunicata (Mollusca: Polvplacophora) mitochondrial genomes. The partial mtDNA sequence of the prosobranch mollusc Littorina saxatilis shows complete congruence with the T. rehtusa gene arrangement with regard to the ribosomal and protein coding genes. This high similarity in gene arrangement is the first to be reported within the protostomes. Sequence analyses of mitochondrial protein coding genes also support a close relationship of the brachiopod with molluscs and annelids, thus supporting the clade Lophotrochozoa. Though being highly informative, sequence analyses of the mitochondrial protein coding genes failed to resolve the branching order within the lophotrochozoa.     

Testing the new animal phylogeny: a phylum level molecular analysis of the animal kingdom

The new animal phylogeny inferred from ribosomal genes some years ago has prompted a number of radical rearrangements of the traditional, morphology based metazoan tree. The two main bilaterian clades, Deuterostomia and Protostomia, find strong support, but the protostomes consist of two sister groups, Ecdysozoa and Lophotrochozoa, not seen in morphology based trees. Although widely accepted, not all recent molecular phylogenetic analyses have supported the tripartite structure of the new animal phylogeny. Furthermore, even if the small ribosomal subunit (SSU) based phylogeny is correct, there is a frustrating lack of resolution of relationships between the phyla that make up the three clades of this tree. To address this issue, we have assembled a dataset including a large number of aligned sequence positions as well as a broad sampling of metazoan phyla. Our dataset consists of sequence data from ribosomal and mitochondrial genes combined with new data from protein coding genes (5139 amino acid and 3524 nucleotide positions in total) from 37 representative taxa sampled across the Metazoa. Our data show strong support for the basic structure of the new animal phylogeny as well as for the Mandibulata including Myriapoda. We also provide some resolution within the Lophotrochozoa, where we confirm support for a monophyletic clade of Echiura, Sipuncula and Annelida and surprising evidence of a close relationship between Brachiopoda and Nemertea.     

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- 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.     

A multigene phylogeny of the Gelidiales including nuclear large-subunit rRNA sequence data

Partial sequences (1032 bp) of the nuclear-encoded large ribosomal RNA gene (LSU) were determined for 16 gelidialean species, and analyzed separately and in combination with plastid rbcL and nuclear SSU gene sequences. The number of informative characters and levels of sequence divergence among taxa are intermediate in LSU sequences as compared to that for rbcL and SSU. Analyses of the separate LSU, and a combined LSU, SSU, and rbcL data sets have identified early-diverging lineages within the Gelidiales including Gelidiella, Pterocladia, Pterocladiella, and a lineage including Gelidium and species classified in other genera. The relationships among most gelidialean taxa are well-resolved and well-supported by analyses of the combined data; however, the relationships of Ptilophora and Capreolia remain unclear. It is speculated that these two lineages have diverged from a common ancestor over an evolutionarily short period of time. This revised version was published online in June 2006 with corrections to the Cover Date.     

The mysterious evolutionary origin for the GNE gene and the root of bilateria

Phylogenomic analyses have revealed several important metazoan clades, such as the Ecdysozoa and the Lophotrochozoa. However, the phylogenetic positions of a few taxa, such as ctenophores, chaetognaths, acoelomorphs, and Xenoturbella, remain contentious. Thus, the findings of qualitative markers or "rare genomic changes" seem ideal to independently test previous phylogenetic hypotheses. We here describe a rare genomic change, the presence of the gene UDP-GlcNAc 2-epimerase/N-acetylmannosamine kinase (GNE). We show that GNE is encoded in the genomes of deuterostomes, acoelomorphs and Xenoturbella, whereas it is absent in protostomes and nonbilaterians. Moreover, the GNE has a complex evolutionary origin involving unique lateral gene transfer events and/or extensive hidden paralogy for each protein domain. However, rather than using GNE as a phylogenetic character, we argue that rare genomic changes such as the one presented here should be used with caution.     

Concatenated SSU and LSU rDNA data confirm the main evolutionary trends within myxosporeans (Myxozoa: Myxosporea) and provide an effective tool for their molecular phylogenetics

Views on myxosporean phylogeny and systematics have recently undergone substantial changes resulting from analyses of SSU rDNA. Here, we further investigate the evolutionary trends within myxosporean lineages by using 35 new sequences of the LSU rDNA. We show a good agreement between the two rRNA genes and confirm the main phylogenetic split between the freshwater and marine lineages. The informative superiority of the LSU data is shown by an increase of the resolution, nodal supports and tree indexes in the LSU rDNA and combined analyses. We determine the most suitable part of LSU for the myxosporean phylogeny by comparing informative content in various regions of the LSU sequences. Based on this comparison, we propose the D5–3′-end part of the LSU rRNA gene as the most informative region which provides in concatenation with the complete SSU a well resolved and robust tree. To allow for simple amplification of the marker, we design specific primer set for this part of LSU rDNA.     

Comparative biochemical study of N-linked glycans from skin of a squid, Todarodes pacificus

For comparative biochemical interest, we analyzed the structures of N-glycans in a squid belonging to the Lophotrochozoa, one of the protostome clades. N-Glycans were prepared from squid skin by hydrazinolysis and re-N-acetylation followed by fluorescent tagging with 2-aminopyridine. The labeled N-glycans were purified, and their structures were determined by the two-dimensional HPLC mapping method combined with glycosidase digestions and mass spectrometry. We found that high mannose-type glycans, paucimannose-type glycans and complex-type glycans with a type-1 structure (Galbeta1-3GlcNAc) were dominant in squid skin. The complex-type glycans detected in the squid were similar to those in vertebrates, but have not yet been found in the Ecdysozoa, which is another protostome clade. However, paucimannose-type glycans are commonly found in the Ecdysozoa. Thus, the N-glycan structures of the squid belonging to the Lophotrochozoa have features common to those in vertebrates and the Ecdysozoa including insects and nematodes.     

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.     

Molecular diversity of the syndinean genus Euduboscquella based on single-cell PCR analysis

The genus Euduboscquella is one of a few described genera within the syndinean dinoflagellates, an enigmatic lineage with abundant diversity in marine environmental clone libraries based on small subunit (SSU) rRNA. The region composed of the SSU through to the partial large subunit (LSU) rRNA was determined from 40 individual tintinnid ciliate loricae infected with Euduboscquella sampled from eight surface water sites in the Northern Hemisphere, producing seven distinct SSU sequences. The corresponding host SSU rRNA region was also amplified from eight host species. The SSU tree of Euduboscquella and syndinean group I sequences from environmental clones had seven well-supported clades and one poorly supported clade across data sets from 57 to 692 total sequences. The genus Euduboscquella consistently formed a supported monophyletic clade within a single subclade of group I sequences. For most parasites with identical SSU sequences, the more variable internal transcribed spacer (ITS) to LSU rRNA regions were polymorphic at 3 to 10 sites. However, in E. cachoni there was variation between ITS to LSU copies at up to 20 sites within an individual, while in a parasite of Tintinnopsis spp., variation between different individuals ranged up to 19 polymorphic sites. However, applying the compensatory base change model to the ITS2 sequences suggested no compensatory changes within or between individuals with the same SSU sequence, while one to four compensatory changes between individuals with similar but not identical SSU sequences were found. Comparisons between host and parasite phylogenies do not suggest a simple pattern of host or parasite specificity.