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

The chloroplast genome of Chrozophora sabulosa Kar. & Kir. and its exploration in the evolutionary position uncertainty of genus Chrozophora

Chrozophora sabulosa Kar. & Kir. is a biennial herbaceous plant that belongs to the Euphorbiaceae family and has medicinal properties. This research aimed to identify the genetic characteristics and phylogenetic position of the Chrozophora genus within the Euphorbiaceae family. The evolutionary position of the Chrozophora genus was previously unknown due to insufficient research. Therefore, to determine the evolutionary link between C. sabulosa and other related species, we conducted a study using the NGS Illumina platform to sequence the C. sabulosa chloroplast (cp.) genome. The study results showed that the genome was 156,488 bp in length. It had a quadripartite structure consisting of two inverted repeats (IRb and IRa) of 24,649-bp, separated by an 87,696-bp LSC region and a 19,494-bp SSC region. The CP genome contained 113 unique genes, including four rRNA genes, 30 tRNA genes, and 79 CDS genes. In the second copy of the inverted repeat, there were 18 duplicated genes. The C. sabulosa lacks the petD, petB, rpl2, and rps16 intron. The analysis of simple sequence repeats (SSRs) revealed 93 SSR loci of 22 types and 78 oligonucleotide repeats of four kinds. The phylogenetic investigation showed that the Chrozophora genus evolved paraphyletically from other members of the Euphorbiaceae family. To support the phylogenetic findings, we selected species from the Euphorbiaceae and Phyllanthaceae families to compare with C. sabulosa for Ks and Ka substitution rates, InDels investigation, IR contraction and expansion, and SNPs analysis. The results of these comparative studies align with the phylogenetic findings. We identified six highly polymorphic regions shared by both families, which could be used as molecular identifiers for the Chrozophora genus (rpl33-rps18, rps18-rpl20, rps15-ycf1, ndhG-ndhI, psaI-ycf4, petA-psbJ). The cp. genome sequence of C. sabulosa reveals the evolution of plastid sequences in Chrozophora species. This is the first time the cp. genome of a Chrozophora genus has been sequenced, serving as a foundation for future sequencing of other species within the Chrozophoreae tribe and facilitating in-depth taxonomic research. The results of this research will also aid in identifying new Chrozophora species.

Genome-wide sequencing reveals geographical variations in the thermal adaptation of an aquaculture species with frequent seedling introductions

Climate change and human activity are essential factors affecting marine biodiversity and aquaculture, and understanding the impacts of human activities on the genetic structure to increasing high temperatures is crucial for sustainable aquaculture and marine biodiversity conservation. As a commercially important bivalve, the Manila clam Ruditapes philippinarum is widely distributed along the coast of China, and it has been frequently introduced from Fujian Province, China, to other regions for aquaculture. In this study, we collected four populations of Manila clams from different areas to evaluate their thermal tolerance by measuring cardiac performance and genetic variations using whole-genome resequencing. The upper thermal limits of the clams showed high variations within and among populations. Different populations displayed divergent genetic compositions, and the admixed population was partly derived from the Zhangzhou population in Fujian Province, implying a complex genomic landscape under the influence of local genetic sources and human introductions. Multiple single nucleotide polymorphisms (SNPs) were associated with the cardiac functional traits, and some of these SNPs can affect the codon usage and the structural stability of the resulting protein. This study shed light on the importance of establishing long-term ecological and genetic monitoring programs at the local level to enhance resilience to future climate change.

Navigating the complexities of multi-domain protein folding

Proteome complexity has expanded tremendously over evolutionary time, enabling biological diversification. Much of this complexity is achieved by combining a limited set of structural units into long polypeptides. This widely used evolutionary strategy poses challenges for folding of the resulting multi-domain proteins. As a consequence, their folding differs from that of small single-domain proteins, which generally fold quickly and reversibly. Co-translational processes and chaperone interactions are important aspects of multi-domain protein folding. In this review, we discuss some of the recent experimental progress toward understanding these processes.

Genomic factors shaping codon usage across the Saccharomycotina subphylum

Codon usage bias, or the unequal use of synonymous codons, is observed across genes, genomes, and between species. The biased use of synonymous codons has been implicated in many cellular functions, such as translation dynamics and transcript stability, but can also be shaped by neutral forces. The Saccharomycotina, the fungal subphylum containing the yeasts Saccharomyces cerevisiae and Candida albicans , has been a model system for studying codon usage. We characterized codon usage across 1,154 strains from 1,051 species to gain insight into the biases, molecular mechanisms, evolution, and genomic features contributing to codon usage patterns across the subphylum. We found evidence of a general preference for A/T-ending codons and correlations between codon usage bias, GC content, and tRNA-ome size. Codon usage bias is also distinct between the 12 orders within the subphylum to such a degree that yeasts can be classified into orders with an accuracy greater than 90% using a machine learning algorithm trained on codon usage. We also characterized the degree to which codon usage bias is impacted by translational selection. Interestingly, the degree of translational selection was influenced by a combination of genome features and assembly metrics that included the number of coding sequences, BUSCO count, and genome length. Our analysis also revealed an extreme bias in codon usage in the Saccharomycodales associated with a lack of predicted arginine tRNAs. The order contains 24 species, and 23 are computationally predicted to lack tRNAs that decode CGN codons, leaving only the AGN codons to encode arginine. Analysis of Saccharomycodales gene expression, tRNA sequences, and codon evolution suggests that extreme avoidance of the CGN codons is associated with a decline in arginine tRNA function. Codon usage bias within the Saccharomycotina is generally consistent with previous investigations in fungi, which show a role for both genomic features and GC bias in shaping codon usage. However, we find cases of extreme codon usage preference and avoidance along yeast lineages, suggesting additional forces may be shaping the evolution of specific codons.

Molecular Evolution of Protein Sequences and Codon Usage in Monkeypox Viruses

The monkeypox virus (mpox virus, MPXV) epidemic in 2022 has posed a significant public health risk. Yet, the evolutionary principles of MPXV remain largely unknown. Here, we examined the evolutionary patterns of protein sequences and codon usage in MPXV. We first demonstrated the signal of positive selection in OPG027, specifically in the Clade I lineage of MPXV. Subsequently, we discovered accelerated protein sequence evolution over time in the variants responsible for the 2022 outbreak. Furthermore, we showed strong epistasis between amino acid substitutions located in different genes. The codon adaptation index (CAI) analysis revealed that MPXV genes tended to use more non-preferred codons compared to human genes, and the CAI decreased over time and diverged between clades, with Clade I > IIa and IIb-A > IIb-B. While the decrease in fatality rate among the three groups aligned with the CAI pattern, it remains unclear whether this correlation was coincidental or if the deoptimization of codon usage in MPXV led to a reduction in fatality rates. This study sheds new light on the mechanisms that govern the evolution of MPXV in human populations.

Quantifying negative selection on synonymous variants

