Mitogenomics at the base of Metazoa

Mito-nuclear discordance within Anthozoa, with notes on unique properties of their mitochondrial genomes

Whole mitochondrial genomes are often used in phylogenetic reconstruction. However, discordant patterns in species relationships between mitochondrial and nuclear phylogenies are commonly observed. Within Anthozoa (Phylum Cnidaria), mitochondrial (mt)-nuclear discordance has not yet been examined using a large and comparable dataset. Here, we used data obtained from target-capture enrichment sequencing to assemble and annotate mt genomes and reconstruct phylogenies for comparisons to phylogenies inferred from hundreds of nuclear loci obtained from the same samples. The datasets comprised 108 hexacorals and 94 octocorals representing all orders and > 50% of extant families. Results indicated rampant discordance between datasets at every taxonomic level. This discordance is not attributable to substitution saturation, but rather likely caused by introgressive hybridization and unique properties of mt genomes, including slow rates of evolution driven by strong purifying selection and substitution rate variation. Strong purifying selection across the mt genomes caution their use in analyses that rely on assumptions of neutrality. Furthermore, unique properties of the mt genomes were noted, including genome rearrangements and the presence of nad5 introns. Specifically, we note the presence of the homing endonuclease in ceriantharians. This large dataset of mitochondrial genomes further demonstrates the utility of off-target reads generated from target-capture data for mt genome assembly and adds to the growing knowledge of anthozoan evolution.

Myxozoans (Cnidaria) do not Retain Key Oxygen-Sensing and Homeostasis Toolkit Genes

For aerobic organisms, both the hypoxia-inducible factor pathway and the mitochondrial genomes are key players in regulating oxygen homeostasis. Recent work has suggested that these mechanisms are not as highly conserved as previously thought, prompting more surveys across animal taxonomic levels, which would permit testing of hypotheses about the ecological conditions facilitating evolutionary loss of such genes. The Phylum Cnidaria is known to harbor wide variation in mitochondrial chromosome morphology, including an extreme example, in the Myxozoa, of mitochondrial genome loss. Because myxozoans are obligate endoparasites, frequently encountering hypoxic environments, we hypothesize that variation in environmental oxygen availability could be a key determinant in the evolution of metabolic gene networks associated with oxygen-sensing, hypoxia-response, and energy production. Here, we surveyed genomes and transcriptomes across 46 cnidarian species for the presence of HIF pathway members, as well as for an assortment of hypoxia, mitochondrial, and stress-response toolkit genes. We find that presence of the HIF pathway, as well as number of genes associated with mitochondria, hypoxia, and stress response, do not vary in parallel to mitochondrial genome morphology. More interestingly, we uncover evidence that myxozoans have lost the canonical HIF pathway repression machinery, potentially altering HIF pathway functionality to work under the specific conditions of their parasitic lifestyles. In addition, relative to other cnidarians, myxozoans show loss of large proportions of genes associated with the mitochondrion and involved in response to hypoxia and general stress. Our results provide additional evidence that the HIF regulatory machinery is evolutionarily labile and that variations in the canonical system have evolved in many animal groups.

Mitochondrial genomes of 10 Mantidae species and their phylogenetic implications

This article aims to present a phylogenetic evaluation of Mantidae based on a mitochondrial genome (mitogenome) data set. The mitogenome of 10 Mantidae species were sequenced using next-generation sequencing. The length of nine the complete mitogenomes ranged from 15,371 bp in Tenodera aridifolia to 16,063 bp in Hierodula longa. Mantidae mitogenomes have 37 genes and control region with two exceptions: five trnR copies in Statilia maculata, and H. zhangi was incomplete missing trnI, trnQ, trnM and a portion of the control region. There was a large noncoding region (LNC) between trnM and nad2 in H. chinensis, H. longa, H. maculata and Titanodula sp. Most of protein-coding genes (PCGs) used the typical start ATN codon and TAA/TAG stop codons. All tRNAs fold into the typical clover-leaf secondary structure except trnS1 which lacks a dihydrouracil (DHU) arm. Nucleotide diversity and Ka/Ks analysis of 13 PCGs showed that atp8 had the highest variability and fastest evolutionary rate. Phylogenetic relationships among 42 Mantidae species were reconstructed using the 13 PCGs and two rRNA genes using Bayesian Inference (BI) and Maximum Likelihood (ML) methods. Of the seven mantid subfamilies included in this analysis, only four had multiple exemplars, and of those only Mantinae and Vatinae formed monophyletic groups in BI and ML trees. Consistent with previous studies, the monophyly of the Hierudulinae and Tenoderinae were not been supported. The present results imply that it is necessary to combine nuclear molecular markers and external characteristic to understand the phylogenetic relationships within Mantidae.

Comparative Analysis of the Mitochondrial Genome of Galatheanthemum sp. MT-2020 (Actiniaria Galatheanthemidae) From a Depth of 9,462 m at the Mariana Trench

The hadal zone, which represents the deepest marine habitat on Earth (6,000–11,000 m), is a harsh environment mainly characterized by extremely high hydrostatic pressure, and this habitat is believed to have a high degree of endemism. The deep-sea anemone family Galatheanthemidae comprises two valid species exclusively from the hadal; however, no other information about this family is currently available. In the present study, a sea anemone was collected from a depth of 9,462 m at the Mariana Trench and was defined as Galatheanthemum sp. MT-2020 (Actiniaria Galatheanthemidae). The mitochondrial genome of Galatheanthemum sp. MT-2020 was circular, was 16,633 bp in length, and contained two ribosomal RNA genes, 13 protein-coding genes and two transfer RNA genes. The order of the genes of Galatheanthemum sp. MT-2020 was identical to that of the majority of the species of the order Actiniaria. The value of the AT-skew was the lowest in the whole mitochondrial genome, with a positive GC skew value for the atp8 gene, while other species, except Antholoba achates, had the negative values of the GC skew. Galatheanthemum sp. MT-2020 was clustered with another abyssal species, Paraphelliactis xishaensis, in the phylogenetic tree, and these species diverged in the early Jurassic approximately 200 Mya from the shallow-sea species. The usage ratio of valine, which is one of the five amino acids with the strongest barophilic properties, in the mitochondrial genomes of the two abyssal species was significantly higher than that in other species with habitats above the depth of 3,000 m. The ω (dN/dS) ratio of the genomes was 2.45-fold higher than that of the shallow-sea species, indicating a slower evolutionary rate. Overall, the present study is the first to provide a complete mitogenome of sea anemones from the hadal and reveal some characteristics that may be associated with adaptation to an extreme environment.

Contrasting modes of mitochondrial genome evolution in sister taxa of wood-eating marine bivalves (Teredinidae and Xylophagaidae)

