Comparative tRNA sequencing and RNA mass spectrometry for surveying tRNA modifications

A Novel tsRNA, m7G-3' tiRNA LysTTT, Promotes Bladder Cancer Malignancy Via Regulating ANXA2 Phosphorylation

Emerging evidence indicates that transfer RNA (tRNA)-derived small RNAs (tsRNAs), originated from tRNA with high abundance RNA modifications, play an important role in many complex physiological and pathological processes. However, the biological functions and regulatory mechanisms of modified tsRNAs in cancer remain poorly understood. Here, it is screened for and confirmed the presence of a novel m7G-modified tsRNA, m7G-3'-tiRNA LysTTT (mtiRL), in a variety of chemical carcinogenesis models by combining small RNA sequencing with an m7G small RNA-modified chip. Moreover, it is found that mtiRL, catalyzed by the tRNA m7G-modifying enzyme mettl1, promotes bladder cancer (BC) malignancy in vitro and in vivo. Mechanistically, mtiRL is found to specifically bind the oncoprotein Annexin A2 (ANXA2) to promote its Tyr24 phosphorylation by enhancing the interactions between ANXA2 and Yes proto-oncogene 1 (Yes1), leading to ANXA2 activation and increased p-ANXA2-Y24 nuclear localization in BC cells. Together, these findings define a critical role for mtiRL and suggest that targeting this novel m7G-modified tsRNA can be an efficient way for to treat BC.

Biosynthesis and function of 7-deazaguanine derivatives in bacteria and phages

SUMMARYDeazaguanine modifications play multifaceted roles in the molecular biology of DNA and tRNA, shaping diverse yet essential biological processes, including the nuanced fine-tuning of translation efficiency and the intricate modulation of codon-anticodon interactions. Beyond their roles in translation, deazaguanine modifications contribute to cellular stress resistance, self-nonself discrimination mechanisms, and host evasion defenses, directly modulating the adaptability of living organisms. Deazaguanine moieties extend beyond nucleic acid modifications, manifesting in the structural diversity of biologically active natural products. Their roles in fundamental cellular processes and their presence in biologically active natural products underscore their versatility and pivotal contributions to the intricate web of molecular interactions within living organisms. Here, we discuss the current understanding of the biosynthesis and multifaceted functions of deazaguanines, shedding light on their diverse and dynamic roles in the molecular landscape of life.

Mapping the tRNA modification landscape of Bartonella henselae Houston I and Bartonella quintana Toulouse

Transfer RNA (tRNA) modifications play a crucial role in maintaining translational fidelity and efficiency, and they may function as regulatory elements in stress response and virulence. Despite their pivotal roles, a comprehensive mapping of tRNA modifications and their associated synthesis genes is still limited, with a predominant focus on free-living bacteria. In this study, we employed a multidisciplinary approach, incorporating comparative genomics, mass spectrometry, and next-generation sequencing, to predict the set of tRNA modification genes responsible for tRNA maturation in two intracellular pathogens-Bartonella henselae Houston I and Bartonella quintana Toulouse, which are causative agents of cat-scratch disease and trench fever, respectively. This analysis presented challenges, particularly because of host RNA contamination, which served as a potential source of error. However, our approach predicted 26 genes responsible for synthesizing 23 distinct tRNA modifications in B. henselae and 22 genes associated with 23 modifications in B. quintana. Notably, akin to other intracellular and symbiotic bacteria, both Bartonella species have undergone substantial reductions in tRNA modification genes, mostly by simplifying the hypermodifications present at positions 34 and 37. Bartonella quintana exhibited the additional loss of four modifications and these were linked to examples of gene decay, providing snapshots of reductive evolution.

The modification landscape of P. aeruginosa tRNAs

RNA modifications have a substantial impact on tRNA function, with modifications in the anticodon loop contributing to translational fidelity and modifications in the tRNA core impacting structural stability. In bacteria, tRNA modifications are crucial for responding to stress and regulating the expression of virulence factors. Although tRNA modifications are well-characterized in a few model organisms, our knowledge of tRNA modifications in human pathogens, such as Pseudomonas aeruginosa, remains limited. Here we leveraged two orthogonal approaches to build a reference landscape of tRNA modifications in E. coli, which enabled us to identify similar modifications in P. aeruginosa. Our analysis revealed a substantial degree of conservation between the two organisms, while also uncovering potential sites of tRNA modification in P. aeruginosa tRNAs that are not present in E. coli. The mutational signature at one of these sites, position 46 of tRNAGln1(UUG) is dependent on the P. aeruginosa homolog of TapT, the enzyme responsible for the 3-(3-amino-3-carboxypropyl) uridine (acp3U) modification. Identifying which modifications are present on different tRNAs will uncover the pathways impacted by the different tRNA modifying enzymes, some of which play roles in determining virulence and pathogenicity.

Mapping the tRNA Modification Landscape of Bartonella henselae Houston I and Bartonella quintana Toulouse

Transfer RNA (tRNA) modifications play a crucial role in maintaining translational fidelity and efficiency, and they may function as regulatory elements in stress response and virulence. Despite their pivotal roles, a comprehensive mapping of tRNA modifications and their associated synthesis genes is still limited, with a predominant focus on free-living bacteria. In this study, we employed a multidisciplinary approach, incorporating comparative genomics, mass spectrometry, and next-generation sequencing, to predict the set of tRNA modification genes responsible for tRNA maturation in two intracellular pathogens- Bartonella henselae Houston I and Bartonella quintana Toulouse, which are causative agents of cat-scratch disease and trench fever, respectively. This analysis presented challenges, particularly because of host RNA contamination, which served as a potential source of error. However, our approach predicted 26 genes responsible for synthesizing 23 distinct tRNA modifications in B. henselae and 22 genes associated with 23 modifications in B. quintana . Notably, akin to other intracellular and symbiotic bacteria, both Bartonella species have undergone substantial reductions in tRNA modification genes, mostly by simplifying the hypermodifications present at positions 34 and 37. B. quintana exhibited the additional loss of four modifications and these were linked to examples of gene decay, providing snapshots of reductive evolution.

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.

