Novel anticodon composition of transfer RNAs in Micrococcus luteus, a bacterium with a high genomic G + C content. Correlation with codon usage

tRNAVal allows four-way decoding with unmodified uridine at the wobble position in Lactobacillus casei

Modifications at the wobble position (position 34) of tRNA facilitate interactions that enable or stabilize non-Watson-Crick base pairs. In bacterial tRNA, 5-hydroxyuridine (ho5U) derivatives xo5U [x: methyl (mo5U), carboxymethyl (cmo5U), and methoxycarbonylmethyl (mcmo5U)] present at the wobble positions of tRNAs are responsible for the recognition of NYN codon families. These modifications of U34 allow base-pairing not only with A and G but also with U, and in some cases, C. mo5U was originally found in Gram-positive bacteria, and cmo5U and mcmo5U were found in Gram-negative bacteria. tRNAs of Mycoplasma species, mitochondria, and chloroplasts adopt four-way decoding in which unmodified U34 recognizes codons ending in A, G, C, and U. Lactobacillus casei, Gram-positive bacteria, and lactic acid bacteria lack the modification enzyme genes for xo5U biosynthesis. Nevertheless, L. casei has only one type of tRNAVal with the anticodon UAC [tRNAVal(UAC)]. However, the genome of L. casei encodes an undetermined tRNA (tRNAUnd) gene, and the sequence corresponding to the anticodon region is GAC. Here, we confirm that U34 in L. casei tRNAVal is unmodified and that there is no tRNAUnd expression in the cells. In addition, in vitro transcribed tRNAUnd was not aminoacylated by L. casei valyl-tRNA synthetase, suggesting that tRNAUnd is not able to accept valine, even if expressed in cells. Correspondingly, native tRNAVal(UAC) with unmodified U34 bound to all four valine codons in the ribosome A site. This suggests that L. casei tRNAVal decodes all valine codons by four-way decoding, similarly to tRNAs from Mycoplasma species, mitochondria, and chloroplasts.

Comparative analysis of codon usage between Gossypium hirsutum and G. barbadense mitochondrial genomes

Gossypium hirsutum and G. barbadense mitochondrial genomes were analyzed to understand the factors shaping codon usage. While most analyses of codon usage suggest minimal to no bias, nucleotide composition, specifically GC content, was significantly correlated with codon usage. In general, both mitochondrial genomes favor codons that end in A or U, with a secondary preference for pyrimidine rich codons. These observations are similar to previous reports of codon usage in cotton nuclear genomes, possibly suggestive of a general bias spanning genomic compartment. Although evidence for codon usage bias is weak for most genes, we identified six genes (i.e. atp8, atp9, sdh3, sdh4, mttB and rpl2) with significant nonrandom codon usage. In general, we find multiple factors that influence cotton mitochondrial genome codon usage, which may include selection in a subset of genes.

The Role of Orthogonality in Genetic Code Expansion

The genetic code defines how information in the genome is translated into protein. Aside from a handful of isolated exceptions, this code is universal. Researchers have developed techniques to artificially expand the genetic code, repurposing codons and translational machinery to incorporate nonstandard amino acids (nsAAs) into proteins. A key challenge for robust genetic code expansion is orthogonality; the engineered machinery used to introduce nsAAs into proteins must co-exist with native translation and gene expression without cross-reactivity or pleiotropy. The issue of orthogonality manifests at several levels, including those of codons, ribosomes, aminoacyl-tRNA synthetases, tRNAs, and elongation factors. In this concept paper, we describe advances in genome recoding, translational engineering and associated challenges rooted in establishing orthogonality needed to expand the genetic code.

An integrated, structure- and energy-based view of the genetic code

The principles of mRNA decoding are conserved among all extant life forms. We present an integrative view of all the interaction networks between mRNA, tRNA and rRNA: the intrinsic stability of codon-anticodon duplex, the conformation of the anticodon hairpin, the presence of modified nucleotides, the occurrence of non-Watson-Crick pairs in the codon-anticodon helix and the interactions with bases of rRNA at the A-site decoding site. We derive a more information-rich, alternative representation of the genetic code, that is circular with an unsymmetrical distribution of codons leading to a clear segregation between GC-rich 4-codon boxes and AU-rich 2:2-codon and 3:1-codon boxes. All tRNA sequence variations can be visualized, within an internal structural and energy framework, for each organism, and each anticodon of the sense codons. The multiplicity and complexity of nucleotide modifications at positions 34 and 37 of the anticodon loop segregate meaningfully, and correlate well with the necessity to stabilize AU-rich codon-anticodon pairs and to avoid miscoding in split codon boxes. The evolution and expansion of the genetic code is viewed as being originally based on GC content with progressive introduction of A/U together with tRNA modifications. The representation we present should help the engineering of the genetic code to include non-natural amino acids.

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.