The bivalve families Teredinidae and Xylophagaidae include voracious consumers of wood in shallow and deep-water marine environments, respectively. The taxa are sister clades whose members consume wood as food with the aid of intracellular cellulolytic endosymbionts housed in their gills. This combination of adaptations is found in no other group of animals and was likely present in the common ancestor of both families. Despite these commonalities, the two families have followed dramatically different evolutionary paths with respect to anatomy, life history and distribution. Here we present 42 new mitochondrial genome sequences from Teredinidae and Xylophagaidae and show that distinct trajectories have also occurred in the evolution and organization of their mitochondrial genomes. Teredinidae display significantly greater rates of amino acid substitution but absolute conservation of protein-coding gene order, whereas Xylophagaidae display significantly less amino acid change but have undergone numerous and diverse changes in genome organization since their divergence from a common ancestor. As with many bivalves, these mitochondrial genomes encode two ribosomal RNAs, 12 protein coding genes, and 22 tRNAs; atp8 was not detected. We further show that their phylogeny, as inferred from amino acid sequences of 12 concatenated mitochondrial protein-coding genes, is largely congruent with those inferred from their nuclear genomes based on 18S and 28S ribosomal RNA sequences. Our results provide a robust phylogenetic framework to explore the tempo and mode of mitochondrial genome evolution and offer directions for future phylogenetic and taxonomic studies of wood-boring bivalves.

Mitogenomics provides new insights into the phylogenetic relationships and evolutionary history of deep-sea sea stars (Asteroidea)

The deep sea (> 200 m) is considered as the largest and most remote biome, which characterized by low temperatures, low oxygen level, scarce food, constant darkness, and high hydrostatic pressure. The sea stars (class Asteroidea) are ecologically important and diverse echinoderms in all of the world’s oceans, occurring from the intertidal to the abyssal zone (to about 6000 m). To date, the phylogeny of the sea stars and the relationships of deep-sea and shallow water groups have not yet been fully resolved. Here, we recovered five mitochondrial genomes of deep-sea asteroids. The A+T content of the mtDNA in deep-sea asteroids were significantly higher than that of the shallow-water groups. The gene orders of the five new mitogenomes were identical to that of other asteroids. The phylogenetic analysis showed that the orders Valvatida, Paxillosida, Forcipulatida are paraphyletic. Velatida was the sister order of all the others and then the cladeValvatida-Spinulosida-Paxillosida-Notomyotida versus Forcipulatida-Brisingida. Deep-sea asteroids were nested in different lineages, instead of a well-supported clade. The tropical Western Pacific was suggested as the original area of asteroids, and the temperate water was initially colonized with asteroids by the migration events from the tropical and cold water. The time-calibrated phylogeny showed that Asteroidea originated during Devonian-Carboniferous boundary and the major lineages of Asteroidea originated during Permian–Triassic boundary. The divergence between the deep-sea and shallow-water asteroids coincided approximately with the Triassic-Jurassic extinction. Total 29 positively selected sites were detected in fifteen mitochondrial genes of five deep-sea lineages, implying a link between deep-sea adaption and mitochondrial molecular biology in asteroids.

Mitochondrial Genome of Nonmodel Marine Metazoans by Next-Generation Sequencing (NGS)

Mitochondrial genomes (mtgenome) represent an important source of information for addressing fundamental evolutionary, phylogeographic, systematic, and ecological questions in marine organisms. In the last two decades the advent of high-throughput next-generation sequencing (NGS) has provided an unprecedented possibility to access large amount of genomic data and, as such, there has been a rapid growth in mtgenome resources and studies. In particular, NGS strategies represent a great advantage for investigating nonmodel marine organisms for which no or limited genomic resources are available. Here, we describe a routinely used standardized protocol to obtain mtgenome of nonmodel marine organisms by NGS. The protocol is composed of five main steps, including DNA extraction, DNA fragmentation, library preparation, high-throughput sequencing, and bioinformatic analyses. Each of the first three steps is followed by size/quality and concentration validations. The advantages of the described protocol rely on the assumption that no a priori information on mtgenome of the studied organism is needed and on its versatility as researchers may choose several kits for DNA extraction and library preparation and adopt different methods for DNA fragmentation depending on their needs, experience, and suppliers.

Comparative analysis of twelve mitogenomes of Caliscelidae (Hemiptera: Fulgoromorpha) and their phylogenetic implications

Here, the complete mitochondrial genomes (mitogenomes) of 12 Caliscelidae species, Augilina tetraina, Augilina triaina, Symplana brevistrata, Symplana lii, Neosymplana vittatum, Pseudosymplanella nigrifasciata, Symplanella brevicephala, Symplanella unipuncta, Augilodes binghami, Cylindratus longicephalus, Caliscelis shandongensis, and Peltonotellus sp., were determined and comparatively analyzed. The genomes varied from 15,424 to 16,746 bp in size, comprising 37 mitochondrial genes and an A+T-rich region. The typical gene content and arrangement were similar to those of most Fulgoroidea species. The nucleotide compositions of the mitogenomes were biased toward A/T. All protein-coding genes (PCGs) started with a canonical ATN or GTG codon and ended with TAN or an incomplete stop codon, single T. Among 13 PCGs in 16 reported Caliscelidae mitogenomes, cox1 and atp8 showed the lowest and highest nucleotide diversity, respectively. All PCGs evolved under purifying selection, with atp8 considered a comparatively fast-evolving gene. Phylogenetic relationships were reconstructed based on 13 PCGs in 16 Caliscelidae species and five outgroups using maximum likelihood and Bayesian inference analyses. All species of Caliscelidae formed a steadily monophyletic group with high support. Peltonotellini was present at the basal position of the phylogenetic tree. Augilini was the sister group to Caliscelini and Peltonotellini.

The complete mitogenome of Lophops carinata (Hemiptera: Fulgoroidea: Lophopidae)

In the present study, the complete mitochondrial genome of Lophops carinata (Hemiptera: Lophopidae) was sequenced for the first time through next-generation sequencing. The complete mitogenome of L. carinata is 15,553 bp in length, with the typical gene content and arrangement usually observed in Hexapods. The mitogenome consisted of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 D-loop. The overall nucleotide composition of the mitogenome was 44.6% A, 14.0% C, 8.3% G, and 33.2% T, with an A + T bias of 77.8%. Phylogenetic analyses from 12 Fulgoroidea species by maximum likelihood were consistent and well supported the basal position of Delphacidae, a close affinity among the families Ricaniidae, Issidae, and Flatidae, and a close relationship between Achilidae and Fulgoridae. And L. carinata belong to a separate lineage, located in the middle of the phylogenetic tree.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

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.