Genome-Wide Profiling of tRNA Using an Unexplored Reverse Transcriptase with High Processivity

Monitoring the dynamic changes of cellular tRNA pools is challenging, due to the extensive post-transcriptional modifications of individual species. The most critical component in tRNAseq is a processive reverse transcriptase (RT) that can read through each modification with high efficiency. Here we show that the recently developed group-II intron RT Induro has the processivity and efficiency necessary to profile tRNA dynamics. Using our Induro-tRNAseq, simpler and more comprehensive than the best methods to date, we show that Induro progressively increases readthrough of tRNA over time and that the mechanism of increase is selective removal of RT stops, without altering the misincorporation frequency. We provide a parallel dataset of the misincorporation profile of Induro relative to the related TGIRT RT to facilitate the prediction of non-annotated modifications. We report an unexpected modification profile among human proline isoacceptors, absent from mouse and lower eukaryotes, that indicates new biology of decoding proline codons.

Integrative Approach to Probe Alternative Redox Mechanisms in RNA Modifications

RNA modifications found in most RNAs, particularly in tRNAs and rRNAs, reveal an abundance of chemical alterations of nucleotides. Over 150 distinct RNA modifications are known, emphasizing a remarkable diversity of chemical moieties in RNA molecules. These modifications play pivotal roles in RNA maturation, structural integrity, and the fidelity and efficiency of translation processes. The catalysts responsible for these modifications are RNA-modifying enzymes that use a striking array of chemistries to directly influence the chemical landscape of RNA. This diversity is further underscored by instances where the same modification is introduced by distinct enzymes that use unique catalytic mechanisms and cofactors across different domains of life. This phenomenon of convergent evolution highlights the biological importance of RNA modification and the vast potential within the chemical repertoire for nucleotide alteration. While shared RNA modifications can hint at conserved enzymatic pathways, a major bottleneck is to identify alternative routes within species that possess a modified RNA but are devoid of known RNA-modifying enzymes. To address this challenge, a combination of bioinformatic and experimental strategies proves invaluable in pinpointing new genes responsible for RNA modifications. This integrative approach not only unveils new chemical insights but also serves as a wellspring of inspiration for biocatalytic applications and drug design. In this Account, we present how comparative genomics and genome mining, combined with biomimetic synthetic chemistry, biochemistry, and anaerobic crystallography, can be judiciously implemented to address unprecedented and alternative chemical mechanisms in the world of RNA modification. We illustrate these integrative methodologies through the study of tRNA and rRNA modifications, dihydrouridine, 5-methyluridine, queuosine, 8-methyladenosine, 5-carboxymethylamino-methyluridine, or 5-taurinomethyluridine, each dependent on a diverse array of redox chemistries, often involving organic compounds, organometallic complexes, and metal coenzymes. We explore how vast genome and tRNA databases empower comparative genomic analyses and enable the identification of novel genes that govern RNA modification. Subsequently, we describe how the isolation of a stable reaction intermediate can guide the synthesis of a biomimetic to unveil new enzymatic pathways. We then discuss the usefulness of a biochemical "shunt" strategy to study catalytic mechanisms and to directly visualize reactive intermediates bound within active sites. While we primarily focus on various RNA-modifying enzymes studied in our laboratory, with a particular emphasis on the discovery of a SAM-independent methylation mechanism, the strategies and rationale presented herein are broadly applicable for the identification of new enzymes and the elucidation of their intricate chemistries. This Account offers a comprehensive glimpse into the evolving landscape of RNA modification research and highlights the pivotal role of integrated approaches to identify novel enzymatic pathways.

Transfer RNAs as dynamic and critical regulators of cancer progression

Transfer RNAs (tRNAs) have been historically viewed as non-dynamic adaptors that decode the genetic code into proteins. Recent work has uncovered dynamic regulatory roles for these fascinating molecules. Advances in tRNA detection methods have revealed that specific tRNAs can become modulated upon DNA copy number and chromatin alterations and can also be perturbed by oncogenic signalling and transcriptional regulators in cancer cells or the tumour microenvironment. Such alterations in the levels of specific tRNAs have been shown to causally impact cancer progression, including metastasis. Moreover, sequencing methods have identified tRNA-derived small RNAs that influence various aspects of cancer progression, such as cell proliferation and invasion, and could serve as diagnostic and prognostic biomarkers or putative therapeutic targets in various cancers. Finally, there is accumulating evidence, including from genetic models, that specific tRNA synthetases - the enzymes responsible for charging tRNAs with amino acids - can either promote or suppress tumour formation. In this Review, we provide an overview of how deregulation of tRNAs influences cancer formation and progression.

The tRNA methyltransferase TrmB is critical for Acinetobacter baumannii stress responses and pulmonary infection

IMPORTANCE

As deficiencies in tRNA modifications have been linked to human diseases such as cancer and diabetes, much research has focused on the modifications' impacts on translational regulation in eukaryotes. However, the significance of tRNA modifications in bacterial physiology remains largely unexplored. In this paper, we demonstrate that the m7G tRNA methyltransferase TrmB is crucial for a top-priority pathogen, Acinetobacter baumannii, to respond to stressors encountered during infection, including oxidative stress, low pH, and iron deprivation. We show that loss of TrmB dramatically attenuates a murine pulmonary infection. Given the current efforts to use another tRNA methyltransferase, TrmD, as an antimicrobial therapeutic target, we propose that TrmB, and other tRNA methyltransferases, may also be viable options for drug development to combat multidrug-resistant A. baumannii.

A tRNA modification in Mycobacterium tuberculosis facilitates optimal intracellular growth

