Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Research article

EST based phylogenomics of Syndermata questions monophyly of Eurotatoria

Alexander Witek1, Holger Herlyn2, Achim Meyer3, Louis Boell4, Gregor Bucher4 and Thomas Hankeln1*

Author Affiliations

1 Institute of Molecular Genetics, Johannes Gutenberg-University Mainz, J. J.-Becherweg 32, D-55099 Mainz, Germany

2 Institute of Anthropology, Johannes Gutenberg-University Mainz, Colonel-Kleinmann-Weg 2, D-55099 Mainz, Germany

3 Institute of Zoology, Johannes Gutenberg-University Mainz, Müllerweg 6, D-55099 Mainz, Germany

4 Johann Friedrich Blumenbach Institute of Anthropology and Zoology, Georg-August-University Göttingen, J. v. Liebig-Weg 11, D-37077 Göttingen, Germany

For all author emails, please log on.

BMC Evolutionary Biology 2008, 8:345  doi:10.1186/1471-2148-8-345

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/8/345


Received:16 June 2008
Accepted:29 December 2008
Published:29 December 2008

© 2008 Witek et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The metazoan taxon Syndermata comprising Rotifera (in the classical sense of Monogononta+Bdelloidea+Seisonidea) and Acanthocephala has raised several hypotheses connected to the phylogeny of these animal groups and the included subtaxa. While the monophyletic origin of Syndermata and Acanthocephala is well established based on morphological and molecular data, the phylogenetic position of Syndermata within Spiralia, the monophyletic origin of Monogononta, Bdelloidea, and Seisonidea and the acanthocephalan sister group are still a matter of debate. The comparison of the alternative hypotheses suggests that testing the phylogenetic validity of Eurotatoria (Monogononta+Bdelloidea) is the key to unravel the phylogenetic relations within Syndermata. The syndermatan phylogeny in turn is a prerequisite for reconstructing the evolution of the acanthocephalan endoparasitism.

Results

Here we present our results from a phylogenomic approach studying i) the phylogenetic position of Syndermata within Spiralia, ii) the monophyletic origin of monogononts and bdelloids and iii) the phylogenetic relations of the latter two taxa to acanthocephalans. For this analysis we have generated EST libraries of Pomphorhynchus laevis, Echinorhynchus truttae (Acanthocephala) and Brachionus plicatilis (Monogononta). By extending these data with database entries of B. plicatilis, Philodina roseola (Bdelloidea) and 25 additional metazoan species, we conducted phylogenetic reconstructions based on 79 ribosomal proteins using maximum likelihood and bayesian approaches. Our findings suggest that the phylogenetic position of Syndermata within Spiralia is close to Platyhelminthes, that Eurotatoria are not monophyletic and that bdelloids are more closely related to acanthocephalans than monogononts.

Conclusion

Mapping morphological character evolution onto molecular phylogeny suggests the (partial or complete) reduction of the corona and the emergence of a retractable anterior end (rostrum, proboscis) before the separation of Acanthocephala. In particular, the evolution of a rostrum might have been a key event leading to the later evolution of the acanthocephalan endoparasitism, given the enormous relevance of the proboscis for anchoring of the adults to the definitive hosts' intestinal wall.

Background