The Mitochondrial Genome of the Sea Anemone Stichodactyla haddoni Reveals Catalytic Introns, Insertion-Like Element, and Unexpected Phylogeny

A hallmark of sea anemone mitochondrial genomes (mitogenomes) is the presence of complex catalytic group I introns. Here, we report the complete mitogenome and corresponding transcriptome of the carpet sea anemone Stichodactyla haddoni (family Stichodactylidae). The mitogenome is vertebrate-like in size, organization, and gene content. Two mitochondrial genes encoding NADH dehydrogenase subunit 5 (ND5) and cytochrome c oxidase subunit I (COI) are interrupted with complex group I introns, and one of the introns (ND5-717) harbors two conventional mitochondrial genes (ND1 and ND3) within its sequence. All the mitochondrial genes, including the group I introns, are expressed at the RNA level. Nonconventional and optional mitochondrial genes are present in the mitogenome of S. haddoni. One of these gene codes for a COI-884 intron homing endonuclease and is organized in-frame with the upstream COI exon. The insertion-like orfA is expressed as RNA and translocated in the mitogenome as compared with other sea anemones. Phylogenetic analyses based on complete nucleotide and derived protein sequences indicate that S. haddoni is embedded within the family Actiniidae, a finding that challenges current taxonomy.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

Mitochondrial Genome Evolution of Placozoans: Gene Rearrangements and Repeat Expansions

Placozoans, nonbilaterian animals with the simplest known metazoan bauplan, are currently classified into 20 haplotypes belonging to three genera, Polyplacotoma, Trichoplax, and Hoilungia. The latter two comprise two and five clades, respectively. In Trichoplax and Hoilungia, previous studies on six haplotypes belonging to four different clades have shown that their mtDNAs are circular chromosomes of 32–43 kb in size, which encode 12 protein-coding genes, 24 tRNAs, and two rRNAs. These mitochondrial genomes (mitogenomes) also show unique features rarely seen in other metazoans, including open reading frames (ORFs) of unknown function, and group I and II introns. Here, we report seven new mitogenomes, covering the five previously described haplotypes H2, H17, H19, H9, and H11, as well as two new haplotypes, H23 (clade III) and H24 (clade VII). The overall gene content is shared between all placozoan mitochondrial genomes, but genome sizes, gene orders, and several exon–intron boundaries vary among clades. Phylogenomic analyses strongly support a tree topology different from previous 16S rRNA analyses, with clade VI as the sister group to all other Hoilungia clades. We found small inverted repeats in all 13 mitochondrial genomes of the Trichoplax and Hoilungia genera and evaluated their distribution patterns among haplotypes. Because Polyplacotoma mediterranea (H0), the sister to the remaining haplotypes, has a small mitochondrial genome with few small inverted repeats and ORFs, we hypothesized that the proliferation of inverted repeats and ORFs substantially contributed to the observed increase in the size and GC content of the Trichoplax and Hoilungia mitochondrial genomes.

Mitogenomics reveals phylogenetic relationships of Arcoida (Mollusca, Bivalvia) and multiple independent expansions and contractions in mitochondrial genome size

Arcoida, comprising about 570 species of blood cockles, is an ecologically and economically important lineage of bivalve molluscs. Current classification of arcoids is largely based on morphology, which shows widespread homoplasy. Despite two recent studies employing multi-locus analyses with broad sampling of Arcoida, evolutionary relationships among major lineages remain controversial. Interestingly, mitochondrial genomes of several ark shell species are 2-3 times larger than those found in most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. These results highlight the need of detailed phylogenetic study to explore evolutionary relationships within Arcoida so that the evolution of mitochondrial genome size can be understood. To this end, we sequenced 17 mitochondrial genomes and compared them with publicly available data, including those from other lineages of Arcoida with emphasis on the subclade Arcoidea species. Our phylogenetic analyses indicate that Noetiidae, Cucullaeidae and Glycymerididae are nested within a polyphyletic Arcidae. Moreover, we find multiple independent expansions and potential contractions of mitochondrial genome size, suggesting that the large mitochondrial genome is not a shared ancestral feature in Arcoida. We also examined tandem repeats and inverted repeats in non-coding regions and investigated the presence of such repeats with relation to genome size variation. Our results suggest that tandem repeats might facilitate intraspecific mitochondrial genome size variation, and that inverted repeats, which could be derived from transposons, might be responsible for mitochondrial genome expansions and contractions. We show that mitochondrial genome size in Arcoida is more dynamic than previously understood and provide insights into evolution of mitochondrial genome size variation in metazoans.

Conserved novel ORFs in the mitochondrial genome of the ctenophore Beroe forskalii

To date, five ctenophore species' mitochondrial genomes have been sequenced, and each contains open reading frames (ORFs) that if translated have no identifiable orthologs. ORFs with no identifiable orthologs are called unidentified reading frames (URFs). If truly protein-coding, ctenophore mitochondrial URFs represent a little understood path in early-diverging metazoan mitochondrial evolution and metabolism. We sequenced and annotated the mitochondrial genomes of three individuals of the beroid ctenophore Beroe forskalii and found that in addition to sharing the same canonical mitochondrial genes as other ctenophores, the B. forskalii mitochondrial genome contains two URFs. These URFs are conserved among the three individuals but not found in other sequenced species. We developed computational tools called pauvre and cuttlery to determine the likelihood that URFs are protein coding. There is evidence that the two URFs are under negative selection, and a novel Bayesian hypothesis test of trinucleotide frequency shows that the URFs are more similar to known coding genes than noncoding intergenic sequence. Protein structure and function prediction of all ctenophore URFs suggests that they all code for transmembrane transport proteins. These findings, along with the presence of URFs in other sequenced ctenophore mitochondrial genomes, suggest that ctenophores may have uncharacterized transmembrane proteins present in their mitochondria.

Mitogenomics suggests a sister relationship of Relicanthus daphneae (Cnidaria: Anthozoa: Hexacorallia: incerti ordinis) with Actiniaria

Relicanthus daphneae (formerly Boloceroides daphneae) was first described in 2006 as a giant sea anemone based on morphology. In 2014, its classification was challenged based on molecular data: using five genes, Relicanthus was resolved sister to zoanthideans, but with mixed support. To better understand the evolutionary relationship of Relicanthus with other early-branching metazoans, we present 15 newly-sequenced sea anemone mitochondrial genomes and a mitogenome-based phylogeny including all major cnidarian groups, sponges, and placozoans. Our phylogenetic reconstruction reveals a moderately supported sister relationship between Relicanthus and the Actiniaria. Morphologically, the cnidae of Relicanthus has apical flaps, the only existing synapomorphy for sea anemones. Based on both molecular and morphological results, we propose a third suborder (Helenmonae) within the Actiniaria to accommodate Relicanthus. Although Relicanthus shares the same gene order and content with other available actiniarian mitogenomes, it is clearly distinct at the nucleotide level from anemones within the existing suborders. The phylogenetic position of Relicanthus could reflect its association with the periphery of isolated hydrothermal vents, which, although patchy and ephemeral, harbor unique chemosynthetic communities that provide a relatively stable food source to higher trophic levels over long evolutionary timescales. The ability to colonize the deep sea and the periphery of new vent systems may be facilitated by Relicanthus’ large and extremely yolky eggs.

