Synthetic biology approach revealed enhancement in haeme oxygenase-1 gene expression by codon pair optimization while reduction by codon deoptimization

Haem oxygenase-1 (HO-1) is a ubiquitously expressed gene involved in cellular homoeostasis, and its imbalance in expression results in various disorders. To alleviate such disorders, HO-1 gene expression needs to be modulated. Codon usage bias results from evolutionary forces acting on any nucleotide sequence and determines the gene expression. Like codon usage bias, codon pair bias also exists, playing a role in gene expression. In the present study, HO-1 gene was recoded by manipulating codon and codon pair bias, and four such constructs were made through codon/codon pair deoptimization and codon/codon pair optimization to reduce and enhance the HO-1 gene expression. Codon usage analysis was done for these constructs for four tissues brain, heart, pancreas and liver. Based on codon usage in different tissues, gene expression of these tissues was determined in terms of the codon adaptation index. Based on the codon adaptation index, minimum free energy, and translation efficiency, constructs were evaluated for enhanced or decreased HO-1 expression. The analysis revealed that for enhancing gene expression, codon pair optimization, while for reducing gene expression, codon deoptimization is efficacious. The recoded constructs developed in the study could be used in gene therapy regimens to cure HO-1 over or underexpression-associated disorders.

A smart RBS library and its prediction model for robust and accurate fine-tuning of gene expression in Bacillus species

Bacillus species, such as Bacillus subtilis and Bacillus licheniformis, are important industrial bacteria. However, there is a lack of standardized and predictable genetic tools for convenient and reproducible assembly of genetic modules in Bacillus species to realize their full potential. In this study, we constructed a Ribosome Binding Site (RBS) library in B. licheniformis, which provides incremental regulation of expression levels over a 104-fold range. Additionally, we developed a model to quantify the resulting translation rates. We successfully demonstrated the robust expression of various target genes using the RBS library and showed that the model accurately predicts the translation rates of arbitrary coding genes. Importantly, we also extended the use of the RBS library and prediction model to B. subtilis, B. thuringiensis, and B. amyloliquefacie. The versatility of the RBS library and its prediction model enables quantification of biological behavior, facilitating reliable forward engineering of gene expression.

Polysome Profiling in Adult Mouse Testes

Polysome profiling is widely used to isolate and analyze polysome fractions, which consist of actively translating mRNAs and ribosomes. Compared to ribosome profiling and translating ribosome affinity purification, polysome profiling is simpler and less time consuming in sample preparation and library constructions. Spermiogenesis, i.e., the post-meiotic phase of male germ cell development, is a highly coordinated developmental process in which transcription and translation are decoupled because of nuclear condensation, resulting in translation regulation as the major mode for the regulation of gene expression in post-meiotic spermatids. To understand the translation regulation during spermiogenesis, an overview of translational state of spermiogenic mRNAs is required. Here, we describe a protocol to identify translating mRNAs using polysome profiling. Briefly, mouse testes are gently homogenized to release polysomes containing translating mRNAs, following polysome-bound mRNAs isolated by sucrose density gradient purification and characterized by RNA-seq. This protocol allows to quickly isolate translating mRNAs from testes and analyze the discrepancy of translational efficiency in mouse testes from different mouse lines. Key features Quickly obtain polysome RNAs from testes. Omit RNase digestion and RNA recovery from gel. High efficiency and robustness compared to ribo-seq. Graphical overview Schematic illustrating the experimental design for polysome profiling in mouse testes. Mouse testes are homogenized and lysed in Sample preparation, and polysome RNAs are enriched by sucrose gradient centrifugation and used to calculate translation efficiency in Sample analysis.

25KDa branched polyethylenimine increases interferon-γ production in natural killer cells via improving translation efficiency

BACKGROUND

Ex vivo cultivation is a promising strategy for increasing the number of NK cells and enhancing their antitumor activity prior to clinical application. Recent studies show that stimulation with 25KDa branched polyethylenimine (25KbPEI) generates NK cells with enhanced antitumor activity. To better understand how 25KbPEI primes NK cells, we explored the mechanism underlying increase in production of IFN-γ.

