Misreading of termination codons in eukaryotes by natural nonsense suppressor tRNAs

Translational stop codon readthrough provides a regulatory mechanism of gene expression that is extensively utilised by positive-sense ssRNA viruses. The misreading of termination codons is achieved by a variety of naturally occurring suppressor tRNAs whose structure and function is the subject of this survey. All of the nonsense suppressors characterised to date (with the exception of selenocysteine tRNA) are normal cellular tRNAs that are primarily needed for reading their cognate sense codons. As a consequence, recognition of stop codons by natural suppressor tRNAs necessitates unconventional base pairings in anticodon–codon interactions. A number of intrinsic features of the suppressor tRNA contributes to the ability to read non-cognate codons. Apart from anticodon–codon affinity, the extent of base modifications within or 3′ of the anticodon may up- or down-regulate the efficiency of suppression. In order to out-compete the polypeptide chain release factor an absolute prerequisite for the action of natural suppressor tRNAs is a suitable nucleotide context, preferentially at the 3′ side of the suppressed stop codon. Three major types of viral readthrough sites, based on similar sequences neighbouring the leaky stop codon, can be defined. It is discussed that not only RNA viruses, but also the eukaryotic host organism might gain some profit from cellular suppressor tRNAs.     

Phylogeny of rodentia (Mammalia) inferred from the nuclear-encoded gene IRBP

The order Rodentia includes nearly half of all living mammalian species. Phylogenetic relationships among 22 species of rodents were investigated by use of a 1.2-kb region from exon 1 of the single-copy nuclear gene IRBP. IRBP has been extensively used for study of interordinal phylogeny in mammals, which allowed inclusion of 50 outgroup species, representing every eutherian order plus seven marsupials. Several clades were strongly supported, regardless of analytical method or inclusion/exclusion of data. These include a monophyletic Muroidea, with a clade including Spalax and Rhizomys as the first divergence; a clade uniting Zapus with Dipus, but excluding Sicista; a monophyletic Myodonta (Muroidea plus Dipodidae); and a clade including Aplodontidae as sister to Sciuridae. One bipartition, separating Hystricognathi and Geomyoidea from the remaining rodents, is strongly supported in all analyses that include third-position sites but almost completely absent from analyses that exclude third-position sites. A combination of nonstationary nucleotide composition and branch length effects may be causing all methods examined (including those using the LogDet distance) to support an incorrect conclusion when third-position sites are analyzed together with first- and second-position sites.     

Identifying conflicting signal in a multigene analysis reveals a highly resolved tree: the phylogeny of Rodentia (Mammalia)

Homoplasy among morphological characters has hindered inference of higher level rodent phylogeny for over 100 years. Initial molecular studies, based primarily on single genes, likewise produced little resolution of the deep relationships among rodent families. Two recent molecular studies (Huchon et al., 2002, Mol. Biol. Evol. 19:1053-1065; Adkins et al., 2003, Mol. Phylogenet. Evol. 26:409-420), using larger samples from the nuclear genome, have produced phylogenies that are generally concordant with each other, but many of the deep superfamilial nodes were still lacking substantial statistical support. Data are presented here for a total of approximately 3,600 base pairs from portions of three different nuclear protein-coding genes, CB1, IRBP, and RAG2, from 19 rodents and three outgroups. Separate analyses, with data partitioned according to both genes and codon position, produced conflicting results. Trees obtained from all partitions of CB1 and RAG2 and those obtained from the first- plus second-position sites of IRBP were generally concordant with each other and the trees from the two recent studies, whereas trees obtained from the third-position sites of IRBP were not. Although the IRBP third-position sites represent only 1/9 of the total data set, combined analyses using either parsimony or likelihood resulted in trees in agreement with the IRBP third-position sites and in disagreement with the remaining 8/9 of the sites from this data set and the two recent multigene studies. In contrast, maximum-likelihood analysis using a site-specific rates model did recover a tree that is highly congruent with the trees in the two recent studies. If the IRBP third-position sites are removed from the current data set, then combined likelihood analyses obtain a tree that is highly congruent with those of the two recent studies. This analysis also provides, for the first time in a study of rodent phylogeny, robust statistical support for every bipartition, with just one exception. This tree divides rodents into two major clades. The first contains Myodonta (Muroidea plus Dipodidae) and the only unresolved trichotomy, from which descend Geomyoidea, Pedetidae, and Castoridae. On the other side of the root is a clade containing Sciuroidea plus Gliridae, and Hystricognathi. Some uncertainty remains on the placement of the root. Trees on which the Hystricognathi are the basal sister group to Myodonta, Geomyoidea, Pedetidae, and Castoridae are also found within a Bayesian 95% credible set, as estimated by Metropolis-coupled Markov chain Monte Carlo sampling.     