What does mitogenomics tell us about the evolutionary history of the Drosophila buzzatii cluster (repleta group)?

The Drosophila repleta group is an array of more than 100 species endemic to the “New World”, many of which are cactophilic. The ability to utilize decaying cactus tissues as breeding and feeding sites is a key aspect that allowed the successful diversification of the repleta group in American deserts and arid lands. Within this group, the Drosophila buzzatii cluster is a South American clade of seven closely related species in different stages of divergence, making them a valuable model system for evolutionary research. Substantial effort has been devoted to elucidating the phylogenetic relationships among members of the D. buzzatii cluster, including molecular phylogenetic studies that have generated ambiguous results where different tree topologies have resulted dependent on the kinds of molecular marker used. Even though mitochondrial DNA regions have become useful markers in evolutionary biology and population genetics, none of the more than twenty Drosophila mitogenomes assembled so far includes this cluster. Here, we report the assembly of six complete mitogenomes of five species: D. antonietae, D. borborema, D. buzzatii, two strains of D. koepferae and D. seriema, with the aim of revisiting phylogenetic relationships and divergence times by means of mitogenomic analyses. Our recovered topology using complete mitogenomes supports the hypothesis of monophyly of the D. buzzatii cluster and shows two main clades, one including D. buzzatii and D. koepferae (both strains), and the other containing the remaining species. These results are in agreement with previous reports based on a few mitochondrial and/or nuclear genes, but conflict with the results of a recent large-scale nuclear phylogeny, indicating that nuclear and mitochondrial genomes depict different evolutionary histories.

Divergence history and hydrothermal vent adaptation of decapod crustaceans: A mitogenomic perspective

Decapod crustaceans, such as alvinocaridid shrimps, bythograeid crabs and galatheid squat lobsters are important fauna in the hydrothermal vents and have well adapted to hydrothermal vent environments. In this study, eighteen mitochondrial genomes (mitogenomes) of hydrothermal vent decapods were used to explore the evolutionary history and their adaptation to the hydrothermal vent habitats. BI and ML algorithms produced consistent phylogeny for Decapoda. The phylogenetic relationship revealed more evolved positions for all the hydrothermal vent groups, indicating they migrated from non-vent environments, instead of the remnants of ancient hydrothermal vent species, which support the extinction/repopulation hypothesis. The divergence time estimation on the Alvinocarididae, Bythograeidae and Galatheoidea nodes are located at 75.20, 56.44 and 47.41–50.43 Ma, respectively, which refers to the Late Cretaceous origin of alvinocaridid shrimps and the Early Tertiary origin of bythograeid crabs and galatheid squat lobsters. These origin stories are thought to associate with the global deep-water anoxic/dysoxic events. Total eleven positively selected sites were detected in the mitochondrial OXPHOS genes of three lineages of hydrothermal vent decapods, suggesting a link between hydrothermal vent adaption and OXPHOS molecular biology in decapods. This study adds to the understanding of the link between mitogenome evolution and ecological adaptation to hydrothermal vent habitats in decapods.

The complete mitochondrial genome of the Ctenophore Beroe cucumis, a mitochondrial genome showing rapid evolutionary rates

We described the complete mitochondrial genome of the Ctenophore Beroe cucumis, which is a circular molecule of 10,487 bp in length. The new mitochondrial genome comprised only 12 genes, making it one of the smallest animals' mtDNA. Both nucleotide substitution and gene order rearrangements exhibited extreme high evolutionary rate in mitogenomes of Ctenophore. The phylogenetic analysis based on mitogenomics failed to reveal the basal position of Ctenophore within metazoan, owing to the extreme evolutionary rate. Based on the available Ctenophora mitogenomes, we found the optimized primers designed by Geller et al. for DNA barcoding suited for the taxon.

Giant group I intron in a mitochondrial genome is removed by RNA back-splicing

Background

The mitochondrial genomes of mushroom corals (Corallimorpharia) are remarkable for harboring two complex group I introns; ND5-717 and COI-884. How these autocatalytic RNA elements interfere with mitochondrial RNA processing is currently not known. Here, we report experimental support for unconventional processing events of ND5-717 containing RNA.

Results

We obtained the complete mitochondrial genome sequences and corresponding mitochondrial transcriptomes of the two distantly related corallimorpharian species Ricordea yuma and Amplexidiscus fenestrafer. All mitochondrial genes were found to be expressed at the RNA-level. Both introns were perfectly removed by autocatalytic splicing, but COI-884 excision appeared more efficient than ND5-717. ND5-717 was organized into giant group I intron elements of 18.1 kb and 19.3 kb in A. fenestrafer and R. yuma, respectively. The intron harbored almost the entire mitochondrial genome embedded within the P8 peripheral segment.

Conclusion

ND5-717 was removed by group I intron splicing from a small primary transcript that contained a permutated intron–exon arrangement. The splicing pathway involved a circular exon-containing RNA intermediate, which is a hallmark of RNA back-splicing. ND5-717 represents the first reported natural group I intron that becomes excised by back-splicing from a permuted precursor RNA. Back-splicing may explain why Corallimorpharia mitochondrial genomes tolerate giant group I introns.

Electronic supplementary material

The online version of this article (10.1186/s12867-019-0134-y) contains supplementary material, which is available to authorized users.

Linear Mitochondrial Genome in Anthozoa (Cnidaria): A Case Study in Ceriantharia

