The many faces of ribosome translocation along the mRNA: reading frame maintenance, ribosome frameshifting and translational bypassing.

In each round of translation elongation, the ribosome translocates along the mRNA by precisely one codon. Translocation is promoted by elongation factor G (EF-G) in bacteria (eEF2 in eukaryotes) and entails a number of precisely-timed large-scale structural rearrangements. As a rule, the movements of the ribosome, tRNAs, mRNA and EF-G are orchestrated to maintain the exact codon-wise step size. However, signals in the mRNA, as well as environmental cues, can change the timing and dynamics of the key rearrangements leading to recoding of the mRNA into production of trans-frame peptides from the same mRNA. In this review, we discuss recent advances on the mechanics of translocation and reading frame maintenance. Furthermore, we describe the mechanisms and biological relevance of non-canonical translocation pathways, such as hungry and programmed frameshifting and translational bypassing, and their link to disease and infection.

Assessing optimal: inequalities in codon optimization algorithms.

BACKGROUND: Custom genes have become a common resource in recombinant biology over the last 20 years due to the plummeting cost of DNA synthesis. These genes are often "optimized" to non-native sequences for overexpression in a non-native host by substituting synonymous codons within the coding DNA sequence (CDS). A handful of studies have compared native and optimized CDSs, reporting different levels of soluble product due to the accumulation of misfolded aggregates, variable activity of enzymes, and (at least one report of) a change in substrate specificity. No study, to the best of our knowledge, has performed a practical comparison of CDSs generated from different codon optimization algorithms or reported the corresponding protein yields. RESULTS: In our efforts to understand what factors constitute an optimized CDS, we identified that there is little consensus among codon-optimization algorithms, a roughly equivalent chance that an algorithm-optimized CDS will increase or diminish recombinant yields as compared to the native DNA, a near ubiquitous use of a codon database that was last updated in 2007, and a high variability of output CDSs by some algorithms. We present a case study, using KRas4B, to demonstrate that a median codon frequency may be a better predictor of soluble yields than the more commonly utilized CAI metric. CONCLUSIONS: We present a method for visualizing, analyzing, and comparing algorithm-optimized DNA sequences for recombinant protein expression. We encourage researchers to consider if DNA optimization is right for their experiments, and work towards improving the reproducibility of published recombinant work by publishing non-native CDSs.

The Standard Genetic Code can Evolve from a Two-Letter GC Code Without Information Loss or Costly Reassignments.

It is widely agreed that the standard genetic code must have been preceded by a simpler code that encoded fewer amino acids. How this simpler code could have expanded into the standard genetic code is not well understood because most changes to the code are costly. Taking inspiration from the recently synthesized six-letter code, we propose a novel hypothesis: the initial genetic code consisted of only two letters, G and C, and then expanded the number of available codons via the introduction of an additional pair of letters, A and U. Various lines of evidence, including the relative prebiotic abundance of the earliest assigned amino acids, the balance of their hydrophobicity, and the higher GC content in genome coding regions, indicate that the original two nucleotides were indeed G and C. This process of code expansion probably started with the third base, continued with the second base, and ended up as the standard genetic code when the second pair of letters was introduced into the first base. The proposed process is consistent with the available empirical evidence, and it uniquely avoids the problem of costly code changes by positing instead that the code expanded its capacity via the creation of new codons with extra letters.

Nucleotide composition of the Zika virus RNA genome and its codon usage.

BACKGROUND: RNA viruses have genomes with a distinct nucleotide composition and codon usage. We present the global characteristics of the RNA genome of Zika virus (ZIKV), an emerging pathogen within the Flavivirus genus. ZIKV was first isolated in 1947 in Uganda, caused a widespread epidemic in South and Central America and the Caribbean in 2015 and has recently been associated with microcephaly in newborns. METHODS: The nearly 11 kb positive-stranded RNA genome of ZIKV was analyzed for its nucleotide composition, also in the context of the folded RNA molecule. Nucleotide trends were investigated along the genome length by skew analyses and we analyzed the codons used for translation of the ZIKV proteins. RESULTS: ZIKV RNA has a biased nucleotide composition in being purine-rich and pyrimidine-poor. This preference for purines is a general characteristic of the mosquito-borne and tick-borne flaviviruses. The virus-specific nucleotide bias is further enriched in the unpaired, single-stranded regions of the structured ZIKV RNA genome, thus further imposing this ZIKV-specific signature. The codons used for translation of the ZIKV proteins is also unusual, but we show that it is the underlying bias in nucleotide composition of the viral RNA that largely dictates these codon preferences. CONCLUSIONS: The ZIKV RNA genome has a biased nucleotide composition that dictates the codon usage of this flavivirus. We discuss the evolutionary scenarios and molecular mechanisms that may be responsible for these distinctive ZIKV RNA genome features.

