Codon optimality is a major determinant of mRNA stability

Codon optimality influences homeostatic gene expression in zebrafish

The ribosome plays a crucial role in translating mRNA into protein; however, the genetic code extends beyond merely specifying amino acids. Upon translation, codons, the three-nucleotide sequences interpreted by ribosomes, have regulatory properties affecting mRNA stability, a phenomenon known as codon optimality. Codon optimality has been previously observed in vertebrates during embryogenesis, where specific codons can influence the stability and degradation rates of mRNA transcripts. In our previous work, we demonstrated that codon optimality impacts mRNA stability in human cell lines. However, the extent to which codon content influences vertebrate gene expression in vivo remained unclear. In this study, we expand on our previous findings by demonstrating that codon optimality has a robust effect on homeostatic mRNA and protein levels in whole zebrafish during normal physiological conditions. Using reporters with nearly identical nucleotide sequences but different codon compositions, all expressed from the same genomic locus, we show that codon composition can significantly influence gene expression. This study provides new insights into the regulatory roles of codon usage in vertebrate gene expression and underscores the importance of considering codon optimality in genetic and translational research. These findings have broad implications for understanding the complexities of gene regulation and could inform the design of synthetic genes and therapeutic strategies targeting mRNA stability.

Advances in the design and delivery of RNA vaccines for infectious diseases

RNA medicines represent a paradigm shift in treatment and prevention of critical diseases of global significance, e.g., infectious diseases. The highly successful messenger RNA (mRNA) vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) were developed at record speed during the coronavirus disease 2019 pandemic. A consequence of this is exceptionally shortened vaccine development times, which in combination with adaptability makes the RNA vaccine technology highly attractive against infectious diseases and for pandemic preparedness. Here, we review state of the art in the design and delivery of RNA vaccines for infectious diseases based on different RNA modalities, including linear mRNA, self-amplifying RNA, trans-amplifying RNA, and circular RNA. We provide an overview of the clinical pipeline of RNA vaccines for infectious diseases, and present analytical procedures, which are paramount for characterizing quality attributes and guaranteeing their quality, and we discuss future perspectives for using RNA vaccines to combat pathogens beyond SARS-CoV-2.

Dengue virus preferentially uses human and mosquito non-optimal codons

Codon optimality refers to the effect that codon composition has on messenger RNA (mRNA) stability and translation level and implies that synonymous codons are not silent from a regulatory point of view. Here, we investigated the adaptation of virus genomes to the host optimality code using mosquito-borne dengue virus (DENV) as a model. We demonstrated that codon optimality exists in mosquito cells and showed that DENV preferentially uses nonoptimal (destabilizing) codons and avoids codons that are defined as optimal (stabilizing) in either human or mosquito cells. Human genes enriched in the codons preferentially and frequently used by DENV are upregulated during infection, and so is the tRNA decoding the nonoptimal and DENV preferentially used codon for arginine. We found that adaptation during single-host passaging in human or mosquito cells results in the selection of synonymous mutations towards DENV's preferred nonoptimal codons that increase virus fitness. Finally, our analyses revealed that hundreds of viruses preferentially use nonoptimal codons, with those infecting a single host displaying an even stronger bias, suggesting that host-pathogen interaction shapes virus-synonymous codon choice.

CbCyp51-Mediated Demethylation Inhibitor Resistance Is Modulated by Codon Bias

Cercospora leaf spot, caused by the fungus Cercospora beticola, is the most destructive foliar disease of sugarbeet worldwide. Resistance to the sterol demethylation inhibitor (DMI) fungicide tetraconazole has been previously correlated with synonymous and nonsynonymous mutations in CbCyp51. Here, we extend these analyses to the DMI fungicides prothioconazole, difenoconazole, and mefentrifluconazole in addition to tetraconazole to confirm whether the synonymous and nonsynonymous mutations at amino acid positions 144 and 170 are associated with resistance to these fungicides. Nearly half of the 593 isolates of C. beticola collected in the Red River Valley of North Dakota and Minnesota in 2021 were resistant to all four DMIs. Another 20% were resistant to tetraconazole and prothioconazole but sensitive to difenoconazole and mefentrifluconazole. A total of 13% of isolates were sensitive to all DMIs tested. We found five CbCyp51 haplotypes and associated them with phenotypes to the four DMIs. The most predominant haplotype (E170_A/L144F_C) correlated with resistance to all four DMIs with up to 97.6% accuracy. The second most common haplotype (E170_A/L144) consisted of isolates associated with resistance phenotypes to tetraconazole and prothioconazole while also exhibiting sensitive phenotypes to difenoconazole and mefentrifluconazole with up to 98.4% accuracy. Quantitative PCR did not identify differences in CbCyp51 expression between haplotypes. This study offers an understanding of the importance of codon usage in fungicide resistance and provides crop management acuity for fungicide application decision-making.

Protein folding and quality control during nuclear transport

The spatial separation of protein synthesis from the compartmental destiny of proteins led to the evolution of transport systems that are efficient and yet highly specific. Co-translational transport has emerged as a strategy to avoid cytosolic aggregation of folding intermediates and the need for energy-consuming unfolding strategies to enable transport through narrow conduits connecting compartments. While translation and compartmental translocation are at times tightly coordinated, we know very little about the temporal coordination of translation, protein folding, and nuclear import. Here, we consider the implications of co-translational engagement of nuclear import machinery. We propose that the dynamic interplay of karyopherins and intrinsically disordered nucleoporins create a favorable protein folding environment for cargo en route to the nuclear compartment while maintaining a barrier function of the nuclear pore complex. Our model is discussed in the context of neurological disorders that are tied to defects in nuclear transport and protein quality control.

Unlocking the potential of Extensin Signal peptide and Elastin-like polypeptide tag fused to Shigella dysenteriae's IpaDSTxB to improve protein expression and purification in Nicotiana tabacum and Medicagosativa

Plants are often seen as a potent tool in the recombinant protein production industry. However, unlike bacterial expression, it is not a popular method due to the low yield and difficulty of protein extraction and purification. Therefore, developing a new high efficient and easy to purify platform is crucial. One of the best approaches to make extraction easier is to utilize the Extensin Signal peptide (EXT) to translocate the recombinant protein to the outside of the cell, along with incorporating an Elastin-like polypeptide tag (ELP) to enhance purification and accumulation rates. In this research, we transiently expressed Shigella dysenteriae's IpaDSTxB fused to both NtEXT and ELP in both Nicotiana tabacum and Medicago sativa. Our results demonstrated that N. tabacum, with an average yield of 6.39 ng/μg TSP, outperforms M. sativa, which had an average yield of 3.58 ng/μg TSP. On the other hand, analyzing NtEXT signal peptide indicated that merging EXT to the constructs facilitates translocation of IpaDSTxB to the apoplast by 78.4% and 65.9% in N. tabacum and M. sativa, respectively. Conversely, the mean level for constructs without EXT was below 25% for both plants. Furthermore, investigation into the orientation of ELP showed that merging it to the C-terminal of IpaDSTxB leads to a higher accumulation rate in both N. tabacum and M. sativa by 1.39 and 1.28 times, respectively. It also facilitates purification rate by over 70% in comparison to 20% of the 6His tag. The results show a highly efficient and easy to purify platform for the expression of heterologous proteins in plant.

