Stable tRNA-based phylogenies using only 76 nucleotides

Investigations into the Diversity and Distribution of tRNA and Phylogenetics of Translation Factors in Amoebozoa-Infecting Nucleocytoviricota

Translation is a sine qua non process for life as we know it. Translation factors (TFs) and tRNAs are rare among viruses but are commonly found in giant viruses of the class Megaviricetes. In this study, we explored the diversity and distribution of tRNAs in giant viruses that were isolated and replicated in amoebae (phylum Amoebozoa), and investigated the evolutionary history of TFs to gain insights into their origins in these viruses. We analyzed the genomes of 77 isolated giant viruses, 52 of which contained at least 1 tRNA. In most of these viruses, tRNA sequences are dispersed throughout the genome, except in Tupanviruses and Yasmineviruses, where most tRNAs are clustered in specific genomic islands. The tRNAs in giant viruses often contain introns, with 73.1% of the genomes exhibiting at least one intronic region in these genes. Codon usage bias (CUB) analysis of various giant viruses revealed at least two distinct patterns of codon preferences among closely related viruses. We did not observe a clear correlation between the presence of tRNAs and CUB in giant viruses. Due to the limited size of these genes, we could not confidently investigate their phylogenetic relationships. However, phylogenetic analysis of TFs found in giant viruses often position these viruses as sister groups or embedded between different eukaryotic taxa with high statistical support. Overall, our findings reinforce the complexity of key components of the translation apparatus in different members of Nucleocytoviricota isolated from different regions of Earth.

Enzymic recognition of amino acids drove the evolution of primordial genetic codes

How genetic information gained its exquisite control over chemical processes needed to build living cells remains an enigma. Today, the aminoacyl-tRNA synthetases (AARS) execute the genetic codes in all living systems. But how did the AARS that emerged over three billion years ago as low-specificity, protozymic forms then spawn the full range of highly-specific enzymes that distinguish between 22 diverse amino acids? A phylogenetic reconstruction of extant AARS genes, enhanced by analysing modular acquisitions, reveals six AARS with distinct bacterial, archaeal, eukaryotic, or organellar clades, resulting in a total of 36 families of AARS catalytic domains. Small structural modules that differentiate one AARS family from another played pivotal roles in discriminating between amino acid side chains, thereby expanding the genetic code and refining its precision. The resulting model shows a tendency for less elaborate enzymes, with simpler catalytic domains, to activate amino acids that were not synthesised until later in the evolution of the code. The most probable evolutionary route for an emergent amino acid type to establish a place in the code was by recruiting older, less specific AARS, rather than adapting contemporary lineages. This process, retrofunctionalisation, differs from previously described mechanisms through which amino acids would enter the code.

Identification and analysis of putative tRNA genes in baculovirus genomes

Transfer RNAs (tRNAs) genes are both coded for and arranged along some viral genomes representing the entire virosphere and seem to play different biological functions during infection, other than transferring the correct amino acid to a growing peptide chain. Baculovirus genome description and annotation has focused mostly on protein-coding genes, microRNA, and homologous regions. Here we carried out a large-scale in silico search for putative tRNA genes in baculovirus genomes. Ninety-six of 257 baculovirus genomes analyzed was found to contain at least one putative tRNA gene. We found great diversity in primary and secondary structure, in location within the genome, in intron presence and size, and in anti-codon identity. In some cases, genes of tRNA-containing genomes were found to have a bias for the codons specified by the tRNAs present in such genomes. Moreover, analysis revealed that most of the putative tRNA genes possessed conserved motifs for tRNA type 2 promoters, including the A-box and B-box motifs with few mismatches from the eukaryotic canonical motifs. From publicly available small RNA deep sequencing datasets of baculovirus-infected insect cells, we found evidence that a putative Autographa californica multiple nucleopolyhedrovirus Gln-tRNA gene was transcribed and modified with the addition of the non-templated 3'-CCA tail found at the end of all tRNAs. Further research is needed to determine the expression and functionality of these viral tRNAs.

tRNA modification profiles in obligate and moderate thermophilic bacilli

Transfer RNAs (tRNAs) are the most ancient RNA molecules in the cell, modification pattern of which is linked to phylogeny. The aim of this study was to determine the tRNA modification profiles of obligate (Anoxybacillus, Geobacillus, Paragebacillus) and moderate (Bacillus, Brevibacillus, Ureibacillus, Paenibacillus) thermophilic aerobic bacilli strains to find out its linkage to phylogenetic variations between species. LC-MS was applied for the quantification of modified nucleosides using both natural and isotopically labeled standards. The presence of m²A and m⁷G modifications at high levels was determined in all species. Relatively high level of i⁶A and m⁵C modification was observed for Paenibacillus and Ureibacillus, respectively. The lowest level of Cm modification was found in Bacillus. The modification ms²i⁶A and m¹G were absent in Brevibacillus and Ureibacillus, respectively, while modifications Am and m²₂G were observed only for Ureibacillus. While both obligate and moderate thermophilic species contain Gm, m¹G and ms²i⁶A modifications, large quantities of them (especially Gm and ms²i⁶A modification) were detected in obligate thermophilic ones (Geobacillus, Paragebacillus and Anoxybacillus). The collective set of modified tRNA bases is genus-specific and linked to the phylogeny of bacilli. In addition, the dataset could be applied to distinguish obligate thermophilic bacilli from moderate ones.

