Auxiliary tRNAs: large-scale analysis of tRNA genes reveals patterns of tRNA repertoire dynamics

Examining the Effects of Temperature on the Evolution of Bacterial tRNA Pools

The genetic code consists of 61 codons coding for 20 amino acids. These codons are recognized by transfer RNAs (tRNAs) that bind to specific codons during protein synthesis. All organisms utilize less than all 61 possible anticodons due to base pair wobble: the ability to have a mismatch with a codon at its third nucleotide. Previous studies observed a correlation between the tRNA pool of bacteria and the temperature of their respective environments. However, it is unclear if these patterns represent biological adaptations to maintain the efficiency and accuracy of protein synthesis in different environments. A mechanistic mathematical model of mRNA translation is used to quantify the expected elongation rates and error rate for each codon based on an organism's tRNA pool. A comparative analysis across a range of bacteria that accounts for covariance due to shared ancestry is performed to quantify the impact of environmental temperature on the evolution of the tRNA pool. We find that thermophiles generally have more anticodons represented in their tRNA pool than mesophiles or psychrophiles. Based on our model, this increased diversity is expected to lead to increased missense errors. The implications of this for protein evolution in thermophiles are discussed.

Genome engineering on size reduction and complexity simplification: A review

BACKGROUND

Genome simplification is an important topic in the field of life sciences that has attracted attention from its conception to the present day. It can help uncover the essential components of the genome and, in turn, shed light on the underlying operating principles of complex biological systems. This has made it a central focus of both basic and applied research in the life sciences. With the recent advancements in related technologies and our increasing knowledge of the genome, now is an opportune time to delve into this topic.

AIM OF REVIEW

Our review investigates the progress of genome simplification from two perspectives: genome size reduction and complexity simplification. In addition, we provide insights into the future development trends of genome simplification.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Reducing genome size requires eliminating non-essential elements as much as possible. This process has been facilitated by advances in genome manipulation and synthesis techniques. However, we still need a better and clearer understanding of living systems to reduce genome complexity. As there is a lack of quantitative and clearly defined standards for this task, we have opted to approach the topic from various perspectives and present our findings accordingly.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Determining the effects of temperature on the evolution of bacterial tRNA pools

The genetic code consists of 61 codon coding for 20 amino acids. These codons are recognized by transfer RNAs (tRNA) that bind to specific codons during protein synthesis. Most organisms utilize less than all 61 possible anticodons due to base pair wobble: the ability to have a mismatch with a codon at its third nucleotide. Previous studies observed a correlation between the tRNA pool of bacteria and the temperature of their respective environments. However, it is unclear if these patterns represent biological adaptations to maintain the efficiency and accuracy of protein synthesis in different environments. A mechanistic mathematical model of mRNA translation is used to quantify the expected elongation rates and error rate for each codon based on an organism's tRNA pool. A comparative analysis across a range of bacteria that accounts for covariance due to shared ancestry is performed to quantify the impact of environmental temperature on the evolution of the tRNA pool. We find that thermophiles generally have more anticodons represented in their tRNA pool than mesophiles or psychrophiles. Based on our model, this increased diversity is expected to lead to increased missense errors. The implications of this for protein evolution in thermophiles are discussed.

The layered costs and benefits of translational redundancy

The rate and accuracy of translation hinges upon multiple components - including transfer RNA (tRNA) pools, tRNA modifying enzymes, and rRNA molecules - many of which are redundant in terms of gene copy number or function. It has been hypothesized that the redundancy evolves under selection, driven by its impacts on growth rate. However, we lack empirical measurements of the fitness costs and benefits of redundancy, and we have poor a understanding of how this redundancy is organized across components. We manipulated redundancy in multiple translation components of Escherichia coli by deleting 28 tRNA genes, 3 tRNA modifying systems, and 4 rRNA operons in various combinations. We find that redundancy in tRNA pools is beneficial when nutrients are plentiful and costly under nutrient limitation. This nutrient-dependent cost of redundant tRNA genes stems from upper limits to translation capacity and growth rate, and therefore varies as a function of the maximum growth rate attainable in a given nutrient niche. The loss of redundancy in rRNA genes and tRNA modifying enzymes had similar nutrient-dependent fitness consequences. Importantly, these effects are also contingent upon interactions across translation components, indicating a layered hierarchy from copy number of tRNA and rRNA genes to their expression and downstream processing. Overall, our results indicate both positive and negative selection on redundancy in translation components, depending on a species' evolutionary history with feasts and famines.

On distinguishing between canonical tRNA genes and tRNA gene fragments in prokaryotes

