Synonymous codon bias is related to gene length in Escherichia coli: selection for translational accuracy?

The levels of synonymous codon bias is shown to be positively correlated to gene length in Escherichia coli genes which are thought to be expressed at similar levels; these are genes whose products are present in multimeric proteins in equimolar amounts. It is argued that the positive correlation could be caused by selection to avoid missense errors during translation. Since the cost of producing a protein is proportional to its length, selection in favor of codons which increase accuracy should be greater in longer genes, and long genes should therefore have higher synonymous codon bias. It is also shown that there is variation in synonymous codon use which is independent of either expression level, gene length, amino acid composition, or chromosomal location. This variation is consistent with selection for translational accuracy but may have other origins.      

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

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

The influence of translational selection on codon usage in fishes from the family Cyprinidae

In this paper, the main factors shaping codon usage in three species of fishes that belong to the family Cyprinidae (namely Brachidanio rerio, Cyprinus carpio, and Carassius auratus) are reported. Correspondence analysis (COA), a commonly used multivariate statistical approach, was used to analyze codon usage bias. Our results show that the main trend is strongly correlated with the GC(3) content at silent sites of each sequence. On the other hand, the second axis discriminates between presumed highly and lowly expressed genes, a result that is confirmed by the distribution of matching expressed sequence tags (ESTs) along that axis. Translational selection appears, therefore, to influence synonymous codon usage in these fishes. The comparison of codon usages of the sequences displaying the extreme values on the second axis indicates that several codons are significantly incremented among the heavily expressed sequences. Interestingly, several of these triplets are not only shared by the three fishes but also by Xenopus laevis, another cold-blooded vertebrate in which translational selection influences codon choices. We postulate that natural selection was operative for codon usage in the last common ancestor of these fishes and Xenopus, and will probably be detected in cold-blooded vertebrates in general. Finally, we raise the possibility that the same phenomena will be found among warm-blooded vertebrates.     

The role of GC-biased gene conversion in shaping the fastest evolving regions of the human genome

GC-biased gene conversion (gBGC) is a recombination-associated evolutionary process that accelerates the fixation of guanine or cytosine alleles, regardless of their effects on fitness. gBGC can increase the overall rate of substitutions, a hallmark of positive selection. Many fast-evolving genes and noncoding sequences in the human genome have GC-biased substitution patterns, suggesting that gBGC-in contrast to adaptive processes-may have driven the human changes in these sequences. To investigate this hypothesis, we developed a substitution model for DNA sequence evolution that quantifies the nonlinear interacting effects of selection and gBGC on substitution rates and patterns. Based on this model, we used a series of lineage-specific likelihood ratio tests to evaluate sequence alignments for evidence of changes in mode of selection, action of gBGC, or both. With a false positive rate of less than 5% for individual tests, we found that the majority (76%) of previously identified human accelerated regions are best explained without gBGC, whereas a substantial minority (19%) are best explained by the action of gBGC alone. Further, more than half (55%) have substitution rates that significantly exceed local estimates of the neutral rate, suggesting that these regions may have been shaped by positive selection rather than by relaxation of constraint. By distinguishing the effects of gBGC, relaxation of constraint, and positive selection we provide an integrated analysis of the evolutionary forces that shaped the fastest evolving regions of the human genome, which facilitates the design of targeted functional studies of adaptation in humans.     

Direct selection of mutants influencing gene conversion in the yeast Schizosaccharomyces pombe