Production and evaluation of three kinds of vaccines against largemouth bass virus, and DNA vaccines show great application prospects

Largemouth bass virus (LMBV) infections has resulted in high mortality and economic losses to the global largemouth bass industry and has seriously restricted the healthy development of the bass aquaculture industry. There are currently no antiviral therapies available for the control of this disease. In this study, we developed three types of vaccine against LMBV; whole virus inactivated vaccine (I), a subunit vaccine composed of the major viral capsid protein MCP (S) as well as an MCP DNA vaccine(D), These were employed using differing immunization and booster strategies spaced 2 weeks apart as follows: II, SS, DD and DS. We found that all vaccine groups induced humoral and cellular immune responses and protected largemouth bass from a lethal LMBV challenge to varying degrees and DD produced the best overall effect. Specifically, the levels of specific IgM in serum in all immunized groups were elevated and significantly higher than those in the control group. Moreover, the expression of humoral immunity (CD4 and IgM) and cellular immunity (MHCI-α) as well as cytokines (IL-1β) was increased, and the activity of immunity-related enzymes ACP, AKP, LZM, and T-SOD in the serum was significantly enhanced. In addition, the relative percent survival of fish following an LMBV lethal challenge 4 weeks after the initial immunizations were high for each group: DD(89.5 %),DS(63.2 %),SS(50 %) and II (44.7 %). These results indicated that the MCP DNA vaccine is the most suitable and promising vaccine candidate for the effective control of LMBV disease.

The regulatory landscape of interacting RNA and protein pools in cellular homeostasis and cancer

Biological systems encompass intricate networks governed by RNA-protein interactions that play pivotal roles in cellular functions. RNA and proteins constituting 1.1% and 18% of the mammalian cell weight, respectively, orchestrate vital processes from genome organization to translation. To date, disentangling the functional fraction of the human genome has presented a major challenge, particularly for noncoding regions, yet recent discoveries have started to unveil a host of regulatory functions for noncoding RNAs (ncRNAs). While ncRNAs exist at different sizes, structures, degrees of evolutionary conservation and abundances within the cell, they partake in diverse roles either alone or in combination. However, certain ncRNA subtypes, including those that have been described or remain to be discovered, are poorly characterized given their heterogeneous nature. RNA activity is in most cases coordinated through interactions with RNA-binding proteins (RBPs). Extensive efforts are being made to accurately reconstruct RNA-RBP regulatory networks, which have provided unprecedented insight into cellular physiology and human disease. In this review, we provide a comprehensive view of RNAs and RBPs, focusing on how their interactions generate functional signals in living cells, particularly in the context of post-transcriptional regulatory processes and cancer.

A spatiotemporally resolved atlas of mRNA decay in the C. elegans embryo reveals differential regulation of mRNA stability across stages and cell types

During embryonic development, cells undergo dynamic changes in gene expression that are required for appropriate cell fate specification. Although both transcription and mRNA degradation contribute to gene expression dynamics, patterns of mRNA decay are less well understood. Here, we directly measure spatiotemporally resolved mRNA decay rates transcriptome-wide throughout C. elegans embryogenesis by transcription inhibition followed by bulk and single-cell RNA sequencing. This allows us to calculate mRNA half-lives within specific cell types and developmental stages, and identify differentially regulated mRNA decay throughout embryonic development. We identify transcript features that are correlated with mRNA stability and find that mRNA decay rates are associated with distinct peaks in gene expression over time. Moreover, we provide evidence that, on average, mRNA is more stable in the germline than in the soma and in later embryonic stages than in earlier stages. This work suggests that differential mRNA decay across cell states and time helps to shape developmental gene expression, and it provides a valuable resource for studies of mRNA turnover regulatory mechanisms.

The DNA damage response and RNA Polymerase II regulator Def1 has posttranscriptional functions in the cytoplasm

Yeast Def1 mediates RNA polymerase II degradation and transcription elongation during stress. Def1 is predominantly cytoplasmic, and DNA damage signals cause its proteolytic processing, liberating its N-terminus to enter the nucleus. Cytoplasmic functions for this abundant protein have not been identified. Proximity-labeling (BioID) experiments indicate that Def1 binds to an array of proteins involved in posttranscriptional control and translation of mRNAs. Deleting DEF1 reduces both mRNA synthesis and decay rates, indicating transcript buffering in the mutant. Directly tethering Def1 to a reporter mRNA suppressed expression, suggesting that Def1 directly regulates mRNAs. Surprisingly, we found that Def1 interacts with polyribosomes, which requires its ubiquitin-binding domain located in its N-terminus. The binding of Def1 to ribosomes requires the ubiquitylation of eS7a (Rsp7A) in the small subunit by the Not4 protein in the Ccr4-Not complex. Not4 ubiquitylation of the ribosome regulates translation quality control and co-translational mRNA decay. The polyglutamine-rich unstructured C-terminus of Def1 is required for its interaction with decay and translation factors, suggesting that Def1 acts as a ubiquitin-dependent scaffold to link translation status to mRNA decay. Thus, we have identified a novel function for this transcription and DNA damage response factor in posttranscriptional regulation in the cytoplasm and establish Def1 as a master regulator of gene expression, functioning during transcription, mRNA decay, and translation.

Transfer RNA levels are tuned to support differentiation during Drosophila neurogenesis

Neural differentiation requires a multifaceted program to alter gene expression along the proliferation to differentiation axis. While critical changes occur at the level of transcription, post-transcriptional mechanisms allow fine-tuning of protein output. We investigated the role of tRNAs in regulating gene expression during neural differentiation by quantifying tRNA abundance in neural progenitor-biased and neuron-biased Drosophila larval brains. We found that tRNA profiles are largely consistent between progenitor-biased and neuron-biased brains but significant variation occurs for 10 cytoplasmic isodecoders (individual tRNA genes) and this establishes differential tRNA levels for 8 anticodon groups. We used these tRNA data to investigate relationships between tRNA abundance, codon optimality-mediated mRNA decay, and translation efficiency in progenitors and neurons. Our data reveal that tRNA levels strongly correlate with codon optimality-mediated mRNA decay within each cell type but generally do not explain differences in stabilizing versus destabilizing codons between cell types. Regarding translation efficiency, we found that tRNA expression in neural progenitors preferentially supports translation of mRNAs whose products are in high demand in progenitors, such as those associated with protein synthesis. In neurons, tRNA expression shifts to disfavor translation of proliferation-related transcripts and preferentially support translation of transcripts tied to neuron-specific functions like axon pathfinding and synapse formation. Overall, our analyses reveal that changes in tRNA levels along the neural differentiation axis support optimal gene expression in progenitors and neurons.

