Rare codon content affects the solubility of recombinant proteins in a codon bias-adjusted Escherichia coli strain

Background  

The expression of heterologous proteins in Escherichia coli is strongly affected by codon bias. This phenomenon occurs when the codon usage of the mRNA coding for the foreign protein differs from that of the bacterium. The ribosome pauses upon encountering a rare codon and may detach from the mRNA, thereby the yield of protein expression is reduced. Several bacterial strains have been engineered to overcome this effect. However, the increased rate of translation may lead to protein misfolding and insolubilization. In order to prove this assumption, the solubility of several recombinant proteins from plants was studied in a codon bias-adjusted E. coli strain.     

Differences in codon bias cannot explain differences in translational power among microbes

Background  

Translational power is the cellular rate of protein synthesis normalized to the biomass invested in translational machinery. Published data suggest a previously unrecognized pattern: translational power is higher among rapidly growing microbes, and lower among slowly growing microbes. One factor known to affect translational power is biased use of synonymous codons. The correlation within an organism between expression level and degree of codon bias among genes of Escherichia coli and other bacteria capable of rapid growth is commonly attributed to selection for high translational power. Conversely, the absence of such a correlation in some slowly growing microbes has been interpreted as the absence of selection for translational power. Because codon bias caused by translational selection varies between rapidly growing and slowly growing microbes, we investigated whether observed differences in translational power among microbes could be explained entirely by differences in the degree of codon bias. Although the data are not available to estimate the effect of codon bias in other species, we developed an empirically-based mathematical model to compare the translation rate of E. coli to the translation rate of a hypothetical strain which differs from E. coli only by lacking codon bias.     

Codon usages in different gene classes of the Escherichia coli genome

A new measure for assessing codon bias of one group of genes with respect to a second group of genes is introduced. In this formulation, codon bias correlations for Escherichia coli genes are evaluated for level of expression, for contrasts along genes, for genes in different 200 kb (or longer) contigs around the genome, for effects of gene size, for variation over different function classes, for codon bias in relation to possible lateral transfer and for dicodon bias for some gene classes. Among the function classes, codon biases of ribosomal proteins are the most deviant from the codon frequencies of the average E. coli gene. Other classes of ‘highly expressed genes’ (e.g. amino acyl tRNA synthetases, chaperonins, modification genes essential to translation activities) show less extreme codon biases. Consistently for genes with experimentally determined expression rates in the exponential growth phase, those of highest molar abundances are more deviant from the average gene codon frequencies and are more similar in codon frequencies to the average ribosomal protein gene. Independent of gene size, the codon biases in the 5′ third of genes deviate by more than a factor of two from those in the middle and 3′ thirds. In this context, there appear to be conflicting selection pressures imposed by the constraints of ribosomal binding, or more generally the early phase of protein synthesis (about the first 50 codons) may be more biased than the complete nascent polypeptide. In partitioning the E. coli genome into 10 equal lengths, pronounced differences in codon site 3 G+C frequencies accumulate. Genes near to oriC have 5% greater codon site 3 G+C frequencies than do genes from the ter region. This difference also is observed between small (100–300 codons) and large (>800 codons) genes. This result contrasts with that for eukaryotic genomes (including human, Caenorhabditis elegans and yeast) where long genes tend to have site 3 more AT rich than short genes. Many of the above results are special for E. coli genes and do not apply to genes of most bacterial genomes. A gene is defined as alien (possibly horizontally transferred) if its codon bias relative to the average gene exceeds a high threshold and the codon bias relative to ribosomal proteins is also appropriately high. These are identified, including four clusters (operons). The bulk of these genes have no known function.     

Construction and performance of heterologous polyketide‐producing K‐12‐ and B‐derived Escherichia coli

