Codon usage determines translation rate in Escherichia coli

Determinants of translation efficiency and accuracy

A given protein sequence can be encoded by an astronomical number of alternative nucleotide sequences. Recent research has revealed that this flexibility provides evolution with multiple ways to tune the efficiency and fidelity of protein translation and folding.

Proper functioning of biological cells requires that the process of protein expression be carried out with high efficiency and fidelity. Given an amino-acid sequence of a protein, multiple degrees of freedom still remain that may allow evolution to tune efficiency and fidelity for each gene under various conditions and cell types. Particularly, the redundancy of the genetic code allows the choice between alternative codons for the same amino acid, which, although ‘synonymous,' may exert dramatic effects on the process of translation. Here we review modern developments in genomics and systems biology that have revolutionized our understanding of the multiple means by which translation is regulated. We suggest new means to model the process of translation in a richer framework that will incorporate information about gene sequences, the tRNA pool of the organism and the thermodynamic stability of the mRNA transcripts. A practical demonstration of a better understanding of the process would be a more accurate prediction of the proteome, given the transcriptome at a diversity of biological conditions.

The codon usage model of the context flanking each cleavage site in the polyprotein of foot-and-mouth disease virus

To investigate the codon usage pattern of the contexts flanking 11 cleavage sites of foot-and-mouth disease virus (FMDV) polyprotein, the codon usage model of the corresponding codon position and the synonymous codon usage in the target contexts of 66 strains were characterized by two simple methods based on the relative synonymous codon usage value. The synonymous codons usage pattern was also compared between this virus and two species of hosts (cattle and domestic pig). It is indicated that FMDV bore a general resemblance to the hosts in terms of the synonymous codon usage pattern. This feature may help FMDV to utilize translational resources of host efficiently. The two amino acid residues constituting each cleavage site contain at least one conserved residue. It was noticed that the codon usage model with the strong bias appeared in some specific positions in the target contexts, and the under-represented synonymous codons, AUA for Ile, CUA for Leu, UUA for Leu and GUA for Val, are preferentially used in these positions. These under-represented synonymous codons likely play role in regulating the translation rate and influencing the secondary structure of the contexts flanking the cleavage sites.

Quantitative transcript analysis of the inducible expression system pSIP: comparison of the overexpression of Lactobacillus spp. β-galactosidases in Lactobacillus plantarum

BACKGROUND

Two sets of overlapping genes, lacLMReu and lacLMAci, encoding heterodimeric β-galactosidases from Lactobacillus reuteri and Lactobacillus acidophilus, respectively, have previously been cloned and expressed using the pSIP vector system and Lactobacillus plantarum WCSF1 as host. Despite the high similarity between these lacLM genes and the use of identical cloning and expression strategies, strains harboring lacLMReu produced about twenty-fold more β-galactosidase than strains containing lacLMAci.

RESULTS

In this study, the plasmid copy numbers (PCN) of expression vectors pEH9R (lacLMReu) and pEH9A (lacLMAci) as well as the transcription levels of both lacLM genes were compared using quantitative PCR methods. Analyses of parallel fermentations of L. plantarum harboring either pEH9R or pEH9A showed that the expression plasmids were present in similar copy numbers. However, transcript levels of lacLM from L. reuteri (pEH9R) were up to 18 times higher than those of lacLM from L. acidophilus (pEH9A). As a control, it was shown that the expression levels of regulatory genes involved in pheromone-induced promoter activation were similar in both strains.

CONCLUSION

The use of identical expression strategies for highly similar genes led to very different mRNA levels. The data indicate that this difference is primarily caused by translational effects that are likely to affect both mRNA synthesis rates and mRNA stability. These translational effects thus seem to be a dominant determinant for the success of gene expression efforts in lactobacilli.

Coupling between codon usage, translation and protein export in Escherichia coli

Proteins destined for export via the Sec-dependent pathway are synthesized with a short N-terminal signal peptide. A requirement for export is that the proteins are in a translocationally competent state. This is a loosely folded state that allows the protein to pass through the SecYEG apparatus and pass into the periplasm. In order to maintain pre-secretory proteins in an export-competent state, there are many factors that slow the folding of the pre-secretory protein in the cytoplasm. These include cytoplasmic chaperones, such as SecB, and the signal recognition particle, which bind the pre-secretory protein and direct it to the cytoplasmic membrane for export. Recently, evidence has been published that non-optimal codons in the signal sequence are important for a time-critical early event to allow the correct folding of pre-secretory proteins. This review details the recent developments in folding of the signal peptide and the pre-secretory protein.

Stochastic sequence-level model of coupled transcription and translation in prokaryotes

Background

In prokaryotes, transcription and translation are dynamically coupled, as the latter starts before the former is complete. Also, from one transcript, several translation events occur in parallel. To study how events in transcription elongation affect translation elongation and fluctuations in protein levels, we propose a delayed stochastic model of prokaryotic transcription and translation at the nucleotide and codon level that includes the promoter open complex formation and alternative pathways to elongation, namely pausing, arrests, editing, pyrophosphorolysis, RNA polymerase traffic, and premature termination. Stepwise translation can start after the ribosome binding site is formed and accounts for variable codon translation rates, ribosome traffic, back-translocation, drop-off, and trans-translation.

Results

First, we show that the model accurately matches measurements of sequence-dependent translation elongation dynamics. Next, we characterize the degree of coupling between fluctuations in RNA and protein levels, and its dependence on the rates of transcription and translation initiation. Finally, modeling sequence-specific transcriptional pauses, we find that these affect protein noise levels.

Conclusions

