Codon usage determines translation rate in Escherichia coli

Ribosomal mutations affecting the translation of genes that use non-optimal codons

Genes that are laterally acquired by a new host species often contain codons that are non-optimal to the tRNA repertoire of the new host, which may lead to insufficient translational levels. Inefficient translation can be overcome by different mechanisms, such as incremental amelioration of the coding sequence, compensatory mutations in the regulatory sequences leading to increased transcription or increase in gene copy number. However, there is also a possibility that ribosomal mutations can improve the expression of such genes. To test this hypothesis, we examined the effects of point mutations in the endogenous ribosomal proteins S12 and S5 in Escherichia coli, which are known to be involved in the decoding of the mRNA, on the efficiency of translation of exogenous genes that use non-optimal codons, in vivo. We show that an S12 mutant in E. coli is able to express exogenous genes, with non-optimal codons, to higher levels than the wild-type, and explore the mechanisms underlying this phenomenon in this mutant. Our results suggest that the transient emergence of mutants that allow efficient expression of exogenous genes with non-optimal codons could also increase the chances of fixation of laterally transferred genes.

Codon-by-codon modulation of translational speed and accuracy via mRNA folding

Secondary structure in mRNAs modulates the speed of protein synthesis codon-by-codon to improve accuracy at important sites while ensuring high speed elsewhere.

Rapid cell growth demands fast protein translational elongation to alleviate ribosome shortage. However, speedy elongation undermines translational accuracy because of a mechanistic tradeoff. Here we provide genomic evidence in budding yeast and mouse embryonic stem cells that the efficiency–accuracy conflict is alleviated by slowing down the elongation at structurally or functionally important residues to ensure their translational accuracies while sacrificing the accuracy for speed at other residues. Our computational analysis in yeast with codon resolution suggests that mRNA secondary structures serve as elongation brakes to control the speed and hence the fidelity of protein translation. The position-specific effect of mRNA folding on translational accuracy is further demonstrated experimentally by swapping synonymous codons in a yeast transgene. Our findings explain why highly expressed genes tend to have strong mRNA folding, slow translational elongation, and conserved protein sequences. The exquisite codon-by-codon translational modulation uncovered here is a testament to the power of natural selection in mitigating efficiency–accuracy conflicts, which are prevalent in biology.

Protein synthesis by ribosomal translation is a vital cellular process, but our understanding of its regulation has been poor. Because the number of ribosomes in the cell is limited, rapid growth relies on fast translational elongation. The accuracy of translation must also be maintained, and in an ideal scenario, both speed and accuracy should be maximized to sustain rapid and productive growth. However, existing data suggest a tradeoff between speed and accuracy, making it impossible to simultaneously maximize both. A potential solution is slowing the elongation at functionally or structurally important sites to ensure their translational accuracies, while sacrificing accuracy for speed at other sites. Here, we show that budding yeast and mouse embryonic stem cells indeed use this strategy. We discover that a codon-by-codon adaptive modulation of translational elongation is accomplished by mRNA secondary structures, which serve as brakes to control the elongation speed and hence translational fidelity. Our findings explain why highly expressed genes tend to have strong mRNA folding, slow translational elongation, and conserved protein sequences. The exquisite translational modulation reflects the power of natural selection in mitigating efficiency–accuracy conflicts, and our study offers a general framework for analyzing similar conflicts, which are widespread in biology.

The effects of codon context on in vivo translation speed

We developed a bacterial genetic system based on translation of the his operon leader peptide gene to determine the relative speed at which the ribosome reads single or multiple codons in vivo. Low frequency effects of so-called "silent" codon changes and codon neighbor (context) effects could be measured using this assay. An advantage of this system is that translation speed is unaffected by the primary sequence of the His leader peptide. We show that the apparent speed at which ribosomes translate synonymous codons can vary substantially even for synonymous codons read by the same tRNA species. Assaying translation through codon pairs for the 5'- and 3'- side positioning of the 64 codons relative to a specific codon revealed that the codon-pair orientation significantly affected in vivo translation speed. Codon pairs with rare arginine codons and successive proline codons were among the slowest codon pairs translated in vivo. This system allowed us to determine the effects of different factors on in vivo translation speed including Shine-Dalgarno sequence, rate of dipeptide bond formation, codon context, and charged tRNA levels.

Redundancy of the genetic code enables translational pausing

The codon redundancy (“degeneracy”) found in protein-coding regions of mRNA also prescribes Translational Pausing (TP). When coupled with the appropriate interpreters, multiple meanings and functions are programmed into the same sequence of configurable switch-settings. This additional layer of Ontological Prescriptive Information (PI_o) purposely slows or speeds up the translation-decoding process within the ribosome. Variable translation rates help prescribe functional folding of the nascent protein. Redundancy of the codon to amino acid mapping, therefore, is anything but superfluous or degenerate. Redundancy programming allows for simultaneous dual prescriptions of TP and amino acid assignments without cross-talk. This allows both functions to be coincident and realizable. We will demonstrate that the TP schema is a bona fide rule-based code, conforming to logical code-like properties. Second, we will demonstrate that this TP code is programmed into the supposedly degenerate redundancy of the codon table. We will show that algorithmic processes play a dominant role in the realization of this multi-dimensional code.

Understanding the influence of codon translation rates on cotranslational protein folding

Protein domains can fold into stable tertiary structures while they are synthesized by the ribosome in a process known as cotranslational folding. If a protein does not fold cotranslationally, however, it has the opportunity to do so post-translationally, that is, after the nascent chain has been fully synthesized and released from the ribosome. The rate at which a ribosome adds an amino acid encoded by a particular codon to the elongating nascent chain can vary significantly and is called the codon translation rate. Recent experiments have illustrated the profound impact that codon translation rates can have on the cotranslational folding process and the acquisition of function by nascent proteins. Synonymous codon mutations in an mRNA molecule change the chemical identity of a codon and its translation rate without changing the sequence of the synthesized protein. This change in codon translation rate can, however, cause a nascent protein to malfunction as a result of cotranslational misfolding. In some situations, such dysfunction can have profound implications; for example, it can alter the substrate specificity of an ABC transporter protein, resulting in patients who are nonresponsive to chemotherapy treatment. Thus, codon translation rates are crucial in coordinating protein folding in a cellular environment and can affect downstream cellular processes that depend on the proper functioning of newly synthesized proteins. As the importance of codon translation rates makes clear, a necessary aspect of fully understanding cotranslational folding lies in considering the kinetics of the process in addition to its thermodynamics. In this Account, we examine the contributions that have been made to elucidating the mechanisms of cotranslational folding by using the theoretical and computational tools of chemical kinetics, molecular simulations, and systems biology. These efforts have extended our ability to understand, model, and predict the influence of codon translation rates on cotranslational protein folding and misfolding. The application of such approaches to this important problem is creating a framework for making quantitative predictions of the impact of synonymous codon substitutions on cotranslational folding that has led to a novel hypothesis regarding the role of fast-translating codons in coordinating cotranslational folding. In addition, it is providing new insights into proteome-wide cotranslational folding behavior and making it possible to identify potential molecular mechanisms by which molecular chaperones can influence such behavior during protein synthesis. As we discuss in this Account, bringing together these theoretical developments with experimental approaches is increasingly helping answer fundamental questions about the nature of nascent protein folding on the ribosome.

