Roles for Synonymous Codon Usage in Protein Biogenesis

Effect of genome composition and codon bias on infectious bronchitis virus evolution and adaptation to target tissues

Background

Infectious bronchitis virus (IBV) is one of the most relevant viruses affecting the poultry industry, and several studies have investigated the factors involved in its biological cycle and evolution. However, very few of those studies focused on the effect of genome composition and the codon bias of different IBV proteins, despite the remarkable increase in available complete genomes. In the present study, all IBV complete genomes were downloaded (n = 383), and several statistics representative of genome composition and codon bias were calculated for each protein-coding sequence, including but not limited to, the nucleotide odds ratio, relative synonymous codon usage and effective number of codons. Additionally, viral codon usage was compared to host codon usage based on a collection of highly expressed genes in IBV target and nontarget tissues.

Results

The results obtained demonstrated a significant difference among structural, non-structural and accessory proteins, especially regarding dinucleotide composition, which appears under strong selective forces. In particular, some dinucleotide pairs, such as CpG, a probable target of the host innate immune response, are underrepresented in genes coding for pp1a, pp1ab, S and N. Although genome composition and dinucleotide bias appear to affect codon usage, additional selective forces may act directly on codon bias. Variability in relative synonymous codon usage and effective number of codons was found for different proteins, with structural proteins and polyproteins being more adapted to the codon bias of host target tissues. In contrast, accessory proteins had a more biased codon usage (i.e., lower number of preferred codons), which might contribute to the regulation of their expression level and timing throughout the cell cycle.

Conclusions

The present study confirms the existence of selective forces acting directly on the genome and not only indirectly through phenotype selection. This evidence might help understanding IBV biology and in developing attenuated strains without affecting the protein phenotype and therefore immunogenicity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07559-5.

Smart-ORF: a single-molecule method for accessing ribosome dynamics in both upstream and main open reading frames

Upstream open reading frame (uORF) translation disrupts scanning 43S flux on mRNA and modulates main open reading frame (mORF) translation efficiency. Current tools, however, have limited access to ribosome dynamics in both upstream and main ORFs of an mRNA. Here, we develop a new two-color in vitro fluorescence assay, Smart-ORF, that monitors individual uORF and mORF translation events in real-time with single-molecule resolution. We demonstrate the utility of Smart-ORF by applying it to uORF-encoded arginine attenuator peptide (AAP)-mediated translational regulation. The method enabled quantification of uORF and mORF initiation efficiencies, 80S dwell time, polysome formation, and the correlation between uORF and mORF translation dynamics. Smart-ORF revealed that AAP-mediated 80S stalling in the uORF stimulates the uORF initiation efficiency and promotes clustering of slower uORF-translating ribosomes. This technology provides a new tool that can reveal previously uncharacterized dynamics of uORF-containing mRNA translation.

The effects of codon bias and optimality on mRNA and protein regulation

The central dogma of molecular biology entails that genetic information is transferred from nucleic acid to proteins. Notwithstanding retro-transcribing genetic elements, DNA is transcribed to RNA which in turn is translated into proteins. Recent advancements have shown that each stage is regulated to control protein abundances for a variety of essential physiological processes. In this regard, mRNA regulation is essential in fine-tuning or calibrating protein abundances. In this review, we would like to discuss one of several mRNA-intrinsic features of mRNA regulation that has been gaining traction of recent-codon bias and optimality. Specifically, we address the effects of codon bias with regard to codon optimality in several biological processes centred on translation, such as mRNA stability and protein folding among others. Finally, we examine how different organisms or cell types, through this system, are able to coordinate physiological pathways to respond to a variety of stress or growth conditions.

Inferring translational heterogeneity from Saccharomyces cerevisiae ribosome profiling

Translation of mRNAs into proteins by the ribosome is the most important step of protein biosynthesis. Accordingly, translation is tightly controlled and heavily regulated to maintain cellular homeostasis. Ribosome profiling (Ribo-seq) has revolutionized the study of translation by revealing many of its underlying mechanisms. However, equally many aspects of translation remain mysterious, in part also due to persisting challenges in the interpretation of data obtained from Ribo-seq experiments. Here, we show that some of the variability observed in Ribo-seq data has biological origins and reflects programmed heterogeneity of translation. Through a comparative analysis of Ribo-seq data from Saccharomyces cerevisiae, we systematically identify short 3-codon sequences that are differentially translated (DT) across mRNAs, that is, identical sequences that are translated sometimes fast and sometimes slowly beyond what can be attributed to variability between experiments. Remarkably, the thus identified DT sequences link to mechanisms known to regulate translation elongation and are enriched in genes important for protein and organelle biosynthesis. Our results thus highlight examples of translational heterogeneity that are encoded in the genomic sequences and tuned to optimizing cellular homeostasis. More generally, our work highlights the power of Ribo-seq to understand the complexities of translation regulation.

Effects of codon usage on gene expression are promoter context dependent

Codon usage bias is a universal feature of all genomes. Although codon usage has been shown to regulate mRNA and protein levels by influencing mRNA decay and transcription in eukaryotes, little or no genome-wide correlations between codon usage and mRNA levels are detected in mammalian cells, raising doubt on the significance of codon usage effect on gene expression. Here we show that gene-specific regulation reduces the genome-wide codon usage and mRNA correlations: Constitutively expressed genes exhibit much higher genome-wide correlations than differentially expressed genes from fungi to human cells. Using Drosophila S2 cells as a model system, we showed that the effect of codon usage on mRNA expression level is promoter-dependent. Regions downstream of the core promoters of differentially expressed genes can repress the codon usage effects on mRNA expression. An element in the Hsp70 promoter was identified to be necessary and sufficient for this inhibitory effect. The promoter-dependent codon usage effects on mRNA levels are regulated at the transcriptional level through modulation of histone modifications, nucleosome densities and premature termination. Together, our results demonstrate that promoters play a major role in determining whether codon usage influences gene expression and further establish the transcription-dependent codon usage effects on gene expression.

Modification of Transfer RNA Levels Affects Cyclin Aggregation and the Correct Duplication of Yeast Cells

