Evolution of linear mitochondrial genomes in medusozoan cnidarians

Characterization and phylogenetic analysis of the complete mitochondrial genome of Cotylorhiza tuberculata assembled using next-generation sequencing

In this study, the complete mitochondrial genome (mitogenome) of Cotylorhiza tuberculata (Scyphozoa; Rhizostomeae; Cepheidae) was assembled by the next-generation sequencing data. The complete mitogenome spanned 16,590 bp and contained 14 protein-coding genes, two transfer RNA genes, and two ribosomal RNA genes. Total AT% content was 67.7%, comprising A 30.22%, C 16.16%, G 17.05%, and T 36.56%. The gene arrangement exhibited consistency with the known mitogenomes of other jellyfish species. Furthermore, the phylogenetic relationship of C. tuberculata was investigated based on analysis of the 13 common protein-coding genes. Results indicated a close relationship between C. tuberculata and both Cassiopea xamachana and Cassiopea andromeda. These findings provide a valuable reference for advancing understanding of the phylogenetic relationships, taxonomic classification, and phylogeography of jellyfish species.

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

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

Dinoflagellate chloroplasts as a model for extreme genome reduction and fragmentation in organelles - The COCOA principle for gene retention

The genomes of peridinin-containing dinoflagellate chloroplasts have a very unusual organisation. These genomes are highly fragmented and greatly reduced, with most of the usual complement of chloroplast genes relocated to the nucleus. Dinoflagellate chloroplasts highlight evolutionary changes that are found to varying extents in a number of other organelle genomes. These include the chloroplast genome of the green alga Boodlea and other Cladophorales, and the mitochondrial genomes of blood-sucking and chewing lice, the parasitic plant Rhopalocnemis phalloides, the red alga Rhodosorus marinus and other members of the Stylonematophyceae, diplonemid flagellates, and some Cnidaria. Consideration of the coding content of the remnant chloroplast genomes indicates that organelles may preferentially retain genes for proteins important in initiating assembly of complexes, and the same is largely true for mitochondria. We propose a new principle, of CO-location for COntrol of Assembly (COCOA), indicating the importance of retaining these genes in the organelle. This adds to, but does not invalidate, the existing hypotheses of the multisubunit completion principle, CO-location for Redox Regulation (CORR) and Control by Epistasy of Synthesis (CES).

The genome of the colonial hydroid Hydractinia reveals that their stem cells use a toolkit of evolutionarily shared genes with all animals

Hydractinia is a colonial marine hydroid that shows remarkable biological properties, including the capacity to regenerate its entire body throughout its lifetime, a process made possible by its adult migratory stem cells, known as i-cells. Here, we provide an in-depth characterization of the genomic structure and gene content of two Hydractinia species, Hydractinia symbiolongicarpus and Hydractinia echinata, placing them in a comparative evolutionary framework with other cnidarian genomes. We also generated and annotated a single-cell transcriptomic atlas for adult male H. symbiolongicarpus and identified cell-type markers for all major cell types, including key i-cell markers. Orthology analyses based on the markers revealed that Hydractinia's i-cells are highly enriched in genes that are widely shared amongst animals, a striking finding given that Hydractinia has a higher proportion of phylum-specific genes than any of the other 41 animals in our orthology analysis. These results indicate that Hydractinia's stem cells and early progenitor cells may use a toolkit shared with all animals, making it a promising model organism for future exploration of stem cell biology and regenerative medicine. The genomic and transcriptomic resources for Hydractinia presented here will enable further studies of their regenerative capacity, colonial morphology, and ability to distinguish self from nonself.

Evolution of myxozoan mitochondrial genomes: insights from myxobolids

BACKGROUND

Myxozoa is a class of cnidarian parasites that encompasses over 2,400 species. Phylogenetic relationships among myxozoans remain highly debated, owing to both a lack of informative morphological characters and a shortage of molecular markers. Mitochondrial (mt) genomes are a common marker in phylogeny and biogeography. However, only five complete myxozoan mt genomes have been sequenced: four belonging to two closely related genera, Enteromyxum and Kudoa, and one from the genus Myxobolus. Interestingly, while cytochrome oxidase genes could be identified in Enteromyxum and Kudoa, no such genes were found in Myxobolus squamalis, and another member of the Myxobolidae (Henneguya salminicola) was found to have lost its entire mt genome. To evaluate the utility of mt genomes to reconstruct myxozoan relationships and to understand if the loss of cytochrome oxidase genes is a characteristic of myxobolids, we sequenced the mt genome of five myxozoans (Myxobolus wulii, M. honghuensis, M. shantungensis, Thelohanellus kitauei and, Sphaeromyxa zaharoni) using Illumina and Oxford Nanopore platforms.

RESULTS

Unlike Enteromyxum, which possesses a partitioned mt genome, the five mt genomes were encoded on single circular chromosomes. An mt plasmid was found in M. wulii, as described previously in Kudoa iwatai. In all new myxozoan genomes, five protein-coding genes (cob, cox1, cox2, nad1, and nad5) and two rRNAs (rnl and rns) were recognized, but no tRNA. We found that Myxobolus and Thelohanellus species shared unidentified reading frames, supporting the view that these mt open reading frames are functional. Our phylogenetic reconstructions based on the five conserved mt genes agree with previously published trees based on the 18S rRNA gene.

CONCLUSIONS

Our results suggest that the loss of cytochrome oxidase genes is not a characteristic of all myxobolids, the ancestral myxozoan mt genome was likely encoded on a single circular chromosome, and mt plasmids exist in a few lineages. Our findings indicate that myxozoan mt sequences are poor markers for reconstructing myxozoan phylogenetic relationships because of their fast-evolutionary rates and the abundance of repeated elements, which complicates assembly.

The mitochondrial genome of Monoserius pennarius (Linnaeus, 1758) from the East China Sea

