Email updates

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

Open Access Research article

Eurotatorian paraphyly: Revisiting phylogenetic relationships based on the complete mitochondrial genome sequence of Rotaria rotatoria (Bdelloidea: Rotifera: Syndermata)

Gi-Sik Min1 and Joong-Ki Park2*

Author Affiliations

1 Department of Biological Sciences, Inha University, Incheon 402-751, Republic of Korea

2 Department of Parasitology, College of Medicine, Chungbuk National University, Cheongju 361-763, Republic of Korea

For all author emails, please log on.

BMC Genomics 2009, 10:533  doi:10.1186/1471-2164-10-533


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


Received:16 July 2009
Accepted:17 November 2009
Published:17 November 2009

© 2009 Min and Park; 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 Syndermata (Rotifera+Acanthocephala) is one of the best model systems for studying the evolutionary origins and persistence of different life styles because it contains a series of lineage-specific life histories: Monogononta (cyclic parthenogenetic and free-living), Bdelloidea (entirely parthenogenetic and mostly benthic dweller), Seisonidea (exclusively bisexual and epizoic or ectoparasitic), and Acanthocephala (sexual and obligatory endoparasitic). Providing phylogenetic resolution to the question of Eurotatoria (Monogononta and Bdelloidea) monophyly versus paraphyly is a key factor for better understanding the evolution of different life styles, yet this matter is not clearly resolved. In this study, we revisited this issue based on comparative analysis of complete mitochondrial genome information for major groups of the Syndermata.

Results

We determined the first complete mitochondrial genome sequences (15,319 bp) of a bdelloid rotifer, Rotaria rotatoria. In order to examine the validity of Eurotatoria (Monogononta and Bdelloidea) monophyly/paraphyly, we performed phylogenetic analysis of amino acid sequences for eleven protein-coding genes sampled from a wide variety of bilaterian representatives. The resulting mitochondrial genome trees, inferred using different algorithms, consistently failed to recover Monogononta and Bdelloidea as monophyletic, but instead identified them as a paraphyletic assemblage. Bdelloidea (as represented by R. rotatoria) shares most common ancestry with Acanthocephala (as represented by L. thecatus) rather than with monogonont B. plicatilis, the other representative of Eurotatoria.

Conclusion

Comparisons of inferred amino acid sequence and gene arrangement patterns with those of other metazoan mtDNAs (including those of acanthocephalan L. thecatus and monogonont B. plicatilis) support the hypothesis that Bdelloidea shares most common ancestry with Acanthocephala rather than with Monogononta. From this finding, we suggest that the obligatory asexuality of bdelloideans may have secondarily derived from some other preexisting condition in earlier lineage of rotifers. Providing a more complete assessment of phylogenetic relationships and inferring patterns of evolution of different types of life styles among Syndermata awaits comparisons requiring mitochondrial genome sequencing of Seisonidea.

Background

The Rotifera (also called Rotatoria) is a group of aquatic micrometazoans, mostly less than a millimeter in size. It includes more than 2,000 described species that usually occur in freshwater and marine environments throughout the world, but some species are found in wet terrestrial habitats, such as moist soil, or mosses and lichens living on fallen trees and rocks [1-3]. This group is usually distinguished from other metazoans by presence of the corona ("wheel organ" = the crown of cilia) located in the cephalic region, which is used for locomotion and food gathering. Due to their great abundance and high reproductive potential, rotifers have been considered to play a significant role in the food webs of certain freshwater environments through their involvement in energy cycling and nutrient transfer [4-6]. Most authorities accept that the Rotifera consists of three classes, each having a unique reproductive strategy [2,7,8]: Bdelloidea (exclusively parthenogenetic), Monogononta (cyclic parthenogenetic with facultative sexual reproduction), and Seisonidea (exclusively bisexual). Of these, the class Bdelloidea is very unique in that its species (several hundred) are exclusively female. This is the largest metazoan asexual group where no sexual reproduction has ever been reported and represents an ancient origin of asexuality with great evolutionary success in the diversification of the species [9,10]. The class Monogononta, representing the largest group of rotifers, comprises more than 1,500 species. They are mostly found in freshwater, brackish and marine waters and are characterized by a unique life cycle, cyclical parthenogenesis, that alternates between two reproductive modes, i.e. amictic and mictic phases of reproduction according to the absence/presence of males, respectively [2,5]. The class Seisonidea consists of three epizoic species (belong to two genera, Seison and Paraseison [11]) on marine crustaceans (Nebalia species) and all of these are known to reproduce by amphimixis [12].

There is little doubt about the close relationship between Rotifera and Acanthocephala, the clade known as Syndermata because they share the feature of the syncytial epidermis [13,14]. This affinity has received broad support from many earlier studies utilizing different sources of phylogenetic information: morphology [15-17], SSU rDNA sequences [18-20], combined analysis of molecular and morphological characters [21,22], and combined analysis of SSU and LSU ribosomal DNA sequences [23]. In contrast to syndermatan monophyly, internal phylogeny within the clade (among Bdelloidea, Monogononta, Seisonidea, and Acanthocephala) has been the subject of relatively vigorous contention [24-27]. Most controversies are related to Eurotatorian monophyly versus paraphyly, although there are differences regarding the phylogenetic position of Seisonidea within the clade. The monophyly/non-monophyly of eurotatorians (Monogononta and Bdelloidea) relative to Acanthocephala has received much attention because this issue plays a key role in understanding the evolution and ecological diversity found among major groups of Syndermata. Lorenzen [15] first recognized the lemnisci (paired projections of the neck epidermis into the body cavity) and proboscis (an invaginable anterior part of the body) as synapomorphic characters for unifying the Bdelloidea and Acanthocephala (inferred the latter as being the highly specialized sister group to the former; see [24] for different view), and this relationship has also been recovered in many subsequent molecular analyses (SSU+LSU+mtDNA cox1 [26]; SSU [19]; SSU+mtDNA 16S [20]; SSU+LSU+histone H3+mtDNA cox1[28]). Based on morphological and molecular perspectives, Garey et al. [19] recognized Acanthocephala as a subtaxon of Rotifera and united these two groups under the superclass 'Lemniscea' (Lemniscea hypothesis-eurotatorian paraphyly). However, eurotatorian monophyly (Monogononta+Bdelloidea) and sister-group relationship of Seisonidea and Acanthocephala have also been suggested from some earlier studies based on morphological evidence (ultrastructural similarity [14,29]; cladistic analysis of morphological dataset [7,30,31]) or from molecular analysis (the partial sequence of nuclear heat-shock protein hsp82 [8]; SSU rDNA [32]; combined analysis of hsp82+SSU data [33]). Furthermore, employing different methods for phylogenetic analysis resulted in inconsistent tree topologies when SSU data were analyzed [34]. Although the most recently published work from EST-based phylogenomic analysis supported Eurotatoria paraphyly [27], but the phylogenetic issue regarding the Eurotatoria monophyly/paraphyly still awaits other types of data for corroboration.

