Codon optimality, bias and usage in translation and mRNA decay

Deep learning suggests that gene expression is encoded in all parts of a co-evolving interacting gene regulatory structure

Understanding the genetic regulatory code governing gene expression is an important challenge in molecular biology. However, how individual coding and non-coding regions of the gene regulatory structure interact and contribute to mRNA expression levels remains unclear. Here we apply deep learning on over 20,000 mRNA datasets to examine the genetic regulatory code controlling mRNA abundance in 7 model organisms ranging from bacteria to Human. In all organisms, we can predict mRNA abundance directly from DNA sequence, with up to 82% of the variation of transcript levels encoded in the gene regulatory structure. By searching for DNA regulatory motifs across the gene regulatory structure, we discover that motif interactions could explain the whole dynamic range of mRNA levels. Co-evolution across coding and non-coding regions suggests that it is not single motifs or regions, but the entire gene regulatory structure and specific combination of regulatory elements that define gene expression levels.

Base editing: advances and therapeutic opportunities

Base editing - the introduction of single-nucleotide variants (SNVs) into DNA or RNA in living cells - is one of the most recent advances in the field of genome editing. As around half of known pathogenic genetic variants are due to SNVs, base editing holds great potential for the treatment of numerous genetic diseases, through either temporary RNA or permanent DNA base alterations. Recent advances in the specificity, efficiency, precision and delivery of DNA and RNA base editors are revealing exciting therapeutic opportunities for these technologies. We expect the correction of single point mutations will be a major focus of future precision medicine.

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

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

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

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

Mutation bias within oncogene families is related to proliferation-specific codon usage

It is well known that in cancer gene families some members are more frequently mutated in tumor samples than their family counterparts. A paradigmatic case of this phenomenon is KRAS from the RAS family. Different explanations have been proposed ranging from differential interaction with other proteins to preferential expression or localization. Interestingly, it has been described that despite the high amino acid identity between RAS family members, KRAS employs an intriguing differential codon usage. Here, we found that this phenomenon is not exclusive to the RAS family. Indeed, in the RAS family and other oncogene families with two or three members, the most prevalently mutated gene in tumor samples employs a differential codon usage that is characteristic of genes involved in proliferation. Prompted by these observations, we chose the RAS family to experimentally demonstrate that the translation efficiency of oncogenes that are preferentially mutated in tumor samples is increased in proliferative cells compared to quiescent cells. These results were further validated by assessing the translation efficiency of KRAS in cell lines that differ in their tRNA expression profile. These differences are related to the cell division rate of the studied cells and thus suggest an important role in context-specific oncogene expression regulation. Altogether, our study demonstrates that dynamic translation programs contribute to shaping the expression profiles of oncogenes. Therefore, we propose this codon bias as a regulatory layer to control cell context-specific expression and explain the differential prevalence of mutations in certain members of oncogene families.

CUBAP: an interactive web portal for analyzing codon usage biases across populations

Synonymous codon usage significantly impacts translational and transcriptional efficiency, gene expression, the secondary structure of both mRNA and proteins, and has been implicated in various diseases. However, population-specific differences in codon usage biases remain largely unexplored. Here, we present a web server, https://cubap.byu.edu, to facilitate analyses of codon usage biases across populations (CUBAP). Using the 1000 Genomes Project, we calculated and visually depict population-specific differences in codon frequencies, codon aversion, identical codon pairing, co-tRNA codon pairing, ramp sequences, and nucleotide composition in 17,634 genes. We found that codon pairing significantly differs between populations in 35.8% of genes, allowing us to successfully predict the place of origin for African and East Asian individuals with 98.8% and 100% accuracy, respectively. We also used CUBAP to identify a significant bias toward decreased CTG pairing in the immunity related GTPase M (IRGM) gene in East Asian and African populations, which may contribute to the decreased association of rs10065172 with Crohn's disease in those populations. CUBAP facilitates in-depth gene-specific and codon-specific visualization that will aid in analyzing candidate genes identified in genome-wide association studies, identifying functional implications of synonymous variants, predicting population-specific impacts of synonymous variants and categorizing genetic biases unique to certain populations.

Principles of mRNA control by human PUM proteins elucidated from multimodal experiments and integrative data analysis

The human PUF-family proteins, PUM1 and PUM2, posttranscriptionally regulate gene expression by binding to a PUM recognition element (PRE) in the 3'-UTR of target mRNAs. Hundreds of PUM1/2 targets have been identified from changes in steady-state RNA levels; however, prior studies could not differentiate between the contributions of changes in transcription and RNA decay rates. We applied metabolic labeling to measure changes in RNA turnover in response to depletion of PUM1/2, showing that human PUM proteins regulate expression almost exclusively by changing RNA stability. We also applied an in vitro selection workflow to precisely identify the binding preferences of PUM1 and PUM2. By integrating our results with prior knowledge, we developed a "rulebook" of key contextual features that differentiate functional versus nonfunctional PREs, allowing us to train machine learning models that accurately predict the functional regulation of RNA targets by the human PUM proteins.

Comparative analysis of protein synthesis rate in COVID-19 with other human coronaviruses