METHODS

Chemical priming was performed on NK-92MI cells by incubating them with 5 μg/ml of 25KbPEI. The production of IFN-γ was evaluated by RT-qPCR, ELISA, and Flow cytometry. By evaluating the effect of pharmacological inhibition of ERK/mTOR-eIF4E signaling pathways on IFN-γ translation, the function of these signaling pathways in IFN-γ translation was examined. To comprehend the level of 25KbPEI activity on immune-related components in NK cells, RNA sequencing and proteomics analyses were conducted.

RESULTS

25KbPEI enhances the production of IFN-γ by NK cells without transcriptional activation. Activation of ERK and mTOR signaling pathways was found to be associated with 25KbPEI-mediated calcium influx in NK cells. The activation of ERK/mTOR signaling was linked to the phosphorylation of 4E-BP1, which resulted in the activation of translation initiation complex and subsequent IFN-γ translation. Analysis of RNA sequencing and proteomics data revealed that the activity of 25KbPEI to improve translation efficiency in NK cells could be extended to additional immune-related molecules.

CONCLUSIONS

This study provides substantial insight into the process by which 25KbPEI primes NK cells. Our data demonstrated that the 25KbPEI mediated activation of ERK/mTOR signaling and subsequent stimulation of eIF4E is the primary mechanism by which the chemical stimulates translation of IFN-γ in NK cells. Video abstract.

Translational activation of ribosome-related genes at initial photoreception is dependent on signals derived from both the nucleus and the chloroplasts in Arabidopsis thaliana

Light is one of the indispensable elements that plants need in order to grow and develop. In particular, it is essential for inducing morphogenesis, such as suppression of hypocotyl elongation and cotyledon expansion, that plants undergo when they first emerge after germination. However, there is a lack of knowledge about the gene expression and, in particular, the translational levels that induce a response upon light exposure. We have investigated the translational expression of nuclear genes in Arabidopsis thaliana seedlings germinated in the dark and then exposed to blue monochromatic light. In this study, ribosome profiling analysis was performed in the blue-light-receptor mutant cry1cry2 and the light-signaling mutant hy5 to understand which signaling pathways are responsible for the changes in gene expression at the translational level after blue-light exposure. The analysis showed that the expression of certain chloroplast- and ribosome-related genes was up-regulated at the translational level in the wild type. However, in both mutants the translational up-regulation of ribosome-related genes was apparently compromised. This suggests that light signaling through photoreceptors and the HY5 transcription factor are responsible for translation of ribosome-related genes. To further understand the effect of photoreception by chloroplasts on nuclear gene expression, chloroplast function was inhibited by adding a photosynthesis inhibitor, 3-(3,4-dichlorophenyl)-1,1-dimethylurea (DCMU), and a carotenoid synthesis inhibitor, norflurazon. The results show that inhibition of chloroplast function did not lead to an increase in the expression of ribosome-related genes at the translational level. These results suggest that signals from both the nucleus and chloroplasts are required to activate translation of ribosome-related genes during blue-light reception.

Click-iT trinucleotide cap analog: Synthesis, mRNA translation, and detection

