The Effects of Mutation and Natural Selection on Codon Bias in the Genes of Drosophila

Codon bias varies widely among the loci of Drosophila melanogaster, and some of this diversity has been explained by variation in the strength of natural selection. A study of correlations between intron and coding region base composition shows that variation in mutation pattern also contributes to codon bias variation. This finding is corroborated by an analysis of variance (ANOVA), which shows a tendency for introns from the same gene to be similar in base composition. The strength of base composition correlations between introns and codon third positions is greater for genes with low codon bias than for genes with high codon bias. This pattern can be explained by an overwhelming effect of natural selection, relative to mutation, in highly biased loci. In particular, this correlation is absent when examining fourfold degenerate sites of highly biased genes. In general, it appears that selection acts more strongly in choosing among fourfold degenerate codons than among twofold degenerate codons. Although the results indicate regional variation in mutational bias, no evidence is found for large scale regions of compositional homogeneity.     

Codon usage in plant genes.

We have examined codon bias in 207 plant gene sequences collected from Genbank and the literature. When this sample was further divided into 53 monocot and 154 dicot genes, the pattern of relative use of synonymous codons was shown to differ between these taxonomic groups, primarily in the use of G + C in the degenerate third base. Maize and soybean codon bias were examined separately and followed the monocot and dicot codon usage patterns respectively. Codon preference in ribulose 1,5 bisphosphate and chlorophyll a/b binding protein, two of the most abundant proteins in leaves was investigated. These highly expressed are more restricted in their codon usage than plant genes in general.     

Working of the genetic code

The genetic code is used differently by different kinds of species. Each type of genome has a particular coding strategy, that is, choices among degenerate bases are consistently similar for all genes therein. This uniformity in the selection between degenerate bases within each taxonomic group has been discovered by applying new methods to the study of coding variability. It is now possible to calculate relative distances between genomes, or genome types, based on use of the codon catalog by the mRNAs therein.     

Preferential codon usage in genes

We present a method which permits comparison of the preferential use of degenerate codons within any gene. The method makes use of the triplet frequencies in the noncoding frames to assess whether a preference is specific to the reading frame. Preference is given a statistical meaning by use of the analysis of variance coupled to Duncan's multiple range test.Preferential use of degenerate codons is gene-specific and independent of gene size. The data suggest that any correlation between codon frequency distribution and tRNA levels is unreliable. In those animal genes examined, codons ending in C or G are preferred; in animal viruses tested, codons ending in U or A are preferred. Similarly, the bacterial genes and the genes of single-stranded DNA phages that we analyzed differed from each other as well as from eukaryotic genes in the third base of the codon.     

Coding strategy differences between constant and variable segments of immunoglobulin genes

Vertebrate immunoglobulin (Ig) mRNAs reveal intraspecies variation in codon usage distinct from that seen with yeast or bacterial genes. Comparison of all available Ig gene sequences shows that %(G + C) in codon position III is consistently lower in variable (V) segments than in constant (C) segments. I find an even lower %(G + C) in the hypervariable domains of V segments. This analysis suggests that base substitution in Ig genes correlates positively with local A + T content.     

Codon usage bias in prokaryotic pyrimidine-ending codons is associated with the degeneracy of the encoded amino acids

Synonymous codons are unevenly distributed among genes, a phenomenon termed codon usage bias. Understanding the patterns of codon bias and the forces shaping them is a major step towards elucidating the adaptive advantage codon choice can confer at the level of individual genes and organisms. Here, we perform a large-scale analysis to assess codon usage bias pattern of pyrimidine-ending codons in highly expressed genes in prokaryotes. We find a bias pattern linked to the degeneracy of the encoded amino acid. Specifically, we show that codon-pairs that encode two- and three-fold degenerate amino acids are biased towards the C-ending codon while codons encoding four-fold degenerate amino acids are biased towards the U-ending codon. This codon usage pattern is widespread in prokaryotes, and its strength is correlated with translational selection both within and between organisms. We show that this bias is associated with an improved correspondence with the tRNA pool, avoidance of mis-incorporation errors during translation and moderate stability of codon-anticodon interaction, all consistent with more efficient translation.     

