On the evolution and expression of Chlamydomonas reinhardtii nucleus-encoded transfer RNA genes

Advances in methods for tRNA sequencing and quantification

In the past decade tRNA sequencing (tRNA-seq) has attracted considerable attention as an important tool for the development of novel approaches to quantify highly modified tRNA species and to propel tRNA research aimed at understanding the cellular physiology and disease and development of tRNA-based therapeutics. Many methods are available to quantify tRNA abundance while accounting for modifications and tRNA charging/acylation. Advances in both library preparation methods and bioinformatic workflows have enabled developments in next-generation sequencing (NGS) workflows. Other approaches forgo NGS applications in favor of hybridization-based approaches. In this review we provide a brief comparative overview of various tRNA quantification approaches, focusing on the advantages and disadvantages of these methods, which together facilitate reliable tRNA quantification.

Adaptation of codon and amino acid use for translational functions in highly expressed cricket genes

BACKGROUND

For multicellular organisms, much remains unknown about the dynamics of synonymous codon and amino acid use in highly expressed genes, including whether their use varies with expression in different tissue types and sexes. Moreover, specific codons and amino acids may have translational functions in highly transcribed genes, that largely depend on their relationships to tRNA gene copies in the genome. However, these relationships and putative functions are poorly understood, particularly in multicellular systems.

RESULTS

Here, we studied codon and amino acid use in highly expressed genes from reproductive and nervous system tissues (male and female gonad, somatic reproductive system, brain and ventral nerve cord, and male accessory glands) in the cricket Gryllus bimaculatus. We report an optimal codon, defined as the codon preferentially used in highly expressed genes, for each of the 18 amino acids with synonymous codons in this organism. The optimal codons were mostly shared among tissue types and both sexes. However, the frequency of optimal codons was highest in gonadal genes. Concordant with translational selection, a majority of the optimal codons had abundant matching tRNA gene copies in the genome, but sometimes obligately required wobble tRNAs. We suggest the latter may comprise a mechanism for slowing translation of abundant transcripts, particularly for cell-cycle genes. Non-optimal codons, defined as those least commonly used in highly transcribed genes, intriguingly often had abundant tRNAs, and had elevated use in a subset of genes with specialized functions (gametic and apoptosis genes), suggesting their use promotes the translational upregulation of particular mRNAs. In terms of amino acids, we found evidence suggesting that amino acid frequency, tRNA gene copy number, and amino acid biosynthetic costs (size/complexity) had all interdependently evolved in this insect model, potentially for translational optimization.

CONCLUSIONS

Collectively, the results suggest a model whereby codon use in highly expressed genes, including optimal, wobble, and non-optimal codons, and their tRNA abundances, as well as amino acid use, have been influenced by adaptation for various functional roles in translation within this cricket. The effects of expression in different tissue types and the two sexes are discussed.

Evidence of multifaceted functions of codon usage in translation within the model beetle Tribolium castaneum

Synonymous codon use is non-random. Codons most used in highly transcribed genes, often called optimal codons, typically have high gene counts of matching tRNA genes (tRNA abundance) and promote accurate and/or efficient translation. Non-optimal codons, those least used in highly expressed genes, may also affect translation. In multicellular organisms, codon optimality may vary among tissues. At present, however, tissue specificity of codon use remains poorly understood. Here, we studied codon usage of genes highly transcribed in germ line (testis and ovary) and somatic tissues (gonadectomized males and females) of the beetle Tribolium castaneum. The results demonstrate that: (i) the majority of optimal codons were organism-wide, the same in all tissues, and had numerous matching tRNA gene copies (Opt-codon_↑tRNAs), consistent with translational selection; (ii) some optimal codons varied among tissues, suggesting tissue-specific tRNA populations; (iii) wobble tRNA were required for translation of certain optimal codons (Opt-codon_wobble), possibly allowing precise translation and/or protein folding; and (iv) remarkably, some non-optimal codons had abundant tRNA genes (Nonopt-codon_↑tRNAs), and genes using those codons were tightly linked to ribosomal and stress-response functions. Thus, Nonopt-codon_↑tRNAs codons may regulate translation of specific genes. Together, the evidence suggests that codon use and tRNA genes regulate multiple translational processes in T. castaneum.