The first example of the synthesis of a new trinucleotide cap analog containing propargyl group such as m^{7,3'-O-propargyl}G(5')PPP(5')A_mpG is reported. The effect of the propargyl group in trinucleotide analog with a standard trinucleotide cap analog (GAG), m⁷G(5')ppp(5')A_mpG was evaluated with respect to their capping efficiency, in vitro T7 RNA polymerase transcription efficiency, and translation activity using cultured A549 lung carcinoma epithelial cells. The new propargyl cap analog is a substrate for T7 RNA polymerase. Notably, the mRNA capped with the propargyl cap is translated ∼ 1.3 times more efficiently than the mRNA capped with the GAG cap. The most characteristic feature of the new propargyl cap analog is that the presence of the propargyl group allows further modification of the mRNA by chemical ligation of an azide-containing fluorescent-labeled substrate to the mRNA via click chemistry.

Conversion of Ribosome-Protected mRNA to DNA and Deep Sequencing for Ribosome Profiling in Haloferax volcanii

The translation of messenger RNA (mRNA) into protein is an essential process for all forms of life. The ability to monitor this process in a quantitative way by ribosome profiling-based approaches has revolutionized our ability to monitor protein synthesis in vivo and to explore and model complex cellular processes. Ribosome profiling is a high-throughput technique that globally analyzes the full set of ribosomes engaged in translation, providing insights into important aspects of the mechanism of protein synthesis and its regulation. This protocol covers the construction of a ribosome profiling library from culture harvesting, footprint isolation via ultracentrifugation, gel-based size fractionation, and footprint sequencing for a model halophilic archaeon, Haloferax volcanii. This approach has revealed the first global view of translation in the archaea.

N6-methyladenosine RNA modification regulates strawberry fruit ripening in an ABA-dependent manner

BACKGROUND

Epigenetic mark such as DNA methylation plays pivotal roles in regulating ripening of both climacteric and non-climacteric fruits. However, it remains unclear whether mRNA m6A methylation, which has been shown to regulate ripening of the tomato, a typical climacteric fruit, is functionally conserved for ripening control among different types of fruits.

RESULTS

Here we show that m6A methylation displays a dramatic change at ripening onset of strawberry, a classical non-climacteric fruit. The m6A modification in coding sequence (CDS) regions appears to be ripening-specific and tends to stabilize the mRNAs, whereas m6A around the stop codons and within the 3' untranslated regions is generally negatively correlated with the abundance of associated mRNAs. We identified thousands of transcripts with m6A hypermethylation in the CDS regions, including those of NCED5, ABAR, and AREB1 in the abscisic acid (ABA) biosynthesis and signaling pathway. We demonstrate that the methyltransferases MTA and MTB are indispensable for normal ripening of strawberry fruit, and MTA-mediated m6A modification promotes mRNA stability of NCED5 and AREB1, while facilitating translation of ABAR.

CONCLUSION

Our findings uncover that m6A methylation regulates ripening of the non-climacteric strawberry fruit by targeting the ABA pathway, which is distinct from that in the climacteric tomato fruit.

Comparative analysis of the relationship between translation efficiency and sequence features of endogenous proteins in multiple organisms

The translation efficiency of protein genes is known to be affected by sequence features. Previous studies have found that various sequence features based on codon usage and mRNA secondary structure contribute to translation efficiency. However, most studies have focused on a specific organism, usually a model organism such as Escherichia coli or Saccharomyces cerevisiae. Here, we investigate whether the relationship between translation efficiency and sequence features is conserved among multiple organisms using publicly available ribosome profiling data and RNA-Seq data. We analyze nine organisms from various taxa: Staphylococcus aureus, five species of Streptomyces, two strains of E. coli, and S. cerevisiae. We reveal that the relationship between translation efficiency and sequence features differs across organisms, partly reflecting their taxonomy. The codon adaptation index shows high correlation in all analyzed organisms. Our study provides an insight into the diversity and commonality of sequence determinants of protein expression in these organisms.

Reduced translation efficiency due to novel splicing variants in 5' untranslated region and identification of novel cis-regulatory elements in canine and human BRCA2

Background

Breast cancer 2, early onset (BRCA2) is a tumor suppressor gene. The protein encoded by this gene plays an important role in homologous recombination (HR)-mediated DNA repair. Deleterious mutations in BRCA2 and downregulation of its expression have been associated with tumorigenesis in dogs and humans. Thus, regulation of BRCA2 expression level is important for maintaining homeostasis in homologous recombination.

Results

In this study, the mechanisms that regulate the expression of BRCA2 were proposed. Novel splicing variants were identified in the 5′ untranslated region (UTR) of canine and human BRCA2 in canine testis, canine ovary, and canine and human cultured cell lines. In cultured cells, the ratio of BRCA2 splicing variants at the 5′ UTR was altered by serum starvation. These novel splicing variants, excluding one of the canine splicing variants, were found to reduce the translational efficiency. Additionally, the DNA sequence in human BRCA2 intron 1 harbored novel cis-regulatory elements. Three silencer and two enhancer cis-regulatory elements were identified in human BRCA2 intron 1.

Conclusions

This study demonstrates that BRCA2 expression level is regulated via 5′ UTR splicing variants and that the BRCA2 intron 1 region harbors cis-regulatory elements.

Ribosome profiling reveals the effects of nitrogen application translational regulation of yield recovery after abrupt drought-flood alternation in rice

Abrupt drought-flood alternation is a frequent meteorological disaster during the summer in Southern China. The study of physiological and translation mechanisms of rice yield recovery after abrupt drought-flood alternation has great potential benefits in field production. Our results showed that yield recovery upon nitrogen (N) application after abrupt drought-flood alternation was due to the increase in effective panicle numbers per plant. The N application resulted in the regulation of physiological and biochemical as well as growth development processes, which led to a rapid growth recovery effect after abrupt drought-flood alternation stress in rice. Using ribosome profiling combined with RNA sequencing (RNA-seq) technology, the interactions between transcription and translation for N application after abrupt drought-flood alternation were analyzed. It was found that a small proportion of response genes were shared at the transcriptional and translational levels, that is, 14% of the expressed genes were upregulated and 6.6% downregulated. Further analysis revealed that the translation efficiency (TE) of the genes was influenced by their sequence characteristics, including their GC content, coding sequence length and normalized minimal free energy. Compared with the number of untranslated upstream open reading frames (uORFs), the increased number of translated uORFs promoted the improvement of TE. The TE of the uORFs for N application was lower than the control without N application after abrupt drought-flood alternation. This study characterizes the translational regulatory pattern in response to N application after abrupt drought-flood alternation stress.

Determinants of translation efficiency in the evolutionarily-divergent protist Trichomonas vaginalis

BACKGROUND

Trichomonas vaginalis, the causative agent of a prevalent urogenital infection in humans, is an evolutionarily divergent protozoan. Protein-coding genes in T. vaginalis are largely controlled by two core promoter elements, producing mRNAs with short 5' UTRs. The specific mechanisms adopted by T. vaginalis to fine-tune the translation efficiency (TE) of mRNAs remain largely unknown.

RESULTS

Using both computational and experimental approaches, this study investigated two key factors influencing TE in T. vaginalis: codon usage and mRNA secondary structure. Statistical dependence between TE and codon adaptation index (CAI) highlighted the impact of codon usage on mRNA translation in T. vaginalis. A genome-wide interrogation revealed that low structural complexity at the 5' end of mRNA followed closely by a highly structured downstream region correlates with TE variation in this organism. To validate these findings, a synthetic library of 15 synonymous iLOV genes was created, representing five mRNA folding profiles and three codon usage profiles. Fluorescence signals produced by the expression of these synonymous iLOV genes in T. vaginalis were consistent with and validated our in silico predictions.

CONCLUSIONS

This study demonstrates the role of codon usage bias and mRNA secondary structure in TE of T. vaginalis mRNAs, contributing to a better understanding of the factors that influence, and possibly regulate, gene expression in this human pathogen.

Ribosome Profiling in Trypanosomatids

Ribosomes are the machinery responsible for reading mRNAs and translating them into proteins. The ribosome profiling approach is based on high-throughput sequencing of ribosome-protected mRNAs. RNAs not harboring ribosomes are removed by nuclease digestion leaving the so-called ribosome "footprints." The purified "footprint" RNA molecules are processed into DNA libraries and their individual abundance is determined by deep sequencing. Ribosome profiling reveals the portion of transcripts which are actually protein-coding and can be used for differential gene expression analysis addressing rates of protein synthesis, and translational control and efficiency.

Changes in prolactin receptor homodimer availability may cause late feathering in chickens

Chicken early (EF) and late feathering (LF) are sex-linked phenotypes conferred by wild-type k⁺ and dominant K alleles on chromosome Z, respectively. Besides prolactin (PRL) receptor (PRLR) and sperm flagellar 2 (SPEF2) genes, the K allele contains a fusion gene in which partially duplicated PRLR (dPRLR) and SPEF2 (dSPEF2) genes are linked in a tail-to-tail manner. The causative dPRLR gene encodes a C-terminal truncated receptor. LF chickens have short or no primaries at hatching; however, their feather growth rate is higher than that of EF chickens. This study aimed to elucidate the molecular basis of the K allele's biphasic effect on feather development. By 3'RACE and RT-PCR analyses, we demonstrated that dSPEF2 gene transcription occurred beyond all coding exons of the dPRLR gene on the opposite strand and that dPRLR mRNA was less abundant than PRLR mRNA. In addition, a 5'UTR splice variant (SPV) of PRL receptor mRNAs was increased in LF chickens. In vitro expression analysis of 5'UTR linked to the luciferase reporter gene revealed higher translation efficiency of SPV. RT-qPCR showed that the dPRLR mRNA level was higher in embryos; conversely, SPV was higher in hatched chickens, as was dSPEF2 mRNA. These findings suggest that the K allele inhibits feather development at the fetal stage by expressing dPRLR to attenuate PRLR function and promotes feather growth after hatching by increasing PRLR through dSPEF2 mRNA expression. Increased SPV may cause greater feather growth than that in EF chickens by increasing the availability of PRLR homodimers and enhancing PRL signaling.

Genome-scale analysis of syngas fermenting acetogenic bacteria reveals the translational regulation for its autotrophic growth

Background

Acetogenic bacteria constitute promising biocatalysts for the conversion of CO₂/H₂ or synthesis gas (H₂/CO/CO₂) into biofuels and value-added biochemicals. These microorganisms are naturally capable of autotrophic growth via unique acetogenesis metabolism. Despite their biosynthetic potential for commercial applications, a systemic understanding of the transcriptional and translational regulation of the acetogenesis metabolism remains unclear.

Results

By integrating genome-scale transcriptomic and translatomic data, we explored the regulatory logic of the acetogenesis to convert CO₂ into biomass and metabolites in Eubacterium limosum. The results indicate that majority of genes associated with autotrophic growth including the Wood-Ljungdahl pathway, the reduction of electron carriers, the energy conservation system, and gluconeogenesis were transcriptionally upregulated. The translation efficiency of genes in cellular respiration and electron bifurcation was also highly enhanced. In contrast, the transcriptionally abundant genes involved in the carbonyl branch of the Wood-Ljungdahl pathway, as well as the ion-translocating complex and ATP synthase complex in the energy conservation system, showed decreased translation efficiency. The translation efficiencies of genes were regulated by 5′UTR secondary structure under the autotrophic growth condition.

Conclusions

The results illustrated that the acetogenic bacteria reallocate protein synthesis, focusing more on the translation of genes for the generation of reduced electron carriers via electron bifurcation, rather than on those for carbon metabolism under autotrophic growth.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5238-0) contains supplementary material, which is available to authorized users.