The species Monoserius pennarius (Linnaeus, 1758), is particularly abundant in the tropical Indo-West Pacific east of Sri Lanka, yet very limited genetic information exists for this species. Here, we report the assembled-linear mitochondrial genome of M. pennarius collected from the East China Sea. The 15,197 bp mitogenome contains 13 protein-coding genes (PCGs), two tRNA genes, and two rRNA genes. Notably, the gene order in this mitogenome differs from that of other hydrozoans within the same taxonomic order. Phylogenetic analysis, based on 13 concatenated mitochondrial PCGs, recovered M. pennarius as a sister of Nemalecium lighti (Hargitt, 1924), outside the other Leptothecata hydrozoans, suggesting paraphyly of Leptothecata. The mitogenome of M. pennarius, serving as the first publicly available for the family Aglaopheniidae, holds foreseeable value for investigating Leptothecata evolution.

Giants among Cnidaria: Large Nuclear Genomes and Rearranged Mitochondrial Genomes in Siphonophores

Siphonophores (Cnidaria: Hydrozoa) are abundant predators found throughout the ocean and are important constituents of the global zooplankton community. They range in length from a few centimeters to tens of meters. They are gelatinous, fragile, and difficult to collect, so many aspects of the biology of these roughly 200 species remain poorly understood. To survey siphonophore genome diversity, we performed Illumina sequencing of 32 species sampled broadly across the phylogeny. Sequencing depth was sufficient to estimate nuclear genome size from k-mer spectra in six specimens, ranging from 0.7 to 2.3 Gb, with heterozygosity estimates between 0.69% and 2.32%. Incremental k-mer counting indicates k-mer peaks can be absent with nearly 20× read coverage, suggesting minimum genome sizes range from 1.4 to 5.6 Gb in the 25 samples without peaks in the k-mer spectra. This work confirms most siphonophore nuclear genomes are large relative to the genomes of other cnidarians, but also identifies several with reduced size that are tractable targets for future siphonophore nuclear genome assembly projects. We also assembled complete mitochondrial genomes for 33 specimens from these new data, indicating a conserved gene order shared among nonsiphonophore hydrozoans, Cystonectae, and some Physonectae, revealing the ancestral mitochondrial gene order of siphonophores. Our results also suggest extensive rearrangement of mitochondrial genomes within other Physonectae and in Calycophorae. Though siphonophores comprise a small fraction of cnidarian species, this survey greatly expands our understanding of cnidarian genome diversity. This study further illustrates both the importance of deep phylogenetic sampling and the utility of k-mer-based genome skimming in understanding the genomic diversity of a clade.

Comparative mitogenomics of Brachiopods reveals conservatism in articulate species and unusualness in inarticulate species

BACKGROUND

Brachiopods are a phylum of marine invertebrates with over 10,000 fossil species. Today, there are fewer than 500 extant species assigned to the class Articulata or Inarticulata and for which knowledge of evolutionary genetics and genomics is still poor. Until now, complete mitogenome sequences of two inarticulate species and four articulate species were available.

METHODS AND RESULTS

The complete mitogenome of the inarticulate brachiopod species Lingula reevii (20,778 bp) was obtained by using next generation sequencing. It contains 12 protein-coding genes (the annotation of atp8 is unsure), two ribosomal RNA genes, 26 transfer RNA genes, and one supernumerary ORF that is also conserved in the inarticulate species Lingula anatina. It is hypothesized that this ORF could represent a Lingula-specific mtORFan gene (without obvious homology to other genes). Comparative mitogenomics indicate the mitochondrial gene order of L. reevii is unique among brachiopods, and that compared to articulate species, inarticulate species exhibit massive mitogenome rearrangements, deviant ATP8 protein sequences and supernumerary ORFs, possibly representing species- or lineage-specific mtORFan genes.

CONCLUSION

The results of this study enrich genetics knowledge of extant brachiopods, which may eventually help to test hypotheses about their decline.

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

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

The complete mitochondrial genome of the mauve stinger jellyfish Pelagia noctiluca Forskål, 1775 (Cnidaria, Scyphozoa, Semaeostomeae) with phylogenetic analysis

This study determined the complete mitochondrial genome of the jellyfish Pelagia noctiluca (Scyphozoa, Semaeostomeae) for the first time. The genome was a linear molecule of 16,390 bp in length and 59.3% AT. It comprised of 13 typical protein-coding genes (cox1-3, nd1-6, nd4L, atp6, atp8, and cytB), two ribosomal RNAs (16S and 12S rRNA), and two tRNAs (trnM and trnW). In addition, we detected two additional open reading frames (polB and ORF314) at one end of the genome. The gene-coding structures were identical to those of other scyphozoans. Based on a molecular phylogeny constructed using 13 protein-coding genes, P. noctiluca has the closest genetic relationship with the genus Chrysaora (Semaeostomeae).

Complete mitochondrial genome of Mastigias papua (Scyphozoa: Rhizostomeae: Mastigiidae) based on next-generation sequencing and phylogenetic analysis

We assembled the complete mitochondrial genome (mitogenome) of Mastigias papua (Scyphozoa: Rhizostomeae: Mastigiidae) by the data generated from the next-generation sequencing platform. The complete mitogenome of M. papua was 16,560 bp in length, containing 14 protein-coding genes, two transfer RNA genes, and two ribosomal RNA genes. The base compositions were A 30.65%, C 15.16%, G 16.34%, and T 37.86%, with a gene arrangement similar to the mitogenomes derived from other representatives of Scyphozoa. Based on the 13 common protein-coding genes of 16 species within Scyphozoa, we constructed the phylogenetic tree and found that M. papua has a close relationship with Cassiopea andromeda and Cassiopea xamachana. All these species belong to an order of jellyfish Rhizostomeae, which have similar morphological characteristics. This is agreement with the conclusion we got by the phylogenetic relationship analysis using molecular data. This research has practical implications for advancing understanding of the phylogenetic relationships, taxonomic classifications, and phylogeography within Scyphozoa.

