Large dataset of octocoral mitochondrial genomes provides new insights into mt-mutS evolution and function

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 the bubble-gum coral Paragorgia papillata (Octocorallia: Coralliidae) from the seamount in the tropical Western Pacific

The complete mitochondrial genome of Paragorgia papillata Li et al. 2021, a deep-sea gorgonian inhabiting at 858 m in Caroline Ridge, was obtained in this study. The length of the mitochondrial genome is 19,018 bp with 14 protein coding genes, one transfer RNA (tRNA-Met) and two ribosomal RNA genes contained in this circular molecule. Phylogenetic analysis indicated that P. papillata and P. coralloides Bayer, 1993 were two closely related species, and a total of 26 mutational sites (four nonsynonymous mutations included) can be detected between their mitochondrial genomes. This exhibits a case that mitochondrial genomes can be applied to differentiate closely related species in gorgonians. The phylogenetic tree constructed with mitochondrial genomes showed that the families in Octocorallia are reciprocally monophyletic, provided that the family names were revised according to the systematic revision of Octocorallia guided by phylogenomics. However, the relationships of the families within each order were different between the previous phylogenomic work and ours. Integrating mitochondrial genomes from a wider array of Octocorallia families is essential for a more accurate comparison of phylogenies derived from nuclear and mitochondrial sequences in future study.

The Mitogenome of the Subarctic Octocoral Alcyonium digitatum Reveals a Putative tRNAPro Gene Nested within MutS

We sequenced and analyzed the complete mitogenome of a Norwegian isolate of the octocoral Alcyonium digitatum using the Ion Torrent sequencing technology. The 18,790 bp circular mitochondrial genome was found to harbor the same set of 17 genes, which encode 14 protein subunits, two structural ribosomal RNAs and one tRNA, as reported in other octocorals. In addition, we detected a new tRNAPro-like gene sequence nested within the MutS protein coding region. This putative tRNA gene feature appears to be conserved among the octocorals but has not been reported previously. The A. digitatum mitogenome was also shown to harbor an optional gene (ORFA) that encodes a putative protein of 191 amino acids with unknown function. A mitogenome-based phylogenetic analysis, presented as a maximum likelihood tree, showed that A. digitatum clustered with high statistical confidence with two other Alcyonium species endemic to the Mediterranean Sea and the Southeast Pacific Ocean.

Genome copy number predicts extreme evolutionary rate variation in plant mitochondrial DNA

Nuclear and organellar genomes can evolve at vastly different rates despite occupying the same cell. In most bilaterian animals, mitochondrial DNA (mtDNA) evolves faster than nuclear DNA, whereas this trend is generally reversed in plants. However, in some exceptional angiosperm clades, mtDNA substitution rates have increased up to 5,000-fold compared with closely related lineages. The mechanisms responsible for this acceleration are generally unknown. Because plants rely on homologous recombination to repair mtDNA damage, we hypothesized that mtDNA copy numbers may predict evolutionary rates, as lower copy numbers may provide fewer templates for such repair mechanisms. In support of this hypothesis, we found that copy number explains 47% of the variation in synonymous substitution rates of mtDNA across 60 diverse seed plant species representing ~300 million years of evolution. Copy number was also negatively correlated with mitogenome size, which may be a cause or consequence of mutation rate variation. Both relationships were unique to mtDNA and not observed in plastid DNA. These results suggest that homologous recombinational repair plays a role in driving mtDNA substitution rates in plants and may explain variation in mtDNA evolution more broadly across eukaryotes. Our findings also contribute to broader questions about the relationships between mutation rates, genome size, selection efficiency, and the drift-barrier hypothesis.

Bamboozled! Resolving deep evolutionary nodes within the phylogeny of bamboo corals (Octocorallia: Scleralcyonacea: Keratoisididae)

Keratoisididae is a globally distributed, and exclusively deep-sea, family of octocorals that contains species and genera that are polyphyletic. An alphanumeric system, based on a three-gene-region phylogeny, is widely used to describe the biodiversity within this family. That phylogeny identified 12 major groups although it did not have enough signal to explore the relationships among groups. Using increased phylogenomic resolution generated from Ultraconserved Elements and exons (i.e. conserved elements), we aim to resolve deeper nodes within the family and investigate the relationships among those predefined groups. In total, 109 libraries of conserved elements were generated from individuals representing both the genetic and morphological diversity of our keratoisidids. In addition, the conserved element data of 12 individuals from previous studies were included. Our taxon sampling included 11 of the 12 keratoisidid groups. We present two phylogenies, constructed from a 75% (231 loci) and 50% (1729 loci) taxon occupancy matrix respectively, using both Maximum Likelihood and Multiple Species Coalescence methods. These trees were congruent at deep nodes. As expected, S1 keratoisidids were recovered as a well-supported sister clade to the rest of the bamboo corals. S1 corals do not share the same mitochondrial gene arrangement found in other members of Keratoisididae. All other bamboo corals were recovered within two major clades. Clade I comprises individuals assigned to alphanumeric groups B1, C1, D1&D2, F1, H1, I4, and J3 while Clade II contains representatives from A1, I1, and M1. By combining genomics with already published morphological data, we provide evidence that group H1 is not monophyletic, and that the division between other groups - D1 and D2, and A1 and M1 - needs to be reconsidered. Overall, there is a lack of robust morphological markers within Keratoisididae, but subtle characters such as sclerite microstructure and ornamentation seem to be shared within groups and warrant further investigation as taxonomically diagnostic characters.

Complete mitochondrial genome of the abyssal coral Abyssoprimnoagemina Cairns, 2015 (Octocorallia, Primnoidae) from the Clarion-Clipperton Zone, Pacific Ocean

