Codon Usage Influences the Local Rate of Translation Elongation to Regulate Co-translational Protein Folding

Unexpected myriad of co-occurring viral strains and species in one of the most abundant and microdiverse viruses on Earth

Viral genetic microdiversity drives adaptation, pathogenicity, and speciation and has critical consequences for the viral-host arms race occurring at the strain and species levels, which ultimately impact microbial community structure and biogeochemical cycles. Despite the fact that most efforts have focused on viral macrodiversity, little is known about the microdiversity of ecologically important viruses on Earth. Recently, single-virus genomics discovered the putatively most abundant ocean virus in temperate and tropical waters: the uncultured dsDNA virus vSAG 37-F6 infecting Pelagibacter, the most abundant marine bacteria. In this study, we report the cooccurrence of up to ≈1,500 different viral strains (>95% nucleotide identity) and ≈30 related species (80-95% nucleotide identity) in a single oceanic sample. Viral microdiversity was maintained over space and time, and most alleles were the result of synonymous mutations without any apparent adaptive benefits to cope with host translation codon bias and efficiency. Gene flow analysis used to delimitate species according to the biological species concept (BSC) revealed the impact of recombination in shaping vSAG 37-F6 virus and Pelagibacter speciation. Data demonstrated that this large viral microdiversity somehow mirrors the host species diversity since ≈50% of the 926 analyzed Pelagibacter genomes were found to belong to independent BSC species that do not significantly engage in gene flow with one another. The host range of this evolutionarily successful virus revealed that a single viral species can infect multiple Pelagibacter BSC species, indicating that this virus crosses not only formal BSC barriers but also biomes since viral ancestors are found in freshwater.

Ribosome dynamics and mRNA turnover, a complex relationship under constant cellular scrutiny

Eukaryotic gene expression is closely regulated by translation and turnover of mRNAs. Recent advances highlight the importance of translation in the control of mRNA degradation, both for aberrant and apparently normal mRNAs. During translation, the information contained in mRNAs is decoded by ribosomes, one codon at a time, and tRNAs, by specifically recognizing codons, translate the nucleotide code into amino acids. Such a decoding step does not process regularly, with various obstacles that can hinder ribosome progression, then leading to ribosome stalling or collisions. The progression of ribosomes is constantly monitored by the cell which has evolved several translation-dependent mRNA surveillance pathways, including nonsense-mediated decay (NMD), no-go decay (NGD), and non-stop decay (NSD), to degrade certain problematic mRNAs and the incomplete protein products. Recent progress in sequencing and ribosome profiling has made it possible to discover new mechanisms controlling ribosome dynamics, with numerous crosstalks between translation and mRNA decay. We discuss here various translation features critical for mRNA decay, with particular focus on current insights from the complexity of the genetic code and also the emerging role for the ribosome as a regulatory hub orchestrating mRNA decay, quality control, and stress signaling. Even if the interplay between mRNA translation and degradation is no longer to be demonstrated, a better understanding of their precise coordination is worthy of further investigation. This article is categorized under: RNA Turnover and Surveillance > Regulation of RNA Stability Translation > Translation Regulation RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Codon optimality in cancer

A key characteristic of cancer cells is their increased proliferative capacity, which requires elevated levels of protein synthesis. The process of protein synthesis involves the translation of codons within the mRNA coding sequence into a string of amino acids to form a polypeptide chain. As most amino acids are encoded by multiple codons, the nucleotide sequence of a coding region can vary dramatically without altering the polypeptide sequence of the encoded protein. Although mutations that do not alter the final amino acid sequence are often thought of as silent/synonymous, these can still have dramatic effects on protein output. Because each codon has a distinct translation elongation rate and can differentially impact mRNA stability, each codon has a different degree of 'optimality' for protein synthesis. Recent data demonstrates that the codon preference of a transcriptome matches the abundance of tRNAs within the cell and that this supply and demand between tRNAs and mRNAs varies between different cell types. The largest observed distinction is between mRNAs encoding proteins associated with proliferation or differentiation. Nevertheless, precisely how codon optimality and tRNA expression levels regulate cell fate decisions and their role in malignancy is not fully understood. This review describes the current mechanistic understanding on codon optimality, its role in malignancy and discusses the potential to target codon optimality therapeutically in the context of cancer.

Comparative analysis of transcriptomic data shows the effects of multiple evolutionary selection processes on codon usage in Marsupenaeus japonicus and Marsupenaeus pulchricaudatus

BACKGROUND

Kuruma shrimp, a major commercial shrimp species in the world, has two cryptic or sibling species, Marsupenaeus japonicus and Marsupenaeus pulchricaudatus. Codon usage analysis would contribute to our understanding of the genetic and evolutionary characteristics of the two Marsupenaeus species. In this study, we analyzed codon usage and related indices using coding sequences (CDSs) from RNA-seq data.

RESULTS

Using CodonW 1.4.2 software, we performed the codon bias analysis of transcriptomes obtained from hepatopancreas tissues, which indicated weak codon bias. Almost all parameters had similar correlations for both species. The gene expression level (FPKM) was negatively correlated with A/T3s. We determined 12 and 14 optimal codons for M. japonicus and M. pulchricaudatus, respectively, and all optimal codons have a C/G-ending. The two Marsupenaeus species had different usage frequencies of codon pairs, which contributed to further analysis of transcriptional differences between them. Orthologous genes that underwent positive selection (ω > 1) had a higher correlation coefficient than that of experienced purifying selection (ω < 1). Parity Rule 2 (PR2) and effective number of codons (ENc) plot analysis showed that the codon usage patterns of both species were influenced by both mutations and selection. Moreover, the average observed ENc value was lower than the expected value for both species, suggesting that factors other than GC may play roles in these phenomena. The results of multispecies clustering based on codon preference were consistent with traditional classification.

CONCLUSIONS

This study provides a relatively comprehensive understanding of the correlations among codon usage bias, gene expression, and selection pressures of CDSs for M. japonicus and M. pulchricaudatus. The genetic evolution was driven by mutations and selection pressure. Moreover, the results point out new insights into the specificities and evolutionary characteristics of the two Marsupenaeus species.

Heterologous production of cellulose- and starch-degrading hydrolases to expand Saccharomyces cerevisiae substrate utilization: Lessons learnt

Selected strains of Saccharomyces cerevisiae are used for commercial bioethanol production from cellulose and starch, but the high cost of exogenous enzymes for substrate hydrolysis remains a challenge. This can be addressed through consolidated bioprocessing (CBP) where S. cerevisiae strains are engineered to express recombinant glycoside hydrolases during fermentation. Looking back at numerous strategies undertaken over the past four decades to improve recombinant protein production in S. cerevisiae, it is evident that various steps in the protein production "pipeline" can be manipulated depending on the protein of interest and its anticipated application. In this review, we briefly introduce some of the strategies and highlight lessons learned with regards to improved transcription, translation, post-translational modification and protein secretion of heterologous hydrolases. We examine how host strain selection and modification, as well as enzyme compatibility, are crucial determinants for overall success. Finally, we discuss how lessons from heterologous hydrolase expression can inform modern synthetic biology and genome editing tools to provide process-ready yeast strains in future. However, it is clear that the successful expression of any particular enzyme is still unpredictable and requires a trial-and-error approach.