Codon reading by tRNAAla with modified uridine in the wobble position

tRNAs reading four-codon families often have a modified uridine, cmo(5)U(34), at the wobble position of the anticodon. Here, we examine the effects on the decoding mechanism of a cmo(5)U modification in tRNA(1B)(Ala), anticodon C(36)G(35)cmo(5)U(34). tRNA(1B)(Ala) reads its cognate codons in a manner that is very similar to that of tRNA(Phe). As Ala codons are GC rich and Phe codons AU rich, this similarity suggests a uniform decoding mechanism that is independent of the GC content of the codon-anticodon duplex or the identity of the tRNA. The presence of cmo(5)U at the wobble position of tRNA(1B)(Ala) permits fairly efficient reading of non-Watson-Crick and nonwobble bases in the third codon position, e.g., the GCC codon. The ribosome accepts the C-cmo(5)U pair as an almost-correct base pair, unlike third-position mismatches, which lead to the incorporation of incorrect amino acids and are efficiently rejected.     

Sense codon in Escherichia coli are translated in context

The nucleotide frequencies 5' and 3' to the sense codons in highly and weakly expressed genes have been investigated by the chi-squares method. A comparison between the experimental and computer-generated random nucleotide sequences (in which each codon is substituted by a random synonymous one) was made. It was shown that the choice of a particular codon among the synonymous ones in a given position of the gene depends on the three nucleotides 3' and 5' adjacent to the codon in highly expressed genes (the triplet 3' and a single nucleotide 5' to the codons in weakly expressed genes). Concrete patterns for the preferable choice of synonymous codons depending on their contexts are presented. It is suggested that these constraints are related to the efficiency of messenger translation. The constraints on the amino acid sequences of encoded proteins also lead to statistically significant bases in nucleotide frequencies around the sense codons. The biological role of these constraints is discussed.     

Clonorchis sinensis: codon usage in nuclear genes

Codon usage in Clonorchis sinensis was analyzed using 12,515 codons from 38 coding sequences. Total GC content was 49.83%, and GC1, GC2 and GC3 contents were 56.32%, 43.15% and 50.00%, respectively. The effective number of codons converged at 51-53 codons. When plotted against total GC content or GC3, codon usage was distributed in relation to GC3 biases. Relative synonymous codon usage for each codon revealed a single major trend, which was highly correlated with GC content at the third position when codons began with A or U at the first two positions. In codons beginning with G or C base at the first two positions, the G or C base rarely occurred at the third position. These results suggest that codon usage is shaped by a bias towards G or C at the third base, and that this is affected by the first and second bases.     

How Mitochondria Redefine the Code

Annotated, complete DNA sequences are available for 213 mitochondrial genomes from 132 species. These provide an extensive sample of evolutionary adjustment of codon usage and meaning spanning the history of this organelle. Because most known coding changes are mitochondrial, such data bear on the general mechanism of codon reassignment. Coding changes have been attributed variously to loss of codons due to changes in directional mutation affecting the genome GC content (Osawa and Jukes 1988), to pressure to reduce the number of mitochondrial tRNAs to minimize the genome size (Anderson and Kurland 1991), and to the existence of transitional coding mechanisms in which translation is ambiguous (Schultz and Yarus 1994a). We find that a succession of such steps explains existing reassignments well. In particular, (1) Genomic variation in the prevalence of a codon's third-position nucleotide predicts relative mitochondrial codon usage well, though GC content does not. This is because A and T, and G and C, are uncorrelated in mitochondrial genomes. (2) Codons predicted to reach zero usage (disappear) do so more often than expected by chance, and codons that do disappear are disproportionately likely to be reassigned. However, codons predicted to disappear are not significantly more likely to be reassigned. Therefore, low codon frequencies can be related to codon reassignment, but appear to be neither necessary nor sufficient for reassignment. (3) Changes in the genetic code are not more likely to accompany smaller numbers of tRNA genes and are not more frequent in smaller genomes. Thus, mitochondrial codons are not reassigned during demonstrable selection for decreased genome size. Instead, the data suggest that both codon disappearance and codon reassignment depend on at least one other event. This mitochondrial event (leading to reassignment) occurs more frequently when a codon has disappeared, and produces only a small subset of possible reassignments. We suggest that coding ambiguity, the extension of a tRNA's decoding capacity beyond its original set of codons, is the second event. Ambiguity can act alone but often acts in concert with codon disappearance, which promotes codon reassignment. Received: 26 October 2000 / Accepted: 19 January 2001     