Comprehensive analysis of the overall codon usage patterns in equine infectious anemia virus.

BACKGROUND: Equine infectious anemia virus (EIAV) is an important animal model for understanding the relationship between viral persistence and the host immune response during lentiviral infections. Comparison and analysis of the codon usage model between EIAV and its hosts is important for the comprehension of viral evolution. In our study, the codon usage pattern of EIAV was analyzed from the available 29 full-length EIAV genomes through multivariate statistical methods. FINDING: Effective number of codons (ENC) suggests that the codon usage among EIAV strains is slightly biased. The ENC-plot analysis demonstrates that mutation pressure plays a substantial role in the codon usage pattern of EIAV, whereas other factors such as geographic distribution and host translation selection also take part in the process of EIAV evolution. Comparative analysis of codon adaptation index (CAI) values among EIAV and its hosts suggests that EIAV utilize the translational resources of horse more efficiently than that of donkey. CONCLUSION: The codon usage bias in EIAV is slight and mutation pressure is the main factor that affects codon usage variation in EIAV. These results suggest that EIAV genomic biases are the result of the co-evolution of genome composition and the ability to evade the host's immune response.

Chemical modifications to mRNA nucleobases impact translation elongation and termination.

Messenger RNAs (mRNAs) serve as blueprints for protein synthesis by the molecular machine the ribosome. The ribosome relies on hydrogen bonding interactions between adaptor aminoacyl-transfer RNA molecules and mRNAs to ensure the rapid and faithful translation of the genetic code into protein. There is a growing body of evidence suggesting that chemical modifications to mRNA nucleosides impact the speed and accuracy of protein synthesis by the ribosome. Modulations in translation rates have downstream effects beyond protein production, influencing protein folding and mRNA stability. Given the prevalence of such modifications in mRNA coding regions, it is imperative to understand the consequences of individual modifications on translation. In this review we present the current state of our knowledge regarding how individual mRNA modifications influence ribosome function. Our comprehensive comparison of the impacts of 16 different mRNA modifications on translation reveals that most modifications can alter the elongation step in the protein synthesis pathway. Additionally, we discuss the context dependence of these effects, highlighting the necessity of further study to uncover the rules that govern how any given chemical modification in an mRNA codon is read by the ribosome.

RET codon 804 mutations in multiple endocrine neoplasia 2: genotype-phenotype correlations and implications in clinical management.

Multiple endocrine neoplasia type 2 (MEN 2) is a genetic syndrome caused by germline mutations in the RET proto-oncogene. These mutations cause changes in either the cysteine-rich extracellular domain or, less commonly, the non-cysteine intracellular domains of the RET protein. The genotype-phenotype correlations of classical cysteine RET mutations have been the subject of several comprehensive reviews. Less is known about the characteristics of the non-cysteine RET mutations. Studies of familial medullary thyroid cancer and MEN 2A kindreds carrying non-cysteine RET mutations have revealed a wide array of phenotypes, variable penetrance, and a diverse clinical course. The observed heterogeneity in disease expression has important diagnostic, therapeutic and prognostic implications. This review summarizes the genotypic and phenotypic characteristics of RET codon 804 mutation, a prototype for the less well-defined non-cysteine RET mutations associated with MEN 2.

Evaluating phylogenetic informativeness and data-type usage for new protein-coding genes across Vertebrata.

As a resource for vertebrate phylogenetics, we developed 75 new protein-coding genes using a combination of expressed sequence tags (ESTs) available in Genbank, and targeted amplification of complementary DNA (cDNA). In addition, we performed three additional analyses in order to assess the utility of our approach. First, we profiled the phylogenetic informativeness of these new markers using the online program PhyDesign. Next, we compared the utility of four different data-types used in phylogenetics: nucleotides (NUCL), amino acids (AA), 1st and 2nd codon positions only (N12), and modified sequences to account for codon degeneracy (DEGEN1; Regier et al., 2010). Lastly, we use these new markers to construct a vertebrate phylogeny and address the uncertain relationship between higher-level mammal groups: monotremes, marsupials, and placentals. Our results show that phylogenetic informativeness of the 75 new markers varies, both in the amount of phylogenetic signal and optimal timescale. When comparing the four data-types, we find that the NUCL data-type, due to the high level of phylogenetic signal, performs the best across all divergence times. The remaining three data-types (AA, N12, DEGEN1) are less subject to homoplasy, but have greatly reduced levels of phylogenetic signal relative to NUCL. Our phylogenetic inference supports the Theria hypothesis of mammalian relationships, with marsupials and placentals being sister groups.

