Scholarly article on topic 'Phylogeny of Annelida (Lophotrochozoa): total-evidence analysis of morphology and six genes'

Phylogeny of Annelida (Lophotrochozoa): total-evidence analysis of morphology and six genes Academic research paper on "Biological sciences"

Share paper
Academic journal
BMC Evol Biol
OECD Field of science

Academic research paper on topic "Phylogeny of Annelida (Lophotrochozoa): total-evidence analysis of morphology and six genes"

BMC Evolutionary Biology Bi.Medce ta

Open Access

Research article

Phylogeny of Annelida (Lophotrochozoa): total-evidence analysis of morphology and six genes

Jan Zrzavy*1,2, Pavel Riha1, Lubomir Pialek1 and Jan Janouskovec1,2

Address: faculty of Science, University of South Bohemia, Branisovska 31, 370 05 Ceske Budejovice, Czech Republic and 2Biology Center, Academy of Sciences, Branisovska 31, 370 05 Ceske Budejovice, Czech Republic

Email: Jan Zrzavy* -; Pavel Riha -; Lubomir Pialek -; Jan Janouskovec - * Corresponding author

Published: 6 August 2009 Received: 13 July 2008

BMC Evolutionary Biology 2009, 9:189 doi:l0.l 186/1471-2148-9-189 Accepted: 6 August 2009

This article is available from: http://www.biomedcentral.eom/l47l-2l48/9/l89 © 2009 Zrzavy 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: Annelida is one of the major protostome phyla, whose deep phylogeny is very poorly understood. Recent molecular phylogenies show that Annelida may include groups once considered separate phyla (Pogonophora, Echiurida, and Sipunculida) and that Clitellata are derived polychaetes. SThe "total-evidence" analyses combining morphological and molecular characters have been published for a few annelid taxa. No attempt has yet been made to analyse simultaneously morphological and molecular information concerning the Annelida as a whole.

Results: Phylogenetic relationships within Annelida were analysed on the basis of 93 morphological characters and sequences of six genes (I8S, 28S, and I6S rRNA, EFI a, H3, COI), altogether, 87 terminals of all annelid "families" and 3,903 informative characters, by Bayesian and maximum-parsimony methods. The analysis of the combined dataset yields the following scheme of relationships: Phyllodocida and Eunicida are monophyletic groups, together probably forming monophyletic Aciculata (incl. Orbiniidae and Parergodrilidae that form a sister group of the Eunicida). The traditional "Scolecida" and "Canalipalpata" are both polyphyletic, forming instead two clades: one including Cirratuliformia and the "sabelloid-spionoid clade" (incl. Sternaspis, Sabellidae-Serpulidae, Sabellariidae, Spionida s.str.), the other ("terebelloid-capitelloid clade") including Terebelliformia, Arenicolidae-Maldanidae, and Capitellidae-Echiurida. The Clitellata and "clitellate-like polychaetes" (Aeolosomatidae, Potamodrilidae, Hrabeiella) form a monophyletic group. The position of the remaining annelid groups is uncertain - the most problematic taxa are the Opheliidae-Scalibregmatidae clade, the Amphinomida-Aberranta clade, Apistobranchus, Chaetopteridae, Myzostomida, the Sipunculida-Dinophilidae clade, and the "core Archiannelida" (= Protodrilidae, Nerillidae, Polygordiidae, Saccocirridae).

Conclusion: The combined ("total-evidence") phylogenetic analysis provides a modified view of annelid evolution, with several higher-level taxa, i.e. Phyllodocida, Eunicida, orbinioid-parergodrilid clade (OPC), Cirratuliformia, sabelloid-spionoid clade (SSC), terebelloid-capitelloid clade (TCC), and "Clitellatomorpha". Two unorthodox clades, the "core Archiannelida" and Sipunculida-Dinophilidae, are proposed. Although the deep-level evolutionary relationships of Annelida remain poorly understood, we propose the monophyly of the Aciculata, sister-group relationships between the Eunicida and OPC, between the Cirratuliformia and SSC, and possibly also between the "Clitellatomorpha" and Oweniidae-Pogonophora clades.


Annelida, the segmented worms (over 16,500 species described), are distributed worldwide from the deepest marine sediments to freshwater and soil habitats. Throughout most of the 20th century they were split into three or four major groups, Polychaeta, Myzostomida, Oligochaeta and Hirudinea. It is now widely recognized that Oligochaeta and Hirudinea form a clade that is referred to as Clitellata (where leeches are only a derived subgroup of oligochaetes [1-3]). Several interstitial groups were classified as the "Archiannelida", another annelid group; however, they are now generally regarded as secondarily simplified, possibly progenetic polychaetes [4,5]. Several more groups have been hypothesized to belong into the Annelida [6], and there is a growing consensus that the Echiurida, Pogonophora (incl. Vestimen-tifera), and Sipunculida are actually modified annelids [79].

A cladistic analysis of Annelida, based on morphological characters, has resulted in a new classification [10,11], with three major clades of the Polychaeta: Scolecida, Acic-ulata (= Amphinomida + Eunicida + Phyllodocida), and Canalipalpata (= Terebellida + Spionida + Sabellida [incl. Pogonophora]). However, several annelid groups were left outside this classification. They include Clitellata, the freshwater and/or terrestrial "clitellate-like" worms (Par-ergodrilidae, Hrabeiella, and Aphanoneura [= Aeolosoma-tidae + Potamodrilus]), some "archiannelids" (Protodrilida and Polygordiidae, both only tentatively regarded as aberrant canalipalpatans), and Psammodrilidae.