Widespread adoption of DNA sequencing has resulted in large numbers of genetic variants, whose contribution to disease is not easily determined. Although many types of variation are known to disrupt cellular processes in predictable ways, for some categories of variants, the effects may not be directly detectable. A particular example is synonymous variants, that is, those single-nucleotide variants that create a codon substitution, such that the produced amino acid sequence is unaffected. Contrary to the original theory suggesting that synonymous variants are benign, there is a growing volume of research showing that, despite their "silent" mechanism of action, some synonymous variation may be deleterious. Here, we studied the extent of the negative selective pressure acting on different classes of synonymous variants by analyzing the relative enrichment of synonymous singleton variants in the human exomes provided by gnomAD. Using a modification of the mutability-adjusted proportion of singletons (MAPS) metric as a measure of purifying selection, we found that some classes of synonymous variants are subject to stronger negative selection than others. For instance, variants that reduce codon optimality undergo stronger selection than optimality-increasing variants. Besides, selection affects synonymous variants implicated in splice-site-loss or splice-site-gain events. To understand what drives this negative selection, we tested a number of predictors in the aim to explain the variability in the selection scores. Our findings provide insights into the effects of synonymous variants at the population level, highlighting the specifics of the role that these variants play in health and disease.

Contiguous identity between entire coding regions of transgenic and native genes rather than special regions is essential for a strong co-suppression

The discovery of co-suppression in plants has greatly boosted the study of gene silencing mechanisms, but its triggering mechanism has remained a mystery. In this study, we explored its possible trigger mechanism by using Fatty acid desaturase 2 (FAD2) and Fatty acid elongase 1 (FAE1) strong co-suppression systems. Analysis of small RNAs in FAD2 co-suppression lines showed that siRNAs distributed throughout the coding region of FAD2 with an accumulated peak. However, mutations of the peak siRNA-matched site and siRNA derived site had not alleviated the co-suppression of its transgenic lines. Synthetic FAD2 (AtFAD2sm), which has synonymous mutations in the entire coding region, failed to trigger any co-suppression. Furthermore, 5' and 3' portions of AtFAD2 and AtFAD2sm were swapped to form two hybrid genes, AtFAD2-3sm and AtFAD2-5sm. 80 % and 92 % of their transgenic lines exhibited co-suppression, respectively. Finally, FAE1s with different degrees of the continuous sequence identity compared with AtFAE1 were tested in their Arabidopsis transgenic lines, and the results showed the co-suppression frequency was reduced as their continuous sequence identity stepped down. This work suggests that contiguous identity between the entire coding regions of transgenic and native genes rather than a special region is essential for a strong co-suppression.

Specific codons control cellular resources and fitness

As cellular engineering progresses from simply overexpressing proteins to imparting complex phenotypes through multigene expression, judicious appropriation of cellular resources is essential. Since codon use is degenerate and biased, codons may control cellular resources at a translational level. We investigate how partitioning transfer RNA (tRNA) resources by incorporating dissimilar codon usage can drastically alter interdependence of expression level and burden on the host. By isolating the effect of individual codons' use during translation elongation while eliminating confounding factors, we show that codon choice can trans-regulate fitness of the host and expression of other heterologous or native genes. We correlate specific codon usage patterns with host fitness and derive a coding scheme for multigene expression called the Codon Health Index (CHI, χ). This empirically derived coding scheme (χ) enables the design of multigene expression systems that avoid catastrophic cellular burden and is robust across several proteins and conditions.

Fitness benefits of a synonymous substitution in an ancient EF-Tu gene depend on the genetic background

Synonymous mutations are changes to DNA sequence, which occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation elongation factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels as well as global polysome abundance on RNA transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.IMPORTANCEThis study explores the degree to which synonymous mutations in essential genes can influence adaptation in bacteria. An experimental system whereby an Escherichia coli strain harboring an engineered translation protein elongation factor-Tu (EF-Tu) was subjected to laboratory evolution. We find that a synonymous mutation acquired on the gene encoding for EF-Tu is conditionally beneficial for bacterial fitness. Our findings provide insight into the importance of the genetic background when a synonymous substitution is favored by natural selection and how such changes have the potential to impact evolution when critical cellular processes are involved.

DNA-Based Nonviral Gene Therapy─Challenging but Promising

Over the past decades, significant progress has been made in utilizing nucleic acids, including DNA and RNA molecules, for therapeutic purposes. For DNA molecules, although various DNA delivery systems have been established, viral vector systems are the go-to choice for large-scale commercial applications. However, viral systems have certain disadvantages such as immune response, limited payload capacity, insertional mutagenesis and pre-existing immunity. In contrast, nonviral systems are less immunogenic, not size limited, safer, and easier for manufacturing compared with viral systems. What's more, nonviral DNA vectors have demonstrated their capacity to mediate specific protein expression in vivo for diverse therapeutic objectives containing a wide range of diseases such as cancer, rare diseases, neurodegenerative diseases, and infectious diseases, yielding promising therapeutic outcomes. However, exogenous plasmid DNA is prone to degrade and has poor immunogenicity in vivo. Thus, various strategies have been developed: (i) designing novel plasmids with special structures, (ii) optimizing plasmid sequences for higher expression, and (iii) developing more efficient nonviral DNA delivery systems. Based on these strategies, many interesting clinical results have been reported. This Review discusses the development of DNA-based nonviral gene therapy, including novel plasmids, nonviral delivery systems, clinical advances, and prospects. These developments hold great potential for enhancing the efficacy and safety of nonviral gene therapy and expanding its applications in the treatment of various diseases.

Proteomics revealed an association between ribosome-associated proteins and amyloid beta deposition in Alzheimer's disease

Most scholars believe that amyloid-beta (Aβ) has the potential to induce apoptosis, stimulate an inflammatory cascade, promote oxidative stress and exacerbate the pathological progression of Alzheimer's disease (AD). Therefore, it is crucial to investigate the deposition of Aβ in AD. At approximately 6 months of age, APP/PS1 double transgenic mice gradually exhibit the development of plaques, as well as spatial and learning impairment. Notably, the hippocampus is specifically affected in the course of AD. Herein, 6-month-old APP/PS1 double transgenic mice were utilized, and the differentially expressed (DE) proteins in the hippocampus were identified and analyzed using 4D label-free quantitative proteomics technology and parallel reaction monitoring (PRM). Compared to wild-type mice, 29 proteins were upregulated and 25 proteins were downregulated in the AD group. Gene Ontology (GO) enrichment analysis of biological processes (BP) indicated that the DE proteins were mainly involved in 'ribosomal large subunit biogenesis'. Molecular function (MF) analysis results were primarily associated with '5.8S rRNA binding' and 'structural constituent of ribosome'. In terms of cellular components (CC), the DE proteins were mainly found in 'polysomal ribosome', 'cytosolic large ribosomal subunit', 'cytosolic ribosome', and 'large ribosomal subunit', among others. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis demonstrated that the results were mainly enriched in the 'Ribosome signaling pathway'. The key target proteins identified were ribosomal protein (Rp)l18, Rpl17, Rpl19, Rpl24, Rpl35, and Rpl6. The PRM verification results were consistent with the findings of the 4D label-free quantitative proteomics analysis. Overall, these findings suggest that Rpl18, Rpl17, Rpl19, Rpl24, Rpl35, and Rpl6 may have potential therapeutic value for the treatment of AD by targeting Aβ.

