The mitochondrial genome of Hydra oligactis (Cnidaria, Hydrozoa) sheds new light on animal mtDNA evolution and cnidarian phylogeny

One Ring does not rule them all: Linear mtDNA in Metazoa

Recent advances in genome sequencing technologies have facilitated the exploration of the architecture of genomes, including mitochondrial genomes (mtDNA). In particular, whole genome sequencing has provided easier access to mitochondrial genomes with unusual organizations, which were difficult to obtain using traditional PCR-based approaches. As a consequence, there has been a steep increase in complete mtDNA sequences, particularly for Metazoa. The popular view of metazoan mtDNA is that of a small gene-dense circular chromosome. This view clashes with discoveries of a number of linear mtDNAs, particularly in non-bilaterian animals. Here, we review the distribution of linear mtDNA in Metazoa, namely in isopods, cnidarians, and sponges. We discuss the multiple origins of linear mitogenomes in these clades, where linearity has been linked to the likely insertion of a linear plasmid in cnidarians and the demosponge Acanthella acuta, while fixation of a heteroplasmy in the anticodon site of a tRNA might be responsible for the monolinear form of the mtDNA in some isopods. We also summarize our current knowledge of mechanisms that maintain the integrity of linear mitochromosomes, where a recurrent theme is the presence of terminal repeats that likely play the role of telomeres. We caution in defining a linear chromosome as complete, particularly when coding sequences and key features of linear DNA are missing. Finally, we encourage authors interested in mitogenome science to utilize all available data for linear mtDNA, including those tagged as "incomplete" or "unverified" in public databases, as they can still provide useful information such as phylogenetic characters and gene order.

Mitochondrial DNA of the Demosponge Acanthella acuta: Linear Architecture and Other Unique Features

While Acanthella acuta Schmidt 1862, a common demosponge found in the Mediterranean Sea and Atlantic Ocean, is morphologically similar to other sponges, its mitochondrial DNA (mtDNA) is unique within the class. In contrast to all other studied demosponges, the mtDNA of A. acuta is inferred to be linear and displays several unusual features such as inverted terminal repeats, group II introns in three mitochondrial genes, and two unique open reading frames (ORFs): one of which (ORF1535) combines a DNA polymerase domain with a DNA-directed RNA polymerase domain, while the second bears no discernible similarity to any reported sequences. The group II intron within the cox2 gene is the first such intron reported in an animal. Our phylogenetic analyses indicate that the cox1 intron is related to similar introns found in other demosponges, while the cox2 intron is likely not of animal origin. The two domains found within ORF1535 do not share a common origin and, along with the cox2 intron, were likely acquired by horizontal gene transfer. The findings of this paper open new avenues of exploration in the understanding of mtDNA linearization within Metazoa.

The complete mitochondrial genome and phylogenetic analysis of hydrozoan jellyfish Eirene ceylonensis (Cnidaria, Hydrozoa, Eirenidae) in the coastal sea of Qinhuangdao, China

Eirene ceylonensis, a hydrozoan jellyfish species with a complex polymorphic life cycle, is widely distributed in the Chinese coastal sea. In this study, we conducted sequencing and analysis of the first complete mitochondrial genome of E. ceylonensis, obtained from the coastal sea of Qinhuangdao, China. The linear mitochondrial genome is 14,997 bp in length with the overall AT content being 72.8%, encoding 13 protein-coding genes (PCGs), two transfer RNA (tRNA) genes (tRNA-Met and tRNA-Trp) and two ribosomal RNA (rRNA) genes (rrnS and rrnL). Phylogenetic analysis of 13 PCGs suggests that the E. ceylonensis is closely related to Laomedea flexuosa. The availability of the complete mitochondrial genome of E. ceylonensis will be useful for studying the evolutionary relationships of hydrozoan jellyfish species.

Revisiting mitogenome evolution in Medusozoa with eight new mitochondrial genomes

Mitogenomics has improved our understanding of medusozoan phylogeny. However, sequenced medusozoan mitogenomes remain scarce, and Medusozoa phylogeny studies often analyze mitogenomic sequences without incorporating mitogenome rearrangements. To better understand medusozoan evolution, we analyzed Medusozoa mitogenome phylogeny by sequencing and assembling eight mitogenomes from three classes (Cubozoa, Hydrozoa, and Scyphozoa). We reconstructed the mitogenome phylogeny using these mitogenomes and 84 other existing cnidarian mitogenomes to study mitochondrial gene rearrangements. All reconstructed mitogenomes had 13 mitochondrial protein-coding genes and two ribosomal genes typical for Medusozoa. Non-cubozoan mitogenomes were all linear and had typical gene orders, while arrangement of genes in the fragmented Cubozoa (Morbakka sp.) mitogenome differed from other Cubozoa mitogenomes. Gene order comparisons and ancestral state reconstruction suggest minimal rearrangements within medusozoan classes except for Hydrozoa. Our findings support a staurozoan ancestral medusozoan gene order, expand the pool of available medusozoan mitogenomes, and enhance our understanding of medusozoan phylogenetic relationships.

Evolution of mitochondrial and nuclear genomes in Pennatulacea

We examine the phylogeny of sea pens using sequences of whole mitochondrial genomes and the nuclear ribosomal cluster generated through low coverage Illumina sequencing. Taxon sampling includes 30 species in 19 genera representing 13 families. Ancestral state reconstruction shows that most sea pen mitochondrial genomes have the ancestral gene order, and that Pennatulacea with diverse gene orders are found in a single clade. The monophyly of Pennatulidae and Protoptilidae are rejected by both the mitochondrial and nuclear dataset, while the mitochondrial dataset further rejects monophyly of Virgulariidae, and the nuclear dataset rejects monophyly of Kophobelemnidae. We show discordance between nuclear ribosomal gene cluster phylogenies and whole mitochondrial genome phylogenies and highlight key Pennatulacea taxa that could be included in cnidarian genome-wide studies to better resolve the sea pen tree of life. We further illustrate how well frequently sequenced markers capture the overall diversity of the mitochondrial genome and the nuclear ribosomal genes in sea pens.

The complete mitogenome sequence of the coral lily (Lilium pumilum) and the Lanzhou lily (Lilium davidii) in China

Background

The mitogenomes of higher plants are conserved. This study was performed to complete the mitogenome of two China Lilium species (Lilium pumilum Redouté and Lilium davidii var. unicolor (Hoog) cotton).

Methods

Genomic DNA was separately extracted from the leaves of L. pumilum and L. davidii in triplicate and used for sequencing. The mitogenome of Allium cepa was used as a reference. Genome assembly, annotation and phylogenetic tree were analyzed.

Results

The mitogenome of L. pumilum and L. davidii was 988,986 bp and 924,401 bp in length, respectively. There were 22 core protein-coding genes (including atp1, atp4, atp6, atp9, ccmB, ccmC, ccmFc, ccmFN1, ccmFN2, cob, cox3, matR, mttB, nad1, nad2, nad3, nad4, nad4L, nad5, nad6, nad7 and nad9), one open reading frame and one ribosomal protein-coding gene (rps12) in the mitogenomes. Compared with the A. cepa mitogenome, the coding sequence of the 24 genes and intergenic spacers in L. pumilum and L. davidii mitogenome contained 1,621 and 1,617 variable sites, respectively. In the phylogenetic tree, L. pumilum and L. davidii were distinct from A. cepa (NC_030100).