From a molecular perspective, the sequence datasets assembled to date have usually been marked by limited numbers of both taxa and characters. Almost all annelid families are now represented by the nuclear small-subunit ribosomal RNA genes ("18S" hereinafter); unfortunately, even 18S studies using the densest taxon sampling [12-14] were unable to recover a monophyletic Annelida or its major subclades. Even if several genes are concatenated to reconstruct annelid phylogeny in recent papers, none of the morphology-based higher taxa (Polychaeta, Scolecida, Aciculata, Canalipalpata) were recovered [15-17]. Recent papers by Struck et al. [7,8] provided the first molecular trees with several resolved higher taxa of the Annelida. They included Aciculata (excl. Amphinomida), Phyllodocida (incl. Orbiniidae), Terebelliformia, Sabellida-Spion-ida, Cirratuliformia, and Amphinomida.

Synthesis of molecular and morphological data from extant and potentially also extinct taxa remains the strongest test of phylogenetic hypotheses and the best summary of the common signal in the diverse data available for phylogenetics [18]. The "total-evidence" analyses have been published for a few annelid taxa, viz., Clitellata [2], Terebelliformia [19], most Canalipalpata [20], Aphroditi-

formia [21], and most Aciculata [22]. So far, no attempt has been made to analyse simultaneously morphological and molecular information on the Annelida as a whole.

In this paper we present the first comprehensive analysis of higher-level phylogenetic relationships in Annelida based on combined morphological and molecular (four nuclear, two mitochondrial genes) data. The purpose is to identify stable and ustable nodes of the combined annelid tree, to make up reliable phylogenetic hypothesis on Annelida, and thus test the morphology-based classification.


The congruence of data partitions

The combined data matrix included 87 terminals and 3,903 cladistically informative characters (93 morphology [= MOR]; 630 cytochrome c oxidase subunit I [= COI]; 604 elongation factor-1a [= EF1a]; 132 histone H3; 763 18S ribosomal RNA; 274 16S ribosomal RNA; 1,407 28S ribosomal RNA). Out of339,561 data matrix cells (= character states), 6% were polymorphic and 30% were missing (unknown or inapplicable). Distribution of the ambiguous character states is highly uneven, as four terminals (Histriobdellidae, Sphaerodoridae, Aberranta, Parapo-drilus) are represented by two data partitions only (MOR + 18S).

Evaluation of the relative quality of data partitions and their performance in the combined maximum-parsimony analysis indicated that the MOR dataset was highly influential in the simultaneous analyses of all data partitions (Figure 1). The partitioned Bremer support (PBS) values in the combined dataset, limited to the 28 composite terminals that were represented by all seven data partitions (one tree, length 18,195, CI 0.35, RI 0.24; Figure 1), revealed that despite the significant incongruence of the morphological and molecular data partitions, the former contributed positively to Bremer support values of (= supported) 15 clades (60%) and is negative for (= contradicted) two clades (8%) only: one within Terebelliformia, the other concerning position of the Orbiniidae within Aciculata. The molecular partitions analysed together are in conflict with two clades as well, both concerning the placement of the Amphinomidae as a sister group of the Eunicida. The molecular partitions are, however, by no means homogeneous. The PBS analysis showed that four data partitions contributed positively to the combined-tree topology (MOR: 2PBS = 100; 28S: 322; EF1a: 232; 16S: 100), while the other three molecular partitions (COI, H3,18S) are in conflict with the combined tree.

Bayesian (Bl) tree

In the combined BI tree, the basal relationships were not fully resolved. However, there were several well-supported (pp > 0.90) clades (Figure 2): (i) the "sabelloid-spionoid


- Chaetopteridae - Oweniidae


!!■! Vestimentifera

■ Opheliidae


■ Ampharetidae ■■■ Trichobranchidae - Arenicolidae -Maldanidae ■ Echiurida


Cirratulidae -Sabellidae

Sabellariidae Spionidae Amphinomidae Lumbrineridae Eunicidae Orbiniidae

Polynoidae „.„, Hesionidae 0 ■■ Phyllodocidae

ESI i-QQi-Glyceridae



18S 16S 28S


Evaluation of the performance of data partitions in the combined maximum-parsimony tree. The analysis was limited to the 28 terminals, represented by all seven data partitions [see Additional File 4]. Bootstrap values are above the branches (only higher than 50%); partitioned Bremer support (PBS) values under the branches: COMB (total Bremer support of all partitions), MOR (PBS of morphology), MOL (total Bremer support of the six molecular partitions). The support grids: black cells: PBS > 0; white cells: PBS < 0; diagonal: PBS = 0. For color code see Figure 2.

clade" ("SSC" hereinafter; incl. Sabellida s.str. and Spionida s.str.); (ii) Cirratuliformia (incl. Histriobdellidae); (iii) Fauveliopsidae-Paraonidae-Cossuridae; (iv) Clitel-lata-Aphanoneura-Hrabeiella ("Clitellatomorpha" hereinafter); (v) Pogonophora; (vi) Opheliidae-Scalibregmatidae; (vii) the "terebelloid-capitelloid clade" or "TCC" (incl. Terebelliformia, Maldanidae-Arenicolidae, and Capitellidae-Echiurida); (viii) Phyllodocida; (ix) Eunicida; (x) the "orbinioid-parergodrilid clade" or "OPC"; (xi) Aberranta-Amphinomida; and (xii) the "core Archiannelida" (Polygordiidae-Saccocirridae and Proto-drilidae-Nerillidae). The traditional Aciculata and Canali-palpata roughly corresponded to the two major superclades (both, however, with rather low posterior probabilities and very short basal branches). The "Acicu-lata" included Phyllodocida, Eunicida, and OPC. The "Canalipalpata" included SSC, Sternaspis, TCC, Cirratuli-

formia, Fauveliopsidae-Paraonidae-Cossuridae, Pogono-phora, and also "Clitellatomorpha" and Opheliidae-Scalibregmatidae. There were, in addition, several groups outside the two major sister superclades: Oweniidae, Chaetopteridae, Magelona, Myzostomida-Protodriloides, Aberranta-Amphinomida, Apistobranchus-Dinophilidae-Sipunculida, and the "core Archiannelida".

Maximum-parsimony (MP) trees