Assessing the involvement of the placental microbiome and virome in preeclampsia using non coding RNA sequencing

OBJECTIVES

Preeclampsia is a dangerous pregnancy complication. The source of preeclampsia is unknown, though the placenta is believed to have a central role in its pathogenesis. An association between maternal infection and preeclampsia has been demonstrated, yet the involvement of the placental microbiome in the etiology of preeclampsia has not been determined. In this study, we examined whether preeclampsia is associated with an imbalanced microorganism composition in the placenta.

METHODS

To this end, we developed a novel method for the identification of bacteria/viruses based on sequencing of small non-coding RNA, which increases the microorganism-to-host ratio, this being a major challenge in microbiome methods. We validated the method on various infected tissues and demonstrated its efficiency in detecting microorganisms in samples with extremely low bacterial/viral biomass. We then applied the method to placenta specimens from preeclamptic and healthy pregnancies. Since the placenta is a remarkably large and heterogeneous organ, we explored the bacterial and viral RNA at each of 15 distinct locations.

RESULTS

Bacterial RNA was detected at all locations and was consistent with previous studies of the placental microbiome, though without significant differences between the preeclampsia and control groups. Nevertheless, the bacterial RNA composition differed significantly between various areas of the placenta. Viral RNA was detected in extremely low quantities, below the threshold of significance, thus viral abundance could not be determined.

CONCLUSIONS

Our results suggest that the bacterial and viral abundance in the placenta may have only limited involvement in the pathogenesis of preeclampsia. The evidence of a heterogenic bacterial RNA composition in the various placental locations warrants further investigation to capture the true nature of the placental microbiome.

Evolutionary characteristics and phylogeny of cotton chloroplast tRNAs

The novel structural variations were identified in cotton chloroplast tRNAs and gene loss events were more obvious than duplications in chloroplast tRNAs. Transfer RNAs (tRNA) have long been believed an evolutionary-conserved molecular family, which play the key roles in the process of protein biosynthesis in plant life activities. In this study, we detected the evolutionary characteristics and phylogeny of chloroplast tRNAs in cotton plants, an economic and fibered important taxon in the world. We firstly annotated the chloroplast tRNAs of 27 Gossypium species to analyze their genetic composition, structural characteristics and evolution. Compared with the traditional view of evolutionary conservation of tRNA, some novel tRNA structural variations were identified in cotton plants. I.g., tRNA^Val-UAC and tRNA^Ile-GAU only contained one intron in the anti-condon loop region of tRNA secondary structure, respectively. In the variable region, some tRNAs contained a circle structure with a few nucleotides. Interestingly, the calculation result of free energy indicated that the variation of novel tRNAs contributed to the stability of tRNA structure. Phylogenetic analysis suggested that chloroplast tRNAs have evolved from multiple common ancestors, and the tRNA^Met seemed to be an ancestral tRNA, which can be duplicated and diversified to produce other tRNAs. The chloroplast tRNAs contained a group I intron in cotton plants, and the evolutionary analysis of introns indicated that group I intron of chloroplast tRNA originated from cyanobacteria. Analysis of gene duplication and loss events showed that gene loss events were more obvious than duplications in Gossypium chloroplast tRNAs. Additionally, we found that the rate of transition was higher than the ones of transversion in cotton chloroplast tRNAs. This study provided new insights into the structural characteristics and evolution of chloroplast tRNAs in cotton plants.

Recoding of stop codons expands the metabolic potential of two novel Asgardarchaeota lineages

Asgardarchaeota have been proposed as the closest living relatives to eukaryotes, and a total of 72 metagenome-assembled genomes (MAGs) representing six primary lineages in this archaeal phylum have thus far been described. These organisms are predicted to be fermentative heterotrophs contributing to carbon cycling in sediment ecosystems. Here, we double the genomic catalogue of Asgardarchaeota by obtaining 71 MAGs from a range of habitats around the globe, including the deep subsurface, brackish shallow lakes, and geothermal spring sediments. Phylogenomic inferences followed by taxonomic rank normalisation confirmed previously established Asgardarchaeota classes and revealed four additional lineages, two of which were consistently recovered as monophyletic classes. We therefore propose the names Candidatus Sifarchaeia class nov. and Ca. Jordarchaeia class nov., derived from the gods Sif and Jord in Norse mythology. Metabolic inference suggests that both classes represent hetero-organotrophic acetogens, which also have the ability to utilise methyl groups such as methylated amines, with acetate as the probable end product in remnants of a methanogen-derived core metabolism. This inferred mode of energy conservation is predicted to be enhanced by genetic code expansions, i.e., stop codon recoding, allowing the incorporation of the rare 21st and 22nd amino acids selenocysteine (Sec) and pyrrolysine (Pyl). We found Sec recoding in Jordarchaeia and all other Asgardarchaeota classes, which likely benefit from increased catalytic activities of Sec-containing enzymes. Pyl recoding, on the other hand, is restricted to Sifarchaeia in the Asgardarchaeota, making it the first reported non-methanogenic archaeal lineage with an inferred complete Pyl machinery, likely providing members of this class with an efficient mechanism for methylamine utilisation. Furthermore, we identified enzymes for the biosynthesis of ester-type lipids, characteristic of bacteria and eukaryotes, in both newly described classes, supporting the hypothesis that mixed ether-ester lipids are a shared feature among Asgardarchaeota.