A dual-reporter system for investigating and optimizing protein translation and folding in E. coli

Strategies for investigating and optimizing the expression and folding of proteins for biotechnological and pharmaceutical purposes are in high demand. Here, we describe a dual-reporter biosensor system that simultaneously assesses in vivo protein translation and protein folding, thereby enabling rapid screening of mutant libraries. We have validated the dual-reporter system on five different proteins and find an excellent correlation between reporter signals and the levels of protein expression and solubility of the proteins. We further demonstrate the applicability of the dual-reporter system as a screening assay for deep mutational scanning experiments. The system enables high throughput selection of protein variants with high expression levels and altered protein stability. Next generation sequencing analysis of the resulting libraries of protein variants show a good correlation between computationally predicted and experimentally determined protein stabilities. We furthermore show that the mutational experimental data obtained using this system may be useful for protein structure calculations.

Genetic and process engineering strategies for enhanced recombinant N-glycoprotein production in bacteria

Background

The production of N-linked glycoproteins in genetically amenable bacterial hosts offers great potential for reduced cost, faster/simpler bioprocesses, greater customisation, and utility for distributed manufacturing of glycoconjugate vaccines and glycoprotein therapeutics. Efforts to optimize production hosts have included heterologous expression of glycosylation enzymes, metabolic engineering, use of alternative secretion pathways, and attenuation of gene expression. However, a major bottleneck to enhance glycosylation efficiency, which limits the utility of the other improvements, is the impact of target protein sequon accessibility during glycosylation.

Results

Here, we explore a series of genetic and process engineering strategies to increase recombinant N-linked glycosylation, mediated by the Campylobacter-derived PglB oligosaccharyltransferase in Escherichia coli. Strategies include increasing membrane residency time of the target protein by modifying the cleavage site of its secretion signal, and modulating protein folding in the periplasm by use of oxygen limitation or strains with compromised oxidoreductase or disulphide-bond isomerase activity. These approaches achieve up to twofold improvement in glycosylation efficiency. Furthermore, we also demonstrate that supplementation with the chemical oxidant cystine enhances the titre of glycoprotein in an oxidoreductase knockout strain by improving total protein production and cell fitness, while at the same time maintaining higher levels of glycosylation efficiency.

Conclusions

In this study, we demonstrate that improved protein glycosylation in the heterologous host could be achieved by mimicking the coordination between protein translocation, folding and glycosylation observed in native host such as Campylobacter jejuni and mammalian cells. Furthermore, it provides insight into strain engineering and bioprocess strategies, to improve glycoprotein yield and titre, and to avoid physiological burden of unfolded protein stress upon cell growth. The process and genetic strategies identified herein will inform further optimisation and scale-up of heterologous recombinant N-glycoprotein production.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01689-x.

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.

Differential regulation of mRNA fate by the human Ccr4-Not complex is driven by coding sequence composition and mRNA localization

BACKGROUND

Regulation of protein output at the level of translation allows for a rapid adaptation to dynamic changes to the cell's requirements. This precise control of gene expression is achieved by complex and interlinked biochemical processes that modulate both the protein synthesis rate and stability of each individual mRNA. A major factor coordinating this regulation is the Ccr4-Not complex. Despite playing a role in most stages of the mRNA life cycle, no attempt has been made to take a global integrated view of how the Ccr4-Not complex affects gene expression.

RESULTS

This study has taken a comprehensive approach to investigate post-transcriptional regulation mediated by the Ccr4-Not complex assessing steady-state mRNA levels, ribosome position, mRNA stability, and protein production transcriptome-wide. Depletion of the scaffold protein CNOT1 results in a global upregulation of mRNA stability and the preferential stabilization of mRNAs enriched for G/C-ending codons. We also uncover that mRNAs targeted to the ER for their translation have reduced translational efficiency when CNOT1 is depleted, specifically downstream of the signal sequence cleavage site. In contrast, translationally upregulated mRNAs are normally localized in p-bodies, contain disorder-promoting amino acids, and encode nuclear localized proteins. Finally, we identify ribosome pause sites that are resolved or induced by the depletion of CNOT1.

CONCLUSIONS

We define the key mRNA features that determine how the human Ccr4-Not complex differentially regulates mRNA fate and protein synthesis through a mechanism linked to codon composition, amino acid usage, and mRNA localization.

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

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

A Code Within a Code: How Codons Fine-Tune Protein Folding in the Cell

The genetic code sets the correspondence between the sequence of a given nucleotide triplet in an mRNA molecule, called a codon, and the amino acid that is added to the growing polypeptide chain during protein synthesis. With four bases (A, G, U, and C), there are 64 possible triplet codons: 61 sense codons (encoding amino acids) and 3 nonsense codons (so-called, stop codons that define termination of translation). In most organisms, there are 20 common/standard amino acids used in protein synthesis; thus, the genetic code is redundant with most amino acids (with the exception of Met and Trp) are being encoded by more than one (synonymous) codon. Synonymous codons were initially presumed to have entirely equivalent functions, however, the finding that synonymous codons are not present at equal frequencies in mRNA suggested that the specific codon choice might have functional implications beyond coding for amino acid. Observation of nonequivalent use of codons in mRNAs implied a possibility of the existence of auxiliary information in the genetic code. Indeed, it has been found that genetic code contains several layers of such additional information and that synonymous codons are strategically placed within mRNAs to ensure a particular translation kinetics facilitating and fine-tuning co-translational protein folding in the cell via step-wise/sequential structuring of distinct regions of the polypeptide chain emerging from the ribosome at different points in time. This review summarizes key findings in the field that have identified the role of synonymous codons and their usage in protein folding in the cell.

Genome-Wide Association and Selective Sweep Studies Reveal the Complex Genetic Architecture of DMI Fungicide Resistance in Cercospora beticola

The rapid and widespread evolution of fungicide resistance remains a challenge for crop disease management. The demethylation inhibitor (DMI) class of fungicides is a widely used chemistry for managing disease, but there has been a gradual decline in efficacy in many crop pathosystems. Reliance on DMI fungicides has increased resistance in populations of the plant pathogenic fungus Cercospora beticola worldwide. To better understand the genetic and evolutionary basis for DMI resistance in C. beticola, a genome-wide association study (GWAS) and selective sweep analysis were conducted for the first time in this species. We performed whole-genome resequencing of 190 C. beticola isolates infecting sugar beet (Beta vulgaris ssp. vulgaris). All isolates were phenotyped for sensitivity to the DMI tetraconazole. Intragenic markers on chromosomes 1, 4, and 9 were significantly associated with DMI fungicide resistance, including a polyketide synthase gene and the gene encoding the DMI target CbCYP51. Haplotype analysis of CbCYP51 identified a synonymous mutation (E170) and nonsynonymous mutations (L144F, I387M, and Y464S) associated with DMI resistance. Genome-wide scans of selection showed that several of the GWAS mutations for fungicide resistance resided in regions that have recently undergone a selective sweep. Using radial plate growth on selected media as a fitness proxy, we did not find a trade-off associated with DMI fungicide resistance. Taken together, we show that population genomic data from a crop pathogen can allow the identification of mutations conferring fungicide resistance and inform about their origins in the pathogen population.