Diverse intrinsic properties shape transcript stability and stabilization in Mycolicibacterium smegmatis

Mycobacteria regulate transcript degradation to facilitate adaptation to environmental stress. However, the mechanisms underlying this regulation are unknown. Here we sought to gain understanding of the mechanisms controlling mRNA stability by investigating the transcript properties associated with variance in transcript stability and stress-induced transcript stabilization. We measured mRNA half-lives transcriptome-wide in Mycolicibacterium smegmatis in log phase growth and hypoxia-induced growth arrest. The transcriptome was globally stabilized in response to hypoxia, but transcripts of essential genes were generally stabilized more than those of non-essential genes. We then developed machine learning models that enabled us to identify the non-linear collective effect of a compendium of transcript properties on transcript stability and stabilization. We identified properties that were more predictive of half-life in log phase as well as properties that were more predictive in hypoxia, and many of these varied between leadered and leaderless transcripts. In summary, we found that transcript properties are differentially associated with transcript stability depending on both the transcript type and the growth condition. Our results reveal the complex interplay between transcript features and microenvironment that shapes transcript stability in mycobacteria.

RNA modifications in physiology and pathology: Progressing towards application in clinical settings

The potential to modify individual nucleotides through chemical means in order to impact the electrostatic charge, hydrophobic properties, and base pairing of RNA molecules is harnessed in the medical application of stable synthetic RNAs like mRNA vaccines and synthetic small RNA molecules. These modifications are used to either increase or decrease the production of therapeutic proteins. Additionally, naturally occurring biochemical alterations of nucleotides play a role in regulating RNA metabolism and function, thereby modulating essential cellular processes. Research elucidating the mechanisms through which RNA modifications govern fundamental cellular functions in multicellular organisms has enhanced our comprehension of how irregular RNA modification profiles can lead to human diseases. Collectively, these fundamental scientific findings have unveiled the molecular and cellular functions of RNA modifications, offering new opportunities for therapeutic intervention and paving the way for a variety of innovative clinical strategies.

The central role of transfer RNAs in mistranslation

Transfer RNAs (tRNA) are essential small non-coding RNAs that enable the translation of genomic information into proteins in all life forms. The principal function of tRNAs is to bring amino acid building blocks to the ribosomes for protein synthesis. In the ribosome, tRNAs interact with messenger RNA (mRNA) to mediate the incorporation of amino acids into a growing polypeptide chain following the rules of the genetic code. Accurate interpretation of the genetic code requires tRNAs to carry amino acids matching their anticodon identity and decode the correct codon on mRNAs. Errors in these steps cause the translation of codons with the wrong amino acids (mistranslation), compromising the accurate flow of information from DNA to proteins. Accumulation of mutant proteins due to mistranslation jeopardizes proteostasis and cellular viability. However, the concept of mistranslation is evolving, with increasing evidence indicating that mistranslation can be used as a mechanism for survival and acclimatization to environmental conditions. In this review, we discuss the central role of tRNAs in modulating translational fidelity through their dynamic and complex interplay with translation factors. We summarize recent discoveries of mistranslating tRNAs and describe the underlying molecular mechanisms and the specific conditions and environments that enable and promote mistranslation.

Pleiotropic effects of PAB1 deletion: Extensive changes in the yeast proteome, transcriptome, and translatome

Cytoplasmic poly(A)-binding protein (PABPC; Pab1 in yeast) is thought to be involved in multiple steps of post-transcriptional control, including translation initiation, translation termination, and mRNA decay. To understand both the direct and indirect roles of PABPC in more detail, we have employed mass spectrometry to assess the abundance of the components of the yeast proteome, as well as RNA-Seq and Ribo-Seq to analyze changes in the abundance and translation of the yeast transcriptome, in cells lacking the PAB1 gene. We find that pab1Δ cells manifest drastic changes in the proteome and transcriptome, as well as defects in translation initiation and termination. Defects in translation initiation and the stabilization of specific classes of mRNAs in pab1Δ cells appear to be partly indirect consequences of reduced levels of specific initiation factors, decapping activators, and components of the deadenylation complex in addition to the general loss of Pab1's direct role in these processes. Cells devoid of Pab1 also manifested a nonsense codon readthrough phenotype indicative of a defect in translation termination. Collectively, our results indicate that, unlike the loss of simpler regulatory proteins, elimination of cellular Pab1 is profoundly pleiotropic and disruptive to numerous aspects of post-transcriptional regulation.

Translation elongation as a rate limiting step of protein production

The impact of synonymous codon choice on protein output has important implications for understanding endogenous gene expression and design of synthetic mRNAs. Synonymous codons are decoded at different speeds, but simple models predict that this should not drive protein output. Instead, translation initiation should be the rate limiting step for production of protein per mRNA, with little impact of codon choice. Previously, we used a neural network model to design a series of synonymous fluorescent reporters and showed that their protein output in yeast spanned a seven-fold range corresponding to their predicted translation elongation speed. Here, we show that this effect is not due primarily to the established impact of slow elongation on mRNA stability, but rather, that slow elongation further decreases the number of proteins made per mRNA. We combine simulations and careful experiments on fluorescent reporters to show that translation is limited on non-optimally encoded transcripts. Using a genome-wide CRISPRi screen, we find that impairing translation initiation attenuates the impact of slow elongation, showing a dynamic balance between rate limiting steps of protein production. Our results show that codon choice can directly limit protein production across the full range of endogenous variability in codon usage.

Decapping activators Edc3 and Scd6 act redundantly with Dhh1 in post-transcriptional repression of starvation-induced pathways

Degradation of many yeast mRNAs involves decapping by the Dcp1:Dcp2 complex. Previous studies on decapping activators Edc3 and Scd6 suggested their limited roles in mRNA decay. RNA-seq analysis of mutants lacking one or both proteins revealed that Scd6 and Edc3 have largely redundant activities in targeting numerous mRNAs for degradation that are masked in the single mutants. These transcripts also are frequently targeted by decapping activators Dhh1 and Pat1, and the collective evidence suggests that Scd6/Edc3 act interchangeably to recruit Dhh1 to Dcp2. Ribosome profiling shows that redundancy between Scd6 and Edc3 and their functional interactions with Dhh1 and Pat1 extend to translational repression of particular transcripts, including a cohort of poorly translated mRNAs displaying interdependent regulation by all four factors. Scd6/Edc3 also participate with Dhh1/Pat1 in post-transcriptional repression of proteins required for respiration and catabolism of alternative carbon sources, which are normally expressed only in limiting glucose. Simultaneously eliminating Scd6/Edc3 increases mitochondrial membrane potential and elevates metabolites of the tricarboxylic acid and glyoxylate cycles typically observed only during growth in low glucose. Thus, Scd6/Edc3 act redundantly, in parallel with Dhh1 and in cooperation with Pat1, to adjust gene expression to nutrient availability by controlling mRNA decapping and decay.

