Tandem termination signals: myth or reality?

Stop codon usage as a window into genome evolution: mutation, selection, biased gene conversion and the TAG paradox

Protein coding genes terminate with one of three stop codons (TAA, TGA, or TAG) that, like synonymous codons, are not employed equally. With TGA and TAG having identical nucleotide content, analysis of their differential usage provides an unusual window into the forces operating on what are ostensibly functionally identical residues. Across genomes and between isochores within the human genome, TGA usage increases with G + C content but, with a common G + C → A + T mutation bias, this cannot be explained by mutation bias-drift equilibrium. Increased usage of TGA in G + C-rich genomes or genomic regions is also unlikely to reflect selection for the optimal stop codon, as TAA appears to be universally optimal, probably because it has the lowest read-through rate. Despite TAA being favored by selection and mutation bias, as with codon usage bias G + C pressure is the prime determinant of between-species TGA usage trends. In species with strong G + C-biased gene conversion (gBGC), such as mammals and birds, the high usage and conservation of TGA is best explained by an A + T → G + C repair bias. How to explain TGA enrichment in other G + C-rich genomes is less clear. Enigmatically, across bacterial and archaeal species and between human isochores TAG usage is mostly unresponsive to G + C pressure. This unresponsiveness we dub the TAG paradox as currently no mutational, selective, or gBGC model provides a well-supported explanation. That TAG does increase with G + C usage across eukaryotes makes the usage elsewhere yet more enigmatic. We suggest resolution of the TAG paradox may provide insights into either an unknown but common selective preference (probably at the DNA/RNA level) or an unrecognized complexity to the action of gBGC.

Why Is the UAG (Amber) Stop Codon Almost Absent in Highly Expressed Bacterial Genes?

The genome hypothesis postulates that genes in a genome tend to conform to their species’ usage of the codon catalog and the GC content of the DNA. Thus, codon frequencies differ across organisms, including the three termination codons in the standard genetic code. Here, we analyze the frequencies of stop codons in a group of highly expressed genes from 196 prokaryotes under strong translational selection. The occurrence of the three translation termination codons is highly biased, with UAA (ochre) being the most prevalent in almost all bacteria. In contrast, UAG (amber) is the least frequent termination codon, e.g., only 321 occurrences (7.4%) in E. coli K-12 substr. W3110. Of the 253 highly expressed genes, only two end with an UAG codon. The strength of the selective bias against UAG in highly expressed genes varies among bacterial genomes, but it is not affected by the GC content of these genomes. In contrast, increased GC content results in a decrease in UAA abundance with a concomitant increase in UGA abundance. We propose that readthrough efficiency and context effects could explain the prevalence of UAA over UAG, particularly in highly expressed genes. Findings from this communication can be utilized for the optimization of gene expression.

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.

Effective Population Size Predicts Local Rates but Not Local Mitigation of Read-through Errors

In correctly predicting that selection efficiency is positively correlated with the effective population size (Ne), the nearly neutral theory provides a coherent understanding of between-species variation in numerous genomic parameters, including heritable error (germline mutation) rates. Does the same theory also explain variation in phenotypic error rates and in abundance of error mitigation mechanisms? Translational read-through provides a model to investigate both issues as it is common, mostly nonadaptive, and has good proxy for rate (TAA being the least leaky stop codon) and potential error mitigation via "fail-safe" 3' additional stop codons (ASCs). Prior theory of translational read-through has suggested that when population sizes are high, weak selection for local mitigation can be effective thus predicting a positive correlation between ASC enrichment and Ne. Contra to prediction, we find that ASC enrichment is not correlated with Ne. ASC enrichment, although highly phylogenetically patchy, is, however, more common both in unicellular species and in genes expressed in unicellular modes in multicellular species. By contrast, Ne does positively correlate with TAA enrichment. These results imply that local phenotypic error rates, not local mitigation rates, are consistent with a drift barrier/nearly neutral model.

In eubacteria, unlike eukaryotes, there is no evidence for selection favouring fail-safe 3' additional stop codons