For parameter values within realistic intervals, transcription and translation are found to be tightly coupled in Escherichia coli, as the noise in protein levels is mostly determined by the underlying noise in RNA levels. Sequence-dependent events in transcription elongation, e.g. pauses, are found to cause tangible effects in the degree of fluctuations in protein levels.

Purification of poly-histidine-tagged proteins

His-tagging is the most widespread and versatile strategy used to purify recombinant proteins for biochemical and structural studies. Recombinant DNA methods are first used to engineer the addition of a short tract of poly-histidine tag (His-tag) to the N terminus or C terminus of a target protein. The His-tag is then exploited to enable purification of the "tagged" protein by immobilised metal affinity chromatography (IMAC). Here, we describe efficient procedures for the isolation of highly purified His-tagged target proteins from an Escherichia coli host using IMAC.

Synonymous but not the same: the causes and consequences of codon bias

Despite their name, synonymous mutations have significant consequences for cellular processes in all taxa. As a result, an understanding of codon bias is central to fields as diverse as molecular evolution and biotechnology. Although recent advances in sequencing and synthetic biology have helped to resolve longstanding questions about codon bias, they have also uncovered striking patterns that suggest new hypotheses about protein synthesis. Ongoing work to quantify the dynamics of initiation and elongation is as important for understanding natural synonymous variation as it is for designing transgenes in applied contexts.

Design and construction of functional AAV vectors

Using the basic principles of molecular biology and laboratory techniques presented in this chapter, researchers should be able to create a wide variety of AAV vectors for both clinical and basic research applications. Basic vector design concepts are covered for both protein coding gene expression and small non-coding RNA gene expression cassettes. AAV plasmid vector backbones (available via AddGene) are described, along with critical sequence details for a variety of modular expression components that can be inserted as needed for specific applications. Protocols are provided for assembling the various DNA components into AAV vector plasmids in Escherichia coli, as well as for transferring these vector sequences into baculovirus genomes for large-scale production of AAV in the insect cell production system.

The N-terminus of GalE induces tmRNA activity in Escherichia coli

BACKGROUND

The tmRNA quality control system recognizes stalled translation complexes and facilitates ribosome recycling in a process termed 'ribosome rescue'. During ribosome rescue, nascent chains are tagged with the tmRNA-encoded SsrA peptide, which targets tagged proteins for degradation. In Escherichia coli, tmRNA rescues ribosomes arrested on truncated messages, as well as ribosomes that are paused during elongation and termination.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we describe a new translational pausing determinant that leads to SsrA peptide tagging of the E. coli GalE protein (UDP-galactose 4-epimerase). GalE chains are tagged at more than 150 sites, primarily within distinct clusters throughout the C-terminal domain. These tagging sites do not correspond to rare codon clusters and synonymous recoding of the galE gene had little effect on tagging. Moreover, tagging was largely unaffected by perturbations that either stabilize or destabilize the galE transcript. Examination of GalE-thioredoxin (TrxA) fusion proteins showed that the GalE C-terminal domain is no longer tagged when fused to an N-terminal TrxA domain. Conversely, the N-terminus of GalE induced tagging within the fused C-terminal TrxA domain.

CONCLUSIONS/SIGNIFICANCE

These findings suggest that translation of the GalE N-terminus induces subsequent tagging of the C-terminal domain. We propose that co-translational maturation of the GalE N-terminal domain influences ribosome pausing and subsequent tmRNA activity.

Analysis of codon usage in bovine viral diarrhea virus

Bovine viral diarrhea virus (BVDV) is a widespread virus in beef and dairy herds. BVDV has been grouped into two genotypes, genotype 1 and genotype 2. In this study, the relative synonymous codon usage (RSCU) values, effective number of codon (ENC) values and nucleotide content were investigated, and a comparative analysis of codon usage patterns for open reading frames (ORFs) of 22 BVDV genomes, including 14 of genotype 1 and 8 of genotype 2, was carried out. A high A+U content and low codon bias were found in BVDV genomes. Depending on the RSCU data, it was found that there was a significant variation in bias of codon usage between the two genotypes, and a geographic factor exists only in genotype-1 of BVDV. The RSCU data have a negative correlation with general average hydrophobicity (GRAVY), aromaticity and nucleotide content. Furthermore, the overall abundance of C and U has no effect on the synonymous codon usage patterns. In contrast, the A and G content showed a significant correlation with the nucleotide content at the third position. In addition, the codon usage patterns of BVDV are similar to those of 22 conserved genes of Bos taurus. Taken together, the genetic characteristics of BVDV possibly result from interactions between natural seclection and mutation pressure.

Slow sites in an exclusion process with limited resources

We introduce slow bottleneck sites into a recent extension of the totally asymmetric exclusion process where hopping rates are allowed to vary dynamically with the availability of resources. In the context of messenger RNA (mRNA) translation in biology, this refers to the availability of amino acid-transfer-RNA (aa-tRNA) complexes which act as the source of amino acids for protein production. We study a simple designer mRNA with a single defect codon in the center. As well as the familiar queuing behavior we also observe a regime within the queuing phase where the queue becomes less severe as the aa-tRNAs become depleted.

Patterns of codon usage bias in Silene latifolia