Steller's sea cow genome suggests this species began going extinct before the arrival of Paleolithic humans

Anthropogenic activity is the top factor directly related to the extinction of several animal species. The last Steller's sea cow (Hydrodamalis gigas) population on the Commander Islands (Russia) was wiped out in the second half of the 18^th century due to sailors and fur traders hunting it for the meat and fat. However, new data suggests that the extinction process of this species began much earlier. Here, we present a nuclear de novo assembled genome of H. gigas with a 25.4× depth coverage. Our results demonstrate that the heterozygosity of the last population of this animal is low and comparable to the last woolly mammoth population that inhabited Wrangel Island 4000 years ago. Besides, as a matter of consideration, our findings also demonstrate that the extinction of this marine mammal starts along the North Pacific coastal line much earlier than the first Paleolithic humans arrived in the Bering sea region.

On the Track of the Missing tRNA Genes: A Source of Non-Canonical Functions?

Cellular tRNAs appear today as a diverse population of informative macromolecules with conserved general elements ensuring essential common functions and different and distinctive features securing specific interactions and activities. Their differential expression and the variety of post-transcriptional modifications they are subject to, lead to the existence of complex repertoires of tRNA populations adjusted to defined cellular states. Despite the tRNA-coding genes redundancy in prokaryote and eukaryote genomes, it is surprising to note the absence of genes coding specific translational-active isoacceptors throughout the phylogeny. Through the analysis of different releases of tRNA databases, this review aims to provide a general summary about those “missing tRNA genes.” This absence refers to both tRNAs that are not encoded in the genome, as well as others that show critical sequence variations that would prevent their activity as canonical translation adaptor molecules. Notably, while a group of genes are universally missing, others are absent in particular kingdoms. Functional information available allows to hypothesize that the exclusion of isodecoding molecules would be linked to: 1) reduce ambiguities of signals that define the specificity of the interactions in which the tRNAs are involved; 2) ensure the adaptation of the translational apparatus to the cellular state; 3) divert particular tRNA variants from ribosomal protein synthesis to other cellular functions. This leads to consider the “missing tRNA genes” as a source of putative non-canonical tRNA functions and to broaden the concept of adapter molecules in ribosomal-dependent protein synthesis.

Structural analysis of leucine, lysine and tryptophan mitochondrial tRNA of nesting turtles Caretta caretta (Testudines: Chelonioidea) in the Colombian Caribbean

The understanding of the functional properties of mitochondrial transfer RNA (mt tRNAs) depend on the knowledge of its structure. tRNA acts as an interface between polynucleotides and polypeptides thus, they are key molecules in protein biosynthesis. The tRNA molecule has a functional design and, given its importance in the translation of mitochondrial genes, it is plausible that modifications of the structure can affect the synthesis of proteins and the functional properties of the mitochondria. In a previous work, the mitochondrial genome of an individual of the nesting Caretta caretta of the Colombian Caribbean was obtained, where specific mutations were identified in the only tRNALeu (CUN), tRNATrp and tRNALys genes. In order to analyze the effect of these mutations on these three mt tRNAs, the prediction of 2D and 3D structures was performed. Genes were sequenced in 11 nesting loggerhead turtles from the Colombian Caribbean. Two-dimensional structures were inferred using the ARWEN program, and three-dimensional structures were obtained with the RNA Composer 3D program. Two polymorphisms were identified in tRNATrp and another one was located in tRNALys, both specific to C. caretta. The thymine substitution in nucleotide position 14 of tRNATrp could constitute an endemic polymorphism of the nesting colony of the Colombian Caribbean. Two 2D and three 3D patterns were obtained for tRNATrp. In the case of tRNALys and tRNALeu 2D and 3D structures were obtained respectively, which showed compliance to canonical structures, with 4 bp in the D-arm, 4-5 bp in the T-arm, and 5 bp in the anticodon arm. Moderate deviations were found, such as a change in the number of nucleotides, elongation in loops or stems and non-Watson-Crick base pairing: adenine-adenine in stem D of tRNATrp, uracil-uracil and adenine-cytosine in the acceptor arm of the tRNALys and cytosine-cytosine in the anticodon stem of the tRNALeu. In addition, distortions or lack of typical interactions in 3D structures gave them unique characteristics. According to the size of the variable region (4-5 nt), the three analyzed tRNAs belong to class I. The interactions in the three studied tRNAs occur mainly between D loop-variable region, and between spacer bases-variable region, which classifies them as tRNA of typology II. The polymorphisms and structural changes described can, apparently, be post-transcriptionally stabilized. It will be crucial to perform studies at the population and functional levels to elucidate the synthetic pathways affected by these genes. This article analyses for the first time the 1D, 2D and 3D structures of the mitochondrial tRNALys, tRNATrp and tRNALeu in the loggerhead turtle.