The animal taxon Rotifera comprises free-living and commensalic microscopic species of aquatic habitats that are traditionally grouped into the three subtaxa Bdelloidea, Monogononta and Seisonidea [1-3]. Bdelloids (about 460 species) inhabit freshwater, are capable of anhydrobiosis and reproduce strictly by parthenogenesis. Monogononts (about 1,570 species) live in limnic, brackish and marine waters and have a lifecycle with alternating phases of parthenogenetic and sexual reproduction. Thirdly, at least two of the hitherto three described species belonging to Seisonidea are epibionts on marine crustaceans of the genus Nebalia [3,4]. Though Bdelloidea, Monogononta and Seisonidea are subsumed as Rotifera or Rotatoria, the eponymous rotatory organ or corona, a seemingly rotating assembly of cilia at the anterior end of the animal, is absent in Seisonidea. For this and other reasons Bdelloidea and Monogononta are often regarded as sistergroups of a taxon named Eurotatoria [2,5,6]. In contrast to Bdelloidea, Monogononta and Seisonidea, the Acanthocephala are obligatory endoparasites with a complicated lifecycle. Their definite hosts are vertebrates, while their intermediate hosts are insects, chilopods and crustaceans (e.g., Meyer [7]). Along with the endoparasitic life cycle, the acanthocephalan subtaxa share a plethora of derived morphological characters (e.g., [8-11]) so that the monophyly of Acanthocephala as a whole has never been debated. Moreover, the grouping of Acanthocephala, Bdelloidea, Monogononta and Seisonidea into the taxon Syndermata is widely accepted due to special features in epidermal and sperm ultrastructure (e.g., syncytial epidermis, spermatozoon with anteriorly inserted cilium; see [5,8,9,12,13]), as well as congruent results from molecular approaches [14-20]. It is further undisputed that Syndermata are part of a more comprehensive monophylum called Gnathifera [9,21,22]. On the other hand, the phylogenetic position of Syndermata beyond Gnathifera as well as the relationships among the syndermatan subtaxa Acanthocephala, Bdelloidea, Monogononta and Seisonidea are still unresolved. So far, five competing hypotheses on the internal phylogeny of Syndermata have been suggested (Fig. 1A–E). The Lemniscea hypothesis goes back to Lorenzen [23] and favors a sister group relationship of bdelloids and acanthocephalans, with the Monogononta and Seisonidea placed basally to the Lemniscea (Fig. 1A). Morphological evidence for such grouping has been inferred from two lateral intrusions in the neck region and a retractable anterior body section in Acanthocephala and Bdelloidea [23]. The Lemniscea hypothesis received additional support from 16S rRNA, 18S rRNA, 28S rRNA, cytochrome c oxidase subunit 1 (cox 1) and histone H3 data [14-17,19]. The second hypothesis suggests a sistergroup relationship of Monogononta and Bdelloidea (Eurotatoria) and of Seisonidea and Acanthocephala (Pararotatoria) and is herein called Eurotatoria+Pararotatoria hypothesis (Fig. 1B). Besides presumed eurotatorian apomorphies such as the already mentioned corona, the Eurotatoria+Pararotatoria hypothesis is based on ultrastructural peculiarities that have been interpreted as synapomorphies of Seisonidea and Acanthocepahala (spermatozoa with "dense bodies" and epidermis with special filaments [5,9,12]). Additional support for the monophyly of Pararotatoria came from partial 18S rRNA data [24] as well as from a combined dataset of 18S rRNA sequences, heat shock gene sequences (hsp82), and morphological characters [25]. The third hypothesis reflects the classical view of monophyletic Rotifera (Monogononta+Bdelloidea+Seisonidea) and Eurotatoria (Monogononta+Bdelloidea) and proposes Acanthocephala as the sistergroup of Rotifera ("classical Rotifera+Acanthocephala hypothesis", see Fig. 1C). This classical concept has been formulated based on specific features of toe morphology, sensory and masticatory apparatus in Rotifera and Eurotatoria, respectively [13,26], and was supported by 18S rRNA data [27]. The fourth hypothesis has been proposed on the basis of hsp82 sequences, and groups Acanthocephala and Eurotatoria with exclusion of Seisonidea [28] ("Eurotatoria+Acanthocephala", see Fig. 1D). Underlying the fourth hypothesis, the absence of acrosomal structures might represent a synapomorphy of Eurotatoria and Acanthocephala [21]. According to the fifth hypothesis, Bdelloidea, Seisonidea, and Acanthocephala form a monophylum for which the name Hemirotifera has been proposed (Fig. 1E). The Hemirotifera hypothesis has been inferred from a combined dataset of molecular (18S rRNA, 28S rRNA, histone H3, cox 1) and morphological characters [29]. This survey of competing hypotheses demonstrates that the question of phylogenetic relationships within Syndermata and therewith of the evolution of the acanthocephalan endoparasitism is closely connected to the more basal question of monophyly of Eurotatoria.

thumbnailFigure 1. Competing phylogenetic hypotheses amongst Syndermata. Cladograms reflecting the competing hypotheses on the phylogenetic relations among Monogononta, Bdelloidea, Acanthocephala and Seisonidea. A Lemniscea hypothesis [23]. B Eurotatoria+Pararotatoria hypothesis [5,9,12]. C Rotifera+Acanthocephala [8,26]. D Eurotatoria+Acanthocephala [28]. E Hemirotifera [29].

In the present study we analyse the phylogenetic position of Syndermata within the spiralian clade as well as the phylogenetic relations among the syndermatan subtaxa. We particularly focus on the question whether Eurotatoria are monophyletic and – if not – whether bdelloids or monogononts are more closely related to acanthocephalans. As ribosomal proteins are favorable tools for metazoan molecular phylogenetic analyses [18,30-32] and easy to obtain from EST libraries, we compiled a phylogenomic dataset comprising 79 ribosomal proteins. To this end, we generated EST libraries for one monogonont (B. plicatilis) and two acanthocephalans (P. laevis and E. truttae; both Echinorhynchida) and sequenced 1,000–2,000 ESTs per library. The new sequences were complemented with ortholog data from public databases for the monogonont B. plicatilis, the bdelloid P. roseola and 25 additional metazoan taxa. Data of Seisonidea have not been included in the present analysis as it is extremely difficult to obtain sufficient material for the preparation of a cDNA library. As a beneficial side effect, the present tree reconstruction cannot be disturbed by the observed long branch leading to representatives of Seisonidea (see, e.g., [19,24]).

Results

Sequence analyses and ribosomal protein alignment

EST sequencing was performed for three syndermatan species and complemented by sequences from public databases (Tab. 1). A dataset containing the coding sequences of 79 ribosomal proteins was extracted, and derived amino acid sequences were concatenated. After cleaning the raw data from ambiguously aligned positions, the final alignment had a length of 11,276 amino acids. Fifteen to 29 species were sampled per protein type, resulting in a taxon coverage ranging from 36 to 100% outside Syndermata (see Additional files 1 and 2 for the complete matrix of taxa and ribosomal proteins used and the amino acid coverage of each ribosomal protein), and 28 to 89% within Syndermata, compared to the vertebrate reference. Given the length of the complete dataset the minimum coverage within the syndermatan sample (see E. truttae in table 2) still represents 3,204 amino acid positions which is a considerable increase in data compared to previous analyses of syndermatan phylogeny (e.g., [17,24]). The total coverage of Acanthocephala, however, is much higher than suggested by the E. truttae data alone, due to the additional 66% sequence coverage in P. laevis (table 2).