Diverse chemical modifications fine-tune the function and metabolism of tRNA. Although tRNA modification is universal in all kingdoms of life, profiles of modifications, their functions, and physiological roles have not been elucidated in most organisms including the human pathogen, Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis. To identify physiologically important modifications, we surveyed the tRNA of Mtb, using tRNA sequencing (tRNA-seq) and genome-mining. Homology searches identified 23 candidate tRNA modifying enzymes that are predicted to create 16 tRNA modifications across all tRNA species. Reverse transcription-derived error signatures in tRNA-seq predicted the sites and presence of nine modifications. Several chemical treatments prior to tRNA-seq expanded the number of predictable modifications. Deletion of Mtb genes encoding two modifying enzymes, TruB and MnmA, eliminated their respective tRNA modifications, validating the presence of modified sites in tRNA species. Furthermore, the absence of mnmA attenuated Mtb growth in macrophages, suggesting that MnmA-dependent tRNA uridine sulfation contributes to Mtb intracellular growth. Our results lay the foundation for unveiling the roles of tRNA modifications in Mtb pathogenesis and developing new therapeutics against tuberculosis.

Restoration of mitochondrial function through activation of hypomodified tRNAs with pathogenic mutations associated with mitochondrial diseases

Mutations in mitochondrial (mt-)tRNAs frequently cause mitochondrial dysfunction. Mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes (MELAS), and myoclonus epilepsy associated with ragged red fibers (MERRF) are major clinical subgroups of mitochondrial diseases caused by pathogenic point mutations in tRNA genes encoded in mtDNA. We previously reported a severe reduction in the frequency of 5-taurinomethyluridine (τm5U) and its 2-thiouridine derivative (τm5s2U) in the anticodons of mutant mt-tRNAs isolated from the cells of patients with MELAS and MERRF, respectively. The hypomodified tRNAs fail to decode cognate codons efficiently, resulting in defective translation of respiratory chain proteins in mitochondria. To restore the mitochondrial activity of MELAS patient cells, we overexpressed MTO1, a τm5U-modifying enzyme, in patient-derived myoblasts. We used a newly developed primer extension method and showed that MTO1 overexpression almost completely restored the τm5U modification of the MELAS mutant mt-tRNALeu(UUR). An increase in mitochondrial protein synthesis and oxygen consumption rate suggested that the mitochondrial function of MELAS patient cells can be activated by restoring the τm5U of the mutant tRNA. In addition, we confirmed that MTO1 expression restored the τm5s2U of the mutant mt-tRNALys in MERRF patient cells. These findings pave the way for epitranscriptomic therapies for mitochondrial diseases.

Discovering RNA modification enzymes using a comparative genomics approach

Identifying RNA modification enzymes is critical for understanding the biogenesis and function of RNA modification. Among several approaches that enable the identification of RNA modification enzymes, comparative genomics has become particularly useful due to the expanding availability of genomic DNA and annotation data. Here, a detailed protocol for carrying out a computational comparative genomics approach for the discovery of RNA modification enzymes is presented. An illustrative example of the utility of this approach in the discovery of AcpA, an acetyltransferase that synthesizes the newly discovered modification, acacp3U is also provided. This computational framework has applications for the identification of genes involved in other cellular processes.

DpCoA tagSeq: Barcoding dpCoA-Capped RNA for Direct Nanopore Sequencing via Maleimide-Thiol Reaction

Recent discoveries of noncanonical RNA caps, such as nicotinamide adenine dinucleotide (NAD+) and 3'-dephospho-coenzyme A (dpCoA), have expanded our knowledge of RNA caps. Although dpCoA has been known to cap RNAs in various species, the identities of its capped RNAs (dpCoA-RNAs) remained unknown. To fill this gap, we developed a method called dpCoA tagSeq, which utilized a thiol-reactive maleimide group to label dpCoA cap with a tag RNA serving as the 5' barcode. The barcoded RNAs were isolated using a complementary DNA strand of the tag RNA prior to direct sequencing by nanopore technology. Our validation experiments with model RNAs showed that dpCoA-RNA was efficiently tagged and captured using this protocol. To confirm that the tagged RNAs are capped by dpCoA and no other thiol-containing molecules, we used a pyrophosphatase NudC to degrade the dpCoA cap to adenosine monophosphate (AMP) moiety before performing the tagSeq protocol. We identified 44 genes that transcribe dpCoA-RNAs in mouse liver, demonstrating the method's effectiveness in identifying and characterizing the capped RNAs. This strategy provides a viable approach to identifying dpCoA-RNAs that allows for further functional investigations of the cap.

A tRNA modification in Mycobacterium tuberculosis facilitates optimal intracellular growth

Diverse chemical modifications fine-tune the function and metabolism of tRNA. Although tRNA modification is universal in all kingdoms of life, profiles of modifications, their functions, and physiological roles have not been elucidated in most organisms including the human pathogen, Mycobacterium tuberculosis ( Mtb ), the causative agent of tuberculosis. To identify physiologically important modifications, we surveyed the tRNA of Mtb , using tRNA sequencing (tRNA-seq) and genome-mining. Homology searches identified 23 candidate tRNA modifying enzymes that are predicted to create 16 tRNA modifications across all tRNA species. Reverse transcription-derived error signatures in tRNA-seq predicted the sites and presence of 9 modifications. Several chemical treatments prior to tRNA-seq expanded the number of predictable modifications. Deletion of Mtb genes encoding two modifying enzymes, TruB and MnmA, eliminated their respective tRNA modifications, validating the presence of modified sites in tRNA species. Furthermore, the absence of mnmA attenuated Mtb growth in macrophages, suggesting that MnmA-dependent tRNA uridine sulfation contributes to Mtb intracellular growth. Our results lay the foundation for unveiling the roles of tRNA modifications in Mtb pathogenesis and developing new therapeutics against tuberculosis.

Navigating the pitfalls of mapping DNA and RNA modifications

Chemical modifications to nucleic acids occur across the kingdoms of life and carry important regulatory information. Reliable high-resolution mapping of these modifications is the foundation of functional and mechanistic studies, and recent methodological advances based on next-generation sequencing and long-read sequencing platforms are critical to achieving this aim. However, mapping technologies may have limitations that sometimes lead to inconsistent results. Some of these limitations are technical in nature and specific to certain types of technology. Here, however, we focus on common (yet not always widely recognized) pitfalls that are shared among frequently used mapping technologies and discuss strategies to help technology developers and users mitigate their effects. Although the emphasis is primarily on DNA modifications, RNA modifications are also discussed.

