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

Translation Rates and Protein Folding

The mRNA coding sequence defines not only the amino acid sequence of the protein, but also the speed at which the ribosomes move along the mRNA while making the protein. The non-uniform local kinetics - denoted as translational rhythm - is similar among mRNAs coding for related protein folds. Deviations from this conserved rhythm can result in protein misfolding. In this review we summarize the experimental evidence demonstrating how local translation rates affect cotranslational protein folding, with the focus on the synonymous codons and patches of charged residues in the nascent peptide as best-studied examples. Alterations in nascent protein conformations due to disturbed translational rhythm can persist off the ribosome, as demonstrated by the effects of synonymous codon variants of several disease-related proteins. Charged amino acid patches in nascent chains also modulate translation and cotranslational protein folding, and can abrogate translation when placed at the N-terminus of the nascent peptide. During cotranslational folding, incomplete nascent chains navigate through a unique conformational landscape in which earlier intermediate states become inaccessible as the nascent peptide grows. Precisely tuned local translation rates, as well as interactions with the ribosome, guide the folding pathway towards the native structure, whereas deviations from the natural translation rhythm may favor pathways leading to trapped misfolded states. Deciphering the 'folding code' of the mRNA will contribute to understanding the diseases caused by protein misfolding and to rational protein design.

Orb2 enables rare-codon-enriched mRNA expression during Drosophila neuron differentiation

Regulation of codon optimality is an increasingly appreciated layer of cell- and tissue-specific protein expression control. Here, we use codon-modified reporters to show that differentiation of Drosophila neural stem cells into neurons enables protein expression from rare-codon-enriched genes. From a candidate screen, we identify the cytoplasmic polyadenylation element binding (CPEB) protein Orb2 as a positive regulator of rare-codon-dependent mRNA stability in neurons. Using RNA sequencing, we reveal that Orb2-upregulated mRNAs in the brain with abundant Orb2 binding sites have a rare-codon bias. From these Orb2-regulated mRNAs, we demonstrate that rare-codon enrichment is important for mRNA stability and social behavior function of the metabotropic glutamate receptor (mGluR). Our findings reveal a molecular mechanism by which neural stem cell differentiation shifts genetic code regulation to enable critical mRNA stability and protein expression.

Alter codon bias of the P. pastoris genome to overcome a bottleneck in codon optimization strategy development and improve protein expression

Apart from its role in translation, codon bias is also an important mechanism to regulate mRNA levels. The traditional frequency-based codon optimization strategy is rather efficient in organisms such as N. crassa, but much less in yeast P. pastoris which is a popular host for heterologous protein expression. This is because that unlike N. crassa, the preferred codons of P. pastoris are actually AU-rich and hence codon optimization for extremely low GC content comes with issues of pre-mature transcriptional termination or low RNA stability in spite of translational advantages. To overcome this bottleneck, we focused on three reporter genes in P. pastoris first and confirmed the great advantage of GC-prone codon optimization on mRNA levels. Then we altered the codon bias profile of P. pastoris by introducing additional rare tRNA gene copies. Prior to that we constructed IPTG-regulated tRNA species to enable chassis cells to switch between different codon bias status. As demonstrated again with reporter genes, protein yield of luc and 0788 was successfully increased by 4-5 folds in chassis cells. In summary, here we provide an alternative codon optimization strategy for genes with unsatisfactory performance under traditional codon frequency-based optimization.

The Function, Regulation, and Mechanism of Protein Turnover in Circadian Systems in Neurospora and Other Species

Circadian clocks drive a large array of physiological and behavioral activities. At the molecular level, circadian clocks are composed of positive and negative elements that form core oscillators generating the basic circadian rhythms. Over the course of the circadian period, circadian negative proteins undergo progressive hyperphosphorylation and eventually degrade, and their stability is finely controlled by complex post-translational pathways, including protein modifications, genetic codon preference, protein-protein interactions, chaperon-dependent conformation maintenance, degradation, etc. The effects of phosphorylation on the stability of circadian clock proteins are crucial for precisely determining protein function and turnover, and it has been proposed that the phosphorylation of core circadian clock proteins is tightly correlated with the circadian period. Nonetheless, recent studies have challenged this view. In this review, we summarize the research progress regarding the function, regulation, and mechanism of protein stability in the circadian clock systems of multiple model organisms, with an emphasis on Neurospora crassa, in which circadian mechanisms have been extensively investigated. Elucidation of the highly complex and dynamic regulation of protein stability in circadian clock networks would greatly benefit the integrated understanding of the function, regulation, and mechanism of protein stability in a wide spectrum of other biological processes.

Mammalian circadian clock proteins form dynamic interacting microbodies distinct from phase separation

Liquid-liquid phase separation (LLPS) underlies diverse biological processes. Because most LLPS studies were performed in vitro using recombinant proteins or in cells that overexpress protein, the physiological relevance of LLPS for endogenous protein is often unclear. PERIOD, the intrinsically disordered domain-rich proteins, are central mammalian circadian clock components and interact with other clock proteins in the core circadian negative feedback loop. Different core clock proteins were previously shown to form large complexes. Circadian clock studies often rely on experiments that overexpress clock proteins. Here, we show that when Per2 transgene was stably expressed in cells, PER2 protein formed nuclear phosphorylation-dependent slow-moving LLPS condensates that recruited other clock proteins. Super-resolution microscopy of endogenous PER2, however, revealed formation of circadian-controlled, rapidly diffusing nuclear microbodies that were resistant to protein concentration changes, hexanediol treatment, and loss of phosphorylation, indicating that they are distinct from the LLPS condensates caused by protein overexpression. Surprisingly, only a small fraction of endogenous PER2 microbodies transiently interact with endogenous BMAL1 and CRY1, a conclusion that was confirmed in cells and in mice tissues, suggesting an enzyme-like mechanism in the circadian negative feedback process. Together, these results demonstrate that the dynamic interactions of core clock proteins are a key feature of mammalian circadian clock mechanism and the importance of examining endogenous proteins in LLPS and circadian clock studies.

