Codon usage and tRNA genes in eukaryotes: correlation of codon usage diversity with translation efficiency and with CG-dinucleotide usage as assessed by multivariate analysis

Antibody production relies on the tRNA inosine wobble modification to meet biased codon demand

Antibodies are produced at high rates to provide immunoprotection, which puts pressure on the B cell translational machinery. Here, we identified a pattern of codon usage conserved across antibody genes. One feature thereof is the hyperutilization of codons that lack genome-encoded Watson-Crick transfer RNAs (tRNAs), instead relying on the posttranscriptional tRNA modification inosine (I34), which expands the decoding capacity of specific tRNAs through wobbling. Antibody-secreting cells had increased I34 levels and were more reliant on I34 for protein production than naïve B cells. Furthermore, antibody I34-dependent codon usage may influence B cell passage through regulatory checkpoints. Our work elucidates the interface between the tRNA pool and protein production in the immune system and has implications for the design and selection of antibodies for vaccines and therapeutics.

Construction of pseudorabies virus variant attenuated vaccine: codon deoptimization of US3 and UL56 genes based on PRV gE/TK deletion strain

Since 2011, pseudorabies based on the pseudorabies virus (PRV) variant has emerged as a serious health issue in pig farms in China. The PRV gE/TK or gE/gI/TK deletion strains protect against emerging PRV variants. However, these variants may cause lethal infections in newborn piglets without PRV antibodies. Previous studies have shown that codon deoptimization of a virulence gene causes virus attenuation. Accordingly, we deoptimized US3-S (US3 gene encoding a short isoform that represents approximately 95% of the total US3 transcription) and UL56 genes (first 10 or all codons) of PRV gE/TK deletion strain (PRVΔTK&gE-AH02) to generate six recombinant PRVs through bacterial artificial chromosome technology. In swine testicular cells, recombinant PRVs with all codon deoptimization of US3-S or UL56 genes were grown to lower titers than the parental virus. Notably, US3-S or UL56 with all codon deoptimization reduced mRNA and protein expressions. Subsequently, the safety and immunogenicity of recombinant PRVs with codon deoptimization of US3-S or UL56 are evaluated as vaccine candidates in mice and piglets. The mice inoculated with recombinant PRVs with codon deoptimization of US3-S or UL56 showed exceptional survival ability without severe clinical signs. All codons deoptimized (US3-S and UL56) significantly decreased virus load and attenuated pathological changes in the brains of the mice. Moreover, the protection efficiency offered by recombinant PRVs with codon deoptimization of US3-S or UL56 showed similar effects to PRVΔTK&gE-AH02. Remarkably, the 1-day-old PRV antibody-negative piglets inoculated with PRVΔTK&gE-US3-ST-CD (a recombinant PRV with all codon deoptimization of US3-S) presented no abnormal clinical symptoms, including fever. The piglets inoculated with PRVΔTK&gE-US3-ST-CD showed a high serum neutralization index against the PRV variant. In conclusion, these results suggest using codon deoptimization to generate innovative live attenuated PRV vaccine candidates.

Dynamic Landscapes of tRNA Transcriptomes and Translatomes in Diverse Mouse Tissues

Although the function of tRNAs in the translational process is well established, it remains controversial whether tRNA abundance is tightly associated with translational efficiency (TE) in mammals. Moreover, how critically the expression of tRNAs contributes to the establishment of tissue-specific proteomes in mammals has not been well addressed. Here, we measured both tRNA expression using demethylase-tRNA sequencing (DM-tRNA-seq) and TE of mRNAs using ribosome-tagging sequencing (RiboTag-seq) in the brain, heart, and testis of mice. Remarkable variation in the expression of tRNA isodecoders was observed among different tissues. When the statistical effect of isodecoder-grouping on reducing variations is considered through permutating the anticodons, we observed an expected reduction in the variation of anticodon expression across all samples, an unexpected smaller variation of anticodon usage bias, and an unexpected larger variation of tRNA isotype expression at amino acid level. Regardless of whether or not they share the same anticodons, the isodecoders encoding the same amino acids are co-expressed across different tissues. Based on the expression of tRNAs and the TE of mRNAs, we find that the tRNA adaptation index (tAI) and TE are significantly correlated in the same tissues but not between tissues; and tRNA expression and the amino acid composition of translating peptides are positively correlated in the same tissues but not between tissues. We therefore hypothesize that the tissue-specific expression of tRNAs might be due to post-transcriptional mechanisms. This study provides a resource for tRNA and translation studies, as well as novel insights into the dynamics of tRNAs and their roles in translational regulation.

Codon Usage Bias Analysis in Macronuclear Genomes of Ciliated Protozoa

Ciliated protozoa (ciliates) are unicellular eukaryotes, several of which are important model organisms for molecular biology research. Analyses of codon usage bias (CUB) of the macronuclear (MAC) genome of ciliates can promote a better understanding of the genetic mode and evolutionary history of these organisms and help optimize codons to improve gene editing efficiency in model ciliates. In this study, the following indices were calculated: the guanine-cytosine (GC) content, the frequency of the nucleotides at the third position of codons (T3, C3, A3, G3), the effective number of codons (ENc), GC content at the 3rd position of synonymous codons (GC3s), and the relative synonymous codon usage (RSCU). Parity rule 2 plot analysis, Neutrality plot analysis, ENc plot analysis, and correlation analysis were employed to explore the main influencing factors of CUB. The results showed that the GC content in the MAC genomes of each of 21 ciliate species, the genomes of which were relatively complete, was lower than 50%, and the base compositions of GC and GC3s were markedly distinct. Synonymous codon analysis revealed that the codons in most of the 21 ciliates ended with A or T and four codons were the general putative optimal codons. Collectively, our results indicated that most of the ciliates investigated preferred using the codons with anof AT-ending and that codon usage bias was affected by gene mutation and natural selection.

In silico methods for predicting functional synonymous variants

Single nucleotide variants (SNVs) contribute to human genomic diversity. Synonymous SNVs are previously considered to be "silent," but mounting evidence has revealed that these variants can cause RNA and protein changes and are implicated in over 85 human diseases and cancers. Recent improvements in computational platforms have led to the development of numerous machine-learning tools, which can be used to advance synonymous SNV research. In this review, we discuss tools that should be used to investigate synonymous variants. We provide supportive examples from seminal studies that demonstrate how these tools have driven new discoveries of functional synonymous SNVs.

Natural selection pressure exerted on "Silent" mutations during the evolution of SARS-CoV-2: Evidence from codon usage and RNA structure

From the first emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) till now, multiple mutations that caused synonymous and nonsynonymous substitutions have accumulated. Among them, synonymous substitutions were regarded as "silent" mutations that received less attention than nonsynonymous substitutions that cause amino acid variations. However, the importance of synonymous substitutions can not be neglected. This research focuses on synonymous substitutions on SARS-CoV-2 and proves that synonymous substitutions were under purifying selection in its evolution. The evidence of purifying selection is provided by comparing the mutation number per site in coding and non-coding regions. We then study the two forces of purifying selection: synonymous codon usage and RNA secondary structure. Results show that the codon usage optimization leads to an adapted codon usage towards humans. Furthermore, our results show that the maintenance of RNA secondary structure causes the purifying of synonymous substitutions in the structural region. These results explain the selection pressure on synonymous substitutions during the evolution of SARS-CoV-2.