Patterns of codon usage bias (CUB) convey useful information about the selection on synonymous codons induced by gene expression and contribute to an understanding of substitution patterns observed at synonymous sites. They can also be informative about the distinctive evolutionary properties of sex chromosomes such as genetic degeneration of the Y chromosome, dosage compensation, and hemizygosity of the X chromosome in males, which can affect the selection on codon usage. Here, we study CUB in Silene latifolia, a species of interest for studying the early stages of sex chromosome evolution. We have obtained a large expressed sequence tag data set containing more than 1,608 sequence fragments by 454 sequencing. Using three different methods, we conservatively define 21 preferred codons. Interestingly, the preferred codons in S. latifolia are almost identical to those in Arabidopsis thaliana, despite their long divergence time (we estimate average nonsynonymous site divergence to be 0.216, and synonymous sites are saturated). The agreement suggests that the nature of selection on codon usage has not changed significantly during the long evolutionary time separating the two species. As in many other organisms, the frequency of preferred codons is negatively correlated with protein length. For the 43 genes with both exon and intron sequences, we find a positive correlation between gene expression levels and GC content at third codon positions, but a strong negative correlation between expression and intron GC content, suggesting that the CUB we detect in S. latifolia is more likely to be due to natural selection than to mutational bias. Using polymorphism data, we detect evidence of ongoing natural selection on CUB, but we find little support for effects of biased gene conversion. An analysis of ten sex-linked genes reveals that the X chromosome has experienced significantly more unpreferred to preferred than preferred to unpreferred substitutions, suggesting that it may be evolving higher CUB. In contrast, numbers of substitutions between preferred and unpreferred codons are similar in both directions in the Y-linked genes, contrary to the expectation of genetic degeneration.

Effect of correlated tRNA abundances on translation errors and evolution of codon usage bias

Despite the fact that tRNA abundances are thought to play a major role in determining translation error rates, their distribution across the genetic code and the resulting implications have received little attention. In general, studies of codon usage bias (CUB) assume that codons with higher tRNA abundance have lower missense error rates. Using a model of protein translation based on tRNA competition and intra-ribosomal kinetics, we show that this assumption can be violated when tRNA abundances are positively correlated across the genetic code. Examining the distribution of tRNA abundances across 73 bacterial genomes from 20 different genera, we find a consistent positive correlation between tRNA abundances across the genetic code. This work challenges one of the fundamental assumptions made in over 30 years of research on CUB that codons with higher tRNA abundances have lower missense error rates and that missense errors are the primary selective force responsible for CUB.

Codon usage bias (CUB) is a ubiquitous and important phenomenon. CUB is thought to be driven primarily due to selection against missense errors. For over 30 years, the standard model of translation errors has implicitly assumed that the relationship between translation errors and tRNA abundances are inversely related. This is based on an implicit and unstated assumption that the distribution of tRNA abundances across the genetic code are uncorrelated. Examining these abundance distributions across 73 bacterial genomes from 20 different genera, we find a consistent positive correlation between tRNA abundances across the genetic code. We further show that codons with higher tRNA abundances are not always “optimal” with respect to reducing the missense error rate and hence cannot explain the observed patterns of CUB.

Limited resources in a driven diffusion process

The advance of particles in many driven diffusion systems depends on the availability of resources in the surrounding environment. In the balance between supply and demand of such resources we are confronted with a regime in which, under limited resource availability, the flow is markedly reduced. In the context of mRNA translation this represents the finite availability of amino acid-tRNA molecules. In this limited resources regime a severe depletion of amino acid tRNAs is also observed. These dramatic effects are vital to our understanding of translation, and are likely to also be important for the many other applications of driven diffusion models.

Silent mutations in sight: co-variations in tRNA abundance as a key to unravel consequences of silent mutations

Mutations that alter the amino acid sequence are known to potentially exert deleterious effects on protein function, whereas substitutions of nucleotides without amino acid change are assumed to be neutral for the protein's functionality. However, cumulative evidence suggests that synonymous substitutions might also induce phenotypic variability by affecting splicing accuracy, translation fidelity, and conformation and function of proteins. tRNA isoacceptors mediate the translation of codons to amino acids, and asymmetric tRNA abundance causes variations in the rate of translation of each single triplet. Consequently, the effect of a silent point mutation in the coding region could be significant due to differential abundances of the cognate tRNA(s), emphasizing the importance of precise assessment of tRNA composition. Here, we provide an overview of the methods used to quantitatively determine the concentrations of tRNA species and discuss synonymous mutations in the context of tRNA composition of the cell, thus providing a new twist on the detrimental impact of the silent mutations.

Relative codon adaptation index, a sensitive measure of codon usage bias

We propose a simple, sensitive measure of synonymous codon usage bias, the Relative Codon Adaptation Index (rCAI), as a way to discriminate better between highly biased and unbiased regions, compared with the widely used Codon Adaptation Index (CAI). CAI is a geometric mean of the relative usage of codons in a gene, and is calculated using the codon usage table trained with a set of highly expressed genes. In contrast, rCAI is computed by subtracting the background codon usage trained with two noncoding frames of highly expressed genes from the codon usage in the coding frame. rCAI has higher signal-to-noise ratio than CAI, considering that noncoding frames would not show codon bias. Translation efficiency and protein abundance correlates comparably or better with rCAI than CAI or other measures such as ‘effective number of codons’ and ‘SCUMBLE offsets’. Within overlapping coding regions, one of the two coding frames dominates in codon usage bias according to rCAI. Presumably, rCAI could substitute CAI in diverse applications.

Translational selection on SHH genes

Codon usage bias has been observed in various organisms. In this study, the correlation between SHH genes expression in some tissues and codon usage features was analyzed by bioinformatics. We found that translational selection may act on compositional features of this set of genes.

The origins of time-delay in template biopolymerization processes

