Codon Optimality Is a Major Determinant of mRNA Stability

BackgroundRecent experiments have shown that codon optimality is a major determinant of mRNA stability in Saccharomyces cerevisiae and that this phenomenon may be conserved in Escherichia coli and some metazoans, although work in Neurospora crassa is not consistent with this model.ResultsWe examined the association between codon optimality and mRNA stability in the fission yeast Schizosaccharomyces pombe. Our analysis revealed the following points. First, we observe a genome-wide association between codon optimality and mRNA stability also in S. pombe, suggesting evolutionary conservation of the phenomenon. Second, in both S. pombe and S. cerevisiae, mRNA synthesis rates are also correlated at the genome-wide analysis with codon optimality, suggesting that the long-appreciated association between codon optimality and mRNA abundance is due to regulation of both mRNA synthesis and degradation. However, when we examined correlation of codon optimality and either mRNA half-lives or synthesis rates controlling for mRNA abundance, codon optimality was still positively correlated with mRNA half-lives in S. cerevisiae, but the association was no longer significant for mRNA half-lives in S. pombe or for synthesis rates in either organism. This illustrates how only the pairwise analysis of multiple correlating variables may limit these types of analyses. Finally, in S. pombe, codon optimality is associated with known DNA/RNA sequence motifs that are associated with mRNA production/stability, suggesting these two features have been under similar selective pressures for optimal gene expression.ConclusionsConsistent with the emerging body of studies, this study suggests that the association between codon optimality and mRNA stability may be a broadly conserved phenomenon. It also suggests that the association can be explained at least in part by independent adaptations of codon optimality and other transcript features for elevated expression during evolution.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3237-6) contains supplementary material, which is available to authorized users.

In this work, we show that Not4 and Not5 from the Ccr4-Not complex modulate translation elongation dynamics and change ribosome A-site dwelling occupancy in a codon-dependent fashion. These codon-specific changes in not5Δ cells are very robust and independent of codon position within the mRNA, the overall mRNA codon composition, or changes of mRNA expression levels. They inversely correlate with codon-specific changes in cells depleted for eIF5A and positively correlate with those in cells depleted for ribosome-recycling factor Rli1. Not5 resides in punctate loci, co-purifies with ribosomes and Rli1, but not with eIF5A, and limits mRNA solubility. Overexpression of wild-type or non-complementing Rli1 and loss of Rps7A ubiquitination enable Not4 E3 ligase-dependent translation of polyarginine stretches. We propose that Not4 and Not5 modulate translation elongation dynamics to produce a soluble proteome by Rps7A ubiquitination, dynamic condensates that limit mRNA solubility and exclude eIF5A, and a moonlighting function of Rli1.

For a successful yet controlled immune response, cells need to specifically destabilize inflammatory mRNAs but prevent premature removal of those still used. The regulatory circuits controlling quality and timing in the global inflammatory mRNA decay are not understood. Here, we show that the mRNA-destabilizing function of the AU-rich element-binding protein tristetraprolin (TTP) is inversely regulated by the p38 MAPK activity profile such that after inflammatory stimulus the TTP-dependent decay is initially limited to few mRNAs. With time, the TTP-dependent decay gradually spreads resulting in cumulative elimination of one third of inflammation-induced unstable mRNAs in macrophages in vitro. We confirmed this sequential decay model in vivo since LPS-treated mice with myeloid TTP ablation exhibited similar cytokine dysregulation profile as macrophages. The mice were hypersensitive to LPS but otherwise healthy with no signs of hyperinflammation seen in conventional TTP knockout mice demonstrating the requirement for myeloid TTP in re-installment but not maintenance of immune homeostasis. These findings reveal a TTP- and p38 MAPK-dominated regulatory mechanism that is vital for balancing acute inflammation by a temporally and qualitatively controlled mRNA decay.