Codon usage of host-specific P genotypes (VP4) in group A rotavirus

Background

Group A rotavirus (RVA) is a common causative agent of acute gastroenteritis in infants and young children worldwide. RVA P genotypes, determined by VP4 sequences, have been confirmed to infect humans and animals. However, their codon usage patterns that are essential to obtain insights into the viral evolution, host adaptability, and genetic characterization remained unclear, especially across animal hosts.

Results

We performed a comprehensive codon usage analysis of eight host-specific RVA P genotypes, including human RVA (P[4] and P[8]), porcine RVA (P[13] and P[23]), and zoonotic RVA (P[1], P[6], P[7] and P[19]), based on 233 VP4 complete coding sequences. Nucleotide composition, relative synonymous codon usage (RSCU), and effective number of codons (ENC) were calculated. Principal component analysis (PCA) based on RSCU values was used to explore the codon usage patterns of different RVA P genotypes. In addition, mutation pressure and natural selection were identified by using ENC-plot, parity rule 2 plot, and neutrality plot analyses. All VP4 sequences preferred using A/U nucleotides (A: 0.354-0.377, U: 0.267-0.314) than G/C nucleotides across genotypes. Similarly, majority of commonly used synonymous codons were likely to end with A/U nucleotides (A: 9/18-12/18, U: 6/18-9/18). In PCA, human, porcine, and zoonotic genotypes clustered separately in terms of RSCU values, indicating the host-specific codon usage patterns; however, porcine and zoonotic genotypes were partly overlapped. Human genotypes, P[4] and P[8], had stronger codon usage bias, as indicated by more over-represented codons and lower ENC, compared to porcine and zoonotic genotypes. Moreover, natural selection was determined to be a predominant driver in shaping the codon usage bias across the eight P genotypes. In addition, mutation pressure contributed to the codon usage bias of human genotypes.

Conclusions

Our study identified a strong codon usage bias of human RVA P genotypes attributable to both natural selection and mutation pressure, whereas similar codon usage bias between porcine and zoonotic genotypes predominantly attributable to natural selection. It further suggests possible cross-species transmission. Therefore, it warrants further surveillance of RVA P genotypes for early identification of zoonotic infection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08730-2.

Vaccinia Virus Attenuation by Codon Deoptimization of the A24R Gene for Vaccine Development

Poxviruses have large DNA genomes, and they are able to infect multiple vertebrate and invertebrate animals, including humans. Despite the eradication of smallpox, poxvirus infections still remain a significant public health concern. Vaccinia virus (VV) is the prototypic member in the poxviridae family and it has been used extensively for different prophylactic applications, including the generation of vaccines against multiple infectious diseases and/or for oncolytic treatment. Many attempts have been pursued to develop novel attenuated forms of VV with improved safety profiles for their implementation as vaccines and/or vaccines vectors. We and others have previously demonstrated how RNA viruses encoding codon-deoptimized viral genes are attenuated, immunogenic and able to protect, upon a single administration, against challenge with parental viruses. In this study, we employed the same experimental approach based on the use of misrepresented codons for the generation of a recombinant (r)VV encoding a codon-deoptimized A24R gene, which is a key component of the viral RNA polymerase. Similar to our previous studies with RNA viruses, the A24R codon-deoptimized rVV (v-A24cd) was highly attenuated in vivo but able to protect, after a single intranasal dose administration, against an otherwise lethal challenge with parental VV. These results indicate that poxviruses can be effectively attenuated by synonymous codon deoptimization and open the possibility of using this methodology alone or in combination with other experimental approaches for the development of attenuated vaccines for the treatment of poxvirus infection, or to generate improved VV-based vectors. Moreover, this approach could be applied to other DNA viruses. IMPORTANCE The family poxviridae includes multiple viruses of medical and veterinary relevance, being vaccinia virus (VV) the prototypic member in the family. VV was used during the smallpox vaccination campaign to eradicate variola virus (VARV), which is considered a credible bioterrorism threat. Because of novel innovations in genetic engineering and vaccine technology, VV has gained popularity as a viral vector for the development of vaccines against several infectious diseases. Several approaches have been used to generate attenuated VV for its implementation as vaccine and/or vaccine vector. Here, we generated a rVV containing a codon-deoptimized A24R gene (v-A24cd), which encodes a key component of the viral RNA polymerase. v-A24cd was stable in culture cells and highly attenuated in vivo but able to protect against a subsequent lethal challenge with parental VV. Our findings support the use of this approach for the development of safe, stable, and protective live-attenuated VV and/or vaccine vectors.

Three-nucleotide periodicity of nucleotide diversity in a population enables the identification of open reading frames

Accurate prediction of open reading frames (ORFs) is important for studying and using genome sequences. Ribosomes move along mRNA strands with a step of three nucleotides and datasets carrying this information can be used to predict ORFs. The ribosome-protected footprints (RPFs) feature a significant 3-nt periodicity on mRNAs and are powerful in predicting translating ORFs, including small ORFs (sORFs), but the application of RPFs is limited because they are too short to be accurately mapped in complex genomes. In this study, we found a significant 3-nt periodicity in the datasets of populational genomic variants in coding sequences, in which the nucleotide diversity increases every three nucleotides. We suggest that this feature can be used to predict ORFs and develop the Python package ‘OrfPP’, which recovers ~83% of the annotated ORFs in the tested genomes on average, independent of the population sizes and the complexity of the genomes. The novel ORFs, including sORFs, identified from single-nucleotide polymorphisms are supported by protein mass spectrometry evidence comparable to that of the annotated ORFs. The application of OrfPP to tetraploid cotton and hexaploid wheat genomes successfully identified 76.17% and 87.43% of the annotated ORFs in the genomes, respectively, as well as 4704 sORFs, including 1182 upstream and 2110 downstream ORFs in cotton and 5025 sORFs, including 232 upstream and 234 downstream ORFs in wheat. Overall, we propose an alternative and supplementary approach for ORF prediction that can extend the studies of sORFs to more complex genomes.

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.

Intragenic codon usage in proteobacteria: Translational selection, IS expansion and genomic shrinkage

This manuscript presents a method to systematically study intragenic variations in codon usage using correspondence analysis and the effective number of codons. The method is applied to >1100 proteobacteria. Codon usage biases (measured as inertia) increases with genome size, the same is true for the percentage of inertia explained by the first axis. It is shown that there is often a relaxed or more uniform codon usage near the gene termini. Ithis is not seen n small genomes, notably those of intracellular organisms like Buchnera aphidicola or Rickettsia prowazekii where translational selection plays less of a role. When genes from E. coli, for which translational selection is well described, are split into low, intermediate and high expression, respectively, it is shown that the intragenic codon usage pattern with more uniform usage at termini exist across all three expression groups. Furthermore, the correspondence analysis reveals a unique pattern in Bordetella pertussis due to IS expansion. This study thus shows that translational selection, genome shrinkage and IS expansion result in characteristic patterns in intragenic codon usage.