Conclusions

L. pumilum and L. davidii mitogenomes have far distances from other plants. This study provided additional information on the species resources of China Lilium.

The linear mitochondrial genome of commensal hydroid Eutima japonica (Cnidaria, Hydrozoa, Eirenidae)

Here, we present the whole mitochondrial genome of commensal hydroid Eutima japonica McCrady 1859 (family Eirinidae); this is the first specimen of the family to have its mitogenome sequenced. The linear mitogenome is 15,315 bp in length and consists of 13 protein-coding genes (PCGs), large and small ribosomal subunits (rRNA), methionine and tryptophan transfer RNA (tRNA) genes (trnM and trnW), and a partial copy of cytochrome oxidase subunit I (cox1) pseudogene, as is typical for the class Hydrozoa. Nucleotide sequences of two cox1 genes at two ends of the linear mitogenome form a part of inverted terminal repeat. The overall genomic structure and gene arrangement of 13 PCGs were identical to the reported mitochondrial genomes of hydrozoans, except for the positions of two tRNA genes. Phylogenetic analysis of E. japonica 13 PCGs and other cnidarians recovers a closest relationship with the derived cluster of two hydrozoans, Laomedea flexuosa and Obelia longissimi within Leptothecata.

The venom proteome of three common scyphozoan jellyfishes (Chrysaora caliparea, Cyanea nozakii and Lychnorhiza malayensis) (Cnidaria: Scyphozoa) from the coastal waters of India

The jellyfish venom stored in nematocysts contains highly toxic compounds comprising of polypeptides, enzymes and other proteins, which form their chemical defence armoury against predators. We have characterized the proteome of crude venom extract from three bloom-forming scyphozoan jellyfish along the south-west coast of India, Chrysaora caliparea, Cyanea nozakii and Lychnorhiza malayensis using a Quadrupole-Time of Flight (Q/TOF) mass spectrometry analysis. The most abundant toxin identified from Chrysaora caliparea and Lychnorhiza malayensis is similar to the pore-forming toxins and metalloproteinases. A protective antioxidant enzyme called peroxiredoxin was found abundantly in Cyanea nozakii. Metalloproteinase identified from the C. caliparea shows similarity with the venom of pit viper (Bothrops pauloensis), while that of L. malayensis was similar to the venom of snakes such as the Bothrops insularis and Bothrops asper. Kininogen-1 is a secreted protein, identified for the first time from the jellyfish L. malayensis. The proteome analysis of Cyanea nozakii, Chrysaora caliparea and Lychnorhiza malayensis contained 20, 12, 8 unique proteins, respectively. Our study characterized the proteome map of crude venom extract from L. malayensis and C. caliparea for the first time, and the venom profile is compared with published information elsewhere. Proteomic data from this study has been made available in the public domain.

The challenge of Coccidae (Hemiptera: Coccoidea) mitochondrial genomes: The case of Saissetia coffeae with novel truncated tRNAs and gene rearrangements

There have been few reports of complete mitochondrial genomes (mitogenomes) of scale insects, and it has been indicated that complex and novel structures in their mitogenomes may lead to difficulties in sequencing, assembly and annotation. Transfer RNAs (tRNAs) usually possess typical cloverleaf secondary structures, and truncated tRNAs are rarely found in insect mitogenomes. Here, we report a complete Saissetia coffeae mitogenome (15,389 bp) with high A+T content (84.7%) sequenced by next-generation sequencing (NGS) methods. Genes in the mitogenome were annotated, and nine tRNAs were not found using MITOS. Most of the detected tRNAs were significantly truncated without the dihydrouridine (DHU) arm or the TΨC (T) arm. In addition, the 9 "lost" tRNAs containing mismatched base pairs were retrieved based on the tRNA annotation workflow for Coccidae described in our study. The gene arrangement in the Saissetia coffeae mitogenome was significantly different from that in other hemipteran insects. Additionally, Bayesian and maximum likelihood trees based on the mitochondrial genes showed a long branch of the Saissetia lineage, indicating significant nonsynonymous substitutions or high evolutionary rates in the Saissetia lineage. We provide a reference mitogenome for the assembly and annotation of the Coccidae mitogenome and offer insights into the evolution of scale insects.

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.

Complete mitochondrial genome sequence of Aequorea coerulescens

The complete mitochondrial genome sequences of hydromedusa Aequorea coerulescens, a dominant giant jellyfish distributing in the Yellow Sea and northern East China Sea, China, is first reported in this research. Its mitochondrial DNA has 14,804 bp in length with a linear structure, containing 13 protein-coding genes (PCGs), 2rRNA (s-rRNA and l-rRNA), and 2tRNA (trnaW-TGA and tranM-AGT). The A + T content of the whole base composition of the genome is 72.8% (A: 30.15%; C: 12.33%; G: 14.87%; T: 42.64%). ATG, ATA, and ATT are start codons in five (ATP8, ATP6, COX3, NAD6, NAD4l), four (NAD2, NAD5, NAD1, COX1) and two PCGs (NAD3, NAD4), respectively. COB and COX2 began with GTG and CAA as start codon, respectively. TAA and ATA were the stop codon of ATP6 and NAD5 as well as NAD2 and NAD4, respectively. However, other PCGs were terminated with different stop codons. The NJ phylogenetic tree among the related 15 jellyfish species showed that A. coerulescens is close to Hydra spp.

Novel diversity in mitochondrial genomes of deep-sea Pennatulacea (Cnidaria: Anthozoa: Octocorallia)

We present the first documented complete mitogenomes of deep-sea Pennatulacea, representing nine genera and eight families. These include one species each of the deep-sea genera Funiculina, Halipteris, Protoptilum and Distichoptilum, four species each of Umbellula and Pennatula, three species of Kophobelemnon and two species of Anthoptilum, as well as one species of the epi- and mesobenthic genus Virgularia. Seventeen circular genomes ranged from 18,513 bp (Halipteris cf. finmarchica) to 19,171 bp (Distichoptilum gracile) and contained all genes standard to octocoral mitochondrial genomes (14 protein-coding genes, two ribosomal RNA genes and one transfer RNA). We found at least three different gene orders in Pennatulacea: the ancestral gene order, the gene order found in bamboo corals (Family Isididae), and a novel gene order. The mitogenome of one species of Umbellula has a bipartite genome (∼13 kbp and ∼5 kbp), with good evidence that both parts are circular.

Molecular and morphological congruence of three new cryptic Neopetrosia spp. in the Caribbean

Neopetrosia proxima (Porifera: Demospongiae: Haplosclerida) is described as a morphologically variable sponge common on shallow reefs of the Caribbean. However, the range of morphological and reproductive variation within putative N. proxima led us to hypothesize that such variability may be indicative of cryptic species rather than plasticity. Using DNA sequences and morphological characters we confirmed the presence of three previously undescribed species of Neopetrosia. Morphological differences of each new congener were best resolved by partial gene sequences of the mitochondrial cytochrome oxidase subunit 1 over nuclear ones (18S rRNA and 28S rRNA). Several new characters for Neopetrosia were revealed by each new species. For example, N. dendrocrevacea sp. nov. and N. cristata sp. nov. showed the presence of grooves on the surface of the sponge body that converge at the oscula, and a more disorganized skeleton than previously defined for the genus. N. sigmafera sp. nov. adds the (1) presence of sigma microscleres, (2) significantly wider/longer oxeas (>200 μm), and (3) the presence of parenchymella larvae. Sampling of conspecifics throughout several locations in the Caribbean revealed larger spicules in habitats closer to the continental shelf than those in remote island locations. Our study highlights the importance of integrating molecular and morphological systematics for the discrimination of new Neopetrosia spp. despite belonging to one of several polyphyletic groups (families, genera) within the current definition of the order Haplosclerida.

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.