Time-delays are common in many physical and biological systems and they give rise to complex dynamic phenomena. The elementary processes involved in template biopolymerization, such as mRNA and protein synthesis, introduce significant time delays. However, there is not currently a systematic mapping between the individual mechanistic parameters and the time delays in these networks. We present here the development of mathematical, time-delay models for protein translation, based on PDE models, which in turn are derived through systematic approximations of first-principles mechanistic models. Theoretical analysis suggests that the key features that determine the time-delays and the agreement between the time-delay and the mechanistic models are ribosome density and distribution, i.e., the number of ribosomes on the mRNA chain relative to their maximum and their distribution along the mRNA chain. Based on analytical considerations and on computational studies, we show that the steady-state and dynamic responses of the time-delay models are in excellent agreement with the detailed mechanistic models, under physiological conditions that correspond to uniform ribosome distribution and for ribosome density up to 70%. The methodology presented here can be used for the development of reduced time-delay models of mRNA synthesis and large genetic networks. The good agreement between the time-delay and the mechanistic models will allow us to use the reduced model and advanced computational methods from nonlinear dynamics in order to perform studies that are not practical using the large-scale mechanistic models.

Global and local depletion of ternary complex limits translational elongation

The translation of genetic information according to the sequence of the mRNA template occurs with high accuracy and fidelity. Critical events in each single step of translation are selection of transfer RNA (tRNA), codon reading and tRNA-regeneration for a new cycle. We developed a model that accurately describes the dynamics of single elongation steps, thus providing a systematic insight into the sensitivity of the mRNA translation rate to dynamic environmental conditions. Alterations in the concentration of the aminoacylated tRNA can transiently stall the ribosomes during translation which results, as suggested by the model, in two outcomes: either stress-induced change in the tRNA availability triggers the premature termination of the translation and ribosomal dissociation, or extensive demand for one tRNA species results in a competition between frameshift to an aberrant open-reading frame and ribosomal drop-off. Using the bacterial Escherichia coli system, we experimentally draw parallels between these two possible mechanisms.

Signal sequence non-optimal codons are required for the correct folding of mature maltose binding protein

Non-optimal codons are generally characterised by a low concentration of isoaccepting tRNA and a slower translation rate compared to optimal codons. In a previous study, we reported a 20-fold reduction in maltose binding protein (MBP) level when the non-optimal codons in the signal sequence were optimised. In this study, we report that the 20-fold reduction is rescued when MBP is expressed at 28 degrees C instead of 37 degrees C, suggesting that the signal sequence optimised MBP protein (MBP-opt) may be misfolded, and is being degraded at 37 degrees C. Consistent with this idea, transient induction of the heat shock proteases prior to MBP expression at 28 degrees C restores the 20-fold difference, demonstrating that the difference in production levels is due to post-translational degradation of MBP-opt by the heat-shock proteases. Analysis of the structure of purified MBP-wt and MBP-opt grown at 28 degrees C showed that although they have similar secondary structure content, MBP-opt is more resistant to thermal unfolding than is MBP-wt. The two proteins also exhibit different tryptic fragment profiles, further confirming that they are folded into conformationally different states. This is the first study to demonstrate that signal sequence non-optimal codons can influence the folding of the mature exported protein.

Increased incidence of rare codon clusters at 5' and 3' gene termini: implications for function

BACKGROUND

The process of translation can be affected by the use of rare versus common codons within the mRNA transcript.

RESULTS

Here, we show that rare codons are enriched at the 5' and 3' termini of genes from E. coli and other prokaryotes. Genes predicted to be secreted show significant enrichment in 5' rare codon clusters, but not 3' rare codon clusters. Surprisingly, no correlation between 5' mRNA structure and rare codon usage was observed.

CONCLUSIONS

Potential functional roles for the enrichment of rare codons at terminal positions are explored.

Divergent stalling sequences sense and control cellular physiology

Recent studies have identified several amino acid sequences that interact with the ribosomal interior components and arrest their own elongation. Whereas stalling of the inducible class depends on specific low-molecular weight compounds, that of the intrinsic class is released when the nascent chain is transported across or inserted into the membrane. The stalled ribosome alters messenger RNA secondary structure and thereby contributes to regulation of the cis-located target gene expression at different levels. The stalling sequences are divergent but likely to utilize non-uniform nature of the peptide bond formation reactions and are recruited relatively recently to different biological systems, possibly including those to be identified in forthcoming studies.

Computational identification of rare codons of Escherichia coli based on codon pairs preference

Background

Codon bias is believed to play an important role in the control of gene expression. In Escherichia coli, some rare codons, which can limit the expression level of exogenous protein, have been defined by gene engineering operations. Previous studies have confirmed the existence of codon pair's preference in many genomes, but the underlying cause of this bias has not been well established. Here we focus on the patterns of rarely-used synonymous codons. A novel method was introduced to identify the rare codons merely by codon pair bias in Escherichia coli.

Results

In Escherichia coli, we defined the "rare codon pairs" by calculating the frequency of occurrence of all codon pairs in coding sequences. Rare codons which are disliked in genes could make great contributions to forming rare codon pairs. Meanwhile our investigation showed that many of these rare codon pairs contain termination codons and the recognized sites of restriction enzymes. Furthermore, a new index (F_rare) was developed. Through comparison with the classical indices we found a significant negative correlation between F_rare and the indices which depend on reference datasets.

Conclusions

Our approach suggests that we can identify rare codons by studying the context in which a codon lies. Also, the frequency of rare codons (F_rare) could be a useful index of codon bias regardless of the lack of expression abundance information.

Prokaryotic expression, purification, and production of polyclonal antibody against novel human serum inhibited related protein I (SI1)