The ins and outs of Bacillus proteases: activities, functions and commercial significance

Because the majority of bacterial species divide by binary fission, and do not have distinguishable somatic and germline cells, they could be considered to be immortal. However, bacteria ‘age’ due to damage to vital cell components such as DNA and proteins. DNA damage can often be repaired using efficient DNA repair mechanisms. However, many proteins have a functional ‘shelf life’; some are short lived, while others are relatively stable. Specific degradation processes are built into the life span of proteins whose activities are required to fulfil a specific function during a prescribed period of time (e.g. cell cycle, differentiation process, stress response). In addition, proteins that are irreparably damaged or that have come to the end of their functional life span need to be removed by quality control proteases. Other proteases are involved in performing a variety of specific functions that can be broadly divided into three categories: processing, regulation and feeding. This review presents a systematic account of the proteases of Bacillus subtilis and their activities. It reviews the proteases found in, or associated with, the cytoplasm, the cell membrane, the cell wall and the external milieu. Where known, the impacts of the deletion of particular proteases are discussed, particularly in relation to industrial applications.

Chloroplast genome features of Moricandia arvensis (Brassicaceae), a C3-C4 intermediate photosynthetic species

Moricandia arvensis, a plant species originating from the Mediterranean, has been classified as a rare C3-C4 intermediate species, and it is a possible bridge during the evolutionary process from C3 to C4 plant photosynthesis in the family Brassicaceae. Understanding the genomic structure, gene order, and gene content of chloroplasts (cp) of such species can provide a glimpse into the evolution of photosynthesis. In the present study, we obtained a well-annotated cp genome of M. arvensis using long PacBio and short Illumina reads with a de novo assembly strategy. The M. arvensis cp genome was a quadripartite circular molecule with the length of 153,312 bp, including two inverted repeats (IR) regions of 26,196 bp, divided by a small single copy (SSC) region of 17,786 bp and a large single copy (LSC) region of 83,134 bp. We detected 112 unigenes in this genome, comprising 79 protein-coding genes, 29 tRNAs, and four rRNAs. Forty-nine long repeat sequences and 51 simple sequence repeat (SSR) loci of 15 repeat types were identified. The analysis of Ks (synonymous) and Ka (non-synonymous) substitution rates indicated that the genes associated with “subunits of ATP synthase” (atpB), “subunits of NADH-dehydrogenase” (ndhG and ndhE), and “self-replication” (rps12 and rpl16) showed relatively higher Ka/Ks values than those of the other genes. The gene content, gene order, and LSC/IR/SSC boundaries and adjacent genes of the M. arvensis cp genome were highly conserved compared to those in related C3 species. Our phylogenetic analysis demonstrated that M. arvensis was clustered into a subclade with cultivated Brassica species and Raphanus sativus, indicating that M. arvensis was not involved in an independent evolutionary origin event. These results will open the way for further studies on the evolutionary process from C3 to C4 photosynthesis and hopefully provide guidance for utilizing M. arvensis as a resource for improvinng photosynthesis efficiency in cultivated Brassica species.

Luciferase-based reporter system for in vitro evaluation of elongation rate and processivity of ribosomes

The elongation step of translation is a key contributor to the abundance, folding and quality of proteins and to the stability of mRNA. However, control over translation elongation has not been thoroughly investigated. In this study, a Renilla-firefly luciferase fusion reporter system was further developed to investigate the in vitro elongation rate and processivity of ribosomes independent of the initiation and termination steps. The reporter mRNA was constructed to contain a single ORF encoding in-frame Renilla luciferase, a specific domain moiety and firefly luciferase. Such a reporter structure enables the quantitative and individual evaluation of the synthesis of a specific domain. As a proof of principle, the synthesis of three protein domains of different lengths and structures was analyzed. Using a cell-free translation assay, both the elongation rate and processivity of ribosomes were shown to vary depending on the domain synthesized. Additionally, a stalling sequence consisting of ten rare arginine codons notably reduced the elongation rate and the processivity of the ribosomes. All these results are consistent with the previously known dynamics of elongation in vivo. Overall, the methodology presented in this report provides a framework for studying aspects that contribute to the elongation step of translation.

Genome-wide role of codon usage on transcription and identification of potential regulators

Codon usage bias is a fundamental feature of all genomes and plays an important role in determining gene expression levels. The codon usage was thought to influence gene expression mainly due to its impact on translation. Recently, however, codon usage was shown to affect transcription of fungal and mammalian genes, indicating the existence of a gene regulatory phenomenon with unknown mechanism. In Neurospora, codon usage biases strongly correlate with mRNA levels genome-wide, and here we show that the correlation between codon usage and RNA levels is maintained in the nucleus. In addition, codon optimality is tightly correlated with both total and nuclear RNA levels, suggesting that codon usage broadly influences mRNA levels through transcription in a translation-independent manner. A large-scale RNA sequencing-based genetic screen in Neurospora identified 18 candidate factors that when deleted decreased the genome-wide correlation between codon usage and RNA levels and reduced the codon usage effect on gene expression. Most of these factors, such as the H3K36 methyltransferase, are chromatin regulators or transcription factors. Together, our results suggest that the transcriptional effect of codon usage is mediated by multiple transcriptional regulatory mechanisms.

Allosteric Activation of SARS-CoV-2 RNA-Dependent RNA Polymerase by Remdesivir Triphosphate and Other Phosphorylated Nucleotides

The catalytic subunit of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA-dependent RNA polymerase (RdRp) Nsp12 has a unique nidovirus RdRp-associated nucleotidyltransferase (NiRAN) domain that transfers nucleoside monophosphates to the Nsp9 protein and the nascent RNA. The NiRAN and RdRp modules form a dynamic interface distant from their catalytic sites, and both activities are essential for viral replication. We report that codon-optimized (for the pause-free translation in bacterial cells) Nsp12 exists in an inactive state in which NiRAN-RdRp interactions are broken, whereas translation by slow ribosomes and incubation with accessory Nsp7/8 subunits or nucleoside triphosphates (NTPs) partially rescue RdRp activity. Our data show that adenosine and remdesivir triphosphates promote the synthesis of A-less RNAs, as does ppGpp, while amino acid substitutions at the NiRAN-RdRp interface augment activation, suggesting that ligand binding to the NiRAN catalytic site modulates RdRp activity. The existence of allosterically linked nucleotidyl transferase sites that utilize the same substrates has important implications for understanding the mechanism of SARS-CoV-2 replication and the design of its inhibitors.