Sequences and structural attributes of mitochondrial genomes have played a critical role in the clarification of relationships among Cnidaria, a key phylum of early-diverging animals. Among the major lineages of Cnidaria, Ceriantharia (“tube anemones”) remains one of the most enigmatic in terms of its phylogenetic position. We sequenced the mitochondrial genomes of two ceriantharians to see whether the complete organellar genome would provide more support for the phylogenetic placement of Ceriantharia. For both Isarachnanthus nocturnus and Pachycerianthus magnus, the mitochondrial gene sequences could not be assembled into a single circular genome. Instead, our analyses suggest that both species have mitochondrial genomes consisting of multiple linear fragments. Linear mitogenomes are characteristic of members of Medusozoa, one of the major lineages of Cnidaria, but are unreported for Anthozoa, which includes the Ceriantharia. The inferred number of fragments and variation in gene order between species is much greater within Ceriantharia than among the lineages of Medusozoa. We identify origins of replication for each of the five putative chromosomes of the Isarachnanthus nocturnus mitogenome and for each of the eight putative chromosomes of the Pachycerianthus magnus mitogenome. At 80,923 bp, I. nocturnus now holds the record for the largest animal mitochondrial genome reported to date. The novelty of the mitogenomic structure in Ceriantharia highlights the distinctiveness of this lineage but, because it appears to be both unique to and diverse within Ceriantharia, it is uninformative about the phylogenetic position of Ceriantharia relative to other Anthozoa. The presence of tRNA^Met and tRNA^Trp in both ceriantharian mitogenomes supports a closer relationship between Ceriantharia and Hexacorallia than between Ceriantharia and any other cnidarian lineage, but phylogenetic analysis of the genes contained in the mitogenomes suggests that Ceriantharia is sister to a clade containing Octocorallia + Hexacorallia indicating a possible suppression of tRNA^Trp in Octocorallia.

Comparative analysis of the complete mitochondrial genomes of five Achilidae species (Hemiptera: Fulgoroidea) and other Fulgoroidea reveals conserved mitochondrial genome organization

In the present study, the complete mitochondrial genomes (mitogenomes) of five Achilidae (Hemiptera: Fulgoroidea), Betatropis formosana, two new species (Magadhaideus luodiana sp. nov and Peltatavertexalis horizontalis sp. nov), Plectoderini sp. and Paracatonidia sp., were sequenced for the first time through next-generation sequencing. The five mitogenomes ranged from 15,214 to 16,216 bp in length, with the typical gene content and arrangement usually observed in Hexapods. The motif "ATGATAA" between atp8 and atp6 was found in all the analyzed species. An overlap "AAGCTTA" between trnW and trnC was observed in the mitogenomes of most Fulgoroidea. The structural and compositional analyses of 26 Fulgoroidea mitogenomes, including the gene rearrangement of five tRNAs (trnW, trnC and trnY; trnT and trnP), the A + T content and AT-skew of the whole mitogenomes, and the nuclear acid and amino acid compositions of the protein-coding genes (PCGs), revealed family-level differences between Delphacidae and other families (Achilidae, Flatidae, Fulgoridae, Issidae and Ricaniidae). Phylogenetic analyses of 13 protein-coding genes from 26 Fulgoroidea species by maximum likelihood and Bayesian Inference were consistent and well supported the basal position of Delphacidae, a close affinity among the families Flatidae, Issidae and Ricaniidae, and a close relationship between Achilidae and Fulgoridae.

Complete mitochondrial genomes of three Oxya grasshoppers (Orthoptera) and their implications for phylogenetic reconstruction

Oxya is a genus of grasshoppers (Orthoptera: Acridoidea) attacking rice and other gramineous plants in Africa and Asia. In the present study, we characterized complete mitochondrial genomes (mitogenomes) of three species, Oxya japonica japonica (15,427 bp), Oxya hainanensis (15,443 bp) and Oxya agavisa robusta (15,552 bp) collected from China. The three mitogenomes contained a typical gene set of metazoan mitogenomes and shared the same gene order with other Acridid grasshoppers, including the rearrangement of tRNA^Asp and tRNA^Lys. Analyses of pairwise genetic distances showed that ATP8 was the least conserved gene, while COI the most conserved. To determine the position of Oxya grasshoppers in the phylogeny of Acrididae, we reconstructed phylogenetic trees among 64 species from across 11 subfamilies using nucleotide sequences of mitogenomes. While the tree confirms traditional classifications of Acrididae at major higher-levels, it suggests a few modifications for classifications at lower-levels.

Deep-water sea anemone with a two-chromosome mitochondrial genome

Mitochondrial genome organization of sea anemones appears conserved among species and families, and is represented by a single circular DNA molecule of 17 to 21 kb. The mitochondrial gene content corresponds to the same 13 protein components of the oxidative phosphorylation (OxPhos) system as in vertebrates. Hallmarks, however, include a highly reduced tRNA gene repertoire and the presence of autocatalytic group I introns. Here we demonstrate that the mitochondrial genome of the deep-water sea anemone Protanthea simplex deviates significantly from that of other known sea anemones. The P. simplex mitochondrial genome contains a heavily scrambled order of genes that are coded on both DNA strands and organized along two circular mito-chromosomes, MCh-I and MCh-II. We found MCh-I to be representative of the prototypic sea anemone mitochondrial genome, encoding 12 OxPhos proteins, two ribosomal RNAs, two transfer RNAs, and a group I intron. In contrast, MCh-II was found to be a laterally transferred plasmid-like DNA carrying the conserved cytochrome oxidase II gene and a second allele of the small subunit ribosomal RNA gene.

Mitogenomics Reveals a Novel Genetic Code in Hemichordata

The diverse array of codon reassignments demonstrate that the genetic code is not universal in nature. Exploring mechanisms underlying codon reassignment is critical for understanding the evolution of the genetic code during translation. Hemichordata, comprising worm-like Enteropneusta and colonial filter-feeding Pterobranchia, is the sister taxon of echinoderms and is more distantly related to chordates. However, only a few hemichordate mitochondrial genomes have been sequenced, hindering our understanding of mitochondrial genome evolution within Deuterostomia. In this study, we sequenced four mitochondrial genomes and two transcriptomes, including representatives of both major hemichordate lineages and analyzed together with public available data. Contrary to the current understanding of the mitochondrial genetic code in hemichordates, our comparative analyses suggest that UAA encodes Tyr instead of a "Stop" codon in the pterobranch lineage Cephalodiscidae. We also predict that AAA encodes Lys in pterobranch and enteropneust mitochondrial genomes, contradicting the previous assumption that hemichordates share the same genetic code with echinoderms for which AAA encodes Asn. Thus, we propose a new mitochondrial genetic code for Cephalodiscus and a revised code for enteropneusts. Moreover, our phylogenetic analyses are largely consistent with previous phylogenomic studies. The only exception is the phylogenetic position of the enteropneust Stereobalanus, whose placement as sister to all other described enteropneusts. With broader taxonomic sampling, we provide evidence that evolution of mitochondrial gene order and genetic codes in Hemichordata are more dynamic than previously thought and these findings provide insights into mitochondrial genome evolution within this clade.

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.

Expression of homing endonuclease gene and insertion-like element in sea anemone mitochondrial genomes: Lesson learned from Anemonia viridis