In the unweighted (uMP) analysis (one tree, length 42,492, CI 0.20, RI 0.25), monophyletic Aciculata were nested within paraphyletic "Canalipalpata" (Figure 3). A few taxa seemed to be misplaced due to long-branch artifacts (e.g. Serpulidae far from the Sabellidae, Pectinariidae outside the Terebelliformia, and Tomopteridae outside the Phyllodocida); moreover, there was an obviously arti-factual "basal" clade including "archiannelids", myzosto-mids, and several long-branch aciculatans and canalipalpatans.

The weighted (wMP, "slow-fast") tree (one tree, weighted length 129,822, CI 0.25, RI 0.27; Figure 4) was nearly identical with the BI tree (both included two major sister superclades, one mostly aciculatan, the other mostly canalipalpatan). They differed predominantly only in position of the Opheliidae-Scalibregmatidae clade, "core archiannelids", Dinophilidae-Sipunculida clade, Apisto-branchus, and Tomopteridae.

Taxon-exclusion tests

The preliminary results suggested that ingroup and out-group taxon selection affected internal relationships. To highlight topological problems in placement of the most "problematic" taxa (Myzostomida, Opheliidae, Scalibreg-matidae, Fauveliopsidae, Serpulidae, Magelona, Owenii-dae, Chaetopteridae, Pectinariidae, Apistobranchus, Polygordiidae, Saccocirridae, Protodrilidae, Protodriloides, Nerillidae, Histriobdellidae, Dinophilidae, Sipunculida, Hesionidae, Paralacydoniidae, Nepthyidae, Pilargidae, Tomopteridae, Aberranta, Amphinomida; see Methods for selection criteria), several experimental analyses were performed.

In the taxon-exclusion tests, after removing the "problematic" species/clades, the unrooted "backbone trees" were almost identical irrespective of the used methods (uMP and BI). Both included either monophyletic or para-phyletic (which naturally cannot be distinguished in an unrooted tree) "Aciculata" and "Canalipalpata", and differed only in the specific position of the Pogonophora-"Clitellatomorpha" clade within the "Canalipalpata" (Figure 5, 6). The important "problematica" grouped consistently as follows: Oweniidae as a sister group of the Pogonophora; Fauveliopsidae, Magelona, and Opheliidae-Scalibregmatidae within the Cirratuliformia-SSC clade;

Apistobranchus as a sister group of OPC; and Aberranta-Amphinomida next to (or within) the Phyllodocida. The other "problematic species" were placed at more or less conflicting positions in BI and uMP tests (see below). The ougroups tended to be placed within the "Canalipalpata", close to the Pogonophora and/or "Clitellatomorpha".


The present analysis points out instability of some basal nodes within the annelid tree. However, most analyses, irrespective of the used methods and parameter sets, converge to a highly compatible set of actually (or, in the unrooted trees, potentially) monophyletic taxa. They are (i) Eunicida-OPC; (ii) Phyllodocida; (iii) SSC; (iv) Cirrat-uliformia; (v) TCC; (vi) "Clitellatomorpha"; and (vii) Pogonophora (Figure 2, 4, 5 &6). Several morphology-based taxa like "Scolecida" (= Orbiniidae-Questa, Opheli-idae-Scalibregmatidae, Arenicolidae-Maldanidae, Capitel-lidae), "Terebellida" (= Terebelliformia and Cirratuliformia), and "Spionida" (= Spionida s.str., Magelona, Chaetopteridae, Apistobranchus) [10,11] are definitely polyphyletic.

The most problematic issue is the tree root position. The all-taxa analyses (that included mollusc and brachiopod outgroups) supported that the root is situated between chaetopterids, magelonids, and/or oweniids and the rest of the Annelida (Figure 2, 3 &4). The experimental analyses (Figures 5, 6) suggested that the non-annelid out-groups were placed next to the Oweniidae, Pogonophora, and/or "Clitellatomorpha".

Phyllodocida, Eunicida, and the orbinioid-parergodrilid clade (OPC)

The Phyllodocida (excl. Myzostomida) is consistently a monophyletic group in all analyses, except for "basal" placement of the long-branch Tomopteridae in some MP trees. In all analyses, there is a strongly supported clade including monophyletic Eunicida (less dinophilids and sometimes also histriobdellids) and OPC, corroborated by presence of the lateral/dorsal cirrus organs and larval akrotrochs. However, the position of the lumbrinerids could be more basal, as they appeared as a sister group of the whole Eunicida-Phyllodocida in some purely molecular trees [7,8,23]. The placement of the orbiniids among the aciculatans provides a strong phylogenetic support for the hypothesis that their aciculae are homologous with those of the Euncida and Phyllodocida [24]. The sister-group relationship between Parergodrilidae and Orbinii-dae-Questa is strongly supported by all analyses [8,24]; both groups share gonoducts with a distal glandular part [25]. No closer orbiniid-spionidan relationships [26] have been recovered by the present analyses. Also Apisto-branchus consistently grouped with OPC in the taxon-exclusion analyses. It is, therefore, possible that also the

Apistobranchus acicula-like chaetae [24] might be actually true aciculae.

"Terebelloid-capitelloid clade" (TCC)

One of the most stable clade covers Terebelliformia, Capitellidae-Echiurida and Arenicolidae-Maldanidae. This clade is weakly supported morphologically, by the specific chaetal arrangement and by the presence of a gular membrane. The placement of the terebelliformians close to a generally palpless group provides a phylogenetic support for the hypothesis that terebelliformian "buccal tentacles" are not homologous with the true polychaete palps [27]. Although classified as a separate phylum until recently, there is a growing consensus that Echiurida is a subgroup of the Annelida [28], as they exhibit segmentation traits during ontogeny. The explicit hypotheses about their specific position were published only by molecular phylogeneticists who discovered close relationships of the Echiurida to the Capitellidae [7-9,13-16]. Also the Ophe-liidae and Scalibregmatidae group either within TCC (present MP, [7,8,15]), or as an isolated clade (present BI). However, in absence of all other problematic annelids, the Opheliidae-Scalibregmatidae clade nests close to the SSC.

