Recoding of translation in turtle mitochondrial genomes: programmed frameshift mutations and evidence of a modified genetic code

Complete mitochondrial genome of the southern painted turtle (Chrysemys dorsalis, Testudines: Emydidae) in Korea

The complete mitochondrial genome of Chrysemys dorsalis in Korea was sequenced and characterized. The mitochondrial genome is 17,258 bp in length and the GC content is 39%. It is constituted of 37 genes, 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a noncoding region. Phylogenetic analysis reveals that C. dorsalis forms a monophyletic group with C. picta turtles but is distinctly separated from them, aligning with previous findings. In Korea, C. dorsalis forms a discrete clade, separate from both native and invasive turtle species. No evidence of genetic disturbance or intermingling is observed. This is the first case of a complete mitochondrial genome from C. dorsalis and provides crucial data for understanding C. dorsalis and managing invasive species effectively, emphasizing the need for continued mitochondrial genome data accumulation.

Complete mitochondrial genome of the razor-backed musk turtle (Sternotherus carinatus, testudines: emydidae) in Korea

Sternotherus carinatus has been considered as a potential invasive species in Korea. However, the mitochondrial genome information of S. carinatus which can be used to control its effect on ecosystem is lacking. In this study, the complete mitochondrial genome of S. carinatus in Korea was sequenced and characterized. The mitochondrial genome consists of 37 genes (13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes) and a noncoding region. Phylogenetic analysis based on the mitochondrial genome sequences showed that S. carinatus from Korea is separated from other turtles which are the invasive species in Korea. Sequence divergence calculations indicated near-zero divergence between S. carinatus populations in Korea, the USA, and China, suggesting limited genetic differentiation. In the context of the broader issue of invasive species disrupting ecosystems, this research contributes to the identification of mitochondrial genomes for various freshwater turtle species, emphasizing the need for extended data collection to discern genetic mixing trends between native and non-native species. This study is a significant step toward managing S. carinatus as a potential invasive species in Korea.

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.

More discussion of minimalist species descriptions and clarifying some misconceptions contained in Meier et al. 2021

This is a response to a preprint version of "A re-analysis of the data in Sharkey et al.'s (2021) minimalist revision reveals that BINs do not deserve names, but BOLD Systems needs a stronger commitment to open science", https://www.biorxiv.org/content/10.1101/2021.04.28.441626v2. Meier et al. strongly criticized Sharkey et al.'s publication in which 403 new species were deliberately minimally described, based primarily on COI barcode sequence data. Here we respond to these criticisms. The following points are made: 1) Sharkey et al. did not equate BINs with species, as demonstrated in several examples in which multiple species were found to be in single BINs. 2) We reiterate that BINs were used as a preliminary sorting tool, just as preliminary morphological identification commonly sorts specimens based on color and size into unit trays; despite BINs and species concepts matching well over 90% of species, this matching does not equate to equality. 3) Consensus barcodes were used only to provide a diagnosis to conform to the rules of the International Code of Zoological Nomenclature just as consensus morphological diagnoses are. The barcode of a holotype is definitive and simply part of its cellular morphology. 4) Minimalist revisions will facilitate and accelerate future taxonomic research, not hinder it. 5) We refute the claim that the BOLD sequences of Plesiocoelusvanachterbergi are pseudogenes and demonstrate that they simply represent a frameshift mutation. 6) We reassert our observation that morphological evidence alone is insufficient to recognize species within species-rich higher taxa and that its usefulness lies in character states that are congruent with molecular data. 7) We show that in the cases in which COI barcodes code for the same amino acids in different putative species, data from morphology, host specificity, and other ecological traits reaffirm their utility as indicators of genetically distinct lineages.

Frameshift and wild-type proteins are often highly similar because the genetic code and genomes were optimized for frameshift tolerance

Frameshift mutations have been considered of significant importance for the molecular evolution of proteins and their coding genes, while frameshift protein sequences encoded in the alternative reading frames of coding genes have been considered to be meaningless. However, functional frameshifts have been found widely existing. It was puzzling how a frameshift protein kept its structure and functionality while substantial changes occurred in its primary amino-acid sequence. This study shows that the similarities among frameshifts and wild types are higher than random similarities and are determined at different levels. Frameshift substitutions are more conservative than random substitutions in the standard genetic code (SGC). The frameshift substitutions score of SGC ranks in the top 2.0-3.5% of alternative genetic codes, showing that SGC is nearly optimal for frameshift tolerance. In many genes and certain genomes, frameshift-resistant codons and codon pairs appear more frequently than expected, suggesting that frameshift tolerance is achieved through not only the optimality of the genetic code but, more importantly, the further optimization of a specific gene or genome through the usages of codons/codon pairs, which sheds light on the role of frameshift mutations in molecular and genomic evolution.

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.

The complete mitochondrial genome of Indian Cuckoo Cuculus micropterus (Aves: Cuculiformes)