Analysis of codon usage of Horseshoe Bat Hepatitis B virus and its host

In the present analysis, codon usage strategies and base distribution of Horseshoe bat hepatitis B virus (HBHBV) were analyzed and compared with its host Rhinolophus sinicus, as no work was yet reported. The magnitude of synonymous codon usage bias (CUB) in the virus and its host was low with higher proportion of the base C. Notably, 21 more frequently used codons, 19 less frequently used codons and 3 underrepresented codons (TCG, ACG and GCG) were found to be similar in both virus and its host coding sequences. Neutrality plot analysis reported greater role of natural selection in HBHBV (67.84%) and R. sinicus (76.90%) over mutation pressure. Base skewness and protein properties also influenced the CUB of genes. Further, codon usage analysis depicted, HBHBV and R. sinicus had many similarities in codon usage patterns that might reflect viral adaptation to its host.

The transcriptional landscape of a rewritten bacterial genome reveals control elements and genome design principles

Sequence rewriting enables low-cost genome synthesis and the design of biological systems with orthogonal genetic codes. The error-free, robust rewriting of nucleotide sequences can be achieved with a complete annotation of gene regulatory elements. Here, we compare transcription in Caulobacter crescentus to transcription from plasmid-borne segments of the synthesized genome of C. ethensis 2.0. This rewritten derivative contains an extensive amount of supposedly neutral mutations, including 123'562 synonymous codon changes. The transcriptional landscape refines 60 promoter annotations, exposes 18 termination elements and links extensive transcription throughout the synthesized genome to the unintentional introduction of sigma factor binding motifs. We reveal translational regulation for 20 CDS and uncover an essential translational regulatory element for the expression of ribosomal protein RplS. The annotation of gene regulatory elements allowed us to formulate design principles that improve design schemes for synthesized DNA, en route to a bright future of iteration-free programming of biological systems.

Transcription, mRNA export and immune evasion shape the codon usage of viruses

The nucleotide composition, dinucleotide composition, and codon usage of many viruses differs from their hosts. These differences arise because viruses are subject to unique mutation and selection pressures that do not apply to host genomes; however, the molecular mechanisms that underlie these evolutionary forces are unclear. Here, we analysed the patterns of codon usage in 1,520 vertebrate-infecting viruses, focusing on parameters known to be under selection and associated with gene regulation. We find that GC content, dinucleotide content, and splicing and m6A modification-related sequence motifs are associated with the type of genetic material (DNA or RNA), strandedness, and replication compartment of viruses. In an experimental follow-up, we find that the effects of GC content on gene expression depend on whether the genetic material is delivered to the cell as DNA or mRNA, whether it is transcribed by endogenous or exogenous RNA polymerase, and whether transcription takes place in the nucleus or cytoplasm. Our results suggest that viral codon usage cannot be explained by a simple adaptation to the codon usage of the host - instead, it reflects the combination of multiple selective and mutational pressures, including the need for efficient transcription, export, and immune evasion.

Codon Bias Can Determine Sorting of a Potassium Channel Protein

Due to the redundancy of the genetic code most amino acids are encoded by multiple synonymous codons. It has been proposed that a biased frequency of synonymous codons can affect the function of proteins by modulating distinct steps in transcription, translation and folding. Here, we use two similar prototype K⁺ channels as model systems to examine whether codon choice has an impact on protein sorting. By monitoring transient expression of GFP-tagged channels in mammalian cells, we find that one of the two channels is sorted in a codon and cell cycle-dependent manner either to mitochondria or the secretory pathway. The data establish that a gene with either rare or frequent codons serves, together with a cell-state-dependent decoding mechanism, as a secondary code for sorting intracellular membrane proteins.

Read between the lines: Diversity of non-translational selection pressures on local codon usage

Protein coding genes can contain specific motifs within their nucleotide sequence that function as a signal for various biological pathways. The presence of such sequence motifs within a gene can have beneficial or detrimental effects on the phenotype and fitness of an organism, and this can lead to the enrichment or avoidance of this sequence motif. The degeneracy of the genetic code allows for the existence of alternative synonymous sequences that exclude or include these motifs, while keeping the encoded amino acid sequence intact. This implies that locally, there can be a selective pressure for preferentially using a codon over its synonymous alternative in order to avoid or enrich a specific sequence motif. This selective pressure could -in addition to mutation, drift and selection for translation efficiency and accuracy- contribute to shape the codon usage bias. In this review, we discuss patterns of avoidance of (or enrichment for) the various biological signals contained in specific nucleotide sequence motifs: transcription and translation initiation and termination signals, mRNA maturation signals, and antiviral immune system targets. Experimental data on the phenotypic or fitness effects of synonymous mutations in these sequence motifs confirm that they can be targets of local selection pressures on codon usage. We also formulate the hypothesis that transposable elements could have a similar impact on codon usage through their preferred integration sequences. Overall, selection on codon usage appears to be a combination of a global selection pressure imposed by the translation machinery, and a patchwork of local selection pressures related to biological signals contained in specific sequence motifs.

Pathway regulation mechanism revealed by cotranslational folding of villin headpiece subdomain HP35

Cotranslational folding is one of the most important features of protein folding in vivo. Although many studies have shown that the folding pathways of cotranslational folding are different from free folding in vitro, the detailed mechanism of folding dynamics is lacking. Here we combine all-atom molecular simulations with an ideal ribosome tunnel model to investigate the cotranslational folding of villin headpiece subdomain HP35. By comparing the folding dynamics between cotranslational folding and free folding, we found that cotranslational folding tends to fold along the pathway that is easier to fold into native state in the latter. In addition, the roles of the ribosome tunnel and sequential folding are analyzed separately. Our results show that the ribosome can prevent the untimely folding of the C segment of HP35 to reduce the non-native interactions, while the translation speed can regulate the amounts of native and non-native interactions and the balance between them. Overall, these results give insights into the general mechanisms of cotranslational protein folding.

Analysis of chromosomes and nucleotides in rice to predict gene expression through codon usage pattern

Amino acids are essential measurements for the potential growth stage because of connecting to protein structures and functions. The objective of this paper was to analyze chromosomes feature at plastid region of rice represented by nucleotide, synonymous codon, and amino acid usage to predict gene expression through codon usage pattern. The results showed that the values of the codon adaption index ranged from 0.733 in chromosome 9 to 0.631 in chromosome 8 with full length of these two chromosomes were 3738 and 1635 respectively. The higher value of guanine and cytosine content was 60% in chromosomes 9 while the lower values was 37% in chromosomes 11. Eight chromosomes (ch1, ch2, ch3, ch5, ch7, ch8, ch10, and ch12) were greater value of modified relative codon bias than threshold (threshold: 0.66) especially in cysteine for ch1, ch2, ch5, ch10, and ch12. While other remaining chromosomes were less than the threshold. Relative synonymous codon usage found that the over-represented of amino acids were asparagine, aspartate, cysteine, glutamate, and phenylalanine across all 12 chromosomes. These results would establish a platform for more and further projects concerning rice breeding and genetics and codon optimization in the amino acids for developing varieties. These results also will help breeders to select desirable genes through the genome for improve target traits.

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

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