CAMAP: Artificial neural networks unveil the role of codon arrangement in modulating MHC-I peptides presentation

MHC-I associated peptides (MAPs) play a central role in the elimination of virus-infected and neoplastic cells by CD8 T cells. However, accurately predicting the MAP repertoire remains difficult, because only a fraction of the transcriptome generates MAPs. In this study, we investigated whether codon arrangement (usage and placement) regulates MAP biogenesis. We developed an artificial neural network called Codon Arrangement MAP Predictor (CAMAP), predicting MAP presentation solely from mRNA sequences flanking the MAP-coding codons (MCCs), while excluding the MCC per se. CAMAP predictions were significantly more accurate when using original codon sequences than shuffled codon sequences which reflect amino acid usage. Furthermore, predictions were independent of mRNA expression and MAP binding affinity to MHC-I molecules and applied to several cell types and species. Combining MAP ligand scores, transcript expression level and CAMAP scores was particularly useful to increase MAP prediction accuracy. Using an in vitro assay, we showed that varying the synonymous codons in the regions flanking the MCCs (without changing the amino acid sequence) resulted in significant modulation of MAP presentation at the cell surface. Taken together, our results demonstrate the role of codon arrangement in the regulation of MAP presentation and support integration of both translational and post-translational events in predictive algorithms to ameliorate modeling of the immunopeptidome.

MHC-I associated peptides (MAPs) are small fragments of intracellular proteins presented at the surface of cells and used by the immune system to detect and eliminate cancerous or virus-infected cells. While it is theoretically possible to predict which portions of the intracellular proteins will be naturally processed by the cells to ultimately reach the surface, current methodologies have prohibitively high false discovery rates. Here we introduce an artificial neural network called Codon Arrangement MAP Predictor (CAMAP) which integrates information from mRNA-to-protein translation to other factors regulating MAP biogenesis (e.g. MAP ligand score and transcript expression levels) to improve MAP prediction accuracy. While most MAP predictive approaches focus on MAP sequences per se, CAMAP’s novelty is to analyze the MAP-flanking mRNA sequences, thereby providing completely independent information for MAP prediction. We show on several datasets that the integration of CAMAP scores with other known factors involved in MAP presentation (i.e. MAP ligand score and mRNA expression) significantly improves MAP prediction accuracy, and further validate CAMAP learned features using an in-vitro assay. These findings may have major implications for the design of vaccines against cancers and viruses, and in times of pandemics could accelerate the identification of relevant MAPs of viral origins.

Compelling Evidence Suggesting the Codon Usage of SARS-CoV-2 Adapts to Human After the Split From RaTG13

SARS-CoV-2 needs to efficiently make use of the resources from hosts in order to survive and propagate. Among the multiple layers of regulatory network, mRNA translation is the rate-limiting step in gene expression. Synonymous codon usage usually conforms with tRNA concentration to allow fast decoding during translation. It is acknowledged that SARS-CoV-2 has adapted to the codon usage of human lungs so that the virus could rapidly proliferate in the lung environment. While this notion seems to nicely explain the adaptation of SARS-CoV-2 to lungs, it is unable to tell why other viruses do not have this advantage. In this study, we retrieve the GTEx RNA-seq data for 30 tissues (belonging to over 17 000 individuals). We calculate the RSCU (relative synonymous codon usage) weighted by gene expression in each human sample, and investigate the correlation of RSCU between the human tissues and SARS-CoV-2 or RaTG13 (the closest coronavirus to SARS-CoV-2). Lung has the highest correlation of RSCU to SARS-CoV-2 among all tissues, suggesting that the lung environment is generally suitable for SARS-CoV-2. Interestingly, for most tissues, SARS-CoV-2 has higher correlations with the human samples compared with the RaTG13-human correlation. This difference is most significant for lungs. In conclusion, the codon usage of SARS-CoV-2 has adapted to human lungs to allow fast decoding and translation. This adaptation probably took place after SARS-CoV-2 split from RaTG13 because RaTG13 is less perfectly correlated with human. This finding depicts the trajectory of adaptive evolution from ancestral sequence to SARS-CoV-2, and also well explains why SARS-CoV-2 rather than other viruses could perfectly adapt to human lung environment.

Variability in Codon Usage in Coronaviruses Is Mainly Driven by Mutational Bias and Selective Constraints on CpG Dinucleotide

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third human-emerged virus of the 21st century from the Coronaviridae family, causing the ongoing coronavirus disease 2019 (COVID-19) pandemic. Due to the high zoonotic potential of coronaviruses, it is critical to unravel their evolutionary history of host species breadth, host-switch potential, adaptation and emergence, to identify viruses posing a pandemic risk in humans. We present here a comprehensive analysis of the composition and codon usage bias of the 82 Orthocoronavirinae members, infecting 47 different avian and mammalian hosts. Our results clearly establish that synonymous codon usage varies widely among viruses, is only weakly dependent on their primary host, and is dominated by mutational bias towards AU-enrichment and by CpG avoidance. Indeed, variation in GC3 explains around 34%, while variation in CpG frequency explains around 14% of total variation in codon usage bias. Further insight on the mutational equilibrium within Orthocoronavirinae revealed that most coronavirus genomes are close to their neutral equilibrium, the exception being the three recently infecting human coronaviruses, which lie further away from the mutational equilibrium than their endemic human coronavirus counterparts. Finally, our results suggest that, while replicating in humans, SARS-CoV-2 is slowly becoming AU-richer, likely until attaining a new mutational equilibrium.

Comparative Analysis of Genomic and Transcriptome Sequences Reveals Divergent Patterns of Codon Bias in Wheat and Its Ancestor Species

The synonymous codons usage shows a characteristic pattern of preference in each organism. This codon usage bias is thought to have evolved for efficient protein synthesis. Synonymous codon usage was studied in genes of the hexaploid wheat Triticum aestivum (AABBDD) and its progenitor species, Triticum urartu (AA), Aegilops tauschii (DD), and Triticum turgidum (AABB). Triticum aestivum exhibited stronger usage bias for G/C-ending codons than did the three progenitor species, and this bias was especially higher compared to T. turgidum and Ae. tauschii. High GC content is a primary factor influencing codon usage in T. aestivum. Neutrality analysis showed a significant positive correlation (p<0.001) between GC12 and GC3 in the four species with regression line slopes near zero (0.16–0.20), suggesting that the effect of mutation on codon usage was only 16–20%. The GC3s values of genes were associated with gene length and distribution density within chromosomes. tRNA abundance data indicated that codon preference corresponded to the relative abundance of isoaccepting tRNAs in the four species. Both mutation and selection have affected synonymous codon usage in hexaploid wheat and its progenitor species. GO enrichment showed that GC biased genes were commonly enriched in physiological processes such as photosynthesis and response to acid chemical. In some certain gene families with important functions, the codon usage of small parts of genes has changed during the evolution process of T. aestivum.

Codon Usage Bias: An Endless Tale