The Indian Cuckoo, Cuculus micropterus, belongs to the family Cuculidae. In this paper, we sequenced and analysized the complete mitochondrial genome of C. micropterus. The complete mitochondrial genome of C. micropterus is 17,541 bp in length, which was submitted to the NCBI database under the accession number MZ048030. It contains 13 protein-coding genes, 22 transfer RNA genes, two ribosome RNA genes, and one non-coding control region. The overall base composition of the mitochondrial DNA is 33.2% for A, 24.2% for T, 29.8% for C, and 12.8% for G, with a GC content of 42.6%. In order to explore the molecular phylogenetics evolution of Cuculidae, the nucleotide sequence data of 13 PCGs of C. micropterus and other 11 Cuculiformes were used for the phylogenetic analysis. The result shows that C. micropterus is closely related to Cuculus canorus bakeri. The study contributes to illuminating the taxonomic status of C. micropterus, and may facilitate further investigation of the evolution of Cuculidae.

Complete mitochondrial genome of the Cumberland slider (Trachemys scripta troostii, Testudienes: Emydidae) in Korea

The complete mitochondrial (mt) genome of Trachemys scripta troostii was sequenced and was characterized, which comprised 37 genes (13 protein-coding genes, 22 transfer RNAs, and 2 ribosomal RNAs) and a non-coding control region. Phylogenetic analysis based on the full mt genome indicated that T. s. troostii was more closely related to T. scripta from Canada than to T. s. elegans from China or T. s. scripta fom China. This is the first complete mt genome from T. s. troostii, which provides data for further study of phylogeny in Emydidae.

Characterization of the complete mitogenomes of Baillon's Crake Porzana pusilla and phylogenetic analysis

The complete mitochondrial genome of the Baillon's Crake Porzana pusilla (Gruiformes: Rallidae) are sequenced and annotated, which contained 37 typical genes. The length of the complete mitochondrial genome is 16,966 bp (GenBank No. MW043485), with As, Ts, Cs, Gs, and AT content of the mitochondrial genome is 32.1, 23.2, 30.9, 13.8, and 55.3%, respectively. All protein-coding genes started with ATN except COX1 and ND5, which start with GTG, and all protein-coding genes end with a complete triplet codon (TAA, AGG, AGA, and TAG), except COX3, which ends with an incomplete T. The ND3 gene of P. pusilla with an extra C nucleotide in 174 site. Phylogenetic analysis revealed that the new sequenced species of P. pusilla was closer to the clade of Porzana fusca and Porzana paykullii, and all three Porzana are clustered into one branch.

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.

Prediction of mitochondrial genome-wide variation through sequencing of mitochondrion-enriched extracts

Although mitochondrial DNA (mtDNA) haplotype variation is often applied for estimating population dynamics and phylogenetic relationships, economical and generalized methods for entire mtDNA genome enrichment prior to high-throughput sequencing are not readily available. This study demonstrates the utility of differential centrifugation to enrich for mitochondrion within cell extracts prior to DNA extraction, short-read sequencing, and assembly using exemplars from eight maternal lineages of the insect species, Ostrinia nubilalis. Compared to controls, enriched extracts showed a significant mean increase of 48.2- and 86.1-fold in mtDNA based on quantitative PCR, and proportion of subsequent short sequence reads that aligned to the O. nubilalis reference mitochondrial genome, respectively. Compared to the reference genome, our de novo assembled O. nubilalis mitochondrial genomes contained 82 intraspecific substitution and insertion/deletion mutations, and provided evidence for correction of mis-annotated 28 C-terminal residues within the NADH dehydrogenase subunit 4. Comparison to a more recent O. nubilalis mtDNA assembly from unenriched short-read data analogously showed 77 variant sites. Twenty-eight variant positions, and a triplet ATT codon (Ile) insertion within ATP synthase subunit 8, were unique within our assemblies. This study provides a generalizable pipeline for whole mitochondrial genome sequence acquisition adaptable to applications across a range of taxa.

The role of selection in the evolution of marine turtles mitogenomes

Sea turtles are the only extant chelonian representatives that inhabit the marine environment. One key to successful colonization of this habitat is the adaptation to different energetic demands. Such energetic requirement is intrinsically related to the mitochondrial ability to generate energy through oxidative phosphorylation (OXPHOS) process. Here, we estimated Testudines phylogenetic relationships from 90 complete chelonian mitochondrial genomes and tested the adaptive evolution of 13 mitochondrial protein-coding genes of sea turtles to determine how natural selection shaped mitochondrial genes of the Chelonioidea clade. Complete mitogenomes showed strong support and resolution, differing at the position of the Chelonioidea clade in comparison to the turtle phylogeny based on nuclear genomic data. Codon models retrieved a relatively increased dN/dS (ω) on three OXPHOS genes for sea turtle lineages. Also, we found evidence of positive selection on at least three codon positions, encoded by NADH dehydrogenase genes (ND4 and ND5). The accelerated evolutionary rates found for sea turtles on COX2, ND1 and CYTB and the molecular footprints of positive selection found on ND4 and ND5 genes may be related to mitochondrial molecular adaptation to stress likely resulted from a more active lifestyle in sea turtles. Our study provides insight into the adaptive evolution of the mtDNA genome in sea turtles and its implications for the molecular mechanism of oxidative phosphorylation.