Context-dependent and -independent selection on synonymous mutations revealed by 1,135 genomes of Arabidopsis thaliana

BACKGROUND

Synonymous mutations do not alter the amino acids and therefore are regarded as neutral for a long time. However, they do change the tRNA adaptation index (tAI) of a particular codon (independent of its context), affecting the tRNA availability during translation. They could also change the isoaccepting relationship with its neighboring synonymous codons in particular context, which again affects the local translation process. Evidence of selection pressure on synonymous mutations has emerged.

RESULTS

The proposed selection patterns on synonymous mutations are never formally and systematically tested in plant species. We fully take advantage of the SNP data from 1,135 A. thaliana lines, and found that the synonymous mutations that increase tAI or the isoaccepting mutations in isoaccepting codon context tend to have higher derived allele frequencies (DAF) compared to other synonymous mutations of the opposite effects.

CONCLUSIONS

Synonymous mutations are not strictly neutral. The synonymous mutations that increase tAI or the isoaccepting mutations in isoaccepting codon context are likely to be positively selected. We propose the concept of context-dependent and -independent selection on synonymous mutations. These concepts broaden our knowledge of the functional consequences of synonymous mutations, and should be appealing to phytologists and evolutionary biologists.

Allosteric activation of SARS-CoV-2 RdRp by remdesivir triphosphate and other phosphorylated nucleotides

The catalytic subunit of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp), Nsp12, has a unique NiRAN domain that transfers nucleoside monophosphates to the Nsp9 protein. The NiRAN and RdRp modules form a dynamic interface distant from their catalytic sites and both activities are essential for viral replication. We report that codon-optimized (for the pause-free translation) Nsp12 exists in inactive state in which NiRAN/RdRp interactions are broken, whereas translation by slow ribosomes and incubation with accessory Nsp7/8 subunits or NTPs partially rescue RdRp activity. Our data show that adenosine and remdesivir triphosphates promote synthesis of A-less RNAs, as does ppGpp, while amino acid substitutions at the NiRAN/RdRp interface augment activation, suggesting that ligand binding to the NiRAN catalytic site modulates RdRp activity. The existence of allosterically-linked nucleotidyl transferase sites that utilize the same substrates has important implications for understanding the mechanism of SARS-CoV-2 replication and design of its inhibitors.

HIGHLIGHTS

Codon-optimization of Nsp12 triggers misfolding and activity lossSlow translation, accessory Nsp7 and Nsp8 subunits, and NTPs rescue Nsp12Non-substrate nucleotides activate RNA chain synthesis, likely via NiRAN domainCrosstalk between two Nsp12 active sites that bind the same ligands.

Defects in translation-dependent quality control pathways lead to convergent molecular and neurodevelopmental pathology

Translation-dependent quality control pathways such as no-go decay (NGD), non-stop decay (NSD), and nonsense-mediated decay (NMD) govern protein synthesis and proteostasis by resolving non-translating ribosomes and preventing the production of potentially toxic peptides derived from faulty and aberrant mRNAs. However, how translation is altered and the in vivo defects that arise in the absence of these pathways are poorly understood. Here, we show that the NGD/NSD factors Pelo and Hbs1l are critical in mice for cerebellar neurogenesis but expendable for survival of these neurons after development. Analysis of mutant mouse embryonic fibroblasts revealed translational pauses, alteration of signaling pathways, and translational reprogramming. Similar effects on signaling pathways, including mTOR activation, the translatome and mouse cerebellar development were observed upon deletion of the NMD factor Upf2. Our data reveal that these quality control pathways that function to mitigate errors at distinct steps in translation can evoke similar cellular responses.

Prediction of gene expression under drought stress in spring wheat using codon usage pattern

Spring wheat (Triticum aestivum) is a staple food providing sources of essential proteins for human. In fact, gene expressions of wheat play an important role in growth and productivity that are affected by drought stress. The objective of this work focused on analysis gene feature on spring wheat represented by nucleotide and gene expressions under drought stress. It was found that the higher codon adaptation index was in both wheat root and L-galactono-1, 4-lactone dehydrogenase. It was also found that guanine and cytosine content were high (55.56%) in wheat root. Whereas, guanine and cytosine content were low (41.28%) in L-galactono-1, 4-lactone dehydrogenase. Moreover, the higher relative synonymous codon usage value was observed in codon CAA (1.20), GAA (1.33), GAT (1.00), and ATG (1.00) in wheat root and thus about 62.95% of the total variation in relative synonymous codon was explained by principal component analysis. Additionally, high averages frequency number of codon were (above 15.76) in Met, Lys, Ala, Gly, Phe, Asp, Glu, His, and Tyr; whereas, low averages were in remaining amino acids and majority (90%) of modified relative codon bias values was between 0.40 and 0.90. Shortly, calculations and analysis of codon usage pattern under drought stress would help for genetic engineering, molecular evolution, and gene prediction in wheat studies for developing varieties that associate with drought tolerance.

Wobble tRNA modification and hydrophilic amino acid patterns dictate protein fate

Regulation of mRNA translation elongation impacts nascent protein synthesis and integrity and plays a critical role in disease establishment. Here, we investigate features linking regulation of codon-dependent translation elongation to protein expression and homeostasis. Using knockdown models of enzymes that catalyze the mcm5s2 wobble uridine tRNA modification (U34-enzymes), we show that gene codon content is necessary but not sufficient to predict protein fate. While translation defects upon perturbation of U34-enzymes are strictly dependent on codon content, the consequences on protein output are determined by other features. Specific hydrophilic motifs cause protein aggregation and degradation upon codon-dependent translation elongation defects. Accordingly, the combination of codon content and the presence of hydrophilic motifs define the proteome whose maintenance relies on U34-tRNA modification. Together, these results uncover the mechanism linking wobble tRNA modification to mRNA translation and aggregation to maintain proteome homeostasis.

Pathogenicity and virulence of hepatitis A virus

Hepatitis A is an acute infection of the liver, which is mostly asymptomatic in children and increases the severity with age. Although in most patients the infection resolves completely, in a few of them it may follow a prolonged or relapsed course or even a fulminant form. The reason for these different outcomes is unknown, but it is generally accepted that host factors such as the immunological status, age and the occurrence of underlaying hepatic diseases are the main determinants of the severity. However, it cannot be ruled out that some virus traits may also contribute to the severe clinical outcomes. In this review, we will analyze which genetic determinants of the virus may determine virulence, in the context of a paradigmatic virus in terms of its genomic, molecular, replicative, and evolutionary features.

