Codon Usage Influences the Local Rate of Translation Elongation to Regulate Co-translational Protein Folding

The Key Parameters that Govern Translation Efficiency

Translation of mRNA into protein is a fundamental yet complex biological process with multiple factors that can potentially affect its efficiency. Here, we study a stochastic model describing the traffic flow of ribosomes along the mRNA and identify the key parameters that govern the overall rate of protein synthesis, sensitivity to initiation rate changes, and efficiency of ribosome usage. By analyzing a continuum limit of the model, we obtain closed-form expressions for stationary currents and ribosomal densities, which agree well with Monte Carlo simulations. Furthermore, we completely characterize the phase transitions in the system, and by applying our theoretical results, we formulate design principles that detail how to tune the key parameters we identified to optimize translation efficiency. Using ribosome profiling data from S. cerevisiae, we show that its translation system is generally consistent with these principles. Our theoretical results have implications for evolutionary biology, as well as for synthetic biology.

The complete chloroplast genome sequence of yellow mustard (Sinapis alba L.) and its phylogenetic relationship to other Brassicaceae species

As a member of the large Brassicaceae family, yellow mustard (Sinapis alba L.) has been used as an important gene pool for the genetic improvement of cash crops in Brassicaceae. Understanding the phylogenetic relationship between Sinapis alba (S. alba) and other Brassicaceae crops can provide guidance on the introgression of its favorable alleles into related species. The chloroplast (cp) genome is an ideal model for assessing genome evolution and the phylogenetic relationships of complex angiosperm families. Herein, we de novo assembled the complete cp genome of S. alba by integrating the PacBio and Illumina sequencing platforms. A 153,760 bp quadripartite cycle without any gap was obtained, including a pair of inverted repeats (IRa and IRb) of 26,221 bp, separated by a large single copy (LSC) region of 83,506 bp and a small single copy (SSC) region of 17,821 bp. A total of 78 protein-coding genes, 30 tRNA genes, and four rRNA genes were identified in this cp genome, as were 89 simple sequence repeat (SSR) loci of 18 types. The codon usage analysis revealed a preferential use of the Leu codon with the A/U ending. The phylogenetic analysis using 82 Brassicaceae species demonstrated that S. alba had a close relationship with important Brassica and Raphanus species; moreover, it likely originated from a separate evolutionary pathway compared with the congeneric Sinapis arvensis. The synonymous (Ks) and non-synonymous (Ks) substitution rate analysis showed that genes encoding "Subunits of cytochrome b/f complex" were under the lowest purifying selection pressure, whereas those associated with "Maturase", "Subunit of acetyl-CoA", and "Subunits of NADH-dehydrogenase" underwent relatively higher purifying selection pressures. Our results provide valuable information for fully utilizing the S. alba cp genome as a potential genetic resource for the genetic improvement of Brassica and Raphanus species.

Elucidation of Codon Usage Signatures across the Domains of Life

Because of the degeneracy of the genetic code, multiple codons are translated into the same amino acid. Despite being “synonymous,” these codons are not equally used. Selective pressures are thought to drive the choice among synonymous codons within a genome, while GC content, which is typically attributed to mutational drift, is the major determinant of variation across species. Here, we find that in addition to GC content, interspecies codon usage signatures can also be detected. More specifically, we show that a single amino acid, arginine, is the major contributor to codon usage bias differences across domains of life. We then exploit this finding and show that domain-specific codon bias signatures can be used to classify a given sequence into its corresponding domain of life with high accuracy. We then wondered whether the inclusion of codon usage codon autocorrelation patterns, which reflects the nonrandom distribution of codon occurrences throughout a transcript, might improve the classification performance of our algorithm. However, we find that autocorrelation patterns are not domain-specific, and surprisingly, are unrelated to tRNA reusage, in contrast to previous reports. Instead, our results suggest that codon autocorrelation patterns are a by-product of codon optimality throughout a sequence, where highly expressed genes display autocorrelated “optimal” codons, whereas lowly expressed genes display autocorrelated “nonoptimal” codons.

Impacts of uORF codon identity and position on translation regulation

Translation regulation plays an important role in eukaryotic gene expression. Upstream open reading frames (uORFs) are potent regulatory elements located in 5′ mRNA transcript leaders. Translation of uORFs usually inhibit the translation of downstream main open reading frames, but some enhance expression. While a minority of uORFs encode conserved functional peptides, the coding regions of most uORFs are not conserved. Thus, the importance of uORF coding sequences on their regulatory functions remains largely unknown. We investigated the impact of an uORF coding region on gene regulation by assaying the functions of thousands of variants in the yeast YAP1 uORF. Varying uORF codons resulted in a wide range of functions, including repressing and enhancing expression of the downstream ORF. The presence of rare codons resulted in the most inhibitory YAP1 uORF variants. Inhibitory functions of such uORFs were abrogated by overexpression of complementary tRNA. Finally, regression analysis of our results indicated that both codon identity and position impact uORF function. Our results support a model in which a uORF coding sequence impacts its regulatory functions by altering the speed of uORF translation.

eRF1 mediates codon usage effects on mRNA translation efficiency through premature termination at rare codons

Codon usage bias is a universal feature of eukaryotic and prokaryotic genomes and plays an important role in regulating gene expression levels. A major role of codon usage is thought to regulate protein expression levels by affecting mRNA translation efficiency, but the underlying mechanism is unclear. By analyzing ribosome profiling results, here we showed that codon usage regulates translation elongation rate and that rare codons are decoded more slowly than common codons in all codon families in Neurospora. Rare codons resulted in ribosome stalling in manners both dependent and independent of protein sequence context and caused premature translation termination. This mechanism was shown to be conserved in Drosophila cells. In both Neurospora and Drosophila cells, codon usage plays an important role in regulating mRNA translation efficiency. We found that the rare codon-dependent premature termination is mediated by the translation termination factor eRF1, which recognizes ribosomes stalled on rare sense codons. Silencing of eRF1 expression resulted in codon usage-dependent changes in protein expression. Together, these results establish a mechanism for how codon usage regulates mRNA translation efficiency.

Coding regions affect mRNA stability in human cells

A new paradigm has emerged that coding regions can regulate mRNA stability in model organisms. Here, due to differences in cognate tRNA abundance, synonymous codons are translated at different speeds, and slow codons then stimulate mRNA decay. To ask if this phenomenon also occurs in humans, we isolated RNA stability effects due to coding regions using the human ORFeome collection. We find that many open reading frame (ORF) characteristics, such as length and secondary structure, fail to provide explanations for how coding regions alter mRNA stability, and, instead, that the ORF relies on translation to impact mRNA stability. Consistent with what has been seen in other organisms, codon use is related to the effects of ORFs on transcript stability. Importantly, we found instability-associated codons have longer A-site dwell times, suggesting for the first time in humans a connection between elongation speed and mRNA decay. Thus, we propose that codon usage alters decoding speeds and so affects human mRNA stability.

Comparative genome analysis reveals a distinct influence of nucleotide composition on virus-host species-specific interaction of prawn-infecting nodavirus