Comparative Analysis of Codon Usage Bias in Six Eimeria Genomes

The codon usage bias (CUB) of genes encoded by different species' genomes varies greatly. The analysis of codon usage patterns enriches our comprehension of genetic and evolutionary characteristics across diverse species. In this study, we performed a genome-wide analysis of CUB and its influencing factors in six sequenced Eimeria species that cause coccidiosis in poultry: Eimeria acervulina, Eimeria necatrix, Eimeria brunetti, Eimeria tenella, Eimeria praecox, and Eimeria maxima. The GC content of protein-coding genes varies between 52.67% and 58.24% among the six Eimeria species. The distribution trend of GC content at different codon positions follows GC1 > GC3 > GC2. Most high-frequency codons tend to end with C/G, except in E. maxima. Additionally, there is a positive correlation between GC3 content and GC3s/C3s, but a significantly negative correlation with A3s. Analysis of the ENC-Plot, neutrality plot, and PR2-bias plot suggests that selection pressure has a stronger influence than mutational pressure on CUB in the six Eimeria genomes. Finally, we identified from 11 to 15 optimal codons, with GCA, CAG, and AGC being the most commonly used optimal codons across these species. This study offers a thorough exploration of the relationships between CUB and selection pressures within the protein-coding genes of Eimeria species. Genetic evolution in these species appears to be influenced by mutations and selection pressures. Additionally, the findings shed light on unique characteristics and evolutionary traits specific to the six Eimeria species.

Global regulation via modulation of ribosome pausing by the ABC-F protein EttA

Having multiple rounds of translation of the same mRNA creates dynamic complexities along with opportunities for regulation related to ribosome pausing and stalling at specific sequences. Yet, mechanisms controlling these critical processes and the principles guiding their evolution remain poorly understood. Through genetic, genomic, physiological, and biochemical approaches, we demonstrate that regulating ribosome pausing at specific amino acid sequences can produce ~2-fold changes in protein expression levels which strongly influence cell growth and therefore evolutionary fitness. We demonstrate, both in vivo and in vitro, that the ABC-F protein EttA directly controls the translation of mRNAs coding for a subset of enzymes in the tricarboxylic acid (TCA) cycle and its glyoxylate shunt, which modulates growth in some chemical environments. EttA also modulates expression of specific proteins involved in metabolically related physiological and stress-response pathways. These regulatory activities are mediated by EttA rescuing ribosomes paused at specific patterns of negatively charged residues within the first 30 amino acids of nascent proteins. We thus establish a unique global regulatory paradigm based on sequence-specific modulation of translational pausing.

Genome-wide quantification of RNA flow across subcellular compartments reveals determinants of the mammalian transcript life cycle

Dissecting the regulatory mechanisms controlling mammalian transcripts from production to degradation requires quantitative measurements of mRNA flow across the cell. We developed subcellular TimeLapse-seq to measure the rates at which RNAs are released from chromatin, exported from the nucleus, loaded onto polysomes, and degraded within the nucleus and cytoplasm in human and mouse cells. These rates varied substantially, yet transcripts from genes with related functions or targeted by the same transcription factors and RNA-binding proteins flowed across subcellular compartments with similar kinetics. Verifying these associations uncovered a link between DDX3X and nuclear export. For hundreds of RNA metabolism genes, most transcripts with retained introns were degraded by the nuclear exosome, while the remaining molecules were exported with stable cytoplasmic lifespans. Transcripts residing on chromatin for longer had extended poly(A) tails, whereas the reverse was observed for cytoplasmic mRNAs. Finally, machine learning identified molecular features that predicted the diverse life cycles of mRNAs.

Hypoxia and inflammation induce synergistic transcriptome turnover in macrophages

Macrophages are effector immune cells that experience substantial changes to oxygenation when transiting through tissues, especially when entering tumors or infected wounds. How hypoxia alters gene expression and macrophage effector function at the post-transcriptional level remains poorly understood. Here, we use TimeLapse-seq to measure how inflammatory activation modifies the hypoxic response in primary macrophages. Nucleoside recoding sequencing allows the derivation of steady-state transcript levels, degradation rates, and transcriptional synthesis rates from the same dataset. We find that hypoxia produces distinct responses from resting and inflammatory macrophages. Hypoxia induces destabilization of mRNA transcripts, though inflammatory macrophages substantially increase mRNA degradation compared to resting macrophages. Increased RNA turnover results in the upregulation of ribosomal protein genes and downregulation of extracellular matrix components in inflammatory macrophages. Pathways regulated by mRNA decay in vitro are differentially regulated in tumor-associated macrophages implying that mixed stimuli could induce post-transcriptional regulation of macrophage function in solid tumors.

Enrichment of rare codons at 5' ends of genes is a spandrel caused by evolutionary sequence turnover and does not improve translation

