Celebrating wobble decoding: Half a century and still much is new

[4Fe-4S]-dependent enzymes in non-redox tRNA thiolation

Post-transcriptional modification of nucleosides in transfer RNAs (tRNAs) is an important process for accurate and efficient translation of the genetic information during protein synthesis in all domains of life. In particular, specific enzymes catalyze the biosynthesis of sulfur-containing nucleosides, such as the derivatives of 2-thiouridine (s2U), 4-thiouridine (s4U), 2-thiocytidine (s2C), and 2-methylthioadenosine (ms2A), within tRNAs. Whereas the mechanism that has prevailed for decades involved persulfide chemistry, more and more tRNA thiolation enzymes have now been shown to contain a [4Fe-4S] cluster. This review summarizes the information over the last ten years concerning the biochemical, spectroscopic and structural characterization of [4Fe-4S]-dependent non-redox tRNA thiolation enzymes.

Maximal Genetic Code Symmetry Is a Physicochemical Purine-Pyrimidine Symmetry Language for Transcription and Translation in the Flow of Genetic Information from DNA to Proteins

Until now, research has not taken into consideration the physicochemical purine-pyrimidine symmetries of the genetic code in the transcription and translation processes of proteinogenesis. Our Supersymmetry Genetic Code table, developed in 2022, is common and unique for all RNA and DNA living species. Its basic structure is a purine-pyrimidine symmetry net with double mirror symmetry. Accordingly, the symmetry of the genetic code directly shows its organisation based on the principle of nucleotide Watson-Crick and codon-anticodon pairing. The maximal purine-pyrimidine symmetries of codons show that each codon has a strictly defined and unchangeable position within the genetic code. We discovered that the physicochemical symmetries of the genetic code play a fundamental role in recognising and differentiating codons from mRNA and the anticodon tRNA and aminoacyl-tRNA synthetases in the transcription and translation processes. These symmetries also support the wobble hypothesis with non-Watson-Crick pairing interactions between the translation process from mRNA to tRNA. The Supersymmetry Genetic Code table shows a specific arrangement of the second base of codons, according to which it is possible that an anticodon from tRNA recognises whether a codon from mRNA belongs to an amino acid with two or four codons, which is very important in the purposeful use of the wobble pairing process. Therefore, we show that canonical and wobble pairings essentially do not lead to misreading and errors during translation, and we point out the role of physicochemical purine-pyrimidine symmetries in decreasing disorder according to error minimisation and preserving the integrity of biological processes during proteinogenesis.

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.

RNA modifying enzymes shape tRNA biogenesis and function

Transfer RNAs (tRNAs) are the most highly modified cellular RNAs, both with respect to the proportion of nucleotides that are modified within the tRNA sequence and with respect to the extraordinary diversity in tRNA modification chemistry. However, the functions of many different tRNA modifications are only beginning to emerge. tRNAs have two general clusters of modifications. The first cluster is within the anticodon stem-loop including several modifications essential for protein translation. The second cluster of modifications is within the tRNA elbow, and roles for these modifications are less clear. In general, tRNA elbow modifications are typically not essential for cell growth, but nonetheless several tRNA elbow modifications have been highly conserved throughout all domains of life. In addition to forming modifications, many tRNA modifying enzymes have been demonstrated or hypothesized to also play an important role in folding tRNA acting as tRNA chaperones. In this review, we summarize the known functions of tRNA modifying enzymes throughout the lifecycle of a tRNA molecule, from transcription to degradation. Thereby, we describe how tRNA modification and folding by tRNA modifying enzymes enhance tRNA maturation, tRNA aminoacylation, and tRNA function during protein synthesis, ultimately impacting cellular phenotypes and disease.

m3C32 tRNA modification controls serine codon-biased mRNA translation, cell cycle, and DNA-damage response

The epitranscriptome includes a diversity of RNA modifications that influence gene expression. N3-methylcytidine (m3C) mainly occurs in the anticodon loop (position C32) of certain tRNAs yet its role is poorly understood. Here, using HAC-Seq, we report comprehensive METTL2A/2B-, METTL6-, and METTL2A/2B/6-dependent m3C profiles in human cells. METTL2A/2B modifies tRNA-arginine and tRNA-threonine members, whereas METTL6 modifies the tRNA-serine family. However, decreased m3C32 on tRNA-Ser-GCT isodecoders is only observed with combined METTL2A/2B/6 deletion. Ribo-Seq reveals altered translation of genes related to cell cycle and DNA repair pathways in METTL2A/2B/6-deficient cells, and these mRNAs are enriched in AGU codons that require tRNA-Ser-GCT for translation. These results, supported by reporter assays, help explain the observed altered cell cycle, slowed proliferation, and increased cisplatin sensitivity phenotypes of METTL2A/2B/6-deficient cells. Thus, we define METTL2A/2B/6-dependent methylomes and uncover a particular requirement of m3C32 tRNA modification for serine codon-biased mRNA translation of cell cycle, and DNA repair genes.

Mistranslating tRNA variants have anticodon- and sex-specific impacts on Drosophila melanogaster

Transfer RNAs (tRNAs) are vital in determining the specificity of translation. Mutations in tRNA genes can result in the misincorporation of amino acids into nascent polypeptides in a process known as mistranslation. Since mistranslation has different impacts, depending on the type of amino acid substitution, our goal here was to compare the impact of different mistranslating tRNASer variants on fly development, lifespan, and behaviour. We established two mistranslating fly lines, one with a tRNASer variant that misincorporates serine at valine codons (V→S) and the other that misincorporates serine at threonine codons (T→S). While both mistranslating tRNAs increased development time and developmental lethality, the severity of the impacts differed depending on amino acid substitution and sex. The V→S variant extended embryonic, larval, and pupal development whereas the T→S only extended larval and pupal development. Females, but not males, containing either mistranslating tRNA presented with significantly more anatomical deformities than controls. Mistranslating females also experienced extended lifespan whereas mistranslating male lifespan was unaffected. In addition, mistranslating flies from both sexes showed improved locomotion as they aged, suggesting delayed neurodegeneration. Therefore, although mistranslation causes detrimental effects, we demonstrate that mistranslation also has positive effects on complex traits such as lifespan and locomotion. This has important implications for human health given the prevalence of tRNA variants in humans.

A tRNA modification pattern that facilitates interpretation of the genetic code

Interpretation of the genetic code from triplets of nucleotides to amino acids is fundamental to life. This interpretation is achieved by cellular tRNAs, each reading a triplet codon through its complementary anticodon (positions 34-36) while delivering the amino acid charged to its 3'-end. This amino acid is then incorporated into the growing polypeptide chain during protein synthesis on the ribosome. The quality and versatility of the interpretation is ensured not only by the codon-anticodon pairing, but also by the post-transcriptional modifications at positions 34 and 37 of each tRNA, corresponding to the wobble nucleotide at the first position of the anticodon and the nucleotide on the 3'-side of the anticodon, respectively. How each codon is read by the matching anticodon, and which modifications are required, cannot be readily predicted from the codon-anticodon pairing alone. Here we provide an easily accessible modification pattern that is integrated into the genetic code table. We focus on the Gram-negative bacterium Escherichia coli as a model, which is one of the few organisms whose entire set of tRNA modifications and modification genes is identified and mapped. This work provides an important reference tool that will facilitate research in protein synthesis, which is at the core of the cellular life.

Functional redundancy in tRNA dihydrouridylation

