Possible multiple origins of replication in primate mitochondria: Alternative role of tRNA sequences

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational ‘hot-spots’ or ‘cold-spots’. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Evolutionary History of Inversions in Directional Mutational Pressures in Crustacean Mitochondrial Genomes: Implications for Evolutionary Studies

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.

Comparative Mitogenomic Analysis of Heptageniid Mayflies (Insecta: Ephemeroptera): Conserved Intergenic Spacer and tRNA Gene Duplication

Large intergenic spacers and tRNA gene duplications have been reported in several insect groups, although little is known about mitogenomes of mayflies. Here, we determined complete mitogenomes of ten heptageniid species and systemically analyzed their mitogenomic features. Both a conserved intergenic spacer (IGS) and trnM duplication were detected in those mitogenomes. The IGS, which was observed in heptageniids, could be further folded into a stable stem-loop structure. The tRNA gene duplication was found in almost all analyzed mitogenomes, and a unique gene block trnI-trnM-trnQ-trnM-ND2 was also discovered. Our analysis demonstrates that the heptageniid gene arrangement pattern can be explained by the tandem duplication-random loss (TDRL) model. Phylogenetic analyses using both Bayesian inference (BI) and maximum likelihood (ML) methods based on the nucleotide and amino acid sequence data recovered the genus Epeorus as monophyletic with strong support. Our results provide a better understanding of mitogenomic evolution in Heptageniidae, as well as novel molecular markers for species identification of mayflies.

Negative CG dinucleotide bias: An explanation based on feedback loops between Arginine codon assignments and theoretical minimal RNA rings

Theoretical minimal RNA rings are candidate primordial genes evolved for non-redundant coding of the genetic code's 22 coding signals (one codon per biogenic amino acid, a start and a stop codon) over the shortest possible length: 29520 22-nucleotide-long RNA rings solve this min-max constraint. Numerous RNA ring properties are reminiscent of natural genes. Here we present analyses showing that all RNA rings lack dinucleotide CG (a mutable, chemically instable dinucleotide coding for Arginine), bearing a resemblance to known CG-depleted genomes. CG in "incomplete" RNA rings (not coding for all coding signals, with only 3-12 nucleotides) gradually decreases towards CG absence in complete, 22-nucleotide-long RNA rings. Presumably, feedback loops during RNA ring growth during evolution (when amino acid assignment fixed the genetic code) assigned Arg to codons lacking CG (AGR) to avoid CG. Hence, as a chemical property of base pairs, CG mutability restructured the genetic code, thereby establishing itself as genetically encoded biological information.

Deamination gradients within codons after 1<->2 position swap predict amino acid hydrophobicity and parallel β-sheet conformational preference

Deaminations C->T and A->G are frequent mutations producing nucleotide content gradients across genomes proportional to singlestrandedness during replication/transcription. Hence, within single codons, deamination risks increase from first to third codon positions, while second codon positions are functionally most crucial. Here genetic codes are analyzed assuming that after anticodons protected codons from deaminations, first and second codon positions swapped (N2N1N3->N1N2N3), with lowest deamination risks for N2 in presumed primitive N2N1N3 codons. N2N1N3, not standard N1N2N3, codon structure minimizes deaminations inversely proportionally to cognate amino acid hydrophobicity and parallel betasheet conformational preference. For N1N2N3, deamination minimization increases with genetic code integration order of cognate amino acids: during the presumed N2N1N3->N1N2N3 codon structure transition, protein synthesis combined direct codon-amino acid interactions for late amino acids and tRNA-based translation for early amino acids. Hence N2N1N3 codons would correspond to tRNA-free translation by spontaneous codon-amino acid affinities, and tRNA-mediated translation presumably caused N2N1N3->N1N2N3 swaps. Results show that rational, not arbitrary rules link codon and amino acid structures. Some analyses detect mitochondrial RNAs and peptides in public data corresponding to systematic position swaps, suggesting occasional swapping polymerase activity.

RNA Rings Strengthen Hairpin Accretion Hypotheses for tRNA Evolution: A Reply to Commentaries by Z.F. Burton and M. Di Giulio

Theoretical minimal RNA ring design ensures coding over the shortest length once for each coding signal (start and stop codons, and each amino acid) and their hairpin configuration. These constraints define 25 RNA rings which surprisingly resemble ancestral tRNA loops, suggesting commonalities between RNA ring design and proto-tRNAs. RNA rings share several other properties with tRNAs, suggesting that primordial RNAs were multifunctional peptide coding sequences and structural RNAs. Two hypotheses, respectively, by M. Di Giulio and Z.F. Burton, derived from cloverleaf structural symmetries suggest that two and three, respectively, stem-loop hairpins agglutinated into tRNAs. Their authors commented that their respective structure-based hypotheses reflect better tRNA structure than RNA rings. Unlike these hypotheses, RNA ring design uses no tRNA-derived information, rendering model predictive power comparisons senseless. Some analyses of RNA ring primary and secondary structures stress RNA ring splicing in their predicted anticodon's midst, indicating ancestrality of split tRNAs, as the two-piece model predicts. Advancement of knowledge, rather than of specific hypotheses, gains foremost by examining independent hypotheses for commonalities, and only secondarily for discordances. RNA rings mimick ancestral biomolecules including tRNAs, and their evolution, and constitute an interesting synthetic system for early prebiotic evolution tests/simulations.

The primordial tRNA acceptor stem code from theoretical minimal RNA ring clusters

