Non-optimal codon usage affects expression, structure and function of clock protein FRQ

Transcription factor CBF-1 is critical for circadian gene expression by modulating WHITE COLLAR complex recruitment to the frq locus

Transcription of the Neurospora crassa circadian clock gene frequency (frq) is an essential process in the negative feedback loop that controls circadian rhythms. WHITE COLLAR 1 (WC-1) and WHITE COLLAR 2 (WC-2) forms the WC complex (WCC) that is the main activator of frq transcription by binding to its promoter. Here, we show that Centromere Binding Factor 1 (CBF-1) is a critical component of the N. crassa circadian clock by regulating frq transcription. Deletion of cbf-1 resulted in long period and low amplitude rhythms, whereas overexpression of CBF-1 abolished the circadian rhythms. Loss of CBF-1 resulted in WC-independent FRQ expression and suppression of WCC activity. As WCC, CBF-1 also binds to the C-box at the frq promoter. Overexpression of CBF-1 impaired WCC binding to the C-box to suppress frq transcription. Together, our results suggest that the proper level of CBF-1 is critical for circadian clock function by suppressing WC-independent FRQ expression and by regulating WCC binding to the frq promoter.

Circadian clocks, which measure time on a scale of approximately 24 hours, are generated by a cell-autonomous circadian oscillator comprised of autoregulatory feedback loops. In the Neurospora crassa circadian oscillator, WHITE COLLAR complex (WCC) actives transcription of the frequency (frq) gene. FRQ inhibits the activity of WCC to close the negative feedback loop. Here, we showed that the transcription factor CBF-1 functions as a repressor to modulate WCC recruitment to the C-box of frq promoter. Our data showed that deletion or overexpression of CBF-1 dampened circadian rhythm due to impaired WCC binding at the frq promoter. As CBF-1 is conserved in eukaryotes, our data provide novel insights into the negative feedback mechanism that controls the biological clocks in other organisms.

Systematic Analysis of Mouse Genome Reveals Distinct Evolutionary and Functional Properties Among Circadian and Ultradian Genes

In living organisms, biological clocks regulate 24 h (circadian) molecular, physiological, and behavioral rhythms to maintain homeostasis and synchrony with predictable environmental changes, in particular with those induced by Earth’s rotation on its axis. Harmonics of these circadian rhythms having periods of 8 and 12 h (ultradian) have been documented in several species. In mouse liver, harmonics of the 24-h period of gene transcription hallmarked genes oscillating with a frequency two or three times faster than circadian periodicity. Many of these harmonic transcripts enriched pathways regulating responses to environmental stress and coinciding preferentially with subjective dawn and dusk. At this time, the evolutionary history of genes with rhythmic expression is still poorly known and the role of length-of-day changes due to Earth’s rotation speed decrease over the last four billion years is totally ignored. We hypothesized that ultradian and stress anticipatory genes would be more evolutionarily conserved than circadian genes and background non-oscillating genes. To investigate this issue, we performed broad computational analyses of genes/proteins oscillating at different frequency ranges across several species and showed that ultradian genes/proteins, especially those oscillating with a 12-h periodicity, are more likely to be of ancient origin and essential in mice. In summary, our results show that genes with ultradian transcriptional patterns are more likely to be phylogenetically conserved and associated with the primeval and inevitable dawn/dusk transitions.

Breast cancer associated germline structural variants harboring small noncoding RNAs impact post-transcriptional gene regulation

Copy Number Variants (CNVs) are a class of structural variations of DNA. Germline CNVs are known to confer disease susceptibility, but their role in breast cancer warrants further investigations. We hypothesized that breast cancer associated germline CNVs contribute to disease risk through gene dosage or other post-transcriptional regulatory mechanisms, possibly through tissue specific expression of CNV-embedded small-noncoding RNAs (CNV-sncRNAs). Our objectives are to identify breast cancer associated CNVs using a genome wide association study (GWAS), identify sncRNA genes embedded within CNVs, confirm breast tissue (tumor and normal) expression of the sncRNAs, correlate their expression with germline copy status and identify pathways influenced by the genes regulated by sncRNAs. We used an association study design and accessed germline CNV data generated on Affymetrix Human SNP 6.0 array in 686 (in-house data) and 495 (TCGA data) subjects served as discovery and validation cohorts. We identified 1812 breast cancer associated CNVs harboring miRNAs (n = 38), piRNAs (n = 9865), snoRNAs (n = 71) and tRNAs (n = 12) genes. A subset of CNV-sncRNAs expressed in breast tissue, also showed correlation with germline copy status. We identified targets potentially regulated by miRNAs and snoRNAs. In summary, we demonstrate the potential impact of embedded CNV-sncRNAs on expression and regulation of down-stream targets.

Codon usage of highly expressed genes affects proteome-wide translation efficiency

Although the genetic code is redundant, synonymous codons for the same amino acid are not used with equal frequencies in genomes, a phenomenon termed "codon usage bias." Previous studies have demonstrated that synonymous changes in a coding sequence can exert significant cis effects on the gene's expression level. However, whether the codon composition of a gene can also affect the translation efficiency of other genes has not been thoroughly explored. To study how codon usage bias influences the cellular economy of translation, we massively converted abundant codons to their rare synonymous counterpart in several highly expressed genes in Escherichia coli This perturbation reduces both the cellular fitness and the translation efficiency of genes that have high initiation rates and are naturally enriched with the manipulated codon, in agreement with theoretical predictions. Interestingly, we could alleviate the observed phenotypes by increasing the supply of the tRNA for the highly demanded codon, thus demonstrating that the codon usage of highly expressed genes was selected in evolution to maintain the efficiency of global protein translation.

SNP-based susceptibility-resistance association and mRNA expression regulation analyses of tlr7 to grass carp Ctenopharyngodon idella reovirus

Eleven single nucleotide polymorphisms (SNP) in Ctenopharyngodon idella toll-like receptor 7 (citlr7) gene, containing two in the 5'-flanking region, three within the single intron and six distributed in the coding sequence (CDS), were identified. A case-control study of 73 susceptible individuals and 67 resistant individuals was conducted to test the SNPs-based susceptibility-resistance association and mRNA expression of citlr7 to grass carp reovirus (GCRV), showing that both 820 A/G and 1726 A/G were significantly correlative sites in genotype (P < 0·05). Multifactor dimensionality reduction (MDR) analysis suggested the exertion of antiviral effects of 820 A/G might rely on SNPs interactions of citlr7 and C. idella toll-like receptor 8 (citlr8). Combining the mortality rate and citlr7 mRNA expression, it was suggested that 1726 GG-genotyped individuals might be more resistant than 1726 A/G genotyped individuals, indicating the selection on synonymous mutations in 1726 A/G might be susceptibility-resistance-type specific. In addition, haplotype analysis uncovered no significantly correlative haplotypes in citlr7. These findings may provide an in-depth insight for the further functional research of citlr7. The potential genetic markers identified may contribute to the molecular and transgenic breeding of C. idella.