Strain-specific parallel evolution drives short-term diversification during Pseudomonas aeruginosa biofilm formation

Generation of genetic diversity is a prerequisite for bacterial evolution and adaptation. Short-term diversification and selection within populations is, however, largely uncharacterised, as existing studies typically focus on fixed substitutions. Here, we use whole-genome deep-sequencing to capture the spectrum of mutations arising during biofilm development for two Pseudomonas aeruginosa strains. This approach identified single nucleotide variants with frequencies from 0.5% to 98.0% and showed that the clinical strain 18A exhibits greater genetic diversification than the type strain PA01, despite its lower per base mutation rate. Mutations were found to be strain specific: the mucoid strain 18A experienced mutations in alginate production genes and a c-di-GMP regulator gene; while PA01 acquired mutations in PilT and PilY1, possibly in response to a rapid expansion of a lytic Pf4 bacteriophage, which may use type IV pili for infection. The Pf4 population diversified with an evolutionary rate of 2.43 × 10(-3) substitutions per site per day, which is comparable to single-stranded RNA viruses. Extensive within-strain parallel evolution, often involving identical nucleotides, was also observed indicating that mutation supply is not limiting, which was contrasted by an almost complete lack of noncoding and synonymous mutations. Taken together, these results suggest that the majority of the P. aeruginosa genome is constrained by negative selection, with strong positive selection acting on an accessory subset of genes that facilitate adaptation to the biofilm lifecycle. Long-term bacterial evolution is known to proceed via few, nonsynonymous, positively selected mutations, and here we show that similar dynamics govern short-term, within-population bacterial diversification.

Coupling between protein level selection and codon usage optimization in the evolution of bacteria and archaea

The relationship between the selection affecting codon usage and selection on protein sequences of orthologous genes in diverse groups of bacteria and archaea was examined by using the Alignable Tight Genome Clusters database of prokaryote genomes. The codon usage bias is generally low, with 57.5% of the gene-specific optimal codon frequencies (F_opt) being below 0.55. This apparent weak selection on codon usage contrasts with the strong purifying selection on amino acid sequences, with 65.8% of the gene-specific dN/dS ratios being below 0.1. For most of the genomes compared, a limited but statistically significant negative correlation between F_opt and dN/dS was observed, which is indicative of a link between selection on protein sequence and selection on codon usage. The strength of the coupling between the protein level selection and codon usage bias showed a strong positive correlation with the genomic GC content. Combined with previous observations on the selection for GC-rich codons in bacteria and archaea with GC-rich genomes, these findings suggest that selection for translational fine-tuning could be an important factor in microbial evolution that drives the evolution of genome GC content away from mutational equilibrium. This type of selection is particularly pronounced in slowly evolving, “high-status” genes. A significantly stronger link between the two aspects of selection is observed in free-living bacteria than in parasitic bacteria and in genes encoding metabolic enzymes and transporters than in informational genes. These differences might reflect the special importance of translational fine-tuning for the adaptability of gene expression to environmental changes. The results of this work establish the coupling between protein level selection and selection for translational optimization as a distinct and potentially important factor in microbial evolution.

Selection affects the evolution of microbial genomes at many levels, including both the structure of proteins and the regulation of their production. Here we demonstrate the coupling between the selection on protein sequences and the optimization of codon usage in a broad range of bacteria and archaea. The strength of this coupling varies over a wide range and strongly and positively correlates with the genomic GC content. The cause(s) of the evolution of high GC content is a long-standing open question, given the universal mutational bias toward AT. We propose that optimization of codon usage could be one of the key factors that determine the evolution of GC-rich genomes. This work establishes the coupling between selection at the level of protein sequence and at the level of codon choice optimization as a distinct aspect of genome evolution.

Codon usage bias of the phosphoprotein gene of spring viraemia of carp virus and high codon adaptation to the host

In this study, we calculated the relative synonymous codon usage (RSCU) value and the effective number of codons (ENC) value to carry out principal component analysis (PCA) and correlation analysis of the codon usage pattern of the phosphoprotein gene (P gene) of spring viraemia of carp virus (SVCV). The synonymous codon usage pattern in P genes is geography-specific, based on PCA analysis. The high correlation between (G + C)1,2 % and (G + C)3 % suggests that mutational pressure rather than natural selection is the main factor that determines the codon usage and base components in P genes. At least 40 out of 59 synonymous codons are similarly selected in all functional genes within five complete SVCV genomes, and the hosts based on the RSCU data. These results not only provide insight into variations in the codon usage pattern of SVCV but also may help in understanding the processes governing the evolution of SVCV.

Investigating the role of site specific synonymous variation in disease association studies

Synonymous codon changes may not always be neutral indicating their significance in disease association studies, which is almost always overlooked. Synonymous substitutions may affect protein-folding rates leading to protein misfolding and aggregation. Genome wide analysis of 2301 mitochondrial genomes is performed to evaluate the significance of synonymous codons in disease association studies. The analysis revealed usage of rare codons at several sites in mitochondrial genes with rare codon usage higher for hydrophobic amino acids. The analysis suggests that variation data in association studies should be analyzed using site-specific codon usage values to infer the potential phenotypic impact of synonymous changes.

Expanding Anfinsen's principle: contributions of synonymous codon selection to rational protein design

Anfinsen’s principle asserts that all information required to specify the structure of a protein is encoded in its amino acid sequence. However, during protein synthesis by the ribosome, the N-terminus of the nascent chain can begin to fold before the C-terminus is available. We tested whether this cotranslational folding can alter the folded structure of an encoded protein in vivo, versus the structure formed when refolded in vitro. We designed a fluorescent protein consisting of three half-domains, where the N- and C-terminal half-domains compete with each other to interact with the central half-domain. The outcome of this competition determines the fluorescence properties of the resulting folded structure. Upon refolding after chemical denaturation, this protein produced equimolar amounts of the N- and C-terminal folded structures, respectively. In contrast, translation in Escherichia coli resulted in a 2-fold enhancement in the formation of the N-terminal folded structure. Rare synonymous codon substitutions at the 5′ end of the C-terminal half-domain further increased selection for the N-terminal folded structure. These results demonstrate that the rate at which a nascent protein emerges from the ribosome can specify the folded structure of a protein.

Dependency of codon usage on protein sequence patterns: a statistical study

BACKGROUND

Codon degeneracy and codon usage by organisms is an interesting and challenging problem. Researchers demonstrated the relation between codon usage and various functions or properties of genes and proteins, such as gene regulation, translation rate, translation efficiency, mRNA stability, splicing, and protein domains. Researchers usually represent segments of proteins responsible for specific functions or structures in a family of proteins as sequence patterns or motifs. We asked the question if organisms use the same codons in pattern segments as compared to the rest of the sequence.

METHODS

We used the likelihood ratio test, Pearson's chi-squared test, and mutual information to compare these two codon usages.

RESULTS