BACKGROUND

Theoretical minimal RNA rings code by design over the shortest length once for each of the 20 amino acids, a start and a stop codon, and form stem-loop hairpins. This defines at most 25 RNA rings of 22 nucleotides. As a group, RNA rings mimick numerous prebiotic and early life biomolecular properties: tRNAs, deamination gradients and replication origins, emergence of codon preferences for the natural circular code, and contents of early protein coding genes. These properties result from the RNA ring's in silico design, based mainly on coding nonredundancy among overlapping translation frames, as the genetic code's codon-amino acid assignments determine. RNA rings resemble ancestral tRNAs, defining RNA ring anticodons and corresponding cognate amino acids. Surprisingly, all examined RNA ring properties coevolve with genetic code integration ranks of RNA ring cognates, as if RNA rings mimick prebiotic and early life evolution.

METHODS

Distances between RNA rings were calculated using different evolutionary models. Associations between these distances and genetic code evolutionary hypotheses detect evolutionary models best describing RNA ring diversification.

RESULTS

Here pseudo-phylogenetic analyses of RNA rings produce clusters corresponding to the primordial code in tRNA acceptor stems, more so when substitution matrices from neutrally evolving pseudogenes are used rather than from functional protein coding genes reflecting selection for conserving amino acid properties.

CONCLUSIONS

Results indicate RNA rings with recent cognates evolved from those with early cognates. Hence RNA rings, as designed by the genetic code's structure, simulate tRNA stem evolution and prebiotic history along neutral chemistry-driven mutation regimes.

The Uroboros Theory of Life's Origin: 22-Nucleotide Theoretical Minimal RNA Rings Reflect Evolution of Genetic Code and tRNA-rRNA Translation Machineries

Theoretical minimal RNA rings attempt to mimick life's primitive RNAs. At most 25 22-nucleotide-long RNA rings code once for each biotic amino acid, a start and a stop codon and form a stem-loop hairpin, resembling consensus tRNAs. We calculated, for each RNA ring's 22 potential splicing positions, similarities of predicted secondary structures with tRNA vs. rRNA secondary structures. Assuming rRNAs partly derived from tRNA accretions, we predict positive associations between relative secondary structure similarities with rRNAs over tRNAs and genetic code integration orders of RNA ring anticodon cognate amino acids. Analyses consider for each secondary structure all nucleotide triplets as potential anticodon. Anticodons for ancient, chemically inert cognate amino acids are most frequent in the 25 RNA rings. For RNA rings with primordial cognate amino acids according to tRNA-homology-derived anticodons, tRNA-homology and coding sequences coincide, these are separate for predicted cognate amino acids that presumably integrated late the genetic code. RNA ring secondary structure similarity with rRNA over tRNA secondary structures associates best with genetic code integration orders of anticodon cognate amino acids when assuming split anticodons (one and two nucleotides at the spliced RNA ring 5' and 3' extremities, respectively), and at predicted anticodon location in the spliced RNA ring's midst. Results confirm RNA ring homologies with tRNAs and CDs, ancestral status of tRNA half genes split at anticodons, the tRNA-rRNA axis of RNA evolution, and that single theoretical minimal RNA rings potentially produce near-complete proto-tRNA sets. Hence genetic code pre-existence determines 25 short circular gene- and tRNA-like RNAs. Accounting for each potential splicing position, each RNA ring potentially translates most amino acids, realistically mimicks evolution of the tRNA-rRNA translation machinery. These RNA rings 'of creation' remind the uroboros' (snake biting its tail) symbolism for creative regeneration.

Mechanisms of gene rearrangement in 13 bothids based on comparison with a newly completed mitogenome of the threespot flounder, Grammatobothus polyophthalmus (Pleuronectiformes: Bothidae)

Background

The mitogenomes of 12 teleost fish of the bothid family (order Pleuronectiformes) indicated that the genomic-scale rearrangements characterized in previous work. A novel mechanism of genomic rearrangement called the Dimer-Mitogenome and Non-Random Loss (DMNL) model was used to account for the rearrangement found in one of these bothids, Crossorhombus azureus.

Results

The 18,170 bp mitogenome of G. polyophthalmus contains 37 genes, two control regions (CRs), and the origin of replication of the L-strand (O_L). This mitogenome is characterized by genomic-scale rearrangements: genes located on the L-strand are grouped in an 8-gene cluster (Q-A-C-Y-S₁-ND6-E-P) that does not include tRNA-N; genes found on the H-strand are grouped together (F-12S … CytB-T) except for tRNA-D that was translocated inside the 8-gene L-strand cluster. Compared to non-rearranged mitogenomes of teleost fishes, gene organization in the mitogenome of G. polyophthalmus and in that of the other 12 bothids characterized thus far is very similar. These rearrangements could be sorted into four types (Type I, II, III and IV), differing in the particular combination of the CR, tRNA-D gene and 8-gene cluster and the shuffling of tRNA-V. The DMNL model was used to account for all but one gene rearrangement found in all 13 bothid mitogenomes. Translocation of tRNA-D most likely occurred after the DMNL process in 10 bothid mitogenomes and could have occurred either before or after DMNL in the three other species. During the DMNL process, the tRNA-N gene was retained rather than the expected tRNA-N′ gene. tRNA-N appears to assist in or act as O_L function when the O_L secondary structure could not be formed from intergenic sequences. A striking finding was that each of the non-transcribed genes has degenerated to a shorter intergenic spacer during the DMNL process. These findings highlight a rare phenomenon in teleost fish.