Mammalian circadian clock proteins form dynamic interacting microbodies distinct from phase separation

Liquid-liquid phase separation (LLPS) underlies diverse biological processes. Because most LLPS studies were performed in vitro or in cells that overexpress protein, the physiological relevance of LLPS is unclear. PERIOD proteins are central mammalian circadian clock components and interact with other clock proteins in the core circadian negative feedback loop. Different core clock proteins were previously shown to form large complexes. Here we show that when transgene was stably expressed, PER2 formed nuclear phosphorylation-dependent LLPS condensates that recruited other clock proteins. Super-resolution microscopy of endogenous PER2, however, revealed formation of circadian-controlled, rapidly diffusing microbodies that were resistant to protein concentration changes, hexanediol treatment, and loss of phosphorylation, indicating that they are distinct from the LLPS condensates caused by overexpression. Surprisingly, only a small fraction of endogenous PER2 microbodies transiently interact with endogenous BMAL1 and CRY1, a conclusion that was confirmed in cells and in mice tissues, suggesting an enzyme-like mechanism in the circadian negative feedback process. Together, these results demonstrate that the dynamic interactions of core clock proteins is a key feature of mammalian circadian clock mechanism and the importance of examining endogenous proteins in LLPS and circadian studies.

Orb2 enables rare-codon-enriched mRNA expression during Drosophila neuron differentiation

Regulation of codon optimality is an increasingly appreciated layer of cell- and tissue-specific protein expression control. Here, we use codon-modified reporters to show that differentiation of Drosophila neural stem cells into neurons enables protein expression from rare-codon-enriched genes. From a candidate screen, we identify the cytoplasmic polyadenylation element binding (CPEB) protein Orb2 as a positive regulator of rare-codon-dependent expression in neurons. Using RNA sequencing, we reveal that Orb2-upregulated mRNAs in the brain with abundant Orb2 binding sites have a rare-codon bias. From these Orb2-regulated mRNAs, we demonstrate that rare-codon enrichment is important for expression control and social behavior function of the metabotropic glutamate receptor (mGluR). Our findings reveal a molecular mechanism by which neural stem cell differentiation shifts genetic code regulation to enable critical mRNA and protein expression.

How soluble misfolded proteins bypass chaperones at the molecular level

Subpopulations of soluble, misfolded proteins can bypass chaperones within cells. The extent of this phenomenon and how it happens at the molecular level are unknown. Through a meta-analysis of the experimental literature we find that in all quantitative protein refolding studies there is always a subpopulation of soluble but misfolded protein that does not fold in the presence of one or more chaperones, and can take days or longer to do so. Thus, some misfolded subpopulations commonly bypass chaperones. Using multi-scale simulation models we observe that the misfolded structures that bypass various chaperones can do so because their structures are highly native like, leading to a situation where chaperones do not distinguish between the folded and near-native-misfolded states. More broadly, these results provide a mechanism by which long-time scale changes in protein structure and function can persist in cells because some misfolded states can bypass components of the proteostasis machinery.

Phosphorylation Promotes the Accumulation of PERIOD Protein Foci

Circadian clock drives the 24-h rhythm in our behavior and physiology. The molecular clock consists of a series of transcriptional/translational feedback loops operated by a number of clock genes. A very recent study reported that the clock protein PERIOD (PER) is organized into discrete foci at the nuclear envelope in fly circadian neurons, which is believed to be important for controlling the subcellular localization of clock genes. Loss of inner nuclear membrane protein lamin B receptor (LBR) leads to disruption of these foci, but how they are regulated is yet unknown. Here, we found that PER foci are likely phase-separated condensates, the formation of which is mediated by intrinsically disordered region in PER. Phosphorylation promotes the accumulation of these foci. Protein phosphatase 2A, which is known to dephosphorylate PER, hampers the accumulation of the foci. On the other hand, the circadian kinase DOUBLETIME (DBT) which phosphorylates PER enhances the accumulation of the foci. LBR likely facilitates PER foci accumulation by destabilizing the catalytic subunit of protein phosphatase 2A, MICROTUBULE STAR (MTS). In conclusion, here, we demonstrate a key role for phosphorylation in promoting the accumulation of PER foci, while LBR modulates this process by impinging on the circadian phosphatase MTS.

Codon modification of Tuba1a alters mRNA levels and causes a severe neurodevelopmental phenotype in mice

The tubulinopathies are an umbrella of rare diseases that result from mutations in tubulin genes and are frequently characterised by severe brain malformations. The characteristics of a given disease reflect the expression pattern of the transcript, the function of a given tubulin gene, and the role microtubules play in a particular cell type. Mouse models have proved to be valuable tools that have provided insight into the molecular and cellular mechanisms that underlie the disease state. In this manuscript we compare two Tuba1a mouse models, both of which express wild-type TUBA1A protein but employ different codon usage. We show that modification of the Tuba1a mRNA sequence results in homozygous lethality and a severe neurodevelopmental phenotype. This is associated with a decrease in the number of post-mitotic neurons, PAX6 positive progenitors, and an increase in the number of apoptotic cells. We attribute this to a decrease in the stability of the modified Tuba1a transcript, and the absence of compensation by the other neurogenic tubulins. Our findings highlight the importance of maintaining the wild-type coding sequence when engineering mouse lines and the impact of synonymous genetic variation.

The role of spatiotemporal organization and dynamics of clock complexes in circadian regulation

Circadian clocks are cell autonomous timekeepers that regulate ∼24-h oscillations in the expression of many genes and control rhythms in nearly all our behavior and physiology. Almost every cell in the human body has a molecular clock and networks of cells containing clock proteins orchestrate daily rhythms in many physiological processes, from sleep-wake cycles to metabolism to immunity. All eukaryotic circadian clocks are based on transcription-translation delayed negative feedback loops in which activation of core clock genes is negatively regulated by their cognate protein products. Our current understanding of circadian clocks has been accumulated from decades of genetic and biochemical experiments, however, what remains poorly understood is how clock proteins, genes, and mRNAs are spatiotemporally organized within live clock cells and how such subcellular organization affects circadian rhythms at the single cell level. Here, we review recent progress in understanding how clock proteins and genes are spatially organized within clock cells over the circadian cycle and the role of such organization in generating circadian rhythms and highlight open questions for future studies.

