Sea anemones possess dynamic mitogenome structures

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.

Confirmation of the first intronic sequence in the bivalvian mitochondrial genome of Macoma balthica (Linnaeus, 1758)

In 2020, the first male-type mitochondrial genome from the clam Macoma balthica was published. Apart from the unusual doubly uniparental inheritance of mtDNA, scientists observed a unique (over 4k bp long) extension in the middle of the cox2 gene. We have attempted to replicate these data by NGS DNA sequencing and explore further the expression of the long cox2 gene. In our study, we report an even longer cox2 gene (over 5.5 kbp) with no stop codon separating conserved cox2 domains, as well as, based on the rtPCR, a lower relative gene expression pattern of the middle part of the gene (5' = 1; mid = 0.46; 3' = 0.89). Lastly, we sequenced the cox2 gene transcript proving the excision of the intronic sequence.

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

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

The mitochondrial genome of the grape powdery mildew pathogen Erysiphe necator is intron rich and exhibits a distinct gene organization

Powdery mildews are notorious fungal plant pathogens but only limited information exists on their genomes. Here we present the mitochondrial genome of the grape powdery mildew fungus Erysiphe necator and a high-quality mitochondrial gene annotation generated through cloning and Sanger sequencing of full-length cDNA clones. The E. necator mitochondrial genome consists of a circular DNA sequence of 188,577 bp that harbors a core set of 14 protein-coding genes that are typically present in fungal mitochondrial genomes, along with genes encoding the small and large ribosomal subunits, a ribosomal protein S3, and 25 mitochondrial-encoded transfer RNAs (mt-tRNAs). Interestingly, it also exhibits a distinct gene organization with atypical bicistronic-like expression of the nad4L/nad5 and atp6/nad3 gene pairs, and contains a large number of 70 introns, making it one of the richest in introns mitochondrial genomes among fungi. Sixty-four intronic ORFs were also found, most of which encoded homing endonucleases of the LAGLIDADG or GIY-YIG families. Further comparative analysis of five E. necator isolates revealed 203 polymorphic sites, but only five were located within exons of the core mitochondrial genes. These results provide insights into the organization of mitochondrial genomes of powdery mildews and represent valuable resources for population genetic and evolutionary studies.

Testing small-scale ecological gradients and intraspecific differentiation for hundreds of kelp forest species using haplotypes from metabarcoding

DNA metabarcoding has been increasingly used to detail distributions of hundreds of species. Most analyses focus on creating molecular operational taxonomic units (MOTUs) from complex mixtures of DNA sequences, but much less common is use of the sequence diversity within these MOTUs. Here we use the diversity of COI haplotypes within MOTUs from a California kelp forest to infer patterns of population abundance, dispersal and population history from 527 species of animals and algae from 106 samples of benthic habitats in Monterey Bay. Using haplotypes as a unit we show fine-grained differences of abundance across locations for 15 species, and marked aggregation from sample to sample for most of the common species of plants and animals. Previous analyses could not distinguish these patterns from artefacts of amplification or sequence bias. Our haplotype data also reveal strong population genetic differentiation over small spatial scales for 48 species of red algae, sponges and Bryozoa. Last, phylogenetic analysis of mismatch frequencies among haplotypes show a wide variety of demographic histories from recent expansions to long, stable population sizes. These analyses show that abundant, small-bodied marine species that are often overlooked in ecological surveys can have strikingly different patterns of ecological and genetic structure leading to population, ecological and perhaps adaptive differences between habitats. MOTU diversity data from the same sequencing efforts that generate species-level analyses can greatly increase the scope and value of metabarcoding studies.

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 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.

The fungal mitochondrial Nad5 pan-genic intron landscape

An intron landscape was prepared for the fungal mitochondrial nad5 gene. A hundred and eighty-eight fungal species were examined and a total of 265 introns were noted to be located in 29 intron insertion sites within the examined nad5 genes. Two hundred and sixty-three introns could be classified as group I types and two group II introns were noted. One additional group II intron module was identified nested within a composite group I intron. Based on features related to RNA secondary structures, introns can be classified into different subtypes and it was observed that intron insertion-sites are biased towards phase 0 and they appear to be specific to an intron type. Intron landscapes could be used as a guide map to predict the location of fungal mtDNA mobile introns, which are composite elements that include a ribozyme component and in some instances open reading frames encoding homing endonucleases or reverse transcriptases and all of these have applications in biotechnology.

Putative group I introns in the eukaryote nuclear internal transcribed spacers

