Role of premature stop codons in bacterial evolution

When the stop codons TGA, TAA, and TAG are found in the second and third reading frames of a protein-encoding gene, they are considered premature stop codons (PSC). Deinococcus radiodurans disproportionately favored TGA more than the other two triplets as a PSC. The TGA triplet was also found more often in noncoding regions and as a stop codon, though the bias was less pronounced. We investigated this phenomenon in 72 bacterial species with widely differing chromosomal GC contents. Although TGA and TAG were compositionally similar, we found a great variation in use of TGA but a very limited range of use of TAG. The frequency of use of TGA in the gene sequences generally increased with the GC content of the chromosome, while the frequency of use of TAG, like that of TAA, was inversely proportional to the GC content of the chromosome. The patterns of use of TAA, TGA and TAG as real stop codons were less biased and less influenced by the GC content of the chromosome. Bacteria with higher chromosomal GC contents often contained fewer PSC trimers in their genes. Phylogenetically related bacteria often exhibited similar PSC ratios. In addition, metabolically versatile bacteria have significantly fewer PSC trimers in their genes. The bias toward TGA but against TAG as a PSC could not be explained either by the preferential usage of specific codons or by the GC contents of individual chromosomes. We proposed that the quantity and the quality of the PSC in the genome might be important in bacterial evolution.     

Variation in Release Factor Abundance Is Not Needed to Explain Trends in Bacterial Stop Codon Usage

In bacteria stop codons are recognized by one of two class I release factors (RF1) recognizing TAG, RF2 recognizing TGA, and TAA being recognized by both. Variation across bacteria in the relative abundance of RF1 and RF2 is thus hypothesized to select for different TGA/TAG usage. This has been supported by correlations between TAG:TGA ratios and RF1:RF2 ratios across multiple bacterial species, potentially also explaining why TAG usage is approximately constant despite extensive variation in GC content. It is, however, possible that stop codon trends are determined by other forces and that RF ratios adapt to stop codon usage, rather than vice versa. Here, we determine which direction of the causal arrow is the more parsimonious. Our results support the notion that RF1/RF2 ratios become adapted to stop codon usage as the same trends, notably the anomalous TAG behavior, are seen in contexts where RF1:RF2 ratios cannot be, or are unlikely to be, causative, that is, at 3′untranslated sites never used for translation termination, in intragenomic analyses, and across archaeal species (that possess only one RF1). We conclude that specifics of RF biology are unlikely to fully explain TGA/TAG relative usage. We discuss why the causal relationships for the evolution of synonymous stop codon usage might be different from those affecting synonymous sense codon usage, noting that transitions between TGA and TAG require two-point mutations one of which is likely to be deleterious.     

alpha-L-iduronidase premature stop codons and potential read-through in mucopolysaccharidosis type I patients

alpha-L-Iduronidase is a glycosyl hydrolase involved in the sequential degradation of the glycosaminoglycans heparan sulphate and dermatan sulphate. A deficiency in alpha-L-iduronidase results in the lysosomal accumulation and urinary secretion of partially degraded glycosaminoglycans and is the cause of the lysosomal storage disorder mucopolysaccharidosis type I (MPS I; Hurler and Scheie syndromes; McKusick 25280). The premature stop codons Q70X and W402X are two of the most common alpha-l-iduronidase gene (IDUA) mutations accounting for up to 70% of MPS I disease alleles in some populations. Here, we have reported a new mutation, making a total of 15 different mutations that can cause premature IDUA stop codons and have investigated the biochemistry of these mutations. Natural stop codon read-through was dependent on the fidelity of the codon when evaluated at Q70X and W402X in CHO-K1 cells, but the three possible stop codons TAA, TAG and TGA, had different effects on mRNA stability and this effect was context dependent. In CHO-K1 cells expressing the Q70X and W402X mutations, the level of gentamicin-enhanced stop codon read-through was slightly less than the increment in activity caused by a lower fidelity stop codon. In this system, gentamicin had more effect on read-through for the TAA and TGA stop codons when compared to the TAG stop codon. In an MPS I patient study, premature TGA stop codons were associated with a slightly attenuated clinical phenotype, when compared to classical Hurler syndrome (e.g. W402X/W402X and Q70X/Q70X genotypes with TAG stop codons). Natural read-through of premature stop codons is a potential explanation for variable clinical phenotype in MPS I patients. Enhanced stop codon read-through is a potential treatment strategy for a large sub-group of MPS I patients.     