Dihydrouridine (D) is a common modified base found predominantly in transfer RNA (tRNA). Despite its prevalence, the mechanisms underlying dihydrouridine biosynthesis, particularly in prokaryotes, have remained elusive. Here, we conducted a comprehensive investigation into D biosynthesis in Bacillus subtilis through a combination of genetic, biochemical, and epitranscriptomic approaches. Our findings reveal that B. subtilis relies on two FMN-dependent Dus-like flavoprotein homologs, namely DusB1 and DusB2, to introduce all D residues into its tRNAs. Notably, DusB1 exhibits multisite enzyme activity, enabling D formation at positions 17, 20, 20a and 47, while DusB2 specifically catalyzes D biosynthesis at positions 20 and 20a, showcasing a functional redundancy among modification enzymes. Extensive tRNA-wide D-mapping demonstrates that this functional redundancy impacts the majority of tRNAs, with DusB2 displaying a higher dihydrouridylation efficiency compared to DusB1. Interestingly, we found that BsDusB2 can function like a BsDusB1 when overexpressed in vivo and under increasing enzyme concentration in vitro. Furthermore, we establish the importance of the D modification for B. subtilis growth at suboptimal temperatures. Our study expands the understanding of D modifications in prokaryotes, highlighting the significance of functional redundancy in this process and its impact on bacterial growth and adaptation.

Exploring the Impact of mRNA Modifications on Translation Efficiency and Immune Tolerance to Self-Antigens

Therapeutic modified mRNAs are being developed for a broad range of human diseases. However, the impact of potential miscoding of modified mRNAs on self-tolerance remains unknown. Additionally, more studies are needed to explore the effects of nucleoside alkylation on translation. While all six tested modifications are tolerated as substrates by T7 RNA polymerase and inhibited mRNA immunogenicity, the translation efficiency varied significantly depending on the type of modification. In contrast to methylation, ethylation at the N1 position of pseudouridine (Ψ) hindered translation, suggesting that the C5-C1' glycosidic bond alone is not a critical element for high translation. Inhibition of mRNA translation was also observed with 5-methoxyuridine modification. However, this inhibition was partially alleviated through the optimization of mRNA coding sequences. BALB/c mice immunized with syngeneic ψ-modified mRNA encoding for Wilms' tumor antigen-1 (WT1) developed a low but significant level of anti-WT1 IgG antibodies compared to those immunized with either unmodified or N1-methyl ψ-modified mRNA. Overall, the data indicate that adding a simple ethyl group (-CH2CH3) at the N1 position of ψ has a major negative effect on translation despite its reduced immunogenicity. Additionally, mRNA containing Ψ may alter translation fidelity at certain codons, which could lead to a breakdown of immune tolerance to self-antigens. This concern should be taken into account during gene replacement therapies, although it could benefit mRNA-based vaccines by generating a diverse repertoire of antigens.

tRNA modification reprogramming contributes to artemisinin resistance in Plasmodium falciparum

Plasmodium falciparum artemisinin (ART) resistance is driven by mutations in kelch-like protein 13 (PfK13). Quiescence, a key aspect of resistance, may also be regulated by a yet unidentified epigenetic pathway. Transfer RNA modification reprogramming and codon bias translation is a conserved epitranscriptomic translational control mechanism that allows cells to rapidly respond to stress. We report a role for this mechanism in ART-resistant parasites by combining tRNA modification, proteomic and codon usage analyses in ring-stage ART-sensitive and ART-resistant parasites in response to drug. Post-drug, ART-resistant parasites differentially hypomodify mcm5s2U on tRNA and possess a subset of proteins, including PfK13, that are regulated by Lys codon-biased translation. Conditional knockdown of the terminal s2U thiouridylase, PfMnmA, in an ART-sensitive parasite background led to increased ART survival, suggesting that hypomodification can alter the parasite ART response. This study describes an epitranscriptomic pathway via tRNA s2U reprogramming that ART-resistant parasites may employ to survive ART-induced stress.

Robust genetic codes enhance protein evolvability

The standard genetic code defines the rules of translation for nearly every life form on Earth. It also determines the amino acid changes accessible via single-nucleotide mutations, thus influencing protein evolvability-the ability of mutation to bring forth adaptive variation in protein function. One of the most striking features of the standard genetic code is its robustness to mutation, yet it remains an open question whether such robustness facilitates or frustrates protein evolvability. To answer this question, we use data from massively parallel sequence-to-function assays to construct and analyze 6 empirical adaptive landscapes under hundreds of thousands of rewired genetic codes, including those of codon compression schemes relevant to protein engineering and synthetic biology. We find that robust genetic codes tend to enhance protein evolvability by rendering smooth adaptive landscapes with few peaks, which are readily accessible from throughout sequence space. However, the standard genetic code is rarely exceptional in this regard, because many alternative codes render smoother landscapes than the standard code. By constructing low-dimensional visualizations of these landscapes, which each comprise more than 16 million mRNA sequences, we show that such alternative codes radically alter the topological features of the network of high-fitness genotypes. Whereas the genetic codes that optimize evolvability depend to some extent on the detailed relationship between amino acid sequence and protein function, we also uncover general design principles for engineering nonstandard genetic codes for enhanced and diminished evolvability, which may facilitate directed protein evolution experiments and the bio-containment of synthetic organisms, respectively.

Epigenetic modifications of DNA and RNA in Alzheimer's disease

Alzheimer's disease (AD) is a complex neurodegenerative disorder and the most common form of dementia. There are two main types of AD: familial and sporadic. Familial AD is linked to mutations in amyloid precursor protein (APP), presenilin-1 (PSEN1), and presenilin-2 (PSEN2). On the other hand, sporadic AD is the more common form of the disease and has genetic, epigenetic, and environmental components that influence disease onset and progression. Investigating the epigenetic mechanisms associated with AD is essential for increasing understanding of pathology and identifying biomarkers for diagnosis and treatment. Chemical covalent modifications on DNA and RNA can epigenetically regulate gene expression at transcriptional and post-transcriptional levels and play protective or pathological roles in AD and other neurodegenerative diseases.

Codon decoding by orthogonal tRNAs interrogates the in vivo preferences of unmodified adenosine in the wobble position

The in vivo codon decoding preferences of tRNAs with an authentic adenosine residue at position 34 of the anticodon, the wobble position, are largely unexplored because very few unmodified A34 tRNA genes exist across the three domains of life. The expanded wobble rules suggest that unmodified adenosine pairs most strongly with uracil, modestly with cytosine, and weakly with guanosine and adenosine. Inosine, a modified adenosine, on the other hand, pairs strongly with both uracil and cytosine and to a lesser extent adenosine. Orthogonal pair directed sense codon reassignment experiments offer a tool with which to interrogate the translational activity of A34 tRNAs because the introduced tRNA can be engineered with any anticodon. Our fluorescence-based screen utilizes the absolute requirement of tyrosine at position 66 of superfolder GFP for autocatalytic fluorophore formation. The introduced orthogonal tRNA competes with the endogenous translation machinery to incorporate tyrosine in response to a codon typically assigned another meaning in the genetic code. We evaluated the codon reassignment efficiencies of 15 of the 16 possible orthogonal tRNAs with A34 anticodons. We examined the Sanger sequencing chromatograms for cDNAs from each of the reverse transcribed tRNAs for evidence of inosine modification. Despite several A34 tRNAs decoding closely-related C-ending codons, partial inosine modification was detected for only three species. These experiments employ a single tRNA body with a single attached amino acid to interrogate the behavior of different anticodons in the background of in vivo E. coli translation and greatly expand the set of experimental measurements of the in vivo function of A34 tRNAs in translation. For the most part, unmodified A34 tRNAs largely pair with only U3 codons as the original wobble rules suggest. In instances with GC pairs in the first two codon positions, unmodified A34 tRNAs decode the C- and G-ending codons as well as the expected U-ending codon. These observations support the "two-out-of-three" and "strong and weak" codon hypotheses.