Ribosome inactivation regulates translation elongation in neurons

Cellular plasticity is crucial for adapting to ever-changing stimuli. As a result, cells consistently reshape their translatome, and, consequently, their proteome. The control of translational activity has been thoroughly examined at the stage of translation initiation. However, the regulation of ribosome speed in cells is widely unknown. In this study, we utilized a timed ribosome runoff approach, along with proteomics and transmission electron microscopy, to investigate global translation kinetics in cells. We found that ribosome speeds vary among various cell types, such as astrocytes, induced pluripotent human stem cells, human neural stem cells, and human and rat neurons. Of all cell types studied, mature cortical neurons exhibit the highest rate of translation. This finding is particularly remarkable because mature cortical neurons express the eukaryotic elongation factor 2 (eEF2) at lower levels than other cell types. Neurons solve this conundrum by inactivating a fraction of their ribosomes. As a result, the increase in eEF2 levels leads to a reduction of inactive ribosomes and an enhancement of active ones. Processes that alter the demand for active ribosomes, like neuronal excitation, cause increased inactivation of redundant ribosomes in an eEF2-dependent manner. Our data suggest a novel regulatory mechanism in which neurons dynamically inactivate ribosomes to facilitate translational remodeling. These findings have important implications for developmental brain disorders characterized by, among other things, aberrant translation.

Synonymous edits in the Escherichia coli genome have substantial and condition-dependent effects on fitness

CRISPR-Cas-based genome editing is widely used in bacteria at scales ranging from construction of individual mutants to massively parallel libraries. This procedure relies on guide RNA-directed cleavage of the genome followed by repair with a template that introduces a desired mutation along with synonymous "immunizing" mutations to prevent re-cleavage of the genome after editing. Because the immunizing mutations do not change the protein sequence, they are often assumed to be neutral. However, synonymous mutations can change mRNA structures in ways that alter levels of the encoded proteins. We have tested the assumption that immunizing mutations are neutral by constructing a library of over 50,000 edits that consist of only synonymous mutations in Escherichia coli. Thousands of edits had substantial effects on fitness during growth of E. coli on acetate, a poor carbon source that is toxic at high concentrations. The percentage of high-impact edits varied considerably between genes and at different positions within genes. We reconstructed clones with high-impact edits and found that 69% indeed had significant effects on growth in acetate. Interestingly, fewer edits affected fitness during growth in glucose, a preferred carbon source, suggesting that changes in protein expression caused by synonymous mutations may be most important when an organism encounters challenging conditions. Finally, we showed that synonymous edits can have widespread effects; a synonymous edit at the 5' end of ptsI altered expression of hundreds of genes. Our results suggest that the synonymous immunizing edits introduced during CRISPR-Cas-based genome editing should not be assumed to be innocuous.

Ribosome profiling: a powerful tool in oncological research

Neoplastic cells need to adapt their gene expression pattern to survive in an ever-changing or unfavorable tumor microenvironment. Protein synthesis (or mRNA translation), an essential part of gene expression, is dysregulated in cancer. The emergence of distinct translatomic technologies has revolutionized oncological studies to elucidate translational regulatory mechanisms. Ribosome profiling can provide adequate information on diverse aspects of translation by aiding in quantitatively analyzing the intensity of translating ribosome-protected fragments. Here, we review the primary currently used translatomics techniques and highlight their advantages and disadvantages as tools for translatomics studies. Subsequently, we clarified the areas in which ribosome profiling could be applied to better understand translational control. Finally, we summarized the latest advances in cancer studies using ribosome profiling to highlight the extensive application of this powerful and promising translatomic tool.

An Improved Method for Eliminating or Creating Intragenic Bacterial Promoters

Recent advances in genomic refactoring have been hindered by the ever-present complication of internal or cryptic transcriptional regulation. Typical approaches to these features have been to randomize or perform mass alterations to the gene sequences thought to contain the regulatory motifs; however, this approach can cause problems by altering translational speeds, introducing long distance DNA-DNA interaction effects, and inducing RNA toxicity. Previously, we developed a rational design approach named COdon Restrained Promoter SilEncing (CORPSE) which takes externally identified promoter sequences and uses position-specific scoring matrices as proxy promoter strengths to make minimal changes to promoter sequences to disable their activity. Additionally, through inverting our system we were also able to modify weak internal promoters to increase their activity. In this chapter, we augment our previous process with the biophysical model Promoter Calculator v1.0 developed by LaFleur et al. to combine promoter identification and activity prediction, with our algorithm to silently modify promoter sequences, to provide more robust promoter elimination and creation.

Synonymous codon usage regulates translation initiation

Nonoptimal synonymous codons repress gene expression, but the underlying mechanisms are poorly understood. We and others have previously shown that nonoptimal codons slow translation elongation speeds and thereby trigger messenger RNA (mRNA) degradation. Nevertheless, transcript levels are often insufficient to explain protein levels, suggesting additional mechanisms by which codon usage regulates gene expression. Using reporters in human and Drosophila cells, we find that transcript levels account for less than half of the variation in protein abundance due to codon usage. This discrepancy is explained by translational differences whereby nonoptimal codons repress translation initiation. Nonoptimal transcripts are also less bound by the translation initiation factors eIF4E and eIF4G1, providing a mechanistic explanation for their reduced initiation rates. Importantly, translational repression can occur without mRNA decay and deadenylation, and it does not depend on the known nonoptimality sensor, CNOT3. Our results reveal a potent mechanism of regulation by codon usage where nonoptimal codons repress further rounds of translation.

Translation Rates and Protein Folding

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

A small stretch of poor codon usage at the beginning of dengue virus open reading frame may act as a translational checkpoint

OBJECTIVE

Rare codons were previously shown to be enriched at the beginning of the dengue virus (DENV) open reading frame. However, the role of rare codons in regulating translation efficiency and replication of DENV remains unclear. The present study aims to clarify the significance of rare codon usage at the beginning of DENV transcripts using the codon adaptation index (CAI).

METHODOLOGY

CAIs of the whole starting regions of DENV transcripts as well as 18-codon sliding windows of the regions were analyzed.

RESULTS

One of the intriguing findings is that those rare codons do not typically result in uniformly low CAI in the starting region with rare codons. However, it shows a notable local drop in CAI around the 50th codon in all dengue serotypes. This suggests that there may be a translational checkpoint at this site and that the rare codon usage upstream to this checkpoint may not be related to translational control.

Quantification of elongation stalls and impact on gene expression in yeast