The genetic code contains information that impacts the efficiency and rate of translation. Translation elongation plays a crucial role in determining the composition of the proteome, errors within a protein contributes towards disease processes. It is important to analyze the novel coronavirus (2019-nCoV) at the codon level to find similarities and variations in hosts to compare with other human coronavirus (CoVs). This requires a comparative and comprehensive study of various human and zoonotic nature CoVs relating to codon usage bias, relative synonymous codon usage (RSCU), proportions of slow codons, and slow di-codons, the effective number of codons (ENC), mutation bias, codon adaptation index (CAI), and codon frequencies. In this work, seven different CoVs were analyzed to determine the protein synthesis rate and the adaptation of these viruses to the host cell. The result reveals that the proportions of slow codons and slow di-codons in human host of 2019-nCoV and SARS-CoV found to be similar and very less compared to the other five coronavirus types, which suggest that the 2019-nCoV and SARS-CoV have faster protein synthesis rate. Zoonotic CoVs have high RSCU and codon adaptation index than human CoVs which implies the high translation rate in zoonotic viruses. All CoVs have more AT% than GC% in genetic codon compositions. The average ENC values of seven CoVs ranged between 38.36 and 49.55, which implies the CoVs are highly conserved and are easily adapted to host cells. The mutation rate of 2019-nCoV is comparatively less than MERS-CoV and NL63 that shows an evidence for genetic diversity. Host-specific codon composition analysis portrays the relation between viral host sequences and the capability of novel virus replication in host cells. Moreover, the analysis provides useful measures for evaluating a virus-host adaptation, transmission potential of novel viruses, and thus contributes to the strategies of anti-viral drug design.

CRISPR/Cas9-mediated generation and analysis of N terminus polymorphic models of β2AR in isogenic hPSC-derived cardiomyocytes

During normal- and patho-physiological situations, the behavior of the beta2-adrenoreceptor (β₂AR) is influenced by polymorphic variants. The functional impact of such polymorphisms has been suggested from data derived from genetic association studies, in vitro experiments with primary cells, and transgenic overexpression models. However, heterogeneous genetic background and non-physiological transgene expression levels confound interpretation, leading to conflicting mechanistic conclusions. To overcome these limitations, we used CRISPR/Cas9 gene editing technology in human pluripotent stem cells (hPSCs) to create a unique suite of four isogenic homozygous variants at amino acid positions 16(G/R) and 27(G/Q), which reside in the N terminus of the β₂AR. By producing cardiomyocytes from these hPSC lines, we determined that at a functional level β₂AR signaling dominated over β₁AR . Examining changes in beat rates and responses to isoprenaline, Gi coupling, cyclic AMP (cAMP) production, downregulation, and desensitization indicated that responses were often heightened for the GE variant, implying differential dominance of both polymorphic location and amino acid substitution. This finding was corroborated, since GE showed hypersensitivity to doxorubicin-induced cardiotoxicity relative to GQ and RQ variants. Thus, understanding the effect of β₂AR polymorphisms on cardiac response to anticancer therapy may provide a route for personalized medicine and facilitate immediate clinical impact.

Formation of HERV-K and HERV-Fc1 Envelope Family Members is Suppressed on Transcriptional and Translational Level

The human genome comprises 8% sequences of retroviral origin, so-called human endogenous retroviruses (HERVs). Most of these proviral sequences are defective, but some possess open reading frames. They can lead to the formation of viral transcripts, when activated by intrinsic and extrinsic factors. HERVs are thought to play a pathological role in inflammatory diseases and cancer. Since the consequences of activated proviral sequences in the human body are largely unexplored, selected envelope proteins of human endogenous retroviruses associated with inflammatory diseases, namely HERV-K18, HERV-K113, and HERV-Fc1, were investigated in the present study. A formation of glycosylated envelope proteins was demonstrated in different mammalian cell lines. Nevertheless, protein maturation seemed to be incomplete as no transport to the plasma membrane was observed. Instead, the proteins remained in the ER where they induced the expression of genes involved in unfolded protein response, such as HSPA5 and sXBP1. Furthermore, low expression levels of native envelope proteins were increased by codon optimization. Cell-free expression systems showed that both the transcriptional and translational level is affected. By generating different codon-optimized variants of HERV-K113 envelope, the influence of single rare t-RNA pools in certain cell lines was demonstrated. The mRNA secondary structure also appears to play an important role in the translation of the tested viral envelope proteins. In summary, the formation of certain HERV proteins is basically possible. However, their complete maturation and thus full biologic activity seems to depend on additional factors that might be disease-specific and await elucidation in the future.

Large-scale in silico mutagenesis experiments reveal optimization of genetic code and codon usage for protein mutational robustness

BACKGROUND

How, and the extent to which, evolution acts on DNA and protein sequences to ensure mutational robustness and evolvability is a long-standing open question in the field of molecular evolution. We addressed this issue through the first structurome-scale computational investigation, in which we estimated the change in folding free energy upon all possible single-site mutations introduced in more than 20,000 protein structures, as well as through available experimental stability and fitness data.

RESULTS

At the amino acid level, we found the protein surface to be more robust against random mutations than the core, this difference being stronger for small proteins. The destabilizing and neutral mutations are more numerous in the core and on the surface, respectively, whereas the stabilizing mutations are about 4% in both regions. At the genetic code level, we observed smallest destabilization for mutations that are due to substitutions of base III in the codon, followed by base I, bases I+III, base II, and other multiple base substitutions. This ranking highly anticorrelates with the codon-anticodon mispairing frequency in the translation process. This suggests that the standard genetic code is optimized to limit the impact of random mutations, but even more so to limit translation errors. At the codon level, both the codon usage and the usage bias appear to optimize mutational robustness and translation accuracy, especially for surface residues.

CONCLUSION

Our results highlight the non-universality of mutational robustness and its multiscale dependence on protein features, the structure of the genetic code, and the codon usage. Our analyses and approach are strongly supported by available experimental mutagenesis data.

The Ongoing Quest to Crack the Genetic Code for Protein Production

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

RNA-binding proteins balance brain function in health and disease