Comparative ribosome profiling uncovers a dominant role for translational control in Toxoplasma gondii

Background

The lytic cycle of the protozoan parasite Toxoplasma gondii, which involves a brief sojourn in the extracellular space, is characterized by defined transcriptional profiles. For an obligate intracellular parasite that is shielded from the cytosolic host immune factors by a parasitophorous vacuole, the brief entry into the extracellular space is likely to exert enormous stress. Due to its role in cellular stress response, we hypothesize that translational control plays an important role in regulating gene expression in Toxoplasma during the lytic cycle. Unlike transcriptional profiles, insights into genome-wide translational profiles of Toxoplasma gondii are lacking.

Methods

We have performed genome-wide ribosome profiling, coupled with high throughput RNA sequencing, in intracellular and extracellular Toxoplasma gondii parasites to investigate translational control during the lytic cycle.

Results

Although differences in transcript abundance were mostly mirrored at the translational level, we observed significant differences in the abundance of ribosome footprints between the two parasite stages. Furthermore, our data suggest that mRNA translation in the parasite is potentially regulated by mRNA secondary structure and upstream open reading frames.

Conclusion

We show that most of the Toxoplasma genes that are dysregulated during the lytic cycle are translationally regulated.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4362-6) contains supplementary material, which is available to authorized users.