"Sabelloid-spionoid clade" (SSC) and Cirratuliformia

In all trees, Sabellidae, Spionida s.str. (= Spionidae + Tro-chochaeta + Poecilochaetus), and Sabellariidae form a clade [8]. The chaetal arrangement in Sabellariidae has recently been dismissed as homologous to the sabellid-serpulid chaetal inversion pattern, and therefore provides no support for a sister-group relationship of sabellids and sabel-lariids with exclusion of the Spionida s.str. [29]. Fauveliopsidae, Sternaspis, Cossuridae and Paraonidae are consistently placed next to the SSC and/or Cirratuliformia [[7,17,20,30]; cf. [8]]. The monophyly of the "Terebellida" (joining Cirratuliformia and Terebelliformia [10,11]) has never been found in the molecular and combined trees.


Phylogenetic position of the Clitellata and most "clitel-late-like" annelids is stable and quite surprising in the present analyses: they always form a clade, further split into Clitellata and the Hrabeiella-Aphanoneura subclade. This clade is supported by hermaphroditism, direct development, loss of parapodia, presence of the unbranched type of ciliary ocelli, metanephridial mantle cell, possibly also by the dorsal pharynx, and by primarily freshwater/ soil habitats [3] (though many of these characters are absent or unknown either in Aphanoneura or in Hrabeiella). The intra-clitellate position of the Aphanone-ura was not supported by reliable morphological synapo-morphies [[31,32]; cf. [33]], and exclusion of aphanoneurans from the Clitellata was also indicated by


Bayesian tree for the combined dataset (morphology + six molecular partitions). Posterior probabilities are shown on the branches. Terminals with just the higher-taxa names (e.g. "FlabeNigeridae") indicate that the sequences from different species of that taxon were concatenated [see Additional File 4].


Unweighted maximum-parsimony (uMP) tree for the combined dataset (morphology + six molecular partitions). Bootstrap/Bremer support values are below the branches. Black/white ellipse hashmarks indicate presence/absence of the clade in purely morphological (left) and molecular (right) trees (uMP), respectively (functionally monotypic taxa are not hashmarked). For color code see Figure 2.


Weighted ("slow-fast") maximum-parsimony (wMP) tree for the combined dataset (morphology + six molecular partitions). Presence of a clade in the 11 trees derived from the reduced datasets (see text) is shown below the branches (in %). Black/white ellipse hashmarks indicate presence/absence of the clade in purely morphological (left) and molecular (right) trees (both wMP), respectively.

Figure 5

Results of the taxon-exclusion maximum-parsimony analyses (morphology + six molecular partitions). The

trees are unrooted (position of possible outgroups, Brachiopoda and Mollusca, indicated). The "backbone tree" is colored, the "problematic taxa" appended in one-by-one manner are in white. Unweighted maximum parsimony; bootstrap values/Bremer support values of the "backbone tree" clades are below the branches. Black/white ellipse hashmarks indicate presence/absence of the clade in purely morphological (left) and molecular (right) "backbone taxa" (uMP), respectively.

Figure 6

Results of the taxon-exclusion Bayesian analyses (morphology + six molecular partitions). The trees are unrooted (position of possible outgroups, Brachiopoda and Mollusca, indicated). The "backbone tree" is colored, the "problematic taxa" appended in one-by-one manner are in white. Posterior probabilities are shown on the branches.

most previous molecular analyses [8,14,31,34-37]. The recent molecular analyses, on the contrary, seem to corroborate close relationships between Clitellata and (some of) the "clitellate-like polychaetes" [7,8,16,31,37,38]. The close relationships of the Clitellata to some eunicidans or dinophilids [15,16], or to TCC [7,8] have never been found in the present analyses.

Pogonophora and "basal Canalipalpata"

In the present analyses, the Pogonophora (= Siboglinidae) is closely related to the "Clitellatomorpha" or, in absence of the latter group, Pogonophora alone groups as a sister group of SSC. Three groups of the traditional Canalipal-pata, namely, Oweniidae, Magelona, and Chaetopteridae, were usually placed among the most basal annelids, or even outside the Annelida, in the present all-species analyses [7-9,15,17]. The taxon-exclusion analyses suggested that magelonids and oweniids might in fact be attracted strongly towards the remote annelid outgroups. In absence of all other "problematic" annelids, Oweniidae alone grouped consistently with the Pogonophora (supported predominantly by the intraepidermal nerve cord [7,20]), and Magelona either with the Paraonidae (uMP) or as a sister group of the Cirratuliformia-SSC superclade (BI). On the contrary, basal placement of the Chaetopteri-dae cannot be excluded [9].

Amphinomida and Aberranta

In the most present analyses, Aberranta is a sister group of the Amphinomida, and the whole clade is one of the basalmost annelid branches, far from the other Aciculata [7,8,15]. However, in the taxon-exclusion tests, the Aber-ranta-Amphinomida clade (weakly supported by median prostomial antennae and parapodial branchiae) is placed as a sister group of the Phyllodocida (BI) or even within it (uMP). It is then possible that the Aberranta-Amphinom-ida clade is, in fact, one of subgroups of the monophyletic Aciculata [10,11], misplaced in most molecular trees (due to their strong attraction towards the outgroups?). Therefore, it seems rather premature to regard Amphinomida as a basal group, "based on the tetraneurous organization of the nervous system" [7]; moreover, the presence of several peripheral longitudinal nerves is not limited to the amphinomids [27].

"Core archiannelids" (Nerillidae, Polygordiidae, Protodrilidae, Protodriloides)