Alternate routes to mnm5s2U synthesis in Gram-positive bacteria

The wobble bases of tRNAs that decode split codons are often heavily modified. In bacteria, tRNAGlu, Gln, Asp contains a variety of xnm5s2U derivatives. The synthesis pathway for these modifications is complex and fully elucidated only in a handful of organisms, including the Gram-negative Escherichia coli K12 model. Despite the ubiquitous presence of mnm5s2U modification, genomic analysis shows the absence of mnmC orthologous genes, suggesting the occurrence of alternate biosynthetic schemes for the conversion of cmnm5s2U to mnm5s2U. Using a combination of comparative genomics and genetic studies, a member of the YtqA subgroup of the radical Sam superfamily was found to be involved in the synthesis of mnm5s2U in both Bacillus subtilis and Streptococcus mutans. This protein, renamed MnmL, is encoded in an operon with the recently discovered MnmM methylase involved in the methylation of the pathway intermediate nm5s2U into mnm5s2U in B. subtilis. Analysis of tRNA modifications of both S. mutans and Streptococcus pneumoniae shows that growth conditions and genetic backgrounds influence the ratios of pathway intermediates owing to regulatory loops that are not yet understood. The MnmLM pathway is widespread along the bacterial tree, with some phyla, such as Bacilli, relying exclusively on these two enzymes. Although mechanistic details of these newly discovered components are not fully resolved, the occurrence of fusion proteins, alternate arrangements of biosynthetic components, and loss of biosynthetic branches provide examples of biosynthetic diversity to retain a conserved tRNA modification in Nature.IMPORTANCEThe xnm5s2U modifications found in several tRNAs at the wobble base position are widespread in bacteria where they have an important role in decoding efficiency and accuracy. This work identifies a novel enzyme (MnmL) that is a member of a subgroup of the very versatile radical SAM superfamily and is involved in the synthesis of mnm5s2U in several Gram-positive bacteria, including human pathogens. This is another novel example of a non-orthologous displacement in the field of tRNA modification synthesis, showing how different solutions evolve to retain U34 tRNA modifications.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Contribution of tRNA sequence and modifications to the decoding preferences of E. coli and M. mycoides tRNAGlyUCC for synonymous glycine codons

tRNA superwobbling, used by certain bacteria and organelles, is an intriguing decoding concept in which a single tRNA isoacceptor is used to decode all synonymous codons of a four-fold degenerate codon box. While Escherichia coli relies on three tRNAGly isoacceptors to decode the four glycine codons (GGN), Mycoplasma mycoides requires only a single tRNAGly. Both organisms express tRNAGly with the anticodon UCC, which are remarkably similar in sequence but different in their decoding ability. By systematically introducing mutations and altering the number and type of tRNA modifications using chemically synthesized tRNAs, we elucidated the contribution of individual nucleotides and chemical groups to decoding by the E. coli and M. mycoides tRNAGly. The tRNA sequence was identified as the key factor for superwobbling, revealing the T-arm sequence as a novel pivotal element. In addition, the presence of tRNA modifications, although not essential for providing superwobbling, was shown to delicately fine-tune and balance the decoding of synonymous codons. This emphasizes that the tRNA sequence and its modifications together form an intricate system of high complexity that is indispensable for accurate and efficient decoding.

The tRNA thiolation-mediated translational control is essential for plant immunity

Plants have evolved sophisticated mechanisms to regulate gene expression to activate immune responses against pathogen infections. However, how the translation system contributes to plant immunity is largely unknown. The evolutionarily conserved thiolation modification of transfer RNA (tRNA) ensures efficient decoding during translation. Here, we show that tRNA thiolation is required for plant immunity in Arabidopsis. We identify a cgb mutant that is hyper-susceptible to the pathogen Pseudomonas syringae. CGB encodes ROL5, a homolog of yeast NCS6 required for tRNA thiolation. ROL5 physically interacts with CTU2, a homolog of yeast NCS2. Mutations in either ROL5 or CTU2 result in loss of tRNA thiolation. Further analyses reveal that both transcriptome and proteome reprogramming during immune responses are compromised in cgb. Notably, the translation of salicylic acid receptor NPR1 is reduced in cgb, resulting in compromised salicylic acid signaling. Our study not only reveals a regulatory mechanism for plant immunity but also uncovers an additional biological function of tRNA thiolation.

Quantitative analysis of tRNA abundance and modifications by nanopore RNA sequencing

Transfer RNAs (tRNAs) play a central role in protein translation. Studying them has been difficult in part because a simple method to simultaneously quantify their abundance and chemical modifications is lacking. Here we introduce Nano-tRNAseq, a nanopore-based approach to sequence native tRNA populations that provides quantitative estimates of both tRNA abundances and modification dynamics in a single experiment. We show that default nanopore sequencing settings discard the vast majority of tRNA reads, leading to poor sequencing yields and biased representations of tRNA abundances based on their transcript length. Re-processing of raw nanopore current intensity signals leads to a 12-fold increase in the number of recovered tRNA reads and enables recapitulation of accurate tRNA abundances. We then apply Nano-tRNAseq to Saccharomyces cerevisiae tRNA populations, revealing crosstalks and interdependencies between different tRNA modification types within the same molecule and changes in tRNA populations in response to oxidative stress.

Alternate routes to mnm 5 s 2 U synthesis in Gram-positive bacteria

The wobble bases of tRNAs that decode split codons are often heavily modified. In Bacteria tRNA Glu, Gln, Asp contain a variety of xnm 5 s 2 U derivatives. The synthesis pathway for these modifications is complex and fully elucidated only in a handful of organisms, including the Gram-negative Escherichia coli K12 model. Despite the ubiquitous presence of mnm 5 s 2 U modification, genomic analysis shows the absence of mnmC orthologous genes, suggesting the occurrence of alternate biosynthetic schemes for the installation of this modification. Using a combination of comparative genomics and genetic studies, a member of the YtqA subgroup of the Radical Sam superfamily was found to be involved in the synthesis of mnm 5 s 2 U in both Bacillus subtilis and Streptococcus mutans . This protein, renamed MnmL, is encoded in an operon with the recently discovered MnmM methylase involved in the methylation of the pathway intermediate nm 5 s 2 U into mnm 5 s 2 U in B. subtilis . Analysis of tRNA modifications of both S. mutans and Streptococcus pneumoniae shows that growth conditions and genetic backgrounds influence the ratios of pathways intermediates in regulatory loops that are not yet understood. The MnmLM pathway is widespread along the bacterial tree, with some phyla, such as Bacilli, relying exclusively on these two enzymes. The occurrence of fusion proteins, alternate arrangements of biosynthetic components, and loss of biosynthetic branches provide examples of biosynthetic diversity to retain a conserved tRNA modification in nature.

Importance

The xnm 5 s 2 U modifications found in several tRNAs at the wobble base position are widespread in Bacteria where they have an important role in decoding efficiency and accuracy. This work identifies a novel enzyme (MnmL) that is a member of a subgroup of the very versatile Radical SAM superfamily and is involved in the synthesis of mnm 5 s 2 U in several Gram-positive bacteria, including human pathogens. This is another novel example of a non-orthologous displacement in the field of tRNA modification synthesis, showing how different solutions evolve to retain U34 tRNA modifications.

Phosphate Limitation Responses in Marine Green Algae Are Linked to Reprogramming of the tRNA Epitranscriptome and Codon Usage Bias