The mitochondrial genomes of sea anemones are dynamic in structure. Invasion by genetic elements, such as self-catalytic group I introns or insertion-like sequences, contribute to sea anemone mitochondrial genome expansion and complexity. By using next generation sequencing we investigated the complete mtDNAs and corresponding transcriptomes of the temperate sea anemone Anemonia viridis and its closer tropical relative Anemonia majano. Two versions of fused homing endonuclease gene (HEG) organization were observed among the Actiniidae sea anemones; in-frame gene fusion and pseudo-gene fusion. We provided support for the pseudo-gene fusion organization in Anemonia species, resulting in a repressed HEG from the COI-884 group I intron. orfA, a putative protein-coding gene with insertion-like features, was present in both Anemonia species. Interestingly, orfA and COI expression were significantly up-regulated upon long-term environmental stress corresponding to low seawater pH conditions. This study provides new insights to the dynamics of sea anemone mitochondrial genome structure and function.

More evolution underground: Accelerated mitochondrial substitution rate in Australian burrowing freshwater crayfishes (Decapoda: Parastacidae)

To further understand the evolutionary history and mitogenomic features of Australia's highly distinctive freshwater crayfish fauna, we utilized a recently described rapid mitogenome sequencing pipeline to generate 24 new crayfish mitogenomes including a diversity of burrowing crayfish species and the first for Astacopsis gouldi, the world's largest freshwater invertebrate. Whole mitogenome-based phylogeny estimates using both Bayesian and Maximum Likelihood methods substantially strengthen existing hypotheses for systematic relationships among Australian freshwater crayfish with evidence of pervasive diversifying selection and accelerated mitochondrial substitution rate among the members of the clade representing strongly burrowing crayfish that may reflect selection pressures for increased energy requirement for adaptation to terrestrial environment and a burrowing lifestyle. Further, gene rearrangements are prevalent in the burrowing crayfish mitogenomes involving both tRNA and protein coding genes. In addition, duplicated control regions were observed in two closely related Engaeus species, together with evidence for concerted evolution. This study significantly adds to the understanding of Australian freshwater crayfish evolutionary relationships and suggests a link between mitogenome evolution and adaptation to terrestrial environments and a burrowing lifestyle in freshwater crayfish.

The mitochondrial genome of a sea anemone Bolocera sp. exhibits novel genetic structures potentially involved in adaptation to the deep-sea environment

The deep sea is one of the most extensive ecosystems on earth. Organisms living there survive in an extremely harsh environment, and their mitochondrial energy metabolism might be a result of evolution. As one of the most important organelles, mitochondria generate energy through energy metabolism and play an important role in almost all biological activities. In this study, the mitogenome of a deep‐sea sea anemone (Bolocera sp.) was sequenced and characterized. Like other metazoans, it contained 13 energy pathway protein‐coding genes and two ribosomal RNAs. However, it also exhibited some unique features: just two transfer RNA genes, two group I introns, two transposon‐like noncanonical open reading frames (ORFs), and a control region‐like (CR‐like) element. All of the mitochondrial genes were coded by the same strand (the H‐strand). The genetic order and orientation were identical to those of most sequenced actiniarians. Phylogenetic analyses showed that this species was closely related to Bolocera tuediae. Positive selection analysis showed that three residues (31 L and 42 N in ATP6, 570 S in ND5) of Bolocera sp. were positively selected sites. By comparing these features with those of shallow sea anemone species, we deduced that these novel gene features may influence the activity of mitochondrial genes. This study may provide some clues regarding the adaptation of Bolocera sp. to the deep‐sea environment.

Deep RNA sequencing reveals the smallest known mitochondrial micro exon in animals: The placozoan cox1 single base pair exon

The phylum Placozoa holds a key position for our understanding of the evolution of mitochondrial genomes in Metazoa. Placozoans possess large mitochondrial genomes which harbor several remarkable characteristics such as a fragmented cox1 gene and trans-splicing cox1 introns. A previous study also suggested the existence of cox1 mRNA editing in Trichoplax adhaerens, yet the only formally described species in the phylum Placozoa. We have analyzed RNA-seq data of the undescribed sister species, Placozoa sp. H2 ("Panama" clone), with special focus on the mitochondrial mRNA. While we did not find support for a previously postulated cox1 mRNA editing mechanism, we surprisingly found two independent transcripts representing intermediate cox1 mRNA splicing stages. Both transcripts consist of partial cox1 exon as well as overlapping intron fragments. The data suggest that the cox1 gene harbors a single base pair (cytosine) micro exon. Furthermore, conserved group I intron structures flank this unique micro exon also in other placozoans. We discuss the evolutionary origin of this micro exon in the context of a self-splicing intron gain in the cox1 gene of the last common ancestor of extant placozoans.

Evolution of group I introns in Porifera: new evidence for intron mobility and implications for DNA barcoding

BACKGROUND

Mitochondrial introns intermit coding regions of genes and feature characteristic secondary structures and splicing mechanisms. In metazoans, mitochondrial introns have only been detected in sponges, cnidarians, placozoans and one annelid species. Within demosponges, group I and group II introns are present in six families. Based on different insertion sites within the cox1 gene and secondary structures, four types of group I and two types of group II introns are known, which can harbor up to three encoding homing endonuclease genes (HEG) of the LAGLIDADG family (group I) and/or reverse transcriptase (group II). However, only little is known about sponge intron mobility, transmission, and origin due to the lack of a comprehensive dataset. We analyzed the largest dataset on sponge mitochondrial group I introns to date: 95 specimens, from 11 different sponge genera which provided novel insights into the evolution of group I introns.

RESULTS

For the first time group I introns were detected in four genera of the sponge family Scleritodermidae (Scleritoderma, Microscleroderma, Aciculites, Setidium). We demonstrated that group I introns in sponges aggregate in the most conserved regions of cox1. We showed that co-occurrence of two introns in cox1 is unique among metazoans, but not uncommon in sponges. However, this combination always associates an active intron with a degenerating one. Earlier hypotheses of HGT were confirmed and for the first time VGT and secondary losses of introns conclusively demonstrated.

CONCLUSION

This study validates the subclass Spirophorina (Tetractinellida) as an intron hotspot in sponges. Our analyses confirm that most sponge group I introns probably originated from fungi. DNA barcoding is discussed and the application of alternative primers suggested.

Animal Mitochondrial DNA as We Do Not Know It: mt-Genome Organization and Evolution in Nonbilaterian Lineages

