tRNA genes rapidly change in evolution to meet novel translational demands

How much is Transcription-associated Mutagenesis Driving tRNA Microevolution?

Transfer RNAs (tRNAs) are among the most highly conserved and frequently transcribed genes. Recent studies have identified transcription-associated mutagenesis (TAM) as a significant contributor to sequence variation around tRNA loci. However, the extent to which TAM drives allelic variation in tRNAs remains unclear, largely due to the confounding effects of strong selection pressures to maintain their structural integrity. This complexity arises because TAM-induced mutations primarily involve nucleotide transitions, which tend to preserve base-pairing stability. To address this dichotomy at the population level, we analyzed tRNA allelic variation in contemporary Caenorhabditis elegans strains. We propose a model of tRNA microevolution driven by TAM and demonstrate that the observed secondary structure characteristics align with our predicted TAM-biased patterns. Furthermore, we developed a continuous Markov substitution model that incorporates TAM-specific mutational biases. This TAM-biased model fits the C. elegans tRNA data more effectively than standard models, such as the general time-reversible (GTR) model. Based on these results, we conclude that TAM plays a significant role in shaping tRNA allelic variation within populations. This finding is consistent with recent experimental studies on tRNA fitness in yeast but challenges prior theoretical and computational analyses that emphasize RNA base-pairing as a primary determinant in genotype-phenotype systems.

Genomic insights into Paspalum vaginatum: Mitochondrial and chloroplast genome mapping, evolutionary insights, and organelle-nucleus communication

Paspalum vaginatum, valued for its salt tolerance, is a vital species in the turfgrass and agricultural industries. Despite its significance, there are still gaps in its genetic composition, particularly in the mitochondrial (mtDNA) and chloroplast (cpDNA) genomes. Our study aimed to fill these knowledge gaps by investigating the evolutionary relationships within the paspalum family and examining the functions of organelle-encoded genes as well as the critical role of reactive oxygen species (ROS) in organelle-nucleus communication. By genome sequencing, assembly, and annotation, we determined 504,515 bp of P. vaginatum mtDNA and 140,483 bp of its cpDNA. Comparative analyses with other Paspalum species and major crops highlight the intricate evolutionary dynamics and varying levels of genetic relatedness observed across different organelle genomes. The complex response of organelle gene expression to salt stress in this study will aid in understanding the molecular mechanisms and evolutionary trajectories of P. vaginatum organelle genomes.

A machine learning approach uncovers principles and determinants of eukaryotic ribosome pausing

Nonuniform local translation speed dictates diverse protein biogenesis outcomes. To unify known and uncover unknown principles governing eukaryotic elongation rate, we developed a machine learning pipeline to analyze RiboSeq datasets. We find that the chemical nature of the incoming amino acid determines how codon optimality influences elongation rate, with hydrophobic residues more dependent on transfer RNA (tRNA) levels than charged residues. Unexpectedly, we find that wobble interactions exert a widespread effect on elongation pausing, with wobble-mediated decoding being slower than Watson-Crick decoding, irrespective of tRNA levels. Applying our ribosome pausing principles to ribosome collisions reveals that disomes arise upon apposition of fast-decoding and slow-decoding signatures. We conclude that codon choice and tRNA pools are evolutionarily constrained to harmonize elongation rate with cotranslational folding while minimizing wobble pairing and deleterious stalling.

Multigenerational inheritance drives symbiotic interactions of the bacterium Bacillus subtilis with its plant host

Bacillus subtilis is a beneficial bacterium that supports plant growth and protects plants from bacterial, fungal, and viral infections. Using a simplified system of B. subtilis and Arabidopsis thaliana interactions, we studied the fitness and transcriptome of bacteria detached from the root over generations of growth in LB medium. We found that bacteria previously associated with the root or exposed to its secretions had greater stress tolerance and were more competitive in root colonization than bacteria not previously exposed to the root. Furthermore, our transcriptome results provide evidence that plant secretions induce a microbial stress response and fundamentally alter signaling by the cyclic nucleotide c-di-AMP, a signature maintained by their descendants. The changes in cellular physiology due to exposure to plant exudates were multigenerational, as they allowed not only the bacterial cells that colonized a new plant but also their descendants to have an advance over naive competitors of the same species, while the overall plasticity of gene expression and rapid adaptation were maintained. These changes were hereditary but not permanent. Our work demonstrates a bacterial memory manifested by multigenerational reversible adaptation to plant hosts in the form of activation of the stressosome, which confers an advantage to symbiotic bacteria during competition.

Establishing chromosomal design-build-test-learn through a synthetic chromosome and its combinatorial reconfiguration

Chromosome-level design-build-test-learn cycles (chrDBTLs) allow systematic combinatorial reconfiguration of chromosomes with ease. Here, we established chrDBTL with a redesigned synthetic Saccharomyces cerevisiae chromosome XV, synXV. We designed and built synXV to harbor strategically inserted features, modified elements, and synonymously recoded genes throughout the chromosome. Based on the recoded chromosome, we developed a method to enable chrDBTL: CRISPR-Cas9-mediated mitotic recombination with endoreduplication (CRIMiRE). CRIMiRE allowed the creation of customized wild-type/synthetic combinations, accelerating genotype-phenotype mapping and synthetic chromosome redesign. We also leveraged synXV as a "build-to-learn" model organism for translation studies by ribosome profiling. We conducted a locus-to-locus comparison of ribosome occupancy between synXV and the wild-type chromosome, providing insight into the effects of codon changes and redesigned features on translation dynamics in vivo. Overall, we established synXV as a versatile reconfigurable system that advances chrDBTL for understanding biological mechanisms and engineering strains.

Comprehensive quantitative modeling of translation efficiency in a genome-reduced bacterium

Translation efficiency has been mainly studied by ribosome profiling, which only provides an incomplete picture of translation kinetics. Here, we integrated the absolute quantifications of tRNAs, mRNAs, RNA half-lives, proteins, and protein half-lives with ribosome densities and derived the initiation and elongation rates for 475 genes (67% of all genes), 73 with high precision, in the bacterium Mycoplasma pneumoniae (Mpn). We found that, although the initiation rate varied over 160-fold among genes, most of the known factors had little impact on translation efficiency. Local codon elongation rates could not be fully explained by the adaptation to tRNA abundances, which varied over 100-fold among tRNA isoacceptors. We provide a comprehensive quantitative view of translation efficiency, which suggests the existence of unidentified mechanisms of translational regulation in Mpn.

Copy number variation in tRNA isodecoder genes impairs mammalian development and balanced translation

The number of tRNA isodecoders has increased dramatically in mammals, but the specific molecular and physiological reasons for this expansion remain elusive. To address this fundamental question we used CRISPR editing to knockout the seven-membered phenylalanine tRNA gene family in mice, both individually and combinatorially. Using ATAC-Seq, RNA-seq, ribo-profiling and proteomics we observed distinct molecular consequences of single tRNA deletions. We show that tRNA-Phe-1-1 is required for neuronal function and its loss is partially compensated by increased expression of other tRNAs but results in mistranslation. In contrast, the other tRNA-Phe isodecoder genes buffer the loss of each of the remaining six tRNA-Phe genes. In the tRNA-Phe gene family, the expression of at least six tRNA-Phe alleles is required for embryonic viability and tRNA-Phe-1-1 is most important for development and survival. Our results reveal that the multi-copy configuration of tRNA genes is required to buffer translation and viability in mammals.