Posttranscriptional gene expression including splicing, RNA transport, translation, and RNA decay provides an important regulatory layer in many if not all molecular pathways. Research in the last decades has positioned RNA-binding proteins (RBPs) right in the center of posttranscriptional gene regulation. Here, we propose interdependent networks of RBPs to regulate complex pathways within the central nervous system (CNS). These are involved in multiple aspects of neuronal development and functioning, including higher cognition. Therefore, it is not sufficient to unravel the individual contribution of a single RBP and its consequences but rather to study and understand the tight interplay between different RBPs. In this review, we summarize recent findings in the field of RBP biology and discuss the complex interplay between different RBPs. Second, we emphasize the underlying dynamics within an RBP network and how this might regulate key processes such as neurogenesis, synaptic transmission, and synaptic plasticity. Importantly, we envision that dysfunction of specific RBPs could lead to perturbation within the RBP network. This would have direct and indirect (compensatory) effects in mRNA binding and translational control leading to global changes in cellular expression programs in general and in synaptic plasticity in particular. Therefore, we focus on RBP dysfunction and how this might cause neuropsychiatric and neurodegenerative disorders. Based on recent findings, we propose that alterations in the entire regulatory RBP network might account for phenotypic dysfunctions observed in complex diseases including neurodegeneration, epilepsy, and autism spectrum disorders.

Translational control in aging and neurodegeneration

Protein metabolism plays central roles in age-related decline and neurodegeneration. While a large body of research has explored age-related changes in protein degradation, alterations in the efficiency and fidelity of protein synthesis with aging are less well understood. Age-associated changes occur in both the protein synthetic machinery (ribosomal proteins and rRNA) and within regulatory factors controlling translation. At the same time, many of the interventions that prolong lifespan do so in part by pre-emptively decreasing protein synthesis rates to allow better harmonization to age-related declines in protein catabolism. Here we review the roles of translation regulation in aging, with a specific focus on factors implicated in age-related neurodegeneration. We discuss how emerging technologies such as ribosome profiling and superior mass spectrometric approaches are illuminating age-dependent mRNA-specific changes in translation rates across tissues to reveal a critical interplay between catabolic and anabolic pathways that likely contribute to functional decline. These new findings point to nodes in posttranscriptional gene regulation that both contribute to aging and offer targets for therapy. This article is categorized under: Translation > Translation Regulation Translation > Ribosome Biogenesis Translation > Translation Mechanisms.

GCN sensitive protein translation in yeast

Levels of protein translation by ribosomes are governed both by features of the translation machinery as well as sequence properties of the mRNAs themselves. We focus here on a striking three-nucleotide periodicity, characterized by overrepresentation of GCN codons and underrepresentation of G at the second position of codons, that is observed in Open Reading Frames (ORFs) of mRNAs. Our examination of mRNA sequences in Saccharomyces cerevisiae revealed that this periodicity is particularly pronounced in the initial codons-the ramp region-of ORFs of genes with high protein expression. It is also found in mRNA sequences immediately following non-standard AUG start sites, located upstream or downstream of the standard annotated start sites of genes. To explore the possible influences of the ramp GCN periodicity on translation efficiency, we tested edited ramps with accentuated or depressed periodicity in two test genes, SKN7 and HMT1. Greater conformance to (GCN)n was found to significantly depress translation, whereas disrupting conformance had neutral or positive effects on translation. Our recent Molecular Dynamics analysis of a subsystem of translocating ribosomes in yeast revealed an interaction surface that H-bonds to the +1 codon that is about to enter the ribosome decoding center A site. The surface, comprised of 16S/18S rRNA C1054 and A1196 (E. coli numbering) and R146 of ribosomal protein Rps3, preferentially interacts with GCN codons, and we hypothesize that modulation of this mRNA-ribosome interaction may underlie GCN-mediated regulation of protein translation. Integration of our expression studies with large-scale reporter studies of ramp sequence variants suggests a model in which the C1054-A1196-R146 (CAR) interaction surface can act as both an accelerator and braking system for ribosome translation.

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

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

The interface between coronaviruses and host cell RNA biology: Novel potential insights for future therapeutic intervention

Coronaviruses, including SARS-Cov-2, are RNA-based pathogens that interface with a large variety of RNA-related cellular processes during infection. These processes include capping, polyadenylation, localization, RNA stability, translation, and regulation by RNA binding proteins or noncoding RNA effectors. The goal of this article is to provide an in-depth perspective on the current state of knowledge of how various coronaviruses interact with, usurp, and/or avoid aspects of these cellular RNA biology machineries. A thorough understanding of how coronaviruses interact with RNA-related posttranscriptional processes in the cell should allow for new insights into aspects of viral pathogenesis as well as identify new potential avenues for the development of anti-coronaviral therapeutics. This article is categorized under: RNA in Disease and Development > RNA in Disease.

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

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

Quantitative tRNA-sequencing uncovers metazoan tissue-specific tRNA regulation

Transfer RNAs (tRNA) are quintessential in deciphering the genetic code; disseminating nucleic acid triplets into correct amino acid identity. While this decoding function is clear, an emerging theme is that tRNA abundance and functionality can powerfully impact protein production rate, folding, activity, and messenger RNA stability. Importantly, however, the expression pattern of tRNAs is obliquely known. Here we present Quantitative Mature tRNA sequencing (QuantM-tRNA seq), a technique to monitor tRNA abundance and sequence variants secondary to RNA modifications. With QuantM-tRNA seq, we assess the tRNA transcriptome in mammalian tissues. We observe dramatic distinctions in isodecoder expression and known tRNA modifications between tissues. Remarkably, despite dramatic changes in tRNA isodecoder gene expression, the overall anticodon pool of each tRNA family is similar across tissues. These findings suggest that while anticodon pools appear to be buffered via an unknown mechanism, underlying transcriptomic and epitranscriptomic differences suggest a more complex tRNA regulatory landscape.

The relative abundance of specific tRNA can impact protein production rate, folding, and messenger RNA stability. Here the authors describe QuantM-tRNA seq — a method to monitor tRNA abundance and sequence variants — and uncover distinctions in isodecoder expression between tissues that are independent of the anticodon pool of each tRNA family.

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

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

Correlation between gene expression levels under drought stress and synonymous codon usage in rice plant by in-silico study