We showed that codon usage, in segments of genes that code for a given pattern or motif in a group of proteins, varied from the rest of the gene. The codon usage in these segments was not random. Amino acids with larger number of codons used more specific codon ratios in these segments. We studied the number of amino acids in the pattern (pattern length). As patterns got longer, there was a slight decrease in the fraction of patterns with significant different codon usage in the pattern region as compared to codon usage in the gene region. We defined a measure of specificity of protein patterns, and studied its relation to the codon usage. The difference in the codon usage between pattern region and gene region, was less for the patterns with higher specificity.

CONCLUSIONS

We provided a hypothesis that there are segments on genes that affect the codon usage and thus influence protein translation speed, and these regions are the regions that code protein pattern regions.

Hepatitis A virus: host interactions, molecular epidemiology and evolution

Infection with hepatitis A virus (HAV) is the commonest viral cause of liver disease and presents an important public health problem worldwide. Several unique HAV properties and molecular mechanisms of its interaction with host were recently discovered and should aid in clarifying the pathogenesis of hepatitis A. Genetic characterization of HAV strains have resulted in the identification of different genotypes and subtypes, which exhibit a characteristic worldwide distribution. Shifts in HAV endemicity occurring in different parts of the world, introduction of genetically diverse strains from geographically distant regions, genotype displacement observed in some countries and population expansion detected in the last decades of the 20th century using phylogenetic analysis are important factors contributing to the complex dynamics of HAV infections worldwide. Strong selection pressures, some of which, like usage of deoptimized codons, are unique to HAV, limit genetic variability of the virus. Analysis of subgenomic regions has been proven useful for outbreak investigations. However, sharing short sequences among epidemiologically unrelated strains indicates that specific identification of HAV strains for molecular surveillance can be achieved only using whole-genome sequences. Here, we present up-to-date information on the HAV molecular epidemiology and evolution, and highlight the most relevant features of the HAV-host interactions.

Arrest peptides: cis-acting modulators of translation

Each peptide bond of a protein is generated at the peptidyl transferase center (PTC) of the ribosome and then moves through the exit tunnel, which accommodates ever-changing segments of ≈ 40 amino acids of newly translated polypeptide. A class of proteins, called ribosome arrest peptides, contains specific sequences of amino acids (arrest sequences) that interact with distinct components of the PTC-exit tunnel region of the ribosome and arrest their own translation continuation, often in a manner regulated by environmental cues. Thus, the ribosome that has translated an arrest sequence is inactivated for peptidyl transfer, translocation, or termination. The stalled ribosome then changes the configuration or localization of mRNA, resulting in specific biological outputs, including regulation of the target gene expression and downstream events of mRNA/polypeptide maturation or localization. Living organisms thus seem to have integrated potentially harmful arrest sequences into elaborate regulatory mechanisms to express genetic information in productive directions.

Control of ribosome traffic by position-dependent choice of synonymous codons

Messenger RNA (mRNA) encodes a sequence of amino acids by using codons. For most amino acids, there are multiple synonymous codons that can encode the amino acid. The translation speed can vary from one codon to another, thus there is room for changing the ribosome speed while keeping the amino acid sequence and hence the resulting protein. Recently, it has been noticed that the choice of the synonymous codon, via the resulting distribution of slow- and fast-translated codons, affects not only on the average speed of one ribosome translating the mRNA but also might have an effect on nearby ribosomes by affecting the appearance of 'traffic jams' where multiple ribosomes collide and form queues. To test this 'context effect' further, we here investigate the effect of the sequence of synonymous codons on the ribosome traffic by using a ribosome traffic model with codon-dependent rates, estimated from experiments. We compare the ribosome traffic on wild-type (WT) sequences and sequences where the synonymous codons were swapped randomly. By simulating translation of 87 genes, we demonstrate that the WT sequences, especially those with a high bias in codon usage, tend to have the ability to reduce ribosome collisions, hence optimizing the cellular investment in the translation apparatus. The magnitude of such reduction of the translation time might have a significant impact on the cellular growth rate and thereby have importance for the survival of the species.

Translation of CGA codon repeats in yeast involves quality control components and ribosomal protein L1

Translation of CGA codon repeats in the yeast Saccharomyces cerevisiae is inefficient, resulting in dose-dependent reduction in expression and in production of an mRNA cleavage product, indicative of a stalled ribosome. Here, we use genetics and translation inhibitors to understand how ribosomes respond to CGA repeats. We find that CGA codon repeats result in a truncated polypeptide that is targeted for degradation by Ltn1, an E3 ubiquitin ligase involved in nonstop decay, although deletion of LTN1 does not improve expression downstream from CGA repeats. Expression downstream from CGA codons at residue 318, but not at residue 4, is improved by deletion of either ASC1 or HEL2, previously implicated in inhibition of translation by polybasic sequences. Thus, translation of CGA repeats likely causes ribosomes to stall and exploits known quality control systems. Expression downstream from CGA repeats at amino acid 4 is improved by paromomycin, an aminoglycoside that relaxes decoding specificity. Paromomycin has no effect if native tRNA(Arg(ICG)) is highly expressed, consistent with the idea that failure to efficiently decode CGA codons might occur in part due to rejection of the cognate tRNA(Arg(ICG)). Furthermore, expression downstream from CGA repeats is improved by inactivation of RPL1B, one of two genes encoding the universally conserved ribosomal protein L1. The effects of rpl1b-Δ and of either paromomycin or tRNA(Arg(ICG)) on CGA decoding are additive, suggesting that the rpl1b-Δ mutant suppresses CGA inhibition by means other than increased acceptance of tRNA(Arg(ICG)). Thus, inefficient decoding of CGA likely involves at least two independent defects in translation.

Translational sensitivity of the Escherichia coli genome to fluctuating tRNA availability

The synthesis of protein from messenger RNA during translation is a highly dynamic process that plays a key role in controlling the efficiency and fidelity of genome-wide protein expression. The availability of aminoacylated transfer RNA (tRNA) is a major factor influencing the speed of ribosomal movement, which depending on codon choices, varies considerably along a transcript. Furthermore, it has been shown experimentally that tRNA availability can vary significantly under different growth and stress conditions, offering the cell a way to adapt translational dynamics across the genome. Existing models of translation have neglected fluctuations of tRNA pools, instead assuming fixed tRNA availabilities over time. This has lead to an incomplete understanding of this process. Here, we show for the entire Escherichia coli genome how and to what extent translational speed profiles, which capture local aspects of translational elongation, respond to measured shifts in tRNA availability. We find that translational profiles across the genome are affected to differing degrees, with genes that are essential or related to fundamental processes such as translation, being more robust than those linked to regulation. Furthermore, we reveal how fluctuating tRNA availability influences profiles of specific sequences known to play a significant role in translational control of gene expression.

Analysis of the use of codon pairs in the HE gene of the ISA virus shows a correlation between bias in HPR codon-pair use and mortality rates caused by the virus

BACKGROUND

Segment 6 of the ISA virus codes for hemoagglutinin-esterase (HE). This segment is highly variable, with more than 26 variants identified. The major variation is observed in what is called the high polymorphism region (HPR). The role of the different HPR zones in the viral cycle or evolution remains unknown. However viruses that present the HPR0 are avirulent, while viruses with important deletions in this region have been responsible for outbreaks with high mortality rates. In this work, using bioinformatic tools, we examined the influence of different HPRs on the adaptation of HE genes to the host translational machinery and the relationship to observed virulence.

METHODS

