Terminating eukaryote translation: domain 1 of release factor eRF1 functions in stop codon recognition

Eukaryote ribosomal translation is terminated when release factor eRF1, in a complex with eRF3, binds to one of the three stop codons. The tertiary structure and dimensions of eRF1 are similar to that of a tRNA, supporting the hypothesis that release factors may act as molecular mimics of tRNAs. To identify the yeast eRF1 stop codon recognition domain (analogous to a tRNA anticodon), a genetic screen was performed to select for mutants with disabled recognition of only one of the three stop codons. Nine out of ten mutations isolated map to conserved residues within the eRF1 N-terminal domain 1. A subset of these mutants, although wild-type for ribosome and eRF3 interaction, differ in their respective abilities to recognize each of the three stop codons, indicating codon-specific discrimination defects. Five of six of these stop codon-specific mutants define yeast domain 1 residues (I32, M48, V68, L123, and H129) that locate at three pockets on the eRF1 domain 1 molecular surface into which a stop codon can be modeled. The genetic screen results and the mutant phenotypes are therefore consistent with a role for domain 1 in stop codon recognition; the topology of this eRF1 domain, together with eRF1-stop codon complex modeling further supports the proposal that this domain may represent the site of stop codon binding itself.     

KRAS codon 12 mutations occur very frequently in pancreatic adenocarcinomas.

DNAs from human pancreatic adenocarcinomas were analyzed for the presence of mutations in codons 12, 13 and 61 of the NRAS, KRAS and HRAS gene. Formalin-fixed and paraffin-embedded tissue was used directly in an in vitro amplification reaction to expand the relevant RAS sequences. The mutations were detected by selective hybridization using mutation-specific synthetic oligonucleotides. In 28 of the 30 patients we found a mutation in codon 12 of the KRAS gene. This result confirms the findings of Almoguera et al. [Cell 53 (1988) 549-554] that KRAS mutations occur frequently in adenocarcinomas of the exocrine pancreas. The mutations are predominantly G-T transversions, in contrast to the KRAS mutations in colon tumors which are mainly G-A transitions. Furthermore, in a portion of the tumors the mutation appears to be homozygous.     

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.     

Identification of eRF1 residues that play critical and complementary roles in stop codon recognition

The initiation and elongation stages of translation are directed by codon-anticodon interactions. In contrast, a release factor protein mediates stop codon recognition prior to polypeptide chain release. Previous studies have identified specific regions of eukaryotic release factor one (eRF1) that are important for decoding each stop codon. The cavity model for eukaryotic stop codon recognition suggests that three binding pockets/cavities located on the surface of eRF1's domain one are key elements in stop codon recognition. Thus, the model predicts that amino acid changes in or near these cavities should influence termination in a stop codon-dependent manner. Previous studies have suggested that the TASNIKS and YCF motifs within eRF1 domain one play important roles in stop codon recognition. These motifs are highly conserved in standard code organisms that use UAA, UAG, and UGA as stop codons, but are more divergent in variant code organisms that have reassigned a subset of stop codons to sense codons. In the current study, we separately introduced TASNIKS and YCF motifs from six variant code organisms into eRF1 of Saccharomyces cerevisiae to determine their effect on stop codon recognition in vivo. We also examined the consequences of additional changes at residues located between the TASNIKS and YCF motifs. Overall, our results indicate that changes near cavities two and three frequently mediated significant effects on stop codon selectivity. In particular, changes in the YCF motif, rather than the TASNIKS motif, correlated most consistently with variant code stop codon selectivity.     

Distinct paths to stop codon reassignment by the variant-code organisms Tetrahymena and Euplotes

The reassignment of stop codons is common among many ciliate species. For example, Tetrahymena species recognize only UGA as a stop codon, while Euplotes species recognize only UAA and UAG as stop codons. Recent studies have shown that domain 1 of the translation termination factor eRF1 mediates stop codon recognition. While it is commonly assumed that changes in domain 1 of ciliate eRF1s are responsible for altered stop codon recognition, this has never been demonstrated in vivo. To carry out such an analysis, we made hybrid proteins that contained eRF1 domain 1 from either Tetrahymena thermophila or Euplotes octocarinatus fused to eRF1 domains 2 and 3 from Saccharomyces cerevisiae. We found that the Tetrahymena hybrid eRF1 efficiently terminated at all three stop codons when expressed in yeast cells, indicating that domain 1 is not the sole determinant of stop codon recognition in Tetrahymena species. In contrast, the Euplotes hybrid facilitated efficient translation termination at UAA and UAG codons but not at the UGA codon. Together, these results indicate that while domain 1 facilitates stop codon recognition, other factors can influence this process. Our findings also indicate that these two ciliate species used distinct approaches to diverge from the universal genetic code.     