Overestimation of nonsynonymous/synonymous rate ratio by reverse-translation of aligned amino acid sequences

In the analysis of protein-coding nucleotide sequences, the ratio of the number of nonsynonymous substitutions to that of synonymous substitutions (d(N)/d(S)) is used as an indicator for the direction and magnitude of natural selection operating at the amino acid sequence level. The d(S) and d(N) values are estimated based on the comparison of homologous codons, which are often identified by converting (reverse-translating) aligned amino acid sequences into codon sequences. In this method, however, homologous codons may be mis-identified when frame-shifts occurred or amino acid sequences were mis-aligned, which may lead to overestimation of the d(N)/d(S) ratio. Here the effect of reverse-translating aligned amino acid sequences on the estimation of d(N)/d(S) ratio was examined through a large-scale analysis of protein-coding nucleotide sequences from vertebrate species. Apparently, 1-9% of codon sites that were identified as homologous with reverse-translation contained non-homologous codons, where the d(N)/d(S) ratio was unduly high. By correcting the d(N)/d(S) ratio for these codon sites, it was inferred that the ratio was 5-43% overestimated with reverse-translation. These results suggest that caution should be exerted in the study of natural selection using the d(N)/d(S) ratio by reverse-translating aligned amino acid sequences.     

Third position base changes in codons 5' and 3' adjacent UGA codons affect UGA suppression in vivo

The base sequence around nonsense codons affects the efficiency of nonsense codon suppression. Published data, comparing different nonsense sites in a mRNA, implicate the two bases downstream of the nonsense codon as major determinants of suppression efficiency. However, the results we report here indicate that the nature of the contiguous upstream codon can also affect nonsense suppression, as can the third (wobble) base of the contiguous downstream codon. These conclusions are drawn from experiments in which the two Ser codons UCU233 and UCG235 in a nonsense mutant form (UGA234) of the trpA gene in Escherichia coli have been replaced with other Ser codons by site-directed mutagenesis. Suppression of these trpA mutants has been studied in the presence of a UGA nonsense suppressor derived from glyT. We speculate that the non-site-specific effects of the two adjacent downstream bases may be largely at the level of the termination process, whereas more site-specific or codon-specific effects may operate primarily on the activity of the suppressor tRNA.     

Internal symmetry in nucleotide sequences of genes encoding the dolichol cycle enzymes

In genes alg5, alg8 and swp1 of Saccharomyces cerevisiae, gpt of Schizosaccharomyces pombe and human gene alg6, encoding the dolichol cycle enzymes, a mirror type internal symmetry was found. The symmetry was detected in both complete nucleotide sequences and sequences of the first, second and third nucleotide bases of codons. In the encoding gene regions the density of single- and double-point centres of the internal symmetry for sequences of the second bases was higher in comparison with the sequences of the first and third bases of codons, whereas in the noncoding regions degrees of symmetry of the first, second and third bases sequences did not differ significantly. A clear positive correlation was revealed in the internal symmetry distribution in the second base sequences of codons in genes, on the one hand, and in the gene encoded amino acid sequences, on the other hand. The maximum internal symmetry of gene segments encoding the functionally important regions of proteins was found at the level of the second base sequences. The obtained results corroborate a hypothesis about the determining role of the second bases of codons in encoding amino acid residues. The investigation of internal symmetry in nucleotide sequences has first shown the existence of internal symmetry at the level of gene primary structure.     

Synonymous Codon Usage Bias Dependent on Local Nucleotide Context in the Class Deinococci

