Poriferan mtDNA and animal phylogeny based on mitochondrial gene arrangements

Variability and evolution of gene order rearrangement in mitochondrial genomes of arthropods (except Hexapoda)

In the species-rich Phylum Arthropoda, the mitochondrial genome is relatively well conserved both in terms of number and order of genes. However, specific clades have a 'typical' gene order that differs from the putative arthropod ancestral arrangement. The aim of this work was to compare the rate of mitochondrial gene rearrangements at inter- and intra-taxonomic levels in the Arthropoda and to postulate the most parsimonious ancestral orders representing the four major arthropod lineages. For this purpose, we performed a comparative genomic analysis of arthropod mitochondrial genomes available in the NCBI database. Using a combination of bioinformatics methods that examined mitochondrial gene rearrangements in 464 species of arthropods from three subphyla (Chelicerata, Myriapoda, and Crustacea [except Hexapoda, previously analyzed]), we observed differences in the rate of rearrangement within major lineages. A higher rate of mitochondrial genome rearrangement was observed in Crustacea and Chelicerata compared to Myriapoda. Likewise, early branching clades exhibit less variability in mitochondrial genome order than late branching clades, within each subphylum. We identified 'hot regions' in the mitochondrial genome of each studied subphylum, and postulated the most likely ancestral gene order in each subphylum and taxonomic order. Our work provides new evidence on the evolutionary dynamics of mitochondrial genome gene order in arthropods and new mitochondrial genome architectures in different taxonomic divisions within each major lineage of arthropods.

Phylomitogenomics bolsters the high-level classification of Demospongiae (phylum Porifera)

Class Demospongiae is the largest in the phylum Porifera (Sponges) and encompasses nearly 8,000 accepted species in three subclasses: Keratosa, Verongimorpha, and Heteroscleromorpha. Subclass Heteroscleromorpha contains ∼90% of demosponge species and is subdivided into 17 orders. The higher level classification of demosponges underwent major revision as the result of nearly three decades of molecular studies. However, because most of the previous molecular work only utilized partial data from a small number of nuclear and mitochondrial (mt) genes, this classification scheme needs to be tested by larger datasets. Here we compiled a mt dataset for 136 demosponge species-including 64 complete or nearly complete and six partial mt-genome sequences determined or assembled for this study-and used it to test phylogenetic relationships among Demospongiae in general and Heteroscleromorpha in particular. We also investigated the phylogenetic position of Myceliospongia araneosa, a highly unusual demosponge without spicules and spongin fibers, currently classified as Demospongiae incertae sedis, for which molecular data were not available. Our results support the previously inferred sister-group relationship between Heteroscleromorpha and Keratosa + Verongimorpha and suggest five main clades within Heteroscleromorpha: Clade C0 composed of order Haplosclerida; Clade C1 composed of Scopalinida, Sphaerocladina, and Spongillida; Clade C2 composed of Axinellida, Biemnida, Bubarida; Clade C3 composed of Tetractinellida; and Clade C4 composed of Agelasida, Clionaida, Desmacellida, Merliida, Suberitida, Poecilosclerida, Polymastiida, and Tethyida. The inferred relationships among these clades were (C0(C1(C2(C3+C4)))). Analysis of molecular data from M. araneosa placed it in the C3 clade as a sister taxon to the highly skeletonized tetractinellids Microscleroderma sp. and Leiodermatium sp. Molecular clock analysis dated divergences among the major clades in Heteroscleromorpha from the Cambrian to the Early Silurian, the origins of most heteroscleromorph orders in the middle Paleozoic, and the most basal splits within these orders around the Paleozoic to Mesozoic transition. Overall, the results of this study are mostly congruent with the accepted classification of Heteroscleromorpha, but add temporal perspective and new resolution to phylogenetic relationships within this subclass.

Full Mitochondrial Genomes Reveal Species Differences between the Venerid Clams Ruditapes philippinarum and R. variegatus

In natural sea areas along the coast of China, venerid clams Ruditapes philippinarum and R. variegatus exhibit similar adult shell forms and are especially difficult to distinguish as spat and juveniles. This study used comparative mitochondrial genome analysis to reveal differences between these species. The results showed that: (1) the mitochondrial genomes of R. philippinarum and R. variegatus share a large number of similar gene clusters arranged in consistent order, yet they also display noncommon genes, with both gene rearrangements and random losses found; (2) the 13 protein-coding genes in R. philippinarum as well as two-fold and four-fold degenerate sites in R. variegatus have an evident AT bias; (3) the Ka/Ks ratio of the mitochondrial ATP8 gene was significantly higher in R. philippinarum than in R. variegatus, and an analysis of selection pressure revealed that the mitochondrial NADH dehydrogenase subunit 2 gene and NADH dehydrogenase subunit 6 gene of R. variegatus were under great selective pressure during its evolution; and finally, (4) the two species clustered into one branch on a phylogenetic tree, further affirming their phylogenetic closeness. Based on these results, we speculate that the species differences between R. variegatus and R. philippinarum are largely attributable to adaptive evolution to the environment. The present findings provide a reference for the development of germplasm identification.

Comparative mitochondrial genome analysis and phylogenetic relationship among lepidopteran species

Lepidoptera has rich species including many agricultural pests and economical insects around the world. The mitochondrial genomes (mitogenomes) were utilized to explore the phylogenetic relationships between difference taxonomic levels in Lepidoptera. However, the knowledge of mitogenomic characteristics and phylogenetic position about superfamily-level in this order is unresolved. In this study, we integrated 794 mitogenomes consisting of 37 genes and a noncoding control region, which covered 26 lepidopteran superfamilies from newly sequenced and publicly available genomes for comparative genomic and phylogenetic analysis. In primitive taxon, putative start codon of cox1 gene was ATA or ATT instead of CGA, but stop codon of that showed four types, namely TAA, TAG, TA and T. The 7-bp overlap between atp8 and atp6 presented as "ATGATAA". Moreover, the most frequently utilized amino acids were leucine (UUA) in 13 PCGs. Phylogenetic analysis showed that the main backbone relationship in Lepidoptera was (Hepialoidea + (Nepticuloidea + (Adeloidea + (Tischerioidea + (Tineoidea + (Yponomeutoidea + (Gracillarioidea + (Papilionoidea + ((Zygaenoidea + Tortricoidea) + (Gelechioidea + (Pyraloidea + ((Geometroidea + Noctuoidea) + (Lasiocampoidea + Bombycoidea))))))))))))).

