Scholarly article on topic 'Mitochondrial genome sequence and gene order of Sipunculus nudus give additional support for an inclusion of Sipuncula into Annelida'

Mitochondrial genome sequence and gene order of Sipunculus nudus give additional support for an inclusion of Sipuncula into Annelida Academic research paper on "Biological sciences"

Share paper
Academic journal
BMC Genomics
OECD Field of science

Academic research paper on topic "Mitochondrial genome sequence and gene order of Sipunculus nudus give additional support for an inclusion of Sipuncula into Annelida"

BMC Genomics KoMe^ t

Open Access

Research article

Mitochondrial genome sequence and gene order of Sipunculus nudus give additional support for an inclusion of Sipuncula into Annelida

Adina Mwinyi1, Achim Meyer2, Christoph Bleidorn3, Bernhard Lieb2, Thomas Bartolomaeus4 and Lars Podsiadlowski*1'4

Address: Institut für Zoologie, Freie Universität Berlin, Koenigin-Luise-Str. 1-3, 14195 Berlin, Germany, 2Institut für Zoologie, Johannes Gutenberg Universität Mainz, Müllerweg 6, 55099 Mainz, Germany, 3Institut für Biochemie und Biologie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, 14476 Golm, Germany and 4Institut für Evoilutionsbiologie & Ökologie, Universität Bonn, An der Immenburg 1, 53121 Bonn, Germany

Email: Adina Mwinyi -; Achim Meyer -; Christoph Bleidorn -; Bernhard Lieb -; Thomas Bartolomaeus -; Lars Podsiadlowski* - * Corresponding author

Published: 16 January 2009 Received: 3 September 2008

BMC Genomics 2009, 10:27 doi:l0.ll86/l47l-2l64-l0-27 Accepted: 16 January 2009

This article is available from: © 2009 Mwinyi et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Background: Mitochondrial genomes are a valuable source of data for analysing phylogenetic relationships. Besides sequence information, mitochondrial gene order may add phylogenetically useful information, too. Sipuncula are unsegmented marine worms, traditionally placed in their own phylum. Recent molecular and morphological findings suggest a close affinity to the segmented Annelida.

Results: The first complete mitochondrial genome of a member of Sipuncula, Sipunculus nudus, is presented. All 37 genes characteristic for metazoan mtDNA were detected and are encoded on the same strand. The mitochondrial gene order (protein-coding and ribosomal RNA genes) resembles that of annelids, but shows several derivations so far found only in Sipuncula. Sequence based phylogenetic analysis of mitochondrial protein-coding genes results in significant bootstrap support for Annelida sensu lato, combining Annelida together with Sipuncula, Echiura, Pogonophora and Myzostomida.

Conclusion: The mitochondrial sequence data support a close relationship of Annelida and Sipuncula. Also the most parsimonious explanation of changes in gene order favours a derivation from the annelid gene order. These results complement findings from recent phylogenetic analyses of nuclear encoded genes as well as a report of a segmental neural patterning in Sipuncula.


Molecular sequence analysis has become the method of choice to address phylogenetic questions. The applied techniques improve continually and the rapidly growing amount of available data helps to broaden our knowledge of phylogenetic relationships within the animal kingdom.

Nevertheless, different molecular datasets often show conflicting phylogenetic signals, so that results relying on just one dataset may be interpreted with caution [1]. Unlike nuclear DNA, the mt-genome of animals is normally rather small and simply structured: haploid, without or only few non-coding segments, repetitive regions

and transposable elements. Derived from endosymbiotic bacteria only a few genes are retained in the mitochondrial genomes of Bilateria: 13 protein subunits (nad1-6, nad4L, cox1-3, cob, atp6/8), 2 ribosomal RNAs (rrnL, rrnS) and 22 tRNAs are found encoded on a circular dou-blestranded DNA molecule sized about 15 kb [2,3]. As such sequencing and annotation of mt-genomes is much easier and faster than analysing nuclear genomes, making mt-genomes one of the commonly used sources of sequence data for phylogenetic analyses. Apart from sequence data other features of the genome may contain phylogenetic information, too. Taxon-specific gene order often remains identical over long periods of time [4-6]. Simultaneously, the intra-taxonomic variances of these characteristic orders are quite distinctive and convergent changes in the positioning of single genes are rather unlikely, due to the vast number of possible combinations [7]. Thus changes in the mitochondrial gene order have proved to be valuable tools in phylogenetic analyses [8-10]. Less often secondary structures of tRNAs or rRNAs show distinct differences between taxa (e.g. loss of a stem/ loop region) and hence may also contribute to a phylogenetic analysis [11].

The taxon Sipuncula (peanut worms) comprises about 150 species, being found in all water depths of different marine habitats. The hemisessile organisms dwell in mud and sand, but settle also in empty mollusc shells or coral reef clefts for instance. Their body shows no segmentation, but a subdivision into a posterior trunk and an anterior introvert that can be fully retracted into the trunk is observeable [12]. Fossils that date back into the later cambrian [13] suggest that sipunculans have undergone little morphologically change over the past 520 Myr. The monophyly of this morphologically uniform taxon is well founded by morphological [14] and molecular data [15]. However, the phylogenetic position within Bilateria was highly disputed. Based on morphological characters, very different phylogenetic positions of Sipuncula were discussed. Early in history an affinity to Echinodermata, especially holothurians was mentioned and later again propagated by Nichols [16], but with little acceptance from other authors. Scheltema [17] proposed a close relationship to molluscs based on the presence of the so calles "molluscan cross" organization of micromeres during spiral cleavage. The usefulness of this character for phylogenetic inference was neglected by Malaskova [18]. Other analyses found Sipuncula to be sister group of Mollusca, Annelida and Arthropoda [19], Articulata (Annelida and Arthropoda) [14], Echiura [20], Mollusca [21], Annelida [22] or Annelida+Echiura [23]. More details about the different hypotheses of sipunculid relationships are reviewed in [24].

In contrast to all these studies, molecular analyses of large datasets from 18S/28S data [25], ESTs [26,27] or mito-

chondrial genome data [28,29] favour an inclusion of Sipuncula into annelids. An implication of this hypothesis is that we have to assume that segmentation has been reduced within Sipuncula [30]. A derivation from segmented ancestors of Sipuncula was recently also supported by a segmental mode of neural patterning in ontogeny [31].

Relationhips within Sipuncula are well investigated [15,24,32-34]. An analysis using combined molecular and morphological data recovered five major clades and supports that Sipunculus is the sister group to all other sipunculids [15].

Up to now mt-genome data from Sipuncula was restricted to a partial mtDNA sequence from Phascolosoma gouldii [29], comprising only about half of the complete genome. Here we describe the first complete mitochondrial genome for another representative of the Sipuncula, Sipunculus nudus. We analyse sequence data in comparison with mitochondrial genomes of various Bilateria to evaluate the phylogenetic position of Sipuncula. In addition we compare gene order among Lophotrochozoa and evaluate the most parsimonious explanation for gene order changes.

