Mitochondrial Genome Evolution in Annelida-A Systematic Study on Conservative and Variable Gene Orders and the Factors Influencing its Evolution.

The mitochondrial genomes of Bilateria are relatively conserved in their protein-coding, rRNA, and tRNA gene complement, but the order of these genes can range from very conserved to very variable depending on the taxon. The supposedly conserved gene order of Annelida has been used to support the placement of some taxa within Annelida. Recently, authors have cast doubts on the conserved nature of the annelid gene order. Various factors may influence gene order variability including, among others, increased substitution rates, base composition differences, structure of noncoding regions, parasitism, living in extreme habitats, short generation times, and biomineralization. However, these analyses were neither done systematically nor based on well-established reference trees. Several focused on only a few of these factors and biological factors were usually explored ad-hoc without rigorous testing or correlation analyses. Herein, we investigated the variability and evolution of the annelid gene order and the factors that potentially influenced its evolution, using a comprehensive and systematic approach. The analyses were based on 170 genomes, including 33 previously unrepresented species. Our analyses included 706 different molecular properties, 20 life-history and ecological traits, and a reference tree corresponding to recent improvements concerning the annelid tree. The results showed that the gene order with and without tRNAs is generally conserved. However, individual taxa exhibit higher degrees of variability. None of the analyzed life-history and ecological traits explained the observed variability across mitochondrial gene orders. In contrast, the combination and interaction of the best-predicting factors for substitution rate and base composition explained up to 30% of the observed variability. Accordingly, correlation analyses of different molecular properties of the mitochondrial genomes showed an intricate network of direct and indirect correlations between the different molecular factors. Hence, gene order evolution seems to be driven by molecular evolutionary aspects rather than by life history or ecology. On the other hand, variability of the gene order does not predict if a taxon is difficult to place in molecular phylogenetic reconstructions using sequence data or not. We also discuss the molecular properties of annelid mitochondrial genomes considering canonical views on gene evolution and potential reasons why the canonical views do not always fit to the observed patterns without making some adjustments. [Annelida; compositional biases; ecology; gene order; life history; macroevolution; mitochondrial genomes; substitution rates.].

Syllidae mitochondrial gene order is unusually variable for Annelida.

Complete mitochondrial genomes of five syllids (Streptosyllis sp., Eusyllis blomstrandi, Myrianida brachycephala, Typosyllis antoni and Typosyllis sp.) have been obtained using Illumina sequencing. Together with two previous studied taxa (Ramisyllis multicaudata and Trypanobia cryptica), the analysed sequences represent most of the main lineages within the family Syllidae (Anoplosyllinae, Eusyllinae, Autolytinae and Syllinae). The genomic features, gene order and phylogenetic relationships are examined. Unusual for annelids, syllid mitochondrial genomes are highly variable in their gene order. Considering genomic features, such as length, skewness, gene content, and codon bias, most similar to the rest of annelids are the genomes of E. blomstrandi and M. brachycephala, while Streptosyllis sp. and the analysed sylline taxa (R. multicaudata, T. cryptica, T. antoni and Typosyllis sp.) are the most dissimilar. Two methionine tRNA's (trnM) have been found in T. antoni and Typosyllis sp. The mt genomes of these latter taxa are the longest with numerous non-coding regions. The 13 protein coding genes, as well as the rRNA's are used to perform phylogenetic analyses that recovered the relationships within the family explored before by previous authors. The gene order in Syllidae shows very different patterns. E. blomstrandi and M. prolifera show a similar pattern to the one found in Pleistoannelida; however this might have changed at least twice within Syllidae: in Streptosyllis sp. and within Syllinae. All analysed Syllinae show different gene orders, thereby illustrating more variability as all other pleistoannelids analysed so far. The information provided herein allows a more accurate reconstruction of the possible evolutionary scenarios in Syllidae.

Evolution of mitochondrial gene order in Annelida.

Annelida is a highly diverse animal group with over 21,000 described species. As part of Lophotrochozoa, the vast majority of annelids are currently classified into two groups: Errantia and Sedentaria, together forming Pleistoannelida. Besides these taxa, Sipuncula, Amphinomidae, Chaetopteridae, Oweniidae and Magelonidae can be found branching at the base of the tree. Comparisons of mitochondrial genomes have been used to investigate phylogenetic relationship within animal taxa. Complete annelid mitochondrial genomes are available for some Sedentaria and Errantia and in most cases exhibit a highly conserved gene order. Only two complete genomes have been published from the basal branching lineages and these are restricted to Sipuncula. We describe the first complete mitochondrial genome sequences for all other basal branching annelid families: Owenia fusiformis (Oweniidae), Magelona mirabilis (Magelonidae), Eurythoe complanata (Amphinomidae), Chaetopterus variopedatus and Phyllochaetopterus sp. (Chaetopteridae). The mitochondrial gene order of all these taxa is substantially different from the pattern found in Pleistoannelida. Additionally, we report the first mitochondrial genomes in Annelida that encode genes on both strands. Our findings demonstrate that the supposedly highly conserved mitochondrial gene order suggested for Annelida is restricted to Pleistoannelida, representing the ground pattern of this group. All investigated basal branching annelid taxa show a completely different arrangement of genes than observed in Pleistoannelida. The gene order of protein coding and ribosomal genes in Magelona mirabilis differs only in two transposition events from a putative lophotrochozoan ground pattern and might be the closest to an ancestral annelid pattern. The mitochondrial genomes of Myzostomida show the conserved pattern of Pleistoannelida, thereby supporting their inclusion in this taxon.

