Alteration of protein function by a silent polymorphism linked to tRNA abundance

A generative language model decodes contextual constraints on codon choice for mRNA design

The genetic code allows multiple synonymous codons to encode the same amino acid, creating a vast sequence space for protein-coding regions. Codon choice can impact mRNA function and protein output, a consideration newly relevant with advances in mRNA technology. Genomes preferentially use some codons, but simple optimization methods that select preferred codons miss complex contextual patterns. We present Trias, an encoder-decoder language model trained on millions of eukaryotic coding sequences. Trias learns codon usage rules directly from sequence data, integrating local and global dependencies to generate species-specific codon sequences that align with biological constraints. Without explicit training on protein expression, Trias generates sequences and scores that correlate strongly with experimental measurements of mRNA stability, ribosome load, and protein output. The model outperforms commercial codon optimization tools in generating sequences resembling high-expression codon sequence variants. By modeling codon usage in context, Trias offers a data-driven framework for synthetic mRNA design and for understanding the molecular and evolutionary principles behind codon choice.

Yeast Genome Mutagenesis With Multi-Mismatch PCR: A Rapid and Efficient Strategy for Site-Directed Mutagenesis in Saccharomyces cerevisiae

Saccharomyces cerevisiae (S. cerevisiae) provides an array of cost-effective and time-efficient methods for diverse genome modifications. Among these techniques, site-directed mutagenesis of target genes is a powerful strategy to elucidate intricate structure-function relationships and create specific mutations. While various PCR-based and CRISPR/Cas9-based methods have been developed for introducing point mutations into the S. cerevisiae genome, they often involve multiple steps. In this study, we presented a rapid and effective site-directed mutagenesis strategy using one-step multi-mismatch PCR, termed Yeast Genome Mutagenesis with Multi-mismatch PCR (YGMMP). YGMMP incorporated multiple synonymous mutations proximal to the target point mutations, along with a selection marker cassette and flanking homologous sequences, into the gene segment spanning from the desired mutation to the gene's terminus through overlap PCR. The resulting PCR product was introduced into yeast cells to facilitate the selection of target variants. As a proof of concept, we applied YGMMP to generate an ADE2 mutant. The results demonstrated that the introduction of five and nine synonymous mutations, in addition to the desired single-point mutation, yielded mutagenesis efficiencies of approximately 20% and 30%, respectively. This rapid, straightforward, and efficient method has the potential to greatly simplify site-specific modifications within the S. cerevisiae genome.

Specific tRNAs promote mRNA decay by recruiting the CCR4-NOT complex to translating ribosomes

The CCR4-NOT complex is a major regulator of eukaryotic messenger RNA (mRNA) stability. Slow decoding during translation promotes association of CCR4-NOT with ribosomes, accelerating mRNA degradation. We applied selective ribosome profiling to further investigate the determinants of CCR4-NOT recruitment to ribosomes in mammalian cells. This revealed that specific arginine codons in the P-site are strong signals for ribosomal recruitment of human CNOT3, a CCR4-NOT subunit. Cryo-electron microscopy and transfer RNA (tRNA) mutagenesis demonstrated that the D-arms of select arginine tRNAs interact with CNOT3 and promote its recruitment whereas other tRNA D-arms sterically clash with CNOT3. These effects link codon content to mRNA stability. Thus, in addition to their canonical decoding function, tRNAs directly engage regulatory complexes during translation, a mechanism we term P-site tRNA-mediated mRNA decay.

Disruption of tRNA biogenesis enhances proteostatic resilience, improves later-life health, and promotes longevity

tRNAs are evolutionarily ancient molecular decoders essential for protein translation. In eukaryotes, tRNAs and other short, noncoding RNAs are transcribed by RNA polymerase (Pol) III, an enzyme that promotes ageing in yeast, worms, and flies. Here, we show that a partial reduction in Pol III activity specifically disrupts tRNA levels. This effect is conserved across worms, flies, and mice, where computational models indicate that it impacts mRNA decoding. In all 3 species, reduced Pol III activity increases proteostatic resilience. In worms, it activates the unfolded protein response (UPR) and direct disruption of tRNA metabolism is sufficient to recapitulate this. In flies, decreasing Pol III's transcriptional initiation on tRNA genes by a loss-of-function in the TFIIIC transcription factor robustly extends lifespan, improves proteostatic resilience and recapitulates the broad-spectrum benefits to late-life health seen following partial Pol III inhibition. We provide evidence that a partial reduction in Pol III activity impacts translation, quantitatively or qualitatively, in both worms and flies, indicating a potential mode of action. Our work demonstrates a conserved and previously unappreciated role of tRNAs in animal ageing.

The Folding Pathway of ABC Transporter CFTR: Effective and Robust

De novo protein folding into a native three-dimensional structure is indispensable for biological function, is instructed by its amino acid sequence, and occurs along a vectorial trajectory. The human proteome contains thousands of membrane-spanning proteins, whose biosynthesis begins on endoplasmic reticulum-associated ribosomes. Nearly half of all membrane proteins traverse the membrane more than once, including therapeutically important protein families such as solute carriers, G-protein-coupled receptors, and ABC transporters. These mediate a variety of functions like signal transduction and solute transport and are often of vital importance for cell function and tissue homeostasis. Missense mutations in multispan membrane proteins can lead to misfolding and cause disease; an example is the ABC transporter Cystic Fibrosis Transmembrane Conductance Regulator (CFTR). Even though our understanding of multispan membrane-protein folding still is rather rudimental, the cumulative knowledge of 20 years of basic research on CFTR folding has led to development of drugs that modulate the misfolded protein. This has provided the prospect of a life without CF to the vast majority of patients. In this review we describe our understanding of the folding pathway of CFTR in cells, which is modular and tolerates many defects, making it effective and robust. We address how modulator drugs affect folding and function of CFTR, and distinguish protein stability from its folding process. Since the domain architecture of (mammalian) ABC transporters are highly conserved, we anticipate that the insights we discuss here for folding of CFTR may lay the groundwork for understanding the general rules of ABC-transporter folding.

The Effects of Codon Usage on Protein Structure and Folding

The rate of protein synthesis is slower than many folding reactions and varies depending on the synonymous codons encoding the protein sequence. Synonymous codon substitutions thus have the potential to regulate cotranslational protein folding mechanisms, and a growing number of proteins have been identified with folding mechanisms sensitive to codon usage. Typically, these proteins have complex folding pathways and kinetically stable native structures. Kinetically stable proteins may fold only once over their lifetime, and thus, codon-mediated regulation of the pioneer round of protein folding can have a lasting impact. Supporting an important role for codon usage in folding, conserved patterns of codon usage appear in homologous gene families, hinting at selection. Despite these exciting developments, there remains few experimental methods capable of quantifying translation elongation rates and cotranslational folding mechanisms in the cell, which challenges the development of a predictive understanding of how biology uses codons to regulate protein folding.

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.

Readthrough-induced misincorporated amino acid ratios guide mutant-specific therapeutic approaches for two CFTR nonsense mutations