With a very few exceptions, metazoan mitochondrial genomes are circular DNA molecules (mostly less than 16 kb in size) that encode 37 genes: 13 protein-coding genes (atp8 is missing in many nematode and flatworm species so far reported), two ribosomal RNA genes, and 22 transfer RNA genes [35,36]. Due to its universality and remarkably stable feature in genome content across various metazoan phyla, comparisons of the mitochondrial genome information (e.g., nucleotide sequence, amino acid sequence and gene order rearrangement) have often proven useful for reconstructing the deep node phylogeny and for assessing the phylogenetic relationships among closely related species [36-39]. In recent years, there has been an unprecedented increase in mitochondrial genomic surveys in relation to phylogenetic comparisons of a variety of animal groups. Mitochondrial genome information has now become available for more than a thousand animal species (See NCBI metazoan mitochondrial genome resources). However, the distribution of completely characterized mitochondrial genome sequences has been strongly biased across the metazoan taxa: the subphylum Vertebrata and the phylum Arthropoda account for more than 80% of the metazoan mitochondrial genome data determined so far, whereas there are still a considerable number of metazoan phyla for which there is limited mitochondrial genome information, and some phyla have never been investigated [40,41]. Complete mitochondrial genome sequencing from poorly investigated groups is needed to supplement a gap in our current understanding of mitochondrial DNA evolution in the metazoa.

Mitochondrial genome information from the Syndermata was reported, for the first time, from the acanthocephalan species Leptorhynchoides thecatus [42] and recently thereafter from the monogonont rotifer Brachionus plicatilis [43]. Despite belonging to the same clade (Syndermata), the mitochondrial genomes reported for these two organisms did not share many characteristics in their organization such as chromosome structure, gene order, codon usage and the secondary structure of tRNA molecules. For example, the mitochondrial genome of B. plicatilis, a representative of the Monogononta, is encoded in two separate mitochondrial chromosomes, each having different gene content and copy number. In contrast, as in most metazonas, all genes of L. thecatus mtDNA are encoded in a single type of circular mitochondrial DNA molecule. The lack of common features in the mitochondrial genomes between these species necessitates additional characterization of mitochondrial genomes from the other major groups of the syndermata. Information from Bdelloidea is expected to supplement our understanding of syndermatan mitochondrial genome evolution and provide utility as a molecular marker in resolving internal phylogeny among the major groups of the Syndermata. To this end, we characterized the first complete mitochondrial genome sequence of the bdelloid species Rotaria rotatoria, and compared its mitochondrial genome information with other syndermatan species in order to investigate phylogenetic issues regarding the monophyly or paraphyly of Eurotatoria.

Results and Discussion

General features of the R. rotatoria mitochondrial genome

The mitochondrial genome of R. rotatoria is 15,319 bp in length (GenBank accession no. GQ304898) and contains 12 protein-coding genes (lacking atp8), two rRNA genes, and as is occasionally found in some other metazoans, all but one (trnC) of the typical 22 tRNAs genes. All genes are encoded in the same direction in a single circular mitochondrial DNA (Fig. 1), unlike that of the monogonont rotifer B. plicatilis, in which mitochondrial genes are encoded in two separated circular chromosomes [43]. The genome organization including the gene order, length of overlapping regions between the genes and length of intergenic spacer regions is shown in Table 1. The nucleotide composition of the entire R. rotatoria mtDNA sequence is 29.5% A, 43.7% T, 17.4% G, and 9.5% C (Table 2). The overall A+T content (73.1%) approximates more or less that of the acanthocephalan L. thecatus (71.5% in A+T content), but is significantly higher than those found in the monogonont rotifer B. plicatilis (A+T contents of 63.9% and 62.9% for mtDNA-I and mtDNA-II, respectively). The average of A+T composition of protein-coding genes (73.7%) is very similar to that of the entire sequence, but there is noticeable compositional difference across the codon positions: the average A+T content of the third position (80.9%) is noticeably higher than those of the first (67.8%) and second codon (72.4%) positions, indicating that there is a strong bias toward T (47.2%) and A (33.7%) in third codon position.

Table 1. The mitochondrial genome organization of Rotaria rotatoria

Table 2. Nucleotide composition of the mitochondrial genome of Rotaria rotatoria

thumbnailFigure 1. Circular representation of the mitochondrial genome of Rotaria rotatoria. All genes are encoded in the same direction and 21 tRNA genes are denoted by the one-letter abbreviation and two leucine and two serine tRNA genes are labeled, according to their anticodon sequence, as L1 (trnL-uag), L2 (trnL-uaa), S1 (trnS-ucu), and S2 (trnS-uga), respectively. The intergenic non-coding regions with notable size (>200 bp) are denoted as NCR1-NCR5, respectively.

Protein-coding genes

Twelve protein-coding genes were identified using NCBI ORF Finder and by comparing their sequence with those of homologous genes reported from the acanthocephalan L. thecatus [42] and monogonont B. plicatilis [43]. Even after an exhaustive search using BLAST, we failed to discover any atp8-like protein sequence. The lack of atp8 in the mtDNA genome is not very rare, and quite common in most nematode and platyhelminth species reported so far (cf. Trichinella spiralis, a nematode species where the atp8 exists [44]).

In many metazoan mtDNAs, the codon usage is strongly biased, especially in the third position of synonymous codons within amino acid families [35,45-47]. In the protein-coding genes of R. rotatoria mtDNA, the higher A+T content is related to amino acid composition with higher abundance of A+T-rich codons. Amino acids encoded by T-rich and A-rich codons (more than two Ts and As in a triplet, respectively) account for 47.3% and 22.2% of the total amino acid composition, respectively, totaling 69.5% of the entire protein sequences (Table 3). Another factor that causes nucleotide bias is differential preference of thymine (T) or adenine (A) in the third position of codons. The average relative nucleotide frequency in the third position is 47.2% (T), 33.7% (A), 13.7% (G), and 5.4% (C), respectively. Indeed, both in twofold- and fourfold-degenerate codon families, the T-ending codons are greatly favored in most synonymous codons, whereas C-ending codons are evidently avoided. The relative frequency of TTT (9.5%) coding phenylalanine is much higher than that of its degenerate codon TTC (1.5%). Another notable example of biased codon usage is that the codon ATT (8.0%) is much more frequently used than ATC (0.7%) in the isoleucine family. The propensity of a preference for T, with an apparent bias against C may reflect mutational bias at the third codon position as postulated in other metazoans [48]. Of 12 protein-coding genes, six (atp6, cob, nad4L, cox2, nad3, and nad2) use ATG as the initiation codon, whereas others initiate with ATT (cox1, nad6, nad1, and nad5), or ATA (nad4 and cox3), respectively (Table 1). The majority of genes end with complete termination codon TAA (cox1, nad1, cob, nad4L, nad4, nad5, cox3, and nad3) or TAG (nad6 and atp6), but cox2 and nad2 terminate with the incomplete stop codon T, which is presumed to become a complete termination codon through post-transcriptional polyadenylation [49].