The complete mitochondrial genome of a parasite at the animal-fungal boundary

BACKGROUND

Sphaerothecum destruens is an obligate intracellular fish parasite which has been identified as a serious threat to freshwater fishes. Taxonomically, S. destruens belongs to the order Dermocystida within the class Ichthyosporea (formerly referred to as Mesomycetozoea), which sits at the animal-fungal boundary. Mitochondrial DNA (mtDNA) sequences can be valuable genetic markers for species detection and are increasingly used in environmental DNA (eDNA) based species detection. Furthermore, mtDNA sequences can be used in epidemiological studies by informing detection, strain identification and geographical spread.

METHODS

We amplified the entire mitochondrial (mt) genome of S. destruens in two overlapping long fragments using primers designed based on the cox1, cob and nad5 partial sequences. The mt-genome architecture of S. destruens was then compared to close relatives to gain insights into its evolution.

RESULTS

The complete mt-genome of Sphaerothecum destruens is 23,939 bp in length and consists of 47 genes including 21 protein-coding genes, 2 rRNA, 22 tRNA and two unidentified open reading frames. The mitochondrial genome of S. destruens is intronless and compact with a few intergenic regions and includes genes that are often missing from animal and fungal mt-genomes, such as, the four ribosomal proteins (small subunit rps13 and 14; large subunit rpl2 and 16), tatC (twin-arginine translocase component C), and ccmC and ccmF (cytochrome c maturation protein ccmC and heme lyase).

CONCLUSIONS

We present the first mt-genome of S. destruens which also represents the first mt-genome for the order Dermocystida. The availability of the mt-genome can assist the detection of S. destruens and closely related parasites in eukaryotic diversity surveys using eDNA and assist epidemiological studies by improving molecular detection and tracking the parasite's spread. Furthermore, as the only representative of the order Dermocystida, its mt-genome can be used in the study of mitochondrial evolution of the unicellular relatives of animals.

Complete Mitogenome of a Leaf-Mining Buprestid Beetle, Trachys auricollis, and Its Phylogenetic Implications

A complete mitogenome of Trachys auricollis is reported, and a mitogenome-based phylogenetic tree of Elateriformia with all protein-coding genes (PCGs), rRNAs, and tRNAs is presented for the first time. The complete mitochondrial genome of T. auricollis is 16,429 bp in size and contains 13 PCGs, two rRNA genes, 22 tRNA genes, and an A + T-rich region. The A + T content of the entire genome is approximately 71.1%, and the AT skew and GC skew are 0.10 and -0.20, respectively. According to the the nonsynonymous substitution rate to synonymous substitution rates (Ka/Ks) of all PCGs, the highest and lowest evolutionary rates were observed for atp8 and cox1, respectively, which is a common finding among animals. The start codons of all PCGs are the typical ATN. Ten PCGs have complete stop codons, but three have incomplete stop codons with T or TA. As calculated based on the relative synonymous codon usage (RSCU) values, UUA(L) is the codon with the highest frequency. Except for trnS1, all 22 tRNA genes exhibit typical cloverleaf structures. The A + T-rich region of T. auricollis is located between rrnS and the trnI-trnG-trnM gene cluster, with six 72-bp tandem repeats. Both maximum likelihood (ML) and Bayesian (BI) trees suggest that Buprestoidea is close to Byrrhoidea and that Buprestoidea and Byrrhoidea are sister groups of Elateroidea, but the position of Psephenidae is undetermined. The inclusion of tRNAs might help to resolve the phylogeny of Coleoptera.

Microbiome characterization by high-throughput transfer RNA sequencing and modification analysis

Advances in high-throughput sequencing have facilitated remarkable insights into the diversity and functioning of naturally occurring microbes; however, current sequencing strategies are insufficient to reveal physiological states of microbial communities associated with protein translation dynamics. Transfer RNAs (tRNAs) are core components of protein synthesis machinery, present in all living cells, and are phylogenetically tractable, which make them ideal targets to gain physiological insights into environmental microbes. Here we report a direct sequencing approach, tRNA-seq, and a software suite, tRNA-seq-tools, to recover sequences, abundance profiles, and post-transcriptional modifications of microbial tRNA transcripts. Our analysis of cecal samples using tRNA-seq distinguishes high-fat- and low-fat-fed mice in a comparable fashion to 16S ribosomal RNA gene amplicons, and reveals taxon- and diet-dependent variations in tRNA modifications. Our results provide taxon-specific in situ insights into the dynamics of tRNA gene expression and post-transcriptional modifications within complex environmental microbiomes.

Rooted tRNAomes and evolution of the genetic code

