Armless mitochondrial tRNAs in Enoplea (Nematoda)

The diverse structural modes of tRNA binding and recognition

tRNAs are short noncoding RNAs responsible for decoding mRNA codon triplets, delivering correct amino acids to the ribosome, and mediating polypeptide chain formation. Due to their key roles during translation, tRNAs have a highly conserved shape and large sets of tRNAs are present in all living organisms. Regardless of sequence variability, all tRNAs fold into a relatively rigid three-dimensional L-shaped structure. The conserved tertiary organization of canonical tRNA arises through the formation of two orthogonal helices, consisting of the acceptor and anticodon domains. Both elements fold independently to stabilize the overall structure of tRNAs through intramolecular interactions between the D- and T-arm. During tRNA maturation, different modifying enzymes posttranscriptionally attach chemical groups to specific nucleotides, which not only affect translation elongation rates but also restrict local folding processes and confer local flexibility when required. The characteristic structural features of tRNAs are also employed by various maturation factors and modification enzymes to assure the selection, recognition, and positioning of specific sites within the substrate tRNAs. The cellular functional repertoire of tRNAs continues to extend well beyond their role in translation, partly, due to the expanding pool of tRNA-derived fragments. Here, we aim to summarize the most recent developments in the field to understand how three-dimensional structure affects the canonical and noncanonical functions of tRNA.

Rediscovery of Bipalium admarginatum de Beauchamp, 1933 (Platyhelminthes, Tricladida, Geoplanidae) in Malaysia, with molecular characterisation including the mitogenome

We present here the first observation of Bipalium admarginatum de Beauchamp, 1933 since its original description 90 years ago. Three specimens were found on Perhentian Kecil Island, off Terengganu State, Malaysia and photographed in the field, and two were collected. This report thus includes the first colour photographs published for this species, from a locality close to the type-locality, Tioman Island (which is ca. 200 km south of the locality in this study, on the east coast of Peninsula Malaysia). We describe the external morphology and colour pattern of the species, which correspond well to the original description, itself based only on two preserved specimens. We performed an in-depth molecular characterisation of the species, including its complete mitochondrial genome, the 18S sequence and elongation 1-alpha (EF1-α) sequence. In addition, EF1-α sequences were also retrieved for 5 additional geoplanid species. No tRNA-Thr could be detected in the mitogenome of B. admarginatum, a lack already reported in several species of geoplanids, but we found a 13 bp sequence that contains the anticodon loop and seems to be conserved among geoplanids and might thus possibly represent a non-canonical undetected tRNA. We discuss the difficulties encountered in trying to reconstruct the cluster of nuclear ribosomal genes, a problem already mentioned for other Triclads. Three phylogenies, based respectively on all mitochondrial proteins, 18S, and EF1-α, were computed; the position of B. admarginatum within the Bipaliinae was confirmed in each tree, as sister-group to various bipaliine species according to the sequences available for each tree. In the mitochondrial proteins tree, which had high support, B. admarginatum was sister to Bipalium kewense and Diversibipalium multilineatum.

Mitochondrial Genomes Assembled from Non-Invasive eDNA Metagenomic Scat Samples in Critically Endangered Mammals

The abundance of many large-bodied vertebrates, both in marine and terrestrial environments, has declined substantially due to global and regional climate stressors that define the Anthropocene. The development of genetic tools that can serve to monitor population’s health non-intrusively and inform strategies for the recovery of these species is crucial. In this study, we formally evaluate whether whole mitochondrial genomes can be assembled from environmental DNA (eDNA) metagenomics scat samples. Mitogenomes of four different large vertebrates, the panda bear (Ailuropoda melanoleuca), the moon bear (Ursus thibetanus), the Java pangolin (Manis javanica), and the the North Atlantic right whale (Eubalaena glacialis) were assembled and circularized using the pipeline GetOrganelle with a coverage ranging from 12x to 480x in 14 out of 18 different eDNA samples. Partial mitochondrial genomes were retrieved from three other eDNA samples. The complete mitochondrial genomes of the studied species were AT-rich and comprised 13 protein coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a putative D-loop/control region. Synteny observed in all assembled mitogenomes was identical to that reported for specimens of the same and other closely related species. This study demonstrates that it is possible to assemble accurate whole mitochondrial chromosomes from eDNA samples (scats) using forthright bench and bioinformatics workflows. The retrieval of mitochondrial genomes from eDNA samples represents a tool to support bioprospecting, bio-monitoring, and other non-intrusive conservation strategies in species considered ‘vulnerable’, ‘endangered’, and/or ‘critically endangered’ by the IUCN Red List of Threatened Species.

tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes

tRNAscan-SE has been widely used for transfer RNA (tRNA) gene prediction for over twenty years, developed just as the first genomes were decoded. With the massive increase in quantity and phylogenetic diversity of genomes, the accurate detection and functional prediction of tRNAs has become more challenging. Utilizing a vastly larger training set, we created nearly one hundred specialized isotype- and clade-specific models, greatly improving tRNAscan-SE’s ability to identify and classify both typical and atypical tRNAs. We employ a new comparative multi-model strategy where predicted tRNAs are scored against a full set of isotype-specific covariance models, allowing functional prediction based on both the anticodon and the highest-scoring isotype model. Comparative model scoring has also enhanced the program's ability to detect tRNA-derived SINEs and other likely pseudogenes. For the first time, tRNAscan-SE also includes fast and highly accurate detection of mitochondrial tRNAs using newly developed models. Overall, tRNA detection sensitivity and specificity is improved for all isotypes, particularly those utilizing specialized models for selenocysteine and the three subtypes of tRNA genes encoding a CAU anticodon. These enhancements will provide researchers with more accurate and detailed tRNA annotation for a wider variety of tRNAs, and may direct attention to tRNAs with novel traits.

Adaptation of the Romanomermis culicivorax CCA-Adding Enzyme to Miniaturized Armless tRNA Substrates

The mitochondrial genome of the nematode Romanomermis culicivorax encodes for miniaturized hairpin-like tRNA molecules that lack D- as well as T-arms, strongly deviating from the consensus cloverleaf. The single tRNA nucleotidyltransferase of this organism is fully active on armless tRNAs, while the human counterpart is not able to add a complete CCA-end. Transplanting single regions of the Romanomermis enzyme into the human counterpart, we identified a beta-turn element of the catalytic core that-when inserted into the human enzyme-confers full CCA-adding activity on armless tRNAs. This region, originally identified to position the 3'-end of the tRNA primer in the catalytic core, dramatically increases the enzyme's substrate affinity. While conventional tRNA substrates bind to the enzyme by interactions with the T-arm, this is not possible in the case of armless tRNAs, and the strong contribution of the beta-turn compensates for an otherwise too weak interaction required for the addition of a complete CCA-terminus. This compensation demonstrates the remarkable evolutionary plasticity of the catalytic core elements of this enzyme to adapt to unconventional tRNA substrates.

Naturally Occurring tRNAs With Non-canonical Structures

Transfer RNA (tRNA) is the central molecule in genetically encoded protein synthesis. Most tRNA species were found to be very similar in structure: the well-known cloverleaf secondary structure and L-shaped tertiary structure. Furthermore, the length of the acceptor arm, T-arm, and anticodon arm were found to be closely conserved. Later research discovered naturally occurring, active tRNAs that did not fit the established 'canonical' tRNA structure. This review discusses the non-canonical structures of some well-characterized natural tRNA species and describes how these structures relate to their role in translation. Additionally, we highlight some newly discovered tRNAs in which the structure-function relationship is not yet fully understood.

Localization and RNA Binding of Mitochondrial Aminoacyl tRNA Synthetases

Mitochondria contain a complete translation machinery that is used to translate its internally transcribed mRNAs. This machinery uses a distinct set of tRNAs that are charged with cognate amino acids inside the organelle. Interestingly, charging is executed by aminoacyl tRNA synthetases (aaRS) that are encoded by the nuclear genome, translated in the cytosol, and need to be imported into the mitochondria. Here, we review import mechanisms of these enzymes with emphasis on those that are localized to both mitochondria and cytosol. Furthermore, we describe RNA recognition features of these enzymes and their interaction with tRNA and non-tRNA molecules. The dual localization of mitochondria-destined aaRSs and their association with various RNA types impose diverse impacts on cellular physiology. Yet, the breadth and significance of these functions are not fully resolved. We highlight here possibilities for future explorations.

Transfer RNA: The molecular demiurge in the origin of biological systems

Herein, we review recent works on the role that the tRNA molecule played in the early origins of biological systems. tRNAs gave origin to the first genes (mRNA), the peptidyl transferase center (PTC), the 16S ribosomal molecule, proto-tRNAs were at the core of a proto-translation system, and the anticodon and operational codes appeared in tRNAs molecules. Metabolic pathways emerged from evolutionary pressures of the decoding systems. The transitions from the RNA world to the ribonucleoprotein world to modern biological systems were driven by two kinds of tRNAs transitions, to wit, tRNAs leading to both mRNA and rRNA.

Accretion history of large ribosomal subunits deduced from theoretical minimal RNA rings is congruent with histories derived from phylogenetic and structural methods