Alzheimer's disease (AD) is a neurodegenerative disorder with few biomarkers even though it impacts a relatively large portion of the population and is predicted to affect significantly more individuals in the future. Neuroimaging has been used in concert with genetic information to improve our understanding in relation to how AD arises and how it can be potentially diagnosed. Additionally, evidence suggests synonymous variants can have a functional impact on gene regulatory mechanisms, including those related to AD. Some synonymous codons are preferred over others leading to a codon bias. The bias can arise with respect to codons that are more or less frequently used in the genome. A bias can also result from optimal and non-optimal codons, which have stronger and weaker codon anti-codon interactions, respectively. Although association tests have been utilized before to identify genes associated with AD, it remains unclear how codon bias plays a role and if it can improve rare variant analysis. In this work, rare variants from whole-genome sequencing from the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort were binned into genes using BioBin. An association analysis of the genes with AD-related neuroimaging biomarker was performed using SKAT-O. While using all synonymous variants we did not identify any genomewide significant associations, using only synonymous variants that affected codon frequency we identified several genes as significantly associated with the imaging phenotype. Additionally, significant associations were found using only rare variants that contains an optimal codon in among minor alleles and a non-optimal codon in the major allele. These results suggest that codon bias may play a role in AD and that it can be used to improve detection power in rare variant association analysis.

Attenuated Codon Optimality Contributes to Neural-Specific mRNA Decay in Drosophila

BackgroundThe regulation of messenger RNA (mRNA) stability has a profound impact on gene expression dynamics during embryogenesis. For example, in animals, maternally deposited mRNAs are degraded after fertilization to enable new developmental trajectories. Regulatory sequences in 3′ untranslated regions (3′UTRs) have long been considered the central determinants of mRNA stability. However, recent work indicates that the coding sequence also possesses regulatory information. Specifically, translation in cis impacts mRNA stability in a codon-dependent manner. However, the strength of this mechanism during embryogenesis, as well as its relationship with other known regulatory elements, such as microRNA, remains unclear.ResultsHere, we show that codon composition is a major predictor of mRNA stability in the early embryo. We show that this mechanism works in combination with other cis-regulatory elements to dictate mRNA stability in zebrafish and Xenopus embryos as well as in mouse and human cells. Furthermore, we show that microRNA targeting efficacy can be affected by substantial enrichment of optimal (stabilizing) or non-optimal (destabilizing) codons. Lastly, we find that one microRNA, miR-430, antagonizes the stabilizing effect of optimal codons during early embryogenesis in zebrafish.ConclusionsBy integrating the contributions of different regulatory mechanisms, our work provides a framework for understanding how combinatorial control of mRNA stability shapes the gene expression landscape.

Pan-cancer analyses of synonymous mutations based on tissue-specific codon optimality

Non-optimal codons have been reported to occur with higher frequency in the signal sequences of secretory proteins in Escherichia coli compared to both the 5’ mature region of secretory and non-secretory genes (134). Non-optimal codons are defined by their low-usage in a genome, low concentrations of cognate tRNA and are associated with translational pausing. This thesis presents data on two secretory proteins, maltose binding protein (MBP) and β-lactamase (Bla), which were chosen as model systems to study the role of non-optimal codons in export. For both proteins, the non-optimal codons were changed to an optimal codon within the synonymous codon family. For the optimised MBP allele (MBP-opt), there was an approximate 20-fold reduction in protein expression compared to MBP-wt. The expression of MBP-opt was partially recovered in strains deficient in multiple cytoplasmic and periplasmic proteases. This suggests that non-optimal codons are important to allow for the correct folding of the pre-protein before secretion, in particular for the signal peptide. In the β-lactamase study, the optimised allele (bla-opt) had four-fold less expression compared to the wild-type allele (bla-wt). To investigate the effect of the distribution of non-optimal codons in the signal sequence for expression, a series of mutations in the signal sequence, in which only a subset of the non-optimal codons were changed to optimal, were made. This revealed that non-optimal codons at the start of the signal sequence (bla-F6) were more important for maintaining expression of β-lactamase. Since bla confers resistance to ampicillin only when fully folded in the periplasm, a revertant study was carried out to select for clones that had higher resistance than those of bla-F6. This study revealed nine such revertants, which all had optimal to non-optimal codon reversions at various positions in the front half of the β-lactamase signal sequence. The results of these studies in both MBP and β- lactamase suggest that non-optimal codons are important for a time-critical early event in the export of proteins to the periplasm. Most prokaryotic sec dependent signal peptides have a positive charge at the Nterminus. Yet studies have shown that the positive charge is not required for export (74, 137, 174). To investigate this, the charge distribution of all sec dependent signal peptides was analysed, which revealed a massive bias for lysine at the second and third amino acid positions (P2 and P3) and a slight bias for arginine at P3. Whilst consistent with the requirement for a positive charge at the N-terminus, it raised the question as to why there was a massive bias for lysine and not arginine at P2? Investigating the codon usage at the various positions revealed that the lysine codon AAA is overrepresented at P2. Other studies have shown that codon and amino acid choice at P2 can affect translation initiation (159) and removal of f-Met by methionylamino peptidase (MAP) (51). To experimentally determine the selection constraints on codon and amino acid usage at P2 in secretory proteins, the signal sequences of MBP and alkaline phosphatase (AP) were fused to the mature region of β-lactamase. Both of these signal sequences had the lysine codon AAA at P2, which was changed to various other codons and the effect on export measured by resistance to ampicillin. The results of these experiments demonstrate that a positive charge is not required for protein export, and suggest that the dominant selective pressure at P2 is for codons that promote high translation initiation efficiencies.