Fluorescent labeling of tRNA for rapid kinetic interaction studies with tRNA-binding proteins

Transfer RNA (tRNA) plays a critical role during translation and interacts with numerous proteins during its biogenesis, functional cycle and degradation. In particular, tRNA is extensively post-transcriptionally modified by various tRNA modifying enzymes which each target a specific nucleotide at different positions within tRNAs to introduce different chemical modifications. Fluorescent assays can be used to study the interaction between a protein and tRNA. Moreover, rapid mixing fluorescence stopped-flow assays provide insights into the kinetics of the tRNA-protein interaction in order to elucidate the tRNA binding mechanism for the given protein. A prerequisite for these studies is a fluorescently labeled molecule, such as fluorescent tRNA, wherein a change in fluorescence occurs upon protein binding. In this chapter, we discuss the utilization of tRNA modifications in order to introduce fluorophores at particular positions within tRNAs. Particularly, we focus on in vitro thiolation of a uridine at position 8 within tRNAs using the tRNA modification enzyme ThiI, followed by labeling of the thiol group with fluorescein. As such, this fluorescently labeled tRNA is primarily unmodified, with the exception of the thiolation modification to which the fluorophore is attached, and can be used as a substrate to study the binding of different tRNA-interacting factors. Herein, we discuss the example of studying the tRNA binding mechanism of the tRNA modifying enzymes TrmB and DusA using internally fluorescein-labeled tRNA.

Middle-out sequence confirmation of CRISPR/Cas9 single guide RNA (sgRNA) using DNA primers and ribonuclease T1 digestion

Accurate sequencing of single guide RNAs (sgRNAs) for CRISPR/Cas9 genome editing is critical for patient safety, as the sgRNA guides the Cas9 nuclease to target site-specific cleavages in DNA. An approach to fully sequence sgRNA using protective DNA primers followed by ribonuclease (RNase) T1 digestion was developed to facilitate the analysis of these larger molecules by hydrophilic interaction liquid chromatography coupled with high-resolution mass spectrometry (HILIC-HRMS). Without RNase digestion, top-down mass spectrometry alone struggles to properly fragment precursor ions in large RNA oligonucleotides to provide confidence in sequence coverage. With RNase T1 digestion of these larger oligonucleotides, however, bottom-up analysis cannot confirm full sequence coverage due to the presence of short, redundant digestion products. By combining primer protection with RNase T1 digestion, digestion products are large enough to prevent redundancy and small enough to provide base resolution by tandem mass spectrometry to allow for full sgRNA sequence coverage. An investigation into the general requirements for adequate primer protection of specific regions of the RNA was conducted, followed by the development of a generic protection and digestion strategy that may be applied to different sgRNA sequences. This middle-out technique has the potential to expedite accurate sequence confirmation of chemically modified sgRNA oligonucleotides.

Functional analysis of tRNA modification enzymes using mutational profiling

Transfer RNA (tRNA) delivers amino acids to the ribosome and functions as an essential adapter molecule for decoding codons on the messenger RNA (mRNA) during protein synthesis. Before attaining their proper activity, tRNAs undergo multiple post-transcriptional modifications with highly diversified roles such as stabilization of the tRNA structure, recognition of aminoacyl tRNA synthetases, precise codon-anticodon recognition, support of viral replication and onset of immune responses. The synthesis of the majority of modified nucleosides is catalyzed by a site-specific tRNA modification enzyme. This chapter provides a detailed protocol for using mutational profiling to analyze the enzymatic function of a tRNA methyltransferase in a high-throughput manner. In a previous study, we took tRNA m1A22 methyltransferase TrmK from Geobacillus stearothermophilus as a model tRNA methyltransferase and applied this protocol to gain mechanistic insights into how TrmK recognizes the substrate tRNAs. In theory, this protocol can be used unaltered for studying enzymes that catalyze modifications at the Watson-Crick face such as 1-methyladenosine (m1A), 3-methylcytosine (m3C), 3-methyluridine (m3U), 1-methylguanosine (m1G), and N2,N2-dimethylguanosine (m22G).

choros: correction of sequence-based biases for accurate quantification of ribosome profiling data

Ribosome profiling quantifies translation genome-wide by sequencing ribosome-protected fragments, or footprints. Its single-codon resolution allows identification of translation regulation, such as ribosome stalls or pauses, on individual genes. However, enzyme preferences during library preparation lead to pervasive sequence artifacts that obscure translation dynamics. Widespread over- and under-representation of ribosome footprints can dominate local footprint densities and skew estimates of elongation rates by up to five fold. To address these biases and uncover true patterns of translation, we present choros, a computational method that models ribosome footprint distributions to provide bias-corrected footprint counts. choros uses negative binomial regression to accurately estimate two sets of parameters: (i) biological contributions from codon-specific translation elongation rates; and (ii) technical contributions from nuclease digestion and ligation efficiencies. We use these parameter estimates to generate bias correction factors that eliminate sequence artifacts. Applying choros to multiple ribosome profiling datasets, we are able to accurately quantify and attenuate ligation biases to provide more faithful measurements of ribosome distribution. We show that a pattern interpreted as pervasive ribosome pausing near the beginning of coding regions is likely to arise from technical biases. Incorporating choros into standard analysis pipelines will improve biological discovery from measurements of translation.

"Transfer" of power: The intersection of DNA virus infection and tRNA biology

