Codon usage determines translation rate in Escherichia coli

Thermal adaptation in plants: understanding the dynamics of translation factors and condensates

Plants, as sessile organisms, face the crucial challenge of adjusting growth and development with ever-changing environmental conditions. Protein synthesis is the fundamental process that enables growth of all organisms. Since elevated temperature presents a substantial threat to protein stability and function, immediate adjustments of protein synthesis rates are necessary to circumvent accumulation of proteotoxic stress and to ensure survival. This review provides an overview of the mechanisms that control translation under high-temperature stress by the modification of components of the translation machinery in plants, and compares them to yeast and metazoa. Recent research also suggests an important role for cytoplasmic biomolecular condensates, named stress granules, in these processes. Current understanding of the role of stress granules in translational regulation and of the molecular processes associated with translation that might occur within stress granules is also discussed.

Differential Selection for Translation Efficiency Shapes Translation Machineries in Bacterial Species

Different bacterial species have dramatically different generation times, from 20-30 min in Escherichia coli to about two weeks in Mycobacterium leprae. The translation machinery in a cell needs to synthesize all proteins for a new cell in each generation. The three subprocesses of translation, i.e., initiation, elongation, and termination, are expected to be under stronger selection pressure to optimize in short-generation bacteria (SGB) such as Vibrio natriegens than in the long-generation Mycobacterium leprae. The initiation efficiency depends on the start codon decoded by the initiation tRNA, the optimal Shine-Dalgarno (SD) decoded by the anti-SD (aSD) sequence on small subunit rRNA, and the secondary structure that may embed the initiation signals and prevent them from being decoded. The elongation efficiency depends on the tRNA pool and codon usage. The termination efficiency in bacteria depends mainly on the nature of the stop codon and the nucleotide immediately downstream of the stop codon. By contrasting SGB with long-generation bacteria (LGB), we predict (1) SGB to have more ribosome RNA operons to produce ribosomes, and more tRNA genes for carrying amino acids to ribosomes, (2) SGB to have a higher percentage of genes using AUG as the start codon and UAA as the stop codon than LGB, (3) SGB to exhibit better codon and anticodon adaptation than LGB, and (4) SGB to have a weaker secondary structure near the translation initiation signals than LGB. These differences between SGB and LGB should be more pronounced in highly expressed genes than the rest of the genes. We present empirical evidence in support of these predictions.

Modifications of the 5' region of the CASPONTM tag's mRNA further enhance soluble recombinant protein production in Escherichia coli

BACKGROUND

Escherichia coli is one of the most commonly used host organisms for the production of biopharmaceuticals, as it allows for cost-efficient and fast recombinant protein expression. However, challenging proteins are often produced with low titres or as inclusion bodies, and the manufacturing process needs to be developed individually for each protein. Recently, we developed the CASPONTM technology, a generic fusion tag-based platform process for high-titer soluble expression including a standardized downstream processing and highly specific enzymatic cleavage of the fusion tag. To assess potential strategies for further improvement of the N-terminally fused CASPONTM tag, we modified the 5'UTR and 5' region of the tag-coding mRNA to optimize the ribosome-mRNA interactions.

RESULTS

In the present work, we found that by modifying the 5'UTR sequence of a pET30acer plasmid-based system, expression of the fusion protein CASPONTM-tumour necrosis factor α was altered in laboratory-scale carbon-limited fed-batch cultivations, but no significant increase in expression titre was achieved. Translation efficiency was highest for a construct carrying an expression enhancer element and additionally possessing a very favourable interaction energy between ribosome and mRNA (∆Gtotal). However, a construct with comparatively low transcriptional efficiency, which lacked the expression enhancer sequence and carried the most favourable ∆Gtotal tested, led to the highest recombinant protein formation alongside the reference pET30a construct. Furthermore, we found, that by introducing synonymous mutations within the nucleotide sequence of the T7AC element of the CASPONTM tag, utilizing a combination of rare and non-rare codons, the free folding energy of the nucleotides at the 5' end (-4 to + 37) of the transcript encoding the CASPONTM tag increased by 6 kcal/mol. Surprisingly, this new T7ACrare variant led to improved recombinant protein titres by 1.3-fold up to 5.3-fold, shown with three industry-relevant proteins in lab-scale carbon limited fed-batch fermentations under industrially relevant conditions.

CONCLUSIONS

This study reveals some of the complex interdependencies between the ribosome and mRNA that govern recombinant protein expression. By modifying the 5'UTR to obtain an optimized interaction energy between the mRNA and the ribosome (ΔGtotal), transcript levels were changed, highlighting the different translation efficiencies of individual transcripts. It was shown that the highest recombinant titre was not obtained by the construct with the most efficient translation but by a construct with a generally high transcript amount coupled with a favourable ΔGtotal. Furthermore, an unexpectedly high potential to enhance expression by introducing silent mutations including multiple rare codons into the 5'end of the CAPONTM tag's mRNA was identified. Although the titres of the fusion proteins were dramatically increased, no formation of inclusion bodies or negative impact on cell growth was observed. We hypothesize that the drastic increase in titre is most likely caused by better ribosomal binding site accessibility. Our study, which demonstrates the influence of changes in ribosome-mRNA interactions on protein expression under industrially relevant production conditions, opens the door to the applicability of the new T7ACrare tag in biopharmaceutical industry using the CASPONTM platform process.

SARS-CoV-2 continuously optimizes its codon usage to adapt to human lung environment

Viruses need to utilize the resources from host cells to reproduce themselves. RNA translation rate, which is largely determined by codon usage, is the rate-limiting step across the life cycle of viruses. Adapting to the codon usage of hosts would help virus better proliferate. We retrieved the time-course mutation profile of millions of world-wide SARS-CoV-2 sequences. For synonymous mutations, we defined whether a mutation elevate or reduce the relative synonymous codon usage (RSCU). We found that if a synonymous mutation in SARS-CoV-2 increases the RSCU (calculated from human lungs), denoted as delta RSCU > 0, then this mutation is positively selected because the allele frequency (AF) of this mutation increases with time, and vice versa. The results suggest that in SARS-CoV-2, the synonymous mutations that increase codon optimality are beneficial to the virus and are favored by natural selection. For the first time, we used the dynamics of allele frequency to demonstrate that SARS-CoV-2 is continuously optimizing its codon usage to adapt to human lungs. Nevertheless, adaptation to other human tissues cannot be excluded. These results warn us that under this global pandemic, synonymous mutations in SARS-CoV-2 should not be automatically ignored since they indeed change the fitness of the virus.

Maximizing Heterologous Expression of Engineered Type I Polyketide Synthases: Investigating Codon Optimization Strategies