Summary In Schizosaccharomyces pombe, a suppressor-active mutation at the anticodon site of the tRNA _Ser ^UCA gene sup3 leads to opal (UGA)-specific suppression. Second-site mutations (rX) in sup3 inactivate the suppressor. The sup3-UGA, rX double mutants are genetically unstable in meiotic selfings, due to the intergenic transfer of information between sup3 and the unlinked genes sup9 and sup12 (Hofer et al. 1979; Munz and Leupold 1981; Munz et al. 1982). These three genes have considerable sequence homology over about 200 base pairs (Hottinger et al. 1982).Mutants showing a decrease or an increase of the meiotic instability at sup3 have been selected. One mutation (rec3-8) increases both the genetic instability and the frequency of intragenic recombination in sup3 by one order of magnitude. It has no effect on the stability of the nonsense alleles arg1-230 (UAA), ade6-704 and ura1-61 (UGA) or on the frequency of crossing-over between sup3 and the closely linked gene cdc8.The existence of a common genetic control over intragenic recombination and genetic instability at sup3 provides a direct way of selecting for rec mutants in homothallic haploid strains of S. pombe carrying a suppressor-inactive allele of sup3. It also supports the hypothesis that the instability of mutant alleles of this gene is due to chromosome mispairing at meiosis allowing sup3 to pair with sup9 or sup12 and then to undergo recombination by gene conversion restoring the suppressor-active allele sup3-UGA from the suppressor-inactive allele sup3-UGA,rX. Two mutations (rec2-5 and rec5-11) have no effect on intragenic recombination, but considerably reduce the meiotic instability of the sup3-UGA,rX alleles. They may suppress illegitimate pairing between sup3 and sup9 or sup12.     

Induction and repair of gene conversion in UV-sensitive mutants of yeast

Summary Photoreactivation effect on UV-induced allelic recombination has been examined using various combinations of leu 1 alleles in UV-sensitive and wild type diploid yeast, Saccharomyces cerevisiae. The frequencies of UV-induced heteroallelic reversion in UV-sensitive strains, presumably lacking dark-repair, are strikingly enhanced compared to those in wild type at the same doses under dark condition. However, these enhanced frequencies of reversion are diminished by photoreactivation almost to the level of those in wild type. The induced frequencies of homoallelic reversion (mutation) of relevant alleles are apparently lower than those of heteroallelic reversion. Phenotypic analysis for linked gene leu 1 on UV-induced heteroallelic revertants has shown that most of the revertants are of the nonreciprocal type recombination (mitotic gene conversion). These results would indicate that most of the dark-repairable damage leading to mitotic gene conversion after UV-light is due to pyrimidine dimers.On leave of absence from Radiation Center of Osaka Prefecture, Shinke-cho Sakai, Osaka, Japan.     

Statistical analysis of pentose phosphate pathway genes from eubacteria and eukarya reveals translational selection as a major force in shaping codon usage pattern

Comparative analysis of metabolic pathways among widely diverse species provides an excellent opportunity to extract information about the functional relation of organisms and pentose phosphate pathway exemplifies one such pathway. A comparative codon usage analysis of the pentose phosphate pathway genes of a diverse group of organisms representing different niches and the related factors affecting codon usage with special reference to the major forces influencing codon usage patterns was carried out. It was observed that organism specific codon usage bias percolates into vital metabolic pathway genes irrespective of their near universality. A clear distinction in the codon usage pattern of gram positive and gram negative bacteria, which is a major classification criterion for bacteria, in terms of pentose phosphate pathway was an important observation of this study. The codon utilization scheme in all the organisms indicates the presence of translation selection as a major force in shaping codon usage. Another key observation was the segregation of the H. sapiens genes as a separate cluster by correspondence analysis, which is primarily attributed to the different codon usage pattern in this genus along with its longer gene lengths. We have also analyzed the amino acid distribution comparison of transketolase protein primary structures among all the organisms and found that there is a certain degree of predictability in the composition profile except in A. fumigatus and H. sapiens, where few exceptions are prominent. In A. fumigatus, a human pathogen responsible for invasive aspergillosis, a significantly different codon usage pattern, which finally translated into its amino acid composition model portraying a unique profile in a key pentose phosphate pathway enzyme transketolase was observed.     

Revisiting the codon adaptation index from a whole-genome perspective: analyzing the relationship between gene expression and codon occurrence in yeast using a variety of models