Accuracy and efficiency are fundamental to mRNA translation. Codon usage bias is widespread across species. Despite the long-standing association between optimized codon usage and improved translation, our understanding of its evolutionary basis and functional effects remains limited. Drosophila is widely used to study codon usage bias, but genome-scale experimental data are scarce. Using high-resolution mass spectrometry data from Drosophila melanogaster, we show that optimal codons have lower translation errors than nonoptimal codons after accounting for these biases. Genomic-scale analysis of ribosome profiling data shows that optimal codons are translated more rapidly than nonoptimal codons. Although we find no long-term selection favoring synonymous mutations in D. melanogaster after diverging from D. simulans, we identify signatures of positive selection driving codon optimization in the D. melanogaster population. These findings expand our understanding of the functional consequences of codon optimization and serve as a foundation for future investigations.

The cell cycle is a temporal program that regulates DNA synthesis and cell division. When we compared the codon usage of cell cycle-regulated genes with that of other genes, we discovered that there is a significant preference for non-optimal codons. Moreover, genes encoding proteins that cycle at the protein level exhibit non-optimal codon preferences. Remarkably, cell cycle-regulated genes expressed in different phases display different codon preferences. Here, we show empirically that transfer RNA (tRNA) expression is indeed highest in the G2 phase of the cell cycle, consistent with the non-optimal codon usage of genes expressed at this time, and lowest toward the end of G1, reflecting the optimal codon usage of G1 genes. Accordingly, protein levels of human glycyl-, threonyl-, and glutamyl-prolyl tRNA synthetases were found to oscillate, peaking in G2/M phase. In light of our findings, we propose that non-optimal (wobbly) matching codons influence protein synthesis during the cell cycle. We describe a new mathematical model that shows how codon usage can give rise to cell-cycle regulation. In summary, our data indicate that cells exploit wobbling to generate cell cycle-dependent dynamics of proteins.

The cytoplasmic abundance of mRNAs is strictly controlled through a balance of production and degradation. Whereas the control of mRNA synthesis through transcription has been well characterized, less is known about the regulation of mRNA turnover, and a consensus model explaining the wide variations in mRNA decay rates remains elusive. Here, we combine non-invasive transcriptome-wide mRNA production and stability measurements with selective and acute perturbations to demonstrate that mRNA degradation is tightly coupled to the regulation of translation, and that a competition between translation initiation and mRNA decay -but not codon optimality or elongation- is the major determinant of mRNA stability in yeast. Our refined measurements also reveal a remarkably dynamic transcriptome with an average mRNA half-life of only 4.8 min - much shorter than previously thought. Furthermore, global mRNA destabilization by inhibition of translation initiation induces a dose-dependent formation of processing bodies in which mRNAs can decay over time.

The cytoplasmic abundance of mRNAs is strictly controlled through a balance of production and degradation. Whereas the control of mRNA synthesis through transcription has been well characterized, less is known about the regulation of mRNA turnover, and a consensus model explaining the wide variations in mRNA decay rates remains elusive. Here, we combine non-invasive transcriptome-wide mRNA production and stability measurements with selective and acute perturbations to demonstrate that mRNA degradation is tightly coupled to the regulation of translation, and that a competition between translation initiation and mRNA decay -but not codon optimality or elongation- is the major determinant of mRNA stability in yeast. Our refined measurements also reveal a remarkably dynamic transcriptome with an average mRNA half-life of only 4.8 min - much shorter than previously thought. Furthermore, global mRNA destabilization by inhibition of translation initiation induces a dose-dependent formation of processing bodies in which mRNAs can decay over time.