Discovery of species-specific interaction between the host and virus has drawn the interest of many researchers to study the evolution of the newly emerged virus. Comparative genome analysis provides insights of the virus functional genome evolution and the underlying mechanisms of virus-host interactions. The analysis of nucleotide composition signified the evolution of nodavirus towards host specialization in a host-specific mutation manner. GC-rich genome of betanodavirus was significantly deficient in UpA and UpU dinucleotides composition, whilst the AU-rich genome of gammanodavirus was deficient in CpG dinucleotide. The capsid of MrNV and PvNV of gammanodavirus retains the highest abundance of adenine and uracil at the second codon position, respectively, which were found to be very distinctive from the other genera. ENC-GC3 plot inferred the influence of natural selection and mutational pressure in shaping the evolution of MrNV RdRp and capsid, respectively. Furthermore, CAI/eCAI analysis predicts a comparable adaptability of MrNV in squid, Sepia officinalis than its natural host, Macrobrachium rosenbergii. Thus, further study is warranted to investigate the capacity of MrNV replication in S. officinalis owing to its high codon adaptation index.

A Single Synonymous Variant (c.354G>A [p.P118P]) in ADAMTS13 Confers Enhanced Specific Activity

Synonymous variants within coding regions may influence protein expression and function. We have previously reported increased protein expression levels ex vivo (~120% in comparison to wild-type) from a synonymous polymorphism variant, c.354G>A [p.P118P], of the ADAMTS13 gene, encoding a plasma protease responsible for von Willebrand Factor (VWF) degradation. In the current study, we investigated the potential mechanism(s) behind the increased protein expression levels from this variant and its effect on ADAMTS13 physico-chemical properties. Cell-free assays showed enhanced translation of the c.354G>A variant and the analysis of codon usage characteristics suggested that introduction of the frequently used codon/codon pair(s) may have been potentially responsible for this effect. Limited proteolysis, however, showed no substantial influence of altered translation on protein conformation. Analysis of post-translational modifications also showed no notable differences but identified three previously unreported glycosylation markers. Despite these similarities, p.P118P variant unexpectedly showed higher specific activity. Structural analysis using modeled interactions indicated that subtle conformational changes arising from altered translation kinetics could affect interactions between an exosite of ADAMTS13 and VWF resulting in altered specific activity. This report highlights how a single synonymous nucleotide variation can impact cellular expression and specific activity in the absence of measurable impact on protein structure.

[Synonymous Codon Usage-a Guide for Co-Translational Protein Folding in the Cell]

In the cell, protein folding begins during protein synthesis/translation and thus is a co-translational process. Co-translational protein folding is tightly linked to translation elongation, which is not a uniform process. While there are many reasons for translation non-uniformity, it is generally believed that non-uniform synonymous codon usage is one of the key factors modulating translation elongation rates. Frequent/optimal codons as a rule are translated more rapidly than infrequently used ones and vice versa. Over 30 years ago, it was hypothesized that changes in synonymous codon usage affecting translation elongation rates could impinge on co-translation protein folding and that many synonymous codons are strategically placed within mRNA to ensure a particular translation kinetics facilitating productive step-by-step co-translational folding of proteins. It was suggested that this particular translation kinetics (and, specifically, translation pause sites) may define the window of opportunity for the protein parts to fold locally, particularly at the critical points where folding is far from equilibrium. It was thus hypothesized that synonymous codons may provide a secondary code for protein folding in the cell. Although, mostly accepted now, this hypothesis appeared to be difficult to prove and many convincing results were obtained only relatively recently. Here, I review the progress in the field and explain, why this simple idea appeared to be so challenging to prove.

A pair of non-optimal codons are necessary for the correct biosynthesis of the Aspergillus nidulans urea transporter, UreA

In both prokaryotic and eukaryotic genomes, synonymous codons are unevenly used. Such differential usage of optimal or non-optimal codons has been suggested to play a role in the control of translation initiation and elongation, as well as at the level of transcription and mRNA stability. In the case of membrane proteins, codon usage has been proposed to assist in the establishment of a pause necessary for the correct targeting of the nascent chains to the translocon. By using as a model UreA, the Aspergillus nidulans urea transporter, we revealed that a pair of non-optimal codons encoding amino acids situated at the boundary between the N-terminus and the first transmembrane segment are necessary for proper biogenesis of the protein at 37°C. These codons presumably regulate the translation rate in a previously undescribed fashion, possibly contributing to the correct interaction of ureA-translating ribosome-nascent chain complexes with the signal recognition particle and/or other factors, while the polypeptide has not yet emerged from the ribosomal tunnel. Our results suggest that the presence of the pair of non-optimal codons would not be functionally important in all cellular conditions. Whether this mechanism would affect other proteins remains to be determined.

Effects of codon optimization on coagulation factor IX translation and structure: Implications for protein and gene therapies

Synonymous codons occur with different frequencies in different organisms, a phenomenon termed codon usage bias. Codon optimization, a common term for a variety of approaches used widely by the biopharmaceutical industry, involves synonymous substitutions to increase protein expression. It had long been presumed that synonymous variants, which, by definition, do not alter the primary amino acid sequence, have no effect on protein structure and function. However, a critical mass of reports suggests that synonymous codon variations may impact protein conformation. To investigate the impact of synonymous codons usage on protein expression and function, we designed an optimized coagulation factor IX (FIX) variant and used multiple methods to compare its properties to the wild-type FIX upon expression in HEK293T cells. We found that the two variants differ in their conformation, even when controlling for the difference in expression levels. Using ribosome profiling, we identified robust changes in the translational kinetics of the two variants and were able to identify a region in the gene that may have a role in altering the conformation of the protein. Our data have direct implications for codon optimization strategies, for production of recombinant proteins and gene therapies.

Parsing the synonymous mutations in the maize genome: isoaccepting mutations are more advantageous in regions with codon co-occurrence bias

BACKGROUND

Synonymous mutations do not change amino acids but do sometimes change the tRNAs (anticodons) that decode a particular codon. An isoaccepting codon is a synonymous codon that shares the same tRNA. If a mutated codon could base pair with the same anticodon as the original, the mutation is termed an isoaccepting mutation. An interesting but less-studied type of codon bias is codon co-occurrence bias. There is a trend to cluster the isoaccepting codons in the genome. The proposed advantage of codon co-occurrence bias is that the tRNA released from the ribosome E site could be quickly recharged and subsequently decode the following isoaccepting codons. This advantage would enhance translation efficiency. In plant species, whether there are signals of positive selection on isoaccepting mutations in the codon co-occurred regions has not been studied.

RESULTS

We termed polymorphic mutations in coding regions using publicly available RNA-seq data in maize (Zea mays). Next, we classified all synonymous mutations into three categories according to the context, i.e., the relationship between the focal codon and the previous codon, as follows: isoaccepting, nonisoaccepting and nonsynonymous. We observed higher fractions of isoaccepting mutations in the isoaccepting context. If we looked at the minor allele frequency (MAF) spectrum, the isoaccepting mutations have a higher MAF in the isoaccepting context than that in other regions, and accordingly, the nonisoaccepting mutations have a higher MAF in the nonisoaccepting context.