Conclusions

This result provides significant evidence to support the existence of dynamic dimeric mitogenomes and the DMNL model as the mechanism of gene rearrangement in bothid mitogenomes, which not only promotes the understanding of mitogenome structural diversity, but also sheds light on mechanisms of mitochondrial genome rearrangement and replication.

Characterization of the Complete Mitochondrial Genome of Harpalus sinicus and Its Implications for Phylogenetic Analyses

In this study, we report the complete mitochondrial genome of Harpalus sinicus (occasionally named as the Chinese ground beetle) which is the first mitochondrial genome for Harpalus. The mitogenome is 16,521 bp in length, comprising 37 genes, and a control region. The A + T content of the mitogenome is as high as 80.6%. A mitochondrial origins of light-strand replication (OL)-like region is found firstly in the insect mitogenome, which can form a stem-loop hairpin structure. Thirteen protein-coding genes (PCGs) share high homology, and all of them are under purifying selection. All tRNA genes (tRNAs) can be folded into the classic cloverleaf secondary structures except tRNA-Ser (GCU), which lacks a dihydrouridine (DHU) stem. The secondary structure of two ribosomal RNA genes (rRNAs) is predicted based on previous insect models. Twelve types of tandem repeats and two stem-loop structures are detected in the control region, and two stem-loop structures may be involved in the initiation of replication and transcription. Additionally, phylogenetic analyses based on mitogenomes suggest that Harpalus is an independent lineage in Carabidae, and is closely related to four genera (Abax, Amara, Stomis, and Pterostichus). In general, this study provides meaningful genetic information for Harpalus sinicus and new insights into the phylogenetic relationships within the Carabidae.

Localized Context-Dependent Effects of the "Ambush" Hypothesis: More Off-Frame Stop Codons Downstream of Shifty Codons

The ambush hypothesis speculates that off-frame stop codons increase translational efficiency after ribosomal frameshifts by stopping early frameshifted translation. Some evidences fit this hypothesis: (1) synonymous codon usages increase with their potential contribution to off-frame stops; (2) the genetic code assigns frequent amino acids to codon families contributing to off-frame stops; (3) positive biases for off-frame stops (AT rich) occur despite adverse nucleotide (GC) biases; and (4) mitochondrial off-frame stop codon densities increase with ribosomal structural instability, potential proxy of frameshift frequencies. In this study, analyses of vertebrate mitogenes and tRNA synthetase genes from all superkingdoms and viruses test a new prediction of the ambush hypothesis: sequences immediately downstream of frameshift-inducing homopolymer codons (AAA, CCC, GGG, and TTT) are off-frame stop rich. Codons immediately downstream of homopolymer codons form more than average off-frame stops, biases are stronger than for corresponding upstream distances and for any other group of synonymous codons. Sequences downstream of that high-density region are off-frame stop depleted. This decrease suggests that off-frame stops, combined with suppressor tRNAs regulate translation of overlapping coding sequences. Results show the predictive power of the ambush hypothesis, from macroevolutionary (genetic code structure) to detailed gene sequence anatomy.

Protein Sequences Recapitulate Genetic Code Evolution

Several hypotheses predict ranks of amino acid assignments to genetic code's codons. Analyses here show that average positions of amino acid species in proteins correspond to assignment ranks, in particular as predicted by Juke's neutral mutation hypothesis for codon assignments. In all tested protein groups, including co- and post-translationally folding proteins, 'recent' amino acids are on average closer to gene 5' extremities than 'ancient' ones. Analyses of pairwise residue contact energies matrices suggest that early amino acids stereochemically selected late ones that stablilize residue interactions within protein cores, presumably producing 5'-late-to-3'-early amino acid protein sequence gradients. The gradient might reduce protein misfolding, also after mutations, extending principles of neutral mutations to protein folding. Presumably, in self-perpetuating and self-correcting systems like the genetic code, initial conditions produce similarities between evolution of the process (the genetic code) and 'ontogeny' of resulting structures (here proteins), producing apparent teleonomy between process and product.

The GC Skew Index: A Measure of Genomic Compositional Asymmetry and the Degree of Replicational Selection

Circular bacterial chromosomes have highly polarized nucleotide composition in the two replichores, and this genomic strand asymmetry can be visualized using GC skew graphs. Here we propose and discuss the GC skew index (GCSI) for the quantification of genomic compositional skew, which combines a normalized measure of fast Fourier transform to capture the shape of the skew graph and Euclidean distance between the two vertices in a cumulative skew graph to represent the degree of skew. We calculated GCSI for all available bacterial genomes, and GCSI correlated well with the visibility of GC skew. This novel index is useful for estimating confidence levels for the prediction of replication origin and terminus by methods based on GC skew and for measuring the strength of replicational selection in a genome.

Genetic Code Optimization for Cotranslational Protein Folding: Codon Directional Asymmetry Correlates with Antiparallel Betasheets, tRNA Synthetase Classes