The evolution of annelids reveals two adaptive routes to the interstitial realm.

Many animals permanently inhabit the marine interstitium, the space between sand grains [1, 2]. Different evolutionary scenarios may explain the existence of interstitial animals [3, 4]. These scenarios include (1) that the interstitial realm is the ancestral habitat of bilaterians [5, 6], (2) that interstitial taxa evolved from larger ancestors by miniaturization, or (3) progenesis [3]. The first view mirrors the former hypothesis that interstitial annelids, called archiannelids, were at the base of the annelid radiation [7]. Based on morphological data, however, progenesis is generally favored for interstitial annelids today [3, 4, 8]. Herein, our phylogenomic approach revealed that interstitial archiannelids are robustly placed into two groups nested within the annelid tree. Evolution of the first group comprising among others Dinophilidae is best explained by progenesis. In contrast, the second group comprising Protodrilida and Polygordiidae appears to have evolved by stepwise miniaturization adapting from coarser to finer sediments. Thus, in addition to progenesis [3, 4], miniaturization, thought to be too slow for an adaptation to the interstitium [3], is an important second route allowing adaptation to interstitial environments. Both progenesis and miniaturization should be considered when investigating evolution of interstitial taxa [1, 3].

Elucidating the phylogenetic position of Gnathostomulida and first mitochondrial genomes of Gnathostomulida, Gastrotricha and Polycladida (Platyhelminthes).

Gnathostomulida is a taxon of small marine worms, which exclusively inhabit the interstitium. The evolution of Gnathostomulida has been discussed for decades. Originally regarded as primitive animals with affinities to flatworms, the phylogenetic position of Gnathostomulida has been debated. Given the lack of an anus a close relationship to Platyhelminthes has been maintained (i.e., Plathelminthomorpha hypothesis). Alternative hypotheses proposed Gnathostomulida as being close to Gastrotricha due to the presence of a monociliary epidermis (i.e., Monokonta/Neotrichozoa hypothesis) or to Syndermata based on the complicated jaw apparatus (i.e., Gnathifera hypothesis). Molecular analyses using only few genes were inconclusive. Recent phylogenomic studies brought some progress by placing Gnathostomulida as sister to Syndermata, but support for this relationship was low and depended on the analytical strategy. Herein we present the first data of complete or nearly complete mitochondrial genomes for two gnathostomulids (Gnathostomula paradoxa &G. armata), one gastrotrich (Lepidodermella squamata) and one polyclad flatworm (Stylochoplana maculata) to address the uncertain phylogenetic affinity of Gnathostomulida. Our analyses found Gnathostomulida as sister to Syndermata (Gnathifera hypothesis). Thorough sensitivity analyses addressing taxon instability, branch length heterogeneity (also known as long branch attraction) and base composition heterogeneity showed that the position of Gnathostomulida is consistent across the different analyses and, hence, independent of potential misleading biases. Moreover, by ameliorating these different biases nodal support values could be increased to maximum values. Thus, our data support the hypothesis that the different jaw apparatuses of Syndermata and Gnathostomulida are indeed homologous structures as proposed by the Gnathifera hypothesis.

Platyzoan paraphyly based on phylogenomic data supports a noncoelomate ancestry of spiralia.