Genetic code evolution reveals the neutral emergence of mutational robustness, and information as an evolutionary constraint

The standard genetic code (SGC) is central to molecular biology and its origin and evolution is a fundamental problem in evolutionary biology, the elucidation of which promises to reveal much about the origins of life. In addition, we propose that study of its origin can also reveal some fundamental and generalizable insights into mechanisms of molecular evolution, utilizing concepts from complexity theory. The first is that beneficial traits may arise by non-adaptive processes, via a process of "neutral emergence". The structure of the SGC is optimized for the property of error minimization, which reduces the deleterious impact of point mutations. Via simulation, it can be shown that genetic codes with error minimization superior to the SGC can emerge in a neutral fashion simply by a process of genetic code expansion via tRNA and aminoacyl-tRNA synthetase duplication, whereby similar amino acids are added to codons related to that of the parent amino acid. This process of neutral emergence has implications beyond that of the genetic code, as it suggests that not all beneficial traits have arisen by the direct action of natural selection; we term these "pseudaptations", and discuss a range of potential examples. Secondly, consideration of genetic code deviations (codon reassignments) reveals that these are mostly associated with a reduction in proteome size. This code malleability implies the existence of a proteomic constraint on the genetic code, proportional to the size of the proteome (P), and that its reduction in size leads to an "unfreezing" of the codon - amino acid mapping that defines the genetic code, consistent with Crick's Frozen Accident theory. The concept of a proteomic constraint may be extended to propose a general informational constraint on genetic fidelity, which may be used to explain variously, differences in mutation rates in genomes with differing proteome sizes, differences in DNA repair capacity and genome GC content between organisms, a selective pressure in the evolution of sexual reproduction, and differences in translational fidelity. Lastly, the utility of the concept of an informational constraint to other diverse fields of research is explored.

The effect of multiple evolutionary selections on synonymous codon usage of genes in the Mycoplasma bovis genome

Mycoplasma bovis is a major pathogen causing arthritis, respiratory disease and mastitis in cattle. A better understanding of its genetic features and evolution might represent evidences of surviving host environments. In this study, multiple factors influencing synonymous codon usage patterns in M. bovis (three strains’ genomes) were analyzed. The overall nucleotide content of genes in the M. bovis genome is AT-rich. Although the G and C contents at the third codon position of genes in the leading strand differ from those in the lagging strand (p<0.05), the 59 synonymous codon usage patterns of genes in the leading strand are highly similar to those in the lagging strand. The over-represented codons and the under-represented codons were identified. A comparison of the synonymous codon usage pattern of M. bovis and cattle (susceptible host) indicated the independent formation of synonymous codon usage of M. bovis. Principal component analysis revealed that (i) strand-specific mutational bias fails to affect the synonymous codon usage pattern in the leading and lagging strands, (ii) mutation pressure from nucleotide content plays a role in shaping the overall codon usage, and (iii) the major trend of synonymous codon usage has a significant correlation with the gene expression level that is estimated by the codon adaptation index. The plot of the effective number of codons against the G+C content at the third codon position also reveals that mutation pressure undoubtedly contributes to the synonymous codon usage pattern of M. bovis. Additionally, the formation of the overall codon usage is determined by certain evolutionary selections for gene function classification (30S protein, 50S protein, transposase, membrane protein, and lipoprotein) and translation elongation region of genes in M. bovis. The information could be helpful in further investigations of evolutionary mechanisms of the Mycoplasma family and heterologous expression of its functionally important proteins.

Molecular reconstruction of a fungal genetic code alteration