Since the genetic code is degenerate, several codons are translated to the same amino acid. Although these triplets were historically considered to be "synonymous" and therefore expected to be used at rather equal frequencies in all genomes, we now know that this is not the case. Indeed, since several coding sequences were obtained in the late '70s and early '80s in the last century, coming from either the same or different species, it was evident that (a) each genome, taken globally, displayed different codon usage patterns, which means that different genomes display a particular global codon usage table when all genes are considered together, and (b) there is a strong intragenomic diversity: in other words, within a given species the codon usage pattern can (and usually do) differ greatly among genes in the same genome. These different patterns were attributed to two main factors: first, the mutational bias characteristic of each genome, which determines that GC- poor species display a general bias towards A/T codons while the reverse is true for GC- rich species. Second, the differences in codon usage among genes from the same species are due to natural selection acting at the level of translation, in such a way that highly expressed genes tend to use codons that match with the most abundant isoacceptor tRNAs. Thus, these genes are translated at a highest rate, which in turn leads to avoid the limiting factor in translation which is the number of available ribosomes per cell. Although these explanations are still valid, new factors are almost constantly postulated to affect codon usage. In this mini review, we shall try to summarize them.

Role of a cryptic tRNA gene operon in survival under translational stress

As compared to eukaryotes, bacteria have a reduced tRNA gene set encoding between 30 and 220 tRNAs. Although in most bacterial phyla tRNA genes are dispersed in the genome, many species from distinct phyla also show genes forming arrays. Here, we show that two types of arrays with distinct evolutionary origins exist. This work focuses on long tRNA gene arrays (L-arrays) that encompass up to 43 genes, which disseminate by horizontal gene transfer and contribute supernumerary tRNA genes to the host. Although in the few cases previously studied these arrays were reported to be poorly transcribed, here we show that the L-array of the model cyanobacterium Anabaena sp. PCC 7120, encoding 23 functional tRNAs, is largely induced upon impairment of the translation machinery. The cellular response to this challenge involves a global reprogramming of the transcriptome in two phases. tRNAs encoded in the array are induced in the second phase of the response, directly contributing to cell survival. Results presented here show that in some bacteria the tRNA gene set may be partitioned between a housekeeping subset, which constantly sustains translation, and an inducible subset that is generally silent but can provide functionality under particular conditions.

Rotavirus A Genome Segments Show Distinct Segregation and Codon Usage Patterns

Reassortment of the Rotavirus A (RVA) 11-segment dsRNA genome may generate new genome constellations that allow RVA to expand its host range or evade immune responses. Reassortment may also produce phylogenetic incongruities and weakly linked evolutionary histories across the 11 segments, obscuring reassortment-specific epistasis and changes in substitution rates. To determine the co-segregation patterns of RVA segments, we generated time-scaled phylogenetic trees for each of the 11 segments of 789 complete RVA genomes isolated from mammalian hosts and compared the segments’ geodesic distances. We found that segments 4 (VP4) and 9 (VP7) occupied significantly different tree spaces from each other and from the rest of the genome. By contrast, segments 10 and 11 (NSP4 and NSP5/6) occupied nearly indistinguishable tree spaces, suggesting strong co-segregation. Host-species barriers appeared to vary by segment, with segment 9 (VP7) presenting the weakest association with host species. Bayesian Skyride plots were generated for each segment to compare relative genetic diversity among segments over time. All segments showed a dramatic decrease in diversity around 2007 coinciding with the introduction of RVA vaccines. To assess selection pressures, codon adaptation indices and relative codon deoptimization indices were calculated with respect to different host genomes. Codon usage varied by segment with segment 11 (NSP5) exhibiting significantly higher adaptation to host genomes. Furthermore, RVA codon usage patterns appeared optimized for expression in humans and birds relative to the other hosts examined, suggesting that translational efficiency is not a barrier in RVA zoonosis.

Whole-Genome Transformation Promotes tRNA Anticodon Suppressor Mutations under Stress

tRNAs are encoded by a large gene family, usually with several isogenic tRNAs interacting with the same codon. Mutations in the anticodon region of other tRNAs can overcome specific tRNA deficiencies. Phylogenetic analysis suggests that such mutations have occurred in evolution, but the driving force is unclear. We show that in yeast suppressor mutations in other tRNAs are able to overcome deficiency of the essential TRT2-encoded tRNAThrCGU at high temperature (40°C). Surprisingly, these tRNA suppressor mutations were obtained after whole-genome transformation with DNA from thermotolerant Kluyveromyces marxianus or Ogataea polymorpha strains but from which the mutations did apparently not originate. We suggest that transient presence of donor DNA in the host facilitates proliferation at high temperature and thus increases the chances for occurrence of spontaneous mutations suppressing defective growth at high temperature. Whole-genome sequence analysis of three transformants revealed only four to five nonsynonymous mutations of which one causing TRT2 anticodon stem stabilization and two anticodon mutations in non-threonyl-tRNAs, tRNALysCUU and tRNAeMetCAU, were causative. Both anticodon mutations suppressed lethality of TRT2 deletion and apparently caused the respective tRNAs to become novel substrates for threonyl-tRNA synthetase. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) data could not detect any significant mistranslation, and reverse transcription-quantitative PCR results contradicted induction of the unfolded protein response. We suggest that stress conditions have been a driving force in evolution for the selection of anticodon-switching mutations in tRNAs as revealed by phylogenetic analysis.IMPORTANCE In this work, we have identified for the first time the causative elements in a eukaryotic organism introduced by applying whole-genome transformation and responsible for the selectable trait of interest, i.e., high temperature tolerance. Surprisingly, the whole-genome transformants contained just a few single nucleotide polymorphisms (SNPs), which were unrelated to the sequence of the donor DNA. In each of three independent transformants, we have identified a SNP in a tRNA, either stabilizing the essential tRNAThrCGU at high temperature or switching the anticodon of tRNALysCUU or tRNAeMetCAU into CGU, which is apparently enough for in vivo recognition by threonyl-tRNA synthetase. LC-MS/MS analysis indeed indicated absence of significant mistranslation. Phylogenetic analysis showed that similar mutations have occurred throughout evolution and we suggest that stress conditions may have been a driving force for their selection. The low number of SNPs introduced by whole-genome transformation may favor its application for improvement of industrial yeast strains.

Genome Size Reduction and Transposon Activity Impact tRNA Gene Diversity While Ensuring Translational Stability in Birds

As a highly diverse vertebrate class, bird species have adapted to various ecological systems. How this phenotypic diversity can be explained genetically is intensively debated and is likely grounded in differences in the genome content. Larger and more complex genomes could allow for greater genetic regulation that results in more phenotypic variety. Surprisingly, avian genomes are much smaller compared to other vertebrates but contain as many protein-coding genes as other vertebrates. This supports the notion that the phenotypic diversity is largely determined by selection on non-coding gene sequences. Transfer RNAs (tRNAs) represent a group of non-coding genes. However, the characteristics of tRNA genes across bird genomes have remained largely unexplored. Here, we exhaustively investigated the evolution and functional consequences of these crucial translational regulators within bird species and across vertebrates. Our dense sampling of 55 avian genomes representing each bird order revealed an average of 169 tRNA genes with at least 31% being actively used. Unlike other vertebrates, avian tRNA genes are reduced in number and complexity but are still in line with vertebrate wobble pairing strategies and mutation-driven codon usage. Our detailed phylogenetic analyses further uncovered that new tRNA genes can emerge through multiplication by transposable elements. Together, this study provides the first comprehensive avian and cross-vertebrate tRNA gene analyses and demonstrates that tRNA gene evolution is flexible albeit constrained within functional boundaries of general mechanisms in protein translation.