Monophyly of the "Archiannelida" has been rejected by all modern authors [4], whereas the monophyly of Proto-drilida (= Saccocirridae + Protodrilidae + Protodriloidides) as well as its spionidan affinities has been proposed recently [4]. The polygordiids were regarded as relatives of protodrilidans or of opheliids [39], or as the canalipalpa-tans of uncertain position [11]. Surprisingly, the molecular trees tend to include a clade covering most of the traditional "archiannelids" [8]. In the present all-taxa BI

tree, there is a well-supported basal annelid clade including Nerillidae-Protodrilidae and Polygordiidae-Saccocirri-dae. Although Nerillidae were considered aciculatans of uncertain position [11,22], closer relationships of the Nerillidae with any Aciculata have never been retrieved in the present analyses. In conclusion, it is an interesting working hypothesis that the "Archiannelida" (less Dinophilidae and possibly also Protodriloides) might form a clade [8], either basal, or close to the SSC [4,11,40-43].

Sipunculida and Dinophilidae

The Sipunculida have been considered a separate phylum by most authors [6]. Although they as adults do not exhibit any signs of segmentation, sipunculids show transitional stages of segmentation during development of their ventral nerve cord [44,45]. Analyses of the mitochondrial genomes [46-51] as well as phylogenomics [9] suggest that the sipunculids are an annelid ingroup. The unique correspondence of podocyte lining the metane-phridia appears to join Sipunculida with the Sabellida and Terebelliformia [33]. The pharyngeal apparatus in the sipunculid pelagosphere larvae is similar to that found in the Cirratuliformia [52]. In the present trees, the sipuncu-lids have always been found as an annelid ingroup, often close to the Dinophilidae. No closer relationships between sipunculids and terebelliformians, sabellidans, oweniids or chaetopterids [7,15,33] were recovered here. The enigmatic Dinophilidae were never found to be closely related to the Dorvilleidae (or Eunicida), i.e. in the placement expected by the morphologists [10,42,43]. Even Parapodrilus, one of the presumably dinophilid-related, "progenetic" dorvilleids [42], groups consistently with other dorvilleids within the Eunicida and shows no affinities to the Dinophilidae [25,35]. In the taxon-exclu-sion analyses, both Dinophilidae and Sipunculida group within the Cirratuliformia (uMP), or they are a sister group of the Pogonophora- "Clitellatomorpha" superc-lade (BI). In conclusion, close affinities between Dinophi-lidae and Sipunculida (supported by the shared loss of chaetae, parapodia, circulatory system, and palps) appear as a possible working hypothesis [8], but the precise position of this clade remains uncertain.


Relationships of the Myzostomida, traditionally regarded as aberrant polychaetes, probably phyllodocidans [6,10,11], are uncertain. They may be either basal proto-stomes [53] or platyzoan relatives [54,55]. Their highly aberrant sequences (the longest branch in the present BI tree) and uncertain homology of many morphological characters do not allow to deduce their relationships precisely. Nevertheless, the numerous annelid-like traits of the Myzostomida [56-58] still could corroborate some polychaete-myzostomid proximity. Recently, Bleidorn et al. [48,49] re-examined this issue by analysis of four nuclear genes and a mitochondrial genome and showed

myzostomids are likely part of the annelid radiation. On the contrary, in the comprehensive phylogenomic analysis of the Metazoa [9], Myzostomida are placed within the Platyzoa [cf. [49]]. In the present analyses, the myzostomids seem to belong into the Annelida (possibly into the Aciculata, as suggested by the taxon-exclusion tests). However, the taxon sampling in this paper does not provide a suitable test for alternative hypotheses due to lack of most of the possible non-annelid relatives of the Myzostomida (Gnathifera, Cycliophora, Platyhelminthes [9,53-55]).


The combined ("total-evidence") phylogenetic analysis provides a modified view of annelid evolution. Several higher-level annelid taxa are suggested by our analyses, i.e. Phyllodocida, Eunicida, the orbinioid-parergodrilid clade, Cirratuliformia, sabelloid-spionoid clade, terebel-loid-capitellid clade (incl. Echiurida), and "Clitellatomor-pha". Although the deep-level evolutionary relationships of Annelida remain poorly understood, we propose sister-group relationships between the Phyllodocida-Aberranta-Amphinomida and Eunicida-OPC clades, between the Cirratuliformia and SSC, and possibly also between the "Clitellatomorpha" and Oweniidae-Pogonophora clades. Two unorthodox clades, "core Archiannelida" and Dinophilidae-Sipunculida, are proposed here.


Taxonomy, datasets, and data combination

The present analysis was performed to include all available "families" of the Polychaeta, plus representatives of the Clitellata, Echiurida, Pogonophora, Myzostomida, Sipunculida, Mollusca, and Brachiopoda (as a rooting outgroup). All but a few annelid nominal "families" were included; the exceptions were enigmatic annelid (?) genera Lobatocerebrum, Diurodrilus [5], and Jennaria, "meso-zoan" Orthonectida (annelid affinities of which were proposed occasionally [59]), Hartmaniellidae (Eunicida), and several, mostly pelagic or parasitic subgroups of the Phyllodocida, for which no molecular data were available at time.

The morphological dataset (MOR) included 93 characters [see Additional Files 1 and 2]. The gene sequences were obtained from GenBank. The dataset included six genes, both protein-coding (EF1 a, histone H3, COI) and ribos-omal (18S, 28S, and 16S), and both nuclear (18S, 28S, EF1 a, H3) and mitochondrial (16S, COI) [see Additional File 3]. The protein-coding sequences were translated to amino acids and then aligned with CLUSTAL W under default settings for gap costs (gap opening penalty 10.00; gap extension penalty 0.20). The alignment of sequences for ribosomal RNA genes was conducted in the on-line version of MAFFT v6 in the E-INS-i mode [60]. Ambigu-

ous positions were excluded by using Gblocks [61]. The saturated positions were not excluded [8].