Transfer RNAs (tRNAs) are at the heart of the molecular biology central dogma, functioning to decode messenger RNAs into proteins. As obligate intracellular parasites, viruses depend on the host translation machinery, including host tRNAs. Thus, the ability of a virus to fine-tune tRNA expression elicits the power to impact the outcome of infection. DNA viruses commonly upregulate the output of RNA polymerase III (Pol III)-dependent transcripts, including tRNAs. Decades after these initial discoveries we know very little about how mature tRNA pools change during viral infection, as tRNA sequencing methodology has only recently reached proficiency. Here, we review perturbation of tRNA biogenesis by DNA virus infection, including an emerging player called tRNA-derived fragments (tRFs). We discuss how tRNA dysregulation shifts the power landscape between the host and virus, highlighting the potential for tRNA-based antivirals as a future therapeutic.

ADAR RNA editing on antisense RNAs results in apparent U-to-C base changes on overlapping sense transcripts

Despite hundreds of RNA modifications described to date, only RNA editing results in a change in the nucleotide sequence of RNA molecules compared to the genome. In mammals, two kinds of RNA editing have been described so far, adenosine to inosine (A-to-I) and cytidine to uridine (C-to-U) editing. Recent improvements in RNA sequencing technologies have led to the discovery of a continuously growing number of editing sites. These methods are powerful but not error-free, making routine validation of newly-described editing sites necessary. During one of these validations on DDX58 mRNA, along with A-to-I RNA editing sites, we encountered putative U-to-C editing. These U-to-C edits were present in several cell lines and appeared regulated in response to specific environmental stimuli. The same findings were also observed for the human long intergenic non-coding RNA p21 (hLincRNA-p21). A more in-depth analysis revealed that putative U-to-C edits result from A-to-I editing on overlapping antisense RNAs that are transcribed from the same loci. Such editing events, occurring on overlapping genes transcribed in opposite directions, have recently been demonstrated to be immunogenic and have been linked with autoimmune and immune-related diseases. Our findings, also confirmed by deep transcriptome data, demonstrate that such loci can be recognized simply through the presence of A-to-I and U-to-C mismatches within the same locus, reflective A-to-I editing both in the sense-oriented transcript and in the cis-natural antisense transcript (cis-NAT), implying that such clusters could be a mark of functionally relevant ADAR1 editing events.

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 (m¹A22) and tRNA (m¹A22) methyltransferase (TrmK) for further analysis. Relative quantification of m¹A22 levels in 59 tRNA transcripts by mutational profiling reveals that TrmK selectively methylates a subset of tRNAs. Using 240 variants of tRNA^Leu 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 tRNA^Leu 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.

Developing New Tools to Fight Human Pathogens: A Journey through the Advances in RNA Technologies

A long scientific journey has led to prominent technological advances in the RNA field, and several new types of molecules have been discovered, from non-coding RNAs (ncRNAs) to riboswitches, small interfering RNAs (siRNAs) and CRISPR systems. Such findings, together with the recognition of the advantages of RNA in terms of its functional performance, have attracted the attention of synthetic biologists to create potent RNA-based tools for biotechnological and medical applications. In this review, we have gathered the knowledge on the connection between RNA metabolism and pathogenesis in Gram-positive and Gram-negative bacteria. We further discuss how RNA techniques have contributed to the building of this knowledge and the development of new tools in synthetic biology for the diagnosis and treatment of diseases caused by pathogenic microorganisms. Infectious diseases are still a world-leading cause of death and morbidity, and RNA-based therapeutics have arisen as an alternative way to achieve success. There are still obstacles to overcome in its application, but much progress has been made in a fast and effective manner, paving the way for the solid establishment of RNA-based therapies in the future.

Sequential action of a tRNA base editor in conversion of cytidine to pseudouridine

Post-transcriptional RNA editing modulates gene expression in a condition-dependent fashion. We recently discovered C-to-Ψ editing in Vibrio cholerae tRNA. Here, we characterize the biogenesis, regulation, and functions of this previously undescribed RNA editing process. We show that an enzyme, TrcP, mediates the editing of C-to-U followed by the conversion of U to Ψ, consecutively. AlphaFold-2 predicts that TrcP consists of two globular domains (cytidine deaminase and pseudouridylase) and a long helical domain. The latter domain tethers tRNA substrates during both the C-to-U editing and pseudouridylation, likely enabling a substrate channeling mechanism for efficient catalysis all the way to the terminal product. C-to-Ψ editing both requires and suppresses other modifications, creating an interdependent network of modifications in the tRNA anticodon loop that facilitates coupling of tRNA modification states to iron availability. Our findings provide mechanistic insights into an RNA editing process that likely promotes environmental adaptation.

Higher-Energy Collisional Dissociation Mass Spectral Networks for the Rapid, Semi-automated Characterization of Known and Unknown Ribonucleoside Modifications

Higher-energy collisional dissociation (HCD) of modified ribonucleosides generates characteristic and highly reproducible nucleoside-specific tandem mass spectra (MS/MS). Here, we demonstrate the capability of HCD spectra in combination with spectral matching for the semi-automated characterization of ribonucleosides. This process involved the generation of an HCD spectral library and the establishment of a mass spectral network for rapid detection with high sensitivity and specificity in a retention time-independent fashion. Systematic spectral matching analysis of the MS/MS spectra of tRNA hydrolysates from different organisms has helped us to uncover evidence for the existence of novel ribonucleoside modifications such as s²Cm and OHyW-14. Such an untargeted label-free approach has the potential to be integrated with other methods, including those that use isotope labeling, to simplify the characterization of unknown modified ribonucleosides. These findings suggest the compilation of a universal spectral network, for the characterization of known and unknown ribonucleosides, could accelerate discoveries in the epitranscriptome.

Arsenite toxicity is regulated by queuine availability and oxidation-induced reprogramming of the human tRNA epitranscriptome