A new codon property, codon directional asymmetry in nucleotide content (CDA), reveals a biologically meaningful genetic code dimension: palindromic codons (first and last nucleotides identical, codon structure XZX) are symmetric (CDA = 0), codons with structures ZXX/XXZ are 5'/3' asymmetric (CDA = - 1/1; CDA = - 0.5/0.5 if Z and X are both purines or both pyrimidines, assigning negative/positive (-/+) signs is an arbitrary convention). Negative/positive CDAs associate with (a) Fujimoto's tetrahedral codon stereo-table; (b) tRNA synthetase class I/II (aminoacylate the 2'/3' hydroxyl group of the tRNA's last ribose, respectively); and (c) high/low antiparallel (not parallel) betasheet conformation parameters. Preliminary results suggest CDA-whole organism associations (body temperature, developmental stability, lifespan). Presumably, CDA impacts spatial kinetics of codon-anticodon interactions, affecting cotranslational protein folding. Some synonymous codons have opposite CDA sign (alanine, leucine, serine, and valine), putatively explaining how synonymous mutations sometimes affect protein function. Correlations between CDA and tRNA synthetase classes are weaker than between CDA and antiparallel betasheet conformation parameters. This effect is stronger for mitochondrial genetic codes, and potentially drives mitochondrial codon-amino acid reassignments. CDA reveals information ruling nucleotide-protein relations embedded in reversed (not reverse-complement) sequences (5'-ZXX-3'/5'-XXZ-3').

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.

Evolution of Nucleotide Punctuation Marks: From Structural to Linear Signals

We present an evolutionary hypothesis assuming that signals marking nucleotide synthesis (DNA replication and RNA transcription) evolved from multi- to unidimensional structures, and were carried over from transcription to translation. This evolutionary scenario presumes that signals combining secondary and primary nucleotide structures are evolutionary transitions. Mitochondrial replication initiation fits this scenario. Some observations reported in the literature corroborate that several signals for nucleotide synthesis function in translation, and vice versa. (a) Polymerase-induced frameshift mutations occur preferentially at translational termination signals (nucleotide deletion is interpreted as termination of nucleotide polymerization, paralleling the role of stop codons in translation). (b) Stem-loop hairpin presence/absence modulates codon-amino acid assignments, showing that translational signals sometimes combine primary and secondary nucleotide structures (here codon and stem-loop). (c) Homopolymer nucleotide triplets (AAA, CCC, GGG, TTT) cause transcriptional and ribosomal frameshifts. Here we find in recently described human mitochondrial RNAs that systematically lack mono-, dinucleotides after each trinucleotide (delRNAs) that delRNA triplets include 2x more homopolymers than mitogenome regions not covered by delRNA. Further analyses of delRNAs show that the natural circular code X (a little-known group of 20 translational signals enabling ribosomal frame retrieval consisting of 20 codons {AAC, AAT, ACC, ATC, ATT, CAG, CTC, CTG, GAA, GAC, GAG, GAT, GCC, GGC, GGT, GTA, GTC, GTT, TAC, TTC} universally overrepresented in coding versus other frames of gene sequences), regulates frameshift in transcription and translation. This dual transcription and translation role confirms for X the hypothesis that translational signals were carried over from transcriptional signals.

Unbiased Mitoproteome Analyses Confirm Non-canonical RNA, Expanded Codon Translations

Proteomic MS/MS mass spectrometry detections are usually biased towards peptides cleaved by experimentally added digestion enzyme(s). Hence peptides resulting from spontaneous degradation and natural proteolysis usually remain undetected. Previous analyses of tryptic human proteome data (cleavage after K, R) detected non-canonical tryptic peptides translated according to tetra- and pentacodons (codons expanded by silent mono- and dinucleotides), and from transcripts systematically (a) deleting mono-, dinucleotides after trinucleotides (delRNAs), (b) exchanging nucleotides according to 23 bijective transformations. Nine symmetric and fourteen asymmetric nucleotide exchanges (X ↔ Y, e.g. A ↔ C; and X → Y → Z → X, e.g. A → C → G → A) produce swinger RNAs. Here unbiased reanalyses of these proteomic data detect preferentially non-canonical tryptic peptides despite assuming random cleavage. Unbiased analyses couldn't reconstruct experimental tryptic digestion if most detected non-canonical peptides were false positives. Detected non-tryptic non-canonical peptides map preferentially on corresponding, previously described non-canonical transcripts, as for tryptic non-canonical peptides. Hence unbiased analyses independently confirm previous trypsin-biased analyses that showed translations of del- and swinger RNA and expanded codons. Accounting for natural proteolysis completes trypsin-biased mitopeptidome analyses, independently confirms non-canonical transcriptions and translations.

Chimeric mitochondrial peptides from contiguous regular and swinger RNA

Previous mass spectrometry analyses described human mitochondrial peptides entirely translated from swinger RNAs, RNAs where polymerization systematically exchanged nucleotides. Exchanges follow one among 23 bijective transformation rules, nine symmetric exchanges (X ↔ Y, e.g. A ↔ C) and fourteen asymmetric exchanges (X → Y → Z → X, e.g. A → C → G → A), multiplying by 24 DNA's protein coding potential. Abrupt switches from regular to swinger polymerization produce chimeric RNAs. Here, human mitochondrial proteomic analyses assuming abrupt switches between regular and swinger transcriptions, detect chimeric peptides, encoded by part regular, part swinger RNA. Contiguous regular- and swinger-encoded residues within single peptides are stronger evidence for translation of swinger RNA than previously detected, entirely swinger-encoded peptides: regular parts are positive controls matched with contiguous swinger parts, increasing confidence in results. Chimeric peptides are 200 × rarer than swinger peptides (3/100,000 versus 6/1000). Among 186 peptides with > 8 residues for each regular and swinger parts, regular parts of eleven chimeric peptides correspond to six among the thirteen recognized, mitochondrial protein-coding genes. Chimeric peptides matching partly regular proteins are rarer and less expressed than chimeric peptides matching non-coding sequences, suggesting targeted degradation of misfolded proteins. Present results strengthen hypotheses that the short mitogenome encodes far more proteins than hitherto assumed. Entirely swinger-encoded proteins could exist.

