The Ayerst Award Lecture 1976. Novel factors in protein biosynthesis.

Comparison of nucleotide sequences surrounding the initiation sites of a number of mRNAs reveals few common features. These may be the presence of in- or out-of-phase nonsense codons and (or) polypurine bases complementary to the 16S RNA of the 30S subunit of ribosomes. Since the bases which precede or follow an initiation site vary in length and composition we have examined whether they play a role as spacers between cistrons or whether they have an active function in the termination and initiation of translation. In vitro we have observed that some sequences 5' terminal to AUG are preferred over others in forming an initiation complex. The same bases have much less effect when present at the 3' terminal end of an AUG codon. When the 5' terminal codon is the termination codon UAA, absolutely no initiation complex can be detected. This suggests that spacing may be needed between a stop and a start codon. Conversely, the hexamer AUGUAA failed to elicit chain termination. This was so in systems that terminated when free UAA was added or when a sense triplet was present between the initiation and termination triplets. These results suggest that ribosomes may recognize the stop triplet. Hence ribosomes may not obey simple A and P site models in the termination reaction.     

The comparative method rules! Codon volatility cannot detect positive Darwinian selection using a single genome sequence

All established methods for detecting positive selection at the molecular level rely on comparisons between nucleotide sequences. An exceptional method that purports to detect selection on the basis of a single genomic sequence has recently been proposed. This method uses a measure called "codon volatility," defined for each codon as the ratio between the number of nonsynonymous codons that differ from the codon under study at a single nucleotide position and the number of sense codons that differ from the codon under study at a single nucleotide position. Here, we examine various properties of codon volatility and its derivatives and use simulation of evolutionary processes to determine whether they can be used to detect selective pressures. Codons for only four amino acids (glycine, leucine, arginine, and serine) show any variation in codon volatility. Thus, codon volatility is mainly a proxy for amino acid usage, rather than for codon usage, with 65% of all synonymous changes and 27% of all nonsynonymous changes being undetectable by this measure. Genes identified by the volatility method as being subject to positive selection tend to have idiosyncratic amino acid compositions (e.g., they are glycine rich or arginine poor). An additional property of codon volatility is the near zero variance of its mean expectation, which translates into overestimated statistical significance estimates, especially in the absence of corrections for multiple comparisons. A comparison with measures of selection inferred through comparative methodology reveals no relationship between the results of the two methods. Finally, we show that codon volatility can increase in the absence of positive Darwinian selection; that is, increased codon volatility is not indicative of positive selection.     

The translational stop signal: Codon with a context, or extended factor recognition element?

Wide ranging studies of the readthrough of translational stop codons within the last 25 years have suggested that the stop codon might be only part of the molecular signature for recognition of the termination signal. Such studies do not distinguish between effects on suppression and effects on termination, and so we have used a number of different approaches to deduce whether the stop signal is a codon with a context or an extended factor recognition element. A data base of natural termination sites from a wide range of organisms (148 organisms, 40000 sequences) shows a very marked bias in the bases surrounding the stop codon in the genes for all organisms examined, with the most dramatic bias in the base following the codon (+4). The nature of this base determines the efficiency of the stop signal in vivo, and in Escherichia coli this is reinforced by overexpressing the stimulatory factor, release factor-3. Strong signals, defined by their high relative rates of selecting the decoding release factors, are enhanced whereas weak signals respond relatively poorly. Site-directed cross-linking from the +1, and bases up to +6 but not beyond make close contact with the bacterial release factor-2. The translational stop signal is deduced to be an extended factor recognition sequence with a core element, rather than simply a factor recognition triplet codon influenced by context.     

Codon usage in mitochondrial genomes: distinguishing context-dependent mutation from translational selection

We analyze the frequencies of synonymous codons in animal mitochondrial genomes, focusing particularly on mammals and fish. The frequencies of bases at 4-fold degenerate sites are found to be strongly influenced by context-dependent mutation, which causes correlations between pairs of neighboring bases. There is a pattern of excess of certain dinucleotides and deficit of others that is consistent across large numbers of species, despite the wide variation of single-nucleotide frequencies among species. In many bacteria, translational selection is an important influence on codon usage. In order to test whether translational selection also plays a role in mitochondria, we need to control for context-dependent mutation. Selection for translational accuracy can be detected by comparison of codon usage in conserved and variable sites in the same genes. We give a test of this type that works in the presence of context-dependent mutation. There is very little evidence for translational accuracy selection in the mitochondrial genes considered here. Selection for translational efficiency might lead to preference for codons that match the limited repertoire of anticodons on the mitochondrial tRNAs. This is difficult to detect because the effect would usually be in the same direction in comparable to codon families and so would not cause an observable difference in codon usage between families. Several lines of evidence suggest that this type of selection is weak in most cases. However, we found several cases where unusual bases occur at the wobble position of the tRNA, and in these cases, some evidence for selection on codon usage was found. We discuss the way that these unusual cases are associated with codon reassignments in the mitochondrial genetic code.     