Mitogenomes Reveal Alternative Initiation Codons and Lineage-Specific Gene Order Conservation in Echinoderms

The mitochondrial genetic code is much more varied than the standard genetic code. The invertebrate mitochondrial code, for instance, comprises six initiation codons, including five alternative start codons. However, only two initiation codons are known in the echinoderm and flatworm mitochondrial code, the canonical ATG and alternative GTG. Here, we analyzed 23 Asteroidea mitogenomes, including ten newly sequenced species and unambiguously identified at least two other start codons, ATT and ATC, both of which also initiate translation of mitochondrial genes in other invertebrates. These findings underscore the diversity of the genetic code and expand upon the suite of initiation codons among echinoderms to avoid erroneous annotations. Our analyses have also uncovered the remarkable conservation of gene order among asteroids, echinoids, and holothuroids, with only an interchange between two gene positions in asteroids over ∼500 Ma of echinoderm evolution.

The first two complete mitogenomes of the order Apodida from deep-sea chemoautotrophic environments: New insights into the gene rearrangement, origin and evolution of the deep-sea sea cucumbers

The deep-sea ecosystem is considered as the largest and most remote biome of the world. It is meaningful and important to elucidate the life origins by exploring the origin and adaptive genetic mechanisms of the large deep-sea organisms. Sea cucumbers (Holothuroidea) are abundant and economically important group of echinoderms, living from the shallow-waters to deep-sea. In this study, we present the mitochondrial genomes of the sea cucumber Chiridota heheva and Chiridota sp. collected from the deep-sea cold seep and hydrothermal vent, respectively. This is the first reported mitochondrial genomes from the order Apodida. The mitochondrial genomes of C. heheva (17,200 bp) and Chiridota sp. (17,199 bp) display novel gene arrangements with the first protein-coding gene rearrangements in the class Holothuroidea. Bases composition analysis showed that the A + T content of deep-sea holothurians were significantly higher than that of the shallow-water groups. We compared the arrangement of genes from the 24 available holothurian mitogenomes and found that the transposition, reverse transposition and tandem-duplication-random-losses (TDRL) may be involved in the evolution of mitochondrial gene arrangements in Holothuroidea. Phylogenetic analysis revealed that the Apodida clustered with Elasipodida, forming two basal deep-sea holothurian clades. The divergence between the deep-sea and shallow-water holothurians was located at 386.93 Mya, during the Late Devonian. Mitochondrial protein-coding genes of deep-sea holothurians underwent relaxed purifying selection. There are 57 positive selected amino acids sites for some mitochondrial genes of the three deep-sea clades, implying they may involve in the adaption of deep-sea sea cucumbers.

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.

Extensive gene rearrangements in the mitogenomes of congeneric annelid species and insights on the evolutionary history of the genus Ophryotrocha

BACKGROUND

Annelids are one the most speciose and ecologically diverse groups of metazoans. Although a significant effort has been recently invested in sequencing genomes of a wide array of metazoans, many orders and families within the phylum Annelida are still represented by a single specimen of a single species. The genus of interstitial annelids Ophryotrocha (Dorvilleidae, Errantia, Annelida) is among these neglected groups, despite its extensive use as model organism in numerous studies on the evolution of life history, physiological and ecological traits. To compensate for the paucity of genomic information in this genus, we here obtained novel complete mitochondrial genomes of six Ophryotrocha species using next generation sequencing. In addition, we investigated the evolution of the reproductive mode in the Ophryotrocha genus using a phylogeny based on two mitochondrial markers (COXI and 16S rDNA) and one nuclear fragment (Histone H3).

RESULTS

Surprisingly, gene order was not conserved among the six Ophryotrocha species investigated, and varied greatly as compared to those found in other annelid species within the class Errantia. The mitogenome phylogeny for the six Ophryotrocha species displayed a separation of gonochoric and hermaphroditic species. However, this separation was not observed in the phylogeny based on the COX1, 16S rDNA, and H3 genes. Parsimony and Bayesian ancestral trait reconstruction indicated that gonochorism was the most parsimonious ancestral reproductive mode in Ophryotrocha spp.

CONCLUSIONS

Our results highlight the remarkably high level of gene order variation among congeneric species, even in annelids. This encourages the need for additional mitogenome sequencing of annelid taxa in order to properly understand its mtDNA evolution, high biodiversity and phylogenetic relationships.

Phylomitogenomic analyses on collembolan higher taxa with enhanced taxon sampling and discussion on method selection

Collembola are a basal group of Hexapoda renowned for both unique morphological characters and significant ecological roles. However, a robust and plausible phylogenetic relationship between its deeply divergent lineages has yet to be achieved. We carried out a mitophylogenomic study based on a so far the most comprehensive mitochondrial genome dataset. Our data matrix contained mitogenomes of 31 species from almost all major families of all four orders, with 16 mitogenomes newly sequenced and annotated. We compared the linear arrangements of genes along mitochondria across species. Then we conducted 13 analyses each under a different combination of character coding, partitioning scheme and heterotachy models, and assessed their performance in phylogenetic inference. Several hypothetical tree topologies were also tested. Mitogenomic structure comparison revealed that most species share the same gene order of putative ancestral pancrustacean pattern, while seven species from Onychiuridae, Poduridae and Symphypleona bear different levels of gene rearrangements, indicating phylogenetic signals. Tomoceroidea was robustly recovered for the first time in the presence of all its families and subfamilies. Monophyly of Onychiuroidea was supported using unpartitioned models alleviating LBA. Paronellidae was revealed polyphyletic with two subfamilies inserted independently into Entomobryidae. Although Entomobryomorpha has not been well supported, more than half of the analyses obtained convincing topologies by placing Tomoceroidea within or near remaining Entomobryomorpha. The relationship between elongate-shaped and spherical-shaped collembolans still remained ambiguous, but Neelipleona tend to occupy the basal position in most trees. This study showed that mitochondrial genomes could provide important information for reconstructing the relationships among Collembola when suitable analytical approaches are implemented. Of all the data refining and model selecting schemes used in this study, the combination of nucleotide sequences, partitioning model and exclusion of third codon positions performed better in generating more reliable tree topology and higher node supports than others.

The complete mitochondrial genome of Amesia sanguiflua (Lepidoptera, Zygaenidae)

Amesia sanguiflua (Lepidoptera, Zygaenidae) is found in northern India, Myanmar, Indochina, Malay Peninsula, Sumatra, Java, and China. In the present study, we sequenced the complete mitochondrial genome of A. sanguiflua. The mitochondrial genome was 15,203 bp in length, containing a typical set of 37 genes (13 protein-coding genes, 2 rRNA genes, 22 tRNA genes) and a 346bp non-coding A+T-rich region. Phylogenetic analysis using mitochondrial genomes of 40 species showed that A. sanguiflua formed a well-supported monophyletic group with other Zygaenidae species.