Codon usage bias (the preferential use of certain synonymous codons (optimal) over others is found at the organism level (intergenomic) within specific genomes (intragenomic) and even in certain genes. Whether it is the result of genetic drift due to GC/AT content and/or natural selection is a topic of intense debate. Preferential codons are mostly found in genes encoding highly-expressed proteins, while lowly-expressed proteins usually contain a high proportion of rare (lowly-represented) codons. While optimal codons are decoded by highly expressed tRNAs, rare codons are usually decoded by lowly-represented tRNAs. Whether rare codons play a role in controlling the expression of lowly- or temporarily-expressed proteins is an open question. In this work we approached this question using two strategies, either by replacing rare glycine codons with optimal counterparts in the gene that encodes the cell cycle protein Cdc13, or by overexpression the tRNA ^Gly that decodes rare codons from the fission yeast, Schizosaccharomyces pombe. While the replacement of synonymous codons severely affected cell growth, increasing tRNA levels affected the aggregation status of Cdc13 and cell division. These lead us to think that rare codons in lowly-expressed cyclin proteins are crucial for cell division, and that the overexpression of tRNA that decodes rare codons affects the expression of proteins containing these rare codons. These codons may be the result of the natural selection of codons in genes that encode lowly-expressed proteins.

Comprehensive analysis of coding sequence architecture features and gene expression in Arachis duranensis

Coding sequence (CDS) architecture affects gene expression levels in organisms. Codon optimization can increase the gene expression level. Therefore, understanding codon usage patterns has important implications for research on genetic engineering and exogenous gene expression. To date, the codon usage patterns of many model plants have been analyzed. However, the relationship between CDS architecture and gene expression in Arachis duranensis remains poorly understood. According to the results of genome sequencing, A. duranensis has many resistant genes that can be used to improve the cultivated peanut. In this study, bioinformatic approaches were used to estimate A. duranensis CDS architectures, including frequency of the optimal codon (Fop), polypeptide length and GC contents at the first (GC1), second (GC2) and third (GC3) codon positions. In addition, Arachis RNA-seq datasets were downloaded from PeanutBase. The relationships between gene expression and CDS architecture were assessed both under normal growth as well as nematode and drought stress conditions. A total of 26 codons with high frequency were identified, which preferentially ended with A or T in A. duranensis CDSs under the above-mentioned three conditions. A similar CDS architecture was found in differentially expressed genes (DEGs) under nematode and drought stresses. The GC1 content differed between DEGs and non-differentially expressed genes (NDEGs) under both drought and nematode stresses. The expression levels of DEGs were affected by different CDS architectures compared with NDEGs under drought stress. In addition, no correlation was found between differential gene expression and CDS architecture neither under nematode nor under drought stress. These results aid the understanding of gene expression in A. duranensis.

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.

Translational Control by Ribosome Pausing in Bacteria: How a Non-uniform Pace of Translation Affects Protein Production and Folding

Protein homeostasis of bacterial cells is maintained by coordinated processes of protein production, folding, and degradation. Translational efficiency of a given mRNA depends on how often the ribosomes initiate synthesis of a new polypeptide and how quickly they read the coding sequence to produce a full-length protein. The pace of ribosomes along the mRNA is not uniform: periods of rapid synthesis are separated by pauses. Here, we summarize recent evidence on how ribosome pausing affects translational efficiency and protein folding. We discuss the factors that slow down translation elongation and affect the quality of the newly synthesized protein. Ribosome pausing emerges as important factor contributing to the regulatory programs that ensure the quality of the proteome and integrate the cellular and environmental cues into regulatory circuits of the cell.

DNA Repair in Huntington's Disease and Spinocerebellar Ataxias: Somatic Instability and Alternative Hypotheses

The use of genome wide association studies (GWAS) in Huntington's disease (HD) research, driven by unbiased human data analysis, has transformed the focus of new targets that could affect age at onset. While there is a significant depth of information on DNA damage repair, with many drugs and drug targets, most of this development has taken place in the context of cancer therapy. DNA damage repair in neurons does not rely on DNA replication correction mechanisms. However, there is a strong connection between DNA repair and neuronal metabolism, mediated by nucleotide salvaging and the poly ADP-ribose (PAR) response, and this connection has been implicated in other age-onset neurodegenerative diseases. Validation of leads including the mismatch repair protein MSH3, and interstrand cross-link repair protein FAN1, suggest the mechanism is driven by somatic CAG instability, which is supported by the protective effect of CAA substitutions in the CAG tract. We currently do not understand: how somatic instability is triggered; the state of DNA damage within expanding alleles in the brain; whether this damage induces mismatch repair and interstrand cross-link pathways; whether instability mediates toxicity, and how this relates to human ageing. We discuss DNA damage pathways uncovered by HD GWAS, known roles of other polyglutamine disease proteins in DNA damage repair, and a panel of hypotheses for pathogenic mechanisms.

De novo transcriptome assembly, functional annotation and characterization of the Atlantic bluefin tuna (Thunnus thynnus) larval stage

In the present work, we assembled and characterized a de novo larval transcriptome of the Atlantic bluefin tuna Thunnus thynnus by taking advantage of publicly available databases with the goal of better understanding its larval development. The assembled transcriptome comprised 37,117 protein-coding transcripts, of which 13,633 full-length (>80% coverage), with an Ex90N50 of 3061 bp and 76% of complete and single-copy core vertebrate genes orthologues. Of these transcripts, 34,980 had a hit against the EggNOG database and 14,983 with the KEGG database. Codon usage bias was identified in processes such as translation and muscle development. By comparing our data with a set of representative fish species, 87.1% of tuna transcripts were included in orthogroups with other species and 5.1% in assembly-specific orthogroups, which were enriched in terms related to muscle and bone development, visual system and ion transport. Following this comparative approach, protein families related to myosin, extracellular matrix and immune system resulted significantly expanded in the Atlantic bluefin tuna. Altogether, these results provide a glimpse of how the Atlantic bluefin tuna might have achieved early physical advantages over competing species in the pelagic environment. The information generated lays the foundation for future research on the more detailed exploration of physiological responses at the molecular level in different larval stages and paves the way to evolutionary studies on the Atlantic bluefin tuna.

Network analysis of synonymous codon usage

MOTIVATION

Most amino acids are encoded by multiple synonymous codons, some of which are used more rarely than others. Analyses of positions of such rare codons in protein sequences revealed that rare codons can impact co-translational protein folding and that positions of some rare codons are evolutionarily conserved. Analyses of their positions in protein 3-dimensional structures, which are richer in biochemical information than sequences alone, might further explain the role of rare codons in protein folding.

RESULTS

We model protein structures as networks and use network centrality to measure the structural position of an amino acid. We first validate that amino acids buried within the structural core are network-central, and those on the surface are not. Then, we study potential differences between network centralities and thus structural positions of amino acids encoded by conserved rare, non-conserved rare and commonly used codons. We find that in 84% of proteins, the three codon categories occupy significantly different structural positions. We examine protein groups showing different codon centrality trends, i.e. different relationships between structural positions of the three codon categories. We see several cases of all proteins from our data with some structural or functional property being in the same group. Also, we see a case of all proteins in some group having the same property. Our work shows that codon usage is linked to the final protein structure and thus possibly to co-translational protein folding.

AVAILABILITY AND IMPLEMENTATION

https://nd.edu/∼cone/CodonUsage/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

What determines eukaryotic translation elongation: recent molecular and quantitative analyses of protein synthesis

Protein synthesis from mRNA is an energy-intensive and tightly controlled cellular process. Translation elongation is a well-coordinated, multifactorial step in translation that undergoes dynamic regulation owing to cellular state and environmental determinants. Recent studies involving genome-wide approaches have uncovered some crucial aspects of translation elongation including the mRNA itself and the nascent polypeptide chain. Additionally, these studies have fuelled quantitative and mathematical modelling of translation elongation. In this review, we provide a comprehensive overview of the key determinants of translation elongation. We discuss consequences of ribosome stalling or collision, and how the cells regulate translation in case of such events. Next, we review theoretical approaches and widely used mathematical models that have become an essential ingredient to interpret complex molecular datasets and study translation dynamics quantitatively. Finally, we review recent advances in live-cell reporter and related analysis techniques, to monitor the translation dynamics of single cells and single-mRNA molecules in real time.

Programmed Trade-offs in Protein Folding Networks

Molecular chaperones as specialized protein quality control enzymes form the core of cellular protein homeostasis. How chaperones selectively interact with their substrate proteins thus allocate their overall limited capacity remains poorly understood. Here, I present an integrated analysis of sequence and structural determinants that define interactions of protein domains as the basic protein folding unit with the Saccharomyces cerevisiae Hsp70 Ssb. Structural homologs of single-domain proteins that differentially interact with Ssb for de novo folding were found to systematically differ in complexity of their folding landscapes, selective use of nonoptimal codons, and presence of short discriminative sequences, thus highlighting pervasive trade-offs in chaperone-assisted protein folding landscapes. However, short discriminative sequences were found to contribute by far the strongest signal toward explaining Ssb interactions. This observation suggested that some chaperone interactions may be directly programmed in the amino acid sequences rather than responding to folding challenges, possibly for regulatory advantages.

GTPBP1 resolves paused ribosomes to maintain neuronal homeostasis

Ribosome-associated quality control pathways respond to defects in translational elongation to recycle arrested ribosomes and degrade aberrant polypeptides and mRNAs. Loss of a tRNA gene leads to ribosomal pausing that is resolved by the translational GTPase GTPBP2, and in its absence causes neuron death. Here, we show that loss of the homologous protein GTPBP1 during tRNA deficiency in the mouse brain also leads to codon-specific ribosome pausing and neurodegeneration, suggesting that these non-redundant GTPases function in the same pathway to mitigate ribosome pausing. As observed in Gtpbp2-/- mice (Ishimura et al., 2016), GCN2-mediated activation of the integrated stress response (ISR) was apparent in the Gtpbp1-/- brain. We observed decreased mTORC1 signaling which increased neuronal death, whereas ISR activation was neuroprotective. Our data demonstrate that GTPBP1 functions as an important quality control mechanism during translation elongation and suggest that translational signaling pathways intricately interact to regulate neuronal homeostasis during defective elongation.

Pairs of amino acids at the P- and A-sites of the ribosome predictably and causally modulate translation-elongation rates

Variation in translation-elongation kinetics along a transcript's coding sequence plays an important role in the maintenance of cellular protein homeostasis by regulating co-translational protein folding, localization, and maturation. Translation-elongation speed is influenced by molecular factors within mRNA and protein sequences. For example, the presence of proline in the ribosome's P- or A-site slows down translation, but the effect of other pairs of amino acids, in the context of all 400 possible pairs, has not been characterized. Here, we study Saccharomyces cerevisiae using a combination of bioinformatics, mutational experiments, and evolutionary analyses, and show that many different pairs of amino acids and their associated tRNA molecules predictably and causally encode translation rate information when these pairs are present in the A- and P-sites of the ribosome independent of other factors known to influence translation speed including mRNA structure, wobble base pairing, tripeptide motifs, positively charged upstream nascent chain residues, and cognate tRNA concentration. The fast-translating pairs of amino acids that we identify are enriched four-fold relative to the slow-translating pairs across Saccharomyces cerevisiae's proteome, while the slow-translating pairs are enriched downstream of domain boundaries. Thus, the chemical identity of amino acid pairs contributes to variability in translation rates, elongation kinetics are causally encoded in the primary structure of proteins, and signatures of evolutionary selection indicate their potential role in co-translational processes.

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.

Analysis of the Codon Usage Pattern of HA and NA Genes of H7N9 Influenza A Virus

Novel H7N9 influenza virus transmitted from birds to human and, since March 2013, it has caused five epidemic waves in China. Although the evolution of H7N9 viruses has been investigated, the evolutionary changes associated with codon usage are still unclear. Herein, the codon usage pattern of two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA), was studied to understand the evolutionary changes in relation to host, epidemic wave, and pathogenicity. Both genes displayed a low codon usage bias, with HA higher than NA. The codon usage was driven by mutation pressure and natural selection, although the main contributing factor was natural selection. Additionally, the codon adaptation index (CAI) and deoptimization (RCDI) illustrated the strong adaptability of H7N9 to Gallus gallus. Similarity index (SiD) analysis showed that Homo sapiens posed a stronger selection pressure than Gallus gallus. Thus, we assume that this may be related to the gradual adaptability of the virus to human. In addition, the host strong selection pressure was validated based on CpG dinucleotide content. In conclusion, this study analyzed the usage of codons of two genes of H7N9 and expanded our understanding of H7N9 host specificity. This aids into the development of control measures against H7N9 influenza virus.