A relationship between GC content and coding-sequence length

Since base composition of translational stop codons (TAG, TAA, and TGA) is biased toward a low G+C content, a differential density for these termination signals is expected in random DNA sequences of different base compositions. The expected length of reading frames (DNA segments of sense codons flanked by in-phase stop codons) in random sequences is thus a function of GC content. The analysis of DNA sequences from several genome databases stratified according to GC content reveals that the longest coding sequences—exons in vertebrates and genes in prokaryotes—are GC-rich, while the shortest ones are GC-poor. Exon lengthening in GC-rich vertebrate regions does not result, however, in longer vertebrate proteins, perhaps because of the lower number of exons in the genes located in these regions. The effects on coding-sequence lengths constitute a new evolutionary meaning for compositional variations in DNA GC content. Correspondence to: J. L. Oliver     

The close proximity ofEscherichia coli genes: Consequences for stop codon and synonymous codon use

It is shown that synonymous codon usage is less biased in favor of those codons preferred by highly expressed genes at the end ofEscherichia coli genes than in the middle. This appears to be due to the close proximity of manyE. coli genes. It is shown that a substantial number of genes overlap either the Shine-Dalgarno sequence or the coding sequence of the next gene on the chromosome and that the codons that overlap have lower synonymous codon bias than those which do not. It is also shown that there is an increase in the frequency of A-ending codons, and a decrease in the frequency of G-ending codons at the end ofE. coli genes that lie close to another gene. It is suggested that these trends in composition could be associated with selection against the formation of mRNA secondary structure near the start of the next gene on the chromosome. Stop codon use is also affected by the close proximity of genes; many genes are forced to use TGA and TAG stop codons because they terminate either within the Shine-Dalgarno or coding sequence of the next gene on the chromosome. The implications these results have for the evolution of synonymous codon use are discussed.     

GC-biased gene conversion and selection affect GC content in the Oryza genus (rice)

Base composition varies among and within eukaryote genomes. Although mutational bias and selection have initially been invoked, more recently GC-biased gene conversion (gBGC) has been proposed to play a central role in shaping nucleotide landscapes, especially in yeast, mammals, and birds. gBGC is a kind of meiotic drive in favor of G and C alleles, associated with recombination. Previous studies have also suggested that gBGC could be at work in grass genomes. However, these studies were carried on third codon positions that can undergo selection on codon usage. As most preferred codons end in G or C in grasses, gBGC and selection can be confounded. Here we investigated further the forces that might drive GC content evolution in the rice genus using both coding and noncoding sequences. We found that recombination rates correlate positively with equilibrium GC content and that selfing species (Oryza sativa and O. glaberrima) have significantly lower equilibrium GC content compared with more outcrossing species. As recombination is less efficient in selfing species, these results suggest that recombination drives GC content. We also detected a positive relationship between expression levels and GC content in third codon positions, suggesting that selection favors codons ending with G or C bases. However, the correlation between GC content and recombination cannot be explained by selection on codon usage alone as it was also observed in noncoding positions. Finally, analyses of polymorphism data ruled out the hypothesis that genomic variation in GC content is due to mutational processes. Our results suggest that both gBGC and selection on codon usage affect GC content in the Oryza genus and likely in other grass species.     

Mutation of the Secys residue 266 in human type 2 selenodeiodinase alters 75Se incorporation without affecting its biochemical properties.

Type 2 deiodinase (D2) is a low Km iodothyronine deiodinase that metabolizes thyroxine (T4) to the active metabolite T3. We have recently shown that the cDNA for the human D2 coding region contains two in-frame selenocysteine (TGA) codons. The 3' TGA is seven codons 5' to a universal stop codon, TAA. The human D2 enzyme, transiently expressed in HEK-293 cells, can be in vivo labeled with 75Se as a doublet of approximately 31 kDa. This doublet is consistent with the possibility that the carboxy-terminal TGA codon can either encode selenocysteine or function as a stop codon. To test this hypothesis we mutagenized the second selenocysteine codon to a cysteine (TGC) or to an unambiguous stop codon (TAA). While the selenium incorporation pattern is different between the wild-type and mutant proteins, the deiodination properties of the enzyme are not affected by mutating the 3'TGA codon. Thus, we conclude that neither this residue nor the remaining seven carboxy-terminal amino acids are critical for the deiodination process.     