Analysis of codon usage in Newcastle disease virus.

In this study, the relative synonymous codon usage (RSCU) values, effective number of codon (ENC) values, nucleotide contents, and dinucleotide were used to investigate codon usage pattern of each protein-coding gene and genome among 31 Newcastle disease virus (NDV) isolates. The result shows that the overall extent of codon usage bias in NDV is low (mean ENC = 56.15 > 40). The good correlation between the (C + G)(12)% and (G + C)(3)% suggests that the mutational pressure, rather than natural selection, is the main factor that determines the codon usage bias and base component in NDV. It is observed that synonymous codon usage pattern in NDV genes is gene function and geography specific, but not host specific. By contrasting synonymous codon usage patterns of different NDV isolates, we suggest that more than one genotype of NDV circulates in waterfowl in USA; and gene length has no significant effect on the variations of synonymous codon usage in these virus genes. CpG under-represented is a characteristic for NDV to fit in its host. These results not only provide an insight into the variation of codon usage pattern among the genomes of NDV, but also may help in understanding the processes governing the evolution of NDV.

Kinetics of stop codon recognition by release factor 1.

Recognition of stop codons by class I release factors is a fundamental step in the termination phase of protein synthesis. Since premature termination is costly to the cell, release factors have to efficiently discriminate between stop and sense codons. To understand the mechanism of discrimination between stop and sense codons, we developed a new, pre-steady state kinetic assay to monitor the interaction of RF1 with the ribosome. Our results show that RF1 associates with similar association rate constants with ribosomes programmed with stop or sense codons. However, dissociation of RF1 from sense codons is as much as 3 orders of magnitude faster than from stop codons. Interestingly, the affinity of RF1 for ribosomes programmed with different sense codons does not correlate with the defects in peptide release. Thus, discrimination against sense codons is achieved with both an increase in the dissociation rates and a decrease in the rate of peptide release. These results suggest that sense codons inhibit conformational changes necessary for RF1 to stably bind to the ribosome and catalyze peptide release.

Amino acid composition of proteins: Selection against the genetic code.

Distribution of amino acids in 68 representative proteins is compared with their distribution among 61 codons of the genetic code. Average amounts of lysine, aspartic acid, glutamic acid, and alanine are above the levels anticipated from the genetic code, and arginine, serine, leucine, cysteine, proline, and histidine are below such levels. Arginine plus lysine account for 11.0 percent of codons and aspartic acid plus glutamic acid account for 11.3 percent; thus the average charge is roughly neutral.

Patterns of codon usage in two ciliates that reassign the genetic code: Tetrahymena thermophila and Paramecium tetraurelia.

We used the recently sequenced genomes of the ciliates Tetrahymena thermophila and Paramecium tetraurelia to analyze the codon usage patterns in both organisms; we have analyzed codon usage bias, Gln codon usage, GC content and the nucleotide contexts of initiation and termination codons in Tetrahymena and Paramecium. We also studied how these trends change along the length of the genes and in a subset of highly expressed genes. Our results corroborate some of the trends previously described in Tetrahymena, but also negate some specific observations. In both genomes we found a strong bias toward codons with low GC content; however, in highly expressed genes this bias is smaller and codons ending in GC tend to be more frequent. We also found that codon bias increases along gene segments and in highly expressed genes and that the context surrounding initiation and termination codons are always AT rich. Our results also suggest differences in the efficiency of translation of the reassigned stop codons between the two species and between the reassigned codons. Finally, we discuss some of the possible causes for such translational efficiency differences.

Analysis of synonymous codon usage in 11 human bocavirus isolates.

Human Bocavirus (HBoV) is a novel virus which can cause respiratory tract disease in infants or children. In this study, the codon usage bias and the base composition variations in the available 11 complete HBoV genome sequences have been investigated. Although, there is a significant variation in codon usage bias among different HBoV genes, codon usage bias in HBoV is a little slight, which is mainly determined by the base compositions on the third codon position and the effective number of codons (ENC) value. The results of correspondence analysis (COA) and Spearman's rank correlation analysis reveals that the G+C compositional constraint is the main factor that determines the codon usage bias in HBoV and the gene's function also contributes to the codon usage in this virus. Moreover, it was found that the hydrophobicity of each protein and the gene length are also critical in affecting these viruses' codon usage, although they were less important than that of the mutational bias and the genes' function. At last, the relative synonymous codon usage (RSCU) of 44 genes from these 11 HBoV isolates is analyzed using a hierarchical cluster method. The result suggests that genes with same function yet from different isolates are classified into the same lineage and it does not depend on geographical location. These conclusions not only can offer an insight into the codon usage patterns and gene classification of HBoV, but also may help in increasing the efficiency of gene delivery/expression systems.