•

Chimeric peptides are translated from contiguous regular and swinger RNA

•

They are 200x rarer than mitochondrial swinger peptides

•

Chimeric peptides integrated in regular mitochondrial proteins are downregulated

•

Contiguous regular parts are matched positive controls for swinger parts

•

The last point validates results beyond other statistical tests for robustness

The Complete Female- and Male-Transmitted Mitochondrial Genome of Meretrix lamarckii

Bivalve mitochondrial genomes show many uncommon features, like additional genes, high rates of gene rearrangement, high A-T content. Moreover, Doubly Uniparental Inheritance (DUI) is a distinctive inheritance mechanism allowing some bivalves to maintain and transmit two separate sex-linked mitochondrial genomes. Many bivalve mitochondrial features, such as gene extensions or additional ORFs, have been proposed to be related to DUI but, up to now, this topic is far from being understood. Several species are known to show this unusual organelle inheritance but, being widespread only among Unionidae and Mytilidae, DUI distribution is unclear. We sequenced and characterized the complete female- (F) and male-transmitted (M) mitochondrial genomes of Meretrix lamarckii, which, in fact, is the second species of the family Veneridae where DUI has been demonstrated so far. The two mitochondrial genomes are comparable in length and show roughly the same gene content and order, except for three additional tRNAs found in the M one. The two sex-linked genomes show an average nucleotide divergence of 16%. A 100-aminoacid insertion in M. lamarckii M-cox2 gene was found; moreover, additional ORFs have been found in both F and M Long Unassigned Regions of M. lamarckii. Even if no direct involvement in DUI process has been demonstrated so far, the finding of cox2 insertions and supernumerary ORFs in M. lamarckii both strengthens this hypothesis and widens the taxonomical distribution of such unusual features. Finally, the analysis of inter-sex genetic variability shows that DUI species form two separate clusters, namely Unionidae and Mytilidae+Veneridae; this dichotomy is probably due to different DUI regimes acting on separate taxa.

Swinger RNA self-hybridization and mitochondrial non-canonical swinger transcription, transcription systematically exchanging nucleotides

Stem-loop hairpins punctuate mitochondrial post-transcriptional processing. Regulation of mitochondrial swinger transcription, transcription producing RNAs matching the mitogenome only assuming systematic exchanges between nucleotides (23 bijective transformations along 9 symmetric exchanges X<>Y, e.g. A<>G, and 14 asymmetric exchanges X>Y>Z>X, e.g. A>G>C>A) remains unknown. Does swinger RNA self-hybridization regulate swinger, as regular, transcription? Groups of 8 swinger transformations share canonical self-hybridization properties within each group, group 0 includes identity (regular) transcription. The human mitogenome has more stem-loop hairpins than randomized sequences for all groups. Group 2 transformations reveal complementarity of the light strand replication origin (OL) loop and a neighboring tRNA gene, detecting the longtime presumed OL/tRNA homology. Non-canonical G=U pairings in hairpins increases with swinger RNA detection. These results confirm biological relevancy of swinger-transformed DNA/RNA, independently of, and in combination with, previously detected swinger DNA/RNA and swinger peptides. Swinger-transformed mitogenomes include unsuspected multilayered information.

Systematically frameshifting by deletion of every 4th or 4th and 5th nucleotides during mitochondrial transcription: RNA self-hybridization regulates delRNA expression

In mitochondria, secondary structures punctuate post-transcriptional RNA processing. Recently described transcripts match the human mitogenome after systematic deletions of every 4th, respectively every 4th and 5th nucleotides, called delRNAs. Here I explore predicted stem-loop hairpin formation by delRNAs, and their associations with delRNA transcription and detected peptides matching their translation. Despite missing 25, respectively 40% of the nucleotides in the original sequence, del-transformed sequences form significantly more secondary structures than corresponding randomly shuffled sequences, indicating biological function, independently of, and in combination with, previously detected delRNA and thereof translated peptides. Self-hybridization decreases delRNA abundances, indicating downregulation. Systematic deletions of the human mitogenome reveal new, unsuspected coding and structural informations.

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.

Systematic exchanges between nucleotides: Genomic swinger repeats and swinger transcription in human mitochondria

Chargaff׳s second parity rule, quasi-equal single strand frequencies for complementary nucleotides, presumably results from insertion of repeats and inverted repeats during sequence genesis. Vertebrate mitogenomes escape this rule because repeats are counterselected: their hybridization produces loop bulges whose deletion is deleterious. Some DNA/RNA sequences match mitogenomes only after assuming one among 23 systematic nucleotide exchanges (swinger DNA/RNA: nine symmetric, e.g. A ↔ C; and 14 asymmetric, e.g. A → C → G → A). Swinger-transformed repeats do not hybridize, escaping selection against deletions due to bulge formation. Blast analyses of the human mitogenome detect swinger repeats for all 23 swinger types, more than in randomized sequences with identical length and nucleotide contents. Mean genomic swinger repeat lengths increase with observed human swinger RNA frequencies: swinger repeat and swinger RNA productions appear linked, perhaps by swinger RNA retrotranscription. Mean swinger repeat lengths are proportional to reading frame retrievability, post-swinger transformation, by the natural circular code. Genomic swinger repeats confirm at genomic level, independently of swinger RNA detection, occurrence of swinger polymerizations. They suggest that repeats, and swinger repeats in particular, contribute to genome genesis.