Codon usage pattern and evolutionary forces of mitochondrial ND genes among orders of class Amphibia

In this study, we used a bioinformatics approach to analyze the nucleotide composition and pattern of synonymous codon usage in mitochondrial ND genes in three amphibian groups, that is, orders Anura, Caudata, and Gymnophiona to identify the commonality and the differences of codon usage as no research work was reported yet. The high value of the effective number of codons revealed that the codon usage bias (CUB) was low in mitochondrial ND genes among the orders. Nucleotide composition analysis suggested that for each gene, the compositional features differed among Anura, Caudata, and Gymnophiona and the GC content was lower than AT content. Furthermore, a highly significant difference (p < .05) for GC content was found in each gene among the orders. The heat map showed contrasting patterns of codon usage among different ND genes. The regression of GC12 on GC3 suggested a narrow range of GC3 distribution and some points were located in the diagonal, indicating both mutation pressure and natural selection might influence the CUB. Moreover, the slope of the regression line was less than 0.5 in all ND genes among orders, indicating natural selection might have played the dominant role whereas mutation pressure had played a minor role in shaping CUB of ND genes across orders.

Characterization of codon usage pattern in SARS-CoV-2

The outbreak of coronavirus disease 2019 (COVID-19) due to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has posed significant threats to international health. The genetic traits as well as evolutionary processes in this novel coronavirus are not fully characterized, and their roles in viral pathogenesis are yet largely unknown. To get a better picture of the codon architecture of this newly emerging coronavirus, in this study we perform bioinformatic analysis, based on publicly available nucleotide sequences of SARS-CoV-2 along with those of other members of human coronaviruses as well as non-human coronaviruses in different hosts, to take a snapshot of the genome-wide codon usage pattern of SARS-CoV-2 and uncover that all over-represented codons end with A/U and this newly emerging coronavirus has a relatively low codon usage bias, which is shaped by both mutation pressure and natural selection. Additionally, there is slight variation in the codon usage pattern among the SARS-CoV-2 isolates from different geo-locations. Furthermore, the overall codon usage pattern of SARS-CoV-2 is generally similar to that of its phylogenetic relatives among non-human betacoronaviruses such as RaTG13. Taken together, we comprehensively analyze the characteristics of codon usage pattern in SARS-CoV-2 via bioinformatic approaches. The information from this research may not only be helpful to get new insights into the evolution of SARS-CoV-2, but also have potential value for developing coronavirus vaccines.

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.