Tandem Stop Codons in Ciliates That Reassign Stop Codons

Tandem stop codons are extra stop codons hypothesized to be present downstream of genes to act as a backup in case of read-through of the real stop codon. Although seemingly absent from Escherichia coli, recent studies have confirmed the presence of such codons in yeast. In this paper we will analyze the genomes of two ciliate species—Paramecium tetraurelia and Tetrahymena thermophila—that reassign the stop codons TAA and TAG to glutamine, for the presence of tandem stop codons. We show that there are more tandem stop codons downstream of both Paramecium and Tetrahymena genes than expected by chance given the base composition of the downstream regions. This excess of tandem stop codons is larger in Tetrahymena and Paramecium than in yeast. We propose that this might be caused by a higher frequency of stop codon read-through in these species than in yeast, possibly because of a leaky termination machinery resulting from stop codon reassignment.     

Characterisation of a non-canonical genetic code in the oxymonad Streblomastix strix

The genetic code is one of the most highly conserved characters in living organisms. Only a small number of genomes have evolved slight variations on the code, and these non-canonical codes are instrumental in understanding the selective pressures maintaining the code. Here, we describe a new case of a non-canonical genetic code from the oxymonad flagellate Streblomastix strix. We have sequenced four protein-coding genes from S.strix and found that the canonical stop codons TAA and TAG encode the amino acid glutamine. These codons are retained in S.strix mRNAs, and the legitimate termination codons of all genes examined were found to be TGA, supporting the prediction that this should be the only true stop codon in this genome. Only four other lineages of eukaryotes are known to have evolved non-canonical nuclear genetic codes, and our phylogenetic analyses of alpha-tubulin, beta-tubulin, elongation factor-1 alpha (EF-1 alpha), heat-shock protein 90 (HSP90), and small subunit rRNA all confirm that the variant code in S.strix evolved independently of any other known variant. The independent origin of each of these codes is particularly interesting because the code found in S.strix, where TAA and TAG encode glutamine, has evolved in three of the four other nuclear lineages with variant codes, but this code has never evolved in a prokaryote or a prokaryote-derived organelle. The distribution of non-canonical codes is probably the result of a combination of differences in translation termination, tRNAs, and tRNA synthetases, such that the eukaryotic machinery preferentially allows changes involving TAA and TAG.     

Unusual mammalian usage of TGA stop codons reveals that sequence conservation need not imply purifying selection

The assumption that conservation of sequence implies the action of purifying selection is central to diverse methodologies to infer functional importance. GC-biased gene conversion (gBGC), a meiotic mismatch repair bias strongly favouring GC over AT, can in principle mimic the action of selection, this being thought to be especially important in mammals. As mutation is GC→AT biased, to demonstrate that gBGC does indeed cause false signals requires evidence that an AT-rich residue is selectively optimal compared to its more GC-rich allele, while showing also that the GC-rich alternative is conserved. We propose that mammalian stop codon evolution provides a robust test case. Although in most taxa TAA is the optimal stop codon, TGA is both abundant and conserved in mammalian genomes. We show that this mammalian exceptionalism is well explained by gBGC mimicking purifying selection and that TAA is the selectively optimal codon. Supportive of gBGC, we observe (i) TGA usage trends are consistent at the focal stop codon and elsewhere (in UTR sequences); (ii) that higher TGA usage and higher TAA→TGA substitution rates are predicted by a high recombination rate; and (iii) across species the difference in TAA <-> TGA substitution rates between GC-rich and GC-poor genes is largest in genomes that possess higher between-gene GC variation. TAA optimality is supported both by enrichment in highly expressed genes and trends associated with effective population size. High TGA usage and high TAA→TGA rates in mammals are thus consistent with gBGC’s predicted ability to “drive” deleterious mutations and supports the hypothesis that sequence conservation need not be indicative of purifying selection. A general trend for GC-rich trinucleotides to reside at frequencies far above their mutational equilibrium in high recombining domains supports the generality of these results.