Table 3. Codon usage for 12 protein-coding genes of the mitochondrial genome of Rotaria rotatoria

Transfer RNA and ribosomal RNA genes

Implementation of tRNAscan algorithm failed to find any tRNA-like secondary structure, but we manually identified them by eye. Of 22 tRNA genes normally found in most other metazoan mtDNAs, 21 tRNA-like nucleotide segments (ranging from 47 to 65 bp in size) except for trnC can be folded into a cloverleaf secondary structure with some mismatches or incomplete configuration of the DHU and/or TΨC arms (Additional file 1). In spite of the exhaustive tRNA search, we were not able to identify any candidates for trnC with high confidence. Out of 21 tRNAs found, only 10 (trnA, trnL1, trnM, trnF, trnP, trnS1, trnS2, trnT, trnY, and trnV) display a typical cloverleaf-like structure equipped with DHU and TΨC arms. Many other tRNAs appear to lack a TΨC arm (6 tRNAs; trnD, trnQ, trnE, trnG, trnH, and trnK) or a DHU arm (4 tRNAs; trnN, trnI, trnL2, and trnW), or both (1 tRNAs; trnR). The rarity of having a typical cloverleaf-like secondary structure was also found in the acanthocephalan species Leptorhynchoides thecatus [42] where almost none of the inferred tRNAs displayed a typical stem-and-loop configuration in the DHU and TΨC arms. This differs from most other metazoan mtDNAs, including the monogonont rotifer Brachionus plicatilis [43]. Indeed, the putative secondary structure of 22 tRNAs found in B. plicatilis contains a typical cloverleaf structure comprising a amino-acyl arm (a stem of seven nucleotide pairs; ntp), a DHU arm (a stem of 2-4 ntp with a 2-8 nt loop), an anticodon arm (a stem of 5 ntp with an anticodon) and a TΨC arm (a stem of 2-4 ntp with a 2-7 nt loop).

Additional file 1. Predicted secondary structures of the 21 mitochondrial tRNAs of R. rotatoria.

Format: EPS Size: 1.4MB Download fileOpen Data

Based on sequence comparison with L. thecatus and B. plicatilis, two ribosomal RNA genes were identified: The rrnL (529 bp) is located between trnY and trnL1 as found in L. thecatus, but its size is remarkably smaller than any other metazoan large ribosomal subunit molecule (mostly larger than 1 kb) including those of other syndermatans reported thus far (e.g., 925 bp and 1,107 bp in L. thecatus and B. plicatilis, respectively). The rrnS (521 bp) is located between trnM and cox2.

Non-coding regions

