Mitochondrial Genome of Savalia savaglia (Cnidaria, Hexacorallia) and Early Metazoan Phylogeny

Genome sizes and repeatome evolution in zoantharians (Cnidaria: Hexacorallia: Zoantharia)

Across eukaryotes, large variations of genome sizes have been observed even between closely related species. Transposable elements as part of the repeated DNA have been proposed and confirmed as one of the most important contributors to genome size variation. However, the evolutionary implications of genome size variation and transposable element dynamics are not well understood. Together with phenotypic traits, they are commonly referred to as the "C-value enigma". The order Zoantharia are benthic cnidarians found from intertidal zones to the deep sea, and some species are particularly abundant in coral reefs. Despite their high ecological relevance, zoantharians have yet to be largely studied from the genomic point of view. This study aims at investigating the role of the repeatome (total content of repeated elements) in genome size variations across the order Zoantharia. To this end, whole-genomes of 32 zoantharian species representing five families were sequenced. Genome sizes were estimated and the abundances of different repeat classes were assessed. In addition, the repeat overlap between species was assessed by a sequence clustering method. The genome sizes in the dataset varied up to 2.4 fold magnitude. Significant correlations between genome size, repeated DNA content and transposable elements, respectively (Pearson's correlation test R2 = 0.47, p = 0.0016; R2 = 0.22, p = 0.05) were found, suggesting their involvement in the dynamics of genome expansion and reduction. In all species, long interspersed nuclear elements and DNA transposons were the most abundant identified elements. These transposable elements also appeared to have had a recent expansion event. This was in contrast to the comparative clustering analysis which revealed species-specific patterns of satellite elements' amplification. In summary, the genome sizes of zoantharians likely result from the complex dynamics of repeated elements. Finally, the majority of repeated elements (up to 70%) could not be annotated to a known repeat class, highlighting the need to further investigate non-model cnidarian genomes. More research is needed to understand how repeated DNA dynamics relate to zoantharian evolution and their biology.

The earliest diverging extant scleractinian corals recovered by mitochondrial genomes

Evolutionary reconstructions of scleractinian corals have a discrepant proportion of zooxanthellate reef-building species in relation to their azooxanthellate deep-sea counterparts. In particular, the earliest diverging "Basal" lineage remains poorly studied compared to "Robust" and "Complex" corals. The lack of data from corals other than reef-building species impairs a broader understanding of scleractinian evolution. Here, based on complete mitogenomes, the early onset of azooxanthellate corals is explored focusing on one of the most morphologically distinct families, Micrabaciidae. Sequenced on both Illumina and Sanger platforms, mitogenomes of four micrabaciids range from 19,048 to 19,542 bp and have gene content and order similar to the majority of scleractinians. Phylogenies containing all mitochondrial genes confirm the monophyly of Micrabaciidae as a sister group to the rest of Scleractinia. This topology not only corroborates the hypothesis of a solitary and azooxanthellate ancestor for the order, but also agrees with the unique skeletal microstructure previously found in the family. Moreover, the early-diverging position of micrabaciids followed by gardineriids reinforces the previously observed macromorphological similarities between micrabaciids and Corallimorpharia as well as its microstructural differences with Gardineriidae. The fact that both families share features with family Kilbuchophylliidae ultimately points towards a Middle Ordovician origin for Scleractinia.

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.

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

Six complete and three partial actiniarian mitochondrial genomes were amplified in two semi-circles using long-range PCR and pyrosequenced in a single run on a 454 GS Junior, doubling the number of complete mitogenomes available within the order. Typical metazoan mtDNA features included circularity, 13 protein-coding genes, 2 ribosomal RNA genes, and length ranging from 17,498 to 19,727 bp. Several typical anthozoan mitochondrial genome features were also observed including the presence of only two transfer RNA genes, elevated A + T richness ranging from 54.9 to 62.4%, large intergenic regions, and group 1 introns interrupting NADH dehydrogenase subunit 5 and cytochrome c oxidase subunit I, the latter of which possesses a homing endonuclease gene. Within the sea anemone Alicia sansibarensis, we report the first mitochondrial gene order rearrangement within the Actiniaria, as well as putative novel non-canonical protein-coding genes. Phylogenetic analyses of all 13 protein-coding and 2 ribosomal genes largely corroborated current hypotheses of sea anemone interrelatedness, with a few lower-level differences.

Metabolomic profiling reveals deep chemical divergence between two morphotypes of the zoanthid Parazoanthus axinellae

Metabolomics has recently proven its usefulness as complementary tool to traditional morphological and genetic analyses for the classification of marine invertebrates. Among the metabolite-rich cnidarian order Zoantharia, Parazoanthus is a polyphyletic genus whose systematics and phylogeny remain controversial. Within this genus, one of the most studied species, Parazoanthus axinellae is prominent in rocky shallow waters of the Mediterranean Sea and the NE Atlantic Ocean. Although different morphotypes can easily be distinguished, only one species is recognized to date. Here, a metabolomic profiling approach has been used to assess the chemical diversity of two main Mediterranean morphotypes, the "slender" and "stocky" forms of P. axinellae. Targeted profiling of their major secondary metabolites revealed a significant chemical divergence between the morphotypes. While zoanthoxanthin alkaloids and ecdysteroids are abundant in both morphs, the "slender" morphotype is characterized by the presence of additional and bioactive 3,5-disubstituted hydantoin derivatives named parazoanthines. The absence of these specific compounds in the "stocky" morphotype was confirmed by spatial and temporal monitoring over an annual cycle. Moreover, specimens of the "slender" morphotype are also the only ones found as epibionts of several sponge species, particularly Cymbaxinella damicornis thus suggesting a putative ecological link.

The "naked coral" hypothesis revisited--evidence for and against scleractinian monophyly

The relationship between Scleractinia and Corallimorpharia, Orders within Anthozoa distinguished by the presence of an aragonite skeleton in the former, is controversial. Although classically considered distinct groups, some phylogenetic analyses have placed the Corallimorpharia within a larger Scleractinia/Corallimorpharia clade, leading to the suggestion that the Corallimorpharia are “naked corals” that arose via skeleton loss during the Cretaceous from a Scleractinian ancestor. Scleractinian paraphyly is, however, contradicted by a number of recent phylogenetic studies based on mt nucleotide (nt) sequence data. Whereas the “naked coral” hypothesis was based on analysis of the sequences of proteins encoded by a relatively small number of mt genomes, here a much-expanded dataset was used to reinvestigate hexacorallian phylogeny. The initial observation was that, whereas analyses based on nt data support scleractinian monophyly, those based on amino acid (aa) data support the “naked coral” hypothesis, irrespective of the method and with very strong support. To better understand the bases of these contrasting results, the effects of systematic errors were examined. Compared to other hexacorallians, the mt genomes of “Robust” corals have a higher (A+T) content, codon usage is far more constrained, and the proteins that they encode have a markedly higher phenylalanine content, leading us to suggest that mt DNA repair may be impaired in this lineage. Thus the “naked coral” topology could be caused by high levels of saturation in these mitochondrial sequences, long-branch effects or model violations. The equivocal results of these extensive analyses highlight the fundamental problems of basing coral phylogeny on mitochondrial sequence data.

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.

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.

Updating our view of organelle genome nucleotide landscape

Organelle genomes show remarkable variation in architecture and coding content, yet their nucleotide composition is relatively unvarying across the eukaryotic domain, with most having a high adenine and thymine (AT) content. Recent studies, however, have uncovered guanine and cytosine (GC)-rich mitochondrial and plastid genomes. These sequences come from a small but eclectic list of species, including certain green plants and animals. Here, I review GC-rich organelle DNAs and the insights they have provided into the evolution of nucleotide landscape. I emphasize that GC-biased mitochondrial and plastid DNAs are more widespread than once thought, sometimes occurring together in the same species, and suggest that the forces biasing their nucleotide content can differ both among and within lineages, and may be associated with specific genome architectural features and life history traits.

Unscrambling genetic information at the RNA level

Genomics aims at unraveling the blueprint of life; however, DNA sequence alone does not always reveal the proteins and structural RNAs encoded by the genome. The reason is that genetic information is often encrypted. Recognizing the logic of encryption, and understanding how living cells decode hidden information--at the level of DNA, RNA or protein--is challenging. RNA-level decryption includes topical RNA editing and more 'macroscopic' transcript rearrangements. The latter events involve the four types of introns recognized to date, notably spliceosomal, group I, group II, and archaeal/tRNA splicing. Intricate variants, such as alternative splicing and trans-splicing, have been reported for each intron type, but the biological significance has not always been confirmed. Novel RNA-level unscrambling processes were recently discovered in mitochondria of dinoflagellates and diplonemids, and potentially euglenids. These processes seem not to rely on known introns, and the corresponding molecular mechanisms remain to be elucidated.

Mitogenome rearrangement in the cold-water scleractinian coral Lophelia pertusa (Cnidaria, Anthozoa) involves a long-term evolving group I intron

Group I introns are genetic insertion elements that invade host genomes in a wide range of organisms. In metazoans, however, group I introns are extremely rare, so far only identified within mitogenomes of hexacorals and some sponges. We sequenced the complete mitogenome of the cold-water scleractinian coral Lophelia pertusa, the dominating deep sea reef-building coral species in the North Atlantic Ocean. The mitogenome (16,150 bp) has the same gene content but organized in a unique gene order compared to that of other known scleractinian corals. A complex group I intron (6460 bp) inserted in the ND5 gene (position 717) was found to host seven essential mitochondrial protein genes and one ribosomal RNA gene. Phylogenetic analysis supports a vertical inheritance pattern of the ND5-717 intron among hexacoral mitogenomes with no examples of intron loss. Structural assessments of the Lophelia intron revealed an unusual organization that lacks the universally conserved ωG at the 3' end, as well as a highly compact RNA core structure with overlapping ribozyme and protein coding capacities. Based on phylogenetic and structural analyses we reconstructed the evolutionary history of ND5-717, from its ancestral protist origin, through intron loss in some early metazoan lineages, and into a compulsory feature with functional implications in hexacorals.

Timing of spawning and early development of Palythoa tuberculosa (Anthozoa, Zoantharia, Sphenopidae) in Okinawa, Japan

The spawning behavior and early embryogenesis of Palythoa tuberculosa (Anthozoa, Zoantharia) were observed in August 2009 off Okinawa Island, Japan. P. tuberculosa released zygotes just after high tide around new moon nights. The mean diameter of zygotes was 365.6 ± s.d.14.8 μm, and zygotes did not contain any symbiotic algae (zooxanthellae). About 2 h after spawning, the first cleavage furrow appeared on one side of the zygotes, although it was uncertain when eggs were fertilized. After second cleavage, the arrangement of blastomeres was pseudospherical. At 9 h after spawning, the embryo became a concave-convex dish shape, then gastrulation occurred and the blastopore was formed. Seven-day old larvae were ellipsoid and about 700 μm long, with an open mouth at one end. Two weeks after spawning, the larvae developed a longitudinal band of long cilia (= ventral ciliate band) that is characteristic of zoanthella larvae. In P. tuberculosa, larvae show a non-radial body plan and then metamorphose to almost-radial (in outward appearance) polyps after settlement. These results may support a hypothesis that a common ancestor of Cnidaria had a bilateral body plan that has been secondarily lost in some extant cnidarians.

Truncation of subunit ND2 disrupts the threefold symmetry of the antiporter-like subunits in complex I from higher metazoans

Three of the conserved, membrane-bound subunits in NADH:ubiquinone oxidoreductase (complex I) are related to one another, and to Mrp sodium-proton antiporters. Recent structural analysis of two prokaryotic complexes I revealed that the three subunits each contain fourteen transmembrane helices that overlay in structural alignments: the translocation of three protons may be coordinated by a lateral helix connecting them together (Efremov, R.G., Baradaran, R. and Sazanov, L.A. (2010). The architecture of respiratory complex I. Nature 465, 441-447). Here, we show that in higher metazoans the threefold symmetry is broken by the loss of three helices from subunit ND2; possible implications for the mechanism of proton translocation are discussed.

Unique mitogenomic features in the scleractinian family pocilloporidae (scleractinia: astrocoeniina)

The complete DNA sequences of three mitochondrial (mt) genomes were obtained from the scleractinian corals, Stylophora pistillata, Pocillopora damicornis, and Madracis mirabilis, and were compared to the published mt genomes to elucidate phylogenetically unique features of the family Pocilloporidae. The entire mt genomes of pocilloporid corals ranged from 16,951 to 17,425 bp with the A+T contents of their sense strands ranging from 68.4% to 70.2%. The gene order of protein-coding genes was identical to those of other scleractinian corals. The novel atp8 gene, first described in confamilial Seriatopora corals, was also confirmed using reverse transcription-polymerase chain reaction (RT-PCR), Northern blot, and sequence analyses in other genera of the Pocilloporidae. The intergenic spacer between atp6 and nad4, containing distinct repeated elements, conserved sequence blocks and domains, and functional structures, possesses typical characteristics of a putative control region for the four coral genera. A duplicated trnW, detected in the region close to the cox1 gene and which shares the highly conserved primary and secondary structures of its original counterpart, was discovered in both Seriatopora and Stylophora. These molecular characteristics are unique and provide the phylogenetic information for future evaluation of the status of the family Pocilloporidae in the evolutionary history of scleractinian corals.

Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species

The mitochondrial genome (mtDNA) of Metazoa is a good model system for evolutionary genomic studies and the availability of more than 1000 sequences provides an almost unique opportunity to decode the mechanisms of genome evolution over a large phylogenetic range. In this paper, we review several structural features of the metazoan mtDNA, such as gene content, genome size, genome architecture and the new parameter of gene strand asymmetry in a phylogenetic framework. The data reviewed here show that: (1) the plasticity of Metazoa mtDNA is higher than previously thought and mainly due to variation in number and location of tRNA genes; (2) an exceptional trend towards stabilization of genomic features occurred in deuterostomes and was exacerbated in vertebrates, where gene content, genome architecture and gene strand asymmetry are almost invariant. Only tunicates exhibit a very high degree of genome variability comparable to that found outside deuterostomes. In order to analyse the genomic evolutionary process at short evolutionary distances, we have also compared mtDNAs of species belonging to the same genus: the variability observed in congeneric species significantly recapitulates the evolutionary dynamics observed at higher taxonomic ranks, especially for taxa showing high levels of genome plasticity and/or fast nucleotide substitution rates. Thus, the analysis of congeneric species promises to be a valuable approach for the assessment of the mtDNA evolutionary trend in poorly or not yet sampled metazoan groups.

The mitochondrial genome of Pocillopora (Cnidaria: Scleractinia) contains two variable regions: the putative D-loop and a novel ORF of unknown function

The complete mitochondrial genomes of two individuals attributed to different morphospecies of the scleractinian coral genus Pocillopora have been sequenced. Both genomes, respectively 17,415 and 17,422 nt long, share the presence of a previously undescribed ORF encoding a putative protein made up of 302 amino acids and of unknown function. Surprisingly, this ORF turns out to be the second most variable region of the mitochondrial genome (1% nucleotide sequence difference between the two individuals) after the putative control region (1.5% sequence difference). Except for the presence of this ORF and for the location of the putative control region, the mitochondrial genome of Pocillopora is organized in a fashion similar to the other scleractinian coral genomes published to date. For the first time in a cnidarian, a putative second origin of replication is described based on its secondary structure similar to the stem-loop structure of O(L), the origin of L-strand replication in vertebrates.