Synonymous mutations that regulate translation speed might play a non-negligible role in liver cancer development

Background

Synonymous mutations do not change the protein sequences. Automatically, they have been regarded as neutral events and are ignored in the mutation-based cancer studies. However, synonymous mutations will change the codon optimality, resulting in altered translational velocity.

Methods

We fully utilized the transcriptome and translatome of liver cancer and normal tissue from ten patients. We profiled the mutation spectrum and examined the effect of synonymous mutations on translational velocity.

Results

Synonymous mutations that increase the codon optimality significantly enhanced the translational velocity, and were enriched in oncogenes. Meanwhile, synonymous mutations decreasing codon optimality slowed down translation, and were enriched in tumor suppressor genes. These synonymous mutations significantly contributed to the translational changes in tumor samples compared to normal samples.

Conclusions

Synonymous mutations might play a role in liver cancer development by altering codon optimality and translational velocity. Synonymous mutations should no longer be ignored in the genome-wide studies.

Adaptation of codon and amino acid use for translational functions in highly expressed cricket genes

BACKGROUND

For multicellular organisms, much remains unknown about the dynamics of synonymous codon and amino acid use in highly expressed genes, including whether their use varies with expression in different tissue types and sexes. Moreover, specific codons and amino acids may have translational functions in highly transcribed genes, that largely depend on their relationships to tRNA gene copies in the genome. However, these relationships and putative functions are poorly understood, particularly in multicellular systems.

RESULTS

Here, we studied codon and amino acid use in highly expressed genes from reproductive and nervous system tissues (male and female gonad, somatic reproductive system, brain and ventral nerve cord, and male accessory glands) in the cricket Gryllus bimaculatus. We report an optimal codon, defined as the codon preferentially used in highly expressed genes, for each of the 18 amino acids with synonymous codons in this organism. The optimal codons were mostly shared among tissue types and both sexes. However, the frequency of optimal codons was highest in gonadal genes. Concordant with translational selection, a majority of the optimal codons had abundant matching tRNA gene copies in the genome, but sometimes obligately required wobble tRNAs. We suggest the latter may comprise a mechanism for slowing translation of abundant transcripts, particularly for cell-cycle genes. Non-optimal codons, defined as those least commonly used in highly transcribed genes, intriguingly often had abundant tRNAs, and had elevated use in a subset of genes with specialized functions (gametic and apoptosis genes), suggesting their use promotes the translational upregulation of particular mRNAs. In terms of amino acids, we found evidence suggesting that amino acid frequency, tRNA gene copy number, and amino acid biosynthetic costs (size/complexity) had all interdependently evolved in this insect model, potentially for translational optimization.

CONCLUSIONS

Collectively, the results suggest a model whereby codon use in highly expressed genes, including optimal, wobble, and non-optimal codons, and their tRNA abundances, as well as amino acid use, have been influenced by adaptation for various functional roles in translation within this cricket. The effects of expression in different tissue types and the two sexes are discussed.

Comprehensive annotation and characterization of planarian tRNA and tRNA-derived fragments (tRFs)

tRNA-derived fragments (tRFs) have recently gained a lot of scientific interest due to their diverse regulatory roles in several cellular processes. However, their function in dynamic biological processes such as development and regeneration remains unexplored. Here, we show that tRFs are dynamically expressed during planarian regeneration, suggesting a possible role for these small RNAs in the regulation of regeneration. In order to characterize planarian tRFs, we first annotated 457 tRNAs in S. mediterranea combining two tRNA prediction algorithms. Annotation of tRNAs facilitated the identification of three main species of tRFs in planarians-the shorter tRF-5s and itRFs, and the abundantly expressed 5'-tsRNAs. Spatial profiling of tRFs in sequential transverse sections of planarians revealed diverse expression patterns of these small RNAs, including those that are enriched in the head and pharyngeal regions. Expression analysis of these tRF species revealed dynamic expression of these small RNAs over the course of regeneration suggesting an important role in planarian anterior and posterior regeneration. Finally, we show that 5'-tsRNA in planaria interact with all three SMEDWI proteins and an involvement of AGO1 in the processing of itRFs. In summary, our findings implicate a novel role for tRFs in planarian regeneration, highlighting their importance in regulating complex systemic processes. Our study adds to the catalog of posttranscriptional regulatory systems in planaria, providing valuable insights on the biogenesis and the function of tRFs in neoblasts and planarian regeneration.

Ribosome states signal RNA quality control

Eukaryotic cells integrate multiple quality control (QC) responses during protein synthesis in the cytoplasm. These QC responses are signaled by slow or stalled elongating ribosomes. Depending on the nature of the delay, the signal may lead to translational repression, messenger RNA decay, ribosome rescue, and/or nascent protein degradation. Here, we discuss how the structure and composition of an elongating ribosome in a troubled state determine the downstream quality control pathway(s) that ensue. We highlight the intersecting pathways involved in RNA decay and the crosstalk that occurs between RNA decay and ribosome rescue.

Analysis of fungal genomes reveals commonalities of intron gain or loss and functions in intron-poor species

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron–exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.

Full-length ribosome density prediction by a multi-input and multi-output model

Translation elongation is regulated by a series of complicated mechanisms in both prokaryotes and eukaryotes. Although recent advance in ribosome profiling techniques has enabled one to capture the genome-wide ribosome footprints along transcripts at codon resolution, the regulatory codes of elongation dynamics are still not fully understood. Most of the existing computational approaches for modeling translation elongation from ribosome profiling data mainly focus on local contextual patterns, while ignoring the continuity of the elongation process and relations between ribosome densities of remote codons. Modeling the translation elongation process in full-length coding sequence (CDS) level has not been studied to the best of our knowledge. In this paper, we developed a deep learning based approach with a multi-input and multi-output framework, named RiboMIMO, for modeling the ribosome density distributions of full-length mRNA CDS regions. Through considering the underlying correlations in translation efficiency among neighboring and remote codons and extracting hidden features from the input full-length coding sequence, RiboMIMO can greatly outperform the state-of-the-art baseline approaches and accurately predict the ribosome density distributions along the whole mRNA CDS regions. In addition, RiboMIMO explores the contributions of individual input codons to the predictions of output ribosome densities, which thus can help reveal important biological factors influencing the translation elongation process. The analyses, based on our interpretable metric named codon impact score, not only identified several patterns consistent with the previously-published literatures, but also for the first time (to the best of our knowledge) revealed that the codons located at a long distance from the ribosomal A site may also have an association on the translation elongation rate. This finding of long-range impact on translation elongation velocity may shed new light on the regulatory mechanisms of protein synthesis. Overall, these results indicated that RiboMIMO can provide a useful tool for studying the regulation of translation elongation in the range of full-length CDS.