Structural Variation of the Turtle Mitochondrial Control Region

The present study describes the most comprehensive comparison of turtle mtD-loop regions to date. The primary structure was compared from DNA sequences accessed from GenBank from 48 species in 13 families of extant turtles, and secondary structures of the mtD-loop region were inferred from thermal stabilities, using the program Mfold, for each superfamiliy of turtles. Both primary and secondary structures were found to be highly variable across the order. The Cryptodira showed conservation in the primary structure at conserved sequence blocks (CSBs), but the Pleurodira displayed limited conservation of primary structural characters, other than the coreTAS, a binding site for the helicase TWINKLE, which was highly conserved in the Central and Right Domains across the order. No secondary structure was associated with a TAS, but an AT-rich fold (secondary structure) near the 3' terminus of the mtD-loop region was detected in all turtle superfamilies. Mapping of character states of structural features of the mtD-loop region revealed that most character states were autapomorphies and inferred a number of homoplasies. The Left Domain of turtles, containing no highly conserved structural elements, likely does not serve a functional role; therefore, the Central Domain in turtles is likely equivalent to the Left Domain of mammals. The AT-rich secondary structural element near the 3' terminus of the mtD-loop region may be conserved across turtles because of a functional role, perhaps containing the Light Strand Promotor, or perhaps interacting with the TWINKLE-coreTAS complex in the Central and Right Domains to regulate mtDNA replication and transcription.

Natural pyrrolysine-biased translation of stop codons in mitochondrial peptides entirely coded by expanded codons

During the noncanonical deletion transcription, k nucleotides are systematically skipped/deleted after each transcribed trinucleotide producing deletion-RNAs (delRNAs). Peptides matching delRNAs either result from (a) canonical translation of delRNAs; or (b) noncanonical translation of regular transcripts along expanded codons. Only along frame "0" (start site) (a) and (b) produce identical peptides. Here, mitochondrial mass spectrometry data analyses assume expanded codon/del-transcription with 3 + k (k from 0 to 12) nucleotides. Detected peptides map preferentially on previously identified delRNAs. More peptides were detected for k (1-12) when del-transcriptional and expanded codon translations start sites coincide (i.e. the 0^th frame) than for frames +1 or +2. Hence, both (a) and (b) produced peptides identified here. Biases for frame 0 decrease for k > 2, reflecting codon/anticodon expansion limits. Further analyses find preferential pyrrolysine insertion at stop codons, suggesting Pyl-specific mitochondrial suppressor tRNAs loaded by Pyl-specific tRNA synthetases with unknown origins. Pyl biases at stops are stronger for regular than expanded codons suggesting that Pyl-tRNAs are less competitive with near-cognate tRNAs in expanded codon contexts. Statistical biases for these findings exclude that detected peptides are experimental and/or bioinformatic artefacts implying both del-transcription and expanded codons translation occur in human mitochondria.

Unorthodox features in two venerid bivalves with doubly uniparental inheritance of mitochondria

In animals, strictly maternal inheritance (SMI) of mitochondria is the rule, but one exception (doubly uniparental inheritance or DUI), marked by the transmission of sex-specific mitogenomes, has been reported in bivalves. Associated with DUI is a frequent modification of the mitochondrial cox2 gene, as well as additional sex-specific mitochondrial genes not involved in oxidative phosphorylation. With the exception of freshwater mussels (for 3 families of the order Unionida), these DUI-associated features have only been shown in few species [within Mytilidae (order Mytilida) and Veneridae (order Venerida)] because of the few complete sex-specific mitogenomes published for these orders. Here, we present the complete sex-specific mtDNAs of two recently-discovered DUI species in two families of the order Venerida, Scrobicularia plana (Semelidae) and Limecola balthica (Tellinidae). These species display the largest differences in genome size between sex-specific mitotypes in DUI species (>10 kb), as well as the highest mtDNA divergences (sometimes reaching >50%). An important in-frame insertion (>3.5 kb) in the male cox2 gene is partly responsible for the differences in genome size. The S. plana cox2 gene is the largest reported so far in the Kingdom Animalia. The mitogenomes may be carrying sex-specific genes, indicating that general mitochondrial features are shared among DUI species.

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.

Chimeric Translation for Mitochondrial Peptides: Regular and Expanded Codons