Cystic fibrosis (CF) is a monogenic disease caused by mutations in the CF transmembrane conductance regulator (CFTR) gene. Premature termination codons (PTCs) represent ∼9% of CF mutations that typically cause severe expression defects of the CFTR anion channel. Despite the prevalence of PTCs as the underlying cause of genetic diseases, understanding the therapeutic susceptibilities of their molecular defects, both at the transcript and protein levels remains partially elucidated. Given that the molecular pathologies depend on the PTC positions in CF, multiple pharmacological interventions are required to suppress the accelerated nonsense-mediated mRNA decay (NMD), to correct the CFTR conformational defect caused by misincorporated amino acids, and to enhance the inefficient stop codon readthrough. The G418-induced readthrough outcome was previously investigated only in reporter models that mimic the impact of the local sequence context on PTC mutations in CFTR. To identify the misincorporated amino acids and their ratios for PTCs in the context of full-length CFTR readthrough, we developed an affinity purification (AP)-tandem mass spectrometry (AP-MS/MS) pipeline. We confirmed the incorporation of Cys, Arg, and Trp residues at the UGA stop codons of G542X, R1162X, and S1196X in CFTR. Notably, we observed that the Cys and Arg incorporation was favored over that of Trp into these CFTR PTCs, suggesting that the transcript sequence beyond the proximity of PTCs and/or other factors can impact the amino acid incorporation and full-length CFTR functional expression. Additionally, establishing the misincorporated amino acid ratios in the readthrough CFTR PTCs aided in maximizing the functional rescue efficiency of PTCs by optimizing CFTR modulator combinations. Collectively, our findings contribute to the understanding of molecular defects underlying various CFTR nonsense mutations and provide a foundation to refine mutation-dependent therapeutic strategies for various CF-causing nonsense mutations.

Evolving precision: rRNA expansion segment 7S modulates translation velocity and accuracy in eukaryal ribosomes

Ribosome-enhanced translational miscoding of the genetic code causes protein dysfunction and loss of cellular fitness. During evolution, open reading frame length increased, necessitating mechanisms for enhanced translation fidelity. Indeed, eukaryal ribosomes are more accurate than bacterial counterparts, despite their virtually identical, conserved active centers. During the evolution of eukaryotic organisms ribosome expansions at the rRNA and protein level occurred, which potentially increases the options for translation regulation and cotranslational events. Here we tested the hypothesis that ribosomal RNA expansions can modulate the core function of the ribosome, faithful protein synthesis. We demonstrate that a short expansion segment present in all eukaryotes' small subunit, ES7S, is crucial for accurate protein synthesis as its presence adjusts codon-specific velocities and guarantees high levels of cognate tRNA selection. Deletion of ES7S in yeast enhances mistranslation and causes protein destabilization and aggregation, dramatically reducing cellular fitness. Removal of ES7S did not alter ribosome architecture but altered the structural dynamics of inter-subunit bridges thus affecting A-tRNA selection. Exchanging the yeast ES7S sequence with the human ES7S increases accuracy whereas shortening causes the opposite effect. Our study demonstrates that ES7S provided eukaryal ribosomes with higher accuracy without perturbing the structurally conserved decoding center.

Translation velocity determines the efficacy of engineered suppressor tRNAs on pathogenic nonsense mutations

Nonsense mutations - the underlying cause of approximately 11% of all genetic diseases - prematurely terminate protein synthesis by mutating a sense codon to a premature stop or termination codon (PTC). An emerging therapeutic strategy to suppress nonsense defects is to engineer sense-codon decoding tRNAs to readthrough and restore translation at PTCs. However, the readthrough efficiency of the engineered suppressor tRNAs (sup-tRNAs) largely varies in a tissue- and sequence context-dependent manner and has not yet yielded optimal clinical efficacy for many nonsense mutations. Here, we systematically analyze the suppression efficacy at various pathogenic nonsense mutations. We discover that the translation velocity of the sequence upstream of PTCs modulates the sup-tRNA readthrough efficacy. The PTCs most refractory to suppression are embedded in a sequence context translated with an abrupt reversal of the translation speed leading to ribosomal collisions. Moreover, modeling translation velocity using Ribo-seq data can accurately predict the suppression efficacy at PTCs. These results reveal previously unknown molecular signatures contributing to genotype-phenotype relationships and treatment-response heterogeneity, and provide the framework for the development of personalized tRNA-based gene therapies.

Codon-optimization in gene therapy: promises, prospects and challenges

Codon optimization has evolved to enhance protein expression efficiency by exploiting the genetic code's redundancy, allowing for multiple codon options for a single amino acid. Initially observed in E. coli, optimal codon usage correlates with high gene expression, which has propelled applications expanding from basic research to biopharmaceuticals and vaccine development. The method is especially valuable for adjusting immune responses in gene therapies and has the potenial to create tissue-specific therapies. However, challenges persist, such as the risk of unintended effects on protein function and the complexity of evaluating optimization effectiveness. Despite these issues, codon optimization is crucial in advancing gene therapeutics. This study provides a comprehensive review of the current metrics for codon-optimization, and its practical usage in research and clinical applications, in the context of gene therapy.

Pan-cellular organelles and suborganelles-from common functions to cellular diversity?

Cell diversification is at the base of increasing multicellular organism complexity in phylogeny achieved during ontogeny. However, there are also functions common to all cells, such as cell division, cell migration, translation, endocytosis, exocytosis, etc. Here we revisit the organelles involved in such common functions, reviewing recent evidence of unexpected differences of proteins at these organelles. For instance, centrosomes or mitochondria differ significantly in their protein composition in different, sometimes closely related, cell types. This has relevance for development and disease. Particularly striking is the high amount and diversity of RNA-binding proteins at these and other organelles, which brings us to review the evidence for RNA at different organelles and suborganelles. We include a discussion about (sub)organelles involved in translation, such as the nucleolus and ribosomes, for which unexpected cell type-specific diversity has also been reported. We propose here that the heterogeneity of these organelles and compartments represents a novel mechanism for regulating cell diversity. One reason is that protein functions can be multiplied by their different contributions in distinct organelles, as also exemplified by proteins with moonlighting function. The specialized organelles still perform pan-cellular functions but in a cell type-specific mode, as discussed here for centrosomes, mitochondria, vesicles, and other organelles. These can serve as regulatory hubs for the storage and transport of specific and functionally important regulators. In this way, they can control cell differentiation, plasticity, and survival. We further include examples highlighting the relevance for disease and propose to examine organelles in many more cell types for their possible differences with functional relevance.

Identification and mechanism determination of the efflux pump subunit amrB gene mutations linked to gentamicin susceptibility in clinical Burkholderia pseudomallei from Malaysian Borneo