Ribosome profiling reveals translational regulation of mammalian cells in response to hypoxic stress

Background

Retinal pigment epithelium (RPE) cells transfer oxygen and nutrients from choroid to the neural retina. Reduced oxygen to RPE perturbs development and functions of blood vessels in retina. Previous efforts of genome-wide studies have been largely focused on transcriptional changes of cells in response to hypoxia. Recently developed ribosome profiling provides an opportunity to study genome-wide translational changes. To gain systemic insights into the transcriptional and translational regulation of cellular in response to hypoxic stress, we used simultaneous RNA sequencing and ribosome profiling on an RPE cells line, ARPE-19, under hypoxia condition.

Results

Both HIF-1α and EPAS1 (HIF-2α) proteins were stabilized in ARPE-19 under hypoxic stress treatment at 1 h, 2 h and 4 h. Analysis of simultaneous RNA sequencing and ribosome profiling data showed genome-wide gene expression changes at both transcriptional and translational levels. Comparative analysis of ribosome profiling and RNA-seq data revealed that hypoxia induced changes of more genes at the translational than the transcriptional levels. Ribosomes densities at 5′ untranslated region (UTR) significantly increased under hypoxic stress. Interestingly, the increase in ribosome densities at 5′ UTR is positively correlated with the presence of upstream open reading frames (uORFs) in the 5′ UTR of mRNAs.