A novel serum inhibited related gene (SI1) has been cloned in our lab by using mRNA differential display analysis of U251 cells in the presence or absence of serum, the expression of SI1 was dramatically inhibited by the addition of serum to serum starved cells. Previous reports suggested the potential significance of SI1 in regulating the cell cycle. In this study, the plasmid construction, protein expression and purification, as well as the generation of anti-SI1 polyclonal antibody are described. A full-length cDNA of Si1 was inserted in a prokaryotic expression plasmid pET28-b(+) and efficiently expressed in E. coli Rosetta (DE3) strain after induction by isopropyl-b-D: -thiogalactoside. The expressed 6His-tagged SI1 fusion protein was purified by Ni(+) affinity column and then used to immunize Balb/C mice, and the anti-SI1 polyclonal antibody was purified by protein A column. To determine the sensitivity and specificity of the antibody against SI1, a cell lysate of pEGFP-N2-SI1 plasmid transiently transfected Hela cell was identified by anti-GFP monoclonal antibody and anti-SI1 polyclonal antibody. Both the GFP-SI1 fusion protein and endogenous SI1 protein in Hela cell can be recognized by the anti-SI1 polyclonal antibody. The anti-SI1 polyclonal antibody will provide a useful tool for further characterization of SI1.

Mechanism of regulating the expression of lambdaN gene by ribosomal protein at translational level

In ribosomal protein S12 mutant or L24 mutant the expression of lambdaN gene was depressed at translational level. To study its mechanism the lambdaN gene region of lambdaN -lacZ gene fusion was trimmed from its 5' end to 3' end with DNA exonuclease III (DNA exoIII) in order to alter the TIR (translational initiation region) and the coding region of lambdaN gene. After DNA sequencing 23 species of different lambdaN-lacZ fused genes were obtained. The beta-galactosidase activities of these deletants in ribosomal protein mutant were compared with that in wild type strain. The result indicated that (i) S12 mutant could affect 305 subunit's binding to the TIR of lambdaN gene messenger and cause the difficulty in forming 30s initiation complex and then decrease the efficiency of translational initiation; (ii) in S12 mutant the coding region of lambdaN gene also affected the expression lambdaN gene; (iii) in L24 mutant the inhibition of lambdaN gene expression was not related to translational initiation and the 5' end of the coding region of lambdaN gene, but related to the 3' end of lambdaN gene.

The combination of ERCC1 and XRCC1 gene polymorphisms better predicts clinical outcome to oxaliplatin-based chemotherapy in metastatic colorectal cancer

PURPOSE

To evaluate the effect of excision repair cross-complementing group 1 (ERCC1) and X-ray cross-complementing group 1 (XRCC1) gene polymorphisms on treatment outcome in patients receiving oxaliplatin-based regimens for metastatic colorectal cancer.

METHODS

Hundred and thirteen patients with a diagnosis of metastatic colorectal cancer were treated with oxaliplatin-based chemotherapy. ERCC1 codon 118C/T and XRCC1 codon 399A/G polymorphisms were tested by real-time polymerase chain reaction (RT-PCR) method in peripheral blood lymphocytes of these patients. Disease control rates and survivals were compared by types of genotypes.

RESULTS

Analyses of the patterns of the polymorphism located at ERCC1 codon 118 showed that 55 (48.67%) patients were homozygous for C/C genotype, 15 (13.27%) were homozygous for the T/T genotype, and 43 (38.06%) were heterozygous for C/T genotype. Analyses of the polymorphism located at XRCC1 codon 399 showed that 61 (53.98%) patients were homozygous for A/A genotype, 13 (11.50%) were homozygous for the G/G genotype, and 39 (34.52%) were heterozygous for A/G genotype. After two cycles of chemotherapy, there was complete response (CR) in 1 patient, partial response (PR) in 24 patients, and stable disease (SD) in 56 patients. Altogether in 81 (71.68%) patients the disease was controlled after chemotherapy. Thirty-two (28.32%) patients showed disease progression. After adjusting for some clinical factors, both the ERCC1 polymorphism and the XRCC1 polymorphism lost their roles in predicting DCR (P = 0.662, P = 0.631) and MST (P = 0.692, P = 0.572). But the combination of ERCC1 and XRCC1 polymorphisms was significantly associated with DCR (P = 0.01) and MST (P = 0.000) independently.

CONCLUSIONS

This is the first study which showed that polymorphisms of ERCC1 and XRCC1, in combination not individually, were independent predictors for DCR and OS. This may contribute to the selection of patients who would benefit from oxaliplatin-based chemotherapy for metastatic colorectal cancer.

Selection for minimization of translational frameshifting errors as a factor in the evolution of codon usage

In a wide range of genomes, it was observed that the usage of synonymous codons is biased toward specific codons and codon patterns. Factors that are implicated in the selection for codon usage include facilitation of fast and accurate translation. There are two types of translational errors: missense errors and processivity errors. There is considerable evidence in support of the hypothesis that codon usage is optimized to minimize missense errors. In contrast, little is known about the relationship between codon usage and frameshifting errors, an important form of processivity errors, which appear to occur at frequencies comparable to the frequencies of missense errors. Based on the recently proposed pause-and-slip model of frameshifting, we developed Frameshifting Robustness Score (FRS). We used this measure to test if the pattern of codon usage indicates optimization against frameshifting errors. We found that the FRS values of protein-coding sequences from four analyzed genomes (the bacteria Bacillus subtilis and Escherichia coli, and the yeasts Saccharomyces cerevisiae and Schizosaccharomyce pombe) were typically higher than expected by chance. Other properties of FRS patterns observed in B. subtilis, S. cerevisiae and S. pombe, such as the tendency of FRS to increase from the 5′- to 3′-end of protein-coding sequences, were also consistent with the hypothesis of optimization against frameshifting errors in translation. For E. coli, the results of different tests were less consistent, suggestive of a much weaker optimization, if any. Collectively, the results fit the concept of selection against mistranslation-induced protein misfolding being one of the factors shaping the evolution of both coding and non-coding sequences.