Frameshifting protein translation occasionally results from insertion of amino acids at isolated mono- or dinucleotide-expanded codons by tRNAs with expanded anticodons. Previous analyses of two different types of human mitochondrial MS proteomic data (Fisher and Waters technologies) detect peptides entirely corresponding to expanded codon translation. Here, these proteomic data are reanalyzed searching for peptides consisting of at least eight consecutive amino acids translated according to regular tricodons, and at least eight adjacent consecutive amino acids translated according to expanded codons. Both datasets include chimerically translated peptides (mono- and dinucleotide expansions, 42 and 37, respectively). The regular tricodon-encoded part of some chimeric peptides corresponds to standard human mitochondrial proteins (mono- and dinucleotide expansions, six (AT6, CytB, ND1, 2xND2, ND5) and one (ND1), respectively). Chimeric translation probably increases the diversity of mitogenome-encoded proteins, putatively producing functional proteins. These might result from translation by tRNAs with expanded anticodons, or from regular tricodon translation of RNAs where transcription/posttranscriptional edition systematically deleted mono- or dinucleotides after each trinucleotide. The pairwise matched combination of adjacent peptide parts translated from regular and expanded codons strengthens the hypothesis that translation of stretches of consecutive expanded codons occurs. Results indicate statistical translation producing distributions of alternative proteins. Genetic engineering should account for potential unexpected, unwanted secondary products.

The Complete Mitochondrial Genome of Platysternon megacephalum peguense and Molecular Phylogenetic Analysis

Platysternon megacephalum is the only living representative species of Platysternidae and only three subspecies remain: P. m. megalorcephalum, P. m. shiui, and P. m. peguense. However, previous reports implied that P. m. peguense has distinct morphological and molecular features. The characterization of the mitogenome has been accepted as an efficient means of phylogenetic and evolutionary analysis. Hence, this study first determined the complete mitogenome of P. m. peguense with the aim to identify the structure and variability of the P. m. peguense mitogenome through comparative analysis. Furthermore, the phylogenetic relationship of the three subspecies was tested. Based on different tRNA gene loss and degeneration of these three subspecies, their rearrangement pathways have been inferred. Phylogenetic analysis showed that P. m. peguense is a sister group to (P. m. megalorcephalum and P. m. shiui). Furthermore, the divergence time estimation of these three subspecies coincided with the uplift of the Tibetan Plateau. This study shows that the genetic distances between P. m. peguense and the other two subspecies are comparable to interspecific genetic distances, for example within Mauremys. In general, this study provides new and meaningful insights into the evolution of the three Platysternidae subspecies.

MACSE v2: Toolkit for the Alignment of Coding Sequences Accounting for Frameshifts and Stop Codons

Multiple sequence alignment is a prerequisite for many evolutionary analyses. Multiple Alignment of Coding Sequences (MACSE) is a multiple sequence alignment program that explicitly accounts for the underlying codon structure of protein-coding nucleotide sequences. Its unique characteristic allows building reliable codon alignments even in the presence of frameshifts. This facilitates downstream analyses such as selection pressure estimation based on the ratio of nonsynonymous to synonymous substitutions. Here, we present MACSE v2, a major update with an improved version of the initial algorithm enriched with a complete toolkit to handle multiple alignments of protein-coding sequences. A graphical interface now provides user-friendly access to the different subprograms.

The first complete mitogenome of Indian star tortoise (Geochelone elegans)

The complete mitochondrial genome of Indian star tortoise (Geochelone elegans) was characterized having 16,446 bp nucleotides encoding 37 genes in circular orientation comprising 13 protein-coding genes, 22 tRNA genes and two rRNA genes. The lengths of 12S and 16S ribosomal RNA were 973 bp and 1600 bp. A non-coding control region (D-Loop) of 966 bp was identified between tRNA^Pro and tRNA^Phe having seven interrupted tandem repeats. A single A + 1 frameshift insertion in the ND3 gene (ND3-174) was also discovered. The complete mitogenome of G. elegans would contribute in deeper understanding of the evolutionary dynamics and conservation effort of vulnerable testudine families.

Adaptive Patterns of Mitogenome Evolution Are Associated with the Loss of Shell Scutes in Turtles

The mitochondrial genome encodes several protein components of the oxidative phosphorylation (OXPHOS) pathway and is critical for aerobic respiration. These proteins have evolved adaptively in many taxa, but linking molecular-level patterns with higher-level attributes (e.g., morphology, physiology) remains a challenge. Turtles are a promising system for exploring mitochondrial genome evolution as different species face distinct respiratory challenges and employ multiple strategies for ensuring efficient respiration. One prominent adaptation to a highly aquatic lifestyle in turtles is the secondary loss of keratenized shell scutes (i.e., soft-shells), which is associated with enhanced swimming ability and, in some species, cutaneous respiration. We used codon models to examine patterns of selection on mitochondrial protein-coding genes along the three turtle lineages that independently evolved soft-shells. We found strong evidence for positive selection along the branches leading to the pig-nosed turtle (Carettochelys insculpta) and the softshells clade (Trionychidae), but only weak evidence for the leatherback (Dermochelys coriacea) branch. Positively selected sites were found to be particularly prevalent in OXPHOS Complex I proteins, especially subunit ND2, along both positively selected lineages, consistent with convergent adaptive evolution. Structural analysis showed that many of the identified sites are within key regions or near residues involved in proton transport, indicating that positive selection may have precipitated substantial changes in mitochondrial function. Overall, our study provides evidence that physiological challenges associated with adaptation to a highly aquatic lifestyle have shaped the evolution of the turtle mitochondrial genome in a lineage-specific manner.