Aims: Escherichia coli has emerged as a viable heterologous host for the production of complex, polyketide natural compounds. In this study, polyketide biosynthesis was compared between different E. coli strains for the purpose of better understanding and improving heterologous production. Methods and Results: Both B and K‐12 E. coli strains were genetically modified to support heterologous polyketide biosynthesis [specifically, 6‐deoxyerythronolide B (6dEB)]. Polyketide production was analysed using a helper plasmid designed to overcome rare codon usage within E. coli. Each strain was analysed for recombinant protein production, precursor consumption, by‐product production, and 6dEB biosynthesis. Of the strains tested for biosynthesis, 6dEB production was greatest for E. coli B strains. When comparing biosynthetic improvements as a function of mRNA stability vs codon bias, increased 6dEB titres were observed when additional rare codon tRNA molecules were provided. Conclusions: Escherichia coli B strains and the use of tRNA supplementation led to improved 6dEB polyketide titres. Significance and Impact of the Study: Given the medicinal potential and growing field of polyketide heterologous biosynthesis, the current study provides insight into host‐specific genetic backgrounds and gene expression parameters aiding polyketide production through E. coli.     

An improved strategy for efficient expression and purification of soluble HIV-1 tat protein in <Emphasis Type="Italic">E.coli</Emphasis>

Although the endogenous function of Tat has been elucidated in the past twenty years, the study of its exogenous activity has been hampered due to the difficulty of large scale preparation of the active Tat protein. To express the full-length Tat protein in E.coli, the tat gene was cloned from an HIV infected patient by overlapping PCR. Rare codon usage analysis showed that rare E.coli codons, especially consecutive rare codons for Arg, account for 14% (14 of 101) rare E.coli codons in the tat gene. The expression of the HIV-1 tat gene was verified to be very poor in strain BL21 (DE3) due to the abundance of rare codons; however, tat gene expression was found to be very efficient in the host strain of Rosetta-gami B (DE3), which was supplemented with six rare tRNAs for Arg, Leu, Ile and Pro. Subsequent purification revealed that the proteins are soluble and unusually, the tagged Tat can form dimers independent of cystine disulfide bonds. The purity, integrity and molecular weight of the Tat protein were demonstrated by MALDI-TOF mass spectrometry. Reporter gene activating assay was further confirmed by investigating the transactivation activity of the recombinant Tat protein. Our improved strategy for efficient high level expression and purification of soluble Tat protein has paved the way to fully investigate its exogenous function.     

Evolutionary Lability of Context-Dependent Codon Bias in Bacteria

In bacteria, synonymous codon usage can be considerably affected by base composition at neighboring sites. Such context-dependent biases may be caused by either selection against specific nucleotide motifs or context-dependent mutation biases. Here we consider the evolutionary conservation of context-dependent codon bias across 11 completely sequenced bacterial genomes. In particular, we focus on two contextual biases previously identified in Escherichia coli; the avoidance of out-of-frame stop codons and AGG motifs. By identifying homologues of E. coli genes, we also investigate the effect of gene expression level in Haemophilus influenzae and Mycoplasma genitalium. We find that while context-dependent codon biases are widespread in bacteria, few are conserved across all species considered. Avoidance of out-of-frame stop codons does not apply to all stop codons or amino acids in E. coli, does not hold for different species, does not increase with gene expression level, and is not relaxed in Mycoplasma spp., in which the canonical stop codon, TGA, is recognized as tryptophan. Avoidance of AGG motifs shows some evolutionary conservation and increases with gene expression level in E. coli, suggestive of the action of selection, but the cause of the bias differs between species. These results demonstrate that strong context-dependent forces, both selective and mutational, operate on synonymous codon usage but that these differ considerably between genomes. Received: 6 May 1999 / Accepted: 29 October 1999     

Codon usage in twelve species of <Emphasis Type="Italic">Drosophila</Emphasis>

Background  

Codon usage bias (CUB), the uneven use of synonymous codons, is a ubiquitous observation in virtually all organisms examined. The pattern of codon usage is generally similar among closely related species, but differs significantly among distantly related organisms, e.g., bacteria, yeast, and Drosophila. Several explanations for CUB have been offered and some have been supported by observations and experiments, although a thorough understanding of the evolutionary forces (random drift, mutation bias, and selection) and their relative importance remains to be determined. The recently available complete genome DNA sequences of twelve phylogenetically defined species of Drosophila offer a hitherto unprecedented opportunity to examine these problems. We report here the patterns of codon usage in the twelve species and offer insights on possible evolutionary forces involved.     

Design Parameters to Control Synthetic Gene Expression in Escherichia coli

Background