Codon harmonization - going beyond the speed limit for protein expression

Codon usage distribution has been soundly used by nature to fine tune protein biogenesis. Alteration of the mRNA structure or sequential scheduling of codons can profoundly affect translation, thus altering protein yield, functionality, solubility, and proper folding. Building on these observations, here, we present an evaluation of different recently designed algorithms of sequence adaptation based on Codon Adaptation Index (CAI) profiling. The first algorithm globally harmonizes synonymous codons in the original sequence in full respect to the heterologous expression host codon usage. The second recodes the sequence in accordance with the native sequence CAI profile. Our data, generated on three model proteins, highlights the importance to consider gene recoding as a parameter itself for recombinant protein expression improvement.

Origins of the Mechanochemical Coupling of Peptide Bond Formation to Protein Synthesis

Mechanical forces acting on the ribosome can alter the speed of protein synthesis, indicating that mechanochemistry can contribute to translation control of gene expression. The naturally occurring sources of these mechanical forces, the mechanism by which they are transmitted 10 nm to the ribosome's catalytic core, and how they influence peptide bond formation rates are largely unknown. Here, we identify a new source of mechanical force acting on the ribosome by using in situ experimental measurements of changes in nascent-chain extension in the exit tunnel in conjunction with all-atom and coarse-grained computer simulations. We demonstrate that when the number of residues composing a nascent chain increases, its unstructured segments outside the ribosome exit tunnel generate piconewtons of force that are fully transmitted to the ribosome's P-site. The route of force transmission is shown to be through the nascent polypetide's backbone, not through the wall of the ribosome's exit tunnel. Utilizing quantum mechanical calculations we find that a consequence of such a pulling force is to decrease the transition state free energy barrier to peptide bond formation, indicating that the elongation of a nascent chain can accelerate translation. Since nascent protein segments can start out as largely unfolded structural ensembles, these results suggest a pulling force is present during protein synthesis that can modulate translation speed. The mechanism of force transmission we have identified and its consequences for peptide bond formation should be relevant regardless of the source of the pulling force.

Codon usage biases co-evolve with transcription termination machinery to suppress premature cleavage and polyadenylation

Codon usage biases are found in all genomes and influence protein expression levels. The codon usage effect on protein expression was thought to be mainly due to its impact on translation. Here, we show that transcription termination is an important driving force for codon usage bias in eukaryotes. Using Neurospora crassa as a model organism, we demonstrated that introduction of rare codons results in premature transcription termination (PTT) within open reading frames and abolishment of full-length mRNA. PTT is a wide-spread phenomenon in Neurospora, and there is a strong negative correlation between codon usage bias and PTT events. Rare codons lead to the formation of putative poly(A) signals and PTT. A similar role for codon usage bias was also observed in mouse cells. Together, these results suggest that codon usage biases co-evolve with the transcription termination machinery to suppress premature termination of transcription and thus allow for optimal gene expression.

The cell free protein synthesis system from the model filamentous fungus Neurospora crassa

Cell-free protein synthesis (CFPS) can be used in many applications to produce polypeptides and to analyze mechanisms of mRNA translation. Here we describe how to make and use a CPFS system from the model filamentous fungus Neurospora crassa. The extensive genetic resources available in this system provide capacities to exploit robust CFPS for understanding translational control. Included are procedures for the growth and harvesting of cells, the preparation of cell-free extracts that serve as the source of the translational machinery in the CFPS and the preparation of synthetic mRNA to program the CFPS. Methods to accomplish cell-free translation and analyze protein synthesis, and to map positions of ribosomes on mRNAs by toeprinting, are described.

Comparative genomic analysis for nucleotide, codon, and amino acid usage patterns of mycoplasmas

The evolutionary factors in influencing the genetic characteristics of nucleotide, synonymous codon, and amino acid usage of 18 mycoplasma species were analyzed. The nucleotide usage at the 1st and 2nd codon position which determines amino acid composition of proteins has a significant correlation with the total nucleotide composition of gene population of these mycoplasma species, however, the nucleotide usage at the 3rd codon position which affects synonymous codon usage patterns has a slight correlation with either the total nucleotide composition or the nucleotide usage at the 1st and 2nd codon position. Other evolutionary factors join in the evolutionary process of mycoplasma apart from mutation pressure caused by nucleotide usage constraint based on the relationships between effective number of codons/codon adaptation index and nucleotide usage at the 3rd codon position. Although nucleotide usage of gene population in mycoplasma dominates in forming the overall codon usage trends, the relative abundance of codon with nucleotide context and amino acid usage pattern show that translation selection involved in translation accuracy and efficiency play an important role in synonymous codon usage patterns. In addition, synonymous codon usage patterns of gene population have a bigger power to represent genetic diversity among different species than amino acid usage. These results suggest that although the mycoplasmas reduce its genome size during the evolutionary process and shape the form, which is opposite to their hosts, of AT usages at high levels, this kind organism still depends on nucleotide usage at the 1st and 2nd codon positions to control syntheses of the requested proteins for surviving in their hosts and nucleotide usage at the 3rd codon position to develop genetic diversity of different mycoplasma species. This systemic analysis with 18 mycoplasma species may provide useful clues for further in vivo genetic studies on the related species.

Factors affecting synonymous codon usage of housekeeping genes in Drosophila melanogaster

Housekeeping genes (HK genes) are required for cell survival and the maintenance of basic cellular functions. The investigation of factors affecting codon usage patterns in HK genes of insects can help in understanding the molecular evolution of insects and aid the development of insect pest management strategies. In this study, we employed bioinformatics approaches to analyze the codon usage bias (CUB) of HK genes in the insect model organism, Drosophila melanogaster. A comparison of CUB between 1107 HK genes and 1084 high tissue specificity genes suggested that HK genes have higher CUB in D. melanogaster. In addition, we found that CUB inversely correlates with the non-synonymous substitution rate of HK genes. Therefore, we attempted to identify the factors that potentially influence the codon usage pattern of HK genes. Our results suggest that mutation pressure and natural selection highly correlate with CUB in the HK genes of D. melanogaster and that two topological properties of HK proteins (proportion of protein interacting length and protein connectivity) also correlate with CUB in the HK genes of D. melanogaster. This study provides insight into CUB in the HK genes of D. melanogaster, and the results can support future investigations of potential applications in agricultural and biomedical field.