Complete mitochondrial genome and evolutionary analysis of Turritopsis dohrnii, the "immortal" jellyfish with a reversible life-cycle

Turritopsis dohrnii (Cnidaria, Hydrozoa, Hydroidolina, Anthoathecata) is the only known metazoan that is capable of reversing its life cycle via morph rejuvenation from the adult medusa stage to the juvenile polyp stage. Here, we present a complete mitochondrial (mt) genome sequence of T. dohrnii, which harbors genes for 13 proteins, two transfer RNAs, and two ribosomal RNAs. The T. dohrnii mt genome is characterized by typical features of species in the Hydroidolina subclass, such as a high A+T content (71.5%), reversed transcriptional orientation for the large rRNA subunit gene, and paucity of CGN codons. An incomplete complementary duplicate of the cox1 gene was found at the 5' end of the T. dohrnii mt chromosome, as were variable repeat regions flanking the chromosome. We identified species-specific variations (nad5, nad6, cob, and cox1 genes) and putative selective constraints (atp8, nad1, nad2, and nad5 genes) in the mt genes of T. dohrnii, and predicted alterations in tertiary structures of respiratory chain proteins (NADH4, NADH5, and COX1 proteins) of T. dohrnii. Based on comparative analyses of available hydrozoan mt genomes, we also determined the taxonomic relationships of T. dohrnii, recovering Filifera IV as a paraphyletic taxon, and assessed intraspecific diversity of various Hydrozoa species.

The Complete Mitochondrial Genome of Aleurocanthus camelliae: Insights into Gene Arrangement and Genome Organization within the Family Aleyrodidae

There are numerous gene rearrangements and transfer RNA gene absences existing in mitochondrial (mt) genomes of Aleyrodidae species. To understand how mt genomes evolved in the family Aleyrodidae, we have sequenced the complete mt genome of Aleurocanthus camelliae and comparatively analyzed all reported whitefly mt genomes. The mt genome of A. camelliae is 15,188 bp long, and consists of 13 protein-coding genes, two rRNA genes, 21 tRNA genes and a putative control region (GenBank: KU761949). The tRNA gene, trnI, has not been observed in this genome. The mt genome has a unique gene order and shares most gene boundaries with Tetraleurodes acaciae. Nineteen of 21 tRNA genes have the conventional cloverleaf shaped secondary structure and two (trnS₁ and trnS₂) lack the dihydrouridine (DHU) arm. Using ARWEN and homologous sequence alignment, we have identified five tRNA genes and revised the annotation for three whitefly mt genomes. This result suggests that most absent genes exist in the genomes and have not been identified, due to be lack of technology and inference sequence. The phylogenetic relationships among 11 whiteflies and Drosophila melanogaster were inferred by maximum likelihood and Bayesian inference methods. Aleurocanthus camelliae and T. acaciae form a sister group, and all three Bemisia tabaci and two Bemisia afer strains gather together. These results are identical to the relationships inferred from gene order. We inferred that gene rearrangement plays an important role in the mt genome evolved from whiteflies.

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.

Two-step evolution of endosymbiosis between hydra and algae

In the Hydra vulgaris group, only 2 of the 25 strains in the collection of the National Institute of Genetics in Japan currently show endosymbiosis with green algae. However, whether the other non-symbiotic strains also have the potential to harbor algae remains unknown. The endosymbiotic potential of non-symbiotic strains that can harbor algae may have been acquired before or during divergence of the strains. With the aim of understanding the evolutionary process of endosymbiosis in the H. vulgaris group, we examined the endosymbiotic potential of non-symbiotic strains of the H. vulgaris group by artificially introducing endosymbiotic algae. We found that 12 of the 23 non-symbiotic strains were able to harbor the algae until reaching the grand-offspring through the asexual reproduction by budding. Moreover, a phylogenetic analysis of mitochondrial genome sequences showed that all the strains with endosymbiotic potential grouped into a single cluster (cluster γ). This cluster contained two strains (J7 and J10) that currently harbor algae; however, these strains were not the closest relatives. These results suggest that evolution of endosymbiosis occurred in two steps; first, endosymbiotic potential was gained once in the ancestor of the cluster γ lineage; second, strains J7 and J10 obtained algae independently after the divergence of the strains. By demonstrating the evolution of the endosymbiotic potential in non-symbiotic H. vulgaris group strains, we have clearly distinguished two evolutionary steps. The step-by-step evolutionary process provides significant insight into the evolution of endosymbiosis in cnidarians.

Cloning a Chymotrypsin-Like 1 (CTRL-1) Protease cDNA from the Jellyfish Nemopilema nomurai

An enzyme in a nematocyst extract of the Nemopilema nomurai jellyfish, caught off the coast of the Republic of Korea, catalyzed the cleavage of chymotrypsin substrate in an amidolytic kinetic assay, and this activity was inhibited by the serine protease inhibitor, phenylmethanesulfonyl fluoride. We isolated the full-length cDNA sequence of this enzyme, which contains 850 nucleotides, with an open reading frame of 801 encoding 266 amino acids. A blast analysis of the deduced amino acid sequence showed 41% identity with human chymotrypsin-like (CTRL) and the CTRL-1 precursor. Therefore, we designated this enzyme N. nomurai CTRL-1. The primary structure of N. nomurai CTRL-1 includes a leader peptide and a highly conserved catalytic triad of His(69), Asp(117), and Ser(216). The disulfide bonds of chymotrypsin and the substrate-binding sites are highly conserved compared with the CTRLs of other species, including mammalian species. Nemopilema nomurai CTRL-1 is evolutionarily more closely related to Actinopterygii than to Scyphozoan (Aurelia aurita) or Hydrozoan (Hydra vulgaris). The N. nomurai CTRL1 was amplified from the genomic DNA with PCR using specific primers designed based on the full-length cDNA, and then sequenced. The N. nomurai CTRL1 gene contains 2434 nucleotides and four distinct exons. The 5' donor splice (GT) and 3' acceptor splice sequences (AG) are wholly conserved. This is the first report of the CTRL1 gene and cDNA structures in the jellyfish N. nomurai.

Insights into the transcriptional and translational mechanisms of linear organellar chromosomes in the box jellyfish Alatina alata (Cnidaria: Medusozoa: Cubozoa)

BACKGROUND

In most animals, the mitochondrial genome is characterized by its small size, organization into a single circular molecule, and a relative conservation of the number of encoded genes. In box jellyfish (Cubozoa, Cnidaria), the mitochondrial genome is organized into 8 linear mito-chromosomes harboring between one and 4 genes each, including 2 extra protein-coding genes: mt-polB and orf314. Such an organization challenges the traditional view of mitochondrial DNA (mtDNA) expression in animals. In this study, we investigate the pattern of mitochondrial gene expression in the box jellyfish Alatina alata, as well as several key nuclear-encoded molecular pathways involved in the processing of mitochondrial gene transcription.