Errors throughout gene expression are likely deleterious, hence genomes are under selection to ameliorate their consequences. Additional stop codons (ASCs) are in-frame nonsense ‘codons’ downstream of the primary stop which may be read by translational machinery should the primary stop have been accidentally read through. Prior evidence in several eukaryotes suggests that ASCs are selected to prevent potentially-deleterious consequences of read-through. We extend this evidence showing that enrichment of ASCs is common but not universal for single cell eukaryotes. By contrast, there is limited evidence as to whether the same is true in other taxa. Here, we provide the first systematic test of the hypothesis that ASCs act as a fail-safe mechanism in eubacteria, a group with high read-through rates. Contra to the predictions of the hypothesis we find: there is paucity, not enrichment, of ASCs downstream; substitutions that degrade stops are more frequent in-frame than out-of-frame in 3’ sequence; highly expressed genes are no more likely to have ASCs than lowly expressed genes; usage of the leakiest primary stop (TGA) in highly expressed genes does not predict ASC enrichment even compared to usage of non-leaky stops (TAA) in lowly expressed genes, beyond downstream codon +1. Any effect at the codon immediately proximal to the primary stop can be accounted for by a preference for a T/U residue immediately following the stop, although if anything, TT- and TC- starting codons are preferred. We conclude that there is no compelling evidence for ASC selection in eubacteria. This presents an unusual case in which the same error could be solved by the same mechanism in eukaryotes and prokaryotes but is not. We discuss two possible explanations: that, owing to the absence of nonsense mediated decay, bacteria may solve read-through via gene truncation and in eukaryotes certain prion states cause raised read-through rates.

In all organisms, gene expression is error-prone. One such error, translational read-through, occurs where the primary stop codon of an expressed gene is missed by the translational machinery. Failure to terminate is likely to be costly, hence genomes are under selection to prevent this from happening. One proposed error-proofing strategy involves in-frame proximal additional stop codons (ASCs) which may act as a ‘fail-safe’ mechanism by providing another opportunity for translation to terminate. There is evidence for ASC enrichment in several eukaryotes. We extend this evidence showing it to be common but not universal in single celled eukaryotes. However, the situation in bacteria is poorly understood, despite bacteria having high read-through rates. Here, we test the fail-safe hypothesis within a broad range of bacteria. To our surprise, we find that not only are ASCs not enriched, but they may even be selected against. This provides evidence for an unusual circumstance where eukaryotes and prokaryotes could solve the same problem the same way but don’t. What are we to make of this? We suggest that if read-through is the problem, ASCs are not necessarily the expected solution. Owing to the absence of nonsense-mediated decay, a process that makes gene truncation in eukaryotes less viable, we propose bacteria may rescue a leaky stop by mutation that creates a new stop upstream. Alternatively, raised read-through rates in some particular conditions in eukaryotes might explain the difference.

Engineered transfer RNAs for suppression of premature termination codons

Premature termination codons (PTCs) are responsible for 10–15% of all inherited disease. PTC suppression during translation offers a promising approach to treat a variety of genetic disorders, yet small molecules that promote PTC read-through have yielded mixed performance in clinical trials. Here we present a high-throughput, cell-based assay to identify anticodon engineered transfer RNAs (ACE-tRNA) which can effectively suppress in-frame PTCs and faithfully encode their cognate amino acid. In total, we identify ACE-tRNA with a high degree of suppression activity targeting the most common human disease-causing nonsense codons. Genome-wide transcriptome ribosome profiling of cells expressing ACE-tRNA at levels which repair PTC indicate that there are limited interactions with translation termination codons. These ACE-tRNAs display high suppression potency in mammalian cells, Xenopus oocytes and mice in vivo, producing PTC repair in multiple genes, including disease causing mutations within cystic fibrosis transmembrane conductance regulator (CFTR).

Premature termination codon suppression therapy could be used to treat a range of genetic disorders. Here the authors present a high-throughput cell-based assay to identify anticodon engineered tRNAs with high suppression activity.

Re-exploration of the Codon Context Effect on Amber Codon-Guided Incorporation of Noncanonical Amino Acids in Escherichia coli by the Blue-White Screening Assay

The effect of codon context on amber codon-guided incorporation of noncanonical amino acids (NAAs) has been previously examined by antibiotic selection. Here, we re-explored this effect by screening a library in which three nucleotides upstream and downstream of the amber codon were randomised, and inserted within the lacZ-α gene. Thousands of clones were obtained and distinguished by the depth of blue colour upon exposure to X-gal. Large-scale sequencing revealed remarkable preferences in nucleotides downstream of the amber codon, and moderate preferences for upstream nucleotides. Nucleotide preference was quantified by a dual-luciferase assay, which verified that the optimum context for NAA incorporation, AATTAGACT, was applicable to different proteins. Our work provides a general guide for engineering amber codons into genes of interest in bacteria.