The glutamate metabotropic receptor 5 (GRM5) gene is associated with beef cattle home range and movement tortuosity

Background

The grazing behaviour of herbivores and their grazing personalities might in part be determined genetically, but there are few studies in beef cattle illustrating this. In this study, we investigated for first time the genetic variation within a candidate ‘grazing gene’, the glutamate metabotropic receptor 5 gene (GRM5), and tested associations between variation in that gene and variation in grazing personality behaviours (GP-behaviours) displayed by free-ranging cows during winter grazing in the steep and rugged rangelands of New Zealand. Mature beef cows (n = 303, from 3 to 10 years of age) were tracked with global positioning system (GPS) and, with 5-minutes (min) relocation frequency, various GP-behaviours were calculated. These included horizontal and vertical distances travelled, mean elevation, elevation range, elevation gain, slope, home range and movement tortuosity, variously calculated using daily relocation trajectories with repeated measurements (i.e., 7 to 24 days (d)) and satellite-derived digital elevation models (DEM). The different GP-behaviours were fitted into mixed models to ascertain their associations with variant sequences and genotypes of GRM5.

Results

We discovered three GRM5 variants (A, B and C) and identified the six possible genotypes in the cattle studied. The mixed models revealed that A was significantly associated with elevation range, home range and movement tortuosity. Similarly, GRM5 genotypes were associated (P < 0.05) to home range and movement tortuosity, while trends suggesting association (P <  0.1) were also revealed for elevation range and horizontal distance travelled. Most GP-behaviour models were improved by correcting for cow age-class as a fixed factor. The analysis of GP-behaviours averaged per cow age-class suggests that grazing personality is fully established as beef cows reached 4 years of age. Home range and movement tortuosity were not only associated with GRM5 variation, but also negatively correlated with each other (r = − 0.27, P < 0.001).

Conclusions

There seems to be a genetically determined trade-off between home range and movement tortuosity that may be useful in beef cattle breeding programmes aiming to improve the grazing distribution and utilisation of steep and rugged rangelands.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00755-7.

Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less functional

Some misfolded protein conformations can bypass proteostasis machinery and remain soluble in vivo. This is an unexpected observation, as cellular quality control mechanisms should remove misfolded proteins. Three questions, then, are: how do long-lived, soluble, misfolded proteins bypass proteostasis? How widespread are such misfolded states? And how long do they persist? We address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of proteins exhibit misfolded subpopulations that bypass molecular chaperones, avoid aggregation, and will not be rapidly degraded, with some misfolded states persisting for months or longer. The surface properties of these misfolded states are native-like, suggesting they will remain soluble, while self-entanglements make them long-lived kinetic traps. In terms of function, we predict that one-third of proteins can misfold into soluble less-functional states. For the heavily entangled protein glycerol-3-phosphate dehydrogenase, limited-proteolysis mass spectrometry experiments interrogating misfolded conformations of the protein are consistent with the structural changes predicted by our simulations. These results therefore provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass proteostasis, and indicate, unexpectedly, this may be a wide-spread phenomenon.

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Background

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression.

Results

Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to prevous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding.

Conclusions

Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08635-0).

Senecavirus A Enhances Its Adaptive Evolution via Synonymous Codon Bias Evolution

Synonymous codon bias in the viral genome affects protein translation and gene expression, suggesting that the synonymous codon mutant plays an essential role in influencing virulence and evolution. However, how the recessive mutant form contributes to virus evolvability remains elusive. In this paper, we characterize how the Senecavirus A (SVA), a picornavirus, utilizes synonymous codon mutations to influence its evolution, resulting in the adaptive evolution of the virus to adverse environments. The phylogenetic tree and Median-joining (MJ)-Network of these SVA lineages worldwide were constructed to reveal SVA three-stage genetic development clusters. Furthermore, we analyzed the codon bias of the SVA genome of selected strains and found that SVA could increase the GC content of the third base of some amino acid synonymous codons to enhance the viral RNA adaptive evolution. Our results highlight the impact of recessive mutation of virus codon bias on the evolution of the SVA and uncover a previously underappreciated evolutionary strategy for SVA. They also underline the importance of understanding the genetic evolution of SVA and how SVA adapts to the adverse effects of external stress.

Distinct responses to rare codons in select Drosophila tissues

Codon usage bias has long been appreciated to influence protein production. Yet, relatively few studies have analyzed the impacts of codon usage on tissue-specific mRNA and protein expression. Here, we use codon-modified reporters to perform an organism-wide screen in Drosophila melanogaster for distinct tissue responses to codon usage bias. These reporters reveal a cliff-like decline of protein expression near the limit of rare codon usage in endogenously expressed Drosophila genes. Near the edge of this limit, however, we find the testis and brain are uniquely capable of expressing rare codon-enriched reporters. We define a new metric of tissue-specific codon usage, the tissue-apparent Codon Adaptation Index (taCAI), to reveal a conserved enrichment for rare codon usage in the endogenously expressed genes of both Drosophila and human testis. We further demonstrate a role for rare codons in an evolutionarily young testis-specific gene, RpL10Aa. Optimizing RpL10Aa codons disrupts female fertility. Our work highlights distinct responses to rarely used codons in select tissues, revealing a critical role for codon bias in tissue biology.

Large-Scale Analysis of SARS-CoV-2 Synonymous Mutations Reveals the Adaptation to the Human Codon Usage During the Virus Evolution

Many large national and transnational studies have been dedicated to the analysis of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) genome, most of which focused on missense and nonsense mutations. However, approximately 30 per cent of the SARS-CoV-2 variants are synonymous, therefore changing the target codon without affecting the corresponding protein sequence.

