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

Valerie de Crecy-Lagard, Robert J. Ross, Marshall Jaroch, Virginie Marchand, Christina Eisenhart, Damien Brégeon, Yuri Motorin, Patrick Limbach

► To cite this version:

Valerie de Crecy-Lagard, Robert J. Ross, Marshall Jaroch, Virginie Marchand, Christina Eisenhart, et al.. Survey and Validation of tRNA Modifications and Their Corresponding Genes in Bacillus subtilis sp Subtilis Strain 168. Biomolecules, 2020, 10 (7), pp.977. 10.3390/biom10070977. hal-02915956

HAL Id: hal-02915956

https://hal.univ-lorraine.fr/hal-02915956

Submitted on 17 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Survey and Validation of tRNA Modifications and Their Corresponding Genes in Bacillus subtilis sp Subtilis Strain 168

Valerie de Crecy-Lagard, Robert J. Ross, Marshall Jaroch, Virginie Marchand, Christina Eisenhart, Damien Brégeon, Yuri Motorin, Patrick Limbach

► To cite this version:

Valerie de Crecy-Lagard, Robert J. Ross, Marshall Jaroch, Virginie Marchand, Christina Eisenhart, et al.. Survey and Validation of tRNA Modifications and Their Corresponding Genes in Bacillus subtilis sp Subtilis Strain 168. Biomolecules, 2020, 10 (7), pp.977. 10.3390/biom10070977. hal-02915956

HAL Id: hal-02915956

https://hal.univ-lorraine.fr/hal-02915956v1

Submitted on 17 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

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

Valérie de Crécy-Lagard<fnref n="1" />,<fnref n="2" />,* ID, Robert L. Ross<fnref n="3" />, Marshall Jaroch<fnref n="1" />, Virginie Marchand<fnref n="4" /> ID, Christina Eisenhart<fnref n="5" />, Damien Brégeon<fnref n="6" />, Yuri Motorin<fnref n="4" /> ID and Patrick A. Limbach<fnref n="5" /> ID

Received: 4 June 2020; Accepted: 27 June 2020; Published: 30 June 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** 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.

**Keywords:** tRNA modifications; model bacteria; Gram-positive; methylation; pseudouridine synthase; YhcT; YjbO

## 1. Introduction

As adaptors between the mRNA and the elongation peptide, tRNAs are the central decoding elements of the translation process. All tRNA molecules harbor post-transcriptional modifications, and studies from the last 50 years have firmly established the diverse and critical roles that these play in the efficiency and accuracy of translation [1,2]. tRNA modifications can influence tRNA stability and quality control [3,4], they can be determinants or anti-determinants for other components of

<footnote n="1">¹ Department of Microbiology and Cell Science, University of Florida, Gainesville, FL 32611, USA; mjaroch@ufl.edu</footnote>

<footnote n="2">² Genetics Institute, University of Florida, Gainesville, FL 32611, USA</footnote>

<footnote n="3">³ Department of Cancer Biology, University of Cincinnati, Cincinnati, OH 45221, USA; rossrb@ucmail.uc.edu</footnote>

<footnote n="4">⁴ UMR7365 IMoPA CNRS-UL and UMS2008 CNRS-UL-INSERM, Université de Lorraine, Biopôle UL, 54000 Nancy, France; virginie.marchand@univ-lorraine.fr (V.M.); motorine5@univ-lorraine.fr (Y.M.)</footnote>

<footnote n="5">⁵ Department of Chemistry, University of Cincinnati, Cincinnati, OH 45221, USA; psihouca@mail.uc.edu (C.E.); limbacpa@ucmail.uc.edu (P.A.L.)</footnote>

<footnote n="6">⁶ IBPS, Biology of Aging and Adaptation, Sorbonne University, 7 Quai Saint Bernard, CEDEX 05, F-75252 Paris, France; damien.bregeon@sorbonne-universite.fr</footnote>

<footnote>* Correspondence: vcrecy@ufl.edu; Tel.: +1-352-392-9416</footnote>
<figure><img src="image_2.png" /></figure>

Article

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

Valérie de Crécy-Lagard<fnref n="1" />,<fnref n="2" />,* ID, Robert L. Ross<fnref n="3" />, Marshall Jaroch<fnref n="1" />, Virginie Marchand<fnref n="4" /> ID, Christina Eisenhart<fnref n="5" />, Damien Brégeon<fnref n="6" />, Yuri Motorin<fnref n="4" /> ID and Patrick A. Limbach<fnref n="5" /> ID

Received: 4 June 2020; Accepted: 27 June 2020; Published: 30 June 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** 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.

**Keywords:** tRNA modifications; model bacteria; Gram-positive; methylation; pseudouridine synthase; YhcT; YjbO

## 1. Introduction

As adaptors between the mRNA and the elongation peptide, tRNAs are the central decoding elements of the translation process. All tRNA molecules harbor post-transcriptional modifications, and studies from the last 50 years have firmly established the diverse and critical roles that these play in the efficiency and accuracy of translation [1,2]. tRNA modifications can influence tRNA stability and quality control [3,4], they can be determinants or anti-determinants for other components of

<footnote n="1">¹ Department of Microbiology and Cell Science, University of Florida, Gainesville, FL 32611, USA; mjaroch@ufl.edu</footnote>

<footnote n="2">² Genetics Institute, University of Florida, Gainesville, FL 32611, USA</footnote>

<footnote n="3">³ Department of Cancer Biology, University of Cincinnati, Cincinnati, OH 45221, USA; rossrb@ucmail.uc.edu</footnote>

<footnote n="4">⁴ UMR7365 IMoPA CNRS-UL and UMS2008 CNRS-UL-INSERM, Université de Lorraine, Biopôle UL, 54000 Nancy, France; virginie.marchand@univ-lorraine.fr (V.M.); motorine5@univ-lorraine.fr (Y.M.)</footnote>

<footnote n="5">⁵ Department of Chemistry, University of Cincinnati, Cincinnati, OH 45221, USA; psihouca@mail.uc.edu (C.E.); limbacpa@ucmail.uc.edu (P.A.L.)</footnote>

<footnote n="6">⁶ IBPS, Biology of Aging and Adaptation, Sorbonne University, 7 Quai Saint Bernard, CEDEX 05, F-75252 Paris, France; damien.bregeon@sorbonne-universite.fr</footnote>

<footnote>* Correspondence: vcrecy@ufl.edu; Tel.: +1-352-392-9416</footnote>
the translation apparatus such as aminoacyl-tRNA synthetases, and they can restrict or increase the decoding capacity of a given tRNA [5,6]. Since these roles are also modification specific, the effect of inactivating a given tRNA modification enzyme can range from cell death to only subtle phenotypes. For example, t⁶A37 is essential in *E. coli*, most certainly because it is a determinant for IleRS [7]. Lysidine (or k²C) is essential in many bacteria because it allows the modified minor tRNAIleCAU to decode AUA codons [8–10]. On the other hand, deleting most modification genes only gives rise to subtle fitness defects or to phenotypes which are observed in very specific conditions [1,11].

While the role of tRNA modifications as regulatory devices was proposed in the early 1990s by B.C. Persson [12], it is only recently that this model has gained traction [13,14]. Due to the degeneracy of the genetic code, the same amino acid can be encoded by different codons read by tRNAs with different modification profiles. Hence, if specific genes are enriched in codons decoded by a tRNA that harbors a given modification, changes in the level of this specific modification will modify the proteome. The identification of these Modification Tunable Transcripts (or MoTT) in a wide range of organisms (reviewed in [15]) has led to models for regulation in response to stress or metabolic imbalance. Other regulatory roles for tRNA modifications include their capacity to act as determinants for RNases [3] that can generate regulatory half-tRNA fragments, even if this has mostly been observed in eukaryotes [16,17].

To discover these regulatory roles, it is essential to have the foundational knowledge about which tRNA modifications are present in a given organism and which enzymes catalyze their synthesis. Unfortunately, to date, this knowledge is limited to only two bacteria. One is the Gram-positive *Mycoplasma capricolum*. Indeed, because of its small genome size and the pioneer work that generated all tRNA sequences, it was the first bacterium with a full set of predicted tRNA modification genes [18]. The other is the model Gram-negative *Escherichia coli* K12. As shown in Figure 1, decades of biochemical and genetic analyses have identified all the genes involved in the discovered tRNA modifications in this organism. These 59 tRNA modification enzymes and factors represent 1.5% of the *E. coli* encoded proteome.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** *Escherichia coli* tRNA modifications and corresponding enzymes. Proteins encoded by essential genes are in blue. Proteins with a * also modify rRNAs. The annotations were captured from the Björk and Hargevall review [11], from Uniprot, and the recent literature. In green are proteins with no ortholog in *B. subtilis*. Uniprot IDs for all proteins listed are given Table S1.</figcaption></figure>
the translation apparatus such as aminoacyl-tRNA synthetases, and they can restrict or increase the decoding capacity of a given tRNA [5,6]. Since these roles are also modification specific, the effect of inactivating a given tRNA modification enzyme can range from cell death to only subtle phenotypes. For example, t⁶A37 is essential in *E. coli*, most certainly because it is a determinant for IleRS [7]. Lysidine (or k²C) is essential in many bacteria because it allows the modified minor tRNAIleCAU to decode AUA codons [8–10]. On the other hand, deleting most modification genes only gives rise to subtle fitness defects or to phenotypes which are observed in very specific conditions [1,11].

While the role of tRNA modifications as regulatory devices was proposed in the early 1990s by B.C. Persson [12], it is only recently that this model has gained traction [13,14]. Due to the degeneracy of the genetic code, the same amino acid can be encoded by different codons read by tRNAs with different modification profiles. Hence, if specific genes are enriched in codons decoded by a tRNA that harbors a given modification, changes in the level of this specific modification will modify the proteome. The identification of these Modification Tunable Transcripts (or MoTT) in a wide range of organisms (reviewed in [15]) has led to models for regulation in response to stress or metabolic imbalance. Other regulatory roles for tRNA modifications include their capacity to act as determinants for RNases [3] that can generate regulatory half-tRNA fragments, even if this has mostly been observed in eukaryotes [16,17].

To discover these regulatory roles, it is essential to have the foundational knowledge about which tRNA modifications are present in a given organism and which enzymes catalyze their synthesis. Unfortunately, to date, this knowledge is limited to only two bacteria. One is the Gram-positive *Mycoplasma capricolum*. Indeed, because of its small genome size and the pioneer work that generated all tRNA sequences, it was the first bacterium with a full set of predicted tRNA modification genes [18]. The other is the model Gram-negative *Escherichia coli* K12. As shown in Figure 1, decades of biochemical and genetic analyses have identified all the genes involved in the discovered tRNA modifications in this organism. These 59 tRNA modification enzymes and factors represent 1.5% of the *E. coli* encoded proteome.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** *Escherichia coli* tRNA modifications and corresponding enzymes. Proteins encoded by essential genes are in blue. Proteins with a * also modify rRNAs. The annotations were captured from the Björk and Hargevall review [11], from Uniprot, and the recent literature. In green are proteins with no ortholog in *B. subtilis*. Uniprot IDs for all proteins listed are given Table S1.</figcaption></figure>
Recently, a combination of mass spectrometry (MS) and tRNA-Seq methods to detect tRNA modifications was performed on *Vibrio cholerae* [19], and even though this bacteria is closely related to *E. coli*, novel modifications were identified. Full analysis of the modified nucleosides present in tRNAs has been performed for three Gram-positive organisms, *Lactococcus lactis* [20], *Streptomyces griseus* [21], and *Mycobacterium bovis* BCG [22], but the links with the corresponding genes were not made in these studies. *Bacillus subtilis* sp. subtilis strain 168 is the model Gram-positive that was sequenced in 1997 [23], and an active community is keeping the genome annotation up to date [24,25]. Even the current knowledge of tRNA modification genes has not been well captured in SubtiWiki, as not all tRNAs have been sequenced. There are some clear gaps in knowledge, which we attempted to fill in the current study by analyzing the modifications found in *B. subtilis* tRNAs, using a combination of MS and tRNA-Seq methods, linking them to the corresponding validated or predicted gene, and experimentally testing a paralog case for the synthesis of the pseudouridines at positions 31 and 32.

## 2. Materials and Methods

### 2.1. Bioinformatic Analyses