The CUG codon is decoded in vivo as serine and not leucine in Candida albicans.

Previous studies have shown that the yeast Candida albicans encodes a unique seryl-tRNA(CAG) that should decode the leucine codon CUG as serine. However, in vitro translation of several different CUG-containing mRNAs in the presence of this unusual seryl-tRNA(CAG) result in an apparent increase in the molecular weight of the encoded polypeptides as judged by SDS-PAGE even though the molecular weight of serine is lower than that of leucine. A possible explanation for this altered electrophoretic mobility is that the CUG codon is decoded as modified serine in vitro. To elucidate the nature of CUG decoding in vivo, a reporter system based on the C. albicans gene (RBP1) encoding rapamycin-binding protein (RBP), coupled to the promoter of the C. albicans TEF3 gene, was utilized. Sequencing and mass-spectrometry analysis of the recombinant RBP expressed in C. albicans demonstrated that the CUG codon was decoded exclusively as serine while the related CUU codon was translated as leucine. A database search revealed that 32 out of the 65 C. albicans gene sequences available have CUG codons in their open reading frames. The CUG-containing genes do not belong to any particular gene family. Thus the amino acid specified by the CUG codon has been reassigned within the mRNAs of C. albicans. We argue here that this unique genetic code change in cellular mRNAs cannot be explained by the 'Codon Reassignment Theory'.     

Theoretical foundations for quantitative paleogenetics

Summary REH theory is extended by deriving the theoretical equations that permit one to analyze the nonrandom molecular divergence of homologous genes and proteins. The nonrandomicities considered are amino acid and base composition, the frequencies with which each of the four nucleotides is replaced by one of the other three, unequal usage of degenerate codons, distribution of fixed base replacements at the three nucleotide positions within codons, and distributions of fixed base replacements among codons. The latter two distributions turn out to dominate the accuracy of genetic distance estimates. The negative binomial density is used to allow for the unequal mutability of different codon sites, and the implications of its two limiting forms, the Poisson and geometric distributions, are considered. It is shown that the fixation intensity — the average number of base replacements per variable codon - is expressible as the simple product of two factors, the first describing the asymmetry of the distribution of base replacements over the gene and the second defining the ratio of the average probability that a codon will fix a mutation to the probability that it will not. Tables are given relating these features to experimentally observable quantities in hemoglobin, hemoglobin, myoglobin, cytochromec, and the parvalbumin group of proteins and to the structure of their corre-sponding genes or mRNAs. The principal results are (1) more accurate methods of estimating parameters of evolutionary interest from experimental gene and protein sequence data, and (2) the fact that change in gene and protein structure has been a much less efficient process than previously believed in the sense of requiring many more base replacements to effect a given structural change than earlier estimation procedures had indicated. This inefficiency is directly traceable to Darwinian selection for the nonrandom gene or protein structures necessary for biological function. The application of these methods is illustrated by detailed consideration of the rabbit -and hemoglobin mRNAs and the proteins for which they code. It is found that these two genes are separated by about 425 fixed base replacements, which is a factor of two greater than earlier estimates. The replacements are distributed over approximately 114 codon sites that were free to accept base mutations during the divergence of these two genes.     

Variable strength of translational selection among 12 Drosophila species

