Variation in the pattern of synonymous and nonsynonymous difference between two fungal genomes

The proportion of synonymous nucleotide differences per synonymous site (p(S)) and the proportion of nonsynonymous differences per nonsynonymous site (p(N)) were computed at 1,993,217 individual codons in 4,133 protein-coding genes between the two yeast species Saccharomyces cerevisiae and Saccharomyces paradoxus. When the modified Nei-Gojobori method was used, significantly more codons with p(N) > p(S) were observed than expected, based on random pairing of observed p(S) and p(N) values. However, this finding was most likely explained by the presence of a strong negative correlation between the number of synonymous differences and the number of nonsynonymous differences at codons with at least one difference. As a result of this correlation, codons with p(N) > p(S) were characterized not only by unusually high p(N) but also by unusually low p(S). On the other hand, the number of codons with p(N)>p(S) (where p(S) is the mean p(S) for all codons) was very similar to the random expectation, and the observed number of 30-codon windows with p(N) > p(S) was significantly lower than the random expectation. These results imply that the occurrence of a certain number of codons or codon windows with p(N) > p(S) is expected given the nature of nucleotide substitution and need not imply the action of positive Darwinian selection.     

Codon specificity of starvation induced misreading

Summary Mistranslated derivatives of the coat protein of the bacteriophage MS2 were isolated from infected cells starved for asparagine. This protein contains a high level of lysine for asparagine substitutions. By peptide analysis and amino acid sequencing we show that there is a six-fold greater frequency of errors at AAU codons than at AAC codons. This ratio is the same as that found in unstarved cells where the overall error frequency is 100-fold less. We also demonstrate that, at least for AAC codons, context affects error frequency.     

Escherichia coli translation initiation factor 3 discriminates the initiation codon in vivo

In a genetic selection designed to isolate Escherichia coli mutations that increase expression of the IS 10 transposase gene ( tnp ), we unexpectedly obtained viable mutants defective in translation initiation factor 3 (IF3). Several lines of evidence led us to conclude that transposase expression, per se , was not increased. Rather, these mutations appear to increase expression of the tnp'–'lacZ gene fusions used in this screen, by increasing translation initiation at downstream, atypical initiation codons. To test this hypothesis we undertook a systematic analysis of start codon requirements and measured the effects of IF3 mutations on initiation from various start codons. Beginning with an efficient translation initiation site, we varied the AUG start codon to all possible codons that differed from AUG by one nucleotide. These potential start codons fall into distinct classes with regard to translation efficiency in vivo : Class I codons (AUG, GUG, and UUG) support efficient translation; Class IIA codons (CUG, AUU, AUC, AUA, and ACG) support translation at levels only 1–3% that of AUG; and Class IIB codons (AGG and AAG) permit levels of translation too low for reliable quantification. Importantly, the IF3 mutations had no effect on translation from Class I codons, but they increased translation from Class II codons 3–5-fold, and this same effect was seen in other gene contexts. Therefore, IF3 is generally able to discriminate between efficient and inefficient codons in vivo , consistent with earlier in vitro observations. We discuss these observations as they relate to IF3 autoregulation and the mechanism of IF3 function.     

On the Evolution of the Standard Genetic Code: Vestiges of Critical Scale Invariance from the RNA World in Current Prokaryote Genomes