Comprehensive analysis of stop codon usage in bacteria and its correlation with release factor abundance

We present a comprehensive analysis of stop codon usage in bacteria by analyzing over eight million coding sequences of 4684 bacterial sequences. Using a newly developed program called "stop codon counter," the frequencies of the three classical stop codons TAA, TAG, and TGA were analyzed, and a publicly available stop codon database was built. Our analysis shows that with increasing genomic GC content the frequency of the TAA codon decreases and that of the TGA codon increases in a reciprocal manner. Interestingly, the release factor 1-specific codon TAG maintains a more or less uniform frequency (∼20%) irrespective of the GC content. The low abundance of TAG is also valid with respect to expression level of the genes ending with different stop codons. In contrast, the highly expressed genes predominantly end with TAA, ensuring termination with either of the two release factors. Using three model bacteria with different stop codon usage (Escherichia coli, Mycobacterium smegmatis, and Bacillus subtilis), we show that the frequency of TAG and TGA codons correlates well with the relative steady state amount of mRNA and protein for release factors RF1 and RF2 during exponential growth. Furthermore, using available microarray data for gene expression, we show that in both fast growing and contrasting biofilm formation conditions, the relative level of RF1 is nicely correlated with the expression level of the genes ending with TAG.

Stops making sense: translational trade-offs and stop codon reassignment

Background

Efficient gene expression involves a trade-off between (i) premature termination of protein synthesis; and (ii) readthrough, where the ribosome fails to dissociate at the terminal stop. Sense codons that are similar in sequence to stop codons are more susceptible to nonsense mutation, and are also likely to be more susceptible to transcriptional or translational errors causing premature termination. We therefore expect this trade-off to be influenced by the number of stop codons in the genetic code. Although genetic codes are highly constrained, stop codon number appears to be their most volatile feature.

Results

In the human genome, codons readily mutable to stops are underrepresented in coding sequences. We construct a simple mathematical model based on the relative likelihoods of premature termination and readthrough. When readthrough occurs, the resultant protein has a tail of amino acid residues incorrectly added to the C-terminus. Our results depend strongly on the number of stop codons in the genetic code. When the code has more stop codons, premature termination is relatively more likely, particularly for longer genes. When the code has fewer stop codons, the length of the tail added by readthrough will, on average, be longer, and thus more deleterious. Comparative analysis of taxa with a range of stop codon numbers suggests that genomes whose code includes more stop codons have shorter coding sequences.

Conclusions

We suggest that the differing trade-offs presented by alternative genetic codes may result in differences in genome structure. More speculatively, multiple stop codons may mitigate readthrough, counteracting the disadvantage of a higher rate of nonsense mutation. This could help explain the puzzling overrepresentation of stop codons in the canonical genetic code and most variants.

Tandem termination signal in plant mRNAs

It was proposed that if some mRNA characteristics resulted in a low efficiency of termination signal, an additional closely located stop codon (tandem stop codons) could be used to prevent the harmful readthrough. However, the role of tandem terminators in higher eukaryotes was not verified and remains hypothetical. In this work the sequence features of Arabidopsis thaliana and Oryza sativa mRNAs were analyzed. It was found that plant mRNAs with UGA terminator were characterized by a higher frequency of nonsense codons in the first triplet position of 3'-UTR that could result from a weak natural selection for "reserve" stop signal. Interestingly, the presence of tandem stop codons positively correlated with a specific amino acid composition in the C-terminal position of the encoded proteins. In particular, C-terminal glycine positively correlated with significantly higher frequencies of reserve terminators at the beginning positions of 3'-UTR in UGA-containing mRNAs. This finding coincides with some earlier observations concerning the role of glycine and its codons in inefficient termination of translation and recoding (e.g., 2A oligopeptide).

Evolution of prokaryotic genes by shift of stop codons