Accretions of tRNAs presumably formed the large complex ribosomal RNA structures. Similarities of tRNA secondary structures with rRNA secondary structures increase with the integration order of their cognate amino acid in the genetic code, indicating tRNA evolution towards rRNA-like structures. Here analyses rank secondary structure subelements of three large ribosomal RNAs (Prokaryota: Archaea: Thermus thermophilus; Bacteria: Escherichia coli; Eukaryota: Saccharomyces cerevisiae) in relation to their similarities with secondary structures formed by presumed proto-tRNAs, represented by 25 theoretical minimal RNA rings. These ranks are compared to those derived from two independent methods (ranks provide a relative evolutionary age to the rRNA substructure), (a) cladistic phylogenetic analyses and (b) 3D-crystallography where core subelements are presumed ancient and peripheral ones recent. Comparisons of rRNA secondary structure subelements with RNA ring secondary structures show congruence between ranks deduced by this method and both (a) and (b) (more with (a) than (b)), especially for RNA rings with predicted ancient cognate amino acid. Reconstruction of accretion histories of large rRNAs will gain from adequately integrating information from independent methods. Theoretical minimal RNA rings, sequences deterministically designed in silico according to specific coding constraints, might produce adequate scales for prebiotic and early life molecular evolution.

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.

Arm-less mitochondrial tRNAs conserved for over 30 millions of years in spiders

BACKGROUND

In recent years, Next Generation Sequencing (NGS) has accelerated the generation of full mitogenomes, providing abundant material for studying different aspects of molecular evolution. Some mitogenomes have been observed to harbor atypical sequences with bizarre secondary structures, which origins and significance could only be fully understood in an evolutionary framework.

RESULTS

Here we report and analyze the mitochondrial sequences and gene arrangements of six closely related spiders in the sister genera Parachtes and Harpactocrates, which belong to the nocturnal, ground dwelling family Dysderidae. Species of both genera have compacted mitogenomes with many overlapping genes and strikingly reduced tRNAs that are among the shortest described within metazoans. Thanks to the conservation of the gene order and the nucleotide identity across close relatives, we were able to predict the secondary structures even on arm-less tRNAs, which would be otherwise unattainable for a single species. They exhibit aberrant secondary structures with the lack of either DHU or TΨC arms and many miss-pairings in the acceptor arm but this degeneracy trend goes even further since at least four tRNAs are arm-less in the six spider species studied.

CONCLUSIONS

The conservation of at least four arm-less tRNA genes in two sister spider genera for about 30 myr suggest that these genes are still encoding fully functional tRNAs though they may be post-transcriptionally edited to be fully functional as previously described in other species. We suggest that the presence of overlapping and truncated tRNA genes may be related and explains why spider mitogenomes are smaller than those of other invertebrates.

Small but large enough: structural properties of armless mitochondrial tRNAs from the nematode Romanomermis culicivorax

As adapter molecules to convert the nucleic acid information into the amino acid sequence, tRNAs play a central role in protein synthesis. To fulfill this function in a reliable way, tRNAs exhibit highly conserved structural features common in all organisms and in all cellular compartments active in translation. However, in mitochondria of metazoans, certain dramatic deviations from the consensus tRNA structure are described, where some tRNAs lack the D- or T-arm without losing their function. In Enoplea, this miniaturization comes to an extreme, and functional mitochondrial tRNAs can lack both arms, leading to a considerable size reduction. Here, we investigate the secondary and tertiary structure of two such armless tRNAs from Romanomermis culicivorax. Despite their high AU content, the transcripts fold into a single and surprisingly stable hairpin structure, deviating from standard tRNAs. The three-dimensional form is boomerang-like and diverges from the standard L-shape. These results indicate that such unconventional miniaturized tRNAs can still fold into a tRNA-like shape, although their length and secondary structure are very unusual. They highlight the remarkable flexibility of the protein synthesis apparatus and suggest that the translational machinery of Enoplea mitochondria may show compensatory adaptations to accommodate these armless tRNAs for efficient translation.

The mammalian mitochondrial epitranscriptome

Correct expression of the mitochondrially-encoded genes is critical for the production of the components of the oxidative phosphorylation machinery. Post-transcriptional modifications of mitochondrial transcripts have been emerging as an important regulatory feature of mitochondrial gene expression. Here we review the current knowledge on how the mammalian mitochondrial epitranscriptome participates in regulating mitochondrial homeostasis. In particular, we focus on the latest breakthroughs made towards understanding the roles of the modified nucleotides in mitochondrially-encoded ribosomal and transfer RNAs, the enzymes responsible for introducing these modifications and on recent transcriptome-wide studies reporting modifications to mitochondrial messenger RNAs. This article is part of a Special Issue entitled: mRNA modifications in gene expression control edited by Dr. Matthias Soller and Dr. Rupert Fray.

Post-Transcriptional Regulation of tRNA Pools To Govern the Central Dogma: A Perspective