Ribosomal pauses are a critical part of cotranslational events including protein folding and localization. However, extended ribosome pauses can lead to ribosome collisions, resulting in the activation of ribosome rescue pathways and turnover of protein and mRNA. While this relationship has been known, there has been little exploration of how ribosomal stalls impact translation duration at a quantitative level. We have taken a method used to measure elongation time and adapted it for use in Saccharomyces cerevisiae to quantify the impact of elongation stalls. We find, in transcripts containing Arg CGA codon repeat-induced stalls, a Hel2-mediated dose-dependent decrease in protein expression and mRNA level and an elongation delay on the order of minutes. In transcripts that contain synonymous substitutions to nonoptimal Leu codons, there is a decrease in protein and mRNA levels, as well as similar elongation delay, but this occurs through a non-Hel2-mediated mechanism. Finally, we find that Dhh1 selectively increases protein expression, mRNA level, and elongation rate. This indicates that distinct poorly translated mRNAs will activate different rescue pathways despite similar elongation stall durations. Taken together, these results provide new quantitative mechanistic insight into the surveillance of translation and the roles of Hel2 and Dhh1 in mediating ribosome pausing events.

Establishing chromosomal design-build-test-learn through a synthetic chromosome and its combinatorial reconfiguration

Chromosome-level design-build-test-learn cycles (chrDBTLs) allow systematic combinatorial reconfiguration of chromosomes with ease. Here, we established chrDBTL with a redesigned synthetic Saccharomyces cerevisiae chromosome XV, synXV. We designed and built synXV to harbor strategically inserted features, modified elements, and synonymously recoded genes throughout the chromosome. Based on the recoded chromosome, we developed a method to enable chrDBTL: CRISPR-Cas9-mediated mitotic recombination with endoreduplication (CRIMiRE). CRIMiRE allowed the creation of customized wild-type/synthetic combinations, accelerating genotype-phenotype mapping and synthetic chromosome redesign. We also leveraged synXV as a "build-to-learn" model organism for translation studies by ribosome profiling. We conducted a locus-to-locus comparison of ribosome occupancy between synXV and the wild-type chromosome, providing insight into the effects of codon changes and redesigned features on translation dynamics in vivo. Overall, we established synXV as a versatile reconfigurable system that advances chrDBTL for understanding biological mechanisms and engineering strains.

Synthetic yeast chromosome XI design provides a testbed for the study of extrachromosomal circular DNA dynamics

We describe construction of the synthetic yeast chromosome XI (synXI) and reveal the effects of redesign at non-coding DNA elements. The 660-kb synthetic yeast genome project (Sc2.0) chromosome was assembled from synthesized DNA fragments before CRISPR-based methods were used in a process of bug discovery, redesign, and chromosome repair, including precise compaction of 200 kb of repeat sequence. Repaired defects were related to poor centromere function and mitochondrial health and were associated with modifications to non-coding regions. As part of the Sc2.0 design, loxPsym sequences for Cre-mediated recombination are inserted between most genes. Using the GAP1 locus from chromosome XI, we show that these sites can facilitate induced extrachromosomal circular DNA (eccDNA) formation, allowing direct study of the effects and propagation of these important molecules. Construction and characterization of synXI contributes to our understanding of non-coding DNA elements, provides a useful tool for eccDNA study, and will inform future synthetic genome design.

Combination of hydrophobicity and codon usage bias determines sorting of model K+ channel protein to either mitochondria or endoplasmic reticulum

When the K+ channel-like protein Kesv from Ectocarpus siliculosus virus 1 is heterologously expressed in mammalian cells, it is sorted to the mitochondria. This targeting can be redirected to the endoplasmic reticulum (ER) by altering the codon usage in distinct regions of the gene or by inserting a triplet of hydrophobic amino acids (AAs) into the protein's C-terminal transmembrane domain (ct-TMD). Systematic variations in the flavor of the inserted AAs and/or its codon usage show that a positive charge in the inserted AA triplet alone serves as strong signal for mitochondria sorting. In cases of neutral AA triplets, mitochondria sorting are favored by a combination of hydrophilic AAs and rarely used codons; sorting to the ER exhibits the inverse dependency. This propensity for ER sorting is particularly high when a common codon follows a rarer one in the AA triplet; mitochondria sorting in contrast is supported by codon uniformity. Since parameters like positive charge, hydrophobic AAs, and common codons are known to facilitate elongation of nascent proteins in the ribosome the data suggest a mechanism in which local changes in elongation velocity and co-translational folding in the ct-TMD influence intracellular protein sorting.

VeloPro: A pipeline integrating Ribo-seq and AlphaFold deciphers association patterns between translation velocity and protein structure features

VeloPro integrates Ribo-seq data and AlphaFold2-predicted 3D protein structure information for characterization of the association patterns between translation velocity and many protein structure features in prokaryotic and eukaryotic organisms across different taxonomical clades such as bacteria, fungi, protozoa, nematode, plants, insect, and mammals. We illustrated that association patterns between translation velocity and protein structure features differ across organisms, partially reflecting their taxonomical relationship.

Kinetics of Translating Ribosomes Determine the Efficiency of Programmed Stop Codon Readthrough

During translation, a stop codon on the mRNA signals the ribosomes to terminate the process. In certain mRNAs, the termination fails due to the recoding of the canonical stop codon, and ribosomes continue translation to generate C-terminally extended protein. This process, termed stop codon readthrough (SCR), regulates several cellular functions. SCR is driven by elements/factors that act immediately downstream of the stop codon. Here, we have analysed the process of SCR using a simple mathematical model to investigate how the kinetics of translating ribosomes influences the efficiency of SCR. Surprisingly, the analysis revealed that the rate of translation inversely regulates the efficiency of SCR. We tested this prediction experimentally in mammalian AGO1 and MTCH2 mRNAs. Reduction in translation either globally by harringtonine or locally by rare codons caused an increase in the efficiency of SCR. Thus, our study has revealed a hitherto unknown mode of regulation of SCR.

Proteostasis regulation through ribosome quality control and no-go-decay

Cell functionality relies on the existing pool of proteins and their folding into functional conformations. This is achieved through the regulation of protein synthesis, which requires error-free mRNAs and ribosomes. Ribosomes are quality control hubs for mRNAs and proteins. Problems during translation elongation slow down the decoding rate, leading to ribosome halting and the eventual collision with the next ribosome. Collided ribosomes form a specific disome structure recognized and solved by ribosome quality control (RQC) mechanisms. RQC pathways orchestrate the degradation of the problematic mRNA by no-go decay and the truncated nascent peptide, the repression of translation initiation, and the recycling of the stalled ribosomes. All these events maintain protein homeostasis and return valuable ribosomes to translation. As such, cell homeostasis and function are maintained at the mRNA level by preventing the production of aberrant or unnecessary proteins. It is becoming evident that the crosstalk between RQC and the protein homeostasis network is vital for cell function, as the absence of RQC components leads to the activation of stress response and neurodegenerative diseases. Here, we review the molecular events of RQC discovered through well-designed stalling reporters. Given the impact of RQC in proteostasis, we discuss the relevance of identifying endogenous mRNA regulated by RQC and their preservation in stress conditions. This article is categorized under: RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms Translation > Regulation.