The bacterium Burkholderia pseudomallei is typically resistant to gentamicin but rare susceptible strains have been isolated in certain regions, such as Thailand and Sarawak, Malaysia. Recently, several amino acid substitutions have been reported in the amrB gene (a subunit of the amrAB-oprA efflux pump gene) that confer gentamicin susceptibility. However, information regarding the mechanism of the substitutions conferring the susceptibility is lacking. To understand the mechanism of amino acid substitution that confers susceptibility, this study identifies the corresponding mutations in clinical gentamicin-susceptible B. pseudomallei isolates from the Malaysian Borneo (n = 46; Sarawak: 5; Sabah: 41). Three phenotypically confirmed gentamicin-susceptible (GENs) strains from Sarawak, Malaysia, were screened for mutations in the amrB gene using gene sequences of gentamicin-resistant (GENr) strains (QEH 56, QEH 57, QEH20, and QEH26) and publicly available sequences (AF072887.1 and BX571965.1) as the comparator. The effect of missense mutations on the stability of the AmrB protein was determined by calculating the average energy change value (ΔΔG). Mutagenesis analysis identified a polymorphism-associated mutation, g.1056 T > G, a possible susceptible-associated in-frame deletion, Delta V412, and a previously confirmed susceptible-associated amino acid substitution, T368R, in each of the three GENs isolates. The contribution of Delta V412 needs further confirmation by experimental mutagenesis analysis. The mechanism by which T368R confers susceptibility, as elucidated by in silico mutagenesis analysis using AmrB-modeled protein structures, is proposed to be due to the location of T368R in a highly conserved region, rather than destabilization of the AmrB protein structure.

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.

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.

Implications of tRNA abundance on translation elongation across bovine tissues

Introduction: Translation is a crucial stage of gene expression. It may also act as an additional layer of regulation that plays an important role in gene expression and function. Highly expressed genes are believed to be codon-biased to support increased protein production, in which quickly translated codons correspond to highly abundant tRNAs. Synonymous SNPs, considered to be silent due to the degeneracy of the genetic code, may shift protein abundance and function through alterations in translational efficiency and suboptimal pairing to lowly abundant tRNAs. Methods: Here, we applied Quantitative Mature tRNA sequencing (QuantM-tRNAseq) and ribosome profiling across bovine tissues in order to investigate the relationship between tRNA expression and slowed translation. Results: Moreover, we have identified genes modulated at transcriptional and/or translational levels underlying tissue-specific biological processes. We have also successfully defined pausing sites that depict the regulatory information encoded within the open reading frame of transcripts, which could be related to translation rate and facilitate proper protein folding. This work offers an atlas of distinctive pausing sites across three bovine tissues, which provides an opportunity to predict codon optimality and understand tissue-specific mechanisms of regulating protein synthesis.

Integrated multiplexed assays of variant effect reveal cis-regulatory 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 and 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.

Seryl-tRNA synthetase promotes translational readthrough by mRNA binding and involvement of the selenocysteine incorporation machinery

Translational readthrough of UGA stop codons by selenocysteine-specific tRNA (tRNASec) enables the synthesis of selenoproteins. Seryl-tRNA synthetase (SerRS) charges tRNASec with serine, which is modified into selenocysteine and delivered to the ribosome by a designated elongation factor (eEFSec in eukaryotes). Here we found that components of the human selenocysteine incorporation machinery (SerRS, tRNASec, and eEFSec) also increased translational readthrough of non-selenocysteine genes, including VEGFA, to create C-terminally extended isoforms. SerRS recognizes target mRNAs through a stem-loop structure that resembles the variable loop of its cognate tRNAs. This function of SerRS depends on both its enzymatic activity and a vertebrate-specific domain. Through eCLIP-seq, we identified additional SerRS-interacting mRNAs as potential readthrough genes. Moreover, SerRS overexpression was sufficient to reverse premature termination caused by a pathogenic nonsense mutation. Our findings expand the repertoire of selenoprotein biosynthesis machinery and suggest an avenue for therapeutic targeting of nonsense mutations using endogenous factors.

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.

Predicting functional consequences of SNPs on mRNA translation via machine learning

The functional impact of single nucleotide polymorphisms (SNPs) on translation has yet to be considered when prioritizing disease-causing SNPs from genome-wide association studies (GWAS). Here we apply machine learning models to genome-wide ribosome profiling data to predict SNP function by forecasting ribosome collisions during mRNA translation. SNPs causing remarkable ribosome occupancy changes are named RibOc-SNPs (Ribosome-Occupancy-SNPs). We found that disease-related SNPs tend to cause notable changes in ribosome occupancy, suggesting translational regulation as an essential pathogenesis step. Nucleotide conversions, such as 'G → T', 'T → G' and 'C → A', are enriched in RibOc-SNPs, with the most significant impact on ribosome occupancy, while 'A → G' (or 'A→ I' RNA editing) and 'G → A' are less deterministic. Among amino acid conversions, 'Glu → stop (codon)' shows the most significant enrichment in RibOc-SNPs. Interestingly, there is selection pressure on stop codons with a lower collision likelihood. RibOc-SNPs are enriched at the 5'-coding sequence regions, implying hot spots of translation initiation regulation. Strikingly, ∼22.1% of the RibOc-SNPs lead to opposite changes in ribosome occupancy on alternative transcript isoforms, suggesting that SNPs can amplify the differences between splicing isoforms by oppositely regulating their translation efficiency.

Engineered tRNAs suppress nonsense mutations in cells and in vivo

Nonsense mutations are the underlying cause of approximately 11% of all inherited genetic diseases1. Nonsense mutations convert a sense codon that is decoded by tRNA into a premature termination codon (PTC), resulting in an abrupt termination of translation. One strategy to suppress nonsense mutations is to use natural tRNAs with altered anticodons to base-pair to the newly emerged PTC and promote translation2-7. However, tRNA-based gene therapy has not yielded an optimal combination of clinical efficacy and safety and there is presently no treatment for individuals with nonsense mutations. Here we introduce a strategy based on altering native tRNAs into  efficient suppressor tRNAs (sup-tRNAs) by individually fine-tuning their sequence to the physico-chemical properties of the amino acid that they carry. Intravenous and intratracheal lipid nanoparticle (LNP) administration of sup-tRNA in mice restored the production of functional proteins with nonsense mutations. LNP-sup-tRNA formulations caused no discernible readthrough at endogenous native stop codons, as determined by ribosome profiling. At clinically important PTCs in the cystic fibrosis transmembrane conductance regulator gene (CFTR), the sup-tRNAs re-established expression and function in cell systems and patient-derived nasal epithelia and restored airway volume homeostasis. These results provide a framework for the development of tRNA-based therapies with a high molecular safety profile and high efficacy in targeted PTC suppression.

Monosomes buffer translational stress to allow for active ribosome elongation