We advocate for a tRNA- rather than an mRNA-centric model for evolution of the genetic code. The mechanism for evolution of cloverleaf tRNA provides a root sequence for radiation of tRNAs and suggests a simplified understanding of code evolution. To analyze code sectoring, rooted tRNAomes were compared for several archaeal and one bacterial species. Rooting of tRNAome trees reveals conserved structures, indicating how the code was shaped during evolution and suggesting a model for evolution of a LUCA tRNAome tree. We propose the polyglycine hypothesis that the initial product of the genetic code may have been short chain polyglycine to stabilize protocells. In order to describe how anticodons were allotted in evolution, the sectoring-degeneracy hypothesis is proposed. Based on sectoring, a simple stepwise model is developed, in which the code sectors from a 1→4→8→∼16 letter code. At initial stages of code evolution, we posit strong positive selection for wobble base ambiguity, supporting convergence to 4-codon sectors and ∼16 letters. In a later stage, ∼5–6 letters, including stops, were added through innovating at the anticodon wobble position. In archaea and bacteria, tRNA wobble adenine is negatively selected, shrinking the maximum size of the primordial genetic code to 48 anticodons. Because 64 codons are recognized in mRNA, tRNA-mRNA coevolution requires tRNA wobble position ambiguity leading to degeneracy of the code.

MeDIP-seq and nCpG analyses illuminate sexually dimorphic methylation of gonadal development genes with high historic methylation in turtle hatchlings with temperature-dependent sex determination

Background

DNA methylation alters gene expression but not DNA sequence and mediates some cases of phenotypic plasticity. Temperature-dependent sex determination (TSD) epitomizes phenotypic plasticity where environmental temperature drives embryonic sexual fate, as occurs commonly in turtles. Importantly, the temperature-specific transcription of two genes underlying gonadal differentiation is known to be induced by differential methylation in TSD fish, turtle and alligator. Yet, how extensive is the link between DNA methylation and TSD remains unclear. Here we test for broad differences in genome-wide DNA methylation between male and female hatchling gonads of the TSD painted turtle Chrysemys picta using methyl DNA immunoprecipitation sequencing, to identify differentially methylated candidates for future study. We also examine the genome-wide nCpG distribution (which affects DNA methylation) in painted turtles and test for historic methylation in genes regulating vertebrate gonadogenesis.

Results

Turtle global methylation was consistent with other vertebrates (57% of the genome, 78% of all CpG dinucleotides). Numerous genes predicted to regulate turtle gonadogenesis exhibited sex-specific methylation and were proximal to methylated repeats. nCpG distribution predicted actual turtle DNA methylation and was bimodal in gene promoters (as other vertebrates) and introns (unlike other vertebrates). Differentially methylated genes, including regulators of sexual development, had lower nCpG content indicative of higher historic methylation.

Conclusions

Ours is the first evidence suggesting that sexually dimorphic DNA methylation is pervasive in turtle gonads (perhaps mediated by repeat methylation) and that it targets numerous regulators of gonadal development, consistent with the hypothesis that it may regulate thermosensitive transcription in TSD vertebrates. However, further research during embryogenesis will help test this hypothesis and the alternative that instead, most differential methylation observed in hatchlings is the by-product of sexual differentiation and not its cause.

Electronic supplementary material

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

Mobile genes in the human microbiome are structured from global to individual scales

Recent work has underscored the importance of the microbiome in human health, largely attributing differences in phenotype to differences in the species present across individuals^1,2,3,4,5. But mobile genes can confer profoundly different phenotypes on different strains of the same species. Little is known about the function and distribution of mobile genes in the human microbiome, and in particular whether the gene pool is globally homogenous or constrained by human population structure. Here, we investigate this question by comparing the mobile genes found in the microbiomes of 81 metropolitan North Americans with that of 172 agrarian Fiji islanders using a combination of single-cell genomics and metagenomics. We find large differences in mobile gene content between the Fijian and North American microbiomes, with functional variation that mirrors known dietary differences such as the excess of plant-based starch degradation genes. Remarkably, differences are also observed between the mobile gene pools of proximal Fijian villages, even though microbiome composition across villages is similar. Finally, we observe high rates of recombination leading to individual-specific mobile elements, suggesting that the abundance of some genes may reflect environmental selection rather than dispersal limitation. Together, these data support the hypothesis that human activities and behaviors provide selective pressures that shape mobile gene pools, and that acquisition of mobile genes is important to colonizing specific human populations.

Evolution of initiator tRNAs and selection of methionine as the initiating amino acid

Transfer RNAs (tRNAs) have been important in shaping biomolecular evolution. Initiator tRNAs (tRNAi), a special class of tRNAs, carry methionine (or its derivative, formyl-methionine) to ribosomes to start an enormously energy consuming but a highly regulated process of protein synthesis. The processes of tRNAi evolution, and selection of methionine as the universal initiating amino acid remain an enigmatic problem. We constructed phylogenetic trees using the whole sequence, the acceptor-TψC arm ('minihelix'), and the anticodon-dihydrouridine arm regions of tRNAi from 158 species belonging to all 3 domains of life. All the trees distinctly assembled into 3 domains of life. Large trees, generated using data for all the tRNAs of a vast number of species, fail to reveal the major evolutionary events and identity of the probable elongator tRNA sequences that could be ancestor of tRNAi. Therefore, we constructed trees using the minihelix or the whole sequence of species specific tRNAs, and iterated our analysis on 50 eubacterial species. We identified tRNA(Pro), tRNA(Glu), or tRNA(Thr) (but surprisingly not elongator tRNA(Met)) as probable ancestors of tRNAi. We then determined the factors imposing selection of methionine as the initiating amino acid. Overall frequency of occurrence of methionine, whose metabolic cost of synthesis is the highest among all amino acids, remains almost unchanged across the 3 domains of life. Our correlation analysis shows that its high metabolic cost is independent of many physicochemical properties of the side chain. Our results indicate that selection of methionine, as the initiating amino acid was possibly a consequence of the evolution of one-carbon metabolism, which plays an important role in regulating translation initiation.