Herein two genetic codes from which the primeval RNA code could have originated the standard genetic code (SGC) are derived. One of them, called extended RNA code type I, consists of all codons of the type RNY (purine-any base-pyrimidine) plus codons obtained by considering the RNA code but in the second (NYR type) and third (YRN type) reading frames. The extended RNA code type II, comprises all codons of the type RNY plus codons that arise from transversions of the RNA code in the first (YNY type) and third (RNR) nucleotide bases. In order to test if putative nucleotide sequences in the RNA World and in both extended RNA codes, share the same scaling and statistical properties to those encountered in current prokaryotes, we used the genomes of four Eubacteria and three Archaeas. For each prokaryote, we obtained their respective genomes obeying the RNA code or the extended RNA codes types I and II. In each case, we estimated the scaling properties of triplet sequences via a renormalization group approach, and we calculated the frequency distributions of distances for each codon. Remarkably, the scaling properties of the distance series of some codons from the RNA code and most codons from both extended RNA codes turned out to be identical or very close to the scaling properties of codons of the SGC. To test for the robustness of these results, we show, via computer simulation experiments, that random mutations of current genomes, at the rates of 10⁻¹⁰ per site per year during three billions of years, were not enough for destroying the observed patterns. Therefore, we conclude that most current prokaryotes may still contain relics of the primeval RNA World and that both extended RNA codes may well represent two plausible evolutionary paths between the RNA code and the current SGC.     

Molecular cloning and characterization of a complete Chinese hamster provirus related to intracisternal A particle genomes.

We report here the nucleotide sequence of a full-length Chinese hamster genomic proviral element, CHIAP34. CHIAP34 is 6,403 bp long with long terminal repeats of 311 bp at each end. The genetic organization of CHIAP34 was determined by comparison with intracisternal A particle (IAP) genetic elements from the mouse and Syrian hamster. Extensive homology at the nucleotide and deduced amino acid sequence levels was observed between CHIAP34 and the mouse and Syrian hamster IAP elements. CHIAP34 may represent a defective Chinese hamster IAP genetic element. The gag gene consists of 837 codons, of which 558 codons are in a single long open reading frame followed by several frameshifts. The pol gene begins with a -1 frameshift and consists of a long open reading frame of 753 codons followed by a short open reading frame of 103 codons. The putative env region contains multiple termination codons in all reading frames. CHIAP34 is representative of the predominant retroviral elements in the Chinese hamster ovary cell genome present at around 80 copies per haploid genome.     

Site-Specific Codon Bias in Bacteria

Sequences of the gapA and ompA genes from 10 genera of enterobacteria have been analyzed. There is strong bias in codon usage, but different synonymous codons are preferred at different sites in the same gene. Site-specific preference for unfavored codons is not confined to the first 100 codons and is usually manifest between two codons utilizing the same tRNA. Statistical analyses, based on conclusions reached in an accompanying paper, show that the use of an unfavored codon at a given site in different genera is not due to common descent and must therefore be caused either by sequence-specific mutation or sequence-specific selection. Reasons are given for thinking that sequence-specific mutation cannot be responsible. We are unable to explain the preference between synonymous codons ending in C or T, but synonymous choice between A and G at third sites is largely explained by avoidance of AG-G (where the hyphen indicates the boundary between codons). We also observed that the preferred codon for proline in Enterobacter cloacea has changed from CCG to CCA.     

The Selective Advantage of Synonymous Codon Usage Bias in Salmonella

The genetic code in mRNA is redundant, with 61 sense codons translated into 20 different amino acids. Individual amino acids are encoded by up to six different codons but within codon families some are used more frequently than others. This phenomenon is referred to as synonymous codon usage bias. The genomes of free-living unicellular organisms such as bacteria have an extreme codon usage bias and the degree of bias differs between genes within the same genome. The strong positive correlation between codon usage bias and gene expression levels in many microorganisms is attributed to selection for translational efficiency. However, this putative selective advantage has never been measured in bacteria and theoretical estimates vary widely. By systematically exchanging optimal codons for synonymous codons in the tuf genes we quantified the selective advantage of biased codon usage in highly expressed genes to be in the range 0.2–4.2 x 10⁻⁴ per codon per generation. These data quantify for the first time the potential for selection on synonymous codon choice to drive genome-wide sequence evolution in bacteria, and in particular to optimize the sequences of highly expressed genes. This quantification may have predictive applications in the design of synthetic genes and for heterologous gene expression in biotechnology.     

毕赤酵母表达系统外源基因表达序列优化的软件设计