Translational efficiency of HE genes and their HPR were estimated analyzing codon-pair bias (CPB), adaptation to host codon use (codon adaptation index-CAI) and the adaptation to available tRNAs (tAI). These values were correlated with reported mortality for the respective ISA virus and the ΔG of RNA folding. tRNA abundance was inferred from tRNA gene numbers identified in the Salmo salar genome using tRNAScan-SE. Statistical correlation between data was performed using a non-parametric test.

RESULTS

We found that HPR0 contains zones with codon pairs of low frequency and low availability of tRNA with respect to salmon codon-pair usage, suggesting that HPR modifies HE translational efficiency. Although calculating tAI was impossible because one third of tRNAs (~60.000) were tRNA-ala, translational efficiency measured by CPB shows that as HPR size increases, the CPB value of the HE gene decreases (P = 2x10⁻⁷, ρ = -0.675, n = 63) and that these values correlate positively with the mortality rates caused by the virus (ρ = 0.829, P = 2x10⁻⁷, n = 11). The mortality associated with different virus isolates or their corresponding HPR sizes were not related with the ΔG of HPR RNA folding, suggesting that the secondary structure of HPR RNA does not modify virulence.

CONCLUSIONS

Our results suggest that HPR size affects the efficiency of gene translation, which modulates the virulence of the virus by a mechanism similar to that observed in production of live attenuated vaccines through deoptimization of codon-pair usage.

The effects of the synonymous codon usage and tRNA abundance on protein folding of the 3C protease of foot-and-mouth disease virus

The 3C protease of foot-and-mouth disease virus (FMDV) has a conserved amino acid sequence and is responsible for most cleavage in the viral polyprotein. The effects of the synonymous codon usage of FMDV 3C gene and tRNA abundance of the hosts on shaping different folding units (α-helix, β-strand and the coil) in the 3C protease were analyzed based on the structural information of the FMDV 3C protease from Protein Data Bank (PDB: 2BHG) and 210 genes of 3C for all serotypes of FMDV. The strong correlation between some codons usage and the specific folding unit in the FMDV 3C protease is found. As for the effect of translation speed caused by tRNA abundance on the formation of folding units, the codon positions with lowly abundant tRNA scatter in the 3C gene and there is the obvious fluctuation of tRNA abundance locating in the transition boundaries from the β-strand to the α-helix and the coil, respectively. However, codon positions with lowly abundant tRNA clustering into these boundaries are not found, suggesting that the adjustment of translation speed is likely also achieved by the single codon position with low tRNA abundance rather than a cluster. The observations can provide the information for insight into the role of the synonymous codon usage in the formation of 3C protease of FMDV and effect of the tRNA abundance of the hosts on this formation of protease.

Estimating selection on synonymous codon usage from noisy experimental data

A key goal in molecular evolution is to extract mechanistic insights from signatures of selection. A case study is codon usage, where despite many recent advances and hypotheses, two longstanding problems remain: the relative contribution of selection and mutation in determining codon frequencies and the relative contribution of translational speed and accuracy to selection. The relevant targets of selection--the rate of translation and of mistranslation of a codon per unit time in the cell--can only be related to mechanistic properties of the translational apparatus if the number of transcripts per cell is known, requiring use of gene expression measurements. Perhaps surprisingly, different gene-expression data sets yield markedly different estimates of selection. We show that this is largely due to measurement noise, notably due to differences between studies rather than instrument error or biological variability. We develop an analytical framework that explicitly models noise in expression in the context of the population-genetic model. Estimates of mutation and selection strength in budding yeast produced by this method are robust to the expression data set used and are substantially higher than estimates using a noise-blind approach. We introduce per-gene selection estimates that correlate well with previous scoring systems, such as the codon adaptation index, while now carrying an evolutionary interpretation. On average, selection for codon usage in budding yeast is weak, yet our estimates show that genes range from virtually unselected to average per-codon selection coefficients above the inverse population size. Our analytical framework may be generally useful for distinguishing biological signals from measurement noise in other applications that depend upon measurements of gene expression.

In vivo biochemistry in bacterial cells using FRAP: insight into the translation cycle

In vivo measurements of the mobility and binding kinetics of cellular components are essential to fully understand the biochemical processes occurring inside cells. Here, we describe a fluorescence recovery after photobleaching-based method that can be easily implemented to the study of reaction-diffusion processes in live bacteria despite their small size. We apply this method to provide new, to our knowledge, quantitative insight into multiple aspects of the bacterial translation cycle by measuring the binding kinetics and the micrometer-scale diffusive properties of the 50S ribosomal subunit in live Caulobacter cells. From our measurements, we infer that 70% of 50S subunits are engaged in translation and display, on average, limited motion on the micrometer scale, consistent with little mixing of transcripts undergoing translation. We also extract the average rate constants for the binding of 50S subunits to 30S initiation complexes during initiation and for their release from mRNAs when translation is completed. From this, we estimate the average time of protein synthesis and the average search time of 50S subunits before they engage in the next initiation event. Additionally, our experiments suggest that so-called free 50S subunits do not diffuse freely; instead their mobility is significantly slowed down, possibly through transient associations with mRNA.

Quasi-cellular systems: stochastic simulation analysis at nanoscale range

BACKGROUND

The wet-lab synthesis of the simplest forms of life (minimal cells) is a challenging aspect in modern synthetic biology. Quasi-cellular systems able to produce proteins directly from DNA can be obtained by encapsulating the cell-free transcription/translation system PURESYSTEM(PS) in liposomes. It is possible to detect the intra-vesicle protein production using DNA encoding for GFP and monitoring the fluorescence emission over time. The entrapment of solutes in small-volume liposomes is a fundamental open problem. Stochastic simulation is a valuable tool in the study of biochemical reaction at nanoscale range. QDC (Quick Direct-Method Controlled), a stochastic simulation software based on the well-known Gillespie's SSA algorithm, was used. A suitable model formally describing the PS reactions network was developed, to predict, from inner species concentrations (very difficult to measure in small-volumes), the resulting fluorescence signal (experimentally observable).

RESULTS

Thanks to suitable features specific of QDC, we successfully formalized the dynamical coupling between the transcription and translation processes that occurs in the real PS, thus bypassing the concurrent-only environment of Gillespie's algorithm. Simulations were firstly performed for large liposomes (2.67µm of diameter) entrapping the PS to synthetize GFP. By varying the initial concentrations of the three main classes of molecules involved in the PS (DNA, enzymes, consumables), we were able to stochastically simulate the time-course of GFP-production. The sigmoid fit of the GFP-production curves allowed us to extract three quantitative parameters which are significantly dependent on the various initial states. Then we extended this study for small-volume liposomes (575 nm of diameter), where it is more complex to infer the intra-vesicle composition, due to the expected anomalous entrapment phenomena. We identified almost two extreme states that are forecasted to give rise to significantly different experimental observables.

CONCLUSIONS

The present work is the first one describing in the detail the stochastic behavior of the PS. Thanks to our results, an experimental approach is now possible, aimed at recording the GFP production kinetics in very small micro-emulsion droplets or liposomes, and inferring, by using the simulation as a reverse-engineering procedure, the internal solutes distribution, and shed light on the still unknown forces driving the entrapment phenomenon.