A Comprehensive tRNA Genomic Survey Unravels the Evolutionary History of tRNA Arrays in Prokaryotes

Considering the importance of tRNAs in the translation machinery, scant attention has been paid to tRNA array units defined as genomic regions containing at least 20 tRNA genes with a minimal tRNA gene density of two tRNA genes per kilobase. Our analysis of Acidithiobacillus ferrivorans CF27 and Acidithiobacillus ferrooxidans ATCC 23270(T) genomes showed that both display a tRNA array unit with syntenic conservation which mainly contributed to the tRNA gene redundancy in these two organisms. Our investigations into the occurrence and distribution of tRNA array units revealed that 1) this tRNA organization is limited to few phyla and mainly found in Gram-positive bacteria; and 2) the presence of tRNA arrays favors the redundancy of tRNA genes, in particular those encoding the core tRNA isoacceptors. Finally, comparative array organization revealed that tRNA arrays were acquired through horizontal gene transfer (from Firmicutes or unknown donor), before being subjected to tRNA rearrangements, deletions, and duplications. In Bacilli, the most parsimonious evolutionary history involved two common ancestors and the acquisition of their arrays arose late in evolution, in the genera branches. Functional roles of the array units in organism lifestyle, selective genetic advantage and translation efficiency, as well as the evolutionary advantages of organisms harboring them were proposed. Our study offers new insight into the structural organization and evolution of tRNA arrays in prokaryotic organisms.

Genomic analysis reveals the biotechnological and industrial potential of levan producing halophilic extremophile, Halomonas smyrnensis AAD6T

Halomonas smyrnensis AAD6T is a gram negative, aerobic, and moderately halophilic bacterium, and is known to produce high levels of levan with many potential uses in foods, feeds, cosmetics, pharmaceutical and chemical industries due to its outstanding properties. Here, the whole-genome analysis was performed to gain more insight about the biological mechanisms, and the whole-genome organization of the bacterium. Industrially crucial genes, including the levansucrase, were detected and the genome-scale metabolic model of H. smyrnensis AAD6T was reconstructed. The bacterium was found to have many potential applications in biotechnology not only being a levan producer, but also because of its capacity to produce Pel exopolysaccharide, polyhydroxyalkanoates, and osmoprotectants. The genomic information presented here will not only provide additional information to enhance our understanding of the genetic and metabolic network of halophilic bacteria, but also accelerate the research on systematical design of engineering strategies for biotechnology applications.

Eukaryotic tRNAs fingerprint invertebrates vis-à-vis vertebrates

During translation, aminoacyl-tRNA synthetases recognize the identities of the tRNAs to charge them with their respective amino acids. The conserved identities of 58,244 eukaryotic tRNAs of 24 invertebrates and 45 vertebrates in genomic tRNA database were analyzed and their novel features extracted. The internal promoter sequences, namely, A-Box and B-Box, were investigated and evidence gathered that the intervention of optional nucleotides at 17a and 17b correlated with the optimal length of the A-Box. The presence of canonical transcription terminator sequences at the immediate vicinity of tRNA genes was ventured. Even though non-canonical introns had been reported in red alga, green alga, and nucleomorph so far, fairly motivating evidence of their existence emerged in tRNA genes of other eukaryotes. Non-canonical introns were seen to interfere with the internal promoters in two cases, questioning their transcription fidelity. In a first of its kind, phylogenetic constructs based on tRNA molecules delineated and built the trees of the vast and diverse invertebrates and vertebrates. Finally, two tRNA models representing the invertebrates and the vertebrates were drawn, by isolating the dominant consensus in the positional fluctuations of nucleotide compositions.

tRNA signatures reveal a polyphyletic origin of SAR11 strains among alphaproteobacteria

Molecular phylogenetics and phylogenomics are subject to noise from horizontal gene transfer (HGT) and bias from convergence in macromolecular compositions. Extensive variation in size, structure and base composition of alphaproteobacterial genomes has complicated their phylogenomics, sparking controversy over the origins and closest relatives of the SAR11 strains. SAR11 are highly abundant, cosmopolitan aquatic Alphaproteobacteria with streamlined, A+T-biased genomes. A dominant view holds that SAR11 are monophyletic and related to both Rickettsiales and the ancestor of mitochondria. Other studies dispute this, finding evidence of a polyphyletic origin of SAR11 with most strains distantly related to Rickettsiales. Although careful evolutionary modeling can reduce bias and noise in phylogenomic inference, entirely different approaches may be useful to extract robust phylogenetic signals from genomes. Here we develop simple phyloclassifiers from bioinformatically derived tRNA Class-Informative Features (CIFs), features predicted to target tRNAs for specific interactions within the tRNA interaction network. Our tRNA CIF-based model robustly and accurately classifies alphaproteobacterial genomes into one of seven undisputed monophyletic orders or families, despite great variability in tRNA gene complement sizes and base compositions. Our model robustly rejects monophyly of SAR11, classifying all but one strain as Rhizobiales with strong statistical support. Yet remarkably, conventional phylogenetic analysis of tRNAs classifies all SAR11 strains identically as Rickettsiales. We attribute this discrepancy to convergence of SAR11 and Rickettsiales tRNA base compositions. Thus, tRNA CIFs appear more robust to compositional convergence than tRNA sequences generally. Our results suggest that tRNA-CIF-based phyloclassification is robust to HGT of components of the tRNA interaction network, such as aminoacyl-tRNA synthetases. We explain why tRNAs are especially advantageous for prediction of traits governing macromolecular interactions from genomic data, and why such traits may be advantageous in the search for robust signals to address difficult problems in classification and phylogeny.