Since their initial discovery, tRNAs have risen from sole adapter molecules during protein synthesis to pivotal modulators of gene expression. Through their many interactions with tRNA-associated protein factors, they play a central role in maintaining cell homeostasis, especially regarding the fine-tuning in response to a rapidly changing cellular environment. Here, we provide a perspective on current tRNA topics with a spotlight on the regulation of post-transcriptional shaping of tRNA molecules. First, we give an update on aberrant structural features that a yet functional fraction of mitochondrial tRNAs can exhibit. Then, we outline several aspects of the regulatory contribution of ribonucleases with a focus on tRNA processing versus tRNA elimination. We close with a comment on the possible consequences for the intracellular examination of nascent tRNA precursors regarding respective processing factors that have been shown to associate with the tRNA transcription machinery in alternative moonlighting functions.

The mitochondrial genome of the oribatid mite Paraleius leontonychus: new insights into tRNA evolution and phylogenetic relationships in acariform mites

Bilaterian mitochondrial (mt) genomes are circular molecules that typically contain 37 genes. To date, only a single complete mitogenome sequence is available for the species-rich sarcoptiform mite order Oribatida. We sequenced the mitogenome of Paraleius leontonychus, another species of this suborder. It is 14,186 bp long and contains 35 genes, including only 20 tRNAs, lacking tRNA^Gly and tRNA^Tyr. Re-annotation of the mitogenome of Steganacarus magnus increased the number of mt tRNAs for this species to 12. As typical for acariform mites, many tRNAs are highly truncated in both oribatid species. The total number of tRNAs and the number of tRNAs with a complete cloverleaf-like structure in P. leontonychus, however, clearly exceeds the numbers previously reported for Sarcoptiformes. This indicates, contrary to what has been previously assumed, that reduction of tRNAs is not a general characteristic for sarcoptiform mites. Compared to other Sarcoptiformes, the two oribatid species have the least rearranged mt genome with respect to the pattern observed in Limulus polyphemus, a basal arachnid species. Phylogenetic analysis of the newly sequenced mt genome and previously published data on other acariform mites confirms paraphyly of the Oribatida and an origin of the Astigmata within the Oribatida.

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.

tRNA Modifications: Impact on Structure and Thermal Adaptation

Transfer RNAs (tRNAs) are central players in translation, functioning as adapter molecules between the informational level of nucleic acids and the functional level of proteins. They show a highly conserved secondary and tertiary structure and the highest density of post-transcriptional modifications among all RNAs. These modifications concentrate in two hotspots-the anticodon loop and the tRNA core region, where the D- and T-loop interact with each other, stabilizing the overall structure of the molecule. These modifications can cause large rearrangements as well as local fine-tuning in the 3D structure of a tRNA. The highly conserved tRNA shape is crucial for the interaction with a variety of proteins and other RNA molecules, but also needs a certain flexibility for a correct interplay. In this context, it was shown that tRNA modifications are important for temperature adaptation in thermophilic as well as psychrophilic organisms, as they modulate rigidity and flexibility of the transcripts, respectively. Here, we give an overview on the impact of modifications on tRNA structure and their importance in thermal adaptation.

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.

Mitochondrial genome diversity in dagger and needle nematodes (Nematoda: Longidoridae)

Dagger and needle nematodes included in the family Longidoridae (viz. Longidorus, Paralongidorus, and Xiphinema) are highly polyphagous plant-parasitic nematodes in wild and cultivated plants and some of them are plant-virus vectors (nepovirus). The mitochondrial (mt) genomes of the dagger and needle nematodes, Xiphinema rivesi, Xiphinema pachtaicum, Longidorus vineacola and Paralongidorus litoralis were sequenced in this study. The four circular mt genomes have an estimated size of 12.6, 12.5, 13.5 and 12.7 kb, respectively. Up to date, the mt genome of X. pachtaicum is the smallest genome found in Nematoda. The four mt genomes contain 12 protein-coding genes (viz. cox1-3, nad1-6, nad4L, atp6 and cob) and two ribosomal RNA genes (rrnL and rrnS), but the atp8 gene was not detected. These mt genomes showed a gene arrangement very different within the Longidoridae species sequenced, with the exception of very closely related species (X. americanum and X. rivesi). The sizes of non-coding regions in the Longidoridae nematodes were very small and were present in a few places in the mt genome. Phylogenetic analysis of all coding genes showed a closer relationship between Longidorus and Paralongidorus and different phylogenetic possibilities for the three Xiphinema species.

Mitogenomics reveals high synteny and long evolutionary histories of sympatric cryptic nematode species

Species with seemingly identical morphology but with distinct genetic differences are abundant in the marine environment and frequently co-occur in the same habitat. Such cryptic species are typically delineated using a limited number of mitochondrial and/or nuclear marker genes, which do not yield information on gene order and gene content of the genomes under consideration. We used next-generation sequencing to study the composition of the mitochondrial genomes of four sympatrically distributed cryptic species of the Litoditis marina species complex (PmI, PmII, PmIII, and PmIV). The ecology, biology, and natural occurrence of these four species are well known, but the evolutionary processes behind this cryptic speciation remain largely unknown. The gene order of the mitochondrial genomes of the four species was conserved, but differences in genome length, gene length, and codon usage were observed. The atp8 gene was lacking in all four species. Phylogenetic analyses confirm that PmI and PmIV are sister species and that PmIII diverged earliest. The most recent common ancestor of the four cryptic species was estimated to have diverged 16 MYA. Synonymous mutations outnumbered nonsynonymous changes in all protein-encoding genes, with the Complex IV genes (coxI-III) experiencing the strongest purifying selection. Our mitogenomic results show that morphologically similar species can have long evolutionary histories and that PmIII has several differences in genetic makeup compared to the three other species, which may explain why it is better adapted to higher temperatures than the other species.