We studied the correlation of synonymous codon usage (SCU) on gene expression levels under drought stress in rice. Sixty genes related to drought stress (with high, intermediate and low expression) were selected from rice meta-analysis data and various codon usage indices such as the effective number of codon usage (ENC), codon adaptation index (CAI) and relative synonymous codon usage (RSCU) were calculated. We found that in genes highly expressing under drought 1) GC content was higher, 2) ENC value was lower, 3) the preferred codons of some amino acids changed and 4) the RSCU ratio of GC-end codons relative to AT-end codons for 18 amino acids increased significantly compared with those in other genes. We introduce ARSCU as the Average ratio of RSCUs of GC-end codons to AT-end codons in each gene that could significantly separate high-expression genes under drought from low-expression genes. ARSCU is calculated using the program ARSCU-Calculator developed by our group to help predicting expression level of rice genes under drought. An index above ARSCU threshold is expected to indicate that the gene under study may belong to the "high expression group under drought". This information may be applied for codon optimization of genes for rice genetic engineering. To validate these findings, we further used 60 other genes (randomly selected subset of 43233 genes studied for their response to drought stress). ARSCU value was able to predict the level of expression at 88.33% of the cases. Using third set of 60 genes selected amongst high expressing genes not related to drought, only 31.65% of the genes showed ARSCU value of higher than the set threshold. This indicates that the phenomenon we described in this report may be unique for drought related genes. To justify the observed correlation between CUB and high expressing genes under drought, possible role of tRNA post transcriptional modification and tRFs was hypothesized as possible underlying biological mechanism.

Splicing-accessible coding 3'UTRs control protein stability and interaction networks

BACKGROUND

3'-Untranslated regions (3'UTRs) play crucial roles in mRNA metabolism, such as by controlling mRNA stability, translation efficiency, and localization. Intriguingly, in some genes the 3'UTR is longer than their coding regions, pointing to additional, unknown functions. Here, we describe a protein-coding function of 3'UTRs upon frameshift-inducing alternative splicing in more than 10% of human and mouse protein-coding genes.

RESULTS

3'UTR-encoded amino acid sequences show an enrichment of PxxP motifs and lead to interactome rewiring. Furthermore, an elevated proline content increases protein disorder and reduces protein stability, thus allowing splicing-controlled regulation of protein half-life. This could also act as a surveillance mechanism for erroneous skipping of penultimate exons resulting in transcripts that escape nonsense mediated decay. The impact of frameshift-inducing alternative splicing on disease development is emphasized by a retinitis pigmentosa-causing mutation leading to translation of a 3'UTR-encoded, proline-rich, destabilized frameshift-protein with altered protein-protein interactions.

CONCLUSIONS

We describe a widespread, evolutionarily conserved mechanism that enriches the mammalian proteome, controls protein expression and protein-protein interactions, and has important implications for the discovery of novel, potentially disease-relevant protein variants.

Growth-Optimized Aminoacyl-tRNA Synthetase Levels Prevent Maximal tRNA Charging

Aminoacyl-tRNA synthetases (aaRSs) serve a dual role in charging tRNAs. Their enzymatic activities both provide protein synthesis flux and reduce uncharged tRNA levels. Although uncharged tRNAs can negatively impact bacterial growth, substantial concentrations of tRNAs remain deacylated even under nutrient-rich conditions. Here, we show that tRNA charging in Bacillus subtilis is not maximized due to optimization of aaRS production during rapid growth, which prioritizes demands in protein synthesis over charging levels. The presence of uncharged tRNAs is alleviated by precisely tuned translation kinetics and the stringent response, both insensitive to aaRS overproduction but sharply responsive to underproduction, allowing for just enough aaRS production atop a "fitness cliff." Notably, we find that the stringent response mitigates fitness defects at all aaRS underproduction levels even without external starvation. Thus, adherence to minimal, flux-satisfying protein production drives limited tRNA charging and provides a basis for the sensitivity and setpoints of an integrated growth-control network.

COSMO: A dynamic programming algorithm for multicriteria codon optimization

Codon optimization in protein-coding sequences (CDSs) is a widely used technique to promote the heterologous expression of target genes. In codon optimization, a combinatorial space of nucleotide sequences that code a given amino acid sequence and take into account user-prescribed forbidden sequence motifs is explored to optimize multiple criteria. Although evolutionary algorithms have been used to tackle such complex codon optimization problems, evolutionary codon optimization tools do not provide guarantees to find the optimal solutions for these multicriteria codon optimization problems.

We have developed a novel multicriteria dynamic programming algorithm, COSMO. By using this algorithm, we can obtain all Pareto-optimal solutions for the multiple features of CDS, which include codon usage, codon context, and the number of hidden stop codons. User-prescribed forbidden sequence motifs are rigorously excluded from the Pareto-optimal solutions. To accelerate CDS design by COSMO, we introduced constraints that reduce the number of Pareto-optimal solutions to be processed in a branch-and-bound manner. We benchmarked COSMO for run-time and the number of generated solutions by adapting selected human genes to yeast codon usage frequencies, and found that the constraints effectively reduce the run-time. In addition to the benchmarking of COSMO, a multi-objective genetic algorithm (MOGA) for CDS design was also benchmarked for the same two aspects and their performances were compared. In this comparison, (i) MOGA identified significantly fewer Pareto-optimal solutions than COSMO, and (ii) the MOGA solutions did not achieve the same mean hypervolume values as those provided by COSMO. These results suggest that generating the whole set of the Pareto-optimal solutions of the codon optimization problems is a difficult task for MOGA.

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.

Design and Experimental Evaluation of a Minimal, Innocuous Watermarking Strategy to Distinguish Near-Identical DNA and RNA Sequences