Introduction: The synthesis of proteins is a fundamental process in the life-span of all cells. The activation of ribosomes on transcripts is the starting signal for elongation and, in turn, the translation of an mRNA. Thereby, most mRNAs circulate between single (monosomes) and multi ribosomal particles (polysomes), a process that defines their translational activity. The interplay between monosomes and polysomes is thought to crucially impact translation rate. How monosomes and polysomes are balanced during stress remains, however, elusive. Methods: Here, we set out to investigate the monosome and polysome levels as well as their kinetics under different translational stress conditions including mTOR inhibition, downregulation of the eukaryotic elongation factor 2 (eEF2) and amino acid depletion. Results: By using a timed ribosome runoff approach in combination with polysome profiling, we found that the used translational stressors show very distinct effects on translation. However, they all had in common that the activity of monosomes was preferentially affected. This adaptation seems to be needed for sufficient translation elongation. Even under harsh conditions such as amino acid starvation, we detected active polysomes while monosomes were mostly inactive. Hence, it is plausible that cells compensate the reduced availability of essential factors during stress by adapting the levels of active monosomes to favor sufficient elongation. Discussion: These results suggest that monosome and polysome levels are balanced under stress conditions. Together, our data argue for the existence of translational plasticity that ensure sufficient protein synthesis under stress conditions, a process that is necessary for cell survival and recovery.

Influence of substitution model selection on protein phylogenetic tree reconstruction

Probabilistic phylogenetic tree reconstruction is traditionally performed under a best-fitting substitution model of molecular evolution previously selected according to diverse statistical criteria. Interestingly, some recent studies proposed that this procedure is unnecessary for phylogenetic tree reconstruction leading to a debate in the field. In contrast to DNA sequences, phylogenetic tree reconstruction from protein sequences is traditionally based on empirical exchangeability matrices that can differ among taxonomic groups and protein families. Considering this aspect, here we investigated the influence of selecting a substitution model of protein evolution on phylogenetic tree reconstruction by the analyses of real and simulated data. We found that phylogenetic tree reconstructions based on a selected best-fitting substitution model of protein evolution are the most accurate, in terms of topology and branch lengths, compared with those derived from substitution models with amino acid replacement matrices far from the selected best-fitting model, especially when the data has large genetic diversity. Indeed, we found that substitution models with similar amino acid replacement matrices produce similar reconstructed phylogenetic trees, suggesting the use of substitution models as similar as possible to a selected best-fitting model when the latter cannot be used. Therefore, we recommend the use of the traditional protocol of selection among substitution models of evolution for protein phylogenetic tree reconstruction.

satmut_utils: a simulation and variant calling package for multiplexed assays of variant effect

The impact of millions of individual genetic variants on molecular phenotypes in coding sequences remains unknown. Multiplexed assays of variant effect (MAVEs) are scalable methods to annotate relevant variants, but existing software lacks standardization, requires cumbersome configuration, and does not scale to large targets. We present satmut_utils as a flexible solution for simulation and variant quantification. We then benchmark MAVE software using simulated and real MAVE data. We finally determine mRNA abundance for thousands of cystathionine beta-synthase variants using two experimental methods. The satmut_utils package enables high-performance analysis of MAVEs and reveals the capability of variants to alter mRNA abundance.

Synonymous and Nonsynonymous Substitutions in Dictyostelium discoideum Ammonium Transporter amtA Are Necessary for Functional Complementation in Saccharomyces cerevisiae

Ammonium transporters are present in all three domains of life. They have undergone extensive horizontal gene transfer (HGT), gene duplication, and functional diversification and therefore offer an excellent paradigm to study protein evolution. We attempted to complement a mep1Δmep2Δmep3Δ strain of Saccharomyces cerevisiae (triple-deletion strain), which otherwise cannot grow on ammonium as a sole nitrogen source at concentrations of <3 mM, with amtA of Dictyostelium discoideum, an orthologue of S. cerevisiae MEP2. We observed that amtA did not complement the triple-deletion strain of S. cerevisiae for growth on low-ammonium medium. We isolated two mutant derivatives of amtA (amtA M1 and amtA M2) from a PCR-generated mutant plasmid library that complemented the triple-deletion strain of S. cerevisiae. amtA M1 bears three nonsynonymous and two synonymous substitutions, which are necessary for its functionality. amtA M2 bears two nonsynonymous substitutions and one synonymous substitution, all of which are necessary for functionality. Interestingly, AmtA M1 transports ammonium but does not confer methylamine toxicity, while AmtA M2 transports ammonium and confers methylamine toxicity, demonstrating functional diversification. Preliminary biochemical analyses indicated that the mutants differ in their conformations as well as their mechanisms of ammonium transport. These intriguing results clearly point out that protein evolution cannot be fathomed by studying nonsynonymous and synonymous substitutions in isolation. The above-described observations have significant implications for various facets of biological processes and are discussed in detail. IMPORTANCE Functional diversification following gene duplication is one of the major driving forces of protein evolution. While the role of nonsynonymous substitutions in the functional diversification of proteins is well recognized, knowledge of the role of synonymous substitutions in protein evolution is in its infancy. Using functional complementation, we isolated two functional alleles of the D. discoideum ammonium transporter gene (amtA), which otherwise does not function in S. cerevisiae as an ammonium transporters. One of them is an ammonium transporter, while the other is an ammonium transporter that also confers methylammonium (ammonium analogue) toxicity, suggesting functional diversification. Surprisingly, both alleles require a combination of synonymous and nonsynonymous substitutions for their functionality. These results bring out a hitherto-unknown pathway of protein evolution and pave the way for not only understanding protein evolution but also interpreting single nucleotide polymorphisms (SNPs).

Not1 and Not4 inversely determine mRNA solubility that sets the dynamics of co-translational events

BACKGROUND

The Ccr4-Not complex is mostly known as the major eukaryotic deadenylase. However, several studies have uncovered roles of the complex, in particular of the Not subunits, unrelated to deadenylation and relevant for translation. In particular, the existence of Not condensates that regulate translation elongation dynamics has been reported. Typical studies that evaluate translation efficiency rely on soluble extracts obtained after the disruption of cells and ribosome profiling. Yet cellular mRNAs in condensates can be actively translated and may not be present in such extracts.

RESULTS

In this work, by analyzing soluble and insoluble mRNA decay intermediates in yeast, we determine that insoluble mRNAs are enriched for ribosomes dwelling at non-optimal codons compared to soluble mRNAs. mRNA decay is higher for soluble RNAs, but the proportion of co-translational degradation relative to the overall mRNA decay is higher for insoluble mRNAs. We show that depletion of Not1 and Not4 inversely impacts mRNA solubilities and, for soluble mRNAs, ribosome dwelling according to codon optimality. Depletion of Not4 solubilizes mRNAs with lower non-optimal codon content and higher expression that are rendered insoluble by Not1 depletion. By contrast, depletion of Not1 solubilizes mitochondrial mRNAs, which are rendered insoluble upon Not4 depletion.

CONCLUSIONS

Our results reveal that mRNA solubility defines the dynamics of co-translation events and is oppositely regulated by Not1 and Not4, a mechanism that we additionally determine may already be set by Not1 promoter association in the nucleus.

Ending a bad start: Triggers and mechanisms of co-translational protein degradation