Exploring the energy landscape of the genetic code

New insights into the arrangement of the genetic code table, based on the analysis of the physico-chemical properties of its molecular constituents, are reported in this paper. It will be demonstrated that the code has a twofold symmetry that is not apparent from the conventional code table, but becomes apparent when the codon-anticodon energies are listed for each triplet. The evolutionary development of the current code based on single base replacement mutations (transitions) from an 'iso-energetic' degenerated subset of 16 of the 64 codons is discussed. The energy landscape of all 64 codons is presented. A detailed analysis of the energy changes due to mutations in the 3rd, 1st or 2nd position of a codon reveals that the modern genetic code is highly robust. Changes come in small discrete steps that can be quantified in relation to the thermal noise of the system. The relation of the individual codon to its neighbours in the rearranged codon table can be completely understood based on thermodynamic considerations.     

An analysis of codon usage in mammals: Selection or mutation bias?

Summary A new statistical test has been developed to detect selection on silent sites. This test compares the codon usage within a gene and thus does not require knowledge of which genes are under the greatest selection, that there exist common trends in codon usage across genes, or that genes have the same mutation pattern. It also controls for mutational biases that might be introduced by the adjacent bases. The test was applied to 62 mammalian sequences, the significant codon usage biases were detected in all three species examined (humans, rats, and mice). However, these biases appear not to be the consequence of selection, but of the first base pair in the codon influencing the mutation pattern at the third position.     

Intragenic localization of chain-terminating triplets arising by transitions and transversions

Mutational changes involving transitions can convert only one sense codon to ochre, two codons to amber, and two codons to UGA. One codon, UGG for tryptophan, can be converted by transitions to either amber or UGA. By transversion changes 15 other codons can be converted to ochre and/or amber and/or UGA. Ten amino acids can never be replaced by chain termination as a result of transition and transversion mutagenesis of single base-pairs. For two systems (bacteriophage T4 lysozyme and Escherichia coli K12 tryptophan synthetase A protein) in which the poly-peptide gene product has been completely sequenced one can construct predictive intra-genic distribution maps for the location of all possible chain-terminating mutations arising as a result of transitions and transversions.     

Mutation-selection models of codon substitution and their use to estimate selective strengths on codon usage

Current models of codon substitution are formulated at the levels of nucleotide substitution and do not explicitly consider the separate effects of mutation and selection. They are thus incapable of inferring whether mutation or selection is responsible for evolution at silent sites. Here we implement a few population genetics models of codon substitution that explicitly consider mutation bias and natural selection at the DNA level. Selection on codon usage is modeled by introducing codon-fitness parameters, which together with mutation-bias parameters, predict optimal codon frequencies for the gene. The selective pressure may be for translational efficiency and accuracy or for fine-tuning translational kinetics to produce correct protein folding. We apply the models to compare mitochondrial and nuclear genes from several mammalian species. Model assumptions concerning codon usage are found to affect the estimation of sequence distances (such as the synonymous rate d(S), the nonsynonymous rate d(N), and the rate at the 4-fold degenerate sites d(4)), as found in previous studies, but the new models produced very similar estimates to some old ones. We also develop a likelihood ratio test to examine the null hypothesis that codon usage is due to mutation bias alone, not influenced by natural selection. Application of the test to the mammalian data led to rejection of the null hypothesis in most genes, suggesting that natural selection may be a driving force in the evolution of synonymous codon usage in mammals. Estimates of selection coefficients nevertheless suggest that selection on codon usage is weak and most mutations are nearly neutral. The sensitivity of the analysis on the assumed mutation model is discussed.     

Codon usage in highly expressed genes of Haemophillus influenzae and Mycobacterium tuberculosis: translational selection versus mutational bias