作为真核生物表达系统,毕赤酵母已成功表达了包括病毒、细菌和真核生物在内的多种外源蛋白。但并非所有的外源序列在毕赤酵母中都能得以高表达,外源序列含有毕赤酵母中低频密码子或是序列中A T含量不合理是限制高表达的重要原因。根据毕赤酵母密码子使用的有关特性,利用遗传算法,建立外源基因表达序列优化的数学模型,并依此设计一套毕赤酵母外源基因表达序列的优化软件。     

Genome variability and capsid structural constraints of hepatitis a virus

The number of synonymous mutations per synonymous site (K(s)), the number of nonsynonymous mutations per nonsynonymous site (K(a)), and the codon usage statistic (N(c)) were calculated for several hepatitis A virus (HAV) isolates. While K(s) was similar to those of poliovirus (PV) and foot-and-mouth disease virus (FMDV), K(a) was 1 order of magnitude lower. The N(c) parameter provides information on codon usage bias and decreases when bias increases. The N(c) value in HAV was about 38, while in PV and FMDV, it was about 53. The emergence of 22 rare codons in front of 8 in PV and 7 in FMDV was detected. Most of the conserved rare codons of the P1 region were strategically located at the carboxy borders of beta barrels and alpha helices, their potential function being the assurance of proper folding of the capsid proteins through a decrease in the translation speed. This strategic location was not observed for amino acids encoded by the conserved rare codons of the 3D region. The percentage of bases with low pairing number values was higher in the latter region, suggesting a role of the conserved rare codons in the maintenance of RNA structure. Many of the rare codons in HAV are among the most frequent in humans, unlike in PV or in FMDV. This fact may be explained by the lack of cellular shutoff in HAV. One hypothesis is that HAV has evolved in order to avoid competition with its host for cellular tRNAs.     

De novo genetic codes and pure translation display

It is appealing to envision engineering translation for the genetically encoded synthesis of new classes of molecules. The complete reassignment of codons to unnatural amino acids at one or two non-adjacent sites per protein has already found wide utility (see other papers in this volume). This has been achieved by suppression at stop codons or rarely used sense codons in crude systems and in vivo. However, competing aminoacyl-tRNAs, aminoacyl-tRNA synthetases, and release factors limit efficiencies and generalization. We maximize flexibility by omitting the competing components and by reconstituting translation from His-tagged initiation and elongation factors. This approach opens up all 64 codons to amino acid reassignment and has allowed incorporation of several adjacent unnatural amino acids for the study of translation mechanism. One potential application is "peptidomimetic evolution" for ligand discovery. Toward this goal, we have demonstrated the display of polypeptides on their mRNAs in a purified translation system, termed "pure translation display."     

Construction of "small-intelligent" focused mutagenesis libraries using well-designed combinatorial degenerate primers

Site-saturation mutagenesis is a powerful tool for protein optimization due to its efficiency and simplicity. A degenerate codon NNN or NNS (K) is often used to encode the 20 standard amino acids, but this will produce redundant codons and cause uneven distribution of amino acids in the constructed library. Here we present a novel "small-intelligent" strategy to construct mutagenesis libraries that have a minimal gene library size without inherent amino acid biases, stop codons, or rare codons of Escherichia coli by coupling well-designed combinatorial degenerate primers with suitable PCR-based mutagenesis methods. The designed primer mixture contains exactly one codon per amino acid and thus allows the construction of small-intelligent mutagenesis libraries with one gene per protein. In addition, the software tool DC-Analyzer was developed to assist in primer design according to the user-defined randomization scheme for library construction. This small-intelligent strategy was successfully applied to the randomization of halohydrin dehalogenases with one or two randomized sites. With the help of DC-Analyzer, the strategy was proven to be as simple as NNS randomization and could serve as a general tool to efficiently randomize target genes at positions of interest.     

Roles of the sequence encoding tobacco etch virus capsid protein in genome amplification: requirements for the translation process and a cis-active element.