CONCLUSION

Our results indicate that in regions with codon co-occurrence bias, natural selection maintains this pattern by suppressing the nonisoaccepting mutations. However, if the consecutive codons are nonisoaccepting, mutations tend to switch these codons to become isoaccepting. Our study demonstrates that the codon co-occurrence bias in the maize genome is selectively maintained by natural selection and that the advantage of this trend could potentially be the rapid recharging and reuse of tRNAs to increase translation efficiency.

Alternative Splicing in Heat Shock Protein Transcripts as a Mechanism of Cell Adaptation in Trichophyton rubrum

Heat shock proteins (HSPs) are involved in critical processes like host tissue invasion, resistance, and pathogenicity in dermatophytes. RNA-Seq analysis of Trichophyton rubrum exposed to undecanoic acid (UDA) revealed intron retention events in HSP transcripts. Because HSPs are modulated in response to various stimuli and as alternative splicing (AS) can result in a broad diversity in the proteome of eukaryotic cells, our objective was to confirm the aforementioned retention events, investigating their consequences and extent. Furthermore, we aimed to determine: (1) the expression profile of HSP genes in an infection-like scenario and (2) the importance of Hsp90 for the keratinolytic potential of T. rubrum. RT and qPCR analyses comparing the exposure to UDA and terbinafine (TRB) confirmed the presence of two mRNA isoforms of the hsp7-like gene, with distinct expression patterns in response to UDA and TRB. The HSP expression profile revealed two upregulated, three downregulated, and four unmodulated transcripts; Hsp90 inhibition by 17-AAG resulted in a significant decrease in keratinolytic potential at 37 °C. Altogether, these results broaden the current knowledge on the importance of HSP-mediated pathways for cell adaptation and other aspects of dermatophyte biology, indicating that HSP network proteins can be potential targets for antifungal therapy.

Compositional features and codon usage pattern of TP63 gene

The tumor protein p63encoded by the gene TP63 acts as a homologue of p53 protein. TP63 gene is the transformation factor with two initiation sites for transcriptional process and is related with stress, signal transduction and cell cycle control. The biasness in the preference of a few codons more frequently over other synonymous codons is the codon usage bias (CUB). Natural selection and mutational pressure are the two prime evolutionary forces acting on CUB. Here, the bioinformatic based analysis was performed to investigate the base distribution and CUB of TP63transcript variants (isoforms) as no work was performed earlier. Analysis of compositional features revealed variation in base content across TP63 gene isoforms and the GC content was more than 50%, indicating GC richness of its isoforms. The mean effective number of codons (ENC), a measure of CUB, was 51.83, i.e. overall CUB of TP63 gene was low. Among 13 isoforms of TP63 gene, nature selected against the CTA codon in 8 isoforms and favored five over-represented (RSCU > 1.6) codons namely CTG, CAG, ATC, AAC and GCC during evolution. Correlation between overall nucleotide composition and its 3rd codon position revealed that both mutational pressure and natural selection moulded its CUB. Further, the correlation between ENC and aromaticity depicted that variation of CUB was related to the degree of aromaticity of p63 protein.

Codon bias confers stability to human mRNAs

Codon bias has been implicated as one of the major factors contributing to mRNA stability in several model organisms. However, the molecular mechanisms of codon bias on mRNA stability remain unclear in humans. Here, we show that human cells possess a mechanism to modulate RNA stability through a unique codon bias. Bioinformatics analysis showed that codons could be clustered into two distinct groups-codons with G or C at the third base position (GC3) and codons with either A or T at the third base position (AT3): the former stabilizing while the latter destabilizing mRNA. Quantification of codon bias showed that increased GC3-content entails proportionately higher GC-content. Through bioinformatics, ribosome profiling, and in vitro analysis, we show that decoupling the effects of codon bias reveals two modes of mRNA regulation, one GC3- and one GC-content dependent. Employing an immunoprecipitation-based strategy, we identify ILF2 and ILF3 as RNA-binding proteins that differentially regulate global mRNA abundances based on codon bias. Our results demonstrate that codon bias is a two-pronged system that governs mRNA abundance.

Characterizing the heat response of Arabidopsis thaliana from the perspective of codon usage bias and translational regulation

mRNA translation is carefully regulated at both the initiation and the elongation step. Under heat stress, it is known that particular genes change their expression and translation levels to respond to the environment. Attention has been paid to the detailed mechanisms of how a few proteins work, and little is done to analyze whether the global evolutionary patterns affect the translational changes. Determinants like codon usage bias and its related evolutionary features are less studied in heat stress experiments, especially for plants. Utilizing the RNA-seq and Ribo-seq data of normal and heat-stressed Arabidopsis thaliana generated from a previous study, we conducted gene-level (global) and codon-resolution (local) translation analyses. We studied how codon usage bias and other evolutionary features could impact the translation patterns in the heat response of the plant. We found that the evolutionary features including codon usage bias, tAI, nitrogen cost, and conservation (identity) could affect the global and local translation efficiency. Under heat stress, the optimal and conserved codons are more likely to alter their local translation elongation speed to modulate the global translation of host genes. Meanwhile, we also verified the widely accepted notions that the secondary structures and proline codons could largely slow down the translation rate. Our results revealed the effect of codon usage bias and other evolutionary patterns on the translation regulation under heat stress. Unveiling the effect of these features on translational regulation of plants might be helpful in understanding the relationship and interaction between plants and the environment.

GWIPS-viz: 2018 update