Elongator and codon bias regulate protein levels in mammalian peripheral neurons

Familial dysautonomia (FD) results from mutation in IKBKAP/ELP1, a gene encoding the scaffolding protein for the Elongator complex. This highly conserved complex is required for the translation of codon-biased genes in lower organisms. Here we investigate whether Elongator serves a similar function in mammalian peripheral neurons, the population devastated in FD. Using codon-biased eGFP sensors, and multiplexing of codon usage with transcriptome and proteome analyses of over 6,000 genes, we identify two categories of genes, as well as specific gene identities that depend on Elongator for normal expression. Moreover, we show that multiple genes in the DNA damage repair pathway are codon-biased, and that with Elongator loss, their misregulation is correlated with elevated levels of DNA damage. These findings link Elongator's function in the translation of codon-biased genes with both the developmental and neurodegenerative phenotypes of FD, and also clarify the increased risk of cancer associated with the disease.

Non-equilibrium coupling of protein structure and function to translation-elongation kinetics

Protein folding research has been dominated by the assumption that thermodynamics determines protein structure and function. And that when the folding process is compromised in vivo the proteostasis machinery-chaperones, deaggregases, the proteasome-work to restore proteins to their soluble, functional form or degrade them to maintain the cellular pool of proteins in a quasi-equilibrium state. During the past decade, however, more and more proteins have been identified for which altering only their speed of synthesis alters their structure and function, the efficiency of the down-stream processes they take part in, and cellular phenotype. Indeed, evidence has emerged that evolutionary selection pressures have encoded translation-rate information into mRNA molecules to coordinate diverse co-translational processes. Thus, non-equilibrium physics can play a fundamental role in influencing nascent protein behavior, mRNA sequence evolution, and disease. Here, we discuss how our understanding of this phenomenon is being advanced by the application of theoretical tools from the physical sciences.

Determinants of translation speed are randomly distributed across transcripts resulting in a universal scaling of protein synthesis times

Ribosome profiling experiments have found greater than 100-fold variation in ribosome density along mRNA transcripts, indicating that individual codon elongation rates can vary to a similar degree. This wide range of elongation times, coupled with differences in codon usage between transcripts, suggests that the average codon translation-rate per gene can vary widely. Yet, ribosome run-off experiments have found that the average codon translation rate for different groups of transcripts in mouse stem cells is constant at 5.6 AA/s. How these seemingly contradictory results can be reconciled is the focus of this study. Here, we combine knowledge of the molecular factors shown to influence translation speed with genomic information from Escherichia coli, Saccharomyces cerevisiae and Homo sapiens to simulate the synthesis of cytosolic proteins in these organisms. The model recapitulates a near constant average translation rate, which we demonstrate arises because the molecular determinants of translation speed are distributed nearly randomly amongst most of the transcripts. Consequently, codon translation rates are also randomly distributed and fast-translating segments of a transcript are likely to be offset by equally probable slow-translating segments, resulting in similar average elongation rates for most transcripts. We also show that the codon usage bias does not significantly affect the near random distribution of codon translation rates because only about 10% of the total transcripts in an organism have high codon usage bias while the rest have little to no bias. Analysis of Ribo-Seq data and an in vivo fluorescent assay supports these conclusions.

RNA tales – how embryos read and discard messages from mom

Following fertilization, embryos develop for a substantial amount of time with a transcriptionally silent genome. Thus, early development is maternally programmed, as it solely relies on RNAs and proteins that are provided by the female gamete. However, these maternal instructions are not sufficient to support later steps of embryogenesis and are therefore gradually replaced by novel products synthesized from the zygotic genome. This switch in the origin of molecular players that drive early development is known as the maternal-to-zygotic transition (MZT). MZT is a universal phenomenon among all metazoans and comprises two interconnected processes: maternal mRNA degradation and the transcriptional awakening of the zygotic genome. The recent adaptation of high-throughput methods for use in embryos has deepened our knowledge of the molecular principles underlying MZT. These mechanisms comprise conserved strategies for RNA regulation that operate in many well-studied cellular contexts but that have adapted differently to early development. In this Review, we will discuss advances in our understanding of post-transcriptional regulatory pathways that drive maternal mRNA clearance during MZT, with an emphasis on recent data in zebrafish embryos on codon-mediated mRNA decay, the contributions of microRNAs (miRNAs) and RNA-binding proteins to this process, and the roles of RNA modifications in the stability control of maternal mRNAs.

Synonymous Codon Usage Controls Various Molecular Aspects

Synonymous sites are generally considered to be functionally neutral. However, there are recent contradictory findings suggesting that synonymous alleles might have functional roles in various molecular aspects. For instance, a recent study demonstrated that synonymous single nucleotide polymorphisms have a similar effect size as nonsynonymous single nucleotide polymorphisms in human disease association studies. Researchers have recognized synonymous codon usage bias (SCUB) in the genomes of almost all species and have investigated whether SCUB is due to random nucleotide compositional bias or to natural selection of any functional exposure generated by synonymous mutations. One of the most prominent observations on the non-neutrality of synonymous codons is the correlation between SCUB and levels of gene expression, such that highly expressed genes tend to have a higher preference toward so-called optimal codons than lowly expressed genes. In relation, it is known that amounts of cognate tRNAs that bind to optimal codons are significantly higher than the amounts of cognate tRNAs that bind to non-optimal codons in genomes. In the present paper, we review various functions that synonymous codons might have other than regulating expression levels.

Unique features of nucleotide and codon usage patterns in mycoplasmas revealed by information entropy