If gene products work well in the networks of foreign cells, their genes may transfer horizontally between unrelated genomes. What factors dictate the ability to integrate into foreign networks? Different RNAs and proteins must interact specifically in order to function well as a system. For example, tRNA functions are determined by the interactions they have with other macromolecules. We have developed ways to predict, from genomic data alone, how tRNAs distinguish themselves to their specific interaction partners. Here, as proof of concept, we built a robust computational model from these bioinformatic predictions in seven lineages of Alphaproteobacteria. We validated our model by classifying hundreds of diverse alphaproteobacterial taxa and tested it on eight strains of SAR11, a phylogenetically controversial group that is highly abundant in the world's oceans. We found that different strains of SAR11 are more distantly related, both to each other and to mitochondria, than widely believed. We explain conflicting results about SAR11 as an artifact of bias created by the variability in base contents of alphaproteobacterial genomes. While this bias affects tRNAs too, our classifier appears unexpectedly robust to it. More broadly, our results suggest that traits governing macromolecular interactions may be more faithfully vertically inherited than the macromolecules themselves.

tRNA genes rapidly change in evolution to meet novel translational demands

Changes in expression patterns may occur when organisms are presented with new environmental challenges, for example following migration or genetic changes. To elucidate the mechanisms by which the translational machinery adapts to such changes, we perturbed the tRNA pool of Saccharomyces cerevisiae by tRNA gene deletion. We then evolved the deletion strain and observed that the genetic adaptation was recurrently based on a strategic mutation that changed the anticodon of other tRNA genes to match that of the deleted one. Strikingly, a systematic search in hundreds of genomes revealed that anticodon mutations occur throughout the tree of life. We further show that the evolution of the tRNA pool also depends on the need to properly couple translation to protein folding. Together, our observations shed light on the evolution of the tRNA pool, demonstrating that mutation in the anticodons of tRNA genes is a common adaptive mechanism when meeting new translational demands.

DOI:http://dx.doi.org/10.7554/eLife.01339.001

Genes contain the blueprints for the proteins that are essential for countless biological functions and processes, and the path that leads from a particular gene to the corresponding protein is long and complex. The genetic information stored in the DNA must first be transcribed to produce a messenger RNA molecule, which then has to be translated to produce a string of amino acids that fold to form a protein. The translation step is performed by a molecular machine called the ribosome, with transfer RNA molecules bringing the amino acids that are needed to make the protein.

The information in messenger RNA is stored as a series of letters, with groups of three letters called codons representing the different amino acids. Since there are four letters—A, C, G and U—it is possible to form 64 different codons. And since there are only 20 amino acids, two or more different codons can specify the same amino acid (for example, AGU and AGC both specify serine), and two or more different transfer RNA molecules can take this amino acid to the ribosome. Moreover, some codons are found more often than others in the messenger RNA molecules, so the genes that encode the related transfer RNA molecules are more common than the genes for other transfer RNA molecules.

Environmental pressures mean that organisms must adapt to survive, with some genes and proteins increasing in importance, and others becoming less important. Clearly the relative numbers of the different transfer RNA molecules will also need to change to reflect these evolutionary changes, but the details of how this happens were not understood.

Now Yona et al. have explored this issue by studying yeast cells that lack a gene for one of the less common transfer RNA molecules (corresponding to the codon AGG, which specifies the amino acid arginine). At first this mutation resulted in slower growth of the yeast cells, but after being allowed to evolve over 200 generations, the rate of growth matched that of a normal strain with all transfer RNA genes. Yona et al. found that the gene for a more common transfer RNA molecule, corresponding to the codon AGA, which also specifies arginine, had mutated to AGG. As a result, the mutated yeast was eventually able to produce proteins as quickly as wild type yeast. Moreover, further experiments showed that the levels of some transfer RNAs are kept deliberately low in order to slow down the production of proteins so as to ensure that the proteins assume their correct structure.

But does the way these cells evolved in the lab resemble what happened in nature? To address this question Yona et al. examined a database of transfer RNA sequences from more than 500 species, and found evidence for the same codon-based switching mechanism in many species across the tree of life.

DOI:http://dx.doi.org/10.7554/eLife.01339.002

Complete mitochondrial DNA sequences of the threadfin cichlid (Petrochromis trewavasae) and the blunthead cichlid (Tropheus moorii) and patterns of mitochondrial genome evolution in cichlid fishes