Type I polyketide synthases (T1PKSs) hold enormous potential as a rational production platform for the biosynthesis of specialty chemicals. However, despite great progress in this field, the heterologous expression of PKSs remains a major challenge. One of the first measures to improve heterologous gene expression can be codon optimization. Although controversial, choosing the wrong codon optimization strategy can have detrimental effects on the protein and product levels. In this study, we analyzed 11 different codon variants of an engineered T1PKS and investigated in a systematic approach their influence on heterologous expression in Corynebacterium glutamicum, Escherichia coli, and Pseudomonas putida. Our best performing codon variants exhibited a minimum 50-fold increase in PKS protein levels, which also enabled the production of an unnatural polyketide in each of these hosts. Furthermore, we developed a free online tool (https://basebuddy.lbl.gov) that offers transparent and highly customizable codon optimization with up-to-date codon usage tables. In this work, we not only highlight the significance of codon optimization but also establish the groundwork for the high-throughput assembly and characterization of PKS pathways in alternative hosts.

Bioinformatic Analysis Reveals the Role of Translation Elongation Efficiency Optimisation in the Evolution of Ralstonia Genus

Translation efficiency modulates gene expression in prokaryotes. The comparative analysis of translation elongation efficiency characteristics of Ralstonia genus bacteria genomes revealed that these characteristics diverge in accordance with the phylogeny of Ralstonia. The first branch of this genus is a group of bacteria commonly found in moist environments such as soil and water that includes the species R. mannitolilytica, R. insidiosa, and R. pickettii, which are also described as nosocomial infection pathogens. In contrast, the second branch is plant pathogenic bacteria consisting of R. solanacearum, R. pseudosolanacearum, and R. syzygii. We found that the soil Ralstonia have a significantly lower number and energy of potential secondary structures in mRNA and an increased role of codon usage bias in the optimization of highly expressed genes' translation elongation efficiency, not only compared to phytopathogenic Ralstonia but also to Cupriavidus necator, which is closely related to the Ralstonia genus. The observed alterations in translation elongation efficiency of orthologous genes are also reflected in the difference of potentially highly expressed gene' sets' content among Ralstonia branches with different lifestyles. Analysis of translation elongation efficiency characteristics can be considered a promising approach for studying complex mechanisms that determine the evolution and adaptation of bacteria in various environments.

BamQuery: a proteogenomic tool to explore the immunopeptidome and prioritize actionable tumor antigens

MHC-I-associated peptides deriving from non-coding genomic regions and mutations can generate tumor-specific antigens, including neoantigens. Quantifying tumor-specific antigens' RNA expression in malignant and benign tissues is critical for discriminating actionable targets. We present BamQuery, a tool attributing an exhaustive RNA expression to MHC-I-associated peptides of any origin from bulk and single-cell RNA-sequencing data. We show that many cryptic and mutated tumor-specific antigens can derive from multiple discrete genomic regions, abundantly expressed in normal tissues. BamQuery can also be used to predict MHC-I-associated peptides immunogenicity and identify actionable tumor-specific antigens de novo.

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression

Differences in codon frequency between genomes, genes, or positions along a gene, modulate transcription and translation efficiency, leading to phenotypic and functional differences. Here, we present a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells, quantifying the phenotypic consequences of codon usage bias at different molecular and cellular levels, with an emphasis on translation elongation. Six synonymous versions of an antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was analyzed by means of quantitative transcriptome and proteome assessment, as proxies for gene expression; cellular fluorescence, as a proxy for single-cell level expression; and real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness. We show that differences in codon usage bias strongly impact the molecular and cellular phenotype: (i) they result in large differences in mRNA levels and protein levels, leading to differences of over 15 times in translation efficiency; (ii) they introduce unpredicted splicing events; (iii) they lead to reproducible phenotypic heterogeneity; and (iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression. In human cells in culture, codon usage bias modulates gene expression by modifying mRNA availability and suitability for translation, leading to differences in protein levels and eventually eliciting functional phenotypic changes.

Three Novel Antisense Overlapping Genes in E. coli O157:H7 EDL933

The abundance of long overlapping genes in prokaryotic genomes is likely to be significantly underestimated. To date, only a few examples of such genes are fully established. Using RNA sequencing and ribosome profiling, we found expression of novel overlapping open reading frames in Escherichia coli O157:H7 EDL933 (EHEC). Indeed, the overlapping candidate genes are equipped with typical structural elements required for transcription and translation, i.e., promoters, transcription start sites, as well as terminators, all of which were experimentally verified. Translationally arrested mutants, unable to produce the overlapping encoded protein, were found to have a growth disadvantage when grown competitively against the wild type. Thus, the phenotypes found imply biological functionality of the genes at the level of proteins produced. The addition of 3 more examples of prokaryotic overlapping genes to the currently limited, yet constantly growing pool of such genes emphasizes the underestimated coding capacity of bacterial genomes. IMPORTANCE The abundance of long overlapping genes in prokaryotic genomes is likely to be significantly underestimated, since such genes are not allowed in genome annotations. However, ribosome profiling catches mRNA in the moment of being template for protein production. Using this technique and subsequent experiments, we verified 3 novel overlapping genes encoded in antisense of known genes. This adds more examples of prokaryotic overlapping genes to the currently limited, yet constantly growing pool of such genes.

Poly-Histidine-Tagged Protein Purification Using Immobilized Metal Affinity Chromatography (IMAC)

His-tagging is the most widespread and versatile strategy used to purify recombinant proteins for biochemical and structural studies. Recombinant DNA methods are first used to engineer the addition of a short tract of poly-histidine tag (His-tag) to the N-terminus or C-terminus of a target protein. The His-tag is then exploited to enable purification of the "tagged" protein by immobilized metal affinity chromatography (IMAC). In this chapter, we describe efficient procedures for the isolation of highly purified His-tagged target proteins from an Escherichia coli host using IMAC in a bind-wash-elute strategy that can be performed under both native and denaturing conditions.

Mechanisms governing codon usage bias and the implications for protein expression in the chloroplast of Chlamydomonas reinhardtii

Chloroplasts possess a considerably reduced genome that is decoded via an almost minimal set of tRNAs. These features make an excellent platform for gaining insights into fundamental mechanisms that govern protein expression. Here, we present a comprehensive and revised perspective of the mechanisms that drive codon selection in the chloroplast of Chlamydomonas reinhardtii and the functional consequences for protein expression. In order to extract this information, we applied several codon usage descriptors to genes with different expression levels. We show that highly expressed genes strongly favor translationally optimal codons, while genes with lower functional importance are rather affected by directional mutational bias. We demonstrate that codon optimality can be deduced from codon-anticodon pairing affinity and, for a small number of amino acids (leucine, arginine, serine, and isoleucine), tRNA concentrations. Finally, we review, analyze, and expand on the impact of codon usage on protein yield, secondary structures of mRNA, translation initiation and termination, and amino acid composition of proteins, as well as cotranslational protein folding. The comprehensive analysis of codon choice provides crucial insights into heterologous gene expression in the chloroplast of C. reinhardtii, which may also be applicable to other chloroplast-containing organisms and bacteria.

Bioinformatic Assessment of Factors Affecting the Correlation between Protein Abundance and Elongation Efficiency in Prokaryotes

Protein abundance is crucial for the majority of genetically regulated cell functions to act properly in prokaryotic organisms. Therefore, developing bioinformatic methods for assessing the efficiency of different stages of gene expression is of great importance for predicting the actual protein abundance. One of these steps is the evaluation of translation elongation efficiency based on mRNA sequence features, such as codon usage bias and mRNA secondary structure properties. In this study, we have evaluated correlation coefficients between experimentally measured protein abundance and predicted elongation efficiency characteristics for 26 prokaryotes, including non-model organisms, belonging to diverse taxonomic groups The algorithm for assessing elongation efficiency takes into account not only codon bias, but also number and energy of secondary structures in mRNA if those demonstrate an impact on predicted elongation efficiency of the ribosomal protein genes. The results show that, for a number of organisms, secondary structures are a better predictor of protein abundance than codon usage bias. The bioinformatic analysis has revealed several factors associated with the value of the correlation coefficient. The first factor is the elongation efficiency optimization type-the organisms whose genomes are optimized for codon usage only have significantly higher correlation coefficients. The second factor is taxonomical identity-bacteria that belong to the class Bacilli tend to have higher correlation coefficients among the analyzed set. The third is growth rate, which is shown to be higher for the organisms with higher correlation coefficients between protein abundance and predicted translation elongation efficiency. The obtained results can be useful for further improvement of methods for protein abundance prediction.

The High Mutational Sensitivity of ccdA Antitoxin Is Linked to Codon Optimality

Deep mutational scanning studies suggest that synonymous mutations are typically silent and that most exposed, nonactive-site residues are tolerant to mutations. Here, we show that the ccdA antitoxin component of the Escherichia coli ccdAB toxin-antitoxin system is unusually sensitive to mutations when studied in the operonic context. A large fraction (∼80%) of single-codon mutations, including many synonymous mutations in the ccdA gene shows inactive phenotype, but they retain native-like binding affinity towards cognate toxin, CcdB. Therefore, the observed phenotypic effects are largely not due to alterations in protein structure/stability, consistent with a large region of CcdA being intrinsically disordered. E. coli codon preference and strength of ribosome-binding associated with translation of downstream ccdB gene are found to be major contributors of the observed ccdA mutant phenotypes. In select cases, proteomics studies reveal altered ratios of CcdA:CcdB protein levels in vivo, suggesting that the ccdA mutations likely alter relative translation efficiencies of the two genes in the operon. We extend these results by studying single-site synonymous mutations that lead to loss of function phenotypes in the relBE operon upon introduction of rarer codons. Thus, in their operonic context, genes are likely to be more sensitive to both synonymous and nonsynonymous point mutations than inferred previously.

Characterization and functional of four mutants of hydroxy fatty acid dehydrogenase from Lactobacillus plantarum p-8

In this study, the hydroxy fatty acid dehydrogenase CLA-DH from Lactobacillus plantarum p-8 and its four mutant variants were expressed in Escherichia coli Rosetta (DE3). UV spectrophotometry was employed to verify the catalytic power of the purified CLA-DH to convert ricinoleic acid into 12-oxo-cis-9-octadecenoic acid in the presence of oxidized nicotinamide adenine dinucleotide (NAD+). The optimum reaction temperature for CLA-DH was 45°C, with a maintained stability between 20°C and 40°C. The optimal pH for CLA-DH catalytic activity was 6.0-7.0, with a maintained stability at a pH range of 6.0-8.0. In addition, Fe3+ promoted enzyme activity, whereas Cu2+, Zn2+, and Fe2+ inhibited enzyme activity (P < 0.05). The Km, Vmax, Kcat, and Kcat/Km of CLA-DH were determined as 2.19 ± 0.34 μM, 2.06 ± 0.28 μM min-1, 2.00 ± 0.27 min-1, and 0.92 ± 0.02 min-1μM-1, respectively. Site-directed mutagenesis and molecular dynamics simulations demonstrated that both Tyr156 and Ser143 residues play significant roles in the catalysis of CLA-DH, and its solubility is affected by Lys160 and Asp63. Moreover, Gas chromatography determined that recombinant CLA-DH could be successfully applied to Conjugated linoleic acids production.

mRNA and tRNA modification states influence ribosome speed and frame maintenance during poly(lysine) peptide synthesis

Ribosome speed is dictated by multiple factors including substrate availability, cellular conditions, and product (peptide) formation. Translation slows during the synthesis of cationic peptide sequences, potentially influencing the expression of thousands of proteins. Available evidence suggests that ionic interactions between positively charged nascent peptides and the negatively charged ribosome exit tunnel impede translation. However, this hypothesis was difficult to test directly because of inability to decouple the contributions of amino acid charge from mRNA sequence and tRNA identity/abundance in cells. Furthermore, it is unclear if other components of the translation system central to ribosome function (e.g., RNA modification) influence the speed and accuracy of positively charged peptide synthesis. In this study, we used a fully reconstituted Escherichia coli translation system to evaluate the effects of peptide charge, mRNA sequence, and RNA modification status on the translation of lysine-rich peptides. Comparison of translation reactions on poly(lysine)-encoding mRNAs conducted with either Lys-tRNA^Lys or Val-tRNA^Lys reveals that that amino acid charge, while important, only partially accounts for slowed translation on these transcripts. We further find that in addition to peptide charge, mRNA sequence and both tRNA and mRNA modification status influence the rates of amino acid addition and the ribosome’s ability to maintain frame (instead of entering the −2, −1, and +1 frames) during poly(lysine) peptide synthesis. Our observations lead us to expand the model for explaining how the ribosome slows during poly(lysine) peptide synthesis and suggest that posttranscriptional RNA modifications can provide cells a mechanism to precisely control ribosome movements along an mRNA.

The genomic basis of host and vector specificity in non-pathogenic trypanosomatids

Trypanosoma theileri, a non-pathogenic parasite of bovines, has a predicted surface protein architecture that likely aids survival in its mammalian host. Their surface proteins are encoded by genes which account for ∼10% of their genome. A non-pathogenic parasite of sheep, Trypanosoma melophagium, is transmitted by the sheep ked and is closely related to T. theileri. To explore host and vector specificity between these species, we sequenced the T. melophagium genome and transcriptome and an annotated draft genome was assembled. T. melophagium was compared to 43 kinetoplastid genomes, including T. theileri. T. melophagium and T. theileri have an AT biased genome, the greatest bias of publicly available trypanosomatids. This trend may result from selection acting to decrease the genomic nucleotide cost. The T. melophagium genome is 6.3Mb smaller than T. theileri and large families of proteins, characteristic of the predicted surface of T. theileri, were found to be absent or greatly reduced in T. melophagium. Instead, T. melophagium has modestly expanded protein families associated with the avoidance of complement-mediated lysis. We propose that the contrasting genomic features of these species is linked to their mode of transmission from their insect vector to their mammalian host. This article has an associated First Person interview with the first author of the paper.

Unexpected myriad of co-occurring viral strains and species in one of the most abundant and microdiverse viruses on Earth

Viral genetic microdiversity drives adaptation, pathogenicity, and speciation and has critical consequences for the viral-host arms race occurring at the strain and species levels, which ultimately impact microbial community structure and biogeochemical cycles. Despite the fact that most efforts have focused on viral macrodiversity, little is known about the microdiversity of ecologically important viruses on Earth. Recently, single-virus genomics discovered the putatively most abundant ocean virus in temperate and tropical waters: the uncultured dsDNA virus vSAG 37-F6 infecting Pelagibacter, the most abundant marine bacteria. In this study, we report the cooccurrence of up to ≈1,500 different viral strains (>95% nucleotide identity) and ≈30 related species (80-95% nucleotide identity) in a single oceanic sample. Viral microdiversity was maintained over space and time, and most alleles were the result of synonymous mutations without any apparent adaptive benefits to cope with host translation codon bias and efficiency. Gene flow analysis used to delimitate species according to the biological species concept (BSC) revealed the impact of recombination in shaping vSAG 37-F6 virus and Pelagibacter speciation. Data demonstrated that this large viral microdiversity somehow mirrors the host species diversity since ≈50% of the 926 analyzed Pelagibacter genomes were found to belong to independent BSC species that do not significantly engage in gene flow with one another. The host range of this evolutionarily successful virus revealed that a single viral species can infect multiple Pelagibacter BSC species, indicating that this virus crosses not only formal BSC barriers but also biomes since viral ancestors are found in freshwater.

Rifampicin Increases Expression of Plant Codon-Optimized Bacillus thuringiensis δ-Endotoxin Genes in Escherichia coli

Transgenic crops expressing Cry δ-endotoxins of Bacillus thuringiensis for insect resistance have been commercialized worldwide with increased crop productivity and spectacular socioeconomic gains. To attain the enhanced level of protein expression, the cry genes have to be extensively modified for RNA stability and translation efficiency in the plant systems. However, such modifications in nucleotide sequences make it difficult to express the cry genes in Escherichia coli because of the presence of E. coli rare codons. Induction of gene expression through the T7 promoter/lac operator system results in high levels of transcription but limits the availability of activated tRNA corresponding to rare codons that leads to translation stalling at ribosomes. In the present study, an Isopropyl ß-D-1-thiogalactopyranoside (IPTG)/rifampicin combination-based approach was adopted to induce transcription of cry genes through T7 promoter/lac operator while simultaneously inhibiting the transcription of host genes through rifampicin. The results show that the IPTG/rifampicin combination leads to high-level expression of four plant codon-optimized cry genes (cry2Aa, cry1F, cry1Ac, and cry1AcF). Northern blot analysis of the cry gene expressing E. coli samples showed that the RNA expression level in the IPTG-induced samples was higher as compared to that in the IPTG/rifampicin-induced samples. Diet overlay insect bioassay of IPTG/rifampicin-induced Cry toxins with Helicoverpa armigera larvae showed bioactivity (measured as LC₅₀) similar to the previous studies. The experiment has proved that recombinant synthetic gene (plant codon-optimized gene) with the combination of Rifampicin which inhibits DNA-dependent bacterial RNA polymerase and reduces the excessive baggage of translational machinery of the bacterial cell triggers the production of synthetic protein. Purification of protein using high pH buffer increases the solubility of the protein. Further, LC₅₀ analysis shows no reduction of protein activity leads to protein stability. Further, purified cry toxin protein can be used for crop protection against pests and a purified form of the synthetic protein can be used for antibody production and perform the immunoassay for the identification of the transgenic plant. The crystallographic structure of synthetic protein could be used for interaction study with another insect to see insecticidal activity.

Codon Usage Bias in Phytoplankton

Non-random usage of synonymous codons, known as “codon bias”, has been described in many organisms, from bacteria to Drosophila, but little is known about it in phytoplankton. This phenomenon is thought to be driven by selection for translational efficiency. As the efficacy of selection is proportional to the effective population size, species with large population sizes, such as phytoplankton, are expected to have strong codon bias. To test this, we measured codon bias in 215 strains from Haptophyta, Chlorophyta, Ochrophyta (except diatoms that were studied previously), Dinophyta, Cryptophyta, Ciliophora, unicellular Rhodophyta and Chlorarachniophyta. Codon bias is modest in most groups, despite the astronomically large population sizes of marine phytoplankton. The strength of the codon bias, measured with the effective number of codons, is the strongest in Haptophyta and the weakest in Chlorarachniophyta. The optimal codons are GC-ending in most cases, but several shifts to AT-ending codons were observed (mainly in Ochrophyta and Ciliophora). As it takes a long time to reach a new equilibrium after such shifts, species having AT-ending codons show a lower frequency of optimal codons compared to other species. Genetic diversity, calculated for species with more than three strains sequenced, is modest, indicating that the effective population sizes are many orders of magnitude lower than the astronomically large census population sizes, which helps to explain the modest codon bias in marine phytoplankton. This study represents the first comparative analysis of codon bias across multiple major phytoplankton groups.

Variation in Release Factor Abundance Is Not Needed to Explain Trends in Bacterial Stop Codon Usage

In bacteria stop codons are recognized by one of two class I release factors (RF1) recognizing TAG, RF2 recognizing TGA, and TAA being recognized by both. Variation across bacteria in the relative abundance of RF1 and RF2 is thus hypothesized to select for different TGA/TAG usage. This has been supported by correlations between TAG:TGA ratios and RF1:RF2 ratios across multiple bacterial species, potentially also explaining why TAG usage is approximately constant despite extensive variation in GC content. It is, however, possible that stop codon trends are determined by other forces and that RF ratios adapt to stop codon usage, rather than vice versa. Here, we determine which direction of the causal arrow is the more parsimonious. Our results support the notion that RF1/RF2 ratios become adapted to stop codon usage as the same trends, notably the anomalous TAG behavior, are seen in contexts where RF1:RF2 ratios cannot be, or are unlikely to be, causative, that is, at 3'untranslated sites never used for translation termination, in intragenomic analyses, and across archaeal species (that possess only one RF1). We conclude that specifics of RF biology are unlikely to fully explain TGA/TAG relative usage. We discuss why the causal relationships for the evolution of synonymous stop codon usage might be different from those affecting synonymous sense codon usage, noting that transitions between TGA and TAG require two-point mutations one of which is likely to be deleterious.

Identification of conserved slow codons that are important for protein expression and function

ABSTRASTDue to the redundancy of the genetic code most amino acids are encoded by several 'synonymous' codons. These codons are used unevenly, and each organism demonstrates its own unique codon usage bias, where the 'preferred' codons are associated with tRNAs that are found in high concentrations. Therefore, for decades, the prevailing view had been that preferred and non-preferred codons are linked to high or slow translation rates, respectively.However, this simplified view is contrasted by the frequent failures of codon-optimization efforts and by evidence of non-preferred (i.e. 'slow') codons having specific roles important for efficient production of functional proteins. One such specific role of slower codons is the regulation of co-translational protein folding, a complex biophysical process that is very challenging to model or to measure.Here, we combined a genome-wide approach with experiments to investigate the role of slow codons in protein production and co-translational folding. We analysed homologous gene groups from divergent bacteria and identified positions of inter-species conservation of bias towards slow codons. We then generated mutants where the conserved slow codons are substituted with 'fast' ones, and experimentally studied the effects of these codon substitutions. Using cellular and biochemical approaches we find that at certain locations, slow-to-fast codon substitutions reduce protein expression, increase protein aggregation, and impair protein function.This report provides an approach for identifying functionally relevant regions with slower codons and demonstrates that such codons are important for protein expression and function.

Codon Optimization, Soluble Expression and Purification of PE_PGRS45 Gene from Mycobacterium tuberculosis and Preparation of Its Polyclonal Antibody Protein

Studies have demonstrated that PE_PGRS45 is constitutively expressed under various environmental conditions (such as nutrient depletion, hypoxia, and low pH) of the in vitro growth conditions examined, indicating that PE_PGRS45 protein is critical to the basic functions of Mycobacterium tuberculosis. However, there are few reports about the biochemical function and pathogenic mechanism of PE_PGRS45 protein. The fact that this M. tuberculosis gene is not easily expressed in E. coli may be mainly due to the high content of G+C and the use of unique codons. Fusion tags are indispensable tools used to improve the soluble expression of recombinant proteins and accelerate the characterization of protein structure and function. In the present study, His6, Trx, and His6-MBP were used as fusion tags, but only MBP-PE_PGRS45 was expressed solubly. The purification using His6-MBP tag-specific binding to the Ni column was easy to separate after the tag cleavage. We used the purified PE_PGRS45 to immunize New Zealand rabbits and obtained anti- PE_PGRS45 serum. We found that the titer of polyclonal antibodies against PE_PGR45 was higher than 1:256000. The result shows that purified PE_PGRS45 can induce New Zealand rabbits to produce high-titer antibodies. In conclusion, the recombinant protein PE_PGRS45 was successfully expressed in E. coli and specific antiserum was prepared, which will be followed by further evaluation of these specific antigens to develop highly sensitive and specific diagnostic tests for tuberculosis.

Engineering gene overlaps to sustain genetic constructs in vivo

Evolution is often an obstacle to the engineering of stable biological systems due to the selection of mutations inactivating costly gene circuits. Gene overlaps induce important constraints on sequences and their evolution. We show that these constraints can be harnessed to increase the stability of costly genes by purging loss-of-function mutations. We combine computational and synthetic biology approaches to rationally design an overlapping reading frame expressing an essential gene within an existing gene to protect. Our algorithm succeeded in creating overlapping reading frames in 80% of E. coli genes. Experimentally, scoring mutations in both genes of such overlapping construct, we found that a significant fraction of mutations impacting the gene to protect have a deleterious effect on the essential gene. Such an overlap thus protects a costly gene from removal by natural selection by associating the benefit of this removal with a larger or even lethal cost. In our synthetic constructs, the overlap converts many of the possible mutants into evolutionary dead-ends, reducing the evolutionary potential of the system and thus increasing its stability over time.

Genomes are translated by triplets of nucleotides on two different strands, allowing for six different reading frames. This permits the existence of gene overlaps, often observed in microbial genomes, where two different proteins are encoded on the same piece of DNA, but in different reading frames. Gene overlaps are classically considered an obstacle for both evolution and genetic engineering, as mutations in overlapping regions likely have pleitrotropic effects on several genes. In 2013, we identified specific evolutionary scenarios where the decrease in evolutionary potential caused by gene overlaps could instead be advantageous and selected for. In this work, we demonstrate the use of gene overlaps in another context where reducing evolutionary potential can be useful: preventing evolution from inactivating synthetic circuits. We show that gene overlaps can be engineered to increase the evolutionary stability of genes that are costly to their hosts, by entangling these costly genes with essential genes.

Compelling Evidence Suggesting the Codon Usage of SARS-CoV-2 Adapts to Human After the Split From RaTG13

SARS-CoV-2 needs to efficiently make use of the resources from hosts in order to survive and propagate. Among the multiple layers of regulatory network, mRNA translation is the rate-limiting step in gene expression. Synonymous codon usage usually conforms with tRNA concentration to allow fast decoding during translation. It is acknowledged that SARS-CoV-2 has adapted to the codon usage of human lungs so that the virus could rapidly proliferate in the lung environment. While this notion seems to nicely explain the adaptation of SARS-CoV-2 to lungs, it is unable to tell why other viruses do not have this advantage. In this study, we retrieve the GTEx RNA-seq data for 30 tissues (belonging to over 17 000 individuals). We calculate the RSCU (relative synonymous codon usage) weighted by gene expression in each human sample, and investigate the correlation of RSCU between the human tissues and SARS-CoV-2 or RaTG13 (the closest coronavirus to SARS-CoV-2). Lung has the highest correlation of RSCU to SARS-CoV-2 among all tissues, suggesting that the lung environment is generally suitable for SARS-CoV-2. Interestingly, for most tissues, SARS-CoV-2 has higher correlations with the human samples compared with the RaTG13-human correlation. This difference is most significant for lungs. In conclusion, the codon usage of SARS-CoV-2 has adapted to human lungs to allow fast decoding and translation. This adaptation probably took place after SARS-CoV-2 split from RaTG13 because RaTG13 is less perfectly correlated with human. This finding depicts the trajectory of adaptive evolution from ancestral sequence to SARS-CoV-2, and also well explains why SARS-CoV-2 rather than other viruses could perfectly adapt to human lung environment.

Bicistronic design as recombinant expression enhancer: characteristics, applications, and structural optimization

The bicistronic design (BCD) is characterized by a short fore-cistron sequence and a second Shine-Dalgarno (SD2) sequence upstream of the target gene. The outstanding performance of this expression cassette in promoting recombinant protein production has attracted attention. Recently, the application of the BCD has been further extended to gene expression control, protein translation monitoring, and membrane protein production. In this review, we summarize the characteristics, molecular mechanisms, applications, and structural optimization of the BCD expression cassette. We also specifically discuss the challenges that the BCD system still faces. This is the first review of the BCD expression strategy, and it is believed that an in-depth understanding of the BCD will help researchers to better utilize and develop it. KEY POINTS: • Summary of the characteristics and molecular mechanisms of the BCD system. • Review of the actual applications of the BCD expression cassette. • Summary of the structural optimization of the BCD system.

Successful expression of the synthetic merBps gene in tobacco

Organomercury is the most toxic biomagnifiable state of mercury, and to date, no natural organomercurial detoxification mechanism is encountered in plants. Bacterial merB gene encoding organomercury lyase show low expression in transgenic plants. For ideal expression, a synthetic merBps gene possessing143 out of 213 codons discrete from native merB gene from Escherichia. coli was fabricated based on codon usage in tobacco. Through Agrobacterium-mediated transformation, the merBps gene got successfully integrated into tobacco. Of several putative merBps transformants selected with 200 μg ml-1 kanamycin, only ∼45% were PCR positive for both nptII and merBps genes. Healthy and vigorously growing shoots of few PCR-positive putative transgenic lines were multiplied and rooted. After transplantation and acclimatization, the resultant plants flowered and fruited in pots. Southern analysis revealed the presence of a single copy of the merBps gene in four lines. RT-PCR and Western investigations established successful transcription and translation of the merBps gene in these transgenic lines, respectively. Fabrication of fully functional organomercury lyase in merBps transgenic lines was established based on the potential of their (i) seeds to germinate; (ii) shoots to grow and multiply; and (iii) leaf disc to remain green, even in the presence of 4 nmole ml-1 phenylmercuryacetate (PMA) while the wild type was susceptible to even 1 nmole ml-1 PMA. These findings confirmed that the synthetic merBps gene could be effectively expressed in plants and exploited for remediation of organomercurial contaminated sites.

Variability in Codon Usage in Coronaviruses Is Mainly Driven by Mutational Bias and Selective Constraints on CpG Dinucleotide

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third human-emerged virus of the 21st century from the Coronaviridae family, causing the ongoing coronavirus disease 2019 (COVID-19) pandemic. Due to the high zoonotic potential of coronaviruses, it is critical to unravel their evolutionary history of host species breadth, host-switch potential, adaptation and emergence, to identify viruses posing a pandemic risk in humans. We present here a comprehensive analysis of the composition and codon usage bias of the 82 Orthocoronavirinae members, infecting 47 different avian and mammalian hosts. Our results clearly establish that synonymous codon usage varies widely among viruses, is only weakly dependent on their primary host, and is dominated by mutational bias towards AU-enrichment and by CpG avoidance. Indeed, variation in GC3 explains around 34%, while variation in CpG frequency explains around 14% of total variation in codon usage bias. Further insight on the mutational equilibrium within Orthocoronavirinae revealed that most coronavirus genomes are close to their neutral equilibrium, the exception being the three recently infecting human coronaviruses, which lie further away from the mutational equilibrium than their endemic human coronavirus counterparts. Finally, our results suggest that, while replicating in humans, SARS-CoV-2 is slowly becoming AU-richer, likely until attaining a new mutational equilibrium.

A Code Within a Code: How Codons Fine-Tune Protein Folding in the Cell

The genetic code sets the correspondence between the sequence of a given nucleotide triplet in an mRNA molecule, called a codon, and the amino acid that is added to the growing polypeptide chain during protein synthesis. With four bases (A, G, U, and C), there are 64 possible triplet codons: 61 sense codons (encoding amino acids) and 3 nonsense codons (so-called, stop codons that define termination of translation). In most organisms, there are 20 common/standard amino acids used in protein synthesis; thus, the genetic code is redundant with most amino acids (with the exception of Met and Trp) are being encoded by more than one (synonymous) codon. Synonymous codons were initially presumed to have entirely equivalent functions, however, the finding that synonymous codons are not present at equal frequencies in mRNA suggested that the specific codon choice might have functional implications beyond coding for amino acid. Observation of nonequivalent use of codons in mRNAs implied a possibility of the existence of auxiliary information in the genetic code. Indeed, it has been found that genetic code contains several layers of such additional information and that synonymous codons are strategically placed within mRNAs to ensure a particular translation kinetics facilitating and fine-tuning co-translational protein folding in the cell via step-wise/sequential structuring of distinct regions of the polypeptide chain emerging from the ribosome at different points in time. This review summarizes key findings in the field that have identified the role of synonymous codons and their usage in protein folding in the cell.

Distinct signatures of codon and codon pair usage in 32 primary tumor types in the novel database CancerCoCoPUTs for cancer-specific codon usage

BACKGROUND

Gene expression is highly variable across tissues of multi-cellular organisms, influencing the codon usage of the tissue-specific transcriptome. Cancer disrupts the gene expression pattern of healthy tissue resulting in altered codon usage preferences. The topic of codon usage changes as they relate to codon demand, and tRNA supply in cancer is of growing interest.

METHODS

We analyzed transcriptome-weighted codon and codon pair usage based on The Cancer Genome Atlas (TCGA) RNA-seq data from 6427 solid tumor samples and 632 normal tissue samples. This dataset represents 32 cancer types affecting 11 distinct tissues. Our analysis focused on tissues that give rise to multiple solid tumor types and cancer types that are present in multiple tissues.

RESULTS

We identified distinct patterns of synonymous codon usage changes for different cancer types affecting the same tissue. For example, a substantial increase in GGT-glycine was observed in invasive ductal carcinoma (IDC), invasive lobular carcinoma (ILC), and mixed invasive ductal and lobular carcinoma (IDLC) of the breast. Change in synonymous codon preference favoring GGT correlated with change in synonymous codon preference against GGC in IDC and IDLC, but not in ILC. Furthermore, we examined the codon usage changes between paired healthy/tumor tissue from the same patient. Using clinical data from TCGA, we conducted a survival analysis of patients based on the degree of change between healthy and tumor-specific codon usage, revealing an association between larger changes and increased mortality. We have also created a database that contains cancer-specific codon and codon pair usage data for cancer types derived from TCGA, which represents a comprehensive tool for codon-usage-oriented cancer research.

CONCLUSIONS

Based on data from TCGA, we have highlighted tumor type-specific signatures of codon and codon pair usage. Paired data revealed variable changes to codon usage patterns, which must be considered when designing personalized cancer treatments. The associated database, CancerCoCoPUTs, represents a comprehensive resource for codon and codon pair usage in cancer and is available at https://dnahive.fda.gov/review/cancercocoputs/ . These findings are important to understand the relationship between tRNA supply and codon demand in cancer states and could help guide the development of new cancer therapeutics.

Particle creation and annihilation in a dynamically disordered totally asymmetric simple exclusion process

We study a single-channel dynamically disordered totally asymmetric simple exclusion process with bulk particle attachment and detachment. The continuum mean-field equations are derived and solved numerically to obtain steady-state phase diagrams and density profiles. The effects of various parameters, namely particle attachment rate, defect binding and unbinding rates, and binding constant, have been investigated. An increase in the attachment rate of particles reduces the number of steady-state phases, whereas a variation in defect binding and unbinding rates shifts the phase boundaries. One of the important consequences of introducing particle nonconserving dynamics is the appearance of shock in the steady state. The shock dynamics have been thoroughly examined and the defect strength is found to have a significant effect on the shock position. The mean-field solutions are validated using extensive Monte Carlo simulations.

Codon usage bias and environmental adaptation in microbial organisms

In each genome, synonymous codons are used with different frequencies; this general phenomenon is known as codon usage bias. It has been previously recognised that codon usage bias could affect the cellular fitness and might be associated with the ecology of microbial organisms. In this exploratory study, we investigated the relationship between codon usage bias, lifestyles (thermophiles vs. mesophiles; pathogenic vs. non-pathogenic; halophilic vs. non-halophilic; aerobic vs. anaerobic and facultative) and habitats (aquatic, terrestrial, host-associated, specialised, multiple) of 615 microbial organisms (544 bacteria and 71 archaea). Principal component analysis revealed that species with given phenotypic traits and living in similar environmental conditions have similar codon preferences, as represented by the relative synonymous codon usage (RSCU) index, and similar spectra of tRNA availability, as gauged by the tRNA gene copy number (tGCN). Moreover, by measuring the average tRNA adaptation index (tAI) for each genome, an index that can be associated with translational efficiency, we observed that organisms able to live in multiple habitats, including facultative organisms, mesophiles and pathogenic bacteria, are characterised by a reduced translational efficiency, consistently with their need to adapt to different environments. Our results show that synonymous codon choices might be under strong translational selection, which modulates the choice of the codons to differently match tRNA availability, depending on the organism's lifestyle needs. To our knowledge, this is the first large-scale study that examines the role of codon bias and translational efficiency in the adaptation of microbial organisms to the environment in which they live.

Downregulating of hemB via synthetic antisense RNAs for improving 5-aminolevulinic acid production in Escherichia coli

Aminolevulinic acid (ALA), a type of natural non-protein amino acid, is a key precursor for the biosynthesis of heme, and it has been broadly applied in medicine, agriculture. Several strategies have been applied to enhance ALA synthesis in bacteria. In the present study, we employed synthetic antisense RNAs (asRNAs) of hemB (encodes ALA dehydratase) to weaken metabolic flux of ALA to porphobilinogen (PBG), and investigated their effect on ALA accumulation. For this purpose, we designed and constructed vectors pET28a-hemA-asRNA and pRSFDuet-hemA-asRNA to simultaneously express 5-ALA synthase (ALAS, encoded by hemA) and PTasRNAs (2 inverted repeat DNA sequences sandwiched with the antisense sequence of hemB), selecting the region ranging from - 57 nt upstream to + 139 nt downstream of the start codon of hemB as a target. The qRT-PCR analysis showed that the mRNA levels of hemB were decreased above 50% of the control levels, suggesting that the anti-hemB asRNA was functioning appropriately. ALA accumulation in the hemB weakened strains were 17.6% higher than that obtained using the control strains while accumulating less PBG. These results indicated that asRNAs can be used as a tool for regulating ALA accumulation in E. coli.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02733-8.

Context-dependent and -independent selection on synonymous mutations revealed by 1,135 genomes of Arabidopsis thaliana

BACKGROUND

Synonymous mutations do not alter the amino acids and therefore are regarded as neutral for a long time. However, they do change the tRNA adaptation index (tAI) of a particular codon (independent of its context), affecting the tRNA availability during translation. They could also change the isoaccepting relationship with its neighboring synonymous codons in particular context, which again affects the local translation process. Evidence of selection pressure on synonymous mutations has emerged.

RESULTS

The proposed selection patterns on synonymous mutations are never formally and systematically tested in plant species. We fully take advantage of the SNP data from 1,135 A. thaliana lines, and found that the synonymous mutations that increase tAI or the isoaccepting mutations in isoaccepting codon context tend to have higher derived allele frequencies (DAF) compared to other synonymous mutations of the opposite effects.

CONCLUSIONS

Synonymous mutations are not strictly neutral. The synonymous mutations that increase tAI or the isoaccepting mutations in isoaccepting codon context are likely to be positively selected. We propose the concept of context-dependent and -independent selection on synonymous mutations. These concepts broaden our knowledge of the functional consequences of synonymous mutations, and should be appealing to phytologists and evolutionary biologists.

Comprehensive annotation and characterization of planarian tRNA and tRNA-derived fragments (tRFs)

tRNA-derived fragments (tRFs) have recently gained a lot of scientific interest due to their diverse regulatory roles in several cellular processes. However, their function in dynamic biological processes such as development and regeneration remains unexplored. Here, we show that tRFs are dynamically expressed during planarian regeneration, suggesting a possible role for these small RNAs in the regulation of regeneration. In order to characterize planarian tRFs, we first annotated 457 tRNAs in S. mediterranea combining two tRNA prediction algorithms. Annotation of tRNAs facilitated the identification of three main species of tRFs in planarians-the shorter tRF-5s and itRFs, and the abundantly expressed 5'-tsRNAs. Spatial profiling of tRFs in sequential transverse sections of planarians revealed diverse expression patterns of these small RNAs, including those that are enriched in the head and pharyngeal regions. Expression analysis of these tRF species revealed dynamic expression of these small RNAs over the course of regeneration suggesting an important role in planarian anterior and posterior regeneration. Finally, we show that 5'-tsRNA in planaria interact with all three SMEDWI proteins and an involvement of AGO1 in the processing of itRFs. In summary, our findings implicate a novel role for tRFs in planarian regeneration, highlighting their importance in regulating complex systemic processes. Our study adds to the catalog of posttranscriptional regulatory systems in planaria, providing valuable insights on the biogenesis and the function of tRFs in neoblasts and planarian regeneration.

Impact of Synonymous Genome Recoding on the HIV Life Cycle

Synonymous mutations within protein coding regions introduce changes in DNA or messenger (m) RNA, without mutating the encoded proteins. Synonymous recoding of virus genomes has facilitated the identification of previously unknown virus biological features. Moreover, large-scale synonymous recoding of the genome of human immunodeficiency virus type 1 (HIV-1) has elucidated new antiviral mechanisms within the innate immune response, and has improved our knowledge of new functional virus genome structures, the relevance of codon usage for the temporal regulation of viral gene expression, and HIV-1 mutational robustness and adaptability. Continuous improvements in our understanding of the impacts of synonymous substitutions on virus phenotype - coupled with the decreased cost of chemically synthesizing DNA and improved methods for assembling DNA fragments - have enhanced our ability to identify potential HIV-1 and host factors and other aspects involved in the infection process. In this review, we address how silent mutagenesis impacts HIV-1 phenotype and replication capacity. We also discuss the general potential of synonymous recoding of the HIV-1 genome to elucidate unknown aspects of the virus life cycle, and to identify new therapeutic targets.

The effects of codon bias and optimality on mRNA and protein regulation

The central dogma of molecular biology entails that genetic information is transferred from nucleic acid to proteins. Notwithstanding retro-transcribing genetic elements, DNA is transcribed to RNA which in turn is translated into proteins. Recent advancements have shown that each stage is regulated to control protein abundances for a variety of essential physiological processes. In this regard, mRNA regulation is essential in fine-tuning or calibrating protein abundances. In this review, we would like to discuss one of several mRNA-intrinsic features of mRNA regulation that has been gaining traction of recent-codon bias and optimality. Specifically, we address the effects of codon bias with regard to codon optimality in several biological processes centred on translation, such as mRNA stability and protein folding among others. Finally, we examine how different organisms or cell types, through this system, are able to coordinate physiological pathways to respond to a variety of stress or growth conditions.

Molecular motor traffic with a slow binding site

We discuss how the presence of a slow binding site in molecular motor traffic gives rise to defect-induced "traffic jams" that have properties different from those of the well-studied boundary-induced jams that originate from an imbalance between initiation and termination. To this end we analyze in detail the stationary distribution of a lattice gas model for traffic of molecular motors with a defect. In particular, we obtain analytically the exact spatial distribution of motors, the probability distribution of the random position of the molecular traffic jam and we report unexpected spatial anticorrelations between local molecular motor densities near the defect.

Cost-efficiency tradeoff is optimized in various cancer types revealed by genome-wide analysis

The tradeoff between cost and efficiency is omnipresent in organisms. Specifically, how the evolutionary force shapes the tradeoff between biosynthetic cost and translation efficiency remains unclear. In the cancer community, whether the adjustment of cost-efficiency tradeoff acts as a strategy to facilitate tumor proliferation and contributes to oncogenesis is uninvestigated. To address this issue, we retrieved the gene expression profile in various cancer types and the matched normal samples from The Cancer Genome Atlas (TCGA). We found that the highly expressed genes in cancers generally have higher tAI/nitro ratios than those in normal samples. This is possibly caused by the higher tAI/nitro ratios observed in oncogenes than tumor suppressor genes (TSG). Furthermore, in the cancer samples, derived mutations in oncogenes usually lead to higher tAI/nitro ratios, while those mutations in TSG lead to lower tAI/nitro. For a special case of kidney cancer, we investigated several crucial genes in tumor samples versus normal samples, and discovered that the changes in tAI/nitro ratios are correlated with the changes in translation level. Our study for the first time revealed the optimization of cost-efficiency tradeoff in cancers. The cost-efficiency dilemma is optimized by the tumor cells, and is possibly beneficial for the translation and production of oncogenes, and eventually contributes to proliferation and oncogenesis. Our findings could provide novel perspectives in depicting the cancer genomes and might help unravel the cancer evolution.

Synonymous but Not Silent: The Codon Usage Code for Gene Expression and Protein Folding

Codon usage bias, the preference for certain synonymous codons, is found in all genomes. Although synonymous mutations were previously thought to be silent, a large body of evidence has demonstrated that codon usage can play major roles in determining gene expression levels and protein structures. Codon usage influences translation elongation speed and regulates translation efficiency and accuracy. Adaptation of codon usage to tRNA expression determines the proteome landscape. In addition, codon usage biases result in nonuniform ribosome decoding rates on mRNAs, which in turn influence the cotranslational protein folding process that is critical for protein function in diverse biological processes. Conserved genome-wide correlations have also been found between codon usage and protein structures. Furthermore, codon usage is a major determinant of mRNA levels through translation-dependent effects on mRNA decay and translation-independent effects on transcriptional and posttranscriptional processes. Here, we discuss the multifaceted roles and mechanisms of codon usage in different gene regulatory processes.

Expression of transgenes enriched in rare codons is enhanced by the MAPK pathway

The ability to translate three nucleotide sequences, or codons, into amino acids to form proteins is conserved across all organisms. All but two amino acids have multiple codons, and the frequency that such synonymous codons occur in genomes ranges from rare to common. Transcripts enriched in rare codons are typically associated with poor translation, but in certain settings can be robustly expressed, suggestive of codon-dependent regulation. Given this, we screened a gain-of-function library for human genes that increase the expression of a GFP^rare reporter encoded by rare codons. This screen identified multiple components of the mitogen activated protein kinase (MAPK) pathway enhancing GFP^rare expression. This effect was reversed with inhibitors of this pathway and confirmed to be both codon-dependent and occur with ectopic transcripts naturally coded with rare codons. Finally, this effect was associated, at least in part, with enhanced translation. We thus identify a potential regulatory module that takes advantage of the redundancy in the genetic code to modulate protein expression.

The birth of a bacterial tRNA gene by large-scale, tandem duplication events

Organisms differ in the types and numbers of tRNA genes that they carry. While the evolutionary mechanisms behind tRNA gene set evolution have been investigated theoretically and computationally, direct observations of tRNA gene set evolution remain rare. Here, we report the evolution of a tRNA gene set in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The growth defect caused by deleting the single-copy tRNA gene, serCGA, is rapidly compensated by large-scale (45–290 kb) duplications in the chromosome. Each duplication encompasses a second, compensatory tRNA gene (serTGA) and is associated with a rise in tRNA-Ser(UGA) in the mature tRNA pool. We postulate that tRNA-Ser(CGA) elimination increases the translational demand for tRNA-Ser(UGA), a pressure relieved by increasing serTGA copy number. This work demonstrates that tRNA gene sets can evolve through duplication of existing tRNA genes, a phenomenon that may contribute to the presence of multiple, identical tRNA gene copies within genomes.

A code within the genetic code: codon usage regulates co-translational protein folding

The genetic code is degenerate, and most amino acids are encoded by two to six synonymous codons. Codon usage bias, the preference for certain synonymous codons, is a universal feature of all genomes examined. Synonymous codon mutations were previously thought to be silent; however, a growing body evidence now shows that codon usage regulates protein structure and gene expression through effects on co-translational protein folding, translation efficiency and accuracy, mRNA stability, and transcription. Codon usage regulates the speed of translation elongation, resulting in non-uniform ribosome decoding rates on mRNAs during translation that is adapted to co-translational protein folding process. Biochemical and genetic evidence demonstrate that codon usage plays an important role in regulating protein folding and function in both prokaryotic and eukaryotic organisms. Certain protein structural types are more sensitive than others to the effects of codon usage on protein folding, and predicted intrinsically disordered domains are more prone to misfolding caused by codon usage changes than other domain types. Bioinformatic analyses revealed that gene codon usage correlates with different protein structures in diverse organisms, indicating the existence of a codon usage code for co-translational protein folding. This review focuses on recent literature on the role and mechanism of codon usage in regulating translation kinetics and co-translational protein folding.

Viral adaption of staphylococcal phage: A genome-based analysis of the selective preference based on codon usage Bias

Given the high therapeutic value of the staphylococcal phage, the genome co-evolution of the phage and the host has gained great attention. Though the genome-wide AT richness in staphylococcal phages has been well-studied with nucleotide usage bias, here we proved that host factor, lifestyle and taxonomy are also important factors in understanding the phage nucleotide usages bias using information entropy formula. Such correlation is especially prominent when it comes to the synonymous codon usages of staphylococcal phages, despite the overall scattered codon usage pattern represented by principal component analysis. This strong relationship is explained by nucleotide skew which testified that the usage biases of nucleotide at different codon positions are acting on synonymous codons. Therefore, our study reveals a hidden relationship of genome evolution with host limitation and phagic phenotype, providing new insight into phage genome evolution at genetic level.

The combined influence of codon composition and tRNA copy number regulates translational efficiency by influencing synonymous nucleotide substitution

Codon usage bias is an important genomic phenomenon, where highly expressed genes use optimal codons for smoother translation with high yield, facilitated by the cognate tRNAs. Here, we presented the tRNA co-adaptation index (co-AI) by correlating tRNA gene copy number and codon composition in Saccharomyces cerevisiae. We observed that this co-AI is positively correlated with protein abundance and translation rate. Considering nucleotide substitutions, co-AI influences synonymous substitutions more than gene expression and protein abundance, the most important determinants of evolutionary rate. Co-AI correlates positively with mRNA secondary structure stability and mRNA half-life, which may lead to protein accumulation under high co-AI. However, the highly expressed proteins encoded by high co-AI genes are assisted by molecular chaperones to attain their proper functional conformation and prevent accumulation.

Tracing Brucella evolutionary dynamics in expanding host ranges through nucleotide, codon and amino acid usages in genomes

The host range of Brucella organisms has expanded from terrestrial and marine mammals to fish and amphibians. The high homology genomes of different Brucella organisms promote us to investigate evolutionary patterns for nucleotide, codon and amino acid usage patterns at gene levels among Brucella species. Although the similar patterns for nucleotide and synonymous codon usages exist in gene population, GC composition at the first codon position has significant correlations to that of the second and third codon positions, respectively, suggesting that nucleotide usages surrounding one codon influence synonymous codon usage patterns. Evolutionary patterns represented by synonymous codon and amino acid usages reflect host factor impacting Brucella speciation. As for genetic variations of important virulent factors involved with different biological functions, genes encoding lipoplysaccharides (LPSs) display more distinctive codon adaptation to Brucella than those of the BvrR/BvrS system and type IV secretion system. By Bayesian analysis, the polygenetic constructions for these genes of virulent factors shared by Brucella species display the purifying/positive selections and partially host factor in mediating genetic variations of these genes. The systemic analyses for nucleotide, synonymous codon and amino acid usages at gene level and genetic variations of important virulent factor genes display that host limitation influences either genetic characterizations at gene level or a particular gene involved in virulent factors of Brucella.Communicated by Ramaswamy H. Sarma.

Codon usage similarity between viral and some host genes suggests a codon-specific translational regulation

The codon usage pattern is a specific characteristic of each species; however, the codon usage of all of the genes in a genome is not uniform. Intriguingly, most viruses have codon usage patterns that are vastly different from the optimal codon usage of their hosts. How viral genes with different codon usage patterns are efficiently expressed during a viral infection is unclear. An analysis of the similarity between viral codon usage and the codon usage of the individual genes of a host genome has never been performed. In this study, we demonstrated that the codon usage of human RNA viruses is similar to that of some human genes, especially those involved in the cell cycle. This finding was substantiated by its concordance with previous reports of an upregulation at the protein level of some of these biological processes. It therefore suggests that some suboptimal viral codon usage patterns may actually be compatible with cellular translational machineries in infected conditions.

Synonymous genome recoding: a tool to explore microbial biology and new therapeutic strategies

Synthetic genome recoding is a new means of generating designed organisms with altered phenotypes. Synonymous mutations introduced into the protein coding region tolerate modifications in DNA or mRNA without modifying the encoded proteins. Synonymous genome-wide recoding has allowed the synthetic generation of different small-genome viruses with modified phenotypes and biological properties. Recently, a decreased cost of chemically synthesizing DNA and improved methods for assembling DNA fragments (e.g. lambda red recombination and CRISPR-based editing) have enabled the construction of an Escherichia coli variant with a 4-Mb synthetic synonymously recoded genome with a reduced number of sense codons (n = 59) encoding the 20 canonical amino acids. Synonymous genome recoding is increasing our knowledge of microbial interactions with innate immune responses, identifying functional genome structures, and strategically ameliorating cis-inhibitory signaling sequences related to splicing, replication (in eukaryotes), and complex microbe functions, unraveling the relevance of codon usage for the temporal regulation of gene expression and the microbe mutant spectrum and adaptability. New biotechnological and therapeutic applications of this methodology can easily be envisaged. In this review, we discuss how synonymous genome recoding may impact our knowledge of microbial biology and the development of new and better therapeutic methodologies.

Widespread non-modular overlapping codes in the coding regions

Messenger RNAs (mRNAs) consist of a coding region (open reading frame (ORF)) and two untranslated regions (UTRs), 5'UTR and 3'UTR. Ribosomes travel along the coding region, translating nucleotide triplets (called codons) to a chain of amino acids. The coding region was long believed to mainly encode the amino acid content of proteins, whereas regulatory signals reside in the UTRs and in other genomic regions. However, in recent years we have learned that the ORF is expansively populated with various regulatory signals, or codes, which are related to all gene expression steps and additional intracellular aspects. In this paper, we review the current knowledge related to overlapping codes inside the coding regions, such as the influence of synonymous codon usage on translation speed (and, in turn, the effect of translation speed on protein folding), ribosomal frameshifting, mRNA stability, methylation, splicing, transcription and more. All these codes come together and overlap in the ORF sequence, ensuring production of the right protein at the right time.

Expression and purification of the full-length N-acetylgalactosaminyltransferase and galactosyltransferase from Campylobacter jejuni in Escherichia coli

The successful enzymatic synthesis of various ganglioside-related oligosaccharides requires many available glycan-processing enzymes. However, the number of available glycan-processing enzymes remains limited. In this study, the full-length CgtA43456 (β-(1→4)-N-acetylgalactosaminyltransferase) and CgtB11168 (β-(1→3)-galactosyltransferase) were successfully produced from Escherichia coli through the optimization of E. coli-preferable codon usage, selection of E. coli strain, and use of the molecular chaperone GroEL-GroES (GroEL/ES). The CgtA43456 enzyme was produced as a soluble form in E. coli C41(DE3) co-expressed with codon-optimized CgtA43456 and GroEL/ES. However, soluble CgtB11168 was well expressed in E. coli C41(DE3) with only the codon-optimized CgtB11168. Rather, when co-expressed with GroEL/ES, total production of CgtB11168 was reduced. Using immobilized-metal affinity chromatography, the CgtA43456 and CgtB11168 proteins were obtained with approximately 75-78 % purity. The purified CgtA43456 showed a specific activity of 21 mU/mg using UDP-N-acetylgalactosamine and GM3 trisaccharide as donor and acceptor, respectively. The purified CgtB11168 catalyzed the transfer of galactose from UDP-Gal to GM2 tetrasaccharide with a specific activity of 16 mU/mg. We propose that they could be used as catalysts for enzymatic synthesis of GM1 ganglioside-related oligosaccharides.