The Mitochondrial Genome of the Sipunculid Phascolopsis gouldii Supports Its Association with Annelida Rather than Mollusca

How Many Sipunculan Species Are Hiding in Our Oceans?

Sipuncula, long considered a separate phylum, are now commonly included in the Annelida based on phylogenomic analyses. The sipunculan body consists of an unsegmented trunk and a retractable introvert, usually with a set of tentacles at its anterior end. Unlike other annelids, they have no chaetae, but the introvert is often adorned with proteinaceous hooks that can be important taxonomic characters. Other external taxonomic characters include the tentacles (number, shape and arrangement), body papillae and, in some cases, hardened shields, as well as length ratios. Many species require dissection for correct identification to reveal internal characteristics, such as introvert retractor muscles, nephridia and contractile vessels. Here we summarize the state of the current knowledge of species diversity in sipunculans. We emphasize molecular studies, conducted over the past two decades, that have revealed multiple complexes of cryptic or pseudocryptic species. It has become obvious that diversity is significantly higher than the current taxonomic scheme accounts for, but formal species descriptions are lagging behind. Although the major branches in the sipunculan phylogeny have become increasingly consolidated, the internal relationships within most branches are still in flux.

Assessment of mitochondrial genomes for heterobranch gastropod phylogenetics

Background

Heterobranchia is a diverse clade of marine, freshwater, and terrestrial gastropod molluscs. It includes such disparate taxa as nudibranchs, sea hares, bubble snails, pulmonate land snails and slugs, and a number of (mostly small-bodied) poorly known snails and slugs collectively referred to as the “lower heterobranchs”. Evolutionary relationships within Heterobranchia have been challenging to resolve and the group has been subject to frequent and significant taxonomic revision. Mitochondrial (mt) genomes can be a useful molecular marker for phylogenetics but, to date, sequences have been available for only a relatively small subset of Heterobranchia.

Results

To assess the utility of mitochondrial genomes for resolving evolutionary relationships within this clade, eleven new mt genomes were sequenced including representatives of several groups of “lower heterobranchs”. Maximum likelihood analyses of concatenated matrices of the thirteen protein coding genes found weak support for most higher-level relationships even after several taxa with extremely high rates of evolution were excluded. Bayesian inference with the CAT + GTR model resulted in a reconstruction that is much more consistent with the current understanding of heterobranch phylogeny. Notably, this analysis recovered Valvatoidea and Orbitestelloidea in a polytomy with a clade including all other heterobranchs, highlighting these taxa as important to understanding early heterobranch evolution. Also, dramatic gene rearrangements were detected within and between multiple clades. However, a single gene order is conserved across the majority of heterobranch clades.

Conclusions

Analysis of mitochondrial genomes in a Bayesian framework with the site heterogeneous CAT + GTR model resulted in a topology largely consistent with the current understanding of heterobranch phylogeny. However, mitochondrial genomes appear to be too variable to serve as good phylogenetic markers for robustly resolving a number of deeper splits within this clade.

Comparative mitochondrial genomic analyses of three chemosynthetic vesicomyid clams from deep-sea habitats

Vesicomyid clams of the subfamily Pliocardinae are among the dominant chemosymbiotic bivalves found in sulfide-rich deep-sea habitats. Plastic morphologies and present molecular data could not resolve taxonomic uncertainties. The complete mitochondrial (mt) genomes will provide more data for comparative studies on molecular phylogeny and systematics of this taxonomically uncertain group, and help to clarify generic classifications. In this study, we analyze the features and evolutionary dynamics of mt genomes from three Archivesica species (Archivesica sp., Ar. gigas and Ar. pacifica) pertaining to subfamily Pliocardinae. Sequence coverage is nearly complete for the three newly sequenced mt genomes, with only the control region and some tRNA genes missing. Gene content, base composition, and codon usage are highly conserved in these pliocardiin species. Comparative analysis revealed the vesicomyid have a relatively lower ratio of Ka/Ks, and all 13 protein-coding genes (PGCs) are under strong purifying selection with a ratio of Ka/Ks far lower than one. Minimal changes in gene arrangement among vesicomyid species are due to the translocation trnaG in Isorropodon fossajaponicum. Additional tRNA genes were detected between trnaG and nad2 in Abyssogena mariana (trnaL3), Ab. phaseoliformis (trnaS3), and Phreagena okutanii (trnaM2), and display high similarity to other pliocardiin sequences at the same location. Single base insertion in multiple sites of this location could result in new tRNA genes, suggesting a possible tRNA arising from nongeneic sequence. Phylogenetic analysis based on 12 PCGs (excluding atp8) supports the monophyly of Pliocardiinae. These nearly complete mitogenomes provide relevant data for further comparative studies on molecular phylogeny and systematics of this taxonomically uncertain group of chemosymbiotic bivalves.

A complete logical approach to resolve the evolution and dynamics of mitochondrial genome in bilaterians

Investigating how recombination might modify gene order during the evolution has become a routine part of mitochondrial genome analysis. A new method of genomic maps analysis based on formal logic is described. The purpose of this method is to 1) use mitochondrial gene order of current taxa as datasets 2) calculate rearrangements between all mitochondrial gene orders and 3) reconstruct phylogenetic relationships according to these calculated rearrangements within a tree under the assumption of maximum parsimony. Unlike existing methods mainly based on the probabilistic approach, the main strength of this new approach is that it calculates all the exact tree solutions with completeness and provides logical consequences as highly robust results. Moreover, this method infers all possible hypothetical ancestors and reconstructs character states for all internal nodes of the trees. We started by testing our method using the deuterostomes as a study case. Then, with sponges as an outgroup, we investigated the evolutionary history of mitochondrial genomes of 47 bilaterian phyla and emphasised the peculiar case of chaetognaths. This pilot work showed that the use of formal logic in a hypothetico-deductive background such as phylogeny (where experimental testing of hypotheses is impossible) is very promising to explore mitochondrial gene order in deuterostomes and should be applied to many other bilaterian clades.