Conclusion

Our results characterized translational profiles of mRNAs for a RPE cell line in response to hypoxia. In particular, uORFs play important roles in the regulation of translation efficiency by affecting ribosomes loading onto mRNAs. This study provides the first attempt to understand translational response of mammalian cells under hypoxic condition.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3996-8) contains supplementary material, which is available to authorized users.

The glpD gene is a novel reporter gene for E. coli that is superior to established reporter genes like lacZ and gusA

Reporter genes facilitate the characterization of promoter activities, transcript stabilities, translational efficiencies, or intracellular localization. Various reporter genes for Escherichia coli have been established, however, most of them have drawbacks like transcript instability or the inability to be used in genetic selections. Therefore, the glpD gene encoding glycerol-3-phosphate dehydrogenase was introduced as a novel reporter gene for E. coli. The enzymatic assay was optimized, and it was verified that growth on glycerol strictly depends on the presence of GlpD. The 5'-UTRs of three E. coli genes were chosen and cloned upstream of the new reporter gene glpD as well as the established reporter genes lacZ and gusA. Protein and transcript levels were quantified and translational efficiencies were calculated. The lacZ transcript was very unstable and its level highly depended on its translation, compromising its use as a reporter. The results obtained with gusA and glpD were similar, however, only glpD can be used for genetic selections. Therefore, glpD was found to be a superior novel reporter gene compared to the established reporter genes lacZ and gusA.

HCVIVdb: The hepatitis-C IRES variation database

Background

Sequence variability in the hepatitis C virus (HCV) genome has led to the development and classification of six genotypes and a number of subtypes. The HCV 5′ untranslated region mainly comprises an internal ribosomal entry site (IRES) responsible for cap-independent synthesis of the viral polyprotein and is conserved among all HCV genotypes.

Description

Considering the possible high impact of variations in HCV IRES on viral protein production and thus virus replication, we decided to collect the available data on known nucleotide variants in the HCV IRES and their impact on IRES function in translation initiation. The HCV IRES variation database (HCVIVdb) is a collection of naturally occurring and engineered mutation entries for the HCV IRES. Each entry contains contextual information pertaining to the entry such as the HCV genotypic background and links to the original publication. Where available, quantitative data on the IRES efficiency in translation have been collated along with details on the reporter system used to generate the data. Data are displayed both in a tabular and graphical formats and allow direct comparison of results from different experiments. Together the data provide a central resource for researchers in the IRES and hepatitis C-oriented fields.