Is sequence conservation a sign of purifying selection and hence functional importance? This analysis of why mammals use and conserve the most error-prone stop codon suggests not, consistent with GC-biased gene conversion’s predicted ability to “drive” deleterious mutations and supporting the hypothesis that sequence conservation need not be indicative of purifying selection.     

Relationships Among Stop Codon Usage Bias, Its Context, Isochores, and Gene Expression Level in Various Eukaryotes

It is well known that stop codons play a critical role in the process of protein synthesis. However, little effort has been made to investigate whether stop codon usage exhibits biases, such as widely seen for synonymous codon usage. Here we systematically investigate stop codon usage bias in various eukaryotes as well as its relationships with its context, GC3 content, gene expression level, and secondary structure. The results show that there is a strong bias for stop codon usage in different eukaryotes, i.e., UAA is overrepresented in the lower eukaryotes, UGA is overrepresented in the higher eukaryotes, and UAG is least used in all eukaryotes. Different conserved patterns for each stop codon in different eukaryotic classes are found based on information content and logo analysis. GC3 contents increase with increasing complexity of organisms. Secondary structure prediction revealed that UAA is generally associated with loop structures, whereas UGA is more uniformly present in loop and stem structures, i.e., UGA is less biased toward having a particular structure. The stop codon usage bias, however, shows no significant relationship with GC3 content and gene expression level in individual eukaryotes. The results indicate that genomic complexity and GC3 content might contribute to stop codon usage bias in different eukaryotes. Our results indicate that stop codons, like synonymous codons, exhibit biases in usage. Additional work will be needed to understand the causes of these biases and their relationship to the mechanism of protein termination. [Reviewing Editor: Dr. Manyuan Long]     

Vertebrate codon bias indicates a highly GC-rich ancestral genome

Two factors are thought to have contributed to the origin of codon usage bias in eukaryotes: 1) genome-wide mutational forces that shape overall GC-content and create context-dependent nucleotide bias, and 2) positive selection for codons that maximize efficient and accurate translation. Particularly in vertebrates, these two explanations contradict each other and cloud the origin of codon bias in the taxon. On the one hand, mutational forces fail to explain GC-richness (~ 60%) of third codon positions, given the GC-poor overall genomic composition among vertebrates (~ 40%). On the other hand, positive selection cannot easily explain strict regularities in codon preferences. Large-scale bioinformatic assessment, of nucleotide composition of coding and non-coding sequences in vertebrates and other taxa, suggests a simple possible resolution for this contradiction. Specifically, we propose that the last common vertebrate ancestor had a GC-rich genome (~ 65% GC). The data suggest that whole-genome mutational bias is the major driving force for generating codon bias. As the bias becomes prominent, it begins to affect translation and can result in positive selection for optimal codons. The positive selection can, in turn, significantly modulate codon preferences.     

GC content dependency of open reading frame prediction via stop codon frequencies

A frequently used approach for detecting potential coding regions is to search for stop codons. In the standard genetic code 3 out of 64 trinucleotides are stop codons. Hence, in random or non-coding DNA one can expect every 21st trinucleotide to have the same sequence as a stop codon. In contrast, the open reading frames (ORFs) of most protein-coding genes are considerably longer. Thus, the stop codon frequency in coding sequences deviates from the background frequency of the corresponding trinucleotides. This has been utilized for gene prediction, in particular, in detecting protein-coding ORFs. Traditional methods based on stop codon frequency are based on the assumption that the GC content is about 50%. However, many genomes show significant deviations from that value. With the presented method we can describe the effects of GC content on the selection of appropriate length thresholds of potentially coding ORFs. Conversely, for a given length threshold, we can calculate the probability of observing it in a random sequence. Thus, we can derive the maximum GC content for which ORF length is practicable as a feature for gene prediction methods and the resulting false positive rates. A rough estimate for an upper limit is a GC content of 80%. This estimate can be made more precise by including further parameters and by taking into account start codons as well. We demonstrate the feasibility of this method by applying it to the genomes of the bacteria Rickettsia prowazekii, Escherichia coli and Caulobacter crescentus, exemplifying the effect of GC content variations according to our predictions. We have adapted the method for predicting coding ORFs by stop codon frequency to the case of GC contents different from 50%. Usually, several methods for gene finding need to be combined. Thus, our results concern a specific part within a package of methods. Interestingly, for genomes with low GC content such as that of R. prowazekii, the presented method provides remarkably good results even when applied alone.     