The Clarion-Clipperton Zone (CCZ) in the tropical East Pacific is a region of interest for deep-sea mining due to its underwater deposits of polymetallic nodules containing economically important metals such as nickel, copper, and cobalt. It is also a region of extensive baseline studies aiming to describe the state of the environment, including the biodiversity of the benthic fauna. An abundant component of the abyssal plain ecosystem consists of sessile fauna which encrusts polymetallic nodules and are vulnerable to potential impacts arising from exploitation activities, particularly removal of substrate. Therefore, this fauna is often considered to have key species whose genetic connectivity should be studied to assess their ecological resilience. One such species is Abyssoprimnoagemina Cairns, 2015, a deep-sea coral from the CCZ whose presence in the Interoceanmetal Joint Organization (IOM) claim area has been confirmed during samplings. In this study, we used next-generation sequencing (NGS) to obtain the 18S nuclear rRNA gene and the complete mitochondrial genome of A.gemina from IOM exploration area. The mitogenome is 18,825 bp long and encodes for 14 protein coding genes, 2 rRNAs, and a single tRNA. The two phylogeny reconstructions derived from these data confirm previous studies and display A.gemina within a highly supported cluster of seven species whose mitogenomes are all colinear and of comparable size. This study also demonstrates the suitability of NGS for DNA barcoding of the benthic megafauna of the CCZ, which could become part of the IOM protocol for the assessment of population diversity and genetic connectivity in its claim area.

The complete mitochondrial genome of soft coral, Eleutherobia rubra (Brundin, 1896) (Cnidaria; Anthozoa; Malacalcyonacea; Alcyoniidae)

The mitogenome of a soft coral, Eleutherobia rubra (Brundin, 1896), was completely sequenced for the first time. The total mitogenome length of E. rubra is 18,724 bp with 14 protein-coding genes, two ribosomal RNA genes, one transfer RNA gene (tRNA-Met), and one non-coding region (NCR). The gene order is also consistent with other Alcyoniidae species. The base composition is 30.1% A, 16.7% C, 19.5% G, and 33.7% T, with a G-C content of 36.2%. This is the first record of the complete mitogenome sequence of the genus Eleutherobia.

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.

Selection in coral mitogenomes, with insights into adaptations in the deep sea

Corals are a dominant benthic fauna that occur across a vast range of depths from just below the ocean's surface to the abyssopelagic zone. However, little is known about the evolutionary mechanisms that enable them to inhabit such a wide range of environments. The mitochondrial (mt) genome, which is involved in energetic pathways, may be subject to selection pressures at greater depths to meet the metabolic demands of that environment. Here, we use a phylogenomic framework combined with codon-based models to evaluate whether mt protein-coding genes (PCGs) associated with cellular energy functions are under positive selection across depth in three groups of corals: Octocorallia, Scleractinia, and Antipatharia. The results demonstrated that mt PCGs of deep- and shallow-water species of all three groups were primarily under strong purifying selection (0.0474 < ω < 0.3123), with the exception of positive selection in atp6 (ω = 1.3263) of deep-sea antipatharians. We also found evidence for positive selection at fifteen sites across cox1, mtMutS, and nad1 in deep-sea octocorals and nad3 of deep-sea antipatharians. These results contribute to our limited understanding of mt adaptations as a function of depth and provide insight into the molecular response of corals to the extreme deep-sea environment.

Ancient endosymbiont-mediated transmission of a selfish gene provides a model for overcoming barriers to gene transfer into animal mitochondrial genomes

In contrast to bilaterian animals, non-bilaterian mitochondrial genomes contain atypical genes, often attributed to horizontal gene transfer (HGT) as an ad hoc explanation. Although prevalent in plants, HGT into animal mitochondrial genomes is rare, lacking suitable explanatory models for their occurrence. HGT of the mismatch DNA repair gene (mtMutS) from giant viruses to octocoral (soft corals and their kin) mitochondrial genomes provides a model for how barriers to HGT to animal mitochondria may be overcome. A review of the available literature suggests that this HGT was mediated by an alveolate endosymbiont infected with a lysogenic phycodnavirus that enabled insertion of the homing endonuclease containing mtMutS into octocoral mitochondrial genomes. We posit that homing endonuclease domains and similar selfish elements play a crucial role in such inter-domain gene transfers. Understanding the role of selfish genetic elements in HGT has the potential to aid development of tools for manipulating animal mitochondrial DNA.

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.

Complete mitochondrial genome of the branching octocoral Paramuricea grayi (Johnson, 1861), phylogenetic relationships and divergence analysis

The Gray's sea fan, Paramuricea grayi (Johnson, 1861), typically inhabits deep littoral and circalittoral habitats of the eastern temperate and tropical Atlantic Ocean. Along the Iberian Peninsula, where P. grayi is a dominant constituent of circalittoral coral gardens, two segregating lineages (yellow and purple morphotypes) were recently identified using single-copy nuclear orthologues. The mitochondrial genomes of 9 P. grayi individuals covering both color morphotypes were assembled from RNA-seq data, using samples collected at three sites in southern (Sagres and Tavira) and western (Cape Espichel) Portugal. The complete circular mitogenome is 18,668 bp in length, has an A + T-rich base composition (62.5%) and contains the 17 genes typically found in Octocorallia: 14 protein-coding genes (atp6, atp8, cob, cox1-3, mt-mutS, nad1-6, and nad4L), the small and large subunit rRNAs (rns and rnl), and one transfer RNA (trnM). The mitogenomes were nearly identical for all specimens, though we identified a noteworthy polymorphism (two SNPs 9 bp apart) in the mt-mutS of one purple individual that is shared with the sister species P. clavata. The mitogenomes of the two species have a pairwise sequence identity of 99.0%, with nad6 and mt-mutS having the highest rates of non-synonymous substitutions.