Additional file 1. Supplementary table 1. List of taxa, number of ribosomal proteins, number of amino acids and the percental coverage per taxon used in our analyses.

Format: DOC Size: 51KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 2. Supplementary table 2. Complete matrix of taxa and ribosomal proteins used in this analysis with the number of amino acids per ribosomal protein. Maximum length of each ribosomal protein is shown in brackets under each protein name.

Format: DOC Size: 534KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Table 1. List of the syndermatan species for which new data have been collected in the present analysis

Table 2. Syndermatan coverage in the dataset

Likelihood mapping analysis determined a strong phylogenetic signal in the data. In detail, 99.1% of the quartets were fully resolved and none of the quartet-trees showed a star-like topology (Fig. 2). Furthermore, we found no evidence for horizontal gene transfer in the bdelloid dataset, applying the test statistics proposed by Gladyshev et al. [33]: As to be expected for ribosomal genes, the so-called Alien Index was < 0 for each protein. We therefore consider our dataset a sound basis for assessing the phylogenetic position of Syndermata within Spiralia, and for answering the question of eurotatorian monophyly.

thumbnailFigure 2. Likelihood mapping of the concatenated alignment. Results from likelihood mapping of the concatenated ribosomal protein sequences from 29 species analysed in the present study. Note that 99.1% of all quartets were fully resolved and none of the quartets produced a star-like tree.

Phylogenetic reconstruction

Maximum likelihood and bayesian phylogenetic inference consistently support a monophyletic origin of Spiralia, although with partly moderate support values (PhyML < 50; Treefinder: 91; PhyloBayes: 0.93) (Fig. 3, 4, 5). Interestingly, only the PhyloBayes tree depicts the widely accepted monophylum Ecdysozoa (e.g., [16]) so that bayesian inference might provide more reliable results than the maximum likelihood approaches employed, given the present data (Fig. 3). Irrespective of this detail, all three tree reconstruction methods yield maximum support for a monophyletic origin of B. plicatilis, P. roseola, E. truttae and P. laevis and, thus, for the monophyletic origin of the four syndermatan species covered by the present dataset (PhyML: 100; Treefinder: 100, PhyloBayes: 1.00). Within the spiralian clade, Syndermata either group with Platyhelminthes (PhyML: 59; Treefinder: 92) (Fig. 4 and 5) or with a clade comprising Platyhelminthes, Bryozoa, Mollusca and Annelida (PhyloBayes: 0.93) (Fig. 3). Remarkably, none of the tree reconstruction methods supports a sister group relationship of Monogononta and Bdelloidea. Instead, Bdelloidea consistently appear more closely related to Acanthocephala than to Monogononta (PhyML: 78; Treefinder: 76; PhyloBayes: 0.83). The paraphyly of Eurotatoria is further corroborated by results from hypothesis testing. Thereafter, a grouping of Bdelloidea+Acanthocephala is much more likely than the alternatives Monogononta+Bdelloidea and Monogononta+Acanthocephala (Tab. 3). Final evidence for the robustness of the present analysis comes from testing for the effect of missing data in the full-length dataset on the results of tree reconstruction [41]. Thus, the internal syndermatan phylogeny did not change when tree reconstruction was carried out on the basis of a shorter dataset (24 ribosomal proteins, 3,535 amino acid positions) in which all ribosomal protein sequences had orthologs in acanthocephalans, bdelloids and monogononts. Support for a grouping of Bdelloidea and Acanthocephala was even higher when underlying this shorter dataset (PhyML: 83; Treefinder: 85; PhyloBayes: 0.92). Taken together our data suggest i) Syndermata being Spiralia with a close phylogenetic relation to Platyhelminthes and ii) the paraphyly of Eurotatoria, with iii) Bdelloidea being more closely related to Acanthocephala than to Monogononta.

thumbnailFigure 3. Phylogenetic tree reconstruction using bayesian inference (Phylobayes). Numbers at internal nodes represent posterior probabilities. Syndermata are shown as a basal spiralian taxon. Moreover, Eurotatoria appear paraphyletic, with Bdelloidea being more closely related to Acanthocephala than to Monogononta.

thumbnailFigure 4. Phylogenetic tree reconstruction using maximum likelihood (PhyML). Numbers at internal nodes represent bootstrap values. Syndermata are shown as a spiralian taxon, with a sistergroup relationship to Platyhelminthes. Moreover, Eurotatoria appear paraphyletic, with Bdelloidea being more closely related to Acanthocephala than to Monogononta.

thumbnailFigure 5. Phylogenetic tree reconstruction using maximum likelihood (Treefinder). Numbers at internal nodes represent expected likelihood weights. Syndermata are shown as a spiralian taxon, with a sistergroup relationship to Platyhelminthes. Moreover, Eurotatoria appear paraphyletic, with Bdelloidea being more closely related to Acanthocephala than to Monogononta.

Table 3. Results from hypotheses testing

Discussion