Previously, Tuller et al. found that the first 30-50 codons of the genes of yeast and other eukaryotes are slightly enriched for rare codons. They argued that this slowed translation, and was adaptive because it queued ribosomes to prevent collisions. Today, the translational speeds of different codons are known, and indeed rare codons are translated slowly. We re-examined this 5' slow translation 'ramp.' We confirm that 5' regions are slightly enriched for rare codons; in addition, they are depleted for downstream Start codons (which are fast), with both effects contributing to slow 5' translation. However, we also find that the 5' (and 3') ends of yeast genes are poorly conserved in evolution, suggesting that they are unstable and turnover relatively rapidly. When a new 5' end forms de novo, it is likely to include codons that would otherwise be rare. Because evolution has had a relatively short time to select against these codons, 5' ends are typically slightly enriched for rare, slow codons. Opposite to the expectation of Tuller et al., we show by direct experiment that genes with slowly translated codons at the 5' end are expressed relatively poorly, and that substituting faster synonymous codons improves expression. Direct experiment shows that slow codons do not prevent downstream ribosome collisions. Further informatic studies suggest that for natural genes, slow 5' ends are correlated with poor gene expression, opposite to the expectation of Tuller et al. Thus, we conclude that slow 5' translation is a 'spandrel'--a non-adaptive consequence of something else, in this case, the turnover of 5' ends in evolution, and it does not improve translation.

Human DDX6 regulates translation and decay of inefficiently translated mRNAs

Recent findings indicate that the translation elongation rate influences mRNA stability. One of the factors that has been implicated in this link between mRNA decay and translation speed is the yeast DEAD-box helicase Dhh1p. Here, we demonstrated that the human ortholog of Dhh1p, DDX6, triggers the deadenylation-dependent decay of inefficiently translated mRNAs in human cells. DDX6 interacts with the ribosome through the Phe-Asp-Phe (FDF) motif in its RecA2 domain. Furthermore, RecA2-mediated interactions and ATPase activity are both required for DDX6 to destabilize inefficiently translated mRNAs. Using ribosome profiling and RNA sequencing, we identified two classes of endogenous mRNAs that are regulated in a DDX6-dependent manner. The identified targets are either translationally regulated or regulated at the steady-state-level and either exhibit signatures of poor overall translation or of locally reduced ribosome translocation rates. Transferring the identified sequence stretches into a reporter mRNA caused translation- and DDX6-dependent degradation of the reporter mRNA. In summary, these results identify DDX6 as a crucial regulator of mRNA translation and decay triggered by slow ribosome movement and provide insights into the mechanism by which DDX6 destabilizes inefficiently translated mRNAs.

Codon usage bias in yeasts and its correlation with gene expression, growth temperature, and protein structure

Codon usage bias (CUB) has been described in viruses, prokaryotes, and eukaryotes and has been linked to several cellular and environmental factors, such as the organism's growth temperature, gene expression levels, and regulation of protein synthesis and folding. Most of the studies in this area have been conducted in bacteria and higher eukaryotes, in some cases with different results. In this study, a comparative analysis of CUB in yeasts isolated from cold and template environments was performed in order to evaluate the correlation of CUB with yeast optimal temperature of growth (OTG), gene expression levels, cellular function, and structure of encoded proteins. Among the main findings, highly expressed ORFs tend to have a more similar CUB within and between yeasts, and a direct correlation between codons ending in C and expression level was generally found. A low correspondence between CUB and OTG was observed, with an inverse correlation for some codons ending in C. The clustering of yeasts based on their CUB partially aligns with their OTG, being more consistent for yeasts with lower OTG. In most yeasts, the abundance of preferred codons was generally lower at the 5' end of ORFs, higher in segments encoding beta strand, lower in segments encoding extracellular and transmembrane regions, and higher in "translation" and "energy metabolism" pathways, especially in highly expressed ORFs. Based on our findings, it is suggested that the abundance and distribution of preferred and non-preferred codons along mRNAs contribute to proper protein folding and functionality by regulating protein synthesis rates, becoming a more important factor under conditions that require faster protein synthesis in yeasts.

Massively parallel identification of sequence motifs triggering ribosome-associated mRNA quality control

Decay of mRNAs can be triggered by ribosome slowdown at stretches of rare codons or positively charged amino acids. However, the full diversity of sequences that trigger co-translational mRNA decay is poorly understood. To comprehensively identify sequence motifs that trigger mRNA decay, we use a massively parallel reporter assay to measure the effect of all possible combinations of codon pairs on mRNA levels in S. cerevisiae. In addition to known mRNA-destabilizing sequences, we identify several dipeptide repeats whose translation reduces mRNA levels. These include combinations of positively charged and bulky residues, as well as proline-glycine and proline-aspartate dipeptide repeats. Genetic deletion of the ribosome collision sensor Hel2 rescues the mRNA effects of these motifs, suggesting that they trigger ribosome slowdown and activate the ribosome-associated quality control (RQC) pathway. Deep mutational scanning of an mRNA-destabilizing dipeptide repeat reveals a complex interplay between the charge, bulkiness, and location of amino acid residues in conferring mRNA instability. Finally, we show that the mRNA effects of codon pairs are predictive of the effects of endogenous sequences. Our work highlights the complexity of sequence motifs driving co-translational mRNA decay in eukaryotes, and presents a high throughput approach to dissect their requirements at the codon level.

SAMD9L acts as an antiviral factor against HIV-1 and primate lentiviruses by restricting viral and cellular translation

Sterile alpha motif domain-containing proteins 9 and 9-like (SAMD9/9L) are associated with life-threatening genetic diseases in humans and are restriction factors of poxviruses. Yet, their cellular function and the extent of their antiviral role are poorly known. Here, we found that interferon-stimulated human SAMD9L restricts HIV-1 in the late phases of replication, at the posttranscriptional and prematuration steps, impacting viral translation and, possibly, endosomal trafficking. Surprisingly, the paralog SAMD9 exerted an opposite effect, enhancing HIV-1. More broadly, we showed that SAMD9L restricts primate lentiviruses, but not a gammaretrovirus (MLV), nor 2 RNA viruses (arenavirus MOPV and rhabdovirus VSV). Using structural modeling and mutagenesis of SAMD9L, we identified a conserved Schlafen-like active site necessary for HIV-1 restriction by human and a rodent SAMD9L. By testing a gain-of-function constitutively active variant from patients with SAMD9L-associated autoinflammatory disease, we determined that SAMD9L pathogenic functions also depend on the Schlafen-like active site. Finally, we found that the constitutively active SAMD9L strongly inhibited HIV, MLV, and, to a lesser extent, MOPV. This suggests that the virus-specific effect of SAMD9L may involve its differential activation/sensing and the virus ability to evade from SAMD9L restriction. Overall, our study identifies SAMD9L as an HIV-1 antiviral factor from the cell autonomous immunity and deciphers host determinants underlying the translational repression. This provides novel links and therapeutic avenues against viral infections and genetic diseases.

A second hotspot for pathogenic exon-skipping variants in CDC45

Biallelic pathogenic variants in CDC45 are associated with Meier-Gorlin syndrome with craniosynostosis (MGORS type 7), which also includes short stature and absent/hypoplastic patellae. Identified variants act through a hypomorphic loss of function mechanism, to reduce CDC45 activity and impact DNA replication initiation. In addition to missense and premature termination variants, several pathogenic synonymous variants have been identified, most of which cause increased exon skipping of exon 4, which encodes an essential part of the RecJ-orthologue's DHH domain. Here we have identified a second cohort of families segregating CDC45 variants, where patients have craniosynostosis and a reduction in height, alongside common facial dysmorphisms, including thin eyebrows, consistent with MGORS7. Skipping of exon 15 is a consequence of two different variants, including a shared synonymous variant that is enriched in individuals of East Asian ancestry, while other variants in trans are predicted to alter key intramolecular interactions in α/β domain II, or cause retention of an intron within the 3'UTR. Our cohort and functional data confirm exon skipping is a relatively common pathogenic mechanism in CDC45, and highlights the need for alternative splicing events, such as exon skipping, to be especially considered for variants initially predicted to be less likely to cause the phenotype, particularly synonymous variants.

Engineering circular RNA for molecular and metabolic reprogramming

The role of messenger RNA (mRNA) in biological systems is extremely versatile. However, it's extremely short half-life poses a fundamental restriction on its application. Moreover, the translation efficiency of mRNA is also limited. On the contrary, circular RNAs, also known as circRNAs, are a common and stable form of RNA found in eukaryotic cells. These molecules are synthesized via back-splicing. Both synthetic circRNAs and certain endogenous circRNAs have the potential to encode proteins, hence suggesting the potential of circRNA as a gene expression machinery. Herein, we aim to summarize all engineering aspects that allow exogenous circular RNA (circRNA) to prolong the time that proteins are expressed from full-length RNA signals. This review presents a systematic engineering approach that have been devised to efficiently assemble circRNAs and evaluate several aspects that have an impact on protein production derived from. We have also reviewed how optimization of the key components of circRNAs, including the topology of vector, 5' and 3' untranslated sections, entrance site of the internal ribosome, and engineered aptamers could be efficiently impacting the translation machinery for molecular and metabolic reprogramming. Collectively, molecular and metabolic reprogramming present a novel way of regulating distinctive cellular features, for instance growth traits to neoplastic cells, and offer new possibilities for therapeutic inventions.

Mechanisms and Delivery of tRNA Therapeutics

Transfer ribonucleic acid (tRNA) therapeutics will provide personalized and mutation specific medicines to treat human genetic diseases for which no cures currently exist. The tRNAs are a family of adaptor molecules that interpret the nucleic acid sequences in our genes into the amino acid sequences of proteins that dictate cell function. Humans encode more than 600 tRNA genes. Interestingly, even healthy individuals contain some mutant tRNAs that make mistakes. Missense suppressor tRNAs insert the wrong amino acid in proteins, and nonsense suppressor tRNAs read through premature stop signals to generate full length proteins. Mutations that underlie many human diseases, including neurodegenerative diseases, cancers, and diverse rare genetic disorders, result from missense or nonsense mutations. Thus, specific tRNA variants can be strategically deployed as therapeutic agents to correct genetic defects. We review the mechanisms of tRNA therapeutic activity, the nature of the therapeutic window for nonsense and missense suppression as well as wild-type tRNA supplementation. We discuss the challenges and promises of delivering tRNAs as synthetic RNAs or as gene therapies. Together, tRNA medicines will provide novel treatments for common and rare genetic diseases in humans.

Orb2 enables rare-codon-enriched mRNA expression during Drosophila neuron differentiation

Regulation of codon optimality is an increasingly appreciated layer of cell- and tissue-specific protein expression control. Here, we use codon-modified reporters to show that differentiation of Drosophila neural stem cells into neurons enables protein expression from rare-codon-enriched genes. From a candidate screen, we identify the cytoplasmic polyadenylation element binding (CPEB) protein Orb2 as a positive regulator of rare-codon-dependent mRNA stability in neurons. Using RNA sequencing, we reveal that Orb2-upregulated mRNAs in the brain with abundant Orb2 binding sites have a rare-codon bias. From these Orb2-regulated mRNAs, we demonstrate that rare-codon enrichment is important for mRNA stability and social behavior function of the metabotropic glutamate receptor (mGluR). Our findings reveal a molecular mechanism by which neural stem cell differentiation shifts genetic code regulation to enable critical mRNA stability and protein expression.

2'-O-methylation at internal sites on mRNA promotes mRNA stability

2'-O-methylation (Nm) is a prominent RNA modification well known in noncoding RNAs and more recently also found at many mRNA internal sites. However, their function and base-resolution stoichiometry remain underexplored. Here, we investigate the transcriptome-wide effect of internal site Nm on mRNA stability. Combining nanopore sequencing with our developed machine learning method, NanoNm, we identify thousands of Nm sites on mRNAs with a single-base resolution. We observe a positive effect of FBL-mediated Nm modification on mRNA stability and expression level. Elevated FBL expression in cancer cells is associated with increased expression levels for 2'-O-methylated mRNAs of cancer pathways, implying the role of FBL in post-transcriptional regulation. Lastly, we find that FBL-mediated 2'-O-methylation connects to widespread 3' UTR shortening, a mechanism that globally increases RNA stability. Collectively, we demonstrate that FBL-mediated Nm modifications at mRNA internal sites regulate gene expression by enhancing mRNA stability.

Decoupled degradation and translation enables noise modulation by poly(A) tails

Poly(A) tails are crucial for mRNA translation and degradation, but the exact relationship between tail length and mRNA kinetics remains unclear. Here, we employ a small library of identical mRNAs that differ only in their poly(A)-tail length to examine their behavior in human embryonic kidney cells. We find that tail length strongly correlates with mRNA degradation rates but is decoupled from translation. Interestingly, an optimal tail length of ∼100 nt displays the highest translation rate, which is identical to the average endogenous tail length measured by nanopore sequencing. Furthermore, poly(A)-tail length variability-a feature of endogenous mRNAs-impacts translation efficiency but not mRNA degradation rates. Stochastic modeling combined with single-cell tracking reveals that poly(A) tails provide cells with an independent handle to tune gene expression fluctuations by decoupling mRNA degradation and translation. Together, this work contributes to the basic understanding of gene expression regulation and has potential applications in nucleic acid therapeutics.

Metabolism-dependent secondary effect of anti-MAPK cancer therapy on DNA repair

Amino acid bioavailability impacts mRNA translation in a codon-dependent manner. Here, we report that the anti-cancer MAPK inhibitors (MAPKi) decrease the intracellular concentration of aspartate and glutamate in melanoma cells. This coincides with the accumulation of ribosomes on codons corresponding to these amino acids and triggers the translation-dependent degradation of mRNAs encoding aspartate- and glutamate-rich proteins, involved in DNA metabolism such as DNA replication and repair. Consequently, cells that survive MAPKi degrade aspartate and glutamate likely to generate energy, which simultaneously decreases their requirement for amino acids due to the downregulation of aspartate- and glutamate-rich proteins involved in cell proliferation. Concomitantly, the downregulation of aspartate- and glutamate-rich proteins involved in DNA repair increases DNA damage loads. Thus, DNA repair defects, and therefore mutations, are at least in part a secondary effect of the metabolic adaptation of cells exposed to MAPKi.

Multiple strategies to improve extracellular secretion and activity of feruloyl esterase

Feruloyl esterase has a wide range of applications, but there are still problems with low enzyme yield and activity, and complex purification steps. Our previous research found Lactobacillus amylovorus feruloyl esterase could be secreted extracellular in Escherichia coli. In this study, multiple strategies were implemented to maximize the extracellular production of feruloyl esterase with improved activity in E. coli. Firstly, codon-optimized feruloyl esterase was obtained based on the preference of E. coli, resulting in 41.97 % increase in extracellular secretion. Furthermore, by cascading T7 promoters, replacing the 5' UTR, randomly mutating the N-terminal sequence, and co-expressing secretory cofactors, the extracellular secretion was increased by 36.46 %, 31.25 %, 20.66 % and 25.75 %, respectively. Moreover, the feruloyl esterase were mutated to improve the substrate affinity and activity. The catalytic efficiency of Fae-Q134T and Fae-Q198A increased by 4.62-fold and 5.42-fold. Combining above strategies, extracellular feruloyl esterase activity was increased from 2013.70 U/L to 10,349.04 U/L. These results indicated that the activity and yield of feruloyl esterase secreted by E. coli were significantly increased, which laid a foundation for its industrial application.

Current insight into the role of mRNA decay pathways in fungal pathogenesis

Pathogenic fungal species can cause superficial and mucosal infections, to potentially fatal systemic or invasive infections in humans. These infections are more common in immunocompromised or critically ill patients and have a significant morbidity and fatality rate. Fungal pathogens utilize several strategies to adapt the host environment resulting in efficient and comprehensive alterations in their cellular metabolism. Fungal virulence is regulated by several factors and post-transcriptional regulation mechanisms involving mRNA molecules are one of them. Post-transcriptional controls have emerged as critical regulatory mechanisms involved in the pathogenesis of fungal species. The untranslated upstream and downstream regions of the mRNA, as well as RNA-binding proteins, regulate morphogenesis and virulence by controlling mRNA degradation and stability. The limited number of available therapeutic drugs, the emergence of multidrug resistance, and high death rates associated with systemic fungal illnesses pose a serious risk to human health. Therefore, new antifungal treatments that specifically target mRNA pathway components can decrease fungal pathogenicity and when combined increase the effectiveness of currently available antifungal drugs. This review summarizes the mRNA degradation pathways and their role in fungal pathogenesis.

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

The purpose and ubiquity of turnover

Turnover-constant component production and destruction-is ubiquitous in biology. Turnover occurs across organisms and scales, including for RNAs, proteins, membranes, macromolecular structures, organelles, cells, hair, feathers, nails, antlers, and teeth. For many systems, turnover might seem wasteful when degraded components are often fully functional. Some components turn over with shockingly high rates and others do not turn over at all, further making this process enigmatic. However, turnover can address fundamental problems by yielding powerful properties, including regeneration, rapid repair onset, clearance of unpredictable damage and errors, maintenance of low constitutive levels of disrepair, prevention of stable hazards, and transitions. I argue that trade-offs between turnover benefits and metabolic costs, combined with constraints on turnover, determine its presence and rates across distinct contexts. I suggest that the limits of turnover help explain aging and that turnover properties and the basis for its levels underlie this fundamental component of life.

A new mRNA structure prediction based approach to identifying improved signal peptides for bone morphogenetic protein 2

BACKGROUND

Signal peptide (SP) engineering has proven able to improve production of many proteins yet is a laborious process that still relies on trial and error. mRNA structure around the translational start site is important in translation initiation and has rarely been considered in this context, with recent improvements in in silico mRNA structure potentially rendering it a useful predictive tool for SP selection. Here we attempt to create a method to systematically screen candidate signal peptide sequences in silico based on both their nucleotide and amino acid sequences. Several recently released computational tools were used to predict signal peptide activity (SignalP), localization target (DeepLoc) and predicted mRNA structure (MXFold2). The method was tested with Bone Morphogenetic Protein 2 (BMP2), an osteogenic growth factor used clinically for bone regeneration. It was hoped more effective BMP2 SPs could improve BMP2-based gene therapies and reduce the cost of recombinant BMP2 production.

RESULTS

Amino acid sequence analysis indicated 2,611 SPs from the TGF-β superfamily were predicted to function when attached to BMP2. mRNA structure prediction indicated structures at the translational start site were likely highly variable. The five sequences with the most accessible translational start sites, a codon optimized BMP2 SP variant and the well-established hIL2 SP sequence were taken forward to in vitro testing. The top five candidates showed non-significant improvements in BMP2 secretion in HEK293T cells. All showed reductions in secretion versus the native sequence in C2C12 cells, with several showing large and significant decreases. None of the tested sequences were able to increase alkaline phosphatase activity above background in C2C12s. The codon optimized control sequence and hIL2 SP showed reasonable activity in HEK293T but very poor activity in C2C12.

CONCLUSIONS

These results support the use of peptide sequence based in silico tools for basic predictions around signal peptide activity in a synthetic biology context. However, mRNA structure prediction requires improvement before it can produce reliable predictions for this application. The poor activity of the codon optimized BMP2 SP variant in C2C12 emphasizes the importance of codon choice, mRNA structure, and cellular context for SP activity.

Pervasive translation of Xrn1-sensitive unstable long noncoding RNAs in yeast

Despite being predicted to lack coding potential, cytoplasmic long noncoding (lnc)RNAs can associate with ribosomes. However, the landscape and biological relevance of lncRNA translation remain poorly studied. In yeast, cytoplasmic Xrn1-sensitive unstable transcripts (XUTs) are targeted by nonsense-mediated mRNA decay (NMD), suggesting a translation-dependent degradation process. Here, we report that XUTs are pervasively translated, which impacts their decay. We show that XUTs globally accumulate upon translation elongation inhibition, but not when initial ribosome loading is impaired. Ribo-seq confirmed ribosomes binding to XUTs and identified ribosome-associated 5'-proximal small ORFs. Mechanistically, the NMD-sensitivity of XUTs mainly depends on the 3'-untranslated region length. Finally, we show that the peptide resulting from the translation of an NMD-sensitive XUT reporter exists in NMD-competent cells. Our work highlights the role of translation in the posttranscriptional metabolism of XUTs. We propose that XUT-derived peptides could be exposed to natural selection, while NMD restricts XUT levels.

The dawn of a new Era: mRNA vaccines in colorectal cancer immunotherapy

Colorectal cancer (CRC) is a typical cancer that accounts for 10% of all new cancer cases annually and nearly 10% of all cancer deaths. Despite significant progress in current classical interventions for CRC, these traditional strategies could be invasive and with numerous adverse effects. The poor prognosis of CRC patients highlights the evident and pressing need for more efficient and targeted treatment. Novel strategies regarding mRNA vaccines for anti-tumor therapy have also been well-developed since the successful application for the prevention of COVID-19. mRNA vaccine technology won the 2023 Nobel Prize in Physiology or Medicine, signaling a new direction in human anti-cancer treatment: mRNA medicine. As a promising new immunotherapy in CRC and other multiple cancer treatments, the mRNA vaccine has higher specificity, better efficacy, and fewer side effects than traditional strategies. The present review outlines the basics of mRNA vaccines and their advantages over other vaccines and informs an available strategy for developing efficient mRNA vaccines for CRC precise treatment. In the future, more exploration of mRNA vaccines for CRC shall be attached, fostering innovation to address existing limitations.

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.

Codon usage bias of human papillomavirus type 33 and 58: A comprehensive analysis

Cervical cancer is closely linked to specific strains of human papillomavirus (HPV), notably HPV-33 and HPV-58, which exhibit a significant prevalence among women in China. Nevertheless, the codon usage bias in HPV-33 and HPV-58 is not well comprehended. The objective of this research is to analyze the codon usage patterns HPV-33 and HPV-58, pinpoint the primary factors that influence codon preference. The overall preference for codon usage in two HPV genotypes is not significant. Both HPV genotypes exhibit a preference for codons that end with A/U. The GC3 content for HPV-33 is 25.43% ± 0.35%, and for HPV-58, it is 29.44% ± 0.57%. Out of the 26 favored codons in HPV-33 and HPV-58 (relative synonymous codon usage (RSCU) > 1), 25 conclude with A/U. Principal component analysis (PCA) shows a tight clustering of the entire genome sequences of HPV-33 and HPV-58, suggesting a similarity in their RSCU preferences. Moreover, an examination of dinucleotide abundance indicated that translation selection influenced the development of a distinctive dinucleotide usage pattern in HPV-33 and HPV-58. Additionally, a combined analysis involving an effective number of codons plot, parity rule 2, and neutrality analysis demonstrated that, for HPV-33 and HPV-58, the primary determinant influencing codon usage preference is natural selection. HPV-33 and HPV-58 exhibit a restricted set of favored codons in common with humans, potentially mitigating competition for translation resources. Our discoveries could provide valuable perspectives on the evolutionary patterns and codon usage preferences of HPV-33 and HPV-58 viruses, contributing to the development and application of relevant HPV subtype vaccines.

Alter codon bias of the P. pastoris genome to overcome a bottleneck in codon optimization strategy development and improve protein expression

Apart from its role in translation, codon bias is also an important mechanism to regulate mRNA levels. The traditional frequency-based codon optimization strategy is rather efficient in organisms such as N. crassa, but much less in yeast P. pastoris which is a popular host for heterologous protein expression. This is because that unlike N. crassa, the preferred codons of P. pastoris are actually AU-rich and hence codon optimization for extremely low GC content comes with issues of pre-mature transcriptional termination or low RNA stability in spite of translational advantages. To overcome this bottleneck, we focused on three reporter genes in P. pastoris first and confirmed the great advantage of GC-prone codon optimization on mRNA levels. Then we altered the codon bias profile of P. pastoris by introducing additional rare tRNA gene copies. Prior to that we constructed IPTG-regulated tRNA species to enable chassis cells to switch between different codon bias status. As demonstrated again with reporter genes, protein yield of luc and 0788 was successfully increased by 4-5 folds in chassis cells. In summary, here we provide an alternative codon optimization strategy for genes with unsatisfactory performance under traditional codon frequency-based optimization.

Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase or decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

Cellular energy regulates mRNA degradation in a codon-specific manner

Codon optimality is a major determinant of mRNA translation and degradation rates. However, whether and through which mechanisms its effects are regulated remains poorly understood. Here we show that codon optimality associates with up to 2-fold change in mRNA stability variations between human tissues, and that its effect is attenuated in tissues with high energy metabolism and amplifies with age. Mathematical modeling and perturbation data through oxygen deprivation and ATP synthesis inhibition reveal that cellular energy variations non-uniformly alter the effect of codon usage. This new mode of codon effect regulation, independent of tRNA regulation, provides a fundamental mechanistic link between cellular energy metabolism and eukaryotic gene expression.

Translation-dependent and -independent mRNA decay occur through mutually exclusive pathways defined by ribosome density during T cell activation

mRNA translation and decay are tightly interconnected processes both in the context of mRNA quality-control pathways and for the degradation of functional mRNAs. Cotranslational mRNA degradation through codon usage, ribosome collisions, and the recruitment of specific proteins to ribosomes is an important determinant of mRNA turnover. However, the extent to which translation-dependent mRNA decay (TDD) and translation-independent mRNA decay (TID) pathways participate in the degradation of mRNAs has not been studied yet. Here we describe a comprehensive analysis of basal and signal-induced TDD and TID in mouse primary CD4+ T cells. Our results indicate that most cellular transcripts are decayed to some extent in a translation-dependent manner. Our analysis further identifies the length of untranslated regions, the density of ribosomes, and GC3 content as important determinants of TDD magnitude. Consistently, all transcripts that undergo changes in ribosome density within their coding sequence upon T cell activation display a corresponding change in their TDD level. Moreover, we reveal a dynamic modulation in the relationship between GC3 content and TDD upon T cell activation, with a reversal in the impact of GC3- and AU3-rich codons. Altogether, our data show a strong and dynamic interconnection between mRNA translation and decay in mammalian primary cells.

Attenuating ribosome load improves protein output from mRNA by limiting translation-dependent mRNA decay

Developing an effective mRNA therapeutic often requires maximizing protein output per delivered mRNA molecule. We previously found that coding sequence (CDS) design can substantially affect protein output, with mRNA variants containing more optimal codons and higher secondary structure yielding the highest protein outputs due to their slow rates of mRNA decay. Here, we demonstrate that CDS-dependent differences in translation initiation and elongation rates lead to differences in translation- and deadenylation-dependent mRNA decay rates, thus explaining the effect of CDS on mRNA half-life. Surprisingly, the most stable and highest-expressing mRNAs in our test set have modest initiation/elongation rates and ribosome loads, leading to minimal translation-dependent mRNA decay. These findings are of potential interest for optimization of protein output from therapeutic mRNAs, which may be achieved by attenuating rather than maximizing ribosome load.

UPF1 regulates mRNA stability by sensing poorly translated coding sequences

Post-transcriptional mRNA regulation shapes gene expression, yet how cis-elements and mRNA translation interface to regulate mRNA stability is poorly understood. We find that the strength of translation initiation, upstream open reading frame (uORF) content, codon optimality, AU-rich elements, microRNA binding sites, and open reading frame (ORF) length function combinatorially to regulate mRNA stability. Machine-learning analysis identifies ORF length as the most important conserved feature regulating mRNA decay. We find that Upf1 binds poorly translated and untranslated ORFs, which are associated with a higher decay rate, including mRNAs with uORFs and those with exposed ORFs after stop codons. Our study emphasizes Upf1's converging role in surveilling mRNAs with exposed ORFs that are poorly translated, such as mRNAs with long ORFs, ORF-like 3' UTRs, and mRNAs containing uORFs. We propose that Upf1 regulation of poorly/untranslated ORFs provides a unifying mechanism of surveillance in regulating mRNA stability and homeostasis in an exon-junction complex (EJC)-independent nonsense-mediated decay (NMD) pathway that we term ORF-mediated decay (OMD).

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.