Group I introns are mobile genetic elements that interrupt genes encoding proteins and RNAs. In the rRNA operon, introns can insert in the small subunit (SSU) and large subunit (LSU) of a wide variety of protists and various prokaryotes, but they were never found in the ITS region. In this study, unusually long ITS regions of fungi and closely related unicellular organisms (Polychytrium aggregatum, Mitosporidium daphniae, Amoeboaphelidium occidentale and Nuclearia simplex) were analysed. While the insertion of repeats is responsible for long ITS in other eukaryotes, the increased size of the sequences analysed herein seems rather due to the presence of introns in ITS-1 or ITS-2. The identified insertions can be folded in secondary structures according to group I intron models, and they cluster within introns in conserved core-based phylogeny. In addition, for Mitosporidium, Amoeboaphelidium and Nuclearia, more conventional ITS-2 structures can be deduced once spacer introns are removed. Sequences of five shark species were also analysed for their structure and included in phylogeny because of unpublished work reporting introns in their ITS, obtaining congruent results. Overall, the data presented herein indicate that spacer regions may contain introns.

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

Background

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

Results

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

Conclusion

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

Electronic supplementary material

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

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

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

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

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

Cancer-specific SNPs originate from low-level heteroplasmic variants in human mitochondrial genomes of a matched cell line pair

Low-level mitochondrial heteroplasmy is a common phenomenon in both normal and cancer cells. Here, we investigate the link between low-level heteroplasmy and mitogenome mutations in a human breast cancer matched cell line by high-throughput sequencing. We identified 23 heteroplasmic sites, of which 15 were common between normal cells (Hs578Bst) and cancer cells (Hs578T). Most sites were clustered within the highly conserved Complex IV and ribosomal RNA genes. Two heteroplasmic variants in normal cells were found as fixed mutations in cancer cells. This indicates a positive selection of these variants in cancer cells. RNA-Seq analysis identified upregulated L-strand specific transcripts in cancer cells, which include three mitochondrial long non-coding RNA molecules. We hypothesize that this is due to two cancer cell-specific mutations in the control region.

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 complete mitochondrial genome of the alvinocaridid shrimp Shinkaicaris leurokolos (Decapoda, Caridea): Insight into the mitochondrial genetic basis of deep-sea hydrothermal vent adaptation in the shrimp

Deep-sea hydrothermal vent is one of the most extreme environments on Earth with low oxygen and high levels of toxins. Decapod species from the family Alvinocarididae have colonized and successfully adapted to this extremely harsh environment. Mitochondria plays a vital role in oxygen usage and energy metabolism, thus it may be under selection in the adaptive evolution of the hydrothermal vent shrimps. In this study, the mitochondrial genome (mitogenome) of alvinocaridid shrimp Shinkaicaris leurokolos (Kikuchi & Hashimoto, 2000) was determined through Illumina sequencing. The mitogenome of S. leurokolos was 15,903bp in length, containing 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. The gene order and orientation were identical to those of sequenced alvinocaridids. It has the longest concatenated sequences of protein-coding genes, tRNAs and shortest pooled rRNAs among the alvinocaridids. The control regions (CRs) of alvinocaridid were significantly longer (p<0.01) than those of the other caridaen. The alignment of the alvinocaridid CRs revealed two conserved sequence blocks (CSBs), and each of the CSBs included a noncanonical open reading frame (ORF), which may be involved in adjusting mitochondrial energy metabolism to adapt to the hydrothermal environment. Phylogenetic analysis supported that the deep-sea hydrothermal vent shrimps may have originated from those living in shallow area. Positive selection analysis reveals the evidence of adaptive change in the mitogenome of Alvinocarididae. Thirty potentially important adaptive residues were identified, which were located in atp6, cox1, cox3, cytb and nad1-5. This study explores the mitochondrial genetic basis of hydrothermal vent adaptation in alvinocaridid for the first time, and provides valuable clues regarding the adaptation.

The Multipartite Mitochondrial Genome of Enteromyxum leei (Myxozoa): Eight Fast-Evolving Megacircles

Myxozoans are a large group of poorly characterized cnidarian parasites. To gain further insight into their evolution, we sequenced the mitochondrial (mt) genome of Enteromyxum leei and reevaluate the mt genome structure of Kudoa iwatai. Although the typical animal mt genome is a compact, 13-25 kb, circular chromosome, the mt genome of E. leei was found to be fragmented into eight circular chromosomes of ∼23 kb, making it the largest described animal mt genome. Each chromosome was found to harbor a large noncoding region (∼15 kb), nearly identical between chromosomes. The protein coding genes show an unusually high rate of sequence evolution and possess little similarity to their cnidarian homologs. Only five protein coding genes could be identified and no tRNA genes. Surprisingly, the mt genome of K. iwatai was also found to be composed of two chromosomes. These observations confirm the remarkable plasticity of myxozoan mt genomes.