The phylogenetic position of Syndermata within Spiralia has been described previously based on molecular data such as 18S rRNA and 16SrRNA [15], 28S rRNA and 18S rRNA [20], ribosomal proteins [18] and on morphological characters like spiral cleavage, filiform sperm without accessory centriole and the subepidermal cerebral ganglion [5]. Likewise, a close phylogenetic relationship of Syndermata and Platyhelminthes within the spiralian clade agrees well with results from previous molecular and morphological approaches on metazoan phylogeny (e.g., [15,18,20,34]). On the other hand, conflicting results from the present tree reconstructions indicate that even the analysis of up to 79 ribosomal proteins from up to 29 species cannot settle the question of the definite phylogenetic relationship of Syndermata and Platyhelminthes. In agreement with Dunn et al. [34] we recommend an enlarged taxon sampling and the incorporation of more closely related groups such as gnathostomulids and micrognathozoans for resolving the question of the sistergroup relationships of Syndermata and Gnathifera, respectively.

In contrast to the still contradictory results regarding the phylogenetic position of Syndermata within Spiralia, our tree reconstructions consistently depict Bdelloidea as more closely related to Acanthocephala than to Monogononta, though with partly moderate support (Fig. 3, 4 and 5). As the moderate support values have been calculated on the basis of the full-length dataset, they might be due to missing data in one or more of the syndermatan lineages sampled. This is at least suggested by the higher statistical support for a clade Bdelloidea+Acanthocephala that could be inferred from the shorter dataset without these missing data (PhyML: 83; Treefinder: 85; PhyloBayes: 0.92). Regardless of differences in the support, analyses of both datasets lead to the same topology among the syndermatan representatives, whichever algorithm was employed. We take this high nodal stability (reflected also by unambiguous results from hypotheses testing) as evidence for reliability of the found grouping of Acanthocephala and Bdelloidea (see [35] for a discussion of nodal stability in the formulation of phylogenetic hypothesis). The future incorporation of EST data from Gnathostomulida and Micrognathozoa, and especially their use as outgroup, will be necessary to yield improved nodal support for the implicit paraphyly of Eurotatoria.

As another note of caution, one has to be aware that the present results as well as consistent results from more limited datasets [24-26] could in principle be influenced by an acceleration of sequence evolution on the branch of the only monogonont sampled, i.e. B. plicatilis. It is thus conceivable that a deviating mode of sequence evolution in B. plicatilis (as described for hsp82 [28]) triggered an attraction of the bdelloid and acanthocephalan branches. On the other hand, this is not very likely as the tree reconstruction methods employed herein (maximum likelihood and bayesian inference) are relatively robust to long-branch attraction. Moreover, Eurotatoria appeared paraphyletic in previous analyses comprising several monogononts [17,19] which cannot be explained by long-branch attraction due to an accelerated sequence evolution in B. plicatilis. We therefore do not believe that a faster sequence evolution along the B. plicatilis branch is causative for the found clustering of Bdelloidea and Acanthocephala.

The present evidence for a paraphyly of Eurotatoria is in apparent conflict with three out of the five competing hypotheses on the intra-syndermatan phylogeny, i.e. the Eurotatoria+Pararotatoria hypothesis (Fig. 1B), the classical Rotifera+Acanthocephala hypothesis (Fig. 1C), and the Eurotatoria+Acanthocephala hypothesis (Fig. 1D). At first sight, the observed grouping of Bdelloidea and Acanthocephala (under exclusion of Monogononta) rather supports the predictions of the Lemniscea hypothesis of Lorenzen ([23]). However, considering previous evidence from morphological [5,9,12] and molecular data [24] as well as from approaches combining both types of data [25,29], it is still possible that Seisonidea represent the true acanthocephalan sister taxon. However, it cannot be ruled out that Seisonidea are the sistergroup of Bdelloidea, Pararotatoria or Monogononta+Acanthocephala (see single trees in [19,25]).

Given the uncertain phylogenetic position of Seisonidea within Syndermata, one has to be cautious when inferring the evolution of morphological characters. On the other hand, the well supported closer relation of Bdelloidea to Acanthocephala, with exclusion of Monogononta ([19], present study), allows for some conclusions regarding the evolution of morphological characters that are not bound to the position of Seisonidea within Syndermata. It is thus very likely that the rotatory organ or corona underwent a (partial or total) reduction before the separation of Acanthocephala. A likewise reduction of a newly emerged character (wings) has for example been described in stick insects (Phasmatodea; [36]). Therefore the reduction of the rotatory organ only a few splits after its emergence at the base or within the syndermatan tree is not as unlikely as it might appear at first sight. Another implication of the grouping of Acanthocephala and Bdelloidea is that a retractable anterior end – whether in the shape of a rostrum in Bdelloidea or as a hooked proboscis in Acanthocephala – probably evolved before the separation of the acanthocephalan stem lineage as well. The reduction of the corona as well as the evolution of a retractable anterior end can easily be explained by different life-styles and patterns of locomotion in the syndermatan subtaxa: free living/free swimming in Monogononta; leech-like creeping/free living in Bdelloidea; leech-like creeping/epibiontic in Seisonidea; reduced motility/endoparasitic in Acanthocephala (see also [24,37-39]). Particularly the early evolution of a retractable anterior end might have represented a key event leading to the later evolution of the acanthocephalan endoparasitism, given the crucial role of the proboscis in the anchoring of adult acanthocephalans to the definitive hosts' intestinal wall [40].

Conclusion