Does proteostasis get lost in translation? Implications for protein aggregation across the lifespan

Protein aggregation is a phenomenon of major relevance in neurodegenerative and neuromuscular disorders, cataracts, diabetes and many other diseases. Research has unveiled that proteins also aggregate in multiple tissues during healthy aging yet, the biological and biomedical relevance of this apparently asymptomatic phenomenon remains to be understood. It is known that proteome homeostasis (proteostasis) is maintained by a balanced protein synthesis rate, high protein synthesis accuracy, efficient protein folding and continual tagging of damaged proteins for degradation, suggesting that protein aggregation during healthy aging may be associated with alterations in both protein synthesis and the proteostasis network (PN) pathways. In particular, dysregulation of protein synthesis and alterations in translation fidelity are hypothesized to lead to the production of misfolded proteins which could explain the occurrence of age-related protein aggregation. Nevertheless, some data on this topic is controversial and the biological mechanisms that lead to widespread protein aggregation remain to be elucidated. We review the recent literature about the age-related decline of proteostasis, highlighting the need to build an integrated view of protein synthesis rate, fidelity and quality control pathways in order to better understand the proteome alterations that occur during aging and in age-related diseases.

Computational discovery and modeling of novel gene expression rules encoded in the mRNA

The transcript is populated with numerous overlapping codes that regulate all steps of gene expression. Deciphering these codes is very challenging due to the large number of variables involved, the non-modular nature of the codes, biases and limitations in current experimental approaches, our limited knowledge in gene expression regulation across the tree of life, and other factors. In recent years, it has been shown that computational modeling and algorithms can significantly accelerate the discovery of novel gene expression codes. Here, we briefly summarize the latest developments and different approaches in the field.

Effect of Protein Structure on Evolution of Cotranslational Folding

Cotranslational folding depends on the folding speed and stability of the nascent protein. It remains difficult, however, to predict which proteins cotranslationally fold. Here, we simulate evolution of model proteins to investigate how native structure influences evolution of cotranslational folding. We developed a model that connects protein folding during and after translation to cellular fitness. Model proteins evolved improved folding speed and stability, with proteins adopting one of two strategies for folding quickly. Low contact order proteins evolve to fold cotranslationally. Such proteins adopt native conformations early on during the translation process, with each subsequently translated residue establishing additional native contacts. On the other hand, high contact order proteins tend not to be stable in their native conformations until the full chain is nearly extruded. We also simulated evolution of slowly translating codons, finding that slower translation speeds at certain positions enhances cotranslational folding. Finally, we investigated real protein structures using a previously published data set that identified evolutionarily conserved rare codons in Escherichia coli genes and associated such codons with cotranslational folding intermediates. We found that protein substructures preceding conserved rare codons tend to have lower contact orders, in line with our finding that lower contact order proteins are more likely to fold cotranslationally. Our work shows how evolutionary selection pressure can cause proteins with local contact topologies to evolve cotranslational folding.

Interplay between Position-Dependent Codon Usage Bias and Hydrogen Bonding at the 5' End of ORFeomes