Tandem Duplication and Random Loss for mitogenome rearrangement in Symphurus (Teleost: Pleuronectiformes)

Background

The mitochondrial genomes (mitogenomes) of flatfishes (Pleuronectiformes) exhibit highly diversified types of large-scale gene rearrangements. We have reported that the mitogenomes of Crossorhombus azureus (Bothidae), Samariscus latus (Samaridae) and Cynoglossus fishes (Cynoglossidae) show different types of gene rearrangements.

Results

In the present study, the complete mitogenomes of two Symphurus species (Cynoglossidae), Symphurus plagiusa and Symphurus orientalis, were determined. The gene order in the S. plagiusa mitogenome is the same as that of a typical vertebrate (without any gene rearrangements). Surprisingly, large-scale gene rearrangements have occurred in S. orientalis. In the rearranged fragment from the control region (CR) to the WANCY tRNA cluster (tRNA cluster of tRNA-W, tRNA-A, tRNA-N, tRNA-C and tRNA-Y) in the S. orientalis mitogenome, tRNA-V and tRNA-M have been translocated to the 3’ end of the 16S rRNA gene, with six large intergenic spacers over 20 bp in length. In addition, an origin for light-strand replication (O_L) structure that is typically located in the WANCY region was absent in both the S. plagiusa and S. orientalis mitogenomes. It is generally recognized that a sequence in the WANCY region that encodes tRNAs forms a hairpin structure (O_L-like structure) and can act as the O_L when the typical locus is lost. Moreover, an additional O_L-like structure was identified near the control region in the S. plagiusa mitogenome.

Conclusions

The positions of the intergenic spacers and the rearranged genes of the S. orientalis mitogenome strongly indicate that the mechanism underlying the rearrangement of this mitogenome was Tandem Duplication and Random Loss. Additionally, two O_L-like regions substituting for the typical locus were found in the S. plagiusa mitogenome. We speculate that the ancestral mitogenomes of S. plagiusa and S. orientalis also had this characteristic, such that if both O_L-like structures functioned during mitochondrial replication, they could initiate duplicate replications of the light strand (L-strand), leading to duplication of the region between the two structures. We consider that this mechanism may account for the gene duplication that occurred during the gene rearrangement process in the evolution of the ancestral mitogenome to the S. orientalis mitogenome.

Electronic supplementary material

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

Mitochondrial swinger replication: DNA replication systematically exchanging nucleotides and short 16S ribosomal DNA swinger inserts

Assuming systematic exchanges between nucleotides (swinger RNAs) resolves genomic 'parenthood' of some orphan mitochondrial transcripts. Twenty-three different systematic nucleotide exchanges (bijective transformations) exist. Similarities between transcription and replication suggest occurrence of swinger DNA. GenBank searches for swinger DNA matching the 23 swinger versions of human and mouse mitogenomes detect only vertebrate mitochondrial swinger DNA for swinger type AT+CG (from five different studies, 149 sequences) matching three human and mouse mitochondrial genes: 12S and 16S ribosomal RNAs, and cytochrome oxidase subunit I. Exchange A<->T+C<->G conserves self-hybridization properties, putatively explaining swinger biases for rDNA, against protein coding genes. Twenty percent of the regular human mitochondrial 16S rDNA consists of short swinger repeats (from 13 exchanges). Swinger repeats could originate from recombinations between regular and swinger DNA: duplicated mitochondrial genes of the parthenogenetic gecko Heteronotia binoei include fewer short A<->T+C<->G swinger repeats than non-duplicated mitochondrial genomes of that species. Presumably, rare recombinations between female and male mitochondrial genes (and in parthenogenetic situations between duplicated genes), favors reverse-mutations of swinger repeat insertions, probably because most inserts affect negatively ribosomal function. Results show that swinger DNA exists, and indicate that swinger polymerization contributes to the genesis of genetic material and polymorphism.

The inversion of the Control Region in three mitogenomes provides further evidence for an asymmetric model of vertebrate mtDNA replication

Mitochondrial genomes are known to have a strong strand-specific compositional bias that is more pronounced at fourfold redundant sites of mtDNA protein-coding genes. This observation suggests that strand asymmetries, to a large extent, are caused by mutational asymmetric mechanisms. In vertebrate mitogenomes, replication and not transcription seems to play a major role in shaping compositional bias. Hence, one can better understand how mtDNA is replicated – a debated issue – through a detailed picture of mitochondrial genome evolution. Here, we analyzed the compositional bias (AT and GC skews) in protein-coding genes of almost 2,500 complete vertebrate mitogenomes. We were able to identify three fish mitogenomes with inverted AT/GC skew coupled with an inversion of the Control Region. These findings suggest that the vertebrate mitochondrial replication mechanism is asymmetric and may invert its polarity, with the leading-strand becoming the lagging-strand and vice-versa, without compromising mtDNA maintenance and expression. The inversion of the strand-specific compositional bias through the inversion of the Control Region is in agreement with the strand-displacement model but it is also compatible with the RITOLS model of mtDNA replication.

Species radiation by DNA replication that systematically exchanges nucleotides?