Based on a dataset comprising sequences from up to 79 ribosomal proteins of up to 29 species, we provide evidence for the paraphyly of Eurotatoria. Irrespective of the tree reconstruction method and dataset used (and additionally supported by hypothesis testing) we found Bdelloidea to be more closely related to Acanthocephala than to Monogononta. Although data for Seisonidea have not been included in the dataset, the present findings allow for the rejection of three (Eurotatoria+Pararotatoria, Eurotatoria+Seisonidea, Eurotatoria+Acanthocephala) out of the presently five competing hypothesis regarding the phylogeny within Syndermata. On the other hand, additional data are needed to determine the actual acanthocephalan sistergroup (Seisonidea or Bdelloidea). Irrespective of these limitations it is very likely that a (partial or complete) reduction of the rotatory organ or corona occurred before the separation of Acanthocephala. Likewise, a retractable anterior end most likely emerged before the separation of the acanthocephalan stem lineage. Considering the importance of the proboscis for the attachment of acanthocephalans to the definite host's intestinal wall, the latter step can be regarded as a key event towards the evolution of acanthocephalan endoparasitism.

Methods

Isolation of RNA and cDNA library construction

Total RNA was extracted from frozen pooled specimen using column-based methods (Qiagen RNeasy Plant Mini Kit, Qiagen, Hilden, Germany). Quality of RNA was visually checked on agarose gels and mRNA was subsequently captured using the NucleoTrap mRNA kit (Macherey-Nagel, Düren, Germany) for B. plicatilis and the polyATract mRNA Isolation System III (Promega, Mannheim, Germany) for P. laevis and E. truttae. cDNA libraries were constructed at the Max Planck Institute for Molecular Genetics in Berlin (P. laevis) and the Institute of Molecular Genetics, University of Mainz (E. truttae, B. plicatilis) by primer extension (P. laevis, B. plicatilis) or LD-PCR (E. truttae), size fractionation and directional cloning applying the Creator SMART cDNA Libraries Kit (Clontech, Heidelberg, Germany) with the vectors pDNR-LIB or a modified pSPORT [41]. Clones containing cDNA inserts were sequenced from the 5' end on ABI 3730 capillary sequencer systems using BigDye chemistry (Applied Biosystems, Darmstadt, Germany).

EST processing

EST processing for P. laevis was accomplished at the Center for Integrative Bioinformatics in Vienna. Sequence chromatograms were first base-called and evaluated using the Phred application [42]. Vector, adaptor, poly-A tract and bacterial sequences were removed employing the software tools Lucy http://www.tigr.org webcite, SeqClean http://compbio.dfci.harvard.edu/tgi/software webcite, and CrossMatch http://www.phrap.org webcite, respectively. Clustering and assembly of the clipped sequences was performed using the TIGCL program package http://compbio.dfci.harvard.edu/tgi/software webcite by performing pairwise comparisons (MGIBlast) and a subsequent clustering step (CAP3). Low quality regions were then removed by Lucy. Finally, contigs were tentatively annotated by aligning them pairwise with the 25 best hits retrieved from NCBI's non-redundant protein database using the BlastX algorithm http://www.ncbi.nlm.nih.gov webcite. Alignment and computation of the resulting match scores, on which the annotation was based, were conducted by GeneWise [43] in order to account for frame shift errors.

ESTs for E. truttae and B. plicatilis were processed semi-automatically: removal of vector parts, polyA tails and bad quality sequence from sequence traces was performed by the SeqMan option of the DNASTAR program suite (Lasergene). Overlapping EST sequences were clustered using SeqMan (Lasergene). For B. plicatilis publicly available data from dbEST and the trace archives were included into the clustering process, and public data for Philodina roseola was clustered the same way. For annotation, EST cluster consensus sequences and EST singletons were subjected to BLASTX comparison against the SWISS-PROT protein database at NCBI http://www.ncbi.nlm.nih.gov/ webcite, using a BLAST client tool (Blastcl3, Blast software package, NCBI) setting the cut-off to 1*e-10. The EST data used in our analyses have been deposited in Genbank under the accession numbers [GenBank: AM849482AM849546 (P. laevis), AM980962AM980984 (E. truttae) and AM980946AM980961 (B. plicatilis)].

Sequence analysis and ribosomal proteins alignment

Ribosomal protein sequences were extracted from the newly obtained and publicly available EST data by their annotation. EST sequence contigs were checked for assembly errors by visual inspection and by comparison with corresponding sequences of related taxa, and translated into amino acid sequences. Gladyshev et al. [33] recently reported evidence for gene aquisition by horizontal gene transfer in two bdelloid species. Although it is unlikely that ribosomal proteins are subject to horizontal gene transfer, but as a precaution, we checked whether our data are influenced by horizontal gene transfer or not. Therefore we performed Blastp searches with our amino acid sequences and calculated the 'Alien Index' as introduced by Gladyshev et al. [33]. This index represents a measure of the orders of magnitude by which the BLAST E-value for the best metazoan hit differs from that for the best non-metazoan hit [33]. Additional ribosomal protein data were retrieved from the alignments compiled by Hausdorf et al. [18]. All ribosomal protein sequences obtained were aligned by the ClustalW algorithm using default parameters [44]. The resulting ribosomal protein alignments were inspected and adjusted manually for obviously misaligned positions using GeneDoc [45]. Questionably aligned positions were eliminated with GBlocks [46] using less stringent parameters. To test for the effect of missing data on present results [47], we assembled an additional dataset (24 ribosomal proteins, 3,535 amino acids) from which ribosomal protein sequences without acanthocephalan, bdelloid and/or monogonont orthologs were removed.