Cells respond to environmental stress by regulating gene expression at the level of both transcription and translation. The ∼50 modified ribonucleotides of the human epitranscriptome contribute to the latter, with mounting evidence that dynamic regulation of transfer RNA (tRNA) wobble modifications leads to selective translation of stress response proteins from codon-biased genes. Here we show that the response of human hepatocellular carcinoma cells to arsenite exposure is regulated by the availability of queuine, a micronutrient and essential precursor to the wobble modification queuosine (Q) on tRNAs reading GUN codons. Among oxidizing and alkylating agents at equitoxic concentrations, arsenite exposure caused an oxidant-specific increase in Q that correlated with up-regulation of proteins from codon-biased genes involved in energy metabolism. Limiting queuine increased arsenite-induced cell death, altered translation, increased reactive oxygen species levels, and caused mitochondrial dysfunction. In addition to demonstrating an epitranscriptomic facet of arsenite toxicity and response, our results highlight the links between environmental exposures, stress tolerance, RNA modifications, and micronutrients.

Genome-wide analysis of the in vivo tRNA structurome reveals RNA structural and modification dynamics under heat stress

RNA structure plays roles in myriad cellular events including transcription, translation, and RNA processing. Genome-wide analyses of RNA secondary structure in vivo by chemical probing have revealed critical structural features of mRNAs and long ncRNAs. Here, we examine the in vivo secondary structure of a small RNA class, tRNAs. Study of tRNA structure is challenging because tRNAs are heavily modified and strongly structured. We introduce "tRNA structure-seq," a new workflow that accurately determines in vivo secondary structures of tRNA. The workflow combines dimethyl sulfate (DMS) probing, ultra-processive RT, and mutational profiling (MaP), which provides mutations opposite DMS and natural modifications thereby allowing multiple modifications to be identified in a single read. We applied tRNA structure-seq to E. coli under control and stress conditions. A leading folding algorithm predicts E. coli tRNA structures with only ∼80% average accuracy from sequence alone. Strikingly, tRNA structure-seq, by providing experimental restraints, improves structure prediction under in vivo conditions to ∼95% accuracy, with more than 14 tRNAs predicted completely correctly. tRNA structure-seq also quantifies the relative levels of tRNAs and their natural modifications at single nucleotide resolution, as validated by LC-MS/MS. Our application of tRNA structure-seq yields insights into tRNA structure in living cells, revealing that it is not immutable but has dynamics, with partial unfolding of secondary and tertiary tRNA structure under heat stress that is correlated with a loss of tRNA abundance. This method is applicable to other small RNAs, including those with natural modifications and highly structured regions.

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.

TRMT6/61A-dependent base methylation of tRNA-derived fragments regulates gene-silencing activity and the unfolded protein response in bladder cancer

RNA modifications are important regulatory elements of RNA functions. However, most genome-wide mapping of RNA modifications has focused on messenger RNAs and transfer RNAs, but such datasets have been lacking for small RNAs. Here we mapped N1-methyladenosine (m1A) in the cellular small RNA space. Benchmarked with synthetic m1A RNAs, our workflow identified specific groups of m1A-containing small RNAs, which are otherwise disproportionally under-represented. In particular, 22-nucleotides long 3' tRNA-fragments are highly enriched for TRMT6/61A-dependent m1A located within the seed region. TRMT6/61A-dependent m1A negatively affects gene silencing by tRF-3s. In urothelial carcinoma of the bladder, where TRMT6/61A is over-expressed, higher m1A modification on tRFs is detected, correlated with a dysregulation of tRF targetome. Lastly, TRMT6/61A regulates tRF-3 targets involved in unfolded protein response. Together, our results reveal a mechanism of regulating gene expression via base modification of small RNA.

Exploring the expanding universe of small RNAs

The world of small noncoding RNAs (sncRNAs) is ever-expanding, from small interfering RNA, microRNA and Piwi-interacting RNA to the recently emerging non-canonical sncRNAs derived from longer structured RNAs (for example, transfer, ribosomal, Y, small nucleolar, small nuclear and vault RNAs), showing distinct biogenesis and functional principles. Here we discuss recent tools for sncRNA identification, caveats in sncRNA expression analysis and emerging methods for direct sequencing of sncRNAs and systematic mapping of RNA modifications that are integral to their function.

tRNA modification dynamics from individual organisms to metaepitranscriptomics of microbiomes

tRNA is the most extensively modified RNA in cells. On average, a bacterial tRNA contains 8 modifications per molecule and a eukaryotic tRNA contains 13 modifications per molecule. Recent studies reveal that tRNA modifications are highly dynamic and respond extensively to environmental conditions. Functions of tRNA modification dynamics include enhanced, on-demand decoding of specific codons in response genes and regulation of tRNA fragment biogenesis. This review summarizes recent advances in the studies of tRNA modification dynamics in biological processes, tRNA modification erasers, and human-associated bacteria. Furthermore, we use the term "metaepitranscriptomics" to describe the potential and approach of tRNA modification studies in natural biological communities such as microbiomes. tRNA is highly modified in cells, and tRNA modifications respond extensively to environmental conditions to enhance translation of specific genes and produce tRNA fragments on demand. We review recent advances in tRNA sequencing methods, tRNA modification dynamics in biological processes, and tRNA modification studies in natural communities such as the microbiomes.

Site-Specific Profiling of 4-Thiouridine Across Transfer RNA Genes in Escherichia coli

The transfer RNA (tRNA) modification 4-thiouridine (s4U) acts as a near-ultraviolet (UVA) radiation sensor in Escherichia coli (E. coli), where it induces a growth delay upon exposure to the UVA radiation (∼310-400 nm). Herein, we report sequencing methodology for site-specific profiling of s4U modification in E. coli tRNAs. Upon the addition of iodoacetamide (IA) or iodoacetyl-PEG2-biotin (BIA), the nucleophilic sulfur of s4U forms a reaction product that is extensively characterized by liquid chromatography-mass spectrometry (LC-MS/MS) analysis. This method is readily applied to the alkylation of natively occurring s4U on E. coli tRNA. Next-generation sequencing of BIA-treated tRNA from E. coli revealed misincorporations at position 8 in 19 of the 20 amino acid tRNA species. Alternatively, tRNA from the ΔthiI strain, which cannot introduce the s4U modification, does not exhibit any misincorporation at the corresponding positions, directly linking the base transitions and the tRNA modification. Independently, the s4U modification on E. coli tRNA was further validated by LC-MS/MS sequencing. Nuclease digestion of wild-type and deletion strains E. coli tRNA with RNase T1 generated smaller s4U/U containing fragments that could be analyzed by MS/MS analysis for modification assignment. Furthermore, RNase T1 digestion of tRNAs treated either with IA or BIA showed the specificity of iodoacetamide reagents toward s4U in the context of complex tRNA modifications. Overall, these results demonstrate the utility of the alkylation of s4U in the site-specific profiling of the modified base in native cellular tRNA.