The GWIPS-viz browser (http://gwips.ucc.ie/) is an on-line genome browser which is tailored for exploring ribosome profiling (Ribo-seq) data. Since its publication in 2014, GWIPS-viz provides Ribo-seq data for an additional 14 genomes bringing the current total to 23. The integration of new Ribo-seq data has been automated thereby increasing the number of available tracks to 1792, a 10-fold increase in the last three years. The increase is particularly substantial for data derived from human sources. Following user requests, we added the functionality to download these tracks in bigWig format. We also incorporated new types of data (e.g. TCP-seq) as well as auxiliary tracks from other sources that help with the interpretation of Ribo-seq data. Improvements in the visualization of the data have been carried out particularly for bacterial genomes where the Ribo-seq data are now shown in a strand specific manner. For higher eukaryotic datasets, we provide characteristics of individual datasets using the RUST program which includes the triplet periodicity, sequencing biases and relative inferred A-site dwell times. This information can be used for assessing the quality of Ribo-seq datasets. To improve the power of the signal, we aggregate Ribo-seq data from several studies into Global aggregate tracks for each genome.

Ribosome profiling of selenoproteins in vivo reveals consequences of pathogenic Secisbp2 missense mutations

Recoding of UGA codons as selenocysteine (Sec) codons in selenoproteins depends on a selenocysteine insertion sequence (SECIS) in the 3'-UTR of mRNAs of eukaryotic selenoproteins. SECIS-binding protein 2 (SECISBP2) increases the efficiency of this process. Pathogenic mutations in SECISBP2 reduce selenoprotein expression and lead to phenotypes associated with the reduction of deiodinase activities and selenoprotein N expression in humans. Two functions have been ascribed to SECISBP2: binding of SECIS elements in selenoprotein mRNAs and facilitation of co-translational Sec insertion. To separately probe both functions, we established here two mouse models carrying two pathogenic missense mutations in Secisbp2 previously identified in patients. We found that the C696R substitution in the RNA-binding domain abrogates SECIS binding and does not support selenoprotein translation above the level of a complete Secisbp2 null mutation. The R543Q missense substitution located in the selenocysteine insertion domain resulted in residual activity and caused reduced selenoprotein translation, as demonstrated by ribosomal profiling to determine the impact on UGA recoding in individual selenoproteins. We found, however, that the R543Q variant is thermally unstable in vitro and completely degraded in the mouse liver in vivo, while being partially functional in the brain. The moderate impairment of selenoprotein expression in neurons led to astrogliosis and transcriptional induction of genes associated with immune responses. We conclude that differential SECISBP2 protein stability in individual cell types may dictate clinical phenotypes to a much greater extent than molecular interactions involving a mutated amino acid in SECISBP2.

Genomic Evidence for Simultaneous Optimization of Transcription and Translation through Codon Variants in the pmoCAB Operon of Type Ia Methanotrophs

Microbial methane oxidation plays a fundamental role in the biogeochemical cycle of Earth’s system. Recent reports have provided evidence for the acquisition of methane monooxygenases by horizontal gene transfer in methane-oxidizing bacteria from different environments, but how evolution has shaped the coding sequences to execute methanotrophy efficiently remains unexplored. In this work, we provide genomic evidence that among the different types of methanotrophs, type Ia methanotrophs possess a unique coding sequence of the pmoCAB operon that is under positive selection for optimal resource allocation and efficient synthesis of transcripts and proteins. This adaptive trait possibly enables type Ia methanotrophs to respond robustly to fluctuating methane availability and explains their global prevalence.

Understanding the interplay between genotype and phenotype is a fundamental goal of functional genomics. Methane oxidation is a microbial phenotype with global-scale significance as part of the carbon biogeochemical cycle and a sink for greenhouse gas. Microorganisms that oxidize methane (methanotrophs) are taxonomically diverse and widespread around the globe. In methanotrophic bacteria, enzymes in the methane oxidation metabolic module (KEGG module M00174, conversion of methane to formaldehyde) are encoded in four operons (pmoCAB, mmoXYZBCD, mxaFI, and xoxF). Recent reports have suggested that methanotrophs in Proteobacteria acquired methane monooxygenases through horizontal gene transfer. Here, we used a genomic meta-analysis to infer the transcriptional and translational advantages of coding sequences from the methane oxidation metabolic modules of different types of methanotrophs. By analyzing isolate and metagenome-assembled genomes from phylogenetically and geographically diverse sources, we detected an anomalous nucleotide composition bias in the coding sequences of particulate methane monooxygenase genes (pmoCAB) from type Ia methanotrophs. We found that this nucleotide bias increases the level of codon bias by decreasing the GC content in the third base of codons, a strategy that contrasts with that of other coding sequences in the module. Further codon usage analyses uncovered that codon variants of the type Ia pmoCAB coding sequences deviate from the genomic signature to match ribosomal protein-coding sequences. Subsequently, computation of transcription and translation metrics revealed that the pmoCAB coding sequences of type Ia methanotrophs optimize the usage of codon variants to maximize translation efficiency and accuracy, while minimizing the synthesis cost of transcripts and proteins.

IMPORTANCE Microbial methane oxidation plays a fundamental role in the biogeochemical cycle of Earth’s system. Recent reports have provided evidence for the acquisition of methane monooxygenases by horizontal gene transfer in methane-oxidizing bacteria from different environments, but how evolution has shaped the coding sequences to execute methanotrophy efficiently remains unexplored. In this work, we provide genomic evidence that among the different types of methanotrophs, type Ia methanotrophs possess a unique coding sequence of the pmoCAB operon that is under positive selection for optimal resource allocation and efficient synthesis of transcripts and proteins. This adaptive trait possibly enables type Ia methanotrophs to respond robustly to fluctuating methane availability and explains their global prevalence.

Oxidative stress damages rRNA inside the ribosome and differentially affects the catalytic center

Intracellular levels of reactive oxygen species (ROS) increase as a consequence of oxidative stress and represent a major source of damage to biomolecules. Due to its high cellular abundance RNA is more frequently the target for oxidative damage than DNA. Nevertheless the functional consequences of damage on stable RNA are poorly understood. Using a genome-wide approach, based on 8-oxo-guanosine immunoprecipitation, we present evidence that the most abundant non-coding RNA in a cell, the ribosomal RNA (rRNA), is target for oxidative nucleobase damage by ROS. Subjecting ribosomes to oxidative stress, we demonstrate that oxidized 23S rRNA inhibits the ribosome during protein biosynthesis. Placing single oxidized nucleobases at specific position within the ribosome's catalytic center by atomic mutagenesis resulted in markedly different functional outcomes. While some active site nucleobases tolerated oxidative damage well, oxidation at others had detrimental effects on protein synthesis by inhibiting different sub-steps of the ribosomal elongation cycle. Our data provide molecular insight into the biological consequences of RNA oxidation in one of the most central cellular enzymes and reveal mechanistic insight on the role of individual active site nucleobases during translation.

Nature and Regulation of Protein Folding on the Ribosome

Co-translational protein folding is an essential process by which cells ensure the safe and efficient production and assembly of new proteins in their functional native states following biosynthesis on the ribosome. In this review, we describe recent progress in probing the changes during protein synthesis of the free energy landscapes that underlie co-translational folding and discuss the critical coupling between these landscapes and the rate of translation that ultimately determines the success or otherwise of the folding process. Recent developments have revealed a variety of mechanisms by which both folding and translation can be modulated or regulated, and we discuss how these effects are utilised by the cell to optimise the outcome of protein biosynthesis.

Translation can affect the antisense activity of RNase H1-dependent oligonucleotides targeting mRNAs

RNase H1-dependent antisense oligonucleotides (ASOs) can degrade complementary RNAs in both the nucleus and the cytoplasm. Since cytoplasmic mRNAs are actively engaged in translation, ASO activity may thus be affected by translating ribosomes that scan the mRNAs. Here we show that mRNAs associated with ribosomes can be cleaved using ASOs and that translation can alter ASO activity. Translation inhibition tends to increase ASO activity when targeting the coding regions of efficiently translated mRNAs, but not nuclear non-coding RNAs or less efficiently translated mRNAs. Increasing the level of RNase H1 protein eliminated the enhancing effects of translation inhibition on ASO activity, suggesting that RNase H1 recruitment to ASO/mRNA heteroduplexes is a rate limiting step and that translating ribosomes can inhibit RNase H1 recruitment. Consistently, ASO activity was not increased by translation inhibition when targeting the 3′ UTRs, independent of the translation efficiency of the mRNAs. Contrarily, the activity of 3′ UTR-targeting ASOs tended to be reduced upon translation inhibition, likely due to decreased accessibility. These results indicate that ASO activity can be affected by the translation process, and the findings also provide important information toward helping better ASO drug design.

Protein Translation and Psychiatric Disorders

Although historically research has focused on transcription as the central governor of protein expression, protein translation is now increasingly being recognized as a major factor for determining protein levels within cells. The central nervous system relies on efficient updating of the protein landscape. Thus, coordinated regulation of mRNA localization, initiation, or termination of translation is essential for proper brain function. In particular, dendritic protein synthesis plays a key role in synaptic plasticity underlying learning and memory as well as cognitive processes. Increasing evidence suggests that impaired mRNA translation is a common feature found in numerous psychiatric disorders. In this review, we describe how malfunction of translation contributes to development of psychiatric diseases, including schizophrenia, major depression, bipolar disorder, and addiction.

Variation and selection on codon usage bias across an entire subphylum

Variation in synonymous codon usage is abundant across multiple levels of organization: between codons of an amino acid, between genes in a genome, and between genomes of different species. It is now well understood that variation in synonymous codon usage is influenced by mutational bias coupled with both natural selection for translational efficiency and genetic drift, but how these processes shape patterns of codon usage bias across entire lineages remains unexplored. To address this question, we used a rich genomic data set of 327 species that covers nearly one third of the known biodiversity of the budding yeast subphylum Saccharomycotina. We found that, while genome-wide relative synonymous codon usage (RSCU) for all codons was highly correlated with the GC content of the third codon position (GC3), the usage of codons for the amino acids proline, arginine, and glycine was inconsistent with the neutral expectation where mutational bias coupled with genetic drift drive codon usage. Examination between genes' effective numbers of codons and their GC3 contents in individual genomes revealed that nearly a quarter of genes (381,174/1,683,203; 23%), as well as most genomes (308/327; 94%), significantly deviate from the neutral expectation. Finally, by evaluating the imprint of translational selection on codon usage, measured as the degree to which genes' adaptiveness to the tRNA pool were correlated with selective pressure, we show that translational selection is widespread in budding yeast genomes (264/327; 81%). These results suggest that the contribution of translational selection and drift to patterns of synonymous codon usage across budding yeasts varies across codons, genes, and genomes; whereas drift is the primary driver of global codon usage across the subphylum, the codon bias of large numbers of genes in the majority of genomes is influenced by translational selection.

Posttranscriptional Gene Regulation of the GABA Receptor to Control Neuronal Inhibition

Behavior and higher cognition rely on the transfer of information between neurons through specialized contact sites termed synapses. Plasticity of neuronal circuits, a prerequisite to respond to environmental changes, is intrinsically coupled with the nerve cell’s ability to form, structurally modulate or remove synapses. Consequently, the synaptic proteome undergoes dynamic alteration on demand in a spatiotemporally restricted manner. Therefore, proper protein localization at synapses is essential for synaptic function. This process is regulated by: (i) protein transport and recruitment; (ii) local protein synthesis; and (iii) synaptic protein degradation. These processes shape the transmission efficiency of excitatory synapses. Whether and how these processes influence synaptic inhibition is, however, widely unknown. Here, we summarize findings on fundamental regulatory processes that can be extrapolated to inhibitory synapses. In particular, we focus on known aspects of posttranscriptional regulation and protein dynamics of the GABA receptor (GABAR). Finally, we propose that local (co)-translational control mechanism might control transmission of inhibitory synapses.

Enterovirus A71 Containing Codon-Deoptimized VP1 and High-Fidelity Polymerase as Next-Generation Vaccine Candidate

Enterovirus A71 (EV-A71) is a major pathogen that causes hand-foot-and-mouth disease (HFMD), which occasionally results in severe neurological complications. In this study, we developed four EV-A71 (rgEV-A71) strains by reverse genetics procedures as possible vaccine candidates. The four rgEV-A71 viruses contained various codon-deoptimized VP1 capsid proteins (VP1-CD) and showed replication rates and antigenicity similar to that of the wild-type virus, while a fifth virus, rg4643C4VP-CD, was unable to form plaques but was still able to be examined by median tissue culture infectious dose (TCID₅₀) titers, which were similar to those of the others, indicating the effect of CD on plaque formation. However, the genome stability showed that there were some mutations which appeared during just one passage of the VP1-CD viruses. Thus, we further constructed VP1-CD rgEV-A71 containing high-fidelity determinants in 3D polymerase (CD-HF), and the number of mutations in CD-HF rgEV-A71 was shown to have decreased. The CD-HF viruses showed less virulence than the parental strain in a mouse infection model. After 14 days postimmunization, antibody titers had increased in mice infected with CD-HF viruses. The mouse antisera showed similar neutralizing antibody titers against various CD-HF viruses and different genotypes of EV-A71. The study demonstrates the proof of concept that VP1 codon deoptimization combined with high-fidelity 3D polymerase decreased EV-A71 mutations and virulence in mice but retained their antigenicity, indicating it is a good candidate for next-generation EV-A71 vaccine development.IMPORTANCE EV-A71 can cause severe neurological diseases with fatality in infants and young children, but there are still no effective drugs to date. Here, we developed a novel vaccine strategy with the combination of CD and HF substitutions to generate the genetically stable reverse genetics virus. We found that CD combined with HF polymerase decreased the virulence but maintained the antigenicity of the virus. This work demonstrated the simultaneous introduction of CD genome sequences and HF substitutions as a potential new strategy to develop attenuated vaccine seed virus. Our work provides insight into the development of a low-virulence candidate vaccine virus through a series of genetic editing of virus sequences while maintaining its antigenicity and genome stability, which will provide an additional strategy for next-generation vaccine development of EV-A71.

A pan-cancer analysis of synonymous mutations

Synonymous mutations have been viewed as silent mutations, since they only affect the DNA and mRNA, but not the amino acid sequence of the resulting protein. Nonetheless, recent studies suggest their significant impact on splicing, RNA stability, RNA folding, translation or co-translational protein folding. Hence, we compile 659194 synonymous mutations found in human cancer and characterize their properties. We provide the user-friendly, comprehensive resource for synonymous mutations in cancer, SynMICdb (http://SynMICdb.dkfz.de), which also contains orthogonal information about gene annotation, recurrence, mutation loads, cancer association, conservation, alternative events, impact on mRNA structure and a SynMICdb score. Notably, synonymous and missense mutations are depleted at the 5'-end of the coding sequence as well as at the ends of internal exons independent of mutational signatures. For patient-derived synonymous mutations in the oncogene KRAS, we indicate that single point mutations can have a relevant impact on expression as well as on mRNA secondary structure.

Synonymous mutations do not alter amino acid sequence but may exert oncogenic effects in other ways. Here, the authors present a catalogue of synonymous mutations in cancer and characterise their properties.

Codon Clarity or Conundrum?

Synonymous variations in protein-coding sequences alter protein expression dynamics, which has important implications for cellular physiology and evolutionary fitness, but disentangling the underlying molecular mechanisms remains challenging.

Leveraging synthetic biology for producing bioactive polyketides and non-ribosomal peptides in bacterial heterologous hosts

Bacteria have historically been a rich source of natural products (e.g. polyketides and non-ribosomal peptides) that possess medically-relevant activities. Despite extensive discovery programs in both industry and academia, a plethora of biosynthetic pathways remain uncharacterized and the corresponding molecular products untested for potential bioactivities. This knowledge gap comes in part from the fact that many putative natural product producers have not been cultured in conventional laboratory settings in which the corresponding products are produced at detectable levels. Next-generation sequencing technologies are further increasing the knowledge gap by obtaining metagenomic sequence information from complex communities where production of the desired compound cannot be isolated in the laboratory. For these reasons, many groups are turning to synthetic biology to produce putative natural products in heterologous hosts. This strategy depends on the ability to heterologously express putative biosynthetic gene clusters and produce relevant quantities of the corresponding products. Actinobacteria remain the most abundant source of natural products and the most promising heterologous hosts for natural product discovery and production. However, researchers are discovering more natural products from other groups of bacteria, such as myxobacteria and cyanobacteria. Therefore, phylogenetically similar heterologous hosts have become promising candidates for synthesizing these novel molecules. The downside of working with these microbes is the lack of well-characterized genetic tools for optimizing expression of gene clusters and product titers. This review examines heterologous expression of natural product gene clusters in terms of the motivations for this research, the traits desired in an ideal host, tools available to the field, and a survey of recent progress.

Translation affects mRNA stability in a codon-dependent manner in human cells

mRNA translation decodes nucleotide into amino acid sequences. However, translation has also been shown to affect mRNA stability depending on codon composition in model organisms, although universality of this mechanism remains unclear. Here, using three independent approaches to measure exogenous and endogenous mRNA decay, we define which codons are associated with stable or unstable mRNAs in human cells. We demonstrate that the regulatory information affecting mRNA stability is encoded in codons and not in nucleotides. Stabilizing codons tend to be associated with higher tRNA levels and higher charged/total tRNA ratios. While mRNAs enriched in destabilizing codons tend to possess shorter poly(A)-tails, the poly(A)-tail is not required for the codon-mediated mRNA stability. This mechanism depends on translation; however, the number of ribosome loads into a mRNA modulates the codon-mediated effects on gene expression. This work provides definitive evidence that translation strongly affects mRNA stability in a codon-dependent manner in human cells.

Proteins are made by joining together building blocks called amino acids into strings. The proteins are ‘translated’ from genetic sequences called mRNA molecules. These sequences can be thought of as series of ‘letters’, which are read in groups of three known as codons. Molecules called tRNAs recognize the codons and add the matching amino acids to the end of the protein. Each tRNA can recognize one or several codons, and the levels of different tRNAs inside the cell vary.

There are 61 codons that code for amino acids, but only 20 amino acids. This means that some codons produce the same amino acid. Despite this, there is evidence to suggest that not all of the codons that produce the same amino acid are exactly equivalent. In bacteria, yeast and zebrafish, some codons seem to make the mRNA molecule more stable, and others make it less stable. This might help the cell to control how many proteins it makes. It was not clear whether the same is true for humans.

To find out, Wu et al. used three separate methods to examine mRNA stability in four types of human cell. Overall, the results revealed that some codons help to stabilize the mRNA, while others make the mRNA molecule break down faster. The effect seems to depend on the supply of tRNAs that have a charged amino acid; mRNA molecules were more likely to self-destruct in cells that contained codons with low levels of the tRNA molecules.

Wu et al. also found that conditions in the cell can alter how strongly the codons affect mRNA stability. For example, a cell that has been infected by a virus reduces translation. Under these conditions, the identity of the codons in the mRNA has less effect on the stability of the mRNA molecule.

Changes to protein production happen in many diseases. Understanding what controls these changes could help to reveal more about our fundamental biology, and what happens when it goes wrong.

Nonsynonymous, synonymous and nonsense mutations in human cancer-related genes undergo stronger purifying selections than expectation

Background

Nonsynonymous mutations change the protein sequences and are frequently subjected to natural selection. The same goes for nonsense mutations that introduce pre-mature stop codons into CDSs (coding sequences). Synonymous mutations, however, are intuitively thought to be functionally silent and evolutionarily neutral. Now researchers know that the optimized synonymous codon usage is advantageous in the speedy mRNA translation process. With the advent of NGS technique, the explosion of NGS data generated from the tumor tissues help researchers identify driver mutations in cancer-related genes, but relatively less attention is paid to the SNP data in healthy human populations when studying cancer.

Methods

Here, we analyzed the publically available human SNPs. We classified these SNPs according to their functional and evolutionary categories. By simply dividing the human genes into cancer-related genes and other genes, we compared the features of nonsynonymous, synonymous and nonsense mutations in these two gene sets from multiple aspects.

Results

We provided lines of evidence that the nonsynonymous, synonymous and nonsense mutations in cancer-related genes undergo stronger purifying selection when compared to the expected pattern in other genes. The lower nonsynonymous to synonymous ratio observed in cancer-related genes suggests the suppression of amino acid substitutions in these genes. The synonymous SNPs, after excluding those in splicing regions, exhibit preferred changes in codon usage and higher codon frequencies in cancer-related genes compared to other genes, indicating the constraint exerted on these mutations. Nonsense mutations are less frequent and located closer to stop codons in cancer-related genes than in other genes, which putatively minimize their deleterious effects.

Conclusion

Our study demonstrated the evolutionary constraint on mutations in CDS of cancer-related genes without the requirement of data from cancer tissues or patients. Our work provides novel perspectives on interpreting the constraint on mutations in cancer-related genes. We reveal extra constraint on synonymous mutations in cancer-related genes which is related to codon usage bias and is in addition to the splicing effect.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5572-x) contains supplementary material, which is available to authorized users.

A User's Guide to Cell-Free Protein Synthesis

Cell-free protein synthesis (CFPS) is a platform technology that provides new opportunities for protein expression, metabolic engineering, therapeutic development, education, and more. The advantages of CFPS over in vivo protein expression include its open system, the elimination of reliance on living cells, and the ability to focus all system energy on production of the protein of interest. Over the last 60 years, the CFPS platform has grown and diversified greatly, and it continues to evolve today. Both new applications and new types of extracts based on a variety of organisms are current areas of development. However, new users interested in CFPS may find it challenging to implement a cell-free platform in their laboratory due to the technical and functional considerations involved in choosing and executing a platform that best suits their needs. Here we hope to reduce this barrier to implementing CFPS by clarifying the similarities and differences amongst cell-free platforms, highlighting the various applications that have been accomplished in each of them, and detailing the main methodological and instrumental requirement for their preparation. Additionally, this review will help to contextualize the landscape of work that has been done using CFPS and showcase the diversity of applications that it enables.

Enhancing co-translational folding of heterologous protein by deleting non-essential ribosomal proteins in Pichia pastoris

BACKGROUND

Translational regulation played an important role in the correct folding of heterologous proteins to form bioactive conformations during biogenesis. Translational pausing coordinates protein translation and co-translational folding. Decelerating translation elongation speed has been shown to improve the soluble protein yield when expressing heterologous proteins in industrial expression hosts. However, rational redesign of translational pausing via synonymous mutations may not be feasible in many cases. Our goal was to develop a general and convenient strategy to improve heterologous protein synthesis in Pichia pastoris without mutating the expressed genes.

RESULTS

Here, a large-scale deletion library of ribosomal protein (RP) genes was constructed for heterologous protein expression in Pichia pastoris, and 59% (16/27) RP deletants have significantly increased heterologous protein yield. This is due to the delay of 60S subunit assembly by deleting non-essential ribosomal protein genes or 60S subunit processing factors, thus globally decreased the translation elongation speed and improved the co-translational folding, without perturbing the relative transcription level and translation initiation.

CONCLUSION

Global decrease in the translation elongation speed by RP deletion enhanced co-translational folding efficiency of nascent chains and decreased protein aggregates to improve heterologous protein yield. A potential expression platform for efficient pharmaceutical proteins and industrial enzymes production was provided without synonymous mutation.

Multiplexing polysome profiling experiments to study translation in Escherichia coli

Polysome profiling is a widely used method to monitor the translation status of mRNAs. Although it is theoretically a simple technique, it is labor intensive. Repetitive polysome fractionation rapidly generates a large number of samples to be handled in the downstream processes of protein elimination, RNA extraction and quantification. Here, we propose a multiplex polysome profiling experiment in which distinct cellular extracts are pooled before loading on the sucrose gradient for fractionation. We used the multiplexing method to study translation in E. coli. Multiplexing polysome profiling experiments provided similar mRNA translation status to that obtained with the non-multiplex method with comparable distribution of mRNA copies between the polysome profiling fractions, similar ribosome occupancy and ribosome density. The multiplexing method was used for parallel characterization of gene translational responses to changing mRNA levels. When the mRNA level of two native genes, cysZ and lacZ was increased by transcription induction, their global translational response was similar, with a higher ribosome load leading to increased ribosome occupancy and ribosome densities. However the pattern and the magnitude of the translational response were gene specific. By reducing the number of polysome profiling experiments, the multiplexing method saved time and effort and reduced cost and technical bias. This method would be useful to study the translational effect of mRNA sequence-dependent parameters that often require testing multiple samples and conditions in parallel.

In vivo characterisation of fluorescent proteins in budding yeast

Fluorescent proteins (FPs) are widely used in many organisms, but are commonly characterised in vitro. However, the in vitro properties may poorly reflect in vivo performance. Therefore, we characterised 27 FPs in vivo using Saccharomyces cerevisiae as model organism. We linked the FPs via a T2A peptide to a control FP, producing equimolar expression of the 2 FPs from 1 plasmid. Using this strategy, we characterised the FPs for brightness, photostability, photochromicity and pH-sensitivity, achieving a comprehensive in vivo characterisation. Many FPs showed different in vivo properties compared to existing in vitro data. Additionally, various FPs were photochromic, which affects readouts due to complex bleaching kinetics. Finally, we codon optimized the best performing FPs for optimal expression in yeast, and found that codon-optimization alters FP characteristics. These FPs improve experimental signal readout, opening new experimental possibilities. Our results may guide future studies in yeast that employ fluorescent proteins.

Quantification and discovery of sequence determinants of protein-per-mRNA amount in 29 human tissues

Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein-to-mRNA (PTR) ratios and a quantification of their effects are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post-transcriptional regulatory elements.

Stable G-quadruplex enabling sequences are selected against by the context-dependent codon bias

The distributions of secondary structural elements appear to differ between coding regions (CDS) of mRNAs compared to the untranslated regions (UTRs), presumably as a mechanism to fine-tune gene expression, including efficiency of translation. However, a systematic and comprehensive analysis of secondary structure avoidance because of potential bias in codon usage is difficult as some of the common secondary structures, such as, hairpins can be formed by numerous sequence combinations. Using G-quadruplex (GQ) as the model secondary structure we studied the impact of codon bias on GQs within the CDS. Because GQs can be predicted using specific consensus sequence motifs, they provide an excellent platform for investigation of the selectivity of such putative structures at the codon level. Using a bioinformatics approach, we calculated the frequencies of putative GQs within the CDS of a variety of species. Our results suggest that the most stable GQs appear to be significantly underrepresented within the CDS, through the use of specific synonymous codon combinations. Furthermore, we identified many peptide sequence motifs in which silent mutations can potentially alter translation via stable GQ formation. This work not only provides a comprehensive analysis on how stable secondary structures appear to be avoided within the CDS of mRNA, but also broadens the current understanding of synonymous codon usage as they relate to the structure-function relationship of RNA.

Ribosome Stoichiometry: From Form to Function

The existence of eukaryotic ribosomes with distinct ribosomal protein (RP) stoichiometry and regulatory roles in protein synthesis has been speculated for over 60 years. Recent advances in mass spectrometry (MS) and high-throughput analysis have begun to identify and characterize distinct ribosome stoichiometry in yeast and mammalian systems. In addition to RP stoichiometry, ribosomes host a vast array of protein modifications, effectively expanding the number of human RPs from 80 to many thousands of distinct proteoforms. Is it possible that these proteoforms combine to function as a 'ribosome code' to tune protein synthesis? We outline the specific benefits that translational regulation by specialized ribosomes can offer and discuss the means and methodologies available to correlate and characterize RP stoichiometry with function. We highlight previous research with a focus on formulating hypotheses that can guide future experiments and crack the ribosome code.

Analysis of synonymous codon usage pattern of genes in unique non-blood-sucking leech Whitmania pigra

Whitmania pigra is a unique, fluid-sucking ectoparasite and an anticoagulant medical leech. The codon usage bias (CUB) is the nonuniform usage of synonymous codons in which some codons are more preferred than others. Here, we performed a comprehensive analysis of CUB of genes in W. pigra, analyzing 140 780 transcripts, 59 553 unigenes, and 20 304 qualified coding sequences (CDSs) from the transcriptomic data of W. pigra. The effective number of codons values suggested that the CUB was low in these genes. We recognized profoundly favored codons in W. pigra that have a G/C-ending. Parity rule two-bias plots suggested that both mutation pressure and natural selection might have influenced the CUB. However, neutrality plots revealed that natural selection might have played a major role while mutation pressure might have played a minor role in shaping the CUB. We applied principal component analysis to relative synonymous codon usage values for divided CDSs based on GC content and codon-ending bases. Codon usage in W. pigra had a general inclination toward C-ending codons and natural selection rather than mutation pressure is the dominant force in the genetic evolution of W. pigra. To our knowledge, this is the first study to describe a complete codon usage analysis of W. pigra; this will increase the understanding of CUB and evolution in W. pigra. The analysis of codon usage patterns in W. pigra aids in understanding its evolution and genetic architecture.

Protein production in Escherichia coli is guided by the trade-off between intracellular substrate availability and energy cost

Background

In vivo protein formation is a crucial part of cellular life. The process needs to adapt to growth conditions and is exploited for the production of technical and pharmaceutical proteins in microbes such as Escherichia coli. Accordingly, the elucidation of basic regulatory mechanisms controlling the in vivo translation machinery is of primary interest, not only to improve heterologous protein production but also to elucidate fundamental regulation regimens of cellular growth.

Results

The current modeling analysis elucidates the impact of diffusion for the stochastic supply of crucial substrates such as the elongation factor EFTu, and tRNA species, all regarded as key elements for ensuring optimum transcriptional elongation. Together with the consideration of cellular ribosome numbers, their impact on the proper functioning of the translation machinery was investigated under different in vivo and in vitro conditions and utilizing the formation of non-native GFP and native EFTu as target proteins. The results show that translational elongation was diffusion limited. However, this effect was much more pronounced for the translation of non-native proteins than for the formation of codon-optimized native proteins.

Conclusions

Cellular ATP requirements constrain the options of improving protein production. In the case of non-native protein sequences, an optimized tRNA supply may be the most economical solution, as cells necessarily have to invest in ATP-costly ribosome synthesis to boost translation and increase growth rates.

Comparative Genomic Analysis Reveals the Distribution, Organization, and Evolution of Metal Resistance Genes in the Genus Acidithiobacillus

Members of the genus Acidithiobacillus, which can adapt to extremely high concentrations of heavy metals, are universally found at acid mine drainage (AMD) sites. Here, we performed a comparative genomic analysis of 37 strains within the genus Acidithiobacillus to answer the untouched questions as to the mechanisms and the evolutionary history of metal resistance genes in Acidithiobacillus spp. The results showed that the evolutionary history of metal resistance genes in Acidithiobacillus spp. involved a combination of gene gains and losses, horizontal gene transfer (HGT), and gene duplication. Phylogenetic analyses revealed that metal resistance genes in Acidithiobacillus spp. were acquired by early HGT events from species that shared habitats with Acidithiobacillus spp., such as Acidihalobacter, Thiobacillus, Acidiferrobacter, and Thiomonas species. Multicopper oxidase genes involved in copper detoxification were lost in iron-oxidizing Acidithiobacillus ferridurans, Acidithiobacillus ferrivorans, and Acidithiobacillus ferrooxidans and were replaced by rusticyanin genes during evolution. In addition, widespread purifying selection and the predicted high expression levels emphasized the indispensable roles of metal resistance genes in the ability of Acidithiobacillus spp. to adapt to harsh environments. Altogether, the results suggested that Acidithiobacillus spp. recruited and consolidated additional novel functionalities during the adaption to challenging environments via HGT, gene duplication, and purifying selection. This study sheds light on the distribution, organization, functionality, and complex evolutionary history of metal resistance genes in Acidithiobacillus spp.IMPORTANCE Horizontal gene transfer (HGT), natural selection, and gene duplication are three main engines that drive the adaptive evolution of microbial genomes. Previous studies indicated that HGT was a main adaptive mechanism in acidophiles to cope with heavy-metal-rich environments. However, evidences of HGT in Acidithiobacillus species in response to challenging metal-rich environments and the mechanisms addressing how metal resistance genes originated and evolved in Acidithiobacillus are still lacking. The findings of this study revealed a fascinating phenomenon of putative cross-phylum HGT, suggesting that Acidithiobacillus spp. recruited and consolidated additional novel functionalities during the adaption to challenging environments via HGT, gene duplication, and purifying selection. Altogether, the insights gained in this study have improved our understanding of the metal resistance strategies of Acidithiobacillus spp.

Unravelling Ribosome Function Through Structural Studies

Ribosomes are biological nanomachine that synthesise all proteins within a cell. It took decades to reveal the architecture of this essential cellular component. To understand the structure -function relationship of this nanomachine needed the utilisisation of different biochemical, biophysical and structural techniques. Structural studies combined with mutagenesis of the different ribosomal complexes comprising various RNAs and proteins enabled us to understand how this machine works inside a cell. Nowadays quite a number of ribosomal structures were published that confirmed biochemical studies on particular steps of protein synthesis by the ribosome . Four major steps were identified: initiation , elongation, termination and recycling. These steps lead us to the important question how the ribosome function can be regulated. Advances in technology for cryo electron microscopy: sample preparations, image recording, developments in algorithms for image analysis and processing significantly helped in revelation of structural details of the ribosome . We now have a library of ribosome structures from prokaryotes to eukaryotes that enable us to understand the complex mechanics of this nanomachine. As this structural library continues to grow, we gradually improve our understanding of this process and how it can be regulated and how the specific ribosomes can be stalled or activated, or completely disabled. This article provides a comprehensive overview of ribosomal structures that represent structural snapshots of the ribosome at its different functional states. Better understanding rises more particular questions that have to be addressed by determination structures of more complexes.Synopsis: Structural biology of the ribosome.

Synonymous Codon Usages as an Evolutionary Dynamic for Chlamydiaceae

The family of Chlamydiaceae contains a group of obligate intracellular bacteria that can infect a wide range of hosts. The evolutionary trend of members in this family is a hot topic, which benefits our understanding of the cross-infection of these pathogens. In this study, 14 whole genomes of 12 Chlamydia species were used to investigate the nucleotide, codon, and amino acid usage bias by synonymous codon usage value and information entropy method. The results showed that all the studied Chlamydia spp. had A/T rich genes with over-represented A or T at the third positions and G or C under-represented at these positions, suggesting that nucleotide usages influenced synonymous codon usages. The overall codon usage trend from synonymous codon usage variations divides the Chlamydia spp. into four separate clusters, while amino acid usage divides the Chlamydia spp. into two clusters with some exceptions, which reflected the genetic diversity of the Chlamydiaceae family members. The overall codon usage pattern represented by the effective number of codons (ENC) was significantly positively correlated to gene GC3 content. A negative correlation exists between ENC and the codon adaptation index for some Chlamydia species. These results suggested that mutation pressure caused by nucleotide composition constraint played an important role in shaping synonymous codon usage patterns. Furthermore, codon usage of T3ss and Pmps gene families adapted to that of the corresponding genome. Taken together, analyses help our understanding of evolutionary interactions between nucleotide, synonymous codon, and amino acid usages in genes of Chlamydiaceae family members.

Presyncodon, a Web Server for Gene Design with the Evolutionary Information of the Expression Hosts

In the natural host, most of the synonymous codons of a gene have been evolutionarily selected and related to protein expression and function. However, for the design of a new gene, most of the existing codon optimization tools select the high-frequency-usage codons and neglect the contribution of the low-frequency-usage codons (rare codons) to the expression of the target gene in the host. In this study, we developed the method Presyncodon, available in a web version, to predict the gene code from a protein sequence, using built-in evolutionary information on a specific expression host. The synonymous codon-usage pattern of a peptide was studied from three genomic datasets (Escherichia coli, Bacillus subtilis, and Saccharomyces cerevisiae). Machine-learning models were constructed to predict a selection of synonymous codons (low- or high-frequency-usage codon) in a gene. This method could be easily and efficiently used to design new genes from protein sequences for optimal expression in three expression hosts (E. coli, B. subtilis, and S. cerevisiae). Presyncodon is free to academic and noncommercial users; accessible at http://www.mobioinfor.cn/presyncodon_www/index.html.