Recodon: coalescent simulation of coding DNA sequences with recombination, migration and demography.

BACKGROUND: Coalescent simulations have proven very useful in many population genetics studies. In order to arrive to meaningful conclusions, it is important that these simulations resemble the process of molecular evolution as much as possible. To date, no single coalescent program is able to simulate codon sequences sampled from populations with recombination, migration and growth. RESULTS: We introduce a new coalescent program, called Recodon, which is able to simulate samples of coding DNA sequences under complex scenarios in which several evolutionary forces can interact simultaneously (namely, recombination, migration and demography). The basic codon model implemented is an extension to the general time-reversible model of nucleotide substitution with a proportion of invariable sites and among-site rate variation. In addition, the program implements non-reversible processes and mixtures of different codon models. CONCLUSION: Recodon is a flexible tool for the simulation of coding DNA sequences under realistic evolutionary models. These simulations can be used to build parameter distributions for testing evolutionary hypotheses using experimental data. Recodon is written in C, can run in parallel, and is freely available from http://darwin.uvigo.es/.

Methods for selecting fixed-effect models for heterogeneous codon evolution, with comments on their application to gene and genome data.

BACKGROUND: Models of codon evolution have proven useful for investigating the strength and direction of natural selection. In some cases, a priori biological knowledge has been used successfully to model heterogeneous evolutionary dynamics among codon sites. These are called fixed-effect models, and they require that all codon sites are assigned to one of several partitions which are permitted to have independent parameters for selection pressure, evolutionary rate, transition to transversion ratio or codon frequencies. For single gene analysis, partitions might be defined according to protein tertiary structure, and for multiple gene analysis partitions might be defined according to a gene's functional category. Given a set of related fixed-effect models, the task of selecting the model that best fits the data is not trivial. RESULTS: In this study, we implement a set of fixed-effect codon models which allow for different levels of heterogeneity among partitions in the substitution process. We describe strategies for selecting among these models by a backward elimination procedure, Akaike information criterion (AIC) or a corrected Akaike information criterion (AICc). We evaluate the performance of these model selection methods via a simulation study, and make several recommendations for real data analysis. Our simulation study indicates that the backward elimination procedure can provide a reliable method for model selection in this setting. We also demonstrate the utility of these models by application to a single-gene dataset partitioned according to tertiary structure (abalone sperm lysin), and a multi-gene dataset partitioned according to the functional category of the gene (flagellar-related proteins of Listeria). CONCLUSION: Fixed-effect models have advantages and disadvantages. Fixed-effect models are desirable when data partitions are known to exhibit significant heterogeneity or when a statistical test of such heterogeneity is desired. They have the disadvantage of requiring a priori knowledge for partitioning sites. We recommend: (i) selection of models by using backward elimination rather than AIC or AICc, (ii) use a stringent cut-off, e.g., p = 0.0001, and (iii) conduct sensitivity analysis of results. With thoughtful application, fixed-effect codon models should provide a useful tool for large scale multi-gene analyses.

ChloroMitoCU: Codon patterns across organelle genomes for functional genomics and evolutionary applications.

Organelle genomes are widely thought to have arisen from reduction events involving cyanobacterial and archaeal genomes, in the case of chloroplasts, or α-proteobacterial genomes, in the case of mitochondria. Heterogeneity in base composition and codon preference has long been the subject of investigation of topics ranging from phylogenetic distortion to the design of overexpression cassettes for transgenic expression. From the overexpression point of view, it is critical to systematically analyze the codon usage patterns of the organelle genomes. In light of the importance of codon usage patterns in the development of hyper-expression organelle transgenics, we present ChloroMitoCU, the first-ever curated, web-based reference catalog of the codon usage patterns in organelle genomes. ChloroMitoCU contains the pre-compiled codon usage patterns of 328 chloroplast genomes (29,960 CDS) and 3,502 mitochondrial genomes (49,066 CDS), enabling genome-wide exploration and comparative analysis of codon usage patterns across species. ChloroMitoCU allows the phylogenetic comparison of codon usage patterns across organelle genomes, the prediction of codon usage patterns based on user-submitted transcripts or assembled organelle genes, and comparative analysis with the pre-compiled patterns across species of interest. ChloroMitoCU can increase our understanding of the biased patterns of codon usage in organelle genomes across multiple clades. ChloroMitoCU can be accessed at: http://chloromitocu.cgu.edu.tw/.