The roles of the capsid protein (CP) and the CP coding sequence of tobacco etch potyvirus (TEV) in genome amplification were analyzed. A series of frameshift-stop codon mutations that interrupted translation of the CP coding sequence at various positions were introduced into the TEV genome. A series of 3' deletion mutants that lacked the CP coding sequence beyond each of the frameshift-stop codon mutations were also produced. In addition, a series of 5' CP deletion mutants were generated. Amplification of genomes containing either frameshift-stop codon insertions after codons 1, 59, 103, and 138 or genomes containing the corresponding 3' deletions of the CP coding sequence was reduced by 100- to 1,000-fold relative to that of the parental genome in inoculated protoplasts. In contrast, a mutant containing a frameshift-stop codon after CP position 189 was amplified to 27% of the level of the parental virus, but the corresponding 3' deletion mutant lacking codons 190 to 261 was nonviable. Deletion mutants lacking CP codons 2 to 100, 2 to 150, 2 to 189, and 2 to 210 were amplified relatively efficiently in protoplasts, but a deletion mutant lacking codons 2 to 230 was nonviable. None of the amplification-defective frameshift-stop codon or deletion mutants was rescued in transgenic cells expressing TEV CP, although the transgenic CP was able to rescue intercellular movement defects of replication-competent CP mutants. Coupled with previous results, these data led to the conclusions that (i) TEV genome amplification requires translation to a position between CP codons 138 and 189 but does not require the CP product and (ii) the TEV CP coding sequence contains a cis-active RNA element between codons 211 and 246. The implications of these findings on mechanisms of RNA replication and genome evolution are discussed.     

Novel ceftazidime-resistance beta-lactamases generated by a codon-based mutagenesis method and selection

Four known and nine new ceftazidime-resistance beta-lactamases were generated by a novel, contaminating codon-based mutagenesis approach. In this method, wild-type codons are spiked with a set of mutant codons during oligonucleotide synthesis, generating random combinatorial libraries of primers that contain few codon replacements per variant. Mutant codons are assembled by tandem addition of a diluted mixture of five Fmoc-dimer amidites to the growing oligo and a mixture of four DMTr-monomer amidites to generate 20 trinucleotides that encode a set of 18 amino acids. Wild-type codons are assembled with conventional chemistry and the whole process takes place in only one synthesis column, making its automation feasible. The random and binomial behavior of this approach was tested in the polylinker region of plasmid pUC19 by the synthesis of three oligonucleotide libraries mutagenized at different rates and cloned as mutagenic cassettes. Additionally, the method was biologically assessed by mutating six contiguous codons that encode amino acids 237-243 (ABL numbering) of the TEM(pUC19) beta-lactamase, which is functionally equivalent to the clinically important TEM-1 beta-lactamase. The best ceftazidime-recognizing variant was a triple mutant, R164H:E240K: R241A, displaying a 333-fold higher resistance than the wild-type enzyme.     

Detecting molecular adaptation at individual codons in the pattern recognition protein, lipopolysaccharide- and beta-1,3-glucan-binding protein of decapods

Pattern recognition proteins play an important role in the innate immune response of invertebrates. Herein we report the evolutionary relationships among Gram-negative bacteria binding proteins (GNBPs) that were previously identified and characterized from a wide array of invertebrates. Our results, together with those obtained in previous studies, indicate that decapod lipopolysaccharide- and beta-1,3-glucan binding protein (LGBP/BGBP) has retained the crucial components for glucanase activity, and shares a common ancestor with GNBPs, as well as with the glucanase proteins of a wide range of invertebrates, rather than with GNBPs of some arthropods. However, experimental evidence of earlier studies suggested a lack of glucanase activity by these proteins, thus implying that during evolutionary time these proteins might have lost their glucan binding protein, but retained their glucan binding activity. The present results have also revealed that although a vast majority of the decapod LGBP/BGBP codons are constrained to purifying selection, certain codons are shown to have a higher rate of nonsynonymous substitutions per nonsynonymous site (dN) than synonymous substitutions per synonymous site (dS), indicating these codons have evolved adaptively (dN/dS>1). Although purifying selection (dN/dS<1) appears to be the major driving force in the evolution of a vast majority of LGBP/BGBP codons in decapods, the findings of several hotspots for nonsynonymous substitutions in this protein indicate host immune selection might play an important role in maintaining diversity among these ecologically diversified decapod species.     