Results and discussion Genome organisation

The complete mt-genome of S. nudus is a circular DNA doublestrand of 15502 bp length. As usual in bilateria, 13 genes coding for different protein subunits and two encoding ribosomal RNA genes were identified. In addition 22 tRNA genes were detected and thus all 37 genes typically present in bilaterian mt genomes, were found (Fig. 1, Table 1). All of these genes are located on the (+)-strand, as is the case in annelid and echiurid mt-genomes. There are two small gene overlaps: one between nad4L and nad4 (7 bp), the other one between trnS (AGN) and nad2 (1 bp). The putative control region is 441 bp in length and flanked by trnF and trnT. Besides the control region 15 other non-coding regions are dispersed over the whole genome, ranging from one to 39 base pairs. The three largest of these are located between trnY and trnE (35 bp), trnH and nad5 (39 bp) and nad5 and trnS (AGN) (21 bp).

The GC-skew [(G-C)/(G+C)] reflects the relative number of cytosin to guanine and is often used to describe the strand-specific bias of the nucleotide composition [35]. In S.nudus the complete (+)-strand genome sequence has a clear bias toward Cytosine (GC-skew -0.296). As all genes are coded on (+)-strand, all single gene sequences exhibit a negative GC-skew, too (Table 1), ranging from -0.23 (rrnS) to -0.45 (nad2). A negative GC-skew is also found in most of the mitochondrial genomes known from annelids, pogonophorans, and myzostomids, with the excep-

Table 1: Genome organisation of Sipunculus nudus. Complete circular mtDNA has a lenght of 15502 bp.

Gene Strand Position Length GC-IAT- Start- Stop- Intergenic bp

(start - end) (nuc.) skew codon codon

cox1 + 1 - 1543 1543 -0.24/-0.07 ATG TAA I2

trnN + 1556 - 1624 69 0

cox2 + 1625 - 2319 695 -0.26/-0.07 ATG TA 0

trnD + 2320 - 2385 66 0

atp8 + 2386 - 2544 159 -0.38/0.12 ATG TAG 2

trnY + 2547 - 2609 63 36

trnE + 2646 - 2712 67 I

trnG + 2714 - 2780 67 0

cox3 + 2781 - 3560 780 -0.27/-0.07 ATG TAA 4

trnQ + 3565 - 3632 68 0

nad6 + 3633 - 4106 474 -0.34/-0.I8 ATG TAG I

cob + 4108 - 5247 1 140 -0.30/-0.06 ATG TAA 7

trnP + 5255 - 5322 68 0

trnS-UCN + 5323 - 5389 67 5

trnC + 5395 - 5455 61 5

trnM + 5461 - 5527 67 0

rrnS (12S) + 5528 - 6373 846 -0.23/0.18 0

trnV + 6374 - 6442 69 0

rrnL (16S) + 6443 - 7929 1487 -0.26/0.08 0

trnL-CUN + 7930 - 7995 66 7

trnA + 8003 - 8070 68 0

trnI + 8071 - 8139 69 0

trnK + 8140 - 8207 68 0

nad3 + 8208 - 8565 358 -0.35/-0.05 ATG T 2

trnF + 8568 - 8631 64 0

Major NCR + 8632 - 9072 441 -0.14/0.09 0

trnT + 9073 - 9142 70 0

nad4L + 9143 - 9424 282 -0.43/-0.06 ATG TAA -7

nad4 + 9418 - 10774 1357 -0.37/-0.02 ATG T 7

trnL-UUR + 10782 - 10846 65 0

nad1 + 10847 - 11789 943 -0.32/-0.09 ATG T 0

trnW + 11790 - 1 1855 66 7

atp6 + 11863 - ■ 12550 688 -0.41/-0.14 ATG T 0

trnR + 12551 - 12619 69 I

trnH + 12689 - 12688 68 39

nad5 + 12728 - 12727 1698 -0.37/0.01 ATA TAA 2I

trnS-AGN + 14447 - 14518 72 -I

nad2 + 14518 - ■ 15502 985 -0.45/-0.13 ATG T 0

* start and stop position of ribosomal RNA and NCR according to adjacent gene boundaries

tion of the annelid Eclysippe vanelli [36]. AT-skew of the complete (+)-strand is close to evenness (-0.013) and single gene AT-skews are distributed around evenness with a range between 0.18 (rrnS) and -0.18 (nad6), see also Table 1. AT content of the complete genome is 54.2%, AT contents of protein-coding and rRNA genes are not much derived from this value, between a minimum of 50,3% (nad3) and a maximum of 59,8% (atp8).

Protein coding genes

All but one of the protein subunits begin with start codon ATG, only nad5 starts with ATA. Both are prevalent in mitochondrial genomes. The commonly found stop codons TAA and TAG are present, as well as the abbreviated forms TA (cox2) and T (nad1-4, atp6). Putative short-

ened stop codons were already found in other species and are thought to be complemented via post-transcriptional polyadenylation [37].

Ribosomal RNA genes and control region

The sizes of the ribosomal RNAs (rrnS: 846 bp; rrnL: 1487 bp) are within the range of their sizes in other animals including molluscs and annelids. The two genes only separated by trnV, a feature often found in animals from vertebrates to arthropods, so therefore this represent an ancestral condition. Among annelids and their kin only echiurans (Urechis caupo) and myzostomids (Myzostoma seymourcollegiorum) differ from that condition in that there is no tRNA gene separating the two ribosomal genes. AT content of ribosomal genes is 50.8% (rrnS) and 53.1%

K I A Ln


Circular map of the mitochondrial genome of Sipun-culus nudus.

found between the poly-T stretch and the GA rich region in arthropods.

Transfer RNAs

All typical 22 tRNAs were detected in the mitochondrial genome of S. nudus, their putative secondary structures are depicted in Fig. 2 and Additional file 1. All but three tRNA genes are capable to be folded in the usual cloverleaf structure, consisting of T^C stem and loop, anticodon stem and loop, DHU stem and loop, and the acceptor stem -tRNA-Ser(AGN) and tRNA-Ser(uCN) have no DHU stem. While tRNA-Ser(AGN) shows this feature in many bilate-rian mt-genomes, the other one must have changed its secondary structure in the lineage leading to Sipuncula and after the split of its sister group. The putative secondary structure of tRNA-Cys shows no T^C, in addition there are two mismatches in the anticodon stem and an unusual anticodon (ACA), weakening this secondary structure hypothesis. But intensive search for an alternative sequence of tRNA-Cys was not successful, so we stuck with this hypothesis although we cannot rule out that this is a non-functional sequence or subject to gene editing. In several other tRNAs there are mismatches in the acceptor or anticodon stem.

(rrnL), so well within the range of AT content of protein-coding genes.

Noncoding regions, putative control region