Marine algae are central to global carbon fixation, and their productivity is dictated largely by resource availability. Reduced nutrient availability is predicted for vast oceanic regions as an outcome of climate change; however, there is much to learn regarding response mechanisms of the tiny picoplankton that thrive in these environments, especially eukaryotic phytoplankton. Here, we investigate responses of the picoeukaryote Micromonas commoda, a green alga found throughout subtropical and tropical oceans. Under shifting phosphate availability scenarios, transcriptomic analyses revealed altered expression of transfer RNA modification enzymes and biased codon usage of transcripts more abundant during phosphate-limiting versus phosphate-replete conditions, consistent with the role of transfer RNA modifications in regulating codon recognition. To associate the observed shift in the expression of the transfer RNA modification enzyme complement with the transfer RNAs encoded by M. commoda, we also determined the transfer RNA repertoire of this alga revealing potential targets of the modification enzymes. Codon usage bias was particularly pronounced in transcripts encoding proteins with direct roles in managing phosphate limitation and photosystem-associated proteins that have ill-characterized putative functions in "light stress." The observed codon usage bias corresponds to a proposed stress response mechanism in which the interplay between stress-induced changes in transfer RNA modifications and skewed codon usage in certain essential response genes drives preferential translation of the encoded proteins. Collectively, we expose a potential underlying mechanism for achieving growth under enhanced nutrient limitation that extends beyond the catalog of up- or downregulated protein-encoding genes to the cell biological controls that underpin acclimation to changing environmental conditions.

Queuosine-tRNA promotes sex-dependent learning and memory formation by maintaining codon-biased translation elongation speed

Queuosine (Q) is a modified nucleoside at the wobble position of specific tRNAs. In mammals, queuosinylation is facilitated by queuine uptake from the gut microbiota and is introduced into tRNA by the QTRT1-QTRT2 enzyme complex. By establishing a Qtrt1 knockout mouse model, we discovered that the loss of Q-tRNA leads to learning and memory deficits. Ribo-Seq analysis in the hippocampus of Qtrt1-deficient mice revealed not only stalling of ribosomes on Q-decoded codons, but also a global imbalance in translation elongation speed between codons that engage in weak and strong interactions with their cognate anticodons. While Q-dependent molecular and behavioral phenotypes were identified in both sexes, female mice were affected more severely than males. Proteomics analysis confirmed deregulation of synaptogenesis and neuronal morphology. Together, our findings provide a link between tRNA modification and brain functions and reveal an unexpected role of protein synthesis in sex-dependent cognitive performance.

yaaJ, the tRNA-Specific Adenosine Deaminase, Is Dispensable in Bacillus subtilis

Post-transcriptional modifications of tRNA are crucial for their core function. The inosine (I; 6-deaminated adenosine) at the first position in the anticodon of tRNAArg(ICG) modulates the decoding capability and is generally considered essential for reading CGU, CGC, and CGA codons in eubacteria. We report here that the Bacillus subtilis yaaJ gene encodes tRNA-specific adenosine deaminase and is non-essential for viability. A β-galactosidase reporter assay revealed that the translational activity of CGN codons was not impaired in the yaaJ-deletion mutant. Furthermore, tRNAArg(CCG) responsible for decoding the CGG codon was dispensable, even in the presence or absence of yaaJ. These results strongly suggest that tRNAArg with either the anticodon ICG or ACG has an intrinsic ability to recognize all four CGN codons, providing a fundamental concept of non-canonical wobbling mediated by adenosine and inosine nucleotides in the anticodon. This is the first example of the four-way wobbling by inosine nucleotide in bacterial cells. On the other hand, the absence of inosine modification induced +1 frameshifting, especially at the CGA codon. Additionally, the yaaJ deletion affected growth and competency. Therefore, the inosine modification is beneficial for translational fidelity and proper growth-phase control, and that is why yaaJ has been actually conserved in B. subtilis.

Broad anti-pathogen potential of DEAD box RNA helicase eIF4A-targeting rocaglates

Inhibition of eukaryotic initiation factor 4A has been proposed as a strategy to fight pathogens. Rocaglates exhibit the highest specificities among eIF4A inhibitors, but their anti-pathogenic potential has not been comprehensively assessed across eukaryotes. In silico analysis of the substitution patterns of six eIF4A1 aa residues critical to rocaglate binding, uncovered 35 variants. Molecular docking of eIF4A:RNA:rocaglate complexes, and in vitro thermal shift assays with select recombinantly expressed eIF4A variants, revealed that sensitivity correlated with low inferred binding energies and high melting temperature shifts. In vitro testing with silvestrol validated predicted resistance in Caenorhabditis elegans and Leishmania amazonensis and predicted sensitivity in Aedes sp., Schistosoma mansoni, Trypanosoma brucei, Plasmodium falciparum, and Toxoplasma gondii. Our analysis further revealed the possibility of targeting important insect, plant, animal, and human pathogens with rocaglates. Finally, our findings might help design novel synthetic rocaglate derivatives or alternative eIF4A inhibitors to fight pathogens.

The thiolation of uridine 34 in tRNA, which controls protein translation, depends on a [4Fe-4S] cluster in the archaeum Methanococcus maripaludis

Thiolation of uridine 34 in the anticodon loop of several tRNAs is conserved in the three domains of life and guarantees fidelity of protein translation. U34-tRNA thiolation is catalyzed by a complex of two proteins in the eukaryotic cytosol (named Ctu1/Ctu2 in humans), but by a single NcsA enzyme in archaea. We report here spectroscopic and biochemical experiments showing that NcsA from Methanococcus maripaludis (MmNcsA) is a dimer that binds a [4Fe-4S] cluster, which is required for catalysis. Moreover, the crystal structure of MmNcsA at 2.8 Å resolution shows that the [4Fe-4S] cluster is coordinated by three conserved cysteines only, in each monomer. Extra electron density on the fourth nonprotein-bonded iron most likely locates the binding site for a hydrogenosulfide ligand, in agreement with the [4Fe-4S] cluster being used to bind and activate the sulfur atom of the sulfur donor. Comparison of the crystal structure of MmNcsA with the AlphaFold model of the human Ctu1/Ctu2 complex shows a very close superposition of the catalytic site residues, including the cysteines that coordinate the [4Fe-4S] cluster in MmNcsA. We thus propose that the same mechanism for U34-tRNA thiolation, mediated by a [4Fe-4S]-dependent enzyme, operates in archaea and eukaryotes.

Identification of a novel 5-aminomethyl-2-thiouridine methyltransferase in tRNA modification

The uridine at the 34th position of tRNA, which is able to base pair with the 3'-end codon on mRNA, is usually modified to influence many aspects of decoding properties during translation. Derivatives of 5-methyluridine (xm5U), which include methylaminomethyl (mnm-) or carboxymethylaminomethyl (cmnm-) groups at C5 of uracil base, are widely conserved at the 34th position of many prokaryotic tRNAs. In Gram-negative bacteria such as Escherichia coli, a bifunctional MnmC is involved in the last two reactions of the biosynthesis of mnm5(s2)U, in which the enzyme first converts cmnm5(s2)U to 5-aminomethyl-(2-thio)uridine (nm5(s2)U) and subsequently installs the methyl group to complete the formation of mnm5(s2)U. Although mnm5s2U has been identified in tRNAs of Gram-positive bacteria and plants as well, their genomes do not contain an mnmC ortholog and the gene(s) responsible for this modification is unknown. We discovered that MnmM, previously known as YtqB, is the methyltransferase that converts nm5s2U to mnm5s2U in Bacillus subtilis through comparative genomics, gene complementation experiments, and in vitro assays. Furthermore, we determined X-ray crystal structures of MnmM complexed with anticodon stem loop of tRNAGln. The structures provide the molecular basis underlying the importance of U33-nm5s2U34-U35 as the key determinant for the specificity of MnmM.