Mitochondrial genome of Paruterina candelabraria (Cestoda: Paruterinidae), with implications for the relationships between the genera Cladotaenia and Paruterina

The taxonomic concept for the family Paruterinidae is controversial, especially concerning the position of the genus Cladotaenia, since the latter genus has been placed sometimes in other families, i.e. in the Taeniidae or in the distinct family Cladotaeniidae; finding a solution based on morphological data is difficult and molecular data on paruterinids and related groups are scarce. In this study, the complete mitochondrial (mt) genome sequence of the type-species of the type-genus of the Paruterinidae, Paruterina candelabraria, was determined and annotated. Gene arrangements are identical with those of Cladotaenia vulturi but differing from those of species of the family Taeniidae by the order change between tRNA-SerUCN and tRNA-LeuCUN. Phylogenetic tree was constructed by Bayesian Inference (BI) analysis using the concatenated amino acid sequences of 12 protein-coding genes. The analysis clearly shows that the Paruterinidae and Taeniidae are sister-groups, and Cladotaenia is a sister taxon of Paruterina. This supports the position of the genus Cladotaenia in the family Paruterinidae and reveals the necessity for sequencing additional taxa of the Paruterinidae for better understanding of phylogenetic relationships within the group.

Genome Rearrangement with ILP

The weighted Genome Sorting Problem (wGSP) is to find a minimum-weight sequence of rearrangement operations that transforms a given gene order into another given gene order using rearrangement operations that are associated with a predefined weight. This paper presents a polynomial sized Integer Linear Program -called GeRe-ILP- for solving the wGSP for the following three types of rearrangement operations: inversion , transposition, and inverse transposition. GeRe-ILP uses variables and constraints for gene orders of length . It is studied experimentally on simulated data how different weighting schemes influence the reconstructed scenarios. The influences of the length of the gene orders and of the size of the reconstructed scenarios on the runtime of GeRe-ILP are studied as well.

Phylogeny, evolution and mitochondrial gene order rearrangement in scale worms (Aphroditiformia, Annelida)

Next-generation sequencing (NGS) has become a powerful tool in phylogenetic and evolutionary studies. Here we applied NGS to recover two ribosomal RNA genes (18S and 28S) from 16 species and 15 mitochondrial genomes from 16 species of scale worms representing six families in the suborder Aphroditiformia (Phyllodocida, Annelida), a complex group of polychaetes characterized by the presence of dorsal elytra or scales. The phylogenetic relationship of the several groups of scale worms remains unresolved due to insufficient taxon sampling and low resolution of individual gene markers. Phylogenetic tree topology based on mitochondrial genomes is comparable with that based on concatenated sequences from two mitochondrial genes (cox1 and 16S) and two ribosomal genes (18S and 28S) genes, but has higher statistical support for several clades. Our analyses show that Aphroditiformia is monophyletic, indicating the presence of elytra is an apomorphic trait. Eulepethidae and Aphroditidae together form the sister group to all other families in this suborder, whereas Acoetidae is sister to Iphionidae. Polynoidae is monophyletic, but within this family the deep-sea subfamilies Branchinotogluminae and Macellicephalinae are paraphyletic. Mitochondrial genomes in most scale-worm families have a conserved gene order, but within Polynoidae there are two novel arrangement patterns in the deep-sea clade. Mitochondrial protein-coding genes in polynoids as a whole have evolved under strong purifying selection, but substitution rates in deep-sea species are much higher than those in shallow-water species, indicating that purifying selection is relaxed in deep-sea polynoids. There are positive selected amino acids for some mitochondrial genes of the deep-sea clade, indicating they may involve in the adaption of deep-sea polynoids. Overall, our study (1) provided more evidence for reconstruction of the phylogeny of Aphroditiformia, (2) provided evidence to refute the assumption that mitochondrial gene order in Errantia is conserved, and (3) indicated that the deep-sea extreme environment may have affected the mitochondrial genome evolution rate and gene order arrangement in Polynoidae.

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.

Sequencing, characterization and phylogenomics of the complete mitochondrial genome of Dactylogyrus lamellatus (Monogenea: Dactylogyridae)

Despite the worldwide distribution and pathogenicity of monogenean parasites belonging to the largest helminth genus, Dactylogyrus, there are no complete Dactylogyrinae (subfamily) mitogenomes published to date. In order to fill this knowledge gap, we have sequenced and characterized the complete mitogenome of Dactylogyrus lamellatus, a common parasite on the gills of grass carp (Ctenopharyngodon idella). The circular mitogenome is 15,187 bp in size, containing the standard 22 tRNA genes, 2 rRNA genes, 12 protein-encoding genes and a long non-coding region (NCR). There are two highly repetitive regions in the NCR. We have used concatenated nucleotide sequences of all 36 genes to perform the phylogenetic analysis using Bayesian inference and maximum likelihood approaches. As expected, the two dactylogyrids, D. lamellatus (Dactylogyrinae) and Tetrancistrum nebulosi (Ancyrocephalinae), were closely related to each other. These two formed a sister group with Capsalidae, and this cluster finally formed a further sister group with Gyrodactylidae. Phylogenetic affinity between Dactylogyrinae and Ancyrocephalinae was further confirmed by the similarity in their gene arrangement. The sequencing of the first Dactylogyrinae, along with a more suitable selection of outgroups, has enabled us to infer a much better phylogenetic resolution than recent mitogenomic studies. However, as many lineages of the class Monogenea remain underrepresented or not represented at all, a much larger number of mitogenome sequences will have to be available in order to infer the evolutionary relationships among the monogeneans fully, and with certainty.

Loss and Gain of Group I Introns in the Mitochondrial Cox1 Gene of the Scleractinia (Cnidaria; Anthozoa)