Fungi of the CTG clade translate the Leu CUG codon as Ser. This genetic code alteration is the only eukaryotic sense-to-sense codon reassignment known to date, is mediated by an ambiguous serine tRNA (tRNACAG(Ser)), exposes unanticipated flexibility of the genetic code and raises major questions about its selection and fixation in this fungal lineage. In particular, the origin of the tRNACAG(Ser) and the evolutionary mechanism of CUG reassignment from Leu to Ser remain poorly understood. In this study, we have traced the origin of the tDNACAG(Ser) gene and studied critical mutations in the tRNACAG(Ser) anticodon-loop that modulated CUG reassignment. Our data show that the tRNACAG(Ser) emerged from insertion of an adenosine in the middle position of the 5'-CGA-3'anticodon of a tRNACGA(Ser) ancestor, producing the 5'-CAG-3' anticodon of the tRNACAG(Ser), without altering its aminoacylation properties. This mutation initiated CUG reassignment while two additional mutations in the anticodon-loop resolved a structural conflict produced by incorporation of the Leu 5'-CAG-3'anticodon in the anticodon-arm of a tRNA(Ser). Expression of the mutant tRNACAG(Ser) in yeast showed that it cannot be expressed at physiological levels and we postulate that such downregulation was essential to maintain Ser misincorporation at sub-lethal levels during the initial stages of CUG reassignment. We demonstrate here that such low level CUG ambiguity is advantageous in specific ecological niches and we propose that misreading tRNAs are targeted for degradation by an unidentified tRNA quality control pathway.

Altered tRNA characteristics and 3' maturation in bacterial symbionts with reduced genomes

Translational efficiency is controlled by tRNAs and other genome-encoded mechanisms. In organelles, translational processes are dramatically altered because of genome shrinkage and horizontal acquisition of gene products. The influence of genome reduction on translation in endosymbionts is largely unknown. Here, we investigate whether divergent lineages of Buchnera aphidicola, the reduced-genome bacterial endosymbiont of aphids, possess altered translational features compared with their free-living relative, Escherichia coli. Our RNAseq data support the hypothesis that translation is less optimal in Buchnera than in E. coli. We observed a specific, convergent, pattern of tRNA loss in Buchnera and other endosymbionts that have undergone genome shrinkage. Furthermore, many modified nucleoside pathways that are important for E. coli translation are lost in Buchnera. Additionally, Buchnera’s A + T compositional bias has resulted in reduced tRNA thermostability, and may have altered aminoacyl-tRNA synthetase recognition sites. Buchnera tRNA genes are shorter than those of E. coli, as the majority no longer has a genome-encoded 3' CCA; however, all the expressed, shortened tRNAs undergo 3′ CCA maturation. Moreover, expression of tRNA isoacceptors was not correlated with the usage of corresponding codons. Overall, our data suggest that endosymbiont genome evolution alters tRNA characteristics that are known to influence translational efficiency in their free-living relative.

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.

Evolution of synonymous codon usage in Neurospora tetrasperma and Neurospora discreta

Neurospora comprises a primary model system for the study of fungal genetics and biology. In spite of this, little is known about genome evolution in Neurospora. For example, the evolution of synonymous codon usage is largely unknown in this genus. In the present investigation, we conducted a comprehensive analysis of synonymous codon usage and its relationship to gene expression and gene length (GL) in Neurospora tetrasperma and Neurospora discreta. For our analysis, we examined codon usage among 2,079 genes per organism and assessed gene expression using large-scale expressed sequenced tag (EST) data sets (279,323 and 453,559 ESTs for N. tetrasperma and N. discreta, respectively). Data on relative synonymous codon usage revealed 24 codons (and two putative codons) that are more frequently used in genes with high than with low expression and thus were defined as optimal codons. Although codon-usage bias was highly correlated with gene expression, it was independent of selectively neutral base composition (introns); thus demonstrating that translational selection drives synonymous codon usage in these genomes. We also report that GL (coding sequences [CDS]) was inversely associated with optimal codon usage at each gene expression level, with highly expressed short genes having the greatest frequency of optimal codons. Optimal codon frequency was moderately higher in N. tetrasperma than in N. discreta, which might be due to variation in selective pressures and/or mating systems.

Genome sequence of the Fleming strain of Micrococcus luteus, a simple free-living actinobacterium