Proteins are versatile molecular machines that control and execute virtually all cellular processes. They are synthesized in a multilayered process requiring transfer of information from DNA to RNA and finally into polypeptide, with many opportunities for error. In addition, nascent proteins must successfully navigate a complex folding-energy landscape, in which their functional native state represents one of many possible outcomes. Consequently, newly synthesized proteins are at increased risk of misfolding and toxic aggregation. To maintain proteostasis-the state of proteome balance-cells employ a plethora of molecular chaperones that guide proteins along a productive folding pathway and quality control factors that direct misfolded species for degradation. Achieving the correct balance between folding and degradation therefore represents a fundamental task for the proteostasis network. While many chaperones act co-translationally, protein quality control is generally considered to be a post-translational process, as the majority of proteins will only achieve their final native state once translation is completed. Nevertheless, it has been observed that proteins can be ubiquitinated during synthesis. The extent and the relevance of co-translational protein degradation, as well as the underlying molecular mechanisms, remain areas of open investigation. Recent studies made seminal advances in elucidating ribosome-associated quality control processes, and how their loss of function can lead to proteostasis failure and disease. Here, we discuss current understanding of the situations leading to the marking of nascent proteins for degradation before synthesis is completed, and the emerging quality controls pathways engaged in this task in eukaryotic cells. We also highlight the methods used to study co-translational quality control.

Features of CFTR mRNA and implications for therapeutics development

Cystic fibrosis (CF) is an autosomal recessive disease impacting ∼100,000 people worldwide. This lethal disorder is caused by mutation of the CF transmembrane conductance regulator (CFTR) gene, which encodes an ATP-binding cassette-class C protein. More than 2,100 variants have been identified throughout the length of CFTR. These defects confer differing levels of severity in mRNA and/or protein synthesis, folding, gating, and turnover. Drug discovery efforts have resulted in recent development of modulator therapies that improve clinical outcomes for people living with CF. However, a significant portion of the CF population has demonstrated either no response and/or adverse reactions to small molecules. Additional therapeutic options are needed to restore underlying genetic defects for all patients, particularly individuals carrying rare or refractory CFTR variants. Concerted focus has been placed on rescuing variants that encode truncated CFTR protein, which also harbor abnormalities in mRNA synthesis and stability. The current mini-review provides an overview of CFTR mRNA features known to elicit functional consequences on final protein conformation and function, including considerations for RNA-directed therapies under investigation. Alternative exon usage in the 5'-untranslated region, polypyrimidine tracts, and other sequence elements that influence splicing are discussed. Additionally, we describe mechanisms of CFTR mRNA decay and post-transcriptional regulation mediated through interactions with the 3'-untranslated region (e.g. poly-uracil sequences, microRNAs). Contributions of synonymous single nucleotide polymorphisms to CFTR transcript utilization are also examined. Comprehensive understanding of CFTR RNA biology will be imperative for optimizing future therapeutic endeavors intended to address presently untreatable forms of CF.

Novel compound heterozygous synonymous and missense variants in the MYO7A gene identified by next-generation sequencing in a Chinese family with nonsyndromic hearing loss

BACKGROUND

Variants in the MYO7A gene are increasingly identified among patients suffering from Usher syndrome type 1B (USH1B). However, such mutations are less commonly detected among patients suffering from nonsyndromic hearing loss (NSHL), including autosomal recessive deafness (DFNB2) and autosomal dominant deafness (DFNA11). This research attempts to clarify the genetic base of DFNB2 in a Chinese family and determine the pathogenicity of the identified mutations.

METHOD

Targeted next-generation sequencing (TGS) of 127 known deafness genes was performed for the 14-year-old proband. Then, Sanger sequencing was performed on the available family members. A minigene splicing assay was performed to verify the impact of the novel MYO7A synonymous variant. After performing targeted next-generation sequencing (TGS) of 127 existing hearing loss-related genes in a 14-year-old proband, Sanger sequencing was carried out on the available family members. Then, to confirm the influence of the novel MYO7A synonymous variants, a minigene splicing assay was performed.

RESULTS

Two heteroallelic mutants of MYO7A (NM_000260.3) were identified: a maternally inherited synonymous variant c.2904G > A (p.Glu968=) in exon 23 and a paternally inherited missense variant c.5994G > T (p.Trp1998Cys) in exon 44. The in vitro minigene expression indicated that c.2904G > A may result in skipping of exon 23 resulting in a truncated protein.

CONCLUSIONS

We reported a novel missense (c.5994G > T) and identified, for the first time, a novel pathogenic synonymous (c.2904G > A) variant within MYO7A in a patient with DFNB2. These findings enrich our understanding of the MYO7A variant spectrum of DFNB2 and can contribute to accurate genetic counseling and diagnosis of NSHL patients.

Consequences of Substitution Model Selection on Protein Ancestral Sequence Reconstruction

The selection of the best-fitting substitution model of molecular evolution is a traditional step for phylogenetic inferences, including ancestral sequence reconstruction (ASR). However, a few recent studies suggested that applying this procedure does not affect the accuracy of phylogenetic tree reconstruction. Here, we revisited this debate topic by analyzing the influence of selection among substitution models of protein evolution, with focus on exchangeability matrices, on the accuracy of ASR using simulated and real data. We found that the selected best-fitting substitution model produces the most accurate ancestral sequences, especially if the data present large genetic diversity. Indeed, ancestral sequences reconstructed under substitution models with similar exchangeability matrices were similar, suggesting that if the selected best-fitting model cannot be used for the reconstruction, applying a model similar to the selected one is preferred. We conclude that selecting among substitution models of protein evolution is recommended for reconstructing accurate ancestral sequences.

Polymorphisms in common antihypertensive targets: Pharmacogenomic implications for the treatment of cardiovascular disease

The idea of personalized medicine came to fruition with sequencing the human genome; however, aside from a few cases, the genetic revolution has yet to materialize. Cardiovascular diseases are the leading cause of death globally, and hypertension is a common prelude to nearly all cardiovascular diseases. Thus, hypertension is an ideal candidate disease to apply tenants of personalized medicine to lessen cardiovascular disease. Herein is a survey that visually depicts the polymorphisms in the top eight antihypertensive targets. Although there are numerous genome-wide association studies regarding cardiovascular disease, few studies look at the effects of receptor polymorphisms on drug treatment. With 17,000+ polymorphisms in the combined target proteins examined, it is expected that some of the clinical variability in the treatment of hypertension is due to polymorphisms in the drug targets. Recent advances in techniques and technology, such as high throughput examination of single mutations, structure prediction, computational power for modeling, and CRISPR models of point mutations, allow for a relatively rapid and comprehensive examination of the effects of known and future polymorphisms on drug affinity and effects. As hypertension is easy to measure and has a plethora of clinically viable ligands, hypertension makes an excellent disease to study pharmacogenomics in the lab and the clinic. If the promises of personalized medicine are to materialize, a concerted effort to examine the effects polymorphisms have on drugs is required. A clinician with the knowledge of a patient's genotype can then prescribe drugs that are optimal for treating that specific patient.

The Ramp Atlas: facilitating tissue and cell-specific ramp sequence analyses through an intuitive web interface

