Multiplexed pyrosequencing of nine sea anemone (Cnidaria: Anthozoa: Hexacorallia: Actiniaria) mitochondrial genomes

New mitochondrial gene order arrangements and evolutionary implications in the class Octocorallia

The complete mitochondrial genomes of octocorals typically range from 18.5 kb to 20.5 kb in length and include 14 protein-coding genes (PCGs), two ribosomal RNA genes and one tRNA. To date, seven different gene orders (A-G) have been described, yet comprehensive investigations of the actual number of arrangements, as well as comparative analyses and evolutionary reconstructions of mitochondrial genome evolution within the whole class Octocorallia, have been often overlooked. Here, we considered the complete mitochondrial genomes available for octocorals and explored their structure and gene order variability. Our results updated the actual number of mitochondrial gene order arrangements so far known for octocorals from 7 to 14 and allowed us to explore and preliminarily discuss the role of some of the structural and functional factors in the mitogenomes. We performed comparative mitogenomic analyses on the existing and novel octocoral gene orders, considering different mitogenomic structural features such as genome size, GC percentage, AT and GC skewness. The mitochondrial gene order history mapped on a recently published nuclear loci phylogeny showed that the most common rearrangement events in octocorals are inversions, inverted transpositions and transpositions. Furthermore, gene order rearrangement events were restricted only to some regions of the tree. Overall, different rearrangement events arose independently and from the ancestral and most common gene order, instead of being derived from other rearranged orders. Finally, our data demonstrate how the study of mitochondrial gene orders can be used to explore the evolution of octocorals and in some cases can be used to assess the phylogenetic placement of certain taxa.

The complete mitochondrial genome of an economic sea anemone (Paracondylactis sinensis) in the East China Sea

Paracondylactis sinensis Carlgren, 1934 (Actiniidae, Actiniaria) is an edible sea anemone in China. Their wild population has intensively decreased in recent years due to overharvesting. In this study, the complete mitochondrial genome of this economic species collected in the coast of Zhejiang, China is sequenced and obtained using high throughput methods. The total length of this circular molecule is 20,786 bp. Thirteen protein coding genes, two ribosomal RNA genes, two transfer RNA (tRNATrp, tRNAMet) genes and a putative ORF are annotated in it. Phylogenetic analysis based on the amino acids of mitochondrial genomes indicates that this species belongs to the family of Actiniidae. This result is consistent with the previous work that identified the edible sea anemone as Paracondylactis sinensis although it has always been recognized as Calliactis sinensis (of family Hormathiidae) in most Chinese reports. Overall, the mitochondrial genome produced in this study assists in clarifying the phylogenetic status of this sea anemone and provides a molecular foundation for future protection and breeding work.

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.

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.

The mitochondrial genome of the orange-striped green sea anemone Diadumene lineata (Actiniaria: Diadumenidae): the first complete sequence in the family Diadumenidae

The complete mitogenome of the orange-striped green sea anemone (Diadumene lineata) has been sequenced and annotated for the first time. The total length of the mitogenome is 17,552 bp with an A+T content of 62.6%. Unlike typical metazoan mitogenome, this mitogenome include 14 protein-coding genes (13 energy pathway protein coding genes, and a heg gene), two tRNAs, two rRNAs, and 19 intergenic regions. The COX1 gene possesses a homing endonuclease gene. This circular genome contains two introns, one in ND5 and another in COX1.This sequence is the first sequenced complete mitogenome in Diadumenidae and provides fundamental data for exploring complicated evolutionary relationships in Actiniaria.

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.

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.

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.

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.

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.

Complete mitochondrial genome of the sea anemone, Anthopleura midori (Actiniaria: Actiniidae)

Here we sequenced and characterized the complete mitochondrial genome (mitogenome) of Anthopleura midori from Weihai, China. This is the first species with sequenced mitogenome of genus Anthopleura. The A. midori mitogenome is 20 039 bp in length (GenBank accession no. KT989511). The overall base composition of the H-strand is 26.81% A, 17.45% C, 21.52% G, and 34.22% T. The base composition clearly showed the A-T skew. Similar to most other sea anemone in family Actiniidae, it contains 13 protein-coding genes, 2 rRNAs, 2 tRNAs, and some non-coding regions. The circular mitogenome codes entirely on the heavy strand and ND5 and COI genes were embedded by an intron. Phylogenetic analyses show that A. midori is related most closely to Phymanthus crucifer.

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.