RNA and DNA syntheses share many properties. Therefore, the existence of 'swinger' RNAs, presumed 'orphan' transcripts matching genomic sequences only if transcription systematically exchanged nucleotides, suggests replication producing swinger DNA. Transcripts occur in many short-lived copies, the few cellular DNA molecules are long-lived. Hence pressures for functional swinger DNAs are greater than for swinger RNAs. Protein coding properties of swinger sequences differ from original sequences, suggesting rarity of corresponding swinger DNA. For genes producing structural RNAs, such as tRNAs and rRNAs, three exchanges (A<->T, C<->G and A<->T+C<->G) conserve self-hybridization properties. All nuclear eukaryote swinger DNA sequences detected in GenBank are for rRNA genes assuming A<->T+C<->G exchanges. In brachyuran crabs, 25 species had A<->T+C<->G swinger 18S rDNA, all matching the reverse-exchanged version of regular 18S rDNA of a related species. In this taxon, swinger replication of 18S rDNA apparently associated with, or even resulted in species radiation. A<->T+C<->G transformation doesn't invert sequence direction, differing from inverted repeats. Swinger repeats (detectable only assuming swinger transformations, A<->T+C<->G swinger repeats most frequent) within regular human rRNAs, independently confirm swinger polymerizations for most swinger types. Swinger replication might be an unsuspected molecular mechanism for ultrafast speciation.

Putative anticodons in mitochondrial tRNA sidearm loops: Pocketknife tRNAs?

The hypothesis that tRNA sidearm loops bear anticodons assumes crossovers between anticodon and sidearms, or translation by expressed aminoacylated tRNA halves forming single stem-loops. Only the latter might require ribosomal adaptations. Drosophila mitochondrial codon usages coevolve with sidearm numbers bearing matching putative anticodons (comparing different codon families in one genome, macroevolution) and when comparing different genomes for single codon families (microevolution). Coevolution between Drosophila and yeast mitochondrial antisense tRNAs and codon usages partly confounds microevolutionary patterns for putative sidearm anticodons. Some tRNA sidearm loops have more than seven nucleotides, putative expanded anticodons potentially matching quadruplet codons (tetracodons, codons expanded by a fourth silent position, forming tetragenes (predicted by alignment analyses of Drosophila mitochondrial genomes)). Tetracodon numbers coevolve with expanded tRNA sidearm loops. Sidearm coevolution with amino acid usages and tetragenes occurs for putative anticodons in 5' and 3' sidearms loops (D and TΨC loops, respectively), are stronger for the D-loop. Results slightly favour isolated stem-loops upon crossover hypotheses. An alternative hypothesis, that patterns observed for sidearm 'anticodons' do not imply translational activity, but recognition signals for tRNA synthetases that aminoacylate tRNAs, is incompatible with tetracodon/tetra-anticodon coevolution. Hence analyses strengthen translational hypotheses for tRNA sidearm anticodons, tetragenes, and antisense tRNAs.

Complete mitogenome sequences of four flatfishes (Pleuronectiformes) reveal a novel gene arrangement of L-strand coding genes

Background

Few mitochondrial gene rearrangements are found in vertebrates and large-scale changes in these genomes occur even less frequently. It is difficult, therefore, to propose a mechanism to account for observed changes in mitogenome structure. Mitochondrial gene rearrangements are usually explained by the recombination model or tandem duplication and random loss model.

Results

In this study, the complete mitochondrial genomes of four flatfishes, Crossorhombus azureus (blue flounder), Grammatobothus krempfi, Pleuronichthys cornutus, and Platichthys stellatus were determined. A striking finding is that eight genes in the C. azureus mitogenome are located in a novel position, differing from that of available vertebrate mitogenomes. Specifically, the ND6 and seven tRNA genes (the Q, A, C, Y, S₁, E, P genes) encoded by the L-strand have been translocated to a position between tRNA-T and tRNA-F though the original order of the genes is maintained.

Conclusions

These special features are used to suggest a mechanism for C. azureus mitogenome rearrangement. First, a dimeric molecule was formed by two monomers linked head-to-tail, then one of the two sets of promoters lost function and the genes controlled by the disabled promoters became pseudogenes, non-coding sequences, and even were lost from the genome. This study provides a new gene-rearrangement model that accounts for the events of gene-rearrangement in a vertebrate mitogenome.

The complete mitochondrial genome of Solemya velum (Mollusca: Bivalvia) and its relationships with conchifera

Background

Bivalve mitochondrial genomes exhibit a wide array of uncommon features, like extensive gene rearrangements, large sizes, and unusual ways of inheritance. Species pertaining to the order Solemyida (subclass Opponobranchia) show many peculiar evolutionary adaptations, f.i. extensive symbiosis with chemoautotrophic bacteria. Despite Opponobranchia are central in bivalve phylogeny, being considered the sister group of all Autobranchia, a complete mitochondrial genome has not been sequenced yet.

Results

In this paper, we characterized the complete mitochondrial genome of the Atlantic awning clam Solemya velum: A-T content, gene arrangement and other features are more similar to putative ancestral mollusks than to other bivalves. Two supranumerary open reading frames are present in a large, otherwise unassigned, region, while the origin of replication could be located in a region upstream to the cox3 gene.

Conclusions

We show that S. velum mitogenome retains most of the ancestral conchiferan features, which is unusual among bivalve mollusks, and we discuss main peculiarities of this first example of an organellar genome coming from the subclass Opponobranchia. Mitochondrial genomes of Solemya (for bivalves) and Haliotis (for gastropods) seem to retain the original condition of mollusks, as most probably exemplified by Katharina.