Codon usage bias in Drosophila melanogaster genes has been attributed to negative selection of those codons whose cellular tRNA abundance restricts rates of mRNA translation. Previous studies, which involved limited numbers of genes, can now be compared against analyses of the entire gene complements of 12 Drosophila species whose genome sequences have become available. Using large numbers (6138) of orthologs represented in all 12 species, we establish that the codon preferences of more closely related species are better correlated. Differences between codon usage biases are attributed, in part, to changes in mutational biases. These biases are apparent from the strong correlation (r = 0.92, P < 0.001) among these genomes' intronic G + C contents and exonic G + C contents at degenerate third codon positions. To perform a cross-species comparison of selection on codon usage, while accounting for changes in mutational biases, we calibrated each genome in turn using the codon usage bias indices of highly expressed ribosomal protein genes. The strength of translational selection was predicted to have varied between species largely according to their phylogeny, with the D. melanogaster group species exhibiting the strongest degree of selection.     

Messenger RNAs of a strongly-expressed late gene of cowpox virus contain 5''-terminal poly(A) sequences.

We have identified and characterized one of the most strongly-expressed genes of cowpox virus (CPV). This is the gene encoding the major protein component of the A-type inclusion bodies produced by this virus. This gene (designated the 160K gene) is transcribed late during the infection. Analyses of its mRNAs showed that these late RNAs, unlike all other characterized late mRNAs of poxviruses, are uniform in length. However, the most remarkable feature of the mRNAs of the 160K gene is the structure of their 5'-termini. Most of these mRNAs have 5'-terminal poly(A) sequences containing 5-21 residues. Furthermore, these 5'-terminal poly(A) sequences are not complementary to the corresponding region of the template strand of the viral DNA. Instead, the nucleotide sequences of the mRNA and the viral DNA diverge at the site of the three As in the sequence 5'-TAAATG-3' containing the gene's initiation codon. Consequently, the poly(A) provides the leader sequences of these mRNAs. These unusual 5'-terminal structures suggest that the late mRNAs of pox-virus genes are generated by a novel process.     

Local absence of secondary structure permits translation of mRNAs that lack ribosome-binding sites

The initiation of translation is a fundamental and highly regulated process in gene expression. Translation initiation in prokaryotic systems usually requires interaction between the ribosome and an mRNA sequence upstream of the initiation codon, the so-called ribosome-binding site (Shine-Dalgarno sequence). However, a large number of genes do not possess Shine-Dalgarno sequences, and it is unknown how start codon recognition occurs in these mRNAs. We have performed genome-wide searches in various groups of prokaryotes in order to identify sequence elements and/or RNA secondary structural motifs that could mediate translation initiation in mRNAs lacking Shine-Dalgarno sequences. We find that mRNAs without a Shine-Dalgarno sequence are generally less structured in their translation initiation region and show a minimum of mRNA folding at the start codon. Using reporter gene constructs in bacteria, we also provide experimental support for local RNA unfoldedness determining start codon recognition in Shine-Dalgarno--independent translation. Consistent with this, we show that AUG start codons reside in single-stranded regions, whereas internal AUG codons are usually in structured regions of the mRNA. Taken together, our bioinformatics analyses and experimental data suggest that local absence of RNA secondary structure is necessary and sufficient to initiate Shine-Dalgarno--independent translation. Thus, our results provide a plausible mechanism for how the correct translation initiation site is recognized in the absence of a ribosome-binding site.     

Impact of translational selection on codon usage bias in the archaeon Methanococcus maripaludis

Patterns of codon usage have been extensively studied among Bacteria and Eukaryotes, but there has been little investigation of species from the third domain of life, the Archaea. Here, we examine the nature of codon usage bias in a methanogenic archaeon, Methanococcus maripaludis. Genome-wide patterns of codon usage are dominated by a strong A + T bias, presumably largely reflecting mutation patterns. Nevertheless, there is variation among genes in the use of a subset of putatively translationally optimal codons, which is strongly correlated with gene expression level. In comparison with Bacteria such as Escherichia coli, the strength of selected codon usage bias in highly expressed genes in M. maripaludis seems surprisingly high given its moderate growth rate. However, the pattern of selected codon usage differs between M. maripaludis and E. coli: in the archaeon, strongly selected codon usage bias is largely restricted to twofold degenerate amino acids (AAs). Weaker bias among the codons for fourfold degenerate AAs is consistent with the small number of tRNA genes in the M. maripaludis genome.     