First complete mitochondrial genome of the South American annual fish Austrolebias charrua (Cyprinodontiformes: Rivulidae): peculiar features among cyprinodontiforms mitogenomes

Background

Among teleosts, the South American genus Austrolebias (Cyprinodontiformes: Rivulidae) includes 42 taxa of annual fishes divided into five different species groups. It is a monophyletic genus, but morphological and molecular data do not resolve the relationship among intrageneric clades and high rates of substitution have been previously described in some mitochondrial genes. In this work, the complete mitogenome of a species of the genus was determined for the first time. We determined its structure, gene order and evolutionary peculiar features, which will allow us to evaluate the performance of mitochondrial genes in the phylogenetic resolution at different taxonomic levels.

Results

Regarding gene content and order, the circular mitogenome of A. charrua (17,271 pb) presents the typical pattern of vertebrate mitogenomes. It contains the full complement of 13 proteins-coding genes, 22 tRNA, 2 rRNA and one non-coding control region. Notably, the tRNA-Cys was only 57 bp in length and lacks the D-loop arm. In three full sibling individuals, heteroplasmatic condition was detected due to a total of 12 variable sites in seven protein-coding genes. Among cyprinodontiforms, the mitogenome of A. charrua exhibits the lowest G+C content (37 %) and GCskew, as well as the highest strand asymmetry with a net difference of T over A at 1st and 3rd codon positions. Considering the 12 coding-genes of the H strand, correspondence analyses of nucleotide composition and codon usage show that A and T at 1st and 3rd codon positions have the highest weight in the first axis, and segregate annual species from the other cyprinodontiforms analyzed. Given the annual life-style, their mitogenomes could be under different selective pressures. All 13 protein-coding genes are under strong purifying selection and we did not find any significant evidence of nucleotide sites showing episodic selection (dN >dS) at annual lineages. When fast evolving third codon positions were removed from alignments, the “supergene” tree recovers our reference species phylogeny as well as the Cytb, ND4L and ND6 genes. Therefore, third codon positions seem to be saturated in the aforementioned coding regions at intergeneric Cyprinodontiformes comparisons.

Conclusions

The complete mitogenome obtained in present work, offers relevant data for further comparative studies on molecular phylogeny and systematics of this taxonomic controversial endemic genus of annual fishes.

Electronic supplementary material

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

Large gene overlaps and tRNA processing in the compact mitochondrial genome of the crustacean Armadillidium vulgare

A faithful expression of the mitochondrial DNA is crucial for cell survival. Animal mitochondrial DNA (mtDNA) presents a highly compact gene organization. The typical 16.5 kbp animal mtDNA encodes 13 proteins, 2 rRNAs and 22 tRNAs. In the backyard pillbug Armadillidium vulgare, the rather small 13.9 kbp mtDNA encodes the same set of proteins and rRNAs as compared to animal kingdom mtDNA, but seems to harbor an incomplete set of tRNA genes. Here, we first confirm the expression of 13 tRNA genes in this mtDNA. Then we show the extensive repair of a truncated tRNA, the expression of tRNA involved in large gene overlaps and of tRNA genes partially or fully integrated within protein-coding genes in either direct or opposite orientation. Under selective pressure, overlaps between genes have been likely favored for strong genome size reduction. Our study underlines the existence of unknown biochemical mechanisms for the complete gene expression of A. vulgare mtDNA, and of co-evolutionary processes to keep overlapping genes functional in a compacted mitochondrial genome.

Towards a comprehensive picture of alloacceptor tRNA remolding in metazoan mitochondrial genomes

Remolding of tRNAs is a well-documented process in mitochondrial genomes that changes the identity of a tRNA. It involves a duplication of a tRNA gene, a mutation that changes the anticodon and the loss of the ancestral tRNA gene. The net effect is a functional tRNA that is more closely related to tRNAs of a different alloacceptor family than to tRNAs with the same anticodon in related species. Beyond being of interest for understanding mitochondrial tRNA function and evolution, tRNA remolding events can lead to artifacts in the annotation of mitogenomes and thus in studies of mitogenomic evolution. Therefore, it is important to identify and catalog these events. Here we describe novel methods to detect tRNA remolding in large-scale data sets and apply them to survey tRNA remolding throughout animal evolution. We identify several novel remolding events in addition to the ones previously mentioned in the literature. A detailed analysis of these remoldings showed that many of them are derived from ancestral events.