The multi-faceted regulation of nuclear tRNA gene transcription

Transfer RNAs are among the most ancient molecules of life on earth. Beyond their crucial role in protein synthesis as carriers of amino acids, they are also important players in a plethora of other biological processes. Many debates in term of biogenesis, regulation and function persist around these fascinating non-coding RNAs. Our review focuses on the first step of their biogenesis in eukaryotes, i.e. their transcription from nuclear genes. Numerous and complementary ways have emerged during evolution to regulate transfer RNA gene transcription. Here, we will summarize the different actors implicated in this process: cis-elements, trans-factors, genomic contexts, epigenetic environments and finally three-dimensional organization of nuclear genomes. © 2019 IUBMB Life, 2019 © 2019 IUBMB Life, 71(8):1099-1108, 2019.

Cross-Kingdom Commonality of a Novel Insertion Signature of RTE-Related Short Retroposons

In multicellular organisms, such as vertebrates and flowering plants, horizontal transfer (HT) of genetic information is thought to be a rare event. However, recent findings unveiled unexpectedly frequent HT of RTE-clade LINEs. To elucidate the molecular footprints of the genomic integration machinery of RTE-related retroposons, the sequence patterns surrounding the insertion sites of plant Au-like SINE families were analyzed in the genomes of a wide variety of flowering plants. A novel and remarkable finding regarding target site duplications (TSDs) for SINEs was they start with thymine approximately one helical pitch (ten nucleotides) downstream of a thymine stretch. This TSD pattern was found in RTE-clade LINEs, which share the 3'-end sequence of these SINEs, in the genome of leguminous plants. These results demonstrably show that Au-like SINEs were mobilized by the enzymatic machinery of RTE-clade LINEs. Further, we discovered the same TSD pattern in animal SINEs from lizard and mammals, in which the RTE-clade LINEs sharing the 3'-end sequence with these animal SINEs showed a distinct TSD pattern. Moreover, a significant correlation was observed between the first nucleotide of TSDs and microsatellite-like sequences found at the 3'-ends of SINEs and LINEs. We propose that RTE-encoded protein could preferentially bind to a DNA region that contains a thymine stretch to cleave a phosphodiester bond downstream of the stretch. Further, determination of cleavage sites and/or efficiency of primer sites for reverse transcription may depend on microsatellite-like repeats in the RNA template. Such a unique mechanism may have enabled retroposons to successfully expand in frontier genomes after HT.

Refining the Ambush Hypothesis: Evidence That GC- and AT-Rich Bacteria Employ Different Frameshift Defence Strategies

Stop codons are frequently selected for beyond their regular termination function for error control. The “ambush hypothesis” proposes out-of-frame stop codons (OSCs) terminating frameshifted translations are selected for. Although early indirect evidence was partially supportive, recent evidence suggests OSC frequencies are not exceptional when considering underlying nucleotide content. However, prior null tests fail to control amino acid/codon usages or possible local mutational biases. We therefore return to the issue using bacterial genomes, considering several tests defining and testing against a null. We employ simulation approaches preserving amino acid order but shuffling synonymous codons or preserving codons while shuffling amino acid order. Additionally, we compare codon usage in amino acid pairs, where one codon can but the next, otherwise identical codon, cannot encode an OSC. OSC frequencies exceed expectations typically in AT-rich genomes, the +1 frame and for TGA/TAA but not TAG. With this complex evidence, simply rejecting or accepting the ambush hypothesis is not warranted. We propose a refined post hoc model, whereby AT-rich genomes have more accidental frameshifts, handled by RF2–RF3 complexes (associated with TGA/TAA) and are mostly +1 (or −2) slips. Supporting this, excesses positively correlate with in silico predicted frameshift probabilities. Thus, we propose a more viable framework, whereby genomes broadly adopt one of the two strategies to combat frameshifts: preventing frameshifting (GC-rich) or permitting frameshifts but minimizing impacts when most are caught early (AT-rich). Our refined framework holds promise yet some features, such as the bias of out-of-frame sense codons, remain unexplained.

Metabolic Labeling and Profiling of Transfer RNAs Using Macroarrays