Modopathies Caused by Mutations in Genes Encoding for Mitochondrial RNA Modifying Enzymes: Molecular Mechanisms and Yeast Disease Models

In eukaryotes, mitochondrial RNAs (mt-tRNAs and mt-rRNAs) are subject to specific nucleotide modifications, which are critical for distinct functions linked to the synthesis of mitochondrial proteins encoded by mitochondrial genes, and thus for oxidative phosphorylation. In recent years, mutations in genes encoding for mt-RNAs modifying enzymes have been identified as being causative of primary mitochondrial diseases, which have been called modopathies. These latter pathologies can be caused by mutations in genes involved in the modification either of tRNAs or of rRNAs, resulting in the absence of/decrease in a specific nucleotide modification and thus on the impairment of the efficiency or the accuracy of the mitochondrial protein synthesis. Most of these mutations are sporadic or private, thus it is fundamental that their pathogenicity is confirmed through the use of a model system. This review will focus on the activity of genes that, when mutated, are associated with modopathies, on the molecular mechanisms through which the enzymes introduce the nucleotide modifications, on the pathological phenotypes associated with mutations in these genes and on the contribution of the yeast Saccharomyces cerevisiae to confirming the pathogenicity of novel mutations and, in some cases, for defining the molecular defects.

Application of mutational profiling: New functional analyses reveal the tRNA recognition mechanism of tRNA m1A22 methyltransferase

Transfer RNAs undergo diverse posttranscriptional modifications to regulate a myriad of cellular events including translation, stress response, and viral replication. These posttranscriptional modifications are synthesized by site-specific modification enzymes. Recent RNA-seq techniques have revealed multiple features of tRNA such as tRNA abundance, tRNA modification, and tRNA structure. Here, we adapt a tRNA-sequencing technique and design a new functional analysis where we perform mutational profiling of tRNA modifications to gain mechanistic insights into how tRNA modification enzymes recognize substrate tRNA. Profiling of Geobacillus stearothermophilus tRNAs and protein orthology analysis predict the existence of natural modifications in 44 tRNA molecular species of G. stearothermophilus. We selected the 1-methyladenosine modification at position 22 (m1A22) and tRNA (m1A22) methyltransferase (TrmK) for further analysis. Relative quantification of m1A22 levels in 59 tRNA transcripts by mutational profiling reveals that TrmK selectively methylates a subset of tRNAs. Using 240 variants of tRNALeu transcripts, we demonstrate the conserved nucleosides including U8, A14, G15, G18, G19, U55, Purine57, and A58 are important for the methyl transfer reaction of TrmK. Additional biochemical experiments reveal that TrmK strictly recognizes U8, A14, G18, and U55 in tRNA. Furthermore, these findings from tRNALeu variants were crossvalidated using variants of three different tRNA species. Finally, a model of the TrmK-tRNA complex structure was constructed based on our findings and previous biochemical and structural studies by others. Collectively, our study expands functional analyses of tRNA modification enzyme in a high-throughput manner where our assay rapidly identifies substrates from a large pool of tRNAs.

Anticodon sequence determines the impact of mistranslating tRNAAla variants

Transfer RNAs (tRNAs) maintain translation fidelity through accurate charging by their cognate aminoacyl-tRNA synthetase and codon:anticodon base pairing with the mRNA at the ribosome. Mistranslation occurs when an amino acid not specified by the genetic message is incorporated into proteins and has applications in biotechnology, therapeutics and is relevant to disease. Since the alanyl-tRNA synthetase uniquely recognizes a G3:U70 base pair in tRNAAla and the anticodon plays no role in charging, tRNAAla variants with anticodon mutations have the potential to mis-incorporate alanine. Here, we characterize the impact of the 60 non-alanine tRNAAla anticodon variants on the growth of Saccharomyces cerevisiae. Overall, 36 tRNAAla anticodon variants decreased growth in single- or multi-copy. Mass spectrometry analysis of the cellular proteome revealed that 52 of 57 anticodon variants, not decoding alanine or stop codons, induced mistranslation when on single-copy plasmids. Variants with G/C-rich anticodons resulted in larger growth deficits than A/U-rich variants. In most instances, synonymous anticodon variants impact growth differently, with anticodons containing U at base 34 being the least impactful. For anticodons generating the same amino acid substitution, reduced growth generally correlated with the abundance of detected mistranslation events. Differences in decoding specificity, even between synonymous anticodons, resulted in each tRNAAla variant mistranslating unique sets of peptides and proteins. We suggest that these differences in decoding specificity are also important in determining the impact of tRNAAla anticodon variants.

Short tRNA anticodon stem and mutant eRF1 allow stop codon reassignment

Cognate tRNAs deliver specific amino acids to translating ribosomes according to the standard genetic code, and three codons with no cognate tRNAs serve as stop codons. Some protists have reassigned all stop codons as sense codons, neglecting this fundamental principle^1-4. Here we analyse the in-frame stop codons in 7,259 predicted protein-coding genes of a previously undescribed trypanosomatid, Blastocrithidia nonstop. We reveal that in this species in-frame stop codons are underrepresented in genes expressed at high levels and that UAA serves as the only termination codon. Whereas new tRNAs^Glu fully cognate to UAG and UAA evolved to reassign these stop codons, the UGA reassignment followed a different path through shortening the anticodon stem of tRNA^Trp_CCA from five to four base pairs (bp). The canonical 5-bp tRNA^Trp recognizes UGG as dictated by the genetic code, whereas its shortened 4-bp variant incorporates tryptophan also into in-frame UGA. Mimicking this evolutionary twist by engineering both variants from B. nonstop, Trypanosoma brucei and Saccharomyces cerevisiae and expressing them in the last two species, we recorded a significantly higher readthrough for all 4-bp variants. Furthermore, a gene encoding B. nonstop release factor 1 acquired a mutation that specifically restricts UGA recognition, robustly potentiating the UGA reassignment. Virtually the same strategy has been adopted by the ciliate Condylostoma magnum. Hence, we describe a previously unknown, universal mechanism that has been exploited in unrelated eukaryotes with reassigned stop codons.

On distinguishing between canonical tRNA genes and tRNA gene fragments in prokaryotes

Automated genome annotation is essential for extracting biological information from sequence data. The identification and annotation of tRNA genes is frequently performed by the software package tRNAscan-SE, the output of which is listed for selected genomes in the Genomic tRNA database (GtRNAdb). Here, we highlight a pervasive error in prokaryotic tRNA gene sets on GtRNAdb: the mis-categorization of partial, non-canonical tRNA genes as standard, canonical tRNA genes. Firstly, we demonstrate the issue using the tRNA gene sets of 20 organisms from the archaeal taxon Thermococcaceae. According to GtRNAdb, these organisms collectively deviate from the expected set of tRNA genes in 15 instances, including the listing of eleven putative canonical tRNA genes. However, after detailed manual annotation, only one of these eleven remains; the others are either partial, non-canonical tRNA genes resulting from the integration of genetic elements or CRISPR-Cas activity (seven instances), or attributable to ambiguities in input sequences (three instances). Secondly, we show that similar examples of the mis-categorization of predicted tRNA sequences occur throughout the prokaryotic sections of GtRNAdb. While both canonical and non-canonical prokaryotic tRNA gene sequences identified by tRNAscan-SE are biologically interesting, the challenge of reliably distinguishing between them remains. We recommend employing a combination of (i) screening input sequences for the genetic elements typically associated with non-canonical tRNA genes, and ambiguities, (ii) activating the tRNAscan-SE automated pseudogene detection function, and (iii) scrutinizing predicted tRNA genes with low isotype scores. These measures greatly reduce manual annotation efforts, and lead to improved prokaryotic tRNA gene set predictions.