Reviewing evidence for systematic transcriptional deletions, nucleotide exchanges, and expanded codons, and peptide clusters in human mitochondria

Polymerization sometimes transforms sequences by (a) systematic deletions of mono-, dinucleotides after trinucleotides, or (b) 23 systematic nucleotide exchanges (9 symmetric, X<>Y, e.g. G<>T, 14 asymmetric, X > Y > Z > X, e.g. A > G > T > A), producing del- and swinger RNAs. Some peptides correspond to del- and swinger RNA translations, also according to tetracodons, codons expanded by a silent nucleotide. Here new analyzes assume different proteolytic patterns, partially alleviating false negative peptide detection biases, expanding noncanonical mitoproteome profiles. Mito-genomic, -transcriptomic and -proteomic evidence for noncanonical transcriptions and translations are reviewed and clusters of del- and swinger peptides (also along tetracodons) are described. Noncanonical peptide clusters indicate regulated expression of cryptically encoded mitochondrial protein coding genes. These candidate noncanonical proteins don't resemble known proteins.

Two mitochondrial genomes in Alcedinidae (Ceryle rudis/Halcyon pileata) and the phylogenetic placement of Coraciiformes

Coraciiformes comprises 209 species belonging to ten families with significant divergence on external morphologies and life styles. The phylogenetic placement of Coraciiformes was still in debate. Here, we determined the complete mitochondrial genomes (mitogenomes) of Crested Kingfisher (Ceryle rudis) and Black-capped Kingfisher (Halcyon pileata). The mitogenomes were 17,355 bp (C. rudis) and 17,612 bp (H. pileata) in length, and both of them contained 37 genes (two rRNA genes, 22 tRNA genes and 13 protein-coding genes) and one control region. The gene organizations and characters of two mitogenomes were similar with those of other mitogenomes in Coraciiformes, however the sizes and nucleotide composition of control regions in different mitogenomes were significantly different. Phylogenetic trees were constructed with both Bayesian and Maximum Likelihood methods based on mitogenome sequences from 11 families of six orders. The trees based on two different data sets supported the basal position of Psittacidae (Psittaciformes), the closest relationship between Cuculiformes (Cuculidae) and Trogoniformes (Trogonidae), and the close relationship between Coraciiformes and Piciformes. The phylogenetic placement of the clade including Cuculiformes and Trogoniformes has not been resolved in present study, which need further investigations with more molecular markers and species. The mitogenome sequences presented here provided valuable data for further taxonomic studies on Coraciiformes and other related groups.

Resolution of the enigmatic phylogenetic relationship of the critically endangered Western Swamp Tortoise Pseudemydura umbrina (Pleurodira: Chelidae) using a complete mitochondrial genome

Pseudemydura umbrina is one of the most endangered turtle species in the world, and the imperative for its conservation is its distinctive morphology and relict status among the Chelidae. We use Illumina sequencing to obtain the complete mitogenome for resolving its uncertain phylogenetic position. A novel nuclear paralogue confounded the assembly, and resolution of the authentic mitogenome required further Sanger sequencing. The P. umbrina mitogenome is 16,414bp comprising 37 genes organized in a conserved pattern for other vertebrates. The nuclear paralogue is 547bp, 97.8% identity to the corresponding mitochondrial sequence. Particular features of the mitogenome include an nd3 174+1A frameshift, loss of DHC loop in tRNA^Ser (AGN), and a light-strand replication initiation site in Wancy region that extends into an adjacent tRNA gene. Phylogenetic analysis showed that P. umbrina is the monotypic sister lineage to the remaining Australasian Chelidae, a lineage probably dating back to the Cretaceous.

A novel mitochondrial genome of Arborophila and new insight into Arborophila evolutionary history

The lineage of the Bar-backed Partridge (Arborophila brunneopectus) was investigated to determine the phylogenetic relationships within Arborophila as the species is centrally distributed within an area covered by the distributions of 22 South-east Asian hill partridge species. The complete mitochondrial genome (mitogenome) of A. brunneopectus was determined and compared with four other hill partridge species mitogenomes. NADH subunit genes are radical in hill partridge mitogenomes and contain the most potential positive selective sites around where variable sites are abundant. Together with 44 other mitogenomes of closely related species, we reconstructed highly resolved phylogenetic trees using maximum likelihood (ML) and Bayesian inference (BI) analyses and calculated the divergence and dispersal history of Arborophila using combined datasets composed of their 13-protein coding sequences. Arborophila is reportedly be the oldest group in Phasianidae whose ancestors probably originated in Asia. A. rufipectus shares a closer relationship with A. ardens and A. brunneopectus compared to A. gingica and A. rufogularis, and such relationships were supported and profiled by NADH dehydrogenase subunit 5 (ND5). The intragenus divergence of all five Arborophila species occurred in the Miocene (16.84~5.69 Mya) when there were periods of climate cooling. We propose that these cooling events in the Miocene forced hill partridges from higher to lower altitudes, which led to geographic isolation and speciation. We demonstrated that the apparently deleterious +1 frameshift mutation in NADH dehydrogenase subunit 3 (ND3) found in all Arborophila is an ancient trait that has been eliminated in some younger lineages, such as Passeriformes. It is unclear of the biological advantages of this elimination for the relevant taxa and this requires further investigation.