Transfer RNAs (tRNA) are abundant short non-coding RNA species that are typically 76 to 90 nucleotides in length. tRNAs are directly responsible for protein synthesis by translating codons in mRNA into amino acid sequences. tRNAs were long considered as house-keeping molecules that lacked regulatory functions. However, a growing body of evidence indicates that cellular tRNA levels fluctuate in correspondence to varying conditions such as cell type, environment, and stress. The fluctuation of tRNA expression directly influences gene translation, favoring or repressing the expression of particular proteins. Ultimately comprehending the dynamic of protein synthesis requires the development of methods able to deliver high-quality tRNA profiles. The method that we present here is named SPOt, which stands for Streamlined Platform for Observing tRNA. SPOt consists of three steps starting with metabolic labeling of cell cultures with radioactive orthophosphate, followed by guanidinium thiocyanate-phenol-chloroform extraction of radioactive total RNAs and finally hybridization on in-house printed macroarrays. tRNA levels are estimated by quantifying the radioactivity intensities at each probe spot. In the protocol presented here we profile tRNAs in Mycobacterium smegmatis mc²155, a nonpathogenic bacterium often used as a model organism to study tuberculosis.

Analysis of Translation Elongation Dynamics in the Context of an Escherichia coli Cell

Understanding the mechanisms behind translation and its rate-limiting steps is crucial for both the development of drug targets and improvement of heterologous protein production with many biotechnological applications, such as in pharmaceutical and biofuel industries. Despite many advances in the knowledge of the ribosome structure and function, there is still much discussion around the determinants of translation elongation with experiments and computational studies pointing in different directions. Here, we use a stochastic framework to simulate the process of translation in the context of an Escherichia coli cell by gathering the available biochemical data into a ribosome kinetics description. Our results from the study of translation in E. coli at different growth rates contradict the increase of mean elongation rate with growth rate established in the literature. We show that both the level of tRNA competition and the type of cognate binding interaction contribute to the modulation of elongation rate, and that optimization of a heterologous transcript for faster elongation rate is achieved by combining the two. We derive an equation that can accurately predict codon elongation rates based on the abundances of free tRNA in the cell, and can be used to assist transcript design. Finally, we show that non-cognate tRNA-ribosome binding has an important weight in translation, and plays an active role in the modulation of mean elongation rate as shown by our amino-acid starvation/surplus studies.

SPOt: A novel and streamlined microarray platform for observing cellular tRNA levels

Recent studies have placed transfer RNA (tRNA), a housekeeping molecule, in the heart of fundamental cellular processes such as embryonic development and tumor progression. Such discoveries were contingent on the concomitant development of methods able to deliver high-quality tRNA profiles. The present study describes the proof of concept obtained in Escherichia coli (E. coli) for an original tRNA analysis platform named SPOt (Streamlined Platform for Observing tRNA). This approach comprises three steps. First, E. coli cultures are spiked with radioactive orthophosphate; second, labeled total RNAs are trizol-extracted; third, RNA samples are hybridized on in-house printed microarrays and spot signals, the proxy for tRNA levels, are quantified by phosphorimaging. Features such as reproducibility and specificity were assessed using several tRNA subpopulations. Dynamic range and sensitivity were evaluated by overexpressing specific tRNA species. SPOt does not require any amplification or post-extraction labeling and can be adapted to any organism. It is modular and easily streamlined with popular techniques such as polysome fractionation to profile tRNAs interacting with ribosomes and actively engaged in translation. The biological relevance of these data is discussed in regards to codon usage, tRNA gene copy number, and position on the genome.

Conserved 3' UTR stem-loop structure in L1 and Alu transposons in human genome: possible role in retrotransposition

BACKGROUND

In the process of retrotransposition LINEs use their own machinery for copying and inserting themselves into new genomic locations, while SINEs are parasitic and require the machinery of LINEs. The exact mechanism of how a LINE-encoded reverse transcriptase (RT) recognizes its own and SINE RNA remains unclear. However it was shown for the stringent-type LINEs that recognition of a stem-loop at the 3'UTR by RT is essential for retrotransposition. For the relaxed-type LINEs it is believed that the poly-A tail is a common recognition element between LINE and SINE RNA. However polyadenylation is a property of any messenger RNA, and how the LINE RT recognizes transposon and non-transposon RNAs remains an open question. It is likely that RNA secondary structures play an important role in RNA recognition by LINE encoded proteins.