RESULTS

Read coverage of DNA-seq data is relatively uniform for all 8 mito-chromosomes, suggesting that each mito-chromosome is present in equimolar proportion in the mitochondrion. Comparison of DNA and RNA-seq based assemblies indicates that mito-chromosomes are transcribed into individual transcripts in which the beginning and ending are highly conserved. Expression levels for mt-polB and orf314 are similar to those of other mitochondrial-encoded genes, which provides further evidence for them having functional roles in the mitochondrion. Survey of the transcriptome suggests recognition of the mitochondrial tRNA-Met by the cytoplasmic aminoacyl-tRNA synthetase counterpart and C-to-U editing of the cytoplasmic tRNA-Trp after import into the mitochondrion. Moreover, several mitochondrial ribosomal proteins appear to be lost.

CONCLUSIONS

This study represents the first survey of mitochondrial gene expression of the linear multi-chromosomal mtDNA in box jellyfish (Cubozoa). Future exploration of small RNAs and the proteome of the mitochondrion will test the hypotheses presented herein.

Phylogenetic analysis of higher-level relationships within Hydroidolina (Cnidaria: Hydrozoa) using mitochondrial genome data and insight into their mitochondrial transcription

Hydrozoans display the most morphological diversity within the phylum Cnidaria. While recent molecular studies have provided some insights into their evolutionary history, sister group relationships remain mostly unresolved, particularly at mid-taxonomic levels. Specifically, within Hydroidolina, the most speciose hydrozoan subclass, the relationships and sometimes integrity of orders are highly unsettled. Here we obtained the near complete mitochondrial sequence of twenty-six hydroidolinan hydrozoan species from a range of sources (DNA and RNA-seq data, long-range PCR). Our analyses confirm previous inference of the evolution of mtDNA in Hydrozoa while introducing a novel genome organization. Using RNA-seq data, we propose a mechanism for the expression of mitochondrial mRNA in Hydroidolina that can be extrapolated to the other medusozoan taxa. Phylogenetic analyses using the full set of mitochondrial gene sequences provide some insights into the order-level relationships within Hydroidolina, including siphonophores as the first diverging clade, a well-supported clade comprised of Leptothecata-Filifera III-IV, and a second clade comprised of Aplanulata-Capitata s.s.-Filifera I-II. Finally, we describe our relatively inexpensive and accessible multiplexing strategy to sequence long-range PCR amplicons that can be adapted to most high-throughput sequencing platforms.

Discordance between morphological and molecular species boundaries among Caribbean species of the reef sponge Callyspongia

Sponges are among the most species-rich and ecologically important taxa on coral reefs, yet documenting their diversity is difficult due to the simplicity and plasticity of their morphological characters. Genetic attempts to identify species are hampered by the slow rate of mitochondrial sequence evolution characteristic of sponges and some other basal metazoans. Here we determine species boundaries of the Caribbean coral reef sponge genus Callyspongia using a multilocus, model-based approach. Based on sequence data from one mitochondrial (COI), one ribosomal (28S), and two single-copy nuclear protein-coding genes, we found evolutionarily distinct lineages were not concordant with current species designations in Callyspongia. While C. fallax,C. tenerrima, and C. plicifera were reciprocally monophyletic, four taxa with different morphologies (C. armigera,C. longissima,C. eschrichtii, and C. vaginalis) formed a monophyletic group and genetic distances among these taxa overlapped distances within them. A model-based method of species delimitation supported collapsing these four into a single evolutionary lineage. Variation in spicule size among these four taxa was partitioned geographically, not by current species designations, indicating that in Callyspongia, these key taxonomic characters are poor indicators of genetic differentiation. Taken together, our results suggest a complex relationship between morphology and species boundaries in sponges.

Octocoral mitochondrial genomes provide insights into the phylogenetic history of gene order rearrangements, order reversals, and cnidarian phylogenetics

We use full mitochondrial genomes to test the robustness of the phylogeny of the Octocorallia, to determine the evolutionary pathway for the five known mitochondrial gene rearrangements in octocorals, and to test the suitability of using mitochondrial genomes for higher taxonomic-level phylogenetic reconstructions. Our phylogeny supports three major divisions within the Octocorallia and show that Paragorgiidae is paraphyletic, with Sibogagorgia forming a sister branch to the Coralliidae. Furthermore, Sibogagorgia cauliflora has what is presumed to be the ancestral gene order in octocorals, but the presence of a pair of inverted repeat sequences suggest that this gene order was not conserved but rather evolved back to this apparent ancestral state. Based on this we recommend the resurrection of the family Sibogagorgiidae to fix the paraphyly of the Paragorgiidae. This is the first study to show that in the Octocorallia, mitochondrial gene orders have evolved back to an ancestral state after going through a gene rearrangement, with at least one of the gene orders evolving independently in different lineages. A number of studies have used gene boundaries to determine the type of mitochondrial gene arrangement present. However, our findings suggest that this method known as gene junction screening may miss evolutionary reversals. Additionally, substitution saturation analysis demonstrates that while whole mitochondrial genomes can be used effectively for phylogenetic analyses within Octocorallia, their utility at higher taxonomic levels within Cnidaria is inadequate. Therefore for phylogenetic reconstruction at taxonomic levels higher than subclass within the Cnidaria, nuclear genes will be required, even when whole mitochondrial genomes are available.

Lock, stock and two different barrels: comparing the genetic composition of morphotypes of the indo-pacific sponge Xestospongia testudinaria

The giant barrel sponge Xestospongiatestudinaria is an ecologically important species that is widely distributed across the Indo-Pacific. Little is known, however, about the precise biogeographic distribution and the amount of morphological and genetic variation in this species. Here we provide the first detailed, fine-scaled (<200 km(2)) study of the morphological and genetic composition of X. testudinaria around Lembeh Island, Indonesia. Two mitochondrial (CO1 and ATP6 genes) and one nuclear (ATP synthase β intron) DNA markers were used to assess genetic variation. We identified four distinct morphotypes of X. testudinaria around Lembeh Island. These morphotypes were genetically differentiated with both mitochondrial and nuclear markers. Our results indicate that giant barrel sponges around Lembeh Island, which were all morphologically identified as X. testudinaria, consist of at least two different lineages that appear to be reproductively isolated. The first lineage is represented by individuals with a digitate surface area, CO1 haplotype C5, and is most abundant around the harbor area of Bitung city. The second lineage is represented by individuals with a predominantly smooth surface area, CO1 haplotype C1 and can be found all around Lembeh Island, though to a lesser extent around the harbor of Bitung city. Our findings of two additional unique genetic lineages suggests the presence of an even broader species complex possibly containing more than two reproductively isolated species. The existence of X. testudinaria as a species complex is a surprising result given the size, abundance and conspicuousness of the sponge.

Mitogenomics at the base of Metazoa