Translation elongation is a process in which amino acids are linked into proteins by ribosomes in cells. Translation elongation rates along the mRNAs are not constant, and are regulated by a series of mechanisms, such as codon rarity and mRNA stability. In this study, we modeled the translation elongation process at a full-length coding sequence level and developed a deep learning based approach to predict the translation elongation rates from mRNA sequences, through extracting the regulatory codes of elongation rates from the contextual sequences. The analyses, based on our interpretable metric named codon impact score, for the first time (to the best of our knowledge), revealed that in addition to the neighboring codons of the ribosomal A sites, the remote codons may also have an important impact on the translation elongation rates. This new finding may stimulate additional experiments and shed light on the regulatory mechanisms of protein synthesis.

Easy Synthesis of Complex Biomolecular Assemblies: Wheat Germ Cell-Free Protein Expression in Structural Biology

Cell-free protein synthesis (CFPS) systems are gaining more importance as universal tools for basic research, applied sciences, and product development with new technologies emerging for their application. Huge progress was made in the field of synthetic biology using CFPS to develop new proteins for technical applications and therapy. Out of the available CFPS systems, wheat germ cell-free protein synthesis (WG-CFPS) merges the highest yields with the use of a eukaryotic ribosome, making it an excellent approach for the synthesis of complex eukaryotic proteins including, for example, protein complexes and membrane proteins. Separating the translation reaction from other cellular processes, CFPS offers a flexible means to adapt translation reactions to protein needs. There is a large demand for such potent, easy-to-use, rapid protein expression systems, which are optimally serving protein requirements to drive biochemical and structural biology research. We summarize here a general workflow for a wheat germ system providing examples from the literature, as well as applications used for our own studies in structural biology. With this review, we want to highlight the tremendous potential of the rapidly evolving and highly versatile CFPS systems, making them more widely used as common tools to recombinantly prepare particularly challenging recombinant eukaryotic proteins.

An Overview on the Development of mRNA-Based Vaccines and Their Formulation Strategies for Improved Antigen Expression In Vivo

The mRNA-based vaccine approach is a promising alternative to traditional vaccines due to its ability for prompt development, high potency, and potential for secure administration and low-cost production. Nonetheless, the application has still been limited by the instability as well as the ineffective delivery of mRNA in vivo. Current technological improvements have now mostly overcome these concerns, and manifold mRNA vaccine plans against various forms of malignancies and infectious ailments have reported inspiring outcomes in both humans and animal models. This article summarizes recent mRNA-based vaccine developments, advances of in vivo mRNA deliveries, reflects challenges and safety concerns, and future perspectives, in developing the mRNA vaccine platform for extensive therapeutic use.

OTTER, a new method quantifying absolute amounts of tRNAs

To maintain optimal proteome, both codon choice of each mRNA and supply of aminoacyl-tRNAs are two principal factors in translation. Recent reports have revealed that the amounts of tRNAs in cells are more dynamic than we had expected. High-throughput methods such as RNA-Seq and microarrays are versatile for comprehensive detection of changes in individual tRNA amounts, but they suffer from inability to assess signal production efficiencies of individual tRNA species. Thus, they are not the perfect choice to measure absolute amounts of tRNAs. Here, we introduce a novel method for this purpose, termed Oligonucleotide-directed Three-prime Terminal Extension of RNA (OTTER), which employs fluorescence-labeling at the 3'-terminus of a tRNA by optimized reverse primer extension and an assessment step of each labeling efficiency by northern blotting. Using this method, we quantified the absolute amounts of the 34 individual and 4 pairs of isoacceptor tRNAs out of the total 42 nuclear-encoded isoacceptors in the yeast Saccharomyces cerevisiae. We found that the amounts of tRNAs in log phase yeast cells grown in a rich glucose medium range from 0.030 to 0.73 pmol/µg RNA. The tRNA amounts seem to be altered at the isoacceptor level by a few folds in response to physiological growing conditions. The data obtained by OTTER are poorly correlated with those by simple RNA-Seq, marginally with those by microarrays and by microscale thermophoresis. However, the OTTER data showed good agreement with the data obtained by 2D-gel analysis of in vivo radiolabeled RNAs. Thus, OTTER is a suitable method for quantifying absolute amounts of tRNAs at the level of isoacceptor resolution.

Massively parallel gene expression variation measurement of a synonymous codon library

BACKGROUND

Cell-to-cell variation in gene expression strongly affects population behavior and is key to multiple biological processes. While codon usage is known to affect ensemble gene expression, how codon usage influences variation in gene expression between single cells is not well understood.

RESULTS

Here, we used a Sort-seq based massively parallel strategy to quantify gene expression variation from a green fluorescent protein (GFP) library containing synonymous codons in Escherichia coli. We found that sequences containing codons with higher tRNA Adaptation Index (TAI) scores, and higher codon adaptation index (CAI) scores, have higher GFP variance. This trend is not observed for codons with high Normalized Translation Efficiency Index (nTE) scores nor from the free energy of folding of the mRNA secondary structure. GFP noise, or squared coefficient of variance (CV2), scales with mean protein abundance for low-abundant proteins but does not change at high mean protein abundance.

CONCLUSIONS

Our results suggest that the main source of noise for high-abundance proteins is likely not originating at translation elongation. Additionally, the drastic change in mean protein abundance with small changes in protein noise seen from our library implies that codon optimization can be performed without concerning gene expression noise for biotechnology applications.

Integrated analyses of translatome and proteome identify the rules of translation selectivity in RPS14-deficient cells

In ribosomopathies, the Diamond-Blackfan anemia (DBA) or 5q- syndrome, ribosomal protein (RP) genes are affected by mutation or deletion, resulting in bone marrow erythroid hypoplasia. Unbalanced production of ribosomal subunits leading to a limited ribosome cellular content regulates translation at the expense of the master erythroid transcription factor GATA1. In RPS14-deficient cells mimicking 5q- syndrome erythroid defects, we show that the transcript length, codon bias of the coding sequence (CDS) and 3’UTR (untranslated region) structure are the key determinants of translation. In these cells, short transcripts with a structured 3’UTR and high codon adaptation index (CAI) showed a decreased translation efficiency. Quantitative analysis of the whole proteome confirmed that the post-transcriptional changes depended on the transcript characteristics that governed the translation efficiency in conditions of low ribosome availability. In addition, proteins involved in normal erythroid differentiation share most determinants of translation selectivity. Our findings thus indicate that impaired erythroid maturation due to 5q- syndrome may proceed from a translational selectivity at the expense of the erythroid differentiation program, and suggest that an interplay between the CDS and UTR may regulate 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.

Genome-wide characterization of PEBP family genes in nine Rosaceae tree species and their expression analysis in P. mume

BACKGROUND

Phosphatidylethanolamine-binding proteins (PEBPs) constitute a common gene family found among animals, plants and microbes. Plant PEBP proteins play an important role in regulating flowering time, plant architecture as well as seed dormancy. Though PEBP family genes have been well studied in Arabidopsis and other model species, less is known about these genes in perennial trees.

RESULTS