The *B. subtilis* genome that was used was NC_000964. The presence and numbers of isoforms were extracted from the GtRNA database (http://gtrnadb.ucsc.edu/GtRNAdb2/genomes/bacteria/Baci_subt_subtilis_168/) [26]. The RNA sequence data were extracted from Modomics (http://modomics.genesilico.pl/sequences/list/?type_field=tRNA&subtype=all&species=4&display=Display&nomenclature=abbrev) [27] from tRNAb (http://trnadb.bioinf.uni-leipzig.de/) [28] and TpsiC databases (http://tpsic.igcz.poznan.pl/info/start/) [29]. The information on all genes involved in tRNA modification in *E. coli* and *B. subtilis* (Tables S1 and S2) were gathered from the literature and further checked using Subtiwiki [24], Uniprot Id mapping tools [30], and BlastP [31] searches at NCBI [32], as well as Pathosystems Resource Integration Center (PATRIC) feature group tools [33]. The SEED database [34] was used to analyze physical clustering. Sequences were visualized using the integrated genome viewer (IGV) (http://software.broadinstitute.org/software/igv/) [35]. Clustergrammer (http://amp.pharm.mssm.edu/clustergrammer/) was used to visualize fitness and expression data as heatmaps [36]. The Colombos v3.0 database (http://colombos.net/) [37] was used to mine the available *B. subtilis* expression data. Using the list of predicted *B. subtilis* tRNA modification genes from Table S2 as input, an expression matrix comprised of 90 experiments (or contrasts) where these genes showed significant expression was generated from the data available in the database, and further analysis was conducted on this dataset. Each contrast is a single RNA-Seq or microarray experiment where RNA transcript levels for a test condition are compared to a reference culture. For each contrast, Colombos assigns a score that corresponds to how many genes show similar changes in expression. If the changes in expression, relative to their respective reference, increase or decrease uniformly for all the input genes, then a higher score is awarded to the contrast. In the case where a contrast's tested condition causes the input genes to show dissimilar and noncoherent expression, a lower score is awarded to the contrast. From the list of 90 contrasts, 47 contrasts were selected, with scores less than 0.5, in order to analyze the conditions that caused the input genes to show dissimilar expression. To illustrate drastic changes in expression, 12 genes were selected that showed unique changes in expression values from 4.5 (red) to -5.0 (blue) across the 47 contrasts.

### 2.2. Strains and Media and Bulk tRNA Preparation

*Bacillus subtilis* subsp. subtilis str. 168 (RefSeq: NC_000964.3) was used as the wild-type strain. The mutants used were BKE11620 (*trpC2 ΔyjbO::kan*) and BKE09210 (*trpC2 ΔyhcT::kan*) [38]. Cells were grown in Luria-Bertani medium (LB), supplemented with kanamycin (5µg.mL⁻¹) at 30 °C, and harvested at an OD (A₆₀₀nm) of 0.6. Bulk RNA was extracted by hot acidic phenol [35] or using chaotropic agents [39].
Recently, a combination of mass spectrometry (MS) and tRNA-Seq methods to detect tRNA modifications was performed on *Vibrio cholerae* [19], and even though this bacteria is closely related to *E. coli*, novel modifications were identified. Full analysis of the modified nucleosides present in tRNAs has been performed for three Gram-positive organisms, *Lactococcus lactis* [20], *Streptomyces griseus* [21], and *Mycobacterium bovis* BCG [22], but the links with the corresponding genes were not made in these studies. *Bacillus subtilis* sp. subtilis strain 168 is the model Gram-positive that was sequenced in 1997 [23], and an active community is keeping the genome annotation up to date [24,25]. Even the current knowledge of tRNA modification genes has not been well captured in SubtiWiki, as not all tRNAs have been sequenced. There are some clear gaps in knowledge, which we attempted to fill in the current study by analyzing the modifications found in *B. subtilis* tRNAs, using a combination of MS and tRNA-Seq methods, linking them to the corresponding validated or predicted gene, and experimentally testing a paralog case for the synthesis of the pseudouridines at positions 31 and 32.

## 2. Materials and Methods

### 2.1. Bioinformatic Analyses

The *B. subtilis* genome that was used was NC_000964. The presence and numbers of isoforms were extracted from the GtRNA database (http://gtrnadb.ucsc.edu/GtRNAdb2/genomes/bacteria/Baci_subt_subtilis_168/) [26]. The RNA sequence data were extracted from Modomics (http://modomics.genesilico.pl/sequences/list/?type_field=tRNA&subtype=all&species=4&display=Display&nomenclature=abbrev) [27] from tRNAb (http://trnadb.bioinf.uni-leipzig.de/) [28] and TpsiC databases (http://tpsic.igcz.poznan.pl/info/start/) [29]. The information on all genes involved in tRNA modification in *E. coli* and *B. subtilis* (Tables S1 and S2) were gathered from the literature and further checked using Subtiwiki [24], Uniprot Id mapping tools [30], and BlastP [31] searches at NCBI [32], as well as Pathosystems Resource Integration Center (PATRIC) feature group tools [33]. The SEED database [34] was used to analyze physical clustering. Sequences were visualized using the integrated genome viewer (IGV) (http://software.broadinstitute.org/software/igv/) [35]. Clustergrammer (http://amp.pharm.mssm.edu/clustergrammer/) was used to visualize fitness and expression data as heatmaps [36]. The Colombos v3.0 database (http://colombos.net/) [37] was used to mine the available *B. subtilis* expression data. Using the list of predicted *B. subtilis* tRNA modification genes from Table S2 as input, an expression matrix comprised of 90 experiments (or contrasts) where these genes showed significant expression was generated from the data available in the database, and further analysis was conducted on this dataset. Each contrast is a single RNA-Seq or microarray experiment where RNA transcript levels for a test condition are compared to a reference culture. For each contrast, Colombos assigns a score that corresponds to how many genes show similar changes in expression. If the changes in expression, relative to their respective reference, increase or decrease uniformly for all the input genes, then a higher score is awarded to the contrast. In the case where a contrast's tested condition causes the input genes to show dissimilar and noncoherent expression, a lower score is awarded to the contrast. From the list of 90 contrasts, 47 contrasts were selected, with scores less than 0.5, in order to analyze the conditions that caused the input genes to show dissimilar expression. To illustrate drastic changes in expression, 12 genes were selected that showed unique changes in expression values from 4.5 (red) to -5.0 (blue) across the 47 contrasts.

### 2.2. Strains and Media and Bulk tRNA Preparation

*Bacillus subtilis* subsp. subtilis str. 168 (RefSeq: NC_000964.3) was used as the wild-type strain. The mutants used were BKE11620 (*trpC2 ΔyjbO::kan*) and BKE09210 (*trpC2 ΔyhcT::kan*) [38]. Cells were grown in Luria-Bertani medium (LB), supplemented with kanamycin (5µg.mL⁻¹) at 30 °C, and harvested at an OD (A₆₀₀nm) of 0.6. Bulk RNA was extracted by hot acidic phenol [35] or using chaotropic agents [39].
## 2.3. *Mass Spectrometry Analyses*

### 2.3.1. LC-MS of Nucleosides

Digestion of the purified tRNA fraction was performed as previously described [39]. Separation was accomplished by reversed phase chromatography using an Acquity UPLC HSS T3, 1.8 µm, 1 mm × 100 mm column, (Waters, Milford, MA, USA) on a Vanquish Flex Quaternary UHPLC system (Thermo Fisher Scientific, San Jose, CA, USA). Mobile phase A consisted of 5.3 mM ammonium acetate (Thermo Fisher Scientific) in LC-MS grade water (Alfa Aesar, Haverhill, MA, USA), pH 5.3. Mobile phase B consisted of a 60:40 mixture of 5.3 mM ammonium acetate and acetonitrile (Honeywell Burdick & Jackson, Morris Plains, NJ, USA) with a gradient of 0% B (from 0 to 1.8 min), 2% B at 3 to 3.5 min, 3% B at 4.1 min, 5% B at 7 min, 25% B at 9 min, 35% B at 15 min, 99% B at 15.5 min (hold for 4.5 min), and 99% B at 20 min, then returning to 0% B at 25.5 min, at a flow rate of 100 µL min⁻¹. The column temperature was set to 40 °C. Accurate mass values and physio-chemical properties were determined using Marvin 17.3.13.0 (https://chemaxon.com).

High resolution accurate mass (HRAM) analysis was performed on an Orbitrap Fusion Lumos Tribrid mass spectrometer (Thermo Fisher Scientific) interfaced with an H-ESI electrospray source in positive polarity mode. Full scan data were acquired at a resolution of 120,000, mass range 220–900 *m/z*, AGC 7.5e4, and IT 100 ms. Data-dependent top speed MS/MS spectra (1 s cycle, CID 40%) were acquired in the ion trap with rapid scan rate, AGC 1.0e5, and an injection time of 300 ms, with scan range mode set to auto *m/z* normal. The other instrumental conditions were as follows: quadrupole isolation of 1 *m/z*; RF 35%; sheath gas, auxiliary gas, and sweep gas of 30, 10, and 0 arbitrary units, respectively; ion transfer tube temperature of 289 °C; vaporizer temperature of 92 °C; and spray voltage of 3500 V. Data were analyzed using Xcalibur 4.0, Compound Discoverer 3.0, and mzVault 2.1 (Thermo Fisher Scientific). Fragment Ion Search (FISh) settings were set at a high mass accuracy of 2.5 ppm, a low mass accuracy of 0.5 Da, and a S/N of 5000.

### 2.3.2. LC-MS Analysis of Oligonucleotides

Separation was accomplished by HILIC chromatography, using a Shodex HILICPak VN-50 2D, 5 µm, 2 mm × 150 mm column, (Showa Denko, New York, NY, USA) on a Vanquish UHPLC system (Thermo Fisher Scientific, San Jose, CA, USA). Mobile phase A consisted of 50 mM ammonium bicarbonate (Sigma Aldrich, St. Louis, MO, USA) in 18 Ohm millipure water, pH 9.0. Mobile phase B consisted of acetonitrile (Honeywell Burdick & Jackson, Morris Plains, NJ, USA) 0.1% NH₄OH, with a gradient of 30% B for 1.5 min, increasing to 35% B at 17 min, 70% B at 30 min, and holding for 2 min, then returning to 30% B at 32.1 min and re-equilibrating to 45 min at a flow rate of 200 µL min⁻¹. The column temperature was set at 20 °C.

RNase T1 digestion products were analyzed using the Orbitrap Fusion Lumos Tribrid mass spectrometer (Thermo Fisher Scientific) interfaced with an H-ESI electrospray source in negative polarity mode. Data-dependent top speed MS/MS spectra (1 s cycle, CID 40%) were acquired in the ion trap with rapid scan rate, AGC 1.0e5, and an injection time of 300 ms. Data-dependent top speed MS/MS spectra (1 s cycle, CID 40%) were acquired in the ion trap with rapid scan rate, AGC 1.0e5, and an injection time of 300 ms. The other instrumental conditions were as follows: quadrupole isolation of 1.6 *m/z*; RF 30%; sheath gas, auxiliary gas, and sweep gas of 40, 8, and 1 arbitrary units, respectively; ion transfer tube temperature of 275 °C; vaporizer temperature of 320 °C; and spray voltage of 2500 V. Data were analyzed using Xcalibur 4.0 and RAMM 1.0 [40].

## 2.4. *Identification of Modification by Next Generation Sequencing Methods*

Deep sequencing-based identification of the *B. subtilis* tRNA modifications was performed using RibOMethSeq protocol, allowing us to map 2'-O-methylations by their protection to cleavage [41,42]. Reverse Transcriptase (RT) misincorporation signatures at RT-arresting nucleotides [43,44] were extracted from RibOMethSeq data and also from the RT-primer extension, using TTO-based ScriptSeq v2
## 2.3. *Mass Spectrometry Analyses*

### 2.3.1. LC-MS of Nucleosides

Digestion of the purified tRNA fraction was performed as previously described [39]. Separation was accomplished by reversed phase chromatography using an Acquity UPLC HSS T3, 1.8 µm, 1 mm × 100 mm column, (Waters, Milford, MA, USA) on a Vanquish Flex Quaternary UHPLC system (Thermo Fisher Scientific, San Jose, CA, USA). Mobile phase A consisted of 5.3 mM ammonium acetate (Thermo Fisher Scientific) in LC-MS grade water (Alfa Aesar, Haverhill, MA, USA), pH 5.3. Mobile phase B consisted of a 60:40 mixture of 5.3 mM ammonium acetate and acetonitrile (Honeywell Burdick & Jackson, Morris Plains, NJ, USA) with a gradient of 0% B (from 0 to 1.8 min), 2% B at 3 to 3.5 min, 3% B at 4.1 min, 5% B at 7 min, 25% B at 9 min, 35% B at 15 min, 99% B at 15.5 min (hold for 4.5 min), and 99% B at 20 min, then returning to 0% B at 25.5 min, at a flow rate of 100 µL min⁻¹. The column temperature was set to 40 °C. Accurate mass values and physio-chemical properties were determined using Marvin 17.3.13.0 (https://chemaxon.com).

High resolution accurate mass (HRAM) analysis was performed on an Orbitrap Fusion Lumos Tribrid mass spectrometer (Thermo Fisher Scientific) interfaced with an H-ESI electrospray source in positive polarity mode. Full scan data were acquired at a resolution of 120,000, mass range 220–900 m/z, AGC 7.5e4, and IT 100 ms. Data-dependent top speed MS/MS spectra (1 s cycle, CID 40%) were acquired in the ion trap with rapid scan rate, AGC 1.0e5, and an injection time of 300 ms, with scan range mode set to auto m/z normal. The other instrumental conditions were as follows: quadrupole isolation of 1 m/z; RF 35%; sheath gas, auxiliary gas, and sweep gas of 30, 10, and 0 arbitrary units, respectively; ion transfer tube temperature of 289 °C; vaporizer temperature of 92 °C; and spray voltage of 3500 V. Data were analyzed using Xcalibur 4.0, Compound Discoverer 3.0, and mzVault 2.1 (Thermo Fisher Scientific). Fragment Ion Search (FISh) settings were set at a high mass accuracy of 2.5 ppm, a low mass accuracy of 0.5 Da, and a S/N of 5000.

### 2.3.2. LC-MS Analysis of Oligonucleotides

Separation was accomplished by HILIC chromatography, using a Shodex HILICPak VN-50 2D, 5 µm, 2 mm × 150 mm column, (Showa Denko, New York, NY, USA) on a Vanquish UHPLC system (Thermo Fisher Scientific, San Jose, CA, USA). Mobile phase A consisted of 50 mM ammonium bicarbonate (Sigma Aldrich, St. Louis, MO, USA) in 18 Ohm millipure water, pH 9.0. Mobile phase B consisted of acetonitrile (Honeywell Burdick & Jackson, Morris Plains, NJ, USA) 0.1% NH₄OH, with a gradient of 30% B for 1.5 min, increasing to 35% B at 17 min, 70% B at 30 min, and holding for 2 min, then returning to 30% B at 32.1 min and re-equilibrating to 45 min at a flow rate of 200 µL min⁻¹. The column temperature was set at 20 °C.

RNase T1 digestion products were analyzed using the Orbitrap Fusion Lumos Tribrid mass spectrometer (Thermo Fisher Scientific) interfaced with an H-ESI electrospray source in negative polarity mode. Data-dependent top speed MS/MS spectra (1 s cycle, CID 40%) were acquired in the ion trap with rapid scan rate, AGC 1.0e5, and an injection time of 300 ms. Data-dependent top speed MS/MS spectra (1 s cycle, CID 40%) were acquired in the ion trap with rapid scan rate, AGC 1.0e5, and an injection time of 300 ms. The other instrumental conditions were as follows: quadrupole isolation of 1.6 m/z; RF 30%; sheath gas, auxiliary gas, and sweep gas of 40, 8, and 1 arbitrary units, respectively; ion transfer tube temperature of 275 °C; vaporizer temperature of 320 °C; and spray voltage of 2500 V. Data were analyzed using Xcalibur 4.0 and RAMM 1.0 [40].

## 2.4. *Identification of Modification by Next Generation Sequencing Methods*

Deep sequencing-based identification of the *B. subtilis* tRNA modifications was performed using RibOMethSeq protocol, allowing us to map 2'-O-methylations by their protection to cleavage [41,42]. Reverse Transcriptase (RT) misincorporation signatures at RT-arresting nucleotides [43,44] were extracted from RibOMethSeq data and also from the RT-primer extension, using TTO-based ScriptSeq v2
(Illumina, San Diego, CA, USA) library preparation kit. RiboMethSeq analysis of tRNAs was performed using alkaline fragmentation of total *B. subtilis* RNA, followed by library preparation and sequencing [42]. RT misincorporation signatures were extracted from Samtools mpileup format and manually validated by inspection of aligned reads in \*.bam file in IGV. Pseudouridine mapping (Psiseq) was performed in the tRNA-enriched fraction of *B. subtilis* RNAs, using a chemical-based protocol derived from CMCT (1-cyclohexyl-(2-morpholinoethyl)carbodiimide metho-p-toluene sulfonate)-chemical mapping protocol, and specifically adapted to tRNAs [45–47].

# 3. Results and Discussion

## 3.1. Identity and Position of tRNA Modifications in the Model Gram-Positive *B. subtilis* 168

A total of 86 tRNA genes specifying 35 different iso-acceptors can be identified in the *Bacillus subtilis* 168 genome (Figure 2). The RNA sequences of 24 of these have been reported (black dots in Figure 2), but 11 tRNA isoacceptors remain to be sequenced (red dots in Figure 2). To fill this knowledge gap, different methods were combined. First, a comprehensive MS analysis of the bulk tRNA extracted from *B. subtilis* cells was performed. These data were complemented by analysis of the RT-arresting nucleotides and 2'-O-methylation by RiboMethSeq, as well as pseudouridines by PsiSeq.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Compilation of sequenced *B. subtilis* tRNA genes and tRNA molecules. Circular representation of the genetic code displaying modified nucleosides based on Grosjean and Westhof [48]. The codon sequence of each amino acid is read in the inside-out direction (1-2-3). The numbers of isoforms are shown by the number in parentheses and the data were extracted from the GtRNA-db 2.0 [26]. A red dot specifies that the tRNA has not been sequenced, while a black dot specifies that the tRNA sequence is available based on the information available in Modomics [27]. The “e” and “i” specify the elongator or initiator tRNA<sup>Met</sup><sub>CAU</sub>, respectively.</figcaption></figure>

Two types of LC-MS/MS analyses were performed. The first consisted of a total nucleoside identification assay, where an aliquot of the bulk tRNA is taken to its nucleoside monomers by enzymatic digestion to identify which modifications are present in the sample within the detection limits. In the second assay, an aliquot of the bulk tRNA is digested with base specific nucleases to
(Illumina, San Diego, CA, USA) library preparation kit. RiboMethSeq analysis of tRNAs was performed using alkaline fragmentation of total *B. subtilis* RNA, followed by library preparation and sequencing [42]. RT misincorporation signatures were extracted from Samtools mpileup format and manually validated by inspection of aligned reads in \*.bam file in IGV. Pseudouridine mapping (Psiseq) was performed in the tRNA-enriched fraction of *B. subtilis* RNAs, using a chemical-based protocol derived from CMCT (1-cyclohexyl-(2-morpholinoethyl)carbodiimide metho-p-toluene sulfonate)-chemical mapping protocol, and specifically adapted to tRNAs [45–47].

# 3. Results and Discussion

## 3.1. Identity and Position of tRNA Modifications in the Model Gram-Positive *B. subtilis* 168

A total of 86 tRNA genes specifying 35 different iso-acceptors can be identified in the *Bacillus subtilis* 168 genome (Figure 2). The RNA sequences of 24 of these have been reported (black dots in Figure 2), but 11 tRNA isoacceptors remain to be sequenced (red dots in Figure 2). To fill this knowledge gap, different methods were combined. First, a comprehensive MS analysis of the bulk tRNA extracted from *B. subtilis* cells was performed. These data were complemented by analysis of the RT-arresting nucleotides and 2'-O-methylation by RiboMethSeq, as well as pseudouridines by PsiSeq.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Compilation of sequenced *B. subtilis* tRNA genes and tRNA molecules. Circular representation of the genetic code displaying modified nucleosides based on Grosjean and Westhof [48]. The codon sequence of each amino acid is read in the inside-out direction (1-2-3). The numbers of isoforms are shown by the number in parentheses and the data were extracted from the GtRNA-db 2.0 [26]. A red dot specifies that the tRNA has not been sequenced, while a black dot specifies that the tRNA sequence is available based on the information available in Modomics [27]. The “e” and “i” specify the elongator or initiator tRNA<sup>Met</sup><sub>CAU</sub>, respectively.</figcaption></figure>

Two types of LC-MS/MS analyses were performed. The first consisted of a total nucleoside identification assay, where an aliquot of the bulk tRNA is taken to its nucleoside monomers by enzymatic digestion to identify which modifications are present in the sample within the detection limits. In the second assay, an aliquot of the bulk tRNA is digested with base specific nucleases to
create oligonucleotides that can be separated chromatographically by length and whose fragmentation ladder allows the placement of the modifications detected in the nucleoside assay into their proper sequence context. For modification mapping by LC-MS/MS, an initial survey scan was conducted using lower resolution mass analyzers (i.e., ion trap and q-TOF), and these data provided a rough map to guide higher resolution analyses conducted on the orbitrap. The mass spectrometry data reported here resulted from final verification using the orbitrap, with a few exceptions that are noted.

Next Generation Sequencing (NGS)-based analysis of the tRNA modifications (and, eventually, rRNA modifications) was performed using three complementary techniques. First, analysis of the RT-arrests and misincorporation profiles was performed using the TTO-based Illumina library preparation kit. However, random priming during the RT-primer extension step makes this kit more suitable for long RNAs (like rRNAs) and not for short tRNAs, thus these data were only used as a complement. The second and major approach used in the study was standard RiboMethSeq, initially developed for 2'-O-methylation detection. We also noticed that, in addition to 2'-O-methylations, moderate RiboMethSeq signals were frequently detected also for pseudouridines, dihydrouridines (D), and m<sup>7</sup>G residues, present in tRNAs. These data were used only as supplementary confirmation of known sites. In addition, the same RiboMethSeq datasets were used for the detection of RT-misincorporations. Direct ligation of 3'- and 5'-adapters to RNA or RNA fragments makes this approach unsuitable for the detection of hard RT-stops, but strong RT misincorporation profiles (e.g., Inosine, m<sup>1</sup>A, m<sup>3</sup>U, m<sup>6,2</sup>A, m<sup>2,2</sup>G, s<sup>4</sup>U, etc.) were still detected, albeit with a somewhat lower efficiency due to the robustness of SuperScript III, which was used in the RT step. The third deep sequencing-based approach used a modified pseudouridine-specific chemical modification protocol and allowed the detection of candidate sites in the tRNAs. A summary of the RiboMethSeq and PsiSeq analyses of the *B. subtilis* tRNAs is given in Table S3.

Thirty-five modified ribonucleosides were identified in the bulk purified tRNA by LC-MS/MS (Table 1). Of these 35, 21 had previously been mapped in *B. subtilis*, as listed in Modomics or in additional studies [49,50]. Several modifications (epoxyQ (oQ), t<sup>6</sup>A, mnm<sup>5</sup>U, and nm<sup>5</sup>U) are intermediates in complex modification pathways [51–53]. Another two (2'-O-methylations Cm and U*m*) were only detected by RiboMethSeq. Several modified ribonucleosides must result from rRNA contamination copurifying with the tRNA during anion exchange fractionation [39]. These include m<sup>4</sup>Cm, m<sup>6,2</sup>A, m<sup>2</sup>G, m<sup>2,2</sup>G, and Am. m<sup>4</sup>Cm is not found in any tRNA but is present at position 1402 of *E. coli* 16S rRNA [54]. This modification is conserved in other Gram-positive bacteria [55], with *B. subtilis* encoding the orthologs of the *E. coli* RmsH/RsmI enzymes involved in m<sup>4</sup>Cm1402 synthesis (BSU_00360 and BSU_00360, respectively). The contaminating rRNAs were abundant enough that they could be analyzed by NGS, and a candidate site for this rRNA modification (C*m*) in *B. subtilis* rRNA was detected by RiboMethSeq, in a highly conserved sequence at position 1410 (equivalent to *E. coli* position 1402) (Figure S1). m<sup>6,2</sup>A is another 16S rRNA modification, found at two adjacent adenosines (A1518 and A1519, *E. coli* numbering), that is very conserved in rRNA but not found in tRNAs, and it is introduced by the KsgA enzyme [56,57] that is present in *B. subtilis* (BSU_00420). Both m<sup>6,2</sup>A residues are visible by their RT-signatures at the corresponding positions A1527 and A1528 in *B. subtilis* 16S rRNA, as shown Figure S1. Two other modifications detected in the MS analysis of bulk nucleosides were also attributed to rRNA contamination. The methylated guanosines m<sup>2</sup>G and m<sup>2,2</sup>G are found in archaeal and eukaryotic tRNAs, where they have structural roles [58–60], but they have never been identified in bacterial tRNAs. Both are found in bacterial rRNAs [55,61]. An ortholog of *E. coli* RsmD involved in m<sup>2</sup>G966 formation in 16S rRNA is present in *B. subtilis* (YlbH/BSU_15010). In *Clostridium acetobutylicum* and *Thermus thermophilus*, this position is actually modified to m<sup>2,2</sup>G [55], so this could also be the case in *B. subtilis*. In addition, an ortholog of *E. coli* RlmKL that is involved in m<sup>2</sup>G2445 in the formation 23S rRNA was identified in *B. subtilis* (YpsC/BSU_22170). We were able to tentatively map the m<sup>2,2</sup>G modification to CpG dimers (Cp[m<sup>2,2</sup>G]p) at the following positions: G285 in 16S rRNA and G2884 in 23S rRNA (Figure S1). Finally, the Am modification could also come from an rRNA
create oligonucleotides that can be separated chromatographically by length and whose fragmentation ladder allows the placement of the modifications detected in the nucleoside assay into their proper sequence context. For modification mapping by LC-MS/MS, an initial survey scan was conducted using lower resolution mass analyzers (i.e., ion trap and q-TOF), and these data provided a rough map to guide higher resolution analyses conducted on the orbitrap. The mass spectrometry data reported here resulted from final verification using the orbitrap, with a few exceptions that are noted.

Next Generation Sequencing (NGS)-based analysis of the tRNA modifications (and, eventually, rRNA modifications) was performed using three complementary techniques. First, analysis of the RT-arrests and misincorporation profiles was performed using the TTO-based Illumina library preparation kit. However, random priming during the RT-primer extension step makes this kit more suitable for long RNAs (like rRNAs) and not for short tRNAs, thus these data were only used as a complement. The second and major approach used in the study was standard RiboMethSeq, initially developed for 2'-O-methylation detection. We also noticed that, in addition to 2'-O-methylations, moderate RiboMethSeq signals were frequently detected also for pseudouridines, dihydrouridines (D), and m<sup>7</sup>G residues, present in tRNAs. These data were used only as supplementary confirmation of known sites. In addition, the same RiboMethSeq datasets were used for the detection of RT-misincorporations. Direct ligation of 3'- and 5'-adapters to RNA or RNA fragments makes this approach unsuitable for the detection of hard RT-stops, but strong RT misincorporation profiles (e.g., Inosine, m<sup>1</sup>A, m<sup>3</sup>U, m<sup>6,2</sup>A, m<sup>2,2</sup>G, s<sup>4</sup>U, etc.) were still detected, albeit with a somewhat lower efficiency due to the robustness of SuperScript III, which was used in the RT step. The third deep sequencing-based approach used a modified pseudouridine-specific chemical modification protocol and allowed the detection of candidate sites in the tRNAs. A summary of the RiboMethSeq and PsiSeq analyses of the *B. subtilis* tRNAs is given in Table S3.

Thirty-five modified ribonucleosides were identified in the bulk purified tRNA by LC-MS/MS (Table 1). Of these 35, 21 had previously been mapped in *B. subtilis*, as listed in Modomics or in additional studies [49,50]. Several modifications (epoxyQ (oQ), t<sup>6</sup>A, mnm<sup>5</sup>U, and nm<sup>5</sup>U) are intermediates in complex modification pathways [51–53]. Another two (2'-O-methylations Cm and U*m*) were only detected by RiboMethSeq. Several modified ribonucleosides must result from rRNA contamination copurifying with the tRNA during anion exchange fractionation [39]. These include m<sup>4</sup>Cm, m<sup>6,2</sup>A, m<sup>2</sup>G, m<sup>2,2</sup>G, and Am. m<sup>4</sup>Cm is not found in any tRNA but is present at position 1402 of *E. coli* 16S rRNA [54]. This modification is conserved in other Gram-positive bacteria [55], with *B. subtilis* encoding the orthologs of the *E. coli* RmsH/RsmI enzymes involved in m<sup>4</sup>Cm1402 synthesis (BSU_00360 and BSU_00360, respectively). The contaminating rRNAs were abundant enough that they could be analyzed by NGS, and a candidate site for this rRNA modification (C*m*) in *B. subtilis* rRNA was detected by RiboMethSeq, in a highly conserved sequence at position 1410 (equivalent to *E. coli* position 1402) (Figure S1). m<sup>6,2</sup>A is another 16S rRNA modification, found at two adjacent adenosines (A1518 and A1519, *E. coli* numbering), that is very conserved in rRNA but not found in tRNAs, and it is introduced by the KsgA enzyme [56,57] that is present in *B. subtilis* (BSU_00420). Both m<sup>6,2</sup>A residues are visible by their RT-signatures at the corresponding positions A1527 and A1528 in *B. subtilis* 16S rRNA, as shown Figure S1. Two other modifications detected in the MS analysis of bulk nucleosides were also attributed to rRNA contamination. The methylated guanosines m<sup>2</sup>G and m<sup>2,2</sup>G are found in archaeal and eukaryotic tRNAs, where they have structural roles [58–60], but they have never been identified in bacterial tRNAs. Both are found in bacterial rRNAs [55,61]. An ortholog of *E. coli* RsmD involved in m<sup>2</sup>G966 formation in 16S rRNA is present in *B. subtilis* (YlbH/BSU_15010). In *Clostridium acetobutylicum* and *Thermus thermophilus*, this position is actually modified to m<sup>2,2</sup>G [55], so this could also be the case in *B. subtilis*. In addition, an ortholog of *E. coli* RlmKL that is involved in m<sup>2</sup>G2445 in the formation 23S rRNA was identified in *B. subtilis* (YpsC/BSU_22170). We were able to tentatively map the m<sup>2,2</sup>G modification to CpG dimers (Cp[m<sup>2,2</sup>G]p) at the following positions: G285 in 16S rRNA and G2884 in 23S rRNA (Figure S1). Finally, the Am modification could also come from an rRNA
modification, since a number of candidate Am sites were found (positions 224 and 1455 in 16S rRNA as well as 1075, 1524, 2220, and 2827 in 23S rRNA) (Figure S1).

Several observed modifications must be artifacts of the sample preparation. m¹I has only been observed in archaeal tRNA [62] but could arise from the chemical or enzymatic deamination of m¹A during the enzymatic digestion to nucleosides [63]. Similarly, oxidation inside the cell or during analysis is most certainly the source of 8-oxoG formation [64].

The base position cannot be determined for simple methylations by oligonucleotide mass mapping; for example, m¹A, m²A, and m⁶A are all detected as only a methylated adenosine. We were able to map an adenosine methylation at position 37 on tRNA MetCAU, as well as ac⁴C at position 34 (Figure S2A). The sequence of the *E. coli* equivalent tRNA suggests that the probable modification is m⁶A. A methylated adenosine was also mapped to position 22 of the tRNA Leu UAA (Figure S2B). This adenosine was identified as m¹A, based on prior literature [65], and this was confirmed by RT-signature analysis (Figure S1 and Table S3). We also mapped the m²A found at position 37 of tRNA Glu UUC, having mnm⁵s²U at U34 (Figure S2C). The presence of m¹A previously reported at position 58 of tRNA Leu CAG [66] was never reproduced (Louis Droogmans, personal communication) and was not detected in our analyses. The unknown C32 modification in tRNA LysTTT reported by Vold et al. [67] was not detected by either the LC-MS or the RiboMethSeq approach.

In addition, we mapped mnm⁵s²U (Figure S2D), as well as mnm⁵U (Figure S2E) and cmnm⁵s²U (Figure S2F), to different tRNAs (Table 1). The modification 5-methoxyluridine (mo⁵U) is found in a number of *E. coli* tRNAs at position 34. Here, we were able to map mo⁵U position 34 of tRNA Leu UAG (Figure S2G). Moreover, 4-thiouridine was mapped to position 8 of tRNA GlnUUG by MS (Figure S2H) and of tRNA Leu UAA by RT (Figure S1). The hypermodification ct⁶A, along with the t⁶A counterpart, was found in tRNA Thr UGU (Figure S2I,J) and tRNA Ile GAU (Figure S2K,L) at position 37. Queuosine was mapped to its normal position 34 in tRNA HisGUG and tRNA Asp GUC. (Figure S2M,N) with spectral complexity hindering the mapping of tRNA Tyr and tRNA Asn. Finally, the modification inosine was mapped to position 34 of the tRNA Arg ACG (Figures S1 and S2O). Due to either low abundance or spectral complexity we were unable to map all the modifications detected through nucleoside mass spectrometric analysis to an oligonucleotide sequence. For example, the modification cmnm⁵s²Um and ms²t⁶A were mapped to the anticodon of tRNA Leu UAA, however the complexity of the MS/MS spectra makes this assignation tentative. The presence of 2'-O-methylation at the modified U*34 in tRNA Leu UAA was also confirmed by RiboMethSeq (Figure S1). Mapped modification position for those which were detected are listed in Table 1.

<table><caption>**Table 1.** Modified ribonucleosides identified in *B. subtilis* tRNAs.</caption><thead><tr><th>Detected in Bulk Nucleoside Analysis and Mapped by RiboMethSeq</th><th>Previously Mapped</th><th>Mapped by MS</th><th>Mapped by NGS</th><th>Pathway Intermediate</th></tr></thead><tbody><tr><td colspan="5"><strong>tRNA modifications</strong></td></tr><tr><td>1-methyladenosine, m<sup>1</sup>A</td><td>M%</td><td>Glu22 UUC<sup>1</sup><br>Tyr22 GUA<sup>1</sup></td><td>Cys22 GCA<br>Gln22 UUG<br>Glu22 UUC<br>Leu22 CAA<br>Leu22 UAA<br>Leu22 GAG<br>Tyr22 GUA<br>Ser22 GGA<br>Ser22 GCT<br>Ser22 UGA</td><td></td></tr><tr><td>1-methylguanosine, m<sup>1</sup>G</td><td>M</td><td>Pro37 UGG<br>Leu37 CAG</td><td></td><td></td></tr><tr><td>2-methyladenosine, m<sup>2</sup>A</td><td></td><td>Glu37 UUC<br>Leu37 UAA</td><td></td><td></td></tr><tr><td>2-methylthio-N6-isopentenyladenosine, ms<sup>2</sup>t<sup>6</sup>A</td><td>M</td><td></td><td></td><td></td></tr><tr><td>2-methylthio-N6-threonylcarbamoyladenosine, ms<sup>2</sup>t<sup>6</sup>A</td><td>M<sup>2</sup></td><td>Ser37 UGA</td><td></td><td></td></tr><tr><td>2'-O-methylcytidine, Cm</td><td></td><td></td><td>Leu34 CAA</td><td></td></tr><tr><td>2'-O-methylguanosine, Gm</td><td>M</td><td></td><td>Phe34 GAA</td><td></td></tr></tbody></table>
modification, since a number of candidate Am sites were found (positions 224 and 1455 in 16S rRNA as well as 1075, 1524, 2220, and 2827 in 23S rRNA) (Figure S1).

Several observed modifications must be artifacts of the sample preparation. m¹I has only been observed in archaeal tRNA [62] but could arise from the chemical or enzymatic deamination of m¹A during the enzymatic digestion to nucleosides [63]. Similarly, oxidation inside the cell or during analysis is most certainly the source of 8-oxoG formation [64].

The base position cannot be determined for simple methylations by oligonucleotide mass mapping; for example, m¹A, m²A, and m⁶A are all detected as only a methylated adenosine. We were able to map an adenosine methylation at position 37 on tRNA MetCAU, as well as ac⁴C at position 34 (Figure S2A). The sequence of the *E. coli* equivalent tRNA suggests that the probable modification is m⁶A. A methylated adenosine was also mapped to position 22 of the tRNA Leu UAA (Figure S2B). This adenosine was identified as m¹A, based on prior literature [65], and this was confirmed by RT-signature analysis (Figure S1 and Table S3). We also mapped the m²A found at position 37 of tRNA Glu UUC, having mnm⁵s²U at U34 (Figure S2C). The presence of m¹A previously reported at position 58 of tRNA Leu CAG [66] was never reproduced (Louis Droogmans, personal communication) and was not detected in our analyses. The unknown C32 modification in tRNA LysTTT reported by Vold et al. [67] was not detected by either the LC-MS or the RiboMethSeq approach.

In addition, we mapped mnm⁵s²U (Figure S2D), as well as mnm⁵U (Figure S2E) and cmnm⁵s²U (Figure S2F), to different tRNAs (Table 1). The modification 5-methoxyluridine (mo⁵U) is found in a number of *E. coli* tRNAs at position 34. Here, we were able to map mo⁵U position 34 of tRNA Leu UAG (Figure S2G). Moreover, 4-thiouridine was mapped to position 8 of tRNA GlnUUG by MS (Figure S2H) and of tRNA Leu UAA by RT (Figure S1). The hypermodification ct⁶A, along with the t⁶A counterpart, was found in tRNA Thr UGU (Figure S2I,J) and tRNA Ile GAU (Figure S2K,L) at position 37. Queuosine was mapped to its normal position 34 in tRNA HisGUG and tRNA Asp GUC. (Figure S2M,N) with spectral complexity hindering the mapping of tRNA Tyr and tRNA Asn. Finally, the modification inosine was mapped to position 34 of the tRNA Arg ACG (Figures S1 and S2O). Due to either low abundance or spectral complexity we were unable to map all the modifications detected through nucleoside mass spectrometric analysis to an oligonucleotide sequence. For example, the modification cmnm⁵s²Um and ms²t⁶A were mapped to the anticodon of tRNA Leu UAA, however the complexity of the MS/MS spectra makes this assignation tentative. The presence of 2'-O-methylation at the modified U*34 in tRNA Leu UAA was also confirmed by RiboMethSeq (Figure S1). Mapped modification position for those which were detected are listed in Table 1.

<table><caption>**Table 1.** Modified ribonucleosides identified in *B. subtilis* tRNAs.</caption><thead><tr><th>Detected in Bulk Nucleoside Analysis and Mapped by RiboMethSeq</th><th>Previously Mapped</th><th>Mapped by MS</th><th>Mapped by NGS</th><th>Pathway Intermediate</th></tr></thead><tbody><tr><td colspan="5"><strong>tRNA modifications</strong></td></tr><tr><td>1-methyladenosine, m¹A</td><td>M%</td><td>Glu22 UUC¹<br>Tyr22 GUA¹</td><td>Cys22 GCA<br>Gln22 UUG<br>Glu22 UUC<br>Leu22 CAA<br>Leu22 UAA<br>Leu22 GAG<br>Tyr22 GUA<br>Ser22 GGA<br>Ser22 GCT<br>Ser22 UGA</td><td></td></tr><tr><td>1-methylguanosine, m¹G</td><td>M</td><td>Pro37 UGG<br>Leu37 CAG</td><td></td><td></td></tr><tr><td>2-methyladenosine, m²A</td><td></td><td>Glu37 UUC<br>Leu37 UAA</td><td></td><td></td></tr><tr><td>2-methylthio-N6-isopentenyladenosine, ms²i⁶A</td><td>M</td><td></td><td></td><td></td></tr><tr><td>2-methylthio-N6-threonylcarbamoyladenosine, ms²t⁶A</td><td>M²</td><td>Ser37 UGA</td><td></td><td></td></tr><tr><td>2'-O-methylcytidine, Cm</td><td></td><td></td><td>Leu34 CAA</td><td></td></tr><tr><td>2'-O-methylguanosine, Gm</td><td>M</td><td></td><td>Phe34 GAA</td><td></td></tr></tbody></table>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th>Detected in Bulk Nucleoside Analysis and Mapped by RiboMethSeq</th><th>Previously Mapped</th><th>Mapped by MS</th><th>Mapped by NGS</th><th>Pathway Intermediate</th></tr></thead><tbody><tr><td>Modified U* with 2'-O-methyluridine, U*m</td><td></td><td></td><td>Lys34 UUU</td><td></td></tr><tr><td>Modified U* with 2'-O-methyluridine, U*m</td><td></td><td></td><td>Gln34 UUG</td><td></td></tr><tr><td></td><td></td><td></td><td>Leu8 UAA</td><td></td></tr><tr><td></td><td></td><td></td><td>Gln8 UUG</td><td></td></tr><tr><td></td><td></td><td></td><td>Tyr8 GUA</td><td></td></tr><tr><td>4-thiouridine, s<sup>4</sup>U</td><td>M</td><td>Gln8 UUG</td><td></td><td></td></tr><tr><td>5-aminomethyluridine, nm<sup>5</sup>U</td><td></td><td></td><td></td><td></td></tr><tr><td>5-carboxymethylaminomethyl-2'-O-methyluridine, cmm<sup>5</sup>Um</td><td>p& [49]</td><td>Leu34 UAA</td><td>Leu34 UAA</td><td>P [49]</td></tr><tr><td>5-carboxymethylaminomethyl-2-thiouridine, cmm<sup>5</sup>s<sup>2</sup>U</td><td>M/[49]</td><td>Gln34 UUG</td><td></td><td></td></tr><tr><td>5-methoxyuridine, mo<sup>5</sup>U</td><td>M</td><td>Leu34 UAG</td><td></td><td></td></tr><tr><td>5-methylaminomethyl-2-thiouridine, mnm<sup>5</sup>s<sup>2</sup>U</td><td>[49]</td><td>Glu34 UUC</td><td></td><td></td></tr><tr><td>5-methylaminomethyluridine, mnm<sup>5</sup>U</td><td></td><td>Gln34 UUG</td><td></td><td></td></tr><tr><td>5-methyluridine, m<sup>5</sup>U</td><td>M</td><td>Gln34 UUG</td><td></td><td>[49]</td></tr><tr><td>7-methylguanosine, m<sup>7</sup>G</td><td>M</td><td></td><td></td><td></td></tr><tr><td>cyclic N6-threonylcarbamoyladenosine, ct<sup>6</sup>A</td><td>M<sup>2</sup></td><td>Leu37 UAA</td><td></td><td></td></tr><tr><td></td><td></td><td>Ile37 GAU</td><td></td><td></td></tr><tr><td></td><td></td><td>Thr37 UGU</td><td></td><td></td></tr><tr><td></td><td></td><td>Asp21-22 GUC</td><td></td><td></td></tr><tr><td>Dihydouridine, D</td><td>M</td><td>Ile21-22 GAU</td><td>X</td><td></td></tr><tr><td></td><td></td><td>Ile21-22 CAU</td><td></td><td></td></tr><tr><td>Epoxyqueuosine, oQ</td><td></td><td></td><td></td><td></td></tr><tr><td>Inosine, I</td><td>M</td><td>Arg34 ACG</td><td>Arg34 ACG</td><td>X</td></tr><tr><td>Lysidine, k<sup>2</sup>C</td><td>M</td><td></td><td></td><td></td></tr><tr><td>N4-acetylcytidine, ac<sup>4</sup>C</td><td>[50]</td><td>Met34 CAU</td><td></td><td></td></tr><tr><td>N6-isopentenyladenosine, i<sup>6</sup>A</td><td>M</td><td></td><td></td><td></td></tr><tr><td>N6-methyladenosine, m<sup>6</sup>A</td><td>M</td><td>Met37 CAU</td><td></td><td></td></tr><tr><td>N6-threonylcarbamoyladenosine, t<sup>6</sup>A</td><td></td><td>Ile37 GAU</td><td></td><td>M<sup>3</sup></td></tr><tr><td></td><td></td><td>Thr37 UGU</td><td></td><td></td></tr><tr><td>Pseudouridine, ψ</td><td>M</td><td>His34 GUG</td><td>many, ψ 31/32/55</td><td></td></tr><tr><td>Queuosine, Q</td><td>M</td><td>Asp34 GUC</td><td></td><td></td></tr><tr><td><strong>rRNA modifications</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>2-methylguanosine, m<sup>2</sup>G</td><td></td><td></td><td></td><td></td></tr><tr><td>N4,2'-O-dimethylcytidine, m<sup>4</sup>C<sub>m</sub></td><td></td><td></td><td>16S rRNA 1410</td><td></td></tr><tr><td>N2,N2-dimethylguanosine, m<sup>2,2</sup>G</td><td></td><td></td><td>putative</td><td></td></tr><tr><td>N6,N6-dimethyladenosine, m<sup>6,6</sup>A</td><td></td><td>Yes</td><td>16S rRNA1527/28</td><td></td></tr><tr><td>2'-O-methyladenosine, Am</td><td></td><td></td><td>putative</td><td></td></tr><tr><td><strong>Chemical or biochemical artifacts</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>1-methylinosine, m<sup>1</sup>I</td><td></td><td></td><td></td><td></td></tr><tr><td>8-oxoguanosine, oxoG</td><td></td><td></td><td></td><td></td></tr></tbody></table>

## 3.2. Compilation of Predicted and Validated *B. subtilis* tRNA Modification Genes

Combining literature searches with sequence similarity searches using the set of known enzymes of *E. coli* as input, and crossing these data with the list of tRNA modification genes in *SubtiWiki* [24], led to a list of tRNA modification genes in *B. subtilis* (Figure 3 and Table S3). The level of confidence for the various predictions was given by a score that ranged from (5) to (1) (see Figure 3 legend for the full description of the scoring system). This analysis led to the identification of 41 predicted tRNA modification genes in *B. subtilis* that will be discussed in more detail below.

<footnote n="9">⁹M, listed in Modomics; &P, predicted. ¹Note 1: m¹A found in low resolution mapping data but not independently verified in high resolution data. ²Note 2: ms²ct⁶A37 previously reported [68] but not detected here. ³Note 3: t⁶A is an intermediate but is also found when the tRNA is prepared in standard conditions and hence this is the molecule reported in Modomics.</footnote>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th>Detected in Bulk Nucleoside Analysis and Mapped by RiboMethSeq</th><th>Previously Mapped</th><th>Mapped by MS</th><th>Mapped by NGS</th><th>Pathway Intermediate</th></tr></thead><tbody><tr><td>Modified U* with 2'-O-methyluridine, U*m</td><td></td><td></td><td>Lys34 UUU</td><td></td></tr><tr><td>Modified U* with 2'-O-methyluridine, U*m</td><td></td><td></td><td>Gln34 UUG</td><td></td></tr><tr><td></td><td></td><td></td><td>Leu8 UAA</td><td></td></tr><tr><td></td><td></td><td></td><td>Gln8 UUG</td><td></td></tr><tr><td></td><td></td><td></td><td>Tyr8 GUA</td><td></td></tr><tr><td>4-thiouridine, s<sup>4</sup>U</td><td>M</td><td>Gln8 UUG</td><td></td><td></td></tr><tr><td>5-aminomethyluridine, nm<sup>5</sup>U</td><td></td><td></td><td></td><td></td></tr><tr><td>5-carboxymethylaminomethyl-2'-O-methyluridine, cmm<sup>5</sup>Um</td><td>p& [49]</td><td>Leu34 UAA</td><td>Leu34 UAA</td><td>P [49]</td></tr><tr><td>5-carboxymethylaminomethyl-2-thiouridine, cmm<sup>5</sup>s<sup>2</sup>U</td><td>M/[49]</td><td>Gln34 UUG</td><td></td><td></td></tr><tr><td>5-methoxyuridine, mo<sup>5</sup>U</td><td>M</td><td>Leu34 UAG</td><td></td><td></td></tr><tr><td>5-methylaminomethyl-2-thiouridine, mnm<sup>5</sup>s<sup>2</sup>U</td><td>[49]</td><td>Glu34 UUC</td><td></td><td></td></tr><tr><td>5-methylaminomethyluridine, mnm<sup>5</sup>U</td><td></td><td>Gln34 UUG</td><td></td><td></td></tr><tr><td>5-methyluridine, m<sup>5</sup>U</td><td>M</td><td>Gln34 UUG</td><td></td><td>[49]</td></tr><tr><td>7-methylguanosine, m<sup>7</sup>G</td><td>M</td><td></td><td></td><td></td></tr><tr><td>cyclic N6-threonylcarbamoyladenosine, ct<sup>6</sup>A</td><td>M<sup>2</sup></td><td>Leu37 UAA</td><td></td><td></td></tr><tr><td></td><td></td><td>Ile37 GAU</td><td></td><td></td></tr><tr><td></td><td></td><td>Thr37 UGU</td><td></td><td></td></tr><tr><td></td><td></td><td>Asp21-22 GUC</td><td></td><td></td></tr><tr><td>Dihydouridine, D</td><td>M</td><td>Ile21-22 GAU</td><td>X</td><td></td></tr><tr><td></td><td></td><td>Ile21-22 CAU</td><td></td><td></td></tr><tr><td>Epoxyqueuosine, oQ</td><td></td><td></td><td></td><td></td></tr><tr><td>Inosine, I</td><td>M</td><td>Arg34 ACG</td><td>Arg34 ACG</td><td>X</td></tr><tr><td>Lysidine, k<sup>2</sup>C</td><td>M</td><td></td><td></td><td></td></tr><tr><td>N4-acetylcytidine, ac<sup>4</sup>C</td><td>[50]</td><td>Met34 CAU</td><td></td><td></td></tr><tr><td>N6-isopentenyladenosine, i<sup>6</sup>A</td><td>M</td><td></td><td></td><td></td></tr><tr><td>N6-methyladenosine, m<sup>6</sup>A</td><td>M</td><td>Met37 CAU</td><td></td><td></td></tr><tr><td>N6-threonylcarbamoyladenosine, t<sup>6</sup>A</td><td></td><td>Ile37 GAU</td><td></td><td>M<sup>3</sup></td></tr><tr><td></td><td></td><td>Thr37 UGU</td><td></td><td></td></tr><tr><td>Pseudouridine, ψ</td><td>M</td><td>His34 GUG</td><td>many, ψ 31/32/55</td><td></td></tr><tr><td>Queuosine, Q</td><td>M</td><td>Asp34 GUC</td><td></td><td></td></tr><tr><td><strong>rRNA modifications</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>2-methylguanosine, m<sup>2</sup>G</td><td></td><td></td><td></td><td></td></tr><tr><td>N4,2'-O-dimethylcytidine, m<sup>4</sup>C<sub>m</sub></td><td></td><td></td><td>16S rRNA 1410</td><td></td></tr><tr><td>N2,N2-dimethylguanosine, m<sup>2,2</sup>G</td><td></td><td></td><td>putative</td><td></td></tr><tr><td>N6,N6-dimethyladenosine, m<sup>6,6</sup>A</td><td></td><td>Yes</td><td>16S rRNA1527/28</td><td></td></tr><tr><td>2'-O-methyladenosine, Am</td><td></td><td></td><td>putative</td><td></td></tr><tr><td><strong>Chemical or biochemical artifacts</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>1-methylinosine, m<sup>1</sup>I</td><td></td><td></td><td></td><td></td></tr><tr><td>8-oxoguanosine, oxoG</td><td></td><td></td><td></td><td></td></tr></tbody></table>

## 3.2. Compilation of Predicted and Validated *B. subtilis* tRNA Modification Genes

Combining literature searches with sequence similarity searches using the set of known enzymes of *E. coli* as input, and crossing these data with the list of tRNA modification genes in *SubtiWiki* [24], led to a list of tRNA modification genes in *B. subtilis* (Figure 3 and Table S3). The level of confidence for the various predictions was given by a score that ranged from (5) to (1) (see Figure 3 legend for the full description of the scoring system). This analysis led to the identification of 41 predicted tRNA modification genes in *B. subtilis* that will be discussed in more detail below.

<footnote n="9">⁹M, listed in Modomics; &P, predicted. ¹Note 1: m¹A found in low resolution mapping data but not independently verified in high resolution data. ²Note 2: ms²ct⁶A37 previously reported [68] but not detected here. ³Note 3: t⁶A is an intermediate but is also found when the tRNA is prepared in standard conditions and hence this is the molecule reported in Modomics.</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** *B. subtilis* tRNA modifications with predicted and validated corresponding enzymes. The four modifications in red were found in the current analysis but were not present in the *B. subtilis* tRNA sequences listed in Modomics. Proteins with stars (*) are encoded by essential genes in *B. subtilis*. Proteins in green are non-orthologous to the enzyme catalyzing the same reaction in *E. coli*. Proteins in blue insert modifications that are not found in *E. coli*. The confidence score for the prediction is as follows: (5), validated in vitro with the *B. subtilis* enzymes and in vivo with the *B. subtilis* deletion mutant; (4), orthologous to a validated enzyme in other species, genetic data in *B. subtilis*; (3), orthologous to a validated enzyme in other species, with mapping data in *B. subtilis*; (2), orthologous to a validated enzyme in other species, without mapping data in *B. subtilis*; (1), similarity with experimentally validated in other species but paralogs.</figcaption></figure>

### 3.3. The Majority of tRNA Modifications Are Introduced by Orthologs in *E. coli* and *B. subtilis* with Some Differences Observed in the Complex Pathways

Many enzymes of the tRNA modification machinery are conserved between *E. coli* and *B. subtilis*. Indeed, out of the 41 proteins predicted to be involved in inserting tRNA modifications into the Gram-positive model, the majority (32/41) are clear orthologs of *E. coli* proteins. These include proteins involved in the synthesis of I34 (TadA), k²C34 (TilS), m¹G37 (TrmD), i⁶A (MiaA) ms²i⁶A37 (MiaB/YmcA), ψ38/39 (TruA), m⁷G46 (TrmB), ψ55 (TruB), and m⁶A37 (TrmN/YabB), for which there is also experimental evidence of the presence of the modifications at these positions in the tRNA (Table 1). For RlmN/YloN, no tRNA sequence containing this modification is found in Modomics, but our analysis revealed a possible m²A at position 37 in tRNAGluUUC (see discussion above and Figure S2C).

For the enzymes involved in ribose methylation at position 34, the *B. subtilis* CspR is an ortholog of the *E. coli* TrmL, but U*m/Cm* are found in *E. coli* tRNAs, while U*m/Gm/Cm* were identified in the *B. subtilis* tRNAs by RiboMethSeq. We confirmed the 2'-O-methylation of Gₘ34 in tRNAPheGAA and identified Cₘ34 in tRNALeuCAA and U*m34 in tRNAGlnUUG and tRNALeuUAA (Figure S1). No signals corresponding to the 2'-O-methylation inserted in *E. coli* by the TrmH and TrmJ enzymes (Gm at position 18 and Um/Cm at position 32, respectively) (Figure 1) were identified in the *B. subtilis* tRNAs (data not shown).

In the cases of most multi-step pathways, such as those involved in s⁴U34, Q34, t⁶A37, and XM⁵s²U34 syntheses, many small variations can be found between the two model organisms; even the enzymes involved in the core machineries of these pathways are conserved. The Q synthesis
<figure><img src="image_6.png" /><figcaption>**Figure 3.** *B. subtilis* tRNA modifications with predicted and validated corresponding enzymes. The four modifications in red were found in the current analysis but were not present in the *B. subtilis* tRNA sequences listed in Modomics. Proteins with stars (*) are encoded by essential genes in *B. subtilis*. Proteins in green are non-orthologous to the enzyme catalyzing the same reaction in *E. coli*. Proteins in blue insert modifications that are not found in *E. coli*. The confidence score for the prediction is as follows: (5), validated in vitro with the *B. subtilis* enzymes and in vivo with the *B. subtilis* deletion mutant; (4), orthologous to a validated enzyme in other species, genetic data in *B. subtilis*; (3), orthologous to a validated enzyme in other species, with mapping data in *B. subtilis*; (2), orthologous to a validated enzyme in other species, without mapping data in *B. subtilis*; (1), similarity with experimentally validated in other species but paralogs.</figcaption></figure>

### 3.3. The Majority of tRNA Modifications Are Introduced by Orthologs in *E. coli* and *B. subtilis* with Some Differences Observed in the Complex Pathways

Many enzymes of the tRNA modification machinery are conserved between *E. coli* and *B. subtilis*. Indeed, out of the 41 proteins predicted to be involved in inserting tRNA modifications into the Gram-positive model, the majority (32/41) are clear orthologs of *E. coli* proteins. These include proteins involved in the synthesis of I34 (TadA), k²C34 (TilS), m¹G37 (TrmD), i⁶A (MiaA) ms²i⁶A37 (MiaB/YmcA), ψ38/39 (TruA), m⁷G46 (TrmB), ψ55 (TruB), and m⁶A37 (TrmN/YabB), for which there is also experimental evidence of the presence of the modifications at these positions in the tRNA (Table 1). For RlmN/YloN, no tRNA sequence containing this modification is found in Modomics, but our analysis revealed a possible m²A at position 37 in tRNAGluUUC (see discussion above and Figure S2C).

For the enzymes involved in ribose methylation at position 34, the *B. subtilis* CspR is an ortholog of the *E. coli* TrmL, but U*m/Cm* are found in *E. coli* tRNAs, while U*m/Gm/Cm* were identified in the *B. subtilis* tRNAs by RiboMethSeq. We confirmed the 2'-O-methylation of Gₘ34 in tRNAPheGAA and identified Cₘ34 in tRNALeuCAA and U*m34 in tRNAGlnUUG and tRNALeuUAA (Figure S1). No signals corresponding to the 2'-O-methylation inserted in *E. coli* by the TrmH and TrmJ enzymes (Gm at position 18 and Um/Cm at position 32, respectively) (Figure 1) were identified in the *B. subtilis* tRNAs (data not shown).

In the cases of most multi-step pathways, such as those involved in s⁴U34, Q34, t⁶A37, and XM⁵s²U34 syntheses, many small variations can be found between the two model organisms; even the enzymes involved in the core machineries of these pathways are conserved. The Q synthesis
pathways are identical in *E. coli* and *B. subtilis*, with every step catalyzed by an orthologous enzyme. The only differences are that an additional enzyme (FolE2) can catalyze the first step in *B. subtilis* [69] and there is no ortholog of the GluQRS/YadB enzyme in *B. subtilis* [70]. The t⁶A37 synthesis core machineries are also very similar in the two model organisms, with differences found both at the beginning and end of the pathways. The enzymes that catalyze the first step are slightly different: TsaC in *E. coli* [71,72] and TsaC2, which carries an extra SUA5 domain, in *B. subtilis* [71,73]. In addition, different t⁶A hypermodifications are found in the two organisms. Whereas m⁶t⁶A is present in *E. coli* [74], ms²t⁶A is found *B. subtilis*, synthesized by the MiaB paralog MtaB [75,76]. An ortholog of the *E. coli* ct⁶A formation enzyme TcdA was validated in *B. subtilis* and both ct⁶A and ms²ct⁶A can be detected in this organism [68], even if we detected only ct⁶A in the current study (Table 1). It is not clear if additional sulfur transfer proteins are required in *B. subtilis*, like they are in *E. coli* (CsdAE) [53]. Thexm⁵U derivatives identified in the *B. subtilis* tRNAs (Table 1) are similar to those present in *E. coli*, as confirmed by the Armengod laboratory [49]. The exception is cm⁵s²U, which might be an intermediate as the full synthesis of these complex modifications varies greatly depending on the media conditions [49]. Orthologs of two of the three *E. coli* enzymes involved in synthesizing this complex modification (MnmE, MnmG) are found in *B. subtilis*. However, no MnmC ortholog is present, an observation already made by the Armengod group, who also showed that another enzyme catalyzing the same activity must be present and that it does not require S-adenosyl-L-methionine [49]. There are also differences in the thiolation enzymes which synthesize the xs²U34 derivatives. If the enzyme that introduces the sulfur in tRNA, MnmA, is similar in both organisms, *E. coli* uses a whole sulfur relay system (TusABCDE) [77], while *B. subtilis* uses only one relay enzyme of the IscS family, YrvO/IscS1 [78,79]. Finally, in *E. coli*, the ribose methylase required to synthesize cmnm⁵Um is the multifunctional TrmL; it is therefore parsimonious to propose that, similarly, CspR catalyzes this reaction in *B. subtilis*, even though this has not been confirmed experimentally.

### 3.4. Modifications and Enzymes Specific to *Bacillus subtilis*

In addition to ms²t⁶A, another modification that is present only in *B. subtilis* and not *E. coli* is m¹A22. This modification is inserted by the TrmK gene, which was described over 10 years ago [80] and is a marker of Gram-positive bacteria, as it is even found in those bacteria with minimal genomes, such as *M. capricolum* [18]. G<sub>m</sub>34 in tRNA<sup>Phe</sup>GAA (Table S3) is also found in *B. subtilis* but not *E. coli* and is reminiscent of eukaryotic tRNA modification profiles [81–83].

Moreover, beyond the MnmC discussed above, where a non-homologous enzyme, still to be discovered, must catalyze the same reaction in *B. subtilis*, several cases of non-orthologous replacement are found between the two model organisms. The first is the well-studied replacement of the SAM dependent *E. coli* methylase TrmA (COG2265) by the *B. subtilis* folate dependent enzyme TrmFO (COG1206) in order to synthesize m⁵U54 [84]. Another case has recently been discovered. In the synthesis of ac⁴C34, the *E. coli* TmcA acetyltransferase of the COG1444 family, which uses acetyl-coA donor for the acetyl group [85], has been replaced by the TmcAL (COG1323) protein in *B. subtilis*, which uses acetate and ATP to create an acetyl-AMP intermediate as the acetyl donor in different mechanisms [50].

The two other non-orthologous displacement examples are slightly more complex. Both *E. coli* and *B. subtilis* harbor the s⁴U8 modification, but the modification enzymes are slightly different [78,86]. The *B. subtilis* Thil is shorter than the *E. coli* homolog, lacks the rhodanese domain, and is not involved in thiamin synthesis. In addition, Gram-positive organisms harbor dedicated cysteine desulfurase for specific thiofactor generation, whereas *E. coli* uses the generalist IscS for all of them. NifZ is the one dedicated cysteine desulfurase that is involved in s⁴U synthesis in *B. subtilis* [86]. Finally, different enzymes modify the ho⁵U precursor in order to give different xmo⁵U34 derivatives. In *E. coli*, the two CmoA and CmoB enzymes synthesize cmo⁵U34 in a complex reaction that uses a S-adenosyl-S-carboxymethyl-l-homocysteine (Cx-SAM) intermediate [87]; the cmo⁵U34 intermediate can be further modified to mcmo⁵U, and then mcmo⁵Um, by CmoM and TrmL, respectively [88].
pathways are identical in *E. coli* and *B. subtilis*, with every step catalyzed by an orthologous enzyme. The only differences are that an additional enzyme (FolE2) can catalyze the first step in *B. subtilis* [69] and there is no ortholog of the GluQRS/YadB enzyme in *B. subtilis* [70]. The t⁶A37 synthesis core machineries are also very similar in the two model organisms, with differences found both at the beginning and end of the pathways. The enzymes that catalyze the first step are slightly different: TsaC in *E. coli* [71,72] and TsaC2, which carries an extra SUA5 domain, in *B. subtilis* [71,73]. In addition, different t⁶A hypermodifications are found in the two organisms. Whereas m⁶t⁶A is present in *E. coli* [74], ms²t⁶A is found *B. subtilis*, synthesized by the MiaB paralog MtaB [75,76]. An ortholog of the *E. coli* ct⁶A formation enzyme TcdA was validated in *B. subtilis* and both ct⁶A and ms²ct⁶A can be detected in this organism [68], even if we detected only ct⁶A in the current study (Table 1). It is not clear if additional sulfur transfer proteins are required in *B. subtilis*, like they are in *E. coli* (CsdAE) [53]. Thexm⁵U derivatives identified in the *B. subtilis* tRNAs (Table 1) are similar to those present in *E. coli*, as confirmed by the Armengod laboratory [49]. The exception is cm⁵s²U, which might be an intermediate as the full synthesis of these complex modifications varies greatly depending on the media conditions [49]. Orthologs of two of the three *E. coli* enzymes involved in synthesizing this complex modification (MnmE, MnmG) are found in *B. subtilis*. However, no MnmC ortholog is present, an observation already made by the Armengod group, who also showed that another enzyme catalyzing the same activity must be present and that it does not require S-adenosyl-L-methionine [49]. There are also differences in the thiolation enzymes which synthesize the xs²U34 derivatives. If the enzyme that introduces the sulfur in tRNA, MnmA, is similar in both organisms, *E. coli* uses a whole sulfur relay system (TusABCDE) [77], while *B. subtilis* uses only one relay enzyme of the IscS family, YrvO/IscS1 [78,79]. Finally, in *E. coli*, the ribose methylase required to synthesize cmnm⁵Um is the multifunctional TrmL; it is therefore parsimonious to propose that, similarly, CspR catalyzes this reaction in *B. subtilis*, even though this has not been confirmed experimentally.

### 3.4. Modifications and Enzymes Specific to *Bacillus subtilis*

In addition to ms²t⁶A, another modification that is present only in *B. subtilis* and not *E. coli* is m¹A22. This modification is inserted by the TrmK gene, which was described over 10 years ago [80] and is a marker of Gram-positive bacteria, as it is even found in those bacteria with minimal genomes, such as *M. capricolum* [18]. G<sub>m</sub>34 in tRNA<sup>Phe</sup>GAA (Table S3) is also found in *B. subtilis* but not *E. coli* and is reminiscent of eukaryotic tRNA modification profiles [81–83].

Moreover, beyond the MnmC discussed above, where a non-homologous enzyme, still to be discovered, must catalyze the same reaction in *B. subtilis*, several cases of non-orthologous replacement are found between the two model organisms. The first is the well-studied replacement of the SAM dependent *E. coli* methylase TrmA (COG2265) by the *B. subtilis* folate dependent enzyme TrmFO (COG1206) in order to synthesize m⁵U54 [84]. Another case has recently been discovered. In the synthesis of ac⁴C34, the *E. coli* TmcA acetyltransferase of the COG1444 family, which uses acetyl-coA donor for the acetyl group [85], has been replaced by the TmcAL (COG1323) protein in *B. subtilis*, which uses acetate and ATP to create an acetyl-AMP intermediate as the acetyl donor in different mechanisms [50].

The two other non-orthologous displacement examples are slightly more complex. Both *E. coli* and *B. subtilis* harbor the s⁴U8 modification, but the modification enzymes are slightly different [78,86]. The *B. subtilis* Thil is shorter than the *E. coli* homolog, lacks the rhodanese domain, and is not involved in thiamin synthesis. In addition, Gram-positive organisms harbor dedicated cysteine desulfurase for specific thiofactor generation, whereas *E. coli* uses the generalist IscS for all of them. NifZ is the one dedicated cysteine desulfurase that is involved in s⁴U synthesis in *B. subtilis* [86]. Finally, different enzymes modify the ho⁵U precursor in order to give different xmo⁵U34 derivatives. In *E. coli*, the two CmoA and CmoB enzymes synthesize cmo⁵U34 in a complex reaction that uses a S-adenosyl-S-carboxymethyl-l-homocysteine (Cx-SAM) intermediate [87]; the cmo⁵U34 intermediate can be further modified to mcmo⁵U, and then mcmo⁵Um, by CmoM and TrmL, respectively [88].
In *B. subtilis*, a classical methyltransferase (TrmR) directly methylates ho⁵U34 to mo⁵U34 [89], and no orthologs of CmoM or TrmL are found. The enzymes involved in ho⁵U synthesis have been identified recently [90]. In *E. coli*, TrhP (formerly known as YegQ), a peptidase U32 family protein, is involved in prephenate-dependent ho⁵U formation, while TrhO (formerly known as YceA), a rhodanese family protein, catalyzes oxygen-dependent ho⁵U formation and bypasses cmo⁵U biogenesis in a subset of tRNAs under aerobic conditions. In *B. subtilis*, the TrhP homologs TrhP1/YrrN and TrhP2/YrrO are also involved in one pathway, but another pathway has also been predicted and the corresponding genes have not yet been identified [91].

### 3.5. Paralog Issues Make Direct Annotation Transfers Problematic for a Few Modifications

There are several cases in which the homologs of validated *E. coli* tRNA modification enzymes are found in *B. subtilis* but, because of the paralog issue, it is not possible to give a high confidence functional annotation without additional experiments. This is the case for the identification of the *B. subtilis* dihydrourine insertion enzymes. In *E. coli*, DusC modifies position 16, DusB modifies position 17, and DusA modifies positions 20 and 20a [92]. The three proteins are members of the same COG0042 family. Dihydouridine is found at four different positions in *B. subtilis* tRNAs: D17, D20, D20a, and D47 (Figure 3). However, only two COG0042 members are encoded in the *B. subtilis* genome: Dus1 (BSU00810) and Dus2 (BSU08030). At this stage, without further biochemical and/or genetic experiments, it is not possible to predict which position each enzyme modifies.

Another case is the prediction of the *B. subtilis* pseudouridine synthase involved in the synthesis of ψ31/32 catalyzed in *E. coli* by RluA, an enzyme that modifies both rRNA and tRNA substrates. Searching the *B. subtilis* genome, using *E. coli* RluA as an input and using BlastP with default settings, identified the two RluA homologs YhcT and YjbO, both of which are COG0564 members. Here, again, it was impossible to predict if these proteins had the same substrate specificity as the *E. coli* RluA, or if they had each evolved to modify just tRNA or rRNA or just one position (31 or 32), or if they had another function in the cell. We therefore analyzed by Pseudo-Seq the bulk tRNA preparations extracted from *B. subtilis* strains lacking each gene, respectively. As shown in Figure 4 and Figure S3, the absence of YjbO led to the absence of ψ at position 31 in tRNA<sup>TrpCCA</sup>, tRNA<sup>CysGCA</sup>, tRNA<sup>ArgCCG</sup>, and tRNA<sup>ArgCCU</sup>, while these were not affected by the absence of YhcT. Similar results were observed for the presence of ψ at position 32 in tRNA<sup>ArgCCU</sup> and tRNA<sup>ProUGG</sup> (Figure S3), suggesting that YhbO was responsible for modifying both the 31 and 32 positions in the tRNAs.

Since YhcT did not seem to be involved in tRNA modification, we tested whether it could be modifying rRNA. Previous studies from the Ofengand laboratory [93,94] mapped ψ modifications in the 23S rRNAs of *E. coli* and *B. subtilis* and identified the corresponding *E. coli* genes. As shown in Table 2, 10 residues are modified to ψ in *E. coli* 23 S rRNA, but only four of those are also modified in *B. subtilis*. It was previously shown that the functional homolog of RluB in *B. subtilis* was YpuL/BSU23160, which was actually closer in sequence to the *E. coli* enzyme RsuA that modifies the small subunit [95]. PsiSeq analysis suggests that YhcT is the functional homolog of *E. coli* RluD and modifies positions 1940, 1944, and 1946 of the 23S rRNA, even if it is closer in sequence to RluA (Figure 4).
In *B. subtilis*, a classical methyltransferase (TrmR) directly methylates ho⁵U34 to mo⁵U34 [89], and no orthologs of CmoM or TrmL are found. The enzymes involved in ho⁵U synthesis have been identified recently [90]. In *E. coli*, TrhP (formerly known as YegQ), a peptidase U32 family protein, is involved in prephenate-dependent ho⁵U formation, while TrhO (formerly known as YceA), a rhodanese family protein, catalyzes oxygen-dependent ho⁵U formation and bypasses cmo⁵U biogenesis in a subset of tRNAs under aerobic conditions. In *B. subtilis*, the TrhP homologs TrhP1/YrrN and TrhP2/YrrO are also involved in one pathway, but another pathway has also been predicted and the corresponding genes have not yet been identified [91].

### 3.5. Paralog Issues Make Direct Annotation Transfers Problematic for a Few Modifications

There are several cases in which the homologs of validated *E. coli* tRNA modification enzymes are found in *B. subtilis* but, because of the paralog issue, it is not possible to give a high confidence functional annotation without additional experiments. This is the case for the identification of the *B. subtilis* dihydrourine insertion enzymes. In *E. coli*, DusC modifies position 16, DusB modifies position 17, and DusA modifies positions 20 and 20a [92]. The three proteins are members of the same COG0042 family. Dihydouridine is found at four different positions in *B. subtilis* tRNAs: D17, D20, D20a, and D47 (Figure 3). However, only two COG0042 members are encoded in the *B. subtilis* genome: Dus1 (BSU00810) and Dus2 (BSU08030). At this stage, without further biochemical and/or genetic experiments, it is not possible to predict which position each enzyme modifies.

Another case is the prediction of the *B. subtilis* pseudouridine synthase involved in the synthesis of ψ31/32 catalyzed in *E. coli* by RluA, an enzyme that modifies both rRNA and tRNA substrates. Searching the *B. subtilis* genome, using *E. coli* RluA as an input and using BlastP with default settings, identified the two RluA homologs YhcT and YjbO, both of which are COG0564 members. Here, again, it was impossible to predict if these proteins had the same substrate specificity as the *E. coli* RluA, or if they had each evolved to modify just tRNA or rRNA or just one position (31 or 32), or if they had another function in the cell. We therefore analyzed by Pseudo-Seq the bulk tRNA preparations extracted from *B. subtilis* strains lacking each gene, respectively. As shown in Figure 4 and Figure S3, the absence of YjbO led to the absence of ψ at position 31 in tRNA<sup>TrpCCA</sup>, tRNA<sup>CysGCA</sup>, tRNA<sup>ArgCCG</sup>, and tRNA<sup>ArgCCU</sup>, while these were not affected by the absence of YhcT. Similar results were observed for the presence of ψ at position 32 in tRNA<sup>ArgCCU</sup> and tRNA<sup>ProUGG</sup> (Figure S3), suggesting that YhbO was responsible for modifying both the 31 and 32 positions in the tRNAs.

Since YhcT did not seem to be involved in tRNA modification, we tested whether it could be modifying rRNA. Previous studies from the Ofengand laboratory [93,94] mapped ψ modifications in the 23S rRNAs of *E. coli* and *B. subtilis* and identified the corresponding *E. coli* genes. As shown in Table 2, 10 residues are modified to ψ in *E. coli* 23 S rRNA, but only four of those are also modified in *B. subtilis*. It was previously shown that the functional homolog of RluB in *B. subtilis* was YpuL/BSU23160, which was actually closer in sequence to the *E. coli* enzyme RsuA that modifies the small subunit [95]. PsiSeq analysis suggests that YhcT is the functional homolog of *E. coli* RluD and modifies positions 1940, 1944, and 1946 of the 23S rRNA, even if it is closer in sequence to RluA (Figure 4).
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Analysis of pseudouridine residues in tRNAs and rRNAs by PsiSeq. Normalized PsiSeq signals are shown for selected regions in *B. subtilis* tRNAs and 23S rRNA. The top row corresponds to the wild-type strain, the middle one for the *ΔyjbO* mutant, and the bottom one for the *ΔyhcT* mutant. Positions of known and detected pseudouridine residues are shown by red arrows. Conventional numbering of tRNA positions may not correspond to real position of nucleotide in the sequence due to missing residues and inclusion of 17a, 20a, 20b, and variable loop nucleotides. Identity of tRNA and its anticodon are indicated on the top of each panel. YjbO is involved in pseudouridylation of position 31 (and also 32) in several tRNAs. YhcT is involved in pseudouridylation of 23S rRNA positions 1940, 1944, and 1946 and is the functional equivalent of *E. coli* RluD.</figcaption></figure>

<table><caption>**Table 2.** Comparison of ψ sites and corresponding modification enzymes in 23S rRNA of *E. coli* and *B. subtilis*.</caption><thead><tr><th>E.c. Mod</th><th>B.s. Mod (n° in Ref [93,94])</th><th>E.c Protein</th><th>B.s Protein</th><th>Evidence</th></tr></thead><tbody><tr><td>Ψ746</td><td>U793 (791)</td><td>RluA</td><td>NA</td><td></td></tr><tr><td>Ψ955</td><td>U1001 (999)</td><td>RluC</td><td>NA</td><td></td></tr><tr><td>Ψ1911</td><td>Ψ1940 (1938)</td><td>RluD</td><td>YhcT/BSU09210</td><td>PsiSeq</td></tr><tr><td>Ψ1915*</td><td>Ψ1944* (1942)</td><td>RluD</td><td>YhcT/BSU09210</td><td>PsiSeq</td></tr><tr><td>Ψ1917</td><td>Ψ1946 (1944)</td><td>RluD</td><td>YhcT/BSU09210</td><td>PsiSeq</td></tr><tr><td>Ψ2457</td><td>U2486(2484)</td><td>RluE</td><td>NA</td><td></td></tr><tr><td>U2492</td><td>Ψ2521 (2520)</td><td>NA</td><td>YlmL/BSU15460</td><td>Only sequence similarity</td></tr><tr><td>Ψ2504</td><td>U2533 (2532)</td><td>RluC</td><td>NA</td><td></td></tr><tr><td>Ψ2580</td><td>U2609 (2608)</td><td>RluC</td><td>NA</td><td></td></tr><tr><td>Ψ2604</td><td>U2633 (2632)</td><td>RluF</td><td>NA</td><td></td></tr><tr><td>Ψ2605</td><td>Ψ2634 (2633)</td><td>RluB</td><td>YpuL/BSU23160</td><td>Experimental [95]</td></tr></tbody></table>

<footnote>* Certainly further modified to m³ψ, as attested by RT signature for this position (Figure S1).</footnote>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Analysis of pseudouridine residues in tRNAs and rRNAs by PsiSeq. Normalized PsiSeq signals are shown for selected regions in *B. subtilis* tRNAs and 23S rRNA. The top row corresponds to the wild-type strain, the middle one for the *ΔyjbO* mutant, and the bottom one for the *ΔyhcT* mutant. Positions of known and detected pseudouridine residues are shown by red arrows. Conventional numbering of tRNA positions may not correspond to real position of nucleotide in the sequence due to missing residues and inclusion of 17a, 20a, 20b, and variable loop nucleotides. Identity of tRNA and its anticodon are indicated on the top of each panel. YjbO is involved in pseudouridylation of position 31 (and also 32) in several tRNAs. YhcT is involved in pseudouridylation of 23S rRNA positions 1940, 1944, and 1946 and is the functional equivalent of *E. coli* RluD.</figcaption></figure>

<table><caption>**Table 2.** Comparison of ψ sites and corresponding modification enzymes in 23S rRNA of *E. coli* and *B. subtilis*.</caption><thead><tr><th>E.c. Mod</th><th>B.s. Mod (n° in Ref [93,94])</th><th>E.c Protein</th><th>B.s Protein</th><th>Evidence</th></tr></thead><tbody><tr><td>Ψ746</td><td>U793 (791)</td><td>RluA</td><td>NA</td><td></td></tr><tr><td>Ψ955</td><td>U1001 (999)</td><td>RluC</td><td>NA</td><td></td></tr><tr><td>Ψ1911</td><td>Ψ1940 (1938)</td><td>RluD</td><td>YhcT/BSU09210</td><td>PsiSeq</td></tr><tr><td>Ψ1915*</td><td>Ψ1944* (1942)</td><td>RluD</td><td>YhcT/BSU09210</td><td>PsiSeq</td></tr><tr><td>Ψ1917</td><td>Ψ1946 (1944)</td><td>RluD</td><td>YhcT/BSU09210</td><td>PsiSeq</td></tr><tr><td>Ψ2457</td><td>U2486(2484)</td><td>RluE</td><td>NA</td><td></td></tr><tr><td>U2492</td><td>Ψ2521 (2520)</td><td>NA</td><td>YlmL/BSU15460</td><td>Only sequence similarity</td></tr><tr><td>Ψ2504</td><td>U2533 (2532)</td><td>RluC</td><td>NA</td><td></td></tr><tr><td>Ψ2580</td><td>U2609 (2608)</td><td>RluC</td><td>NA</td><td></td></tr><tr><td>Ψ2604</td><td>U2633 (2632)</td><td>RluF</td><td>NA</td><td></td></tr><tr><td>Ψ2605</td><td>Ψ2634 (2633)</td><td>RluB</td><td>YpuL/BSU23160</td><td>Experimental [95]</td></tr></tbody></table>

<footnote>* Certainly further modified to m³ψ, as attested by RT signature for this position (Figure S1).</footnote>
This leaves one pseudouridine synthase to be identified for the modification of position 2521 (2520) in the 23S rRNA in *B. subtilis*. The equivalent position in *E. coli* (U2492) is not modified. A plausible candidate is the third COG0564 *B. subtilis* member, YlmL (BSU15460), a prediction that requires experimental validation. As a side note of this analysis, no ortholog of *E. coli* RluF was identified in *B. subtilis*, and in *E. coli*, this enzyme also modifies position 35 of tRNA<sup>Tyr</sup>QΨA [96]. We therefore predict that the corresponding *B. subtilis* tRNA is not modified with Ψ at this same position, and this is confirmed by the PsiSeq analysis (Table S3).

### 3.6. System Analysis of tRNA Modification Genes in *B. subtilis*

As *B. subtilis* is such a prominent model organism, many whole genome datasets are available. These include mutant collections, phenotypic screens, TnSeq datasets, and gene expression data [24]. Many of the tRNA modification genes are poorly annotated and/or were only recently discovered, so we revisited several published genome-wide studies that used 42 *B. subtilis* tRNA modification genes as input.

Seven out of the 42 *B. subtilis* modification genes are essential (Figure 3) [38]. *FolE/mtrA* is essential because of its role in the synthesis of the tetrahydrofolate cofactor, but not its role in Q synthesis, as it can be replaced by *folE2/yciA* only under zinc limitation [69]. *TilS* is essential in most bacteria [10], with the rare exception of organisms such as *Mycoplasma mobile* that have changed the minor tRNA<sup>Ile</sup> anticodon from CAU to UAU, making the lysidine modification dispensable [97]. The essentiality of *trmD* in many bacteria has already made it an antibacterial target [98]. Interestingly, two t<sup>6</sup>A37 synthesis genes are essential (*tsaB* and *tsaD*), whereas two others (*tsaC* and *tsaE*) are not, while all four orthologs are essential in *E. coli* [7]. This suggests that the essentiality of *tsaBD* in *B. subtilis* could be due to another role in the cell besides t<sup>6</sup>A formation in the tRNA. Finally, the *iscS1* and *mnmA* involved in s<sup>2</sup>U34 formation are essential in *B. subtilis* but not in *E. coli* [99], for reasons that are not totally clear.

A phenotype screen testing several conditions (growth, C or N sources, and temperature sensitivity) of two ordered, barcoded, erythromycin-resistance- and kanamycin-resistance-marked single-gene deletion libraries is available [38]. Fitness scores are given for every gene in every condition tested. We extracted from these data the scores for the 35 predicted tRNA modification genes that are not essential (Table S4), and we visualized these data using clustergrammer, after transforming them into Z-scores. While transforming the fitness data into Z-scores may have exaggerated the differences between the ErmR and KanR datasets, it was noted in their original publication [40] that discrepancies between datasets could be caused by the downstream effects of the resistance cassettes that were used to prepare each library of mutants. Since we observed differences in the fitness data between the ErmR and KanR datasets for our list of genes, only the ErmR dataset was used to prepare Figure 5. Any phenotype observations from these datasets should be experimentally validated.

As shown in Figure 5, genes in same modification pathways, such as *tsaC/tsaE* (t<sup>6</sup>A37), *tgt/queA* (Q), *queDCEF* (Q precursor), and cmnm<sup>5</sup>U34 (*mnmE/G*), have similar or the same fitness profiles, suggesting a biological relevance that will need to be further studied.

We also looked at published *B. subtilis* Tn-Seq studies and found that three tRNA modification genes have been identified as being involved in sporulation (*mnmE*, *mnmG*, and *miaA*) [100], and *mnmG* was also important for germination [101].
This leaves one pseudouridine synthase to be identified for the modification of position 2521 (2520) in the 23S rRNA in *B. subtilis*. The equivalent position in *E. coli* (U2492) is not modified. A plausible candidate is the third COG0564 *B. subtilis* member, YlmL (BSU15460), a prediction that requires experimental validation. As a side note of this analysis, no ortholog of *E. coli* RluF was identified in *B. subtilis*, and in *E. coli*, this enzyme also modifies position 35 of tRNA<sup>Tyr</sup>QΨA [96]. We therefore predict that the corresponding *B. subtilis* tRNA is not modified with Ψ at this same position, and this is confirmed by the PsiSeq analysis (Table S3).

### 3.6. System Analysis of tRNA Modification Genes in *B. subtilis*

As *B. subtilis* is such a prominent model organism, many whole genome datasets are available. These include mutant collections, phenotypic screens, TnSeq datasets, and gene expression data [24]. Many of the tRNA modification genes are poorly annotated and/or were only recently discovered, so we revisited several published genome-wide studies that used 42 *B. subtilis* tRNA modification genes as input.

Seven out of the 42 *B. subtilis* modification genes are essential (Figure 3) [38]. *FolE/mtrA* is essential because of its role in the synthesis of the tetrahydrofolate cofactor, but not its role in Q synthesis, as it can be replaced by *folE2/yciA* only under zinc limitation [69]. *TilS* is essential in most bacteria [10], with the rare exception of organisms such as *Mycoplasma mobile* that have changed the minor tRNA<sup>Ile</sup> anticodon from CAU to UAU, making the lysidine modification dispensable [97]. The essentiality of *trmD* in many bacteria has already made it an antibacterial target [98]. Interestingly, two t<sup>6</sup>A37 synthesis genes are essential (*tsaB* and *tsaD*), whereas two others (*tsaC* and *tsaE*) are not, while all four orthologs are essential in *E. coli* [7]. This suggests that the essentiality of *tsaBD* in *B. subtilis* could be due to another role in the cell besides t<sup>6</sup>A formation in the tRNA. Finally, the *iscS1* and *mnmA* involved in s<sup>2</sup>U34 formation are essential in *B. subtilis* but not in *E. coli* [99], for reasons that are not totally clear.

A phenotype screen testing several conditions (growth, C or N sources, and temperature sensitivity) of two ordered, barcoded, erythromycin-resistance- and kanamycin-resistance-marked single-gene deletion libraries is available [38]. Fitness scores are given for every gene in every condition tested. We extracted from these data the scores for the 35 predicted tRNA modification genes that are not essential (Table S4), and we visualized these data using clustergrammer, after transforming them into Z-scores. While transforming the fitness data into Z-scores may have exaggerated the differences between the ErmR and KanR datasets, it was noted in their original publication [40] that discrepancies between datasets could be caused by the downstream effects of the resistance cassettes that were used to prepare each library of mutants. Since we observed differences in the fitness data between the ErmR and KanR datasets for our list of genes, only the ErmR dataset was used to prepare Figure 5. Any phenotype observations from these datasets should be experimentally validated.

As shown in Figure 5, genes in same modification pathways, such as *tsaC/tsaE* (t<sup>6</sup>A37), *tgt/queA* (Q), *queDCEF* (Q precursor), and *cmnm<sup>5</sup>*U34 (*mnmE/G*), have similar or the same fitness profiles, suggesting a biological relevance that will need to be further studied.

We also looked at published *B. subtilis* Tn-Seq studies and found that three tRNA modification genes have been identified as being involved in sporulation (*mnmE*, *mnmG*, and *miaA*) [100], and *mnmG* was also important for germination [101].
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Fitness Data Similarity Matrix for predicted tRNA modification genes. The fitness scores for the ErmR mutant set are displayed as a similarity matrix, with the deleted gene shown on each of the axes of the matrix. The position of the tRNA modification that each gene mediates is denoted by the colored box next to each gene's name. Fitness scores with similar values are shown in red while dissimilar scores are shown in blue. The full dataset is given in Table S4.</figcaption></figure>

Using the list of locus tags as input, we analyzed the expression of 42 *B. subtilis* RNA modification genes in all transcriptomic data available in the Colombos database (41 tRNA modification and one involved in rRNA modification) [37] (Table S5). Some conditions lead to reductions or increases in most tRNA modification genes and are logically conditions where all translation genes would be expected to be affected, such as those that affect growth rate or transcription inhibition. To further explore the conditions where the response was specific to just a few genes, we focused on the conditions that led to differences in expression for just a subset of genes, as described in the Methods section. As shown in Figure 6, we found that twelve genes involved in the synthesis of Anticodon Stem Loop modifications (ψ31/32, Q34, ψ38/39, i⁶A37, t⁶A37, k²C34, s²U34) were up or down regulated in specific stress conditions, opening possible avenues for further exploration.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Fitness Data Similarity Matrix for predicted tRNA modification genes. The fitness scores for the ErmR mutant set are displayed as a similarity matrix, with the deleted gene shown on each of the axes of the matrix. The position of the tRNA modification that each gene mediates is denoted by the colored box next to each gene's name. Fitness scores with similar values are shown in red while dissimilar scores are shown in blue. The full dataset is given in Table S4.</figcaption></figure>

Using the list of locus tags as input, we analyzed the expression of 42 *B. subtilis* RNA modification genes in all transcriptomic data available in the Colombos database (41 tRNA modification and one involved in rRNA modification) [37] (Table S5). Some conditions lead to reductions or increases in most tRNA modification genes and are logically conditions where all translation genes would be expected to be affected, such as those that affect growth rate or transcription inhibition. To further explore the conditions where the response was specific to just a few genes, we focused on the conditions that led to differences in expression for just a subset of genes, as described in the Methods section. As shown in Figure 6, we found that twelve genes involved in the synthesis of Anticodon Stem Loop modifications (ψ31/32, Q34, ψ38/39, i⁶A37, t⁶A37, k²C34, s²U34) were up or down regulated in specific stress conditions, opening possible avenues for further exploration.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Most dramatic observed expression changes for tRNA modification genes. A subset of genes with dramatic shifts in expression levels across 47 contrasts were selected to create a separate heatmap. Genes with increased expression levels are shown in red and genes with decreased expression levels are shown in blue. The gene names are shown on the left axis and the contrasts are listed on the top axis. See Table S5 for the complete data and the lists of conditions.</figcaption></figure>

## 4. Conclusions

This study strengthens the efficiency of using an integrated strategy to identify tRNA modifications and the corresponding genes in a given organism, as recently illustrated with *V. cholerae* [19]. Mass spectrometry, tRNA-Seq, and comparative genomic methods do not give accurate or complete results on their own, but when they are used in combination, one can generate a high confidence tRNA modification map. Such integrated approaches involving both LC-MS analysis and the deep sequencing-based mapping of RNA modifications will be particularly useful for the analysis of still unexplored bacterial, archaeal, and eukaryotic organisms.

Even though we still have not been able to sequence all tRNA molecules of *B. subtilis*, and hence some modifications might remain to be discovered and/or mapped, we do believe that the set of tRNA modification genes in *B. subtilis* is nearly complete, with only two predicted genes missing (the functional homolog of *mnmC* and an alternate ho⁵U synthesis gene) (Figure 3 and Table S3). The model Gram-positive *B. subtilis* harbors slightly fewer modifications (2/3) than the model Gram-negative *E. coli* but contains twice the amount found in the minimal Gram-positive *M. capricolum* (Table 3).

<table><caption>**Table 3.** Comparison of numbers of tRNA modifications and associated genes in three model organisms.</caption><thead><tr><th>Organism</th><th>Modifications at Specific Positions in tRNAs</th><th>Modification Genes</th></tr></thead><tbody><tr><td><i>Bacillus subtilis</i> subsp. subtilis str. 168</td><td>35</td><td>41 + 2 missing</td></tr><tr><td><i>Mycoplasma capricolum</i></td><td>17</td><td>22</td></tr><tr><td><i>Escherichia coli</i> K12 MG1655</td><td>45</td><td>59</td></tr></tbody></table>

As shown in Figure S4, the levels of specific tRNA modifications can vary greatly between rare modifications such as k²C34 or abundant modifications such as i⁶A37. Further studies will be required to follow how modification levels change with growth conditions or stresses, but the transcriptomics and phenotype data do suggest that the levels of specific modifications are affected by
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Most dramatic observed expression changes for tRNA modification genes. A subset of genes with dramatic shifts in expression levels across 47 contrasts were selected to create a separate heatmap. Genes with increased expression levels are shown in red and genes with decreased expression levels are shown in blue. The gene names are shown on the left axis and the contrasts are listed on the top axis. See Table S5 for the complete data and the lists of conditions.</figcaption></figure>

## 4. Conclusions

This study strengthens the efficiency of using an integrated strategy to identify tRNA modifications and the corresponding genes in a given organism, as recently illustrated with *V. cholerae* [19]. Mass spectrometry, tRNA-Seq, and comparative genomic methods do not give accurate or complete results on their own, but when they are used in combination, one can generate a high confidence tRNA modification map. Such integrated approaches involving both LC-MS analysis and the deep sequencing-based mapping of RNA modifications will be particularly useful for the analysis of still unexplored bacterial, archaeal, and eukaryotic organisms.

Even though we still have not been able to sequence all tRNA molecules of *B. subtilis*, and hence some modifications might remain to be discovered and/or mapped, we do believe that the set of tRNA modification genes in *B. subtilis* is nearly complete, with only two predicted genes missing (the functional homolog of *mnmC* and an alternate *ho*⁵*U* synthesis gene) (Figure 3 and Table S3). The model Gram-positive *B. subtilis* harbors slightly fewer modifications (2/3) than the model Gram-negative *E. coli* but contains twice the amount found in the minimal Gram-positive *M. capricolum* (Table 3).

<table><caption>**Table 3.** Comparison of numbers of tRNA modifications and associated genes in three model organisms.</caption><thead><tr><th>Organism</th><th>Modifications at Specific Positions in tRNAs</th><th>Modification Genes</th></tr></thead><tbody><tr><td><i>Bacillus subtilis</i> subsp. subtilis str. 168</td><td>35</td><td>41 + 2 missing</td></tr><tr><td><i>Mycoplasma capricolum</i></td><td>17</td><td>22</td></tr><tr><td><i>Escherichia coli</i> K12 MG1655</td><td>45</td><td>59</td></tr></tbody></table>

As shown in Figure S4, the levels of specific tRNA modifications can vary greatly between rare modifications such as *k*²C34 or abundant modifications such as *i*⁶A37. Further studies will be required to follow how modification levels change with growth conditions or stresses, but the transcriptomics and phenotype data do suggest that the levels of specific modifications are affected by
growth conditions and that cellular processes that integrate tRNA modifications could be part of the regulatory loop in *B. subtilis*.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2218-273X/10/7/977/s1, Table S1: Validated *E. coli* tRNA modification genes, Table S2: Predicted and validated *B. subtilis* modification genes, Table S3: Summary of B. subtilis tRNA modification mapping by RiboMethSeq and PsiSeq, Table S4: Fitness data for tRNA modification genes in *B. subtilis*, Table S5: Expression data tRNA modification genes in *B. subtilis*, Figure S1: Mapping of tRNA and rRNA modification by RiboMethSeq and derived RT-signatures, Figure S2: Tandem mass spectrum of different tRNAs at different positions, Figure S3: Bar plots for normalized PsiScore values for six *B. subtilis* tRNAs containing ψ31 and ψ32 residues, Figure S4: Relative ion abundances of tRNA modifications found in *B. subtilis* 168.

**Author Contributions:** Conceptualization, V.d.C.-L. and R.L.R.; Data curation, V.d.C.-L. and Y.M.; Formal analysis, V.d.C.-L., R.L.R. and Y.M.; Funding acquisition, V.d.C.-L., Y.M. and P.A.L.; Investigation, R.L.R., M.J., V.M., C.E. and D.B.; Methodology, R.L.R.; Project administration, V.d.C.-L. and Y.M.; Supervision, V.d.C.-L. and P.A.L.; Visualization, R.L.R., M.J. and Y.M.; Writing—original draft, V.d.C.-L., R.L.R., M.J., D.B. and Y.M.; Writing—review & editing, V.d.C.-L., C.E., Y.M. and P.A.L. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the National Institutes of Health (R01 GM70641 to V dC-L), EpiARN FRCR project from Grand Est Region to Y.M. and the National Science Foundation (NSF CHE 150357 to P.A.L.).

**Acknowledgments:** We thank Djemel Handame, and Louis Droogmans for insightful discussions. We acknowledge the early work in tRNA modification mapping by LC-MS/MS conducted by Collin Wetzel. VdC-L would like to thank the Eduardo Rocha laboratory (Pasteur Institute) for having hosted her during her 2017 sabbatical when this analysis was initiated.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

**References**

1. El Yacoubi, B.; Bailly, M.; de Crécy-Lagard, V. Biosynthesis and function of posttranscriptional modifications of Transfer RNAs. *Annu. Rev. Genet.* **2012**, *46*, 69–95. [CrossRef] [PubMed]

2. Agris, P.F.; Eruysal, E.R.; Narendran, A.; Väre, V.Y.P.; Vangaveti, S.; Ranganathan, S.V. Celebrating wobble decoding: Half a century and still much is new. *RNA Biol.* **2018**, *15*, 537–553. [CrossRef]

3. Kimura, S.; Waldor, M.K. The RNA degradosome promotes tRNA quality control through clearance of hypomodified tRNA. *Proc. Natl. Acad. Sci. USA* **2019**, *116*, 1394–1403. [CrossRef] [PubMed]

4. Alexandrov, A.; Chernyakov, I.; Gu, W.; Hiley, S.L.; Hughes, T.R.; Grayhack, E.J.; Phizicky, E.M. Rapid tRNA decay can result from lack of nonessential modifications. *Mol. Cell* **2006**, *21*, 87–96. [CrossRef] [PubMed]

5. Koh, C.S.; Sarin, L.P. Transfer RNA modification and infection-Implications for pathogenicity and host responses. *Biochim. Biophys. Acta-Gene Regul. Mech.* **2018**, *1861*, 419–432. [CrossRef] [PubMed]

6. Lorenz, C.; Lünse, C.; Mörl, M. tRNA modifications: Impact on structure and thermal adaptation. *Biomolecules* **2017**, *7*, 35. [CrossRef]

7. Thiaville, P.C.; El Yacoubi, B.; Köhrer, C.; Thiaville, J.J.; Deutsch, C.; Iwata-Reuyl, D.; Bacusmo, J.M.; Armengaud, J.; Bessho, Y.; Wetzel, C.; et al. Essentiality of threonylcarbamoyladenosine (t6A), a universal tRNA modification, in bacteria. *Mol. Microbiol.* **2015**, *98*, 1199–1221. [CrossRef]

8. Soma, A.; Ikeuchi, Y.; Kanemasa, S.; Kobayashi, K.; Ogasawara, N.; Ote, T.; Kato, J.; Watanabe, K.; Sekine, Y.; Suzuki, T.; et al. An RNA-modifying enzyme that governs both the codon and amino acid specificities of isoleucine tRNA. *Mol. Cell.* **2003**, *12*, 689–698. [CrossRef]

9. Ikeuchi, Y.; Soma, A.; Ote, T.; Kato, J.; Sekine, Y.; Suzuki, T. Molecular mechanism of lysidine synthesis that determines tRNA identity and codon recognition. *Mol. Cell* **2005**, *19*, 235–246. [CrossRef]

10. Fabret, C.; Dervyn, E.; Dalmais, B.; Guillot, A.; Marck, C.; Grosjean, H.; Noirot, P. Life without the essential bacterial tRNA Ile2-lysidine synthetase TisS: A case of tRNA gene recruitment in *Bacillus subtilis*. *Mol. Microbiol.* **2011**, *80*, 1062–1074. [CrossRef]

11. Björk, G.R.; Hagervall, T.G. Transfer RNA modification: Presence, synthesis, and function. *EcoSal Plus* **2014**, *6*. [CrossRef] [PubMed]

12. Persson, B.C. Modification of tRNA as a regulatory device. *Mol. Microbiol.* **1993**, *8*, 1011–1016. [CrossRef] [PubMed]
growth conditions and that cellular processes that integrate tRNA modifications could be part of the regulatory loop in *B. subtilis*.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2218-273X/10/7/977/s1, Table S1: Validated *E. coli* tRNA modification genes, Table S2: Predicted and validated *B. subtilis* modification genes, Table S3: Summary of B. subtilis tRNA modification mapping by RiboMethSeq and PsiSeq, Table S4: Fitness data for tRNA modification genes in *B. subtilis*, Table S5: Expression data tRNA modification genes in *B. subtilis*, Figure S1: Mapping of tRNA and rRNA modification by RiboMethSeq and derived RT-signatures, Figure S2: Tandem mass spectrum of different tRNAs at different positions, Figure S3: Bar plots for normalized PsiScore values for six *B. subtilis* tRNAs containing ψ31 and ψ32 residues, Figure S4: Relative ion abundances of tRNA modifications found in *B. subtilis* 168.

**Author Contributions:** Conceptualization, V.d.C.-L. and R.L.R.; Data curation, V.d.C.-L. and Y.M.; Formal analysis, V.d.C.-L., R.L.R. and Y.M.; Funding acquisition, V.d.C.-L., Y.M. and P.A.L.; Investigation, R.L.R., M.J., V.M., C.E. and D.B.; Methodology, R.L.R.; Project administration, V.d.C.-L. and Y.M.; Supervision, V.d.C.-L. and P.A.L.; Visualization, R.L.R., M.J. and Y.M.; Writing—original draft, V.d.C.-L., R.L.R., M.J., D.B. and Y.M.; Writing—review & editing, V.d.C.-L., C.E., Y.M. and P.A.L. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the National Institutes of Health (R01 GM70641 to V dC-L), EpiARN FRCR project from Grand Est Region to Y.M. and the National Science Foundation (NSF CHE 150357 to P.A.L.).

**Acknowledgments:** We thank Djemel Handame, and Louis Droogmans for insightful discussions. We acknowledge the early work in tRNA modification mapping by LC-MS/MS conducted by Collin Wetzel. VdC-L would like to thank the Eduardo Rocha laboratory (Pasteur Institute) for having hosted her during her 2017 sabbatical when this analysis was initiated.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. El Yacoubi, B.; Bailly, M.; de Crécy-Lagard, V. Biosynthesis and function of posttranscriptional modifications of Transfer RNAs. *Annu. Rev. Genet.* **2012**, *46*, 69–95. [CrossRef] [PubMed]

2. Agris, P.F.; Eruysal, E.R.; Narendran, A.; Väre, V.Y.P.; Vangaveti, S.; Ranganathan, S.V. Celebrating wobble decoding: Half a century and still much is new. *RNA Biol.* **2018**, *15*, 537–553. [CrossRef]

3. Kimura, S.; Waldor, M.K. The RNA degradosome promotes tRNA quality control through clearance of hypomodified tRNA. *Proc. Natl. Acad. Sci. USA* **2019**, *116*, 1394–1403. [CrossRef] [PubMed]

4. Alexandrov, A.; Chernyakov, I.; Gu, W.; Hiley, S.L.; Hughes, T.R.; Grayhack, E.J.; Phizicky, E.M. Rapid tRNA decay can result from lack of nonessential modifications. *Mol. Cell* **2006**, *21*, 87–96. [CrossRef] [PubMed]

5. Koh, C.S.; Sarin, L.P. Transfer RNA modification and infection-Implications for pathogenicity and host responses. *Biochim. Biophys. Acta-Gene Regul. Mech.* **2018**, *1861*, 419–432. [CrossRef] [PubMed]

6. Lorenz, C.; Lünse, C.; Mörl, M. tRNA modifications: Impact on structure and thermal adaptation. *Biomolecules* **2017**, *7*, 35. [CrossRef]

7. Thiaville, P.C.; El Yacoubi, B.; Köhrer, C.; Thiaville, J.J.; Deutsch, C.; Iwata-Reuyl, D.; Bacusmo, J.M.; Armengaud, J.; Bessho, Y.; Wetzel, C.; et al. Essentiality of threonylcarbamoyladenosine (t6A), a universal tRNA modification, in bacteria. *Mol. Microbiol.* **2015**, *98*, 1199–1221. [CrossRef]

8. Soma, A.; Ikeuchi, Y.; Kanemasa, S.; Kobayashi, K.; Ogasawara, N.; Ote, T.; Kato, J.; Watanabe, K.; Sekine, Y.; Suzuki, T.; et al. An RNA-modifying enzyme that governs both the codon and amino acid specificities of isoleucine tRNA. *Mol. Cell.* **2003**, *12*, 689–698. [CrossRef]

9. Ikeuchi, Y.; Soma, A.; Ote, T.; Kato, J.; Sekine, Y.; Suzuki, T. Molecular mechanism of lysidine synthesis that determines tRNA identity and codon recognition. *Mol. Cell* **2005**, *19*, 235–246. [CrossRef]

10. Fabret, C.; Dervyn, E.; Dalmais, B.; Guillot, A.; Marck, C.; Grosjean, H.; Noirot, P. Life without the essential bacterial tRNA Ile2-lysidine synthetase TisS: A case of tRNA gene recruitment in *Bacillus subtilis*. *Mol. Microbiol.* **2011**, *80*, 1062–1074. [CrossRef]

11. Björk, G.R.; Hagervall, T.G. Transfer RNA modification: Presence, synthesis, and function. *EcoSal Plus* **2014**, *6*. [CrossRef] [PubMed]

12. Persson, B.C. Modification of tRNA as a regulatory device. *Mol. Microbiol.* **1993**, *8*, 1011–1016. [CrossRef] [PubMed]
13. Gu, C.; Begley, T.J.; Dedon, P.C. tRNA modifications regulate translation during cellular stress. *FEBS Lett.* **2014**, *588*, 4287-4296. [CrossRef] [PubMed]

14. Chan, C.; Pham, P.; Dedon, P.C.; Begley, T.J. Lifestyle modifications: Coordinating the tRNA epitranscriptome with codon bias to adapt translation during stress responses. *Genome Biol.* **2018**, *19*, 228. [CrossRef] [PubMed]

15. Pollo-Oliveira, L.; de Crécy-Lagard, V. Can protein expression be regulated by modulation of tRNA Modification profiles? *Biochemistry* **2019**, *58*, 355-362. [CrossRef]

16. Huber, S.M.; Leonardi, A.; Dedon, P.C.; Begley, T.J. The versatile roles of the tRNA epitranscriptome during cellular responses to toxic exposures and environmental stress. *Toxics* **2019**, *7*, 17. [CrossRef]

17. Raina, M.; Ibba, M. tRNAs as regulators of biological processes. *Front. Genet.* **2014**, *5*, 171. [CrossRef]

18. De Crécy-Lagard, V.; Marck, C.; Brochier-Armanet, C.; Grosjean, H. Comparative RNomics and Modomics in Mollicutes: Prediction of gene function and evolutionary implications. *IUBMB Life* **2007**, *59*, 634-658. [CrossRef]

19. Kimura, S.; Dedon, P.C.; Waldor, M.K. Comparative tRNA sequencing and RNA mass spectrometry for surveying tRNA modifications. *Nat. Chem. Biol.* **2020**, in press. [CrossRef]

20. Puri, P.; Wetzel, C.; Saffert, P.; Gaston, K.W.; Russell, S.P.; Cordero Varela, J.A.; van der Vlies, P.; Zhang, G.; Limbach, P.A.; Ignatova, Z.; et al. Systematic identification of tRNAome and its dynamics in *Lactococcus lactis*. *Mol. Microbiol.* **2014**, *93*, 944-956. [CrossRef]

21. Cao, X.; Limbach, P.A. Enhanced detection of post-transcriptional modifications using a mass-exclusion list strategy for RNA modification mapping by LC-MS/MS. *Anal. Chem.* **2015**, *87*, 8433-8440. [CrossRef] [PubMed]

22. Chionh, Y.H.; McBee, M.; Babu, I.R.; Hia, F.; Lin, W.; Zhao, W.; Cao, J.; Dziergowska, A.; Malkiewicz, A.; Begley, T.J.; et al. tRNA-mediated codon-biased translation in mycobacterial hypoxic persistence. *Nat. Commun.* **2016**, *7*, 13302. [CrossRef] [PubMed]

23. Kunst, F.; Ogasawara, N.; Moszer, I.; Albertini, A.M.; Alloni, G.; Azevedo, V.; Bertero, M.G.; Bessieres, P.; Bolotin, A.; Borchert, S.; et al. The complete genome sequence of the gram-positive bacterium *Bacillus subtilis*. *Nature* **1997**, *390*, 249-256. [CrossRef]

24. Flórez, L.A.; Roppel, S.F.; Schmeisky, A.G.; Lammers, C.R.; Stülke, J. A community-curated consensual annotation that is continuously updated: The *Bacillus subtilis* centered wiki SubtiWiki. *Database* **2009**, *2009*, bap012. [CrossRef] [PubMed]

25. Belda, E.; Sekowska, A.; Le Fèvre, F.; Morgat, A.; Mornico, D.; Ouzounis, C.; Vallenet, D.; Médigue, C.; Danchin, A. An updated metabolic view of the *Bacillus subtilis* 168 genome. *Microbiology* **2013**, *159*, 757-770. [CrossRef] [PubMed]

26. Chan, P.P.; Lowe, T.M. GtRNAdb 2.0: An expanded database of transfer RNA genes identified in complete and draft genomes. *Nucleic Acids Res.* **2016**, *44*, 184-189. [CrossRef]

27. Boccaletto, P.; MacHnicka, M.A.; Purta, E.; Pitkowski, P.; Baginski, B.; Wirecki, T.K.; de Crécy-Lagard, V.; Ross, R.; Limbach, P.A.; Kotter, A.; et al. MODOMICS: A database of RNA modification pathways. 2017 update. *Nucleic Acids Res.* **2018**, *46*, 303-307. [CrossRef]

28. Jühling, F.; Mörl, M.; Hartmann, R.K.; Sprinzl, M.; Stadler, P.F.; Pütz, J. tRNAdb 2009: Compilation of tRNA sequences and tRNA genes. *Nucleic Acids Res.* **2009**, *37*, 159-162. [CrossRef]

29. Sajek, M.P.; Woźniak, T.; Sprinzl, M.; Jaruzelska, J.; Barciszewski, J. T-psi-C: User friendly database of tRNA sequences and structures. *Nucleic Acids Res.* **2020**, *48*, 256-260. [CrossRef]

30. The UniProt Consortium. UniProt: The universal protein knowledgebase. *Nucleic Acids Res.* **2017**, *45*, 158-169. [CrossRef]

31. Altschul, S.F.; Madden, T.L.; Schaffer, A.A.; Zhang, J.H.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. *FASEB J.* **1997**, *25*, 1326. [CrossRef] [PubMed]

32. NCBI Resource Coordinators. Database resources of the National Center for Biotechnology Information. *Nucleic Acids Res.* **2018**, *46*, 8-13. [CrossRef] [PubMed]

33. Wattam, A.R.; Davis, J.J.; Assaf, R.; Boisvert, S.; Brettin, T.; Bun, C.; Conrad, N.; Dietrich, E.M.; Disz, T.; Gabbard, J.L.; et al. Improvements to PATRIC, the all-bacterial bioinformatics database and analysis resource center. *Nucleic Acids Res.* **2017**, *45*, 535-542. [CrossRef] [PubMed]
13. Gu, C.; Begley, T.J.; Dedon, P.C. tRNA modifications regulate translation during cellular stress. *FEBS Lett.* **2014**, *588*, 4287–4296. [CrossRef] [PubMed]

14. Chan, C.; Pham, P.; Dedon, P.C.; Begley, T.J. Lifestyle modifications: Coordinating the tRNA epitranscriptome with codon bias to adapt translation during stress responses. *Genome Biol.* **2018**, *19*, 228. [CrossRef] [PubMed]

15. Pollo-Oliveira, L.; de Crécy-Lagard, V. Can protein expression be regulated by modulation of tRNA Modification profiles? *Biochemistry* **2019**, *58*, 355–362. [CrossRef]

16. Huber, S.M.; Leonardi, A.; Dedon, P.C.; Begley, T.J. The versatile roles of the tRNA epitranscriptome during cellular responses to toxic exposures and environmental stress. *Toxics* **2019**, *7*, 17. [CrossRef]

17. Raina, M.; Ibba, M. tRNAs as regulators of biological processes. *Front. Genet.* **2014**, *5*, 171. [CrossRef]

18. De Crécy-Lagard, V.; Marck, C.; Brochier-Armanet, C.; Grosjean, H. Comparative RNomics and Modomics in Mollicutes: Prediction of gene function and evolutionary implications. *IUBMB Life* **2007**, *59*, 634–658. [CrossRef]

19. Kimura, S.; Dedon, P.C.; Waldor, M.K. Comparative tRNA sequencing and RNA mass spectrometry for surveying tRNA modifications. *Nat. Chem. Biol.* **2020**, in press. [CrossRef]

20. Puri, P.; Wetzel, C.; Saffert, P.; Gaston, K.W.; Russell, S.P.; Cordero Varela, J.A.; van der Vlies, P.; Zhang, G.; Limbach, P.A.; Ignatova, Z.; et al. Systematic identification of tRNAome and its dynamics in *Lactococcus lactis*. *Mol. Microbiol.* **2014**, *93*, 944–956. [CrossRef]

21. Cao, X.; Limbach, P.A. Enhanced detection of post-transcriptional modifications using a mass-exclusion list strategy for RNA modification mapping by LC-MS/MS. *Anal. Chem.* **2015**, *87*, 8433–8440. [CrossRef] [PubMed]

22. Chionh, Y.H.; McBee, M.; Babu, I.R.; Hia, F.; Lin, W.; Zhao, W.; Cao, J.; Dziergowska, A.; Malkiewicz, A.; Begley, T.J.; et al. tRNA-mediated codon-biased translation in mycobacterial hypoxic persistence. *Nat. Commun.* **2016**, *7*, 13302. [CrossRef] [PubMed]

23. Kunst, F.; Ogasawara, N.; Moszer, I.; Albertini, A.M.; Alloni, G.; Azevedo, V.; Bertero, M.G.; Bessieres, P.; Bolotin, A.; Borchert, S.; et al. The complete genome sequence of the gram-positive bacterium *Bacillus subtilis*. *Nature* **1997**, *390*, 249–256. [CrossRef]

24. Flórez, L.A.; Roppel, S.F.; Schmeisky, A.G.; Lammers, C.R.; Stülke, J. A community-curated consensual annotation that is continuously updated: The *Bacillus subtilis* centered wiki SubtiWiki. *Database* **2009**, *2009*, bap012. [CrossRef] [PubMed]

25. Belda, E.; Sekowska, A.; Le Fèvre, F.; Morgat, A.; Mornico, D.; Ouzounis, C.; Vallenet, D.; Médigue, C.; Danchin, A. An updated metabolic view of the *Bacillus subtilis* 168 genome. *Microbiology* **2013**, *159*, 757–770. [CrossRef] [PubMed]

26. Chan, P.P.; Lowe, T.M. GtRNAdb 2.0: An expanded database of transfer RNA genes identified in complete and draft genomes. *Nucleic Acids Res.* **2016**, *44*, 184–189. [CrossRef]

27. Boccaletto, P.; MacHnicka, M.A.; Purta, E.; Pitkowski, P.; Baginski, B.; Wirecki, T.K.; de Crécy-Lagard, V.; Ross, R.; Limbach, P.A.; Kotter, A.; et al. MODOMICS: A database of RNA modification pathways. 2017 update. *Nucleic Acids Res.* **2018**, *46*, 303–307. [CrossRef]

28. Jühling, F.; Mörl, M.; Hartmann, R.K.; Sprinzl, M.; Stadler, P.F.; Pütz, J. tRNAdb 2009: Compilation of tRNA sequences and tRNA genes. *Nucleic Acids Res.* **2009**, *37*, 159–162. [CrossRef]

29. Sajek, M.P.; Woźniak, T.; Sprinzl, M.; Jaruzelska, J.; Barciszewski, J. T-psi-C: User friendly database of tRNA sequences and structures. *Nucleic Acids Res.* **2020**, *48*, 256–260. [CrossRef]

30. The UniProt Consortium. UniProt: The universal protein knowledgebase. *Nucleic Acids Res.* **2017**, *45*, 158–169. [CrossRef]

31. Altschul, S.F.; Madden, T.L.; Schaffer, A.A.; Zhang, J.H.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. *FASEB J.* **1997**, *25*, 1326. [CrossRef] [PubMed]

32. NCBI Resource Coordinators. Database resources of the National Center for Biotechnology Information. *Nucleic Acids Res.* **2018**, *46*, 8–13. [CrossRef] [PubMed]

33. Wattam, A.R.; Davis, J.J.; Assaf, R.; Boisvert, S.; Brettin, T.; Bun, C.; Conrad, N.; Dietrich, E.M.; Disz, T.; Gabbard, J.L.; et al. Improvements to PATRIC, the all-bacterial bioinformatics database and analysis resource center. *Nucleic Acids Res.* **2017**, *45*, 535–542. [CrossRef] [PubMed]
34. Overbeek, R.; Olson, R.; Pusch, G.D.; Olsen, G.J.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Parrello, B.; Shukla, M.; et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). *Nucleic Acids Res.* **2014**, *42*, 206–214. [CrossRef] [PubMed]

35. Robinson, J.T.; Thorvaldsdóttir, H.; Winckler, W.; Guttman, M.; Lander, E.S.; Getz, G.; Mesirov, J.P. Integrative genomics viewer. *Nat. Biotechnol.* **2011**, *29*, 24–26. [CrossRef]

36. Fernandez, N.F.; Gundersen, G.W.; Rahman, A.; Grimes, M.L.; Rikova, K.; Hornbeck, P.; Ma'ayan, A. Clustergrammer, a web-based heatmap visualization and analysis tool for high-dimensional biological data. *Sci. Data* **2017**, *4*, 170151. [CrossRef]

37. Moretto, M.; Sonego, P.; Dierckxsens, N.; Brilli, M.; Bianco, L.; Ledezma-Tejeida, D.; Gama-Castro, S.; Galardini, M.; Romualdi, C.; Laukens, K.; et al. COLOMBOS v3.0: Leveraging gene expression compendia for cross-species analyses. *Nucleic Acids Res.* **2016**, *44*, 620–623. [CrossRef]

38. Koo, B.-M.; Kritikos, G.; Farelli, J.D.; Todor, H.; Tong, K.; Kimsey, H.; Wapinski, I.; Galardini, M.; Cabal, A.; Peters, J.M.; et al. Construction and analysis of two genome-scale deletion libraries for *Bacillus subtilis*. *Cell Syst.* **2017**, *4*, 291–305. [CrossRef]

39. Ross, R.; Cao, X.; Yu, N.; Limbach, P.A. Sequence mapping of transfer RNA chemical modifications by liquid chromatography tandem mass spectrometry. *Methods* **2016**, *107*, 73–78. [CrossRef]

40. Yu, N.; Lobue, P.A.; Cao, X.; Limbach, P.A. RNAModMapper: RNA modification mapping software for analysis of liquid chromatography tandem mass spectrometry data. *Anal. Chem.* **2017**, *89*, 10744–10752. [CrossRef]

41. Marchand, V.; Blanloeil-Oillo, F.; Helm, M.; Motorin, Y. Illumina-based RiboMethSeq approach for mapping of 22032-O-Me residues in RNA. *Nucleic Acids Res.* **2016**, *44*, e135. [CrossRef]

42. Marchand, V.; Pichot, F.; Thüring, K.; Ayadi, L.; Freund, I.; Dalpke, A.; Helm, M.; Motorin, Y. Next-generation sequencing-based ribomethseq protocol for analysis of tRNA 2'-O-methylation. *Biomolecules* **2017**, *7*, 13. [CrossRef]

43. Werner, S.; Schmidt, L.; Marchand, V.; Kemmer, T.; Falschlunger, C.; Sednev, M.V.; Bec, G.; Ennifar, E.; Höbartner, C.; Micura, R.; et al. Machine learning of reverse transcription signatures of variegated polymerases allows mapping and discrimination of methylated purines in limited transcriptomes. *Nucleic Acids Res.* **2020**, *48*, 3734–3746. [CrossRef] [PubMed]

44. Motorin, Y.; Marchand, V. Detection and analysis of RNA ribose 2'-O-methylations: Challenges and solutions. *Genes* **2018**, *9*, 642. [CrossRef] [PubMed]

45. Motorin, Y.; Muller, S.; Behm-Ansmant, I.; Branlant, C. Identification of modified residues in RNAs by reverse transcription-based methods. *Methods Enzymol.* **2007**, *425*, 21–53. [PubMed]

46. Behm-Ansmant, I.; Helm, M.; Motorin, Y. Use of specific chemical reagents for detection of modified nucleotides in RNA. *J. Nucleic Acids* **2011**, *2011*, 408053. [CrossRef] [PubMed]

47. De Brouwer, A.P.M.; Abou Jamra, R.; Körtel, N.; Soyris, C.; Polla, D.L.; Safra, M.; Zisso, A.; Powell, C.A.; Rebelo-Guiomar, P.; Dinges, N.; et al. Variants in PUS7 cause intellectual disability with speech delay, microcephaly, short stature, and aggressive behavior. *Am. J. Hum. Genet.* **2018**, *103*, 1045–1052. [CrossRef] [PubMed]

48. Grosjean, H.; Westhof, E. An integrated, structure- and energy-based view of the genetic code. *Nucleic Acids Res.* **2016**, *44*, 8020–8040. [CrossRef]

49. Moukadiri, I.; Villarroya, M.; Benítez-Páez, A.; Armengod, M.-E. *Bacillus subtilis* exhibits MnmC-like tRNA modification activities. *RNA Biol.* **2018**, *15*, 1167–1173. [CrossRef]

50. Taniguchi, T.; Miyauchi, K.; Sakaguchi, Y.; Yamashita, S.; Soma, A.; Tomita, K.; Suzuki, T. Acetate-dependent tRNA acetylation required for decoding fidelity in protein synthesis. *Nat. Chem. Biol.* **2018**, *14*, 1010–1020. [CrossRef]

51. Miles, Z.D.; McCarty, R.M.; Molnar, G.; Bandarian, V. Discovery of epoxyqueuosine (oQ) reductase reveals parallels between halorespiration and tRNA modification. *Proc. Natl. Acad. Sci. USA* **2011**, *108*, 7368–7372. [CrossRef] [PubMed]

52. Moukadiri, I.; Prado, S.; Piera, J.; Velázquez-Campoy, A.; Björk, G.R.; Armengod, M.E. Evolutionarily conserved proteins MnmE and GidA catalyze the formation of two methyluridine derivatives at tRNA wobble positions. *Nucleic Acids Res.* **2009**, *37*, 7177–7193. [CrossRef] [PubMed]

53. Miyauchi, K.; Kimura, S.; Suzuki, T. A cyclic form of N6-threonylcarbamoyladenosine as a widely distributed tRNA hypermodification. *Nat. Chem. Biol.* **2013**, *9*, 105–111. [CrossRef] [PubMed]
34. Overbeek, R.; Olson, R.; Pusch, G.D.; Olsen, G.J.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Parrello, B.; Shukla, M.; et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). *Nucleic Acids Res.* **2014**, *42*, 206–214. [CrossRef] [PubMed]

35. Robinson, J.T.; Thorvaldsdóttir, H.; Winckler, W.; Guttman, M.; Lander, E.S.; Getz, G.; Mesirov, J.P. Integrative genomics viewer. *Nat. Biotechnol.* **2011**, *29*, 24–26. [CrossRef]

36. Fernandez, N.F.; Gundersen, G.W.; Rahman, A.; Grimes, M.L.; Rikova, K.; Hornbeck, P.; Ma'ayan, A. Clustergrammer, a web-based heatmap visualization and analysis tool for high-dimensional biological data. *Sci. Data* **2017**, *4*, 170151. [CrossRef]

37. Moretto, M.; Sonego, P.; Dierckxsens, N.; Brilli, M.; Bianco, L.; Ledezma-Tejeida, D.; Gama-Castro, S.; Galardini, M.; Romualdi, C.; Laukens, K.; et al. COLOMBOS v3.0: Leveraging gene expression compendia for cross-species analyses. *Nucleic Acids Res.* **2016**, *44*, 620–623. [CrossRef]

38. Koo, B.-M.; Kritikos, G.; Farelli, J.D.; Todor, H.; Tong, K.; Kimsey, H.; Wapinski, I.; Galardini, M.; Cabal, A.; Peters, J.M.; et al. Construction and analysis of two genome-scale deletion libraries for *Bacillus subtilis*. *Cell Syst.* **2017**, *4*, 291–305. [CrossRef]

39. Ross, R.; Cao, X.; Yu, N.; Limbach, P.A. Sequence mapping of transfer RNA chemical modifications by liquid chromatography tandem mass spectrometry. *Methods* **2016**, *107*, 73–78. [CrossRef]

40. Yu, N.; Lobue, P.A.; Cao, X.; Limbach, P.A. RNAModMapper: RNA modification mapping software for analysis of liquid chromatography tandem mass spectrometry data. *Anal. Chem.* **2017**, *89*, 10744–10752. [CrossRef]

41. Marchand, V.; Blanloeil-Oillo, F.; Helm, M.; Motorin, Y. Illumina-based RibOMethSeq approach for mapping of 22032-O-Me residues in RNA. *Nucleic Acids Res.* **2016**, *44*, e135. [CrossRef]

42. Marchand, V.; Pichot, F.; Thüring, K.; Ayadi, L.; Freund, I.; Dalpke, A.; Helm, M.; Motorin, Y. Next-generation sequencing-based ribomethseq protocol for analysis of tRNA 2'-O-methylation. *Biomolecules* **2017**, *7*, 13. [CrossRef]

43. Werner, S.; Schmidt, L.; Marchand, V.; Kemmer, T.; Falschlunger, C.; Sednev, M.V.; Bec, G.; Ennifar, E.; Höbartner, C.; Micura, R.; et al. Machine learning of reverse transcription signatures of variegated polymerases allows mapping and discrimination of methylated purines in limited transcriptomes. *Nucleic Acids Res.* **2020**, *48*, 3734–3746. [CrossRef] [PubMed]

44. Motorin, Y.; Marchand, V. Detection and analysis of RNA ribose 2'-O-methylations: Challenges and solutions. *Genes* **2018**, *9*, 642. [CrossRef] [PubMed]

45. Motorin, Y.; Muller, S.; Behm-Ansmant, I.; Branlant, C. Identification of modified residues in RNAs by reverse transcription-based methods. *Methods Enzymol.* **2007**, *425*, 21–53. [PubMed]

46. Behm-Ansmant, I.; Helm, M.; Motorin, Y. Use of specific chemical reagents for detection of modified nucleotides in RNA. *J. Nucleic Acids* **2011**, *2011*, 408053. [CrossRef] [PubMed]

47. De Brouwer, A.P.M.; Abou Jamra, R.; Körtel, N.; Soyris, C.; Polla, D.L.; Safra, M.; Zisso, A.; Powell, C.A.; Rebelo-Guiomar, P.; Dinges, N.; et al. Variants in PUS7 cause intellectual disability with speech delay, microcephaly, short stature, and aggressive behavior. *Am. J. Hum. Genet.* **2018**, *103*, 1045–1052. [CrossRef] [PubMed]

48. Grosjean, H.; Westhof, E. An integrated, structure- and energy-based view of the genetic code. *Nucleic Acids Res.* **2016**, *44*, 8020–8040. [CrossRef]

49. Moukadiri, I.; Villarroya, M.; Benítez-Páez, A.; Armengod, M.-E. *Bacillus subtilis* exhibits MnmC-like tRNA modification activities. *RNA Biol.* **2018**, *15*, 1167–1173. [CrossRef]

50. Taniguchi, T.; Miyauchi, K.; Sakaguchi, Y.; Yamashita, S.; Soma, A.; Tomita, K.; Suzuki, T. Acetate-dependent tRNA acetylation required for decoding fidelity in protein synthesis. *Nat. Chem. Biol.* **2018**, *14*, 1010–1020. [CrossRef]

51. Miles, Z.D.; McCarty, R.M.; Molnar, G.; Bandarian, V. Discovery of epoxyqueuosine (oQ) reductase reveals parallels between halorespiration and tRNA modification. *Proc. Natl. Acad. Sci. USA* **2011**, *108*, 7368–7372. [CrossRef] [PubMed]

52. Moukadiri, I.; Prado, S.; Piera, J.; Velázquez-Campoy, A.; Björk, G.R.; Armengod, M.E. Evolutionarily conserved proteins MnmE and GidA catalyze the formation of two methyluridine derivatives at tRNA wobble positions. *Nucleic Acids Res.* **2009**, *37*, 7177–7193. [CrossRef] [PubMed]

53. Miyauchi, K.; Kimura, S.; Suzuki, T. A cyclic form of N6-threonylcarbamoyladenosine as a widely distributed tRNA hypermodification. *Nat. Chem. Biol.* **2013**, *9*, 105–111. [CrossRef] [PubMed]
54. Kimura, S.; Suzuki, T. Fine-tuning of the ribosomal decoding center by conserved methyl-modifications in the *Escherichia coli* 16S rRNA. *Nucleic Acids Res.* **2010**, *38*, 1341–1352. [CrossRef] [PubMed]

55. Emmerechts, G.; Barbe, S.; Herdewijn, P.; Anne, J.; Rozenski, J. Post-transcriptional modification mapping in the *Clostridium acetobutylicum* 16S rRNA by mass spectrometry and reverse transcriptase assays. *Nucleic Acids Res.* **2007**, *35*, 3494–3503. [CrossRef]

56. Desai, P.M.; Culver, G.M.; Rife, J.P. Site-directed mutants of 16S rRNA reveal important RNA domains for KsgA function and 30S subunit assembly. *Biochemistry* **2011**, *50*, 854–863. [CrossRef]

57. Poldermans, B.; Bakker, H.; Van Knippenberg, P.H. Studies on the function of two adjacent N6, N6-dimethyladenosines near the 3' end of 16 S ribosomal RNA of *Escherichia coli*. IV. The effect of the methylgroups on ribosomal subunit interaction. *Nucleic Acids Res.* **1980**, *8*, 143–152. [CrossRef]

58. Grosjean, H.; Edqvist, J.; Straby, K.B.; Gieg'e, R. Enzymatic formation of modified nucleosides in tRNA: Dependence on tRNA architecture. *J. Mol. Biol.* **1996**, *255*, 67–85. [CrossRef]

59. Armengaud, J.; Urbonavicius, J.; Fernandez, B.; Chaussinand, G.; Bujnicki, J.M.; Grosjean, H. N2-methylation of guanosine at position 10 in tRNA is catalyzed by a THUMP domain-containing, S-Adenosylmethionine-dependent methyltransferase, conserved in Archaea and Eukaryota. *J. Biol. Chem.* **2004**, *279*, 37142–37152. [CrossRef]

60. Pallan, P.S.; Kreutz, C.; Bosio, S.; Micura, R.; Egli, M. Effects of N(2),N(2)-dimethylguanosine on RNA structure and stability: Crystal structure of an RNA duplex with tandem m(2) (2)G:A pairs. *RNA* **2008**, *14*, 2452–2453. [CrossRef]

61. Golovina, A.Y.; Dzama, M.M.; Osterman, I.A.; Sergiev, P.V.; Serebryakova, M.V.; Bogdanov, A.A.; Dontsova, O.A. The last rRNA methyltransferase of *E. coli* revealed: The *yhiR* gene encodes adenine-N6 methyltransferase specific for modification of A2030 of 23S ribosomal RNA. *RNA* **2012**, *18*, 1725–1734. [CrossRef] [PubMed]

62. Grosjean, H.; Constantinesco, F.; Foiret, D.; Benachenhou, N. A novel enzymatic pathway leading to 1-methylinosine modification in *Haloferax volcanii* tRNA. *Nucleic Acids Res.* **1995**, *23*, 4312–4319. [CrossRef] [PubMed]

63. Jora, M.; Lobue, P.A.; Ross, R.L.; Williams, B.; Addepalli, B. Detection of ribonucleoside modifications by liquid chromatography coupled with mass spectrometry. *Biochim. Biophys. Acta-Gene Regul. Mech.* **2019**, *1862*, 280–290. [CrossRef]

64. Singh, R.; Farmer, P.B. Liquid chromatography-electrospray ionization-mass spectrometry: The future of DNA adduct detection. *Carcinogenesis* **2005**, *27*, 178–196. [CrossRef]

65. Arnold, H.H.; Keith, G. The nucleotide sequence of phenylalanine tRNA from *Bacillus subtilis*. *Nucleic Acids Res.* **1977**, *4*, 2821–2829. [CrossRef] [PubMed]

66. Matsumoto, T.; Ohta, T.; Kumagai, I.; Oshima, T.; Murao, K.; Hasegawa, T.; Ishikura, H.; Watanabe, K. A thermostable Gm-methylase recognizes the tertiary structure of tRNA 1. *J. Biochem.* **1987**, *101*, 1191–1198. [CrossRef]

67. Vold, B.S.; Keith, D.E., Jr.; Buck, M.; McCloskey, J.A.; Pang, H. Lysine tRNAs from *Bacillus subtilis* 168: Structural analysis. *Nucleic Acids Res.* **1982**, *10*, 3125–3132. [CrossRef]

68. Kang, B.-I.; Miyauchi, K.; Matuszewski, M.; D'Almeida, G.S.; Rubio, M.A.T.; Alfonzo, J.D.; Inoue, K.; Sakaguchi, Y.; Suzuki, T.; Sochacka, E.; et al. Identification of 2-methylthio cyclic N6-threonylcarbamoyladenosine (ms2ct6A) as a novel RNA modification at position 37 of tRNAs. *Nucleic Acids Res.* **2017**, *45*, 2124–2136. [CrossRef]

69. Sankaran, B.; Bonnett, S.A.; Shah, K.; Gabriel, S.; Reddy, R.; Schimmel, P.; Rodionov, D.A.; de Crécy-Lagard, V.; Helmann, J.D.; Iwata-Reuyl, D.; et al. Zinc-independent folate biosynthesis: Genetic, biochemical, and structural investigations reveal new metal dependence for GTP cyclohydrolase IB. *J. Bacteriol.* **2009**, *191*, 6936–6949. [CrossRef]

70. Dubois, D.Y.; Blaise, M.; Becker, H.D.; Campanacci, V.; Keith, G.; Giegé, R.; Cambillau, C.; Lapointe, J.; Kern, D. An aminoacyl-tRNA synthetase-like protein encoded by the *Escherichia coli* yadB gene glutamylates specifically tRNAAsp. *Proc. Natl. Acad. Sci. USA* **2004**, *101*, 7530–7535. [CrossRef]

71. El Yacoubi, B.; Lyons, B.; Cruz, Y.; Reddy, R.; Nordin, B.; Agnelli, F.; Williamson, J.R.; Schimmel, P.; Swairjo, M.A.; de Crécy-Lagard, V. The universal YrdC/Sua5 family is required for the formation of threonylcarbamoyladenosine in tRNA. *Nucleic Acids Res.* **2009**, *37*, 2894–2909. [CrossRef]
54. Kimura, S.; Suzuki, T. Fine-tuning of the ribosomal decoding center by conserved methyl-modifications in the *Escherichia coli* 16S rRNA. *Nucleic Acids Res.* **2010**, *38*, 1341–1352. [CrossRef] [PubMed]

55. Emmerechts, G.; Barbe, S.; Herdewijn, P.; Anne, J.; Rozenski, J. Post-transcriptional modification mapping in the *Clostridium acetobutylicum* 16S rRNA by mass spectrometry and reverse transcriptase assays. *Nucleic Acids Res.* **2007**, *35*, 3494–3503. [CrossRef]

56. Desai, P.M.; Culver, G.M.; Rife, J.P. Site-directed mutants of 16S rRNA reveal important RNA domains for KsgA function and 30S subunit assembly. *Biochemistry* **2011**, *50*, 854–863. [CrossRef]

57. Poldermans, B.; Bakker, H.; Van Knippenberg, P.H. Studies on the function of two adjacent N6, N6-dimethyladenosines near the 3' end of 16 S ribosomal RNA of *Escherichia coli*. IV. The effect of the methylgroups on ribosomal subunit interaction. *Nucleic Acids Res.* **1980**, *8*, 143–152. [CrossRef]

58. Grosjean, H.; Edqvist, J.; Straby, K.B.; Gieg'e, R. Enzymatic formation of modified nucleosides in tRNA: Dependence on tRNA architecture. *J. Mol. Biol.* **1996**, *255*, 67–85. [CrossRef]

59. Armengaud, J.; Urbonavicius, J.; Fernandez, B.; Chaussinand, G.; Bujnicki, J.M.; Grosjean, H. N2-methylation of guanosine at position 10 in tRNA is catalyzed by a THUMP domain-containing, S-Adenosylmethionine-dependent methyltransferase, conserved in Archaea and Eukaryota. *J. Biol. Chem.* **2004**, *279*, 37142–37152. [CrossRef]

60. Pallan, P.S.; Kreutz, C.; Bosio, S.; Micura, R.; Egli, M. Effects of N(2),N(2)-dimethylguanosine on RNA structure and stability: Crystal structure of an RNA duplex with tandem m(2) (2)G:A pairs. *RNA* **2008**, *14*, 2452–2453. [CrossRef]

61. Golovina, A.Y.; Dzama, M.M.; Osterman, I.A.; Sergiev, P.V.; Serebryakova, M.V.; Bogdanov, A.A.; Dontsova, O.A. The last rRNA methyltransferase of *E. coli* revealed: The *yhiR* gene encodes adenine-N6 methyltransferase specific for modification of A2030 of 23S ribosomal RNA. *RNA* **2012**, *18*, 1725–1734. [CrossRef] [PubMed]

62. Grosjean, H.; Constantinesco, F.; Foiret, D.; Benachenhou, N. A novel enzymatic pathway leading to 1-methylinosine modification in *Haloferax volcanii* tRNA. *Nucleic Acids Res.* **1995**, *23*, 4312–4319. [CrossRef] [PubMed]

63. Jora, M.; Lobue, P.A.; Ross, R.L.; Williams, B.; Addepalli, B. Detection of ribonucleoside modifications by liquid chromatography coupled with mass spectrometry. *Biochim. Biophys. Acta-Gene Regul. Mech.* **2019**, *1862*, 280–290. [CrossRef]

64. Singh, R.; Farmer, P.B. Liquid chromatography-electrospray ionization-mass spectrometry: The future of DNA adduct detection. *Carcinogenesis* **2005**, *27*, 178–196. [CrossRef]

65. Arnold, H.H.; Keith, G. The nucleotide sequence of phenylalanine tRNA from *Bacillus subtilis*. *Nucleic Acids Res.* **1977**, *4*, 2821–2829. [CrossRef] [PubMed]

66. Matsumoto, T.; Ohta, T.; Kumagai, I.; Oshima, T.; Murao, K.; Hasegawa, T.; Ishikura, H.; Watanabe, K. A thermostable Gm-methylase recognizes the tertiary structure of tRNA 1. *J. Biochem.* **1987**, *101*, 1191–1198. [CrossRef]

67. Vold, B.S.; Keith, D.E., Jr.; Buck, M.; McCloskey, J.A.; Pang, H. Lysine tRNAs from *Bacillus subtilis* 16S: Structural analysis. *Nucleic Acids Res.* **1982**, *10*, 3125–3132. [CrossRef]

68. Kang, B.-I.; Miyauchi, K.; Matuszewski, M.; D'Almeida, G.S.; Rubio, M.A.T.; Alfonzo, J.D.; Inoue, K.; Sakaguchi, Y.; Suzuki, T.; Sochacka, E.; et al. Identification of 2-methylthio cyclic N6-threonylcarbamoyladenosine (ms2ct6A) as a novel RNA modification at position 37 of tRNAs. *Nucleic Acids Res.* **2017**, *45*, 2124–2136. [CrossRef]

69. Sankaran, B.; Bonnett, S.A.; Shah, K.; Gabriel, S.; Reddy, R.; Schimmel, P.; Rodionov, D.A.; de Crécy-Lagard, V.; Helmann, J.D.; Iwata-Reuyl, D.; et al. Zinc-independent folate biosynthesis: Genetic, biochemical, and structural investigations reveal new metal dependence for GTP cyclohydrolase IB. *J. Bacteriol.* **2009**, *191*, 6936–6949. [CrossRef]

70. Dubois, D.Y.; Blaise, M.; Becker, H.D.; Campanacci, V.; Keith, G.; Giegé, R.; Cambillau, C.; Lapointe, J.; Kern, D. An aminoacyl-tRNA synthetase-like protein encoded by the *Escherichia coli* yadB gene glutamylates specifically tRNAAsp. *Proc. Natl. Acad. Sci. USA* **2004**, *101*, 7530–7535. [CrossRef]

71. El Yacoubi, B.; Lyons, B.; Cruz, Y.; Reddy, R.; Nordin, B.; Agnelli, F.; Williamson, J.R.; Schimmel, P.; Swairjo, M.A.; de Crécy-Lagard, V. The universal YrdC/Sua5 family is required for the formation of threonylcarbamoyladenosine in tRNA. *Nucleic Acids Res.* **2009**, *37*, 2894–2909. [CrossRef]
72. Deutsch, C.; El Yacoubi, B.; De Crécy-Lagard, V.; Iwata-Reuyl, D. Biosynthesis of threonylcarbamoyl adenosine (t⁶A), a universal tRNA nucleoside. *J. Biol. Chem.* **2012**, *287*, 13666–13673. [CrossRef]

73. Lauhon, C.T. Mechanism of N6-threonylcarbamoyladenosine (t⁶A) biosynthesis: Isolation and characterization of the Intermediate threonylcarbamoyl-AMP. *Biochemistry* **2012**, *51*, 8950–8963. [CrossRef] [PubMed]

74. Kimura, S.; Miyauchi, K.; Ikeuchi, Y.; Thiaville, P.C.; de Crécy-Lagard, V.; Suzuki, T. Discovery of the beta-barrel-type RNA methyltransferase. *Nucleic Acids Res.* **2014**, *42*, 9350–9365. [CrossRef] [PubMed]

75. Anton, B.P.; Russell, S.P.; Vertrees, J.; Kasif, S.; Raleigh, E.A.; Limbach, P.A.; Roberts, R.J. Functional characterization of the YmcB and YqeV tRNA methylthiotransferases of *Bacillus subtilis*. *Nucleic Acids Res.* **2010**, *38*, 6195–6205. [CrossRef] [PubMed]

76. Arragain, S.; Handelman, S.K.; Forouhar, F.; Wei, F.Y.; Tomizawa, K.; Hunt, J.F.; Douki, T.; Fontecave, M.; Mulliez, E.; Atta, M. Identification of eukaryotic and prokaryotic methylthiotransferase for biosynthesis of 2-methylthio-N6-threonylcarbamoyladenosine in tRNA. *J. Biol. Chem.* **2010**, *285*, 28425–28433. [CrossRef] [PubMed]

77. Ikeuchi, Y.; Shigi, N.; Kato, J.; Nishimura, A.; Suzuki, T. Mechanistic insights into sulfur relay by multiple sulfur mediators involved in thiouridine biosynthesis at tRNA wobble positions. *Mol. Cell.* **2006**, *21*, 97–108. [CrossRef] [PubMed]

78. Zheng, C.; Black, K.A.; Dos Santos, P.C. Diverse mechanisms of sulfur decoration in bacterial tRNA and their cellular functions. *Biomolecules* **2017**, *7*, 33. [CrossRef]

79. Black, K.A.; dos Santos, P.C. Abbreviated pathway for biosynthesis of 2-thiouridine in *Bacillus subtilis*. *J. Bacteriol.* **2015**, *196*, 1952–1962. [CrossRef]

80. Roovers, M.; Kaminska, K.H.; Tkaczuk, K.L.; Gigot, D.; Droogmans, L.; Bujnicki, J.M. The YqfN protein of *Bacillus subtilis* is the tRNA: m¹A22 methyltransferase (TrmK). *Nucleic Acids Res.* **2008**, *36*, 3252–3262. [CrossRef]

81. Guy, M.P.; Podyma, B.M.; Preston, M.A.; Shaheen, H.H.; Krivos, K.L.; Limbach, P.A.; Hopper, A.K.; Phizicky, E.M. Yeast Trm7 interacts with distinct proteins for critical modifications of the tRNAPhe anticodon loop. *RNA* **2012**, *18*, 1921–1933. [CrossRef] [PubMed]

82. Guy, M.P.; Phizicky, E.M. Conservation of an intricate circuit for crucial modifications of the tRNAPhe anticodon loop in eukaryotes. *RNA* **2015**, *21*, 61–74. [CrossRef] [PubMed]

83. Angelova, M.T.; Dimitrova, D.G.; Da Silva, B.; Marchand, V.; Jacquier, C.; Achour, C.; Brazane, M.; Goyenvalle, C.; Bourguignon-Igel, V.; Shehzada, S.; et al. tRNA 2'-O-methylation by a duo of TRM7/FTSJ1 proteins modulates small RNA silencing in Drosophila. *Nucleic Acids Res.* **2020**, *48*, 2050–2072. [CrossRef]

84. Urbonavicius, J.; Skouloubris, S.; Myllykallio, H.; Grosjean, H. Identification of a novel gene encoding a flavin-dependent tRNA: m5U methyltransferase in bacteria–evolutionary implications. *Nucleic Acids Res.* **2005**, *33*, 3955–3964. [CrossRef]

85. Ikeuchi, Y.; Kitachara, K.; Suzuki, T. The RNA acetyltransferase driven by ATP hydrolysis synthesizes N4-acetylcytidine of tRNA anticodon. *EMBO J.* **2008**, *27*, 2194–2203. [CrossRef] [PubMed]

86. Rajakovich, L.J.; Tomlinson, J.; Dos Santos, P.C. Functional analysis of *Bacillus subtilis* genes involved in the biosynthesis of 4-thiouridine in tRNA. *J. Bacteriol.* **2012**, *194*, 4933–4940. [CrossRef] [PubMed]

87. Kim, J.; Xiao, H.; Bonanno, J.B.; Kalyanaraman, C.; Brown, S.; Tang, X.; Al-Obaidi, N.F.; Patskovsky, Y.; Babbitt, P.C.; Jacobson, M.P.; et al. Structure-guided discovery of the metabolite carboxy-SAM that modulates tRNA function. *Nature* **2013**, *498*, 123–126. [CrossRef]

88. Sakai, Y.; Miyauchi, K.; Kimura, S.; Suzuki, T. Biogenesis and growth phase-dependent alteration of 5-methoxycarbonylmethoxyuridine in tRNA anticodons. *Nucleic Acids Res.* **2016**, *44*, 509–523. [CrossRef]

89. Ryu, H.; Grove, T.L.; Almo, S.C.; Kim, J. Identification of a novel tRNA wobble uridine modifying activity in the biosynthesis of 5-methoxyuridine. *Nucleic Acids Res.* **2018**, *46*, 9160–9169. [CrossRef]

90. Sakai, Y.; Kimura, S.; Suzuki, T. Dual pathways of tRNA hydroxylation ensure efficient translation by expanding decoding capability. *Nat. Commun.* **2019**, *10*, 2858. [CrossRef] [PubMed]

91. Lauhon, C.T. Identification and characterization of genes required for 5-hydroxyuridine synthesis in *Bacillus subtilis* and *Escherichia coli*. *J. Bacteriol.* **2019**, *201*, e00433-19. [CrossRef] [PubMed]

92. Bou-Nader, C.; Montémont, H.; Guérineau, V.; Jean-Jean, O.; Brégeon, D.; Hamdane, D. Unveiling structural and functional divergences of bacterial tRNA dihydrouridine synthases: Perspectives on the evolution scenario. *Nucleic Acids Res.* **2018**, *46*, 1386–1394. [CrossRef] [PubMed]
72. Deutsch, C.; El Yacoubi, B.; De Crécy-Lagard, V.; Iwata-Reuyl, D. Biosynthesis of threonylcarbamoyl adenosine (t⁶A), a universal tRNA nucleoside. *J. Biol. Chem.* **2012**, *287*, 13666–13673. [CrossRef]

73. Lauhon, C.T. Mechanism of N6-threonylcarbamoyladenosine (t⁶A) biosynthesis: Isolation and characterization of the Intermediate threonylcarbamoyl-AMP. *Biochemistry* **2012**, *51*, 8950–8963. [CrossRef] [PubMed]

74. Kimura, S.; Miyauchi, K.; Ikeuchi, Y.; Thiaville, P.C.; de Crécy-Lagard, V.; Suzuki, T. Discovery of the beta-barrel-type RNA methyltransferase. *Nucleic Acids Res.* **2014**, *42*, 9350–9365. [CrossRef] [PubMed]

75. Anton, B.P.; Russell, S.P.; Vertrees, J.; Kasif, S.; Raleigh, E.A.; Limbach, P.A.; Roberts, R.J. Functional characterization of the YmcB and YqeV tRNA methylthiotransferases of *Bacillus subtilis*. *Nucleic Acids Res.* **2010**, *38*, 6195–6205. [CrossRef] [PubMed]

76. Arragain, S.; Handelman, S.K.; Forouhar, F.; Wei, F.Y.; Tomizawa, K.; Hunt, J.F.; Douki, T.; Fontecave, M.; Mulliez, E.; Atta, M. Identification of eukaryotic and prokaryotic methylthiotransferase for biosynthesis of 2-methylthio-N6-threonylcarbamoyladenosine in tRNA. *J. Biol. Chem.* **2010**, *285*, 28425–28433. [CrossRef] [PubMed]

77. Ikeuchi, Y.; Shigi, N.; Kato, J.; Nishimura, A.; Suzuki, T. Mechanistic insights into sulfur relay by multiple sulfur mediators involved in thiouridine biosynthesis at tRNA wobble positions. *Mol. Cell.* **2006**, *21*, 97–108. [CrossRef] [PubMed]

78. Zheng, C.; Black, K.A.; Dos Santos, P.C. Diverse mechanisms of sulfur decoration in bacterial tRNA and their cellular functions. *Biomolecules* **2017**, *7*, 33. [CrossRef]

79. Black, K.A.; dos Santos, P.C. Abbreviated pathway for biosynthesis of 2-thiouridine in *Bacillus subtilis*. *J. Bacteriol.* **2015**, *196*, 1952–1962. [CrossRef]

80. Roovers, M.; Kaminska, K.H.; Tkaczuk, K.L.; Gigot, D.; Droogmans, L.; Bujnicki, J.M. The YqfN protein of *Bacillus subtilis* is the tRNA: m¹A22 methyltransferase (TrmK). *Nucleic Acids Res.* **2008**, *36*, 3252–3262. [CrossRef]

81. Guy, M.P.; Podyma, B.M.; Preston, M.A.; Shaheen, H.H.; Krivos, K.L.; Limbach, P.A.; Hopper, A.K.; Phizicky, E.M. Yeast Trm7 interacts with distinct proteins for critical modifications of the tRNAPhe anticodon loop. *RNA* **2012**, *18*, 1921–1933. [CrossRef] [PubMed]

82. Guy, M.P.; Phizicky, E.M. Conservation of an intricate circuit for crucial modifications of the tRNAPhe anticodon loop in eukaryotes. *RNA* **2015**, *21*, 61–74. [CrossRef] [PubMed]

83. Angelova, M.T.; Dimitrova, D.G.; Da Silva, B.; Marchand, V.; Jacquier, C.; Achour, C.; Brazane, M.; Goyenvalle, C.; Bourguignon-Igel, V.; Shehzada, S.; et al. tRNA 2'-O-methylation by a duo of TRM7/FTSJ1 proteins modulates small RNA silencing in Drosophila. *Nucleic Acids Res.* **2020**, *48*, 2050–2072. [CrossRef]

84. Urbonavicius, J.; Skouloubris, S.; Myllykallio, H.; Grosjean, H. Identification of a novel gene encoding a flavin-dependent tRNA: m5U methyltransferase in bacteria–evolutionary implications. *Nucleic Acids Res.* **2005**, *33*, 3955–3964. [CrossRef]

85. Ikeuchi, Y.; Kitachara, K.; Suzuki, T. The RNA acetyltransferase driven by ATP hydrolysis synthesizes N4-acetylcytidine of tRNA anticodon. *EMBO J.* **2008**, *27*, 2194–2203. [CrossRef] [PubMed]

86. Rajakovich, L.J.; Tomlinson, J.; Dos Santos, P.C. Functional analysis of *Bacillus subtilis* genes involved in the biosynthesis of 4-thiouridine in tRNA. *J. Bacteriol.* **2012**, *194*, 4933–4940. [CrossRef] [PubMed]

87. Kim, J.; Xiao, H.; Bonanno, J.B.; Kalyanaraman, C.; Brown, S.; Tang, X.; Al-Obaidi, N.F.; Patskovsky, Y.; Babbitt, P.C.; Jacobson, M.P.; et al. Structure-guided discovery of the metabolite carboxy-SAM that modulates tRNA function. *Nature* **2013**, *498*, 123–126. [CrossRef]

88. Sakai, Y.; Miyauchi, K.; Kimura, S.; Suzuki, T. Biogenesis and growth phase-dependent alteration of 5-methoxycarbonylmethoxyuridine in tRNA anticodons. *Nucleic Acids Res.* **2016**, *44*, 509–523. [CrossRef]

89. Ryu, H.; Grove, T.L.; Almo, S.C.; Kim, J. Identification of a novel tRNA wobble uridine modifying activity in the biosynthesis of 5-methoxyuridine. *Nucleic Acids Res.* **2018**, *46*, 9160–9169. [CrossRef]

90. Sakai, Y.; Kimura, S.; Suzuki, T. Dual pathways of tRNA hydroxylation ensure efficient translation by expanding decoding capability. *Nat. Commun.* **2019**, *10*, 2858. [CrossRef] [PubMed]

91. Lauhon, C.T. Identification and characterization of genes required for 5-hydroxyuridine synthesis in *Bacillus subtilis* and *Escherichia coli*. *J. Bacteriol.* **2019**, *201*, e00433-19. [CrossRef] [PubMed]

92. Bou-Nader, C.; Montémont, H.; Guérineau, V.; Jean-Jean, O.; Brégeon, D.; Hamdane, D. Unveiling structural and functional divergences of bacterial tRNA dihydrouridine synthases: Perspectives on the evolution scenario. *Nucleic Acids Res.* **2018**, *46*, 1386–1394. [CrossRef] [PubMed]
93. Ofengand, J.; Bakin, A. Mapping to nucleotide resolution of pseudouridine residues in large subunit ribosomal RNAs from representative eukaryotes, prokaryotes, archaebacteria, mitochondria and chloroplasts. *J. Mol. Biol.* **1997**, *266*, 246–268. [CrossRef] [PubMed]

94. Del Campo, M.; Kaya, Y.; Ofengand, J. Identification and site of action of the remaining four putative pseudouridine synthases in *Escherichia coli*. *RNA* **2001**, *7*, 1603–1615. [PubMed]

95. Niu, L.; Ofengand, J. Cloning and characterization of the 23S RNA pseudouridine 2633 synthase from *Bacillus subtilis*. *Biochemistry* **1999**, *38*, 629–635. [CrossRef] [PubMed]

96. Addepalli, B.; Limbach, P.A. Pseudouridine in the anticodon of *Escherichia coli* tRNATyr(QΨA) Is catalyzed by the dual specificity enzyme RluF. *J. Biol. Chem.* **2016**, *291*, 22327–22337. [CrossRef]

97. Taniguchi, T.; Miyauchi, K.; Nakane, D.; Miyata, M.; Muto, A.; Nishimura, S.; Suzuki, T. Decoding system for the AUA codon by tRNAle with the UAU anticodon in *Mycoplasma mobile*. *Nucleic Acids Res.* **2013**, *41*, 2621–2631. [CrossRef]

98. Hill, P.J.; Abibi, A.; Albert, R.; Andrews, B.; Gagnon, M.M.; Gao, N.; Grebe, T.; Hajec, L.I.; Huang, J.; Livchak, S.; et al. Selective inhibitors of bacterial t-RNA-(N1G37) methyltransferase (TrmD) that demonstrate novel ordering of the lid domain. *J. Med. Chem.* **2013**, *56*, 7278–7288. [CrossRef]

99. Kambampati, R.; Lauhon, C.T. MnmA and IscS are required for in vitro 2-thiouridine biosynthesis in *Escherichia coli*. *Biochemistry* **2003**, *42*, 1109–1117. [CrossRef]

100. Meeske, A.J.; Rodrigues, C.D.A.; Brady, J.; Lim, H.C.; Bernhardt, T.G.; Rudner, D.Z. High-throughput genetic screens identify a large and diverse collection of new sporulation genes in *Bacillus subtilis*. *PLoS Biol.* **2016**, *14*, e1002341. [CrossRef]

101. Sayer, C.V.; Barat, B.; Popham, D.L. Identification of L-Valine-initiated-germination-active genes in *Bacillus subtilis* using Tn-seq. *PLoS ONE* **2019**, *14*, e0218220. [CrossRef] [PubMed]

<figure><img src="image_10.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
93. Ofengand, J.; Bakin, A. Mapping to nucleotide resolution of pseudouridine residues in large subunit ribosomal RNAs from representative eukaryotes, prokaryotes, archaebacteria, mitochondria and chloroplasts. *J. Mol. Biol.* **1997**, *266*, 246–268. [CrossRef] [PubMed]

94. Del Campo, M.; Kaya, Y.; Ofengand, J. Identification and site of action of the remaining four putative pseudouridine synthases in *Escherichia coli*. *RNA* **2001**, *7*, 1603–1615. [PubMed]

95. Niu, L.; Ofengand, J. Cloning and characterization of the 23S RNA pseudouridine 2633 synthase from *Bacillus subtilis*. *Biochemistry* **1999**, *38*, 629–635. [CrossRef] [PubMed]

96. Addepalli, B.; Limbach, P.A. Pseudouridine in the anticodon of *Escherichia coli* tRNATyr(QΨA) Is catalyzed by the dual specificity enzyme RluF. *J. Biol. Chem.* **2016**, *291*, 22327–22337. [CrossRef]

97. Taniguchi, T.; Miyauchi, K.; Nakane, D.; Miyata, M.; Muto, A.; Nishimura, S.; Suzuki, T. Decoding system for the AUA codon by tRNAle with the UAU anticodon in *Mycoplasma mobile*. *Nucleic Acids Res.* **2013**, *41*, 2621–2631. [CrossRef]

98. Hill, P.J.; Abibi, A.; Albert, R.; Andrews, B.; Gagnon, M.M.; Gao, N.; Grebe, T.; Hajec, L.I.; Huang, J.; Livchak, S.; et al. Selective inhibitors of bacterial t-RNA-(N1G37) methyltransferase (TrmD) that demonstrate novel ordering of the lid domain. *J. Med. Chem.* **2013**, *56*, 7278–7288. [CrossRef]

99. Kambampati, R.; Lauhon, C.T. MnmA and IscS are required for in vitro 2-thiouridine biosynthesis in *Escherichia coli*. *Biochemistry* **2003**, *42*, 1109–1117. [CrossRef]

100. Meeske, A.J.; Rodrigues, C.D.A.; Brady, J.; Lim, H.C.; Bernhardt, T.G.; Rudner, D.Z. High-throughput genetic screens identify a large and diverse collection of new sporulation genes in *Bacillus subtilis*. *PLoS Biol.* **2016**, *14*, e1002341. [CrossRef]

101. Sayer, C.V.; Barat, B.; Popham, D.L. Identification of L-Valine-initiated-germination-active genes in *Bacillus subtilis* using Tn-seq. *PLoS ONE* **2019**, *14*, e0218220. [CrossRef] [PubMed]

<figure><img src="image_10.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).