RESULTS

Here we selected a set of L1 and Alu elements from the human genome and investigated their sequences for the presence of position-specific stem-loop structures. We found highly conserved stem-loop positions at the 3'UTR. Comparative structural analyses of a human L1 3'UTR stem-loop showed a similarity to 3'UTR stem-loops of the stringent-type LINEs, which were experimentally shown to be recognized by LINE RT. The consensus stem-loop structure consists of 5-7 bp loop, 8-10 bp stem with a bulge at a distance of 4-6 bp from the loop. The results show that a stem loop with a bulge exists at the 3'-end of Alu. We also found conserved stem-loop positions at 5'UTR and at the end of ORF2 and discuss their possible role.

CONCLUSIONS

Here we presented an evidence for the presence of a highly conserved 3'UTR stem-loop structure in L1 and Alu retrotransposons in the human genome. Both stem-loops show structural similarity to the stem-loops of the stringent-type LINEs experimentally confirmed as essential for retrotransposition. Here we hypothesize that both L1 and Alu RNA are recognized by L1 RT via the 3'-end RNA stem-loop structure. Other conserved stem-loop positions in L1 suggest their possible functions in protein-RNA interactions but to date no experimental evidence has been reported.

Evolution of codon usage in the smallest photosynthetic eukaryotes and their giant viruses

Prasinoviruses are among the largest viruses (>200 kb) and encode several hundreds of protein coding genes, including most genes of the DNA replication machinery and several genes involved in transcription and translation, as well as transfer RNAs (tRNAs). They can infect and lyse small eukaryotic planktonic marine green algae, thereby affecting global algal population dynamics. Here, we investigate the causes of codon usage bias (CUB) in one prasinovirus, OtV5, and its host Ostreococcus tauri, during a viral infection using microarray expression data. We show that 1) CUB in the host and in the viral genes increases with expression levels and 2) optimal codons use those tRNAs encoded by the most abundant host tRNA genes, supporting the notion of translational optimization by natural selection. We find evidence that viral tRNA genes complement the host tRNA pool for those viral amino acids whose host tRNAs are in short supply. We further discuss the coevolution of CUB in hosts and prasinoviruses by comparing optimal codons in three evolutionary diverged host–virus-specific pairs whose complete genome sequences are known.

RNA-Mediated Gene Duplication and Retroposons: Retrogenes, LINEs, SINEs, and Sequence Specificity

A substantial number of “retrogenes” that are derived from the mRNA of various intron-containing genes have been reported. A class of mammalian retroposons, long interspersed element-1 (LINE1, L1), has been shown to be involved in the reverse transcription of retrogenes (or processed pseudogenes) and non-autonomous short interspersed elements (SINEs). The 3′-end sequences of various SINEs originated from a corresponding LINE. As the 3′-untranslated regions of several LINEs are essential for retroposition, these LINEs presumably require “stringent” recognition of the 3′-end sequence of the RNA template. However, the 3′-ends of mammalian L1s do not exhibit any similarity to SINEs, except for the presence of 3′-poly(A) repeats. Since the 3′-poly(A) repeats of L1 and Alu SINE are critical for their retroposition, L1 probably recognizes the poly(A) repeats, thereby mobilizing not only Alu SINE but also cytosolic mRNA. Many flowering plants only harbor L1-clade LINEs and a significant number of SINEs with poly(A) repeats, but no homology to the LINEs. Moreover, processed pseudogenes have also been found in flowering plants. I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized a specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

PlantRNA, a database for tRNAs of photosynthetic eukaryotes