Yaoyang Chuang, Marcelo Kitahara, Hironobu Fukami, Dianne Tracey, David J. Miller, and Chaolun Allen Chen (2017) Group I introns encoding a homing endonuclease gene (HEG) that is potentially capable of sponsoring mobility are present in the cytochrome oxidase subunit 1 (cox1) gene of some Hexacorallia, including a number of scleractinians assigned to the "robust" coral clade. In an e ort to infer the evolutionary history of this cox1 group I intron, DNA sequences were determined for 12 representative "basal" and "complex" corals and for 11 members of the Corallimorpharia, a sister order of the Scleractinia. Comparisons of insertion sites, secondary structures, and amino acid sequences of the HEG implied a common origin for cox1 introns of corallimorpharians, and basal and complex corals, but cox1 introns of robust corals were highly divergent, most likely reflecting independent acquisition. Phylogenetic analyses with a calibrated molecular clock suggested that cox1 introns of scleractinians and corallimorpharians have persisted at the same insertion site as that in the common ancestor 552 million years ago (mya). This ancestral intron was probably lost in complex corals around 213 to 190 mya at the junction between the Trassic and Jurassic. The coral cox1 gene remained intronless until new introns, probably from sponges or fungi, reinvaded different positions of the cox1 gene in robust corals around 135 mya in the Cretaceous, and then it subsequently began to lose them around 65.5 mya in some robust coral lineages coincident with the later Maastrichtian extinction at the Cretaceous-Tertiary boundary.

The complete mitochondrial genome of the deep-sea sponge Poecillastra laminaris (Astrophorida, Vulcanellidae)

The complete mitochondrial genome of a New Zealand specimen of the deep-sea sponge Poecillastra laminaris (Sollas, 1886) (Astrophorida, Vulcanellidae), from the Colville Ridge, New Zealand, was sequenced using the 454 Life Science pyrosequencing system. To identify homologous mitochondrial sequences, the 454 reads were mapped to the complete mitochondrial genome sequence of Geodia neptuni (GeneBank No. NC_006990). The P. laminaris genome is 18,413 bp in length and includes 14 protein-coding genes, 24 transfer RNA genes and 2 ribosomal RNA genes. Gene order resembled that of other demosponges. The base composition of the genome is A (29.1%), T (35.2%), C (14.0%) and G (21.7%). This is the second published mitogenome for a sponge of the order Astrophorida and will be useful in future phylogenetic analysis of deep-sea sponges.

The mitochondrial genome of Paragyrodactylus variegatus (Platyhelminthes: Monogenea): differences in major non-coding region and gene order compared to Gyrodactylus

BACKGROUND

Paragyrodactylus Gvosdev and Martechov, 1953, a viviparous genus of ectoparasite within the Gyrodactylidae, contains three nominal species all of which infect Asian river loaches. The group is suspected to be a basal lineage within Gyrodactylus Nordmann, 1832 sensu lato although this remains unclear. Further molecular study, beyond characterization of the standard Internal Transcribed Spacer region, is needed to clarify the evolutionary relationships within the family and the placement of this genus.

METHODS

The mitochondrial genome of Paragyrodactylus variegatus You, King, Ye and Cone, 2014 was amplified in six parts from a single worm, sequenced using primer walking, annotated and analyzed using bioinformatic tools.

RESULTS

The mitochondrial genome of P. variegatus is 14,517 bp, containing 12 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a major non-coding region (NCR). The overall A + T content of the mitochondrial genome is 76.3%, which is higher than all reported mitochondrial genomes of monogeneans. All of the 22 tRNAs have the typical cloverleaf secondary structure, except tRNACys, tRNASer1 and tRNASer2 that lack the dihydrouridine (DHU) arm. There are six domains (domain III is absent) and three domains in the inferred secondary structures of the large ribosomal subunit (rrnL) and small ribosomal subunit (rrnS), respectively. The NCR includes six 40 bp tandem repeat units and has the double identical poly-T stretches, stem-loop structure and some surrounding structure elements. The gene order (tRNAGln, tRNAMet and NCR) differs in arrangement compared to the mitochondrial genomes reported from Gyrodactylus spp.

CONCLUSION

The Duplication and Random Loss Model and Recombination Model together are the most plausible explanations for the variation in gene order. Both morphological characters and characteristics of the mitochondrial genome support Paragyrodactylus as a distinct genus from Gyrodactylus. Considering their specific distribution and known hosts, we believe that Paragyrodactylus is a relict freshwater lineage of viviparous monogenean isolated in the high plateaus of central Asia on closely related river loaches.

Phylogenetic and biological significance of evolutionary elements from metazoan mitochondrial genomes

The evolutionary history of living species is usually inferred through the phylogenetic analysis of molecular and morphological information using various mathematical models. New challenges in phylogenetic analysis are centered mostly on the search for accurate and efficient methods to handle the huge amounts of sequence data generated from newer genome sequencing. The next major challenge is the determination of relationships between the evolution of structural elements and their functional implementation, which is largely ignored in previous analyses. Here, we described the discovery of structural elements in metazoan mitochondrial genomes, termed key K-strings, that can serve as a basis for phylogenetic tree construction. Although comprising only a small fraction (0.73%) of all K-strings, these key K-strings are pivotal to the tree construction because they allow for a significant reduction in the computational time required to construct phylogenetic trees, and more importantly, they make significant improvement to the results of phylogenetic inference. The trees constructed from the key K-strings were consistent overall to our current view of metazoan phylogeny and exhibited a more rational topology than the trees constructed by using other conventional methods. Surprisingly, the key K-strings tended to accumulate in the conserved regions of the original sequences, which were most likely due to strong selection pressure. Furthermore, the special structural features of the key K-strings should have some potential applications in the study of the structures and functions relationship of proteins and in the determination of evolutionary trajectory of species. The novelty and potential importance of key K-strings lead us to believe that they are essential evolutionary elements. As such, they may play important roles in the process of species evolution and their physical existence. Further studies could lead to discoveries regarding the relationship between evolution and processes of speciation.

The complete maternally and paternally inherited mitochondrial genomes of the endangered freshwater mussel Solenaia carinatus (Bivalvia: Unionidae) and implications for Unionidae taxonomy

Doubly uniparental inheritance (DUI) is an exception to the typical maternal inheritance of mitochondrial (mt) DNA in Metazoa, and found only in some bivalves. In species with DUI, there are two highly divergent gender-associated mt genomes: maternal (F) and paternal (M), which transmit independently and show different tissue localization. Solenaia carinatus is an endangered freshwater mussel species exclusive to Poyang Lake basin, China. Anthropogenic events in the watershed greatly threaten the survival of this species. Nevertheless, the taxonomy of S. carinatus based on shell morphology is confusing, and the subfamilial placement of the genus Solenaia remains unclear. In order to clarify the taxonomic status and discuss the phylogenetic implications of family Unionidae, the entire F and M mt genomes of S. carinatus were sequenced and compared with the mt genomes of diverse freshwater mussel species. The complete F and M mt genomes of S. carinatus are 16716 bp and 17102 bp in size, respectively. The F and M mt genomes of S. carinatus diverge by about 40% in nucleotide sequence and 48% in amino acid sequence. Compared to F counterparts, the M genome shows a more compact structure. Different gene arrangements are found in these two gender-associated mt genomes. Among these, the F genome cox2-rrnS gene order is considered to be a genome-level synapomorphy for female lineage of the subfamily Gonideinae. From maternal and paternal mtDNA perspectives, the phylogenetic analyses of Unionoida indicate that S. carinatus belongs to Gonideinae. The F and M clades in freshwater mussels are reciprocal monophyly. The phylogenetic trees advocate the classification of sampled Unionidae species into four subfamilies: Gonideinae, Ambleminae, Anodontinae, and Unioninae, which is supported by the morphological characteristics of glochidia.