Characterization of mouse muc6 and evidence of conservation of the gel-forming mucin gene cluster between human and mouse

Using degenerate primers designed from conserved cysteine-rich domains of gel-forming mucins, we cloned two new mouse mucin cDNAs. Blast searching showed that they belong to the same new gene assigned to chromosome 7 band F5. This gene is clustered with the three secreted large gel-forming mucins Muc2, Muc5ac, and Muc5b in a region that exhibits synteny with human chromosome 11p15. Computer analysis and sequence alignments with mucin genes predict that the new gene is composed of 33 exons and spans 30 kb from the initiation ATG codon to the Stop codon. Sequence similarities, domain organization of the deduced peptide, and expression analysis allow us to conclude that this newly cloned mouse gene is Muc6, i.e., the mouse ortholog of human MUC6. Like those of their human homologs, the genomic order and arrangement of the four mucins within the cluster of mucin genes are conserved.     

mRNAs have greater negative folding free energies than shuffled or codon choice randomized sequences

An examination of 51 mRNA sequences in GenBank has revealed that calculated mRNA folding is more stable than expected by chance. Free energy minimization calculations of native mRNA sequences are more negative than randomized mRNA sequences with the same base composition and length. Randomization of the coding region of genes yields folding free energies of less negative magnitude than the original native mRNA sequence. Randomization of codon choice, while still preserving original base composition, also results in less stable mRNAs. This suggests that a bias in the selection of codons favors the potential formation of mRNA structures which contribute to folding stability.     

Alpha-thalassemia due to the deletion of nucleotides -2 and -3 preceding the AUG initiation codon affects translation efficiency both in vitro and in vivo.

We previously hypothesized that a 2 nucleotide deletion, causing a A-greater than C change at position -3 preceding the ATG initiation codon of alpha globin gene, reduced translation efficiency of alpha globin mRNA and was responsible for a form of alpha + thalassemia displayed by an Algerian patient. We presently show that this deletion leads to a 30-45% reduction in translation efficiency of synthetic alpha globin mRNA in rabbit reticulocyte lysate. In other experiments, we constructed alpha/G gamma hybrid globin genes in which the 3' end of normal or mutated alpha globin genes downstream to the ATG initiation codon was substituted by the 3' part of a G gamma globin gene. COS cells transfected with either of these 2 hybrid genes were shown to synthesize a similar amount of alpha/G gamma hybrid mRNAs but 50% less G gamma globin when transfected with the alpha/G gamma hybrid gene carrying the deletion. These results definitively establish that the 2 nucleotide deletion reduces translation efficiency by 30-50%. This contrasts with the 93% reduction induced by a similar A-greater than C change at position -3 in the different nucleotide context preceding the ATG codon of the rat preproinsulin gene.     

Preferential codon usage in prokaryotic genes: the optimal codon-anticodon interaction energy and the selective codon usage in efficiently expressed genes

By considering the nucleotide sequence of several highly expressed coding regions in bacteriophage MS2 and mRNAs from Escherichia coli, it is possible to deduce some rules which govern the selection of the most appropriate synonymous codons NNU or NNC read by tRNAs having GNN, QNN or INN as anticodon. The rules fit with the general hypothesis that an efficient in-phase translation is facilitated by proper choice of degenerate codewords promoting a codon-anticodon interaction with intermediate strength (optimal energy) over those with very strong or very weak interaction energy. Moreover, codons corresponding to minor tRNAs are clearly avoided in these efficiently expressed genes. These correlations are clearcut in the normal reading frame but not in the corresponding frameshift sequences +1 and +2. We hypothesize that both the optimization of codon-anticodon interaction energy and the adaptation of the population to codon frequency or vice versa in highly expressed mRNAs of E. coli are part of a strategy that optimizes the efficiency of translation. Conversely, codon usage in weakly expressed genes such as repressor genes follows exactly the opposite rules. It may be concluded that, in addition to the need for coding an amino acid sequence, the energetic consideration for codon-anticodon pairing, as well as the adaptation of codons to the tRNA population, may have been important evolutionary constraints on the selection of the optimal nucleotide sequence.     