Selection intensity on preferred codons correlates with overall codon usage bias in Caenorhabditis remanei

Adaptive codon usage provides evidence of natural selection in one of its most subtle forms: a fitness benefit of one synonymous codon relative to another. Codon usage bias is evident in the coding sequences of a broad array of taxa, reflecting selection for translational efficiency and/or accuracy as well as mutational biases. Here, we quantify the magnitude of selection acting on alternative codons in genes of the nematode Caenorhabditis remanei, an outcrossing relative of the model organism C. elegans, by fitting the expected mutation-selection-drift equilibrium frequency distribution of preferred and unpreferred codon variants to the empirical distribution. This method estimates the intensity of selection on synonymous codons in genes with high codon bias as N(e)s = 0.17, a value significantly greater than zero. In addition, we demonstrate for the first time that estimates of ongoing selection on codon usage among genes, inferred from nucleotide polymorphism data, correlate strongly with long-term patterns of codon usage bias, as measured by the frequency of optimal codons in a gene. From the pattern of polymorphisms in introns, we also infer that these findings do not result from the operation of biased gene conversion toward G or C nucleotides. We therefore conclude that coincident patterns of current and ancient selection are responsible for shaping biased codon usage in the C. remanei genome.     

Codon usage <Emphasis Type="BoldItalic">vis-a-vis</Emphasis> start and stop codon context analysis of three dicot species

To understand the variation in genomic composition and its effect on codon usage, we performed the comparative analysis of codon usage and nucleotide usage in the genes of three dicots, Glycine max, Arabidopsis thaliana and Medicago truncatula. The dicot genes were found to be A/T rich and have predominantly A-ending and/or T-ending codons. GC3s directly mimic the usage pattern of global GC content. Relative synonymous codon usage analysis suggests that the high usage frequency of A/T over G/C mononucleotide containing codons in AT-rich dicot genome is due to compositional constraint as a factor of codon usage bias. Odds ratio analysis identified the dinucleotides TpG, TpC, GpA, CpA and CpT as over-represented, where, CpG and TpA as under-represented dinucleotides. The results of (NcExp?NcObs)/NcExp plot suggests that selection pressure other than mutation played a significant role in influencing the pattern of codon usage in these dicots. PR2 analysis revealed the significant role of selection pressure on codon usage. Analysis of varience on codon usage at start and stop site showed variation in codon selection in these sites. This study provides evidence that the dicot genes were subjected to compositional selection pressure.     

Across bacterial phyla, distantly-related genomes with similar genomic GC content have similar patterns of amino acid usage

The GC content of bacterial genomes ranges from 16% to 75% and wide ranges of genomic GC content are observed within many bacterial phyla, including both gram negative and gram positive phyla. Thus, divergent genomic GC content has evolved repeatedly in widely separated bacterial taxa. Since genomic GC content influences codon usage, we examined codon usage patterns and predicted protein amino acid content as a function of genomic GC content within eight different phyla or classes of bacteria. We found that similar patterns of codon usage and protein amino acid content have evolved independently in all eight groups of bacteria. For example, in each group, use of amino acids encoded by GC-rich codons increased by approximately 1% for each 10% increase in genomic GC content, while the use of amino acids encoded by AT-rich codons decreased by a similar amount. This consistency within every phylum and class studied led us to conclude that GC content appears to be the primary determinant of the codon and amino acid usage patterns observed in bacterial genomes. These results also indicate that selection for translational efficiency of highly expressed genes is constrained by the genomic parameters associated with the GC content of the host genome.     

Codon usage bias covaries with expression breadth and the rate of synonymous evolution in humans, but this is not evidence for selection.