Unraveling the base of metazoan evolution is of crucial importance for rooting the metazoan Tree of Life. This subject has attracted substantial attention for more than a century and recently fueled a burst of modern phylogenetic studies. Conflicting scenarios from different studies and incongruent results from nuclear versus mitochondrial markers challenge current molecular phylogenetic approaches. Here we analyze the presently most comprehensive data sets of mitochondrial genomes from non-bilaterian animals to illuminate the phylogenetic relationships among early branching metazoan phyla. The results of our analyses illustrate the value of mitogenomics and support previously known topologies between animal phyla but also identify several problematic taxa, which are sensitive to long branch artifacts or missing data.

The complete mitochondrial genome of Hydra vulgaris (Hydroida: Hydridae)

The complete mitochondrial genome of Hydra vulgaris (Hydroida: Hydridae) is composed of two linear DNA molecules. The mitochondrial DNA (mtDNA) molecule 1 is 8010 bp long and contains six protein-coding genes, large subunit rRNA, methionine and tryptophan tRNAs, two pseudogenes consisting respectively of a partial copy of COI, and terminal sequences at two ends of the linear mtDNA, while the mtDNA molecule 2 is 7576 bp long and contains seven protein-coding genes, small subunit rRNA, methionine tRNA, a pseudogene consisting of a partial copy of COI and terminal sequences at two ends of the linear mtDNA. COI gene begins with GTG as start codon, whereas other 12 protein-coding genes start with a typical ATG initiation codon. In addition, all protein-coding genes are terminated with TAA as stop codon.

Exceptionally high cumulative percentage of NUMTs originating from linear mitochondrial DNA molecules in the Hydra magnipapillata genome

Background

In contrast to most animal genomes, mitochondrial genomes in species belonging to the phylum Cnidaria show distinct variations in genome structure, including the mtDNA structure (linear or circular) and the presence or absence of introns in protein-coding genes. Therefore, the analysis of nuclear insertions of mitochondrial sequences (NUMTs) in cnidarians allows us to compare the NUMT content in animals with different mitochondrial genome structures.

Results

NUMT identification in the Hydra magnipapillata, Nematostella vectensis and Acropora digitifera genomes showed that the NUMT density in the H. magnipapillata genome clearly exceeds that in other two cnidarians with circular mitochondrial genomes. We found that H. magnipapillata is an exceptional ancestral metazoan with a high NUMT cumulative percentage but a large genome, and its mitochondrial genome linearisation might be responsible for the NUMT enrichment. We also detected the co-transposition of exonic and intronic fragments within NUMTs in N. vectensis and provided direct evidence that mitochondrial sequences can be transposed into the nuclear genome through DNA-mediated fragment transfer. In addition, NUMT expression analyses showed that NUMTs are co-expressed with adjacent protein-coding genes, suggesting the relevance of their biological function.

Conclusions

Taken together, our results provide valuable information for understanding the impact of mitochondrial genome structure on the interaction of mitochondrial molecules and nuclear genomes.

Characterization and localization of mitochondrial DNA-encoded tRNAs and nuclear DNA-encoded tRNAs in the sea anemone Metridium senile

The mitochondrial (mt) genome of the sea anemone Metridium senile contains genes for only two transfer RNAs (tRNAs), tRNAf-Met and tRNATrp. Experiments were conducted to seek evidence for the occurrence of functional tRNAs corresponding to these genes and for the participation of nuclear DNA-encoded tRNAs in mt-protein synthesis. RNA sequences corresponding to the two mt-tRNA genes were located in mitochondria and it was shown that 3'-CC (and possibly A, but no other nucleotide) is added post-transcriptionally to the 3' end of at least 50 % of mt-tRNAf-Met molecules and to a small fraction of the mt-tRNATrp molecules. Using specific oligonucleotide primers based on expected nuclear DNA-encoded tRNAs in a series of RACE experiments, we located the nuclear genes for tRNAGln, tRNAIle, tRNAi-Met, tRNAVal and tRNAThr. Data from Northern blot analyses indicated that mtDNA-encoded tRNAf-Met is limited to mitochondria but that nuclear DNA-encoded tRNAVal and tRNAi-Met are present in the cytoplasm and in mitochondria. These data provide direct evidence that in M. senile, mature, functional tRNAs are transcribed from the mtDNA-encoded tRNAf-Met and tRNATrp genes, and are consistent with the interpretation that both nuclear DNA-encoded tRNAVal and tRNAi-Met are utilized in mitochondrial and cytosolic protein synthesis.

Complete mitochondrial genomes of the Japanese pink coral (Corallium elatius) and the Mediterranean red coral (Corallium rubrum): a reevaluation of the phylogeny of the family Coralliidae based on molecular data

Precious corals are soft corals belonging to the family Coralliidae (Anthozoa: Octocorallia: Alcyonacea) and class Anthozoa, whose skeletal axes are used for jewelry. The family Coralliidae includes ca. 40 species and was originally thought to comprise of the single genus Corallium. In 2003, Corallium was split into two genera, Corallium and Paracorallium, and seven species were moved to this newly identified genus on the bases of morphological features. Previously, we determined the complete mitochondrial genome sequence of two precious corals Paracorallium japonicum and Corallium konojoi, in order to clarify their systematic positions. The two genomes showed high nucleotide sequence identity, but their gene order arrangements were not identical. Here, we determined three complete mitochondrial genome sequences from the one specimen of Mediterranean Corallium rubrum and two specimens of Corallium elatius coming from Kagoshima (South Japan). The circular mitochondrial genomes of C. rubrum and C. elatius are 18,915bp and 18,969-18,970bp in length, respectively, and encode 14 typical octocorallian protein-coding genes (nad1-6, nad4L, cox1-3, cob, atp6, atp8, and mtMutS, which is an octocoral-specific mismatch repair gene homologue), two ribosomal RNA genes (rns and rnl), and one transfer RNA (trnM). The overall nucleotide differences between C. konojoi and each C. elatius haplotype (T2007 and I2011) are only 10 and 11 nucleotides, respectively; this degree of similarity indicates that C. elatius and C. konojoi are very closely related species. Notably, the C. rubrum mitochondrial genome shows more nucleotide sequence identity to P. japonicum (99.5%) than to its congeneric species C. konojoi (95.3%) and C. elatius (95.3%). Moreover, the gene order arrangement of C. rubrum was the same as that of P. japonicum, while that of C. elatius was the same as C. konojoi. Phylogenetic analysis based on three mitochondrial genes from 24 scleraxonian species shows that the family Coralliidae is separated into two distinct groups, recovering Corallium as a paraphyletic genus. Our results indicate that the currently accepted generic classification of Coralliidae should be reconsidered.

The complete mitochondrial genome of Chinese green hydra, Hydra sinensis (Hydroida: Hydridae)

The complete mitochondrial genome of Chinese green hydra, Hydra sinensis (Hydroida: Hydridae) is a linear molecule of 16,189 bp in length, containing 13 protein-coding genes, small and large subunit ribosomal RNAs, methionine and tryptophan transfer RNAs, a pseudogene consisting of a partial copy of COI and terminal sequences at two ends of the linear mitochondrial DNA. The A + T content of the overall base composition of H-strand is 77.2% (T: 41.7%; C: 10.9%; A: 35.5%; and G: 11.9%). COI and ND1 genes begin with GTG as start codon, while other 11 protein-coding genes start with a typical ATG initiation codon. COII, ATP8, ATP6, COIII, ND5, ND6, ND3, ND1, ND4 and COI genes are terminated with TAA as stop codon, ND4L ends with TAG, ND2 ends with TA and Cyt b ends with T.