The combined datasets were completed by introducing question marks for the absent data partitions. To minimize the number of missing entries in the dataset, composite terminals were constructed with individual partitions from different species of a higher taxon, usually a nominal "family" [see Additional File 4].

Phylogenetic analyses

Bayesian phylogenetic analysis was conducted with a Metropolis-coupled Markov chain Monte Carlo algorithm [62] as implemented in MrBayes v3.1.2 [63]. MrModel-Test v2.2 [64], a simplified version of ModelTest 3.06 [65], and PAUP* v4.0b10 [66] were used to estimate GTR+A substitution model as the best-fitting for all molecular character sets (based on both AIC and hLRT criteria). Morphological characters were treated with the standard discrete model assuming gamma-shaped rate variation and variable coding bias. Model parameters were unlinked across partitions. Two independent runs of combined analysis with 10 Markov chains each were conducted for 10,000,000 generations with a sample frequency of 100 (heating 0.1). The first 61,000 trees from each run were discarded as burn-in; convergence between the two runs was estimated using diagnostics criteria produced by the "sump" command in MrBayes (PSRF [TL] = 1.001). The remaining 78,000 trees were used for reconstruction of a 50% majority-rule consensus tree. Testing the influence of burnin value on the consensual tree revealed a high stability of the tree topology within the whole investigated range (burn-in 61,000 to 90,000).

The maximum-parsimony (MP) analysis was applied to MOR, molecular, and combined data matrices (NONA version 2.0 [67]: heuristics, option "hold10000 mult*100 hold/100", unconstrained search strategy with TBR branch swapping). Bremer indices of branch support and bootstrap support values were calculated by NONA (options "bsupport10000" and "mult*100 hold/10" with 1,000 replications, respectively).

In addition, two "experimental" MP analyses (both in NONA: option "hold10000 mult*100 hold/100", unconstrained search strategy with TBR branch swapping) were performed.

(1) The modified "slow-fast" method [68] was used to remove characters (both morphological and molecular) that were supposedly responsible for stochastic information noise. New datasets were constructed, in which only characters with no observed variability within one of the 11 clades (that were present in both all-taxa BI and MP

trees) were included. The analysed clades included Cirrat-uliformia (5 spp.), SSC (6), "Clitellatomorpha" (7), TCC (10), Phyllodocida (19), Eunicida (6), Pogonophora (4), Amphinomida (2), Opheliidae-Scalibregmatidae (2), Sipunculida (2), and Mollusca (3). All reduced datasets were then concatenated to a new, weighted dataset, in which each original character received a weight from 11 to zero. Moreover, all characters with more than 80% of ambiguous (unknown, inapplicable, or polymorphic) character states were excluded. The final weighted dataset included 17,268 informative "characters" (601 morphological and 16,667 molecular).

(2) To test preliminary results suggesting that ingroup and outgroup taxon selection affected internal relationships, we ran several taxon-exclusion analyses. The "problematic" species/clades (Myzostomida, Opheliidae, Scalibreg-matidae, Fauveliopsidae, Serpulidae, Magelona, Oweniidae, Chaetopteridae, Pectinariidae, Apistobranchus, Polygordiidae, Saccocirridae, Protodrilidae, Protodriloides, Nerillidae, Histriobdellidae, Dinophilidae, Sipunculida, Hesionidae, Paralacydoniidae, Nepthyidae, Pilargidae, Tomopteridae, Aberranta, Amphinomida) were identified as follows: (i) their position was highly unstable in comparison of the all-taxa uMP, wMP and BI trees; and/or (ii) their presence/absence caused important topological changes; and/or (iii) they were consistenly placed basally or even outside the Annelida in both MP and BI trees. Topological effects of including/removing the species identified as "problematic" were tested by constructing an unrooted "backbone tree" to which individual "problematic" species were appended in one-by-one manner.

Authors' contributions

JZ conceived this study, performed the MP analyses, and wrote the manuscript. PR and JJ edited and aligned sequences and performed some analyses. LP performed most of the Bayesian analyses. All authors read and approved the final manuscript.

Additional material

Additional file 1

List of morphological characters. The data provided include an annotated list of 93 morphological characters used in the combined phyloge-netic analyses. Click here for file


Additional file 2

Morphological dataset. The data provided include a data matrix including 89 terminals (including outgroups, Psammodrilidae and Spinther) and 93 morphological characters. Click here for file


Additional file 3

Concatenated molecular dataset. The data provided include an alignment of the six concatenated molecular data partitions (COI, EF1 a, H3, 18S, 16S, 28S). Click here for file


Additional file 4

List of taxa used in the analyses. The data provided include a list of taxa used in the analyses with GenBank accession numbers. Click here for file



Our study was supported by research project No. MSM600766580I

(Czech Ministry of Education). Part of this work was carried out by using

the resources of the Computational Biology Service Unit from Cornell University (which is partially funded by Microsoft Corp.).


1. Erseus C, Kallersjö M: I8S rDNA phylogeny of Clitellata (Annelida). Zool Scr 2004, 33:187-196.

2. Marotta R, Ferraguti M, Erseus C, Gustavsson LM: Combined-data phylogenetics and character evolution of Clitellata (Annelida) using I8S rDNA and morphology. Zool J Linn Soc 2008, 154:1-26.

3. Rousset V, Plaisance L, Erseus C, Siddall ME, Rouse GW: Evolution of habitat preference in Clitellata (Annelida). Biol J Linn Soc 2008, 95:447-464.

4. Struck TH: Progenetic species in polychaetes (Annelida) and problems assessing their phylogenetic affiliation. Integr Comp Biol 2006, 46:558-568.

5. Worsaae K, Rouse GW: Is Diurodrilus an annelid? J Morphol 2008, 269:1426-1455.

6. Nielsen C: Animal Evolution Oxford: Oxford University Press; 2001.

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

8. Struck TH, Nesnidal MP, Purschke G, Halanych KM: Detecting possibly saturated positions in 18S and 28S sequences and their influence on phylogenetic reconstruction of Annelida (Lophotrochozoa). Mol Phylogen Evol 2008, 48:628-645.

9. Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, Sorensen MV, Haddock SHD, Schmidt-Rhaesa A, Okusu A, Kristensen RM, Wheeler WC, Martindale MQ, Giribet G: Broad phylogenomic sampling improves resolution of the animal tree of life. Nature 2008, 452:745-749.

10. Rouse GW, Fauchald K: Cladistics and polychaetes. Zool Scr 1997, 26:139-204.

11. Rouse GW, Pleijel F: Polychaetes London: Oxford University Press; 2001.

12. Bleidorn C, Vogt L, Bartolomaeus T: A contribution to sedentary polychaete phylogeny using I8S rRNA sequence data. J Zool Syst Evol Res 2003, 41:186-195.

13. Bleidorn C, Vogt L, Bartolomaeus T: New insights into polychaete phylogeny (Annelida) inferred from I8S rDNA sequences. Mol Phylogen Evol 2003, 29:279-288.

14. Hall KA, Hutchings PA, Colgan DJ: Further phylogenetic studies of the Polychaeta using I8S rDNA sequence data. J Mar Biol Ass UK 2004, 84:949-960.

15. Colgan DJ, Hutchings PA, Braune M: A multigene framework for polychaete phylogenetic studies. Organ Divers Evol 2006, 6:220-235.

16. Rousset V, Pleijel F, Rouse GW, Erseus C, Siddall ME: A molecular phylogeny of annelids. Cladistics 2007, 23:41-63.

17. Colgan DJ, Hutchings PA, Beacham E: Multi-gene analysis of the phylogenetic relationships among the Mollusca, Annelida, and Arthropoda. Zool Sci 2008, 47:338-351.

18. Gatesy J, O'Grady P, Baker RH: Corroboration among data sets in simultaneous analysis: hidden support for phylogenetic relationships among higher level artiodactyl taxa. Cladistics 1999, 15:271-313.

19. Rousset V, Rouse GW, Feral JP, Desbruyeres D, Pleijel F: Molecular and morphological evidence of Alvinellidae relationships (Terebelliformia, Polychaeta, Annelida). Zool Scr 2003, 32:185-197.

20. Rousset V, Rouse GW, Siddall ME, Tillier A, Pleijel F: The phylogenetic position of Siboglinidae (Annelida) inferred from I8S rRNA, 28S rRNA and morphological data. Cladistics 2004, 20:518-533.

21. Wiklund H, Nygren A, Pleijel F, Sundberg P: Phylogeny of Aphro-ditiformia (Polychaeta) based on molecular and morphological data. Mol Phylogen Evol 2005, 37:494-502.

22. Worsaae K, Nygren A, Rouse GW, Giribet G, Persson J, Sundberg P, Pleijel F: Phylogenetic position of Nerillidae and Aberranta (Polychaeta, Annelida), analysed by direct optimization of combined molecular and morphological data. Zool Scr 2005, 34:313-328.

23. Struck TH, Purschke G, Halanych KM: Phylogeny of Eunicida (Annelida) and exploring data congruence using a partition addition bootstrap alteration (PABA) approach. Syst Biol 2006, 55(1): 1-20.

24. Hoffmann S, Hausen H: Chaetal arrangement in Orbiniidae (Annelida, Polychaeta) and its significance for systematics. Zoomorphology 2007, 126:215-227.

25. Struck TH, Westheide W, Purschke G: Progenesis in Eunicida ("Polychaeta," Annelida) - separate evolutionary events? Evidence from molecular data. Mol Phylogen Evol 2002, 25:190-199.

26. Purschke G, Jördens J: Male genital organs in the eulittoral meiofaunal polychaete Stygocapitella subterranea (Annelida, Parergodrilidae): ultrastructure, functional and phyloge-netic significance. Zoomorphology 2008, 126:283-297.

27. Orrhage L, Müller MCM: Morphology of the nervous system of Polychaeta (Annelida). Hydrobiologia 2005, 535:79-111.

28. Hessling R, Westheide W: Are Echiura derived from a segmented ancestor? Immunohistochemical analysis of the nervous system in developmental stages of Bonellia viridis. J Morphol 2002, 252:100-113.

29. Kieselbach D, Hausen H: Chaetal arrangement provides no support for a close relationship of Sabellidae and Sabellariidae (Annelida). J Morphol 2008, 269:104-117.

30. Burnette AB, Struck TH, Halanych KM: Holopelagic Poeobius meseres ("Poeobiidae", Annelida) is derived from benthic fla-belligerid worms. Biol Bull 2005, 208:213-220.

31. Struck TH, Purschke G: The sister group relationship of Aeolo-somatidae and Potamodrilidae (Annelida: "Polychaeta") - a molecular phylogenetic approach based on I8S rDNA and cytochrome oxidase I. Zool Anz 2005, 243:281-293.

32. Hessling R, Purschke G: Immunohistochemical (cLSM) and ultrastructural analysis of the central nervous system and sense organs in Aeolosoma hemprichi (Annelida, Aeolosoma-tidae). Zoomorphology 2000, 120:65-78.

33. Bartolomaeus T, Quast B: Structure and development of nephridia in Annelida and related taxa. Hydrobiologia 2005, 535:139-165.

34. Martin P: On the origin of the Hirudinea and the demise of the Oligochaeta. Proc R Soc Lond B 2001, 268:1089-1098.

35. Struck TH, Halanych KM, Purschke G: Dinophilidae (Annelida) is most likely not a progenetic Eunicida: evidence from 18S and 28S rDNA. Mol Phylogen Evol 2005, 37:619-623.

36. Struck T, Hessling R, Purschke G: The phylogenetic position of the Aeolosomatidae and Parergodrilidae, two enigmatic oli-gochaete-like taxa of the ,Polychaeta", based on molecular data from I8S rDNA sequences. J Zool Syst Evol Res 2002, 40:155-163.