The putative control region is found between nad3/trnF on one side and trnT/nad4L/nad4 on the other side. While gene order (or protein-coding and rRNA genes) in Annelida is more or less conserved there is a great variation in the position of the control region: (a) Species from Clitel-lata, Maldanidae and Terebellidae have a major non-coding region between atp6/trnR and trnH/nad5; (b) in Orbinia it is located between nad4/trnC/trnL2 and trnLl/ trnM/rrnS; (c) in Platynereis it is found between cox2/trnG and trnV/atp8 [8,28,36,38]. Such great variability is not found in other taxa like Arthropoda or Vertebrata, where also the control region is found in the same position in different species, when gene order of the rest of the mt-genome is conserved.

In Sipunculus nudus the major non-coding region has a size of 441 bp and is clearly more AT rich (66.1%) than the rest of the genome (53.9%). Structural elements know from arthropod mitochondrial control regions [39] are present also in S. nudus: (1) a poly-TA(A) stretch of 50 bp including a tenfold TA repeat; (2) a poly-T stretch flanked by purin bases; (3) a GA-rich block of 16 bases length. Although we examined the complete non-coding region intensively by software and by eye, no large stem-loop structure was identified. Such a structure is normally

Mitochondrial gene order

Fig. 3 shows a comparison of lophotrochozoan mito-chondrial gene orders and the ground pattern of Bilateria (as mentioned in [40]). We restrict the discussion of gene order to the protein-coding and rRNA genes, as tRNA genes change their relative position much faster than the former, as seen in gene order comparisons of e.g. annelids [8] or crustaceans [41]. The annelids, pogonophorans and myzostomids do not differ from each other in the relative positions of protein-coding and rRNA genes. Compared to the ground pattern of Bilateria several genes have a different relative position: (1) nad6/cob are found right after cox3, (2) atp6 is found between cob and nad5, (3) nad5 and nad4L/nad4 have interchanged positions, and (4) nad3 is found between nadl and nad2 (numbers refer also to hypothesized events in Fig. 4). Mollusca (Conus textile [42], Ilyanassa obsoleta [43]) and Brachiopoda (Terebratu-lina retusa [44]) show a different pattern, with derived positions for three gene blocks: rrnS/rrnL/nadl , cox3/nad3 and nad6/cob. The translocation of nad6/cob may be explained as a commonly derived feature of Lophotrocho-zoa, or a subtaxon of it including Mollusca, Phoronida, Brachiopoda, Nemertea, Annelida s. l. (including Pogono-phora, Echiura and Myzostomida) and Sipuncula (compare Fig. 4). The other translocation events found in annelids and their kin (2.-4.) seem to be restricted to that group. The gene order so far known from Nemertea (Cephalothrix rufifrons, partial genome [45]) can be easily derived with one change (translocation of nad6) from the pattern of the brachiopod Terebratulina and the gene order


Secondary structure of tRNAs identified in the mitochondrial genome of S. nudus. The best found putative secondary structure of tRNA-Cys (box) seems to be strongly derived, probably non-functional or subject to gene editing. This figure shows only part of the full image, please see also Additional file I.

Phoronida (Phoronis psammophila)

(+) (-)

coxl cox2

nad5 nad4 q! cob nad6 nadl rrnL rrnS

cox3 atp6 nad3 nad2

Some molluscs (Polyplacophora: Katharina; Gastropoda; Haliotis; Cephalopoda; Octopus)

(+) (-)

coxl cox2 atp6 cox3 nad3 nad2

nad5 nad4 ü! cob nad6 nadl rrnL rrnS

Cephalothrix rufifrons (Nemertea)

(+) -\c°xl\-

rrnL nadl cob nad4 nad5 nad6 cox3 nad3 nad2

Terebratulina retusa (Brachiopoda), some molluscs (Gastropoda: Conus, Ilyanassa, Thais, Lophiotoma)

(+) (-)

Annelida sensu strictu ("Polychaeta", "Oligochaeta", Hirudinea)

(+) coxl cox2 cox3 nad6 cob atp6 nad5 Í nad4 rrnS rrnL nadl nad3 nad2

Riftia pachyptila (Pogonophora: Vestimentifera)

(+) coxl cox2 cox3 nad6 cob atp6 nad4 rrnS rrnL nadl nad3 nad2

Galathealinum brachiosum (Pogonophora: Perviata)

(+) (-)

coxl cox2 cox3 nad6 cob rrnL nadl nad3 nad2

Myzostoma seymourcollegiorum (Myzostomida)

(+) (-)


Sipunculus nudus (Sipuncula)

(+) (-)

Phascolopsis gouldii (Sipuncula)

(+) (-)

Urechis caupo (Echiura)

(+) (-)


Comparison of mitochondrial gene order (protein-coding genes and ribosomal RNAs only) of several lopho-trochozoan taxa compared and the putative bilaterian ground pattern (according to [40]). Genome segments from the bilaterian ground pattern are colour coded for a better visualization of differences between gene orders. For complete species names and accession numbers see Table 3.

nad6 cob

nad6 cob

-■— Phoronida (Phoronis*)

-Mollusca (Katharina, Octopus)

■-Nemertea (Cephalothrix*)

-Brachiopoda (Terebratulina)

Annelida s. str.

-Pogonophora (Riftia, Galathealinum)

-Myzostomida (Myzostoma)

i— - - - Echiura (Urechis) 8 9

Sipuncula (Sipunculus, Phascolopsis)

1 2 3 4 5 6 7 8 9

Species nad 6+ atp6 nad5 nad3 COX1- nad2 rrnS- rmS+r aip6+

cob atp8 rrnL rnL nad 5

Platynereis dumerilü „Polychaeta" X X X X O 0 0 0 0

Orbinia latreilii „Polychaeta" X X X X 0 0 0 0 0

Clymenella torquata „Polychaeta" X X X X 0 0 0 0 0

Lumbricus terrestris „Oligochaeta" X X X X 0 0 0 0 0

Helobdella robusta* Hirudinea X X 7 X 0 0 0 ? 0

Galathealinum br. * Perviata X 7 7 X 0 0 0 ? ?

Riftia pachyptila* Vestimentifera X X ? X 0 0 0 0 ?

Urechis caupo Echiura X X X X X X X 0 0

Phascolopsis gouldii* Sipuncula X (X) ? (X) 0 0 0 ? X

Sipunculus nudus Sipuncula X (X) (x) 0 0 0 0 X X

Myzostomum seym. * Myzostomida X X X ? 0 0 0 0 0

*: partial genome data


Cladogram for changes in gene order of lophotrochozoan taxa (only changes in protein-coding and rRNA genes were analysed). The translocation of a gene or a gene block is treated as an apomorphic feature (small box) with numbers according to translocated genes in the table below. "x" indicates derived gene positions, circles stand for an unvaried order. "(x)" symbolizes that although the position of the gene is now different there is evidence that it. Questionmarks indicate missing sequence data or putative secondary events complicating the interpretation. Changes not mentioned in the table: (10) translocation of cox3/nad3; (II) translocation of rrnS/rmUnadi; (12) translocation of nad6; (I 3) large inversion of a segment spanning from rrnS to nad5; (I4) translocation of atp6. See text for further details.