Cnidarian phylogenetic relationships as revealed by mitogenomics

BACKGROUND

Cnidaria (corals, sea anemones, hydroids, jellyfish) is a phylum of relatively simple aquatic animals characterized by the presence of the cnidocyst: a cell containing a giant capsular organelle with an eversible tubule (cnida). Species within Cnidaria have life cycles that involve one or both of the two distinct body forms, a typically benthic polyp, which may or may not be colonial, and a typically pelagic mostly solitary medusa. The currently accepted taxonomic scheme subdivides Cnidaria into two main assemblages: Anthozoa (Hexacorallia + Octocorallia) - cnidarians with a reproductive polyp and the absence of a medusa stage - and Medusozoa (Cubozoa, Hydrozoa, Scyphozoa, Staurozoa) - cnidarians that usually possess a reproductive medusa stage. Hypothesized relationships among these taxa greatly impact interpretations of cnidarian character evolution.

RESULTS

We expanded the sampling of cnidarian mitochondrial genomes, particularly from Medusozoa, to reevaluate phylogenetic relationships within Cnidaria. Our phylogenetic analyses based on a mitochogenomic dataset support many prior hypotheses, including monophyly of Hexacorallia, Octocorallia, Medusozoa, Cubozoa, Staurozoa, Hydrozoa, Carybdeida, Chirodropida, and Hydroidolina, but reject the monophyly of Anthozoa, indicating that the Octocorallia + Medusozoa relationship is not the result of sampling bias, as proposed earlier. Further, our analyses contradict Scyphozoa [Discomedusae + Coronatae], Acraspeda [Cubozoa + Scyphozoa], as well as the hypothesis that Staurozoa is the sister group to all the other medusozoans.

CONCLUSIONS

Cnidarian mitochondrial genomic data contain phylogenetic signal informative for understanding the evolutionary history of this phylum. Mitogenome-based phylogenies, which reject the monophyly of Anthozoa, provide further evidence for the polyp-first hypothesis. By rejecting the traditional Acraspeda and Scyphozoa hypotheses, these analyses suggest that the shared morphological characters in these groups are plesiomorphies, originated in the branch leading to Medusozoa. The expansion of mitogenomic data along with improvements in phylogenetic inference methods and use of additional nuclear markers will further enhance our understanding of the phylogenetic relationships and character evolution within Cnidaria.

Molecular and morphological data support reclassification of the octocoral genus Isidoides

The rare octocoral genus Isidoides Nutting, 1910 was originally placed in the Gorgonellidae (now the Ellisellidae), even though it showed a remarkable similarity to the Isidae (now the Isididae). Isidoides was not classified in the Isididae mostly because the type specimen lacked skeletal nodes, a defining characteristic of that family. The genus was later assigned to the Chrysogorgiidae based on sclerite morphology. Specimens were recently collected in the south-western Pacific, providing material for genetic analysis and detailed characterisation of the morphology, and allowing us to consider the systematic placement of this taxon within the suborder Calcaxonia. A previously reported phylogeny allowed us to reject monophyly with the Chrysogorgiidae, and infer a close relationship with the Isididae subfamily Keratoisidinae. While scanning for molecular variation across mitochondrial genes, we discovered a novel gene order that is, based on available data, unique among metazoans. Despite these new data, the systematic placement of Isidoides is still unclear, as (1) the phylogenetic relationships among Isididae subfamilies remain poorly resolved, (2) genetic distances between mitochondrial mtMutS sequences from Isidoides and Keratoisidinae are characteristic of intra-familial distances, and (3) mitochondrial gene rearrangements may occur among confamilial genera. For these reasons, and because a revision of the Isididae is beyond the scope of this contribution, we amend the familial placement of Isidoides to incertae sedis.

Phylogenetic placement of Hydra and relationships within Aplanulata (Cnidaria: Hydrozoa)

The model organism Hydra belongs to the hydrozoan clade Aplanulata. Despite being a popular model system for development, little is known about the phylogenetic placement of this taxon or the relationships of its closest relatives. Previous studies have been conflicting regarding sister group relationships and have been unable to resolve deep nodes within the clade. In addition, there are several putative Aplanulata taxa that have never been sampled for molecular data or analyzed using multiple markers. Here, we combine the fast-evolving cytochrome oxidase 1 (CO1) mitochondrial marker with mitochondrial 16S, nuclear small ribosomal subunit (18S, SSU) and large ribosomal subunit (28S, LSU) sequences to examine relationships within the clade Aplanulata. We further discuss the relative contribution of four different molecular markers to resolving phylogenetic relationships within Aplanulata. Lastly, we report morphological synapomorphies for some of the major Aplanulata genera and families, and suggest new taxonomic classifications for two species of Aplanulata, Fukaurahydra anthoformis and Corymorpha intermedia, based on a preponderance of molecular and morphological data that justify the designation of these species to different genera.

Mitochondrial genome of the freshwater jellyfish Craspedacusta sowerbyi and phylogenetics of Medusozoa

The 17,922 base pairs (bp) nucleotide sequence of the linear mitochondrial DNA (mtDNA) molecule of the freshwater jellyfish Craspedacusta sowerbyi (Hydrozoa, Trachylina, Limnomedusae) has been determined. This sequence exhibits surprisingly low A+T content (57.1%), containing genes for 13 energy pathway proteins, a small and a large subunit rRNAs, and methionine and tryptophan tRNAs. Mitochondrial ancestral medusozoan gene order (AMGO) was found in the C. sowerbyi, as those found in Cubaia aphrodite (Hydrozoa, Trachylina, Limnomedusae), discomedusan Scyphozoa and Staurozoa. The genes of C. sowerbyi mtDNA are arranged in two clusters with opposite transcriptional polarities, whereby transcription proceeds toward the ends of the DNA molecule. Identical inverted terminal repeats (ITRs) flank the ends of the mitochondrial DNA molecule, a characteristic typical of medusozoans. In addition, two open reading frames (ORFs) of 354 and 1611 bp in length were found downstream of the large subunit rRNA gene, similar to the two ORFs of ORF314 and polB discovered in the linear mtDNA of C. aphrodite, discomedusan Scyphozoa and Staurozoa. Phylogenetic analyses of C. sowerbyi and other cnidarians were carried out based on both nucleotide and inferred amino acid sequences of the 13 mitochondrial energy pathway genes. Our working hypothesis supports the monophyletic Medusozoa being a sister group to Octocorallia (Cnidaria, Anthozoa). Within Medusozoa, the phylogenetic analysis suggests that Staurozoa may be the earliest diverging class and the sister group of all other medusozoans. Cubozoa and coronate Scyphozoa form a clade that is the sister group of Hydrozoa plus discomedusan Scyphozoa. Hydrozoa is the sister group of discomedusan Scyphozoa. Semaeostomeae is a paraphyletic clade with Rhizostomeae, while Limnomedusae (Trachylina) is the sister group of hydroidolinans and may be the earliest diverging lineage among Hydrozoa.

Mitochondrial DNA of Clathrina clathrus (Calcarea, Calcinea): six linear chromosomes, fragmented rRNAs, tRNA editing, and a novel genetic code