Combinational effect of mutational bias and translational selection for translation efficiency in tomato (Solanum lycopersicum) cv. Micro-Tom

We conducted a comprehensive analysis of codon usage bias (CUB) based on the available non-redundant full-length cDNA (nrFLcDNA) and expressed sequence tags (ESTs) data of cultivar Micro-Tom and evaluated the associations of observed CUB and measurements of transcriptional and translational effectiveness. The analysis presented in our study suggests a correlation, which is negative but highly correlated between Axis 1 and GC3s (r=-0.827, P<0.01), indicating that mutational bias has a significant and dominant repressive role to the choices of GC3. We also observed a strong positive correlation between codon adaptation index (CAI) and translational adaptation index (tAIg) (0.407, P<0.01), which demonstrates the facilitation of efficient translation by the optimal codon usage patterns of the highly expressed genes. We believe that the complete set of optimal codon usage patterns detected in this study will serve as a model to enhance the transgenesis in the studied cultivar of Solanum lycopersicum.

Molecular cloning, recombinant expression, and antimicrobial activity of EC-hepcidin3, a new four-cysteine hepcidin isoform from Epinephelus coioides

Hepcidin, a cysteine-rich antimicrobial peptide, is widespread in fish and shows multiple activities, including antimicrobial, antivirus, and antitumor. Here, a new four-cysteine hepcidin isoform gene, EC-hepcidin3, was cloned from the marine-cultured orange-spotted grouper (Epinephelus coioides). The complete cDNA sequence consisted of 603 bases with an open reading frame (ORF) of 270 bases. The genomic DNA sequence was composed of two introns and three exons, and its 312-bp upstream region had multiple putative transcription factor binding sites. Soluble recombinant protein EC-proHep3 containing a His-tag at the C-terminus was obtained from expression plasmid pET-28a/EC-proHep3 in Escherichia coli Rosetta. It was purified by immobilized metal affinity chromatography (IMAC), and it showed antibacterial activity in vitro. Kinetic studies indicated that recombinant EC-proHep3 has strong, rapid activity against Staphylococcus aureus and Pseudomonas stutzeri. The results indicate that EC-hepcidin3 might be an effective component in the innate immune system of groupers.

Sounds of silence: synonymous nucleotides as a key to biological regulation and complexity

Messenger RNA is a key component of an intricate regulatory network of its own. It accommodates numerous nucleotide signals that overlap protein coding sequences and are responsible for multiple levels of regulation and generation of biological complexity. A wealth of structural and regulatory information, which mRNA carries in addition to the encoded amino acid sequence, raises the question of how these signals and overlapping codes are delineated along non-synonymous and synonymous positions in protein coding regions, especially in eukaryotes. Silent or synonymous codon positions, which do not determine amino acid sequences of the encoded proteins, define mRNA secondary structure and stability and affect the rate of translation, folding and post-translational modifications of nascent polypeptides. The RNA level selection is acting on synonymous sites in both prokaryotes and eukaryotes and is more common than previously thought. Selection pressure on the coding gene regions follows three-nucleotide periodic pattern of nucleotide base-pairing in mRNA, which is imposed by the genetic code. Synonymous positions of the coding regions have a higher level of hybridization potential relative to non-synonymous positions, and are multifunctional in their regulatory and structural roles. Recent experimental evidence and analysis of mRNA structure and interspecies conservation suggest that there is an evolutionary tradeoff between selective pressure acting at the RNA and protein levels. Here we provide a comprehensive overview of the studies that define the role of silent positions in regulating RNA structure and processing that exert downstream effects on proteins and their functions.

Good codons, bad transcript: large reductions in gene expression and fitness arising from synonymous mutations in a key enzyme

Biased codon usage in protein-coding genes is pervasive, whereby amino acids are largely encoded by a specific subset of possible codons. Within individual genes, codon bias is stronger at evolutionarily conserved residues, favoring codons recognized by abundant tRNAs. Although this observation suggests an overall pattern of selection for translation speed and/or accuracy, other work indicates that transcript structure or binding motifs drive codon usage. However, our understanding of codon bias evolution is constrained by limited experimental data on the fitness effects of altering codons in functional genes. To bridge this gap, we generated synonymous variants of a key enzyme-coding gene in Methylobacterium extorquens. We found that mutant gene expression, enzyme production, enzyme activity, and fitness were all significantly lower than wild-type. Surprisingly, encoding the gene using only rare codons decreased fitness by 40%, whereas an allele coded entirely by frequent codons decreased fitness by more than 90%. Increasing gene expression restored mutant fitness to varying degrees, demonstrating that the fitness disadvantage of synonymous mutants arose from a lack of beneficial protein rather than costs of protein production. Protein production was negatively correlated with the frequency of motifs with high affinity for the anti-Shine-Dalgarno sequence, suggesting ribosome pausing as the dominant cause of low mutant fitness. Together, our data support the idea that, although a particular set of codons are favored on average across a genome, in an individual gene selection can either act for or against codons depending on their local context.

Effects of gene length on the dynamics of gene expression

In Escherichia coli, the nucleotide length of a gene is bound to affect its expression dynamics. From simulations of a stochastic model of gene expression at the nucleotide and codon levels we show that, within realistic parameter values, the nucleotide length affects RNA and protein mean levels, as well as the expected transient time for RNA and protein numbers to change, following a signal. Fluctuations in RNA and protein numbers are found to be minimized for a small range of lengths, which matches the means of the distributions of lengths found in E. coli of both essential and non-essential genes. The variance of the length distribution for essential genes is found to be smaller than for non-essential genes, implying that these distributions are far from random. Finally, gene lengths are shown to affect the kinetics of a genetic switch, namely, the correlation between temporal proteins numbers, the stability of the two noisy attractors of the switch, and how biased is the choice of noisy attractor. The stability increases with gene length due to increased 'memory' about the previous states of the switch. We argue that, by affecting the dynamics of gene expression and of genetic circuits, gene lengths are subject to selection.

Computational codon optimization of synthetic gene for protein expression

Background

The construction of customized nucleic acid sequences allows us to have greater flexibility in gene design for recombinant protein expression. Among the various parameters considered for such DNA sequence design, individual codon usage (ICU) has been implicated as one of the most crucial factors affecting mRNA translational efficiency. However, previous works have also reported the significant influence of codon pair usage, also known as codon context (CC), on the level of protein expression.

Results

In this study, we have developed novel computational procedures for evaluating the relative importance of optimizing ICU and CC for enhancing protein expression. By formulating appropriate mathematical expressions to quantify the ICU and CC fitness of a coding sequence, optimization procedures based on genetic algorithm were employed to maximize its ICU and/or CC fitness. Surprisingly, the in silico validation of the resultant optimized DNA sequences for Escherichia coli, Lactococcus lactis, Pichia pastoris and Saccharomyces cerevisiae suggests that CC is a more relevant design criterion than the commonly considered ICU.

Conclusions

The proposed CC optimization framework can complement and enhance the capabilities of current gene design tools, with potential applications to heterologous protein production and even vaccine development in synthetic biotechnology.

Speeding with control: codon usage, tRNAs, and ribosomes