Precise mapping and dynamics of tRNA-derived fragments (tRFs) in the development of Triops cancriformis (tadpole shrimp)

Background

In a deep sequencing analysis of small RNAs prepared from a living fossil, the tadpole shrimp Triops cancriformis, a 32-nt small RNA was specifically detected in the adult stage. A nucleotide sequence comparison between the 32-nt small RNA and predicted tRNA sequences in the draft nuclear genomic DNA showed that the small RNA was derived from tRNA^Gly(GCC). To determine the overall features of the tRNA-derived fragments (tRFs) of T. cancriformis, the small RNA sequences in each of the six developmental stages (egg, 1st − 4th instar larvae, and adult) were compared with the mitochondrial and nuclear tRNA sequences.

Results

We found that the tRFs were derived from mitochondrial and nuclear tRNAs corresponding to 16 and 39 anticodons, respectively. The total read number of nuclear tRFs was approximately 400 times larger than the number of mitochondrial tRFs. Interestingly, the main regions in each parental tRNA from which these tRFs were derived differed, depending on the parental anticodon. Mitochondrial tRF^Ser(GCU)s were abundantly produced from the 5’ half regions of the parental tRNA, whereas mitochondrial tRF^Val(UAC)s were mainly produced from the 3’ end regions. Highly abundant nuclear tRFs, tRF^Gly(GCC)s, tRF^Gly(CCC)s, tRF^Glu(CUC)s, and tRF^Lys(CUU)s were derived from the 5’ half regions of the parental tRNAs. Further analysis of the tRF read counts in the individual developmental stages suggested that the expression of mitochondrial and nuclear tRFs differed during the six stages. Based on these data, we precisely summarized the positions of the tRFs in their parental tRNAs and their expression changes during development.

Conclusions

Our results reveal the entire dynamics of the tRFs from both the nuclear and mitochondrial genomes of T. cancriformis and indicate that the majority of tRFs in the cell are derived from nuclear tRNAs. This study provides the first examples of developmentally expressed mitochondrial tRFs.

Electronic supplementary material

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

The ancestor of modern Holozoa acquired the CCA-adding enzyme from Alphaproteobacteria by horizontal gene transfer

Transfer RNAs (tRNAs) require the absolutely conserved sequence motif CCA at their 3′-ends, representing the site of aminoacylation. In the majority of organisms, this trinucleotide sequence is not encoded in the genome and thus has to be added post-transcriptionally by the CCA-adding enzyme, a specialized nucleotidyltransferase. In eukaryotic genomes this ubiquitous and highly conserved enzyme family is usually represented by a single gene copy. Analysis of published sequence data allows us to pin down the unusual evolution of eukaryotic CCA-adding enzymes. We show that the CCA-adding enzymes of animals originated from a horizontal gene transfer event in the stem lineage of Holozoa, i.e. Metazoa (animals) and their unicellular relatives, the Choanozoa. The tRNA nucleotidyltransferase, acquired from an α-proteobacterium, replaced the ancestral enzyme in Metazoa. However, in Choanoflagellata, the group of Choanozoa that is closest to Metazoa, both the ancestral and the horizontally transferred CCA-adding enzymes have survived. Furthermore, our data refute a mitochondrial origin of the animal tRNA nucleotidyltransferases.

The identity of the discriminator base has an impact on CCA addition

CCA-adding enzymes synthesize and maintain the C-C-A sequence at the tRNA 3'-end, generating the attachment site for amino acids. While tRNAs are the most prominent substrates for this polymerase, CCA additions on non-tRNA transcripts are described as well. To identify general features for substrate requirement, a pool of randomized transcripts was incubated with the human CCA-adding enzyme. Most of the RNAs accepted for CCA addition carry an acceptor stem-like terminal structure, consistent with tRNA as the main substrate group for this enzyme. While these RNAs show no sequence conservation, the position upstream of the CCA end was in most cases represented by an adenosine residue. In tRNA, this position is described as discriminator base, an important identity element for correct aminoacylation. Mutational analysis of the impact of the discriminator identity on CCA addition revealed that purine bases (with a preference for adenosine) are strongly favoured over pyrimidines. Furthermore, depending on the tRNA context, a cytosine discriminator can cause a dramatic number of misincorporations during CCA addition. The data correlate with a high frequency of adenosine residues at the discriminator position observed in vivo. Originally identified as a prominent identity element for aminoacylation, this position represents a likewise important element for efficient and accurate CCA addition.

The complete mitochondrial genome of the facultative entomopathogenic nematode Oscheius chongmingensis (Rhabditida: Rhabditidae)