Micrococcus luteus (NCTC2665, "Fleming strain") has one of the smallest genomes of free-living actinobacteria sequenced to date, comprising a single circular chromosome of 2,501,097 bp (G+C content, 73%) predicted to encode 2,403 proteins. The genome shows extensive synteny with that of the closely related organism, Kocuria rhizophila, from which it was taxonomically separated relatively recently. Despite its small size, the genome harbors 73 insertion sequence (IS) elements, almost all of which are closely related to elements found in other actinobacteria. An IS element is inserted into the rrs gene of one of only two rrn operons found in M. luteus. The genome encodes only four sigma factors and 14 response regulators, a finding indicative of adaptation to a rather strict ecological niche (mammalian skin). The high sensitivity of M. luteus to beta-lactam antibiotics may result from the presence of a reduced set of penicillin-binding proteins and the absence of a wblC gene, which plays an important role in the antibiotic resistance in other actinobacteria. Consistent with the restricted range of compounds it can use as a sole source of carbon for energy and growth, M. luteus has a minimal complement of genes concerned with carbohydrate transport and metabolism and its inability to utilize glucose as a sole carbon source may be due to the apparent absence of a gene encoding glucokinase. Uniquely among characterized bacteria, M. luteus appears to be able to metabolize glycogen only via trehalose and to make trehalose only via glycogen. It has very few genes associated with secondary metabolism. In contrast to most other actinobacteria, M. luteus encodes only one resuscitation-promoting factor (Rpf) required for emergence from dormancy, and its complement of other dormancy-related proteins is also much reduced. M. luteus is capable of long-chain alkene biosynthesis, which is of interest for advanced biofuel production; a three-gene cluster essential for this metabolism has been identified in the genome.

Evolving genetic code

In 1985, we reported that a bacterium, Mycoplasma capricolum, used a deviant genetic code, namely UGA, a "universal" stop codon, was read as tryptophan. This finding, together with the deviant nuclear genetic codes in not a few organisms and a number of mitochondria, shows that the genetic code is not universal, and is in a state of evolution. To account for the changes in codon meanings, we proposed the codon capture theory stating that all the code changes are non-disruptive without accompanied changes of amino acid sequences of proteins. Supporting evidence for the theory is presented in this review. A possible evolutionary process from the ancient to the present-day genetic code is also discussed.

Gender-specific selection on codon usage in plant genomes

Background

Currently, there is little data available regarding the role of gender-specific gene expression on synonymous codon usage (translational selection) in most organisms, and particularly plants. Using gender-specific EST libraries (with > 4000 ESTs) from Zea mays and Triticum aestivum, we assessed whether gender-specific gene expression per se and gender-specific gene expression level are associated with selection on codon usage.

Results

We found clear evidence of a greater bias in codon usage for genes expressed in female than in male organs and gametes, based on the variation in GC content at third codon positions and the frequency of species-preferred codons. This finding holds true for both highly and for lowly expressed genes. In addition, we found that highly expressed genes have greater codon bias than lowly expressed genes for both female- and male-specific genes. Moreover, in both species, genes with female-specific expression show a greater usage of species-specific preferred codons for each of the 18 amino acids having synonymous codons. A supplemental analysis of Brassica napus suggests that bias in codon usage could also be higher in genes expressed in male gametophytic tissues than in heterogeneous (flower) tissues.

Conclusion

This study reports gender-specific bias in codon usage in plants. The findings reported here, based on the analysis of 1 497 876 codons, are not caused either by differences in the biological functions of the genes or by differences in protein lengths, nor are they likely attributable to mutational bias. The data are best explained by gender-specific translational selection. Plausible explanations for these findings and the relevance to these and other organisms are discussed.

Differences in codon bias cannot explain differences in translational power among microbes

Background

Translational power is the cellular rate of protein synthesis normalized to the biomass invested in translational machinery. Published data suggest a previously unrecognized pattern: translational power is higher among rapidly growing microbes, and lower among slowly growing microbes. One factor known to affect translational power is biased use of synonymous codons. The correlation within an organism between expression level and degree of codon bias among genes of Escherichia coli and other bacteria capable of rapid growth is commonly attributed to selection for high translational power. Conversely, the absence of such a correlation in some slowly growing microbes has been interpreted as the absence of selection for translational power. Because codon bias caused by translational selection varies between rapidly growing and slowly growing microbes, we investigated whether observed differences in translational power among microbes could be explained entirely by differences in the degree of codon bias. Although the data are not available to estimate the effect of codon bias in other species, we developed an empirically-based mathematical model to compare the translation rate of E. coli to the translation rate of a hypothetical strain which differs from E. coli only by lacking codon bias.