Accurate annotation of protein-coding genes in mitochondrial genomes

Mitochondrial genome sequences are available in large number and new sequences become published nowadays with increasing pace. Fast, automatic, consistent, and high quality annotations are a prerequisite for downstream analyses. Therefore, we present an automated pipeline for fast de novo annotation of mitochondrial protein-coding genes. The annotation is based on enhanced phylogeny-aware hidden Markov models (HMMs). The pipeline builds taxon-specific enhanced multiple sequence alignments (MSA) of already annotated sequences and corresponding HMMs using an approximation of the phylogeny. The MSAs are enhanced by fixing unannotated frameshifts, purging of wrong sequences, and removal of non-conserved columns from both ends. A comparison with reference annotations highlights the high quality of the results. The frameshift correction method predicts a large number of frameshifts, many of which are unknown. A detailed analysis of the frameshifts in nad3 of the Archosauria-Testudines group has been conducted.

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.

Translation of mitochondrial swinger RNAs according to tri-, tetra- and pentacodons

Transcriptomes and proteomes include RNA and protein fragments not matching regular transcription/translation. Some 'non-canonical' mitochondrial transcripts match mitogenomes after assuming one among 23 systematic exchanges between nucleotides, producing swinger RNAs (nine symmetric, X↔Y, example C↔T; 14 asymmetric, X→Y→Z→X, example A→T→G→A) in GenBank's EST database. Here, reanalyzes of (a) public human mitochondrial transcriptome data (Illumina: RNA-seq) allowed to detect mitochondrial swinger RNAs for all 23 exchanges and (b) independent public human mitochondrial trypsinized proteomic mass spectrometry data allowed to detect peptides predicted from translation of parts of swinger-transformed mitogenomes covered by detected swinger reads. RNA-seq and previous EST swinger transcript data converge. Swinger RNA translation frequently inserts various amino acids at stop codons. Swinger RNA-peptide associations exist also for peptides matching systematically frameshifting translation, peptides entirely coded by tetra- and pentacodons (regular codons expanded by silent mononucleotides at 4th, and silent dinucleotides at 4th and 5th position(s), respectively). Swinger peptides differ from regular mitochondrial proteins: not membrane embedded, reflect warmer, anaerobic, low resource conditions, reminding a free-living ancestor. Tetra- and pentacoded peptides associate with low, high GC contents, respectively, suggesting expanded codon translations associate with thermic stresses. Results confirm experimentally predicted swinger, tetra- and pentacoded mitochondrial peptides, increasing mitogenomic coding density.

Mitochondrial genome diversity among six laboratory zebrafish (Danio rerio) strains

The mitochondrial genome (mtDNA) is distinct from the nuclear genome and is known to play a significant role in several disease phenotypes, such as longevity and fertility. Here we characterize the complete mitochondrial genomes (∼16 590 bp), and the extent of within and between strain variation for 27 adult zebrafish, representing five commonly used laboratory strains (AB, TL, HL, WIK, and SJD) and one line acquired from a local pet shop. These data were subsequently analyzed to determine the phylogenetic relationships between strains and ascertain if positive selection might be operating on any mtDNA genes. Relationships between strains are not entirely consistent with those observed previously using nuclear DNA. Further there is a substantial body of variation within current zebrafish lines, with 172 variants described across lines. Of these, 27 changes are non-synonymous and there is nominal evidence for positive selection in the mtDNA sequences at some of these sites. We further identify novel frameshift mutations in eight genes, which are all predicted to have functional consequences. Our study provides the first information on mtDNA diversity in zebrafish, identifies multiple non-synonymous substitutions and other mutations expected to have functional effects, and represents an important first step in establishing zebrafish as a model for investigating the phenotypic effects of mtDNA mutations.

Codon expansion and systematic transcriptional deletions produce tetra-, pentacoded mitochondrial peptides

Genes include occasionally isolated codons with a fourth (and fifth) silent nucleotide(s). Assuming tetracodons, translated hypothetical peptides align with regular GenBank proteins; predicted tetracodons coevolve with predicted tRNAs with expanded anticodons in each mammal, Drosophila and Lepidosauria mitogenomes, GC contents and with lepidosaurian body temperatures, suggesting that expanded codons are an adaptation of translation to high temperature. Hypothetically, continuous stretches of tetra- and pentacodons code for peptides. Both systematic nucleotide deletions during transcription, and translation by tRNAs with expanded anticodons could produce these peptides. Reanalyses of human nanoLc mass spectrometry peptidome data detect numerous tetra- and pentapeptides translated from the human mitogenome. These map preferentially on (BLAST-detected) human RNAs matching the human mitogenome, assuming systematic mono- and dinucleotide deletions after each third nucleotide (delRNAs). Translation by expanded anticodons is incompatible with silent nucleotides in the midst rather than at codon 3' extremity. More than 1/3 of detected tetra- and pentapeptides assume silent positions at codon extremity, suggesting that both mechanisms, regular translation of delRNAs and translation of regular RNAs by expanded anticodons, produce this peptide subgroup. Results show that systematically deleting polymerization occurs, and confirm serial translation of expanded codons. Non-canonical transcriptions and translations considerably expand the coding potential of DNA and RNA sequences.