Production of proteins as therapeutic agents, research reagents and molecular tools frequently depends on expression in heterologous hosts. Synthetic genes are increasingly used for protein production because sequence information is easier to obtain than the corresponding physical DNA. Protein-coding sequences are commonly re-designed to enhance expression, but there are no experimentally supported design principles.

Principal Findings

To identify sequence features that affect protein expression we synthesized and expressed in E. coli two sets of 40 genes encoding two commercially valuable proteins, a DNA polymerase and a single chain antibody. Genes differing only in synonymous codon usage expressed protein at levels ranging from undetectable to 30% of cellular protein. Using partial least squares regression we tested the correlation of protein production levels with parameters that have been reported to affect expression. We found that the amount of protein produced in E. coli was strongly dependent on the codons used to encode a subset of amino acids. Favorable codons were predominantly those read by tRNAs that are most highly charged during amino acid starvation, not codons that are most abundant in highly expressed E. coli proteins. Finally we confirmed the validity of our models by designing, synthesizing and testing new genes using codon biases predicted to perform well.

Conclusion

The systematic analysis of gene design parameters shown in this study has allowed us to identify codon usage within a gene as a critical determinant of achievable protein expression levels in E. coli. We propose a biochemical basis for this, as well as design algorithms to ensure high protein production from synthetic genes. Replication of this methodology should allow similar design algorithms to be empirically derived for any expression system.     

Impact of translational selection on codon usage bias in the archaeon Methanococcus maripaludis

Patterns of codon usage have been extensively studied among Bacteria and Eukaryotes, but there has been little investigation of species from the third domain of life, the Archaea. Here, we examine the nature of codon usage bias in a methanogenic archaeon, Methanococcus maripaludis. Genome-wide patterns of codon usage are dominated by a strong A + T bias, presumably largely reflecting mutation patterns. Nevertheless, there is variation among genes in the use of a subset of putatively translationally optimal codons, which is strongly correlated with gene expression level. In comparison with Bacteria such as Escherichia coli, the strength of selected codon usage bias in highly expressed genes in M. maripaludis seems surprisingly high given its moderate growth rate. However, the pattern of selected codon usage differs between M. maripaludis and E. coli: in the archaeon, strongly selected codon usage bias is largely restricted to twofold degenerate amino acids (AAs). Weaker bias among the codons for fourfold degenerate AAs is consistent with the small number of tRNA genes in the M. maripaludis genome.     

Comparative Analysis of Codon Usage Bias and Codon Context Patterns between Dipteran and Hymenopteran Sequenced Genomes

Background

Codon bias is a phenomenon of non-uniform usage of codons whereas codon context generally refers to sequential pair of codons in a gene. Although genome sequencing of multiple species of dipteran and hymenopteran insects have been completed only a few of these species have been analyzed for codon usage bias.

Methods and Principal Findings

Here, we use bioinformatics approaches to analyze codon usage bias and codon context patterns in a genome-wide manner among 15 dipteran and 7 hymenopteran insect species. Results show that GAA is the most frequent codon in the dipteran species whereas GAG is the most frequent codon in the hymenopteran species. Data reveals that codons ending with C or G are frequently used in the dipteran genomes whereas codons ending with A or T are frequently used in the hymenopteran genomes. Synonymous codon usage orders (SCUO) vary within genomes in a pattern that seems to be distinct for each species. Based on comparison of 30 one-to-one orthologous genes among 17 species, the fruit fly Drosophila willistoni shows the least codon usage bias whereas the honey bee (Apis mellifera) shows the highest bias. Analysis of codon context patterns of these insects shows that specific codons are frequently used as the 3′- and 5′-context of start and stop codons, respectively.

Conclusions

Codon bias pattern is distinct between dipteran and hymenopteran insects. While codon bias is favored by high GC content of dipteran genomes, high AT content of genes favors biased usage of synonymous codons in the hymenopteran insects. Also, codon context patterns vary among these species largely according to their phylogeny.     

Genomic transcriptional response to loss of chromosomal supercoiling in <Emphasis Type="Italic">Escherichia coli</Emphasis>

Background  

The chromosome of Escherichia coli is maintained in a negatively supercoiled state, and supercoiling levels are affected by growth phase and a variety of environmental stimuli. In turn, supercoiling influences local DNA structure and can affect gene expression. We used microarrays representing nearly the entire genome of Escherichia coli MG1655 to examine the dynamics of chromosome structure.     

