Eight new mtDNA sequences of glass sponges reveal an extensive usage of +1 frameshifting in mitochondrial translation

Phylomitogenomics bolsters the high-level classification of Demospongiae (phylum Porifera)

Class Demospongiae is the largest in the phylum Porifera (Sponges) and encompasses nearly 8,000 accepted species in three subclasses: Keratosa, Verongimorpha, and Heteroscleromorpha. Subclass Heteroscleromorpha contains ∼90% of demosponge species and is subdivided into 17 orders. The higher level classification of demosponges underwent major revision as the result of nearly three decades of molecular studies. However, because most of the previous molecular work only utilized partial data from a small number of nuclear and mitochondrial (mt) genes, this classification scheme needs to be tested by larger datasets. Here we compiled a mt dataset for 136 demosponge species-including 64 complete or nearly complete and six partial mt-genome sequences determined or assembled for this study-and used it to test phylogenetic relationships among Demospongiae in general and Heteroscleromorpha in particular. We also investigated the phylogenetic position of Myceliospongia araneosa, a highly unusual demosponge without spicules and spongin fibers, currently classified as Demospongiae incertae sedis, for which molecular data were not available. Our results support the previously inferred sister-group relationship between Heteroscleromorpha and Keratosa + Verongimorpha and suggest five main clades within Heteroscleromorpha: Clade C0 composed of order Haplosclerida; Clade C1 composed of Scopalinida, Sphaerocladina, and Spongillida; Clade C2 composed of Axinellida, Biemnida, Bubarida; Clade C3 composed of Tetractinellida; and Clade C4 composed of Agelasida, Clionaida, Desmacellida, Merliida, Suberitida, Poecilosclerida, Polymastiida, and Tethyida. The inferred relationships among these clades were (C0(C1(C2(C3+C4)))). Analysis of molecular data from M. araneosa placed it in the C3 clade as a sister taxon to the highly skeletonized tetractinellids Microscleroderma sp. and Leiodermatium sp. Molecular clock analysis dated divergences among the major clades in Heteroscleromorpha from the Cambrian to the Early Silurian, the origins of most heteroscleromorph orders in the middle Paleozoic, and the most basal splits within these orders around the Paleozoic to Mesozoic transition. Overall, the results of this study are mostly congruent with the accepted classification of Heteroscleromorpha, but add temporal perspective and new resolution to phylogenetic relationships within this subclass.

Genome skimming elucidates the evolutionary history of Octopoda

Phylogenies for Octopoda have, until now, been based on morphological characters or a few genes. Here we provide the complete mitogenomes and the nuclear 18S and 28S ribosomal genes of twenty Octopoda specimens, comprising 18 species of Cirrata and Incirrata, representing 13 genera and all five putative families of Cirrata (Cirroctopodidae, Cirroteuthidae, Grimpoteuthidae, Opisthoteuthidae and Stauroteuthidae) and six families of Incirrata (Amphitretidae, Argonautidae, Bathypolypodidae, Eledonidae, Enteroctopodidae, and Megaleledonidae) which were assembled using genome skimming. Phylogenetic trees were built using Maximum Likelihood and Bayesian Inference with several alignment matrices. All mitochondrial genomes had the 'typical' genome composition and gene order previously reported for octopodiforms, except Bathypolypus ergasticus, which appears to lack ND5, two tRNA genes that flank ND5 and two other tRNA genes. Argonautoidea was revealed as sister to Octopodidae by the mitochondrial protein-coding gene dataset, however, it was recovered as sister to all other incirrate octopods with strong support in an analysis using nuclear rRNA genes. Within Cirrata, our study supports two existing classifications suggesting neither is likely in conflict with the true evolutionary history of the suborder. Genome skimming is useful in the analysis of phylogenetic relationships within Octopoda; inclusion of both mitochondrial and nuclear data may be key.

Low genetic diversity and predation threaten a rediscovered marine sponge

Discovered in 1819 in the tropical waters off Singapore, the magnificent Neptune's cup sponge Cliona patera (Hardwicke, 1820) was harvested for museums and collectors until it was presumed extinct worldwide for over a century since 1907. Recently in 2011, seven living individuals were rediscovered in Singapore with six relocated to a marine protected area in an effort to better monitor and protect the population, as well as to enhance external fertilisation success. To determine genetic diversity within the population, we sequenced the complete mitochondrial genomes and nuclear ribosomal DNA of these six individuals and found extremely limited variability in their genes. The low genetic diversity of this rediscovered population is confirmed by comparisons with close relatives of C. patera and could compromise the population's ability to recover from environmental and anthropogenic pressures associated with the highly urbanised coastlines of Singapore. This lack of resilience is compounded by severe predation which has been shrinking sponge sizes by up to 5.6% every month. Recovery of this highly endangered population may require ex situ approaches and crossbreeding with other populations, which are also rare.