Sponges (phylum Porifera) are a large and ancient group of morphologically simple but ecologically important aquatic animals. Although their body plan and lifestyle are relatively uniform, sponges show extensive molecular and genetic diversity. In particular, mitochondrial genomes from three of the four previously studied classes of Porifera (Demospongiae, Hexactinellida, and Homoscleromorpha) have distinct gene contents, genome organizations, and evolutionary rates. Here, we report the mitochondrial genome of Clathrina clathrus (Calcinea, Clathrinidae), a representative of the fourth poriferan class, the Calcarea, which proves to be the most unusual. Clathrina clathrus mitochondrial DNA (mtDNA) consists of six linear chromosomes 7.6-9.4 kb in size and encodes at least 37 genes: 13 protein codings, 2 ribosomal RNAs (rRNAs), and 24 transfer RNAs (tRNAs). Protein genes include atp9, which has now been found in all major sponge lineages, but no atp8. Our analyses further reveal the presence of a novel genetic code that involves unique reassignments of the UAG codons from termination to tyrosine and of the CGN codons from arginine to glycine. Clathrina clathrus mitochondrial rRNAs are encoded in three (srRNA) and ≥6 (lrRNA) fragments distributed out of order and on several chromosomes. The encoded tRNAs contain multiple mismatches in the aminoacyl acceptor stems that are repaired posttranscriptionally by 3'-end RNA editing. Although our analysis does not resolve the phylogenetic position of calcareous sponges, likely due to their high rates of mitochondrial sequence evolution, it confirms mtDNA as a promising marker for population studies in this group. The combination of unusual mitochondrial features in C. clathrus redefines the extremes of mtDNA evolution in animals and further argues against the idea of a "typical animal mtDNA."

Genetic aspects of mitochondrial genome evolution

Many years of extensive studies of metazoan mitochondrial genomes have established differences in gene arrangements and genetic codes as valuable phylogenetic markers. Understanding the underlying mechanisms of replication, transcription and the role of the control regions which cause e.g. different gene orders is important to assess the phylogenetic signal of such events. This review summarises and discusses, for the Metazoa, the general aspects of mitochondrial transcription and replication with respect to control regions as well as several proposed models of gene rearrangements. As whole genome sequencing projects accumulate, more and more observations about mitochondrial gene transfer to the nucleus are reported. Thus occurrence and phylogenetic aspects concerning nuclear mitochondrial-like sequences (NUMTS) is another aspect of this review.

The mitochondrial genome of Paraminabea aldersladei (Cnidaria: Anthozoa: Octocorallia) supports intramolecular recombination as the primary mechanism of gene rearrangement in octocoral mitochondrial genomes

Sequencing of the complete mitochondrial genome of the soft coral Paraminabea aldersladei (Alcyoniidae) revealed a unique gene order, the fifth mt gene arrangement now known within the cnidarian subclass Octocorallia. At 19,886 bp, the mt genome of P. aldersladei is the second largest known for octocorals; its gene content and nucleotide composition are, however, identical to most other octocorals, and the additional length is due to the presence of two large, noncoding intergenic regions. Relative to the presumed ancestral octocoral gene order, in P. aldersladei a block of three protein-coding genes (nad6–nad3–nad4l) has been translocated and inverted. Mapping the distribution of mt gene arrangements onto a taxonomically comprehensive phylogeny of Octocorallia suggests that all of the known octocoral gene orders have evolved by successive inversions of one or more evolutionarily conserved blocks of protein-coding genes. This mode of genome evolution is unique among Metazoa, and contrasts strongly with that observed in Hexacorallia, in which extreme gene shuffling has occurred among taxonomic orders. Two of the four conserved gene blocks found in Octocorallia are, however, also conserved in the linear mt genomes of Medusozoa and in one group of Demospongiae. We speculate that the rate and mechanism of gene rearrangement in octocorals may be influenced by the presence in their mt genomes of mtMutS, a putatively active DNA mismatch repair protein that may also play a role in mediating intramolecular recombination.

Widespread atypical mitochondrial DNA structure in isopods (Crustacea, Peracarida) related to a constitutive heteroplasmy in terrestrial species

Metazoan mitochondrial DNA (mtDNA) is generally composed of circular monomeric molecules. However, a few exceptions do exist and among them two terrestrial isopods Armadillidium vulgare and Porcellionides pruinosus have an atypical mtDNA composed of linear monomers associated with circular "head-to-head" dimers: a very unusual structure for animal mtDNA genome. To assess the distribution of this atypical mtDNA among isopods, we performed RFLP and Southern blot analyses on mtDNA of 16 terrestrial (Oniscidea family) and two aquatic isopod species: the marine Sphaeroma serratum (suborder Flabellifera, sister group of Oniscidea) and the freshwater Asellus aquaticus (Asellota, early derived taxon of isopod). The atypical mtDNA structure was observed in 15 terrestrial isopod species and A. aquaticus, suggesting a wide distribution of atypical mtDNA among isopods. However, a typical metazoan mtDNA structure was detected in the marine isopod S. serratum and the Oniscidea Ligia oceanica . Our results suggest two possible scenarios: an early origin of the atypical mtDNA in isopods followed by reversion to the typical ancestral mtDNA structure for several species, or a convergent appearance of the atypical mtDNA structure in two isopod suborders. We compare this distribution of the atypical mtDNA structure with the presence of a heteroplasmy also observed in the mtDNA of several terrestrial isopod species. We discuss if this transmitted heteroplasmy is vectored by the atypical mtDNA and its impact on the maintenance of the atypical mtDNA in isopods.

The Oxytricha trifallax mitochondrial genome

The Oxytricha trifallax mitochondrial genome contains the largest sequenced ciliate mitochondrial chromosome (~70 kb) plus a ~5-kb linear plasmid bearing mitochondrial telomeres. We identify two new ciliate split genes (rps3 and nad2) as well as four new mitochondrial genes (ribosomal small subunit protein genes: rps- 2, 7, 8, 10), previously undetected in ciliates due to their extreme divergence. The increased size of the Oxytricha mitochondrial genome relative to other ciliates is primarily a consequence of terminal expansions, rather than the retention of ancestral mitochondrial genes. Successive segmental duplications, visible in one of the two Oxytricha mitochondrial subterminal regions, appear to have contributed to the genome expansion. Consistent with pseudogene formation and decay, the subtermini possess shorter, more loosely packed open reading frames than the remainder of the genome. The mitochondrial plasmid shares a 251-bp region with 82% identity to the mitochondrial chromosome, suggesting that it most likely integrated into the chromosome at least once. This region on the chromosome is also close to the end of the most terminal member of a series of duplications, hinting at a possible association between the plasmid and the duplications. The presence of mitochondrial telomeres on the mitochondrial plasmid suggests that such plasmids may be a vehicle for lateral transfer of telomeric sequences between mitochondrial genomes. We conjecture that the extreme divergence observed in ciliate mitochondrial genomes may be due, in part, to repeated invasions by relatively error-prone DNA polymerase-bearing mobile elements.

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.

First complete mitochondrial genome sequence from a box jellyfish reveals a highly fragmented linear architecture and insights into telomere evolution