Automated genome annotation is essential for extracting biological information from sequence data. The identification and annotation of tRNA genes is frequently performed by the software package tRNAscan-SE, the output of which is listed for selected genomes in the Genomic tRNA database (GtRNAdb). Here, we highlight a pervasive error in prokaryotic tRNA gene sets on GtRNAdb: the mis-categorization of partial, non-canonical tRNA genes as standard, canonical tRNA genes. Firstly, we demonstrate the issue using the tRNA gene sets of 20 organisms from the archaeal taxon Thermococcaceae. According to GtRNAdb, these organisms collectively deviate from the expected set of tRNA genes in 15 instances, including the listing of eleven putative canonical tRNA genes. However, after detailed manual annotation, only one of these eleven remains; the others are either partial, non-canonical tRNA genes resulting from the integration of genetic elements or CRISPR-Cas activity (seven instances), or attributable to ambiguities in input sequences (three instances). Secondly, we show that similar examples of the mis-categorization of predicted tRNA sequences occur throughout the prokaryotic sections of GtRNAdb. While both canonical and non-canonical prokaryotic tRNA gene sequences identified by tRNAscan-SE are biologically interesting, the challenge of reliably distinguishing between them remains. We recommend employing a combination of (i) screening input sequences for the genetic elements typically associated with non-canonical tRNA genes, and ambiguities, (ii) activating the tRNAscan-SE automated pseudogene detection function, and (iii) scrutinizing predicted tRNA genes with low isotype scores. These measures greatly reduce manual annotation efforts, and lead to improved prokaryotic tRNA gene set predictions.

Non-Redundant tRNA Reference Sequences for Deep Sequencing Analysis of tRNA Abundance and Epitranscriptomic RNA Modifications

Analysis of RNA by deep-sequencing approaches has found widespread application in modern biology. In addition to measurements of RNA abundance under various physiological conditions, such techniques are now widely used for mapping and quantification of RNA modifications. Transfer RNA (tRNA) molecules are among the frequent targets of such investigation, since they contain multiple modified residues. However, the major challenge in tRNA examination is related to a large number of duplicated and point-mutated genes encoding those RNA molecules. Moreover, the existence of multiple isoacceptors/isodecoders complicates both the analysis and read mapping. Existing databases for tRNA sequencing provide near exhaustive listings of tRNA genes, but the use of such highly redundant reference sequences in RNA-seq analyses leads to a large number of ambiguously mapped sequencing reads. Here we describe a relatively simple computational strategy for semi-automatic collapsing of highly redundant tRNA datasets into a non-redundant collection of reference tRNA sequences. The relevance of the approach was validated by analysis of experimentally obtained tRNA-sequencing datasets for different prokaryotic and eukaryotic model organisms. The data demonstrate that non-redundant tRNA reference sequences allow improving unambiguous mapping of deep sequencing data.

Construction of anti-codon table of the plant kingdom and evolution of tRNA selenocysteine (tRNASec)

Background

The tRNAs act as a bridge between the coding mRNA and incoming amino acids during protein translation. The anti-codon of tRNA recognizes the codon of the mRNA and deliver the amino acid into the protein translation chain. However, we did not know about the exact abundance of anti-codons in the genome and whether the frequency of abundance remains same across the plant lineage or not.

Results

Therefore, we analysed the tRNAnome of 128 plant species and reported an anti-codon table of the plant kingdom. We found that CAU anti-codon of tRNA^Met has highest (5.039%) whereas GCG anti-codon of tRNA^Arg has lowest (0.004%) abundance. However, when we compared the anti-codon frequencies according to the tRNA isotypes, we found tRNA^Leu (7.808%) has highest abundance followed by tRNA^Ser (7.668%) and tRNA^Gly (7.523%). Similarly, suppressor tRNA (0.036%) has lowest abundance followed by tRNA^Sec (0.066%) and tRNA^His (2.109). The genome of Ipomoea nil, Papaver somniferum, and Zea mays encoded the highest number of anti-codons (isoacceptor) at 59 each whereas the genome of Ostreococcus tauri was found to encode only 18 isoacceptors. The tRNA^Sec genes undergone losses more frequently than duplication and we found that tRNA^Sec showed anti-codon switch during the course of evolution.

Conclusion

The anti-codon table of the plant tRNA will enable us to understand the synonymous codon usage of the plant kingdom and can be very helpful to understand which codon is preferred over other during the translation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07216-3.

The birth of a bacterial tRNA gene by large-scale, tandem duplication events