of Phoronida (Phoronis psammophila, partial genome [46]) from that of the mollusc Katharina tunicata with only one event (translocation of atp6). Much more variation is seen within Mollusca [6,47] and Brachiopoda [48-50] (not shown).

Compared to the Annelida and their kin, the mitochondrial gene order of Sipunculus nudus differs clearly: (a) atp6 and nad5 are found between nadl and nad2. This may be interpreted as two events restricted to the sipunculid lineage and independently achieved from the bilaterian or lophotrochozoan ground pattern. But another explanation would be a singular event translocating the block atp6/nad5 compared to the annelid ground pattern (No. 8 in Fig. 4); (b) rrnS/rrnL found a different position, between cob and nad3 - this is as well different from the situation in Brachiopoda and Mollusca, so probably another event in the lineage leading to Sipuncula (No. 9 in Fig. 4); (c) nad3 is found right after rrnL and adjacent to nad4L/nad4. This is different from its position in annelids, pogonophorans, myzostomids and echiuran taxa and is more similar to the bilaterian ground pattern. Visualized in Fig. 4 the most parsimonious explanation of sipunculid gene order is that Sipuncula share two events with annelids, but lack the translocation of nad3. In addition two events have to be assumed in the lineage of Sipunula (rrnS/rrnL and atp6/nad5, corresponding to 8 and 9 in Fig. 4). Derivation of the Sipunculus gene order directly from the bilaterian ground pattern would demand four translocation events (nad6/cob, rrnS/rrnL, atp6, nad5) from which only one is shared with other lophotrochozoan taxa (nad6/cob). So this hypothesis is in demand of three additional events instead of two for the "annelid" hypothesis. Derivation of the sipunculid gene order from the brachio-pod/mollusc pattern is in demand of five additional events. Therefore the most parsimonious explanation of gene order changes would be that Sipuncula is sister group to a group comprising Annelida s.str., Myzostomida, Echiura and Pogonophora.

At first sight gene order of the echiurid Urechis caupo [51] is completely different from that of annelids and Sipunculus, but the position of atp6 between cob and nad5 and that of nad3 adjacent to nadl clearly hint to the derived features postulated for the annelid ground pattern (see b and c in the discussion of annelid gene order above). As well adjacency of nad6 to cox3 is found in all annelids and Sipu-nuculus. So the gene order of Urechis may be derived from the annelid ground pattern, with additional translocations of three genome segments: (a) cox1/cox2/atp8, (b) rrnS/rrnL and (c) nad2.

Phylogenetic analysis of mitochondrial sequences

The phylogenetic analysis was performed with a concatenated amino acid alignment of 11 protein-coding genes

(exept atp8 and nad4L) from 74 species. Fig. 5 shows the best tree of the Maximum Likelihood analysis with RaxML (mtREV+G+I). A close relationship of Sipunculus and Phas-colopsis and thus monophyletic Sipuncula is well supported (ML bt: 100%). Sipuncula appears to be close related to the classic "Annelida", Echiura and Pogono-phora - this assemblage has a bootstrap support of 93%. This assemblage is also other recovered in recent molecular analyses of 18S/28S rRNA and EFla [25] or EST data [26]. The internal relationships of these taxa are not well resolved by our analysis. With high bootstrap support Clitellata (98%) and Pogonophora (100%) appear mono-phyletic, while their sister group relationship found only weak support (bootstrap: 75%). Sister group to the Sipun-cula/Annelida/Echiura/Pogonophora taxon is Myzostom-ida (ML bt: 91%), this relationship is also supported by morphological characters and mitochondrial gene order as recently detailed elsewhere [8]. The position of this "Annelida sensu lato" among other Lophotrochozan subtaxa is not well resolved in our analysis.

Probably due to long branch effects, Ecdysozoa and Lophotrochoza appear not to be monophyletic in our analysis. While the former miss Nematoda, the latter miss Platyhelminthes, Ectoprocta, Rotifera, Acanthocephala and some molluscs. All these taxa are associated with long branches and form a probably artificial clade, which was never recovered in analyses with molecular data from nuclear genes or morphological data. Apart from this the most "problematic" taxon are Mollusca, with some taxa (Lottia, Argopecten, Venerupis, Siphonodentalium) found clustering with the above mentioned nematode-platy-helminth assemblage, others (Katharina, Haliotis, Nautilus, Octopus) clustering with Nemertea, Phoronida and Ento-procta, while Lampsilis appears as sister taxon to Brachiop-oda.

For further evaluation the interrelationships of Annelida sensu lato, we performed additional phylogenetic analyses with a smaller taxon set comprising 30 species (all species from the lophotrochozoan branch of the larger taxon set). ML analyses were done comparing mtREV (RaxML) and mtART (Treefinder) models; in addition a Bayesian analysis was performed with mtREV model (MrBayes). Myzostomida, Sipuncula and other Annelida formed a monophyletic group (Fig. 6) supported by ML bootstrapping (mtREV: 92%, mtART: 98%), but not by BI, where support is below 0.95 (Bayesian posterior probabilities). Sipuncula and Annelida together form a clade well supported by all three analyses, while Annelida without Sipuncula found best support only in BI, while the ML analyses do not significantly support this group, leaving open if there is a basal split between Sipuncula and the rest of the annelids. In the best ML-mtART tree Platynereis is found as sister to Sipuncula tree, but with bootstrap


Best tree from the Maximum Likelihood analysis, inferred from the mitochondrial amino acid data set of 11 protein coding genes (RaxML 7.00, mtREV, G+I, single gene partitions). Numbers beneath nodes are ML bootstrap percentages, bold branches indicate bootstrap percentages >85%. See Table 3 for complete species names and accession numbers. Asterisks indicate taxa with incomplete mt-genome information.


Best tree from the Maximum Likelihood analysis (RAxML 7.00, mtREV, G+I, single gene partitions) of the reduced taxon set (30 lophotrochozoan species). Numbers beneath nodes indicate support (from left to right or up to down, respectively): (1) through RaxML bootstrapping (1000 pseudoreplicates) (2) ML analysis with Treefinder (1000 pseudor-eplicates), model mtART+G+I, (3) Bayesian posterior probabilities (model mtREV+G+I). Triple asterisks indicate maximum support from all three analyses (100/100/1.0). See Table 3 for complete species names and accession numbers. Single asterisks indicate taxa with incomplete mt-genome information. Scalebar depicts substitutions per site in the best RAxML tree.