By performing a large-scale analysis of sequencing data generated from almost 400,000 SARS-CoV-2 samples, we show that silent mutations increasing the similarity of viral codons to the human ones tend to fixate in the viral genome overtime. This indicates that SARS-CoV-2 codon usage is adapting to the human host, likely improving its effectiveness in using the human aminoacyl-tRNA set through the accumulation of deceitfully neutral silent mutations.

One-Sentence Summary. Synonymous SARS-CoV-2 mutations related to the activity of different mutational processes may positively impact viral evolution by increasing its adaptation to the human codon usage.

Strong negative correlation between codon usage bias and protein structural disorder impedes protein expression after codon optimization

As a common phenomenon existing in almost all genomes, codon usage bias has been studied for a long time. Codon optimization is a frequently used strategy to accelerate protein synthesis rate. Besides regulating protein translation speed, codon usage bias has also be reported to affect co-translation folding and transcription. P. pastoris is a well-developed expression system, whose efficiency is tightly correlated with commercial value. However, few studies focus on the role of codon usage bias in affecting protein expression in P. pastoris. Besides, many genes in P. pastoris genome show significant negative correlation between codon usage bias and protein structural disorder tendency. It's not known whether this feature is important for their expression. In order to answer these questions, we picked 4P. pastoris gene candidates with strong negative correlation between codon usage bias and protein structural disorder. We then performed full-length codon optimization which completely eliminated the correlation. Protein and RNA assays were then used to compare protein and mRNA levels before and after codon optimization. As a result, codon optimization failed to elevate their protein expression levels, and even resulted in a decrease. As represented by the trypsin sensitivity assays, codon optimization also altered the protein structure of 0616 and 0788. Besides protein, codon optimization also affected mRNA levels. Shown by in vitro and in vivo RNA degradation assays, the mRNA stability of 0616, 0788 and 0135 were also altered by codon optimization. For each gene, the detailed effect may be related with its specific sequence and protein structure. Our results suggest that codon usage bias is an important factor to regulate gene expression level, as well as mRNA and protein stabilities in P. pastoris. "Extreme" codon optimization in genes with strong negative correlation between codon usage bias and protein structural disorder tendency may not be favored. Compromised strategies should be tried if expression is not successful. Besides, codon optimization may affect protein structural conformation more severely in structural disordered proteins.

Codon usage and protein length-dependent feedback from translation elongation regulates translation initiation and elongation speed

Essential cellular functions require efficient production of many large proteins but synthesis of large proteins encounters many obstacles in cells. Translational control is mostly known to be regulated at the initiation step. Whether translation elongation process can feedback to regulate initiation efficiency is unclear. Codon usage bias, a universal feature of all genomes, plays an important role in determining gene expression levels. Here, we discovered that there is a conserved but codon usage-dependent genome-wide negative correlation between protein abundance and CDS length. The codon usage effects on protein expression and ribosome flux on mRNAs are influenced by CDS length; optimal codon usage preferentially promotes production of large proteins. Translation of mRNAs with long CDS and non-optimal codon usage preferentially induces phosphorylation of initiation factor eIF2α, which inhibits translation initiation efficiency. Deletion of the eIF2α kinase CPC-3 (GCN2 homolog) in Neurospora preferentially up-regulates large proteins encoded by non-optimal codons. Surprisingly, CPC-3 also inhibits translation elongation rate in a codon usage and CDS length-dependent manner, resulting in slow elongation rates for long CDS mRNAs. Together, these results revealed a codon usage and CDS length-dependent feedback mechanism from translation elongation to regulate both translation initiation and elongation kinetics.

Clock proteins regulate spatiotemporal organization of clock genes to control circadian rhythms

Circadian clocks regulate ∼24-h oscillations in gene expression, behavior, and physiology. While the genetic and molecular mechanisms of circadian rhythms are well characterized, what remains poorly understood are the intracellular dynamics of circadian clock components and how they affect circadian rhythms. Here, we elucidate how spatiotemporal organization and dynamics of core clock proteins and genes affect circadian rhythms in Drosophila clock neurons. Using high-resolution imaging and DNA-fluorescence in situ hybridization techniques, we demonstrate that Drosophila clock proteins (PERIOD and CLOCK) are organized into a few discrete foci at the nuclear envelope during the circadian repression phase and play an important role in the subnuclear localization of core clock genes to control circadian rhythms. Specifically, we show that core clock genes, period and timeless, are positioned close to the nuclear periphery by the PERIOD protein specifically during the repression phase, suggesting that subnuclear localization of core clock genes might play a key role in their rhythmic gene expression. Finally, we show that loss of Lamin B receptor, a nuclear envelope protein, leads to disruption of PER foci and per gene peripheral localization and results in circadian rhythm defects. These results demonstrate that clock proteins play a hitherto unexpected role in the subnuclear reorganization of core clock genes to control circadian rhythms, revealing how clocks function at the subcellular level. Our results further suggest that clock protein foci might regulate dynamic clustering and spatial reorganization of clock-regulated genes over the repression phase to control circadian rhythms in behavior and physiology.

The codon usage code for co-translational folding of viral capsids

Codon bias is common to all organisms and is the result of mutation, drift, and selection. Selection for the efficiency and accuracy of translation is well recognized as a factor shaping the codon usage. In contrast, fewer studies report the control of the rate of translation as an additional selective pressure influencing the codon usage of an organism. Experimental molecular evolution using RNA virus populations is a powerful tool for the identification of mechanisms underlying the codon bias. Indeed, the role of deoptimized codons on the cotranslational folding has been proven in the capsids of two fecal-orally transmitted picornaviruses, poliovirus, and the hepatitis A virus, emphasizing the role of the frequency of codons in determining the phenotype. However, most studies on virus codon usage rely only on computational analyses, and experimental studies should be encouraged to clearly define the role of selection on codon evolution.

Codon usage bias regulates gene expression and protein conformation in yeast expression system P. pastoris

Background

Protein synthesis is one of the extremely important anabolic pathways in the yeast expression system Pichia pastoris. Codon optimization is a commonly adopted strategy for improved protein expression, although unexpected failures did appear sometimes waiting for further exploration. Recently codon bias has been studied to regulate protein folding and activity in many other organisms.