Codon usage bias exerts control over a wide variety of molecular processes. The positioning of synonymous codons within coding sequences (CDSs) dictates protein expression by mechanisms such as local translation efficiency, mRNA Gibbs free energy, and protein cotranslational folding. In this work, we explore how codon usage affects the position-dependent content of hydrogen bonding, which in turn influences energy requirements for unwinding double-stranded DNA (dsDNA). We categorized codons according to their hydrogen bond content and found differential effects on hydrogen bonding encoded by codon variants. The specific positional disposition of codon variants within CDSs creates a ramp of hydrogen bonding at the 5' end of the ORFeome in Escherichia coli CDSs occupying the first position of operons are subjected to selective pressure that reduces their hydrogen bonding compared to internal CDSs, and highly transcribed CDSs demand a lower maximum capacity of hydrogen bonds per codon, suggesting that the energetic requirement for unwinding the dsDNA in highly transcribed CDSs has evolved to be minimized in E. coli Subsequent analysis of over 14,000 ORFeomes showed a pervasive ramp of hydrogen bonding at the 5' end in Bacteria and Archaea that positively correlates with the probability of mRNA secondary structure formation. Both the ramp and the correlation were not found in Fungi The position-dependent hydrogen bonding might be part of the mechanism that contributes to the coordination between transcription and translation in Bacteria and Archaea A Web-based application to analyze the position-dependent hydrogen bonding of ORFeomes has been developed and is publicly available (https://juanvillada.shinyapps.io/hbonds/).IMPORTANCE Redundancy of the genetic code creates a vast space of alternatives to encode a protein. Synonymous codons exert control over a variety of molecular and physiological processes of cells mainly through influencing protein biosynthesis. Recent findings have shown that synonymous codon choice affects transcription by controlling mRNA abundance, mRNA stability, transcription termination, and transcript biosynthesis cost. In this work, by analyzing thousands of Bacteria, Archaea, and Fungi genomes, we extend recent findings by showing that synonymous codon choice, corresponding to the number of hydrogen bonds in a codon, can also have an effect on the energetic requirements for unwinding double-stranded DNA in a position-dependent fashion. This report offers new perspectives on the mechanism behind the transcription-translation coordination and complements previous hypotheses on the resource allocation strategies used by Bacteria and Archaea to manage energy efficiency in gene expression.

Synonymous Mutations of Porcine Igf1r Extracellular Domain Affect Differentiation and Mineralization in MC3T3-E1 Cells

Owing to the wide application of miniature pigs in biomedicine, the formation mechanism of its short stature must be elucidated. The insulin-like growth factor 1 receptor (IGF-1R), which receives signals through the extracellular domain (ECD) binding with ligands, is crucial in regulating cell growth and bone matrix mineralization. In this study, two haplotypes of Igf1r with four synonymous mutations in the coding sequences of IGF-1R ECD between large pigs (LP) and Bama pigs (BM) were stably expressed in the Igf1r-knockout MC3T3-E1 cells and named as MC3T3-LP cells (LP group) and MC3T3-BM cells (BM group), respectively. IGF-1R expression was lower in the BM group than in the LP group both in terms of transcription and translation levels, and IGF-1R expression inhibited cell proliferation. In addition, IGF-1R expression in the BM group promoted early-stage differentiation but delayed late-stage differentiation, which not only suppressed the expression of bone-related factors but also reduced alkaline phosphatase activity and calcium deposition. Moreover, different haplotypes of Igf1r changed the stability and conformation of the protein, further affecting the binding with IGF-1. Our data indicated that the four synonymous mutations of IGF1R ECD encoded by affect gene transcription and translation, thereby further leading to differences in the downstream pathways and functional changes of osteoblasts.

Effects of genetic variations in Acads gene on the risk of chronic obstructive pulmonary disease

Short-chain acyl-CoA dehydrogenase (SCAD), encoded by the Acads gene, functions in the mitochondrial β-oxidation of saturated short-chain fatty acids. SCAD deficiency results in mitochondrial dysfunction, which is one underlying biological mechanism of chronic obstructive pulmonary disease (COPD) pathogenesis. In this case-control study, we aimed to examine the effects of Acads gene polymorphisms on the susceptibility to COPD. A total of 16 tagging single-nucleotide polymorphisms (SNPs) in Acads gene region was identified and genotyped in 646 unrelated ethnic Chinese Han individuals including 279 patients with COPD and 367 healthy controls, their allelic and genotypic associations with COPD were determined by different genetic models. Furthermore, we estimated the linkage disequilibrium and haplotypes from these tested variants and determined the effects of haplotypes on the risk of COPD. The allelic and genotypic frequencies of SNPs rs2239686 and rs487915 in Acads gene were significantly different between COPD patients and controls, no statistically significant results were observed for other SNPs. Minor alleles A of rs2239686 and T of rs487915 were associated with a decreased pulmonary function and an increased COPD risk in a dominant manner. Functional analysis indicated that the risk allele A of rs2239686 could increase Acads expressions and the intracellular reactive oxygen species content. Haplotype analysis revealed that the haplotypes CTCCT in block 2 (rs3794216-rs3794215-rs34491494-rs558314-rs7312316) as well as GC in block 3 (rs2239686-rs487915) were protective against COPD, while haplotypes CTCGC in block 2 and AT in block 3 exhibited significant associations with the increased susceptibility to COPD. Our results suggest that Acads gene could potentially be a risk factor of COPD and thus its genetic variants might be as genetic biomarkers to predict the COPD susceptibility.

Comparative analysis of codon usage between Gossypium hirsutum and G. barbadense mitochondrial genomes

Gossypium hirsutum and G. barbadense mitochondrial genomes were analyzed to understand the factors shaping codon usage. While most analyses of codon usage suggest minimal to no bias, nucleotide composition, specifically GC content, was significantly correlated with codon usage. In general, both mitochondrial genomes favor codons that end in A or U, with a secondary preference for pyrimidine rich codons. These observations are similar to previous reports of codon usage in cotton nuclear genomes, possibly suggestive of a general bias spanning genomic compartment. Although evidence for codon usage bias is weak for most genes, we identified six genes (i.e. atp8, atp9, sdh3, sdh4, mttB and rpl2) with significant nonrandom codon usage. In general, we find multiple factors that influence cotton mitochondrial genome codon usage, which may include selection in a subset of genes.

Porcine IGF-1R synonymous mutations in the intracellular domain affect cell proliferation and alter kinase activity

Miniature pigs are regarded as ideal organ donors for xenotransplantation into humans. Elucidating the formation mechanism of miniature pigs is important. The insulin-like growth factor 1 receptor (IGF-1R) is crucial in the regulation of cell proliferation and organismal growth. According to our previous research, the IGF-1R expression levels between large and miniature pigs showed different profiles in liver and muscle tissues. Here, five synonymous mutations of IGF-1R in the coding sequence (CDS) of intracellular domain (ICD) between large and miniature pigs were analysed by constructing expression vectors of two haplotypes and named pcDNA3.1-LP (with the CDS of IGF-1R ICD of Large White pigs, LP group) and pcDNA3.1-BM (with the CDS of IGF-1R ICD of Bama Xiang pigs, BM group). The IGF-1R of the BM group was expressed lower than that of the LP group in transcription, translation and autophosphorylation levels. The IGF-1R of the BM group also down-regulated the protein levels of p-AKT/p-ERK than that of the LP group. PK-15 and C2C12 cell proliferation were detected to further understand the function of the haplotype. Results showed that the proliferation viability of PK-15 and C2C12 cells weakened in the BM group. Moreover, the mRNA and protein stabilities of the BM group were higher than those of the LP group. Our data indicated that two haplotypes of IGF-1R CDS in ICD between large and miniature pigs altered IGF-1R expression and down-regulated AKT and ERK signalling pathways at translation levels, resulting in an inhibitory effect on PK-15 and C2C12 cell proliferation.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum

BACKGROUND

Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts.

RESULTS

We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes.

CONCLUSIONS

These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Analysis of computational codon usage models and their association with translationally slow codons

Improved computational modeling of protein translation rates, including better prediction of where translational slowdowns along an mRNA sequence may occur, is critical for understanding co-translational folding. Because codons within a synonymous codon group are translated at different rates, many computational translation models rely on analyzing synonymous codons. Some models rely on genome-wide codon usage bias (CUB), believing that globally rare and common codons are the most informative of slow and fast translation, respectively. Others use the CUB observed only in highly expressed genes, which should be under selective pressure to be translated efficiently (and whose CUB may therefore be more indicative of translation rates). No prior work has analyzed these models for their ability to predict translational slowdowns. Here, we evaluate five models for their association with slowly translated positions as denoted by two independent ribosome footprint (RFP) count experiments from S. cerevisiae, because RFP data is often considered as a “ground truth” for translation rates across mRNA sequences. We show that all five considered models strongly associate with the RFP data and therefore have potential for estimating translational slowdowns. However, we also show that there is a weak correlation between RFP counts for the same genes originating from independent experiments, even when their experimental conditions are similar. This raises concerns about the efficacy of using current RFP experimental data for estimating translation rates and highlights a potential advantage of using computational models to understand translation rates instead.

Widespread non-modular overlapping codes in the coding regions

Messenger RNAs (mRNAs) consist of a coding region (open reading frame (ORF)) and two untranslated regions (UTRs), 5'UTR and 3'UTR. Ribosomes travel along the coding region, translating nucleotide triplets (called codons) to a chain of amino acids. The coding region was long believed to mainly encode the amino acid content of proteins, whereas regulatory signals reside in the UTRs and in other genomic regions. However, in recent years we have learned that the ORF is expansively populated with various regulatory signals, or codes, which are related to all gene expression steps and additional intracellular aspects. In this paper, we review the current knowledge related to overlapping codes inside the coding regions, such as the influence of synonymous codon usage on translation speed (and, in turn, the effect of translation speed on protein folding), ribosomal frameshifting, mRNA stability, methylation, splicing, transcription and more. All these codes come together and overlap in the ORF sequence, ensuring production of the right protein at the right time.

From SARS and MERS CoVs to SARS-CoV-2: Moving toward more biased codon usage in viral structural and nonstructural genes

Background

Severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) is an emerging disease with fatal outcomes. In this study, a fundamental knowledge gap question is to be resolved by evaluating the differences in biological and pathogenic aspects of SARS‐CoV‐2 and the changes in SARS‐CoV‐2 in comparison with the two prior major COV epidemics, SARS and Middle East respiratory syndrome (MERS) coronaviruses.