support below 50%. Well supported subtaxa of annelids are Pogonophora (s.lato), Clitellata, Pogonophora+Clit-ellata, Orbiniidae (Scoloplos+Orbinia). Topologies obtained in the three analyses differ in the position of Urechis (Echiura), which is found as sister to Maldanidae+Ter-ebelliformia in the best ML tree with mtREV model (bootstrap support 65%), as sister to Orbiniidae in the best tree with mtART model (bootstrap support below 50%) and as sister to Pogonophora+Clitellata in BI (BPP below 0.95).

In addition we performed an AU test as implemented in CONSEL to statistically test the hypothesis of a sister group relationship between Sipuncula and Mollusca. We were able to significantly reject (p < 0.001) this hypothesis compared to the best ML-tree (mtREV).


Annelida, in traditional phylogenetic systems the sister group to Arthropoda, are nowadays included in the taxon Lophotrochozoa by almost all large scale analyses [26,27,52-54]. In this view more and more molecular studies no longer support the monophyly of the classical Annelida ("polychaetes" and clitellates). As well as the unsegmented Pogonophora, Echiura, and Myzostomida the Sipuncula have also been under suspect to be included in what was called Annelida sensu lato [8,25,26,28]. The complete mitochondrial genomic sequence of Sipunculus nudus presented in this paper, adds an important piece of evidence to answer the question of sipunculid position in the metazoan tree of life. Our sequence data and gene order analysis clearly support an affinity of Sipuncula to

Annelida s. l. (including Pogonophora, Echiura and Myzostomida) rather than to Mollusca or any other phylum. It still remains an open question if Sipuncula and the whole Annelida s. l. are sister groups (as the most parsimonious explanation of gene order data suggests), or if Myzostomids form the sister group to Sipuncula and the remaining Annelida (as sequence based analyses favour). In sequence-based analyses the myzostomid is the annelid taxon with the longest branch, suggesting a more rapid evolution of mitochondrial sequence in this taxon. Therefore analyses placing Myzostomids outside the Annelida are probably misleading due to higher substitution rates in myzostomids.


Animals, DNA purification

A specimen of S. nudus was collected in Concarneau, France and conserved in 100% ethanol. Using the DNeasy® Blood & Tissue kit (Qiagen, Hilden, Germany) we followed the instructions given to extract DNA from animal tissues and used approximately 1 x 1 cm of the body wall from one individual.

PCR and purification of DNA fragments

EST sequence fragments for the genes nadl, nad3, rrnL, cob, coxl, cox2 and cox3 were used to design the first species specific primer pairs [27]. The complete mitochondrial genome of S. nudus was amplified in PCR fragments generated with species specific primer pairs from EST information (see Table 2). All PCRs were done with Eppendorf Mastercycler or Eppendorf Mastercycler Gradient thermocyclers. PCRs were carried out in 50 |l volumes

Table 2: Primer pairs and corresponding annealing temperatures used for successful amplification of mitochondrial genome fragments

from Sipunculus nudus

Primer names Primer sequence (5'-3') Annealing temperature Approxim. size of PCR product








(41.75 |l water, 5 |l 10x buffer, 0.25 |l Taq polymerase (5 U/|l), 1 ml dNTP mixture, 1 |l template DNA, 1 |l primer mixture (10 |M each)) using the Eppendorf 5-prime kit (Eppendorf, Germany). The cycling conditions were as follows: 94°C for 2 min (initial denaturation); 40 cycles of 94°C for 30 sec (denaturation); primer-specific temperature (see Table 2) for 1 min (annealing), 68°C for 1 min (elongation), was followed by 68 °C for 2 min (final elongation). After 40 cycles the samples were stored at 4°C and visualised on a 1% ethidium bromide-stained TBE agarose gel, respectively. DNA fragments expected to be larger than 3 kb, were amplified in 25 | l volumes (16.75 |l water, 2.5 |l buffer, 0.25 |l Takara LA Taq polymerase, 4 |l dNTP mixture, 1 |l template DNA, 0.5 |l primer mixture (10 |M each)) under the following long PCR conditions (Takara LA kit): 94°C for 2 min (initial denaturation); 40 cycles of 94°C for 30 sec (denatura-tion), primer-specific temperature for 1 min (annealing) and 72°C for 10 min (elongation). After the final elongation step (68°C for 2 min), samples were treated as described above. PCR products were purified with minispin columns provided in the Nucleo Spin Extract II kit (Macherey & Nagel) and the Blue Matrix PCR/DNA clean up DNA Purification kit (EurX, Gdansk, Poland). Dependent on the band intensity on the agarose gel, DNA was eluted in 30-60 |l elution buffer and stored at -20 °C. Slightly contaminated samples were cut from a 1% ethid-ium bromide-stained TAE agarose gel and purified with the QIAquick Gel Extraction kit (Qiagen) afterwards.


If the DNA amount, obtained by PCR, turned out to be insufficient for sequencing, the respective fragment was cloned in a pGEM-T Easy Vector (Promega). Ligation was carried out in 5 |l volumes instead of the double amount, proposed in the protocol. In each case 2 |l of the sample were used for transformation in 50 |l competent E. coli XL Gold (Stratagene) cells. Colonies, containing recombinant plasmids, were detected via blue-white screen on LB selection plates, charged with IPTG, ampicillin and X-gal. To check whether the desired insert had been really transferred to the picked out colonies, a minimum amount of each colony (approximately half of it) was utilized as DNA template in a colony PCR. PCRs were run in 50 | l volumes (ingredients, amounts and conditions as above named), using M13F and M13R vector primers. Products were checked on 1% TBE agarose gels and - if they contained an insert of the anticipated size - transferred to LB/ampicillin medium. After proliferation over night, samples were purified according to the guidelines of the Quantum Prep-Kit (Bio Rad) and finally stored at -20°C.

Sequencing and gene annotation

The amplified fragments were set up in 10 | l reaction volumes (2.5 |l DNA, 2.5 |l water, 1 |l primer (10 |M), 4 |l

DCTS master mix) and sequencing PCR reactions were carried out according to the following procedure: 96° C for 20 sec (denaturation); primer-specific temperature for 20 sec (annealing); 60°C for 2 min (elongation). After 30 cycles the samples were sequenced with a CEQ™8000 capillary sequencer (Beckmann-Coulter) and the appropriate CEQ DCTS Quick Start kit (Beckmann-Coulter).

While the first checking of the sequences was carried out with the CEQ 8000 software (Beckman-Coulter), the actual sequence assemblage was done with BioEdit, version 7.0.5 [55]. Protein coding and ribosomal RNA genes, encoded in the mtDNA, were identified by BLAST (blastn, tblastx) searches on NCBI databases and by aligning the different sipunculid fragments with the mt genome of the echiurid Urechis caupo. To revise the final consensus sequence of S. nudus, further mt-genome data of relatively closely related taxa were retrieved from the OGRe database [56]. The species used for sequence comparison were: Platynereis dumerilii (Annelida), Clymenella torquata (Annelida), Orbinia latreillii (Annelida), Lumbricus terres-tris (Annelida), Terebratalia transversa (Brachiopoda), Ter-ebratulina retusa (Brachiopoda), Laqueus rubellus (Brachiopoda), Urechis caupo (echiura), Epiperipatus biol-leyi (Onychophora), and Flustrellidra hispida (Bryozoa), see Table 3 for accession numbers. Transfer RNA genes and their putative secondary structures, were determined with the tRNAscan-SE [57] and ARWEN [58] and for the missing ones by eye inspection of candidate regions. The genome sequence was deposited in NCBI database [GenBank: FT422961].