Microevolutionary dynamics of eccDNA in Chinese hamster ovary cells grown in fed-batch cultures under control and lactate-stressed conditions

Chinese hamster ovary (CHO) cell lines are widely used to manufacture biopharmaceuticals. However, CHO cells are not an optimal expression host due to the intrinsic plasticity of the CHO genome. Genome plasticity can lead to chromosomal rearrangements, transgene exclusion, and phenotypic drift. A poorly understood genomic element of CHO cell line instability is extrachromosomal circular DNA (eccDNA) in gene expression and regulation. EccDNA can facilitate ultra-high gene expression and are found within many eukaryotes including humans, yeast, and plants. EccDNA confers genetic heterogeneity, providing selective advantages to individual cells in response to dynamic environments. In CHO cell cultures, maintaining genetic homogeneity is critical to ensuring consistent productivity and product quality. Understanding eccDNA structure, function, and microevolutionary dynamics under various culture conditions could reveal potential engineering targets for cell line optimization. In this study, eccDNA sequences were investigated at the beginning and end of two-week fed-batch cultures in an ambr®250 bioreactor under control and lactate-stressed conditions. This work characterized structure and function of eccDNA in a CHO-K1 clone. Gene annotation identified 1551 unique eccDNA genes including cancer driver genes and genes involved in protein production. Furthermore, RNA-seq data is integrated to identify transcriptionally active eccDNA genes.

Anticodon sequence determines the impact of mistranslating tRNAAla variants

Transfer RNAs (tRNAs) maintain translation fidelity through accurate charging by their cognate aminoacyl-tRNA synthetase and codon:anticodon base pairing with the mRNA at the ribosome. Mistranslation occurs when an amino acid not specified by the genetic message is incorporated into proteins and has applications in biotechnology, therapeutics and is relevant to disease. Since the alanyl-tRNA synthetase uniquely recognizes a G3:U70 base pair in tRNAAla and the anticodon plays no role in charging, tRNAAla variants with anticodon mutations have the potential to mis-incorporate alanine. Here, we characterize the impact of the 60 non-alanine tRNAAla anticodon variants on the growth of Saccharomyces cerevisiae. Overall, 36 tRNAAla anticodon variants decreased growth in single- or multi-copy. Mass spectrometry analysis of the cellular proteome revealed that 52 of 57 anticodon variants, not decoding alanine or stop codons, induced mistranslation when on single-copy plasmids. Variants with G/C-rich anticodons resulted in larger growth deficits than A/U-rich variants. In most instances, synonymous anticodon variants impact growth differently, with anticodons containing U at base 34 being the least impactful. For anticodons generating the same amino acid substitution, reduced growth generally correlated with the abundance of detected mistranslation events. Differences in decoding specificity, even between synonymous anticodons, resulted in each tRNAAla variant mistranslating unique sets of peptides and proteins. We suggest that these differences in decoding specificity are also important in determining the impact of tRNAAla anticodon variants.

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.

The Evolution of tRNA Copy Number and Repertoire in Cellular Life

tRNAs are universal decoders that bridge the gap between transcriptome and proteome. They can also be processed into small RNA fragments with regulatory functions. In this work, we show that tRNA copy number is largely controlled by genome size in all cellular organisms, in contrast to what is observed for protein-coding genes that stop expanding between ~20,000 and ~35,000 loci per haploid genome in eukaryotes, regardless of genome size. Our analyses indicate that after the bacteria/archaea split, the tRNA gene pool experienced the evolution of increased anticodon diversity in the archaeal lineage, along with a tRNA gene size increase and mature tRNA size decrease. The evolution and diversification of eukaryotes from archaeal ancestors involved further expansion of the tRNA anticodon repertoire, additional increase in tRNA gene size and decrease in mature tRNA length, along with an explosion of the tRNA gene copy number that emerged coupled with accelerated genome size expansion. Our findings support the notion that macroscopic eukaryotes with a high diversity of cell types, such as land plants and vertebrates, independently evolved a high diversity of tRNA anticodons along with high gene redundancy caused by the expansion of the tRNA copy number. The results presented here suggest that the evolution of tRNA genes played important roles in the early split between bacteria and archaea, and in eukaryogenesis and the later emergence of complex eukaryotes, with potential implications in protein translation and gene regulation through tRNA-derived RNA fragments.

tRNA derived small RNAs—Small players with big roles

In the past 2 decades, small non-coding RNAs derived from tRNA (tsRNAs or tRNA derived fragments; tRFs) have emerged as new powerful players in the field of small RNA mediated regulation of gene expression, translation, and epigenetic control. tRFs have been identified from evolutionarily divergent organisms from Archaea, the higher plants, to humans. Recent studies have confirmed their roles in cancers and other metabolic disorders in humans and experimental models. They have been implicated in biotic and abiotic stress responses in plants as well. In this review, we summarize the current knowledge on tRFs including types of tRFs, their biogenesis, and mechanisms of action. The review also highlights recent studies involving differential expression profiling of tRFs and elucidation of specific functions of individual tRFs from various species. We also discuss potential considerations while designing experiments involving tRFs identification and characterization and list the available bioinformatics tools for this purpose.

Preferred synonymous codons are translated more accurately: Proteomic evidence, among-species variation, and mechanistic basis

A commonly stated cause of unequal uses of synonymous codons is their differential translational accuracies. A key prediction of this long-standing translational accuracy hypothesis (TAH) of codon usage bias is higher translational accuracies of more frequently used synonymous codons, which, however, has had no direct evidence beyond case studies. Analyzing proteomic data from Escherichia coli, we present direct, global evidence for this prediction. The experimentally measured codon-specific translational accuracies validate a sequence-based proxy; this proxy provides support for the TAH from the vast majority of over 1000 taxa surveyed in all domains of life. We find that the relative translational accuracies of synonymous codons vary substantially among taxa and are strongly correlated with the amounts of cognate transfer RNAs (tRNAs) relative to those of near-cognate tRNAs. These and other observations suggest a model in which selections for translational efficiency and accuracy drive codon usage bias and its coevolution with the tRNA pool.

The Ramp Atlas: facilitating tissue and cell-specific ramp sequence analyses through an intuitive web interface