Biases in the codon usage and base compositions at three codon sites in different genes of A+T-rich Gram-negative bacterium Haemophillus influenzae and G+C-rich Gram-positive bacterium Mycobacterium tuberculosis have been examined to address the following questions: (1) whether the synonymous codon usage in organisms having highly skewed base compositions is totally dictated by the mutational bias as reported previously (Sharp, P.M., Devine, K.M., 1989. Codon usage and gene expression level in Dictyostelium discoideum: highly expressed genes do `prefer' optimal codons. Nucleic Acids Res. 17, 5029–5039), or is also controlled by translational selection; (2) whether preference of G in the first codon positions by highly expressed genes, as reported in Escherichia coli (Gutierrez, G., Marquez, L., Marin, A., 1996. Preference for guanosine at first codon position in highly expressed Escherichia coli genes. A relationship with translational efficiency. Nucleic Acids Res. 24, 2525–2527), is true in other bacteria; and (3) whether the usage of bases in three codon positions is species-specific. Result presented here show that even in organisms with high mutational bias, translational selection plays an important role in dictating the synonymous codon usage, though the set of optimal codons is chosen in accordance with the mutational pressure. The frequencies of G-starting codons are positively correlated to the level of expression of genes, as estimated by their Codon Adaptation Index (CAI) values, in M. tuberculosis as well as in H. influenzae in spite of having an A+T-rich genome. The present study on the codon preferences of two organisms with oppositely skewed base compositions thus suggests that the preference of G-starting codons by highly expressed genes might be a general feature of bacteria, irrespective of their overall G+C contents. The ranges of variations in the frequencies of individual bases at the first and second codon positions of genes of both H. influenzae and M. tuberculosis are similar to those of E. coli, implying that though the composition of all three codon positions is governed by a selection-mutation balance, the mutational pressure has little influence in the choice of bases at the first two codon positions, even in organisms with highly biased base compositions.     

Inferring natural selection operating on conservative and radical substitution at single amino acid sites

Natural selection operating on amino acid substitution at single amino acid sites can be detected by comparing the rates of synonymous (r(S)) and nonsynonymous (r(N)) nucleotide substitution at single codon sites. Amino acid substitutions can be classified as conservative or radical according to whether they retain the properties of the substituted amino acid. Here methods for comparing the rates of conservative (r(C)) and radical (r(R)) nonsynonymous substitution with r(S) at single codon sites were developed to detect natural selection operating on these substitutions at single amino acid sites. A method for comparing r(C) and r(R) at single codon sites was also developed to detect biases toward these substitutions at single amino acid sites. Charge was used as the property of the amino acids. In a computer simulation, false-positive rates of these methods were always < 5%, unless termination sites were included in the computation of the numbers of sites and estimates of transition/transversion rate ratio were highly biased. The frequency of detection of natural selection operating on conservative substitution was almost independent of the presence of natural selection operating on radical substitution, and vice versa. Natural selection operating specifically on conservative and radical substitution was detected more efficiently by comparing r(S) with r(C) and r(S) with r(R) than by comparing r(S) with r(N). These methods also appeared to be robust against the occurrence of recombination during evolution. In an analysis of class I human leukocyte antigen, negative selection operating on conservative substitution, but not positive selection operating on radical substitution, was observed at some of the codon sites with r(R) > r(C), suggesting that r(R) > r(C) may not necessarily be an indicator of positive selection operating on radical substitution.     

Mutagenesis of the three bases preceding the start codon of the beta-galactosidase mRNA and its effect on translation in Escherichia coli.

The effect on the translation efficiency of various mutations in the three bases (the -1 triplet) that precede the AUG start codon of the beta-galactosidase mRNA in Escherichia coli was studied. Of the 39 mutants examined, the level of expression varies over a 20-fold range. The most favorable combinations of bases in the -1 triplet are UAU and CUU. The expression levels in the mutants with UUC, UCA or AGG as the -1 triplet are 20-fold lower than those with UAU or CUU. In general, a U residue immediately preceding the start codon is more favorable for expression than any other base; furthermore, an A residue at the -2 position enhances the translation efficiency in most instances. In both cases, however, the degree of enhancement depends on its context, i.e. the neighboring bases. Although the rules derived from this study are complex, the results show that mutations in any of the three bases preceding the start codon can strongly affect the translational efficiency of the beta-galactosidase mRNA.     

Advantages of a mechanistic codon substitution model for evolutionary analysis of protein-coding sequences

BACKGROUND: A mechanistic codon substitution model, in which each codon substitution rate is proportional to the product of a codon mutation rate and the average fixation probability depending on the type of amino acid replacement, has advantages over nucleotide, amino acid, and empirical codon substitution models in evolutionary analysis of protein-coding sequences. It can approximate a wide range of codon substitution processes. If no selection pressure on amino acids is taken into account, it will become equivalent to a nucleotide substitution model. If mutation rates are assumed not to depend on the codon type, then it will become essentially equivalent to an amino acid substitution model. Mutation at the nucleotide level and selection at the amino acid level can be separately evaluated. RESULTS: The present scheme for single nucleotide mutations is equivalent to the general time-reversible model, but multiple nucleotide changes in infinitesimal time are allowed. Selective constraints on the respective types of amino acid replacements are tailored to each gene in a linear function of a given estimate of selective constraints. Their good estimates are those calculated by maximizing the respective likelihoods of empirical amino acid or codon substitution frequency matrices. Akaike and Bayesian information criteria indicate that the present model performs far better than the other substitution models for all five phylogenetic trees of highly-divergent to highly-homologous sequences of chloroplast, mitochondrial, and nuclear genes. It is also shown that multiple nucleotide changes in infinitesimal time are significant in long branches, although they may be caused by compensatory substitutions or other mechanisms. The variation of selective constraint over sites fits the datasets significantly better than variable mutation rates, except for 10 slow-evolving nuclear genes of 10 mammals. An critical finding for phylogenetic analysis is that assuming variable mutation rates over sites lead to the overestimation of branch lengths.     

The Mechanisms of Codon Reassignments in Mitochondrial Genetic Codes

Many cases of nonstandard genetic codes are known in mitochondrial genomes. We carry out analysis of phylogeny and codon usage of organisms for which the complete mitochondrial genome is available, and we determine the most likely mechanism for codon reassignment in each case. Reassignment events can be classified according to the gain-loss framework. The "gain" represents the appearance of a new tRNA for the reassigned codon or the change of an existing tRNA such that it gains the ability to pair with the codon. The "loss" represents the deletion of a tRNA or the change in a tRNA so that it no longer translates the codon. One possible mechanism is codon disappearance (CD), where the codon disappears from the genome prior to the gain and loss events. In the alternative mechanisms the codon does not disappear. In the unassigned codon mechanism, the loss occurs first, whereas in the ambiguous intermediate mechanism, the gain occurs first. Codon usage analysis gives clear evidence of cases where the codon disappeared at the point of the reassignment and also cases where it did not disappear. CD is the probable explanation for stop to sense reassignments and a small number of reassignments of sense codons. However, the majority of sense-to-sense reassignments cannot be explained by CD. In the latter cases, by analysis of the presence or absence of tRNAs in the genome and of the changes in tRNA sequences, it is sometimes possible to distinguish between the unassigned codon and the ambiguous intermediate mechanisms. We emphasize that not all reassignments follow the same scenario and that it is necessary to consider the details of each case carefully.     

Identification of a novel in-frame translational stop codon in human intestine apoB mRNA

Human apolipoprotein (apo) B exists in plasma as two isoproteins designated apoB-100 and apoB-48. ApoB-100 (512 kDa) and apoB-48 (250 kDa) are synthesized by the liver and intestine respectively. Analysis of apoB cDNA clones isolated from a human intestinal cDNA library revealed that the intestinal apoB mRNA contains a new in-frame translational stop codon. This premature stop codon is generated by a single base substitution of a 'C' to 'T' at nucleotide 6538 which converts the codon 'CAA' coding for the amino acid glutamine residue 2153 to an in-frame stop codon 'TAA'. The generation of a stop codon in the intestinal apoB mRNA appears to be tissue specific since it has not been reported in cDNA clones isolated from human liver cDNA libraries which code for the 4536 amino acid apoB-100. A potential polyadenylation signal sequence 'AATAAA' was also identified 390 bases downstream from the new stop codon. The new stop codon in the human intestinal apoB mRNA provides a potential mechanism for the biosynthesis of intestinal apoB-48.     

Aminoglycoside antibiotics mediate context-dependent suppression of termination codons in a mammalian translation system

The translation machinery recognizes codons that enter the ribosomal A site with remarkable accuracy to ensure that polypeptide synthesis proceeds with a minimum of errors. When a termination codon enters the A site of a eukaryotic ribosome, it is recognized by the release factor eRF1. It has been suggested that the recognition of translation termination signals in these organisms is not limited to a simple trinucleotide codon, but is instead recognized by an extended tetranucleotide termination signal comprised of the stop codon and the first nucleotide that follows. Interestingly, pharmacological agents such as aminoglycoside antibiotics can reduce the efficiency of translation termination by a mechanism that alters this ribosomal proofreading process. This leads to the misincorporation of an amino acid through the pairing of a near-cognate aminoacyl tRNA with the stop codon. To determine whether the sequence context surrounding a stop codon can influence aminoglycoside-mediated suppression of translation termination signals, we developed a series of readthrough constructs that contained different tetranucleotide termination signals, as well as differences in the three bases upstream and downstream of the stop codon. Our results demonstrate that the sequences surrounding a stop codon can play an important role in determining its susceptibility to suppression by aminoglycosides. Furthermore, these distal sequences were found to influence the level of suppression in remarkably distinct ways. These results suggest that the mRNA context influences the suppression of stop codons in response to subtle differences in the conformation of the ribosomal decoding site that result from aminoglycoside binding.     

Amino acid sequence of guinea pig liver transglutaminase from its cDNA sequence

Transglutaminases (EC 2.3.2.13) catalyze the formation of epsilon-(gamma-glutamyl)lysine cross-links and the substitution of a variety of primary amines for the gamma-carboxamide groups of protein-bound glutaminyl residues. These enzymes are involved in many biological phenomena. In this paper, the complete amino acid sequence of guinea pig liver transglutaminase, a typical tissue-type nonzymogenic transglutaminase, was predicted by the cloning and sequence analysis of DNA complementary to its mRNA. The cDNA clones carrying the sequences for the 5'- and 3'-end regions of mRNA were obtained by use of the sequence of the partial-length cDNA of guinea pig liver transglutaminase [Ikura, K., Nasu, T., Yokota, H., Sasaki, R., & Chiba, H. (1987) Agric. Biol. Chem. 51, 957-961]. A total of 3695 bases were identified from sequence data of four overlapping cDNA clones. Northern blot analysis of guinea pig liver poly(A+) RNA showed a single species of mRNA with 3.7-3.8 kilobases, indicating that almost all of the mRNA sequence was analyzed. The composite cDNA sequence contained 68 bases of a 5'-untranslated region, 2073 bases of an open reading frame that encoded 691 amino acids, a stop codon (TAA), 1544 bases of a 3'-noncoding region, and a part of a poly(A) tail (7 bases). The molecular weight of guinea pig liver transglutaminase was calculated to be 76,620 from the amino acid sequence deduced, excluding the initiator Met. This enzyme contained no carbohydrate [Folk, J. E., & Chung, S. I. (1973) Adv. Enzymol. Relat. Areas Mol. Biol. 38, 109-191], but six potential Asn-linked glycosylation sites were found in the sequence deduced.(ABSTRACT TRUNCATED AT 250 WORDS)     

Identification of a ribosome-binding site for a leader peptide encoded by Rous sarcoma virus RNA.

A new method for identifying ribosome-binding sites was developed to determine whether AUG codons in the 5'-terminal RNA sequence of Rous sarcoma virus were used to initiate protein synthesis. We found that when translation is inhibited, the major ribosome-binding site on Rous sarcoma virus RNA is at the 5'-proximal AUG codon, even though the primary translational product from this RNA, Pr76gag, is encoded behind the fourth AUG codon 331 bases downstream from the observed initiation site. These results suggest that ribosomes can initiate translation on Rous sarcoma virus RNA at more than one site, thereby producing a seven-amino-acid peptide, as well as the gag gene polyprotein precursor of Mr 76,000.     

Nondefective rotavirus mutants with an NSP1 gene which has a deletion of 500 nucleotides, including a cysteine-rich zinc finger motif-encoding region (nucleotides 156 to 248), or which has a nonsense codon at nucleotides 153-155.

We isolated two nondefective bovine rotavirus mutants (A5-10 and A5-16 clones) which have nonsense mutations in the early portion of the open reading frame of the NSP1 gene. In the NSP1 gene (1,587 bases long) of A5-10, a nonsense codon is present at nucleotides 153 to 155 just upstream of the coding region (nucleotides 156 to 230) of a cysteine-rich Zn finger motif. A5-16 gene 5 (1,087 bases long) was found to have a large deletion of 500 bases corresponding to nucleotides 142 to 641 of a parent A5-10 NSP1 gene and to have a nonsense codon at nucleotides 183 to 185, which resulted from the deletion. Expression of gene 5-specific NSP1 could not be detected in MA-104 cells infected with the A5-10 or A5-16 clone or in an in vitro translation system using the plasmids with gene 5 cDNA from A5-10 or A5-16. Nevertheless, both A5-10 and A5-16 replicated well in cultured cells, although the plaque size of A5-16 was extremely small.