Highly expressed genes in many bacteria and small eukaryotes often have a strong compositional bias, in terms of codon usage. Two widely used numerical indices, the codon adaptation index (CAI) and the codon usage, use this bias to predict the expression level of genes. When these indices were first introduced, they were based on fairly simple assumptions about which genes are most highly expressed: the CAI was originally based on the codon composition of a set of only 24 highly expressed genes, and the codon usage on assumptions about which functional classes of genes are highly expressed in fast-growing bacteria. Given the recent advent of genome-wide expression data, we should be able to improve on these assumptions. Here, we measure, in yeast, the degree to which consideration of the current genome-wide expression data sets improves the performance of both numerical indices. Indeed, we find that by changing the parameterization of each model its correlation with actual expression levels can be somewhat improved, although both indices are fairly insensitive to the exact way they are parameterized. This insensitivity indicates a consistent codon bias amongst highly expressed genes. We also attempt direct linear regression of codon composition against genome-wide expression levels (and protein abundance data). This has some similarity with the CAI formalism and yields an alternative model for the prediction of expression levels based on the coding sequences of genes. More information is available at http://bioinfo.mbb.yale.edu/expression/codons.     

Quantitative correlation between mRNA secondary structure around the region downstream of the initiation codon and translational efficiency in Escherichia coli

Translational efficiency in Escherichia coli is known to be strongly influenced by the secondary structure around the ribosome‐binding site and the initiation codon in the translational‐initiation region of the mRNA. Several quantitative studies have reported that translational efficiency is attributable to effects on ribosome accessibility predominantly caused by the secondary structure surrounding the ribosome‐binding site. However, the influence of mRNA secondary structure around regions downstream of the initiation codon on translational efficiency after ribosome‐binding step has not been quantitatively studied. Here, we quantitatively analyzed the relationship between secondary structure of mRNA surrounding the region downstream of the initiation codon, referred to as the downstream region (DR), and protein expression levels. Modified hairpin structures containing the initiation codon were constructed by site‐directed mutagenesis, and their effects on expression were analyzed in vivo. The minimal folding free energy (ΔG) of a local hairpin structure was found to be linearly correlated with the relative expression level over a range of fourfold change. These results demonstrate that expression level can be quantitatively controlled by changing the stability of the secondary structure surrounding the DR. Biotechnol. Bioeng. 2009; 104: 611–616 © 2009 Wiley Periodicals, Inc.     

The detection of mitotic and meiotic aneuploidy in yeast using a gene dosage selection system

Summary A system is described in which spontaneous and chemically-induced mitotic and meiotic hyperploidy can be assayed in the same diploid culture of Saccharomyces cerevisiae. Monitoring gene dosage changes at two loci on chromosome VIII, the test utilizes a leaky temperature-sensitive allele arg4-8 and low level copper resistance conferred by the single copy allele cup1 ^s. An extra chromosome VIII provides simultaneous increased dosage for both genes, resulting in colonies that are both prototrophic for arginine at 30° C and copper resistant. During mitotic cell divisions in diploids, spontaneous chromosome VIII hyperploids (trisomes and tetrasomes) occur at a frequency of 6.4×10^-6 per viable cell. Among ascospores, the spontaneous chromosome VIII disome frequency is 5.5×10^-6 per viable spore. The tubulin-binding reagent methyl benzimidazol-2-yl carbamate (MBC) elicits enhanced levels of mitotic and meiotic aneuploidy relative to control levels. The system represents a novel model for examining chromosome behavior during mitosis and meiosis and provides a sensitive and quantifiable procedure for examining chemically induced aneuploidy.     

Programmed translational frameshifting and initiation at an AUU codon in gene expression of bacterial insertion sequence IS911.