Results

Here the codon bias profile of P. pastoris genome was examined first and a direct correlation between codon translation efficiency and usage frequency was identified. By manipulating the codon choices of both endogenous and heterologous signal peptides, secretion abilities of N-terminal signal peptides were shown to be tolerant towards codon changes. Then two gene candidates with different levels of structural disorder were studied, and full-length codon optimization was found to affect their expression profiles differentially. Finally, more evidences were provided to support possible protein conformation change brought by codon optimization in structurally disordered proteins.

Conclusion

Our results suggest that codon bias regulates gene expression by modulating several factors including transcription and translation efficiency, protein folding and activity. Because of sequences difference, the extent of affection may be gene specific. For some genes, special codon optimization strategy should be adopted to ensure appropriate expression and conformation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01580-9.

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

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

The Making and Breaking of RNAs: Dynamics of Rhythmic RNA Expression in Mammals

The identification and characterization of rhythmically expressed mRNAs have been an active area of research over the past 20 years, as these mRNAs are believed to produce the daily rhythms in a wide range of biological processes. Circadian transcriptome studies have used mature mRNA as a primary readout and focused largely on rhythmic RNA synthesis as a regulatory mechanism underlying rhythmic mRNA expression. However, RNA synthesis, RNA degradation, or a combination of both must be rhythmic to drive rhythmic RNA profiles, and it is still unclear to what extent rhythmic synthesis leads to rhythmic RNA profiles. In addition, circadian RNA expression is also often tissue specific. Although a handful of genes cycle in all or most tissues, others are rhythmic only in certain tissues, even though the same core clock mechanism is believed to control the rhythmic RNA profiles in all tissues. This review focuses on the dynamics of rhythmic RNA synthesis and degradation and discusses how these steps collectively determine the rhythmicity, phase, and amplitude of RNA accumulation. In particular, we highlight a possible role of RNA degradation in driving tissue-specific RNA rhythms. By unifying findings from experimental and theoretical studies, we will provide a comprehensive overview of how rhythmic gene expression can be achieved and how each regulatory step contributes to tissue-specific circadian transcriptome output in mammals.

Nonoptimal Codon Usage Is Critical for Protein Structure and Function of the Master General Amino Acid Control Regulator CPC-1

Under amino acid starvation conditions, eukaryotic organisms activate a general amino acid control response. In Neurospora crassa, Cross Pathway Control Protein 1 (CPC-1), the ortholog of the Saccharomyces cerevisiae bZIP transcription factor GCN4, functions as the master regulator of the general amino acid control response. Codon usage biases are a universal feature of eukaryotic genomes and are critical for regulation of gene expression. Although codon usage has also been implicated in the regulation of protein structure and function, genetic evidence supporting this conclusion is very limited. Here, we show that Neurospora cpc-1 has a nonoptimal NNU-rich codon usage profile that contrasts with the strong NNC codon preference in the genome. Although substitution of the cpc-1 NNU codons with synonymous NNC codons elevated CPC-1 expression in Neurospora, it altered the CPC-1 degradation rate and abolished its amino acid starvation-induced protein stabilization. The codon-manipulated CPC-1 protein also exhibited different sensitivity to limited protease digestion. Furthermore, CPC-1 functions in rescuing the cell growth of the cpc-1 deletion mutant and activation of the expression of its target genes were impaired by the synonymous codon changes. Together, these results reveal the critical role of codon usage in regulation of CPC-1 expression and function and establish a genetic example of the importance of codon usage in protein folding.IMPORTANCE The general amino acid control response is critical for adaptation of organisms to amino acid starvation conditions. The preference to use certain synonymous codons is a universal feature of all genomes. Synonymous codon changes were previously thought to be silent mutations. In this study, we showed that the Neurospora cpc-1 gene has an unusual codon usage profile compared to other genes in the genome. We found that codon optimization of the cpc-1 gene without changing its amino acid sequence resulted in elevated CPC-1 expression, an altered protein degradation rate, and impaired protein functions due to changes in protein structure. Together, these results reveal the critical role of synonymous codon usage in regulation of CPC-1 expression and function and establish a genetic example of the importance of codon usage in protein structure.

The Ongoing Quest to Crack the Genetic Code for Protein Production

Understanding the genetic design principles that determine protein production remains a major challenge. Although the key principles of gene expression were discovered 50 years ago, additional factors are still being uncovered. Both protein-coding and non-coding sequences harbor elements that collectively influence the efficiency of protein production by modulating transcription, mRNA decay, and translation. The influences of many contributing elements are intertwined, which complicates a full understanding of the individual factors. In natural genes, a functional balance between these factors has been obtained in the course of evolution, whereas for genetic-engineering projects, our incomplete understanding still limits optimal design of synthetic genes. However, notable advances have recently been made, supported by high-throughput analysis of synthetic gene libraries as well as by state-of-the-art biomolecular techniques. We discuss here how these advances further strengthen understanding of the gene expression process and how they can be harnessed to optimize protein production.

A Crosstalk on Codon Usage in Genes Associated with Leukemia

Leukemia is the outcome of aggregation of damaged white blood cells. Several genes were reported to be associated with the pathogenesis of leukemia. These genes were computationally analyzed to decipher their codon usage bias (CUB) and to identify the prime factors influencing the codon usage profile as no work was reported yet. The mean values of synonymous codon usage order (SCUO) parameter indicated low CUB of the genes. Significant positive association of SCUO with overall GC and positional GCs might signal the presence of mutational pressure. However, neutrality plot suggested the dominant role of natural selection across the genes. Along with natural selection, the role of mutation pressure was also prominent and that might be responsible for lower CUB (SCUO = 0.19) of genes. Low translational speed might permit accuracy in the process. A strong inverse relationship of translational rate was observed with CUB of genes and folding energy.

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

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

A comparative analysis of the mobility of 45 proteins in the synaptic bouton