Revisiting mitogenome evolution in Medusozoa with eight new mitochondrial genomes

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

Description of a common stauromedusa on the Pacific Coast of the United States and Canada, Haliclystus sanjuanensis new species (Cnidaria: Staurozoa)

Haliclystus "sanjuanensis" nomen nudum is the most common staurozoan on the west coast of the United States and Canada. This species was described in the M.S. Thesis by Gellermann (1926) and although that name has been in use nearly continuously since that time, no published description exists. Furthermore, the most popular operative name for this species has varied between several related species names over time, resulting in confusion. Herein, we provide a detailed description and synonymy of Haliclystus sanjuanensis n. sp., whose distribution is verified from Unalaska Island in the Aleutians (53.4° N, 166.8° W) in the northwest, to Santa Barbara County, California, just north of Point Conception (34.5° N, 120.5° W), in the south. Haliclystus sanjuanensis n. sp. is compared with the twelve other described species of Haliclystus and illustrations of both macroscopic and microscopic anatomy are provided. Haliclystus sanjuanensis n. sp. is unique among species of Haliclystus in the arrangement of the bright-white nematocyst spots in its calyx and the pattern of dark stripes running the length of the stalk and up the outside of the calyx.

The genome of the colonial hydroid Hydractinia reveals their stem cells utilize a toolkit of evolutionarily shared genes with all animals

Hydractinia is a colonial marine hydroid that exhibits remarkable biological properties, including the capacity to regenerate its entire body throughout its lifetime, a process made possible by its adult migratory stem cells, known as i-cells. Here, we provide an in-depth characterization of the genomic structure and gene content of two Hydractinia species, H. symbiolongicarpus and H. echinata, placing them in a comparative evolutionary framework with other cnidarian genomes. We also generated and annotated a single-cell transcriptomic atlas for adult male H. symbiolongicarpus and identified cell type markers for all major cell types, including key i-cell markers. Orthology analyses based on the markers revealed that Hydractinia's i-cells are highly enriched in genes that are widely shared amongst animals, a striking finding given that Hydractinia has a higher proportion of phylum-specific genes than any of the other 41 animals in our orthology analysis. These results indicate that Hydractinia's stem cells and early progenitor cells may use a toolkit shared with all animals, making it a promising model organism for future exploration of stem cell biology and regenerative medicine. The genomic and transcriptomic resources for Hydractinia presented here will enable further studies of their regenerative capacity, colonial morphology, and ability to distinguish self from non-self.

Mitochondrial genome comparison reveals the evolution of cnidarians

Cnidarians are the most primitive metazoans, but their evolutionary relationships are poorly understood, although recent studies present several phylogenetic hypotheses. Here, we collected 266 complete cnidarian mitochondrial genomes and re-evaluated the phylogenetic relationships between the major lineages. We described the gene rearrangement patterns of Cnidaria. Anthozoans had significantly greater mitochondrial genome size and lower A + T content than medusozoans. Most of the protein-coding genes in anthozoans such as COX 13, ATP6, and CYTB displayed a faster rate of evolution based on selection analysis. There were 19 distinct patterns of mitochondrial gene order, including 16 unique gene orders in anthozoans and 3 mtDNA gene orders pattern in medusozoans, were identified among cnidarians. The gene order arrangement suggested that a linearized mtDNA structure may be more conducive to Medusozoan mtDNA stability. Based on phylogenetic analyses, the monophyly of the Anthozoa was strongly supported compared to previous mitochondrial genome-based analyses rather than octocorals forming a sister group relationship with medusozoans. In addition, Staurozoa were more closely related to Anthozoa than to Medusozoa. In conclusion, these results largely support the traditional phylogenetic view of the relationships of cnidarians and provide new insights into the evolutionary processes for studying the most ancient animal radiations.

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.

The metazoan landscape of mitochondrial DNA gene order and content is shaped by selection and affects mitochondrial transcription

Mitochondrial DNA (mtDNA) harbors essential genes in most metazoans, yet the regulatory impact of the multiple evolutionary mtDNA rearrangements has been overlooked. Here, by analyzing mtDNAs from ~8000 metazoans we found high gene content conservation (especially of protein and rRNA genes), and codon preferences for mtDNA-encoded tRNAs across most metazoans. In contrast, mtDNA gene order (MGO) was selectively constrained within but not between phyla, yet certain gene stretches (ATP8-ATP6, ND4-ND4L) were highly conserved across metazoans. Since certain metazoans with different MGOs diverge in mtDNA transcription, we hypothesized that evolutionary mtDNA rearrangements affected mtDNA transcriptional patterns. As a first step to test this hypothesis, we analyzed available RNA-seq data from 53 metazoans. Since polycistron mtDNA transcripts constitute a small fraction of the steady-state RNA, we enriched for polycistronic boundaries by calculating RNA-seq read densities across junctions between gene couples encoded either by the same strand (SSJ) or by different strands (DSJ). We found that organisms whose mtDNA is organized in alternating reverse-strand/forward-strand gene blocks (mostly arthropods), displayed significantly reduced DSJ read counts, in contrast to organisms whose mtDNA genes are preferentially encoded by one strand (all chordates). Our findings suggest that mtDNA rearrangements are selectively constrained and likely impact mtDNA regulation.

Myxozoans (Cnidaria) do not Retain Key Oxygen-Sensing and Homeostasis Toolkit Genes