The proteins expressed by insertion sequence IS911, a member of the widespread IS3 family of elements, have been analyzed. The results indicate that three major species are produced from two consecutive reading frames. A protein of Mr 11,500, ORFA, is synthesized from an upstream reading frame. A larger protein, ORFAB, uses the same initiation codon and is produced by a -1 programmed translational frameshift between orfA and a downstream frame, orfB, whose amino acid sequence shows significant homology with retroviral integrase proteins. The orfB frame is also expressed independently in two alternative forms: the first uses a rare AUU initiation codon in the orfB phase whereas the second appears to initiate in the orfA phase and is produced by a -1 frameshift mechanism similar to that used in ORFAB expression. A specific IS911 integration reaction using a minimal active junction composed of 51 base-pairs of the right inverted repeat and a flanking phase lambda sequence resembling a second end in inverted orientation has been developed to analyze the functions of these proteins by transcomplementation in vivo. The orfA and orfB frames are shown to be essential and production of ORFAB is shown to stimulate integration in this system, suggesting that this fusion protein is the IS911 transposase.     

The RBE of neutrons for induced mitotic gene conversion in "error-prone repair" defective yeast

Mitotic gene conversion was induced in the diploid yeast strain D7.rad6 which lacks "error-prone repair" and thus does not mutate. Neutrons (14.5 MeV), 60Co gamma rays, and 150 kVp X rays delivered under oxic or anoxic conditions were compared for their ability to induce gene conversion. Doses were chosen to minimize cell killing. A lack of induced mutation in this strain at the ilv1-92 allele was confirmed. Gene conversion of the trp5-27/trp5-12 alleles was induced with a linear dose response, and the yield of convertants per gray was significantly enhanced over yields reported previously for a wild-type stain. The relative biological effectiveness (RBE) of neutrons relative to low-LET radiations was found to be about 2.2 for either oxic or anoxic radiation in contrast to wild-type where the oxic RBE was 1.7 and the anoxic RBE 2.7. Absence of the rad6 function was therefore associated with an altered RBE for the conversional end point. The oxygen enhancement ratio (OER) for gene conversion was found to be about 1.7 for all radiations in contrast to the wild type where the OER for neutrons was 1.7, but for low-LET radiations it was 2.7. As repair of ionizing damage in the rad6 strain did not lead to mutation, owing to the loss of "error-prone repair," the changes in yield, RBE, and OER were consistent with the hypothesis that some of the lesions processed by wild type to generate mutations could, in the rad6 strain, lead instead to gene conversion.     

The relative biological effectiveness of 14.5-MeV neutrons for the induction of gene conversion and mutation in yeast

The relative biological effectiveness (RBE) and oxygen enhancement ratio (OER) were determined in the yeast Saccharomyces cerevisiae for the induction of gene conversion (the product of recombinational repair) and mutation (the product of error prone repair) by 14.5-MeV neutrons in comparison with 60Co gamma rays and 150 KVp X rays. Neutron irradiation in oxic or anoxic conditions induced significantly higher yields of convertants and mutants than sparsely ionizing radiations under the same conditions. RBEs for both gene conversion and mutation under anoxia were significantly higher than under oxic conditions. RBEs for mutant induction under anoxia were lower than the RBEs for gene conversion under the same conditions. The data support the hypothesis that the production of lesions leading to the genetic consequences of gene conversion and mutation differ in their dependence upon LET and the presence of oxygen during irradiation, and therefore the two DNA repair processes which produce these end points recognize, at least in part, different classes of damage.     

The relative importance of meiotic gene conversion,selection and mutation pressure,in population genetics and evolution

Disparity in the direction of meiotic gene conversion can change allele frequencies, favouring one allele of a pair in heterozygotes. Equilibrium allele frequencies for large diploid populations are examined by means of equations relating them to meiotic gene conversion, selection and mutation for deleterious recessives, deleterious dominants, and deleterious alleles with no dominance. Using observed conversion parameters from various fungi,Zea mays andDrosophila, it is shown that conversion is generally much more important than mutation pressure and may be of greater or lesser importance than selection, depending on dominance and the strength of selection and conversion forces for the alleles involved.