Organisms differ in the types and numbers of tRNA genes that they carry. While the evolutionary mechanisms behind tRNA gene set evolution have been investigated theoretically and computationally, direct observations of tRNA gene set evolution remain rare. Here, we report the evolution of a tRNA gene set in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The growth defect caused by deleting the single-copy tRNA gene, serCGA, is rapidly compensated by large-scale (45–290 kb) duplications in the chromosome. Each duplication encompasses a second, compensatory tRNA gene (serTGA) and is associated with a rise in tRNA-Ser(UGA) in the mature tRNA pool. We postulate that tRNA-Ser(CGA) elimination increases the translational demand for tRNA-Ser(UGA), a pressure relieved by increasing serTGA copy number. This work demonstrates that tRNA gene sets can evolve through duplication of existing tRNA genes, a phenomenon that may contribute to the presence of multiple, identical tRNA gene copies within genomes.

Identification and Characterization of Genes Required for 5-Hydroxyuridine Synthesis in Bacillus subtilis and Escherichia coli tRNA

In bacteria, tRNAs that decode 4-fold degenerate family codons and have uridine at position 34 of the anticodon are typically modified with either 5-methoxyuridine (mo⁵U) or 5-methoxycarbonylmethoxyuridine (mcmo⁵U). These modifications are critical for extended recognition of some codons at the wobble position. Whereas the alkylation steps of these modifications have been described, genes required for the hydroxylation of U34 to give 5-hydroxyuridine (ho⁵U) remain unknown. Here, a number of genes in Escherichia coli and Bacillus subtilis are identified that are required for wild-type (wt) levels of ho⁵U. The yrrMNO operon is identified in B. subtilis as important for the biosynthesis of ho⁵U. Both yrrN and yrrO are homologs to peptidase U32 family genes, which includes the rlhA gene required for ho⁵C synthesis in E. coli Deletion of either yrrN or yrrO, or both, gives a 50% reduction in mo⁵U tRNA levels. In E. coli, yegQ was found to be the only one of four peptidase U32 genes involved in ho⁵U synthesis. Interestingly, this mutant shows the same 50% reduction in (m)cmo⁵U as that observed for mo⁵U in the B. subtilis mutants. By analyzing the genomic context of yegQ homologs, the ferredoxin YfhL is shown to be required for ho⁵U synthesis in E. coli to the same extent as yegQ Additional genes required for Fe-S biosynthesis and biosynthesis of prephenate give the same 50% reduction in modification. Together, these data suggest that ho⁵U biosynthesis in bacteria is similar to that of ho⁵C, but additional genes and substrates are required for complete modification.IMPORTANCE Modified nucleotides in tRNA serve to optimize both its structure and function for accurate translation of the genetic code. The biosynthesis of these modifications has been fertile ground for uncovering unique biochemistry and metabolism in cells. In this work, genes that are required for a novel anaerobic hydroxylation of uridine at the wobble position of some tRNAs are identified in both Bacillus subtilis and Escherichia coli These genes code for Fe-S cluster proteins, and their deletion reduces the levels of the hydroxyuridine by 50% in both organisms. Additional genes required for Fe-S cluster and prephenate biosynthesis and a previously described ferredoxin gene all display a similar reduction in hydroxyuridine levels, suggesting that still other genes are required for the modification.

Wobbling Forth and Drifting Back: The Evolutionary History and Impact of Bacterial tRNA Modifications

Along with tRNAs, enzymes that modify anticodon bases are a key aspect of translation across the tree of life. tRNA modifications extend wobble pairing, allowing specific (“target”) tRNAs to recognize multiple codons and cover for other (“nontarget”) tRNAs, often improving translation efficiency and accuracy. However, the detailed evolutionary history and impact of tRNA modifying enzymes has not been analyzed. Using ancestral reconstruction of five tRNA modifications across 1093 bacteria, we show that most modifications were ancestral to eubacteria, but were repeatedly lost in many lineages. Most modification losses coincided with evolutionary shifts in nontarget tRNAs, often driven by increased bias in genomic GC and associated codon use, or by genome reduction. In turn, the loss of tRNA modifications stabilized otherwise highly dynamic tRNA gene repertoires. Our work thus traces the complex history of bacterial tRNA modifications, providing the first clear evidence for their role in the evolution of bacterial translation.

A-to-I editing on tRNAs: biochemical, biological and evolutionary implications

Inosine on transfer RNAs (tRNAs) are post-transcriptionally formed by a deamination mechanism of adenosines at positions 34, 37 and 57 of certain tRNAs. Despite its ubiquitous nature, the biological role of inosine in tRNAs remains poorly understood. Recent developments in the study of nucleotide modifications are beginning to indicate that the dynamics of such modifications are used in the control of specific genetic programs. Likewise, the essentiality of inosine-modified tRNAs in genome evolution and animal biology is becoming apparent. Here we review our current understanding on the role of inosine in tRNAs, the enzymes that catalyze the modification and the evolutionary link between such enzymes and other deaminases.