Many proteins involved in synaptic transmission are well known, and their features, as their abundance or spatial distribution, have been analyzed in systematic studies. This has not been the case, however, for their mobility. To solve this, we analyzed the motion of 45 GFP‐tagged synaptic proteins expressed in cultured hippocampal neurons, using fluorescence recovery after photobleaching, particle tracking, and modeling. We compared synaptic vesicle proteins, endo‐ and exocytosis cofactors, cytoskeleton components, and trafficking proteins. We found that movement was influenced by the protein association with synaptic vesicles, especially for membrane proteins. Surprisingly, protein mobility also correlated significantly with parameters as the protein lifetimes, or the nucleotide composition of their mRNAs. We then analyzed protein movement thoroughly, taking into account the spatial characteristics of the system. This resulted in a first visualization of overall protein motion in the synapse, which should enable future modeling studies of synaptic physiology.

Codon usage trend in genes associated with obesity

Obesity is not only a social menace but also an economic burden as it reduces productivity and increases health care cost. We used bioinformatic tools to analyze the CUB of obesity associated genes and compared with housekeeping genes (control) to explore the similarities and differences between two data sets as no work was reported yet. The mean effective number of codons (ENC) in genes associated with obesity and housekeeping gene was 50.45 and 52.03 respectively, indicating low CUB. The relative synonymous codon usage (RSCU) suggested that codons namely CTG and GTG were over-represented in both obesity and housekeeping genes while under-represented codons were TCG, TTA, CTA, CCG, CAA, CGT, ATA, ACG, GTA and GCG in obesity genes and TCG, TTA, CCG, ATA, ACG, GTA, and GCG in housekeeping genes. t test analysis suggested that 11 codons namely TTA (Leu), TTG (Leu), CCG (Pro), CAC (His), CAA (Gln), CAG (Gln), CGT (Arg), AGA (Arg), ATA (Ile), ATT (Ile) and GCG (Ala) were significantly differed (p < 0.05 or p < 0.01) between obesity and housekeeping genes. Highly significant correlation was observed between GC12 and GC3 in obesity and housekeeping genes i.e. r = 0.580** and r = 0.498** (p < 0.01) respectively indicating the effect of directional mutation pressure present in all codon positions.

Emerging Roles of Translational Control in Circadian Timekeeping

A large part of mammalian physiology and behaviour shows regular daily variations. This temporal organisation is driven by the activity of an endogenous circadian clock, whose molecular basis consists of diurnal waves in gene expression. Circadian transcription is the major driver of these rhythms, yet post-transcriptional mechanisms, some of which occur in response to systemic cues and in a tissue-specific fashion, have central roles in ultimately establishing the oscillatory gene expression programme as well. Regulatory control that occurs at the level of translation is emerging as an important player in the generation and modulation of protein accumulation rhythms. As a mechanism, translation lies at a privileged position to integrate genetically encoded rhythmic signals with other, external and internal stimuli, including nutrient-derived cues. In this review, we summarise our current knowledge of how diurnal control of translation affects both bulk protein levels and gene-specific protein biosynthesis. We discuss mechanisms of regulation, in particular with regard to the complex interplay between circadian cycles and feeding/fasting cycles, as well as emerging roles for upstream open reading frames in clock control.

Genome-wide codon usage pattern analysis reveals the correlation between codon usage bias and gene expression in Cuscuta australis

The protein-coding genes and pseudogenes of Cuscuta australis had the diverse contribution to the formation and evolution of parasitism. The codon usage pattern analysis of these two type genes could be used to understand the gene transcription and translation. In this study, we systematically analyzed the codon usage patterns of protein-coding sequences and pseudogenes sequences in C. australis. The results showed that the high frequency codons of protein coding sequences and pseudogenes had the same A/U bias in the third position. However, these two sequences had converse bias at the third base in optimal codons: the protein coding sequences preferred G/C-ending codons while pseudogene sequences preferred A/U-ending codons. Neutrality plot and effective number of codons plot revealed that natural selection played a more important role than mutation pressure in two sequences codon usage bias. Furthermore, the gene expression level had a significant positive correlation with codon usage bias in C. australis. Highly-expressed protein coding genes exhibited a higher codon bias than lowly-expressed genes. Meanwhile, the high-expression genes tended to use G/C-ending synonymous codons. This result further verified the optimal codons usage bias and its correlation with the gene expression in C. australis.

Cotranslational folding allows misfolding-prone proteins to circumvent deep kinetic traps

Many large proteins suffer from slow or inefficient folding in vitro. It has long been known that this problem can be alleviated in vivo if proteins start folding cotranslationally. However, the molecular mechanisms underlying this improvement have not been well established. To address this question, we use an all-atom simulation-based algorithm to compute the folding properties of various large protein domains as a function of nascent chain length. We find that for certain proteins, there exists a narrow window of lengths that confers both thermodynamic stability and fast folding kinetics. Beyond these lengths, folding is drastically slowed by nonnative interactions involving C-terminal residues. Thus, cotranslational folding is predicted to be beneficial because it allows proteins to take advantage of this optimal window of lengths and thus avoid kinetic traps. Interestingly, many of these proteins' sequences contain conserved rare codons that may slow down synthesis at this optimal window, suggesting that synthesis rates may be evolutionarily tuned to optimize folding. Using kinetic modeling, we show that under certain conditions, such a slowdown indeed improves cotranslational folding efficiency by giving these nascent chains more time to fold. In contrast, other proteins are predicted not to benefit from cotranslational folding due to a lack of significant nonnative interactions, and indeed these proteins' sequences lack conserved C-terminal rare codons. Together, these results shed light on the factors that promote proper protein folding in the cell and how biomolecular self-assembly may be optimized evolutionarily.

eRF1 mediates codon usage effects on mRNA translation efficiency through premature termination at rare codons