The Mitochondrial Genomes of Aquila fasciata and Buteo lagopus (Aves, Accipitriformes): Sequence, Structure and Phylogenetic Analyses

The family Accipitridae is one of the largest groups of non-passerine birds, including 68 genera and 243 species globally distributed. In the present study, we determined the complete mitochondrial sequences of two species of accipitrid, namely Aquila fasciata and Buteo lagopus, and conducted a comparative mitogenome analysis across the family. The mitogenome length of A. fasciata and B. lagopus are 18,513 and 18,559 bp with an A + T content of 54.2% and 55.0%, respectively. For both the two accipitrid birds mtDNAs, obvious positive AT-skew and negative GC-skew biases were detected for all 12 PCGs encoded by the H strand, whereas the reverse was found in MT-ND6 encoded by the L strand. One extra nucleotide'C'is present at the position 174 of MT-ND3 gene of A. fasciata, which is not observed at that of B. lagopus. Six conserved sequence boxes in the Domain II, named boxes F, E, D, C, CSBa, and CSBb, respectively, were recognized in the CRs of A. fasciata and B. lagopus. Rates and patterns of mitochondrial gene evolution within Accipitridae were also estimated. The highest dN/dS was detected for the MT-ATP8 gene (0.32493) among Accipitridae, while the lowest for the MT-CO1 gene (0.01415). Mitophylogenetic analysis supported the robust monophyly of Accipitriformes, and Cathartidae was basal to the balance of the order. Moreover, we performed phylogenetic analyses using two other data sets (two mitochondrial loci, and combined nuclear and mitochondrial loci). Our results indicate that the subfamily Aquilinae and all currently polytypic genera of this subfamily are monophyletic. These two novel mtDNA data will be useful in refining the phylogenetic relationships and evolutionary processes of Accipitriformes.

The mitochondrial genome of the Saunders's gull Chroicocephalus saundersi (Charadriiformes: Laridae) and a higher phylogeny of shorebirds (Charadriiformes)

The complete mitogenome of Chroicocephalus saundersi was characterized and compared with the 6 published Charadriiformes mitogenomes. The mitogenome of C. saundersi is a closed circular molecule 16,739 bp in size, and contains 37 genes and a control region. The AT and GC skews are positive and negative, respectively, and in agreement with those of the other Charadriiformes mitogenomes. The mitogenome of C. saundersi contains 3 start codons (ATG, GTG, and ATT), 4 stop codons (TAA, TAG, AGG, and AGA), and an incomplete stop codon (T-) in 13 PCGs. A codon usage analysis of all available Charadriiformes mitogenomes showed that the ATG (78%) and TAA (50.5%) were the most common start codon and stop codon, respectively. An unusual start codon, ATT, is commonly found in the ND3s of Charadriiformes mitogenomes, whereas the more common start codons, ATC and ATA, are rarely found. In all the Laridae species, one extra cytosine was inserted at position 174 in ND3. The control region of C. saundersi is 1180-bp long, with a nucleotide composition of 30.2% A, 28.6% T, 27.3% C, and 14.0% G. Variable numbers of tandem repeats (VNTRs) with nine copies of the 10 bp repeat sequence (AACAACAAAC) are found within the CSB domain of the control region. The ML/BI analyses, based on the amino acids of the 13 mitochondrial PCGs, strongly support the monophyly of the order Charadriiformes, with the suborder Lari considered sister to the Scolopaci, which is in turn a sister group to the suborder Charadrii.

The complete mitochondrial genome sequence of the little egret (Egretta garzetta)

Many phylogenetic questions in the Ciconiiformes remain unresolved and complete mitogenome data are urgently needed for further molecular investigation. In this work, we determined the complete mitogenome sequence of the little egret (Egretta garzetta). The genome was 17,361 bp in length and the gene organization was typical of other avian mtDNA. In protein-coding genes (PCGs), a C insertion was found in ND3, and COIII and ND4 terminated with incomplete stop codons (T). tRNA-Val and tRNA-Ser (AGY) were unable to fold into canonical cloverleaf secondary structures because they had lost the DHU arms. Long repetitive sequences consisting of five types of tandem repeats were found at the 3' end of Domain III in the control region. A phylogenetic analysis of 11 species of Ciconiiformes was done using complete mitogenome data and 12 PCGs. The tree topologies obtained with these two strategies were identical, which strongly confirmed the monophyly of Ardeidae, Threskiorothidae and Ciconiidae. The phylogenetic analysis also revealed that Egretta was more closely related to Ardea than to Nycticorax in the Ardeidae, and Platalea was more closely related to Threskiornis than to Nipponia in the Threskiornithidae. These findings contribute to our understanding of the phylogenetic relationships of Ciconiiformes based on complete mitogenome data.