Phylogenetic analysis

The amino acid alignments of the protein-coding genes (except the two short and highly variable genes atp8 and nad4L) were concatenated. Sequence data from 74 species were included in the large analyses (see Table 3 for all species names and accession numbers). The tree was rooted with two representatives of Cnidaria. Maximum likelihood analysis was performed with RAxML, ver. 7.00 [59,60]. mtREV+G+I was chosen as model for aminoacid substitutions. The complete dataset was partitioned, so that model parameters and amino acid frequencies were optimized for each single gene alignment. 100 bootstrap replicates were performed to infer the support of clades from the best tree. A second set of analyses were done with a reduced dataset of 30 species. This dataset was analyzed with RAxML as described above (model mtREV+G+I, partitioned according to the 12 single gene sequences), with 1000 bootstrap replicates. Secondly we did a Bayesian analysis with MrBayes ver. 3.1.2 [61]. In BI the mtREV+G+I model was used and 1.000.000 generations were run with 8 chains in parallel. Trees were sampled every 1000 generations, while the first 200 trees were discarded as burn-in (according to the likelihood plot). In addition we performed a ML analysis using the

Table 3: Species, systematic position and accession number of mitochondrial genome sequences used in the phylogenetic analysis and/ or for of gene order comparisons

Species Taxonomic position Accession no.

Sipunculus nudus Sipuncula FJ422961

Phascolopsis gouldii* Sipuncula AF374337

Urechis caupo Echiura NC 006379

Myzostoma seymourcollegiorum* Myzostomida EF506562

Lumbricus terrestris Annelida - Clitellata NC 00I677

Perionyx excavatus Annelida - Clitellata NC 00963I

Helobdella robusta* Annelida - Clitellata AFI78678

Platynereis dumerilii Annelida - "Polychaeta" NC 000931

Orbinia latreillii Annelida - "Polychaeta" NC 007933

Eclysippe vanelli* Annelida - "Polychaeta" EU239687

Clymenella torquata Annelida - "Polychaeta" NC 00632I

Pista cristata Annelida - "Polychaeta" NC 0I I0I I

Terebellides stroemi Annelida - "Polychaeta" NC 0I I0I4

Scoloplos armiger* Annelida - "Polychaeta" DQ5I7436

Galathealinum brachiosum* Annelida - Pogonophora AF178679

Riftia pachyptila* Annelida - Pogonophora AY741662

Epiperipatus biolleyi Onychophora NC 009082

Limulus polyphemus Chelicerata - Xiphosura NC 003057

Heptathela hangzhouensis Chelicerata - Araneae NC 005924

Pseudocellus pearsei Chelicerata - Ricinulei NC 009985

Lithobius forficatus Myriapoda - Chilopoda NC 002629

Petrobius brevistylis Hexapoda - Archaeognatha NC 007689

Locusta migratoria Hexapoda - Orthoptera NC 00I7I2

Artemia franciscana Crustacea - Anostraca NC 00I620

Triops cancriformis Crustacea - Phyllopoda NC 004465

Penaeus monodon Crustacea - Decapoda NC 002I84

Priapulus caudatus Priapulida NC 008557

Cephalothrix rufifrons* Nemertea EFI40788

Phoronis psammophila* Phoronida AY368231

Terebratulina retusa Brachiopoda NC 000941

Laqueus rubellus Brachiopoda NC 002322

Terebratalia transversa Brachiopoda NC 003086

Katharina tunicata Mollusca - Polyplacophora NC 00I636

Lottia digitalis Mollusca - Gastropoda NC 007782

Haliotis rubra Mollusca - Gastropoda NC 005940

Conus textile Mollusca - Gastropoda NC 008797

Ilyanassa obsoloeta Mollusca - Gastropoda NC 007781

Thais clavigera Mollusca - Gastropoda NC 0I0090

Lophiotoma cerithiformis Mollusca - Gastropoda NC 008098

Nautilus macromphalus Mollusca - Cephalopoda NC 007980

Octopus ocellatus Mollusca - Cephalopoda NC 007896

Venerupis phllippinarum Mollusca - Bivalvia NC 003354

Argopecten irradians Mollusca - Bivalvia NC 009687

Lampsilis ornata Mollusca - Bivalvia NC 005335

Siphonodentalium lobatum Mollusca - Scaphopoda NC 005840

Loxocorone allax Entoprocta NC 0I043I

Loxosomella aloxiata Entoprocta NC 0I0432

Flustrellidra hispida Bryozoa/Ectoprocta NC 008I92

Paraspadella gotoi Chaetognatha NC 006083

Spadella cephaloptera Chaetognatha NC 006386

Brachionus plicatilis Rotifera NC 0I0484

Leptorhynchoides thecatus Acanthocephala NC 006892

Anisakis simplex Nematoda NC 007934

Agamermis sp. Nematoda NC 008231

Onchocercas volvulus Nematoda NC 00I86I

Caenorhabditis elegans Nematoda NC 00I328

Trichinella spiralis Nematoda NC 002681

Xiphinema americanum Nematoda NC 005928

Table 3: Species, systematic position and accession number of mitochondrial genome sequences used in the phylogenetic analysis and/ or for of gene order comparisons (Continued)

Paratomella rubra* Acoela AY228758

Microstomum lineare* Platyhelminthes - "Turbellaria" AY228756

Fasciola hepatica Platyhelminthes - "Trematoda" NC 002546

Paragoniums westermanni Platyhelminthes - "Trematoda" NC 002354

Gyrodactylus salaris Platyhelminthes - "Trematoda" NC 008815

Microcotyle sebastis Platyhelminthes - "Trematoda" NC 009055

Schistosoma haematobium Platyhelminthes - "Trematoda" NC 008074

Trichobilharzia regenti Platyhelminthes - NC 009680

Hymenolepis diminuta Platyhelminthes - Cestoda NC 002767

Taenia asiatica Platyhelminthes - Cestoda NC 004826

Echinococcus granulosus Platyhelminthes - Cestoda NC 008075

Balanoglossus carnosus Enteropneusta NC 001887

Saccoglossus kowalevskii Enteropneusta NC 007438

Florometra serratissima Echinodermata - Crinoidea NC 001878

Paracentrotus lividus Echinodermata - Echinoidea NC 001572

Xenoturbella bocki Xenoturbellida NC 008556

Acropora tenuis Cnidaria - Anthozoa NC 003522