Ribo-seq enlightens codon usage bias.

Codon usage is biased between lowly and highly expressed genes in a genome-specific manner. This universal bias has been well assessed in some unicellular species, but remains problematic to assess in more complex species. We propose a new method to compute codon usage bias based on genome wide translational data. A new technique based on sequencing of ribosome protected mRNA fragments (Ribo-seq) allowed us to rank genes and compute codon usage bias with high precision for a great variety of species, including mammals. Genes ranking using Ribo-Seq data confirms the influence of the tRNA pool on codon usage bias and shows a decreasing bias in multicellular species. Ribo-Seq analysis also makes possible to detect preferred codons without information on genes function.

Analysis of Ribosome Stalling and Translation Elongation Dynamics by Deep Learning.

Ribosome stalling is manifested by the local accumulation of ribosomes at specific codon positions of mRNAs. Here, we present ROSE, a deep learning framework to analyze high-throughput ribosome profiling data and estimate the probability of a ribosome stalling event occurring at each genomic location. Extensive validation tests on independent data demonstrated that ROSE possessed higher prediction accuracy than conventional prediction models, with an increase in the area under the receiver operating characteristic curve by up to 18.4%. In addition, genome-wide statistical analyses showed that ROSE predictions can be well correlated with diverse putative regulatory factors of ribosome stalling. Moreover, the genome-wide ribosome stalling landscapes of both human and yeast computed by ROSE recovered the functional interplays between ribosome stalling and cotranslational events in protein biogenesis, including protein targeting by the signal recognition particles and protein secondary structure formation. Overall, our study provides a novel method to complement the ribosome profiling techniques and further decipher the complex regulatory mechanisms underlying translation elongation dynamics encoded in the mRNA sequence.

Expression of the regulatory protein LndI for landomycin E production in Streptomyces globisporus 1912 is controlled by the availability of tRNA for the rare UUA codon.

The gene lndI encodes the activator of landomycin biosynthesis. The utilization of LndI-EGFP fusions led us to investigate the temporal pattern of this gene expression and demonstrated the delay between lndI transcription and translation. The TTA codon in lndI is thought to be the reason for this delay. The replacement of TTA with CTC cancelled the pause between lndI transcription and the translation. The wild-type of the lndI gene is not expressed in the Streptomyces coelicolor bldA- mutant strain, indicating the importance of the bldA tRNA in its mRNA translation.

Mutations in carboxy-terminal part of E2 including PKR/eIF2alpha phosphorylation homology domain and interferon sensitivity determining region of nonstructural 5A of hepatitis C virus 1b: their correlation with response to interferon monotherapy and viral load.

AIM: To study the amino acid substitutions in the carboxy (C)-terminal part of E2 protein and in the interferon (IFN) sensitivity determining region (ISDR) and their correlation with response to IFN and viral load in 85 hepatitis C virus (HCV)-1b-infected patients treated with IFN. METHODS: The C-terminal part of E2 (codons 617-711) including PKR/eIF2alpha phosphorylation homology domain (PePHD) and ISDR was sequenced in 85 HCV-1b-infected patients treated by IFN monotherapy. RESULTS: The amino acid substitutions in PePHD detected only in 4 of 85 patients were not correlated either with response to IFN or with viral load. The presence of substitutions in a N-terminal variable region (codons 617-641) in the C-terminal part of E2 was significantly correlated with both small viral load (33.9% vs 13.8%, P = 0.0394) and sustained response to IFN (25.0% vs 6.9%, P = 0.0429). Four or more substitutions in ISDR were significantly correlated with both small viral load (78.6% vs 16.2%, P < 0.0001) and sustained response to IFN (85.7% vs 2.9%, P < 0.0001). In multivariate analysis, ISDR in nonstructural (NS) 5A (OR = 0.39, P < 0.0001) and N-terminal variable region (OR = 0.51, P = 0.039) was selected as the independent predictors for small viral load, and ISDR (OR = 39.0, P < 0.0001) was selected as the only independent predictor for sustained response. CONCLUSION: The N-terminal variable region in the C-terminal part of E2 correlates with both response to IFN monotherapy and viral load and is one of the factors independently associated with a small viral load.