Phylogenetic reconstruction

The content of phylogenetic information of the alignments was estimated by the likelihood mapping approach as implemented in Tree-Puzzle 5.2 [48,49], testing all 23,751 possible quartets with exact parameter estimation.

Bayesian inference analyses based on the site-heterogenous CAT model (which allows the amino-acid replacement pattern to vary across a protein alignment; [50]) were performed using PhyloBayes v2.1c [51]. Two independent chains were run simultaneously for 11,210 points each. Chain equilibrium was estimated by plotting the log-likelihood and the alpha parameter as a function of the generation number. The first 500 points were subsequently discarded as burn-in. According to the divergence of bipartition frequencies, both chains reached convergence (maximal difference <0.08, mean difference <0.003), supported by the fact that both chains produced the same consensus tree topology. Taking every 10th sampled tree, a 50% majority rule consensus tree was finally computed using both chains.

ProtTest [52] was used to assess the appropriate model of sequence evolution for maximum likelihood-based tree reconstruction. As ribosomal proteins are likely to evolve similarly, the model was determined for the concatenated dataset, instead of for each single protein. Analyses were then conducted using PhyML [53] and Treefinder [54,55] with the rtREV+I+G+F substitution model [56] and 500 bootstrap replicates. Confidence values for the edges of the maximum likelihood tree (Treefinder) were computed by applying expected likelihood weights (ELWs) [57] to all local rearrangements of tree topology around an edge (1,000 replications). Trees produced in the course of the analysis were further edited using TreeView [58].

To test predefined phylogenetic hypotheses, we used constrained trees and the 'resolve multifurcations' option of Treefinder to obtain the maximum likelihood tree for a specified hypothesis. Thereafter we investigated whether the maximum likelihood trees for these hypotheses are part of the confidence set of trees applying the expected likelihood weights method [57].

Authors' contributions

AW participated in the design and coordination of the study, carried out the sequence analyses and alignments, participated in the phylogenetic reconstructions and wrote the manuscript. HH provided acanthocephalans, participated in the design of the study and wrote the manuscript. AM participated in the phylogenetic reconstructions. LB and GB provided monogononts. TH conceived the study, participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We gratefully acknowledge the work of Michael Kube and Richard Reinhardt (Max Planck Institute for Molecular Genetics, Berlin) in cDNA library construction and sequencing, Ingo Ebersberger, Sascha Strauß and Arndt von Haeseler (Max F. Perutz Laboratories, Center for Integrative Bioinformatics, Vienna) in processing ESTs in the framework of the DFG priority programme SPP1174 "Deep Metazoan Phylogeny", David Mark Welch and two anonymous referees for their helpful comments on the earlier version of this manuscript and Ana Rogulja-Ortmann for language editing. The work was funded by the Deutsche Forschungsgemeinschaft (DFG grant Ha 2103/4-1, SPP1174).