Genes for highly abundant proteins in Escherichia coli avoid 5' codons that promote ribosomal initiation

In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria.

Comparative analysis of codon usage patterns and phylogenetic implications of five mitochondrial genomes of the genus Japanagallia Ishihara, 1955 (Hemiptera, Cicadellidae, Megophthalminae)

Japanagallia is a genus of Cicadomorpha in the family of leafhoppers that are plant piercing-sucking insects, and it is difficult to distinguish by morphological characteristics. So far, only one complete mitochondrial genome data has been reported for the genus Japanagallia. Therefore, in order to better understand this group, we assembled and annotated the complete mitochondrial genomes of five Japanagallia species, and analyzed their codon usage patterns. Nucleotide composition analysis showed that AT content was higher than GC content, and the protein-coding sequences preferred to end with A/T at the third codon position. Relative synonymous codon usage analysis revealed most over-represented codon ends with A or T. Parity plot analysis revealed the codon usage bias of mitochondrial genes was influenced by both natural selection and mutation pressure. In the neutrality plot, the slopes of regression lines were < 0.5, suggesting that natural selection was playing a major role while mutation pressure was of minor importance. The effective number of codons showed that the codon usage bias between genes and genomes was low. Correspondence analysis revealed that the codon usage pattern differed among 13 protein-coding genes. Phylogenetic analyses based on three datasets using two methods (maximum likelihood and Bayesian inference), restored the Megophthalminae monophyly with high support values (bootstrap support values (BS) = 100, Bayesian posterior probability (PP) = 1). In the obtained topology, the seven Japanagallia species were clustered into a monophyletic group and formed a sister group with Durgade. In conclusion, our study can provide a reference for the future research on organism evolution, identification and phylogeny relationships of Japanagallia species.

Active role of the protein translation machinery in protecting against stress tolerance in Synechococcus elongatus PCC7942

In vivo protein synthesis is crucial for all domains of life. It is accomplished through translational machinery, and a key step is the translocation of tRNA-mRNA by elongation factor G (EF-G). Genome-based analysis revealed two EF-G encoding genes (S0885 and S2082) in the freshwater cyanobacterium model Synechococcus elongatus PCC7942. S0885 is the essential EF-G gene for photosynthesis. We generated a strain of S. elongatus PCC7942 that overexpressed S0885 (OX-S0885) to identify EF-G functionality. RT-PCR and Western blot analyses revealed increased transcriptional and translational levels in OX-S0885 at 10.5-13.5 and 2.0-3.0 fold, respectively. Overexpression of S0885 led to an increase in specific growth rate. Additionally, polysome-to-monosome ratio (P/M) and RNA-to-protein ratio (R/P) were elevated in OX-S0885 compared with the empty vector. Interestingly, R/P in OX-S0885 was retained at more than 70% under oxidative stress while R/P in the empty vector was severely depleted, suggesting the maintenance of translation. Thus, S0885 appeared to be the important target of oxidative stress because it was protected by the stress response system to maintain its function. These results suggest that cyanobacterial EF-G has a primary function in translation and an unrelated activity during stress conditions. These findings support the substantial role of EF-G in the formation and maintenance of cellular protein formation, and in the protection of the global translational mechanism under oxidative stress condition.

A beneficial synonymous substitution in EF-Tu is contingent on genetic background

Synonymous mutations are changes to DNA sequence that occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation Elongation Factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels, as well as the polysome abundance on global transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.

Effect of mRNA/tRNA mutations on translation speed: Implications for human diseases

Recent discoveries establish tRNAs as central regulators of mRNA translation dynamics, and therefore cotranslational folding and function of the encoded protein. The tRNA pool, whose composition and abundance change in a cell- and tissue-dependent manner, is the main factor which determines mRNA translation velocity. In this review, we discuss a group of pathogenic mutations, in the coding sequences of either protein-coding genes or in tRNA genes, that alter mRNA translation dynamics. We also summarize advances in tRNA biology that have uncovered how variations in tRNA levels on account of genetic mutations affect protein folding and function, and thereby contribute to phenotypic diversity in clinical manifestations.

BamQuery: a proteogenomic tool to explore the immunopeptidome and prioritize actionable tumor antigens

MHC-I-associated peptides deriving from non-coding genomic regions and mutations can generate tumor-specific antigens, including neoantigens. Quantifying tumor-specific antigens' RNA expression in malignant and benign tissues is critical for discriminating actionable targets. We present BamQuery, a tool attributing an exhaustive RNA expression to MHC-I-associated peptides of any origin from bulk and single-cell RNA-sequencing data. We show that many cryptic and mutated tumor-specific antigens can derive from multiple discrete genomic regions, abundantly expressed in normal tissues. BamQuery can also be used to predict MHC-I-associated peptides immunogenicity and identify actionable tumor-specific antigens de novo.

Optimization and Deoptimization of Codons in SARS-CoV-2 and Related Implications for Vaccine Development

The spread of coronavirus disease 2019 (COVID-19), caused by severe respiratory syndrome coronavirus 2 (SARS-CoV-2), has progressed into a global pandemic. To date, thousands of genetic variants have been identified among SARS-CoV-2 isolates collected from patients. Sequence analysis reveals that the codon adaptation index (CAI) values of viral sequences have decreased over time but with occasional fluctuations. Through evolution modeling, it is found that this phenomenon may result from the virus's mutation preference during transmission. Using dual-luciferase assays, it is further discovered that the deoptimization of codons in the viral sequence may weaken protein expression during virus evolution, indicating that codon usage may play an important role in virus fitness. Finally, given the importance of codon usage in protein expression and particularly for mRNA vaccines, it is designed several codon-optimized Omicron BA.2.12.1, BA.4/5, and XBB.1.5 spike mRNA vaccine candidates and experimentally validated their high levels of expression. This study highlights the importance of codon usage in virus evolution and provides guidelines for codon optimization in mRNA and DNA vaccine development.

Heterologous protein production in filamentous fungi

Filamentous fungi are able to produce a wide range of valuable proteins and enzymes for many industrial applications. Recent advances in fungal genomics and experimental technologies are rapidly changing the approaches for the development and use of filamentous fungi as hosts for the production of both homologous and heterologous proteins. In this review, we highlight the benefits and challenges of using filamentous fungi for the production of heterologous proteins. We review various techniques commonly employed to improve the heterologous protein production in filamentous fungi, such as strong and inducible promoters, codon optimization, more efficient signal peptides for secretion, carrier proteins, engineering of glycosylation sites, regulation of the unfolded protein response and endoplasmic reticulum associated protein degradation, optimization of the intracellular transport process, regulation of unconventional protein secretion, and construction of protease-deficient strains. KEY POINTS: • This review updates the knowledge on heterologous protein production in filamentous fungi. • Several fungal cell factories and potential candidates are discussed. • Insights into improving heterologous gene expression are given.