Complete mitochondrial genomes of the yellow-bellied slider turtle Trachemys scripta scripta and anoxia tolerant red-eared slider Trachemys scripta elegans

The complete mitochondrial genomes of the yellow-bellied slider (Trachemys scripta scripta) and anoxia tolerant red-eared slider (Trachemys scripta elegans) turtles were sequenced to analyze gene arrangement. The complete mt genomes of T. s. scripta and elegans were circular molecules of 16,791 bp and 16,810 bp in length, respectively, and included an A + 1 frameshift insertion in ND3 and ND4L genes. The AT content of the overall base composition of scripta and elegans was 61.2%. Nucleotide sequence divergence of the mt-genome (p distance) between scripta and elegans was 0.4%. A detailed comparison between the mitochondrial genomes of the two subspecies is shown.

Complete mitochondrial genome of Cuora trifasciata (Chinese three-striped box turtle), and a comparative analysis with other box turtles

Cuora trifasciata has become one of the most critically endangered species in the world. The complete mitochondrial genome of C. trifasciata (Chinese three-striped box turtle) was determined in this study. Its mitochondrial genome is a 16,575-bp-long circular molecule that consists of 37 genes that are typically found in other vertebrates. And the basic characteristics of the C. trifasciata mitochondrial genome were also determined. Moreover, a comparison of C. trifasciata with Cuora cyclornata, Cuora pani and Cuora aurocapitata indicated that the four mitogenomics differed in length, codons, overlaps, 13 protein-coding genes (PCGs), ND3, rRNA genes, control region, and other aspects. Phylogenetic analysis with Bayesian inference and maximum likelihood based on 12 protein-coding genes of the genus Cuora indicated the phylogenetic position of C. trifasciata within Cuora. The phylogenetic analysis also showed that C. trifasciata from Vietnam and China formed separate monophyletic clades with different Cuora species. The results of nucleotide base compositions, protein-coding genes and phylogenetic analysis showed that C. trifasciata from these two countries may represent different Cuora species.

The complete mitogenome of the Black-tailed gull Larus crassirostris (Charadriiformes: Laridae)

The complete mitogenome (KM507782) of the Black-tailed gull L. crassirostris was 16,746 bp long, with A 30.6% (5128 bp), T 24.3% (4076 bp), C 30.9% (5176 bp), and G 14.1% (2366 bp). Total length of 13 protein-coding genes was 11,396 bp and 12 of them except for ND6 were encoded in heavy strand. The three initiation codons ATG, GTG and ATT were used in the protein-coding genes. The ATG was a common initiation codon which was found in most of the protein-coding genes, whereas GTG was used in COX1 and ND5 and ATT in only ND3, respectively. The total length of 22 tRNA genes was 1550 bp, ranging from 66 bp (tRNA(Ser(AGY))) to 74 bp (tRNA(Leu(UUR)) and tRNA(Ser(UCN))).

The complete mitochondrial genome of western painted turtle, Chrysemys picta bellii (Chrysemys, Emydidae)

In this study, we report the complete mitochondrial genome sequence of western painted turtle, Chrysemys picta bellii. The genome is found to be 16,875 bp in length and has a base composition of A (34.4%), G (13.0%), C (26.0%), and T (26.6%). Similar to other turtles, it contains a typically conserved structure including 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 control region (D-loop). Most of the genes are encoded on H-strand, except for the eight tRNA and ND6 genes. All protein-coding genes start with an ATN codon except for COX1 and ND4, which initiate with GTG instead, and terminate with the typical stop codon (TAA/TAG) or a single T (T-) or an unexpected codon of AGG. The complete mitochondrial genome sequence provided here would be useful for further phylogenetic analysis and conservation genetic studies in C. p. bellii.

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

Three previously studied mitochondrial genomes of glass sponges (phylum Porifera, class Hexactinellida) contained single nucleotide insertions in protein coding genes inferred as sites of +1 translational frameshifting. To investigate the distribution and evolution of these sites and to help elucidate the mechanism of frameshifting, we determined eight new complete or nearly complete mtDNA sequences from glass sponges and examined individual mitochondrial genes from three others. We found nine new instances of single nucleotide insertions in these sequences and analyzed them both comparatively and phylogenetically. The base insertions appear to have been gained and lost repeatedly in hexactinellid mt protein genes, suggesting no functional significance for the frameshifting sites. A high degree of sequence conservation, the presence of unusual tRNAs, and a distinct pattern of codon usage suggest the "out-of-frame pairing" model of translational frameshifting. Additionally, we provide evidence that relaxed selection pressure on glass sponge mtDNA - possibly a result of their low growth rates and deep-water lifestyle - has allowed frameshift insertions to be tolerated for hundreds of millions of years. Our study provides the first example of a phylogenetically diverse and extensive usage of translational frameshifting in animal mitochondrial coding sequences.