Results

Our reanalysis of data from the scientific literature suggests that translational power can differ by a factor of 5 or more between E. coli and slowly growing microbial species. Using empirical codon-specific in vivo translation rates for 29 codons, and several scenarios for extrapolating from these data to estimates over all codons, we find that codon bias cannot account for more than a doubling of the translation rate in E. coli, even with unrealistic simplifying assumptions that exaggerate the effect of codon bias. With more realistic assumptions, our best estimate is that codon bias accelerates translation in E. coli by no more than 60% in comparison to microbes with very little codon bias.

Conclusions

While codon bias confers a substantial benefit of faster translation and hence greater translational power, the magnitude of this effect is insufficient to explain observed differences in translational power among bacterial and archaeal species, particularly the differences between slowly growing and rapidly growing species. Hence, large differences in translational power suggest that the translational apparatus itself differs among microbes in ways that influence translational performance.

Codon usage bias from tRNA's point of view: redundancy, specialization, and efficient decoding for translation optimization

The selection-mutation-drift theory of codon usage plays a major role in the theory of molecular evolution by explaining the co-evolution of codon usage bias and tRNA content in the framework of translation optimization. Because most studies have focused only on codon usage, we analyzed the tRNA gene pool of 102 bacterial species. We show that as minimal generation times get shorter, the genomes contain more tRNA genes, but fewer anticodon species. Surprisingly, despite the wide G+C variation of bacterial genomes these anticodons are the same in most genomes. This suggests an optimization of the translation machinery to use a small subset of optimal codons and anticodons in fast-growing bacteria and in highly expressed genes. As a result, the overrepresented codons in highly expressed genes tend to be the same in very different genomes to match the same most-frequent anticodons. This is particularly important in fast-growing bacteria, which have higher codon usage bias in these genes. Three models were tested to understand the choice of codons recognized by the same anticodons, all providing significant fit, but under different classes of genes and genomes. Thus, co-evolution of tRNA gene composition and codon usage bias in genomes seen from tRNA's point of view agrees with the selection-mutation-drift theory. However, it suggests a much more universal trend in the evolution of anticodon and codon choice than previously thought. It also provides new evidence that a selective force for the optimization of the translation machinery is the maximization of growth.

Driving change: the evolution of alternative genetic codes

Pioneering studies in the 1960s that elucidated the genetic code suggested that all extant forms of life use the same genetic code. This early presumption has subsequently been challenged by the discovery of deviations of the universal genetic code in prokaryotes, eukaryotic nuclear genomes and mitochondrial genomes. These studies have revealed that the genetic code is still evolving despite strong negative forces working against the fixation of mutations that result in codon reassignment. Recent data from in vitro, in vivo and in silico comparative genomics studies are revealing significant, previously overlooked links between modified nucleosides in tRNAs, genetic code ambiguity, genome base composition, codon usage and codon reassignment.

Unnatural RNA display libraries

Combinatorial peptide and protein libraries have now been developed to accommodate unnatural amino acids in a genetically encoded format via in vitro nonsense and sense suppression. General translation features and specific regioselective and stereoselective properties of the ribosome endow these libraries with a broad chemical diversity. Alternatively, amino acid residues can be chemically derivatized post-translationally to add preferred functionality to the encoded peptide. All of these efforts are advancing combinatorial peptide and protein libraries for enhanced ligands against biological targets of interest.

Patterns of codon usage bias in three dicot and four monocot plant species

Codon usage in nuclear genes of four monocot and three dicot species was analyzed to find general patterns in codon choice of plant species. Codon bias was correlated with GC content at the third codon position. GC contents were higher in monocot species than in dicot species at all codon positions. The high GC contents of monocot species might be the result of relatively strong mutational bias that occurred in the lineage of the Poaceae species. In both dicot and monocot species, the effective number of codons (ENCs) for most genes was similar to that for the expected ENCs based on the GC content at the third codon positions. G and C ending codons were detected as the "preferred" codons in monocot species, as in Drosophila. Also, many "preferred" codons are the same in dicot species. Pyrimidine (C and T) is used more frequently than purine (G and A) in four-fold degenerate codon groups.