Translation and mRNA Stability Control

Messenger RNA (mRNA) stability and translational efficiency are two crucial aspects of the post-transcriptional process that profoundly impact protein production in a cell. While it is widely known that ribosomes produce proteins, studies during the past decade have surprisingly revealed that ribosomes also control mRNA stability in a codon-dependent manner, a process referred to as codon optimality. Therefore, codons, the three-nucleotide words read by the ribosome, have a potent effect on mRNA stability and provide cis-regulatory information that extends beyond the amino acids they encode. While the codon optimality molecular mechanism is still unclear, the translation elongation rate appears to trigger mRNA decay. Thus, transfer RNAs emerge as potential master gene regulators affecting mRNA stability. Furthermore, while few factors related to codon optimality have been identified in yeast, the orthologous genes in vertebrates do not necessary share the same functions. Here, we discuss codon optimality findings and gene regulation layers related to codon composition in different eukaryotic species.

A viral pan-end RNA element and host complex define a SARS-CoV-2 regulon

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, generates multiple protein-coding, subgenomic RNAs (sgRNAs) from a longer genomic RNA, all bearing identical termini with poorly understood roles in regulating viral gene expression. Insulin and interferon-gamma, two host-derived, stress-related agents, and virus spike protein, induce binding of glutamyl-prolyl-tRNA synthetase (EPRS1), within an unconventional, tetra-aminoacyl-tRNA synthetase complex, to the sgRNA 3'-end thereby enhancing sgRNA expression. We identify an EPRS1-binding sarbecoviral pan-end activating RNA (SPEAR) element in the 3'-end of viral RNAs driving agonist-induction. Translation of another co-terminal 3'-end feature, ORF10, is necessary for SPEAR-mediated induction, independent of Orf10 protein expression. The SPEAR element enhances viral programmed ribosomal frameshifting, thereby expanding its functionality. By co-opting noncanonical activities of a family of essential host proteins, the virus establishes a post-transcriptional regulon stimulating global viral RNA translation. A SPEAR-targeting strategy markedly reduces SARS-CoV-2 titer, suggesting a pan-sarbecoviral therapeutic modality.

Deoptimization of FMDV P1 Region Results in Robust Serotype-Independent Viral Attenuation

Foot-and-mouth disease (FMD), caused by the FMD virus (FMDV), is a highly contagious disease of cloven-hoofed livestock that can have severe economic impacts. Control and prevention strategies, including the development of improved vaccines, are urgently needed to effectively control FMD outbreaks in endemic settings. Previously, we employed two distinct strategies (codon pair bias deoptimization (CPD) and codon bias deoptimization (CD)) to deoptimize various regions of the FMDV serotype A subtype A12 genome, which resulted in the development of an attenuated virus in vitro and in vivo, inducing varying levels of humoral responses. In the current study, we examined the versatility of the system by using CPD applied to the P1 capsid coding region of FMDV serotype A subtype, A24, and another serotype, Asia1. Viruses carrying recoded P1 (A24-P1Deopt or Asia1-P1Deopt) exhibited different degrees of attenuation (i.e., delayed viral growth kinetics and replication) in cultured cells. Studies in vivo using a mouse model of FMD demonstrated that inoculation with the A24-P1Deopt and Asia1-P1Deopt strains elicited a strong humoral immune response capable of offering protection against challenge with homologous wildtype (WT) viruses. However, different results were obtained in pigs. While clear attenuation was detected for both the A24-P1Deopt and Asia1-P1Deopt strains, only a limited induction of adaptive immunity and protection against challenge was detected, depending on the inoculated dose and serotype deoptimized. Our work demonstrates that while CPD of the P1 coding region attenuates viral strains of multiple FMDV serotypes/subtypes, a thorough assessment of virulence and induction of adaptive immunity in the natural host is required in each case in order to finely adjust the degree of deoptimization required for attenuation without affecting the induction of protective adaptive immune responses.

Design and application of artificial rare L-lysine codons in Corynebacterium glutamicum

Background: L-lysine is widely used in the feed, food, and pharmaceutical industries, and screening for high L-lysine-producing strains has become a key goal for the industry. Methods: We constructed the rare L-lysine codon AAA by corresponding tRNA promoter replacement in C. glutamicum. Additionally, a screening marker related to the intracellular L-lysine content was constructed by converting all L-lysine codons of enhanced green fluorescent protein (EGFP) into the artificial rare codon AAA. The artificial EGFP was then ligated into pEC-XK99E and transformed into competent Corynebacterium glutamicum 23604 cells with the rare L-lysine codon. After atmospheric and room-temperature plasma mutation and induction culture, 55 mutants (0.01% of total cells) with stronger fluorescence were sorted using flow cytometry, and further screened by fermentation in a 96-deep-well plate and 500 mL shaker. Results: The fermentation results showed that the L-lysine production was increased by up to 9.7% in the mutant strains with higher fluorescence intensities, and that the highest screening positive rate was 69%, compared with that in the wild-type strain. Conclusion: The application of artificially constructed rare codons in this study represents an efficient, accurate, and simple method for screening other amino acid-producing microorganisms.

Structural Characteristics and Phylogenetic Analysis of the Mitochondrial Genomes of Four Krisna Species (Hemiptera: Cicadellidae: Iassinae)

Krisna species are insects that have piercing-sucking mouthparts and belong to the Krisnini tribe in the Iassinae subfamily of leafhoppers in the Cicadellidae family. In this study, we sequenced and compared the mitochondrial genomes (mitogenomes) of four Krisna species. The results showed that all four mitogenomes were composed of cyclic double-stranded molecules and contained 13 protein-coding genes (PCGs) and 22 and 2 genes coding for tRNAs and rRNAs, respectively. Those mitogenomes exhibited similar base composition, gene size, and codon usage patterns for the protein-coding genes. The analysis of the nonsynonymous substitution rate (Ka)/synonymous substitution rate (Ks) showed that evolution occurred the fastest in ND4 and the slowest in COI. 13 PCGs that underwent purification selection were suitable for studying phylogenetic relationships within Krisna. ND2, ND6, and ATP6 had highly variable nucleotide diversity, whereas COI and ND1 exhibited the lowest diversity. Genes or gene regions with high nucleotide diversity can provide potential marker candidates for population genetics and species delimitation in Krisna. Analyses of parity and neutral plots showed that both natural selection and mutation pressure affected the codon usage bias. In the phylogenetic analysis, all subfamilies were restored to a monophyletic group; the Krisnini tribe is monophyletic, and the Krisna genus is paraphyletic. Our study provides novel insights into the significance of the background nucleotide composition and codon usage patterns in the CDSs of the 13 mitochondrial PCGs of the Krisna genome, which could enable the identification of a different gene organization and may be used for accurate phylogenetic analysis of Krisna species.

Quality control of protein synthesis in the early elongation stage