Mitochondrial Genomes in Perkinsus Decode Conserved Frameshifts in All Genes

Mitochondrial genomes of apicomplexans, dinoflagellates, and chrompodellids that collectively make up the Myzozoa, encode only three proteins (Cytochrome b [COB], Cytochrome c oxidase subunit 1 [COX1], Cytochrome c oxidase subunit 3 [COX3]), contain fragmented ribosomal RNAs, and display extensive recombination, RNA trans-splicing, and RNA-editing. The early-diverging Perkinsozoa is the final major myzozoan lineage whose mitochondrial genomes remained poorly characterized. Previous reports of Perkinsus genes indicated independent acquisition of non-canonical features, namely the occurrence of multiple frameshifts. To determine both ancestral myzozoan and novel perkinsozoan mitochondrial genome features, we sequenced and assembled mitochondrial genomes of four Perkinsus species. These data show a simple ancestral genome with the common reduced coding capacity but disposition for rearrangement. We identified 75 frameshifts across the four species that occur as distinct types and that are highly conserved in gene location. A decoding mechanism apparently employs unused codons at the frameshift sites that advance translation either +1 or +2 frames to the next used codon. The locations of frameshifts are seemingly positioned to regulate protein folding of the nascent protein as it emerges from the ribosome. The cox3 gene is distinct in containing only one frameshift and showing strong selection against residues that are otherwise frequently encoded at the frameshift positions in cox1 and cob. All genes lack cysteine codons implying a reduction to 19 amino acids in these genomes. Furthermore, mitochondrion-encoded rRNA fragment complements are incomplete in Perkinsus spp. but some are found in the nuclear DNA suggesting import into the organelle. Perkinsus demonstrates further remarkable trajectories of organelle genome evolution including pervasive integration of frameshift translation into genome expression.

The mitogenome of common snipe, Gallinago gallinago gallinago Linnaeus, 1758 and evolutionary implications for the family Scolopacidae

The mitochondrial genome (mitogenome) of Gallinago gallinago gallinago Linnaeus, 1758 was determined by the high-throughput data. The assembled mitogenome was 16,919 bp in length, with a 58.7% A + T content and GC skew of −0.3850. Among 13 PCGs, an unusual start codon (GTG) was identified for the COX1 gene, and incomplete stop codons (T-) were found in the COX3, ND2 and ND4 genes. The function of a cytosine insertion at site 174 in the ND3 gene and its phylogenetic significance are worthy of further scrutiny. In the control region (CR), thirteen 15-bp simple sequence repeats were found in G. g. gallinago. Phylogenetic analysis indicated that Gallinago was clustered at the basal position of the Scolopax clade and that the monophyly of Gallinago was also recovered. The mitogenome data of G. g. gallinago provides useful resources for further studying the evolution of Scolopacidae.

Multiple origins of a frameshift insertion in a mitochondrial gene in birds and turtles

BACKGROUND

During evolutionary history, molecular mechanisms have emerged to cope with deleterious mutations. Frameshift insertions in protein-coding sequences are extremely rare because they disrupt the reading frame. There are a few known examples of their correction through translational frameshifting, a process that enables ribosomes to skip nucleotides during translation to regain proper reading frame. Corrective frameshifting has been proposed to act on the single base pair insertion at position 174 of the mitochondrial NADH dehydrogenase subunit 3 gene (ND3) that has been observed in several turtles and birds. However, the relatively sparse taxonomic representation has hampered our understanding of the evolution of this insertion in vertebrates.

RESULTS

Here, we analyzed 87,707 ND3 sequences from 10,309 vertebrate taxa to reveal the evolutionary history of this insertion and its common genomic characteristics. We confirmed that the insertion only appears in turtles and birds and reconstructed that it evolved independently in both groups with complex patterns of gains and losses. The insertion was observed in almost all bird orders but was absent in all members of the diverse Passeriformes. We found strong conservation in the nucleotides surrounding the insertion in both turtles and birds, which implies that the insertion enforces structural constraints that could be involved in its correction.

CONCLUSIONS

Our study demonstrates that frameshifts can be widespread and can be retained for millions of years if they are embedded in a conserved sequence theme.

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.