De novo origin of coding sequence remains an obscure issue in molecular evolution. One of the possible paths for addition (subtraction) of DNA segments to (from) a gene is stop codon shift. Single nucleotide substitutions can destroy the existing stop codon, leading to uninterrupted translation up to the next stop codon in the gene's reading frame, or create a premature stop codon via a nonsense mutation. Furthermore, short indels-caused frameshifts near gene's end may lead to premature stop codons or to translation past the existing stop codon. Here, we describe the evolution of the length of coding sequence of prokaryotic genes by change of positions of stop codons. We observed cases of addition of regions of 3'UTR to genes due to mutations at the existing stop codon, and cases of subtraction of C-terminal coding segments due to nonsense mutations upstream of the stop codon. Many of the observed stop codon shifts cannot be attributed to sequencing errors or rare deleterious variants segregating within bacterial populations. The additions of regions of 3'UTR tend to occur in those genes in which they are facilitated by nearby downstream in-frame triplets which may serve as new stop codons. Conversely, subtractions of coding sequence often give rise to in-frame stop codons located nearby. The amino acid composition of the added region is significantly biased, compared to the overall amino acid composition of the genes. Our results show that in prokaryotes, shift of stop codon is an underappreciated contributor to functional evolution of gene length.

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.

Comparison of characteristics and function of translation termination signals between and within prokaryotic and eukaryotic organisms

Six diverse prokaryotic and five eukaryotic genomes were compared to deduce whether the protein synthesis termination signal has common determinants within and across both kingdoms. Four of the six prokaryotic and all of the eukaryotic genomes investigated demonstrated a similar pattern of nucleotide bias both 5′ and 3′ of the stop codon. A preferred core signal of 4 nt was evident, encompassing the stop codon and the following nucleotide. Codons decoded by hyper-modified tRNAs were over-represented in the region 5′ to the stop codon in genes from both kingdoms. The origin of the 3′ bias was more variable particularly among the prokaryotic organisms. In both kingdoms, genes with the highest expression index exhibited a strong bias but genes with the lowest expression showed none. Absence of bias in parasitic prokaryotes may reflect an absence of pressure to evolve more efficient translation. Experiments were undertaken to determine if a correlation existed between bias in signal abundance and termination efficiency. In Escherichia coli signal abundance correlated with termination efficiency for UAA and UGA stop codons, but not in mammalian cells. Termination signals that were highly inefficient could be made more efficient by increasing the concentration of the cognate decoding release factor.

Transfer-messenger RNA unfolds as it transits the ribosome

In bacteria, translation of mRNAs lacking stop codons produces truncated polypeptides and traps ribosomes in unproductive complexes. Potentially harmful truncated proteins are tagged with short peptides encoded by the mRNA-like domain of tmRNA and targeted for digestion by housekeeping proteases. We show that altered Escherichia coli transfer-messenger RNAs (tmRNAs) produce in vivo fusion proteins with peptide tags that extend far beyond the conventional termination signal of the wild-type tmRNA. Regions of tmRNA capable of serving as templates for protein synthesis include helix 5, as well as pseudoknots 2, 3, and 4. The removal of all six in-frame stop codons negatively affects tmRNA processing, thereby preventing translation of the 3' portion of the tRNA-like domain. These findings provide evidence that trans-translation can be accompanied by the unfolding of significant portions of the tmRNA molecule. Many of these conformational changes are likely to be required during trans-translation to maintain the ribosomal subunits in close proximity to the tmRNA for monitoring its transit.

Conservation of tandem stop codons in yeasts

This study shows that a statistical excess of stop codons has evolved at the third codon downstream of the real stop codon UAA in yeasts. Comparative analysis indicates that stop codons at this location are considerably more conserved than sense codons, suggesting that these tandem stop codons are maintained by selection.

Background

It has been long thought that the stop codon in a gene is followed by another stop codon that acts as a backup if the real one is read through by a near-cognate tRNA. The existence of such 'tandem stop codons', however, remains elusive.

Results

Here we show that a statistical excess of stop codons has evolved at the third codon downstream of the real stop codon UAA in yeasts. Comparative analysis indicates that stop codons at this location are considerably more conserved than sense codons, suggesting that these tandem stop codons are maintained by selection. We evaluated the influence of expression levels of genes and other biological factors on the distribution of tandem stop codons. Our results suggest that expression level is an important factor influencing the presence of tandem stop codons.

Conclusions

Our study demonstrates the existence of tandem stop codons, which represent one of many meaningful genomic features that are driven by relatively weak selective forces.