A total of 19 intergenic non-coding sequences, ranging from 1 to 806 bp in size (3,092 bp in total accounting for 20.2% of entire genome length), were detected. Of these, five intergenic non-coding regions (NCR) are prominent by their significant lengths (≥ 200 bp), ranging from 236 bp to 806 bp (535 bp-NCR1 [between nad3 and trnT], 806 bp-NCR2 [between trnT and trnP], 760 bp-NCR3 [between trnR and trnS2], 316 bp-NCR4 [between trnS2 and trnV], and 236 bp-NCR5 [between trnV and trnS1, respectively). Despite an exhaustive search using the NCBI ORF Finder, we were not able to find ORF-like candidates of significant length from these spacer regions. The relatively larger genome size of the R. rotatoria mtDNA is attributed to the abundance of conspicuously long, unassigned spacer regions, compared to that of acanthocephalan L. thecatus (13,888 bp [42]). The region located between trnR and trnS2 (NCR3) contains 24 tandemly repeated units of a 18-nt sequences (5'-NRRWYBTYRNGHRRYYYY-3'), each having potential to be folded into a stem-loop structure (a stem of 5-6 ntp with a 6-8 nt loop). Such repeats, known as variable numbers of tandem repeats (VNTRs), have been reported from a variety of animal mtDNAs and they have been utilized as molecular markers for detecting subdivision of the population in ecological studies (See [50] for a review).

Mitochondrial molecular phylogeny

Phylogenetic position of Platyzoa within Lophotrochozoa

The Platyzoa refers to an assemblage of some acoelomate and pseudocoelomate animal taxa including Platyhelminthes, Rotifera, Acanthocephala, Gastrotricha, and Gnathostomulida [51]. Phylogenetic affinity of these groups with other lophotrochozoans has not yet been confirmed. Results of recent molecular analyses depicted Platyzoa as a sister to Trochozoa (morphology+SSU dataset [21]), or as a derived subclade of Lophotrochozoa (LSU+SSU dataset [23]. This disagreement still awaits unambiguous resolution. In this study, we performed phylogenetic analyses for the amino acid dataset of major representatives of Bilateria with two cnidarian sequences (Aurelia aurita and Acropora tenuis) as outgroups. Resulting trees from both ML and BI methods showed that the protostomes and deuterostomes were divided into two separate clades with very high support (Fig. 2). Within the protostome clade, monophyly of Ecdysozoa and monophyly of Lophotrochozoa were also found, but the dichotomy between these two groups was collapsed when the nematode sequences (Caenorhabditis elegans, Trichinella spiralis) were included in the analysis. In all cases, nematode sequences were nested within the Lophotrochozoan clade rather than grouped with ecdysozoan members (Additional file 2), and their position varied considerably according to different phylogenetic methods. Their unusual position appears to be artificial and due to the long-branch attraction (LBA) problem as documented in many previous papers [52-54]. Therefore, in order to reduce this factor in the tree reconstruction, we excluded nematode sequences from subsequent analyses. After the exclusion of nematodes, the concatenated amino acid sequence dataset consisting of 1,556 homologous positions was used for the subsequent analyses. Note that monophyly of Lophotrochozoa and monophyly of Ecdysozoa were recovered after removal of nematode sequences from the analyses, and support values for clades were noticeably improved in all analyses (see table description of Additional file 3 for the phylogenetic position of Chaetognatha within the metozoa). The monophyletic grouping of ecdysozoan members received 100% BP (bootstrap percentage) in ML, and 1.00 BPP (Bayesian posterior probability). The Lophotrochozoan clade was also well supported with 1.00 BPP in BI, but received a very weak bootstrap support value of 56% in ML analysis (Fig. 2). Although the Lophotrochozoan clade is not strongly supported, the current amino acid sequence data support the Ecdysozoa/Lophotrochozoa split. It is particularly noteworthy that both phylogenetic methods (BI and ML) recovered Platyhelminthes, Rotifera, and Acanthocephala (all platyzoan members included in this study) as monophyletic lineages in the Lophotrochozoan clade and nodal support for Platyzoa was very strong (100% in both BI and ML analysis). This result is in agreement with an earlier analysis based on a combined dataset of LSU and SSU sequences that recognized the Platyzoa as a derived monophyletic assemblage within the Lophotrochozoa [23]. As yet, the availability of complete mitochondrial genome information is limited to Platyhelminthes, Acanthocephala, and Rotifera, and further evidence from Gnathostomulida and Gastrotricha is required to confirm the monophyly of Platyzoa and its derived position within the Lophotrochozon clade (sensu [51]).

Additional file 2. Phylogenetic tree inferred from Bayesian analysis for metazoan mitochondrial amino acid sequences before the exclusion of nematode species from the analysis. Numbers above/below branches are Bayesian posterior probability (BPP) and local rearrangement-expected likelihood weight (LR-ELW) edge support estimated from Bayesian and maximum likelihood analyses, respectively (BPP/LR-ELW). The branches that are supported with values of ≤50% or not consistent between Bayesian and maximum likelihood methods in their positions are represented by "-".

Format: EPS Size: 1.1MB Download fileOpen Data

Additional file 3. Taxa/species used for phylogenetic analysis in this study. The Chaetognatha was traditionally considered a member of deuterostomes, but recent molecular analysis using SSU rDNA sequence suggested that they are not deuterostomes [77,78]. Although phylogenetic position of Chaetognatha with other metazoan phyla still remained enigmatic, recent molecular phylogenetic surveys including comparative analysis of the complete mitochondrial genome often suggest its affinity to protostomes (for more details, see [21,79-81]). Therefore, we tentatively included two chaetognath species (Spadella cephaloptera and Paraspadella gotoi) in the Protostome clade for phylogenetic analysis.

Format: DOC Size: 69KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

thumbnailFigure 2. Mitochondrial gene tree from Bayesian analysis showing the phylogenetic relationships among 35 metazoan species. The tree topologies from Bayesian analysis and maximum likelihood are very similar. Numbers above/below branches are Bayesian posterior probability (BPP) and local rearrangement-expected likelihood weight (LR-ELW) edge support, and bootstrap (BP) values from maximum likelihood analysis, respectively (BPP/LR-ELW/BP). The branches that are supported with values of ≤50% or not consistent in their positions between Bayesian and maximum likelihood methods are represented by "-".

Phylogenetic implications for Syndermata: Eurotatoria paraphyly

Phylogenetic relationships among Syndermata, in particular with regard to Eurotatoria monophyly/paraphyly are still vigorously debated. The answer to this question is prerequisite for correct understanding of the evolution of different types of life styles found among major groups of Syndermata [27]. In this study, we revisited this issue based on phylogenetic analysis of amino acid sequence for eleven protein-coding genes (except atp6 and atp8) of the complete mitochondrial genome. The resulting trees both from Bayesian and maximum likelihood analyses recovered syndermatan members (Monogononta, Bdelloidea, and Acanthocephala) as a monophyletic group, which is in turn grouped with Platyhelminthes (Fig. 2). The phylogenetic affinity of the Syndermata with Platyhelminthes in this study is concordant with some earlier publications asserting their sister-relationship [23,55]. It is particularly significant to note that both BI and ML methods placed Bdelloidea and Acanthocephala as sister taxa and this relationship was very strongly supported by 1.00 BPP in BI analysis, but received relatively weaker support (65% BP) in ML analysis. Moreover, none of the phylogenetic methods recovered the two eurotatorian species, Rotaria rotatoria (Bdelloidea) and Brachionus plicatilis (Monogononta) as monophyletic, but instead always depicted them as a paraphyletic assemblage. Bdelloidea represented by R. rotatoria appears to share most recent common ancestry with Acanthocephala rather than with the monogonont B. plicatilis. The sister-group position of the Monogononta to the Bdelloidea+Acanthocephala, received strong support from BI (0.99 BPP) and moderate support from ML (70% BP), respectively. This result indicates a firm support for the 'Eurotatorian paraphyly' hypothesis. In all cases of the analysis, monophyly of Bdelloidea + Acanthocephala relative to Monogononta is always very robust. In order to assess whether our data rejects the alternative hypothesis (Eurotatoria monophyly), we performed statistical comparison of the likelihood scores for the best ML tree without constraint versus the best tree with the topological constraint of Eurotatoria monophyly using theTreefinder program. The result of topology test failed to detect a significant difference between these two competing hypotheses under various criteria (ELW, SH and AU; Table 4). Although Eurotatoria monophyly was not significantly worse interpretation of these data by tree topology comparison test, the optimal phylogenetic trees from different methods unanimously depicted B. plicatilis (Monogononta) sister to two other syndermatans (bdelloidean R. rotatoria and acanthocephalan L. thecatus), regardless of different options applied in the data analysis (including or excluding inferred gaps in the analysis). In addition, this relation (Eurotatoria paraphyly) is further corroborated by sharing greater similarity in gene arrangement between R. rotatoria and L. thecatus (next section).

Table 4. Results of tree topology test using different criteria: expected likelihood weights (ELW), Shimodaira-Hasegawa (SH), and approximately unbiased (AU)

A close relationship between Bdelloidea and Acanthocephala is also evident when gene arrangement patterns of R. rotatoria mtDNA are compared with other metazoan animal groups. Gene arrangement of R. rotatoria mtDNA is very different from those found in other metazoans, but shows the highest similarity to the acanthocephalan L. thecatus. The gene order comparison between R. rotatoria and L. thecatus identified the gene order of trnN-nad2-cox1-trnG, trnY-rrnL-trnL1-nad6, and nad4L-nad4-trnH-nad5 shared between the representatives of the Bdelloidea and Acanthocephala, respectively (Fig. 3). The gene order of nad4L-nad4-trnH-nad5 can also be found across many various metazoan mtDNAs, such as Bryozoa [56], Brachiura (Arthropoda [57]), Gastropoda (Mollusca [58]), Vertebrata (Teleostei [59] and Amphibia [60]), and Priapulida [54]. However, this conserved gene order is not shared by the B. plicatilis mtDNA and we suggest it is due to the idiosyncratic features of B. plicatilis mtDNA in which 36 genes are encoded separately into bipartite mitochondrial genomes [43]. In contrast, note that the gene order similarity found between R. rotatoria and L. thecatus (trnN-nad2-cox1-trnG and trnY-rrnL-trnL1-nad6) is very unique among the metazoans and this shared gene rearrangement lends another line of strong support for their common ancestry. The phylogenetic conclusion of Eurotatoria paraphyly suggested from the present study is also consistent with the most recent EST based phylogenomic analysis [27].

thumbnailFigure 3. Linearized comparison of the mitochondrial gene arrangement of three syndermatan species. Gene and genome size are not to scale. All genes are transcribed in the same direction (from left to right). The tRNAs are labeled by single-letter abbreviations. Gene clusters shared between L. thecatus (Acanthocephala) and R. rotatoria (Bdelloidea) are represented by shadowed areas.

In general, the evolutionary potential for adaptive radiation and its subsequent diversification of extant biota are often associated with preferential adoption of different life styles during evolutionary history. In most cases, a robust phylogeny is a prerequisite to correctly understand the evolution of life styles among the groups. The Syndermata (Rotifera+Acanthocephala) is one of the best model systems for studying the evolutionary origin and persistence of different life styles because it shows a variety of lineage-specific reproduction modes. Of these, Bdelloidea is particularly among the most unique in that it contains exclusively several hundred species and all individuals are females, the largest metazoan group with no sexual reproduction that has ever been documented [9]. This group is also considered to represent an ancient origin of asexuality with a great evolutionary success in their diversification of asexual species [9,10]. There are some additional cytological [61] and molecular evidence [62] suggesting that the evolution and maintenance of asexuality in Bdelloidea is of relatively ancient origin, dating back to more than tens of millions of years ago. Inferring the history of the obligatory asexual class Bdelloidea within the syndermatan phylogenetic framework is a central topic in evolutionary biology of asexuality because these organisms have long been regarded as perfect candidates to test the advantage of sexual reproduction [63]. The phylogenetic conclusion obtained from this comparative mitochondrial genome study provides essential clues in extrapolating the directionality of long-term asexual evolution of bdelloidean lineage. The sister group position of Monogononta to the Bdelloidea/Acanthocephala clade suggests that the obligatory asexuality of bdelloideans may have been derived secondarily from some other preexisting condition in earlier lineage of rotifers. However, because no mitochondrial genome information is available from Seisonidea, we are not able to postulate what is the most likely ancestral condition that gave rise to bdelloidean asexuality. There are only a few complete mitochondrial DNA sequences available from the Syndermata: B. plicatilis (Monogononta [43]), R. rotatoria (Bdelloidea; this study) and L. thecatus (Acanthocephala [42]). More information, especially from Seisonidea needs to be obtained to confirm the phylogenetic position of Seisonidea within Syndermata. This further phylolgenetic information would be a very substantial ingredient for precisely estimating the origin of asexuality in Bdelloidea and at the same time, contribute to a better understanding of mitochondrial genome evolution within the syndermatan phylogenetic framework.

Conclusion

In this study, we revisited phylogenetic relationships of Eurotatoria (monophyly versus paraphyly) based on comparative analysis of the complete mitochondrial genome information for major groups of the Syndermata. For this purpose, we determined the first complete mitochondrial genome sequence (15,319 bp) of a bdelloid rotifer, Rotaria rotatoria. Comparisons of inferred amino acid sequence and gene arrangement pattern with those of other syndermatan mtDNAs (from the acanthocephalan L. thecatus and monogonont B. plicatilis) supports the hypothesis that the Bdelloidea shares a most recent common ancestor with Acanthocephala, and not with Monogononta. From this finding, we suggest that the obligatory asexuality of bdelloideans may have secondarily derived from some preexisting condition in earlier lineages of rotifers. Definitively, testing this question of phylogenetic relationships and evolution of different types of life style among the major members of Syndermata awaits further investigation of mitochondrial genome sequencing from Seisonidea.

Methods

Sampling and molecular techniques

Rotifers were collected from a small pond at Inha University campus in South Korea (37°26'58.08"N, 126°39'20.09"E) and washed several times with distilled water. Total genomic DNA was extracted from pooled worms using a QIAamp tissue kit (QIAGEN Inc.) and used as a template for PCR amplification. Four small fragments of Rotaria rotatoria mtDNA (ranging in size from ~296 to 712 bp) were initially PCR-amplified using corresponding primer sets for each of four gene regions (cox1 [LCO1490/HCO2198], lrRNA [16SA/16SB], cob [Cytb-Uni5-2/Cytb-Uni3-2], and cox2 [CO2-Uni5/CO2-Uni3]; see Table 5). PCR reactions were conducted in a total volume of 50 μl reaction mixture containing 0.13 μg/μl of genomic DNA, 10× PCR buffer, 10 mM dNTP mixture, 10 pmole each primer, 25 mM MgCl2, and 2.5 units Taq polymerase (Roche Co.) using the following cycling conditions: 1 cycle (92°C for 2 min), 35 cycles (92°C for 1 min, 45°C to 55°C for 30 sec, 72°C for 1 min 30 sec), and 1 cycle (72°C for 10 min). The nucleotide sequence determined from these four fragments was then used to design R. rotatoria mtDNA-specific primers for long PCR amplification (see also Table 5 for details of the primer information). Four overlapping fragments (ranging in size from 1,559 bp to 5,481 bp) covering the entire mitochondrial genome of R. rotatoria were amplified using the Expand Long Template PCR System (Roche, USA) under the following conditions: 1 cycle of 2 min at 94°C (initial denaturation), 35 cycles of denaturation-primer annealing-elongation (10s at 95°C, 1 min at 55°C, and 10 min at 68°C), and 1 cycle of the final extension (10 min at 68°C). A negative control (no template) was included for every PCR run to detect any potential contamination of the PCR products. The amplified PCR products were isolated on a 1.0% agarose gel containing crystal violet, excised and extracted according to the TOPO XL gel-purifying protocol (Invitrogen Co.). Purified DNAs were cloned into E. coli competent cells using TOPO XL PCR Cloning kit, as recommended by manufacturer. Sequencing reactions of the target fragments were performed in both directions by 'primer walking' using a Big Dye Terminator Cycle-Sequencing Kit (Applied Biosystems), and overlapping fragments were assembled to complete the sequence of the entire genome.

Table 5. Primer sequence information used in this study

Gene Annotation

Twelve protein-coding genes and two ribosomal RNA genes of R. rotatoria were identified by sequence comparison with the mtDNA sequences of the acanthocephalan L. thecatus and monogonont rotifer B. plicatilis and by using the NCBI ORF (open reading frame) Finder. We attempted to find tRNAs using tRNAscan-SE 1.21, but no tRNA-like structure was detected from this search. Therefore, tRNA genes were identified by searching for anticodon consensus motif sequences (TxxxR; xxx = anticodon) and by recognizing potential secondary structures by eye. In this process, we were aided by a web-based automatic annotation program for organellar genomes (DOGMA [64]).

Phylogenetic analysis

Eleven mitochondrial protein-coding genes (except atp6 and atp8, which are highly variable in their lengths among the groups) from 35 species (6 ecdysozoans, 22 lophotrochozoans, 5 deuterostomes, and two cnidarian outgroups) representing major clades of the Bilateria, including those of R. rotatoria mtDNA, were used in phylogenetic analysis (see Additional file 3 for details of taxon sampling). A multiple alignment of the amino acid sequences for each protein-coding gene was performed using ClustalX [65] with default options. Implementation of ClustalX for multiple sequence alignment does not always guarantee an unambiguous result due to the length and sequence variation among the species. This becomes more problematic when taxon sampling includes a wide range of animal taxa that evolve at very different rates. To improve reliability, we selected conserved blocks from aligned amino acid sequences for each of protein-coding genes using the Gblocks program [66] and a concatenated dataset was then prepared for the subsequent phylogenetic analyses. To reconstruct mitochondrial gene phylogeny, maximum likelihood (ML) and Bayesian inference (BI) were conducted for the concatenated amino acid sequence dataset. For ML analysis, the best-fit model of our amino acid sequence datasets was estimated using the Akaike Information Criterion (AIC) using ProtTest version 2.0 [67]. Maximum likelihood analysis was performed in Treefinder October version [68] using the MtArt matrix [69], which was selected as the best-fit model of amino acid substitution from ProtTest. Nodal support of the resulting ML trees was estimated by nonparametric bootstrap analysis with 500 random replications using Treefinder. We compared the likelihood scores of the competing hypotheses (the best tree versus alternative hypothesis) using various criteria (Expected-Likelihood Weight, ELW [70]; Shimodaira-Hasegawa [71]; Approximately Unbiased, AU [72] tests) implemented in Treefinder. The current version of the MrBayes program does not include the MtArt model of protein sequence, and thus we used the MtRev model with the likelihood parameter setting to "ngammacat = 4", "rates = invgamma" as an alternative best-fit model for Bayesian analysis. The analysis was run for 106 generations, sampled every 100 generations with four Markov Chain Monte Carlo (MCMC) chains using MrBayes 3.1.2 [73]. Bayesian posterior probability (BPP) values were estimated after the initial 200 saved trees (the first 2 × 105 generations) were discarded as burn-in. We also conducted maximum parsimony (MP) and neighbour-joining (NJ) analyses and nodal support in the resulting tree was estimated by nonparametric bootstrap analysis with 1,000 random replications in PAUP* 4.0b10 version [74].

Abbreviations

atp6, and atp8: genes for ATP synthase subunits 6 and 8; BI: Bayesian inference; bp: base pair; BI: Bayesian inference; BP: bootstrap percentage; BPP: Bayesian posterior probability; cob: gene for cytochrome oxidase b; cox1-cox3: genes for cytochrome oxidase c subunit 1-3; dNTP: deoxyribonucleotide triphosphate; EST: expressed sequence tag; kb: kilo base; LSU: large subunit nuclear ribosomal DNA; ML: maximum likelihood; MP: maximum parsimony; mtDNA: mitochondrial DNA; nad1-6: and nad4L: genes for NADH dehydrogenase subunits 1-6 and 4L; NCR: non-coding region; NJ: neighbor joining; nt: nucleotide; ntp: nucleotide pair; ORF: open reading frame; PCR: polymerase chain reaction; rrnS: and rrnL: genes for small and large mitochondrial ribosomal RNA subunits; SSU: small subunit nuclear ribosomal DNA; tRNA: transfer RNA.

Authors' contributions

GSM participated in the design of the study and carried out most of the molecular work. JKP designed the study, performed phylogenetic analyses, interpreted the results and drafted the manuscript. Both authors read and approved the final manuscript.

Acknowledgements

We are indebted to two anonymous reviewers and Steven Nadler for their insightful comments that helped to improve this manuscript. We also thank Mi-Hyun Park, Jina Baek (Inha University), Seokha Kang (Chungbuk National University) for their assistance in performing laboratory works. This work was supported by National Research Foundation of Korea Grant funded by the Korean Government (KRF-2008-313-C00813) and the research grant of the Chungbuk National University in 2008.

References

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

    Zootaxa 2007, 1564:1-104. OpenURL

  2. Wallace RL: Rotifera. In Encyclopedia of reproduction. Volume 4. Edited by Knobil E, Neil JD. San Diego: Academic Press; 1998::118-129. OpenURL

  3. Wallace RL, Snell TW, Ricci C: Rotifera: Volume 1: Biology, Ecology and Systematics. In Guides to the Identification of the Microinvertebrates of the Continental Waters of the World. Volume 23. Edited by Segers H. Ghent: Kenobi Productions; Leiden: Backhuys Publishers; 2006. OpenURL

  4. Herzig A: The analysis of planktonic rotifer populations: A plea for long-term investigations.

    Hydrobiologia 1987, 147:163-180. Publisher Full Text OpenURL

  5. Wallace RL, Snell TW: Phylum Rotifera. In Ecology and Classification of North American Freshwater Invertebrates. Edited by Thorp JH, Covich AP. New York: Academic Press; 1991:187-248. OpenURL

  6. Rublee PA: Rotifers in arctic North America with particular reference to their role in microplankton community structure and response to ecosystem perturbations in Alaskan Arctic LTER lakes.

    Hydrobiologia 1998, 387/388:153-160. Publisher Full Text OpenURL

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

    Hydrobiologia 1998, 387/388:101-107. Publisher Full Text OpenURL

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

    Invertebr Biol 2000, 119:17-26. OpenURL

  9. Mark Welch DB, Meselson M: Evidence for the evolution of bdelloid rotifers without sexual reproduction or genetic exchange.

    Science 2000, 288:1211-1215. PubMed Abstract | Publisher Full Text OpenURL

  10. Fontaneto D, Herniou EA, Boschetti C, Caprioli M, Melone G, Ricci C, Barraclough TG: Independently evolving species in asexual bdelloid rotifers.

    PLoS Biol 2007, 5(4):e87. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

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

    Zool Stud 2005, 44:34-43. OpenURL

  12. Ricci C, Melone G, Sotgia C: Old and new data on Seisonidea (Rotifera).

    Hydrobiologia 1993, 255/256:495-511. Publisher Full Text OpenURL

  13. Ahlrichs WH:

    Zur Ultrastruktur und Phylogenie von Seison nebaliae Grube,1859 und Seison annulatus Claus, 1876 - Hypothesen zu phylogenetischen Verwandtschaftsverhältnissen innerhalb der Bilateria. Cuvillier, Göttingen. 1995. OpenURL

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

    Zoomorphology 1997, 117:41-48. Publisher Full Text OpenURL

  15. 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: Oxford University Press; 1985:210-223. OpenURL

  16. Nielsen C, Scharff N, Eibye-Jacobsen D: Cladistic analyses of the animal kingdom.

    Biol J Linnean Soc 1996, 57:385-410. Publisher Full Text OpenURL

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

    Invertebr Biol 1996, 115:104-112. Publisher Full Text OpenURL

  18. Winnepenninckx B, Backeljau T, Mackey LY, Brooks JM, De Wachter R, Kumar S, Garey JR: 18S rRNA data indicate that the aschelminthes are polyphyletic in origin and consist of at least three distinct clades.

    Mol Biol Evol 1995, 12:1132-1137. PubMed Abstract | Publisher Full Text OpenURL

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

    J Mol Evol 1996, 43:287-292. PubMed Abstract | Publisher Full Text OpenURL

  20. Garey JR, Schmidt-Rhaesa A: The essential role of "minor" phyla in molecular studies of animal evolution.

    Amer Zool 1998, 38:907-917. OpenURL

  21. Giribet G, Distel DL, Polz M, Sterrer W, Wheeler WC: Triploblastic relationships with emphasis on the acoelomates and the position of Gnathostomulida, Cyclophora, Platyhelminthes, and Chaetognatha: A combined approach of 18S rDNA sequences and morphology.

    Syst Biol 2000, 49:539-562. PubMed Abstract | Publisher Full Text OpenURL

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

    Mol Evol Phylogenet 2006, 40:585-608. Publisher Full Text OpenURL

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

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

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

    Hydrobiologia 1998, 387/388:93-96. Publisher Full Text OpenURL

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

    Hydrobiologia 1998, 387/388:83-91. Publisher Full Text OpenURL

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

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

  27. Witek A, Herlyn H, Meyer A, Boell L, Bucher G, Hankeln T: EST based phylogenomics of Syndermata questions monophyly of Eurotatoria.

    BMC Evol Biol 2008, 8:345. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

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

    Cladistics 2004, 20:1-13. Publisher Full Text OpenURL

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

    Zoomorphology 1998, 118:255-261. Publisher Full Text OpenURL

  30. Wallace RL, Colburn RA: Phylogenetic relationships within phylum Rotifera: orders and genus Notholca.

    Hydrobiologia 1989, 186/187:311-318. Publisher Full Text OpenURL

  31. Sørensen MV: On the evolution and morphology of the rotiferan trophi, with a cladistic analysis of Rotifera.

    J Zool Syst Evol Res 2002, 40:129-154. Publisher Full Text OpenURL

  32. García-Varela M, Pérez-Ponce de León G, de la Torre P, Cummings MP, Sarma SS, Laclette JP: Phylogenetic Relationships of Acanthocephala Based on Analysis of 18S Ribosomal RNA Gene Sequences.

    J Mol Evol 2000, 50:532-540. PubMed Abstract | Publisher Full Text OpenURL

  33. Mark Welch DB: Bayesian and maximum likelihood analyses of rotifer-acnathocephalan relationships.

    Hydrobiologia 2005, 546:47-54. Publisher Full Text OpenURL

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

    Integr Comp Biol 2002, 42:668-677. Publisher Full Text OpenURL

  35. Wolstenholme DR: Animal mitochondrial DNA: structure and evolution.

    Int Rev Cytol 1992, 141:173-216. PubMed Abstract | Publisher Full Text OpenURL

  36. Boore JL: Animal mitochondrial genomes.

    Nucleic Acids Res 1999, 27:1767-1780. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Lavrov DV, Lang BF: Poriferan mtDNA and animal phylogeny based on mitochondrial gene arrangements.

    Syst Biol 2005, 54:651-659. PubMed Abstract | Publisher Full Text OpenURL

  38. Dellaporta SL, Xu A, Sagasser S, Jakob W, Moreno MA, Buss SW, Schierwater B: Mitochondrial genome of Trichoplax adhaerens supports Placozoa as the basal lower metazoan phylum.

    Proc Natl Acad Sci USA 2006, 103:8751-8756. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Littlewood DTJ, Lockyer AE, Webster BL, Johnston DA, Le TH: The complete mitochondrial genomes of Schistosoma haematobium and Schistosoma spindale and the evolutionary history of mitochondrial genome changes among parasitic flatworms.

    Mol Phylogenet Evol 2006, 39:452-467. PubMed Abstract | Publisher Full Text OpenURL

  40. Gissi C, Iannelli F, Pesole G: Evolution of the mitochondrial genome of metazoan as exemplified by comparison of congeneric species.

    Heredity 2008, 101:301-320. PubMed Abstract | Publisher Full Text OpenURL

  41. Shen X, Ma X, Ren J, Zhao F: A close phylogenetic relationship between Sipuncula and Annelida evidenced from the complete mitochondrial genome sequence of Phascolosoma esculenta.

    BMC Genomics 2009, 10:136. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  42. Steinauer ML, Nickol BB, Broughton R, Ortí G: First sequenced mitochondrial genome from the phylum Acanthocephala (Leptorhynchoides thecatus) and its phylogenetic position within metazoan.

    J Mol Evol 2005, 60:706-715. PubMed Abstract | Publisher Full Text OpenURL

  43. Suga K, Mark Welch DB, Tanaka Y, Sakakura Y, Hagiwara A: Two circular chromosomes of unequal copy number make up the mitochondrial genome of the rotifer Brachionus plicatilis.

    Mol Biol Evol 2008, 25:1129-1137. PubMed Abstract | Publisher Full Text OpenURL

  44. Lavrov DV, Brown WM: Trichinella spiralis mtDNA: a nematode mitochondrial genome that encodes a putative ATP8 and normally structured tRNAs and has a gene arrangement relatable to those of coelomate metazoans.

    Genetics 2001, 157:621-637. PubMed Abstract | PubMed Central Full Text OpenURL

  45. Boore JL, Brown WM: Mitochondrial genomes of Galathealinum, Helobdella, and Platynereis: sequence and gene arrangement comparisons indicate that Pogonophora is not a phylum and Annelida and Arthropoda are not sister taxa.

    Mol Biol Evol 2000, 17:87-106. PubMed Abstract | Publisher Full Text OpenURL

  46. Scouras A, Smith MJ: A novel mitochondrial gene order in the crinoid echinoderm Florometra serratissima.

    Mol Biol Evol 2001, 18:61-73. PubMed Abstract | Publisher Full Text OpenURL

  47. Herbeck JT, Novembre J: Codon usage patterns in cytochrome oxidase I across multiple insect Orders.

    J Mol Evol 2003, 56:691-701. PubMed Abstract | Publisher Full Text OpenURL

  48. Helfenbein KG, Brown WM, Boore JL: The complete mitochondrial genome of the articulate brachiopod Terebratalia transversa.

    Mol Biol Evol 2001, 18:1734-1744. PubMed Abstract | Publisher Full Text OpenURL

  49. Ojala D, Montoya J, Attardi G: tRNA punctuation model of RNA processing in human mitochondria.

    Nature 1981, 290:470-474. PubMed Abstract | Publisher Full Text OpenURL

  50. Lunt DH, Whipple LE, Hyman BC: Mitochondrial DNA variable number tandem repeats (VNTRs): utility and problems in molecular ecology.

    Mol Ecol 1998, 7:1441-1455. PubMed Abstract | Publisher Full Text OpenURL

  51. Cavalier-Smith T: A revised six-kingdom system of life.

    Biol Rev 1998, 73:203-266. PubMed Abstract | Publisher Full Text OpenURL

  52. Aguinaldo AM, Turbeville JM, Linford LS, Rivera MC, Garey JR, Raff RA, Lake JA: Evidence for a clade of nematodes, arthropods and other moulting animals.

    Nature 1997, 387:489-493. PubMed Abstract | Publisher Full Text OpenURL

  53. Philippe H, Lartillot N, Brinkmann H: Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia.

    Mol Biol Evol 2005, 22:1246-1253. PubMed Abstract | Publisher Full Text OpenURL

  54. Webster BL, Copley RR, Jenner RA, Mackenzie-Dodds JA, Bourlat SJ, Rota-Stabelli O, Littlewood DTJ, Telford MJ: Mitogenomics and phylogenomics reveal priapulid worms as extant models of the ancestral Ecdysozoan.

    Evol Dev 2006, 8:502-510. PubMed Abstract | Publisher Full Text OpenURL

  55. Littlewood DTJ, Telford MJ, Clough KA, Rohde K: Gnathostomulida--an enigmatic metazoan phylum from both morphological and molecular perspectives.

    Mol Phylogenet Evol 1998, 9:72-79. PubMed Abstract | Publisher Full Text OpenURL

  56. Jang KH, Hwang UW: Complete mitochondrial genome of Bugula neritina (Bryozoa, Gymnolaemata, Cheilostomata): phylogenetic position of Bryozoa and phylogeny of lophophorates within the Lophotrochozoa.

    BMC Genomics 2009, 10:167. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  57. Lavrov DV, Brown WM, Boore JL: Phylogenetic position of the Pentastomida and (pan)crustacean relationships.

    Proc R Soc Lond B 2004, 271:537-544. Publisher Full Text OpenURL

  58. Simison WB, Lindberg DR, Boore JL: Rolling circle amplification of metazoan mitochondrial genomes.

    Mol Phylogenet Evol 2006, 39:562-567. PubMed Abstract | Publisher Full Text OpenURL

  59. Wang X, Wang J, He S, Mayden RL: The complete mitochondrial genome of the Chinese hook snout carp Opsariichthys bidens (Actinopterygii: Cypriniformes) and an alternative pattern of mitogenomic evolution in vertebrate.

    Gene 2007, 399:11-19. PubMed Abstract | Publisher Full Text OpenURL

  60. Mueller RL, Macey JR, Jaekel M, Wake DB, Boore JL: Morphological homoplasy, life history evolution, and historical biogeography of plethodontid salamanders inferred from complete mitochondrial genomes.

    Proc Natl Acad Sci USA 2004, 101:13820-13825. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  61. Mark Welch JL, Mark Welch DB, Meselson M: Cytogenetic evidence for asexual evolution of bdelloid rotifers.

    Proc Natl Acad Sci USA 2004, 101:1618-1621. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Mark Welch DB, Cummings MP, Hillis DM, Meselson M: Divergent gene copies in the asexual class Bdelloidea (Rotifera) separated before the bdelloidean radiation or within bdelloid families.

    Proc Natl Acad Sci USA 2004, 101:1622-1625. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Wuethrich B: Evolution of sex: The asexual life.

    Science 1998, 281:1981. Publisher Full Text OpenURL

  64. Wyman SK, Jansen RK, Boore JL: Automatic annotation of organellar genomes with DOGMA.

    Bioinformatics 2004, 20:3252-3255. PubMed Abstract | Publisher Full Text OpenURL

  65. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.

    Nucleic Acids Res 1997, 24:4876-4883. Publisher Full Text OpenURL

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

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

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

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

  68. Jobb G: TREEFINDER.

    version of October 2008. Munich, Germany 2008.

    Distributed by the author at http://www.treefinder.de webcite

    OpenURL

  69. Abascal F, Zardoya R, Posada D: MtArt: A New Model of Amino Acid Replacement for Arthropoda.

    Mol Evol Biol 2007, 24:1-5. Publisher Full Text OpenURL

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

    Proc R Soc Lond B 2002, 269:137-142. Publisher Full Text OpenURL

  71. Shimodaira H, Hasagawa M: Multiple comparisons of log-likelihoods with applications to phylogenetic inference.

    Mol Biol Evol 1999, 16:1114-1116. OpenURL

  72. Shimodaira H: An approximately unbiased test of phylogenetic tree selection.

    Syst Biol 2002, 51:492-508. PubMed Abstract | Publisher Full Text OpenURL

  73. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models.

    Bioinformatics 2003, 19:1572-1574. PubMed Abstract | Publisher Full Text OpenURL

  74. Swofford DL: PAUP*: Phylogenetic analysis using parsimony (*and other methods) version 4.0b10. Sunderland, Massachusetts: Sinauer Associates; 2002. PubMed Abstract | Publisher Full Text OpenURL

  75. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates.

    Mol Mar Biol Biotechnol 1994, 3:294-299. PubMed Abstract OpenURL

  76. Xiong B, Kocher TD: Comparison of mitochondrial DNA sequences of seven morphospecies of black flies (Diptera: Simuliidae).

    Genome 1991, 34:306-311. PubMed Abstract OpenURL

  77. Telford MJ, Holland PWH: The phylogenetic affinities of the chaetognaths: A molecular analysis.

    Mol Biol Evol 1993, 10:660-676. PubMed Abstract | Publisher Full Text OpenURL

  78. Wada H, Satoh N: Details of the evolutionary history from invertebrates to vertebrates, as deduced from the sequences of 18S rDNA.

    Proc Natl Acad Sci USA 1994, 91:1801-1804. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  79. Helfenbein KG, Fourcade HM, Vanjani RG, Boore JL: The mitochondrial genome of Paraspadella gotoi is highly reduced and reveals that chaetognaths are a sister group to protostomes.

    Proc Natl Acad Sci USA 2004, 101:10639-10643. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  80. Papillon D, Perez Y, Caubit X, Le Parco Y: Identification of chaetognaths as protostomes is supported by the analysis of their mitochondrial genome.

    Mol Biol Evol 2004, 21:2122-2129. PubMed Abstract | Publisher Full Text OpenURL

  81. Halanych KM: The new view of animal phylogeny.

    Annu Rev Ecol Syst 2004, 35:229-256. Publisher Full Text OpenURL