Reconstruction of family-level phylogenetic relationships within Demospongiae (Porifera) using nuclear encoded housekeeping genes

BACKGROUND

Demosponges are challenging for phylogenetic systematics because of their plastic and relatively simple morphologies and many deep divergences between major clades. To improve understanding of the phylogenetic relationships within Demospongiae, we sequenced and analyzed seven nuclear housekeeping genes involved in a variety of cellular functions from a diverse group of sponges.

METHODOLOGY/PRINCIPAL FINDINGS

We generated data from each of the four sponge classes (i.e., Calcarea, Demospongiae, Hexactinellida, and Homoscleromorpha), but focused on family-level relationships within demosponges. With data for 21 newly sampled families, our Maximum Likelihood and Bayesian-based approaches recovered previously phylogenetically defined taxa: Keratosa(p), Myxospongiae(p), Spongillida(p), Haploscleromorpha(p) (the marine haplosclerids) and Democlavia(p). We found conflicting results concerning the relationships of Keratosa(p) and Myxospongiae(p) to the remaining demosponges, but our results strongly supported a clade of Haploscleromorpha(p)+Spongillida(p)+Democlavia(p). In contrast to hypotheses based on mitochondrial genome and ribosomal data, nuclear housekeeping gene data suggested that freshwater sponges (Spongillida(p)) are sister to Haploscleromorpha(p) rather than part of Democlavia(p). Within Keratosa(p), we found equivocal results as to the monophyly of Dictyoceratida. Within Myxospongiae(p), Chondrosida and Verongida were monophyletic. A well-supported clade within Democlavia(p), Tetractinellida(p), composed of all sampled members of Astrophorina and Spirophorina (including the only lithistid in our analysis), was consistently revealed as the sister group to all other members of Democlavia(p). Within Tetractinellida(p), we did not recover monophyletic Astrophorina or Spirophorina. Our results also reaffirmed the monophyly of order Poecilosclerida (excluding Desmacellidae and Raspailiidae), and polyphyly of Hadromerida and Halichondrida.

CONCLUSIONS/SIGNIFICANCE

These results, using an independent nuclear gene set, confirmed many hypotheses based on ribosomal and/or mitochondrial genes, and they also identified clades with low statistical support or clades that conflicted with traditional morphological classification. Our results will serve as a basis for future exploration of these outstanding questions using more taxon- and gene-rich datasets.

Comparative analysis of complete mitochondrial genome sequences confirms independent origins of plant-parasitic nematodes

Background

The nematode infraorder Tylenchomorpha (Class Chromadorea) includes plant parasites that are of agricultural and economic importance, as well as insect-associates and fungal feeding species. Among tylenchomorph plant parasites, members of the superfamily Tylenchoidea, such as root-knot nematodes, have great impact on agriculture. Of the five superfamilies within Tylenchomorpha, one (Aphelenchoidea) includes mainly fungal-feeding species, but also some damaging plant pathogens, including certain Bursaphelenchus spp. The evolutionary relationships of tylenchoid and aphelenchoid nematodes have been disputed based on classical morphological features and molecular data. For example, similarities in the structure of the stomatostylet suggested a common evolutionary origin. In contrast, phylogenetic hypotheses based on nuclear SSU ribosomal DNA sequences have revealed paraphyly of Aphelenchoidea, with, for example, fungal-feeding Aphelenchus spp. within Tylenchomorpha, but Bursaphelenchus and Aphelenchoides spp. more closely related to infraorder Panagrolaimomorpha. We investigated phylogenetic relationships of plant-parasitic tylenchoid and aphelenchoid species in the context of other chromadorean nematodes based on comparative analysis of complete mitochondrial genome data, including two newly sequenced genomes from Bursaphelenchus xylophilus (Aphelenchoidea) and Pratylenchus vulnus (Tylenchoidea).

Results

The complete mitochondrial genomes of B. xylophilus and P. vulnus are 14,778 bp and 21,656 bp, respectively, and identical to all other chromadorean nematode mtDNAs in that they contain 36 genes (lacking atp8) encoded in the same direction. Their mitochondrial protein-coding genes are biased toward use of amino acids encoded by T-rich codons, resulting in high A+T richness. Phylogenetic analyses of both nucleotide and amino acid sequence datasets using maximum likelihood and Bayesian methods did not support B. xylophilus as most closely related to Tylenchomorpha (Tylenchoidea). Instead, B. xylophilus, was nested within a strongly supported clade consisting of species from infraorders Rhabditomorpha, Panagrolaimomorpha, Diplogasteromorpha, and Ascaridomorpha. The clade containing sampled Tylenchoidea (P. vulnus, H. glycines, and R. similis) was sister to all analyzed chromadoreans. Comparison of gene arrangement data was also consistent with the phylogenetic relationships as inferred from sequence data. Alternative tree topologies depicting a monophyletic grouping of B. xylophilus (Aphelenchoidea) plus Tylenchoidea, Tylenchoidea plus Diplogasteromorpha (Pristionchus pacificus), or B. xylophilus plus Diplogasteromorpha were significantly worse interpretations of the mtDNA data.

Conclusions

Phylogenetic trees inferred from nucleotide and amino acid sequences of mtDNA coding genes are in agreement that B. xylophilus (the single representative of Aphelenchoidea) is not closely related to Tylenchoidea, indicating that these two groups of plant parasites do not share an exclusive most recent common ancestor, and that certain morphological similarities between these stylet-bearing nematodes must result from convergent evolution. In addition, the exceptionally large mtDNA genome size of P. vulnus, which is the largest among chromadorean nematode mtDNAs sequenced to date, results from lengthy repeated segments in non-coding regions.