Investigating the relationship between genome structure, composition, and ecology in prokaryotes

Our thesis is that the DNA composition and structure of genomes are selected in part by mutation bias (GC pressure) and in part by ecology. To illustrate this point, we compare and contrast the oligonucleotide composition and the mosaic structure in 36 complete genomes and in 27 long genomic sequences from archaea and eubacteria. We report the following findings (1) High-GC-content genomes show a large underrepresentation of short distances between G(n) and C(n) homopolymers with respect to distances between A(n) and T(n) homopolymers; we discuss selection versus mutation bias hypotheses. (2) The oligonucleotide compositions of the genomes of Neisseria (meningitidis and gonorrhoea), Helicobacter pylori and Rhodobacter capsulatus are more biased than the other sequenced genomes. (3) The genomes of free-living species or nonchronic pathogens show more mosaic-like structure than genomes of chronic pathogens or intracellular symbionts. (4) Genome mosaicity of intracellular parasites has a maximum corresponding to the average gene length; in the genomes of free-living and nonchronic pathogens the maximum occurs at larger length scales. This suggests that free-living species can incorporate large pieces of DNA from the environment, whereas for intracellular parasites there are recombination events between homologous genes. We discuss the consequences in terms of evolution of genome size. (5) Intracellular symbionts and obligate pathogens show small, but not zero, amount of chromosome mosaicity, suggesting that recombination events occur in these species.

Studies of codon usage and tRNA genes of 18 unicellular organisms and quantification of Bacillus subtilis tRNAs: gene expression level and species-specific diversity of codon usage based on multivariate analysis

We examined codon usage in Bacillus subtilis genes by multivariate analysis, quantified its cellular levels of individual tRNAs, and found a clear constraint of tRNA contents on synonymous codon choice. Individual tRNA levels were proportional to the copy number of the respective tRNA genes. This indicates that the tRNA gene copy number is an important factor to determine in cellular tRNA levels, which is common with Escherichia coli and yeast Saccharomyces cerevisiae. Codon usage in 18 unicellular organisms whose genomes have been sequenced completely was analyzed and compared with the composition of tRNA genes. The 18 organisms are as follows: yeast S. cerevisiae, Aquifex aeolicus, Archaeoglobus fulgidus, B. subtilis, Borrelia burgdorferi, Chlamydia trachomatis, E. coli, Haemophilus influenzae, Helicobacterpylori, Methanococcusjannaschii, Methanobacterium thermoautotrophicum, Mycobacterium tuberculosis, Mycoplasma genitalium, Mycoplasma pneumoniae, Pyrococcus horikoshii, Rickettsia prowazekii, Synechocystis sp., and Treponema pallidum. Codons preferred in highly expressed genes were related to the codons optimal for the translation process, which were predicted by the composition of isoaccepting tRNA genes. Genes with specific codon usage are discussed in connection with their evolutionary origins and functions. The origin and terminus of replication could be predicted on the basis of codon usage when the usage was analyzed relative to the transcription direction of individual genes.

Exploiting unassigned codons in Micrococcus luteus for tRNA-based amino acid mutagenesis

An alternative to suppression of stop codons for the biosynthetic insertion of non-natural amino acids has been developed. Micrococcus luteus , a Gram-positive bacterium, is incapable of translating at least two codons. One of these unused codons was inserted in a gene to act as a nonsense site. An aminoacylated tRNA was synthesized which was complementary to this codon. The gene containing the missing codon was expressed in vitro in a M.luteus transcription/translation system. Read-through of the missing codon occurred only when the complementary tRNA was included. The results demonstrate that M.luteus can be used for incorporation of amino acids via synthetically prepared aminoacylated tRNAs. The use of a M. luteus translation system provides a method for incorporation of non-natural amino acids which avoids the use of stop codons.

Non-universal decoding of the leucine codon CUG in several Candida species