Methods

The genome composition, nucleotide analysis, codon usage indices, relative synonymous codons usage, and effective number of codons (ENc) were analyzed in the four structural genes; Spike (S), Envelope (E), membrane (M), and Nucleocapsid (N) genes, and two of the most important nonstructural genes comprising RNA‐dependent RNA polymerase and main protease (Mpro) of SARS‐CoV‐2, Beta‐CoV from pangolins, bat SARS, MERS, and SARS CoVs.

Results

SARS‐CoV‐2 prefers pyrimidine rich codons to purines. Most high‐frequency codons were ending with A or T, while the low frequency and rare codons were ending with G or C. SARS‐CoV‐2 structural proteins showed 5 to 20 lower ENc values, compared with SARS, bat SARS, and MERS CoVs. This implies higher codon bias and higher gene expression efficiency of SARS‐CoV‐2 structural proteins. SARS‐CoV‐2 encoded the highest number of over‐biased and negatively biased codons. Pangolin Beta‐CoV showed little differences with SARS‐CoV‐2 ENc values, compared with SARS, bat SARS, and MERS CoV.

Conclusion

Extreme bias and lower ENc values of SARS‐CoV‐2, especially in Spike, Envelope, and Mpro genes, are suggestive for higher gene expression efficiency, compared with SARS, bat SARS, and MERS CoVs.

Synonymous mutation rs2515641 affects CYP2E1 mRNA and protein expression and susceptibility to drug-induced liver injury

Aim: To evaluate whether the synonymous mutant rs2515641 could affect cytochrome P450 2E1 (CYP2E1) expression and the response to acetaminophen (APAP) or triptolide (TP) treatment. Materials & methods: HepG2 cells were transfected with lentiviral vector containing either CYP2E1-1263C or CYP2E1-1263T. Some of these recombinant cells were then treated with APAP or TP. CYP2E1 gene expression was detected by PCR and western blot. Results: CYP2E1 gene expression decreased significantly both in mRNA and protein level after rs2515641 mutation, indicating that this polymorphism can affect both transcription and translation. Furthermore, rs2515641 mutation dramatically changes the response of CYP2E1 expression to APAP or TP treatment. Conclusion: Rs2515641 significantly changes CYP2E1 expression and function, which would be expected to affect drug disposition and response.

Synonymous codon substitutions perturb cotranslational protein folding in vivo and impair cell fitness

In the cell, proteins are synthesized from N to C terminus and begin to fold during translation. Cotranslational folding mechanisms are therefore linked to elongation rate, which varies as a function of synonymous codon usage. However, synonymous codon substitutions can affect many distinct cellular processes, which has complicated attempts to deconvolve the extent to which synonymous codon usage can promote or frustrate proper protein folding in vivo. Although previous studies have shown that some synonymous changes can lead to different final structures, other substitutions will likely be more subtle, perturbing predominantly the protein folding pathway without radically altering the final structure. Here we show that synonymous codon substitutions encoding a single essential enzyme lead to dramatically slower cell growth. These mutations do not prevent active enzyme formation; instead, they predominantly alter the protein folding mechanism, leading to enhanced degradation in vivo. These results support a model in which synonymous codon substitutions can impair cell fitness by significantly perturbing cotranslational protein folding mechanisms, despite the chaperoning provided by the cellular protein homeostasis network.

Associations of noise kurtosis, genetic variations in NOX3 and lifestyle factors with noise-induced hearing loss

BACKGROUND

Noise-induced hearing loss (NIHL) is a complex disease caused by environmental and genetic risk factors. This study was to explore the association of noise kurtosis, triphosphopyridine nucleotide oxidase 3 (NOX3) and lifestyles with NIHL.

METHODS

This case-control study included 307 patients with NIHL and 307 matched control individuals from Zhejiang province of China. General characteristics, noise exposure data, the exfoliated cells of the oral mucosa, and lifestyle details of individuals were collected. The kompetitive allele specific polymerase chain reaction (KASP) method was used to analyze the genotypes of three single nucleotide polymorphisms (SNPs) of NOX3.

RESULTS

People who exposed to complex noise had a higher risk of NIHL than those exposed to steady noise (adjusted: OR = 1.806, P = 0.002). The GT genotype of additive model and TT + GT genotype of dominant model in NOX3 rs12195525 decreased the risk of NIHL (adjusted: OR = 0.618, P = 0.043; OR = 0.622, P = 0.036). Smoking and exposure to high video volume increased the risk of NIHL (adjusted: OR = 1.486, P = 0.038; OR = 1.611, P = 0.014). Oppositely, regular physical exercise decreased the risk of NIHL (adjusted: OR = 0.598, P = 0.004). A positive interaction was found between complex noise and lifestyles including high video volume exposure and no physical exercise in the additive models (RERI = 1.088, P < 0.001; RERI = 1.054, P = 0.024). A positive interaction was also found between NOX3 rs12195525 GG genotype and lifestyles including smoking and high video volume exposure in the additive models (RERI = 1.042, P = 0.005; RERI = 0.774, P = 0.044).

CONCLUSIONS

Noise temporal structure, NOX3 rs12195525 polymorphism, and the three lifestyles of smoking, video volume, and physical exercise were related to the NIHL. There were the interactions between noise temporal structure and the lifestyle of video volume or physical exercise, as well as between NOX3 and the lifestyle of smoking or video volume. These results provide a theoretical basis for the prevention and genetic testing of NIHL.