Partial Mitochondrial Genome Sequences of Two Abyssal Sponges (Porifera: Hexactinellida), Bathydorus laniger and Docosaccus maculatus

We announce the nearly complete mitochondrial genome sequences of two hexactinellid sponges, Bathydorus laniger and Docosaccus maculatus. A contiguous region of over 15,000 bp was sequenced from each genome. An uncommon structural element was identified as a series of repetitive elements with sequences matching cob in the genome of D. maculatus.

An integrative systematic framework helps to reconstruct skeletal evolution of glass sponges (Porifera, Hexactinellida)

Background

Glass sponges (Class Hexactinellida) are important components of deep-sea ecosystems and are of interest from geological and materials science perspectives. The reconstruction of their phylogeny with molecular data has only recently begun and shows a better agreement with morphology-based systematics than is typical for other sponge groups, likely because of a greater number of informative morphological characters. However, inconsistencies remain that have far-reaching implications for hypotheses about the evolution of their major skeletal construction types (body plans). Furthermore, less than half of all described extant genera have been sampled for molecular systematics, and several taxa important for understanding skeletal evolution are still missing. Increased taxon sampling for molecular phylogenetics of this group is therefore urgently needed. However, due to their remote habitat and often poorly preserved museum material, sequencing all 126 currently recognized extant genera will be difficult to achieve. Utilizing morphological data to incorporate unsequenced taxa into an integrative systematics framework therefore holds great promise, but it is unclear which methodological approach best suits this task.

Results

Here, we increase the taxon sampling of four previously established molecular markers (18S, 28S, and 16S ribosomal DNA, as well as cytochrome oxidase subunit I) by 12 genera, for the first time including representatives of the order Aulocalycoida and the type genus of Dactylocalycidae, taxa that are key to understanding hexactinellid body plan evolution. Phylogenetic analyses suggest that Aulocalycoida is diphyletic and provide further support for the paraphyly of order Hexactinosida; hence these orders are abolished from the Linnean classification. We further assembled morphological character matrices to integrate so far unsequenced genera into phylogenetic analyses in maximum parsimony (MP), maximum likelihood (ML), Bayesian, and morphology-based binning frameworks. We find that of these four approaches, total-evidence analysis using MP gave the most plausible results concerning congruence with existing phylogenetic and taxonomic hypotheses, whereas the other methods, especially ML and binning, performed more poorly. We use our total-evidence phylogeny of all extant glass sponge genera for ancestral state reconstruction of morphological characters in MP and ML frameworks, gaining new insights into the evolution of major hexactinellid body plans and other characters such as different spicule types.

Conclusions

Our study demonstrates how a comprehensive, albeit in some parts provisional, phylogeny of a larger taxon can be achieved with an integrative approach utilizing molecular and morphological data, and how this can be used as a basis for understanding phenotypic evolution. The datasets and associated trees presented here are intended as a resource and starting point for future work on glass sponge evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12983-017-0191-3) contains supplementary material, which is available to authorized users.

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.

Ribosomal frameshifting and transcriptional slippage: From genetic steganography and cryptography to adventitious use

Genetic decoding is not ‘frozen’ as was earlier thought, but dynamic. One facet of this is frameshifting that often results in synthesis of a C-terminal region encoded by a new frame. Ribosomal frameshifting is utilized for the synthesis of additional products, for regulatory purposes and for translational ‘correction’ of problem or ‘savior’ indels. Utilization for synthesis of additional products occurs prominently in the decoding of mobile chromosomal element and viral genomes. One class of regulatory frameshifting of stable chromosomal genes governs cellular polyamine levels from yeasts to humans. In many cases of productively utilized frameshifting, the proportion of ribosomes that frameshift at a shift-prone site is enhanced by specific nascent peptide or mRNA context features. Such mRNA signals, which can be 5′ or 3′ of the shift site or both, can act by pairing with ribosomal RNA or as stem loops or pseudoknots even with one component being 4 kb 3′ from the shift site. Transcriptional realignment at slippage-prone sequences also generates productively utilized products encoded trans-frame with respect to the genomic sequence. This too can be enhanced by nucleic acid structure. Together with dynamic codon redefinition, frameshifting is one of the forms of recoding that enriches gene expression.

Extensive Mitochondrial mRNA Editing and Unusual Mitochondrial Genome Organization in Calcaronean Sponges