For aerobic organisms, both the hypoxia-inducible factor pathway and the mitochondrial genomes are key players in regulating oxygen homeostasis. Recent work has suggested that these mechanisms are not as highly conserved as previously thought, prompting more surveys across animal taxonomic levels, which would permit testing of hypotheses about the ecological conditions facilitating evolutionary loss of such genes. The Phylum Cnidaria is known to harbor wide variation in mitochondrial chromosome morphology, including an extreme example, in the Myxozoa, of mitochondrial genome loss. Because myxozoans are obligate endoparasites, frequently encountering hypoxic environments, we hypothesize that variation in environmental oxygen availability could be a key determinant in the evolution of metabolic gene networks associated with oxygen-sensing, hypoxia-response, and energy production. Here, we surveyed genomes and transcriptomes across 46 cnidarian species for the presence of HIF pathway members, as well as for an assortment of hypoxia, mitochondrial, and stress-response toolkit genes. We find that presence of the HIF pathway, as well as number of genes associated with mitochondria, hypoxia, and stress response, do not vary in parallel to mitochondrial genome morphology. More interestingly, we uncover evidence that myxozoans have lost the canonical HIF pathway repression machinery, potentially altering HIF pathway functionality to work under the specific conditions of their parasitic lifestyles. In addition, relative to other cnidarians, myxozoans show loss of large proportions of genes associated with the mitochondrion and involved in response to hypoxia and general stress. Our results provide additional evidence that the HIF regulatory machinery is evolutionarily labile and that variations in the canonical system have evolved in many animal groups.

Evolution of mitochondrial and nuclear genomes in Pennatulacea

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

The Phylogenetic Position of the Enigmatic, Polypodium hydriforme (Cnidaria, Polypodiozoa): Insights from Mitochondrial Genomes

Polypodium hydriforme is an enigmatic parasite that belongs to the phylum Cnidaria. Its taxonomic position has been debated: whereas it was previously suggested to be part of Medusozoa, recent phylogenomic analyses based on nuclear genes support the view that P. hydriforme and Myxozoa form a clade called Endocnidozoa. Medusozoans have linear mitochondrial (mt) chromosomes, whereas myxozoans, as most metazoan species, have circular chromosomes. In this work, we determined the structure of the mt genome of P. hydriforme, using Illumina and Oxford Nanopore Technologies reads, and showed that it is circular. This suggests that P. hydriforme is not nested within Medusozoa, as this would entail linearization followed by recirculation. Instead, our results support the view that P. hydriforme is a sister clade to Myxozoa, and mt linearization in the lineage leading to medusozoans occurred after the divergence of Myxozoa + P. hydriforme. Detailed analyses of the assembled P. hydriforme mt genome show that: (1) it is encoded on a single circular chromosome with an estimated size of ∼93,000 base pairs, making it one of the largest metazoan mt genomes; (2) around 78% of the genome encompasses a noncoding region composed of several repeat types; (3) similar to Myxozoa, no mt tRNAs were identified; (4) the codon TGA is a stop codon and does not encode for tryptophan as in other cnidarians; (5) similar to myxozoan mt genomes, it is extremely fast evolving.

Caryophylliids (Anthozoa, Scleractinia) and mitochondrial gene order: insights from mitochondrial and nuclear phylogenomics

Molecularly, the family Caryophylliidae is polyphyletic and different sets of genetic data converge towards a consensus that a taxonomic review of this family is necessary. Overall, the order of genes in the mitochondrial genome (mitogenome) together with DNA sequences have been used to successfully untangle evolutionary relationships in several groups of organisms. Published mitogenomes of two caryophylliid genera (Desmophyllum and Solenosmilia) present a transposition of the gene block containing cob, nad2, and nad6, which is located between nad5 5' exon and trnW, while that of Polycyathus chaishanensis presents the same gene order as the majority of scleractinian corals. In molecular-based evolutionary reconstructions, caryophylliids that have the mitochondrial gene rearrangement were recovered as a monophyletic lineage ("true" caryophylliids), while members of the genus Polycyathus were placed in a different position. In this study, additional mitogenomes of this family were assembled and included in evolutionary reconstructions of Scleractinia in order to improve our understanding on whether the mitogenome gene rearrangement is limited to and, therefore, could be a synapomorphy of the actual members of Caryophylliidae. Specimens of Caryophyllia scobinosa, Premocyathus sp., Heterocyathus sulcatus, and Trochocyathus caryophylloides, as well as Desmophyllum pertusum and Solenosmilia variabilis from the Southwest Atlantic were sequenced using Illumina platforms. Then, mitochondrial genomes were assembled and annotated, and nuclear datasets were recovered in-silico from assembled contigs using a previously published set of baits. Evolutionary reconstructions were performed using mitochondrial and nuclear datasets and based on Maximum Likelihood and Bayesian Inference. Obtained mitogenomes are circular and range between 15,816 and 18,225 bp in size and from 30.76% to 36.63% in GC content. The gene rearrangement is only seen in C. scobinosa, D. pertusum, Premocyathus sp., and S. variabilis, which were recovered as a monophyletic clade in both mitochondrial and nuclear phylogenies. On the other hand, the "caryophylliids" with the canonical mitogenome gene order were not recovered within this clade. Differences in features of the skeleton of "true" caryophylliids in comparison to traditional members of the family were observed and offer further support that the gene rearrangement might be seen as a synapomorphy of family Caryophylliidae.

Fragmentation in mitochondrial genomes in relation to elevated sequence divergence and extreme rearrangements

Background

A single circular mitochondrial (mt) genome is a common feature across most metazoans. The mt-genome includes protein-coding genes involved in oxidative phosphorylation, as well as RNAs necessary for translation of mt-RNAs, whose order and number are highly conserved across animal clades, with few known exceptions of alternative mt-gene order or mt-genome architectures. One such exception consists of the fragmented mitochondrial genome, a type of genome architecture where mt-genes are split across two or more mt-chromosomes. However, the origins of mt-genome fragmentation and its effects on mt-genome evolution are unknown. Here, we investigate these origin and potential mechanisms underlying mt-genome fragmentation, focusing on a genus of booklice, Liposcelis, which exhibits elevated sequence divergence, frequent rearrangement of mt-gene order, and fragmentation of the mt genome, and compare them to other Metazoan clades.

Results