Genome-wide prediction of stop codon readthrough during translation in the yeast Saccharomyces cerevisiae

In-frame stop codons normally signal termination during mRNA translation, but they can be read as 'sense' (readthrough) depending on their context, comprising the 6 nt preceding and following the stop codon. To identify novel contexts directing readthrough, under-represented 5' and 3' stop codon contexts from Saccharomyces cerevisiae were identified by genome-wide survey in silico. In contrast with the nucleotide bias 3' of the stop codon, codon bias in the two codon positions 5' of the termination codon showed no correlation with known effects on stop codon readthrough. However, individually, poor 5' and 3' context elements were equally as effective in promoting stop codon readthrough in vivo, readthrough which in both cases responded identically to changes in release factor concentration. A novel method analysing specific nucleotide combinations in the 3' context region revealed positions +1,2,3,5 and +1,2,3,6 after the stop codon were most predictive of termination efficiency. Downstream of yeast open reading frames (ORFs), further in-frame stop codons were significantly over-represented at the +1, +2 and +3 codon positions after the ORF, acting to limit readthrough. Thus selection against stop codon readthrough is a dominant force acting on 3', but not on 5', nucleotides, with detectable selection on nucleotides as far downstream as +6 nucleotides. The approaches described can be employed to define potential readthrough contexts for any genome.

Reinitiation and Recycling are Distinct Processes Occurring Downstream of Translation Termination in Yeast

The circularisation model of the polysome suggests that ribosome recycling is facilitated by 5'-3' interactions mediated by the cap-binding complex eIF4F and the poly(A)-binding protein, Pab1. Alternatively, downstream of a short upstream open reading frame (uORF) in the 5' untranslated region of a gene, posttermination ribosomes can maintain the competence to (re)initiate translation. Our data show that recycling and reinitiation must be distinct processes in Saccharomyces cerevisiae. The role of the 3'UTR in recycling was assessed by restricting ribosome movement along the mRNA using a poly(G) stretch or the mammalian iron regulatory protein bound to the iron responsive element. We find that although 3'UTR structure can influence translation, the main pathway of ribosome recycling does not depend on scanning-like movement through the 3'UTR. Changes in termination kinetics or disruption of the Pab1-eIF4F interaction do not affect recycling, yet the maintenance of normal in vivo mRNP structure is important to this process. Using bicistronic ACT1-LUC constructs, elongating yeast ribosomes were found to maintain the competence to (re)initiate over only short distances. Thus, as the first ORF to be translated is progressively truncated, reinitiation downstream of an uORF of 105nt is found to be just detectable, and increases markedly in efficiency as uORF length is reduced to 15nt. Experiments using a strain mutated in the Cca1 nucleotidyltransferase suggest that the uORF length-dependence of changes in reinitiation competence is affected by peptide elongation kinetics, but that ORF length per se may also be relevant.

Polyamines enhance synthesis of the RNA polymerase sigma 38 subunit by suppression of an amber termination codon in the open reading frame

The mechanisms by which polyamines stimulate synthesis of the RNA polymerase sigma(38) subunit in Escherichia coli were studied. Polyamine stimulation was observed only in strains in which the 33rd codon of RpoS mRNA is a UAG termination codon instead of a CAG codon for glutamine in wild-type E. coli. Readthrough of the termination codon by Gln-tRNA(supE) was stimulated by polyamines. This stimulation was found to be caused by an increase in both the level of suppressor tRNA(supE) and the binding affinity of Gln-tRNA(supE) for ribosomes. The stimulatory effect was observed with a UAG termination codon but not with UGA and UAA codons. Readthrough of the UAG termination codon at the 270th amino acid position of RpoS mRNA was also stimulated by polyamines, indicating that polyamines stimulate readthrough of a UAG codon regardless of its location within the RpoS mRNA. When cell viability of an E. coli strain having a termination codon in the 33rd position of RpoS mRNA was compared using cells cultured with or without putrescine, it was higher in cells cultured with putrescine than in cells cultured without putrescine. The level of sigma(38) subunit in the cells cultured with putrescine was higher than that in cells cultured without putrescine on days 2, 4, and 8, but the level of sigma(70) subunit was almost the same in cells cultured with or without putrescine. These results confirm that elevated expression of the rpoS gene is important for cell viability at late stationary phase.