Mitochondrial RNA processing in absence of tRNA punctuations in octocorals

BACKGROUND

Mitogenome diversity is staggering among early branching animals with respect to size, gene density, content and order, and number of tRNA genes, especially in cnidarians. This last point is of special interest as tRNA cleavage drives the maturation of mitochondrial mRNAs and is a primary mechanism for mt-RNA processing in animals. Mitochondrial RNA processing in non-bilaterian metazoans, some of which possess a single tRNA gene in their mitogenomes, is essentially unstudied despite its importance in understanding the evolution of mitochondrial transcription in animals.

RESULTS

We characterized the mature mitochondrial mRNA transcripts in a species of the octocoral genus Sinularia (Alcyoniidae: Octocorallia), and defined precise boundaries of transcription units using different molecular methods. Most mt-mRNAs were polycistronic units containing two or three genes and 5' and/or 3' untranslated regions of varied length. The octocoral specific, mtDNA-encoded mismatch repair gene, the mtMutS, was found to undergo alternative polyadenylation, and exhibited differential expression of alternate transcripts suggesting a unique regulatory mechanism for this gene. In addition, a long noncoding RNA complementary to the ATP6 gene (lncATP6) potentially involved in antisense regulation was detected.

CONCLUSIONS

Mt-mRNA processing in octocorals possessing a single mt-tRNA is complex. Considering the variety of mitogenome arrangements known in cnidarians, and in general among non-bilaterian metazoans, our findings provide a first glimpse into the complex mtDNA transcription, mt-mRNA processing, and regulation among early branching animals and represent a first step towards understanding its functional and evolutionary implications.

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.

Loss and Gain of Group I Introns in the Mitochondrial Cox1 Gene of the Scleractinia (Cnidaria; Anthozoa)

Yaoyang Chuang, Marcelo Kitahara, Hironobu Fukami, Dianne Tracey, David J. Miller, and Chaolun Allen Chen (2017) Group I introns encoding a homing endonuclease gene (HEG) that is potentially capable of sponsoring mobility are present in the cytochrome oxidase subunit 1 (cox1) gene of some Hexacorallia, including a number of scleractinians assigned to the "robust" coral clade. In an e ort to infer the evolutionary history of this cox1 group I intron, DNA sequences were determined for 12 representative "basal" and "complex" corals and for 11 members of the Corallimorpharia, a sister order of the Scleractinia. Comparisons of insertion sites, secondary structures, and amino acid sequences of the HEG implied a common origin for cox1 introns of corallimorpharians, and basal and complex corals, but cox1 introns of robust corals were highly divergent, most likely reflecting independent acquisition. Phylogenetic analyses with a calibrated molecular clock suggested that cox1 introns of scleractinians and corallimorpharians have persisted at the same insertion site as that in the common ancestor 552 million years ago (mya). This ancestral intron was probably lost in complex corals around 213 to 190 mya at the junction between the Trassic and Jurassic. The coral cox1 gene remained intronless until new introns, probably from sponges or fungi, reinvaded different positions of the cox1 gene in robust corals around 135 mya in the Cretaceous, and then it subsequently began to lose them around 65.5 mya in some robust coral lineages coincident with the later Maastrichtian extinction at the Cretaceous-Tertiary boundary.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Evolutionary and biogeographical implications of degraded LAGLIDADG endonuclease functionality and group I intron occurrence in stony corals (Scleractinia) and mushroom corals (Corallimorpharia)

Group I introns and homing endonuclease genes (HEGs) are mobile genetic elements, capable of invading target sequences in intron-less genomes. LAGLIDADG HEGs are the largest family of endonucleases, playing a key role in the mobility of group I introns in a process known as ‘homing’. Group I introns and HEGs are rare in metazoans, and can be mainly found inserted in the COXI gene of some sponges and cnidarians, including stony corals (Scleractinia) and mushroom corals (Corallimorpharia). Vertical and horizontal intron transfer mechanisms have been proposed as explanations for intron occurrence in cnidarians. However, the central role of LAGLIDADG motifs in intron mobility mechanisms remains poorly understood. To resolve questions regarding the evolutionary origin and distribution of group I introns and HEGs in Scleractinia and Corallimorpharia, we examined intron/HEGs sequences within a comprehensive phylogenetic framework. Analyses of LAGLIDADG motif conservation showed a high degree of degradation in complex Scleractinia and Corallimorpharia. Moreover, the two motifs lack the respective acidic residues necessary for metal-ion binding and catalysis, potentially impairing horizontal intron mobility. In contrast, both motifs are highly conserved within robust Scleractinia, indicating a fully functional endonuclease capable of promoting horizontal intron transference. A higher rate of non-synonymous substitutions (K_a) detected in the HEGs of complex Scleractinia and Corallimorpharia suggests degradation of the HEG, whereas lower K_a rates in robust Scleractinia are consistent with a scenario of purifying selection. Molecular-clock analyses and ancestral inference of intron type indicated an earlier intron insertion in complex Scleractinia and Corallimorpharia in comparison to robust Scleractinia. These findings suggest that the lack of horizontal intron transfers in the former two groups is related to an age-dependent degradation of the endonuclease activity. Moreover, they also explain the peculiar geographical patterns of introns in stony and mushroom corals.

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.