Effect of exonic splicing regulation on synonymous codon usage in alternatively spliced exons of Dscam

Background

Synonymous codon usage is typically biased towards translationally superior codons in many organisms. In Drosophila, genomic data indicates that translationally optimal codons and splice optimal codons are mostly mutually exclusive, and adaptation to translational efficiency is reduced in the intron-exon boundary regions where potential exonic splicing enhancers (ESEs) reside. In contrast to genomic scale analyses on large datasets, a refined study on a well-controlled set of samples can be effective in demonstrating the effects of particular splice-related factors. Down syndrome cell adhesion molecule (Dscam) has the largest number of alternatively spliced exons (ASEs) known to date, and the splicing frequency of each ASE is accessible from the relative abundance of the transcript. Thus, these ASEs comprise a unique model system for studying the effect of splicing regulation on synonymous codon usage.

Results

Codon Bias Indices (CBI) in the 3' boundary regions were reduced compared to the rest of the exonic regions among 48 and 33 ASEs of exon 6 and 9 clusters, respectively. These regional differences in CBI were affected by splicing frequency and distance from adjacent exons. Synonymous divergence levels between the 3' boundary region and the remaining exonic region of exon 6 ASEs were similar. Additionally, another sensitive comparison of paralogous exonic regions in recently retrotransposed processed genes and their parental genes revealed that, in the former, the differences in CBI between what were formerly the central regions and the boundary regions gradually became smaller over time.

Conclusion

Analyses of the multiple ASEs of Dscam allowed direct tests of the effect of splice-related factors on synonymous codon usage and provided clear evidence that synonymous codon usage bias is restricted by exonic splicing signals near the intron-exon boundary. A similar synonymous divergence level between the different exonic regions suggests that the intensity of splice-related selection is generally weak and comparable to that of translational selection. Finally, the leveling off of differences in codon bias over time in retrotransposed genes meets the direct prediction of the tradeoff model that invokes conflict between translational superiority and splicing regulation, and strengthens the conclusions obtained from Dscam.

Generic algorithm to predict the speed of translational elongation: implications for protein biogenesis

Synonymous codon usage and variations in the level of isoaccepting tRNAs exert a powerful selective force on translation fidelity. We have developed an algorithm to evaluate the relative rate of translation which allows large-scale comparisons of the non-uniform translation rate on the protein biogenesis. Using the complete genomes of Escherichia coli and Bacillus subtilis we show that stretches of codons pairing to minor tRNAs form putative sites to locally attenuate translation; thereby the tendency is to cluster in near proximity whereas long contiguous stretches of slow-translating triplets are avoided. The presence of slow-translating segments positively correlates with the protein length irrespective of the protein abundance. The slow-translating clusters are predominantly located down-stream of the domain boundaries presumably to fine-tune translational accuracy with the folding fidelity of multidomain proteins. Translation attenuation patterns at highly structurally and functionally conserved domains are preserved across the species suggesting a concerted selective pressure on the codon selection and species-specific tRNA abundance in these regions.

Biased codon usage in signal peptides: a role in protein export

The signal peptide of proteins exported via the general secretory pathway encodes structural features that enable the targeting and export of the protein to the periplasm. Recent studies have shown biased codon usage at the second amino acid position and a high usage of non-optimal codons within the signal peptide. Altering these biases in codon usage can have deleterious effects on protein folding and export. We propose that these codon-usage biases act in concert to optimize the export process through modulating ribosome spacing on the transcript. This highlights a new aspect of protein export and implies that codon usage in the signal peptide encodes signals that are important for protein targeting and export to the periplasm.

Evolution of the Caenorhabditis elegans genome

A fundamental problem in genome biology is to elucidate the evolutionary forces responsible for generating nonrandom patterns of genome organization. As the first metazoan to benefit from full-genome sequencing, Caenorhabditis elegans has been at the forefront of research in this area. Studies of genomic patterns, and their evolutionary underpinnings, continue to be augmented by the recent push to obtain additional full-genome sequences of related Caenorhabditis taxa. In the near future, we expect to see major advances with the onset of whole-genome resequencing of multiple wild individuals of the same species. In this review, we synthesize many of the important insights to date in our understanding of genome organization and function that derive from the evolutionary principles made explicit by theoretical population genetics and molecular evolution and highlight fertile areas for future research on unanswered questions in C. elegans genome evolution. We call attention to the need for C. elegans researchers to generate and critically assess nonadaptive hypotheses for genomic and developmental patterns, in addition to adaptive scenarios. We also emphasize the potential importance of evolution in the gonochoristic (female and male) ancestors of the androdioecious (hermaphrodite and male) C. elegans as the source for many of its genomic and developmental patterns.

A gripping tale of ribosomal frameshifting: extragenic suppressors of frameshift mutations spotlight P-site realignment