PlantRNA database (http://plantrna.ibmp.cnrs.fr/) compiles transfer RNA (tRNA) gene sequences retrieved from fully annotated plant nuclear, plastidial and mitochondrial genomes. The set of annotated tRNA gene sequences has been manually curated for maximum quality and confidence. The novelty of this database resides in the inclusion of biological information relevant to the function of all the tRNAs entered in the library. This includes 5'- and 3'-flanking sequences, A and B box sequences, region of transcription initiation and poly(T) transcription termination stretches, tRNA intron sequences, aminoacyl-tRNA synthetases and enzymes responsible for tRNA maturation and modification. Finally, data on mitochondrial import of nuclear-encoded tRNAs as well as the bibliome for the respective tRNAs and tRNA-binding proteins are also included. The current annotation concerns complete genomes from 11 organisms: five flowering plants (Arabidopsis thaliana, Oryza sativa, Populus trichocarpa, Medicago truncatula and Brachypodium distachyon), a moss (Physcomitrella patens), two green algae (Chlamydomonas reinhardtii and Ostreococcus tauri), one glaucophyte (Cyanophora paradoxa), one brown alga (Ectocarpus siliculosus) and a pennate diatom (Phaeodactylum tricornutum). The database will be regularly updated and implemented with new plant genome annotations so as to provide extensive information on tRNA biology to the research community.

Co-evolution of mitochondrial tRNA import and codon usage determines translational efficiency in the green alga Chlamydomonas

Mitochondria from diverse phyla, including protozoa, fungi, higher plants, and humans, import tRNAs from the cytosol in order to ensure proper mitochondrial translation. Despite the broad occurrence of this process, our understanding of tRNA import mechanisms is fragmentary, and crucial questions about their regulation remain unanswered. In the unicellular green alga Chlamydomonas, a precise correlation was found between the mitochondrial codon usage and the nature and amount of imported tRNAs. This led to the hypothesis that tRNA import might be a dynamic process able to adapt to the mitochondrial genome content. By manipulating the Chlamydomonas mitochondrial genome, we introduced point mutations in order to modify its codon usage. We find that the codon usage modification results in reduced levels of mitochondrial translation as well as in subsequent decreased levels and activities of respiratory complexes. These effects are linked to the consequential limitations of the pool of tRNAs in mitochondria. This indicates that tRNA mitochondrial import cannot be rapidly regulated in response to a novel genetic context and thus does not appear to be a dynamic process. It rather suggests that the steady-state levels of imported tRNAs in mitochondria result from a co-evolutive adaptation between the tRNA import mechanism and the requirements of the mitochondrial translation machinery.

Mitochondria are endosymbiotic organelles involved in diverse fundamental cellular processes. They contain their own genome that encodes a few but essential proteins (e.g. proteins of the respiratory chain complexes). Their synthesis requires functional mitochondrial translational machinery that necessitates a full set of transfer RNAs (tRNAs). As mitochondrial genomes of various organisms do not code for the complete set of tRNA genes, nucleus-encoded tRNAs have to be imported into mitochondria to compensate. Mitochondrial import of tRNAs is highly specific and tailored to the mitochondrial needs. Because transformation of the mitochondrial genome is possible in Chlamydomonas, we used this green alga as model to know if a fine regulation of the tRNA import process is possible so that the tRNA population can rapidly adapt to codon usage changes in mitochondria. Here we provide evidence that the regulation of tRNA mitochondrial import process is not dynamic but is rather the result of a co-evolutive process between the import and the mitochondrial codon bias in order to optimize the mitochondrial translation efficiency.

Parallel relaxation of stringent RNA recognition in plant and mammalian L1 retrotransposons

L1 elements are mammalian non-long terminal repeat retrotransposons, or long interspersed elements (LINEs), that significantly influence the dynamics and fluidity of the genome. A series of observations suggest that plant L1-clade LINEs, just as mammalian L1s, mobilize both short interspersed elements (SINEs) and certain messenger RNA by recognizing the 3'-poly(A) tail of RNA. However, one L1 lineage in monocots was shown to possess a conserved 3'-end sequence with a solid RNA structure also observed in maize and sorghum SINEs. This strongly suggests that plant LINEs require a particular 3'-end sequence during initiation of reverse transcription. As one L1-clade LINE was also found to share the 3'-end sequence with a SINE in a green algal genome, I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized the specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

Mutation bias is the driving force of codon usage in the Gallus gallus genome

Synonymous codons are used with different frequencies both among species and among genes within the same genome and are controlled by neutral processes (such as mutation and drift) as well as by selection. Up to now, a systematic examination of the codon usage for the chicken genome has not been performed. Here, we carried out a whole genome analysis of the chicken genome by the use of the relative synonymous codon usage (RSCU) method and identified 11 putative optimal codons, all of them ending with uracil (U), which is significantly departing from the pattern observed in other eukaryotes. Optimal codons in the chicken genome are most likely the ones corresponding to highly expressed transfer RNA (tRNAs) or tRNA gene copy numbers in the cell. Codon bias, measured as the frequency of optimal codons (Fop), is negatively correlated with the G + C content, recombination rate, but positively correlated with gene expression, protein length, gene length and intron length. The positive correlation between codon bias and protein, gene and intron length is quite different from other multi-cellular organism, as this trend has been only found in unicellular organisms. Our data displayed that regional G + C content explains a large proportion of the variance of codon bias in chicken. Stepwise selection model analyses indicate that G + C content of coding sequence is the most important factor for codon bias. It appears that variation in the G + C content of CDSs accounts for over 60% of the variation of codon bias. This study suggests that both mutation bias and selection contribute to codon bias. However, mutation bias is the driving force of the codon usage in the Gallus gallus genome. Our data also provide evidence that the negative correlation between codon bias and recombination rates in G. gallus is determined mostly by recombination-dependent mutational patterns.

RNA polymerase III transcription control elements: themes and variations

Eukaryotic genomes are punctuated by a multitude of tiny genetic elements, that share the property of being recognized and transcribed by the RNA polymerase (Pol) III machinery to produce a variety of small, abundant non-protein-coding (nc) RNAs (tRNAs, 5S rRNA, U6 snRNA and many others). The highly selective, efficient and localized action of Pol III at its minute genomic targets is made possible by a handful of cis-acting regulatory elements, located within the transcribed region (where they are bound by the multisubunit assembly factor TFIIIC) and/or upstream of the transcription start site. Most of them participate directly or indirectly in the ultimate recruitment of TFIIIB, a key multiprotein initiation factor able to direct, once assembled, multiple transcription cycles by Pol III. But the peculiar efficiency and selectivity of Pol III transcription also depends on its ability to recognize very simple and precisely positioned termination signals. Studies in the last few years have significantly expanded the set of known Pol III-associated loci in genomes and, concomitantly, have revealed unexpected features of Pol III cis-regulatory elements in terms of variety, function, genomic location and potential contribution to transcriptome complexity. Here we review, in a historical perspective, well established and newly acquired knowledge about Pol III transcription control elements, with the aim of providing a useful reference for future studies of the Pol III system, which we anticipate will be numerous and intriguing for years to come.

Coadaptation of isoacceptor tRNA genes and codon usage bias for translation efficiency in Aedes aegypti and Anopheles gambiae

The transfer RNAs (tRNAs) are essential components of translational machinery. We determined that tRNA isoacceptors (tRNAs with different anticodons but incorporating the same amino acid in protein synthesis) show differential copy number abundance, genomic distribution patterns and sequence evolution between Aedes aegypti and Anopheles gambiae mosquitoes. The tRNA-Ala genes are present in unusually high copy number in the Ae. aegypti genome but not in An. gambiae. Many of the tRNA-Ala genes of Ae. aegypti are flanked by a highly conserved sequence that is not observed in An. gambiae. The relative abundance of tRNA isoacceptor genes is correlated with preferred (or optimal) and nonpreferred (or rare) codons for ∼2-4% of the predicted protein coding genes in both species. The majority (∼74-85%) of these genes are related to pathways involved with translation, energy metabolism and carbohydrate metabolism. Our results suggest that these genes and the related pathways may be under translational selection in these mosquitoes.

A global picture of tRNA genes in plant genomes

Although transfer RNA (tRNA) has a fundamental role in cell life, little is known about tRNA gene organization and expression on a genome-wide scale in eukaryotes, particularly plants. Here, we analyse the content and distribution of tRNA genes in five flowering plants and one green alga. The tRNA gene content is homogenous in plants, and is mostly correlated with genome size. The number of tRNA pseudogenes and organellar-like tRNA genes present in nuclear genomes varies greatly from one plant species to another. These pseudogenes or organellar-like genes appear to be generated or inserted randomly during evolution. Interestingly, we identified a new family of tRNA-related short interspersed nuclear elements (SINEs) in the Populus trichocarpa nuclear genome. In higher plants, intron-containing tRNA genes are rare, and correspond to genes coding for tRNA(Tyr) and tRNA(Mete) . By contrast, in green algae, more than half of the tRNA genes contain an intron. This suggests divergent means of intron acquisition and the splicing process between green algae and land plants. Numerous tRNAs are co-transcribed in Chlamydomonas, but they are mostly transcribed as a single unit in flowering plants. The only exceptions are tRNA(Gly) -snoRNA and tRNA(Mete) -snoRNA cotranscripts in dicots and monocots, respectively. The internal or external motifs required for efficient transcription of tRNA genes by RNA polymerase III are well conserved among angiosperms. A brief analysis of the mitochondrial and plastidial tRNA gene populations is also provided.

Effect of correlated tRNA abundances on translation errors and evolution of codon usage bias

Despite the fact that tRNA abundances are thought to play a major role in determining translation error rates, their distribution across the genetic code and the resulting implications have received little attention. In general, studies of codon usage bias (CUB) assume that codons with higher tRNA abundance have lower missense error rates. Using a model of protein translation based on tRNA competition and intra-ribosomal kinetics, we show that this assumption can be violated when tRNA abundances are positively correlated across the genetic code. Examining the distribution of tRNA abundances across 73 bacterial genomes from 20 different genera, we find a consistent positive correlation between tRNA abundances across the genetic code. This work challenges one of the fundamental assumptions made in over 30 years of research on CUB that codons with higher tRNA abundances have lower missense error rates and that missense errors are the primary selective force responsible for CUB.

Codon usage bias (CUB) is a ubiquitous and important phenomenon. CUB is thought to be driven primarily due to selection against missense errors. For over 30 years, the standard model of translation errors has implicitly assumed that the relationship between translation errors and tRNA abundances are inversely related. This is based on an implicit and unstated assumption that the distribution of tRNA abundances across the genetic code are uncorrelated. Examining these abundance distributions across 73 bacterial genomes from 20 different genera, we find a consistent positive correlation between tRNA abundances across the genetic code. We further show that codons with higher tRNA abundances are not always “optimal” with respect to reducing the missense error rate and hence cannot explain the observed patterns of CUB.

Computational analysis of tRNA identity

I review recent developments in computational analysis of tRNA identity. I suggest that the tRNA-protein interaction network is hierarchically organized, and coevolutionarily flexible. Its functional specificity of recognition and discrimination persists despite generic structural constraints and perturbative evolutionary forces. This flexibility comes from its arbitrary nature as a self-recognizing shape code. A revisualization of predicted Proteobacterial tRNA identity highlights open research problems. tRNA identity elements and their coevolution with proteins must be mapped structurally over the Tree of Life. These traits can also resolve deep roots in the Tree. I show that histidylation identity elements phylogenetically reposition Pelagibacter ubique within alpha-Proteobacteria.

Steady-state levels of imported tRNAs in Chlamydomonas mitochondria are correlated with both cytosolic and mitochondrial codon usages

The mitochondrial genome of Chlamydomonas reinhardtii only encodes three expressed tRNA genes, thus most mitochondrial tRNAs are likely imported. The sharing of tRNAs between chloroplasts and mitochondria has been speculated in this organism. We first demonstrate that no plastidial tRNA is present in mitochondria and that the mitochondrial translation mainly relies on the import of nucleus-encoded tRNA species. Then, using northern analysis, we show that the extent of mitochondrial localization for the 49 tRNA isoacceptor families encoded by the C. reinhardtii nuclear genome is highly variable. Until now the reasons for such variability were unknown. By comparing cytosolic and mitochondrial codon usage with the sub-cellular distribution of tRNAs, we provide unprecedented evidence that the steady-state level of a mitochondrial tRNA is linked not only to the frequency of the cognate codon in mitochondria but also to its frequency in the cytosol, then allowing optimal mitochondrial translation.