Two isoleucyl tRNAs that decode synonymous codons divergently regulate breast cancer metastatic growth by controlling translation of proliferation-regulating genes

The human genome contains 61 codons encoding 20 amino acids. Synonymous codons representing a given amino acid are decoded by a set of transfer RNAs (tRNAs) called isoacceptors. We report the surprising observation that two isoacceptor tRNAs that decode synonymous codons become modulated in opposing directions during breast cancer progression. Specifically, tRNAIleUAU became upregulated, whereas tRNAIleGAU became repressed as breast cancer cells attained enhanced metastatic capacity. Functionally, tRNAIleUAU promoted and tRNAIleGAU suppressed metastatic colonization in mouse xenograft models. These tRNAs mediated opposing effects on codon-dependent translation of growth-promoting genes, consistent with genomic enrichment or depletion of their cognate codons in mitotic genes. Our findings uncover a specific isoacceptor tRNA pair that act in opposition, divergently impacting growth-regulating genes and a disease phenotype. Degeneracy of the genetic code can thus be biologically exploited by human cancer cells via tRNA isoacceptor shifts that causally facilitate the transition toward a growth-promoting state.

Structural basis for sequence-independent substrate selection by eukaryotic wobble base tRNA deaminase ADAT2/3

The essential deamination of adenosine A34 to inosine at the wobble base is the individual tRNA modification with the greatest effects on mRNA decoding, empowering a single tRNA to translate three different codons. To date, many aspects of how eukaryotic deaminases specifically select their multiple substrates remain unclear. Here, using cryo-EM, we present the structure of a eukaryotic ADAT2/3 deaminase bound to a full-length tRNA, revealing that the enzyme distorts the anticodon loop, but in contrast to the bacterial enzymes, selects its substrate via sequence-independent contacts of eukaryote-acquired flexible or intrinsically unfolded motifs distal from the conserved catalytic core. A gating mechanism for substrate entry to the active site is identified. Our multi-step tRNA recognition model yields insights into how RNA editing by A34 deamination evolved, shaped the genetic code, and directly impacts the eukaryotic proteome.

tRNA methylation resolves codon usage bias at the limit of cell viability

Codon usage of each genome is closely correlated with the abundance of tRNA isoacceptors. How codon usage bias is resolved by tRNA post-transcriptional modifications is largely unknown. Here we demonstrate that the N¹-methylation of guanosine at position 37 (m¹G37) on the 3′-side of the anticodon, while not directly responsible for reading of codons, is a neutralizer that resolves differential decoding of proline codons. A genome-wide suppressor screen of a non-viable Escherichia coli strain, lacking m¹G37, identifies proS suppressor mutations, indicating a coupling of methylation with tRNA prolyl-aminoacylation that sets the limit of cell viability. Using these suppressors, where prolyl-aminoacylation is decoupled from tRNA methylation, we show that m¹G37 neutralizes differential translation of proline codons by the major isoacceptor. Lack of m¹G37 inactivates this neutralization and exposes the need for a minor isoacceptor for cell viability. This work has medical implications for bacterial species that exclusively use the major isoacceptor for survival.

Masuda et al. show that loss of m¹G37 from the 3′ side of the tRNA anticodon renders a modified wobble nucleotide of the anticodon insufficient to decode a set of rare codons, providing a functional underpinning for the “modification circuit” between position 37 and the wobble position of the tRNA anticodon.

The absence of the queuosine tRNA modification leads to pleiotropic phenotypes revealing perturbations of metal and oxidative stress homeostasis in Escherichia coli K12

Queuosine (Q) is a conserved hypermodification of the wobble base of tRNA containing GUN anticodons but the physiological consequences of Q deficiency are poorly understood in bacteria. This work combines transcriptomic, proteomic and physiological studies to characterize a Q-deficient Escherichia coli K12 MG1655 mutant. The absence of Q led to an increased resistance to nickel and cobalt, and to an increased sensitivity to cadmium, compared to the wild-type (WT) strain. Transcriptomic analysis of the WT and Q-deficient strains, grown in the presence and absence of nickel, revealed that the nickel transporter genes (nikABCDE) are downregulated in the Q- mutant, even when nickel is not added. This mutant is therefore primed to resist to high nickel levels. Downstream analysis of the transcriptomic data suggested that the absence of Q triggers an atypical oxidative stress response, confirmed by the detection of slightly elevated reactive oxygen species (ROS) levels in the mutant, increased sensitivity to hydrogen peroxide and paraquat, and a subtle growth phenotype in a strain prone to accumulation of ROS.

Structural Significance of Conformational Preferences and Ribose-Ring-Puckering of Hyper Modified Nucleotide 5'-Monophosphate 2-Methylthio Cyclic N6-Threonylcarbamoyladenosine (p-ms2ct6A) Present at 37th Position in Anticodon Loop of tRNALys

Structural significance of conformational preferences and ribose ring puckering of newly discovered hyper modified nucleotide, 5'-monophosphate 2-methylthio cyclic N⁶-threonylcarbamoyladenosine (p-ms²ct⁶A) have been investigated using quantum chemical semi-empirical RM1 and molecular dynamics simulation techniques. Automated geometry optimization of most stable structure of p-ms²ct⁶A has also been carried out with the help of abinitio (HF SCF, DFT) as well as semi empirical quantum chemical (RM1, AM1, PM3, and PM6) methods. Most stable structure of p-ms²ct⁶A is stabilized by intramolecular interactions between N(3)…HC(2'), N(1)…HC(16), O(13)…HC(15), and O(13)…HO(14). The torsion angles alpha (α) and beta (β) show the significant characteristic patterns with the involvement of intramolecular hydrogen bonding to provide stability to the p-ms²ct⁶A. Further, molecular dynamics simulations of p-ms²ct⁶A revealed the role of ribose sugar ring puckering i.e. C2'-endo and C3'-endo on the structural dynamics of ms²ct⁶A side chain. The modified nucleotide p-ms²ct⁶A periodically prefers both the C2'-endo and C3'-endo sugar with 'anti' and 'syn' conformations. This property of p-ms²ct⁶A could be useful to recognize the starting ANN codons. All atom explicit MD simulation of anticodon loop (ACL) of tRNA^Lys of Bacillus subtilis containing ms²ct⁶A at 37th position showed the U-turn feature, base stacking ability with other adjacent bases and hydrogen bonding interactions similar to the isolated base p-ms²ct⁶A. The ribose sugar puckering contributes to the orientation of the side chain conformation of p-ms²ct⁶A. Thus, the present study could be helpful to understand the structure-function relationship of the hypermodified nucleoside, ms²ct⁶A in recognition of the proper codons AAA/AAG during protein biosynthesis.

2-Selenouridine, a Modified Nucleoside of Bacterial tRNAs, Its Reactivity in the Presence of Oxidizing and Reducing Reagents

The 5-substituted 2-selenouridines are natural components of the bacterial tRNA epitranscriptome. Because selenium-containing biomolecules are redox-active entities, the oxidation susceptibility of 2-selenouridine (Se2U) was studied in the presence of hydrogen peroxide under various conditions and compared with previously reported data for 2-thiouridine (S2U). It was found that Se2U is more susceptible to oxidation and converted in the first step to the corresponding diselenide (Se2U)₂, an unstable intermediate that decomposes to uridine and selenium. The reversibility of the oxidized state of Se2U was demonstrated by the efficient reduction of (Se2U)₂ to Se2U in the presence of common reducing agents. Thus, the 2-selenouridine component of tRNA may have antioxidant potential in cells because of its ability to react with both cellular ROS components and reducing agents. Interestingly, in the course of the reactions studied, we found that (Se2U)₂ reacts with Se2U to form new ‘oligomeric nucleosides′ as linear and cyclic byproducts.