Efficient expression of a cnidarian peptide-gated ion channel in mammalian cells

Hydra Na⁺ channels (HyNaCs) are peptide-gated ion channels of the DEG/ENaC gene family that are directly activated by neuropeptides of the Hydra nervous system. They have previously been successfully characterized in Xenopus oocytes. To establish their expression in mammalian cells, we transiently expressed heteromeric HyNaC2/3/5 in human HEK 293 and monkey COS-7 cells. We found that the expression of HyNaC2/3/5 using native cDNAs was inefficient and that codon optimization strongly increased protein expression and current amplitude in patch-clamp experiments. We used the improved expression of codon-optimized channel subunits to perform Ca²⁺ imaging and to demonstrate their glycosylation pattern. In summary, we established efficient expression of a cnidarian ion channel in mammalian cell lines.

Mechanisms of Attenuation by Genetic Recoding of Viruses

The development of safe and effective vaccines against viruses is central to disease control. With advancements in DNA synthesis technology, the production of synthetic viral genomes has fueled many research efforts that aim to generate attenuated viruses by introducing synonymous mutations. Elucidation of the mechanisms underlying virus attenuation through synonymous mutagenesis is revealing interesting new biology that can be exploited for vaccine development. Here, we review recent advancements in this field of synthetic virology and focus on the molecular mechanisms of attenuation by genetic recoding of viruses. We highlight the action of the zinc finger antiviral protein (ZAP) and RNase L, two proteins involved in the inhibition of viruses enriched for CpG and UpA dinucleotides, that are often the products of virus recoding algorithms. Additionally, we discuss current challenges in the field as well as studies that may illuminate how other host functions, such as translation, are potentially involved in the attenuation of recoded viruses.

Hijacking tRNAs From Translation: Regulatory Functions of tRNAs in Mammalian Cell Physiology

Transfer tRNAs (tRNAs) are small non-coding RNAs that are highly conserved in all kingdoms of life. Originally discovered as the molecules that deliver amino acids to the growing polypeptide chain during protein synthesis, tRNAs have been believed for a long time to play exclusive role in translation. However, recent studies have identified key roles for tRNAs and tRNA-derived small RNAs in multiple other processes, including regulation of transcription and translation, posttranslational modifications, stress response, and disease. These emerging roles suggest that tRNAs may be central players in the complex machinery of biological regulatory pathways. Here we overview these non-canonical roles of tRNA in normal physiology and disease, focusing largely on eukaryotic and mammalian systems.

Intra-genome variability in the dinucleotide composition of SARS-CoV-2

CpG dinucleotides are under-represented in the genomes of single-stranded RNA viruses, and SARS-CoV-2 is no exception to this. Artificial modification of CpG frequency is a valid approach for live attenuated vaccine development; if this is to be applied to SARS-CoV-2, we must first understand the role CpG motifs play in regulating SARS-CoV-2 replication. Accordingly, the CpG composition of the SARS-CoV-2 genome was characterised. CpG suppression among coronaviruses does not differ between virus genera but does vary with host species and primary replication site (a proxy for tissue tropism), supporting the hypothesis that viral CpG content may influence cross-species transmission. Although SARS-CoV-2 exhibits overall strong CpG suppression, this varies considerably across the genome, and the Envelope (E) open reading frame (ORF) and ORF10 demonstrate an absence of CpG suppression. Across the Coronaviridae, E genes display remarkably high variation in CpG composition, with those of SARS and SARS-CoV-2 having much higher CpG content than other coronaviruses isolated from humans. This is an ancestrally derived trait reflecting their bat origins. Conservation of CpG motifs in these regions suggests that they have a functionality which over-rides the need to suppress CpG; an observation relevant to future strategies towards a rationally attenuated SARS-CoV-2 vaccine.

Natural Selection Shapes Codon Usage in the Human Genome

Synonymous codon usage has been identified as a determinant of translational efficiency and mRNA stability in model organisms and human cell lines. However, whether natural selection shapes human codon content to optimize translation efficiency is unclear. Furthermore, aside from those that affect splicing, synonymous mutations are typically ignored as potential contributors to disease. Using genetic sequencing data from nearly 200,000 individuals, we uncover clear evidence that natural selection optimizes codon content in the human genome. In deriving intolerance metrics to quantify gene-level constraint on synonymous variation, we discover that dosage-sensitive genes, DNA-damage-response genes, and cell-cycle-regulated genes are particularly intolerant to synonymous variation. Notably, we illustrate that reductions in codon optimality in BRCA1 can attenuate its function. Our results reveal that synonymous mutations most likely play an underappreciated role in human variation.

Designing a novel mRNA vaccine against SARS-CoV-2: An immunoinformatics approach

SARS-CoV-2 is the deadly virus behind COVID-19, the disease that went on to ravage the world and caused the biggest pandemic 21st century has witnessed so far. On the face of ongoing death and destruction, the urgent need for the discovery of a vaccine against the virus is paramount. This study resorted to the emerging discipline of immunoinformatics in order to design a multi-epitope mRNA vaccine against the spike glycoprotein of SARS-CoV-2. Various immunoinformatics tools were utilized to predict T and B lymphocyte epitopes. The epitopes were channeled through a filtering pipeline comprised of antigenicity, toxicity, allergenicity, and cytokine inducibility evaluation with the goal of selecting epitopes capable of generating both T and B cell-mediated immune responses. Molecular docking simulation between the epitopes and their corresponding MHC molecules was carried out. 13 epitopes, a highly immunogenic adjuvant, elements for proper sub-cellular trafficking, a secretion booster, and appropriate linkers were combined for constructing the vaccine. The vaccine was found to be antigenic, almost neutral at physiological pH, non-toxic, non-allergenic, capable of generating a robust immune response and had a decent worldwide population coverage. Based on these parameters, this design can be considered a promising choice for a vaccine against SARS-CoV-2.

•

SARS-CoV-2, the causative agent for COVID-19, has caused the biggest pandemic 21st century has witnessed so far.

•

Here, immunoinformatics was used for designing a muti-epitope mRNA vaccine against the spike glycoprotein of SARS-CoV-2.

•

Through stringent selection of T and B cell epitopes and other necessary elements, a vaccine was constructed in silico.

•

Proposed mechanism of its synthesis, delivery and action has also been presented.

•

The vaccine was found to be immunogenic, almost neutral, non-toxic, non-allergenic and has a decent worldwide coverage.

A Lassa Fever Live-Attenuated Vaccine Based on Codon Deoptimization of the Viral Glycoprotein Gene