In the early stage of bacterial translation, peptidyl-tRNAs frequently dissociate from the ribosome (pep-tRNA drop-off) and are recycled by peptidyl-tRNA hydrolase. Here, we establish a highly sensitive method for profiling of pep-tRNAs using mass spectrometry, and successfully detect a large number of nascent peptides from pep-tRNAs accumulated in Escherichia coli pthts strain. Based on molecular mass analysis, we found about 20% of the peptides bear single amino-acid substitutions of the N-terminal sequences of E. coli ORFs. Detailed analysis of individual pep-tRNAs and reporter assay revealed that most of the substitutions take place at the C-terminal drop-off site and that the miscoded pep-tRNAs rarely participate in the next round of elongation but dissociate from the ribosome. These findings suggest that pep-tRNA drop-off is an active mechanism by which the ribosome rejects miscoded pep-tRNAs in the early elongation, thereby contributing to quality control of protein synthesis after peptide bond formation.

Schlafen 12 restricts HIV-1 latency reversal by a codon-usage dependent post-transcriptional block in CD4+ T cells

Latency is a major barrier towards virus elimination in HIV-1-infected individuals. Yet, the mechanisms that contribute to the maintenance of HIV-1 latency are incompletely understood. Here we describe the Schlafen 12 protein (SLFN12) as an HIV-1 restriction factor that establishes a post-transcriptional block in HIV-1-infected cells and thereby inhibits HIV-1 replication and virus reactivation from latently infected cells. The inhibitory activity is dependent on the HIV-1 codon usage and on the SLFN12 RNase active sites. Within HIV-1-infected individuals, SLFN12 expression in PBMCs correlated with HIV-1 plasma viral loads and proviral loads suggesting a link with the general activation of the immune system. Using an RNA FISH-Flow HIV-1 reactivation assay, we demonstrate that SLFN12 expression is enriched in infected cells positive for HIV-1 transcripts but negative for HIV-1 proteins. Thus, codon-usage dependent translation inhibition of HIV-1 proteins participates in HIV-1 latency and can restrict the amount of virus release after latency reversal.

Insight on genetic features prevalent in five Ipomoea species using comparative codon pattern analysis reveals differences in major codons and reduced GC content at the 5’ end of CDS

Ipomoea plants possess important commercial, medicinal, and ornamental value. Molecular and morphological studies have confirmed that most species of this genus exhibit similar phenotypes but complex phylogenetic relationships. To date, limited information is available on these evolutionary relationships. In this study, systematic analysis of diverse species from Ipomoea was used to elucidate the relationships in this genus. To this end, we employed the concept of codon usage bias (CUB) to analyze the codon usage bias of five Ipomoea species such as effective number of codons (ENC) and GC content at the third synonym codon position (GC3s). Three types of plots including ENC-GC3s, parity rule 2 (PR2) and neutrality plots were employed to discover the factors determining CUB, and the frequency of hydrogen bonds and nucleotide were calculated to dissect changes in GC content at the 5'-end of the coding sequence. Our results showed little distinctness in CUB among the five species, with a reduction of hydrogen bonds content at the 5'-end (with similar changes in cytosines). In addition, optimal codons of Ipomoea aquatica ended with G or C, different from those of the other four species, which ended in A or T. These results may be useful for exploring the evolutionary relationships among this group, and for understanding the reasons for the variation among Ipomoea species.

Ultradeep characterisation of translational sequence determinants refutes rare-codon hypothesis and unveils quadruplet base pairing of initiator tRNA and transcript

Translation is a key determinant of gene expression and an important biotechnological engineering target. In bacteria, 5'-untranslated region (5'-UTR) and coding sequence (CDS) are well-known mRNA parts controlling translation and thus cellular protein levels. However, the complex interaction of 5'-UTR and CDS has so far only been studied for few sequences leading to non-generalisable and partly contradictory conclusions. Herein, we systematically assess the dynamic translation from over 1.2 million 5'-UTR-CDS pairs in Escherichia coli to investigate their collective effect using a new method for ultradeep sequence-function mapping. This allows us to disentangle and precisely quantify effects of various sequence determinants of translation. We find that 5'-UTR and CDS individually account for 53% and 20% of variance in translation, respectively, and show conclusively that, contrary to a common hypothesis, tRNA abundance does not explain expression changes between CDSs with different synonymous codons. Moreover, the obtained large-scale data provide clear experimental evidence for a base-pairing interaction between initiator tRNA and mRNA beyond the anticodon-codon interaction, an effect that is often masked for individual sequences and therefore inaccessible to low-throughput approaches. Our study highlights the indispensability of ultradeep sequence-function mapping to accurately determine the contribution of parts and phenomena involved in gene regulation.

Quantification of elongation stalls and impact on gene expression in yeast

Ribosomal pauses are a critical part of co-translational events including protein folding and localization. However, extended ribosome pauses can lead to ribosome collisions, resulting in the activation of ribosome rescue pathways and turnover of protein and mRNA. While this relationship has been known, the specific threshold between permissible pausing versus activation of rescue pathways has not been quantified. We have taken a method used to measure elongation time and adapted it for use in S. cerevisiae to quantify the impact of elongation stalls. We find, in transcripts containing Arg CGA codon repeat-induced stalls, a Hel2-mediated dose-dependent decrease in protein expression and mRNA level and an elongation delay on the order of minutes. In transcripts that contain synonymous substitutions to non-optimal Leu codons, there is a decrease in protein and mRNA levels, as well as similar elongation delay, but this occurs through a non-Hel2-mediated mechanism. Finally, we find that Dhh1 selectively increases protein expression, mRNA level, and elongation rate. This indicates that distinct poorly translated codons in an mRNA will activate different rescue pathways despite similar elongation stall durations. Taken together, these results provide new quantitative mechanistic insight into the surveillance of translation and the roles of Hel2 and Dhh1 in mediating ribosome pausing events.

Effects of synonymous mutations on kinetic properties and structure of firefly luciferase: Molecular dynamics simulation, molecular docking, RNA folding, and experimental study

Although synonymous mutations have long been thought to lack striking results, a growing body of research shows these mutations have highly variable effects. In this study, the impact of synonymous mutations in the development of thermostable luciferase was investigated using a combination of experimental and theoretical approaches. Using bioinformatics analysis, the codon usage features in the Lampyridae family's luciferases were studied and four synonymous mutations of Arg in luciferase were created. An exciting result was that the analysis of kinetic parameters showed a slight increase in the thermal stability of the mutant luciferase. AutoDock Vina, %MinMax algorithm, and UNAFold Server were used to perform molecular docking, folding rate, and RNA folding, respectively. Here, it was assumed that in the region (Arg337) with a moderate propensity for coil, synonymous mutation altered the rate of translation, which in turn may lead to a slight change in the structure of the enzyme. According to the molecular dynamics simulation data, local minor global flexibility is observed in the context of the protein conformation. A plausible explanation is that this flexibility may strengthen hydrophobic interactions due to its sensitivity to a molecular collision. Accordingly, thermostability originated mainly from hydrophobic interaction.