Codon usage and tRNA abundance are critical parameters for gene synthesis. However, the forces determining codon usage bias within genomes and between organisms, as well as the functional roles of biased codon compositions, remain poorly understood. Similarly, the composition and dynamics of mature tRNA populations in cells in terms of isoacceptor abundances, and the prevalence and function of base modifications are not well understood. As we begin to decipher some of the rules that govern codon usage and tRNA abundances, it is becoming clear that these parameters are a way to not only increase gene expression, but also regulate the speed of ribosomal translation, the efficiency of protein folding, and the coordinated expression of functionally related gene families. Here, we discuss the importance of codon-anticodon interactions in translation regulation and highlight the contribution of non-random codon distributions and post-transcriptional base modifications to this regulation.

tRNAomics: tRNA gene copy number variation and codon use provide bioinformatic evidence of a new anticodon:codon wobble pair in a eukaryote

tRNA genes are interspersed throughout eukaryotic DNA, contributing to genome architecture and evolution in addition to translation of the transcriptome. Codon use correlates with tRNA gene copy number in noncomplex organisms including yeasts. Synonymous codons impact translation with various outcomes, dependent on relative tRNA abundances. Availability of whole-genome sequences allowed us to examine tRNA gene copy number variation (tgCNV) and codon use in four Schizosaccharomyces species and Saccharomyces cerevisiae. tRNA gene numbers vary from 171 to 322 in the four Schizosaccharomyces despite very high similarity in other features of their genomes. In addition, we performed whole-genome sequencing of several related laboratory strains of Schizosaccharomyces pombe and found tgCNV at a cluster of tRNA genes. We examined for the first time effects of wobble rules on correlation of tRNA gene number and codon use and showed improvement for S. cerevisiae and three of the Schizosaccharomyces species. In contrast, correlation in Schizosaccharomyces japonicus is poor due to markedly divergent tRNA gene content, and much worsened by the wobble rules. In japonicus, some tRNA iso-acceptor genes are absent and others are greatly reduced relative to the other yeasts, while genes for synonymous wobble iso-acceptors are amplified, indicating wobble use not apparent in any other eukaryote. We identified a subset of japonicus-specific wobbles that improves correlation of codon use and tRNA gene content in japonicus. We conclude that tgCNV is high among Schizo species and occurs in related laboratory strains of S. pombe (and expectedly other species), and tRNAome-codon analyses can provide insight into species-specific wobble decoding.

A max-plus model of ribosome dynamics during mRNA translation

We examine the dynamics of the translation stage of cellular protein production, in which ribosomes move uni-directionally along an mRNA strand, building amino acid chains as they go. We describe the system using a timed event graph-a class of Petri net useful for studying discrete events, which have to satisfy constraints. We use max-plus algebra to describe a deterministic version of the model, where the constraints represent steric effects which prevent more than one ribosome reading a given codon at a given time and delays associated with the availability of the different tRNAs. We calculate the protein production rate and density of ribosomes on the mRNA and find exact agreement between these analytical results and numerical simulations of the deterministic model, even in the case of heterogeneous mRNAs.

Effect of codon message on xylanase thermal activity

Because the genetic codon is known for degeneracy, its effect on enzyme thermal property is seldom investigated. A dataset was constructed for GH10 xylanase coding sequences and optimal temperatures for activity (T(opt)). Codon contents and relative synonymous codon usages were calculated and respectively correlated with the enzyme T(opt) values, which were used to describe the xylanase thermophilic tendencies without dividing them into two thermophilic and mesophilic groups. After analyses of codon content and relative synonymous codon usages were checked by the Bonferroni correction, we found five codons, with three (AUA, AGA, and AGG) correlating positively and two (CGU and AGC) correlating negatively with the T(opt) value. The three positive codons are purine-rich codons, and the two negative codons have A-ends. The two negative codons are pyridine-rich codons, and one has a C-end. Comparable with the codon C- and A-ending features, C- and A-content within mRNA correlated negatively and positively with the T(opt) value, respectively. Thereby, codons have effects on enzyme thermal property. When the issue is analyzed at the residual level, the effect of codon message is lost. The codons relating to enzyme thermal property are selected by thermophilic force at nucleotide level.

Hepatitis A virus evolution and the potential emergence of new variants escaping the presently available vaccines

Hepatitis A is the most common infection of the liver worldwide and is fecal-orally transmitted. Its incidence tends to decrease with improvements in hygiene conditions but at the same time its severity increases. Hepatitis A virus is the causative agent of acute hepatitis in humans and belongs to the Hepatovirus genus in the Picornaviridae family, and it has very unique characteristics. This article reviews some molecular and biological properties that allow the virus to live in a very quiescent way and to build an extremely stable capsid that is able to persist in and out of the body. Additionally, the relationship between the genomic composition and the structural and antigenic properties of the capsid is discussed, and the potential emergence of antigenic variants is evaluated from an evolutionary perspective.

Codon usage bias in prokaryotic pyrimidine-ending codons is associated with the degeneracy of the encoded amino acids

Synonymous codons are unevenly distributed among genes, a phenomenon termed codon usage bias. Understanding the patterns of codon bias and the forces shaping them is a major step towards elucidating the adaptive advantage codon choice can confer at the level of individual genes and organisms. Here, we perform a large-scale analysis to assess codon usage bias pattern of pyrimidine-ending codons in highly expressed genes in prokaryotes. We find a bias pattern linked to the degeneracy of the encoded amino acid. Specifically, we show that codon-pairs that encode two- and three-fold degenerate amino acids are biased towards the C-ending codon while codons encoding four-fold degenerate amino acids are biased towards the U-ending codon. This codon usage pattern is widespread in prokaryotes, and its strength is correlated with translational selection both within and between organisms. We show that this bias is associated with an improved correspondence with the tRNA pool, avoidance of mis-incorporation errors during translation and moderate stability of codon–anticodon interaction, all consistent with more efficient translation.

Recent advances in microbial production of δ-aminolevulinic acid and vitamin B12

δ-aminolevulinate (ALA) is an important intermediate involved in tetrapyrrole synthesis (precursor for vitamin B12, chlorophyll and heme) in vivo. It has been widely applied in agriculture and medicine. On account of many disadvantages of its chemical synthesis, microbial production of ALA has been received much attention as an alternative because of less expensive raw materials, low pollution, and high productivity. Vitamin B12, one of ALA derivatives, which plays a vital role in prevention of anaemia has also attracted intensive works. In this review, recent advances on the production of ALA and vitamin B12 with novel approaches such as whole-cell enzyme-transformation and metabolic engineering are described. Furthermore, the direction for future research and perspective are also summarized.

Large-scale analysis of conserved rare codon clusters suggests an involvement in co-translational molecular recognition events

MOTIVATION

An increasing amount of evidence from experimental and computational analysis suggests that rare codon clusters are functionally important for protein activity. Most of the studies on rare codon clusters were performed on a limited number of proteins or protein families. In the present study, we present the Sherlocc program and how it can be used for large scale protein family analysis of evolutionarily conserved rare codon clusters and their relation to protein function and structure. This large-scale analysis was performed using the whole Pfam database covering over 70% of the known protein sequence universe. Our program Sherlocc, detects statistically relevant conserved rare codon clusters and produces a user-friendly HTML output.

RESULTS