We found this genus Liposcelis exhibits very low conservation of mt-gene order across species, relative to other metazoans. Levels of gene order rearrangement were, however, unrelated to whether or not mt-genomes were fragmented or intact, suggesting mitochondrial genome fragmentation is not affecting mt-gene order directly. We further investigated possible mechanisms underpinning these patterns and revealed very high conservation of non-coding sequences at the edges of multiple recombination regions across populations of one particular Liposcelis species, supportive of a hypothesis that mt-fragmentation arises from recombination errors between mt-genome copies. We propose these errors may arise as a consequence of a heightened mutation rate in clades exhibiting mt-fragmentation. Consistent with this, we observed a striking pattern across three Metazoan phyla (Arthropoda, Nematoda, Cnidaria) characterised by members exhibiting high levels of mt-gene order rearrangement and cases of mt-fragmentation, whereby the mt-genomes of species more closely related to species with fragmented mt-genomes diverge more rapidly despite experiencing strong purifying selection.

Conclusions

We showed that contrary to expectations, mt-genome fragmentation is not correlated with the increase in mt-genome rearrangements. Furthermore, we present evidence that fragmentation of the mt-genome may be part of a general relaxation of a natural selection on the mt-genome, thus providing new insights into the origins of mt-genome fragmentation and evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01218-7.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

The mitochondrial genome of Nemalecium lighti (Hydrozoa, Leptothecata)

The hydrozoan species Nemalecium lighti (Hargitt, 1924) is widely distributed in tropical marine waters around the world. Here we report the complete linear mitochondrial genome of N. lighti from Sint Eustatius (Lesser Antilles). The mitochondrial genome with a length of 14,320 bp encodes for 13 protein-coding genes, two tRNA genes, and two rRNA genes. Gene arrangement differs from that found in other species of the same taxonomic order and a phylogenetic analysis shows that based on mitochondrial genes, N. lighti clusters outside of the Leptothecata, rendering the order paraphyletic.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

Complete mitochondrial genome of the hydrozoan jellyfish Blackfordia virginica Mayer, 1910 (Cnidaria; Hydrozoa; Leptothecata) with phylogenetic analysis

In this study, we analyzed the complete mitochondrial genome of the hydrozoan jellyfish Blackfordia virginica. The genome was a linear form (15,109 bp long, 73.6% AT), including 13 protein-coding genes (cox2, atp8, atp6, cox3, nad2, nad5, nad6, nad3, nad4L, nad1, nad4, cytB, and cox1), 2 tRNAs (tRNA-Met and tRNA-Trp), and 2 rRNAs (12S and 16S RNA). The genome structure of the B. virginica was completely identical to mitochondrial genomes of other hydrozoans that belonged to Leptothecata and Anthoathecata. Molecular phylogenetic analysis within hydrozoan species showed that B. virginica was the closest to the hydrozoan Laomedea flexuosa.

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

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

Advances and Discoveries in Myxozoan Genomics

Myxozoans are highly diverse and globally distributed cnidarian endoparasites in freshwater and marine habitats. They have adopted a heteroxenous life cycle, including invertebrate and fish hosts, and have been associated with diseases in aquaculture and wild fish stocks. Despite their importance, genomic resources of myxozoans have proven difficult to obtain due to their miniaturized and derived genome character and close associations with fish tissues. The first 'omic' datasets have now become the main resource for a better understanding of host-parasite interactions, virulence, and diversity, but also the evolutionary history of myxozoans. In this review, we discuss recent genomic advances in the field and outline outstanding questions to be answered with continuous and improved efforts of generating myxozoan genomic data.

Mitochondrial Genome Evolution of Placozoans: Gene Rearrangements and Repeat Expansions

Placozoans, nonbilaterian animals with the simplest known metazoan bauplan, are currently classified into 20 haplotypes belonging to three genera, Polyplacotoma, Trichoplax, and Hoilungia. The latter two comprise two and five clades, respectively. In Trichoplax and Hoilungia, previous studies on six haplotypes belonging to four different clades have shown that their mtDNAs are circular chromosomes of 32–43 kb in size, which encode 12 protein-coding genes, 24 tRNAs, and two rRNAs. These mitochondrial genomes (mitogenomes) also show unique features rarely seen in other metazoans, including open reading frames (ORFs) of unknown function, and group I and II introns. Here, we report seven new mitogenomes, covering the five previously described haplotypes H2, H17, H19, H9, and H11, as well as two new haplotypes, H23 (clade III) and H24 (clade VII). The overall gene content is shared between all placozoan mitochondrial genomes, but genome sizes, gene orders, and several exon–intron boundaries vary among clades. Phylogenomic analyses strongly support a tree topology different from previous 16S rRNA analyses, with clade VI as the sister group to all other Hoilungia clades. We found small inverted repeats in all 13 mitochondrial genomes of the Trichoplax and Hoilungia genera and evaluated their distribution patterns among haplotypes. Because Polyplacotoma mediterranea (H0), the sister to the remaining haplotypes, has a small mitochondrial genome with few small inverted repeats and ORFs, we hypothesized that the proliferation of inverted repeats and ORFs substantially contributed to the observed increase in the size and GC content of the Trichoplax and Hoilungia mitochondrial genomes.

Revisiting Ceriantharian (Anthozoa) Mitochondrial Genomes: Casting Doubts about Their Structure and Size

Recently, Stampar et al. (2019. Linear mitochondrial genome in Anthozoa (Cnidaria): a case study in. Sci Rep. 9(1):6094.) uncovered highly atypical mitochondrial genome structures in the cnidarian species Pachycerianthus magnus and Isarachnanthus nocturnus (Anthozoa, Ceriantharia). These two mitochondrial DNAs assembled as linear fragmented genomes, comprising eight and five chromosomes, respectively—architectures unlike any other anthozoan mitogenome described to date. What’s more, they have cumulative lengths of 77.8 (P. magnus) and 80.9 kb (I. nocturnus), making them the largest animal mitochondrial DNAs on record, a finding which garnered significant attention by various news media. Here, I take a closer look at the work of Stampar et al. and question their key results. I provide evidence that the currently available mitogenome sequences for I. nocturnus and P. magnus, including their structures, sizes, and chromosome numbers, should be treated with caution. More work must be done on these genomes before one can say with any certainty that they are linear, fragmented, or the largest animal mitogenomes observed to date.