Ramp sequences occur when the average translational efficiency of codons near the 5′ end of highly expressed genes is significantly lower than the rest of the gene sequence, which counterintuitively increases translational efficiency by decreasing downstream ribosomal collisions. Here, we show that the relative codon adaptiveness within different tissues changes the existence of a ramp sequence without altering the underlying genetic code. We present the first comprehensive analysis of tissue and cell type-specific ramp sequences and report 3108 genes with ramp sequences that change between tissues and cell types, which corresponds with increased gene expression within those tissues and cells. The Ramp Atlas (https://ramps.byu.edu/) allows researchers to query precomputed ramp sequences in 18 388 genes across 62 tissues and 66 cell types and calculate tissue-specific ramp sequences from user-uploaded FASTA files through an intuitive web interface. We used The Ramp Atlas to identify seven SARS-CoV-2 genes and seven human SARS-CoV-2 entry factor genes with tissue-specific ramp sequences that may help explain viral proliferation within those tissues. We anticipate that The Ramp Atlas will facilitate personalized and creative tissue-specific ramp sequence analyses for both human and viral genes that will increase our ability to utilize this often-overlooked regulatory region.

Evolutionary repair reveals an unexpected role of the tRNA modification m1G37 in aminoacylation

The tRNA modification m¹G37, introduced by the tRNA methyltransferase TrmD, is thought to be essential for growth in bacteria because it suppresses translational frameshift errors at proline codons. However, because bacteria can tolerate high levels of mistranslation, it is unclear why loss of m¹G37 is not tolerated. Here, we addressed this question through experimental evolution of trmD mutant strains of Escherichia coli. Surprisingly, trmD mutant strains were viable even if the m¹G37 modification was completely abolished, and showed rapid recovery of growth rate, mainly via duplication or mutation of the proline-tRNA ligase gene proS. Growth assays and in vitro aminoacylation assays showed that G37-unmodified tRNA^Pro is aminoacylated less efficiently than m¹G37-modified tRNA^Pro, and that growth of trmD mutant strains can be largely restored by single mutations in proS that restore aminoacylation of G37-unmodified tRNA^Pro. These results show that inefficient aminoacylation of tRNA^Pro is the main reason for growth defects observed in trmD mutant strains and that proS may act as a gatekeeper of translational accuracy, preventing the use of error-prone unmodified tRNA^Pro in translation. Our work shows the utility of experimental evolution for uncovering the hidden functions of essential genes and has implications for the development of antibiotics targeting TrmD.

Genomic Changes and Genetic Divergence of Vibrio alginolyticus Under Phage Infection Stress Revealed by Whole-Genome Sequencing and Resequencing

Bacteriophages (phages) and their bacterial hosts were the most abundant and genetically highly diverse organisms on the earth. In this study, a series of phage-resistant mutant (PRM) strains derived from Vibrio alginolyticus were isolated and Infrequent-restriction-site PCR (IRS-PCR) was used to investigate the genetic diversity of the PRM strains. Phenotypic variations of eight PRM strains were analyzed using profiles of utilizing carbon sources and chemical sensitivity. Genetic variations of eight PRM strains and coevolved V. alginolyticus populations with phages were analyzed by whole-genome sequencing and resequencing, respectively. The results indicated that eight genetically discrepant PRM stains exhibited abundant and abundant phenotypic variations. Eight PRM strains and coevolved V. alginolyticus populations (VE1, VE2, and VE3) contained numerous single nucleotide variations (SNVs) and insertions/indels (InDels) and exhibited obvious genetic divergence. Most of the SNVs and InDels in coding genes were related to the synthesis of flagellar, extracellular polysaccharide (EPS), which often served as the receptors of phage invasion. The PRM strains and the coevolved cell populations also contained frequent mutations in tRNA and rRNA genes. Two out of three coevolved populations (VE1 and VE2) contained a large mutation segment severely deconstructing gene nrdA, which was predictably responsible for the booming of mutation rate in the genome. In summary, numerous mutations and genetic divergence were detected in the genomes of V. alginolyticus PRM strains and in coevolved cell populations of V. alginolyticus under phage infection stress. The phage infection stress may provide an important force driving genomic evolution of V. alginolyticus.

Regulation of gene expression via translational buffering

Translation of an mRNA represents a critical step during the expression of protein-coding genes. As mechanisms governing post-transcriptional regulation of gene expression are progressively unveiled, it is becoming apparent that transcriptional programs are not fully reflected in the proteome. Herein, we highlight a previously underappreciated post-transcriptional mode of regulation of gene expression termed translational buffering. In principle, translational buffering opposes the impact of alterations in mRNA levels on the proteome. We further describe three types of translational buffering: compensation, which maintains protein levels e.g. across species or individuals; equilibration, which retains pathway stoichiometry; and offsetting, which acts as a reversible mechanism that maintains the levels of selected subsets of proteins constant despite genetic alteration and/or stress-induced changes in corresponding mRNA levels. While mechanisms underlying compensation and equilibration have been reviewed elsewhere, the principal focus of this review is on the less-well understood mechanism of translational offsetting. Finally, we discuss potential roles of translational buffering in homeostasis and disease.

Evolutionary responses to codon usage of horizontally transferred genes in Pseudomonas aeruginosa: gene retention, amelioration and compensatory evolution

Prokaryote genome evolution is characterized by the frequent gain of genes through horizontal gene transfer (HGT). For a gene, being horizontally transferred can represent a strong change in its genomic and physiological context. If the codon usage of a transferred gene deviates from that of the receiving organism, the fitness benefits it provides can be reduced due to a mismatch with the expression machinery. Consequently, transferred genes with a deviating codon usage can be selected against or elicit evolutionary responses that enhance their integration, such as gene amelioration and compensatory evolution. Within bacterial species, the extent and relative importance of these different mechanisms has never been considered altogether. In this study, a phylogeny-based method was used to investigate the occurrence of these different evolutionary responses in Pseudomonas aeruginosa. Selection on codon usage of genes acquired through HGT was observed over evolutionary time, with the overall codon usage converging towards that of the core genome. Gene amelioration, through the accumulation of synonymous mutations after HGT, did not seem to systematically affect transferred genes. This pattern therefore seemed to be mainly driven by selective retention of transferred genes with an initial codon usage similar to that of the core genes. Additionally, variation in the copy number of tRNA genes was often associated with the acquisition of genes for which the observed variation could enhance their expression. This provides evidence that compensatory evolution might be an important mechanism for the integration of horizontally transferred genes.

Homopeptide and homocodon levels across fungi are coupled to GC/AT-bias and intrinsic disorder, with unique behaviours for some amino acids

Homopeptides (runs of one amino-acid type) are evolutionarily important since they are prone to expand/contract during DNA replication, recombination and repair. To gain insight into the genomic/proteomic traits driving their variation, we analyzed how homopeptides and homocodons (which are pure codon repeats) vary across 405 Dikarya, and probed their linkage to genome GC/AT bias and other factors. We find that amino-acid homopeptide frequencies vary diversely between clades, with the AT-rich Saccharomycotina trending distinctly. As organisms evolve, homocodon and homopeptide numbers are majorly coupled to GC/AT-bias, exhibiting a bi-furcated correlation with degree of AT- or GC-bias. Mid-GC/AT genomes tend to have markedly fewer simply because they are mid-GC/AT. Despite these trends, homopeptides tend to be GC-biased relative to other parts of coding sequences, even in AT-rich organisms, indicating they absorb AT bias less or are inherently more GC-rich. The most frequent and most variable homopeptide amino acids favour intrinsic disorder, and there are an opposing correlation and anti-correlation versus homopeptide levels for intrinsic disorder and structured-domain content respectively. Specific homopeptides show unique behaviours that we suggest are linked to inherent slippage probabilities during DNA replication and recombination, such as poly-glutamine, which is an evolutionarily very variable homopeptide with a codon repertoire unbiased for GC/AT, and poly-lysine whose homocodons are overwhelmingly made from the codon AAG.

tRNAs as a Driving Force of Genome Evolution in Yeast

Transfer RNAs (tRNAs) are widely known for their roles in the decoding of the linear mRNA information into amino acid sequences of proteins. They are also multifunctional platforms in the translation process and have other roles beyond translation, including sensing amino acid abundance, interacting with the general stress response machinery, and modulating cellular adaptation, survival, and death. In this mini-review, we focus on the emerging role of tRNA genes in the organization and modification of the genomic architecture of yeast and the role of tRNA misexpression and decoding infidelity in genome stability, evolution, and adaption. We discuss published work showing how quickly tRNA genes can mutate to meet novel translational demands, how tRNAs speed up genome evolution, and how tRNA genes can be sites of genomic instability. We highlight recent works showing that loss of tRNA decoding fidelity and small alterations in tRNA expression have unexpected and profound impacts on genome stability. By dissecting these recent evidence, we hope to lay the groundwork that prompts future investigations on the mechanistic interplay between tRNAs and genome modification that likely triggers genome evolution.

Whole-Genome Transformation Promotes tRNA Anticodon Suppressor Mutations under Stress

tRNAs are encoded by a large gene family, usually with several isogenic tRNAs interacting with the same codon. Mutations in the anticodon region of other tRNAs can overcome specific tRNA deficiencies. Phylogenetic analysis suggests that such mutations have occurred in evolution, but the driving force is unclear. We show that in yeast suppressor mutations in other tRNAs are able to overcome deficiency of the essential TRT2-encoded tRNAThrCGU at high temperature (40°C). Surprisingly, these tRNA suppressor mutations were obtained after whole-genome transformation with DNA from thermotolerant Kluyveromyces marxianus or Ogataea polymorpha strains but from which the mutations did apparently not originate. We suggest that transient presence of donor DNA in the host facilitates proliferation at high temperature and thus increases the chances for occurrence of spontaneous mutations suppressing defective growth at high temperature. Whole-genome sequence analysis of three transformants revealed only four to five nonsynonymous mutations of which one causing TRT2 anticodon stem stabilization and two anticodon mutations in non-threonyl-tRNAs, tRNALysCUU and tRNAeMetCAU, were causative. Both anticodon mutations suppressed lethality of TRT2 deletion and apparently caused the respective tRNAs to become novel substrates for threonyl-tRNA synthetase. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) data could not detect any significant mistranslation, and reverse transcription-quantitative PCR results contradicted induction of the unfolded protein response. We suggest that stress conditions have been a driving force in evolution for the selection of anticodon-switching mutations in tRNAs as revealed by phylogenetic analysis.IMPORTANCE In this work, we have identified for the first time the causative elements in a eukaryotic organism introduced by applying whole-genome transformation and responsible for the selectable trait of interest, i.e., high temperature tolerance. Surprisingly, the whole-genome transformants contained just a few single nucleotide polymorphisms (SNPs), which were unrelated to the sequence of the donor DNA. In each of three independent transformants, we have identified a SNP in a tRNA, either stabilizing the essential tRNAThrCGU at high temperature or switching the anticodon of tRNALysCUU or tRNAeMetCAU into CGU, which is apparently enough for in vivo recognition by threonyl-tRNA synthetase. LC-MS/MS analysis indeed indicated absence of significant mistranslation. Phylogenetic analysis showed that similar mutations have occurred throughout evolution and we suggest that stress conditions may have been a driving force for their selection. The low number of SNPs introduced by whole-genome transformation may favor its application for improvement of industrial yeast strains.