Lassa virus (LASV) infects several hundred thousand people in Western Africa, resulting in many lethal Lassa fever (LF) cases. Licensed LF vaccines are not available, and anti-LF therapy is limited to off-label use of the nucleoside analog ribavirin with uncertain efficacy. We describe the generation of a novel live-attenuated LASV vaccine candidate. This vaccine candidate is based on mutating wild-type (WT) LASV in a key region of the viral genome, the glycoprotein precursor (GPC) gene. These mutations do not change the encoded GPC but interfere with its production in host cells. This mutated LASV (rLASV-GPC/CD) behaves like WT LASV (rLASV-WT) in cell culture, but in contrast to rLASV-WT, does not cause disease in inoculated guinea pigs. Guinea pigs immunized with rLASV-GPC/CD were protected against an otherwise lethal exposure to WT LASV. Our results support the testing of this candidate vaccine in nonhuman primate models ofLF.

Lassa virus (LASV) is endemic in Western Africa and is estimated to infect hundreds of thousands of individuals annually. A considerable number of these infections result in Lassa fever (LF), which is associated with significant morbidity and a case-fatality rate as high as 69% among hospitalized confirmed patients. U.S. Food and Drug Administration-approved LF vaccines are not available. Current antiviral treatment is limited to off-label use of a nucleoside analogue, ribavirin, that is only partially effective and associated with significant side effects. We generated and characterized a recombinant LASV expressing a codon-deoptimized (CD) glycoprotein precursor gene (GPC), rLASV-GPC/CD. Comparison of growth kinetics and peak titers showed that rLASV-GPC/CD is slightly attenuated in cell culture compared to wild-type (WT) recombinant LASV (rLASV-WT). However, rLASV-GPC/CD is highly attenuated in strain 13 and Hartley guinea pigs, as reflected by the absence of detectable clinical signs in animals inoculated with rLASV-GPC/CD. Importantly, a single subcutaneous dose of rLASV-GPC/CD provides complete protection against an otherwise lethal exposure to LASV. Our results demonstrate the feasibility of implementing a CD approach for developing a safe and effective LASV live-attenuated vaccine candidate. Moreover, rLASV-GPC/CD might provide investigators with a tool to safely study LASV outside maximum (biosafety level 4) containment, which could accelerate the elucidation of basic aspects of the molecular and cell biology of LASV and the development of novel LASV medical countermeasures.

The origin of chromosomal histones in a 30S ribosomal protein

Histones are genes that regulate chromatin structure. They are present in both eukaryotes and archaea, and form nucleosomes with DNA, but their exact evolutionary origins have hitherto remained a mystery. A longstanding hypothesis is that they have precursors in ribosomal proteins with whom they share much in common in terms of size and chemistry. By examining the proteome of the Asgard archaeon, Lokiarchaeum, the most conserved of all the histones, H4, is found to plausibly be homologous with one of its 30S ribosomal proteins, RPS6. This is based on both sequence identity and statistical analysis. The N-terminal tail, containing key sites involved in post-translational modifications, is notably present in the precursor gene. Although other archaeal groups possess similar homologs, these are not as close to H4 as the one found in Lokiarchaeum. The other core histones, H2A, H2B and H3, appear to have also evolved from the same ribosomal protein. Parts of H4 are also similar to another ribosomal protein, namely RPS15, suggesting that the ancestral precursor could have resembled both. Eukaryotic histones, in addition, appear to have an independent origin to that of their archaeal counterparts that evolved from similar, but still different, 30S subunit proteins, some of which are much more like histones in terms of their physical structure. The nucleosome may, therefore, be not only of archaeal but also of ribosomal origin.

Genetic evolution and codon usage analysis of NKX-2.5 gene governing heart development in some mammals

NKX-2.5 gene is responsible for cardiac development and its targeted disruption apprehends cardiac development at the linear heart tube stage. Bioinformatic analysis was employed to investigate the codon usage pattern and dN/dS of mammalian NKX-2.5 gene. The relative synonymous codon usage analysis revealed variation in codon usage and two synonymous codons namely ATA (Ile) and GTA (Val) were absent in NKX-2.5 gene across selected mammalian species suggesting that these two codons were possibly selected against during evolution. Parity rule 2 analysis of two and four fold amino acids showed CT bias whereas six-fold amino acids revealed GA bias. Neutrality analysis suggests that selection played a prominent role while mutation had a minor role. The dN/dS analysis suggests synonymous substitution played a significant role and it negatively correlated with p-distance of the gene. Purifying natural selection played a dominant role in the genetic evolution of NKX-2.5 gene in mammals.

Codon Usage Differences among Genes Expressed in Different Tissues of Drosophila melanogaster

Codon usage patterns are affected by both mutational biases and translational selection. The frequency at which each codon is used in the genome is directly linked to the cellular concentrations of their corresponding tRNAs. Transfer RNA abundances—as well as the abundances of other potentially relevant factors, such as RNA-binding proteins—may vary across different tissues, making it possible that genes expressed in different tissues are subject to different translational selection regimes, and thus differ in their patterns of codon usage. These differences, however, are poorly understood, having been studied only in Arabidopsis, rice and human, with controversial results in human. Drosophila melanogaster is a suitable model organism to study tissue-specific codon adaptation given its large effective population size. Here, we compare 2,046 genes, each expressed specifically in one tissue of D. melanogaster. We show that genes expressed in different tissues exhibit significant differences in their patterns of codon usage, and that these differences are only partially due to differences in GC content, expression levels, or protein lengths. Remarkably, these differences are stronger when analyses are restricted to highly expressed genes. Our results strongly suggest that genes expressed in different tissues are subject to different regimes of translational selection.

Characterization of bovine (Bos taurus) imprinted genes from genomic to amino acid attributes by data mining approaches

Genomic imprinting results in monoallelic expression of genes in mammals and flowering plants. Understanding the function of imprinted genes improves our knowledge of the regulatory processes in the genome. In this study, we have employed classification and clustering algorithms with attribute weighting to specify the unique attributes of both imprinted (monoallelic) and biallelic expressed genes. We have obtained characteristics of 22 known monoallelically expressed (imprinted) and 8 biallelic expressed genes that have been experimentally validated alongside 208 randomly selected genes in bovine (Bos taurus). Attribute weighting methods and various supervised and unsupervised algorithms in machine learning were applied. Unique characteristics were discovered and used to distinguish mono and biallelic expressed genes from each other in bovine. To obtain the accuracy of classification, 10-fold cross-validation with concerning each combination of attribute weighting (feature selection) and machine learning algorithms, was used. Our approach was able to accurately predict mono and biallelic genes using the genomics and proteomics attributes.

Widespread selection for extremely high and low levels of secondary structure in coding sequences across all domains of life