Animal mitochondrial DNA (mtDNA) is commonly described as a small, circular molecule that is conserved in size, gene content, and organization. Data collected in the last decade have challenged this view by revealing considerable diversity in animal mitochondrial genome organization. Much of this diversity has been found in nonbilaterian animals (phyla Cnidaria, Ctenophora, Placozoa, and Porifera), which, from a phylogenetic perspective, form the main branches of the animal tree along with Bilateria. Within these groups, mt-genomes are characterized by varying numbers of both linear and circular chromosomes, extra genes (e.g. atp9, polB, tatC), large variation in the number of encoded mitochondrial transfer RNAs (tRNAs) (0-25), at least seven different genetic codes, presence/absence of introns, tRNA and mRNA editing, fragmented ribosomal RNA genes, translational frameshifting, highly variable substitution rates, and a large range of genome sizes. This newly discovered diversity allows a better understanding of the evolutionary plasticity and conservation of animal mtDNA and provides insights into the molecular and evolutionary mechanisms shaping mitochondrial genomes.

Mitochondrial Genomes of Anopheles (Kerteszia) (Diptera: Culicidae) From the Atlantic Forest, Brazil

Mitochondrial genome sequences are widely used as molecular markers for phylogenetic studies of mosquito species complexes, such as the Anopheles albitarsis complex. Except for a few studies that employed a limited number of nuclear or mitochondrial loci to address the genetic structure and species status of Anopheles cruzii, Anopheles bellator, and Anopheles homunculus, little is known about genetic markers that can be employed in studies focusing on Kerteszia species. The complete mitochondrial genomes of seven specimens of An. bellator, An. cruzii, An. homunculus, and Anopheles laneanus were sequenced using long-range polymerase chain reaction and Illumina sequencing. The mitochondrial genomes varied from 15,446 to 15,738 bp in length and contained 37 genes (13 protein-encoding genes, 2 rRNA genes [12S rRNA and 16S rRNA] and 22 tRNA genes), and the AT-rich control region, as all do other Anopheles mitochondrial genomes sequenced to date. Specimens from four populations of An. cruzii showed differences in codon composition.

Comparative study of Hippo pathway genes in cellular conveyor belts of a ctenophore and a cnidarian

Background

The Hippo pathway regulates growth rate and organ size in fly and mouse, notably through control of cell proliferation. Molecular interactions at the heart of this pathway are known to have originated in the unicellular ancestry of metazoans. They notably involve a cascade of phosphorylations triggered by the kinase Hippo, with subsequent nuclear to cytoplasmic shift of Yorkie localisation, preventing its binding to the transcription factor Scalloped, thereby silencing proliferation genes. There are few comparative expression data of Hippo pathway genes in non-model animal species and notably none in non-bilaterian phyla.

Results

All core Hippo pathway genes could be retrieved from the ctenophore Pleurobrachia pileus and the hydrozoan cnidarian Clytia hemisphaerica, with the important exception of Yorkie in ctenophore. Expression study of the Hippo, Salvador and Scalloped genes in tentacle “cellular conveyor belts” of these two organisms revealed striking differences. In P. pileus, their transcripts were detected in areas where undifferentiated progenitors intensely proliferate and where expression of cyclins B and D was also seen. In C. hemisphaerica, these three genes and Yorkie are expressed not only in the proliferating but also in the differentiation zone of the tentacle bulb and in mature tentacle cells. However, using an antibody designed against the C. hemiphaerica Yorkie protein, we show in two distinct cell lineages of the medusa that Yorkie localisation is predominantly nuclear in areas of active cell proliferation and mainly cytoplasmic elsewhere.

Conclusions

This is the first evidence of nucleocytoplasmic Yorkie shift in association with the arrest of cell proliferation in a cnidarian, strongly evoking the cell division-promoting role of this protein and its inhibition by the activated Hippo pathway in bilaterian models. Our results furthermore highlight important differences in terms of deployment and regulation of Hippo pathway genes between cnidarians and ctenophores.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0041-y) contains supplementary material, which is available to authorized users.

Fifteen new earthworm mitogenomes shed new light on phylogeny within the Pheretima complex

The Pheretima complex within the Megascolecidae family is a major earthworm group. Recently, the systematic status of the Pheretima complex based on morphology was challenged by molecular studies. In this study, we carry out the first comparative mitogenomic study in oligochaetes. The mitogenomes of 15 earthworm species were sequenced and compared with other 9 available earthworm mitogenomes, with the main aim to explore their phylogenetic relationships and test different analytical approaches on phylogeny reconstruction. The general earthworm mitogenomic features revealed to be conservative: all genes encoded on the same strand, all the protein coding loci shared the same initiation codon (ATG), and tRNA genes showed conserved structures. The Drawida japonica mitogenome displayed the highest A + T content, reversed AT/GC-skews and the highest genetic diversity. Genetic distances among protein coding genes displayed their maximum and minimum interspecific values in the ATP8 and CO1 genes, respectively. The 22 tRNAs showed variable substitution patterns between the considered earthworm mitogenomes. The inclusion of rRNAs positively increased phylogenetic support. Furthermore, we tested different trimming tools for alignment improvement. Our analyses rejected reciprocal monophyly among Amynthas and Metaphire and indicated that the two genera should be systematically classified into one.

The complete mitochondrial genome sequence of the tubeworm Lamellibrachia satsuma and structural conservation in the mitochondrial genome control regions of Order Sabellida

The control region of the mitochondrial genomes shows high variation in conserved sequence organizations, which follow distinct evolutionary patterns in different species or taxa. In this study, we sequenced the complete mitochondrial genome of Lamellibrachia satsuma from the cold-seep region of Kagoshima Bay, as a part of whole genome study and extensively studied the structural features and patterns of the control region sequences. We obtained 15,037 bp of mitochondrial genome using Illumina sequencing and identified the non-coding AT-rich region or control region (354 bp, AT=83.9%) located between trnH and trnR. We found 7 conserved sequence blocks (CSB), scattered throughout the control region of L. satsuma and other taxa of Annelida. The poly-TA stretches, which commonly form the stem of multiple stem-loop structures, are most conserved in the CSB-I and CSB-II regions. The mitochondrial genome of L. satsuma encodes a unique repetitive sequence in the control region, which forms a unique secondary structure in comparison to Lamellibrachia luymesi. Phylogenetic analyses of all protein-coding genes indicate that L. satsuma forms a monophyletic clade with L. luymesi along with other tubeworms found in cold-seep regions (genera: Lamellibrachia, Escarpia, and Seepiophila). In general, the control region sequences of Annelida could be aligned with certainty within each genus, and to some extent within the family, but with a higher rate of variation in conserved regions.

The ctenophore lineage is older than sponges? That cannot be right! Or can it?