Mitochondrial genome of the nonphotosynthetic mycoheterotrophic plant Hypopitys monotropa, its structure, gene expression and RNA editing

Heterotrophic plants—plants that have lost the ability to photosynthesize—are characterized by a number of changes at all levels of organization. Heterotrophic plants are divided into two large categories—parasitic and mycoheterotrophic (MHT). The question of to what extent such changes are similar in these two categories is still open. The plastid genomes of nonphotosynthetic plants are well characterized, and they exhibit similar patterns of reduction in the two groups. In contrast, little is known about the mitochondrial genomes of MHT plants. We report the structure of the mitochondrial genome of Hypopitys monotropa, a MHT member of Ericaceae, and the expression of its genes. In contrast to its highly reduced plastid genome, the mitochondrial genome of H. monotropa is larger than that of its photosynthetic relative Vaccinium macrocarpon, and its complete size is ~810 Kb. We observed an unusually long repeat-rich structure of the genome that suggests the existence of linear fragments. Despite this unique feature, the gene content of the H. monotropa mitogenome is typical of flowering plants. No acceleration of substitution rates is observed in mitochondrial genes, in contrast to previous observations in parasitic non-photosynthetic plants. Transcriptome sequencing revealed the trans-splicing of several genes and RNA editing in 33 of 38 genes. Notably, we did not find any traces of horizontal gene transfer from fungi, in contrast to plant parasites, which extensively integrate genetic material from their hosts.

Jellyfish genomes reveal distinct homeobox gene clusters and conservation of small RNA processing

The phylum Cnidaria represents a close outgroup to Bilateria and includes familiar animals including sea anemones, corals, hydroids, and jellyfish. Here we report genome sequencing and assembly for true jellyfish Sanderia malayensis and Rhopilema esculentum. The homeobox gene clusters are characterised by interdigitation of Hox, NK, and Hox-like genes revealing an alternate pathway of ANTP class gene dispersal and an intact three gene ParaHox cluster. The mitochondrial genomes are linear but, unlike in Hydra, we do not detect nuclear copies, suggesting that linear plastid genomes are not necessarily prone to integration. Genes for sesquiterpenoid hormone production, typical for arthropods, are also now found in cnidarians. Somatic and germline cells both express piwi-interacting RNAs in jellyfish revealing a conserved cnidarian feature, and evidence for tissue-specific microRNA arm switching as found in Bilateria is detected. Jellyfish genomes reveal a mosaic of conserved and divergent genomic characters evolved from a shared ancestral genetic architecture.

Recent advances in understanding mitochondrial genome diversity

Ever since its discovery, the double-stranded DNA contained in the mitochondria of eukaryotes has fascinated researchers because of its bacterial endosymbiotic origin, crucial role in encoding subunits of the respiratory complexes, compact nature, and specific inheritance mechanisms. In the last few years, high-throughput sequencing techniques have accelerated the sequencing of mitochondrial genomes (mitogenomes) and uncovered the great diversity of organizations, gene contents, and modes of replication and transcription found in living eukaryotes. Some early divergent lineages of unicellular eukaryotes retain certain synteny and gene content resembling those observed in the genomes of alphaproteobacteria (the inferred closest living group of mitochondria), whereas others adapted to anaerobic environments have drastically reduced or even lost the mitogenome. In the three main multicellular lineages of eukaryotes, mitogenomes have pursued diverse evolutionary trajectories in which different types of molecules (circular versus linear and single versus multipartite), gene structures (with or without self-splicing introns), gene contents, gene orders, genetic codes, and transfer RNA editing mechanisms have been selected. Whereas animals have evolved a rather compact mitochondrial genome between 11 and 50 Kb in length with a highly conserved gene content in bilaterians, plants exhibit large mitochondrial genomes of 66 Kb to 11.3 Mb with large intergenic repetitions prone to recombination, and fungal mitogenomes have intermediate sizes of 12 to 236 Kb.

Conserved novel ORFs in the mitochondrial genome of the ctenophore Beroe forskalii

To date, five ctenophore species' mitochondrial genomes have been sequenced, and each contains open reading frames (ORFs) that if translated have no identifiable orthologs. ORFs with no identifiable orthologs are called unidentified reading frames (URFs). If truly protein-coding, ctenophore mitochondrial URFs represent a little understood path in early-diverging metazoan mitochondrial evolution and metabolism. We sequenced and annotated the mitochondrial genomes of three individuals of the beroid ctenophore Beroe forskalii and found that in addition to sharing the same canonical mitochondrial genes as other ctenophores, the B. forskalii mitochondrial genome contains two URFs. These URFs are conserved among the three individuals but not found in other sequenced species. We developed computational tools called pauvre and cuttlery to determine the likelihood that URFs are protein coding. There is evidence that the two URFs are under negative selection, and a novel Bayesian hypothesis test of trinucleotide frequency shows that the URFs are more similar to known coding genes than noncoding intergenic sequence. Protein structure and function prediction of all ctenophore URFs suggests that they all code for transmembrane transport proteins. These findings, along with the presence of URFs in other sequenced ctenophore mitochondrial genomes, suggest that ctenophores may have uncharacterized transmembrane proteins present in their mitochondria.

Exceptional Enlargement of the Mitochondrial Genome Results from Distinct Causes in Different Rain Frogs (Anura: Brevicipitidae: Breviceps)