Association of SNPs in GnRH gene with sperm quality traits of Chinese water buffalo

Hypothalamic gonadotropin-releasing hormone (GnRH) controls the activity of hypothalamic-pituitary-gonadal axis and plays a key role in the reproductive performance of animals. In this study, five single nucleotide polymorphisms (SNPs), namely g.991T > C, g.1041T > C g.3424T > C, g.3462C > A and g.3463Inde A, were detected in the GnRH gene of 162 water buffaloes by Sanger sequencing. Each SNP was associated with more than two sperm quality traits of ejaculate volume, sperm concentration, post-thaw sperm motility and sperm abnormality. g.3424T > C and g.3462C > A were related to these four traits and had a remarkable effect on ejaculate volume. The three other SNPs were related to sperm concentration, post-thaw sperm motility and sperm abnormality. Moreover, six haplotypes (H1: TCCAI, H2: CTTC-, H3: TCCCI, H4: CTTA-, H5: CCTA- and H6: CTCC-) composed of five SNPs comprising seven different combined genotypes were generated by linkage disequilibrium analysis. Statistics followed by one-way ANOVA indicated that water buffaloes with the haplotype combination H1H1 had the highest genotypic frequency, and those with the H4H4 haplotype combination had the highest ejaculate volume. The sperm concentration of those with haplotype combination H1H5 was higher than that of the other genotypes. In summary, our study showed a remarkable association between the SNPs of GnRH and sperm quality traits of Chinese water buffalo.

Selection constraints determine preference for A/U-ending codons in Taxus contorta

Unequal utilization of synonymous codons is a well-known phenomenon among living organisms. This phenomenon plays a major role in the enhancement of the accuracy and efficiency of translation. Gymnosperms are rarely paid attention in this aspect. Understanding the degree of and determining the forces influencing codon usage bias (CUB) in Taxus contorta, an endangered Himalayan gymnosperm, will prove useful in interpreting the evolutionary characteristics of this species. Using RNAseq data, 93 790 assembled transcripts were clustered into 32 701 unigenes. Around 13 061 full-length sequences were utilized for the analysis of CUB. Compositional properties showed that GC-content ranged from 28.76% to 65.22%, with an average value of 44.28%, suggesting an AT-rich genome. The mean effective number of codons (ENC) value revealed that CUB is not strong in T. contorta. The preferred codons tended to be A/U ending, whereas the avoided codons tended to be G/C ending. A P2 index of 0.54 and a Mutation Responsive Index (MRI) value of -0.02 in addition to the results revealed by the neutrality, ENC, and parity plots showed that natural selection is a predominating factor governing CUB. Mutational pressure, gene length, hydropathiciy, aromaticity, and nucleotide composition influence CUB weakly.

Analysis of preferred codon usage in the coronavirus N genes and their implications for genome evolution and vaccine design

•

The nucleotide variations among the N genes of 13 different coronaviruses (CoVs) were interpreted.

•

Overall, 18 amino acids observed with varying preferred codons.

•

The effective number of codon values ranged from 40.43 to 53.85, revealing a slight codon bias.

•

A highly significant correlation between GC3s and ENc values was observed in porcine epidemic diarrhea CoV, followed by Middle East respiratory syndrome CoV.

The nucleocapsid (N) protein of a coronavirus plays a crucial role in virus assembly and in its RNA transcription. It is important to characterize a virus at the nucleotide level to discover the virus’s genomic sequence variations and similarities relative to other viruses that could have an impact on the functions of its genes and proteins. This entails a comprehensive and comparative analysis of the viral genomes of interest for preferred nucleotides, codon bias, nucleotide changes at the 3^rd position (NT3s), synonymous codon usage and relative synonymous codon usage. In this study, the variations in the N proteins among 13 different coronaviruses (CoVs) were analysed at the nucleotide and amino acid levels in an attempt to reveal how these viruses adapt to their hosts relative to their preferred codon usage in the N genes. The results revealed that, overall, eighteen amino acids had different preferred codons and eight of these were over-biased. The N genes had a higher AT% over GC% and the values of their effective number of codons ranged from 40.43 to 53.85, indicating a slight codon bias. Neutrality plots and correlation analyses showed a very high level of GC3s/GC correlation in porcine epidemic diarrhea CoV (pedCoV), followed by Middle East respiratory syndrome-CoV (MERS CoV), porcine delta CoV (dCoV), bat CoV (bCoV) and feline CoV (fCoV) with r values 0.81, 0.68, -0.47, 0.98 and 0.58, respectively. These data implied a high rate of evolution of the CoV genomes and a strong influence of mutation on evolutionary selection in the CoV N genes. This type of genetic analysis would be useful for evaluating a virus’s host adaptation, evolution and is thus of value to vaccine design strategies.

Progress toward Eukaryotic Semisynthetic Organisms: Translation of Unnatural Codons

We have created a bacterial semisynthetic organism (SSO) that retains an unnatural base pair (UBP) in its DNA, transcribes it into mRNA and tRNA with cognate unnatural codons and anticodons, and after the tRNA is charged with a noncanonical amino acid synthesizes proteins containing the noncanonical amino acid. Here, we report the first progress toward the creation of eukaryotic SSOs. After demonstrating proof-of-concept with human HEK293 cells, we show that a variety of different unnatural codon-anticodon pairs can efficiently mediate the synthesis of unnatural proteins in CHO cells. Interestingly, we find that there are both similarities and significant differences between how the prokaryotic and eukaryotic ribosomes recognize the UBP, with the eukaryotic ribosome appearing more tolerant. The results represent the first progress toward eukaryotic SSOs and, in fact, suggest that such SSOs might be able to retain more unnatural information than their bacterial counterparts.

DeepShape: estimating isoform-level ribosome abundance and distribution with Ribo-seq data

Background

Ribosome profiling brings insight to the process of translation. A basic step in profile construction at transcript level is to map Ribo-seq data to transcripts, and then assign a huge number of multiple-mapped reads to similar isoforms. Existing methods either discard the multiple mapped-reads, or allocate them randomly, or assign them proportionally according to transcript abundance estimated from RNA-seq data.

Results

Here we present DeepShape, an RNA-seq free computational method to estimate ribosome abundance of isoforms, and simultaneously compute their ribosome profiles using a deep learning model. Our simulation results demonstrate that DeepShape can provide more accurate estimations on both ribosome abundance and profiles when compared to state-of-the-art methods. We applied DeepShape to a set of Ribo-seq data from PC3 human prostate cancer cells with and without PP242 treatment. In the four cell invasion/metastasis genes that are translationally regulated by PP242 treatment, different isoforms show very different characteristics of translational efficiency and regulation patterns. Transcript level ribosome distributions were analyzed by “Codon Residence Index (CRI)” proposed in this study to investigate the relative speed that a ribosome moves on a codon compared to its synonymous codons. We observe consistent CRI patterns in PC3 cells. We found that the translation of several codons could be regulated by PP242 treatment.