Animal mitochondrial DNAs (mtDNAs) are typically single circular chromosomes, with the exception of those from medusozoan cnidarians (jellyfish and hydroids), which are linear and sometimes fragmented. Most medusozoans have linear monomeric or linear bipartite mitochondrial genomes, but preliminary data have suggested that box jellyfish (cubozoans) have mtDNAs that consist of many linear chromosomes. Here, we present the complete mtDNA sequence from the winged box jellyfish Alatina moseri (the first from a cubozoan). This genome contains unprecedented levels of fragmentation: 18 unique genes distributed over eight 2.9- to 4.6-kb linear chromosomes. The telomeres are identical within and between chromosomes, and recombination between subtelomeric sequences has led to many genes initiating or terminating with sequences from other genes (the most extreme case being 150 nt of a ribosomal RNA containing the 5′ end of nad2), providing evidence for a gene conversion–based model of telomere evolution. The silent-site nucleotide variation within the A. moseri mtDNA is among the highest observed from a eukaryotic genome and may be associated with elevated rates of recombination.

Evolution of linear mitochondrial genomes in medusozoan cnidarians

In nearly all animals, mitochondrial DNA (mtDNA) consists of a single circular molecule that encodes several subunits of the protein complexes involved in oxidative phosphorylation as well as part of the machinery for their expression. By contrast, mtDNA in species belonging to Medusozoa (one of the two major lineages in the phylum Cnidaria) comprises one to several linear molecules. Many questions remain on the ubiquity of linear mtDNA in medusozoans and the mechanisms responsible for its evolution, replication, and transcription. To address some of these questions, we determined the sequences of nearly complete linear mtDNA from 24 species representing all four medusozoan classes: Cubozoa, Hydrozoa, Scyphozoa, and Staurozoa. All newly determined medusozoan mitochondrial genomes harbor the 17 genes typical for cnidarians and map as linear molecules with a high degree of gene order conservation relative to the anthozoans. In addition, two open reading frames (ORFs), polB and ORF314, are identified in cubozoan, schyphozoan, staurozoan, and trachyline hydrozoan mtDNA. polB belongs to the B-type DNA polymerase gene family, while the product of ORF314 may act as a terminal protein that binds telomeres. We posit that these two ORFs are remnants of a linear plasmid that invaded the mitochondrial genomes of the last common ancestor of Medusozoa and are responsible for its linearity. Hydroidolinan hydrozoans have lost the two ORFs and instead have duplicated cox1 at each end of their mitochondrial chromosome(s). Fragmentation of mtDNA occurred independently in Cubozoa and Hydridae (Hydrozoa, Hydroidolina). Our broad sampling allows us to reconstruct the evolutionary history of linear mtDNA in medusozoans.

Extreme mitochondrial evolution in the ctenophore Mnemiopsis leidyi: Insight from mtDNA and the nuclear genome

Recent advances in sequencing technology have led to a rapid accumulation of mitochondrial DNA (mtDNA) sequences, which now represent the wide spectrum of animal diversity. However, one animal phylum--Ctenophora--has, to date, remained completely unsampled. Ctenophores, a small group of marine animals, are of interest due to their unusual biology, controversial phylogenetic position, and devastating impact as invasive species. Using data from the Mnemiopsis leidyi genome sequencing project, we Polymerase Chain Reaction (PCR) amplified and analyzed its complete mitochondrial (mt-) genome. At just over 10 kb, the mt-genome of M. leidyi is the smallest animal mtDNA ever reported and is among the most derived. It has lost at least 25 genes, including atp6 and all tRNA genes. We show that atp6 has been relocated to the nuclear genome and has acquired introns and a mitochondrial targeting presequence, while tRNA genes have been genuinely lost, along with nuclear-encoded mt-aminoacyl tRNA synthetases. The mt-genome of M. leidyi also displays extremely high rates of sequence evolution, which likely led to the degeneration of both protein and rRNA genes. In particular, encoded rRNA molecules possess little similarity with their homologs in other organisms and have highly reduced secondary structures. At the same time, nuclear encoded mt-ribosomal proteins have undergone expansions, likely to compensate for the reductions in mt-rRNA. The unusual features identified in M. leidyi mtDNA make this organism an interesting system for the study of various aspects of mitochondrial biology, particularly protein and tRNA import and mt-ribosome structures, and add to its value as an emerging model species. Furthermore, the fast-evolving M. leidyi mtDNA should be a convenient molecular marker for species- and population-level studies.

The complete mitochondrial genome of Calicogorgia granulosa (Anthozoa: Octocorallia): potential gene novelty in unidentified ORFs formed by repeat expansion and segmental duplication

Mitochondrial genomes of many nonbilaterian animals show high diversity of genome size and gene content, revealing many intergenic regions (IGRs), diverse repeats and additional genes. Here we present a new complete mitogenome of the cnidarian sea fan species, Calicogorgia granulosa (Anthozoa: Octocorallia) and its novel genomic features. The 20,246 bp of the complete mitogenome, which is the largest among the nine octocorals sequenced to date, contains 13 protein coding genes, 2 rRNAs and a tRNA within its circular form of mitochondrial DNA. We found an identical segmental duplication (S1 and S2, 913 bp) composed of an ORF (672 bp) coding for a hypothetical protein within which Direct Variant Repeat (DVR) expansions reside in-frame to the coding sequence. Additionally, the duplicated segmental DNA showed no variation in nucleotide sequences both between S1 and S2 and across multiple individual samples. Upon these observations, we discuss plausible causes for the intramitochondrial segmental duplication and the absence of sequence variation, and a need for further investigation of the novel ORF as well. In conclusion the present mitogenome of C. granulosa adds more information to our understanding of the diversity and evolution of mitogenomes of nonbilaterian animals.

Complete mitochondrial genomes of two Japanese precious corals, Paracorallium japonicum and Corallium konojoi (Cnidaria, Octocorallia, Coralliidae): notable differences in gene arrangement

Precious coral are taxonomically a group of corals that belong to the family Coralliidae within the order Alcyonacea, subclass Octocorallia, and class Anthozoa, whose skeletal axes are used for jewelry. They are distributed in the Mediterranean Sea and in waters adjacent to Japan, Taiwan, Midway Island and the Hawaiian Islands. The genus Corallium of the family Coralliidae was recently divided into two genera, Corallium and Paracorallium, based on morphological observations, but insufficient molecular evidence to support this classification has been presented to date. We determined for the first time the complete mitochondrial genome sequence of two precious corals P. japonicum and C. konojoi, in order to clarify their systematic positions. The circular mitochondrial genomes of P. japonicum and C. konojoi are 18,913bp and 18,969bp in length, respectively, and encode 13 typical energy pathway protein coding genes (nad1-6, nad4L, cox1-3, cob, atp6 and atp8), two ribosomal RNA genes (rns and rnl), a transfer RNA (trnM) and a mismatch repair gene homologue msh1. The two genomes have an overall nucleotide sequence identity of 97.5%, which is comparable to that between Acanella eburnea and Keratoisidinae sp. belonging to Octocorallia. Surprisingly, however, their gene arrangements were not identical. Phylogenetic analyses using seven complete mitochondrial genome sequences belonging to species in the subclass Octocorallia indicated that within the subclass, at least three gene order rearrangement events occurred during evolution. Our results support the validity of the morphological classification that separated the family Coralliidae into two genera, Corallium and Paracorallium.