The mitochondrial (mt) genome of the bushveld rain frog (Breviceps adspersus, Brevicipitidae, Afrobatrachia) is the largest (28.8 kbp) among the vertebrates investigated to date. The major cause of genome size enlargement in this species is the duplication of multiple genomic regions. To investigate the evolutionary lineage, timing, and process of mt genome enlargement, we sequenced the complete mtDNAs of two congeneric rain frogs, B. mossambicus and B. poweri. The mt genomic organization, gene content, and gene arrangements of these two rain frogs are very similar to each other but differ from those of B. adspersus. The B. mossambicus mt genome (22.5 kbp) does not differ significantly from that of most other afrobatrachians. In contrast, the B. poweri mtDNA (28.1 kbp) is considerably larger: currently the second largest among vertebrates, after B. adspersus. The main causes of genome enlargement differ among Breviceps species. Unusual elongation (12.5 kbp) of the control region (CR), a single major noncoding region of the vertebrate mt genome, is responsible for the extremely large mt genome in B. poweri. Based on the current Breviceps phylogeny and estimated divergence age, it can be concluded that the genome enlargements occurred independently in each species lineage within relatively short periods. Furthermore, a high nucleotide substitution rate and relaxation of selective pressures, which are considered to be involved in changes in genome size, were also detected in afrobatrachian lineages. Our results suggest that these factors were not direct causes but may have indirectly affected mt genome enlargements in Breviceps.

Higher Order Organization of the mtDNA: Beyond Mitochondrial Transcription Factor A

The higher order organization of eukaryotic and prokaryotic genomes is pivotal in the regulation of gene expression. Specifically, chromatin accessibility in eukaryotes and nucleoid accessibility in bacteria are regulated by a cohort of proteins to alter gene expression in response to diverse physiological conditions. By contrast, prior studies have suggested that the mitochondrial genome (mtDNA) is coated solely by mitochondrial transcription factor A (TFAM), whose increased cellular concentration was proposed to be the major determinant of mtDNA packaging in the mitochondrial nucleoid. Nevertheless, recent analysis of DNase-seq and ATAC-seq experiments from multiple human and mouse samples suggest gradual increase in mtDNA occupancy during the course of embryonic development to generate a conserved footprinting pattern which correlate with sites that have low TFAM occupancy in vivo (ChIP-seq) and tend to adopt G-quadruplex structures. These findings, along with recent identification of mtDNA binding by known modulators of chromatin accessibility such as MOF, suggest that mtDNA higher order organization is generated by cross talk with the nuclear regulatory system, may have a role in mtDNA regulation, and is more complex than once thought.

Regenerative Capacity of the Upside-down Jellyfish Cassiopea xamachana

This study provides the first observation that umbrellar tissue can lead to the formation of virtually all body structures in jellyfish of the order Rhizostomeae. The regeneration process was observed in two specimens of the upside-down jellyfish Cassiopea xamachana Bigelow, 1892, one housed at the Vienna Zoo, Austria and the other in a laboratory at the University of São Paulo, Brazil. The process was triggered by an injury and ended with the formation of two new sets of body structures. Our observation offers evidence that C. xamachana has a hidden regenerative capacity exceeding that previously recorded.

Mammalian mitochondrial translation - revealing consequences of divergent evolution

Mitochondria are ubiquitous organelles present in the cytoplasm of all nucleated eukaryotic cells. These organelles are described as arising from a common ancestor but a comparison of numerous aspects of mitochondria between different organisms provides remarkable examples of divergent evolution. In humans, these organelles are of dual genetic origin, comprising ∼1500 nuclear-encoded proteins and thirteen that are encoded by the mitochondrial genome. Of the various functions that these organelles perform, it is only oxidative phosphorylation, which provides ATP as a source of chemical energy, that is dependent on synthesis of these thirteen mitochondrially encoded proteins. A prerequisite for this process of translation are the mitoribosomes. The recent revolution in cryo-electron microscopy has generated high-resolution mitoribosome structures and has undoubtedly revealed some of the most distinctive molecular aspects of the mitoribosomes from different organisms. However, we still lack a complete understanding of the mechanistic aspects of this process and many of the factors involved in post-transcriptional gene expression in mitochondria. This review reflects on the current knowledge and illustrates some of the striking differences that have been identified between mitochondria from a range of organisms.

Population Genetics of Paramecium Mitochondrial Genomes: Recombination, Mutation Spectrum, and Efficacy of Selection

The evolution of mitochondrial genomes and their population-genetic environment among unicellular eukaryotes are understudied. Ciliate mitochondrial genomes exhibit a unique combination of characteristics, including a linear organization and the presence of multiple genes with no known function or detectable homologs in other eukaryotes. Here we study the variation of ciliate mitochondrial genomes both within and across 13 highly diverged Paramecium species, including multiple species from the P. aurelia species complex, with four outgroup species: P. caudatum, P. multimicronucleatum, and two strains that may represent novel related species. We observe extraordinary conservation of gene order and protein-coding content in Paramecium mitochondria across species. In contrast, significant differences are observed in tRNA content and copy number, which is highly conserved in species belonging to the P. aurelia complex but variable among and even within the other Paramecium species. There is an increase in GC content from ∼20% to ∼40% on the branch leading to the P. aurelia complex. Patterns of polymorphism in population-genomic data and mutation-accumulation experiments suggest that the increase in GC content is primarily due to changes in the mutation spectra in the P. aurelia species. Finally, we find no evidence of recombination in Paramecium mitochondria and find that the mitochondrial genome appears to experience either similar or stronger efficacy of purifying selection than the nucleus.

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

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

Computational Characterization of the mtORF of Pocilloporid Corals: Insights into Protein Structure and Function in Stylophora Lineages from Contrasting Environments