Re-evaluating the phylogeny of Sipuncula through transcriptomics

Sipunculans (also known as peanut worms) are an ancient group of exclusively marine worms with a global distribution and a fossil record that dates back to the Early Cambrian. The systematics of sipunculans, now considered a distinct subclade of Annelida, has been studied for decades using morphological and molecular characters, and has reached the limits of Sanger-based approaches. Here, we reevaluate their family-level phylogeny by comparative transcriptomic analysis of eight species representing all known families within Sipuncula. Two data matrices with alternative gene occupancy levels (large matrix with 675 genes and 62% missing data; reduced matrix with 141 genes and 23% missing data) were analysed using concatenation and gene-tree methods, yielding congruent results and resolving each internal node with maximum support. We thus corroborate prior phylogenetic work based on molecular data, resolve outstanding issues with respect to the familial relationships of Aspidosiphonidae, Antillesomatidae and Phascolosomatidae, and highlight the next area of focus for sipunculan systematics.

The complete mitochondrial genome of Solemya velum (Mollusca: Bivalvia) and its relationships with conchifera

Background

Bivalve mitochondrial genomes exhibit a wide array of uncommon features, like extensive gene rearrangements, large sizes, and unusual ways of inheritance. Species pertaining to the order Solemyida (subclass Opponobranchia) show many peculiar evolutionary adaptations, f.i. extensive symbiosis with chemoautotrophic bacteria. Despite Opponobranchia are central in bivalve phylogeny, being considered the sister group of all Autobranchia, a complete mitochondrial genome has not been sequenced yet.

Results

In this paper, we characterized the complete mitochondrial genome of the Atlantic awning clam Solemya velum: A-T content, gene arrangement and other features are more similar to putative ancestral mollusks than to other bivalves. Two supranumerary open reading frames are present in a large, otherwise unassigned, region, while the origin of replication could be located in a region upstream to the cox3 gene.

Conclusions

We show that S. velum mitogenome retains most of the ancestral conchiferan features, which is unusual among bivalve mollusks, and we discuss main peculiarities of this first example of an organellar genome coming from the subclass Opponobranchia. Mitochondrial genomes of Solemya (for bivalves) and Haliotis (for gastropods) seem to retain the original condition of mollusks, as most probably exemplified by Katharina.

Mitogenomics does not resolve deep molluscan relationships (yet?)

The origin of molluscs among lophotrochozoan metazoans is unresolved and interclass relationships are contradictory between morphology-based, multi-locus, and recent phylogenomic analyses. Within the "Deep Metazoan Phylogeny" framework, all available molluscan mitochondrial genomes were compiled, covering 6 of 8 classes. Genomes were reannotated, and 13 protein coding genes (PCGs) were analyzed in various taxon settings, under multiple masking and coding regimes. Maximum Likelihood based methods were used for phylogenetic reconstructions. In all cases, molluscs result mixed up with lophotrochozoan outgroups, and most molluscan classes with more than single representatives available are non-monophyletic. We discuss systematic errors such as long branch attraction to cause aberrant, basal positions of fast evolving ingroups such as scaphopods, patellogastropods and, in particular, the gastropod subgroup Heterobranchia. Mitochondrial sequences analyzed either as amino acids or nucleotides may perform well in some (Cephalopoda) but not in other palaeozoic molluscan groups; they are not suitable to reconstruct deep (Cambrian) molluscan evolution. Supposedly "rare" mitochondrial genome level features have long been promoted as phylogenetically informative. In our newly annotated data set, features such as genome size, transcription on one or both strands, and certain coupled pairs of PCGs show a homoplastic, but obviously non-random distribution. Apparently congruent (but not unambiguous) signal for non-trivial subclades, e.g. for a clade composed of pteriomorph and heterodont bivalves, needs confirmation from a more comprehensive bivalve sampling. We found that larger clusters not only of PCGs but also of rRNAs and even tRNAs can bear local phylogenetic signal; adding trnG-trnE to the end of the ancestral cluster trnM-trnC-trnY-trnW-trnQ might be synapomorphic for Mollusca. Mitochondrial gene arrangement and other genome level features explored and reviewed herein thus failed as golden bullets, but are promising as additional characters or evidence supporting deep molluscan clades revealed by other data sets. A representative and dense sampling of molluscan subgroups may contribute to resolve contentious interclass relationships in the future, and is vital for exploring the evolution of especially diverse mitochondrial genomes in molluscs.

Spiral cleavage and early embryology of a loxosomatid entoproct and the usefulness of spiralian apical cross patterns for phylogenetic inferences

Background

Among the four major bilaterian clades, Deuterostomia, Acoelomorpha, Ecdysozoa, and Lophotrochozoa, the latter shows an astonishing diversity of bodyplans. While the largest lophotrochozoan assemblage, the Spiralia, which at least comprises Annelida, Mollusca, Entoprocta, Platyhelminthes, and Nemertea, show a spiral cleavage pattern, Ectoprocta, Brachiopoda and Phoronida (the Lophophorata) cleave radially. Despite a vast amount of recent molecular phylogenetic analyses, the interrelationships of lophotrochozoan phyla remain largely unresolved. Thereby, Entoprocta play a key role, because they have frequently been assigned to the Ectoprocta, despite their differently cleaving embryos. However, developmental data on entoprocts employing modern methods are virtually non-existent and the data available rely exclusively on sketch drawings, thus calling for thorough re-investigation.

Results