We determined the complete mitochondrial genome of the facultative entomopathogenic nematode Oscheius chongmingensis. The mitogenome length was 15,413 bp and similar to other Rhabditids contains genes for 2 rRNAs, 22 tRNAs, and 12 proteins (ATPase subunit 8 is missing). Predicted tRNAs indicated the secondary structure typical for chromadorean nematodes. Gene order is similar to that observed in the genus Caenorhabditis. The control AT-rich region is considerably large (2061 bp, 84% of AT), positioned in between tRNA(Ala) and tRNA(Pro) and has several microsatellite-like (AT)n elements.

tRNA biology in mitochondria

Mitochondria are the powerhouses of eukaryotic cells. They are considered as semi-autonomous because they have retained genomes inherited from their prokaryotic ancestor and host fully functional gene expression machineries. These organelles have attracted considerable attention because they combine bacterial-like traits with novel features that evolved in the host cell. Among them, mitochondria use many specific pathways to obtain complete and functional sets of tRNAs as required for translation. In some instances, tRNA genes have been partially or entirely transferred to the nucleus and mitochondria require precise import systems to attain their pool of tRNAs. Still, tRNA genes have also often been maintained in mitochondria. Their genetic arrangement is more diverse than previously envisaged. The expression and maturation of mitochondrial tRNAs often use specific enzymes that evolved during eukaryote history. For instance many mitochondria use a eukaryote-specific RNase P enzyme devoid of RNA. The structure itself of mitochondrial encoded tRNAs is also very diverse, as e.g., in Metazoan, where tRNAs often show non canonical or truncated structures. As a result, the translational machinery in mitochondria evolved adapted strategies to accommodate the peculiarities of these tRNAs, in particular simplified identity rules for their aminoacylation. Here, we review the specific features of tRNA biology in mitochondria from model species representing the major eukaryotic groups, with an emphasis on recent research on tRNA import, maturation and aminoacylation.

The ribosome as a missing link in the evolution of life

Many steps in the evolution of cellular life are still mysterious. We suggest that the ribosome may represent one important missing link between compositional (or metabolism-first), RNA-world (or genes-first) and cellular (last universal common ancestor) approaches to the evolution of cells. We present evidence that the entire set of transfer RNAs for all twenty amino acids are encoded in both the 16S and 23S rRNAs of Escherichia coli K12; that nucleotide sequences that could encode key fragments of ribosomal proteins, polymerases, ligases, synthetases, and phosphatases are to be found in each of the six possible reading frames of the 16S and 23S rRNAs; and that every sequence of bases in rRNA has information encoding more than one of these functions in addition to acting as a structural component of the ribosome. Ribosomal RNA, in short, is not just a structural scaffold for proteins, but the vestigial remnant of a primordial genome that may have encoded a self-organizing, self-replicating, auto-catalytic intermediary between macromolecules and cellular life.

Chemical synthesis of the 5-taurinomethyl(-2-thio)uridine modified anticodon arm of the human mitochondrial tRNA(Leu(UUR)) and tRNA(Lys)

5-Taurinomethyluridine (τm(5)U) and 5-taurinomethyl-2-thiouridine (τm(5)s(2)U) are located at the wobble position of human mitochondrial (hmt) tRNA(Leu(UUR)) and tRNA(Lys), respectively. Both hypermodified units restrict decoding of the third codon letter to A and G. Pathogenic mutations in the genes encoding hmt-tRNA(Leu(UUR)) and hmt-tRNA(Lys) are responsible for the loss of the discussed modifications and, as a consequence, for the occurrence of severe mitochondrial dysfunctions (MELAS, MERRF). Synthetic oligoribonucleotides bearing modified nucleosides are a versatile tool for studying mechanisms of genetic message translation and accompanying pathologies at nucleoside resolution. In this paper, we present site-specific chemical incorporation of τm(5)U and τm(5)s(2)U into 17-mers related to the sequence of the anticodon arms hmt-tRNA(Leu(UUR)) and hmt-tRNA(Lys), respectively employing phosphoramidite chemistry on CPG support. Selected protecting groups for the sulfonic acid (4-(tert-butyldiphenylsilanyloxy)-2,2-dimethylbutyl) and the exoamine function (-C(O)CF3) are compatible with the blockage of the canonical monomeric units. The synthesis of τm(5)s(2)U-modified RNA fragment was performed under conditions eliminating the formation of side products of 2-thiocarbonyl group oxidation and/or oxidative desulphurization. The structure of the final oligomers was confirmed by mass spectroscopy and enzymatic cleavage data.

Losing the stem-loop structure from metazoan mitochondrial tRNAs and co-evolution of interacting factors