Extracellular N6 -isopentenyladenosine (i6A) addition induces cotranscriptional i6A incorporation into ribosomal RNAs

N⁶ -isopentenyladenosine (i⁶A), a modified adenosine monomer, is known to induce cell death upon its addition to the culture medium. However, the molecular fate of extracellularly added i⁶A has yet to be identified. Here we show that i⁶A addition to cell culture medium results in i⁶A incorporation into cellular RNA in several cell lines, including the 5-fluorouracil (5-FU)-resistant human oral squamous cell carcinoma cell line FR2-SAS and its parental 5-FU-sensitive cell line SAS. i⁶A was predominantly incorporated into 18S and 28S rRNAs, and i⁶A incorporation into total RNA was mostly suppressed by treating these cell lines with an RNA polymerase I (Pol I) inhibitor. i⁶A was incorporated into RNA even upon inactivation of TRIT1, the only cellular i⁶A-modifying enzyme. These results indicate that upon cellular uptake of i⁶A, it is anabolized to be used for Pol I transcription. Interestingly, at lower i⁶A concentrations, the cytotoxic effect of i⁶A was substantially more pronounced in FR2-SAS cells than in SAS cells. Moreover, in FR2-SAS cells, i⁶A treatment decreased the rate of cellular protein synthesis and increased intracellular protein aggregation, and these effects were more pronounced than in SAS cells. Our work provides insights into the molecular fate of extracellularly applied i⁶A in the context of intracellular nucleic acid anabolism and suggests investigation of i⁶A as a candidate for a chemotherapy agent against 5-FU-resistant cancer cells.

Characterization and phylogenetic analysis of the complete mitochondrial genome sequence of Diospyros oleifera, the first representative from the family Ebenaceae

Plant mitochondrial genomes are a valuable source of genetic information for a better understanding of phylogenetic relationships. However, no mitochondrial genome of any species in Ebenaceae has been reported. In this study, we reported the first mitochondrial genome of an Ebenaceae model plant Diospyros oleifera. The mitogenome was 493,958 bp in length, contained 39 protein-coding genes, 27 transfer RNA genes, and 3 ribosomal RNA genes. The rps2 and rps11 genes were missing in the D. oleifera mt genome, while the rps10 gene was identified. The length of the repetitive sequence in the D. oleifera mt genome was 31 kb, accounting for 6.33%. A clear bias in RNA-editing sites were found in the D. oleifera mt genome. We also detected 28 chloroplast-derived fragments significantly associated with D. oleifera mt genes, indicating intracellular tRNA genes transferred frequently from chloroplasts to mitochondria in D. oleifera. Phylogenetic analysis based on the mt genomes of D. oleifera and 27 other taxa reflected the exact evolutionary and taxonomic status of D. oleifera. Ka/Ks analysis revealed that 95.16% of the protein-coding genes in the D. oleifera mt genome had undergone negative selections. But, the rearrangement of mitochondrial genes has been widely occur among D. oleifera and these observed species. These results will lay the foundation for identifying further evolutionary relationships within Ebenaceae.

Measuring thermodynamic preferences to form non-native conformations in nucleic acids using ultraviolet melting

Thermodynamic preferences to form non-native conformations are crucial for understanding how nucleic acids fold and function. However, they are difficult to measure experimentally because this requires accurately determining the population of minor low-abundance (<10%) conformations in a sea of other conformations. Here, we show that melting experiments enable facile measurements of thermodynamic preferences to adopt nonnative conformations in DNA and RNA. The key to this "delta-melt" approach is to use chemical modifications to render specific minor non-native conformations the major state. The validity and robustness of delta-melt is established for four different non-native conformations under various physiological conditions and sequence contexts through independent measurements of thermodynamic preferences using NMR. Delta-melt is faster relative to NMR, simple, and cost-effective and enables thermodynamic preferences to be measured for exceptionally low-populated conformations. Using delta-melt, we obtained rare insights into conformational cooperativity, obtaining evidence for significant cooperativity (1.0 to 2.5 kcal/mol) when simultaneously forming two adjacent Hoogsteen base pairs. We also measured the thermodynamic preferences to form G-C⁺ and A-T Hoogsteen and A-T base open states for nearly all 16 trinucleotide sequence contexts and found distinct sequence-specific variations on the order of 2 to 3 kcal/mol. This rich landscape of sequence-specific non-native minor conformations in the DNA double helix may help shape the sequence specificity of DNA biochemistry. Thus, melting experiments can now be used to access thermodynamic information regarding regions of the free energy landscape of biomolecules beyond the native folded and unfolded conformations.

Enhanced Codon-Anticodon Interaction at In-Frame UAG Stop Codon Improves the Efficiency of Non-Natural Amino Acid Mutagenesis

The introduction of non-natural amino acids into proteins through the stop codon readthrough methodology has been used to design proteins for diverse applications. However, this method suffers from low yields of the modified protein, as the suppressor tRNA that recognizes the stop codon is unable to compete effectively with release factor 1 (RF1), which terminates translation. We reasoned that a suppressor tRNA with improved interaction with the UAG stop codon on the mRNA will be able to compete more effectively with RF1. To test this idea, we inserted two 2,6-diaminopurine (D) units in the tRNA anticodon stem loop, including one at the third position of the tRNA anticodon. The modified suppressor tRNA could potentially form additional H-bonds between the N²-exocyclic amine of D and the C2 carbonyl group of uracil, thereby enhancing mRNA-tRNA interaction and/or altering tRNA conformation. The stronger interaction at the codon-anticodon interface resulted in improved UAG decoding efficiency and a higher yield of the modified protein containing a non-natural amino acid at multiple sites. Our findings are consistent with the importance of hydrogen bonding and tRNA conformation at the tRNA-mRNA duplex interface during in-frame UAG suppression, which improves protein translation at multiple UAG stop sites. This work provides valuable inputs toward improved non-natural amino acid mutagenesis for creating designer proteins.

Bacterial tRNA landscape revisited

The updated Wobble Hypothesis reasonably explains why some 40 tRNA species are sufficient to decode the 61 amino acid codons of the Universal Genetic Code. However, we still have no clue why eubacteria lack tRNA isoacceptors with ANN anticodons, whereas eukaryotes universally lack eight GNN anticodons, only one of which is also absent in bacteria. Direct tRNA sequencing could resolve the patterns of nucleoside modification that had been driving the divergent evolution in prokaryotes and eukaryotes, but this task will require the development of AI-supported base-callers that can recognize modified nucleosides without any subsequent analytical verification. Our knowledge of the bacterial tRNA landscape is moreover broadened by the recent discovery of antisense tRNAs and tRNA-derived fragments that should be examined in their roles for gene expression, translation, bacterial physiology or metabolism.

Escherichia coli tRNA 2-Selenouridine Synthase (SelU): Elucidation of Substrate Specificity to Understand the Role of S-Geranyl-tRNA in the Conversion of 2-Thio- into 2-Selenouridines in Bacterial tRNA