By applying fluorescence staining in combination with confocal microscopy and 3D-imaging techniques, we analyzed early embryonic development of a basal loxosomatid entoproct. We found that cleavage is asynchronous, equal, and spiral. An apical rosette, typical for most spiralian embryos, is formed. We also identified two cross-like cellular arrangements that bear similarities to both, a "molluscan-like" as well as an "annelid-like" cross, respectively.

Conclusions

A broad comparison of cleavage types and apical cross patterns across Lophotrochozoa shows high plasticity of these character sets and we therefore argue that these developmental traits should be treated and interpreted carefully when used for phylogenetic inferences.

Detecting the symplesiomorphy trap: a multigene phylogenetic analysis of terebelliform annelids

BACKGROUND

For phylogenetic reconstructions, conflict in signal is a potential problem for tree reconstruction. For instance, molecular data from different cellular components, such as the mitochondrion and nucleus, may be inconsistent with each other. Mammalian studies provide one such case of conflict where mitochondrial data, which display compositional biases, support the Marsupionta hypothesis, but nuclear data confirm the Theria hypothesis. Most observations of compositional biases in tree reconstruction have focused on lineages with different composition than the majority of the lineages under analysis. However in some situations, the position of taxa that lack compositional bias may be influenced rather than the position of taxa that possess compositional bias. This situation is due to apparent symplesiomorphic characters and known as "the symplesiomorphy trap".

RESULTS

Herein, we report an example of the sympleisomorphy trap and how to detect it. Worms within Terebelliformia (sensu Rouse & Pleijel 2001) are mainly tube-dwelling annelids comprising five 'families': Alvinellidae, Ampharetidae, Terebellidae, Trichobranchidae and Pectinariidae. Using mitochondrial genomic data, as well as data from the nuclear 18S, 28S rDNA and elongation factor-1α genes, we revealed incongruence between mitochondrial and nuclear data regarding the placement of Trichobranchidae. Mitochondrial data favored a sister relationship between Terebellidae and Trichobranchidae, but nuclear data placed Trichobranchidae as sister to an Ampharetidae/Alvinellidae clade. Both positions have been proposed based on morphological data.

CONCLUSIONS

Our investigation revealed that mitochondrial data of Ampharetidae and Alvinellidae exhibited strong compositional biases. However, these biases resulted in a misplacement of Trichobranchidae, rather than Alvinellidae and Ampharetidae. Herein, we document that Trichobranchidae was apparently caught in the symplesiomorphy trap suggesting that in certain situations even homologies can be misleading.

Evidence for inter- and intra-clade recombinations in rabies virus

Homologous recombination is considered rare in negative-strand RNA viruses and has not been reported for rabies virus. In this study, full-length genomes of 44 rabies virus strains were analyzed for potential homologous recombination events. Phylogenetic analysis classified these strains into three clades. By applying six different recombination detection methods, one inter-clade and one intra-clade potential recombination events were identified with high confidence values. Software-predicted recombination break points of the two events were all located within the polymerase gene. This report presents the first evidence suggesting the possibility of homologous recombination in rabies virus, which could provide valuable insights for understanding the diversity and evolution of rabies virus as well as other negative-strand RNA viruses.

Species identification of the Northern shrimp (Pandalus borealis) by polymerase chain reaction-restriction fragment length polymorphism and proteomic analysis

Genomic and proteomic techniques for species identification of meat and seafood products are being widely used. In this study, a genomic approach was used to differentiate Pandalus borealis (the Northern shrimp), which belongs to the superfamily Pandaloidea, from 30 crustaceans consisting of 19 commercially relevant prawns/shrimps species that belong to the superfamily Penaeoidea, which include the families Penaeidae and Solenoceridae, and 11 other crustacean species, including prawns, shrimps, lobsters, and crabs. For this purpose, a polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method was designed based on the amplification of the 16S rRNA/tRNA(Val)/12S rRNA mitochondrial regions using the primers 16S-CruF and 16S-CruR. The 966-bp PCR products were produced and cleaved with the restriction enzymes AluI, TaqI, and HinfI, which provided species-specific restriction patterns. In addition, a proteomic approach, based on matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) and electrospray ionization-ion trap (ESI-IT) mass spectrometry, was used to identify and characterize new P. borealis-specific peptides that could be useful as potential markers of this species in protein-based detection methods. To our knowledge, this is the first time a molecular method has been successfully applied to identify a wide range of prawn and shrimp species, including P. borealis, for either whole individuals or processed products. However, validation of the methods proposed here is required by applying them to a larger sample of individuals from different populations and geographic origins in order to avoid mainly false-negative results.

Sipunculans and segmentation

Comparative molecular, developmental and morphogenetic analyses show that the three major segmented animal groups-Lophotrochozoa, Ecdysozoa and Vertebrata-use a wide range of ontogenetic pathways to establish metameric body organization. Even in the life history of a single specimen, different mechanisms may act on the level of gene expression, cell proliferation, tissue differentiation and organ system formation in individual segments. Accordingly, in some polychaete annelids the first three pairs of segmental peripheral neurons arise synchronously, while the metameric commissures of the ventral nervous system form in anterior-posterior progression. Contrary to traditional belief, loss of segmentation may have occurred more often than commonly assumed, as exemplified in the sipunculans, which show remnants of segmentation in larval stages but are unsegmented as adults. The developmental plasticity and potential evolutionary lability of segmentation nourishes the controversy of a segmented bilaterian ancestor versus multiple independent evolution of segmentation in respective metazoan lineages.

The mitochondrial genome of the soybean cyst nematode, Heterodera glycines