Codon usage bias is a universal feature of eukaryotic and prokaryotic genomes and plays an important role in regulating gene expression levels. A major role of codon usage is thought to regulate protein expression levels by affecting mRNA translation efficiency, but the underlying mechanism is unclear. By analyzing ribosome profiling results, here we showed that codon usage regulates translation elongation rate and that rare codons are decoded more slowly than common codons in all codon families in Neurospora. Rare codons resulted in ribosome stalling in manners both dependent and independent of protein sequence context and caused premature translation termination. This mechanism was shown to be conserved in Drosophila cells. In both Neurospora and Drosophila cells, codon usage plays an important role in regulating mRNA translation efficiency. We found that the rare codon-dependent premature termination is mediated by the translation termination factor eRF1, which recognizes ribosomes stalled on rare sense codons. Silencing of eRF1 expression resulted in codon usage-dependent changes in protein expression. Together, these results establish a mechanism for how codon usage regulates mRNA translation efficiency.

A pair of non-optimal codons are necessary for the correct biosynthesis of the Aspergillus nidulans urea transporter, UreA

In both prokaryotic and eukaryotic genomes, synonymous codons are unevenly used. Such differential usage of optimal or non-optimal codons has been suggested to play a role in the control of translation initiation and elongation, as well as at the level of transcription and mRNA stability. In the case of membrane proteins, codon usage has been proposed to assist in the establishment of a pause necessary for the correct targeting of the nascent chains to the translocon. By using as a model UreA, the Aspergillus nidulans urea transporter, we revealed that a pair of non-optimal codons encoding amino acids situated at the boundary between the N-terminus and the first transmembrane segment are necessary for proper biogenesis of the protein at 37°C. These codons presumably regulate the translation rate in a previously undescribed fashion, possibly contributing to the correct interaction of ureA-translating ribosome-nascent chain complexes with the signal recognition particle and/or other factors, while the polypeptide has not yet emerged from the ribosomal tunnel. Our results suggest that the presence of the pair of non-optimal codons would not be functionally important in all cellular conditions. Whether this mechanism would affect other proteins remains to be determined.

Honey bee predisposition of resistance to ubiquitous mite infestations

Host-parasite co-evolution history is lacking when parasites switch to novel hosts. This was the case for Western honey bees (Apis mellifera) when the ectoparasitic mite, Varroa destructor, switched hosts from Eastern honey bees (Apis cerana). This mite has since become the most severe biological threat to A. mellifera worldwide. However, some A. mellifera populations are known to survive infestations, largely by suppressing mite population growth. One known mechanism is suppressed mite reproduction (SMR), but the underlying genetics are poorly understood. Here, we take advantage of haploid drones, originating from one queen from the Netherlands that developed Varroa-resistance, whole exome sequencing and elastic-net regression to identify genetic variants associated with SMR in resistant honeybees. An eight variants model predicted 88% of the phenotypes correctly and identified six risk and two protective variants. Reproducing and non-reproducing mites could not be distinguished using DNA microsatellites, which is in agreement with the hypothesis that it is not the parasite but the host that adapted itself. Our results suggest that the brood pheromone-dependent mite oogenesis is disrupted in resistant hosts. The identified genetic markers have a considerable potential to contribute to a sustainable global apiculture.

CodSeqGen: A tool for generating synonymous coding sequences with desired GC-contents

Identification of regulatory elements is essential for understanding the mechanism behind regulating gene expression. These regulatory elements—located in or near gene—bind to proteins called transcription factors to initiate the transcription process. Their occurrences are influenced by the GC-content or nucleotide composition. For generating synthetic coding sequences with pre-specified amino acid sequence and desired GC-content, there exist two stochastic methods, multinomial and maximum entropy. Both methods rely on the probability of choosing the codon synonymous for usage in regard to a specific amino acid. In spite the latter exhibited unbiased manner, the produced sequences are not exactly obeying the GC-content constraint. In this paper, we present an algorithmic solution to produce coding sequences that follow exactly a primary amino acid sequence and a desired GC-content. The proposed tool, namely CodSeqGen, depends on random selection for smaller subsets to be traversed using the backtracking approach.

•

Generating random coding sequence is computationally expensive.

•

For proteins of length n, we need to test 6ⁿ coding sequences to find those satisfying GC-content constraints.

•

CodSeqGen uses smart backtracking, and produces synonymous coding sequences with the exact desired GC-contents.

•

Experiments show CodSeqGen produces sequences with GC-contents that satisfy the constraints while NullSeq does not.

O-GlcNAcylation of PERIOD regulates its interaction with CLOCK and timing of circadian transcriptional repression

Circadian clocks coordinate time-of-day-specific metabolic and physiological processes to maximize organismal performance and fitness. In addition to light and temperature, which are regarded as strong zeitgebers for circadian clock entrainment, metabolic input has now emerged as an important signal for clock entrainment and modulation. Circadian clock proteins have been identified to be substrates of O-GlcNAcylation, a nutrient sensitive post-translational modification (PTM), and the interplay between clock protein O-GlcNAcylation and other PTMs is now recognized as an important mechanism by which metabolic input regulates circadian physiology. To better understand the role of O-GlcNAcylation in modulating clock protein function within the molecular oscillator, we used mass spectrometry proteomics to identify O-GlcNAcylation sites of PERIOD (PER), a repressor of the circadian transcriptome and a critical biochemical timer of the Drosophila clock. In vivo functional characterization of PER O-GlcNAcylation sites indicates that O-GlcNAcylation at PER(S942) reduces interactions between PER and CLOCK (CLK), the key transcriptional activator of clock-controlled genes. Since we observe a correlation between clock-controlled daytime feeding activity and higher level of PER O-GlcNAcylation, we propose that PER(S942) O-GlcNAcylation during the day functions to prevent premature initiation of circadian repression phase. This is consistent with the period-shortening behavioral phenotype of per(S942A) flies. Taken together, our results support that clock-controlled feeding activity provides metabolic signals to reinforce light entrainment to regulate circadian physiology at the post-translational level. The interplay between O-GlcNAcylation and other PTMs to regulate circadian physiology is expected to be complex and extensive, and reach far beyond the molecular oscillator.

The codon sequences predict protein lifetimes and other parameters of the protein life cycle in the mouse brain