Codon composition, GC content and local RNA secondary structures can have a profound effect on gene expression, and mutations affecting these parameters, even though they do not alter the protein sequence, are not neutral in terms of selection. Although evidence exists that, in some cases, selection favours more stable RNA secondary structures, we currently lack a concrete idea of how many genes are affected within a species, and whether this is a universal phenomenon in nature. We searched for signs of structural selection in a global manner, analysing a set of 1 million coding sequences from 73 species representing all domains of life, as well as viruses, by means of our newly developed software PACKEIS. We show that codon composition and amino acid identity are main determinants of RNA secondary structure. In addition, we show that the arrangement of synonymous codons within coding sequences is non-random, yielding extremely high, but also extremely low, RNA structuredness significantly more often than expected by chance. Taken together, we demonstrate that selection for high and low levels of secondary structure is a widespread phenomenon. Our results provide another line of evidence that synonymous mutations are less neutral than commonly thought, which is of importance for many evolutionary models.

Nonsynonymous, synonymous and nonsense mutations in human cancer-related genes undergo stronger purifying selections than expectation

Background

Nonsynonymous mutations change the protein sequences and are frequently subjected to natural selection. The same goes for nonsense mutations that introduce pre-mature stop codons into CDSs (coding sequences). Synonymous mutations, however, are intuitively thought to be functionally silent and evolutionarily neutral. Now researchers know that the optimized synonymous codon usage is advantageous in the speedy mRNA translation process. With the advent of NGS technique, the explosion of NGS data generated from the tumor tissues help researchers identify driver mutations in cancer-related genes, but relatively less attention is paid to the SNP data in healthy human populations when studying cancer.

Methods

Here, we analyzed the publically available human SNPs. We classified these SNPs according to their functional and evolutionary categories. By simply dividing the human genes into cancer-related genes and other genes, we compared the features of nonsynonymous, synonymous and nonsense mutations in these two gene sets from multiple aspects.

Results

We provided lines of evidence that the nonsynonymous, synonymous and nonsense mutations in cancer-related genes undergo stronger purifying selection when compared to the expected pattern in other genes. The lower nonsynonymous to synonymous ratio observed in cancer-related genes suggests the suppression of amino acid substitutions in these genes. The synonymous SNPs, after excluding those in splicing regions, exhibit preferred changes in codon usage and higher codon frequencies in cancer-related genes compared to other genes, indicating the constraint exerted on these mutations. Nonsense mutations are less frequent and located closer to stop codons in cancer-related genes than in other genes, which putatively minimize their deleterious effects.

Conclusion

Our study demonstrated the evolutionary constraint on mutations in CDS of cancer-related genes without the requirement of data from cancer tissues or patients. Our work provides novel perspectives on interpreting the constraint on mutations in cancer-related genes. We reveal extra constraint on synonymous mutations in cancer-related genes which is related to codon usage bias and is in addition to the splicing effect.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5572-x) contains supplementary material, which is available to authorized users.

mRNA as novel technology for passive immunotherapy

While active immunization elicits a lasting immune response by the body, passive immunotherapy transiently equips the body with exogenously generated immunological effectors in the form of either target-specific antibodies or lymphocytes functionalized with target-specific receptors. In either case, administration or expression of recombinant proteins plays a fundamental role. mRNA prepared by in vitro transcription (IVT) is increasingly appreciated as a drug substance for delivery of recombinant proteins. With its biological role as transient carrier of genetic information translated into protein in the cytoplasm, therapeutic application of mRNA combines several advantages. For example, compared to transfected DNA, mRNA harbors inherent safety features. It is not associated with the risk of inducing genomic changes and potential adverse effects are only temporary due to its transient nature. Compared to the administration of recombinant proteins produced in bioreactors, mRNA allows supplying proteins that are difficult to manufacture and offers extended pharmacokinetics for short-lived proteins. Based on great progress in understanding and manipulating mRNA properties, efficacy data in various models have now demonstrated that IVT mRNA constitutes a potent and flexible platform technology. Starting with an introduction into passive immunotherapy, this review summarizes the current status of IVT mRNA technology and its application to such immunological interventions.

Accounting for Programmed Ribosomal Frameshifting in the Computation of Codon Usage Bias Indices

Experimental evidence shows that synonymous mutations can have important consequences on genetic fitness. Many organisms display codon usage bias (CUB), where synonymous codons that are translated into the same amino acid appear with distinct frequency. Within genomes, CUB is thought to arise from selection for translational efficiency and accuracy, termed the translational efficiency hypothesis (TEH). Indeed, CUB indices correlate with protein expression levels, which is widely interpreted as evidence for translational selection. However, these tests neglect -1 programmed ribosomal frameshifting (-1 PRF), an important translational disruption effect found across all organisms of the tree of life. Genes that contain -1 PRF signals should cost more to express than genes without. Thus, CUB indices that do not consider -1 PRF may overestimate genes’ true adaptation to translational efficiency and accuracy constraints. Here, we first investigate whether -1 PRF signals do indeed carry such translational cost. We then propose two corrections for CUB indices for genes containing -1 PRF signals. We retest the TEH in Saccharomyces cerevisiae under these corrections. We find that the correlation between corrected CUB index and protein expression remains intact for most levels of uniform -1 PRF efficiencies, and tends to increase when these efficiencies decline with protein expression. We conclude that the TEH is strengthened and that -1 PRF events constitute a promising and useful tool to examine the relationships between CUB and selection for translation efficiency and accuracy.

Cells alter their tRNA abundance to selectively regulate protein synthesis during stress conditions

Decoding the information in mRNA during protein synthesis relies on tRNA adaptors, the abundance of which can affect the decoding rate and translation efficiency. To determine whether cells alter tRNA abundance to selectively regulate protein expression, we quantified changes in the abundance of individual tRNAs at different time points in response to diverse stress conditions in Saccharomyces cerevisiae. We found that the tRNA pool was dynamic and rearranged in a manner that facilitated selective translation of stress-related transcripts. Through genomic analysis of multiple data sets, stochastic simulations, and experiments with designed sequences of proteins with identical amino acids but altered codon usage, we showed that changes in tRNA abundance affected protein expression independently of factors such as mRNA abundance. We suggest that cells alter their tRNA abundance to selectively affect the translation rates of specific transcripts to increase the amounts of required proteins under diverse stress conditions.

The adaptation of codon usage of +ssRNA viruses to their hosts

Viruses depend on their host's cellular structure to survive. Most of them do not have tRNAs, their translation relies on hosts' tRNA pools. Over the course of evolution, viruses needed to optimally exploit cellular processes of their host. Thus, codon usage of a virus should coevolve with its host to efficiently and rapidly replicate. Some viruses can invade a broad spectrum of hosts (BSTVs), while others can invade a narrow spectrum only (NSTVs). Consequently, we test the hypothesis that similarity of codon usage preference and the degree of matching between BSTVs and their hosts will be lower than that of NSTVs, which only need to coevolve with few hosts. We compare the patterns of codon usage in 255 virus genomes to test this hypothesis. Our results show that NSTVs have a higher degree of matching to their hosts' tRNA pools than BSTVs. Further, analysis of the effective number of codons (ENC) infers that codon usage bias of NSTVs is relatively stronger than that of BSTVs. Thus, codon usage of NSTVs tends to better match their host than that of BSTVs. This supports the hypothesis that viruses adapt to the expression system of their host(s).

Biological meaning of "habitable zone" in nucleotide composition space