Ras protein abundance correlates with Ras isoform mutation patterns in cancer

Activating mutations of Ras genes are often observed in cancer. The protein products of the three Ras genes are almost identical. However, for reasons that remain unclear, KRAS is far more frequently mutated than the other Ras isoforms in cancer and RASopathies. We have quantified HRAS, NRAS, KRAS4A and KRAS4B protein abundance across a large panel of cell lines and healthy tissues. We observe consistent patterns of KRAS > NRAS»HRAS protein expression in cells that correlate with the rank order of Ras mutation frequencies in cancer. Our data provide support for the model of a sweet-spot of Ras dosage mediating isoform-specific contributions to cancer and development. We suggest that in most cases, being the most abundant Ras isoform correlates with occupying the sweet-spot and that HRAS and NRAS expression is usually insufficient to promote oncogenesis when mutated. However, our results challenge the notion that rare codons mechanistically underpin the predominance of KRAS mutant cancers. Finally, direct measurement of mutant versus wildtype KRAS protein abundance revealed a frequent imbalance that may suggest additional non-gene duplication mechanisms for optimizing oncogenic Ras dosage.

How synonymous mutations alter enzyme structure and function over long timescales

The specific activity of enzymes can be altered over long timescales in cells by synonymous mutations that alter a messenger RNA molecule's sequence but not the encoded protein's primary structure. How this happens at the molecular level is unknown. Here, we use multiscale modelling of three Escherichia coli enzymes (type III chloramphenicol acetyltransferase, D-alanine-D-alanine ligase B and dihydrofolate reductase) to understand experimentally measured changes in specific activity due to synonymous mutations. The modelling involves coarse-grained simulations of protein synthesis and post-translational behaviour, all-atom simulations to test robustness and quantum mechanics/molecular mechanics calculations to characterize enzymatic function. We show that changes in codon translation rates induced by synonymous mutations cause shifts in co-translational and post-translational folding pathways that kinetically partition molecules into subpopulations that very slowly interconvert to the native, functional state. Structurally, these states resemble the native state, with localized misfolding near the active sites of the enzymes. These long-lived states exhibit reduced catalytic activity, as shown by their increased activation energies for the reactions they catalyse.

Recoding of Nonsense Mutation as a Pharmacological Strategy

Approximately 11% of genetic human diseases are caused by nonsense mutations that introduce a premature termination codon (PTC) into the coding sequence. The PTC results in the production of a potentially harmful shortened polypeptide and activation of a nonsense-mediated decay (NMD) pathway. The NMD pathway reduces the burden of unproductive protein synthesis by lowering the level of PTC mRNA. There is an endogenous rescue mechanism that produces a full-length protein from a PTC mRNA. Nonsense suppression therapies aim to increase readthrough, suppress NMD, or are a combination of both strategies. Therefore, treatment with translational readthrough-inducing drugs (TRIDs) and NMD inhibitors may increase the effectiveness of PTC suppression. Here we discuss the mechanism of PTC readthrough and the development of novel approaches to PTC suppression. We also discuss the toxicity and bioavailability of therapeutics used to stimulate PTC readthrough.

Codon adaptation by synonymous mutations impacts the functional properties of the estrogen receptor-alpha protein in breast cancer cells

Oestrogen receptor-alpha (ERα) positivity is intimately associated with the development of hormone-dependent breast cancers. A major challenge in the treatment of these cancers is to understand and overcome the mechanisms of endocrine resistance. Recently, two distinct translation programmes using specific transfer RNA (tRNA) repertoires and codon usage frequencies were evidenced during cell proliferation and differentiation. Considering the phenotype switch of cancer cells to more proliferating and less-differentiated states, we can speculate that the changes in the tRNA pool and codon usage that likely occur make the ERα coding sequence no longer adapted, impacting translational rate, co-translational folding and the resulting functional properties of the protein. To verify this hypothesis, we generated an ERα synonymous coding sequence whose codon usage was optimized to the frequencies observed in genes expressed specifically in proliferating cells and then investigated the functional properties of the encoded receptor. We demonstrate that such a codon adaptation restores ERα activities to levels observed in differentiated cells, including: (a) an enhanced contribution exerted by transactivation function 1 (AF1) in ERα transcriptional activity; (b) enhanced interactions with nuclear receptor corepressor 1 and 2 [NCoR1 and NCoR2 (also known as SMRT) respectively], promoting repressive capability; and (c) reduced interactions with SRC proto-oncogene, non-receptor tyrosine kinase (Src) and phosphoinositide 3-kinase (PI3K) p85 kinases, inhibiting MAPK and AKT signalling pathway.

Identification of SNPs and expression patterns of ALB, AHSG and GC genes and their association with growth traits in Hu sheep

Growth traits are economically important traits in sheep breeding. This study was conducted to evaluate the polymorphisms of ALB, AHSG and GC genes and their association with growth traits in Hu sheep. We measured and recorded the body weight (BW), body height (BH), body length (BL) and feed conversion ratio (FCR) of 1418 male Hu sheep raised in the same environment from 80 to 180 days of age. The total of four SNPs in the ALB, AHSG and GC genes were identified by direct sequencing technology. The results of association analysis showed that two loci (g.8699 A>T and g.13458 T>C) of ALB gene significantly affect average daily gain (ADG; P < 0.05). The genotypes of SNP g.2454 T>C in AHSG gene were significantly associated with ADG and FCR (P < 0.05). There were significant associations between GC g.19484 A>C and BW, BH and BL (P < 0.05). The results of qRT-PCR showed that ALB, AHSG, and GC genes were extremely significantly higher in H_BW sheep compared with those in the L_BW sheep (P < 0.01). These results revealed that ALB-1 g.8699 A>T, ALB-2 g.13458 T>C, AHSG g.2454 T>C and GC g.19484 A>C loci are potential molecular markers for Hu sheep breeding.

mRNA-Based Therapeutics in Cancer Treatment

Over the past two decades, significant technological innovations have led to messenger RNA (mRNA) becoming a promising option for developing prophylactic and therapeutic vaccines, protein replacement therapies, and genome engineering. The success of the two COVID-19 mRNA vaccines has sparked new enthusiasm for other medical applications, particularly in cancer treatment. In vitro-transcribed (IVT) mRNAs are structurally designed to resemble naturally occurring mature mRNA. Delivery of IVT mRNA via delivery platforms such as lipid nanoparticles allows host cells to produce many copies of encoded proteins, which can serve as antigens to stimulate immune responses or as additional beneficial proteins for supplements. mRNA-based cancer therapeutics include mRNA cancer vaccines, mRNA encoding cytokines, chimeric antigen receptors, tumor suppressors, and other combination therapies. To better understand the current development and research status of mRNA therapies for cancer treatment, this review focused on the molecular design, delivery systems, and clinical indications of mRNA therapies in cancer.