The homeostasis of the proteome depends on the tight regulation of the mRNA and protein abundances, of the translation rates, and of the protein lifetimes. Results from several studies on prokaryotes or eukaryotic cell cultures have suggested that protein homeostasis is connected to, and perhaps regulated by, the protein and the codon sequences. However, this has been little investigated for mammals in vivo. Moreover, the link between the coding sequences and one critical parameter, the protein lifetime, has remained largely unexplored, both in vivo and in vitro. We tested this in the mouse brain, and found that the percentages of amino acids and codons in the sequences could predict all of the homeostasis parameters with a precision approaching experimental measurements. A key predictive element was the wobble nucleotide. G-/C-ending codons correlated with higher protein lifetimes, protein abundances, mRNA abundances and translation rates than A-/U-ending codons. Modifying the proportions of G-/C-ending codons could tune these parameters in cell cultures, in a proof-of-principle experiment. We suggest that the coding sequences are strongly linked to protein homeostasis in vivo, albeit it still remains to be determined whether this relation is causal in nature.

CK1α Collaborates with DOUBLETIME to Regulate PERIOD Function in the Drosophila Circadian Clock

The animal circadian timing system interprets environmental time cues and internal metabolic status to orchestrate circadian rhythms of physiology, allowing animals to perform necessary tasks in a time-of-day-dependent manner. Normal progression of circadian rhythms is dependent on the daily cycling of core transcriptional factors that make up cell-autonomous molecular oscillators. In Drosophila, PERIOD (PER), TIMELESS (TIM), CLOCK (CLK), and CYCLE (CYC) are core clock proteins that function in a transcriptional-translational feedback mechanism to regulate the circadian transcriptome. Posttranslational modifications of core clock proteins provide precise temporal control over when they are active as regulators of clock-controlled genes. In particular, phosphorylation is a key regulatory mechanism that dictates the subcellular localization, stability, and transcriptional activity of clock proteins. Previously, casein kinase 1α (CK1α) has been identified as a kinase that phosphorylates mammalian PER1 and modulates its stability, but the mechanisms by which it modulates PER protein stability is still unclear. Using Drosophila as a model, we show that CK1α has an overall function of speeding up PER metabolism and is required to maintain the 24 h period of circadian rhythms. Our results indicate that CK1α collaborates with the key clock kinase DOUBLETIME (DBT) in both the cytoplasm and the nucleus to regulate the timing of PER-dependent repression of the circadian transcriptome. Specifically, we observe that CK1α promotes PER nuclear localization by antagonizing the activity of DBT to inhibit PER nuclear translocation. Furthermore, CK1α enhances DBT-dependent PER phosphorylation and degradation once PER moves into the nucleus.SIGNIFICANCE STATEMENT Circadian clocks are endogenous timers that integrate environmental signals to impose temporal control over organismal physiology over the 24 h day/night cycle. To maintain the 24 h period length of circadian clocks and to ensure that circadian rhythms are in synchrony with the external environment, key proteins that make up the molecular oscillator are extensively regulated by phosphorylation to ensure that they perform proper time-of-day-specific functions. Casein kinase 1α (CK1α) has previously been identified as a kinase that phosphorylates mammalian PERIOD (PER) proteins to promote their degradation, but the mechanism by which it modulates PER stability is unclear. In this study, we characterize the mechanisms by which CK1α interacts with DOUBLETIME (DBT) to achieve the overall function of speeding up PER metabolism and to ensure proper time-keeping.

Codon usage regulates human KRAS expression at both transcriptional and translational levels

KRAS and HRAS are highly homologous oncogenic Ras GTPase family members that are mutated in a wide spectrum of human cancers. Despite having high amino acid identity, KRAS and HRAS have very different codon usage biases: the HRAS gene contains many common codons, and KRAS is enriched for rare codons. Rare codons in KRAS suppress its protein expression, which has been shown to affect both normal and cancer biology in mammals. Here, using HRAS or KRAS expression in different human cell lines and in vitro transcription and translation assays, we show that KRAS rare codons inhibit both translation efficiency and transcription and that the contribution of these two processes varies among different cell lines. We observed that codon usage regulates mRNA translation efficiency such that WT KRAS mRNA is poorly translated. On the other hand, common codons increased transcriptional rates by promoting activating histone modifications and recruitment of transcriptional coactivators. Finally, we found that codon usage also influences KRAS protein conformation, likely because of its effect on co-translational protein folding. Together, our results reveal that codon usage has multidimensional effects on protein expression, ranging from effects on transcription to protein folding in human cells.

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.

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.

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.

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.

Genetic Code Optimization for Cotranslational Protein Folding: Codon Directional Asymmetry Correlates with Antiparallel Betasheets, tRNA Synthetase Classes

A new codon property, codon directional asymmetry in nucleotide content (CDA), reveals a biologically meaningful genetic code dimension: palindromic codons (first and last nucleotides identical, codon structure XZX) are symmetric (CDA = 0), codons with structures ZXX/XXZ are 5'/3' asymmetric (CDA = - 1/1; CDA = - 0.5/0.5 if Z and X are both purines or both pyrimidines, assigning negative/positive (-/+) signs is an arbitrary convention). Negative/positive CDAs associate with (a) Fujimoto's tetrahedral codon stereo-table; (b) tRNA synthetase class I/II (aminoacylate the 2'/3' hydroxyl group of the tRNA's last ribose, respectively); and (c) high/low antiparallel (not parallel) betasheet conformation parameters. Preliminary results suggest CDA-whole organism associations (body temperature, developmental stability, lifespan). Presumably, CDA impacts spatial kinetics of codon-anticodon interactions, affecting cotranslational protein folding. Some synonymous codons have opposite CDA sign (alanine, leucine, serine, and valine), putatively explaining how synonymous mutations sometimes affect protein function. Correlations between CDA and tRNA synthetase classes are weaker than between CDA and antiparallel betasheet conformation parameters. This effect is stronger for mitochondrial genetic codes, and potentially drives mitochondrial codon-amino acid reassignments. CDA reveals information ruling nucleotide-protein relations embedded in reversed (not reverse-complement) sequences (5'-ZXX-3'/5'-XXZ-3').

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.