In numerous species, from bacteria to Drosophila, evidence suggests that selection acts even on synonymous codon usage: codon bias is greater in more abundantly expressed genes, the rate of synonymous evolution is lower in genes with greater codon bias, and there is consistency between genes in the same species in which codons are preferred. In contrast, in mammals, while nonequal use of alternative codons is observed, the bias is attributed to the background variance in nucleotide concentrations, reflected in the similar nucleotide composition of flanking noncoding and exonic third sites. However, a systematic examination of the covariants of codon usage controlling for background nucleotide content has yet to be performed. Here we present a new method to measure codon bias that corrects for background nucleotide content and apply this to 2396 human genes. Nearly all (99%) exhibit a higher amount of codon bias than expected by chance. The patterns associated with selectively driven codon bias are weakly recovered: Broadly expressed genes have a higher level of bias than do tissue-specific genes, the bias is higher for genes with lower rates of synonymous substitutions, and certain codons are repeatedly preferred. However, while these patterns are suggestive, the first two patterns appear to be methodological artifacts. The last pattern reflects in part biases in usage of nucleotide pairs. We conclude that we find no evidence for selection on codon usage in humans.     

The Usage of Codons Which are Similar to Stop Codons in the Genomes of Xylella fastidiosa and Xanthomonas citri

During the evolution of living organisms, a natural selection event occurs toward the optimization of their genomes regarding the usage of codons. During this process which is known as codon bias, a set of preferred codons is naturally defined in the genome of a given organism, since there are 61 possible codons (plus 3 stop codons) to 20 amino acids. Such event leads to optimization of metabolic cellular processes such as translational efficiency, RNA stability and energy saving. Although we know why, we do not know how exactly a set of preferred codons for each amino acid is defined for a given genome considering that the usage frequency of each synonymous codons is peculiar to each organism. In order to help answering this question, we analyzed the usage frequency of codons which are similar to stop codons, since a minor mutation on these codons may lead to a stop codon into the open reading frame compromising the protein expression as a result. We found a reduced use of those codons in Xanthomomas axonopodis pv. citri which presents an optimized genome regarding codon usage. On the other hand, such codons are more often used in Xylella fastidiosa, which does not seem to have established codon preferences as previously shown. Our results support that a set of preferred codons is not randomly selected and propose new ideas to the field warranting further experiments in this regard.     

Patterns of codon usage in two ciliates that reassign the genetic code: Tetrahymena thermophila and Paramecium tetraurelia

We used the recently sequenced genomes of the ciliates Tetrahymena thermophila and Paramecium tetraurelia to analyze the codon usage patterns in both organisms; we have analyzed codon usage bias, Gln codon usage, GC content and the nucleotide contexts of initiation and termination codons in Tetrahymena and Paramecium. We also studied how these trends change along the length of the genes and in a subset of highly expressed genes. Our results corroborate some of the trends previously described in Tetrahymena, but also negate some specific observations. In both genomes we found a strong bias toward codons with low GC content; however, in highly expressed genes this bias is smaller and codons ending in GC tend to be more frequent. We also found that codon bias increases along gene segments and in highly expressed genes and that the context surrounding initiation and termination codons are always AT rich. Our results also suggest differences in the efficiency of translation of the reassigned stop codons between the two species and between the reassigned codons. Finally, we discuss some of the possible causes for such translational efficiency differences.     

Use of silent mutations in cDNA encoding human glutathione transferase M2-2 for optimized expression in Escherichia coli.

Heterologous expression of human glutathione transferase M2-2 (GST M2-2) using Escherichia coli was improved 140-fold by mutating the cDNA expressing the enzyme. Expression of GST M2-2 from this cDNA clone, pKHXhGM2, generated approximately 190 mg protein per liter of bacterial culture, corresponding to approximately 12% of the total amount of soluble protein. The high-level-expressing cDNA was generated by oligonucleotide-directed mutagenesis introducing alternative silent mutations into the third nucleotide of codons 2, 4-7, and 10-14 in the 5' end of the cDNA coding region. The choice of alternative codons was restricted to those naturally occurring in highly biased genes in E. coli. Furthermore, the wild-type TAG stop codon at the 3' end was replaced with the two stop codons TAA and TGA in tandem to increase translation termination efficiency. The resulting partially randomized cDNA library was assayed for high-level expression using immunoscreening. Sequence similarities between the constructed high-level-expressing GST M2-2 cDNA and a similarly designed cDNA encoding the closely related human GST M1-1 suggest that the codons in the region immediately following the start codon are influential in achieving high-level expression. Pyrimidines seem to be more favorable than purines in the third position of codons in optimizing the expression of these enzymes in E. coli.