Protein Complexes Form a Basis for Complex Hybrid Incompatibility

Proteins are the workhorses of the cell and execute many of their functions by interacting with other proteins forming protein complexes. Multi-protein complexes are an admixture of subunits, change their interaction partners, and modulate their functions and cellular physiology in response to environmental changes. When two species mate, the hybrid offspring are usually inviable or sterile because of large-scale differences in the genetic makeup between the two parents causing incompatible genetic interactions. Such reciprocal-sign epistasis between inter-specific alleles is not limited to incompatible interactions between just one gene pair; and, usually involves multiple genes. Many of these multi-locus incompatibilities show visible defects, only in the presence of all the interactions, making it hard to characterize. Understanding the dynamics of protein-protein interactions (PPIs) leading to multi-protein complexes is better suited to characterize multi-locus incompatibilities, compared to studying them with traditional approaches of genetics and molecular biology. The advances in omics technologies, which includes genomics, transcriptomics, and proteomics can help achieve this end. This is especially relevant when studying non-model organisms. Here, we discuss the recent progress in the understanding of hybrid genetic incompatibility; omics technologies, and how together they have helped in characterizing protein complexes and in turn multi-locus incompatibilities. We also review advances in bioinformatic techniques suitable for this purpose and propose directions for leveraging the knowledge gained from model-organisms to identify genetic incompatibilities in non-model organisms.

Partition of tRNAGly isoacceptors between protein and cell-wall peptidoglycan synthesis in Staphylococcus aureus

The sequence of tRNAs is submitted to evolutionary constraints imposed by their multiple interactions with aminoacyl-tRNA synthetases, translation elongation factor Tu in complex with GTP (EF-Tu•GTP), and the ribosome, each being essential for accurate and effective decoding of messenger RNAs. In Staphylococcus aureus, an additional constraint is imposed by the participation of tRNA^Gly isoacceptors in the addition of a pentaglycine side chain to cell-wall peptidoglycan precursors by transferases FmhB, FemA and FemB. Three tRNA^Gly isoacceptors poorly interacting with EF-Tu•GTP and the ribosome were previously identified. Here, we show that these ‘non-proteogenic’ tRNAs are preferentially recognized by FmhB based on kinetic analyses and on synthesis of stable aminoacyl-tRNA analogues acting as inhibitors. Synthesis of chimeric tRNAs and of helices mimicking the tRNA acceptor arms revealed that this discrimination involves identity determinants exclusively present in the D and T stems and loops of non-proteogenic tRNAs, which belong to an evolutionary lineage only present in the staphylococci. EF-Tu•GTP competitively inhibited FmhB by sequestration of ‘proteogenic’ aminoacyl-tRNAs in vitro. Together, these results indicate that competition for the Gly-tRNA^Gly pool is restricted by both limited recognition of non-proteogenic tRNAs by EF-Tu•GTP and limited recognition of proteogenic tRNAs by FmhB.

Sphingomonas changnyeongensis sp. nov. isolated from the Hapcheon–Changnyeong barrage area in the Nakdong river