Dynamic changes in tRNA modifications and abundance during T cell activation

The tRNA pool determines the efficiency, throughput, and accuracy of translation. Previous studies have identified dynamic changes in the tRNA (transfer RNA) supply and mRNA (messenger RNA) demand during cancerous proliferation. Yet dynamic changes may also occur during physiologically normal proliferation, and these are less well characterized. We examined the tRNA and mRNA pools of T cells during their vigorous proliferation and differentiation upon triggering their antigen receptor. We observed a global signature of switch in demand for codons at the early proliferation phase of the response, accompanied by corresponding changes in tRNA expression levels. In the later phase, upon differentiation, the response of the tRNA pool relaxed back to the basal level, potentially restraining excessive proliferation. Sequencing of tRNAs allowed us to evaluate their diverse base-modifications. We found that two types of tRNA modifications, wybutosine and ms²t6A, are reduced dramatically during T cell activation. These modifications occur in the anticodon loops of two tRNAs that decode "slippery codons," which are prone to ribosomal frameshifting. Attenuation of these frameshift-protective modifications is expected to increase the potential for proteome-wide frameshifting during T cell proliferation. Indeed, human cell lines deleted of a wybutosine writer showed increased ribosomal frameshifting, as detected with an HIV gag-pol frameshifting site reporter. These results may explain HIV's specific tropism toward proliferating T cells since it requires ribosomal frameshift exactly on the corresponding codon for infection. The changes in tRNA expression and modifications uncover a layer of translation regulation during T cell proliferation and expose a potential tradeoff between cellular growth and translation fidelity.

Endogenous tRNA-derived small RNA (tRF3-Thr-AGT) inhibits ZBP1/NLRP3 pathway-mediated cell pyroptosis to attenuate acute pancreatitis (AP)

Endogenous transfer RNA‐derived small RNAs (tsRNAs) are newly identified RNAs that are closely associated with the pathogenesis of multiple diseases, but the involvement of tsRNAs in regulating acute pancreatitis (AP) development has not been reported. In this study, we screened out a novel tsRNA, tRF3‐Thr‐AGT, that was aberrantly downregulated in the acinar cell line AR42J treated with sodium taurocholate (STC) and the pancreatic tissues of STC‐induced AP rat models. In addition, STC treatment suppressed cell viability, induced pyroptotic cell death and cellular inflammation in AP models in vitro and in vivo. Overexpression of tRF3‐Thr‐AGT partially reversed STC‐induced detrimental effects on the AR42J cells. Next, Z‐DNA‐binding protein 1 (ZBP1) was identified as the downstream target of tRF3‐Thr‐AGT. Interestingly, upregulation of tRF3‐Thr‐AGT suppressed NOD‐like receptor protein 3 (NLRP3)‐mediated pyroptotic cell death in STC‐treated AR42J cells via degrading ZBP1. Moreover, the effects of tRF3‐Thr‐AGT overexpression on cell viability and inflammation in AR42J cells were abrogated by upregulating ZBP1 and NLRP3. Collectively, our data indicated that tRF3‐Thr‐AGT suppressed ZBP1 expressions to restrain NLRP3‐mediated pyroptotic cell death and inflammation in AP models. This study, for the first time, identified the role and potential underlying mechanisms by which tRF3‐Thr‐AGT regulated AP pathogenesis.

Adaptor Molecules Epitranscriptome Reprograms Bacterial Pathogenicity

The strong decoration of tRNAs with post-transcriptional modifications provides an unprecedented adaptability of this class of non-coding RNAs leading to the regulation of bacterial growth and pathogenicity. Accumulating data indicate that tRNA post-transcriptional modifications possess a central role in both the formation of bacterial cell wall and the modulation of transcription and translation fidelity, but also in the expression of virulence factors. Evolutionary conserved modifications in tRNA nucleosides ensure the proper folding and stability redounding to a totally functional molecule. However, environmental factors including stress conditions can cause various alterations in tRNA modifications, disturbing the pathogen homeostasis. Post-transcriptional modifications adjacent to the anticodon stem-loop, for instance, have been tightly linked to bacterial infectivity. Currently, advances in high throughput methodologies have facilitated the identification and functional investigation of such tRNA modifications offering a broader pool of putative alternative molecular targets and therapeutic avenues against bacterial infections. Herein, we focus on tRNA epitranscriptome shaping regarding modifications with a key role in bacterial infectivity including opportunistic pathogens of the human microbiome.

RNA Modifications in Pathogenic Bacteria: Impact on Host Adaptation and Virulence

RNA modifications are involved in numerous biological processes and are present in all RNA classes. These modifications can be constitutive or modulated in response to adaptive processes. RNA modifications play multiple functions since they can impact RNA base-pairings, recognition by proteins, decoding, as well as RNA structure and stability. However, their roles in stress, environmental adaptation and during infections caused by pathogenic bacteria have just started to be appreciated. With the development of modern technologies in mass spectrometry and deep sequencing, recent examples of modifications regulating host-pathogen interactions have been demonstrated. They show how RNA modifications can regulate immune responses, antibiotic resistance, expression of virulence genes, and bacterial persistence. Here, we illustrate some of these findings, and highlight the strategies used to characterize RNA modifications, and their potential for new therapeutic applications.

Micro-flow hydrophilic interaction liquid chromatography coupled with triple quadrupole mass spectrometry detects modified nucleosides in the transfer RNA pool of cyanobacteria