To understand the evolution of PEBP genes and their functional roles in flowering control, we identified 56 PEBP members belonging to three gene clades (MFT-like, FT-like, and TFL1-like) and five lineages (FT, BFT, CEN, TFL1, and MFT) across nine Rosaceae perennial species. Structural analysis revealed highly conserved gene structure and protein motifs among Rosaceae PEBP proteins. Codon usage analysis showed slightly biased codon usage across five gene lineages. With selection pressure analysis, we detected strong purifying selection constraining divergence within most lineages, while positive selection driving the divergence of FT-like and TFL1-like genes from the MFT-like gene clade. Spatial and temporal expression analyses revealed the essential role of FT in regulating floral bud breaking and blooming in P. mume. By employing a weighted gene co-expression network approach, we inferred a putative FT regulatory module required for dormancy release and blooming in P. mume.

CONCLUSIONS

We have characterized the PEBP family genes in nine Rosaceae species and examined their phylogeny, genomic syntenic relationship, duplication pattern, and expression profiles during flowering process. These results revealed the evolutionary history of PEBP genes and their functions in regulating floral bud development and blooming among Rosaceae tree species.

Potential role of the X circular code in the regulation of gene expression

The X circular code is a set of 20 trinucleotides (codons) that has been identified in the protein-coding genes of most organisms (bacteria, archaea, eukaryotes, plasmids, viruses). It has been shown previously that the X circular code has the important mathematical property of being an error-correcting code. Thus, motifs of the X circular code, i.e. a series of codons belonging to X and called X motifs, allow identification and maintenance of the reading frame in genes. X motifs are significantly enriched in protein-coding genes, but have also been identified in many transfer RNA (tRNA) genes and in important functional regions of the ribosomal RNA (rRNA), notably in the peptidyl transferase center and the decoding center. Here, we investigate the potential role of X motifs as functional elements of protein-coding genes. First, we identify the codons of the X circular code which are frequent or rare in each domain of life (archaea, bacteria, eukaryota) and show that, for the amino acids with the highest codon bias, the preferred codon is often an X codon. We also observe a correlation between the 20 X codons and the optimal codons/dicodons that have been shown to influence translation efficiency. Then, we examined recently published experimental results concerning gene expression levels in diverse organisms. The approach used is the analysis of X motifs according to their density d_s(X), i.e. the number of X motifs per kilobase in a gene sequence s. Surprisingly, this simple parameter identifies several unexpected relations between the X circular code and gene expression. For example, the X motifs are significantly enriched in the minimal gene set belonging to the three domains of life, and in codon-optimized genes. Furthermore, the density of X motifs generally correlates with experimental measures of translation efficiency and mRNA stability. Taken together, these results lead us to propose that the X motifs may represent a genetic signal contributing to the maintenance of the correct reading frame and the optimization and regulation of gene expression.

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.

Druggable Transient Pockets in Protein Kinases

Drug discovery using small molecule inhibitors is reaching a stalemate due to low selectivity, adverse off-target effects and inevitable failures in clinical trials. Conventional chemical screening methods may miss potent small molecules because of their use of simple but outdated kits composed of recombinant enzyme proteins. Non-canonical inhibitors targeting a hidden pocket in a protein have received considerable research attention. Kii and colleagues identified an inhibitor targeting a transient pocket in the kinase DYRK1A during its folding process and termed it FINDY. FINDY exhibits a unique inhibitory profile; that is, FINDY does not inhibit the fully folded form of DYRK1A, indicating that the FINDY-binding pocket is hidden in the folded form. This intriguing pocket opens during the folding process and then closes upon completion of folding. In this review, we discuss previously established kinase inhibitors and their inhibitory mechanisms in comparison with FINDY. We also compare the inhibitory mechanisms with the growing concept of “cryptic inhibitor-binding sites.” These sites are buried on the inhibitor-unbound surface but become apparent when the inhibitor is bound. In addition, an alternative method based on cell-free protein synthesis of protein kinases may allow the discovery of small molecules that occupy these mysterious binding sites. Transitional folding intermediates would become alternative targets in drug discovery, enabling the efficient development of potent kinase inhibitors.

Use of Synonymous Deoptimization to Derive Modified Live Attenuated Strains of Foot and Mouth Disease Virus

Foot-and-mouth disease (FMD) is one of the most economically important viral diseases that can affect livestock. In the last 70 years, use of an inactivated whole antigen vaccine has contributed to the eradication of disease from many developed nations. However, recent outbreaks in Europe and Eastern Asia demonstrated that infection can spread as wildfire causing economic and social devastation. Therefore, it is essential to develop new control strategies that could confer early protection and rapidly stop disease spread. Live attenuated vaccines (LAV) are one of the best choices to obtain a strong early and long-lasting protection against viral diseases. In proof of concept studies, we previously demonstrated that “synonymous codon deoptimization” could be applied to the P1 capsid coding region of the viral genome to derive attenuated FMDV serotype A12 strains. Here, we demonstrate that a similar approach can be extended to the highly conserved non-structural P2 and P3 coding regions, providing a backbone for multiple serotype FMDV LAV development. Engineered codon deoptimized P2, P3 or P2, and P3 combined regions were included into the A₂₄Cruzeiro infectious clone optimized for vaccine production, resulting in viable progeny that exhibited different degrees of attenuation in cell culture, in mice, and in the natural host (swine). Derived strains were thoroughly characterized in vitro and in vivo. Our work demonstrates that overall, the entire FMDV genome tolerates codon deoptimization, highlighting the potential of using this technology to derive novel improved LAV candidates.

Cost-efficiency tradeoff is optimized in various cancer types revealed by genome-wide analysis

The tradeoff between cost and efficiency is omnipresent in organisms. Specifically, how the evolutionary force shapes the tradeoff between biosynthetic cost and translation efficiency remains unclear. In the cancer community, whether the adjustment of cost-efficiency tradeoff acts as a strategy to facilitate tumor proliferation and contributes to oncogenesis is uninvestigated. To address this issue, we retrieved the gene expression profile in various cancer types and the matched normal samples from The Cancer Genome Atlas (TCGA). We found that the highly expressed genes in cancers generally have higher tAI/nitro ratios than those in normal samples. This is possibly caused by the higher tAI/nitro ratios observed in oncogenes than tumor suppressor genes (TSG). Furthermore, in the cancer samples, derived mutations in oncogenes usually lead to higher tAI/nitro ratios, while those mutations in TSG lead to lower tAI/nitro. For a special case of kidney cancer, we investigated several crucial genes in tumor samples versus normal samples, and discovered that the changes in tAI/nitro ratios are correlated with the changes in translation level. Our study for the first time revealed the optimization of cost-efficiency tradeoff in cancers. The cost-efficiency dilemma is optimized by the tumor cells, and is possibly beneficial for the translation and production of oncogenes, and eventually contributes to proliferation and oncogenesis. Our findings could provide novel perspectives in depicting the cancer genomes and might help unravel the cancer evolution.