The bacterial enzyme tRNA 2-selenouridine synthase (SelU) is responsible for the conversion of 5-substituted 2-thiouridine (R5S2U), present in the anticodon of some bacterial tRNAs, into 5-substituted 2-selenouridine (R5Se2U). We have already demonstrated using synthetic RNAs that transformation S2U→Se2U is a two-step process, in which the S2U-RNA is geranylated and the resulting geS2U-RNA is selenated. Currently, the question is how SelU recognizes its substrates and what the cellular pathway of R5S2U→R5Se2U conversion is in natural tRNA. In the study presented here, we characterized the SelU substrate requirements, identified SelU-associated tRNAs and their specific modifications in the wobble position. Finally, we explained the sequence of steps in the selenation of tRNA. The S2U position within the RNA chain, the flanking sequence of the modification, and the length of the RNA substrate, all have a key influence on the recognition by SelU. MST data on the affinity of SelU to individual RNAs confirmed the presumed process. SelU binds the R5S2U-tRNA and then catalyzes its geranylation to the R5geS2U-tRNA, which remains bound to the enzyme and is selenated in the next step of the transformation. Finally, the R5Se2U-tRNA leaves the enzyme and participates in the translation process. The enzyme does not directly catalyze the R5S2U-tRNA selenation and the R5geS2U-tRNA is the intermediate product in the linear sequence of reactions.

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

A versatile tRNA modification-sensitive northern blot method with enhanced performance

The 22 mitochondrial and ∼45 cytosolic tRNAs in human cells contain several dozen different post-transcriptional modified nucleotides such that each carries a unique constellation that complements its function. Many tRNA modifications are linked to altered gene expression, and deficiencies due to mutations in tRNA modification enzymes (TMEs) are responsible for numerous diseases. Easily accessible methods to detect tRNA hypomodifications can facilitate progress in advancing such molecular studies. Our laboratory developed a northern blot method that can quantify relative levels of base modifications on multiple specific tRNAs ∼10 yr ago, which has been used to characterize four different TME deficiencies and is likely further extendable. The assay method depends on differential annealing efficiency of a DNA-oligo probe to the modified versus unmodified tRNA. The signal of this probe is then normalized by a second probe elsewhere on the same tRNA. This positive hybridization in the absence of modification (PHAM) assay has proven useful for i6A37, t6A37, m3C32, and m2,2G26 in multiple laboratories. Yet, over the years we have observed idiosyncratic inconsistency and variability in the assay. Here we document these for some tRNAs and probes and illustrate principles and practices for improved reliability and uniformity in performance. We provide an overview of the method and illustrate benefits of the improved conditions. This is followed by data that demonstrate quantitative validation of PHAM using a TME deletion control, and that nearby modifications can falsely alter the calculated apparent modification efficiency. Finally, we include a calculator tool for matching probe and hybridization conditions.

Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA

Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.

Epitranscriptomic Reprogramming Is Required to Prevent Stress and Damage from Acetaminophen

Epitranscriptomic marks, in the form of enzyme catalyzed RNA modifications, play important gene regulatory roles in response to environmental and physiological conditions. However, little is known with respect to how acute toxic doses of pharmaceuticals influence the epitranscriptome. Here we define how acetaminophen (APAP) induces epitranscriptomic reprogramming and how the writer Alkylation Repair Homolog 8 (Alkbh8) plays a key gene regulatory role in the response. Alkbh8 modifies tRNA selenocysteine (tRNA^Sec) to translationally regulate the production of glutathione peroxidases (Gpx’s) and other selenoproteins, with Gpx enzymes known to play protective roles during APAP toxicity. We demonstrate that APAP increases toxicity and markers of damage, and decreases selenoprotein levels in Alkbh8 deficient mouse livers, when compared to wildtype. APAP also promotes large scale reprogramming of many RNA marks comprising the liver tRNA epitranscriptome including: 5-methoxycarbonylmethyluridine (mcm⁵U), isopentenyladenosine (i⁶A), pseudouridine (Ψ), and 1-methyladenosine (m¹A) modifications linked to tRNA^Sec and many other tRNA’s. Alkbh8 deficiency also leads to wide-spread epitranscriptomic dysregulation in response to APAP, demonstrating that a single writer defect can promote downstream changes to a large spectrum of RNA modifications. Our study highlights the importance of RNA modifications and translational responses to APAP, identifies writers as key modulators of stress responses in vivo and supports the idea that the epitranscriptome may play important roles in responses to pharmaceuticals.

“Superwobbling” and tRNA-34 Wobble and tRNA-37 Anticodon Loop Modifications in Evolution and Devolution of the Genetic Code

The genetic code evolved around the reading of the tRNA anticodon on the primitive ribosome, and tRNA-34 wobble and tRNA-37 modifications coevolved with the code. We posit that EF-Tu, the closing mechanism of the 30S ribosomal subunit, methylation of wobble U34 at the 5-carbon and suppression of wobbling at the tRNA-36 position were partly redundant and overlapping functions that coevolved to establish the code. The genetic code devolved in evolution of mitochondria to reduce the size of the tRNAome (all of the tRNAs of an organism or organelle). “Superwobbling” or four-way wobbling describes a major mechanism for shrinking the mitochondrial tRNAome. In superwobbling, unmodified wobble tRNA-U34 can recognize all four codon wobble bases (A, G, C and U), allowing a single unmodified tRNA-U34 to read a 4-codon box. During code evolution, to suppress superwobbling in 2-codon sectors, U34 modification by methylation at the 5-carbon position appears essential. As expected, at the base of code evolution, tRNA-37 modifications mostly related to the identity of the adjacent tRNA-36 base. TRNA-37 modifications help maintain the translation frame during elongation.

Co-occurrence of m.15992A>G and m.15077G>A Is Associated With a High Penetrance of Maternally Inherited Hypertension in a Chinese Pedigree

BACKGROUND

Mitochondrial DNA (mtDNA) pathogenic variants have been identified to be associated with maternally inherited essential hypertension (MIEH). However, the distinctive clinical features and molecular pathogenesis of MIEH are not fully understood.

METHODS

In this study, we collected a Chinese MIEH family with extraordinary higher penetrance of essential hypertension (88.89%) and early ages of onset (31-40 years old), and performed clinical and genetic characterization for this family. The complete mitochondrial genome of the proband was sequenced and analyzed.

RESULTS

The maternally related members in this family were presented with severe increased blood pressure, left ventricular remodeling, and metabolic abnormalities. Through sequencing the entire mtDNA of the proband and performing systematic analysis of the mtDNA variants with a phylogenic approach, we identified a potentially pathogenic tRNA variant (m.15992A>G in the MT-TP gene) that may account for the MIEH in this family. One nonsynonymous variant (m.15077G>A in the MT-CYB gene) was identified to play a synergistic role with m.15992A>G to cause a high penetrance of MIEH.

CONCLUSIONS

Our results, together with previous findings, have indicated that tRNA pathogenic variants in the mtDNA could act important roles in the pathogenesis of MIEH through reducing mitochondrial translation and disturbing mitochondrial function.

Efficient and Precise Protein Synthesis in a Cell-Free System Using a Set of In Vitro Transcribed tRNAs with Nucleotide Modifications

Reconstitution of a complicated system with a minimal set of components is essential for understanding the mechanisms of how the input is reflected in the output, which is fundamental for further engineering of the corresponding system. We have recently developed a reconstituted cell-free protein synthesis system equipped only with 21 in vitro transcribed tRNAs, one of the minimal systems for understanding the genetic code decoding mechanisms. Introduction of several nucleotide modifications to the transcribed tRNAs showed improvement of both protein synthesis efficiency and its fidelity, suggesting various combinations of tRNAs and their modifications can be evaluated in the developed system. In this chapter, we describe how to prepare this minimal system. Methods for preparing the transcribed tRNAs, their modifications, and the protein production using the set of prepared tRNAs are shown.