Mutants of translation components which compensate for both -1 and +1 frameshift mutations showed the first evidence for framing malleability. Those compensatory mutants isolated in bacteria and yeast with altered tRNA or protein factors are reviewed here and are considered to primarily cause altered P-site realignment and not altered translocation. Though the first sequenced tRNA mutant which suppressed a +1 frameshift mutation had an extra base in its anticodon loop and led to a textbook "yardstick" model in which the number of anticodon bases determines codon size, this model has long been discounted, although not by all. Accordingly, the reviewed data suggest that reading frame maintenance and translocation are two distinct features of the ribosome. None of the -1 tRNA suppressors have anticodon loops with fewer than the standard seven nucleotides. Many of the tRNA mutants potentially affect tRNA bending and/or stability and can be used for functional assays, and one has the conserved C74 of the 3' CCA substituted. The effect of tRNA modification deficiencies on framing has been particularly informative. The properties of some mutants suggest the use of alternative tRNA anticodon loop stack conformations by individual tRNAs in one translation cycle. The mutant proteins range from defective release factors with delayed decoding of A-site stop codons facilitating P-site frameshifting to altered EF-Tu/EF1alpha to mutant ribosomal large- and small-subunit proteins L9 and S9. Their study is revealing how mRNA slippage is restrained except where it is programmed to occur and be utilized.

Chapter 1. Methods to study no-go mRNA decay in Saccharomyces cerevisiae

In eukaryotic cells, conserved mRNA surveillance systems target and degrade aberrant mRNAs, eliminating translation errors that occur during protein synthesis and thereby imposing quality control of gene expression. Two such cytoplasmic quality control systems, nonsense-mediated mRNA decay and nonstop mRNA decay, have evolved to target mRNAs with aberrancies in translation. A third novel quality control system has been identified for yeast mRNAs with defects in translation elongation due to strong translation pause sites. This subset of mRNAs with ribosome pause sites is recognized and targeted for degradation by an endonucleolytic cleavage in a process referred to as no-go mRNA decay (NGD). The methods described herein are designed to aid in the study of NGD in Saccharomyces cerevisiae. They include procedures to create an efficient translation elongation pause, assay decay characteristics of NGD substrates, and characterize NGD-dependent endonucleolytic cleavage of mRNA. The logic of the design and methods described can be modulated and used for the identification and analysis of novel RNA quality control pathways in other organisms.

A novel bicistronic vector for overexpressing Mycobacterium tuberculosis proteins in Escherichia coli

A putative DNA glycosylase encoded by the Rv3297 gene (MtuNei2) has been identified in Mycobacterium tuberculosis. Our efforts to express this gene in Escherichia coli either by supplementing tRNAs for rare codons or optimizing the gene with preferred codons for E. coli resulted in little or no expression. On the other hand, high-level expression was observed using a bicistronic expression vector in which the target gene was translationally coupled to an upstream leader sequence. Further comparison of the predicted mRNA secondary structures supported the hypothesis that mRNA secondary structure(s) surrounding the translation initiation region (TIR), rather than codon usage, played the dominant role in influencing translation efficiency, although manipulation of codon usage or tRNA supplementation did further enhance expression in the bicistronic vector. Addition of a cleavable N-terminal tag also facilitated gene expression in E. coli, possibly through a similar mechanism. However, since cleavage of N-terminal tags is determined by the amino acid at the P(1)' position downstream of the protease recognition sequence and results in the addition of an extra amino acid in front of the N-terminus of the protein, this strategy is not particularly amenable to Fpg/Nei family DNA glycosylases which carry the catalytic proline residue at the P(1)' position and require a free N-terminus. On the other hand, the bicistronic vector constructed here is potentially valuable particularly when expressing proteins from G/C rich organisms and when the proteins carry proline residues at the N-terminus in their native form. Thus the bicistronic expression system can be used to improve translation efficiency of mRNAs and achieve high-level expression of mycobacterial genes in E. coli.

Differential selection on gene translation efficiency between the filamentous fungus Ashbya gossypii and yeasts

Background

The filamentous fungus Ashbya gossypii grows into a multicellular mycelium that is distinct from the unicellular morphology of its closely related yeast species. It has been proposed that genes important for cell cycle regulation play central roles for such phenotypic differences. Because A. gossypii shares an almost identical set of cell cycle genes with the typical yeast Saccharomyces cerevisiae, the differences might occur at the level of orthologous gene regulation. Codon usage patterns were compared to identify orthologous genes with different gene regulation between A. gossypii and nine closely related yeast species.

Results

Here we identified 3,151 orthologous genes between A. gossypii and nine yeast species. Two groups of genes with significant differences in codon usage (gene translation efficiency) were identified between A. gossypii and yeasts. 333 genes (Group I) and 552 genes (Group II) have significantly higher translation efficiency in A. gossypii and yeasts, respectively. Functional enrichment and pathway analysis show that Group I genes are significantly enriched with cell cycle functions whereas Group II genes are biased toward metabolic functions.

Conclusion

Because translation efficiency of a gene is closely related to its functional importance, the observed functional distributions of orthologous genes with different translation efficiency might account for phenotypic differentiation between A. gossypii and yeast species. The results shed light on the mechanisms for pseudohyphal growth in pathogenic yeast species.

Programmed drug-dependent ribosome stalling

The ribosome has the intrinsic capacity to monitor the sequence and structure of the nascent peptide. This fundamental property of the ribosome is often exploited in regulation of gene expression, in particular, for activation of expression of genes conferring resistance to ribosome-targeting antibiotics. Induction of expression of these genes is controlled by the programmed stalling of the ribosome at a regulatory open reading frame located upstream of the resistance cistron. Formation of the stalled translation complex depends on the presence of an antibiotic in the ribosome exit tunnel and the sequence of the nascent peptide. In this review, we summarize our current understanding of the molecular mechanisms of drug- and nascent peptide-dependent ribosome stalling.

Computationally Optimised DNA Assembly of synthetic genes