The construction of powerful cell factories requires intensive and extensive remodelling of microbial genomes. Considering the rapidly increasing number of these synthetic biology endeavors, there is an increasing need for DNA watermarking strategies that enable the discrimination between synthetic and native gene copies. While it is well documented that codon usage can affect translation, and most likely mRNA stability in eukaryotes, remarkably few quantitative studies explore the impact of watermarking on transcription, protein expression, and physiology in the popular model and industrial yeast Saccharomyces cerevisiae. The present study, using S. cerevisiae as eukaryotic paradigm, designed, implemented, and experimentally validated a systematic strategy to watermark DNA with minimal alteration of yeast physiology. The 13 genes encoding proteins involved in the major pathway for sugar utilization (i.e., glycolysis and alcoholic fermentation) were simultaneously watermarked in a yeast strain using the previously published pathway swapping strategy. Carefully swapping codons of these naturally codon optimized, highly expressed genes, did not affect yeast physiology and did not alter transcript abundance, protein abundance, and protein activity besides a mild effect on Gpm1. The markerQuant bioinformatics method could reliably discriminate native from watermarked genes and transcripts. Furthermore, presence of watermarks enabled selective CRISPR/Cas genome editing, specifically targeting the native gene copy while leaving the synthetic, watermarked variant intact. This study offers a validated strategy to simply watermark genes in S. cerevisiae.

TLR7 dosage polymorphism shapes interferogenesis and HIV-1 acute viremia in women

Type I IFN (IFN-I) production by plasmacytoid DCs (pDCs) occurs during acute HIV-1 infection in response to TLR7 stimulation, but the role of pDC-derived IFN-I in controlling or promoting HIV-1 infection is ambiguous. We report here a sex-biased interferogenic phenotype for a frequent single-nucleotide polymorphism of human TLR7, rs179008, displaying an impact on key parameters of acute HIV-1 infection. We show allele rs179008 T to determine lower TLR7 protein abundance in cells from women, specifically - likely by diminishing TLR7 mRNA translation efficiency through codon usage. The hypomorphic TLR7 phenotype is mirrored by decreased TLR7-driven IFN-I production by female pDCs. Among women from the French ANRS PRIMO cohort of acute HIV-1 patients, carriage of allele rs179008 T associated with lower viremia, cell-associated HIV-1 DNA, and CXCL10 (IP-10) plasma concentrations. RNA viral load was decreased by 0.85 log10 (95% CI, -1.51 to -0.18) among T/T homozygotes, who also exhibited a lower frequency of acute symptoms. TLR7 emerges as an important control locus for acute HIV-1 viremia, and the clinical phenotype for allele rs179008 T, carried by 30%-50% of European women, supports a beneficial effect of toning down TLR7-driven IFN-I production by pDCs during acute HIV-1 infection.

Tracing Brucella evolutionary dynamics in expanding host ranges through nucleotide, codon and amino acid usages in genomes

The host range of Brucella organisms has expanded from terrestrial and marine mammals to fish and amphibians. The high homology genomes of different Brucella organisms promote us to investigate evolutionary patterns for nucleotide, codon and amino acid usage patterns at gene levels among Brucella species. Although the similar patterns for nucleotide and synonymous codon usages exist in gene population, GC composition at the first codon position has significant correlations to that of the second and third codon positions, respectively, suggesting that nucleotide usages surrounding one codon influence synonymous codon usage patterns. Evolutionary patterns represented by synonymous codon and amino acid usages reflect host factor impacting Brucella speciation. As for genetic variations of important virulent factors involved with different biological functions, genes encoding lipoplysaccharides (LPSs) display more distinctive codon adaptation to Brucella than those of the BvrR/BvrS system and type IV secretion system. By Bayesian analysis, the polygenetic constructions for these genes of virulent factors shared by Brucella species display the purifying/positive selections and partially host factor in mediating genetic variations of these genes. The systemic analyses for nucleotide, synonymous codon and amino acid usages at gene level and genetic variations of important virulent factor genes display that host limitation influences either genetic characterizations at gene level or a particular gene involved in virulent factors of Brucella.Communicated by Ramaswamy H. Sarma.

Emerging Roles of tRNAs in RNA Virus Infections

Viruses rely on the host cell translation machinery for efficient synthesis of their own proteins. Emerging evidence highlights different roles for host transfer RNAs (tRNAs) in the process of virus replication. For instance, different RNA viruses manipulate host tRNA pools to favor viral protein translation. Interestingly, specific host tRNAs are used as reverse transcription primers and are packaged into retroviral virions. Recent data also demonstrate the formation of tRNA-derived fragments (tRFs) upon infection to facilitate viral replication. Here, we comprehensively discuss how RNA viruses exploit distinct aspects of the host tRNA biology for their benefit. In light of the recent advances in the field, we propose that host tRNA-related pathways and mechanisms represent promising cellular targets for the development of novel antiviral strategies.

On the origin and continuing evolution of SARS-CoV-2

The SARS-CoV-2 epidemic started in late December 2019 in Wuhan, China, and has since impacted a large portion of China and raised major global concern. Herein, we investigated the extent of molecular divergence between SARS-CoV-2 and other related coronaviruses. Although we found only 4% variability in genomic nucleotides between SARS-CoV-2 and a bat SARS-related coronavirus (SARSr-CoV; RaTG13), the difference at neutral sites was 17%, suggesting the divergence between the two viruses is much larger than previously estimated. Our results suggest that the development of new variations in functional sites in the receptor-binding domain (RBD) of the spike seen in SARS-CoV-2 and viruses from pangolin SARSr-CoVs are likely caused by natural selection besides recombination. Population genetic analyses of 103 SARS-CoV-2 genomes indicated that these viruses had two major lineages (designated L and S), that are well defined by two different SNPs that show nearly complete linkage across the viral strains sequenced to date. We found that L lineage was more prevalent than the S lineage within the limited patient samples we examined. The implication of these evolutionary changes on disease etiology remains unclear. These findings strongly underscores the urgent need for further comprehensive studies that combine viral genomic data, with epidemiological studies of coronavirus disease 2019 (COVID-19).