To study the evolution of mutation biased synonymous codon usage, we examined nucleotide co-occurrence patterns in the Deinococcus radiodurans, D. geothermalis, and Thermus thermophilus genomes for nucleotide replacement dependent on the surrounding nucleotide context. Nucleotides on the third codon site were found to be strongly correlated with nucleotide sites at most six nucleotides away in all three species, where abundance patterns were dependent on whether two nucleotides share the same purine(R)/pyrimidine(Y) status. In the class Deinococci adjacent third site nucleotides were strongly correlated, where NNR|NNR and NNY|NNY codon pairs were overabundant while NNR|NNY and NNY|NNR codon pairs were underabundant. By far the largest deviations in all three species occur for NN(YR)|(YR)NN codon pairs. In the Thermus species, the NNY|YNN and NNR|RNN codon pairs were overabundant versus the underabundant NNY|RNN and NNR|YNN codon pairs, whereas in the Deinococcus species the opposite over-/underabundance relationship held for adjacent (GC) bases. We also observed a weaker overabundance of NNR|NRN and NNY|NYN codon pairs versus the underabundant NNR|NYN and NNY|NRN codon pairs. The perfect purine/pyrimidine symmetry of each of these cases, plus the lack of significant deviations for nucleotide pairs on other length scales up to 20 codons apart demonstrates that a pervasive pattern of nucleotide replacement dependent on local nucleotide context, and not codon bias, has occurred in these species. This nucleotide replacement has led to modified synonymous codon usage within the class Deinococci that affects which codons are positioned at particular codon sites dependent on the local nucleotide context.     

Structural Constraints on RNA Virus Evolution

The recently discovered hepatitis G virus (HGV) or GB virus C (GBV-C) is widely distributed in human populations, and homologues such as HGV/GBV-CCPZ and GBV-A are found in a variety of different primate species. Both epidemiological and phylogenetic analyses support the hypothesis that GB viruses coevolved with their primate hosts, although their degree of sequence similarity appears incompatible with the high rate of sequence change of HGV/GBV-C over short observation periods. Comparison of complete coding sequences (8,500 bases) of different genotypes of HGV/GBV-C showed an excess of invariant synonymous sites (at 23% of all codons) compared with the frequency expected by chance (10%). To investigate the hypothesis that RNA secondary-structure formation through internal base pairing limited sequence variability at these sites, an algorithm was developed to detect covariant sites among HGV/GBV-C sequences of different genotypes. At least 35 covariant sites that were spatially associated with potential stem-loop structures were detected, whose positions correlated with positions in the genome that showed reductions in synonymous variability. Although the functional roles of the predicted secondary structures remain unclear, the restriction of sequence change imposed by secondary-structure formation provides a mechanism for differences in net rate of accumulation of nucleotide substitutions at different sites. However, the resulting disparity between short- and long-term rates of sequence change of HGV/GBV-C violates the assumptions of the "molecular clock." This places a major restriction on the use of nucleotide or amino acid sequence comparisons to calculate times of divergence of other viruses evolving under the same structural constraints as GB viruses.     

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

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

Unusually biased nucleotide sequences on sense strands of Flavobacterium sp. genes produce nonstop frames on the corresponding antisense strands.

From investigation of eight Flavobacterium sp. genes encoding enzyme proteins, it was found that six genes had nonstop frames (NSFs) on the antisense strands, and base sequences of the genes are mainly composed of repeating triplet sequence(s), 5'-GNC-3' (where G and C are guanine and cytosine, and N is either of the four bases), in the reading frames. Thus, we concluded that the biased nucleotide sequences on the sense strands produce NSFs on the corresponding antisense strands. Furthermore, from the precise alignments of both nucleotide and amino acid sequences of two related Flavobacterium sp. genes, nyIB and nyIB', it was found that base replacements might have occurred symmetrically in the codons. That is, transversions between G and C were observed at high frequencies at the first and third positions of codons, but not at the second positions. At the first position, AG base transitions were observed much more than similar CT transitions, whereas CT transitions were found at the third positions at a relatively high frequency. These suggest that symmetrical base replacements in codons might be the main contribution to evolution in Flavobacterium sp. genes.     

Heterogeneous Selective Pressure Acting on Influenza B Victoria- and Yamagata-Like Hemagglutinins