It has been reported that CUG, a universal leucine codon, is read as serine in an asporogenic yeast, Candida cylindracea. The distribution of this non-universal genetic code in various yeast species was studied using an in vitro translation assay system with a synthetic messenger RNA containing CUG codons in-frame. It was found that CUG is used as a serine codon in six out of the fourteen species examined, while it is used for leucine in the remaining eight. The tRNA species responsible for the translation of codon CUG as serine was detected in all the six species in which CUG is translated as serine. The grouping according to the CUG codon assignments in these yeast species shows a good correlation with physiological classification by the chain lengths of the isoprenoid moiety of ubiquinone and the cell-wall sugar contained in the yeasts. The six Candida species examined in which CUG is used as serine belong to one distinct group in Hemiascomycetes.

Directional mutation pressure, mutator mutations, and dynamics of molecular evolution

Using a general form of the directional mutation theory, this paper analyzes the effect of mutations in mutator genes on the G+C content of DNA, the frequency of substitution mutations, and evolutionary changes (cumulative mutations) under various degrees of selective constraints. Directional mutation theory predicts that when the mutational bias between A/T and G/C nucleotide pairs is equilibrated with the base composition of a neutral set of DNA nucleotides, the mutation frequency per gene will be much lower than the frequency immediately after the mutator mutation takes place. This prediction explains the wide variation of the DNA G+C content among unicellular organisms and possibly also the wide intragenomic heterogeneity of third codon positions for the genes of multicellular eukaryotes. The present analyses lead to several predictions that are not consistent with a number of the frequently held assumptions in the field of molecular evolution, including belief in a constant rate of evolution, symmetric branching of phylogenetic trees, the generality of higher mutation frequency for neutral sets of nucleotides, the notion that mutator mutations are generally deleterious because of their high mutation rates, and teleological explanations of DNA base composition.

Recent evidence for evolution of the genetic code

The genetic code, formerly thought to be frozen, is now known to be in a state of evolution. This was first shown in 1979 by Barrell et al. (G. Barrell, A. T. Bankier, and J. Drouin, Nature [London] 282:189-194, 1979), who found that the universal codons AUA (isoleucine) and UGA (stop) coded for methionine and tryptophan, respectively, in human mitochondria. Subsequent studies have shown that UGA codes for tryptophan in Mycoplasma spp. and in all nonplant mitochondria that have been examined. Universal stop codons UAA and UAG code for glutamine in ciliated protozoa (except Euplotes octacarinatus) and in a green alga, Acetabularia. E. octacarinatus uses UAA for stop and UGA for cysteine. Candida species, which are yeasts, use CUG (leucine) for serine. Other departures from the universal code, all in nonplant mitochondria, are CUN (leucine) for threonine (in yeasts), AAA (lysine) for asparagine (in platyhelminths and echinoderms), UAA (stop) for tyrosine (in planaria), and AGR (arginine) for serine (in several animal orders) and for stop (in vertebrates). We propose that the changes are typically preceded by loss of a codon from all coding sequences in an organism or organelle, often as a result of directional mutation pressure, accompanied by loss of the tRNA that translates the codon. The codon reappears later by conversion of another codon and emergence of a tRNA that translates the reappeared codon with a different assignment. Changes in release factors also contribute to these revised assignments. We also discuss the use of UGA (stop) as a selenocysteine codon and the early history of the code.

Evolution of tRNAs and tRNA genes in Acholeplasma laidlawii

The genes for 22 tRNA species from Acholeplasma laidawii, belonging to the class Mollicutes (Mycoplasmas), have been cloned and sequenced. Sixteen genes are organized in 3 clusters consisting of eleven, three and two tRNA genes, respectively, and the other 6 genes exist as a single gene. The arrangement of tRNA genes in the 11-gene, the 3-gene and the 2-gene clusters reveals extensive similarity to several parts of the 21-tRNA or 16-tRNA gene cluster in Bacillus subtilis. The 11-gene cluster is also similar to the tRNA gene clusters found in other mycoplasma species, the 9-tRNA gene cluster in M.capricolum and in M.mycoides, and the 10-tRNA gene cluster in Spiroplasma meliferm. The results suggest that the tRNA genes in mycoplasmas have evolved from large tRNA gene clusters in the ancestral Gram-positive bacterial genome common to mycoplasmas and B.subtilis. The anticodon sequences including base modifications of 15 tRNA species from A.laidlawii were determined. The anticodon composition and codon-recognition patterns of A.laidlawii resemble those of Bacillus subtilis rather than those of other mycoplasma species.