Ramp sequences occur when the average translational efficiency of codons near the 5′ end of highly expressed genes is significantly lower than the rest of the gene sequence, which counterintuitively increases translational efficiency by decreasing downstream ribosomal collisions. Here, we show that the relative codon adaptiveness within different tissues changes the existence of a ramp sequence without altering the underlying genetic code. We present the first comprehensive analysis of tissue and cell type-specific ramp sequences and report 3108 genes with ramp sequences that change between tissues and cell types, which corresponds with increased gene expression within those tissues and cells. The Ramp Atlas (https://ramps.byu.edu/) allows researchers to query precomputed ramp sequences in 18 388 genes across 62 tissues and 66 cell types and calculate tissue-specific ramp sequences from user-uploaded FASTA files through an intuitive web interface. We used The Ramp Atlas to identify seven SARS-CoV-2 genes and seven human SARS-CoV-2 entry factor genes with tissue-specific ramp sequences that may help explain viral proliferation within those tissues. We anticipate that The Ramp Atlas will facilitate personalized and creative tissue-specific ramp sequence analyses for both human and viral genes that will increase our ability to utilize this often-overlooked regulatory region.

High frequency of complex CFTR alleles associated with c.1521_1523delCTT (F508del) in Russian cystic fibrosis patients

Cystic fibrosis (CF, MIM# 219,700) is an autosomal recessive disease caused by pathogenic variants within the CFTR gene. It was shown that genetic variants located in cis can affect disease severity or treatment response because of additive or epistatic effects. Studies on the prevalence of complex alleles in Russian CF patients have just begun.

Aim

To evaluate frequencies and genetic background of complex alleles carrying c.1521_1523delCTT (F508del) and c.1399C>T (L467F), c.2562T>G (T854=) or c.4389G>A (Q1463=) in cis; to determine clinical consequences of complex allele c.[1399C>T;1521_1523delCTT] ([L467;F508del]) in Russian CF patients.

Methods

Sequencing of coding regions of CFTR gene and analysis of polymorphic markers in CF patients carrying F508del variant. Comparing of clinical features in two groups patients having genotypes [L467F;F508del];[F508del] (group 1) and [F508del];[F508del] (group 2).

Results

Frequency of [L467F;F508del] allele linked to 2–2–21–6–17–13 haplotype was 4.42%, of [F508del;T854=;Q1463=] allele linked to haplotype 1–2–21–6–17–13 – 2.2% in F508del chromosomes. No differences in disease severity in patients carrying complex allele [L467F;F508del] and patients homozygous for F508del was found.

Conclusion

The frequency of complex alleles associated with F508del was at least 6.6% in Russian CF patients, which should be taken into account for the decision on optimal treatment options with CFTR modulators.

WT1, NR0B1, NR5A1, LHX9, ZFP92, ZNF275, INSL3, and NRIP1 Genetic Variants in Patients with Premature Ovarian Insufficiency in a Mexican Cohort

Premature ovarian insufficiency (POI) is one of the main causes of female premature infertility. POI is a genetically heterogeneous disorder with a complex molecular etiology; as such, the genetic causes remain unknown in the majority of patients. Therefore, this study aimed to identify mutations and characterize the associated molecular contribution of gonadogenesis-determinant genes to POI. Genomic assays, including PCR-SSCP and Sanger sequencing, followed by in silico analyses were used to investigate the underpinnings of ovarian deficiency in 11 women affected by POI. Large deletions and nucleotide insertions and duplications were excluded by PCR. Thirteen genetic variants were identified in the WT1 (c.213G>T, c.609T>C, c.873A>G, c.1122G>A), NR0B1 (c.353C>T, c.425G>A), NR5A1 (c.437G>C, IVS4-20C>T), LHX9 (IVS2-12G>C, IVS3+13C>T, c.741T>C), ZNF275 (c.969C>T), and NRIP1 (c.3403C>T) genes. Seven novel genetic variants and five unpublished substitutions were identified. No genetic aberrations were detected in the ZFP92 and INSL3 genes. Each variant was genotyped using PCR-SSCP in 100 POI-free subjects, and their allelic frequencies were similar to the patients. These analyses indicated that allelic variation in the WT1, NR0B1, NR5A1, LHX9, ZFP92, ZNF275, INSL3, and NRIP1 genes may be a non-disease-causing change or may not contribute significantly to the genetics underlying POI disorders. Findings support the polygenic nature of this clinical disorder, with the SNVs identified representing only a probable contribution to the variability of the human genome.

A splice site mutation in the TSEN2 causes a new syndrome with craniofacial and central nervous system malformations, and atypical hemolytic uremic syndrome

Recessive mutations in the genes encoding the four subunits of the tRNA splicing endonuclease complex (TSEN54, TSEN34, TSEN15, and TSEN2) cause various forms of pontocerebellar hypoplasia, a disorder characterized by hypoplasia of the cerebellum and the pons, microcephaly, dysmorphisms, and other variable clinical features. Here, we report an intronic recessive founder variant in the gene TSEN2 that results in abnormal splicing of the mRNA of this gene, in six individuals from four consanguineous families affected with microcephaly, multiple craniofacial malformations, radiological abnormalities of the central nervous system, and cognitive retardation of variable severity. Remarkably, unlike patients with previously described mutations in the components of the TSEN complex, all the individuals that we report developed atypical hemolytic uremic syndrome (aHUS) with thrombotic microangiopathy, microangiopathic hemolytic anemia, thrombocytopenia, proteinuria, severe hypertension, and end-stage kidney disease (ESKD) early in life. Bulk RNA sequencing of peripheral blood cells of four affected individuals revealed abnormal tRNA transcripts, indicating an alteration of the tRNA biogenesis. Morpholino-mediated skipping of exon 10 of tsen2 in zebrafish produced phenotypes similar to human patients. Thus, we have identified a novel syndrome accompanied by aHUS suggesting the existence of a link between tRNA biology and vascular endothelium homeostasis, which we propose to name with the acronym TRACK syndrome (TSEN2 Related Atypical hemolytic uremic syndrome, Craniofacial malformations, Kidney failure).

mRNAs sequestered in stress granules recover nearly completely for translation

Stress granules (SGs) are membrane-less condensates composed of RNA and protein that assemble in response to stress stimuli and disassemble when stress is lifted. Both assembly and disassembly are tightly controlled processes, yet, it remains elusive whether mRNAs in SGs completely recover for translation following stress relief. Using RNA-seq of translating fractions in human cell line, we found that higher fraction of the m⁶A-modified mRNAs recovered for translation compared to unmodified mRNAs, i.e. 95% vs 84%, respectively. Considering structural mRNA analysis, we found that the m⁶A modification enhances structuring at nucleotides in its close vicinity. Our results suggest that SG-sequestered mRNAs disassemble nearly completely from SGs and the m⁶A modification may display some advantage to the mRNAs in their recovery for translation likely by m⁶A-driven structural stabilization.

Availability of Arg, but Not tRNA, Is a Rate-Limiting Factor for Intracellular Arginylation

Protein arginylation, mediated by arginyltransferase ATE1, is a posttranslational modification of emerging biological importance that consists of transfer of the amino acid Arg from tRNA to protein and peptide targets. ATE1 can bind tRNA and exhibits specificity toward particular tRNA types, but its dependence on the availability of the major components of the arginylation reaction has never been explored. Here we investigated key intracellular factors that can potentially regulate arginylation in vivo, including several tRNA types that show strong binding to ATE1, as well as availability of free Arg, in an attempt to identify intracellular rate limiting steps for this enzyme. Our results demonstrate that, while modulation of tRNA levels in cells does not lead to any changes in intracellular arginylation efficiency, availability of free Arg is a potentially rate-limiting factor that facilitates arginylation if added to the cultured cells. Our results broadly outline global pathways that may be involved in the regulation of arginylation in vivo.

Coordination of -1 programmed ribosomal frameshifting by transcript and nascent chain features revealed by deep mutational scanning

Programmed ribosomal frameshifting (PRF) is a translational recoding mechanism that enables the synthesis of multiple polypeptides from a single transcript. During translation of the alphavirus structural polyprotein, the efficiency of -1PRF is coordinated by a 'slippery' sequence in the transcript, an adjacent RNA stem-loop, and a conformational transition in the nascent polypeptide chain. To characterize each of these effectors, we measured the effects of 4530 mutations on -1PRF by deep mutational scanning. While most mutations within the slip-site and stem-loop reduce the efficiency of -1PRF, the effects of mutations upstream of the slip-site are far more variable. We identify several regions where modifications of the amino acid sequence of the nascent polypeptide impact the efficiency of -1PRF. Molecular dynamics simulations of polyprotein biogenesis suggest the effects of these mutations primarily arise from their impacts on the mechanical forces that are generated by the translocon-mediated cotranslational folding of the nascent polypeptide chain. Finally, we provide evidence suggesting that the coupling between cotranslational folding and -1PRF depends on the translation kinetics upstream of the slip-site. These findings demonstrate how -1PRF is coordinated by features within both the transcript and nascent chain.

Detection of a novel single nucleotide polymorphism in IGF2 gene with a negative impact on egg production and body weight in Japanese quail (Coturnix japonica)

Background

Insulin-like growth factor 2 (IGF2) is one of three hormones that share high structural similarity to insulin. It is involved in several insulin-like growth-regulating and mitogenic activities. This study was conducted to genotype the coding regions of the IGF2 gene in Japanese quail (Coturnix japonica) using PCR-SSCP-sequencing, and to assess the possible association of the polymorphism of these regions with the main egg production traits. A total of 240 female birds with an equal number of three Japanese quail populations (Brown or BR, Black or BL, and White or WT) were included in this study.

Results

All the genotyped regions exerted no heterogeneity in their sequences with one exception detected in the exon 2. In this locus, a novel single nucleotide polymorphism (SNP) was detected in which “A” was substituted with “G” at 81 position with a silent effect (p.F79=SNP) on IGF2 protein. Association analyses indicated a significant (P < 0.05) relation of this SNP with egg number (EN) and bird weight (BW) in the analyzed populations, in which the birds with AG genotype had lower EN and BW than those with AA genotype. The p.F79=SNP was largely detected in the WT line than the other two lines.

Conclusion

The detected p.F79=SNP has a highly negative effect on EN and BW in Japanese quail. Thus, the implementation of the variations of the IGF2 gene can be a useful marker in the marker-assisted selection of Japanese quail. This is the first report to describe IGF2 gene variations in Japanese quail, which strongly suggests raising the birds from the BR line with AA genotype when breeders desire to raise Japanese quail for large-scale egg production.

mRNA codon optimization with quantum computers

Reverse translation of polypeptide sequences to expressible mRNA constructs is a NP-hard combinatorial optimization problem. Each amino acid in the protein sequence can be represented by as many as six codons, and the process of selecting the combination that maximizes probability of expression is termed codon optimization. This work investigates the potential impact of leveraging quantum computing technology for codon optimization. A Quantum Annealer (QA) is compared to a standard genetic algorithm (GA) programmed with the same objective function. The QA is found to be competitive in identifying optimal solutions. The utility of gate-based systems is also evaluated using a simulator resulting in the finding that while current generations of devices lack the hardware requirements, in terms of both qubit count and connectivity, to solve realistic problems, future generation devices may be highly efficient.

Assembly defects of human tRNA splicing endonuclease contribute to impaired pre-tRNA processing in pontocerebellar hypoplasia

Introns of human transfer RNA precursors (pre-tRNAs) are excised by the tRNA splicing endonuclease TSEN in complex with the RNA kinase CLP1. Mutations in TSEN/CLP1 occur in patients with pontocerebellar hypoplasia (PCH), however, their role in the disease is unclear. Here, we show that intron excision is catalyzed by tetrameric TSEN assembled from inactive heterodimers independently of CLP1. Splice site recognition involves the mature domain and the anticodon-intron base pair of pre-tRNAs. The 2.1-Å resolution X-ray crystal structure of a TSEN15-34 heterodimer and differential scanning fluorimetry analyses show that PCH mutations cause thermal destabilization. While endonuclease activity in recombinant mutant TSEN is unaltered, we observe assembly defects and reduced pre-tRNA cleavage activity resulting in an imbalanced pre-tRNA pool in PCH patient-derived fibroblasts. Our work defines the molecular principles of intron excision in humans and provides evidence that modulation of TSEN stability may contribute to PCH phenotypes.

tRNA overexpression rescues peripheral neuropathy caused by mutations in tRNA synthetase

[Figure: see text].

Not4 and Not5 modulate translation elongation by Rps7A ubiquitination, Rli1 moonlighting, and condensates that exclude eIF5A

In this work, we show that Not4 and Not5 from the Ccr4-Not complex modulate translation elongation dynamics and change ribosome A-site dwelling occupancy in a codon-dependent fashion. These codon-specific changes in not5Δ cells are very robust and independent of codon position within the mRNA, the overall mRNA codon composition, or changes of mRNA expression levels. They inversely correlate with codon-specific changes in cells depleted for eIF5A and positively correlate with those in cells depleted for ribosome-recycling factor Rli1. Not5 resides in punctate loci, co-purifies with ribosomes and Rli1, but not with eIF5A, and limits mRNA solubility. Overexpression of wild-type or non-complementing Rli1 and loss of Rps7A ubiquitination enable Not4 E3 ligase-dependent translation of polyarginine stretches. We propose that Not4 and Not5 modulate translation elongation dynamics to produce a soluble proteome by Rps7A ubiquitination, dynamic condensates that limit mRNA solubility and exclude eIF5A, and a moonlighting function of Rli1.

A synonymous variant in MYO15A enriched in the Ashkenazi Jewish population causes autosomal recessive hearing loss due to abnormal splicing

Nonsyndromic hearing loss is genetically heterogeneous. Despite comprehensive genetic testing, many cases remain unsolved because the clinical significance of identified variants is uncertain or because biallelic pathogenic variants are not identified for presumed autosomal recessive cases. Common synonymous variants are often disregarded. Determining the pathogenicity of synonymous variants may improve genetic diagnosis. We report a synonymous variant c.9861 C > T/p.(Gly3287=) in MYO15A in homozygosity or compound heterozygosity with another pathogenic or likely pathogenic MYO15A variant in 10 unrelated families with nonsyndromic sensorineural hearing loss. Biallelic variants in MYO15A were identified in 21 affected and were absent in 22 unaffected siblings. A mini-gene assay confirms that the synonymous variant leads to abnormal splicing. The variant is enriched in the Ashkenazi Jewish population. Individuals carrying biallelic variants involving c.9861 C > T often exhibit progressive post-lingual hearing loss distinct from the congenital profound deafness typically associated with biallelic loss-of-function MYO15A variants. This study establishes the pathogenicity of the c.9861 C > T variant in MYO15A and expands the phenotypic spectrum of MYO15A-related hearing loss. Our work also highlights the importance of multicenter collaboration and data sharing to establish the pathogenicity of a relatively common synonymous variant for improved diagnosis and management of hearing loss.

The genomic loci of specific human tRNA genes exhibit ageing-related DNA hypermethylation

The epigenome has been shown to deteriorate with age, potentially impacting on ageing-related disease. tRNA, while arising from only ˜46 kb (<0.002% genome), is the second most abundant cellular transcript. tRNAs also control metabolic processes known to affect ageing, through core translational and additional regulatory roles. Here, we interrogate the DNA methylation state of the genomic loci of human tRNA. We identify a genomic enrichment for age-related DNA hypermethylation at tRNA loci. Analysis in 4,350 MeDIP-seq peripheral-blood DNA methylomes (16-82 years), identifies 44 and 21 hypermethylating specific tRNAs at study-and genome-wide significance, respectively, contrasting with none hypomethylating. Validation and replication (450k array and independent targeted Bisuphite-sequencing) supported the hypermethylation of this functional unit. Tissue-specificity is a significant driver, although the strongest consistent signals, also independent of major cell-type change, occur in tRNA-iMet-CAT-1-4 and tRNA-Ser-AGA-2-6. This study presents a comprehensive evaluation of the genomic DNA methylation state of human tRNA genes and reveals a discreet hypermethylation with advancing age.

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

Background

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

Methods

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

Results

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

Conclusions

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

On the Track of the Missing tRNA Genes: A Source of Non-Canonical Functions?

Cellular tRNAs appear today as a diverse population of informative macromolecules with conserved general elements ensuring essential common functions and different and distinctive features securing specific interactions and activities. Their differential expression and the variety of post-transcriptional modifications they are subject to, lead to the existence of complex repertoires of tRNA populations adjusted to defined cellular states. Despite the tRNA-coding genes redundancy in prokaryote and eukaryote genomes, it is surprising to note the absence of genes coding specific translational-active isoacceptors throughout the phylogeny. Through the analysis of different releases of tRNA databases, this review aims to provide a general summary about those “missing tRNA genes.” This absence refers to both tRNAs that are not encoded in the genome, as well as others that show critical sequence variations that would prevent their activity as canonical translation adaptor molecules. Notably, while a group of genes are universally missing, others are absent in particular kingdoms. Functional information available allows to hypothesize that the exclusion of isodecoding molecules would be linked to: 1) reduce ambiguities of signals that define the specificity of the interactions in which the tRNAs are involved; 2) ensure the adaptation of the translational apparatus to the cellular state; 3) divert particular tRNA variants from ribosomal protein synthesis to other cellular functions. This leads to consider the “missing tRNA genes” as a source of putative non-canonical tRNA functions and to broaden the concept of adapter molecules in ribosomal-dependent protein synthesis.