Statistically significant rare codon clusters were detected in a multitude of Pfam protein families. The most statistically significant rare codon clusters were predominantly identified in N-terminal Pfam families. Many of the longest rare codon clusters are found in membrane-related proteins which are required to interact with other proteins as part of their function, for example in targeting or insertion. We identified some cases where rare codon clusters can play a regulating role in the folding of catalytically important domains. Our results support the existence of a widespread functional role for rare codon clusters across species. Finally, we developed an online filter-based search interface that provides access to Sherlocc results for all Pfam families.

AVAILABILITY

The Sherlocc program and search interface are open access and are available at http://bcb.med.usherbrooke.ca

Balanced codon usage optimizes eukaryotic translational efficiency

Cellular efficiency in protein translation is an important fitness determinant in rapidly growing organisms. It is widely believed that synonymous codons are translated with unequal speeds and that translational efficiency is maximized by the exclusive use of rapidly translated codons. Here we estimate the in vivo translational speeds of all sense codons from the budding yeast Saccharomyces cerevisiae. Surprisingly, preferentially used codons are not translated faster than unpreferred ones. We hypothesize that this phenomenon is a result of codon usage in proportion to cognate tRNA concentrations, the optimal strategy in enhancing translational efficiency under tRNA shortage. Our predicted codon–tRNA balance is indeed observed from all model eukaryotes examined, and its impact on translational efficiency is further validated experimentally. Our study reveals a previously unsuspected mechanism by which unequal codon usage increases translational efficiency, demonstrates widespread natural selection for translational efficiency, and offers new strategies to improve synthetic biology.

Although an amino acid can be encoded by multiple synonymous codons, these codons are not used equally frequently in a genome. Biased codon usage is believed to improve translational efficiency because it is thought that preferentially used codons are translated faster than unpreferred ones. Surprisingly, we find similar translational speeds among synonymous codons. We show that translational efficiency is optimized by a previously unknown mechanism that relies on proportional use of codons according to their cognate tRNA concentrations. Our results provide important molecular details of protein translation, answer why codon usage is unequal, demonstrate widespread natural selection for translational efficiency, and can guide designs of synthetic genomes and cells with efficient translation systems.

The anti-Shine-Dalgarno sequence drives translational pausing and codon choice in bacteria

Protein synthesis by ribosomes takes place on a linear substrate but at variable speeds. Transient pausing of ribosomes can impact a variety of co-translational processes, including protein targeting and folding¹. These pauses are influenced by the sequence of the mRNA². Thus redundancy in the genetic code allows the same protein to be translated at different rates. However, our knowledge of both the position and the mechanism of translational pausing in vivo is highly limited. Here we present a genome-wide analysis of translational pausing in bacteria using ribosome profiling–deep sequencing of ribosome-protected mRNA fragments^3-5. This approach enables high-resolution measurement of ribosome density profiles along most transcripts at unperturbed, endogenous expression levels. Unexpectedly, we found that codons decoded by rare tRNAs do not lead to slow translation under nutrient-rich conditions. Instead, Shine-Dalgarno-(SD)⁶ like features within coding sequences cause pervasive translational pausing. Using an orthogonal ribosome^7,8 possessing an altered anti-SD sequence, we demonstrated that pausing is due to hybridization between mRNA and the 16S rRNA of the translating ribosome. In protein coding sequences, internal SD sequences are disfavoured, which leads to biased usage, avoiding codons and codon pairs that resemble canonical SD sites. Our results indicate that internal SD-like sequences are a major determinant of translation rates and a global driving force for the coding of bacterial genomes.

Genetic code redundancy and its influence on the encoded polypeptides

The genetic code is said to be redundant in that the same amino acid residue can be encoded by multiple, so-called synonymous, codons. If all properties of synonymous codons were entirely equivalent, one would expect that they would be equally distributed along protein coding sequences. However, many studies over the last three decades have demonstrated that their distribution is not entirely random. It has been postulated that certain codons may be translated by the ribosome faster than others and thus their non-random distribution dictates how fast the ribosome moves along particular segments of the mRNA. The reasons behind such segmental variability in the rates of protein synthesis, and thus polypeptide emergence from the ribosome, have been explored by theoretical and experimental approaches. Predictions of the relative rates at which particular codons are translated and their impact on the nascent chain have not arrived at unequivocal conclusions. This is probably due, at least in part, to variation in the basis for classification of codons as “fast” or “slow”, as well as variability in the number and types of genes and proteins analyzed. Recent methodological advances have allowed nucleotide-resolution studies of ribosome residency times in entire transcriptomes, which confirm the non-uniform movement of ribosomes along mRNAs and shed light on the actual determinants of rate control. Moreover, experiments have begun to emerge that systematically examine the influence of variations in ribosomal movement and the fate of the emerging polypeptide chain.

Engineering genes for predictable protein expression

The DNA sequence used to encode a polypeptide can have dramatic effects on its expression. Lack of readily available tools has until recently inhibited meaningful experimental investigation of this phenomenon. Advances in synthetic biology and the application of modern engineering approaches now provide the tools for systematic analysis of the sequence variables affecting heterologous expression of recombinant proteins. We here discuss how these new tools are being applied and how they circumvent the constraints of previous approaches, highlighting some of the surprising and promising results emerging from the developing field of gene engineering.

Synthesis of genetically engineered protein polymers (recombinamers) as an example of advanced self-assembled smart materials

In this chapter, we describe two methods for bio-producing recombinant repetitive polypeptide polymers for use in biomedical devices. These polymers, known as elastin-like recombinamers (ELRs), are derived from the repetition of selected amino acid domains of extracellular matrix proteins with the aim of recreating their mechanical and physiological features. The proteinaceous nature of ELRs allows us to make use of the natural biosynthetic machinery of heterologous hosts to express advanced and large polymers or "recombinamers." Despite the essentially unlimited possibilities for designing recombinamers, the production of synthetic genes to encode them should allow us to overcome the difficulties surrounding bioproduction of these non-natural monotonous DNA and protein sequences. The aim of this work is to supply the biotechnologist with fine-tuning methods to biosynthesize advanced self-assembled smart materials.

mRNA pseudoknot structures can act as ribosomal roadblocks

Several viruses utilize programmed ribosomal frameshifting mediated by mRNA pseudoknots in combination with a slippery sequence to produce a well defined stochiometric ratio of the upstream encoded to the downstream-encoded protein. A correlation between the mechanical strength of mRNA pseudoknots and frameshifting efficiency has previously been found; however, the physical mechanism behind frameshifting still remains to be fully understood. In this study, we utilized synthetic sequences predicted to form mRNA pseudoknot-like structures. Surprisingly, the structures predicted to be strongest lead only to limited frameshifting. Two-dimensional gel electrophoresis of pulse labelled proteins revealed that a significant fraction of the ribosomes were frameshifted but unable to pass the pseudoknot-like structures. Hence, pseudoknots can act as ribosomal roadblocks, prohibiting a significant fraction of the frameshifted ribosomes from reaching the downstream stop codon. The stronger the pseudoknot the larger the frameshifting efficiency and the larger its roadblocking effect. The maximal amount of full-length frameshifted product is produced from a structure where those two effects are balanced. Taking ribosomal roadblocking into account is a prerequisite for formulating correct frameshifting hypotheses.