Optimal codon identities in bacteria: implications from the conflicting results of two different methods

A correlation method was recently adopted to identify selection-favored 'optimal' codons from 675 bacterial genomes. Surprisingly, the identities of these optimal codons were found to track the bacterial GC content, leading to a conclusion that selection would generally shape the codon usages to the same direction as the overall mutation does. Raising several concerns, here we report a thorough comparative study on 203 well-selected bacterial species, which strongly suggest that the previous conclusion is likely an illusion. Firstly, the previous study did not preclude species that are suffering weak or no selection pressures on their codon usages. For these species, as showed in this study, the optimal codon identities are prone to be incorrect and follow GC content. Secondly, the previous study only adopted the correlation method, without considering another method to test the reliability of inferred optimal codons. Actually by definition, optimal codons can also be identified by simply comparing codon usages between high- and low-expression genes. After using both methods to identify optimal codons for the selected species, we obtained highly conflicting results, suggesting at least one method is misleading. Further we found a critical problem of correlation method at the step of calculating gene bias level. Due to a failure of accurately defining the background mutation, the problem would result in wrong optimal codon identities. In other words, partial mutational effects on codon choices were mistakenly regarded as selective influences, leading to incorrect and biased optimal codon identities. Finally, considering the translational dynamics, optimal codons identified by comparison method can be well-explained by tRNA compositions, whereas optimal codons identified by correlation method can not be. For all above reasons, we conclude that real optimal codons actually do not track the genomic GC content, and correlation method is misleading in identifying optimal codons and better be avoided.     

Stop codon recognition in ciliates: Euplotes release factor does not respond to reassigned UGA codon

In eukaryotes, the polypeptide release factor 1 (eRF1) is involved in translation termination at all three stop codons. However, the mechanism for decoding stop codons remains unknown. A direct interaction of eRF1 with the stop codons has been postulated. Recent studies focus on eRF1 from ciliates in which some stop codons are reassigned to sense codons. Using an in vitro assay based on mammalian ribosomes, we show that eRF1 from the ciliate Euplotes aediculatus responds to UAA and UAG as stop codons and lacks the capacity to decipher the UGA codon, which encodes cysteine in this organism. This result strongly suggests that in ciliates with variant genetic codes eRF1 does not recognize the reassigned codons. Recent hypotheses describing stop codon discrimination by eRF1 are not fully consistent with the set of eRF1 sequences available so far and require direct experimental testing.     

Likelihood models of somatic mutation and codon substitution in cancer genes

The role of somatic mutation in cancer is well established and several genes have been identified that are frequent targets. This has enabled large-scale screening studies of the spectrum of somatic mutations in cancers of particular organs. Cancer gene mutation databases compile the results of many studies and can provide insight into the importance of specific amino acid sequences and functional domains in cancer, as well as elucidate aspects of the mutation process. Past studies of the spectrum of cancer mutations (in particular genes) have examined overall frequencies of mutation (at specific nucleotides) and of missense, nonsense, and silent substitution (at specific codons) both in the sequence as a whole and in a specific functional domain. Existing methods ignore features of the genetic code that allow some codons to mutate to missense, or stop, codons more readily than others (i.e., by one nucleotide change, vs. two or three). A new codon-based method to estimate the relative rate of substitution (fixation of a somatic mutation in a cancer cell lineage) of nonsense vs. missense mutations in different functional domains and in different tumor tissues is presented. Models that account for several potential influences on rates of somatic mutation and substitution in cancer progenitor cells and allow biases of mutation rates for particular dinucleotide sequences (CGs and dipyrimidines), transition vs. transversion bias, and variable rates of silent substitution across functional domains (useful in detecting investigator sampling bias) are considered. Likelihood-ratio tests are used to choose among models, using cancer gene mutation data. The method is applied to analyze published data on the spectrum of p53 mutations in cancers. A novel finding is that the ratio of the probability of nonsense to missense substitution is much lower in the DNA-binding and transactivation domains (ratios near 1) than in structural domains such as the linker, tetramerization (oligomerization), and proline-rich domains (ratios exceeding 100 in some tissues), implying that the specific amino acid sequence may be less critical in structural domains (e.g., amino acid changes less often lead to cancer). The transition vs. transversion bias and effect of CpG dinucleotides on mutation rates in p53 varied greatly across cancers of different organs, likely reflecting effects of different endogenous and exogenous factors influencing mutation in specific organs.     

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.     

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.     