As a consequence of immune pressure, influenza virus hemagglutinin presents some of its amino acids under positive selection. Several authors have reported the existence of influenza A hemagglutinin codons under positive selective pressure (PSP). In this framework, the present work objectives were to demonstrate the presence of PSP and evaluate its effects on Victoria- and Yamagata-like influenza B viruses. Methodology adopted consisted in estimating the acceptance rate of nonsynonymous substitutions (ω = d_N/d_S) that describe the strength of selective pressure and identifying codons that may be positively selected, applying a set of continuous-time Markov chain codon-substitution models. Two groups of HA1 sequences (140 from Yamagata and 60 from Victoria lineage) were used. All the model maximum-likelihood estimates were obtained using codeml software application (PAML 3.15). The hypothesis of no existence of sites under PSP was rejected for both lineages (p < 0.001), using likelihood ratio tests. These results demonstrate the presence of positive selection acting on hemagglutinin of both Yamagata- and Victoria-like influenza B viruses. Several different sites were identified to be under PSP on Yamagata and Victoria hemagglutinins. Sites found with a posterior probability > 0.95 were codons 197 and 199 in both lineages, codon 75 in the Yamagata lineage, and codon 129 in the Victoria lineage. The detected amino acids are located at or near antigenic sites in influenza A virus H3 hemagglutinin. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

不具有3-碱基周期性的编码序列初探

对120个较短编码序列(＜1 200 bp)的Fourier频谱进行分析表明,3-碱基周期性在短编码序列中并不是绝对存在的.统计分析提示,编码序列有无3-碱基周期性与序列的碱基组成和分布、所编码蛋白质氨基酸的选用和顺序以及同义密码子的使用都有一定的关系.一般地,非周期-3序列中A+U含量高于G+C含量,周期-3序列的情况则相反;非周期-3序列中碱基在密码子三个位点上的分布比周期-3序列中的分布均匀;非周期-3序列密码子和氨基酸的使用偏向没有周期-3序列的大.在利用Fourier分析方法预测DNA序列中的基因和外显子时,应充分考虑到这些现象.     

Phylogenetic inference: how much evolutionary history is knowable?

In order to reconstruct phylogenetic trees from extremely dissimilar sequences it is necessary to estimate accurately the extent of sequence divergence. In this paper a new method of sequence analysis, Markov triple analysis, is developed for determining the relative frequencies of nucleotide substitutions within the three branches of a three-taxon dendrogram. Assuming that nucleotide sites are independently and identically distributed and assuming a Markov model for nucleotide (or protein) evolution, it is shown that the unique Markov matrices can be reconstructed given only the joint probability distribution relating three taxa. (In the much simpler case involving only two taxa and two character states, Markov matrices can also be reconstructed, provided symmetry assumptions are placed on the elements of the matrices.) The method is illustrated using sequence data from the combined first and second codon positions derived from complete human, mouse, and cow mitochondrial sequences.      

Codon reading patterns in Drosophila melanogaster mitochondria based on their tRNA sequences: a unique wobble rule in animal mitochondria.

Mitochondrial (mt) tRNA(Trp), tRNA(Ile), tRNA(Met), tRNA(Ser)GCU, tRNA(Asn)and tRNA(Lys)were purified from Drosophila melanogaster (fruit fly) and their nucleotide sequences were determined. tRNA(Lys)corresponding to both AAA and AAG lysine codons was found to contain the anticodon CUU, C34 at the wobble position being unmodified. tRNA(Met)corresponding to both AUA and AUG methionine codons was found to contain 5-formylcytidine (f(5)C) at the wobble position, although the extent of modification is partial. These results suggest that both C and f(5)C as the wobble bases at the anticodon first position (position 34) can recognize A at the codon third position (position 3) in the fruit fly mt translation system. tRNA(Ser)GCU corresponding to AGU, AGC and AGA serine codons was found to contain unmodified G at the anticodon wobble position, suggesting the utilization of an unconventional G34-A3 base pair during translation. When these tRNA anticodon sequences are compared with those of other animal counterparts, it is concluded that either unmodified C or G at the wobble position can recognize A at the codon third position and that modification from A to t(6)A at position 37, 3'-adjacent to the anticodon, seems to be important for tRNA possessing C34 to recognize A3 in the mRNA in the fruit fly mt translation system.     

Comparative analysis of UV-induced mutability of ten different codon units in position 211 of the Escherichia coli trpA gene

Summary To study the influence of the local base, composition upon UV-induced mutability the reversion frequencies of ten trpA mutants with known codon sequences at position 211 were compared. Comparison of mutant strains reverting by the same base substitution type but with different dipyrimidinic, sequences reveals the mutagenic character of the pyrimidine-pyrimidine (6-4) photoproduct. The codons GAC and GAT, both reverting by AT-GC transitions and AT-CG transversions in the middle position and harboring the dipyrimidinic sequence CT in the opposite strand, differ, in their reversion frequencies sixfold. This difference can only be due to the influence of the different bases in the third codon position upon the mutability of adjacent second bases.