Currently, the comparison between GC usage pattern at the 3rd codon position and codon usage index is commonly used to estimate the roles of evolutionary forces in shaping synonymous codon usages, however, this kind of analysis often losses the information about the role of A/T usage bias in shaping synonymous codon usage bias. To overcome this limitation and better understand the interplay between nucleotide and codon usages for the evolution of bacteria at gene levels, in this study, we employed the information entropy method with some modification to estimate roles of nucleotide compositions in the overall codon usage bias for 18 mycoplasma species in combination with Davies-Bouldin index. At gene levels, the overall nucleotide usage bias represents A content as the highest, followed by T, G and C for mycoplasmas, resulting in a low GC content. This feature is universal across these species derived from different hosts, suggesting that the hosts have the limited impact on nucleotide usage bias of mycoplasmas. Information entropy and Davies-Bouldin index can better reveal that the nucleotide usage bias at the 3rd codon position is essential in shaping the overall nucleotide bias for all given mycoplasmas except M. pneumoniae M129. Davies-Bouldin index revealed that the 1st and 2nd codon position play more important role in synonymous codon usage bias than that of the 3rd position at gene levels. To our knowledge, this is the first comprehensive investigation into cooperation between nucleotide and codon usages for mycoplasma and extends our knowledge of the mechanisms that contribute to codon usage and evolution of this microorganism.

%MinMax: A versatile tool for calculating and comparing synonymous codon usage and its impact on protein folding