We sequenced the entire coding region of the mitochondrial genome of Heterodera glycines. The sequence obtained comprised 14.9 kb, with PCR evidence indicating that the entire genome comprised a single, circular molecule of approximately 21-22 kb. The genome is the most T-rich nematode mitochondrial genome reported to date, with T representing over half of all nucleotides on the coding strand. The genome also contains the highest number of poly(T) tracts so far reported (to our knowledge), with 60 poly(T) tracts ≥ 12 Ts. All genes are transcribed from the same mitochondrial strand. The organization of the mitochondrial genome of H. glycines shows a number of similarities compared with Radopholus similis, but fewer similarities when compared with Meloidogyne javanica. Very few gene boundaries are shared with Globodera pallida or Globodera rostochiensis. Partial mitochondrial genome sequences were also obtained for Heterodera cardiolata (5.3 kb) and Punctodera chalcoensis (6.8 kb), and these had identical organizations compared with H. glycines. We found PCR evidence of a minicircular mitochondrial genome in P. chalcoensis, but at low levels and lacking a noncoding region. Such circularised genome fragments may be present at low levels in a range of nematodes, with multipartite mitochondrial genomes representing a shift to a condition in which these subgenomic circles predominate.

Determination of the complete mitochondrial DNA sequence of Octopus minor

In this study, we have determined the complete nucleotide sequence of the mitochondrial genome of Octopus minor. It is 15,974 nucleotide pairs and encodes 13 proteins, two ribosomal RNAs and 22 tRNAs of the mitochondrion's own protein synthesizing system. Seven of thirteen proteins are encoded by the H-strand, while the other six proteins, as well as the two ribosomal RNAs are encoded by the L-strand. The nucleotide composition of the proteins showed a nucleotide bias against G encoded by the H-strand, while they showed a nucleotide bias against A and C encoded by the L-strand. Two of the 13 protein coding genes of O. minor began with the unorthodox translation initiation codon ATA and all others use the standard ATG. In addition, six of thirteen mt proteins of O. minor have unambiguous termination codons. There are four cases where tRNA genes appear to overlap. The long noncoding region (LNCR) of O. minor was 930 nucleotides and no repeated sequences were found in this LNCR. The gene arrangements of O. minor showed remarkable similarity to that of O. ocellatus and O. vulgaris. Phylogenetic analysis demonstrated that O. minor appears as sister taxan to the monophyletic group combined by O. ocellatus and O. vulgaris, suggesting a relative distant genetic relationship between O. minor and the other two octopus species.

Cellular and muscular growth patterns during sipunculan development

Sipuncula is a lophotrochozoan taxon with annelid affinities, albeit lacking segmentation of the adult body. Here, we present data on cell proliferation and myogenesis during development of three sipunculan species, Phascolosoma agassizii, Thysanocardia nigra, and Themiste pyroides. The first anlagen of the circular body wall muscles appear simultaneously and not subsequently as in the annelids. At the same time, the rudiments of four longitudinal retractor muscles appear. This supports the notion that four introvert retractors were part of the ancestral sipunculan bodyplan. The longitudinal muscle fibers form a pattern of densely arranged fibers around the retractor muscles, indicating that the latter evolved from modified longitudinal body wall muscles. For a short time interval, the distribution of S-phase mitotic cells shows a metameric pattern in the developing ventral nerve cord during the pelagosphera stage. This pattern disappears close to metamorphic competence. Our findings are congruent with data on sipunculan neurogenesis, as well as with recent molecular analyses that place Sipuncula within Annelida, and thus strongly support a segmental ancestry of Sipuncula.

Genetic characterization of ticks from southwestern Romania by sequences of mitochondrial cox1 and nad5 genes

In the present study, samples representing three hard tick species and one soft tick species, namely Dermacentor marginatus, Haemaphysalis punctata, Ixodes ricinus and Argas persicus from southwestern Romania, and one hard tick, Haemaphysalis longicornis, from China were characterized genetically by a portion of mitochondrial cytochrome c oxidase subunit 1 gene (pcox1) and a portion of nicotinamide adenine dinucleotide dehydrogenase subunit 5 gene (pnad5). The pcox1 and pnad5 were amplified separately from individual ticks by PCR, sequenced and analyzed. The length of pcox1 and pnad5 sequences of all samples was 732 and 519 bp, respectively. The intra-specific sequence variation in De. marginatus was 0.1-1.0% for pcox1 and 0.2-1.2% for pnad5, whereas in Ha. punctata it was 0.4-1.9% for pcox1 and 0.4-1.0% for pnad5. For the tick species examined in the present study, sequence comparison revealed that the inter-specific sequence differences were higher: 15.9-27.6% for pcox1 and 20.3-42.4% for pnad5. This suggests that the cox1 and nad5 sequences could provide useful genetic markers for the specific identification and genetic characterization of ticks in Romania and elsewhere.

Expression of FoxA and GATA transcription factors correlates with regionalized gut development in two lophotrochozoan marine worms: Chaetopterus (Annelida) and Themiste lageniformis (Sipuncula)

Background

A through gut is present in almost all metazoans, and most likely represents an ancient innovation that enabled bilaterian animals to exploit a wide range of habitats. Molecular developmental studies indicate that Fox and GATA regulatory genes specify tissue regions along the gut tube in a broad diversity of taxa, although little is known about gut regionalization within the Lophotrochozoa. In this study, we isolated FoxA and GATA456 orthologs and used whole mount in situ hybridization during larval gut formation in two marine worms: the segmented, polychaete annelid Chaetopterus, which develops a planktotrophic larva with a tripartite gut, and the non-segmented sipunculan Themiste lageniformis, which develops a lecithotrophic larva with a U-shaped gut.

Results