Conclusion

The collation of over 1900 mutations enables systematic analysis of the HCV IRES. The database is mainly dedicated to detailed comparative and functional analysis of all the HCV IRES domains, which can further lead to the development of site-specific drug designs and provide a guide for future experiments. HCVIVdb is available at http://www.hcvivdb.org.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-016-0804-6) contains supplementary material, which is available to authorized users.

Design, synthesis and biological evaluation of dinucleotide mRNA cap analog containing propargyl moiety

The first example of the synthesis of new dinucleotide cap analog containing propargyl group such as m(7,3'-O-propargyl)G[5']ppp[5']G is reported. The effect of propargyl cap analog with standard cap was evaluated with respect to their capping efficiency, in vitro T7 RNA polymerase transcription efficiency, and translation activity using cultured HeLa cells. It is noteworthy that propargyl cap analog outperforms standard cap by 3.1 fold in terms of translational properties. The propargyl cap analog forms a more stable complex with translation initiation factor eIF4E based on the molecular modeling studies.

The Code of Silence: Widespread Associations Between Synonymous Codon Biases and Gene Function

Some mutations in gene coding regions exchange one synonymous codon for another, and thus do not alter the amino acid sequence of the encoded protein. Even though they are often called 'silent,' these mutations may exhibit a plethora of effects on the living cell. Therefore, they are often selected during evolution, causing synonymous codon usage biases in genomes. Comparative analyses of bacterial, archaeal, fungal, and human cancer genomes have found many links between a gene's biological role and the accrual of synonymous mutations during evolution. In particular, highly expressed genes in certain functional categories are enriched with optimal codons, which are decoded by the abundant tRNAs, thus enhancing the speed and accuracy of the translating ribosome. The set of genes exhibiting codon adaptation differs between genomes, and these differences show robust associations to organismal phenotypes. In addition to selection for translation efficiency, other distinct codon bias patterns have been found in: amino acid starvation genes, cyclically expressed genes, tissue-specific genes in animals and plants, oxidative stress response genes, cellular differentiation genes, and oncogenes. In addition, genomes of organisms harboring tRNA modifications exhibit particular codon preferences. The evolutionary trace of codon bias patterns across orthologous genes may be examined to learn about a gene's relevance to various phenotypes, or, more generally, its function in the cell.

Directed evolution study unveiling key sequence factors that affect translation efficiency in Escherichia coli

Synonymous mutations in protein coding genes significantly impact translation efficiency. We synthesized a pair of genes encoding green fluorescent protein that were separated by 160 synonymous mutations to investigate key factors that affect translation efficiency. One sequence was optimized for Escherichia coli (GFP(Eco)) and the other for Bacillus subtilis (GFP(Bsu)). When the genes were expressed in E. coli, GFP(Eco) fluoresced 12-fold stronger than GFP(Bsu), confirming the suboptimal nature of the GFP(Bsu) gene. We then employed directed evolution to improve the expression of GFP(Bsu). Random mutagenesis and DNA shuffling was used to generate mutant libraries, which were screened for fluorescence. A variant showing 6-fold fluorescence enhancement was identified, which contained a single mutation (G10A) in a rare codon for Gly-4. However, the substitution generated another type of rare codon, AGA, for Arg, suggesting that the improvement was caused by a factor other than the rare codon. We next applied saturation mutagenesis to Gly-4. The darkest variant contained a GGG codon (GFP(Bsu)-G) for Gly-4. Taking the location of the mutation into account, we hypothesized that destabilization of the mRNA secondary structure around the initiation codon improved the expression. We then randomized the nucleotide triplet in 5'-untranslated region (5'UTR) of GFP(Bsu), which is complementary to the Gly-4 codon. A variant showing 6-fold fluorescence enhancement was identified, which exhibited a destabilized secondary structure. When this 5'UTR sequence was combined with GFP(Bsu)-G, 22-fold fluorescent improvement was achieved. Collectively, the stability of the mRNA secondary structure around the initiation codon predominantly affected the translation efficiency.