The novel bacterial strain C33T was isolated from a freshwater sample collected from the Hapcheon–Changnyeong barrage. The Gram-negative, motile, yellow-pigmented strain C33T was characterized as a rod-shaped and strictly aerobic bacterium. A 16S-rRNA phylogenetic analysis revealed that this strain was most closely related to Sphingomonas changbaiensis V2M44T, Sphingomonas tabacisoli X1-8T, and Sphingomonas flavalba ZLT-5T with 97.1, 97.0, and 95.0 % 16S-rRNA sequence similarities, respectively. The genomic DNA GC content of strain C33T was estimated at 65.0 mol%. The average nucleotide identity of strain C33T relative to S. changbaiensis V2M44T and S. flavalba ZLT-5T was found to be 77.0 and 75.6%, with average amino-acid identities of 69.9, and 66.7%, and the digital DNA–DNA hybridization values of 21.3 and 17.7 %, respectively. The cells grew at 19–37 °C and pH 6–9 with 0–0.5 % (w/v) NaCl (optimum: 28 °C, pH 6.5, and 0 % NaCl). The major component identified in the polyamine pattern was sym-homospermidine, and the main ubiquinone was Q-10. The predominant polar lipids characterized were diphophatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, phosphatidylmonomethylethanolamine, phosphatidyldimethylethanolamine, and sphingoglycolipid. Iso-C15 : 0, C15 : 0 anteiso, and summed feature 3 (C16 : 1 ω6c and/or C16 : 1 ω7c) were found to be the primary cellular fatty acids in strain C33T. Based on these genotypic and phenotypic characteristics, strain C33T was classified as a novel species of the genus Sphingomonas ; and the name Sphingomonas changnyeongensis sp. nov. is proposed (=KACC 21511T=JCM 33880T).

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.

Alanine tRNAs Translate Environment Into Behavior in Caenorhabditis elegans

Caenorhabditis elegans nematodes produce and maintain imprints of attractive chemosensory cues to which they are exposed early in life. Early odor-exposure increases adult chemo-attraction to the same cues. Imprinting is transiently or stably inherited, depending on the number of exposed generations. We show here that the Alanine tRNA (UGC) plays a central role in regulating C. elegans chemo-attraction. Naive worms fed on tRNAAla (UGC) purified from odor-experienced worms, acquire odor-specific imprints. Chemo-attractive responses require the tRNA-modifying Elongator complex sub-units 1 (elpc-1) and 3 (elpc-3) genes. elpc-3 deletions impair chemo-attraction, which is fully restored by wild-type tRNAAla (UGC) feeding. A stably inherited decrease of odor-specific responses ensues from early odor-exposition of elpc-1 deletion mutants. tRNAAla (UGC) may adopt various chemical forms to mediate the cross-talk between innately-programmed and environment-directed chemo-attractive behavior.

Interchangeable parts: The evolutionarily dynamic tRNA population in plant mitochondria

Transfer RNAs (tRNAs) remain one of the very few classes of genes still encoded in the mitochondrial genome. These key components of the protein translation system must interact with a large enzymatic network of nuclear-encoded gene products to maintain mitochondrial function. Plants have an evolutionarily dynamic mitochondrial tRNA population, including ongoing tRNA gene loss and replacement by both horizontal gene transfer from diverse sources and import of nuclear-expressed tRNAs from the cytosol. Thus, plant mitochondria represent an excellent model for understanding how anciently divergent genes can act as "interchangeable parts" during the evolution of complex molecular systems. In particular, understanding the integration of the mitochondrial translation system with elements of the corresponding machinery used in cytosolic protein synthesis is a key area for eukaryotic cellular evolution. Here, we review the increasingly detailed phylogenetic data about the evolutionary history of mitochondrial tRNA gene loss, transfer, and functional replacement that has created extreme variation in mitochondrial tRNA populations across plant species. We describe emerging tRNA-seq methods with promise for refining our understanding of the expression and subcellular localization of tRNAs. Finally, we summarize current evidence and identify open questions related to coevolutionary changes in nuclear-encoded enzymes that have accompanied turnover in mitochondrial tRNA populations.

Elucidation of Codon Usage Signatures across the Domains of Life

Because of the degeneracy of the genetic code, multiple codons are translated into the same amino acid. Despite being “synonymous,” these codons are not equally used. Selective pressures are thought to drive the choice among synonymous codons within a genome, while GC content, which is typically attributed to mutational drift, is the major determinant of variation across species. Here, we find that in addition to GC content, interspecies codon usage signatures can also be detected. More specifically, we show that a single amino acid, arginine, is the major contributor to codon usage bias differences across domains of life. We then exploit this finding and show that domain-specific codon bias signatures can be used to classify a given sequence into its corresponding domain of life with high accuracy. We then wondered whether the inclusion of codon usage codon autocorrelation patterns, which reflects the nonrandom distribution of codon occurrences throughout a transcript, might improve the classification performance of our algorithm. However, we find that autocorrelation patterns are not domain-specific, and surprisingly, are unrelated to tRNA reusage, in contrast to previous reports. Instead, our results suggest that codon autocorrelation patterns are a by-product of codon optimality throughout a sequence, where highly expressed genes display autocorrelated “optimal” codons, whereas lowly expressed genes display autocorrelated “nonoptimal” codons.

Hidden cases of tRNA gene duplication and remolding in mitochondrial genomes of amphipods

The evolution of tRNA genes in mitochondrial (mt) genomes is a complex process that includes duplications, degenerations, and transpositions, as well as a specific process of identity change through mutations in the anticodon (tRNA gene remolding or tRNA gene recruitment). Using amphipod-specific tRNA models for annotation, we show that tRNA duplications are more common in the mt genomes of amphipods than what was revealed by previous annotations. Seventeen cases of tRNA gene duplications were detected in the mt genomes of amphipods, and ten of them were tRNA genes that underwent remolding. The additional tRNA gene findings were verified using phylogenetic analysis and genetic distance analysis. The majority of remolded tRNA genes (seven out of ten cases) were found in the mt genomes of endemic amphipod species from Lake Baikal. All additional mt tRNA genes arose independently in the Baikalian amphipods, indicating the unusual plasticity of tRNA gene evolution in these species assemblages. The possible reasons for the unusual abundance of additional tRNA genes in the mt genomes of Baikalian amphipods are discussed. The amphipod-specific tRNA models developed for MiTFi refine existing predictions of tRNA genes in amphipods and reveal additional cases of duplicated tRNA genes overlooked by using less specific Metazoa-wide models. The application of these models for mt tRNA gene prediction will be useful for the correct annotation of mt genomes of amphipods and probably other crustaceans.

Evolution of intron splicing towards optimized gene expression is based on various Cis- and Trans-molecular mechanisms

Splicing expands, reshapes, and regulates the transcriptome of eukaryotic organisms. Despite its importance, key questions remain unanswered, including the following: Can splicing evolve when organisms adapt to new challenges? How does evolution optimize inefficiency of introns’ splicing and of the splicing machinery? To explore these questions, we evolved yeast cells that were engineered to contain an inefficiently spliced intron inside a gene whose protein product was under selection for an increased expression level. We identified a combination of mutations in Cis (within the gene of interest) and in Trans (in mRNA-maturation machinery). Surprisingly, the mutations in Cis resided outside of known intronic functional sites and improved the intron’s splicing efficiency potentially by easing tight mRNA structures. One of these mutations hampered a protein’s domain that was not under selection, demonstrating the evolutionary flexibility of multi-domain proteins as one domain functionality was improved at the expense of the other domain. The Trans adaptations resided in two proteins, Npl3 and Gbp2, that bind pre-mRNAs and are central to their maturation. Interestingly, these mutations either increased or decreased the affinity of these proteins to mRNA, presumably allowing faster spliceosome recruitment or increased time before degradation of the pre-mRNAs, respectively. Altogether, our work reveals various mechanistic pathways toward optimizations of intron splicing to ultimately adapt gene expression patterns to novel demands.