References

  1. Westheide W, Rieger R: Spezielle Zoologie, Teil 1. Einzeller und wirbellose Tiere. Stuttgart, Jena, NewYork: Verlag Gustav Fischer; 1996. PubMed Abstract OpenURL

  2. Ricci C: Are lemnisci and proboscis present in Bdelloidea?

    Hydrobiologia 1998, 387–388(4):93-96. OpenURL

  3. Segers H: Annotated checklist of the rotifers (Phylum Rotifera), with notes on nomenclature, taxonomy and distribution.

    Zootaxa 2007, (1564):3-104. OpenURL

  4. Sørensen MV, Segers H, Funch P: On a new Seison Grube, 1861 from coastal waters of Kenya, with a reappraisal of the classification of the Seisonida (Rotifera).

    Zoological Studies 2005, 44(1):34-43. OpenURL

  5. Ahlrichs WH: Ultrastruktur und Phylogenie von Seison nebaliae (Grube 1859) und Seison annulatus (Claus 1876). In Hypothesen und Verwandschaftsverhältnisse innerhalb der Bilateria. Göttingen: Cuvillier; 1995. OpenURL

  6. Wallace RL, Colburn RA: Phylogenetic-Relationships within Phylum Rotifera – Orders and Genus Notholca.

    Hydrobiologia 1989, 186:311-318. OpenURL

  7. Meyer A: Acanthocephala. In Bronns Klassen und Ordnungen des Tierreichs 4 (2/2). Edited by Bronn HG. Leipzig: Akademische Verlagsgesellschaft mbH; 1933:582. OpenURL

  8. Wallace RL, Ricci C, Melone G: A cladistic analysis of pseudocoelomate (aschelminth) morphology.

    Invertebrate Biology 1996, 115(2):104-112. OpenURL

  9. Ahlrichs WH: Epidermal ultrastructure of Seison nebaliae and Seison annulatus, and a comparison of epidermal structures within the Gnathifera.

    Zoomorphology 1997, 117(1):41-48. OpenURL

  10. Herlyn H, Ehlers U: Organisation of the praesoma in Acanthocephalus anguillae (Acanthocephala, Palaeacanthocephala) with special reference to the muscular system.

    Zoomorphology 2001, 121(1):13-18. OpenURL

  11. Herlyn H, Rohrig H: Ultrastructure and overall organization of ligament sac, uterine bell, uterus and vagina in Paratenuisentis ambiguus (Acanthocephala, Eoacanthocephala) – the character evolution within the Acanthocephala.

    Acta Zoologica 2003, 84(3):239-247. OpenURL

  12. Ahlrichs WH: Spermatogenesis and ultrastructure of the spermatozoa of Seison nebaliae (Syndermata).

    Zoomorphology 1998, 118(4):255-261. OpenURL

  13. Melone G, Ricci C, Segers H, Wallace RL: Phylogenetic relationsships of phylum Rotifera with emphasis on the families of Bdelloidea.

    Hydrobiologia 1998, 387–388:101-107. OpenURL

  14. Garey JR, Near TJ, Nonnemacher MR, Nadler SA: Molecular evidence for Acanthocephala as a subtaxon of Rotifera.

    Journal of Molecular Evolution 1996, 43(3):287-292. PubMed Abstract | Publisher Full Text OpenURL

  15. Garey JR, Schmidt-Rhaesa A, Near TJ, Nadler SA: The evolutionary relationships of rotifers and acanthocephalans.

    Hydrobiologia 1998, 388:83-91. OpenURL

  16. Giribet G, Distel DL, Polz M, Sterrer W, Wheeler WC: 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.

    Systematic Biology 2000, 49(3):539-562. PubMed Abstract OpenURL

  17. Giribet G, Sørensen MV, Funch P, Kristensen RM, Sterrer W: Investigations into the phylogenetic position of Micrognathozoa using four molecular loci.

    Cladistics-the International Journal of the Willi Hennig Society 2004, 20(1):1-13. OpenURL

  18. Hausdorf B, Helmkampf M, Meyer A, Witek A, Herlyn H, Bruchhaus I, Hankeln T, Struck TH, Lieb B: Spiralian phylogenomics supports the resurrection of Bryozoa comprising Ectoprocta and Entoprocta.

    Mol Biol Evol 2007, 24(12):2723-2729. PubMed Abstract | Publisher Full Text OpenURL

  19. Garcia-Varela M, Nadler SA: Phylogenetic relationships among Syndermata inferred from nuclear and mitochondrial gene sequences.

    Mol Phylogenet Evol 2006, 40(1):61-72. PubMed Abstract | Publisher Full Text OpenURL

  20. Passamaneck Y, Halanych KM: Lophotrochozoan phylogeny assessed with LSU and SSU data: Evidence of lophophorate polyphyly.

    Mol Phylogenet Evol 2006, 40(1):20-28. PubMed Abstract | Publisher Full Text OpenURL

  21. Sørensen MV, Funch P, Willerslev E, Hansen AJ, Olesen J: On the phylogeny of the Metazoa in the light of Cycliophora and Micrognathozoa.

    Zoologischer Anzeiger 2000, 239:297-318. OpenURL

  22. Herlyn H, Ehlers U: Ultrastructure and function of the pharynx of Gnathostomula paradoxa (Gnathostomulida).

    Zoomorphology 1997, 117(3):135-145. OpenURL

  23. Lorenzen S: Phylogenetic aspects of pseudocoelomate evolution. In The origins and relationships of lower invertebrates. Edited by Conway Morris S, George JD, Gibson R, Platt HM. Oxford: Clarendon Press; 1985:210-223. OpenURL

  24. Herlyn H, Piskurek O, Schmitz J, Ehlers U, Zischler H: The syndermatan phylogeny and the evolution of acanthocephalan endoparasitism as inferred from 18S rDNA sequences.

    Mol Phylogenet Evol 2003, 26(1):155-164. PubMed Abstract | Publisher Full Text OpenURL

  25. Zrzavy J: The interrelationships of metazoan parasites: a review of phylum- and higher-level hypotheses from recent morphological and molecular phylogenetic analyses.

    Folia Parasitologica 2001, 48(2):81-103. PubMed Abstract OpenURL

  26. Nielsen C: Animal Evolution. In Interrelationships of the Living phyla. Oxford: Oxford University Press; 1995. OpenURL

  27. Garcia-Varela M, de Leon GPP, de la Torre P, Cummings MP, Sarma SSS, Laclette JP: Phylogenetic relationships of Acanthocephala based on analysis of 18S ribosomal RNA gene sequences.

    Journal of Molecular Evolution 2000, 50(6):532-540. PubMed Abstract | Publisher Full Text OpenURL

  28. Mark Welch DB: Evidence from a protein-coding gene that acanthocephalans are rotifers.

    Invertebrate Biology 2000, 119(1):17-26. OpenURL

  29. Sørensen MV, Giribet G: A modern approach to rotiferan phylogeny: Combining morphological and molecular data.

    Mol Phylogenet Evol 2006, 40(2):585-608. PubMed Abstract | Publisher Full Text OpenURL

  30. Philippe H, Snell EA, Bapteste E, Lopez P, Holland PW, Casane D: Phylogenomics of eukaryotes: impact of missing data on large alignments.

    Mol Biol Evol 2004, 21(9):1740-1752. PubMed Abstract | Publisher Full Text OpenURL

  31. Hughes J, Longhorn SJ, Papadopoulou A, Theodorides K, de Riva A, Mejia-Chang M, Foster PG, Vogler AP: Dense taxonomic EST sampling and its applications for molecular systematics of the Coleoptera (beetles).

    Mol Biol Evol 2006, 23(2):268-278. PubMed Abstract | Publisher Full Text OpenURL

  32. Struck TH, Fisse F: Phylogenetic position of Nemertea derived from phylogenomic data.

    Mol Biol Evol 2008. PubMed Abstract | Publisher Full Text OpenURL

  33. Gladyshev EA, Meselson M, Arkhipova IR: Massive horizontal gene transfer in bdelloid rotifers.

    Science 2008, 320(5880):1210-1213. PubMed Abstract | Publisher Full Text OpenURL

  34. Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, et al.: Broad phylogenomic sampling improves resolution of the animal tree of life.

    Nature 2008, 452(7188):745-749. PubMed Abstract | Publisher Full Text OpenURL

  35. Giribet G: Stability in phylogenetic formulations and its relationship to nodal support.

    Systematic Biology 2003, 52(4):554-564. PubMed Abstract OpenURL

  36. Whiting MF, Bradler S, Maxwell T: Loss and recovery of wings in stick insects.

    Nature 2003, 421(6920):264-267. PubMed Abstract | Publisher Full Text OpenURL

  37. Near TJ, Garey JR, Nadler SA: Phylogenetic relationships of the Acanthocephala inferred from 18S ribosomal DNA sequences.

    Mol Phylogenet Evol 1998, 10(3):287-298. PubMed Abstract | Publisher Full Text OpenURL

  38. Herlyn H: Zur Ultrastruktur, Morphologie und Phylogenie der Acanthocephala. Berlin: Logos Verlag; 2000.

  39. Near TJ: Acanthocephalan phylogeny and the evolution of parasitism.

    Integrative and Comparative Biology 2002, 42(3):668-677. OpenURL

  40. Taraschewski H: Host-parasite interactions in acanthocephala: a morphological approach.

    Advances in Parasitology 2000, 46:1-179. PubMed Abstract OpenURL

  41. Wellenreuther R, Schupp I, Poustka A, Wiemann S: SMART amplification combined with cDNA size fractionation in order to obtain large full-length clones.

    BMC Genomics 2004, 5(1):36. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment.

    Genome Res 1998, 8(3):175-185. PubMed Abstract | Publisher Full Text OpenURL

  43. Birney E, Clamp M, Durbin R: GeneWise and Genomewise.

    Genome Res 2004, 14(5):988-995. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Thompson JD, Higgins DG, Gibson TJ: Clustal-W – Improving the Sensitivity of Progressive Multiple Sequence Alignment through Sequence Weighting, Position-Specific Gap Penalties and Weight Matrix Choice.

    Nucleic Acids Research 1994, 22(22):4673-4680. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Nicholas KB, N HB Jr, Deerfield DWI: GeneDoc: Analysis and Visualization of Genetic Variation.

    EMBNEWNEWS 1997, 4:14. OpenURL

  46. Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis.

    Mol Biol Evol 2000, 17(4):540-552. PubMed Abstract | Publisher Full Text OpenURL

  47. Hartmann S, Vision TJ: Using ESTs for phylogenomics: Can one accurately infer a phylogenetic tree from a gappy alignment?

    Bmc Evolutionary Biology 2008., 8 PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing.

    Bioinformatics 2002, 18(3):502-504. PubMed Abstract | Publisher Full Text OpenURL

  49. Strimmer K, von Haeseler A: Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignment.

    Proc Natl Acad Sci USA 1997, 94(13):6815-6819. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Lartillot N, Philippe H: A Bayesian mixture model for across-site heterogeneities in the amino-acid replacement process.

    Mol Biol Evol 2004, 21(6):1095-1109. PubMed Abstract | Publisher Full Text OpenURL

  51. Blanquart S, Lartillot N: A Bayesian compound stochastic process for modeling nonstationary and nonhomogeneous sequence evolution.

    Mol Biol Evol 2006, 23(11):2058-2071. PubMed Abstract | Publisher Full Text OpenURL

  52. Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution.

    Bioinformatics 2005, 21(9):2104-2105. PubMed Abstract | Publisher Full Text OpenURL

  53. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood.

    Syst Biol 2003, 52(5):696-704. PubMed Abstract OpenURL

  54. Jobb G: TREEFINDER.

    2007.

    Available from Munich, distributed by the author at http://www.treefinder.de webcite

  55. Jobb G, von Haeseler A, Strimmer K: TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics.

    Bmc Evolutionary Biology 2004., 4 PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Dimmic MW, Rest JS, Mindell DP, Goldstein RA: rtREV: an amino acid substitution matrix for inference of retrovirus and reverse transcriptase phylogeny.

    J Mol Evol 2002, 55(1):65-73. PubMed Abstract | Publisher Full Text OpenURL

  57. Strimmer K, Rambaut A: Inferring confidence sets of possibly misspecified gene trees.

    Proc Biol Sci 2002, 269(1487):137-142. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  58. Page RD: TreeView: an application to display phylogenetic trees on personal computers.

    Comput Appl Biosci 1996, 12(4):357-358. PubMed Abstract OpenURL