37. Jördens J, Struck T, Purschke G: Phylogenetic inference regarding Parergodrilidae and Hrabeiella periglandulata ("Polycha-eta", Annelida) based on I8S rDNA, 28S rDNA and COI sequences. J Zool Syst Evol Res 2004, 42:270-280.

38. Purschke G: Is Hrabeiella periglandulata (Annelida, "Polycha-eta") the sister group of Clitellata? Evidence from an ultrastructural analysis of the dorsal pharynx in H. periglan-dulata and Enchytraeus minutus (Annelida, Clitellata). Zoomor-phology 2003, 122:55-66.

39. Wilkens V, Purschke G: Central nervous system and sense organs, with special reference to photoreceptor-like sensory elements, in Polygordius appendiculatus (Annelida), an interstitial polychaete with uncertain phylogenetic affinities. Inver-tebr Biol 2009, 128:46-64.

40. Rouse GW: Trochophore concepts: ciliary bands and the evolution of larvae in spiralian Metazoa. Biol J Linn Soc 1999, 66:411-464.

41. Rouse GW: Bias? What bias? The evolution of downstream larval feeding in animals. Zool Scr 2000, 29:213-236.

42. Eibye-Jacobsen D, Kristensen RM: A new genus and species of Dorvilleidae (Annelida, Polychaeta) from Bermuda, with a phylogenetic analysis of Dorvilleidae, Iphitimidae and Dinophilidae. Zool Scr 1994, 23:107-131.

43. Müller MCM, Westheide W: Comparative analysis of the nervous systems in presumptive progenetic dinophilid and dorvil-leid polychaetes (Annelida) by immunohistochemistry and cLSM. Acta Zool 2002, 83:33-48.

44. Wanninger A, Koop D, Bromham L, Noonan E, Degnan BM: Nervous and muscle system development in Phascolion strombus (Sipuncula). Dev Gene Evol 2005, 215:509-518.

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

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

47. Jennings RM, Halanych KM: Mitochondrial genomes of Clymenella torquata (Maldanidae) and Riftia pachyptila (Sibo-glinidae): evidence for conserved gene order in Annelida. Mol Biol Evol 2005, 22:210-222.

48. Bleidorn C, Eeckhaut I, Podsiadlowski L, Schult N, McHugh Damhnait, Halanych 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.

49. Bleidorn C, Podsiadlowski L, Min Z, Eeckhaut I, Hartmann S, Halanych KM, Tiedemann R: On the phylogenetic position of Myzostomida: can 77 genes get it wrong? BMC Evol Biol 2009, 9:150.

50. Mwinyi A, Meyer A, Bleidorn C, Lieb B, Bartolomaeus T, Podsiadlowski L: Mitochondrial genome sequence and gene order of Sipunculus nudus give additional support for an inclusion of Sipuncula into Annelida. BMC Genom 2009, 10:27.

51. Shen X, Ma XY, Ren JF, Zhao FQ: A close phylogenetic relationship between Sipuncula and Annelida evidenced from the complete mitochondrial genome sequence of Phascolosoma esculenta. BMC Genom 2009, 10:136.

52. Tzetlin AB, Purschke G: Fine structure of the pharyngeal apparatus of the pelagosphera larva in Phascolosoma agassizii (Sipuncula) and its phylogenetic significance. Zoomorphology 2006, 125:109-117.

53. Passamaneck Y, Halanych KM: Lophotrochozoan phylogeny assessed with LSU and SSU data: evidence of lophophorate polyphyly. Mol Phylogen Evol 2006, 40:20-28.

54. Eeckhaut I, McHugh D, Mardulyn P, Tiedemann P, Monteyne D, Jangoux M, Milinkovitch MC: Myzostomida: a link between tro-chozoans and flatworms? Proc R Soc Lond B 2000, 267:1383-1392.

55. Zrzavy J, Hypsa V, Tietz DF: Myzostomida are not annelids: molecular and morphological support for a clade of animals with anterior sperm flagella. Cladistics 2001, 17(2): 170-198.

56. Eeckhaut I, Lanterbecq F: Myzostomida: a review of the phylog-eny and ultrastructure. Hydrobiologia 2005, 535:253-275.

57. Eeckhaut I, Fievez L, Müller MCM: Larval development of Myzos-toma cirriferum (Myzostomida). J Morphol 2003, 258:269-283.

58. Müller MCM, Westheide W: Structure of the nervous system of Myzostoma cirriferum (Annelida) as revealed by immunohistochemistry and cLSM analyses. J Morphol 2000, 245(2):87-98.

59. Slyusarev GS, Kristensen RM: structure of the ciliated cells and ciliary rootlets of Intoshia variabili Fine (Orthonectida). Zoomorphology 2003, 122:33-39.

60. Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform 2008, 9:286-298 [].

61. Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol 2000, 17:540-552 [ Gblocks.html].

62. Altekar G, Dwarkadas S, Huelsenbeck JP, Ronquist F: Parallel Metropolis-coupled Markov Chain Monte Carlo for Bayesian phylogenetic inference. Bioinformatics 2004, 20:407-415.

63. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19:1572-1574.

64. Nylander JAA: MrModelTest v2. [Computer Software]. Uppsala: Evolutionary Biology Centre, Uppsala University; 2004.

65. Posada D, Crandall KA: ModelTest: testing the model of DNA substitution. Bioinformatics 1998, 14:817-818.

66. Swofford DL: Paup*. Phylogenetic Analysis Using Parsimony (*and Other Methods), Version 4. Sunderland, Massachusetts: Sinauer Associates; 2003.

67. Goloboff P: NONA, Version 2.0. Computer software and manual. 1999 [].

68. Brinkmann H, Philippe H: Archaea sister group of bacteria? Indications from tree reconstruction artifacts in ancient phylog-enies. Mol Biol Evol 1999, 16:817-825.

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 ^