On the origin and continuing evolution of SARS-CoV-2

The SARS-CoV-2 epidemic started in late December 2019 in Wuhan, China, and has since impacted a large portion of China and raised major global concern. Herein, we investigated the extent of molecular divergence between SARS-CoV-2 and other related coronaviruses. Although we found only 4% variability in genomic nucleotides between SARS-CoV-2 and a bat SARS-related coronavirus (SARSr-CoV; RaTG13), the difference at neutral sites was 17%, suggesting the divergence between the two viruses is much larger than previously estimated. Our results suggest that the development of new variations in functional sites in the receptor-binding domain (RBD) of the spike seen in SARS-CoV-2 and viruses from pangolin SARSr-CoVs are likely caused by natural selection besides recombination. Population genetic analyses of 103 SARS-CoV-2 genomes indicated that these viruses had two major lineages (designated L and S), that are well defined by two different SNPs that show nearly complete linkage across the viral strains sequenced to date. We found that L lineage was more prevalent than the S lineage within the limited patient samples we examined. The implication of these evolutionary changes on disease etiology remains unclear. These findings strongly underscores the urgent need for further comprehensive studies that combine viral genomic data, with epidemiological studies of coronavirus disease 2019 (COVID-19).

Emerging Roles of Translational Control in Circadian Timekeeping

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

Evidence of multifaceted functions of codon usage in translation within the model beetle Tribolium castaneum

Synonymous codon use is non-random. Codons most used in highly transcribed genes, often called optimal codons, typically have high gene counts of matching tRNA genes (tRNA abundance) and promote accurate and/or efficient translation. Non-optimal codons, those least used in highly expressed genes, may also affect translation. In multicellular organisms, codon optimality may vary among tissues. At present, however, tissue specificity of codon use remains poorly understood. Here, we studied codon usage of genes highly transcribed in germ line (testis and ovary) and somatic tissues (gonadectomized males and females) of the beetle Tribolium castaneum. The results demonstrate that: (i) the majority of optimal codons were organism-wide, the same in all tissues, and had numerous matching tRNA gene copies (Opt-codon_↑tRNAs), consistent with translational selection; (ii) some optimal codons varied among tissues, suggesting tissue-specific tRNA populations; (iii) wobble tRNA were required for translation of certain optimal codons (Opt-codon_wobble), possibly allowing precise translation and/or protein folding; and (iv) remarkably, some non-optimal codons had abundant tRNA genes (Nonopt-codon_↑tRNAs), and genes using those codons were tightly linked to ribosomal and stress-response functions. Thus, Nonopt-codon_↑tRNAs codons may regulate translation of specific genes. Together, the evidence suggests that codon use and tRNA genes regulate multiple translational processes in T. castaneum.

Impact of C-terminal amino acid composition on protein expression in bacteria

The C-terminal sequence of a protein is involved in processes such as efficiency of translation termination and protein degradation. However, the general relationship between features of this C-terminal sequence and levels of protein expression remains unknown. Here, we identified C-terminal amino acid biases that are ubiquitous across the bacterial taxonomy (1,582 genomes). We showed that the frequency is higher for positively charged amino acids (lysine, arginine), while hydrophobic amino acids and threonine are lower. We then studied the impact of C-terminal composition on protein levels in a library of Mycoplasma pneumoniae mutants, covering all possible combinations of the two last codons. We found that charged and polar residues, in particular lysine, led to higher expression, while hydrophobic and aromatic residues led to lower expression, with a difference in protein levels up to fourfold. We further showed that modulation of protein degradation rate could be one of the main mechanisms driving these differences. Our results demonstrate that the identity of the last amino acids has a strong influence on protein expression levels.

Comparative Analysis of the Complete Plastid Genome of Five Bupleurum Species and New Insights into DNA Barcoding and Phylogenetic Relationship

Bupleurum L. (Apiaceae) is a perennial and herbal genus, most species of which have high medicinal value. However, few studies have been performed using plastome data in this genus, and the phylogenetic relationships have always been controversial. In this study, the plastid genomes of Bupleurum chinense and Bupleurum commelynoideum were sequenced, and their gene content, order, and structure were counted and analyzed. The only three published Bupleurum species (B. boissieuanum, B. falcatum, and B. latissimum) and other fifteen allied species were selected to conduct a series of comparative and phylogenetic analyses. The genomes of B. chinense and B. commelynoideum were 155,869 and 155,629 bp in length, respectively, both of which had a typical quadripartite structure. The genome length, structure, guanine and cytosine (GC) content, and gene distribution were highly similar to the other three Bupleurum species. The five Bupleurum species had nearly the same codon usages, and eight regions (petN-psbM, rbcL-accD, ccsA-ndhD, trnK(UUU)-rps16, rpl32-trnL(UAG)-ccsA, petA-psbJ, ndhF-rpl32, and trnP(UGG)-psaJ-rpl33) were found to possess relatively higher nucleotide diversity, which may be the promising DNA barcodes in Bupleurum. Phylogenetic analysis revealed that all Bupleurum species clustered into a monophyletic clade with high bootstrap support and diverged after the Chamaesium clade. Overall, our study provides new insights into DNA barcoding and phylogenetic relationship between Bupleurum and its related genera, and will facilitate the population genomics, conservation genetics, and phylogenetics of Bupleurum in Apiaceae.

Codon usage bias affects α-amylase mRNA level by altering RNA stability and cytosine methylation patterns in Escherichia coli