An improved method for determining codon variability in a gene and its application to the rate of fixation of mutations in evolution

If one has the amino acid sequences of a set of homologous proteins as well as their phylogenetic relationships, one can easily determine the minimum number of mutations (nucleotide replacements) which must have been fixed in each codon since their common ancestor. It is found that for 29 species of cytochrome c the data fit the assumption that there is a group of approximately 32 invariant codons and that the remainder compose two Poisson-distributed groups of size 65 and 16 codons, the latter smaller group fixing mutations at about 3.2 times the rate of the larger. It is further found that the size of the invariant group increases as the range of species is narrowed. Extrapolation suggests that less than 10% of the codons in a given mammalian cytochrome c gene are capable of accepting a mutation. This is consistent with the view that at any one point in time only a very restricted number of positions can fix mutations but that as mutations are fixed the positions capable of accepting mutations also change so that examination of a wide range of species reveals a wide range of altered positions. We define this restricted group as the concomitantly variable codons. Given this restriction, the fixation rates for mutations in concomitantly variable codons in cytochrome c and fibrinopeptide A are not very different, a result which should be the case if most of these mutations are in fact selectively neutral as Kimura suggests.Paper number 1382 from the Laboratory of Genetics. Work performed in part at the University of Iowa, Department of Preventive Medicine and Environmental Health and Department of Statistics, Iowa City, Iowa. Computing supported by the Graduate College, University of Iowa.     

An estimate on the effect of point mutation and natural selection on the rate of amino acid replacement in proteins

We outline a method for estimating quantitatively the influence of point mutations and selection on the frequencies of codons and amino acids. We show how the mutation rate, i.e., the rate of amino acid replacement due to point mutation, can be affected by the codon usage as well as by the rates of the involved base exchanges. A comparison of the mutation rates calculated from reliable values of codon usage and base exchange probabilities with those that would be expected on the basis of chance reveals a notable suppression of replacements leading to tryptophan, glutamate, lysine, and methionine, and particularly of those leading to the termination codons. If selection constraints are neglected and only mutations are taken into account, the best agreement between expected and observed frequencies of both codons and amino acids is obtained for alpha = 1.13-1.15, where (Formula: see text). The "selection values" of codons and amino acids derived by our method show a pattern that partially deviates from others in the literature. For example, the selection pressure on methionine and cysteine turns out to be much more pronounced than expected if only the discrepancies between their observed and expected occurrences in proteins are considered. To estimate to what extent randomly occurring amino acid replacements are accepted by selection, we constructed an "acceptability matrix" from the well-established matrix of accepted point mutations. On the basis of this matrix "acceptability values" of the amino acids can be defined that correlate with their selection values. We also examine the significance of mutations and selection of amino acids with respect to their physicochemical properties and functions in proteins. The conservatism of amino acid replacements with respect to certain properties such as polarity can be brought about by the mutational process alone, whereas the conservatism with respect to other relevant properties--among them all measures of bulkiness--obviously is the result of additional selectional constraints on the evolution of protein structures.     

The involvement of base 1054 in 16S rRNA for UGA stop codon dependent translational termination.

The deletion of the highly conserved cytidine nucleotide at position 1054 in E. coli 16S rRNA has been characterized to confer an UGA stop codon specific suppression activity which suggested a functional participation of small subunit rRNA in translational termination. Based on this structure-function correlation we constructed the three point mutations at site 1054, changing the wild-type C residue to an A, G or U base. The mutations were expressed from a complete plasmid encoded rRNA operon (rrnB) using a conditional expression system with the lambda PL-promoter. All three altered 16S rRNA molecules were expressed and incorporated into 70S ribosomal particles. Structural analysis of the protein and 16S rRNA moieties of the mutant ribosomes showed no differences when compared to wild-type particles. The phenotypic analysis revealed that only the 1054G base change led to a significantly reduced generation time of transformed cells, which could be correlated with the inability of the mutant ribosomes to specifically stop at UGA stop codons in vivo. The response towards UAA and UAG termination codons was not altered. Furthermore, in vitro RF-2 termination factor binding experiments indicated that the association behaviour of mutant ribosomes was not changed, enforcing the view that the UGA stop codon suppression is a direct consequence of the rRNA mutation. Taken together, these results argue for a direct participation of that 16S rRNA motif in UGA dependent translational termination and furthermore, suggest that termination factor binding and stop codon recognition are two separate steps of the termination event.     