An experimental evolution study involving an inefficiently spliced intron reveals that the splicing machinery, introns, and RNA quality control factors evolve in Cis and in Trans when cells optimize their transcriptome to new challenges.

Indole Derivatives Maintain the Status Quo Between Beneficial Biofilms and Their Plant Hosts

Biofilms formed by bacteria on plant roots play an important role in maintaining an optimal rhizosphere environment that supports plant growth and fitness. Bacillus subtilis is a potent plant growth promoter, forming biofilms that play a key role in protecting the host from fungal and bacterial infections. In this work, we demonstrate that the development of B. subtilis biofilms is antagonized by specific indole derivatives that accumulate during symbiotic interactions with plant hosts. Indole derivatives are more potent signals when the plant polysaccharide xylan serves as a carbon source, a mechanism to sustain beneficial biofilms at a biomass that can be supported by the plant. Moreover, B. subtilis biofilms formed by mutants resistant to indole derivatives become deleterious to the plants due to their capacity to consume and recycle plant polysaccharides. These results demonstrate how a dynamic metabolite-based dialogue can promote homeostasis between plant hosts and their beneficial biofilm communities.

Evolthon: A community endeavor to evolve lab evolution

In experimental evolution, scientists evolve organisms in the lab, typically by challenging them to new environmental conditions. How best to evolve a desired trait? Should the challenge be applied abruptly, gradually, periodically, sporadically? Should one apply chemical mutagenesis, and do strains with high innate mutation rate evolve faster? What are ideal population sizes of evolving populations? There are endless strategies, beyond those that can be exposed by individual labs. We therefore arranged a community challenge, Evolthon, in which students and scientists from different labs were asked to evolve Escherichia coli or Saccharomyces cerevisiae for an abiotic stress—low temperature. About 30 participants from around the world explored diverse environmental and genetic regimes of evolution. After a period of evolution in each lab, all strains of each species were competed with one another. In yeast, the most successful strategies were those that used mating, underscoring the importance of sex in evolution. In bacteria, the fittest strain used a strategy based on exploration of different mutation rates. Different strategies displayed variable levels of performance and stability across additional challenges and conditions. This study therefore uncovers principles of effective experimental evolutionary regimens and might prove useful also for biotechnological developments of new strains and for understanding natural strategies in evolutionary arms races between species. Evolthon constitutes a model for community-based scientific exploration that encourages creativity and cooperation.

This Community Page article describes Evolthon; a first-of-its-kind community-based effort, involving about 30 participant labs around the world, aiming to explore the best strategy for evolving microorganisms to cope with an environmental challenge.

The stop-and-go traffic regulating protein biogenesis: How translation kinetics controls proteostasis

Generating a functional proteome requires the ribosome to carefully regulate disparate co-translational processes that determine the fate of nascent polypeptides. With protein synthesis being energetically expensive, the ribosome must balance the costs of efficiently making a protein with those of properly folding it. Emerging as a primary means of regulating this trade-off is the nonuniform rate of translation elongation that defines translation kinetics. The varying speeds with which the ribosome progresses along a transcript have been implicated in several aspects of protein biogenesis, including co-translational protein folding and translational fidelity, as well as gene expression by mediating mRNA decay and protein quality control pathways. The optimal translation kinetics required to efficiently execute these processes can be distinct. Thus, the ribosome is tasked with tightly regulating translation kinetics to balance these processes while maintaining adaptability for changing cellular conditions. In this review, we first discuss the regulatory role of translation elongation in protein biogenesis and what factors influence elongation kinetics. We then describe how changes in translation kinetics signal downstream pathways that dictate the fate of nascent polypeptides. By regulating these pathways, the kinetics of translation elongation has emerged as a critical tool for driving gene expression and maintaining proteostasis through varied mechanisms, including nascent chain folding and binding different ribosome-associated machinery. Indeed, a growing number of examples demonstrate the important role of local changes in elongation kinetics in modulating the pathophysiology of human disease.

Cells alter their tRNA abundance to selectively regulate protein synthesis during stress conditions

Decoding the information in mRNA during protein synthesis relies on tRNA adaptors, the abundance of which can affect the decoding rate and translation efficiency. To determine whether cells alter tRNA abundance to selectively regulate protein expression, we quantified changes in the abundance of individual tRNAs at different time points in response to diverse stress conditions in Saccharomyces cerevisiae. We found that the tRNA pool was dynamic and rearranged in a manner that facilitated selective translation of stress-related transcripts. Through genomic analysis of multiple data sets, stochastic simulations, and experiments with designed sequences of proteins with identical amino acids but altered codon usage, we showed that changes in tRNA abundance affected protein expression independently of factors such as mRNA abundance. We suggest that cells alter their tRNA abundance to selectively affect the translation rates of specific transcripts to increase the amounts of required proteins under diverse stress conditions.

Repertoires of tRNAs: The Couplers of Genomics and Proteomics

The pool of transfer RNA (tRNA) molecules in cells allows the ribosome to decode genetic information. This repertoire of molecular decoders is positioned in the crossroad of the genome, the transcriptome, and the proteome. Omics and systems biology now allow scientists to explore the entire repertoire of tRNAs of many organisms, revealing basic exciting biology. The tRNA gene set of hundreds of species is now characterized, in addition to the tRNA genes of organelles and viruses. Genes encoding tRNAs for certain anticodon types appear in dozens of copies in a genome, while others are universally absent from any genome. Transcriptome measurement of tRNAs is challenging, but in recent years new technologies have allowed researchers to determine the dynamic expression patterns of tRNAs. These advances reveal that availability of ready-to-translate tRNA molecules is highly controlled by several transcriptional and posttranscriptional regulatory processes. This regulation shapes the proteome according to the cellular state. The tRNA pool profoundly impacts many aspects of cellular and organismal life, including protein expression level, translation accuracy, adequacy of folding, and even mRNA stability. As a result, the shape of the tRNA pool affects organismal health and may participate in causing conditions such as cancer and neurological conditions.

Evolutionary instability of CUG-Leu in the genetic code of budding yeasts

The genetic code used in nuclear genes is almost universal, but here we report that it changed three times in parallel during the evolution of budding yeasts. All three changes were reassignments of the codon CUG, which is translated as serine (in 2 yeast clades), alanine (1 clade), or the ‘universal’ leucine (2 clades). The newly discovered Ser2 clade is in the final stages of a genetic code transition. Most species in this clade have genes for both a novel tRNA^Ser(CAG) and an ancestral tRNA^Leu(CAG) to read CUG, but only tRNA^Ser(CAG) is used in standard growth conditions. The coexistence of these alloacceptor tRNA genes indicates that the genetic code transition occurred via an ambiguous translation phase. We propose that the three parallel reassignments of CUG were not driven by natural selection in favor of their effects on the proteome, but by selection to eliminate the ancestral tRNA^Leu(CAG).

The genetic code for amino acids is nearly universal, and among eukaryotic nuclear genomes the only known reassignments are of codon CUG in yeasts. Here, the authors identify a third independent CUG transition in budding yeasts that is still ongoing with alternative tRNAs present in the genome.