The great screen anomaly--a new frontier in product discovery through functional metagenomics

Functional metagenomics, the study of the collective genome of a microbial community by expressing it in a foreign host, is an emerging field in biotechnology. Over the past years, the possibility of novel product discovery through metagenomics has developed rapidly. Thus, metagenomics has been heralded as a promising mining strategy of resources for the biotechnological and pharmaceutical industry. However, in spite of innovative work in the field of functional genomics in recent years, yields from function-based metagenomics studies still fall short of producing significant amounts of new products that are valuable for biotechnological processes. Thus, a new set of strategies is required with respect to fostering gene expression in comparison to the traditional work. These new strategies should address a major issue, that is, how to successfully express a set of unknown genes of unknown origin in a foreign host in high throughput. This article is an opinionating review of functional metagenomic screening of natural microbial communities, with a focus on the optimization of new product discovery. It first summarizes current major bottlenecks in functional metagenomics and then provides an overview of the general metagenomic assessment strategies, with a focus on the challenges that are met in the screening for, and selection of, target genes in metagenomic libraries. To identify possible screening limitations, strategies to achieve optimal gene expression are reviewed, examining the molecular events all the way from the transcription level through to the secretion of the target gene product.

Nascentome analysis uncovers futile protein synthesis in Escherichia coli

Although co-translational biological processes attract much attention, no general and easy method has been available to detect cellular nascent polypeptide chains, which we propose to call collectively a “nascentome.” We developed a method to selectively detect polypeptide portions of cellular polypeptidyl-tRNAs and used it to study the generality of the quality control reactions that rescue dead-end translation complexes. To detect nascent polypeptides, having their growing ends covalently attached to a tRNA, cellular extracts are separated by SDS-PAGE in two dimensions, first with the peptidyl-tRNA ester bonds preserved and subsequently after their in-gel cleavage. Pulse-labeled nascent polypeptides of Escherichia coli form a characteristic line below the main diagonal line, because each of them had contained a tRNA of nearly uniform size in the first-dimension electrophoresis but not in the second-dimension. The detection of nascent polypeptides, separately from any translation-completed polypeptides or degradation products thereof, allows us to follow their fates to gain deeper insights into protein biogenesis and quality control pathways. It was revealed that polypeptidyl-tRNAs were significantly stabilized in E. coli upon dysfunction of the tmRNA-ArfA ribosome-rescuing system, whose function had only been studied previously using model constructs. Our results suggest that E. coli cells are intrinsically producing aberrant translation products, which are normally eliminated by the ribosome-rescuing mechanisms.

Wobble base-pairing slows in vivo translation elongation in metazoans

In the universal genetic code, most amino acids can be encoded by multiple trinucleotide codons, and the choice among available codons can influence position-specific translation elongation rates. By using sequence-based ribosome profiling, we obtained transcriptome-wide profiles of in vivo ribosome occupancy as a function of codon identity in Caenorhabditis elegans and human cells. Particularly striking in these profiles was a universal trend of higher ribosome occupancy for codons translated via G:U wobble base-pairing compared with synonymous codons that pair with the same tRNA family using G:C base-pairing. These data support a model in which ribosomal translocation is slowed at wobble codon positions.

Mutation bias is the driving force of codon usage in the Gallus gallus genome

Synonymous codons are used with different frequencies both among species and among genes within the same genome and are controlled by neutral processes (such as mutation and drift) as well as by selection. Up to now, a systematic examination of the codon usage for the chicken genome has not been performed. Here, we carried out a whole genome analysis of the chicken genome by the use of the relative synonymous codon usage (RSCU) method and identified 11 putative optimal codons, all of them ending with uracil (U), which is significantly departing from the pattern observed in other eukaryotes. Optimal codons in the chicken genome are most likely the ones corresponding to highly expressed transfer RNA (tRNAs) or tRNA gene copy numbers in the cell. Codon bias, measured as the frequency of optimal codons (Fop), is negatively correlated with the G + C content, recombination rate, but positively correlated with gene expression, protein length, gene length and intron length. The positive correlation between codon bias and protein, gene and intron length is quite different from other multi-cellular organism, as this trend has been only found in unicellular organisms. Our data displayed that regional G + C content explains a large proportion of the variance of codon bias in chicken. Stepwise selection model analyses indicate that G + C content of coding sequence is the most important factor for codon bias. It appears that variation in the G + C content of CDSs accounts for over 60% of the variation of codon bias. This study suggests that both mutation bias and selection contribute to codon bias. However, mutation bias is the driving force of the codon usage in the Gallus gallus genome. Our data also provide evidence that the negative correlation between codon bias and recombination rates in G. gallus is determined mostly by recombination-dependent mutational patterns.

The dynamics of supply and demand in mRNA translation

We study the elongation stage of mRNA translation in eukaryotes and find that, in contrast to the assumptions of previous models, both the supply and the demand for tRNA resources are important for determining elongation rates. We find that increasing the initiation rate of translation can lead to the depletion of some species of aa-tRNA, which in turn can lead to slow codons and queueing. Particularly striking “competition” effects are observed in simulations of multiple species of mRNA which are reliant on the same pool of tRNA resources. These simulations are based on a recent model of elongation which we use to study the translation of mRNA sequences from the Saccharomyces cerevisiae genome. This model includes the dynamics of the use and recharging of amino acid tRNA complexes, and we show via Monte Carlo simulation that this has a dramatic effect on the protein production behaviour of the system.

In this paper we show that the rate at which proteins are produced can be controlled at the elongation stage of mRNA translation. Regulation of translation initiation has been a focus of much study, but the subsequent effect of changes in the initiation rate on the overall translation rate, and the role of slow and fast codon usage in mRNA sequences is still not fully understood. We consider a model of elongation in which the dynamics of tRNA use and recharging are considered for real mRNA sequences. We find that the balance between the demand for, and supply of tRNAs is crucial in determining translation rates. Particularly interesting “competition” effects are observed when the simultaneous translation of multiple mRNA is considered. We show indeed that, via the choice of slow or fast codons, it is in principle possible to control how variation of the supply and demand for tRNA resources changes the rate of protein production from different mRNAs.

Comparative analysis of ovine adenovirus 287 and human adenovirus 2 and 5 based on their codon usage

Ovine adenovirus 287 (OAdV287) emerges as one of the most promising gene vectors resulting from its unique biological characteristics. To obtain a more detailed knowledge about the codon usage of OAdV287, a comparative study based on the codon usage of OAdV287 and the prototypes of human adenovirus serotypes 2 and 5 (HAdV2/5) was carried out. Some commonly used indices measuring the codon usage patterns, including effective number of codons, relative synonymous codon usage, and statistical methods, were adopted. Overall, OAdV287 had a more biased and conservative codon usage pattern than that of HAdV2/5. Both mutation pressure and natural selection played important roles in shaping the codon usage patterns of these three adenoviruses. All the preference codons of OAdV287 had A/U ends and were totally different from those of sheep and humans; however, the preference codons of HAdV2/5 mostly had G/C ends and were mostly coincident with those of sheep and humans. The codon usage analysis in this study supplies some clues for further comprehending the unique biological characteristics of OAdV287 as gene vectors.