Recent phylogenetic analyses resulting from collection of whole genome data suggest that ctenophores, or comb jellies, are sister to all other animals. Even before publication, this result prompted discussion among researchers. Here, I counter common criticisms raised about this result and show that assumptions placing sponges as the basal-most extant animal lineage are based on limited evidence and questionable premises. For example, the idea that sponges are simple and the reported similarity of sponge choanocytes to Choanflagellata do not provide useful characters for determining the positions of sponges within the animal tree. Intertwined with discussion of basal metazoan phylogeny is consideration of the evolution of neuronal systems. Recent data show that neural systems of ctenophores are vastly different from those of other animals and use different sets of cellular and genetic mechanisms. Thus, neural systems appear to have at least two independent origins regardless of whether ctenophores or sponges are the earliest branching extant animal lineage.

Complete mitochondrial genome of four pheretimoid earthworms (Clitellata: Oligochaeta) and their phylogenetic reconstruction

Among oligochaetes, the Pheretima complex within the Megascolecidae is a major earthworm group. Recently, however, the systematics of the Pheretima complex based on morphology are challenged by molecular studies. Since little comparative analysis of earthworm complete mitochondrial genomes has been reported yet, we sequenced mitogenomes of four pheretimoid earthworm species to explore their phylogenetic relationships. The general earthworm genomic features are also found in four earthworms: all genes transcribed from the same strand, the same initiation codon ATG for each PCGs, and conserved structures of RNA genes. Interestingly we find an extra potential tRNA-leucine (CUN) in Amynthas longisiphonus. The earthworm mitochondrial ATP8 exhibits the highest evolutionary rate, while the gene CO1 evolves slowest. Phylogenetic analysis based on protein-coding genes (PCGs) strongly supports the monophyly of the Clitellata, Hirudinea, Oligochaeta, Megascolecidae and Pheretima complex. Our analysis, however, reveals non-monophyly within the genara Amynthas and Metaphire. Thus the generic divisions based on morphology in the Pheretima complex should be reconsidered.

The relevance of the storage of subunit c of ATP synthase in different forms and models of Batten disease (NCLs)

The discoveries of specific protein storage in the NCLs, particularly of subunit c of ATP synthase in most, and the sphingolipid activator proteins, SAPs or saposins A and D in CLN1, CLN10 and an unassigned form are reviewed. The subunit c stored in the relevant NCLs is the complete mature molecule including an unusual modification found only in animal species, trimethylation of its lysine-43. Because of its strongly hydrophobic and lipid-like properties subunit c is easily overlooked or incorrectly described. This is becoming more of a problem as subunit c is not detected in standard proteomic investigations. Methods are reviewed that allow its unequivocal characterisation. Subunit c storage and cellular storage body accumulation do not cause the neuropathology characteristic of these diseases. The function of the trimethyl group on lysine-43 of subunit c is considered, along with some indications of where its normal turnover may be disrupted in the NCLs.

Employing Phylogenomics to Resolve the Relationships among Cnidarians, Ctenophores, Sponges, Placozoans, and Bilaterians

Despite an explosion in the amount of sequence data, phylogenomics has failed to settle controversy regarding some critical nodes on the animal tree of life. Understanding relationships among Bilateria, Ctenophora, Cnidaria, Placozoa, and Porifera is essential for studying how complex traits such as neurons, muscles, and gastrulation have evolved. Recent studies have cast doubt on the historical viewpoint that sponges are sister to all other animal lineages with recent studies recovering ctenophores as sister. However, the ctenophore-sister hypothesis has been criticized as unrealistic and caused by systematic error. We review past phylogenomic studies and potential causes of systematic error in an effort to identify areas that can be improved in future studies. Increased sampling of taxa, less missing data, and a priori removal of sequences and taxa that may cause systematic error in phylogenomic inference will likely be the most fruitful areas of focus when assembling future datasets. Ultimately, we foresee metazoan relationships being resolved with higher support in the near future, and we caution against dismissing novel hypotheses merely because they conflict with historical viewpoints of animal evolution.

Error, signal, and the placement of Ctenophora sister to all other animals

Elucidating relationships among early animal lineages has been difficult, and recent phylogenomic analyses place Ctenophora sister to all other extant animals, contrary to the traditional view of Porifera as the earliest-branching animal lineage. To date, phylogenetic support for either ctenophores or sponges as sister to other animals has been limited and inconsistent among studies. Lack of agreement among phylogenomic analyses using different data and methods obscures how complex traits, such as epithelia, neurons, and muscles evolved. A consensus view of animal evolution will not be accepted until datasets and methods converge on a single hypothesis of early metazoan relationships and putative sources of systematic error (e.g., long-branch attraction, compositional bias, poor model choice) are assessed. Here, we investigate possible causes of systematic error by expanding taxon sampling with eight novel transcriptomes, strictly enforcing orthology inference criteria, and progressively examining potential causes of systematic error while using both maximum-likelihood with robust data partitioning and Bayesian inference with a site-heterogeneous model. We identified ribosomal protein genes as possessing a conflicting signal compared with other genes, which caused some past studies to infer ctenophores and cnidarians as sister. Importantly, biases resulting from elevated compositional heterogeneity or elevated substitution rates are ruled out. Placement of ctenophores as sister to all other animals, and sponge monophyly, are strongly supported under multiple analyses, herein.

Conservation of complete trimethylation of lysine-43 in the rotor ring of c-subunits of metazoan adenosine triphosphate (ATP) synthases

The rotors of ATP synthases turn about 100 times every second. One essential component of the rotor is a ring of hydrophobic c-subunits in the membrane domain of the enzyme. The rotation of these c-rings is driven by a transmembrane proton-motive force, and they turn against a surface provided by another membrane protein, known as subunit a. Together, the rotating c-ring and the static subunit a provide a pathway for protons through the membrane in which the c-ring and subunit a are embedded. Vertebrate and invertebrate c-subunits are well conserved. In the structure of the bovine F1-ATPase-c-ring subcomplex, the 75 amino acid c-subunit is folded into two transmembrane α-helices linked by a short loop. Each bovine rotor-ring consists of eight c-subunits with the N- and C-terminal α-helices forming concentric inner and outer rings, with the loop regions exposed to the phospholipid head-group region on the matrix side of the inner membrane. Lysine-43 is in the loop region and its ε-amino group is completely trimethylated. The role of this modification is unknown. If the trimethylated lysine-43 plays some important role in the functioning, assembly or degradation of the c-ring, it would be expected to persist throughout vertebrates and possibly invertebrates also. Therefore, we have carried out a proteomic analysis of c-subunits across representative species from different classes of vertebrates and from invertebrate phyla. In the twenty-nine metazoan species that have been examined, the complete methylation of lysine-43 is conserved, and it is likely to be conserved throughout the more than two million extant metazoan species. In unicellular eukaryotes and prokaryotes, when the lysine is conserved it is unmethylated, and the stoichiometries of c-subunits vary from 9-15. One possible role for the trimethylated residue is to provide a site for the specific binding of cardiolipin, an essential component of ATP synthases in mitochondria.