Synonymous codon usage in Cryptosporidium parvum: identification of two distinct trends among genes

The usage of alternative synonymous codons in the apicomplexan Cryptosporidium parvum has been investigated. A data set of 54 genes was analysed. Overall, A- and U-ending codons predominate, as expected in an A+T-rich genome. Two trends of codon usage variation among genes were identified using correspondence analysis. The primary trend is in the extent of usage of a subset of presumably translationally optimal codons, that are used at significantly higher frequencies in genes expected to be expressed at high levels. Fifteen of the 18 codons identified as optimal are more G+C-rich than the otherwise common codons, so that codon selection associated with translation opposes the general mutation bias. Among 40 genes with lower frequencies of these optimal codons, a secondary trend in G+C content was identified. In these genes, G+C content at synonymously variable third positions of codons is correlated with that in 5' and 3' flanking sequences, indicative of regional variation in G+C content, perhaps reflecting regional variation in mutational biases.     

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.     

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.     

The effect of context on synonymous codon usage in genes with low codon usage bias.

The effect of neighbouring bases on the usage of synonymous codons in genes with low codon usage bias in yeast and E. coli is examined. The codon adaptation index is employed to identify a group of genes in each organism with low codon usage bias, which are likely to be weakly expressed. A similar pattern is found in complementary sequences with respect to synonymous usage of A vs G or of U vs C. It is suggested that this may reflect an effect of context on mutation rates in weakly expressed genes.     

Stop codons in bacteria are not selectively equivalent

ABSTRACT: BACKGROUND: The evolution and genomic stop codon frequencies have not been rigorously studied with the exception of coding of non-canonical amino acids. Here we study the rate of evolution and frequency distribution of stop codons in bacterial genomes. RESULTS: We show that in bacteria stop codons evolve slower than synonymous sites, suggesting the action of weak negative selection. However, the frequency of stop codons relative to genomic nucleotide content indicated that this selection regime is not straightforward. The frequency of TAA and TGA stop codons is GC-content dependent, with TAA decreasing and TGA increasing with GC-content, while TAG frequency is independent of GC-content. Applying a formal, analytical model to these data we found that the relationship between stop codon frequencies and nucleotide content cannot be explained by mutational biases or selection on nucleotide content. However, with weak nucleotide content-dependent selection on TAG, -0.5 < Nes < 1.5, the model fits all of the data and recapitulates the relationship between TAG and nucleotide content. For biologically plausible rates of mutations we show that, in bacteria, TAG stop codon is universally associated with lower fitness, with TAA being the optimal for G-content < 16% while for G-content > 16% TGA has a higher fitness than TAG. CONCLUSIONS: Our data indicate that TAG codon is universally suboptimal in the bacterial lineage, such that TAA is likely to be the preferred stop codon for low GC content while the TGA is the preferred stop codon for high GC content. The optimization of stop codon usage may therefore be useful in genome engineering or gene expression optimization applications.ReviewersThis article was reviewed by Michail Gelfand, Arcady Mushegian and Shamil Sunyaev. For the full reviews, please go to the Reviewers' Comments section.     

Combined point mutations in codon 12 and 13 of KRAS oncogene in prostate carcinomas

Prostate cancer is a common malignancy that develops by structural mutation(s) and/or other genetic alterations in specific genes.The G to T transversions in codon 12 and C to T transitions in codon 13 of KRAS proto-oncogene are predominant point mutations that occur in about 20% of different cancers in human. In the current study it was aimed to investigate the prevalence and predictive significance of KRAS mutations in patients with prostate carcinomas. In a total of 30 fresh tumoural tissue specimens were investigated in patients with prostate carcinoma. All tumoural specimens were histo-pathologically diagnosed and genotyped for codon 12, 13 KRAS point mutations by reverse hybridisation and direct sequencing methods. KRAS mutations were found in 12 (40%) samples with 29 samples deriving from adenocarcinomas and 1 sample was small cell prostate carcinoma. In 1 (3.44%) sample codon 12 was found to be mutated and in 2 (6.8%) samples codon 13 and in 9 (31%) samples combined codon 12 and 13 were found to be mutated particularly in higher grade of tumoural tissues. Our study, based on representative collection of human prostate tumours, indicates that combined mutations in codons 12 and 13 KRAS are relatively infrequent and most commonly occur in prostate carcinomas.