Codon usage of highly expressed genes affects proteome-wide translation efficiency

Although the genetic code is redundant, synonymous codons for the same amino acid are not used with equal frequencies in genomes, a phenomenon termed "codon usage bias." Previous studies have demonstrated that synonymous changes in a coding sequence can exert significant cis effects on the gene's expression level. However, whether the codon composition of a gene can also affect the translation efficiency of other genes has not been thoroughly explored. To study how codon usage bias influences the cellular economy of translation, we massively converted abundant codons to their rare synonymous counterpart in several highly expressed genes in Escherichia coli This perturbation reduces both the cellular fitness and the translation efficiency of genes that have high initiation rates and are naturally enriched with the manipulated codon, in agreement with theoretical predictions. Interestingly, we could alleviate the observed phenotypes by increasing the supply of the tRNA for the highly demanded codon, thus demonstrating that the codon usage of highly expressed genes was selected in evolution to maintain the efficiency of global protein translation.

Pausing of Chloroplast Ribosomes Is Induced by Multiple Features and Is Linked to the Assembly of Photosynthetic Complexes

Many mRNAs contain pause sites that briefly interrupt the progress of translation. Specific features that induce ribosome pausing have been described; however, their individual contributions to pause-site formation, and the overall biological significance of ribosome pausing, remain largely unclear. We have taken advantage of the compact genome of chloroplasts to carry out a plastid genome-wide survey of pause sites, as a basis for studying the impact of pausing on posttranslational processes. Based on ribosomal profiling of Arabidopsis (Arabidopsis thaliana) chloroplast mRNAs, we demonstrate that a combination of factors-mRNA secondary structure, internal Shine-Dalgarno sequences, and positively charged amino acids in the nascent peptide chain-explains 95% of the major pause sites on plastid mRNAs, whereas codon usage has little impact. The distribution of the pause sites is nonrandom and conforms to distinct patterns in the vicinity of sequences coding for transmembrane domains, which depend on their orientation within the membrane as well as being next to sequences coding for cofactor binding sites. We found strong indications that the mechanisms causing ribosomal pausing and at least some of the ribosomes pause sites are conserved between distantly related plant species. In addition, the positions of features that cause pausing are well conserved in photoautotrophic plants, but less so in their nonphotosynthetic, parasitic relatives, implying that the synthesis and assembly of photosynthetic multiprotein complexes requires localized ribosome pausing.

Bug mapping and fitness testing of chemically synthesized chromosome X

Debugging a genome sequence is imperative for successfully building a synthetic genome. As part of the effort to build a designer eukaryotic genome, yeast synthetic chromosome X (synX), designed as 707,459 base pairs, was synthesized chemically. SynX exhibited good fitness under a wide variety of conditions. A highly efficient mapping strategy called pooled PCRTag mapping (PoPM), which can be generalized to any watermarked synthetic chromosome, was developed to identify genetic alterations that affect cell fitness ("bugs"). A series of bugs were corrected that included a large region bearing complex amplifications, a growth defect mapping to a recoded sequence in FIP1, and a loxPsym site affecting promoter function of ATP2 PoPM is a powerful tool for synthetic yeast genome debugging and an efficient strategy for phenotype-genotype mapping.

Evolving Mistranslating tRNAs Through a Phenotypically Ambivalent Intermediate in Saccharomyces cerevisiae

The genetic code converts information from nucleic acid into protein. The genetic code was thought to be immutable, yet many examples in nature indicate that variations to the code provide a selective advantage. We used a sensitive selection system involving suppression of a deleterious allele (tti2-L187P) in Saccharomyces cerevisiae to detect mistranslation and identify mechanisms that allow genetic code evolution. Though tRNA^Ser containing a proline anticodon (UGG) is toxic, using our selection system we identified four tRNA^Ser_UGG variants, each with a single mutation, that mistranslate at a tolerable level. Mistranslating tRNA^Leu_UGG variants were also obtained, demonstrating the generality of the approach. We characterized two of the tRNA^Ser_UGG variants. One contained a G26A mutation, which reduced cell growth to 70% of the wild-type rate, induced a heat shock response, and was lost in the absence of selection. The reduced toxicity of tRNA^Ser_UGG-G26A is likely through increased turnover of the tRNA, as lack of methylation at G26 leads to degradation via the rapid tRNA decay pathway. The second tRNA^Ser_UGG variant, with a G9A mutation, had minimal effect on cell growth, was relatively stable in cells, and gave rise to less of a heat shock response. In vitro, the G9A mutation decreases aminoacylation and affects folding of the tRNA. Notably, the G26A and G9A mutations were phenotypically neutral in the context of an otherwise wild-type tRNA^Ser These experiments reveal a model for genetic code evolution in which tRNA anticodon mutations and mistranslation evolve through phenotypically ambivalent intermediates that reduce tRNA function.

Selfing in Haploid Plants and Efficacy of Selection: Codon Usage Bias in the Model Moss Physcomitrella patens

A long-term reduction in effective population size will lead to major shift in genome evolution. In particular, when effective population size is small, genetic drift becomes dominant over natural selection. The onset of self-fertilization is one evolutionary event considerably reducing effective size of populations. Theory predicts that this reduction should be more dramatic in organisms capable for haploid than for diploid selfing. Although theoretically well-grounded, this assertion received mixed experimental support. Here, we test this hypothesis by analyzing synonymous codon usage bias of genes in the model moss Physcomitrella patens frequently undergoing haploid selfing. In line with population genetic theory, we found that the effect of natural selection on synonymous codon usage bias is very weak. Our conclusion is supported by four independent lines of evidence: 1) Very weak or nonsignificant correlation between gene expression and codon usage bias, 2) no increased codon usage bias in more broadly expressed genes, 3) no evidence that codon usage bias would constrain synonymous and nonsynonymous divergence, and 4) predominant role of genetic drift on synonymous codon usage predicted by a model-based analysis. These findings show striking similarity to those observed in AT-rich genomes with weak selection for optimal codon usage and GC content overall. Our finding is in contrast to a previous study reporting adaptive codon usage bias in the moss P. patens.

APOBEC3B cytidine deaminase targets the non-transcribed strand of tRNA genes in yeast

Variations in mutation rates across the genome have been demonstrated both in model organisms and in cancers. This phenomenon is largely driven by the damage specificity of diverse mutagens and the differences in DNA repair efficiency in given genomic contexts. Here, we demonstrate that the single-strand DNA-specific cytidine deaminase APOBEC3B (A3B) damages tRNA genes at a 1000-fold higher efficiency than other non-tRNA genomic regions in budding yeast. We found that A3B-induced lesions in tRNA genes were predominantly located on the non-transcribed strand, while no transcriptional strand bias was observed in protein coding genes. Furthermore, tRNA gene mutations were exacerbated in cells where RNaseH expression was completely abolished (Δrnh1Δrnh35). These data suggest a transcription-dependent mechanism for A3B-induced tRNA gene hypermutation. Interestingly, in strains proficient in DNA repair, only 1% of the abasic sites formed upon excision of A3B-deaminated cytosines were not repaired leading to mutations in tRNA genes, while 18% of these lesions failed to be repaired in the remainder of the genome. A3B-induced mutagenesis in tRNA genes was found to be efficiently suppressed by the redundant activities of both base excision repair (BER) and the error-free DNA damage bypass pathway. On the other hand, deficiencies in BER did not have a profound effect on A3B-induced mutations in CAN1, the reporter for protein coding genes. We hypothesize that differences in the mechanisms underlying ssDNA formation at tRNA genes and other genomic loci are the key determinants of the choice of the repair pathways and consequently the efficiency of DNA damage repair in these regions. Overall, our results indicate that tRNA genes are highly susceptible to ssDNA-specific DNA damaging agents. However, increased DNA repair efficacy in tRNA genes can prevent their hypermutation and maintain both genome and proteome homeostasis.