Genome-wide polymorphism and signatures of selection in the symbiotic sea anemone Aiptasia

Background

Coral reef ecosystems are declining in response to global climate change and anthropogenic impacts. Yet patterns of standing genetic variation within cnidarian species, a major determinant of adaptive potential, are virtually unknown at genome-scale resolution. We explore patterns of genome-wide polymorphism and identify candidate loci under selection in the sea anemone Aiptasia, an important laboratory model system for studying the symbiosis between corals and dinoflagellate algae of the genus Symbiodinium.

Results

Low coverage genome sequencing revealed large genetic distances among globally widespread lineages, novel candidate targets of selection, and considerably higher heterozygosity than previously reported for Aiptasia. More than 670,000 single nucleotide polymorphisms were identified among 10 Aiptasia individuals including two pairs of genetic clones. Evolutionary relationships based on genome-wide polymorphism supported the current paradigm of a genetically distinct population from the US South Atlantic that harbors diverse Symbiodinium clades. However, anemones from the US South Atlantic demonstrated a striking lack of shared derived polymorphism. Heterozygosity was an important feature shaping nucleotide diversity patterns: at any given SNP site, more than a third of individuals genotyped were heterozygotes, and heterozygosity within individual genomes ranged from 0.37–0.58 %. Analysis of nonsynonymous and synonymous sites suggested that highly heterozygous regions are evolving under relaxed purifying selection compared to the rest of the Aiptasia genome. Genes previously identified as having elevated evolutionary rates in Aiptasia compared to other cnidarians were found in our study to be under strong purifying selection within Aiptasia. Candidate targets of selection, including lectins and genes involved in Rho GTPase signalling, were identified based on unusual signatures of nucleotide diversity, Tajima’s D, and heterozygosity compared to genome-wide averages.

Conclusions

This study represents the first genome-wide analysis of Tajima’s D in a cnidarian. Our results shed light on patterns of intraspecific genome-wide polymorphism in a model for studies of coral-algae symbiosis and present genetic targets for future research on evolutionary and cellular processes in early-diverging metazoans.

Electronic supplementary material

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

Mitochondrial group I and group II introns in the sponge orders Agelasida and Axinellida

Background

Self-splicing introns are present in the mitochondria of members of most eukaryotic lineages. They are divided into Group I and Group II introns, according to their secondary structure and splicing mechanism. Being rare in animals, self-splicing introns were only described in a few sponges, cnidarians, placozoans and one annelid species. In sponges, three types of mitochondrial Group I introns were previously described in two demosponge families (Tetillidae, and Aplysinellidae) and in the homoscleromorph family Plakinidae. These three introns differ in their insertion site, secondary structure and in the sequence of the LAGLIDADG gene they encode. Notably, no group II introns have been previously described in sponges.

Results

We report here the presence of mitochondrial introns in the cytochrome oxidase subunit 1 (COI) gene of three additional sponge species from three different families: Agelas oroides (Agelasidae, Agelasida), Cymbaxinella^pverrucosa (Hymerhabdiidae, Agelasida) and Axinella polypoides (Axinellidae, Axinellida). We show, for the first time, that sponges can also harbour Group II introns in their COI gene, whose presence in animals’ mitochondria has so far been described in only two phyla, Placozoa and Annelida. Surprisingly, two different Group II introns were discovered in the COI gene of C. verrucosa. Phylogenetic analysis indicates that the Group II introns present in C. verrucosa are related to red algae (Rhodophyta) introns.

Conclusions

The differences found among intron secondary structures and the phylogenetic inferences support the hypothesis that the introns originated from independent horizontal gene transfer events. Our results thus suggest that self-splicing introns are more diverse in the mitochondrial genome of sponges than previously anticipated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0556-1) contains supplementary material, which is available to authorized users.

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.