Most amino acids can be encoded by more than one synonymous codon, but these are rarely used with equal frequency. In many coding sequences the usage patterns of rare versus common synonymous codons is nonrandom and under selection. Moreover, synonymous substitutions that alter these patterns can have a substantial impact on the folding efficiency of the encoded protein. This has ignited broad interest in exploring synonymous codon usage patterns. For many protein chemists, biophysicists and structural biologists, the primary motivation for codon analysis is identifying and preserving usage patterns most likely to impact high-yield production of functional proteins. Here we describe the core functions and new features of %MinMax, a codon usage calculator freely available as a web-based portal and downloadable script (http://www.codons.org). %MinMax evaluates the relative usage frequencies of the synonymous codons used to encode a protein sequence of interest and compares these results to a rigorous null model. Crucially, for analyzing codon usage in common host organisms %MinMax requires only the coding sequence as input; with a user-input codon frequency table, %MinMax can be used to evaluate synonymous codon usage patterns for any coding sequence from any fully sequenced genome. %MinMax makes no assumptions regarding the impact of transfer ribonucleic acid concentrations or other molecular-level interactions on translation rates, yet its output is sufficient to predict the effects of synonymous codon substitutions on cotranslational folding mechanisms. A simple calculation included within %MinMax can be used to harmonize codon usage frequencies for heterologous gene expression.

Codon usage regulates protein structure and function by affecting translation elongation speed in Drosophila cells

Codon usage biases are found in all eukaryotic and prokaryotic genomes and have been proposed to regulate different aspects of translation process. Codon optimality has been shown to regulate translation elongation speed in fungal systems, but its effect on translation elongation speed in animal systems is not clear. In this study, we used a Drosophila cell-free translation system to directly compare the velocity of mRNA translation elongation. Our results demonstrate that optimal synonymous codons speed up translation elongation while non-optimal codons slow down translation. In addition, codon usage regulates ribosome movement and stalling on mRNA during translation. Finally, we show that codon usage affects protein structure and function in vitro and in Drosophila cells. Together, these results suggest that the effect of codon usage on translation elongation speed is a conserved mechanism from fungi to animals that can affect protein folding in eukaryotic organisms.

Codon optimality, bias and usage in translation and mRNA decay

The advent of ribosome profiling and other tools to probe mRNA translation has revealed that codon bias - the uneven use of synonymous codons in the transcriptome - serves as a secondary genetic code: a code that guides the efficiency of protein production, the fidelity of translation and the metabolism of mRNAs. Recent advancements in our understanding of mRNA decay have revealed a tight coupling between ribosome dynamics and the stability of mRNA transcripts; this coupling integrates codon bias into the concept of codon optimality, or the effects that specific codons and tRNA concentrations have on the efficiency and fidelity of the translation machinery. In this Review, we first discuss the evidence for codon-dependent effects on translation, beginning with the basic mechanisms through which translation perturbation can affect translation efficiency, protein folding and transcript stability. We then discuss how codon effects are leveraged by the cell to tailor the proteome to maintain homeostasis, execute specific gene expression programmes of growth or differentiation and optimize the efficiency of protein production.

Evidence of evolutionary selection for cotranslational folding

Recent experiments and simulations have demonstrated that proteins can fold on the ribosome. However, the extent and generality of fitness effects resulting from cotranslational folding remain open questions. Here we report a genome-wide analysis that uncovers evidence of evolutionary selection for cotranslational folding. We describe a robust statistical approach to identify loci within genes that are both significantly enriched in slowly translated codons and evolutionarily conserved. Surprisingly, we find that domain boundaries can explain only a small fraction of these conserved loci. Instead, we propose that regions enriched in slowly translated codons are associated with cotranslational folding intermediates, which may be smaller than a single domain. We show that the intermediates predicted by a native-centric model of cotranslational folding account for the majority of these loci across more than 500 Escherichia coli proteins. By making a direct connection to protein folding, this analysis provides strong evidence that many synonymous substitutions have been selected to optimize translation rates at specific locations within genes. More generally, our results indicate that kinetics, and not just thermodynamics, can significantly alter the efficiency of self-assembly in a biological context.

Reduced Protein Expression in a Virus Attenuated by Codon Deoptimization

A general means of viral attenuation involves the extensive recoding of synonymous codons in the viral genome. The mechanistic underpinnings of this approach remain unclear, however. Using quantitative proteomics and RNA sequencing, we explore the molecular basis of attenuation in a strain of bacteriophage T7 whose major capsid gene was engineered to carry 182 suboptimal codons. We do not detect transcriptional effects from recoding. Proteomic observations reveal that translation is halved for the recoded major capsid gene, and a more modest reduction applies to several coexpressed downstream genes. We observe no changes in protein abundances of other coexpressed genes that are encoded upstream. Viral burst size, like capsid protein abundance, is also decreased by half. Together, these observations suggest that, in this virus, reduced translation of an essential polycistronic transcript and diminished virion assembly form the molecular basis of attenuation.

Predicting synonymous codon usage and optimizing the heterologous gene for expression in E. coli

Of the 20 common amino acids, 18 are encoded by multiple synonymous codons. These synonymous codons are not redundant; in fact, all of codons contribute substantially to protein expression, structure and function. In this study, the codon usage pattern of genes in the E. coli was learned from the sequenced genomes of E. coli. A machine learning based method, Presyncodon was proposed to predict synonymous codon selection in E. coli based on the learned codon usage patterns of the residue in the context of the specific fragment. The predicting results indicate that Presycoden could be used to predict synonymous codon selection of the gene in the E. coli with the high accuracy. Two reporter genes (egfp and mApple) were designed with a combination of low- and high-frequency-usage codons by the method. The fluorescence intensity of eGFP and mApple expressed by the (egfp and mApple) designed by this method was about 2.3- or 1.7- folds greater than that from the genes with only high-frequency-usage codons in E. coli. Therefore, both low- and high-frequency-usage codons make positive contributions to the functional expression of the heterologous proteins. This method could be used to design synthetic genes for heterologous gene expression in biotechnology.

Increasing Protein Production Rates Can Decrease the Rate at Which Functional Protein Is Produced and Their Steady-State Levels

The rate at which soluble, functional protein is produced by the ribosome has recently been found to vary in complex and unexplained ways as various translation-associated rates are altered through synonymous codon substitutions. To understand this phenomenon, here, we combine a well-established ribosome-traffic model with a master-equation model of cotranslational domain folding to explore the scenarios that are possible for the protein production rate, J, and the functional-nascent protein production rate, F, as the rates of various translation processes are altered for five different E. coli proteins. We find that while J monotonically increases as the rates of translation-initiation, -elongation, and -termination increase, F can either increase or decrease. We show that F's nonmonotonic behavior arises within the model from two opposing trends: the tendency for increased translation rates to produce more total protein but less cotranslationally folded protein. We further demonstrate that under certain conditions these nonmonotonic changes in F can result in nonmonotonic variations in post-translational, steady-state levels of functional protein. These results provide a potential explanation for recent experimental observations in which the specific activity of enzymatic proteins decreased with increased synthesis rates. Additionally our model has the potential to be used to rationally design transcripts to maximize the production of functional nascent protein by simultaneously optimizing translation initiation, elongation, and termination rates.

DNA Replication Is Required for Circadian Clock Function by Regulating Rhythmic Nucleosome Composition

Although the coupling between circadian and cell cycles allows circadian clocks to gate cell division and DNA replication in many organisms, circadian clocks were thought to function independently of cell cycle. Here, we show that DNA replication is required for circadian clock function in Neurospora. Genetic and pharmacological inhibition of DNA replication abolished both overt and molecular rhythmicities by repressing frequency (frq) gene transcription. DNA replication is essential for the rhythmic changes of nucleosome composition at the frq promoter. The FACT complex, known to be involved in histone disassembly/reassembly, is required for clock function and is recruited to the frq promoter in a replication-dependent manner to promote replacement of histone H2A.Z by H2A. Finally, deletion of H2A.Z uncoupled the dependence of the circadian clock on DNA replication. Together, these results establish circadian clock and cell cycle as interdependent coupled oscillators and identify DNA replication as a critical process in the circadian mechanism.

Up-Frameshift Protein UPF1 Regulates Neurospora crassa Circadian and Diurnal Growth Rhythms

Nonsense-mediated RNA decay (NMD) is a crucial post-transcriptional regulatory mechanism that recognizes and eliminates aberrantly processed transcripts, and mediates the expression of normal gene transcripts. In this study, we report that in the filamentous fungus Neurospora crassa, the NMD factors play a conserved role in regulating the surveillance of NMD targets including premature termination codon (PTC)-containing transcripts and normal transcripts. The circadian rhythms in all of the knockout strains of upf1-3 genes, which encode the Up-frameshift proteins, were aberrant. The upf1 knockout strain displays a shortened circadian period, which can be restored by constantly expressing exogenous Up-frameshift protein 1 (UPF1). UPF1 regulates the circadian clock by modulating the splicing of the core clock gene frequency (frq) through spliceosome and spliceosome-related arginine/serine-rich splicing factors, which partly account for the short periods in the upf1 knockout strain. We also demonstrated that the clock genes including White Collar (WC)-1, WC-2, and FRQ are involved in controlling the diurnal growth rhythm, and UPF1 may affect the growth rhythms by mediating the FRQ protein levels in the daytime. These findings suggest that the NMD factors play important roles in regulating the circadian clock and diurnal growth rhythms in Neurospora.

Alteration of protein function by a silent polymorphism linked to tRNA abundance

Synonymous single nucleotide polymorphisms (sSNPs) are considered neutral for protein function, as by definition they exchange only codons, not amino acids. We identified an sSNP that modifies the local translation speed of the cystic fibrosis transmembrane conductance regulator (CFTR), leading to detrimental changes to protein stability and function. This sSNP introduces a codon pairing to a low-abundance tRNA that is particularly rare in human bronchial epithelia, but not in other human tissues, suggesting tissue-specific effects of this sSNP. Up-regulation of the tRNA cognate to the mutated codon counteracts the effects of the sSNP and rescues protein conformation and function. Our results highlight the wide-ranging impact of sSNPs, which invert the programmed local speed of mRNA translation and provide direct evidence for the central role of cellular tRNA levels in mediating the actions of sSNPs in a tissue-specific manner.

Synonymous single nucleotide polymorphisms (sSNPs) occur at high frequency in the human genome and are associated with ~50 diseases in humans; the responsible molecular mechanisms remain enigmatic. Here, we investigate the impact of the common sSNP, T2562G, on cystic fibrosis transmembrane conductance regulator (CFTR). Although this sSNP, by itself, does not cause cystic fibrosis (CF), it is prevalent in patients with CFTR-related disorders. T2562G sSNP modifies the local translation speed at the Thr854 codon, leading to changes in CFTR stability and channel function. This sSNP introduces a codon pairing to a low-abundance tRNA, which is particularly rare in human bronchial epithelia, but not in other human tissues, suggesting a tissue-specific effect of this sSNP. Enhancement of the cellular concentration of the tRNA cognate to the mutant ACG codon rescues the stability and conduction defects of T2562G-CFTR. These findings reveal an unanticipated mechanism—inverting the programmed local speed of mRNA translation in a tRNA-dependent manner—for sSNP-associated diseases.

Codon usage affects the structure and function of the Drosophila circadian clock protein PERIOD

Fu et al. show that Drosophila period (dper) codon usage is important for circadian clock function. Codon optimization of dper resulted in conformational changes of dPER protein, altered dPER phosphorylation profile and stability, and impaired dPER function in the circadian negative feedback loop, which manifests into changes in molecular rhythmicity and abnormal circadian behavioral output.

Codon usage bias is a universal feature of all genomes, but its in vivo biological functions in animal systems are not clear. To investigate the in vivo role of codon usage in animals, we took advantage of the sensitivity and robustness of the Drosophila circadian system. By codon-optimizing parts of Drosophila period (dper), a core clock gene that encodes a critical component of the circadian oscillator, we showed that dper codon usage is important for circadian clock function. Codon optimization of dper resulted in conformational changes of the dPER protein, altered dPER phosphorylation profile and stability, and impaired dPER function in the circadian negative feedback loop, which manifests into changes in molecular rhythmicity and abnormal circadian behavioral output. This study provides an in vivo example that demonstrates the role of codon usage in determining protein structure and function in an animal system. These results suggest a universal mechanism in eukaryotes that uses a codon usage “code” within genetic codons to regulate cotranslational protein folding.

Widespread position-specific conservation of synonymous rare codons within coding sequences

Synonymous rare codons are considered to be sub-optimal for gene expression because they are translated more slowly than common codons. Yet surprisingly, many protein coding sequences include large clusters of synonymous rare codons. Rare codons at the 5’ terminus of coding sequences have been shown to increase translational efficiency. Although a general functional role for synonymous rare codons farther within coding sequences has not yet been established, several recent reports have identified rare-to-common synonymous codon substitutions that impair folding of the encoded protein. Here we test the hypothesis that although the usage frequencies of synonymous codons change from organism to organism, codon rarity will be conserved at specific positions in a set of homologous coding sequences, for example to tune translation rate without altering a protein sequence. Such conservation of rarity–rather than specific codon identity–could coordinate co-translational folding of the encoded protein. We demonstrate that many rare codon cluster positions are indeed conserved within homologous coding sequences across diverse eukaryotic, bacterial, and archaeal species, suggesting they result from positive selection and have a functional role. Most conserved rare codon clusters occur within rather than between conserved protein domains, challenging the view that their primary function is to facilitate co-translational folding after synthesis of an autonomous structural unit. Instead, many conserved rare codon clusters separate smaller protein structural motifs within structural domains. These smaller motifs typically fold faster than an entire domain, on a time scale more consistent with translation rate modulation by synonymous codon usage. While proteins with conserved rare codon clusters are structurally and functionally diverse, they are enriched in functions associated with organism growth and development, suggesting an important role for synonymous codon usage in organism physiology. The identification of conserved rare codon clusters advances our understanding of distinct, functional roles for otherwise synonymous codons and enables experimental testing of the impact of synonymous codon usage on the production of functional proteins.

Proteins are long linear polymers that must fold into complex three-dimensional shapes in order to carry out their cellular functions. Every protein is synthesized by the ribosome, which decodes each trinucleotide codon in an mRNA coding sequence in order to select the amino acid residue that will occupy each position in the protein sequence. Most amino acids can be encoded by more than one codon, but these synonymous codons are not used with equal frequency. Rare codons are associated with generally slower rates for protein synthesis, and for this reason have traditionally been considered mildly deleterious for efficient protein production. However, because synonymous codon substitutions do not change the sequence of the encoded protein, the majority view is that they merely reflect genomic ‘background noise’. To the contrary, here we show that the positions of many synonymous rare codons are conserved in mRNA sequences that encode structurally similar proteins from a diverse range of organisms. These results suggest that rare codons have a functional role related to the production of functional proteins, potentially to regulate the rate of protein synthesis and the earliest steps of protein folding, while synthesis is still underway.

Effect of Nascent Peptide Steric Bulk on Elongation Kinetics in the Ribosome Exit Tunnel

All proteins are synthesized by the ribosome, a macromolecular complex that accomplishes the life-sustaining tasks of faithfully decoding mRNA and catalyzing peptide bond formation at the peptidyl transferase center (PTC). The ribosome has evolved an exit tunnel to host the elongating new peptide, protect it from proteolytic digestion, and guide its emergence. It is here that the nascent chain begins to fold. This folding process depends on the rate of translation at the PTC. We report here that besides PTC events, translation kinetics depend on steric constraints on nascent peptide side chains and that confined movements of cramped side chains within and through the tunnel fine-tune elongation rates.

Determinants of Helix Formation for a Kv1.3 Transmembrane Segment inside the Ribosome Exit Tunnel

Proteins begin to fold in the ribosome, and misfolding has pathological consequences. Among the earliest folding events in biogenesis is the formation of a helix, an elementary structure that is ubiquitously present and required for correct protein folding in all proteomes. The determinants underlying helix formation in the confined space of the ribosome exit tunnel are relatively unknown. We chose the second transmembrane segment, S2, of a voltage-gated potassium channel, Kv1.3, as a model to probe this issue. Since the N terminus of S2 is initially in an extended conformation in the folding vestibule of the ribosome yet ultimately emerges at the exit port as a helix, S2 is ideally suited for delineating sequential events and folding determinants of helix formation inside the ribosome. We show that S2's extended N terminus inside the tunnel is converted into a helix by a single, distant mutation in the nascent peptide. This transition depends on nascent peptide sequence at specific tunnel locations. Co-translational secondary folding of nascent chains inside the ribosome has profound physiological consequences that bear on correct membrane insertion, tertiary folding, oligomerization, and biochemical modification of the newborn protein during biogenesis.

Making Time: Conservation of Biological Clocks from Fungi to Animals

The capacity for biological timekeeping arose at least three times through evolution, in prokaryotic cyanobacteria, in cells that evolved into higher plants, and within the group of organisms that eventually became the fungi and the animals. Neurospora is a tractable model system for understanding the molecular bases of circadian rhythms in the last of these groups, and is perhaps the most intensively studied circadian cell type. Rhythmic processes described in fungi include growth rate, stress responses, developmental capacity, and sporulation, as well as much of metabolism; fungi use clocks to anticipate daily environmental changes. A negative feedback loop comprises the core of the circadian system in fungi and animals. In Neurospora, the best studied fungal model, it is driven by two transcription factors, WC-1 and WC-2, that form the White Collar Complex (WCC). WCC elicits expression of the frq gene. FRQ complexes with other proteins, physically interacts with the WCC, and reduces its activity; the kinetics of these processes is strongly influenced by progressive phosphorylation of FRQ. When FRQ becomes sufficiently phosphorylated that it loses the ability to influence WCC activity, the circadian cycle starts again. Environmental cycles of light and temperature influence frq and FRQ expression and thereby reset the internal circadian clocks. The molecular basis of circadian output is also becoming understood. Taken together, molecular explanations are emerging for all the canonical circadian properties, providing a molecular and regulatory framework that may be extended to many members of the fungal and animal kingdoms, including humans.

Transcriptome Analysis of Core Dinoflagellates Reveals a Universal Bias towards "GC" Rich Codons

Although dinoflagellates are a potential source of pharmaceuticals and natural products, the mechanisms for regulating and producing these compounds are largely unknown because of extensive post-transcriptional control of gene expression. One well-documented mechanism for controlling gene expression during translation is codon bias, whereby specific codons slow or even terminate protein synthesis. Approximately 10,000 annotatable genes from fifteen “core” dinoflagellate transcriptomes along a range of overall guanine and cytosine (GC) content were used for codonW analysis to determine the relative synonymous codon usage (RSCU) and the GC content at each codon position. GC bias in the analyzed dataset and at the third codon position varied from 51% and 54% to 66% and 88%, respectively. Codons poor in GC were observed to be universally absent, but bias was most pronounced for codons ending in uracil followed by adenine (UA). GC bias at the third codon position was able to explain low abundance codons as well as the low effective number of codons. Thus, we propose that a bias towards codons rich in GC bases is a universal feature of core dinoflagellates, possibly relating to their unique chromosome structure, and not likely a major mechanism for controlling gene expression.

When mRNA translation meets decay

Messenger RNA (mRNA) translation and mRNA degradation are important determinants of protein output, and they are interconnected. Previously, it was thought that translation of an mRNA, as a rule, prevents its degradation. mRNA surveillance mechanisms, which degrade mRNAs as a consequence of their translation, were considered to be exceptions to this rule. Recently, however, it has become clear that many mRNAs are degraded co-translationally, and it has emerged that codon choice, by influencing the rate of ribosome elongation, affects the rate of mRNA decay. In this review, we discuss the links between translation and mRNA stability, with an emphasis on emerging data suggesting that codon optimality may regulate mRNA degradation.

Bicodon bias can determine the role of synonymous SNPs in human diseases

Background

For a long time synonymous single nucleotide polymorphisms were considered as silent mutations. However, nowadays it is well known that they can affect protein conformation and function, leading to altered disease susceptibilities, differential prognosis and/or drug responses, among other clinically relevant genetic traits. This occurs through different mechanisms: by disrupting the splicing signals of precursor mRNAs, affecting regulatory binding-sites of transcription factors and miRNAs, or by modifying the secondary structure of mRNAs.

Results

In this paper we considered 22 human genetic diseases or traits, linked to 35 synonymous single nucleotide polymorphisms in 27 different genes. We performed a local sequence context analysis in terms of the ribosomal pause propensity affected by synonymous single nucleotide polymorphisms. We found that synonymous mutations related to the above mentioned mechanisms presented small pause propensity changes, whereas synonymous mutations that were not related to those mechanisms presented large pause propensity changes. On the other hand, we did not observe large variations in the codon usage of codons associated with these mutations. Furthermore, we showed that the changes in the pause propensity associated with benign sSNPs are significantly lower than the pause propensity changes related to sSNPs associated to diseases.

Conclusions

These results suggest that the genetic diseases or traits related to synonymous mutations with large pause propensity changes, could be the consequence of another mechanism underlying non-silent synonymous mutations. Namely, alternative protein configuration related, in turn, to alterations in the ribosome-mediated translational attenuation program encoded by pairs of consecutive codons, not codons. These findings shed light on the latter mechanism based on the perturbation of the co-translational folding process.

Electronic supplementary material

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

Synonymous Codons: Choose Wisely for Expression

The genetic code, which defines the amino acid sequence of a protein, also contains information that influences the rate and efficiency of translation. Neither the mechanisms nor functions of codon-mediated regulation were well understood. The prevailing model was that the slow translation of codons decoded by rare tRNAs reduces efficiency. Recent genome-wide analyses have clarified several issues. Specific codons and codon combinations modulate ribosome speed and facilitate protein folding. However, tRNA availability is not the sole determinant of rate; rather, interactions between adjacent codons and wobble base pairing are key. One mechanism linking translation efficiency and codon use is that slower decoding is coupled to reduced mRNA stability. Changes in tRNA supply mediate biological regulationfor instance,, changes in tRNA amounts facilitate cancer metastasis.

Modulation of the Proteostasis Machinery to Overcome Stress Caused by Diminished Levels of t6A-Modified tRNAs in Drosophila

Transfer RNAs (tRNAs) harbor a subset of post-transcriptional modifications required for structural stability or decoding function. N⁶-threonylcarbamoyladenosine (t⁶A) is a universally conserved modification found at position 37 in tRNA that pair A-starting codons (ANN) and is required for proper translation initiation and to prevent frame shift during elongation. In its absence, the synthesis of aberrant proteins is likely, evidenced by the formation of protein aggregates. In this work, our aim was to study the relationship between t⁶A-modified tRNAs and protein synthesis homeostasis machinery using Drosophila melanogaster. We used the Gal4/UAS system to knockdown genes required for t⁶A synthesis in a tissue and time specific manner and in vivo reporters of unfolded protein response (UPR) activation. Our results suggest that t⁶A-modified tRNAs, synthetized by the threonyl-carbamoyl transferase complex (TCTC), are required for organismal growth and imaginal cell survival, and is most likely to support proper protein synthesis.

Distribution bias analysis of germline and somatic single-nucleotide variations that impact protein functional site and neighboring amino acids

Single nucleotide variations (SNVs) can result in loss or gain of protein functional sites. We analyzed the effects of SNVs on enzyme active sites, ligand binding sites, and various types of post translational modification (PTM) sites. We found that, for most types of protein functional sites, the SNV pattern differs between germline and somatic mutations as well as between synonymous and non-synonymous mutations. From a total of 51,138 protein functional site affecting SNVs (pfsSNVs), a pan-cancer analysis revealed 142 somatic pfsSNVs in five or more cancer types. By leveraging patient information for somatic pfsSNVs, we identified 17 loss of functional site SNVs and 60 gain of functional site SNVs which are significantly enriched in patients with specific cancer types. Of the key pfsSNVs identified in our analysis above, we highlight 132 key pfsSNVs within 17 genes that are found in well-established cancer associated gene lists. For illustrating how key pfsSNVs can be prioritized further, we provide a use case where we performed survival analysis showing that a loss of phosphorylation site pfsSNV at position 105 in MEF2A is significantly associated with decreased pancreatic cancer patient survival rate. These 132 pfsSNVs can be used in developing genetic testing pipelines.

Systems Biology-Derived Discoveries of Intrinsic Clocks

A systems approach to studying biology uses a variety of mathematical, computational, and engineering tools to holistically understand and model properties of cells, tissues, and organisms. Building from early biochemical, genetic, and physiological studies, systems biology became established through the development of genome-wide methods, high-throughput procedures, modern computational processing power, and bioinformatics. Here, we highlight a variety of systems approaches to the study of biological rhythms that occur with a 24-h period-circadian rhythms. We review how systems methods have helped to elucidate complex behaviors of the circadian clock including temperature compensation, rhythmicity, and robustness. Finally, we explain the contribution of systems biology to the transcription-translation feedback loop and posttranslational oscillator models of circadian rhythms and describe new technologies and "-omics" approaches to understand circadian timekeeping and neurophysiology.

Differential Regulation of the Melanoma Proteome by eIF4A1 and eIF4E

Small molecules and antisense oligonucleotides that inhibit the translation initiation factors eIF4A1 and eIF4E have been explored as broad-based therapeutic agents for cancer treatment, based on the frequent upregulation of these two subunits of the eIF4F cap-binding complex in many cancer cells. Here, we provide support for these therapeutic approaches with mechanistic studies of eIF4F-driven tumor progression in a preclinical model of melanoma. Silencing eIF4A1 or eIF4E decreases melanoma proliferation and invasion. There were common effects on the level of cell-cycle proteins that could explain the antiproliferative effects in vitro Using clinical specimens, we correlate the common cell-cycle targets of eIF4A1 and eIF4E with patient survival. Finally, comparative proteomic and transcriptomic analyses reveal extensive mechanistic divergence in response to eIF4A1 or eIF4E silencing. Current models indicate that eIF4A1 and eIF4E function together through the 5'UTR to increase translation of oncogenes. In contrast, our data demonstrate that the common effects of eIF4A1 and eIF4E on translation are mediated by the coding region and 3'UTR. Moreover, their divergent effects occur through the 5'UTR. Overall, our work shows that it will be important to evaluate subunit-specific inhibitors of eIF4F in different disease contexts to fully understand their anticancer actions. Cancer Res; 77(3); 613-22. ©2016 AACR.

Analysis of codon usage patterns in Ginkgo biloba reveals codon usage tendency from A/U-ending to G/C-ending

As one of the most ancient tree species, the codon usage pattern analysis of Ginkgo biloba is a useful way to understand its evolutionary and genetic mechanisms. Several studies have been conducted on angiosperms, but seldom on gymnosperms. Based on RNA-Seq data of the G. biloba transcriptome, amount to 17,579 unigenes longer than 300 bp were selected and analyzed from 68,547 candidates. The codon usage pattern tended towards more frequently use of A/U-ending codons, which showed an obvious gradient progressing from gymnosperms to dicots to monocots. Meanwhile, analysis of high/low-expression unigenes revealed that high-expression unigenes tended to use G/C-ending codons together with more codon usage bias. Variation of unigenes with different functions suggested that unigenes involving in environment adaptation use G/C-ending codons more frequently with more usage bias, and these results were consistent with the conclusion that the formation of G. biloba codon usage bias was dominated by natural selection.

Improving virus production through quasispecies genomic selection and molecular breeding

Virus production still is a challenging issue in antigen manufacture, particularly with slow-growing viruses. Deep-sequencing of genomic regions indicative of efficient replication may be used to identify high-fitness minority individuals suppressed by the ensemble of mutants in a virus quasispecies. Molecular breeding of quasispecies containing colonizer individuals, under regimes allowing more than one replicative cycle, is a strategy to select the fittest competitors among the colonizers. A slow-growing cell culture-adapted hepatitis A virus strain was employed as a model for this strategy. Using genomic selection in two regions predictive of efficient translation, the internal ribosome entry site and the VP1-coding region, high-fitness minority colonizer individuals were identified in a population adapted to conditions of artificially-induced cellular transcription shut-off. Molecular breeding of this population with a second one, also adapted to transcription shut-off and showing an overall colonizer phenotype, allowed the selection of a fast-growing population of great biotechnological potential.

The Need for Integrated Approaches in Metabolic Engineering

This review highlights state-of-the-art procedures for heterologous small-molecule biosynthesis, the associated bottlenecks, and new strategies that have the potential to accelerate future accomplishments in metabolic engineering. We emphasize that a combination of different approaches over multiple time and size scales must be considered for successful pathway engineering in a heterologous host. We have classified these optimization procedures based on the "system" that is being manipulated: transcriptome, translatome, proteome, or reactome. By bridging multiple disciplines, including molecular biology, biochemistry, biophysics, and computational sciences, we can create an integral framework for the discovery and implementation of novel biosynthetic production routes.

Codon optimality controls differential mRNA translation during amino acid starvation

It is common wisdom that codon usage bias has evolved in the selection for efficient translation, in which highly expressed genes are encoded predominantly by optimal codons. However, a growing body of evidence suggests regulatory roles for non-optimal codons in translation dynamics. Here we report that in mammalian cells, non-optimal codons play a critical role in promoting selective mRNA translation during amino acid starvation. During starvation, in contrast to genes encoding ribosomal proteins whose translation is highly sensitive to amino acid deprivation, translation of genes involved in the cellular protein degradation pathways remains unaffected. We found that these two gene groups bear different codon composition, with non-optimal codons being highly enriched in genes encoding the ubiquitin-proteasome system. Supporting the selective tRNA charging model originally proposed in Escherichia coli, we demonstrated that tRNA isoacceptors decoding rare codons are maintained in translating ribosomes under amino acid starvation. Finally, using luciferase reporters fused with endogenous gene-derived sequences, we show that codon optimality contributes to differential mRNA translation in response to amino acid starvation. These results highlight the physiological significance of codon usage bias in cellular adaptation to stress.

Single Synonymous Mutations in KRAS Cause Transformed Phenotypes in NIH3T3 Cells

Synonymous mutations in the KRAS gene are clustered at G12, G13, and G60 in human cancers. We constructed 9 stable NIH3T3 cell lines expressing KRAS, each with one of these synonymous mutations. Compared to the negative control cell line expressing the wild type human KRAS gene, all the synonymous mutant lines expressed more KRAS protein, grew more rapidly and to higher densities, and were more invasive in multiple assays. Three of the cell lines showed dramatic loss of contact inhibition, were more refractile under phase contrast, and their refractility was greatly reduced by treatment with trametinib. Codon usage at these glycines is highly conserved in KRAS compared to HRAS, indicating selective pressure. These transformed phenotypes suggest that synonymous mutations found in driver genes such as KRAS may play a role in human cancers.

Codon usage is an important determinant of gene expression levels largely through its effects on transcription

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.