Polymerization of non-complementary RNA: systematic symmetric nucleotide exchanges mainly involving uracil produce mitochondrial RNA transcripts coding for cryptic overlapping genes

Usual DNA→RNA transcription exchanges T→U. Assuming different systematic symmetric nucleotide exchanges during translation, some GenBank RNAs match exactly human mitochondrial sequences (exchange rules listed in decreasing transcript frequencies): C↔U, A↔U, A↔U+C↔G (two nucleotide pairs exchanged), G↔U, A↔G, C↔G, none for A↔C, A↔G+C↔U, and A↔C+G↔U. Most unusual transcripts involve exchanging uracil. Independent measures of rates of rare replicational enzymatic DNA nucleotide misinsertions predict frequencies of RNA transcripts systematically exchanging the corresponding misinserted nucleotides. Exchange transcripts self-hybridize less than other gene regions, self-hybridization increases with length, suggesting endoribonuclease-limited elongation. Blast detects stop codon depleted putative protein coding overlapping genes within exchange-transcribed mitochondrial genes. These align with existing GenBank proteins (mainly metazoan origins, prokaryotic and viral origins underrepresented). These GenBank proteins frequently interact with RNA/DNA, are membrane transporters, or are typical of mitochondrial metabolism. Nucleotide exchange transcript frequencies increase with overlapping gene densities and stop densities, indicating finely tuned counterbalancing regulation of expression of systematic symmetric nucleotide exchange-encrypted proteins. Such expression necessitates combined activities of suppressor tRNAs matching stops, and nucleotide exchange transcription. Two independent properties confirm predicted exchanged overlap coding genes: discrepancy of third codon nucleotide contents from replicational deamination gradients, and codon usage according to circular code predictions. Predictions from both properties converge, especially for frequent nucleotide exchange types. Nucleotide exchanging transcription apparently increases coding densities of protein coding genes without lengthening genomes, revealing unsuspected functional DNA coding potential.

Systematic asymmetric nucleotide exchanges produce human mitochondrial RNAs cryptically encoding for overlapping protein coding genes

GenBank's EST database includes RNAs matching exactly human mitochondrial sequences assuming systematic asymmetric nucleotide exchange-transcription along exchange rules: A→G→C→U/T→A (12 ESTs), A→U/T→C→G→A (4 ESTs), C→G→U/T→C (3 ESTs), and A→C→G→U/T→A (1 EST), no RNAs correspond to other potential asymmetric exchange rules. Hypothetical polypeptides translated from nucleotide-exchanged human mitochondrial protein coding genes align with numerous GenBank proteins, predicted secondary structures resemble their putative GenBank homologue's. Two independent methods designed to detect overlapping genes (one based on nucleotide contents analyses in relation to replicative deamination gradients at third codon positions, and circular code analyses of codon contents based on frame redundancy), confirm nucleotide-exchange-encrypted overlapping genes. Methods converge on which genes are most probably active, and which not, and this for the various exchange rules. Mean EST lengths produced by different nucleotide exchanges are proportional to (a) extents that various bioinformatics analyses confirm the protein coding status of putative overlapping genes; (b) known kinetic chemistry parameters of the corresponding nucleotide substitutions by the human mitochondrial DNA polymerase gamma (nucleotide DNA misinsertion rates); (c) stop codon densities in predicted overlapping genes (stop codon readthrough and exchanging polymerization regulate gene expression by counterbalancing each other). Numerous rarely expressed proteins seem encoded within regular mitochondrial genes through asymmetric nucleotide exchange, avoiding lengthening genomes. Intersecting evidence between several independent approaches confirms the working hypothesis status of gene encryption by systematic nucleotide exchanges.

Coding constraints modulate chemically spontaneous mutational replication gradients in mitochondrial genomes

Distances from heavy and light strand replication origins determine duration mitochondrial DNA remains singlestranded during replication. Hydrolytic deaminations from A->G and C->T occur more on single- than doublestranded DNA. Corresponding replicational nucleotide gradients exist across mitochondrial genomes, most at 3rd, least 2^nd codon positions. DNA singlestrandedness during RNA transcription causes gradients mainly in long-lived species with relatively slow metabolism (high transcription/replication ratios). Third codon nucleotide contents, evolutionary results of mutation cumulation, follow replicational, not transcriptional gradients in Homo; observed human mutations follow transcriptional gradients. Synonymous third codon position transitions potentially alter adaptive off frame information. No mutational gradients occur at synonymous positions forming off frame stops (these adaptively stop early accidental frameshifted protein synthesis), nor in regions coding for putative overlapping genes according to an overlapping genetic code reassigning stop codons to amino acids. Deviation of 3rd codon nucleotide contents from deamination gradients increases with coding importance of main frame 3rd codon positions in overlapping genes (greatest if these are 2^nd position in overlapping genes). Third codon position deamination gradients calculated separately for each codon family are strongest where synonymous transitions are rarely pathogenic; weakest where transitions are frequently pathogenic. Synonymous mutations affect translational accuracy, such as error compensation of misloaded tRNAs by codon-anticodon mismatches (prevents amino acid misinsertion despite tRNA misacylation), a potential cause of pathogenic mutations at synonymous codon positions. Indeed, codon-family-specific gradients are inversely proportional to error compensation associated with gradient-promoted transitions. Deamination gradients reflect spontaneous chemical reactions in singlestranded DNA, but functional coding constraints modulate gradients.