Factors That Shape Eukaryotic tRNAomes: Processing, Modification and Anticodon-Codon Use

Transfer RNAs (tRNAs) contain sequence diversity beyond their anticodons and the large variety of nucleotide modifications found in all kingdoms of life. Some modifications stabilize structure and fit in the ribosome whereas those to the anticodon loop modulate messenger RNA (mRNA) decoding activity more directly. The identities of tRNAs with some universal anticodon loop modifications vary among distant and parallel species, likely to accommodate fine tuning for their translation systems. This plasticity in positions 34 (wobble) and 37 is reflected in codon use bias. Here, we review convergent evidence that suggest that expansion of the eukaryotic tRNAome was supported by its dedicated RNA polymerase III transcription system and coupling to the precursor-tRNA chaperone, La protein. We also review aspects of eukaryotic tRNAome evolution involving G34/A34 anticodon-sparing, relation to A34 modification to inosine, biased codon use and regulatory information in the redundancy (synonymous) component of the genetic code. We then review interdependent anticodon loop modifications involving position 37 in eukaryotes. This includes the eukaryote-specific tRNA modification, 3-methylcytidine-32 (m³C₃₂) and the responsible gene, TRM140 and homologs which were duplicated and subspecialized for isoacceptor-specific substrates and dependence on i⁶A₃₇ or t⁶A₃₇. The genetics of tRNA function is relevant to health directly and as disease modifiers.

Gene Architectures that Minimize Cost of Gene Expression

Gene expression burdens cells by consuming resources and energy. While numerous studies have investigated regulation of expression level, little is known about gene design elements that govern expression costs. Here, we ask how cells minimize production costs while maintaining a given protein expression level and whether there are gene architectures that optimize this process. We measured fitness of ∼14,000 E. coli strains, each expressing a reporter gene with a unique 5' architecture. By comparing cost-effective and ineffective architectures, we found that cost per protein molecule could be minimized by lowering transcription levels, regulating translation speeds, and utilizing amino acids that are cheap to synthesize and that are less hydrophobic. We then examined natural E. coli genes and found that highly expressed genes have evolved more forcefully to minimize costs associated with their expression. Our study thus elucidates gene design elements that improve the economy of protein expression in natural and heterologous systems.

Emerging Themes in Regulation of Global mRNA Turnover in cis

mRNA is the molecule that conveys genetic information from DNA to the translation apparatus. mRNAs in all organisms display a wide range of stability, and mechanisms have evolved to selectively and differentially regulate individual mRNA stability in response to intracellular and extracellular cues. In recent years, three seemingly distinct aspects of RNA biology-mRNA N⁶-methyladenosine (m6A) modification, alternative 3' end processing and polyadenylation (APA), and mRNA codon usage-have been linked to mRNA turnover, and all three aspects function to regulate global mRNA stability in cis. Here, we discuss the discovery and molecular dissection of these mechanisms in relation to how they impact the intrinsic decay rate of mRNA in eukaryotes, leading to transcriptome reprogramming.

Design, synthesis, and testing toward a 57-codon genome

Recoding--the repurposing of genetic codons--is a powerful strategy for enhancing genomes with functions not commonly found in nature. Here, we report computational design, synthesis, and progress toward assembly of a 3.97-megabase, 57-codon Escherichia coli genome in which all 62,214 instances of seven codons were replaced with synonymous alternatives across all protein-coding genes. We have validated 63% of recoded genes by individually testing 55 segments of 50 kilobases each. We observed that 91% of tested essential genes retained functionality with limited fitness effect. We demonstrate identification and correction of lethal design exceptions, only 13 of which were found in 2229 genes. This work underscores the feasibility of rewriting genomes and establishes a framework for large-scale design, assembly, troubleshooting, and phenotypic analysis of synthetic organisms.

Orthologs, turn-over, and remolding of tRNAs in primates and fruit flies

Background

Transfer RNAs (tRNAs) are ubiquitous in all living organism. They implement the genetic code so that most genomes contain distinct tRNAs for almost all 61 codons. They behave similar to mobile elements and proliferate in genomes spawning both local and non-local copies. Most tRNA families are therefore typically present as multicopy genes. The members of the individual tRNA families evolve under concerted or rapid birth-death evolution, so that paralogous copies maintain almost identical sequences over long evolutionary time-scales. To a good approximation these are functionally equivalent. Individual tRNA copies thus are evolutionary unstable and easily turn into pseudogenes and disappear. This leads to a rapid turnover of tRNAs and often large differences in the tRNA complements of closely related species. Since tRNA paralogs are not distinguished by sequence, common methods cannot not be used to establish orthology between tRNA genes.

Results

In this contribution we introduce a general framework to distinguish orthologs and paralogs in gene families that are subject to concerted evolution. It is based on the use of uniquely aligned adjacent sequence elements as anchors to establish syntenic conservation of sequence intervals. In practice, anchors and intervals can be extracted from genome-wide multiple sequence alignments. Syntenic clusters of concertedly evolving genes of different families can then be subdivided by list alignments, leading to usually small clusters of candidate co-orthologs. On the basis of recent advances in phylogenetic combinatorics, these candidate clusters can be further processed by cograph editing to recover their duplication histories. We developed a workflow that can be conceptualized as stepwise refinement of a graph of homologous genes. We apply this analysis strategy with different types of synteny anchors to investigate the evolution of tRNAs in primates and fruit flies. We identified a large number of tRNA remolding events concentrated at the tips of the phylogeny. With one notable exception all phylogenetically old tRNA remoldings do not change the isoacceptor class.

Conclusions

Gene families evolving under concerted evolution are not amenable to classical phylogenetic analyses since paralogs maintain identical, species-specific sequences, precluding the estimation of correct gene trees from sequence differences. This leaves conservation of syntenic arrangements with respect to “anchor elements” that are not subject to concerted evolution as the only viable source of phylogenetic information. We have demonstrated here that a purely synteny-based analysis of tRNA gene histories is indeed feasible. Although the choice of synteny anchors influences the resolution in particular when tight gene clusters are present, and the quality of sequence alignments, genome assemblies, and genome rearrangements limits the scope of the analysis, largely coherent results can be obtained for tRNAs. In particular, we conclude that a large fraction of the tRNAs are recent copies. This proliferation is compensated by rapid pseudogenization as exemplified by many very recent alloacceptor remoldings.

Electronic supplementary material

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