Modification of Transfer RNA Levels Affects Cyclin Aggregation and the Correct Duplication of Yeast Cells

Codon usage bias (the preferential use of certain synonymous codons (optimal) over others is found at the organism level (intergenomic) within specific genomes (intragenomic) and even in certain genes. Whether it is the result of genetic drift due to GC/AT content and/or natural selection is a topic of intense debate. Preferential codons are mostly found in genes encoding highly-expressed proteins, while lowly-expressed proteins usually contain a high proportion of rare (lowly-represented) codons. While optimal codons are decoded by highly expressed tRNAs, rare codons are usually decoded by lowly-represented tRNAs. Whether rare codons play a role in controlling the expression of lowly- or temporarily-expressed proteins is an open question. In this work we approached this question using two strategies, either by replacing rare glycine codons with optimal counterparts in the gene that encodes the cell cycle protein Cdc13, or by overexpression the tRNA ^Gly that decodes rare codons from the fission yeast, Schizosaccharomyces pombe. While the replacement of synonymous codons severely affected cell growth, increasing tRNA levels affected the aggregation status of Cdc13 and cell division. These lead us to think that rare codons in lowly-expressed cyclin proteins are crucial for cell division, and that the overexpression of tRNA that decodes rare codons affects the expression of proteins containing these rare codons. These codons may be the result of the natural selection of codons in genes that encode lowly-expressed proteins.

Positive epistasis between disease-causing missense mutations and silent polymorphism with effect on mRNA translation velocity

Epistasis refers to the dependence of a mutation on other mutation(s) and the genetic context in general. In the context of human disorders, epistasis complicates the spectrum of disease symptoms and has been proposed as a major contributor to variations in disease outcome. The nonadditive relationship between mutations and the lack of complete understanding of the underlying physiological effects limit our ability to predict phenotypic outcome. Here, we report positive epistasis between intragenic mutations in the cystic fibrosis transmembrane conductance regulator (CFTR)-the gene responsible for cystic fibrosis (CF) pathology. We identified a synonymous single-nucleotide polymorphism (sSNP) that is invariant for the CFTR amino acid sequence but inverts translation speed at the affected codon. This sSNP in cis exhibits positive epistatic effects on some CF disease-causing missense mutations. Individually, both mutations alter CFTR structure and function, yet when combined, they lead to enhanced protein expression and activity. The most robust effect was observed when the sSNP was present in combination with missense mutations that, along with the primary amino acid change, also alter the speed of translation at the affected codon. Functional studies revealed that synergistic alteration in ribosomal velocity is the underlying mechanism; alteration of translation speed likely increases the time window for establishing crucial domain-domain interactions that are otherwise perturbed by each individual mutation.

Synonymous but Not Silent: The Codon Usage Code for Gene Expression and Protein Folding

Codon usage bias, the preference for certain synonymous codons, is found in all genomes. Although synonymous mutations were previously thought to be silent, a large body of evidence has demonstrated that codon usage can play major roles in determining gene expression levels and protein structures. Codon usage influences translation elongation speed and regulates translation efficiency and accuracy. Adaptation of codon usage to tRNA expression determines the proteome landscape. In addition, codon usage biases result in nonuniform ribosome decoding rates on mRNAs, which in turn influence the cotranslational protein folding process that is critical for protein function in diverse biological processes. Conserved genome-wide correlations have also been found between codon usage and protein structures. Furthermore, codon usage is a major determinant of mRNA levels through translation-dependent effects on mRNA decay and translation-independent effects on transcriptional and posttranscriptional processes. Here, we discuss the multifaceted roles and mechanisms of codon usage in different gene regulatory processes.