Role of GC-biased mutation pressure on synonymous codon choice in Micrococcus luteus, a bacterium with a high genomic GC-content.

The GC (G + C, or G or C)-contents of codon silent positions in all two-codon sets and three codons AUY/A (IIe), and in most of the family boxes of Micrococcus luteus (genomic GC-content: 74%) are 95% to 100% in both the highly and weakly expressed genes. In some family boxes, there is a decrease in NNC codons and an increase in NNG codons from the highly expressed to weakly expressed genes without apparent involvement of NNU and NNA codons. From these observations, we conclude that the selective use of synonymous codons in M. luteus may be largely determined by GC-biased mutation pressure and that in the highly expressed genes tRNAs would act as a weak selection pressure in some family boxes. Available data suggest that the effect of selection pressure by tRNAs on the synonymous codon choice becomes more apparent in the highly expressed genes in eubacteria with intermediate GC-contents such as Escherichia coli and Bacillus subtilis, and that the U/C ratio of the codon third positions in NNU/C-type two-codon sets in the weakly expressed genes would represent the approximate magnitude of directional mutation pressure throughout eubacteria.     

Synonymous codon preferences in bacteriophage T4: a distinctive use of transfer RNAs from T4 and from its host Escherichia coli.

Codon usage data of bacteriophage T4 genes were compiled and synonymous codon preferences were investigated in comparison with tRNA availabilities in an infected cell. Since the genome of T4 is highly AT rich and its codon usage pattern is significantly different from that of its host Escherichia coli, certain codons of T4 genes need to be translated by appropriate host transfer RNAs present in minor amounts. To avoid this predicament, T4 phage seems to direct the synthesis of its own tRNA molecules and these phage tRNAs are suggested to supplement the host tRNA population with isoacceptors that are normally present in minor amounts. A positive correlation was found in that the frequency of E. coli optimal codons in T4 genes increases as the number of protein monomers per phage particle increases. A negative correlation was also found between the number of protein monomers per phage and the frequency of "T4 optimal codons", which are defined as those codons that are efficiently recognized by T4 tRNAs. From these observations it was proposed that tRNAs from the host are predominantly used for translation of highly expressed T4 genes while tRNAs from T4 tend to be used for translation of weakly expressed T4 genes. This distinctive tRNA-usage in T4 may be an optimization of translational efficiency, and an adjustment of T4-encoded tRNAs to the synonymous codon preferences, which are largely influenced by the high genomic AT-content, would have occurred during evolution.     

Growth rate dependence of transfer RNA abundance in Escherichia coli.

We have tested the predictions of a model that accounts for the codon preferences of bacteria in terms of a growth maximization strategy. According to this model the tRNA species cognate to minor and major codons should be regulated differently under different growth conditions: the isoacceptors cognate to major codons should increase at fast growth rates while those cognate to minor codons should decrease at fast growth rates. We have used a quantitative Northern blotting technique to measure the abundance of the methionine and the leucine isoacceptor families over growth rates ranging from 0.5 to 2.1 doublings per hour. Five tRNA species that are cognate to major codons (tRNA(eMet), tRNA(1fMet), tRNA(2fMet), tRNA(1Leu) and tRNA(3Leu) increase both as a relative fraction of total tRNA and in absolute concentration with increasing growth rates. Three tRNA species that are cognate to minor codons (tRNA(2Leu), tRNA(4Leu) and tRNA(5Leu) decrease as a relative fraction of total RNA and in absolute concentration with increasing growth rates. These data suggest that the abundances of groups of tRNA species are regulated in different ways, and that they are not regulated simply according to isoacceptor specificity. In particular, the data support the growth optimization model for codon bias.