Gene synthesis is hampered by two obstacles: improper assembly of oligonucleotides; oligonucleotide defects incurred during chemical synthesis. To overcome the first problem, we describe the employment of a Computationally Optimised DNA Assembly (CODA) algorithm that uses the degeneracy of the genetic code to design overlapping oligonucleotides with thermodynamic properties for self-assembly into a single, linear, DNA product. To address the second problem, we describe a hierarchical assembly strategy that reduces the incorporation of defective oligonucleotides into full-length gene constructs. The CODA algorithm and these biological methods enable fast, simple and reliable assemblies of sequence-correct full-length genes.

Synonymous mutations and ribosome stalling can lead to altered folding pathways and distinct minima

How can we understand a case in which a given amino acid sequence folds into structurally and functionally distinct molecules? Synonymous single-nucleotide polymorphisms in the MDR1 (multidrug resistance 1 or ABCB1) gene involving frequent-to-rare codon substitutions lead to identical protein sequences. Remarkably, these alternative sequences give a protein product with similar but different structures and functions. Here, we propose that long-enough ribosomal pause time scales may lead to alternate folding pathways and distinct minima on the folding free energy surface. While the conformational and functional differences between the native and alternate states may be minor, the MDR1 case illustrates that the barriers may nevertheless constitute sufficiently high hurdles in physiological time scales, leading to kinetically trapped states with altered structures and functions. Different folding pathways leading to conformationally similar trapped states may be due to swapping of (fairly symmetric) segments. Domain swapping is more likely in the no-pause case in which the chain elongates and folds simultaneously; on the other hand, sufficiently long pause times between such segments may be expected to lessen the chances of swapping events. Here, we review the literature in this light.

Ribosome collisions and translation efficiency: optimization by codon usage and mRNA destabilization

Individual mRNAs are translated by multiple ribosomes that initiate translation with an interval of a few seconds. The ribosome speed is codon dependent, and ribosome queuing has been suggested to explain specific data for translation of some mRNAs in vivo. By modeling the stochastic translation process as a traffic problem, we here analyze conditions and consequences of collisions and queuing. The model allowed us to determine the on-rate (0.8 to 1.1 initiations/s) and the time (1 s) the preceding ribosome occludes initiation for Escherichia coli lacZ mRNA in vivo. We find that ribosome collisions and queues are inevitable consequences of a stochastic translation mechanism that reduce the translation efficiency substantially on natural mRNAs. The cells minimize collisions by having its mRNAs being unstable and by a highly selected codon usage in the start of the mRNA. The cost of mRNA breakdown is offset by the concomitant increase in translation efficiency.

Amino acid starvation and colicin D treatment induce A-site mRNA cleavage in Escherichia coli

Escherichia coli possesses a unique RNase activity that cleaves stop codons in the ribosomal aminoacyl-tRNA binding site (A-site) during inefficient translation termination. This A-site mRNA cleavage allows recycling of arrested ribosomes by facilitating recruitment of the tmRNA*SmpB ribosome rescue system. To test whether A-site nuclease activity also cleaves sense codons, we induced ribosome pausing at each of the six arginine codons using three strategies; rare codon usage, arginine starvation, and inactivation of arginine tRNAs with colicin D. In each instance, ribosome pausing induced mRNA cleavage within the target arginine codons, and resulted in tmRNA-mediated SsrA-peptide tagging of the nascent polypeptide. A-site mRNA cleavage did not require the stringent factor ppGpp, or bacterial toxins such as RelE, which mediates a similar nuclease activity. However, the efficiency of A-site cleavage was modulated by the identity of the two codons immediately upstream (5' side) of the A-site codon. Starvation for histidine and tryptophan also induced A-site cleavage at histidine and tryptophan codons, respectively. Thus, A-site mRNA cleavage is a general response to ribosome pausing, capable of cleaving a variety of sense and stop codons. The induction of A-site cleavage during amino acid starvation suggests this nuclease activity may help to regulate protein synthesis during nutritional stress.

Genome landscapes and bacteriophage codon usage

Across all kingdoms of biological life, protein-coding genes exhibit unequal usage of synonymous codons. Although alternative theories abound, translational selection has been accepted as an important mechanism that shapes the patterns of codon usage in prokaryotes and simple eukaryotes. Here we analyze patterns of codon usage across 74 diverse bacteriophages that infect E. coli, P. aeruginosa, and L. lactis as their primary host. We use the concept of a “genome landscape,” which helps reveal non-trivial, long-range patterns in codon usage across a genome. We develop a series of randomization tests that allow us to interrogate the significance of one aspect of codon usage, such as GC content, while controlling for another aspect, such as adaptation to host-preferred codons. We find that 33 phage genomes exhibit highly non-random patterns in their GC3-content, use of host-preferred codons, or both. We show that the head and tail proteins of these phages exhibit significant bias towards host-preferred codons, relative to the non-structural phage proteins. Our results support the hypothesis of translational selection on viral genes for host-preferred codons, over a broad range of bacteriophages.

Any protein can be encoded by multiple, synonymous spellings. But organisms typically prefer one spelling over another—a phenomenon known as codon bias. Codon bias is generally understood to result from selection for synonymous spellings that increase the rate and accuracy of protein translation. In this work, we have examined the complete genomes of all sequenced viruses that infect the bacteria E. coli, P. aeruginosa, and L. lactis, and have found that many of these viral genomes also exhibit codon bias. Moreover, the degree of codon bias varies across the viral genome, as visualized using a technique called a “genome landscape.” By comparing the observed genomes to randomly drawn genomes, we demonstrate that the regions of high codon bias in these viral genomes often coincide with regions encoding structural proteins. Thus, the proteins that a virus needs to produce in high copy number utilize the same encoding as its host organism does for highly expressed proteins. Our results extend the translational theory of codon bias to the viral kingdom: parts of the viral genome are selected to obey the preferences of its host.