Codon usage bias exists in almost every organism and is reported to regulate protein translation efficiency and folding. Besides translation, the preliminary role of codon usage bias on gene transcription has also been revealed in some eukaryotes such as Neurospora crassa. In this study, we took as an example the α-amylase-coding gene (amyA) and examined the role of codon usage bias in regulating gene expression in the typical prokaryote Escherichia coli. We confirmed the higher translation efficiency on codon-optimized amyA RNAs and found that the RNA level itself was also affected by codon optimization. The decreased RNA level was caused at least in part by altered mRNA stability at the post-transcriptional level. Codon optimization also altered the number of cytosine methylation sites. Examination on dcm knockouts suggested that cytosine methylation may be a minor mechanism adopted by codon bias to regulate gene RNA levels. More studies are required to verify the global effect of codon usage and to reveal its detailed mechanism on transcription.

Widespread non-modular overlapping codes in the coding regions

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

Analyses of genetics and pathogenesis of Salmonella enterica QH with narrow spectrum of antibiotic resistance isolated from yak

Salmonella is an important pathogen for public health due to food poisoning and acute infectious intestinal disease by zoonotic trait. We isolated Salmonella enterica QH which represents the normal growth condition in Luria-Bertani culture and displays a wide range of susceptibility for multiple antibiotics. To further investigate genetic and pathogenic traits of S. enterica QH, the sequencing genome of S. enterica QH and oral Salmonella infection in mice were performed in this study. Compared with other Salmonella strains, several large sequences containing prophages and genomic islands were inserted into S. enterica QH genome. Furthermore, nucleotide and synonymous codon usage patterns display mutation pressure and natural selection serving as drivers for the evolutionary trend of S. enterica QH at gene level. The unique codon usage pattern of S. enterica QH probably contributes to adaptation to environmental/host niches and to pathogenicity. In an early oral S. enterica QH infection, the levels of CD4⁺ and CD8⁺ lymphocytes significantly reduce in peripheral blood of mice, but the increasing transcription levels of some cytokines (IFN-β1, IFN-γ and CXCL10) might have pleiotypic immune effects against S. enterica QH infection. Of note, IL10 displays significant enhancement at levels of transcription and translation, suggesting that immunosuppressive effects mediated by IL10 may function as an early oral S. enterica QH infection. The systemic investigations, including genomic and genetic characterizations and biological traits of S. enterica QH in vivo and in vitro may reflect the basic lifestyle of S. enterica QH, requiring intestine colonization, undergoing environmental stresses and performing dissemination.

Ribosome pausing, a dangerous necessity for co-translational events

In recent years translation elongation has emerged as an important contributor to the regulation of gene expression. There are multiple quality control checkpoints along the way of producing mature proteins and targeting them to the right cellular compartment, or associating them correctly with their partners. Ribosomes pause to allow co-translational protein folding, protein targeting or protein interactions, and the pausing is dictated by a combination of the mRNA sequence and structure, the tRNA availability and the nascent peptide. However, ribosome pausing can also lead to ribosome collisions and co-translational degradation of both mRNA and nascent chain. Understanding how the translating ribosome tunes the different maturation steps that nascent proteins must undergo, what the timing of these maturation events is, and how degradation can be avoided when pausing is needed, is now possible by the emergence of methods to follow ribosome dynamics in vivo. This review summarizes some of the recent studies that have advanced our knowledge about co-translational events using the power of ribosome profiling, and some of the questions that have emerged from these studies.

The human methyltransferase ZCCHC4 catalyses N6-methyladenosine modification of 28S ribosomal RNA

RNA methylations are essential both for RNA structure and function, and are introduced by a number of distinct methyltransferases (MTases). In recent years, N⁶-methyladenosine (m⁶A) modification of eukaryotic mRNA has been subject to intense studies, and it has been demonstrated that m⁶A is a reversible modification that regulates several aspects of mRNA function. However, m⁶A is also found in other RNAs, such as mammalian 18S and 28S ribosomal RNAs (rRNAs), but the responsible MTases have remained elusive. 28S rRNA carries a single m⁶A modification, found at position A4220 (alternatively referred to as A4190) within a stem–loop structure, and here we show that the MTase ZCCHC4 is the enzyme responsible for introducing this modification. Accordingly, we found that ZCCHC4 localises to nucleoli, the site of ribosome assembly, and that proteins involved in RNA metabolism are overrepresented in the ZCCHC4 interactome. Interestingly, the absence of m⁶A4220 perturbs codon-specific translation dynamics and shifts gene expression at the translational level. In summary, we establish ZCCHC4 as the enzyme responsible for m⁶A modification of human 28S rRNA, and demonstrate its functional significance in mRNA translation.

Significance and roles of synonymous codon usage in the evolutionary process of Proteus

Proteus spp. bacteria frequently serve as opportunistic pathogens that can infect many animals and show positive survival and existence in various natural environments. The evolutionary pattern of Proteus spp. is an unknown topic, which benefits understanding the different evolutionary dynamics for excellent bacterial adaptation to various environments. Here, the eight whole genomes of different Proteus species were analyzed for the interplay between nucleotide usage and synonymous codon usage. Although the orthologous average nucleotide identity and average nucleotide identity display the genetic diversity of these Proteus species at the genome level, the principal component analysis further shows that these species sustain the specific genetic niche at the aspect of synonymous codon usage patterns. Interestingly, although these Proteus species have A/T rich genes with underrepresented G (guanine) or C (cytosine) at the third codon positions and overrepresented A or T at these positions, some synonymous codons with A or T end are obviously suppressed in usage. The overall codon usage pattern reflected by the effective number of codons (ENC) has a significantly positive correlation with GC3 content (GC content at the third codon position), and ENC has a significantly negative correlation with the adaptation index for these species. These results suggest that the mutation pressure caused by nucleotide composition constraint serves as a dominant evolutionary dynamic driving evolutionary trend of Proteus spp., along with other selections related to natural selection, replication and fine-tune translation, and so on. Taken together, the analyses help to understand the evolutionary interplay between nucleotide and codon usage at the gene level of Proteus.

Mass-spectrometry-based draft of the Arabidopsis proteome