Organisms generally display two contrasting properties: large biodiversity and a uniform state of "life". In this study, we focused on the question of how genome sequences describe "life" where a large number of biomolecules are harmonized. We analyzed the whole genome sequence of 2664 organisms, paying attention to the nucleotide composition which is an intensive parameter from the genome sequence. The results showed that all organisms were plotted in narrow regions of the nucleotide composition space of the first and second letters of the codon. Since all genome sequences overlap irrespective of the living environment, it can be called a "habitable zone". The habitable zone deviates by 500 times the standard deviation from the nucleotide composition expected from the random sequence, indicating that unexpectedly rare sequences are realized. Furthermore, we found that the habitable zones at the first and second letters of the codon serve as the background mechanisms for the functional network of biological systems. The habitable zone at the second letter of the codon controls the formation of transmembrane regions and the habitable zone at the first letter controls the formation of molecular recognition unit. These analyses showed that the habitable zone of the nucleotide composition space and the exquisite arrangement of amino acids in the codon table are conjugated to form biological systems. Finally, we discussed the evolution of the higher order of genome sequences.

Efficient mAb production in CHO cells with optimized signal peptide, codon, and UTR

Antibody drugs have been used to treat a number of diseases successfully. Producing antibodies with high yield and quality is necessary for clinical applications of antibodies. For a candidate molecule, optimization of a vector to produce sufficient yield and an accurate primary structure is indispensable in the early stage of the production process development. It is especially important to maintain the fidelity of N-terminal sequence. In order to produce antibodies with a high yield and accurate N-terminal, the expression vector was systematically optimized in this study. First, the heavy chain and light chain were co-expressed in Chinese hamster ovary (CHO) cells with different signal peptides. Mass spectrometry (MS) revealed that signal peptides Esp-K, Bsp-H, and 8Hsp-H were accurately deleted from mature antibodies. Further, the yield was doubled by codon optimization and increased by 50% with the presence of untranslated regions (UTR). The combination of UTR with optimal codon and signal peptide to form an expression vector resulted in yield improvement of 150% and correct N-terminal sequences. Moreover, the main product peak was above 98% as assessed by size-exclusion chromatography (SEC). Additionally, the bioactivity of products made from optimized transient gene expression (TGE) was almost identical to the standard sample. The production efficiency and product quality from the identified TGE optimization strategy was further demonstrated through application to two other antibodies. The expression level of SGE (stable gene expression) can also be improved effectively with this optimization strategy. In conclusion, vector optimization via combination of optimized signal peptide, codon, and UTR is an alternative approach for efficient antibody production with high fidelity N-terminal sequence in CHO cells.

Essentiality, conservation, evolutionary pressure and codon bias in bacterial genomes

Essential genes constitute the core of genes which cannot be mutated too much nor lost along the evolutionary history of a species. Natural selection is expected to be stricter on essential genes and on conserved (highly shared) genes, than on genes that are either nonessential or peculiar to a single or a few species. In order to further assess this expectation, we study here how essentiality of a gene is connected with its degree of conservation among several unrelated bacterial species, each one characterised by its own codon usage bias. Confirming previous results on E. coli, we show the existence of a universal exponential relation between gene essentiality and conservation in bacteria. Moreover, we show that, within each bacterial genome, there are at least two groups of functionally distinct genes, characterised by different levels of conservation and codon bias: i) a core of essential genes, mainly related to cellular information processing; ii) a set of less conserved nonessential genes with prevalent functions related to metabolism. In particular, the genes in the first group are more retained among species, are subject to a stronger purifying conservative selection and display a more limited repertoire of synonymous codons. The core of essential genes is close to the minimal bacterial genome, which is in the focus of recent studies in synthetic biology, though we confirm that orthologs of genes that are essential in one species are not necessarily essential in other species. We also list a set of highly shared genes which, reasonably, could constitute a reservoir of targets for new anti-microbial drugs.

Codon-Resolution Analysis Reveals a Direct and Context-Dependent Impact of Individual Synonymous Mutations on mRNA Level

Codon usage bias (CUB) refers to the observation that synonymous codons are not used equally frequently in a genome. CUB is stronger in more highly expressed genes, a phenomenon commonly explained by stronger natural selection on translational accuracy and/or efficiency among these genes. Nevertheless, this phenomenon could also occur if CUB regulates gene expression at the mRNA level, a hypothesis that has not been tested until recently. Here, we attempt to quantify the impact of synonymous mutations on mRNA level in yeast using 3,556 synonymous variants of a heterologous gene encoding green fluorescent protein (GFP) and 523 synonymous variants of an endogenous gene TDH3. We found that mRNA level was positively correlated with CUB among these synonymous variants, demonstrating a direct role of CUB in regulating transcript concentration, likely via regulating mRNA degradation rate, as our additional experiments suggested. More importantly, we quantified the effects of individual synonymous mutations on mRNA level and found them dependent on 1) CUB and 2) mRNA secondary structure, both in proximal sequence contexts. Our study reveals the pleiotropic effects of synonymous codon usage and provides an additional explanation for the well-known correlation between CUB and gene expression level.

Dissimilar substitution rates between two strands of DNA influence codon usage pattern in some human genes

We illustrated the descriptive aspects of codon usage of some important human genes and their expression potential in E. coli. By comparing the results of various codon usage parameters, effects that are due to selection and mutational pressures have been deciphered. The variation in GC3s explains a significant proportion of the variation in codon usage patterns. The codons CGC, CGG, CTG and GCG showed strong positive correlation with GC3, which suggested that codon usage had been influenced by GC bias. We also found that ACC (Thr, RSCU-1.77), GCC (Ala, RSCU-1.67), CCC (Pro, RSCU-1.54), TCC (Ser, RSCU-1.47) were frequently used which signified that C was common at 2nd and 3rd codon positions. Correspondence analysis revealed that F1 axis had significant correlation with various GC contents suggesting that compositional properties under mutation pressure might affect codon usage bias. Nc-GC3 plot analysis suggested that both mutation pressure and natural selection might affect the codon usage bias which is also supported by neutrality plot analysis. The dinucleotide CT, TG and AG were significantly over-represented and CG, TA, AT, TT, and GT were underrepresented due to high rate of spontaneous mutation resulting from cytosine deamination.

Synonymous Codon Usage Controls Various Molecular Aspects

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

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.

Recombination, meiotic expression and human codon usage

Synonymous codon usage (SCU) varies widely among human genes. In particular, genes involved in different functional categories display a distinct codon usage, which was interpreted as evidence that SCU is adaptively constrained to optimize translation efficiency in distinct cellular states. We demonstrate here that SCU is not driven by constraints on tRNA abundance, but by large-scale variation in GC-content, caused by meiotic recombination, via the non-adaptive process of GC-biased gene conversion (gBGC). Expression in meiotic cells is associated with a strong decrease in recombination within genes. Differences in SCU among functional categories reflect differences in levels of meiotic transcription, which is linked to variation in recombination and therefore in gBGC. Overall, the gBGC model explains 70% of the variance in SCU among genes. We argue that the strong heterogeneity of SCU induced by gBGC in mammalian genomes precludes any optimization of the tRNA pool to the demand in codon usage.

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

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

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