Complete mitochondrial genome of the free-living earwig, Challia fletcheri (Dermaptera: Pygidicranidae) and phylogeny of Polyneoptera

The insect order Dermaptera, belonging to Polyneoptera, includes ∼2,000 extant species, but no dermapteran mitochondrial genome has been sequenced. We sequenced the complete mitochondrial genome of the free-living earwig, Challia fletcheri, compared its genomic features to other available mitochondrial sequences from polyneopterous insects. In addition, the Dermaptera, together with the other known polyneopteran mitochondrial genome sequences (protein coding, ribosomal RNA, and transfer RNA genes), were employed to understand the phylogeny of Polyneoptera, one of the least resolved insect phylogenies, with emphasis on the placement of Dermaptera. The complete mitochondrial genome of C. fletcheri presents the following several unusual features: the longest size in insects is 20,456 bp; it harbors the largest tandem repeat units (TRU) among insects; it displays T- and G-skewness on the major strand and A- and C-skewness on the minor strand, which is a reversal of the general pattern found in most insect mitochondrial genomes, and it possesses a unique gene arrangement characterized by a series of gene translocations and/or inversions. The reversal pattern of skewness is explained in terms of inversion of replication origin. All phylogenetic analyses consistently placed Dermaptera as the sister to Plecoptera, leaving them as the most basal lineage of Polyneoptera or sister to Ephemeroptera, and placed Odonata consistently as the most basal lineage of the Pterygota.

A comparative study of nemertean complete mitochondrial genomes, including two new ones for Nectonemertes cf. mirabilis and Zygeupolia rubens, may elucidate the fundamental pattern for the phylum Nemertea

BACKGROUND

The mitochondrial genome is important for studying genome evolution as well as reconstructing the phylogeny of organisms. Complete mitochondrial genome sequences have been reported for more than 2200 metazoans, mainly vertebrates and arthropods. To date, from a total of about 1275 described nemertean species, only three complete and two partial mitochondrial DNA sequences from nemerteans have been published. Here, we report the entire mitochondrial genomes for two more nemertean species: Nectonemertes cf. mirabilis and Zygeupolia rubens.

RESULTS

The sizes of the entire mitochondrial genomes are 15365 bp for N. cf. mirabilis and 15513 bp for Z. rubens. Each circular genome contains 37 genes and an AT-rich non-coding region, and overall nucleotide composition is AT-rich. In both species, there is significant strand asymmetry in the distribution of nucleotides, with the coding strand being richer in T than A and in G than C. The AT-rich non-coding regions of the two genomes have some repeat sequences and stem-loop structures, both of which may be associated with the initiation of replication or transcription. The 22 tRNAs show variable substitution patterns in nemerteans, with higher sequence conservation in genes located on the H strand. Gene arrangement of N. cf. mirabilis is identical to that of Paranemertes cf. peregrina, both of which are Hoplonemertea, while that of Z. rubens is the same as in Lineus viridis, both of which are Heteronemertea. Comparison of the gene arrangements and phylogenomic analysis based on concatenated nucleotide sequences of the 12 mitochondrial protein-coding genes revealed that species with closer relationships share more identical gene blocks.

CONCLUSION

The two new mitochondrial genomes share many features, including gene contents, with other known nemertean mitochondrial genomes. The tRNA families display a composite substitution pathway. Gene order comparison to the proposed ground pattern of Bilateria and some lophotrochozoans suggests that the nemertean ancestral mitochondrial gene order most closely resembles the heteronemertean type. Phylogenetic analysis proposes a sister-group relationship between Hetero- and Hoplonemertea, which supports one of two recent alternative hypotheses of nemertean phylogeny.

The complete mitochondrial genome sequence of Oncicola luehei (Acanthocephala: Archiacanthocephala) and its phylogenetic position within Syndermata

In the present study, we determined the complete mitochondrial genome sequence of Oncicola luehei (14,281bp), the first archiacanthocephalan representative and the second complete sequence from the phylum Acanthocephala. The complete genome contains 36 genes including 12 protein coding genes, 22 transfer RNA (tRNA) genes and 2 ribosomal RNA genes (rrnL and rrnS) as reported for other syndermatan species. All genes are encoded on the same strand. The overall nucleotide composition of O. luehei mtDNA is 37.7% T, 29.6% G, 22.5% A, and 10.2% C. The overall A+T content (60.2%) is much lower, compared to other syndermatan species reported so far, due to the high frequency (18.3%) of valine encoded by GTN in its protein-coding genes. Results from phylogenetic analyses of amino acid sequences for 10 protein-coding genes from 41 representatives of major metazoan groups including O. luehei supported monophyly of the phylum Acanthocephala and of the clade Syndermata (Acanthocephala+Rotifera), and the paraphyly of the clade Eurotatoria (classes Bdelloidea+Monogononta from phylum Rotifera). Considering the position of the acanthocephalan species within Syndermata, it is inferred that obligatory parasitism characteristic of acanthocephalans was acquired after the common ancestor of acanthocephalans diverged from its sister group, Bdelloidea. Additional comparison of complete mtDNA sequences from unsampled acanthocephalan lineages, especially classes Polyacanthocephala and Eoacanthocephala, is required to test if mtDNA provides reliable information for the evolutionary relationships and pattern of life history diversification found in the syndermatan groups.

The complete mitochondrial genome sequence of Euphausia pacifica (Malacostraca: Euphausiacea) reveals a novel gene order and unusual tandem repeats

Euphausiid krill are dominant organisms in the zooplankton population and play a central role in marine ecosystems. Euphausia pacifica (Malacostraca: Euphausiacea) is one of the most important and dominant crustaceans in the North Pacific Ocean. In this paper, we described the gene content, organization, and codon usage of the E. pacifica mitochondrial genome. The mitochondrial genome of E. pacifica is 16 898 bp in length and contains a standard set of 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Translocation of three transfer RNAs (trnL1, trnL2, and trnW) was found in the E. pacifica mitochondrial genome when comparing with the pancrustacean ground pattern. The rate of Ka/Ks in 13 protein-coding genes among three krill is much less than 1, which indicates a strong purifying selection within this group. The largest noncoding region in the E. pacifica mitochondrial genome contains one section with tandem repeats (4.7 × 154 bp), which are the largest tandem repeats found in malacostracan mitochondrial genomes so far. All analyses based on nucleotide and amino acid data strongly support the monophyly of Stomatopoda, Penaeidae, Caridea, Brachyura, and Euphausiacea. The Bayesian analysis of nucleotide and amino acid datasets strongly supports the close relationship between Euphausiacea and Decapoda, which confirms traditional findings. The maximum likelihood analysis based on amino acid data strongly supports the close relationship between Euphausiacea and Penaeidae, which destroys the monophyly of Decapoda.