The cichlid fishes of the East African Great Lakes represent a model especially suited to study adaptive radiation and speciation. With several African cichlid genome projects being in progress, a promising set of closely related genomes is emerging, which is expected to serve as a valuable data base to solve questions on genotype-phenotype relations. The mitochondrial (mt) genomes presented here are the first results of the assembly and annotation process for two closely related but eco-morphologically highly distinct Lake Tanganyika cichlids, Petrochromis trewavasae and Tropheus moorii. The genomic sequences comprise 16,588 bp (P. trewavasae) and 16,590 bp (T. moorii), and exhibit the typical mitochondrial structure, with 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a non-coding control region. Analyses confirmed that the two species are very closely related with an overall sequence similarity of 96%. We analyzed the newly generated sequences in the phylogenetic context of 21 published labroid fish mitochondrial genomes. Consistent with other vertebrates, the D-loop region was found to evolve faster than protein-coding genes, which in turn are followed by the rRNAs; the tRNAs vary greatly in the rate of sequence evolution, but on average evolve the slowest. Within the group of coding genes, ND6 evolves most rapidly. Codon usage is similar among examined cichlid tribes and labroid families; although a slight shift in usage patterns down the gene tree could be observed. Despite having a clearly different nucleotide composition, ND6 showed a similar codon usage. C-terminal ends of Cox1 exhibit variations, where the varying number of amino acids is related to the structure of the obtained phylogenetic tree. This variation may be of functional relevance for Cox1 synthesis.

Life without tRNAArg-adenosine deaminase TadA: evolutionary consequences of decoding the four CGN codons as arginine in Mycoplasmas and other Mollicutes

In most bacteria, two tRNAs decode the four arginine CGN codons. One tRNA harboring a wobble inosine (tRNA(Arg)ICG) reads the CGU, CGC and CGA codons, whereas a second tRNA harboring a wobble cytidine (tRNA(Arg)CCG) reads the remaining CGG codon. The reduced genomes of Mycoplasmas and other Mollicutes lack the gene encoding tRNA(Arg)CCG. This raises the question of how these organisms decode CGG codons. Examination of 36 Mollicute genomes for genes encoding tRNA(Arg) and the TadA enzyme, responsible for wobble inosine formation, suggested an evolutionary scenario where tadA gene mutations first occurred. This allowed the temporary accumulation of non-deaminated tRNA(Arg)ACG, capable of reading all CGN codons. This hypothesis was verified in Mycoplasma capricolum, which contains a small fraction of tRNA(Arg)ACG with a non-deaminated wobble adenosine. Subsets of Mollicutes continued to evolve by losing both the mutated tRNA(Arg)CCG and tadA, and then acquired a new tRNA(Arg)UCG. This permitted further tRNA(Arg)ACG mutations with tRNA(Arg)GCG or its disappearance, leaving a single tRNA(Arg)UCG to decode the four CGN codons. The key point of our model is that the A-to-I deamination activity had to be controlled before the loss of the tadA gene, allowing the stepwise evolution of Mollicutes toward an alternative decoding strategy.

Ancient origin of the divergent forms of leucyl-tRNA synthetases in the Halobacteriales

Background

Horizontal gene transfer (HGT) has greatly impacted the genealogical history of many lineages, particularly for prokaryotes, with genes frequently moving in and out of a line of descent. Many genes that were acquired by a lineage in the past likely originated from ancestral relatives that have since gone extinct. During the course of evolution, HGT has played an essential role in the origin and dissemination of genetic and metabolic novelty.

Results

Three divergent forms of leucyl-tRNA synthetase (LeuRS) exist in the archaeal order Halobacteriales, commonly known as haloarchaea. Few haloarchaeal genomes have the typical archaeal form of this enzyme and phylogenetic analysis indicates it clusters within the Euryarchaeota as expected. The majority of sequenced halobacterial genomes possess a bacterial form of LeuRS. Phylogenetic reconstruction puts this larger group of haloarchaea at the base of the bacterial domain. The most parsimonious explanation is that an ancient transfer of LeuRS took place from an organism related to the ancestor of the bacterial domain to the haloarchaea. The bacterial form of LeuRS further underwent gene duplications and/or gene transfers within the haloarchaea, with some genomes possessing two distinct types of bacterial LeuRS. The cognate tRNA^Leu also reveals two distinct clusters for the haloarchaea; however, these tRNA^Leu clusters do not coincide with the groupings found in the LeuRS tree, revealing that LeuRS evolved independently of its cognate tRNA.

Conclusions

The study of leucyl-tRNA synthetase in haloarchaea illustrates the importance of gene transfer originating in lineages that went extinct since the transfer occurred. The haloarchaeal LeuRS and tRNA^Leu did not co-evolve.

Thermal, autonomous replicator made from transfer RNA

Evolving systems rely on the storage and replication of genetic information. Here we present an autonomous, purely thermally driven replication mechanism. A pool of hairpin molecules, derived from transfer RNA replicates the succession of a two-letter code. Energy is first stored thermally in metastable hairpins. Thereafter, energy is released by a highly specific and exponential replication with a duplication time of 30 s, which is much faster than the tendency to produce false positives in the absence of template. Our experiments propose a physical rather than a chemical scenario for the autonomous replication of protein encoding information in a disequilibrium setting.

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.