Challenges Associated with Heterologous Expression of <Emphasis Type="Italic">Flavobacterium psychrophilum</Emphasis> Proteins in <Emphasis Type="Italic">Escherichia coli</Emphasis>

A two-parameter statistical model was used to predict the solubility of 96 putative virulence-associated proteins of Flavobacterium psychrophilum (CSF259-93) upon over expression in Escherichia coli. This analysis indicated that 88.5% of the F. psychrophilum proteins would be expressed as insoluble aggregates (inclusion bodies). These solubility predictions were verified experimentally by colony filtration blot for six different F. psychrophilum proteins. A comprehensive analysis of codon usage identified over a dozen codons that are used frequently in F. psychrophilum, but that are rarely used in E. coli. Expression of F. psychrophilum proteins in E. coli was often associated with production of minor molecular weight products, presumably because of the codon usage bias between these two organisms. Expression of recombinant protein in the presence of rare tRNA genes resulted in marginal improvements in the expressed products. Consequently, Vibrio parahaemolyticus was developed as an alternative expression host because its codon usage is similar to F. psychrophilum. A full-length recombinant F. psychrophilum hemolysin was successfully expressed and purified from V. parahaemolyticus in soluble form, whereas this protein was insoluble upon expression in E. coli. We show that V. parahaemolyticus can be used as an alternate heterologous expression system that can remedy challenges associated with expression and production of F. psychrophilum recombinant proteins.     

Use of Modified BL21(DE3) Escherichia coli Cells for High-Level Expression of Recombinant Peanut Allergens Affected by Poor Codon Usage

We previously cloned a panel of peanut allergens by phage display technology. Examination of the codons used in these sequences indicated that most of the cDNAs contain an excess of the least used codons in Escherichia coli, namely AGG/AGA, that correspond to a minor tRNA, the product of the dnaY gene. To achieve high-level expression of the peanut allergens, the cDNAs were subcloned into an expression vector of the pET series (Novagen) in order to produce (His)₁₀-tagged fusion proteins in conventional E. coli BL21(DE3) cells. The peanut allergens Ara h 1, Ara h 2, and Ara h 6 with an AGG/AGA codon content of 8–10% were only marginally expressed, whereas the peanut profilin Ara h 5, with an AGG/AGA codon content of only 0.8%, was efficiently expressed in these cells. Hence, by using modified BL21(DE3) E. coli cells, namely BL21-CodonPlus(DE3)-RIL cells (Stratagene) with extra copies of E. coli argU, ileY, and leuW tRNA genes, it was possible to attain high-level expression of the proteins affected by rare codon usage. IPTG-induced expression of several recombinant peanut allergens, such as Ara h 1, Ara h 2, and Ara h 6, was greatly increased in these special cells compared to the expression yield achieved by conventional E. coli hosts. The purification of the soluble and the insoluble fraction of Ara h 2 was performed by metal-affinity chromatography and yielded a total of about 30 mg (His)₁₀-tagged recombinant protein per liter of culture of transformed BL21(DE3)CodonPlus-RIL cells. This is over 100 times more than achieved by production of Ara h 2 in conventional BL21(DE3) cells.     

The close proximity ofEscherichia coli genes: Consequences for stop codon and synonymous codon use

It is shown that synonymous codon usage is less biased in favor of those codons preferred by highly expressed genes at the end ofEscherichia coli genes than in the middle. This appears to be due to the close proximity of manyE. coli genes. It is shown that a substantial number of genes overlap either the Shine-Dalgarno sequence or the coding sequence of the next gene on the chromosome and that the codons that overlap have lower synonymous codon bias than those which do not. It is also shown that there is an increase in the frequency of A-ending codons, and a decrease in the frequency of G-ending codons at the end ofE. coli genes that lie close to another gene. It is suggested that these trends in composition could be associated with selection against the formation of mRNA secondary structure near the start of the next gene on the chromosome. Stop codon use is also affected by the close proximity of genes; many genes are forced to use TGA and TAG stop codons because they terminate either within the Shine-Dalgarno or coding sequence of the next gene on the chromosome. The implications these results have for the evolution of synonymous codon use are discussed.     