In protein-coding genes, synonymous codon usage and amino acid composition correlate to expression in some eukaryotes, and may result from translational selection. Here, we studied large-scale RNA-seq data from three divergent arthropod models, including cricket (Gryllus bimaculatus), milkweed bug (Oncopeltus fasciatus), and the amphipod crustacean Parhyale hawaiensis, and tested for optimization of codon and amino acid usage relative to expression level. We report strong signals of AT3 optimal codons (those favored in highly expressed genes) in G. bimaculatus and O. fasciatus, whereas weaker signs of GC3 optimal codons were found in P. hawaiensis, suggesting selection on codon usage in all three organisms. Further, in G. bimaculatus and O. fasciatus, high expression was associated with lowered frequency of amino acids with large size/complexity (S/C) scores in favor of those with intermediate S/C values; thus, selection may favor smaller amino acids while retaining those of moderate size for protein stability or conformation. In P. hawaiensis, highly transcribed genes had elevated frequency of amino acids with large and small S/C scores, suggesting a complex dynamic in this crustacean. In all species, the highly transcribed genes appeared to favor short proteins, high optimal codon usage, specific amino acids, and were preferentially involved in cell-cycling and protein synthesis. Together, based on examination of 1,680,067, 1,667,783, and 1,326,896 codon sites in G. bimaculatus, O. fasciatus, and P. hawaiensis, respectively, we conclude that translational selection shapes codon and amino acid usage in these three Pancrustacean arthropods.

Tremella fuciformis is a popular edible fungus with fruiting bodies that can be produced in large quantities at low costs, while it is easy to transform and cultivate as yeast. This makes it an attractive potential bioreactor. Enhanced heterologous gene expression through codon optimization would be useful, but until now codon usage preferences in T. fuciformis remain unknown. To precisely determine the preferred codon usage of T. fuciformis we sequenced the genome of strain Tr26 resulting in a 24.2 Mb draft genome with 10,040 predicted genes. 3288 of the derived predicted proteins matched the UniProtKB/Swiss-Prot databases with 40% or more similarity. Corresponding gene models of this subset were subsequently optimized through repetitive comparison of alternative start codons and selection of best length matching gene models. For experimental confirmation of gene models, 96 random clones from an existing T. fuciformis cDNA library were sequenced, generating 80 complete CDSs. Calculated optimal codons for the 3288 predicted and the 80 cloned CDSs were highly similar, indicating sufficient accuracy of predicted gene models for codon usage analysis. T. fuciformis showed a strong preference for C and then G at the third base pair position of used codons, while average GC content of predicted genes was slightly higher than the total genome sequence average. Most optimal codons ended in C or G except for one, and an increased frequency of C ending codons was observed in genes with higher expression levels. Surprisingly, the preferred codon usage in T. fuciformis strongly differed from T. mesenterica and C. neoformans. Instead, optimal codon usage was similar to more distant related species such as Ustilago maydis and Neurospora crassa. Despite much higher overall sequence homology between T. fuciformis and T. mesenterica, only 7 out of 21 optimal codons were equal, whereas T. fuciformis shared up to 20 out of 21 optimal codons with other species. Clearly, codon usage in Tremella can differ largely and should be estimated for individual species. The precise identification of optimal and high expression related codons is therefore an important step in the development of T. fuciformis as a bioreactor system.

Codon usage biases are found in all eukaryotic and prokaryotic genomes, and preferred codons are more frequently used in highly expressed genes. The effects of codon usage on gene expression were previously thought to be mainly mediated by its impacts on translation. Here, we show that codon usage strongly correlates with both protein and mRNA levels genome-wide in the filamentous fungus Neurospora Gene codon optimization also results in strong up-regulation of protein and RNA levels, suggesting that codon usage is an important determinant of gene expression. Surprisingly, we found that the impact of codon usage on gene expression results mainly from effects on transcription and is largely independent of mRNA translation and mRNA stability. Furthermore, we show that histone H3 lysine 9 trimethylation is one of the mechanisms responsible for the codon usage-mediated transcriptional silencing of some genes with nonoptimal codons. Together, these results uncovered an unexpected important role of codon usage in ORF sequences in determining transcription levels and suggest that codon biases are an adaptation of protein coding sequences to both transcription and translation machineries. Therefore, synonymous codons not only specify protein sequences and translation dynamics, but also help determine gene expression levels.