Plants are essential for life and are extremely diverse organisms with unique molecular capabilities¹. Here we present a quantitative atlas of the transcriptomes, proteomes and phosphoproteomes of 30 tissues of the model plant Arabidopsis thaliana. Our analysis provides initial answers to how many genes exist as proteins (more than 18,000), where they are expressed, in which approximate quantities (a dynamic range of more than six orders of magnitude) and to what extent they are phosphorylated (over 43,000 sites). We present examples of how the data may be used, such as to discover proteins that are translated from short open-reading frames, to uncover sequence motifs that are involved in the regulation of protein production, and to identify tissue-specific protein complexes or phosphorylation-mediated signalling events. Interactive access to this resource for the plant community is provided by the ProteomicsDB and ATHENA databases, which include powerful bioinformatics tools to explore and characterize Arabidopsis proteins, their modifications and interactions.

Bacterial RNA in extracellular vesicles: A new regulator of host-pathogen interactions?

Extracellular vesicles (EVs) are released by cells from all kingdoms and represent one form of cell-cell interaction. This universal system of communication blurs cells type boundaries, offering an new avenue for pathogens to infect their hosts. EVs carry with them an arsenal of virulence factors that have been the focus of numerous studies. During the last years, the RNA content of EVs has also gained increasing attention, particularly in the context of infection. Secreted RNA in EVs (evRNA) from several bacterial pathogens have been characterised but the exact mechanisms promoting pathogenicity remain elusive. In this review, we evaluate the potential of such transcripts to directly interact with targets in infected cells and, by this, represent a novel angle of host cell manipulation during bacterial infection. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.

The Genetic Landscape of Dystrophin Mutations in Italy: A Nationwide Study

Dystrophinopathies are inherited diseases caused by mutations in the dystrophin (DMD) gene for which testing is mandatory for genetic diagnosis, reproductive choices and eligibility for personalized trials. We genotyped the DMD gene in our Italian cohort of 1902 patients (BMD n = 740, 39%; DMD n =1162, 61%) within a nationwide study involving 11 diagnostic centers in a 10-year window (2008–2017). In DMD patients, we found deletions in 57%, duplications in 11% and small mutations in 32%. In BMD, we found deletions in 78%, duplications in 9% and small mutations in 13%. In BMD, there are a higher number of deletions, and small mutations are more frequent than duplications. Among small mutations that are generally frequent in both phenotypes, 44% of DMD and 36% of BMD are nonsense, thus, eligible for stop codon read-through therapy; 63% of all out-of-frame deletions are eligible for single exon skipping. Patients were also assigned to Italian regions and showed interesting regional differences in mutation distribution. The full genetic characterization in this large, nationwide cohort has allowed us to draw several correlations between DMD/BMD genotype landscapes and mutation frequency, mutation types, mutation locations along the gene, exon/intron architecture, and relevant protein domain, with effects on population genetic characteristics and new personalized therapies.

Tailoring mRNA Vaccine to Balance Innate/Adaptive Immune Response

mRNA vaccine platforms present numerous advantages, such as versatility, rapid production, and induction of cellular and humoral responses. Moreover, mRNAs have inherent adjuvant properties due to their complex interaction with pattern recognition receptors (PRRs). This recognition can be either beneficial in activating antigen-presenting cells (APCs) or detrimental by indirectly blocking mRNA translation. To decipher this Janus effect, we describe the different innate response mechanisms triggered by mRNA molecules and how each element from the 5' cap to the poly-A tail interferes with innate/adaptive immune responses. Then, we emphasize the importance of some critical steps such as production, purification, and formulation as key events to further improve the quality of immune responses and balance innate and adaptive immunity.

Translational efficiency across healthy and tumor tissues is proliferation-related

Different tissues express genes with particular codon usage and anticodon tRNA repertoires. However, the codon-anticodon co-adaptation in humans is not completely understood, nor is its effect on tissue-specific protein levels. Here, we first validated the accuracy of small RNA-seq for tRNA quantification across five human cell lines. We then analyzed the tRNA abundance of more than 8,000 tumor samples from TCGA, together with their paired mRNA-seq and proteomics data, to determine the Supply-to-Demand Adaptation. We thereby elucidate that the dynamic adaptation of the tRNA pool is largely related to the proliferative state across tissues. The distribution of such tRNA pools over the whole cellular translatome affects the subsequent translational efficiency, which functionally determines a condition-specific expression program both in healthy and tumor states. Furthermore, the aberrant translational efficiency of some codons in cancer, exemplified by ProCCA and GlyGGT, is associated with poor patient survival. The regulation of these tRNA profiles is partly explained by the tRNA gene copy numbers and their promoter DNA methylation.

Codon and amino acid content are associated with mRNA stability in mammalian cells

Messenger RNA (mRNA) degradation plays a critical role in regulating transcript levels in the cell and is a major control point for modulating gene expression. In yeast and other model organisms, codon identity is a powerful determinant of transcript stability, contributing broadly to impact half-lives. General principles governing mRNA stability are poorly understood in mammalian systems. Importantly, however, the degradation machinery is highly conserved, thus it seems logical that mammalian transcript half-lives would also be strongly influenced by coding determinants. Herein we characterize the contribution of coding sequence towards mRNA decay in human and Chinese Hamster Ovary cells. In agreement with previous studies, we observed that synonymous codon usage impacts mRNA stability in mammalian cells. Surprisingly, however, we also observe that the amino acid content of a gene is an additional determinant correlating with transcript stability. The impact of codon and amino acid identity on mRNA decay appears to be associated with underlying tRNA and intracellular amino acid concentrations. Accordingly, genes of similar physiological function appear to coordinate their mRNA stabilities in part through codon and amino acid content. Together, these results raise the possibility that intracellular tRNA and amino acid levels interplay to mediate coupling between translational elongation and mRNA degradation rate in mammals.

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

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