Monophyly of clade III nematodes is not supported by phylogenetic analysis of complete mitochondrial genome sequences

Background

The orders Ascaridida, Oxyurida, and Spirurida represent major components of zooparasitic nematode diversity, including many species of veterinary and medical importance. Phylum-wide nematode phylogenetic hypotheses have mainly been based on nuclear rDNA sequences, but more recently complete mitochondrial (mtDNA) gene sequences have provided another source of molecular information to evaluate relationships. Although there is much agreement between nuclear rDNA and mtDNA phylogenies, relationships among certain major clades are different. In this study we report that mtDNA sequences do not support the monophyly of Ascaridida, Oxyurida and Spirurida (clade III) in contrast to results for nuclear rDNA. Results from mtDNA genomes show promise as an additional independently evolving genome for developing phylogenetic hypotheses for nematodes, although substantially increased taxon sampling is needed for enhanced comparative value with nuclear rDNA. Ultimately, topological incongruence (and congruence) between nuclear rDNA and mtDNA phylogenetic hypotheses will need to be tested relative to additional independent loci that provide appropriate levels of resolution.

Results

For this comparative phylogenetic study, we determined the complete mitochondrial genome sequences of three nematode species, Cucullanus robustus (13,972 bp) representing Ascaridida, Wellcomia siamensis (14,128 bp) representing Oxyurida, and Heliconema longissimum (13,610 bp) representing Spirurida. These new sequences were used along with 33 published nematode mitochondrial genomes to investigate phylogenetic relationships among chromadorean orders. Phylogenetic analyses of both nucleotide and amino acid sequence datasets support the hypothesis that Ascaridida is nested within Rhabditida. The position of Oxyurida within Chromadorea varies among analyses; in most analyses this order is sister to the Ascaridida plus Rhabditida clade, with representative Spirurida forming a distinct clade, however, in one case Oxyurida is sister to Spirurida. Ascaridida, Oxyurida, and Spirurida (the sampled clade III taxa) do not form a monophyletic group based on complete mitochondrial DNA sequences. Tree topology tests revealed that constraining clade III taxa to be monophyletic, given the mtDNA datasets analyzed, was a significantly worse result.

Conclusion

The phylogenetic hypotheses from comparative analysis of the complete mitochondrial genome data (analysis of nucleotide and amino acid datasets, and nucleotide data excluding 3^rd positions) indicates that nematodes representing Ascaridida, Oxyurida and Spirurida do not share an exclusive most recent common ancestor, in contrast to published results based on nuclear ribosomal DNA. Overall, mtDNA genome data provides reliable support for nematode relationships that often corroborates findings based on nuclear rDNA. It is anticipated that additional taxonomic sampling will provide a wealth of information on mitochondrial genome evolution and sequence data for developing phylogenetic hypotheses for the phylum Nematoda.

The enigmatic mitochondrial genome of Rhabdopleura compacta (Pterobranchia) reveals insights into selection of an efficient tRNA system and supports monophyly of Ambulacraria

Background

The Hemichordata comprises solitary-living Enteropneusta and colonial-living Pterobranchia, sharing morphological features with both Chordata and Echinodermata. Despite their key role for understanding deuterostome evolution, hemichordate phylogeny is controversial and only few molecular data are available for phylogenetic analysis. Furthermore, mitochondrial sequences are completely lacking for pterobranchs. Therefore, we determined and analyzed the complete mitochondrial genome of the pterobranch Rhabdopleura compacta to elucidate deuterostome evolution. Thereby, we also gained important insights in mitochondrial tRNA evolution.

Results

The mitochondrial DNA of Rhabdopleura compacta corresponds in size and gene content to typical mitochondrial genomes of metazoans, but shows the strongest known strand-specific mutational bias in the nucleotide composition among deuterostomes with a very GT-rich main-coding strand. The order of the protein-coding genes in R. compacta is similar to that of the deuterostome ground pattern. However, the protein-coding genes have been highly affected by a strand-specific mutational pressure showing unusual codon frequency and amino acid composition. This composition caused extremely long branches in phylogenetic analyses. The unusual codon frequency points to a selection pressure on the tRNA translation system to codon-anticodon sequences of highest versatility instead of showing adaptations in anticodon sequences to the most frequent codons. Furthermore, an assignment of the codon AGG to Lysine has been detected in the mitochondrial genome of R. compacta, which is otherwise observed only in the mitogenomes of some arthropods. The genomes of these arthropods do not have such a strong strand-specific bias as found in R. compacta but possess an identical mutation in the anticodon sequence of the tRNA_Lys.

Conclusion

A strong reversed asymmetrical mutational constraint in the mitochondrial genome of Rhabdopleura compacta may have arisen by an inversion of the replication direction and adaptation to this bias in the protein sequences leading to an enigmatic mitochondrial genome. Although, phylogenetic analyses of protein coding sequences are hampered, features of the tRNA system of R. compacta support the monophyly of Ambulacraria. The identical reassignment of AGG to Lysine in two distinct groups may have occurred by convergent evolution in the anticodon sequence of the tRNA_Lys.

The phylogenetic position of Acoela as revealed by the complete mitochondrial genome of Symsagittifera roscoffensis

Background

Acoels are simply organized unsegmented worms, lacking hindgut and anus. Several publications over recent years challenge the long-held view that acoels are early offshoots of the flatworms. Instead a basal position as sister group to all other bilaterian animals was suggested, mainly based on molecular evidence. This led to the view that features of acoels might reflect those of the last common ancestor of Bilateria, and resulted in several evo-devo studies trying to interpret bilaterian evolution using acoels as a proxy model for the "Urbilateria".

Results

We describe the first complete mitochondrial genome sequence of a member of the Acoela, Symsagittifera roscoffensis. Gene content and circular organization of the mitochondrial genome does not significantly differ from other bilaterian animals. However, gene order shows no similarity to any other mitochondrial genome within the Metazoa. Phylogenetic analyses of concatenated alignments of amino acid sequences from protein coding genes support a position of Acoela and Nemertodermatida as the sister group to all other Bilateria. Our data provided no support for a sister group relationship between Xenoturbellida and Acoela or Acoelomorpha. The phylogenetic position of Xenoturbella bocki as sister group to or part of the deuterostomes was also unstable.

Conclusions