Metridium senile Cnidaria - Anthozoa NC 000933

*: incomplete genome sequence

mtART+G+I model with Treefinder [62] and "edge support" analysis, again with a partitioned dataset (= independently optimizing model parameters for the 12 genes).

the manuscript, all other authors helped in interpretation of data and discussion of results.

Additional material

For comparison of the hypothesis that sipunculids might be closely related with molluscs and our best tree, we used a constraint for a ML-analysis (Sipuncula + Mollusca) of the sequence dataset using RaxML [59] with parameters described above. We computed per-site log-likelihoods with RAxML for both topologies (best tree and constrained topology) and conducted an au-test as implemented in CONSEL [63].

Additional File 1

Full version of figure2. Secondary structure of tRNAs identified in the mitochondrial genome of S. nudus. The best found putative secondary structure of tRNA-Cys (box) seems to be strongly derived, probably nonfunctional or subject to gene editing. Click here for file



atp6 and 8: genes encoding ATPase subunit 6 and 8; bp: base pairs; bt: bootstrap; cox 1-3: genes encoding cytochrome oxidase subunits I-III; cob: gene encoding cytochrome b; BI: Bayesian Inference; ML: Maximum Likelihood; mtDNA: mitochondrial DNA; mt-genome: mitochondrial genome; nad1-6 and nad4L: genes encoding NADH dehydroenase subunits 1-6 and 4L; PCR: polymerase chain reaction; rRNA: ribosomal RNA; rrnL: large rRNA subunit (16S); rrnS: small rRNA subunit (12S); tRNA: transfer RNA; trnX: tRNA gene (X is replaced by one letter amino acid code).

Authors' contributions

LP and TB conceived and supervised this study. AcM and BL constructed the EST libarary. AdM did all PCR experiments and sequencing of the mitochondrial genome. AdM and LP annotated the genome. AdM, LP and CB performed phylogenetic analysis of sequence data, LP analyzed gene order data. AdM and LP wrote the main part of


We are grateful to Dominik Kieselbach for providing us with a S. nudus specimen. Financial support for this study came from Deutsche Forschungsgemeinschaft (DFG): DFG Ba 1520/10-1,2 to TB and LP; DFG BL 787/2-1 to CB, both from priority programme "Deep Metazoan Phylogeny" SPP 1 174.


1. Rokas A, Williams BL, King N, Carroll SB: Genome-scale approaches to resolving incongruence in molecular phyloge-nies. Nature 2003, 425:798-804.

2. Boore JL: Animal mitochondrial genomes. Nucleic Acids Res 1999, 27:1767-1780.

3. Wolstenholme DR: Animal mitochondrial DNA: structure and evolution. Int Rev Cytol 1992, 141:173-216.

4. Shao R, Aoki Y, Mitani H, Tabuchi N, Barker SC, Fukunaga M: The mitochondrial genomes of soft ticks have an arrangement of genes that has remained unchanged for over 400 million years. Insect Mol Biol 2004, 13:219-224.

5. Boore JL, Collins TM, Stanton D, Daehler LL, Brown WM: Deducing the pattern of arthropod phylogeny from mitochondrial DNA rearrangements. Nature London 1995, 376:163-165.

6. Valles Y, Boore JL: Lophotrochozoan mitochondrial genomes. Integrative and Comparative Biology 2006, 46:544-557.

7. Dowton M, Castro LR, Austin AD: Mitochondrial gene rearrangements as phylogenetic characters in the invertebrates: The examination of genome 'morphology'. Invertebrate Systematics 2002, 16:345-356.

8. Bleidorn C, Eeckhaut I, Podsiadlowski L, Schult N, McHugh D, Halan-ych KM, Milinkovitch MC, Tiedemann R: Mitochondrial genome and nuclear sequence data support myzostomida as part of the annelid radiation. Mol Biol Evol 2007, 24:1690-1701.

9. Boore JL, Lavrov DV, Brown WM: Gene translocation links insects and crustaceans. Nature 1998, 392:667-668.

10. Lavrov DV, Brown WM, Boore JL: Phylogenetic position of the Pentastomida and (pan)crustacean relationships. Proc R Soc Lond B Biol Sci 2004, 271:537-544.

11. Haen KM, Lang BF, Pomponi SA, Lavrov DV: Glass sponges and bilaterian animals share derived mitochondrial genomic features: a common ancestry or parallel evolution? Mol Biol Evol 2007, 24:1518-1527.

12. Cutler EB: The Sipuncula their systematics, biology and evolution Ithaca, N.Y: Comstock Pub. Associates; 1994.

13. Huang DY, Chen JY, Vannier J, Salinas JIS: Early Cambrian sipun-culan worms from southwest China. Proceedings of the Royal Society of London Series B-Biological Sciences 2004, 271:1671-1676.

14. Ax P: Multicellular Animals. The Phylogenetic System of the Metazoa Berlin, Heidelberg, New York: Springer; 2000.

15. Schulze A, Cutler EB, Giribet G: Phylogeny of sipunculan worms: A combined analysis of four gene regions and morphology. Molecular Phylogenetics and Evolution 2007, 42:171-192.

16. Nichols D: The origin of echinoderms. Symposia of the Zoological Society London 1967, 20:209-229.

17. Scheltema AH: Aplacophora As Progenetic Aculiferans and the Coelomate Origin of Mollusks As the Sister Taxon of Sipuncula. Biological Bulletin 1993, 184:57-78.

18. Maslakova SA, Martindale MQ, Norenburg JL: Fundamental properties of the spiralian developmental program are displayed by the basal nemertean Carinoma tremaphoros (Palaeone-mertea, Nemertea). Developmental Biology 2004, 267:342-360.

19. Nielsen C: Animal Evolution. Interrelationships of the living phyla Oxford: Oxford University Press; 2001.

20. Meglitsch PA, Schram FR: Invertebrate zoology New York; Oxford: Oxford University Press; 1991.

21. Brusca RC, Brusca GJ: Invertebrates Sunderland, Mass: Sinauer Associates; 2003.

22. Erber A, Riemer D, Bovenschulte M, Weber K: Molecular phylogeny of metazoan intermediate filament proteins. Journal of Molecular Evolution 1998, 47:751-762.

23. Eernisse DJ, Albert JS, Anderson FE: Annelida and Arthropoda Are Not Sister Taxa - A Phylogenetic Analysis of Spiralian Metazoan Morphology. Systematic Biology 1992, 41:305-330.

24. Schulze A, Cutler EB, Giribet G: Reconstructing the phylogeny of the Sipuncula. Hydrobiologia 2005, 535:277-296.

25. Struck TH, Schult N, Kusen T, Hickman E, Bleidorn C, McHugh D, Halanych KM: Annelid phylogeny and the status of Sipuncula and Echiura. BMC Evol Biol 2007, 7:57.

26. Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, et al.: Broad phylog-enomic sampling improves resolution of the animal tree of life. Nature 2008, 452:745-7U5.

27. 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. Molecular Biology and Evolution 2007, 24:2723-2729.

28. Bleidorn C, Podsiadlowski L, Bartolomaeus T: The complete mito-chondrial genome of the orbiniid polychaete Orbinia latreillii (Annelida, Orbiniidae) - A novel gene order for Annelida and implications for annelid phylogeny. Gene 2006, 370:96-103.

29. Boore JL, Staton JL: The mitochondrial genome of the Sipuncu-lid Phascolopsis gouldii supports its association with Annelida rather than Mollusca. Mol Biol Evol 2002, 19:127-137.

30. Bleidorn C: The role of character loss in phylogenetic reconstruction as exemplified for the Annelida. Journal of Zoological Systematics and Evolutionary Research 2007, 45:299-307.

31. Kristof A, Wollesen T, Wanninger A: Segmental mode of neural patterning in Sipuncula. Curr Biol 2008, 18:1129-1132.

32. Cutler EB, Gibbs PE: A Phylogenetic Analysis of Higher Taxa in the Phylum Sipuncula. Systematic Zoology 1985, 34:162-173.

33. Schulze A, Cutler EB, Giribet G: Molecular and morphological evolution in sipunculan worms. Integrative and Comparative Biology 2005, 45:1070.

34. Staton JL: Phylogenetic analysis of the mitochondrial cytochrome c oxidase subunit 1 gene from 13 sipunculan genera: intra- and interphylum relationships. Invertebrate Biology 2003, 122:252-264.

35. Perna NT, Kocher TD: Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol 1995, 41:353-358.

36. Zhong M, Struck TH, Halanych KM: Phylogenetic information from three mitochondrial genomes of Terebelliformia (Annelida) worms and duplication of the methionine tRNA. Gene 2008, 416:11-21.

37. Ojala D, Montoya J, Attardi G: tRNA punctuation model of RNA processing in human mitochondria. Nature 1981, 290:470-474.

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

39. Zhang DX, Hewitt GM: Insect mitochondrial control region: A review of its structure, evolution and usefulness in evolutionary studies. Biochemical Systematics and Ecology 1997, 25:99-120.

40. Lavrov DV, Lang BF: Poriferan mtDNA and animal phylogeny based on mitochondrial gene arrangements. Systematic Biology 2005, 54:651-659.

41. Kilpert F, Podsiadlowski L: The complete mitochondrial genome of the common sea slater, Ligia oceanica (Crustacea, Isopoda) bears a novel gene order and unusual control region features. BMC Genomics 2006, 7:241.

42. Bandyopadhyay PK, Stevenson BJ, Ownby JP, Cady MT, Watkins M, Olivera BM: The mitochondrial genome of Conus textile, coxI-coxII intergenic sequences and Conoidean evolution. Molecular Phylogenetics and Evolution 2008, 46:215-223.

43. Simison WB, Lindberg DR, Boore JL: Rolling circle amplification of metazoan mitochondrial genomes. Molecular Phylogenetics and Evolution 2006, 39:562-567.

44. Stechmann A, Schlegel M: Analysis of the complete mitochon-drial DNA sequence of the brachiopod Terebratulina retusa places Brachiopoda within the protostomes. Proceedings of the Royal Society Biological Sciences Series B 1999, 266:2043-2052.

45. Turbeville JM, Smith DM: The partial mitochondrial genome of the Cephalothrix rufifrons (Nemertea, Palaeonemertea): Characterization and implications for the phylogenetic position of Nemertea. Mol Phylogenet Evol 2007, 43:1056-1065.

46. Helfenbein KG, Boore JL: The mitochondrial genome of Pho-ronis architecta-comparisons demonstrate that phoronids are lophotrochozoan protostomes. Mol Biol Evol 2004, 21:153-157.

47. Yokobori SI, Lindsay DJ, Tsuchiya K, Yaniagishi A, Maruyama T, Shirna T: Mitochondrial genome structure and evolution in the living fossil vampire squid, Vampyroteuthis infernalis, and extant cephalopods. Molecular Phylogenetics and Evolution 2007, 44:898-910.

48. Endo K, Noguchi Y, Ueshima R, Jacobs HT: Novel repetitive structures, deviant protein-encoding sequences and unidentified ORFs in the mitochondrial genome of the brachiopod Lin-gula anatina. J Mol Evol 2005, 61:36-53.

49. Helfenbein KG, Brown WM, Boore JL: The complete mitochondrial genome of the articulate brachiopod Terebratalia transversa. Mol Biol Evol 2001, 18:1734-1744.

50. Noguchi Y, Endo K, Tajima F, Ueshima R: The mitochondrial genome of the brachiopod Laqueus rubellus. Genetics 2000, 155:245-259.

51. Boore JL: Complete mitochondrial genome sequence of Urechis caupo, a representative of the phylum Echiura. BMC

Genomics 2004, 5:67.

52. Halanych KM, Bacheller JD, Aguinaldo AM, Liva SM, Hillis DM, Lake JA: Evidence from I8S ribosomal DNA that the lopho-phorates are protostome animals. Science 1995, 267:1641-1643.

53. Halanych KM: The new view of animal phylogeny. Annual Review of Ecology Evolution and Systematics 2004, 35:229-256.

54. Mallatt J, Winchell CJ: Testing the new animal phylogeny: First use of combined large-subunit and small-subunit rRNA gene

sequences to classify the protostomes. Molecular Biology and Evolution 2002, 19:289-301.

55. Hall TA: BioEdit a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 1999, 41:95-98.

56. Jameson D, Gibson AP, Hudelot C, Higgs PG: OGRe: a relational database for comparative analysis of mitochondrial genomes. Nucleic Acids Res 2003, 31:202-206.

57. Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res 1997, 25:955-964.

58. Laslett D, Canback B: ARWEN: a program to detect tRNA genes in metazoan mitochondrial nucleotide sequences. Bio-informatics 2008, 24:172-175.

59. Stamatakis A: RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 2006, 22:2688-2690.

60. Stamatakis A, Hoover P, Rougemont J: A rapid bootstrap algorithm for the RAxML web servers. Systematic Biology 2008, 57:758-771.

61. Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 2001, 17:754-755.

62. Jobb G: Treefinder. Version of Feb. 2007. Munich, Germany. 2007 []. no journal

63. Shimodaira H: An approximately unbiased test of phylogenetic tree selection. Syst Biol 2002, 51:492-508.

Publish with BioMecl Central and every scientist can read your work free of charge

"BioMed Central will be the most significant development for disseminating the results of biomedical research in our lifetime." Sir Paul Nurse, Cancer Research UK

Your research papers will be:

• available free of charge to the entire biomedical community

• peer reviewed and published immediately upon acceptance

• cited in PubMed and archived on PubMed Central

• yours — you keep the copyright

Submit your manuscript here: i } BioMedcentral ^