FoxA and GATA456 transcripts are predominantly restricted to gut tissue, and together show regional expression spanning most of the alimentary canal in each of these lophotrochozoans, although neither FoxA nor GATA456 is expressed in the posterior intestine of Chaetopterus. In both species, FoxA is expressed at the blastula stage, transiently in presumptive endoderm before formation of a definitive gut tube, and throughout early larval development in discrete foregut and hindgut domains. GATA456 genes are expressed during endoderm formation, and in endoderm and mesoderm associated with the midgut in each species. Several species-specific differences were detected, including an overlap of FoxA and GATA456 expression in the intestinal system of Themiste, which is instead complimentary in Chaetopterus. Other differences include additional discrete expression domains of FoxA in ectodermal trunk cells in Themiste but not Chaetopterus, and expression of GATA456 in anterior ectoderm and midgut cells unique to Chaetopterus.

Conclusions

This study of gene expression in a sipunculan contributes new comparative developmental insights from lophotrochozoans, and shows that FoxA and GATA456 transcription factors are part of an ancient patterning mechanism that was deployed during early evolution of the metazoan through gut. The common utilization of FoxA and GATA456 throughout gut formation by species with contrasting life history modes indicates that both genes are core components of a gut-specific gene regulatory network in spiralians. Despite a highly conserved pattern of early development, and probably similar ontogenic origins of gut tissue, there are molecular differences in gut regionalization between lophotrochozoan species.

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

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 (18S, 28S, and 16S rRNA, EF1alpha, 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.

Identification of complete mitochondrial genome of the tufted deer

The tufted deer Elaphodus cephalophus are endangered animals in the world and little is understood about their mitochondrial (mt) genome. In our study, the mt genome of the tufted deer is identified--which is about 16 kb in length and contains 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and a non-coding sequence (control region). The distinguishing feature is that GTG is the start codon of the NADH4L gene and the cyt b gene has a subterminal AAA followed by the stop codon TAG. According to 12 H strand protein-coding genes and phylogenetic analysis, Elaphodus may have a sister relationship with another deer group Muntiacus.

Shaping the things to come: ontogeny of lophotrochozoan neuromuscular systems and the tetraneuralia concept

Despite the large variation in adult bodyplan phenotypes, a worm-shaped morphology is considered plesiomorphic for both Lophotrochozoa and Bilateria. Although almost all larval and adult lophotrochozoan worms have serially arranged ring muscles in their body wall, a comparison of their ontogeny reveals no less than six different developmental pathways that lead to this homogenous arrangement of ring muscles. However, in all taxa, with the exception of chaetodermomorph molluscs and the segmented annelids, ring muscle development starts with synchronous formation of certain pioneer myocytes, which is thus considered basal for Lophotrochozoa. Recent studies on spiralian neurogenesis revealed remnants of ancestral segmentation in echiurans and sipunculans, thus confirming molecular phylogenetic studies that propose a close relationship of these three taxa. Larval entoprocts exhibit a mosaic of larval and adult molluscan characters and, among other apomorphies, share with polyplacophoran Mollusca a complex larval apical organ and a tetraneurous nervous system, strongly suggesting a monophyletic assemblage of Entoprocta and Mollusca. The term Tetraneuralia is proposed herein for this lophotrochozoan clade. Overall, formation of the lophotrochozoan neuromuscular bodyplan appears as a highly dynamic process on both the ontogenetic and the evolutionary timescales, highlighting the importance of insights into these processes for reconstructing ancestral bodyplan features and phylogenetic relationships.

A close phylogenetic relationship between Sipuncula and Annelida evidenced from the complete mitochondrial genome sequence of Phascolosoma esculenta

Background

There are many advantages to the application of complete mitochondrial (mt) genomes in the accurate reconstruction of phylogenetic relationships in Metazoa. Although over one thousand metazoan genomes have been sequenced, the taxonomic sampling is highly biased, left with many phyla without a single representative of complete mitochondrial genome. Sipuncula (peanut worms or star worms) is a small taxon of worm-like marine organisms with an uncertain phylogenetic position. In this report, we present the mitochondrial genome sequence of Phascolosoma esculenta, the first complete mitochondrial genome of the phylum.

Results

The mitochondrial genome of P.esculenta is 15,494 bp in length. The coding strand consists of 32.1% A, 21.5% C, 13.0% G, and 33.4% T bases (AT = 65.5%; AT skew = -0.019; GC skew = -0.248). It contains thirteen protein-coding genes (PCGs) with 3,709 codons in total, twenty-two transfer RNA genes, two ribosomal RNA genes and a non-coding AT-rich region (AT = 74.2%). All of the 37 identified genes are transcribed from the same DNA strand. Compared with the typical set of metazoan mt genomes, sipunculid lacks trnR but has an additional trnM. Maximum Likelihood and Bayesian analyses of the protein sequences show that Myzostomida, Sipuncula and Annelida (including echiurans and pogonophorans) form a monophyletic group, which supports a closer relationship between Sipuncula and Annelida than with Mollusca, Brachiopoda, and some other lophotrochozoan groups.

Conclusion

This is the first report of a complete mitochondrial genome as a representative within the phylum Sipuncula. It shares many more similar features with the four known annelid and one echiuran mtDNAs. Firstly, sipunculans and annelids share quite similar gene order in the mitochondrial genome, with all 37 genes located on the same strand; secondly, phylogenetic analyses based on the concatenated protein sequences also strongly support the sipunculan + annelid clade (including echiurans and pogonophorans). Hence annelid "key-characters" including segmentation may be more labile than previously assumed.

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

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.

Lophotrochozoa internal phylogeny: new insights from an up-to-date analysis of nuclear ribosomal genes

Resolving the relationships among animal phyla is a key biological problem that remains to be solved. Morphology is unable to determine the relationships among most phyla and although molecular data have unveiled a new evolutionary scenario, they have their own limitations. Nuclear ribosomal genes (18S and 28S rDNA) have been used effectively for many years. However, they are considered of limited use for resolving deep divergences such as the origin of the bilaterians, due to certain drawbacks such as the long-branch attraction (LBA) problem. Here, we attempt to overcome these pitfalls by combining several methods suggested in previous studies and routinely used in contemporary standard phylogenetic analyses but that have not yet been applied to any bilaterian phylogeny based on these genes. The methods used include maximum likelihood and Bayesian inference, the application of models with rate heterogeneity across sites, wide taxon sampling and compartmentalized analyses for each problematic clade. The results obtained show that the combination of the above-mentioned methodologies minimizes the LBA effect, and a new Lophotrochozoa phylogeny emerges. Also, the Acoela and Nemertodermatida are confirmed with maximum support as the first branching bilaterians. Ribosomal RNA genes are thus a reliable source for the study of deep divergences in the metazoan tree, provided that the data are treated carefully.

Beyond linear sequence comparisons: the use of genome-level characters for phylogenetic reconstruction

The first whole genomes to be compared for phylogenetic inference were those of mitochondria, which provided the first sets of genome-level characters for phylogenetic reconstruction. Most powerful among these characters has been the comparisons of the relative arrangements of genes, which has convincingly resolved numerous branch points, including those that had remained recalcitrant even to very large molecular sequence comparisons. Now the world faces a tsunami of complete nuclear genome sequences. In addition to the tremendous amount of DNA sequence that is becoming available for comparison, there is also a potential for many more genome-level characters to be developed, including the relative positions of introns, the domain structures of proteins, gene family membership, the presence of particular biochemical pathways, aspects of DNA replication or transcription, and many others. These characters can be especially convincing owing to their low likelihood of reverting to a primitive condition or occurring independently in separate lineages, thereby reducing the occurrence of homoplasy. The comparisons of organelle genomes pioneered the way for using such features for phylogenetic reconstructions, and it is almost certainly true, as ever more genomic sequence becomes available, that further use of genome-level characters will play a big role in outlining the relationships among major animal groups.

Segmental mode of neural patterning in sipuncula

Recent molecular phylogenetic analyses suggest a close relationship between two worm-shaped phyla, the nonsegmented Sipuncula (peanut worms) and the segmented Annelida (e.g., earthworms and polychaetes) [1-5]. The striking differences in their bodyplans are exemplified by the annelids' paired, ladder-like ventral nervous system, which contains segmentally arranged ganglia, and the sipunculans' single ventral nerve cord (VNC), which is devoid of any segmental structures [6, 7]. Investigating central nervous system (CNS) formation with serotonin and FMRFamide labeling in a representative sipunculan, Phascolosoma agassizii, we found that neurogenesis initially follows a segmental pattern similar to that of annelids. Starting out with paired FMRFamidergic and serotonergic axons, four pairs of associated serotonergic perikarya and interconnecting commissures form one after another in an anterior-posterior progression. In late-stage larvae, the two serotonergic axons of the VNCs fuse, the commissures disappear, and one additional pair of perikarya is formed. These cells (ten in total) migrate toward one another, eventually forming two clusters of five cells each. These neural-remodeling processes result in the single nonmetameric CNS of the adult sipunculan. Our data confirm the segmental ancestry of Sipuncula and render Phascolosoma a textbook example for the Haeckelian hypothesis of ontogenetic recapitulation of the evolutionary history of a species [8].

Phylogenetic information from three mitochondrial genomes of Terebelliformia (Annelida) worms and duplication of the methionine tRNA

Mitochondrial genomes have been useful for inferring animal phylogeny across a wide range of clades, however they are still poorly sampled in some animal taxa, limiting our knowledge of mtDNA evolution. For example, despite being one of the most diverse animal phyla, only 5 complete annelid mitochrondial genomes have been published. To address this paucity of information, we obtained complete mitochondrial genomic sequences from Pista cristata (Terebellidae) and Terebellides stroemi (Trichobranchidae) as well as one nearly complete mitochondrial genome from Eclysippe vanelli (Ampharetidae). These taxa are within Terebelliformia (Annelida), which include spaghetti worms, icecream cone worms and their relatives. In contrast to the 37 genes found in most bilaterian metazoans, we recover 38 genes in the mitochondrial genomes of T. stroemi and P. cristata due to the presence of a second methionine tRNA (trnM). Interestingly, the two trnMs are located next to each other and are possibly a synapomorphy of these two taxa. The E. vanelli partial mitochondrial genome lacks this additional trnM at the same position, but it may be present in the region not sampled. Compared to other annelids, gene orders of these three mitochondrial genomes are generally conserved except for the atp6-mSSU region. Phylogenetic analyses reveal that mtDNA data strongly supports a Trichobranchidae/Terebellidae clade.

Annelid phylogeny and the status of Sipuncula and Echiura

BACKGROUND

Annelida comprises an ancient and ecologically important animal phylum with over 16,500 described species and members are the dominant macrofauna of the deep sea. Traditionally, two major groups are distinguished: Clitellata (including earthworms, leeches) and "Polychaeta" (mostly marine worms). Recent analyses of molecular data suggest that Annelida may include other taxa once considered separate phyla (i.e., Echiura, and Sipuncula) and that Clitellata are derived annelids, thus rendering "Polychaeta" paraphyletic; however, this contradicts classification schemes of annelids developed from recent analyses of morphological characters. Given that deep-level evolutionary relationships of Annelida are poorly understood, we have analyzed comprehensive datasets based on nuclear and mitochondrial genes, and have applied rigorous testing of alternative hypotheses so that we can move towards the robust reconstruction of annelid history needed to interpret animal body plan evolution.

RESULTS

Sipuncula, Echiura, Siboglinidae, and Clitellata are all nested within polychaete annelids according to phylogenetic analyses of three nuclear genes (18S rRNA, 28S rRNA, EF1alpha; 4552 nucleotide positions analyzed) for 81 taxa, and 11 nuclear and mitochondrial genes for 10 taxa (additional: 12S rRNA, 16S rRNA, ATP8, COX1-3, CYTB, NAD6; 11,454 nucleotide positions analyzed). For the first time, these findings are substantiated using approximately unbiased tests and non-scaled bootstrap probability tests that compare alternative hypotheses. For echiurans, the polychaete group Capitellidae is corroborated as the sister taxon; while the exact placement of Sipuncula within Annelida is still uncertain, our analyses suggest an affiliation with terebellimorphs. Siboglinids are in a clade with other sabellimorphs, and clitellates fall within a polychaete clade with aeolosomatids as their possible sister group. None of our analyses support the major polychaete clades reflected in the current classification scheme of annelids, and hypothesis testing significantly rejects monophyly of Scolecida, Palpata, Canalipalpata, and Aciculata.

CONCLUSION

Using multiple genes and explicit hypothesis testing, we show that Echiura, Siboglinidae, and Clitellata are derived annelids with polychaete sister taxa, and that Sipuncula should be included within annelids. The traditional composition of Annelida greatly underestimates the morphological diversity of this group, and inclusion of Sipuncula and Echiura implies that patterns of segmentation within annelids have been evolutionarily labile. Relationships within Annelida based on our analyses of multiple genes challenge the current classification scheme, and some alternative hypotheses are provided.

The complete mitochondrial genome of the black coral Chrysopathes formosa (Cnidaria:Anthozoa:Antipatharia) supports classification of antipatharians within the subclass Hexacorallia

Black corals comprise a globally distributed shallow- and deep-water taxon whose phylogenetic position within the Anthozoa has been debated. We sequenced the complete mitochondrial genome of the antipatharian Chrysopathes formosa to further evaluate its phylogenetic relationships. The circular mitochondrial genome (18,398 bp) consists of 13 energy pathway protein-coding genes and two ribosomal RNAs, but only two transfer RNA genes (trnM and trnW), as well as a group I intron within the nad5 gene that contains the only copies of nad1 and nad3. No novel genes were found in the antipatharian mitochondrial genome. Gene order and genome content are most similar to those of the sea anemone Metridium senile (subclass Hexacorallia), with differences being the relative location of three contiguous genes (cox2-nad4-nad6) and absence (from the antipatharian) of a group I intron within the cox1 gene. Phylogenetic analyses of multiple protein-coding genes support classifying the Antipatharia within the subclass Hexacorallia and not the subclass Ceriantipatharia; however, the sister-taxon relationships of black corals within Hexacorallia remain inconclusive.

Mitochondrial genomes of parasitic arthropods: implications for studies of population genetics and evolution

Over 39000 species of arthropods parasitize humans, domestic animals and wildlife. Despite their medical, veterinary and economic importance, most aspects of the population genetics and evolution of the vast majority of parasitic arthropods are poorly understood. Mitochondrial genomes are a rich source of markers for studies of population genetics and evolution. These markers include (1) nucleotide sequences of each of the 37 mitochondrial genes and non-coding regions; (2) concatenated nucleotide sequences of 2 or more genes; and (3) genomic features, such as gene duplications, gene rearrangements, and changes in gene content and secondary structures of RNAs. To date, the mitochondrial genomes of over 700 species of multi-cellular animals have been sequenced entirely, however, only 24 of these species are parasitic arthropods. Of the mitochondrial genome markers, only the nucleotide sequences of 4 mitochondrial genes,cox1,cob,rrnSandrrnL, have been well explored in population genetic and evolutionary studies of parasitic arthropods whereas the sequences of the other 33 genes, and various genomic features have not. We review current knowledge of the mitochondrial genomes of parasitic arthropods, summarize applications of mitochondrial genes and genomic features in population genetic and evolutionary studies, and highlight prospects for future research.

The development of the larval nervous system, musculature and ciliary bands of Pomatoceros lamarckii (Annelida): heterochrony in polychaetes

Background

To understand the evolution of animals it is essential to have taxon sampling across a representative spread of the animal kingdom. With the recent rearrangement of most of the Bilateria into three major clades (Ecdysozoa, Lophotrochozoa and Deuterostomia) it has become clear that the Lophotrochozoa are relatively poorly represented in our knowledge of animal development, compared to the Ecdysozoa and Deuterostomia. We aim to contribute towards redressing this balance with data on the development of the muscular, nervous and ciliary systems of the annelid Pomatoceros lamarckii (Serpulidae). We compare our data with other lophotrochozoans.

Results

P. lamarckii develops locomotory and feeding structures that enable it to become a swimming, planktotrophic larva within 24 hours. Formation of the trochophore includes development of a prototroch, metatroch and neurotroch, development of apical and posterior nervous elements at similar times, and development of musculature around the ciliary bands and digestive tract prior to development of any body wall muscles. The adult nervous and muscular systems are essentially preformed in the late larva. Interestingly, the muscular systems of the larvae and juvenile worms do not include the circular muscles of the body wall, which are considered to be plesiomorphic for annelids, although the possibility that circular muscles develop after these stages cannot be ruled out at this point.

Conclusion

A comparison between polychaetes shows variability in the timing (heterochrony) of development of body wall muscles and elements of the nervous system. These heterochronies are one route for evolution of different life history strategies, such as adaptations to feeding requirements.

The complete sequences and gene organisation of the mitochondrial genomes of the heterodont bivalves Acanthocardia tuberculata and Hiatella arctica--and the first record for a putative Atpase subunit 8 gene in marine bivalves

BACKGROUND

Mitochondrial (mt) gene arrangement is highly variable among molluscs and especially among bivalves. Of the 30 complete molluscan mt-genomes published to date, only one is of a heterodont bivalve, although this is the most diverse taxon in terms of species numbers. We determined the complete sequence of the mitochondrial genomes of Acanthocardia tuberculata and Hiatella arctica, (Mollusca, Bivalvia, Heterodonta) and describe their gene contents and genome organisations to assess the variability of these features among the Bivalvia and their value for phylogenetic inference.

RESULTS

The size of the mt-genome in Acanthocardia tuberculata is 16.104 basepairs (bp), and in Hiatella arctica 18.244 bp. The Acanthocardia mt-genome contains 12 of the typical protein coding genes, lacking the Atpase subunit 8 (atp8) gene, as all published marine bivalves. In contrast, a complete atp8 gene is present in Hiatella arctica. In addition, we found a putative truncated atp8 gene when re-annotating the mt-genome of Venerupis philippinarum. Both mt-genomes reported here encode all genes on the same strand and have an additional trnM. In Acanthocardia several large non-coding regions are present. One of these contains 3.5 nearly identical copies of a 167 bp motive. In Hiatella, the 3' end of the NADH dehydrogenase subunit (nad)6 gene is duplicated together with the adjacent non-coding region. The gene arrangement of Hiatella is markedly different from all other known molluscan mt-genomes, that of Acanthocardia shows few identities with the Venerupis philippinarum. Phylogenetic analyses on amino acid and nucleotide levels robustly support the Heterodonta and the sister group relationship of Acanthocardia and Venerupis. Monophyletic Bivalvia are resolved only by a Bayesian inference of the nucleotide data set. In all other analyses the two unionid species, being to only ones with genes located on both strands, do not group with the remaining bivalves.

CONCLUSION

The two mt-genomes reported here add to and underline the high variability of gene order and presence of duplications in bivalve and molluscan taxa. Some genomic traits like the loss of the atp8 gene or the encoding of all genes on the same strand are homoplastic among the Bivalvia. These characters, gene order, and the nucleotide sequence data show considerable potential of resolving phylogenetic patterns at lower taxonomic levels.

The relationship between the rate of molecular evolution and the rate of genome rearrangement in animal mitochondrial genomes

Evolution of mitochondrial genes is far from clock-like. The substitution rate varies considerably between species, and there are many species that have a significantly increased rate with respect to their close relatives. There is also considerable variation among species in the rate of gene order rearrangement. Using a set of 55 complete arthropod mitochondrial genomes, we estimate the evolutionary distance from the common ancestor to each species using protein sequences, tRNA sequences, and breakpoint distances (a measure of the degree of genome rearrangement). All these distance measures are correlated. We use relative rate tests to compare pairs of related species in several animal phyla. In the majority of cases, the species with the more highly rearranged genome also has a significantly higher rate of sequence evolution. Species with higher amino acid substitution rates in mitochondria also have more variable amino acid composition in response to mutation pressure. We discuss the possible causes of variation in rates of sequence evolution and gene rearrangement among species and the possible reasons for the observed correlation between the two rates.

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

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

The complete mitochondrial genome of the orbiniid polychaete Orbinia latreillii (Annelida, Orbiniidae)--A novel gene order for Annelida and implications for annelid phylogeny

Relationships of annelid subtaxa are controversially discussed and additional markers are necessarily needed to get further insights into their evolution. Due to their high content of information, mitochondrial genomes have been proven very useful in phylogenetic analyses. Whereas many complete mitochondrial genomes of arthropods are available, lophotrochozoan taxa are only scarcely represented and this is especially true for annelids. Here we present the complete mitochondrial genome of the orbiniid polychaete Orbinia latreillii. The circular genome is 15,558 bp in size and contains the same 37 genes as found in most other metazoans. As in the case for all studied annelids all genes are transcribed from the same strand. Compared with the known data from other annelids at least five gene translocations must be hypothesized for O. latreillii. A comparison of the available data shows that gene translocations within Annelida seem to be less frequent than in molluscs, but more frequent as previously assumed. Phylogenetic analyses of mitochondrial DNA sequence data and amino acid data support an inclusion of Sipuncula within Annelida and a closer relationship to orbiniids is recovered for this taxon.

Complete DNA sequence of the mitochondrial genome of the sea-slug, Aplysia californica: conservation of the gene order in Euthyneura

We have sequenced and characterized the complete mitochondrial genome of the sea slug, Aplysia californica, an important model organism in experimental biology and a representative of Anaspidea (Opisthobranchia, Gastropoda). The mitochondrial genome of Aplysia is in the small end of the observed sizes of animal mitochondrial genomes (14,117 bp, NCBI Accession No. NC_005827). The Aplysia genome, like most other mitochondrial genomes, encodes genes for 2 ribosomal subunit RNAs (small and large rRNAs), 22 tRNAs, and 13 protein subunits (cytochrome c oxidase subunits 1-3, cytochrome b apoenzyme, ATP synthase subunits 6 and 8, and NADH dehydrogenase subunits 1-6 and 4L). The gene order is virtually identical between opisthobranchs and pulmonates, with the majority of differences arising from tRNA translocations. In contrast, the gene order from representatives of basal gastropods and other molluscan classes is significantly different from opisthobranchs and pulmonates. The Aplysia genome was compared to all other published molluscan mitochondrial genomes and phylogenetic analyses were carried out using a concatenated protein alignment. Phylogenetic analyses using maximum likelihood based analyses of the well aligned regions of the protein sequences support both monophyly of Euthyneura (a group including both the pulmonates and opisthobranchs) and Opisthobranchia (as a more derived group). The Aplysia mitochondrial genome sequenced here will serve as an important platform in both comparative and neurobiological studies using this model organism.