One of the unusual features of DNA-containing organelles in general and mitochondria in particular is the frequent occurrence of RNA editing [1]. The term "RNA editing" refers to a variety of mechanistically unrelated biochemical processes that alter RNA sequence during or after transcription [2]. The editing can be insertional, deletional, or substitutional and has been found in all major types of RNAs [3, 4]. Although mitochondrial mRNA editing is widespread in some eukaryotic lineages [5-7], it is rare in animals, with reported cases limited both in their scope and in phylogenetic distribution [8-11] (see also [12]). While analyzing genomic data from calcaronean sponges Sycon ciliatum and Leucosolenia complicata, we were perplexed by the lack of recognizable mitochondrial coding sequences. Comparison of genomic and transcriptomic data from these species revealed the presence of mitochondrial cryptogenes whose transcripts undergo extensive editing. This editing consisted of single or double uridylate (U) insertions in pre-existing short poly(U) tracts. Subsequent analysis revealed the presence of similar editing in Sycon coactum and the loss of editing in Petrobiona massiliana, a hypercalcified calcaronean sponge. In addition, mitochondrial genomes of at least some calcaronean sponges were found to have a highly unusual architecture, with nearly all genes located on individual and likely linear chromosomes. Phylogenetic analysis of mitochondrial coding sequences revealed accelerated rates of sequence evolution in this group. The latter observation presents a challenge for the mutational-hazard hypothesis [13], which posits that mRNA editing should not occur in lineages with an elevated mutation rate.

The mitochondrial genome of the deep-sea glass sponge Lophophysema eversa (Porifera, Hexacinellida, Hyalonematidae)

We reported a nearly complete mitochondrial genome (mitogenome) from the glass sponge Lophophysema eversa, the second mitogenome in the order Amphidiscosida and the ninth in the class Hexactinellida. It is 20,651 base pairs in length and contains 39 genes including 13 protein-coding genes, 2 ribosomal RNA subunit genes and 24 tRNA genes. The gene content and order of L. eversa are identical to those of Tabachnickia sp., the other species with a sequenced mitogenome in Amphidiscosida, except with two additional tRNAs and three tRNA translocations. The cob gene has a +1 translational frameshift. These results will contribute to a better understanding of the phylogeny of glass sponges.

Mitochondrial Genome Sequence of the Glass Sponge Oopsacas minuta

We report the complete mitochondrial genome sequence of the Mediterranean glass sponge Oopsacas minuta. This 19-kb mitochondrial genome has 24 noncoding genes (22 tRNAs and 2 rRNAs) and 14 protein-encoding genes coding for 11 subunits of respiratory chain complexes and 3 ATP synthase subunits.

Cytonuclear Interactions in the Evolution of Animal Mitochondrial tRNA Metabolism

The evolution of mitochondrial information processing pathways, including replication, transcription and translation, is characterized by the gradual replacement of mitochondrial-encoded proteins with nuclear-encoded counterparts of diverse evolutionary origins. Although the ancestral enzymes involved in mitochondrial transcription and replication have been replaced early in eukaryotic evolution, mitochondrial translation is still carried out by an apparatus largely inherited from the α-proteobacterial ancestor. However, variation in the complement of mitochondrial-encoded molecules involved in translation, including transfer RNAs (tRNAs), provides evidence for the ongoing evolution of mitochondrial protein synthesis. Here, we investigate the evolution of the mitochondrial translational machinery using recent genomic and transcriptomic data from animals that have experienced the loss of mt-tRNAs, including phyla Cnidaria and Ctenophora, as well as some representatives of all four classes of Porifera. We focus on four sets of mitochondrial enzymes that directly interact with tRNAs: Aminoacyl-tRNA synthetases, glutamyl-tRNA amidotransferase, tRNA^Ile lysidine synthetase, and RNase P. Our results support the observation that the fate of nuclear-encoded mitochondrial proteins is influenced by the evolution of molecules encoded in mitochondrial DNA, but in a more complex manner than appreciated previously. The data also suggest that relaxed selection on mitochondrial translation rather than coevolution between mitochondrial and nuclear subunits is responsible for elevated rates of evolution in mitochondrial translational proteins.

The mitochondrial genome of Ixobrychus flavicollis (Ardeidae: Ixobrychus)

Ixobrychus flavicollis is the largest bittern in Ixobrychus which has very wide geographical range. In our results, the complete mitochondrial genome of I. flavicollis is 17,016 bp in length. Thirteen protein-coding genes, two rRNA genes, twenty-two tRNA genes and one non-coding region (D-loop) are present in this mitochondrial genome. No special arrangement when compared with other typical mitochondrial of vertebrates. The percentage of A and T (56.7%) was slightly higher than G and C.