Conclusion

In summary, we demonstrate that DeepShape can serve as a powerful tool for Ribo-seq data analysis.

Codon Usage in the Iflaviridae Family Is Not Diverse Though the Family Members Are Isolated from Diverse Host Taxa

All iflavirus members belong to the unique genus, Iflavirus, of the family, Iflaviridae. The host taxa and sequence identities of these viruses are diverse. A codon usage bias, maintained by a balance between selection, mutation, and genetic drift, exists in a wide variety of organisms. We characterized the codon usage patterns of 44 iflavirus genomes that were isolated from the classes, Insecta, Arachnida, Mammalia, and Malacostraca. Iflaviruses lack a strong codon usage bias when they are evaluated using an effective number of codons. The odds ratios of the majority of dinucleotides are within the normal range. However, the dinucleotides at the 1st-2nd codon positions are more biased than those at the 2nd-3rd codon positions. Plots of effective numbers of codons, relative neutrality analysis, and PR2 bias analysis all indicate that selection pressure dominates mutations in shaping codon usage patterns in the family, Iflaviridae. When these viruses were grouped into their host taxa, we found that the indices, including the nucleotide composition, effective number of codons, relative synonymous codon usage, and the influencing factors behind the codon usage patterns, all show that there are non-significant differences between the six host-taxa-groups. Our results disagree with our assumption that diverse viruses should possess diverse codon usage patterns, suggesting that the nucleotide composition and codon usage in the family, Iflaviridae, are not host taxa-specific signatures.

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.

Effects of codon optimization on coagulation factor IX translation and structure: Implications for protein and gene therapies

Synonymous codons occur with different frequencies in different organisms, a phenomenon termed codon usage bias. Codon optimization, a common term for a variety of approaches used widely by the biopharmaceutical industry, involves synonymous substitutions to increase protein expression. It had long been presumed that synonymous variants, which, by definition, do not alter the primary amino acid sequence, have no effect on protein structure and function. However, a critical mass of reports suggests that synonymous codon variations may impact protein conformation. To investigate the impact of synonymous codons usage on protein expression and function, we designed an optimized coagulation factor IX (FIX) variant and used multiple methods to compare its properties to the wild-type FIX upon expression in HEK293T cells. We found that the two variants differ in their conformation, even when controlling for the difference in expression levels. Using ribosome profiling, we identified robust changes in the translational kinetics of the two variants and were able to identify a region in the gene that may have a role in altering the conformation of the protein. Our data have direct implications for codon optimization strategies, for production of recombinant proteins and gene therapies.

SDHA gain-of-function engages inflammatory mitochondrial retrograde signaling via KEAP1-Nrf2

Whether screening the metabolic activity of immune cells facilitates discovery of molecular pathology remains unknown. Here we prospectively screened the extracellular acidification rate as a measure of glycolysis and the oxygen consumption rate as a measure of mitochondrial respiration in B cells from patients with primary antibody deficiency. The highest oxygen consumption rate values were detected in three study participants with persistent polyclonal B cell lymphocytosis (PPBL). Exome sequencing identified germline mutations in SDHA, which encodes succinate dehydrogenase subunit A, in all three patients with PPBL. SDHA gain-of-function led to an accumulation of fumarate in PPBL B cells, which engaged the KEAP1-Nrf2 system to drive the transcription of genes encoding inflammatory cytokines. In a single patient trial, blocking the activity of the cytokine interleukin-6 in vivo prevented systemic inflammation and ameliorated clinical disease. Overall, our study has identified pathological mitochondrial retrograde signaling as a disease modifier in primary antibody deficiency.

Protein surface charge effect on 3D domain swapping in cells for c-type cytochromes

Many c-type cytochromes (cyts) can form domain-swapped oligomers. The positively charged Hydrogenobacter thermophilus (HT) cytochrome (cyt) c₅₅₂ forms domain-swapped oligomers during expression in the Escherichia coli (E. coli) expression system, but the factors influencing the oligomerization remain unrevealed. Here, we found that the dimer of the negatively charged Shewanella violacea (SV) cyt c₅ exhibits a domain-swapped structure, in which the N-terminal helix is exchanged between protomers, similar to the structures of the HT cyt c₅₅₂ and Pseudomonas aeruginosa (PA) cyt c₅₅₁ domain-swapped dimers. Positively charged horse cyt c and HT cyt c₅₅₂ domain swapped during expression in E. coli, whereas negatively charged PA cyt c₅₅₁ and SV cyt c₅ did not. Oligomers were formed during expression in E. coli for HT cyt c₅₅₂ attached to either a co- or post-translational signal peptide for transportation through the cytoplasm membrane, but not for PA cyt c₅₅₁ attached to either signal peptide. HT cyt c₅₅₂ formed oligomers in E. coli in the presence and absence of rare codons. More oligomers were obtained from the in vitro folding of horse cyt c and HT cyt c₅₅₂ by the addition of negatively charged liposomes during folding, whereas the amount of oligomers for the in vitro folding of PA cyt c₅₅₁ and SV cyt c₅ did not change significantly by the addition. These results indicate that the protein surface charge affects the oligomerization of c-type cyts in cells; positively charged c-type cyts assemble on a negatively charged membrane, inducing formation of domain-swapped oligomers during folding.

Decomposition of time-dependent fluorescence signals reveals codon-specific kinetics of protein synthesis

During protein synthesis, the nascent peptide chain traverses the peptide exit tunnel of the ribosome. We monitor the co-translational movement of the nascent peptide using a fluorescent probe attached to the N-terminus of the nascent chain. Due to fluorophore quenching, the time-dependent fluorescence signal emitted by an individual peptide is determined by co-translational events, such as secondary structure formation and peptide-tunnel interactions. To obtain information on these individual events, the measured ensemble fluorescence signal has to be decomposed into position-dependent intensities. Here, we describe mRNA translation as a Markov process with specific fluorescence intensities assigned to the different states of the process. Combining the computed stochastic time evolution of the translation process with a sequence of observed ensemble fluorescence time courses, we compute the unknown position-specific intensities and obtain detailed information on the kinetics of the translation process. In particular, we find that translation of poly(U) mRNAs dramatically slows down at the fourth UUU codon. The method presented here detects subtle differences in the position-specific fluorescence intensities and thus provides a novel approach to study translation kinetics in ensemble experiments.