Codon Usage in Plastid Genes Is Correlated with Context, Position Within the Gene, and Amino Acid Content

Highly expressed plastid genes display codon adaptation, which is defined as a bias toward a set of codons which are complementary to abundant tRNAs. This type of adaptation is similar to what is observed in highly expressed Escherichia coli genes and is probably the result of selection to increase translation efficiency. In the current work, the codon adaptation of plastid genes is studied with regard to three specific features that have been observed in E. coli and which may influence translation efficiency. These features are (1) a relatively low codon adaptation at the 5′ end of highly expressed genes, (2) an influence of neighboring codons on codon usage at a particular site (codon context), and (3) a correlation between the level of codon adaptation of a gene and its amino acid content. All three features are found in plastid genes. First, highly expressed plastid genes have a noticeable decrease in codon adaptation over the first 10–20 codons. Second, for the twofold degenerate NNY codon groups, highly expressed genes have an overall bias toward the NNC codon, but this is not observed when the 3′ neighboring base is a G. At these sites highly expressed genes are biased toward NNT instead of NNC. Third, plastid genes that have higher codon adaptations also tend to have an increased usage of amino acids with a high G + C content at the first two codon positions and GNN codons in particular. The correlation between codon adaptation and amino acid content exists separately for both cytosolic and membrane proteins and is not related to any obvious functional property. It is suggested that at certain sites selection discriminates between nonsynonymous codons based on translational, not functional, differences, with the result that the amino acid sequence of highly expressed proteins is partially influenced by selection for increased translation efficiency. Received: 21 July 1999 / Accepted: 5 November 1999     

Enhanced expression of PCV2 capsid protein in Escherichia coli and Lactococcus lactis by codon optimization

Capsid protein (Cap) of porcine circovirus type 2 (PCV2) encoded by orf2 is a main structural protein with strong immunoreactivity. However, capsid protein is expressed poorly in prokaryotic organisms because of differences in codon usage. In this study, we introduce 24 synonymous mutations into orf2 by mutagenic primers and overlap extension polymerase chain reaction (OE-PCR) technique. Fourteen rare codons of orf2 were replaced with preferable codons used in Escherichia coli cells. Moreover, the nuclear localization signal (NLS) region rich in rare codon clusters at the 5′ end was deleted. The codon-optimized genes demonstrated higher levels of expression compared with wild-type genes. The influence of rare codons on the gene expression was eliminated by mutation. Western blot analysis confirmed the immunoreactivity of the proteins expressed by mutated genes. Further testing demonstrated that the mutated genes were also expressed successfully in Lactococcus lactis NZ9000. The immunologically active Cap proteins produced by recombinant strains have the potential applications for serological diagnostic assays and vaccine development against PCV2-associated diseases.     

Effect of codon adaptation on codon-level and gene-level translation efficiency in vivo

Background

There is a significant difference between synonymous codon usage in many organisms, and it is known that codons used more frequently generally showed efficient decoding rate. At the gene level, however, there are conflicting reports on the existence of a correlation between codon adaptation and translation efficiency, even in the same organism.

Results

To resolve this issue, we cultured Escherichia coli under conditions designed to maintain constant levels of mRNA and protein and subjected the cells to ribosome profiling (RP) and mRNA-seq analyses. We showed that the RP results correlated more closely with protein levels generated under similar culture conditions than with the mRNA abundance from the mRNA-seq. Our result indicated that RP/mRNA ratio could be used as a measure of translation efficiency at gene level. On the other hand, the RP data showed that codon-specific ribosome density at the decoding site negatively correlated with codon usage, consistent with the hypothesis that preferred codons display lower ribosome densities due to their faster decoding rate. However, highly codon-adapted genes showed higher ribosome densities at the gene level, indicating that the efficiency of translation initiation, rather than higher elongation efficiency of preferred codons, exerted a greater effect on ribosome density and thus translation efficiency.

Conclusions

These findings indicate that evolutionary pressure on highly expressed genes influenced both codon bias and translation initiation efficiency and therefore explains contradictory findings that codon usage bias correlates with translation efficiency of native genes, but not with the artificially created gene pool, which was not subjected to evolution pressure.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1115) contains supplementary material, which is available to authorized users.