Based on molecular data three major clades have been recognized within Bilateria: Deuterostomia, Ecdysozoa, and Spiralia. Within Spiralia, small-sized and simply organized animals such as flatworms, gastrotrichs, and gnathostomulids have recently been grouped together as Platyzoa. However, the representation of putative platyzoans was low in the respective molecular phylogenetic studies, in terms of both, taxon number and sequence data. Furthermore, increased substitution rates in platyzoan taxa raised the possibility that monophyletic Platyzoa represents an artifact due to long-branch attraction. In order to overcome such problems, we employed a phylogenomic approach, thereby substantially increasing 1) the number of sampled species within Platyzoa and 2) species-specific sequence coverage in data sets of up to 82,162 amino acid positions. Using established and new measures (long-branch score), we disentangled phylogenetic signal from misleading effects such as long-branch attraction. In doing so, our phylogenomic analyses did not recover a monophyletic origin of platyzoan taxa that, instead, appeared paraphyletic with respect to the other spiralians. Platyhelminthes and Gastrotricha formed a monophylum, which we name Rouphozoa. To the exclusion of Gnathifera, Rouphozoa and all other spiralians represent a monophyletic group, which we name Platytrochozoa. Platyzoan paraphyly suggests that the last common ancestor of Spiralia was a simple-bodied organism lacking coelomic cavities, segmentation, and complex brain structures, and that more complex animals such as annelids evolved from such a simply organized ancestor. This conclusion contradicts alternative evolutionary scenarios proposing an annelid-like ancestor of Bilateria and Spiralia and several independent events of secondary reduction.

A comprehensive analysis of bilaterian mitochondrial genomes and phylogeny.

About 2800 mitochondrial genomes of Metazoa are present in NCBI RefSeq today, two thirds belonging to vertebrates. Metazoan phylogeny was recently challenged by large scale EST approaches (phylogenomics), stabilizing classical nodes while simultaneously supporting new sister group hypotheses. The use of mitochondrial data in deep phylogeny analyses was often criticized because of high substitution rates on nucleotides, large differences in amino acid substitution rate between taxa, and biases in nucleotide frequencies. Nevertheless, mitochondrial genome data might still be promising as it allows for a larger taxon sampling, while presenting a smaller amount of sequence information. We present the most comprehensive analysis of bilaterian relationships based on mitochondrial genome data. The analyzed data set comprises more than 650 mitochondrial genomes that have been chosen to represent a profound sample of the phylogenetic as well as sequence diversity. The results are based on high quality amino acid alignments obtained from a complete reannotation of the mitogenomic sequences from NCBI RefSeq database. However, the results failed to give support for many otherwise undisputed high-ranking taxa, like Mollusca, Hexapoda, Arthropoda, and suffer from extreme long branches of Nematoda, Platyhelminthes, and some other taxa. In order to identify the sources of misleading phylogenetic signals, we discuss several problems associated with mitochondrial genome data sets, e.g. the nucleotide and amino acid landscapes and a strong correlation of gene rearrangements with long branches.

Mitochondrial genomes to the rescue--Diurodrilidae in the myzostomid trap.

Diurodrilidae is a taxon of Lophotrochozoa comprising about six, exclusively interstitial species, which are up to 500µm long and dorsoventrally flattened. Traditionally, Diurodrilidae had been regarded as an annelid family. However, recently Diurodrilidae had been excluded from Annelida and been placed in closer relationship to platyzoan taxa based on both morphological and nuclear rRNA data. Since both, Diurodrilidae and platyzoan taxa, exhibit long branches in the molecular analyses, the close relationship might be due to a long branch attraction artifact. The annelid taxon Myzostomida had been trapped in a similar long branch attraction artifact with platyzoan taxa using nuclear rRNA data, but determination of the nearly complete mitochondrial genome of myzostomids revealed their annelid affinity. Therefore, we determined the nearly complete mitochondrial genome of Diurodrilus subterraneus as well as new nuclear rRNA data for D. subterraneus and some platyzoan taxa. All our analyses of nuclear rRNA and mitochondrial sequence and gene order data presented herein clearly place Diurodrilidae within Annelida and with strong nodal support values in some analyses. Therefore, the previously suggested exclusion of Diurodrilidae from Annelida and its close relationship with platyzoan taxa can be attributed to a long branch artifact. Morphological data do not unambiguously support a platyzoan affinity of Diurodrilidae, but instead would also be in line with a progenetic origin of Diurodrilidae within Annelida.

The impact of mitochondrial genome analyses on the understanding of deuterostome phylogeny.

Deuterostomia, one of the three major lineages of Bilateria, comprises many well-known animals such as vertebrates, sea squirts, sea stars and sea urchins. Whereas monophyly of Deuterostomia and several subtaxa is well supported, the relationships of these to each other and, hence, deuterostome relationships are still uncertain. To address these issues in deuterostome phylogeny we analyzed datasets comprising more than 300 complete deuterostome mitochondrial genomes. Based on sequence information, the results revealed support for several relationships such as a basal position of Xenoturbella within Deuterostomia or for taxa like Craniota or Ambulacraria, but yielded also problems in some taxa, e.g. Tunicata, Pterobranchia and Ophiuroidea, due to long-branch artifacts. However, within tunicates the relationships are well supported. Variation in the genetic code was also informative and, e.g., supported the taxon Ambulacraria including Pterobranchia.

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.