Our phylogenetic analysis supports the view that acoels and nemertodermatids are the earliest divergent extant lineage of Bilateria. As such they remain a valid source for seeking primitive characters present in the last common ancestor of Bilateria. Gene order of mitochondrial genomes seems to be very variable among Acoela and Nemertodermatida and the groundplan for the metazoan mitochondrial genome remains elusive. More data are needed to interpret mitochondrial genome evolution at the base of Bilateria.

Sessile snails, dynamic genomes: gene rearrangements within the mitochondrial genome of a family of caenogastropod molluscs

BACKGROUND

Widespread sampling of vertebrates, which comprise the majority of published animal mitochondrial genomes, has led to the view that mitochondrial gene rearrangements are relatively rare, and that gene orders are typically stable across major taxonomic groups. In contrast, more limited sampling within the Phylum Mollusca has revealed an unusually high number of gene order arrangements. Here we provide evidence that the lability of the molluscan mitochondrial genome extends to the family level by describing extensive gene order changes that have occurred within the Vermetidae, a family of sessile marine gastropods that radiated from a basal caenogastropod stock during the Cenozoic Era.

RESULTS

Major mitochondrial gene rearrangements have occurred within this family at a scale unexpected for such an evolutionarily young group and unprecedented for any caenogastropod examined to date. We determined the complete mitochondrial genomes of four species (Dendropoma maximum, D. gregarium, Eualetes tulipa, and Thylacodes squamigerus) and the partial mitochondrial genomes of two others (Vermetus erectus and Thylaeodus sp.). Each of the six vermetid gastropods assayed possessed a unique gene order. In addition to the typical mitochondrial genome complement of 37 genes, additional tRNA genes were evident in D. gregarium (trnK) and Thylacodes squamigerus (trnV, trnLUUR). Three pseudogenes and additional tRNAs found within the genome of Thylacodes squamigerus provide evidence of a past duplication event in this taxon. Likewise, high sequence similarities between isoaccepting leucine tRNAs in Thylacodes, Eualetes, and Thylaeodus suggest that tRNA remolding has been rife within this family. While vermetids exhibit gene arrangements diagnostic of this family, they also share arrangements with littorinimorph caenogastropods, with which they have been linked based on sperm morphology and primary sequence-based phylogenies.

CONCLUSIONS

We have uncovered major changes in gene order within a family of caenogastropod molluscs that are indicative of a highly dynamic mitochondrial genome. Studies of mitochondrial genomes at such low taxonomic levels should help to illuminate the dynamics of gene order change, since the telltale vestiges of gene duplication, translocation, and remolding have not yet been erased entirely. Likewise, gene order characters may improve phylogenetic hypotheses at finer taxonomic levels than once anticipated and aid in investigating the conditions under which sequence-based phylogenies lack resolution or prove misleading.

Ecdysozoan mitogenomics: evidence for a common origin of the legged invertebrates, the Panarthropoda

Ecdysozoa is the recently recognized clade of molting animals that comprises the vast majority of extant animal species and the most important invertebrate model organisms--the fruit fly and the nematode worm. Evolutionary relationships within the ecdysozoans remain, however, unresolved, impairing the correct interpretation of comparative genomic studies. In particular, the affinities of the three Panarthropoda phyla (Arthropoda, Onychophora, and Tardigrada) and the position of Myriapoda within Arthropoda (Mandibulata vs. Myriochelata hypothesis) are among the most contentious issues in animal phylogenetics. To elucidate these relationships, we have determined and analyzed complete or nearly complete mitochondrial genome sequences of two Tardigrada, Hypsibius dujardini and Thulinia sp. (the first genomes to date for this phylum); one Priapulida, Halicryptus spinulosus; and two Onychophora, Peripatoides sp. and Epiperipatus biolleyi; and a partial mitochondrial genome sequence of the Onychophora Euperipatoides kanagrensis. Tardigrada mitochondrial genomes resemble those of the arthropods in term of the gene order and strand asymmetry, whereas Onychophora genomes are characterized by numerous gene order rearrangements and strand asymmetry variations. In addition, Onychophora genomes are extremely enriched in A and T nucleotides, whereas Priapulida and Tardigrada are more balanced. Phylogenetic analyses based on concatenated amino acid coding sequences support a monophyletic origin of the Ecdysozoa and the position of Priapulida as the sister group of a monophyletic Panarthropoda (Tardigrada plus Onychophora plus Arthropoda). The position of Tardigrada is more problematic, most likely because of long branch attraction (LBA). However, experiments designed to reduce LBA suggest that the most likely placement of Tardigrada is as a sister group of Onychophora. The same analyses also recover monophyly of traditionally recognized arthropod lineages such as Arachnida and of the highly debated clade Mandibulata.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

The mitochondrial genome structure of Xenoturbella bocki (phylum Xenoturbellida) is ancestral within the deuterostomes

Background

Mitochondrial genome comparisons contribute in multiple ways when inferring animal relationships. As well as primary sequence data, rare genomic changes such as gene order, shared gene boundaries and genetic code changes, which are unlikely to have arisen through convergent evolution, are useful tools in resolving deep phylogenies. Xenoturbella bocki is a morphologically simple benthic marine worm recently found to belong among the deuterostomes. Here we present analyses comparing the Xenoturbella bocki mitochondrial gene order, genetic code and control region to those of other metazoan groups.

Results

The complete mitochondrial genome sequence of Xenoturbella bocki was determined. The gene order is most similar to that of the chordates and the hemichordates, indicating that this conserved mitochondrial gene order might be ancestral to the deuterostome clade. Using data from all phyla of deuterostomes, we infer the ancestral mitochondrial gene order for this clade. Using inversion and breakpoint analyses of metazoan mitochondrial genomes, we test conflicting hypotheses for the phylogenetic placement of Xenoturbella and find a closer affinity to the hemichordates than to other metazoan groups. Comparative analyses of the control region reveal similarities in the transcription initiation and termination sites and origin of replication of Xenoturbella with those of the vertebrates. Phylogenetic analyses of the mitochondrial sequence indicate a weakly supported placement as a basal deuterostome, a result that may be the effect of compositional bias.

Conclusion

The mitochondrial genome of Xenoturbella bocki has a very conserved gene arrangement in the deuterostome group, strikingly similar to that of the hemichordates and the chordates, and thus to the ancestral deuterostome gene order. Similarity to the hemichordates in particular is suggested by inversion and breakpoint analysis. Finally, while phylogenetic analyses of the mitochondrial sequences support a basal deuterostome placement, support for this decreases with the use of more sophisticated models of sequence evolution.

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.