More than a decade ago, a new mitochondrial Open Reading Frame (mtORF) was discovered in corals of the family Pocilloporidae and has been used since then as an effective barcode for these corals. Recently, mtORF sequencing revealed the existence of two differentiated Stylophora lineages occurring in sympatry along the environmental gradient of the Red Sea (18.5°C to 33.9°C). In the endemic Red Sea lineage RS_LinB, the mtORF and the heat shock protein gene hsp70 uncovered similar phylogeographic patterns strongly correlated with environmental variations. This suggests that the mtORF too might be involved in thermal adaptation. Here, we used computational analyses to explore the features and putative function of this mtORF. In particular, we tested the likelihood that this gene encodes a functional protein and whether it may play a role in adaptation. Analyses of full mitogenomes showed that the mtORF originated in the common ancestor of Madracis and other pocilloporids, and that it encodes a transmembrane protein differing in length and domain architecture among genera. Homology-based annotation and the relative conservation of metal-binding sites revealed traces of an ancient hydrolase catalytic activity. Furthermore, signals of pervasive purifying selection, lack of stop codons in 1830 sequences analyzed, and a codon-usage bias similar to that of other mitochondrial genes indicate that the protein is functional, i.e., not a pseudogene. Other features, such as intrinsically disordered regions, tandem repeats, and signals of positive selection particularly in StylophoraRS_LinB populations, are consistent with a role of the mtORF in adaptive responses to environmental changes.

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.

The genome of the jellyfish Clytia hemisphaerica and the evolution of the cnidarian life-cycle

Jellyfish (medusae) are a distinctive life-cycle stage of medusozoan cnidarians. They are major marine predators, with integrated neurosensory, muscular and organ systems. The genetic foundations of this complex form are largely unknown. We report the draft genome of the hydrozoan jellyfish Clytia hemisphaerica and use multiple transcriptomes to determine gene use across life-cycle stages. Medusa, planula larva and polyp are each characterized by distinct transcriptome signatures reflecting abrupt life-cycle transitions and all deploy a mixture of phylogenetically old and new genes. Medusa-specific transcription factors, including many with bilaterian orthologues, associate with diverse neurosensory structures. Compared to Clytia, the polyp-only hydrozoan Hydra has lost many of the medusa-expressed transcription factors, despite similar overall rates of gene content evolution and sequence evolution. Absence of expression and gene loss among Clytia orthologues of genes patterning the anthozoan aboral pole, secondary axis and endomesoderm support simplification of planulae and polyps in Hydrozoa, including loss of bilateral symmetry. Consequently, although the polyp and planula are generally considered the ancestral cnidarian forms, in Clytia the medusa maximally deploys the ancestral cnidarian-bilaterian transcription factor gene complement.

The Mitochondrial Endonuclease M20 Participates in the Down-Regulation of Mitochondrial DNA in Pollen Cells

Maintaining the appropriate number of mitochondrial DNA (mtDNA) molecules is crucial for supporting mitochondrial metabolism and function in both plant and animal cells. For example, a substantial decrease in mtDNA levels occurs as a key part of pollen development. The molecular mechanisms regulating mtDNA copy number are largely unclear, particularly with regard to those that reduce mtDNA levels. Here, we identified and purified a 20-kD endonuclease, M20, from maize (Zea mays) pollen mitochondria. We found M20 to be an His-Asn-His/Asn (H-N-H/N) nuclease that degrades linear and circular DNA in the presence of Mg²⁺ or Mn²⁺ Arabidopsis (Arabidopsis thaliana) AtM20, which shared high sequence similarity with maize M20, localized to the mitochondria, had a similar H-N-H/N structure, and degraded both linear and circular DNA. AtM20 transcript levels increased during pollen development, in parallel with a rapid reduction in mtDNA. Clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 genome-editing techniques were used to generate knockout lines of AtM20 (atm20), which exhibited a significant delay in the reduction in mtDNA levels in pollen vegetative cells but normal mtDNA levels in somatic cells. The delayed reduction in pollen mtDNA levels was rescued by the transgenic expression of AtM20 in atm20 plants. This study thus uncovers an endonucleolytic DNase in plant mitochondria and its crucial role in reducing mtDNA levels, pointing to the complex mechanism regulating mtDNA levels in plants.

Unraveling chloroplast transcriptomes with ChloroSeq, an organelle RNA-Seq bioinformatics pipeline

Online sequence repositories are teeming with RNA sequencing (RNA-Seq) data from a wide range of eukaryotes. Although most of these data sets contain large numbers of organelle-derived reads, researchers tend to ignore these data, focusing instead on the nuclear-derived transcripts. Consequently, GenBank contains massive amounts of organelle RNA-Seq data that are just waiting to be downloaded and analyzed. Recently, a team of scientists designed an open-source bioinformatics program called ChloroSeq, which systemically analyzes an organelle transcriptome using RNA-Seq. The ChloroSeq pipeline uses RNA-Seq alignment data to deliver detailed analyses of organelle transcriptomes, which can be fed into statistical software for further analysis and for generating graphical representations of the data. In addition to providing data on expression levels via coverage statistics, ChloroSeq can examine splicing efficiency and RNA editing profiles. Ultimately, ChloroSeq provides a well-needed avenue for researchers of all stripes to start exploring organelle transcription and could be a key step toward a more thorough understanding of organelle gene expression.

Recent Asian origin of chytrid fungi causing global amphibian declines

Globalized infectious diseases are causing species declines worldwide, but their source often remains elusive. We used whole-genome sequencing to solve the spatiotemporal origins of the most devastating panzootic to date, caused by the fungus Batrachochytrium dendrobatidis, a proximate driver of global amphibian declines. We traced the source of B. dendrobatidis to the Korean peninsula, where one lineage, BdASIA-1, exhibits the genetic hallmarks of an ancestral population that seeded the panzootic. We date the emergence of this pathogen to the early 20th century, coinciding with the global expansion of commercial trade in amphibians, and we show that intercontinental transmission is ongoing. Our findings point to East Asia as a geographic hotspot for B. dendrobatidis biodiversity and the original source of these lineages that now parasitize amphibians worldwide.