Codon replacement in the PGK1 gene of Saccharomyces cerevisiae: experimental approach to study the role of biased codon usage in gene expression.

The coding sequences of genes in the yeast Saccharomyces cerevisiae show a preference for 25 of the 61 possible coding triplets. The degree of this biased codon usage in each gene is positively correlated to its expression level. Highly expressed genes use these 25 major codons almost exclusively. As an experimental approach to studying biased codon usage and its possible role in modulating gene expression, systematic codon replacements were carried out in the highly expressed PGK1 gene. The expression of phosphoglycerate kinase (PGK) was studied both on a high-copy-number plasmid and as a single copy gene integrated into the chromosome. Replacing an increasing number (up to 39% of all codons) of major codons with synonymous minor ones at the 5' end of the coding sequence caused a dramatic decline of the expression level. The PGK protein levels dropped 10-fold. The steady-state mRNA levels also declined, but to a lesser extent (threefold). Our data indicate that this reduction in mRNA levels was due to destabilization caused by impaired translation elongation at the minor codons. By preventing translation of the PGK mRNAs by the introduction of a stop codon 3' and adjacent to the start codon, the steady-state mRNA levels decreased dramatically. We conclude that efficient mRNA translation is required for maintaining mRNA stability in S. cerevisiae. These findings have important implications for the study of the expression of heterologous genes in yeast cells.     

Translation efficiencies of synonymous codons are not always correlated with codon usage in tobacco chloroplasts

Codon usage in chloroplasts is different from that in prokaryotic and eukaryotic nuclear genomes. However, no experimental approach has been made to analyse the translation efficiency of individual codons in chloroplasts. We devised an in vitro assay for translation efficiencies using synthetic mRNAs, and measured the translation efficiencies of five synonymous codon groups in tobacco chloroplasts. Among four alanine codons (GCN, where N is U, C, A or G), GCU was the most efficient for translation, whereas the chloroplast genome lacks tRNA genes corresponding to GCU. Phenylalanine and tyrosine are each encoded by two codons (UUU/C and UAU/C, respectively). Phenylalanine UUC and tyrosine UAC were translated more than twice as efficiently than UUU and UAU, respectively, contrary to their codon usage, whereas translation efficiencies of synonymous codons for alanine, aspartic acid and asparagine were parallel to their codon usage. These observations indicate that translation efficiencies of individual codons are not always correlated with codon usage in vitro in chloroplasts. This raises an important issue for foreign gene expression in chloroplasts.     

Compilation and comparison of the sequence context around the AUG startcodons in Saccharomyces cerevisiae mRNAs.

The nucleotide sequence of the translation initiation regions of 96 Saccharomyces cerevisiae mRNAs was compiled and compared. The entire 5' untranslated sequence of most mRNAs is very rich in A-residues. G-residues are underrepresented in the untranslated region. The AUG startcodon context appeared to be distinctly different from that of animal mRNAs, although an A-residue at -3 also occurs very frequently (81 percent) in yeast mRNAs. The prevailing codon 3' adjacent to the AUG is the UCU serine codon. All these features are more extreme in the highly expressed genes. Fifty percent of all highly expressed genes use the UCU serine codon as second triplet. In this group G-residues are completely absent in the 7 bases preceding the startcodon and an A-residue occurs at position -1 and -3 at a frequency of 89 percent and 100 percent, respectively. The abundance of A-residues throughout the leader suggests that unstructured mRNA is required for efficient translation initiation in yeast. The consensus sequence for the AUG context in highly expressed genes can be summarized as follows: (Sequence: see text).