Post-transcriptional modification of nucleosides is observed in almost all elements of RNA. Modified nucleosides finely tune the structure of RNA molecules and affect vital functions, such as the modified wobble position 34 of transfer RNAs expanding the reading preference of anticodons to codons. Recent investigations have revealed that the modification species and their frequencies in an RNA element are not stable but vary with specific cellular factors including metabolites and particular proteins (writers, readers, and erasers). To understand the link between dynamic RNA modifications and biological processes, sensitive and reliable methods for determining modified nucleosides are required. In this study, micro-flow (8 μL/min) hydrophilic interaction liquid chromatography was coupled with triple quadrupole mass spectrometry for the simultaneous determination of adenosine, uridine, cytidine, guanosine, and 20 modified nucleosides. The method was calibrated using 0.1-1000 nM standards (∼0.03-300 ng/mL) and successfully applied to the determination of transfer RNA modifications in the model cyanobacterium Synechococcus elongatus PCC 7942. A protocol for the isolation of a clean transfer RNA pool was optimized, requiring only 25 ng for the identification and quantification of transfer RNA modifications. This micro-flow liquid chromatography-tandem mass spectrometry method constitutes the first step toward monitoring dynamic ribonucleoside modifications in a limited RNA sample.

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.

The expanding world of tRNA modifications and their disease relevance

Transfer RNA (tRNA) is an adapter molecule that links a specific codon in mRNA with its corresponding amino acid during protein synthesis. tRNAs are enzymatically modified post-transcriptionally. A wide variety of tRNA modifications are found in the tRNA anticodon, which are crucial for precise codon recognition and reading frame maintenance, thereby ensuring accurate and efficient protein synthesis. In addition, tRNA-body regions are also frequently modified and thus stabilized in the cell. Over the past two decades, 16 novel tRNA modifications were discovered in various organisms, and the chemical space of tRNA modification continues to expand. Recent studies have revealed that tRNA modifications can be dynamically altered in response to levels of cellular metabolites and environmental stresses. Importantly, we now understand that deficiencies in tRNA modification can have pathological consequences, which are termed 'RNA modopathies'. Dysregulation of tRNA modification is involved in mitochondrial diseases, neurological disorders and cancer.

Instrumental analysis of RNA modifications

Organisms from all domains of life invest a substantial amount of energy for the introduction of RNA modifications into nearly all transcripts studied to date. Instrumental analysis of RNA can focus on the modified residues and reveal the function of these epitranscriptomic marks. Here, we will review recent advances and breakthroughs achieved by NMR spectroscopy, sequencing, and mass spectrometry of the epitranscriptome.

Distinct evolutionary pathways for the synthesis and function of tRNA modifications

Transfer ribonucleicacids (RNAs) (tRNAs) are essential adaptor molecules for translation. The functions and stability of tRNAs are modulated by their post-transcriptional modifications (tRNA modifications). Each domain of life has a specific set of modifications that include ones shared in multiple domains and ones specific to a domain. In some cases, different tRNA modifications across domains have similar functions to each other. Recent studies uncovered that distinct enzymes synthesize the same modification in different organisms, suggesting that such modifications are acquired through independent evolution. In this short review, I outline the mechanisms by which various modifications contribute to tRNA function, including modulation of decoding and tRNA stability, using recent findings. I also focus on modifications that are synthesized by distinct biosynthetic pathways.

Functions of Bacterial tRNA Modifications: From Ubiquity to Diversity

Modified nucleotides in tRNA are critical components of the translation apparatus, but their importance in the process of translational regulation had until recently been greatly overlooked. Two breakthroughs have recently allowed a fuller understanding of the importance of tRNA modifications in bacterial physiology. One is the identification of the full set of tRNA modification genes in model organisms such as Escherichia coli K12. The second is the improvement of available analytical tools to monitor tRNA modification patterns. The role of tRNA modifications varies greatly with the specific modification within a given tRNA and with the organism studied. The absence of these modifications or reductions can lead to cell death or pleiotropic phenotypes or may have no apparent visible effect. By linking translation through their decoding functions to metabolism through their biosynthetic pathways, tRNA modifications are emerging as important components of the bacterial regulatory toolbox.

Probing the diversity and regulation of tRNA modifications

Transfer RNAs (tRNAs) are non-coding RNAs essential for protein synthesis. tRNAs are heavily decorated with a variety of post-transcriptional modifications (tRNA modifications). Recent methodological advances provide new tools for rapid profiling of tRNA modifications and have led to discoveries of novel modifications and their regulation. Here, we provide an overview of the techniques for investigating tRNA modifications and of the expanding knowledge of their chemistry and regulation.

Survey and Validation of tRNA Modifications and Their Corresponding Genes in Bacillus subtilis sp Subtilis Strain 168

Extensive knowledge of both the nature and position of tRNA modifications in all cellular tRNAs has been limited to two bacteria, Escherichia coli and Mycoplasma capricolum. Bacillus subtilis sp subtilis strain 168 is the model Gram-positive bacteria and the list of the genes involved in tRNA modifications in this organism is far from complete. Mass spectrometry analysis of bulk tRNA extracted from B. subtilis, combined with next generation sequencing technologies and comparative genomic analyses, led to the identification of 41 tRNA modification genes with associated confidence scores. Many differences were found in this model Gram-positive bacteria when compared to E. coli. In general, B. subtilis tRNAs are less modified than those in E. coli, even if some modifications, such as m¹A22 or ms²t⁶A, are only found in the model Gram-positive bacteria. Many examples of non-orthologous displacements and of variations in the most complex pathways are described. Paralog issues make uncertain direct annotation transfer from E. coli to B. subtilis based on homology only without further experimental validation. This difficulty was shown with the identification of the B. subtilis enzyme that introduces ψ at positions 31/32 of the tRNAs. This work presents the most up to date list of tRNA modification genes in B. subtilis, identifies the gaps in knowledge, and lays the foundation for further work to decipher the physiological role of tRNA modifications in this important model organism and other bacteria.