Conventional tRNAs have highly conserved sequences, four-armed cloverleaf secondary structures, and L-shaped tertiary structures. However, metazoan mitochondrial tRNAs contain several exceptional structures. Almost all tRNAs^Ser for AGY/N codons lack the D-arm. Furthermore, in some nematodes, no four-armed cloverleaf-type tRNAs are present: two tRNAs^Ser without the D-arm and 20 tRNAs without the T-arm are found. Previously, we showed that in nematode mitochondria, an extra elongation factor Tu (EF-Tu) has evolved to support interaction with tRNAs lacking the T-arm, which interact with C-terminal domain 3 in conventional EF-Tu. Recent mitochondrial genome analyses have suggested that in metazoan lineages other than nematodes, tRNAs without the T-arm are present. Furthermore, even more simplified tRNAs are predicted in some lineages. In this review, we discuss mitochondrial tRNAs with divergent structures, as well as protein factors, including EF-Tu, that support the function of truncated metazoan mitochondrial tRNAs.

Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing

Little is known about the variations of nematode mitogenomes (mtDNA). Sequencing a complete mtDNA using a PCR approach remains a challenge due to frequent genome reorganizations and low sequence similarities between divergent nematode lineages. Here, a genome skimming approach based on HiSeq sequencing (shotgun) was used to assemble de novo the first complete mtDNA sequence of a root-knot nematode (Meloidogyne graminicola). An AT-rich genome (84.3%) of 20,030 bp was obtained with a mean sequencing depth superior to 300. Thirty-six genes were identified with a semi-automated approach. A comparison with a gene map of the M. javanica mitochondrial genome indicates that the gene order is conserved within this nematode lineage. However, deep genome rearrangements were observed when comparing with other species of the superfamily Hoplolaimoidea. Repeat elements of 111 bp and 94 bp were found in a long non-coding region of 7.5 kb, as similarly reported in M. javanica and M. hapla. This study points out the power of next generation sequencing to produce complete mitochondrial genomes, even without a reference sequence, and possibly opening new avenues for species/race identification, phylogenetics and population genetics of nematodes.

Idiosyncrasies in decoding mitochondrial genomes

Mitochondria originate from the α-proteobacterial domain of life. Since this unique event occurred, mitochondrial genomes of protozoans, fungi, plants and metazoans have highly derived and diverged away from the common ancestral DNA. These resulting genomes highly differ from one another, but all present-day mitochondrial DNAs have a very reduced coding capacity. Strikingly however, ATP production coupled to electron transport and translation of mitochondrial proteins are the two common functions retained in all mitochondrial DNAs. Paradoxically, most components essential for these two functions are now expressed from nuclear genes. Understanding how mitochondrial translation evolved in various eukaryotic models is essential to acquire new knowledge of mitochondrial genome expression. In this review, we provide a thorough analysis of the idiosyncrasies of mitochondrial translation as they occur between organisms. We address this by looking at mitochondrial codon usage and tRNA content. Then, we look at the aminoacyl-tRNA-forming enzymes in terms of peculiarities, dual origin, and alternate function(s). Finally we give examples of the atypical structural properties of mitochondrial tRNAs found in some organisms and the resulting adaptive tRNA-protein partnership.

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.

Biological evidence for the world's smallest tRNAs

Due to their function as adapters in translation, tRNA molecules share a common structural organization in all kingdoms and organelles with ribosomal protein biosynthesis. A typical tRNA has a cloverleaf-like secondary structure, consisting of acceptor stem, D-arm, anticodon arm, a variable region, and T-arm, with an average length of 73 nucleotides. In several mitochondrial genomes, however, tRNA genes encode transcripts that show a considerable deviation of this standard, having reduced D- or T-arms or even completely lack one of these elements, resulting in tRNAs as small as 66 nts. An extreme case of such truncations is found in the mitochondria of Enoplea. Here, several tRNA genes are annotated that lack both the D- and the T-arm, suggesting even shorter transcripts with a length of only 42 nts. However, direct evidence for these exceptional tRNAs, which were predicted by purely computational means, has been lacking so far. Here, we demonstrate that several of these miniaturized armless tRNAs consisting only of acceptor- and anticodon-arms are indeed transcribed and correctly processed by non-encoded CCA addition in the mermithid Romanomermis culicivorax. This is the first direct evidence for the existence and functionality of the smallest tRNAs ever identified so far. It opens new possibilities towards exploration/assessment of minimal structural motifs defining a functional tRNA and their evolution.

Genetic aspects of mitochondrial genome evolution

Many years of extensive studies of metazoan mitochondrial genomes have established differences in gene arrangements and genetic codes as valuable phylogenetic markers. Understanding the underlying mechanisms of replication, transcription and the role of the control regions which cause e.g. different gene orders is important to assess the phylogenetic signal of such events. This review summarises and discusses, for the Metazoa, the general aspects of mitochondrial transcription and replication with respect to control regions as well as several proposed models of gene rearrangements. As whole genome sequencing projects accumulate, more and more observations about mitochondrial gene transfer to the nucleus are reported. Thus occurrence and phylogenetic aspects concerning nuclear mitochondrial-like sequences (NUMTS) is another aspect of this review.