A tRNA modification in Mycobacterium tuberculosis facilitates optimal intracellular growth

Version of Record

The authors declare this version of their article to be the Version of Record.

About eLife's process

Download
Cite
Share
CommentOpen annotations (there are currently 0 annotations on this page).

516 views

Version of Record published: September 27, 2023 (This version)
Reviewed preprint version 2: August 8, 2023 (Go to version)
Reviewed preprint version 1: April 12, 2023 (Go to version)
Sent for peer review: February 27, 2023
Preprint posted: February 20, 2023 (Go to version)

1. Of interest
Evolution towards simplicity in bacterial small heat shock protein system

Piotr Karaś, Klaudia Kochanowicz ... Krzysztof Liberek

Research Article Dec 8, 2023
Further reading

Abstract
eLife assessment
Introduction
Results
Discussion
Materials and methods
Data availability
References
Peer review
Author response
Article and author information
Metrics

Abstract

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

eLife assessment

This is a valuable addition to the literature as it helps us understand the role of tRNA modifying enzymes in Mycobacterium tuberculosis. By knocking out one of the enzymes, the authors convincingly demonstrate the importance of tRNA-modifying enzymes for intra-host growth of tubercle bacteria. Some of the claims regarding modification as well as the role in virulence could be strengthened through further bioinformatics and phylogenetic analyses as well as experimental approaches. The work will be of interest to microbiologists.

https://doi.org/10.7554/eLife.87146.3.sa0

About eLife assessments

Introduction

tRNA is an adaptor molecule that enables protein synthesis by converting the triplet genetic code in mRNA into amino acids. The fidelity of base pairing between mRNA codons and tRNA anticodons is monitored within ribosomes and is critical for properly incorporating the amino acids bound to the 3′ ends of tRNAs into growing polypeptides. For optimal translation, the abundances, and properties of tRNA isoacceptors are fine tuned by diverse mechanisms, including chemical modifications (Huang and Hopper, 2016; Björk and Hagervall, 2014; Shepherd and Ibba, 2015; Torrent et al., 2018). Dysregulation of tRNA abundance and/or structure leads to defective decoding and results in ribosome pausing and collisions, protein misfolding, stress responses and can have detrimental or lethal effects on the cell (Orellana et al., 2022; Suzuki, 2021; Liu et al., 2022; Goodarzi et al., 2016; Delaunay et al., 2022; Nedialkova and Leidel, 2015).

Chemical modifications of tRNA (tRNA modifications) are found in all kingdoms of life and fine-tune tRNA properties including mRNA decoding efficiency, recognition by aminoacyl-tRNA synthetases, half-life, and structural stability (Björk and Hagervall, 2014; Nedialkova and Leidel, 2015; Kimura and Waldor, 2019; Giegé and Eriani, 2023). Modifications are prevalent in the anticodon loop, particularly at the first letter of the anticodon. Modifications of the anticodon loop directly modulate codon recognition, whereas modifications in the tRNA body region primarily stabilize tRNA tertiary structure, protecting them from degradation in the cell. tRNA modifications are generated by dedicated site-specific enzymes referred to as tRNA modifying enzymes. tRNA modifications have been extensively characterized in a few model organisms (Björk and Hagervall, 2014; de Crécy-Lagard and Jaroch, 2021), but their profiles, regulation, and functions in non-model organisms, including bacterial pathogens, are understudied (de Crécy-Lagard and Jaroch, 2021).

tRNA sequencing (tRNA-seq) allows for the rapid and systematic prediction of many tRNA modification sites (Zhang et al., 2022; Zheng et al., 2015). We recently developed a comparative tRNA-seq protocol to profile tRNA modifications in organisms with uncharted tRNA modification profiles; in Vibrio cholerae, this approach led to the discovery of a new RNA modification and RNA editing process (Kimura et al., 2020). tRNA-seq enables rapid prediction of modified sites through detection of reverse transcription-derived signatures, such as nucleotide misincorporation and early termination, both of which occur more frequently at modified sites. Furthermore, several chemical treatments of tRNA can convert modifications that are not recognizable as reverse transcription-derived signatures into detectable signals, expanding the repertoire of modifications that can be distinguished by tRNA-seq (Motorin and Helm, 2019; Finet et al., 2022; Draycott et al., 2022; Dai et al., 2023).

Mycobacterium tuberculosis (Mtb), the agent of tuberculosis (TB), is a global pathogen that caused >10.5 million cases and over 1.5 million deaths worldwide in 2020 (W.H.O, 2021). Several studies have uncovered roles for non-Mtb mycobacterial tRNA modifications in stress responses, adaptation to environmental changes, and persister cell formation (Chionh et al., 2016). Mycobacterium bovis BCG, an organism closely related to Mtb, responds to hypoxia by reprogramming 40 ribonucleoside modifications in tRNA to facilitate translation of a subset of proteins that promote survival in hypoxic conditions (Chionh et al., 2016). To date, studies of the profiles and functions of Mtb tRNA modifications have been limited.

Here, we conducted tRNA-seq in Mtb. We assigned modifications to many of the reverse transcription-derived signatures identified, using information on the presence of homologs of known modifying enzymes in Mtb. Chemical treatments of tRNAs carried out prior to tRNA-seq increased the detectability of certain modifications. We constructed two Mtb deletion mutant strains, with deletions of mnmA and truB, and confirmed that the absence of these modification enzymes eliminated the predicted signals in tRNA-seq data. Furthermore, while deletion of mnmA in Mtb did not affect the pathogen’s growth in in vitro laboratory growth conditions, the mnmA knockout strain’s growth was attenuated in a macrophage infection model. Our findings suggest that tRNA modifications warrant further study as we unravel the complexity of Mtb infections, as they may serve as targets for new therapeutics.

Results

In silico prediction of Mtb tRNA modifying enzymes

To predict tRNA modifications in Mtb, we used Basic Local Alignment Search Tool (BLAST) (Altschul et al., 1990) to identify homologs of all tRNA modification enzymes registered in Modomics in the Mtb genome (Boccaletto et al., 2022; Supplementary file 1). With a stringent threshold (E-value <1 × 10⁻¹⁰), 31 Mtb genes homologous to genes encoding known RNA modification enzymes were identified (Supplementary file 2). Twenty-three of these genes are predicted to synthesize 16 tRNA modifications in Mtb (Supplementary files 2 and 3), including miaA and miaB for 2-methylthio-6-isopentenyl-adenosine (ms²i⁶A), tsaD, tsaB, tsaE, and sua5 for N⁶-threonylcarbamoyladenosine (t⁶A), mnmA for 2-thiouridine (s²U), truB, truA, rluA, and pus9 for pseudouridine (Ψ), trm2 for 5-methyluridine (m⁵U), trmD for 1-methylguanosine (m¹G), trmI for 1-methyladenosine (m¹A), trmL for 2′-O-methylcytidine (Cm) or 2′-O-methyluridine (Um), trmB for 7-methylguanosine (m⁷G), two trmH for 2′-O-methylguanosine (Gm), trmR for 5-methoxyuridine (mo⁵U), trcM for 5-methylcytidine (m⁵C), dusB for dihydrouridine (D), tadA for inosine (I), and tilS for lysidine (k²C). While eight additional genes met the threshold for homology to modification enzymes, they exhibited greater similarity to a Glutamyl-tRNA synthase; Rv2992c, ribosome associated GTPases (Der; Rv1713 and Era; Rv2364c), molybdopterin biosynthesis proteins (MoeW; Rv2338c and MoeB2; Rv3116), thiosulfate sulfur transferases (CysA; Rv3117 and SseA; Rv3283) and a riboflavin biosynthesis protein (RibG; Rv1409), and likely do not correspond to tRNA modification enzymes (Supplementary file 2).

We mined data from genome-wide Tn-seq and CRISPRi screens (DeJesus et al., 2017; Bosch et al., 2021) to assess the impacts of Mtb tRNA modifications on its growth. Five genes encoding Mtb tRNA modifying enzymes, trmD, tilS, tadA, tsaC2, and miaA, were reported to be essential for Mtb growth in both Tn-seq (DeJesus et al., 2017) and CRISPRi screens (Bosch et al., 2021; Supplementary file 3), suggesting that the modifications they produce are critical for tRNA functions. Indeed, E. coli trmD, tilS, tadA, and tsaC2 are also essential for growth and critical for codon decoding, aminoacylation, and reading frame maintenance (Masuda et al., 2022; Soma et al., 2003; Wolf et al., 2002; El Yacoubi et al., 2009). The Mtb modifications synthesized by these enzymes likely have similar impacts on tRNA functions. Unexpectedly, one modifying enzyme, MiaA, is non-essential in E. coli grown in nutrient-rich medium, but apparently essential in Mtb, suggesting that i⁶A, the modification introduced by MiaA, may have more profound roles in Mtb translation than in E. coli.

Several of the putative Mtb tRNA modifying enzymes are conserved across all three domains of life (e.g., TruB and TsaD) (Figure 1). By contrast, some enzymes were limited to bacterial species closely related to Mtb, possibly suggesting their species-specific physiological roles, including pathogenesis. For example, TrmI (Rv2118c) homologs are widely present in Archaea and Eukaryotes, but are sparsely distributed in bacteria, where they are primarily limited to Actinomycetia, including Mtb, and several thermophilic bacterial species (e.g., Thermus thermophilus). TrmI synthesizes m¹A at position 58 in eukaryotes (Anderson et al., 1998) and at position 57/58 in archaea (Guelorget et al., 2010); indeed, TrmI in Mtb has been proven to generate m¹A at position 58 (Varshney et al., 2004). Pus9 (Rv3300c) has many homologs exhibiting weak similarity across eukaryotes and bacteria. However, some species, including Pseudomonas aeruginosa, Acinetobacter baumannii, Neisseria gonorrhoeae, and several species of Actinomycetia, encode proteins showing strong homology to Mtb Pus9, suggesting that this set of pseudouridylases have a distinctive property such as substrate specificity.

Figure 1

Download asset Open asset

Phylogenetic distribution of *Mtb* tRNA modifying enzyme homologs.

Heatmap of log10 E-values from BLAST search results. BLAST searches were conducted against 120 manually picked organisms using *Mtb* tRNA modifying enzymes as queries. When one organism has multiple hits, the lowest log10(Eval) values among hits are shown. iTol (Letunic and Bork, 2021) was used to depict the results.

Figure 1—source data 1 log10 E-values from BLAST searching for the homologs of Mtb tRNA modifying enzymes in 120 organisms.: https://cdn.elifesciences.org/articles/87146/elife-87146-fig1-data1-v1.xlsx
Download elife-87146-fig1-data1-v1.xlsx

tRNA modifying enzymes predicted in Mtb were observed in mycobacteria species, including Mycolicibacterium smegmatis and Mycobacteroides abscessus (Figure 1), suggesting that tRNA modification patterns are similar among mycobacterium species.

Profiling Mtb tRNA modification sites by tRNA-seq

To begin profiling detectable Mtb tRNA modifications, we sequenced tRNAs isolated from wild-type Mtb strain H37Rv grown in 7H9 medium. In this protocol, tRNAs are first reversed transcribed to cDNA (Kimura et al., 2020). During cDNA synthesis, chemical modifications on tRNA nucleotides disrupt Watson–Crick base pairing and increase the frequency of reverse transcriptase errors, leading to incorporation of the incorrect nucleotide or early termination of cDNA synthesis (Kellner et al., 2010). These reverse transcription-derived ‘signatures’ typically correspond to modified sites (Zhang et al., 2022; Zheng et al., 2015) and are depicted in the heatmap in Figure 2.

Figure 2 with 1 supplement see all

Download asset Open asset

Heatmap of misincorporation and early termination frequency in sequencing of tRNAs from wild-ype *Mtb*.

Heatmaps show misincorporation (A) and termination (B) frequencies at all positions across tRNAs (read 5′ to 3′). Predicted modifications are labeled based on similarity to known modifications in other organisms and the presence of the tRNA modifying enzyme homologs (Supplementary file 3). The positions with more than 10% misincorporation in *Mtb* but not in *E. coli* are depicted in white in A. Representative data of two independent experiments with similar results are shown. (C) *M. tuberculosis* tRNA modifications predicted in this study. Schematic tRNA secondary structure with sites of modifications identified either by the presence of modifying enzymes and/or tRNA sequencing (tRNA-seq). Modifications and tRNA species that are not observed in *E. coli* are shown in red. Modifications and positions that are predicted by both RT-derived signature and the presence of the homologs of tRNA modifying enzymes are shown in yellow (without chemical treatment) and green (with chemical treatment), whereas modifications that are only predicted by the presence of the homologs are shown in light blue. Genes reported to be essential in *Mtb* are shown in bold.

Figure 2—source data 1 Misincorporation and early termination frequencies in sequencing of tRNAs from wild-type Mtb.: https://cdn.elifesciences.org/articles/87146/elife-87146-fig2-data1-v1.xlsx
Download elife-87146-fig2-data1-v1.xlsx

Comparison of the reverse transcription-derived signatures observed in Mtb to E. coli, where tRNA modifications have been well characterized (Kimura et al., 2020) enables the prediction of the presence of common modifications, including ms²i⁶A, m¹G, I, and k²C (Figure 2 and Figure 2—figure supplement 1). These predictions are strongly supported by the set of tRNA modification enzymes identified in the Mtb genome (Figure 1 and Supplementary file 3), including miaA and miaB (ms²i⁶A), trmD (m¹G), tadA (I), and tilS (k²C). Some tRNA modifications were observed in Mtb but are not present in E. coli. In other actinobacteria, A58 and A59 are likely modified to m¹A (Schwartz et al., 2018), and since trmI, the methylase that generates this modification is present in Mtb (Varshney et al., 2004), most Mtb tRNAs likely contain this modification as well. Nucleoside variation in the sequence of tRNA genes can account for some of the variations in modified sites between E. coli and Mtb. For example, in Mtb, termination signatures derived from G at position 37 were detected in tRNA-Arg2, -Gln1, and -Gln2, whereas this position in these tRNAs in E. coli are modified A, such as m²A and m⁶A, which are silent in tRNA-seq. Since m¹G induces strong termination during reverse transcription, these positions are likely modified to m¹G, as observed in the Bacillus subtilis tRNA-Arg2 gene position 37 G (Jühling et al., 2009).

tRNA samples were also treated with several chemical treatments prior to sequencing, to expand the set of tRNA modifications detectable by tRNA-seq. These treatments included iodoacetamide (IAA), for detection of sulfur modifications, 1-cyclohexyl-(2-morpholinoethyl) carbodiimide (CMC) for detection of Ψ, and alkali for detection of D and m⁷G. The chemical treatment protocols were first carried out with E. coli, to validate the methods.

IAA is a thiol-reactive compound that covalently attaches carboxyamidomethyl to thiolated uridines via nucleophilic substitution (Herzog et al., 2017) and modified s⁴U is detected as C instead of U. In IAA-treated samples, positions 8 and 9, corresponding to s⁴U in many tRNA species, had high misincorporation frequencies, confirming that IAA treatment modifies s⁴U, leading to elevated misincorporation (Figure 3—figure supplement 1). Furthermore, we observed higher misincorporation and termination signals at the positions corresponding to other sulfur modification, s²C, s²U and their derivatives, such as position 32 in tRNA-Arg3, -Arg5, -Ser3, and Arg4, and 34 in tRNA-Glu, -Gln1, and -Lys (Figure 3 and Figure 3—figure supplement 1), revealing that IAA treatment facilitates the detection of not only s⁴U but also additional sulfur modifications, which are weakly detected without the IAA treatment.

Figure 3 with 1 supplement see all

Download asset Open asset

Iodoacetamide (IAA) treatment promotes detection of sulfur modifications on tRNAs by enhancing termination signals.

Heatmaps of the termination signals of *E. coli* tRNAs treated with (A) or without (B) IAA. Known modification sites, including sulfur modifications (s⁴U, s²C, s²U in white) are shown. (C) Termination frequency at s⁴U (n = 48), s²C (n = 5), and s²U (n = 3) sites of tRNAs treated with or without IAA. The experiment was performed once.

Figure 3—source data 1 Termination frequencies of E. coli tRNAs treated with or without iodoacetamide (IAA).: https://cdn.elifesciences.org/articles/87146/elife-87146-fig3-data1-v1.xlsx
Download elife-87146-fig3-data1-v1.xlsx

Next, we applied IAA treatment to Mtb tRNA-seq. IAA treatment increased termination signals from position 34 in tRNA-Glu1, -Gln1, and -Lys, which contain s²U derivatives in E. coli (Figure 4A). The 2-thiouridine modification is carried out by MnmA in E. coli (Kambampati and Lauhon, 2003), and a homolog, Rv3024c, of this enzyme was identified in the Mtb genome (Supplementary file 3; Kapopoulou et al., 2011). We used double-stranded DNA-based recombineering to delete Rv3024c in Mtb, yielding strain MtbΔmnmA (van Kessel and Hatfull, 2007). Sequencing of tRNA isolated from MtbΔmnmA with prior IAA treatment showed reduced termination signals from position 34 in tRNA-Glu1, -Gln1, and -Lys in treated samples, indicating that Rv3024c plays a critical role in the modification responsible for increased termination frequency derived from position 34 in these tRNAs (Figure 4B, C). Together, these observations strongly suggest that Rv3024c encodes an MnmA-like enzyme that sulfurates position 34 uridines in three Mtb tRNA isoacceptors.

Figure 4

Download asset Open asset

Heatmap and plot of early termination frequency from sequencing of tRNAs from wild-type and *MtbΔmnmA* with and without RNA alkylation.

Heatmap of early termination frequencies across tRNA molecules and positions for wild-type (WT) (A) and *MtbΔmnmA* (B). Sulfur modification is shown in white. (C) Plot of termination frequencies at position 37 in WT *Mtb* and *MtbΔmnmA* for lysine_UUU, glutamate_UUG, and glutamine_UUC isoacceptors (n = 3). IAA: iodoacetamide. The experiment was performed once.

Figure 4—source data 1 Termination frequencies from sequencing of tRNAs from wild-type and MtbΔmnmA with and without RNA alkylation.: https://cdn.elifesciences.org/articles/87146/elife-87146-fig4-data1-v1.xlsx
Download elife-87146-fig4-data1-v1.xlsx

As shown previously (Carlile et al., 2014), CMC treatment increased both misincorporation and termination signatures at a subset of Ψs in E. coli tRNAs (Figure 5, and Figure 5—figure supplement 1 and Figure 5—figure supplement 2). In addition, CMC-treated samples showed increased frequencies of both misincorporation and termination at positions 16, 17, 20, and 20A corresponding to D, and m⁷G at position 46. Both modifications are known to undergo base elimination in mild alkali conditions (Marchand, 2021). Since these signals were also observed in the reaction condition in which CMC was not added, these signals are likely attributable to the alkali treatment that is common to both conditions.

Figure 5 with 4 supplements see all

Download asset Open asset

1-Cyclohexyl-(2-morpholinoethyl) carbodiimide (CMC) and alkali treatment facilitate detection of ψ, D, and m⁷G modifications in E. *coli.*

Heatmaps of the misincorporation signals of *E. coli* tRNAs treated with (A) or without (B) CMC. In both conditions, tRNAs are incubated in alkali condition. Known Ψ, D, and m⁷G sites are shown. Ψ is shown in white. (C) Misincorporation frequency at known Ψ (n = 80), D (n = 76), and m⁷G (n = 25) sites of tRNAs treated with CMC+/alkali+ or CMC−/alkali+, and tRNAs without treatment. The experiment was performed once.

Figure 5—source data 1 Misincorporation frequencies from sequencing of E. coli tRNAs treated with or without 1-cyclohexyl-(2-morpholinoethyl) carbodiimide (CMC).: https://cdn.elifesciences.org/articles/87146/elife-87146-fig5-data1-v1.xlsx
Download elife-87146-fig5-data1-v1.xlsx

Reverse transcription-derived signatures derived from alkali-treated D showed a distinctive pattern. With this treatment, when two Ds are at consecutive positions, for example, D16 and D17, termination signals were elevated at the following position, that is, position 18 (Figure 5—figure supplement 1, Figure 5—figure supplement 3, and Figure 5—figure supplement 4). Furthermore, alkali treatment also led to higher misincorporation frequencies at singlet Ds (Figure 5, Figure 5—figure supplement 3, and Figure 5—figure supplement 4). Thus, termination and misincorporation signatures enabled the prediction of known E. coli tRNA sites modified to D.

CMC/alkali treatment facilitated the identification of additional modifications in Mtb tRNAs. Regardless of CMC treatment, alkali-treated samples showed increased misincorporation and termination frequencies derived from U located at positions 16, 17, 20, and 20A (Figure 6 and Figure 6—figure supplement 1). As observed in E. coli, termination signals at positions 18 and 21 likely correspond to consecutive Ds at positions 16 and 17, and 20 and 20A, respectively. Mtb Rv0823c is a homolog of dihydrouridylase DusB (Supplementary file 3), which likely accounts for the synthesis of D at these positions. Furthermore, alkali treatment increased the misincorporation frequencies at G at position 46 (Figure 6—figure supplement 2). Since Rv0208c is a homolog of TrmB, which synthesizes m⁷G at position 46 in E. coli, multiple Mtb tRNA species likely contain m⁷G at position 46 (Figure 6 and Figure 6—figure supplement 2).

Figure 6 with 3 supplements see all

Download asset Open asset

Heatmap of early termination frequency from sequencing of tRNAs isolated from wild-type (WT) and *MtbΔtruB* following 1-cyclohexyl-(2-morpholinoethyl) carbodiimide (CMC) treatment.

Heatmap of early termination frequencies across tRNA molecules and positions for WT (left) and *MtbΔtruB* (right). Termination signals derived from position 55 are shown in white. The experiment was performed once.

Figure 6—source data 1 Termination frequencies across tRNA molecules and positions for wild-type (WT) and MtbΔtruB.: https://cdn.elifesciences.org/articles/87146/elife-87146-fig6-data1-v1.xlsx
Download elife-87146-fig6-data1-v1.xlsx

CMC treatment also increased the termination frequency at several sites. Termination signatures derived from position 55, which is exclusively uridine in all tRNA species, increased in most tRNA species, suggesting that Mtb tRNAs contain pseudouridines at this position (Jühling et al., 2009). Rv2793c is an Mtb homolog of E. coli TruB and deletion of Rv2793c reduced the termination frequencies at this position in tRNAs isolated from MtbΔtruB (Figure 6 and Figure 6-figure supplement 3; van Kessel and Hatfull, 2007). Together, these observations suggest that Rv2793c encodes a TruB-like enzyme that modifies position 55 uridines to Ψ across tRNA species. Furthermore, the presence of a TruA homolog in Mtb (Rv3455c) suggests that in multiple tRNA species U at positions 38–40 can be modified to Ψ. Indeed, the termination signatures derived from positions 38 and 39 increased depending on CMC treatment, strongly suggesting that these positions are modified to Ψ (Figure 6 and Figure 6—figure supplement 3).

In total, among 16 tRNA Mtb tRNA modifications predicted by the presence of tRNA modifying enzymes (Figure 1 and Supplementary file 3), 9 species of modifications were detected based on reverse-transcription-derived signatures (Figure 2C).

Growth of MtbΔmnmA is attenuated in a macrophage infection model

To address whether Mtb tRNA modifications impact the pathogen’s growth in the host environment, we used the MtbTnDB transposon insertion sequencing (Tn-seq) database (Zhang et al., 2013; Jinich et al., 2021) to determine if transposon insertions in genes encoding tRNA modifying enzymes have been associated with in vivo growth defects. Transposon insertions in mnmA were reported to attenuate Mtb growth in mice infected with a library of Mtb transposon mutants, suggesting that mnmA facilitates Mtb growth in vivo. We found that the growth of WT and MtbΔmnmA were similar in 7H9 medium (Figure 7), suggesting the absence of s²U modification at position 34 does not impair Mtb growth in culture. In contrast, the MtbΔmnmA mutant was significantly impaired for growth in a macrophage infection model (Jinich et al., 2021; Figure 7). Defective growth of the MtbΔmnmA mutant was also observed in macrophages treated with all-trans retinoic acid (ATRA), which promotes macrophage control of Mtb infection (Babunovic et al., 2022). These observations strongly suggest that modification of U to s²U by MnmA facilitates Mtb growth in macrophages.

Figure 7

Download asset Open asset

*MtbΔmnmA* is attenuated in a macrophage infection model.

(A) Wild-type and *MtbΔmnmA* do not display growth differences in 7H9 medium. (B) Auto-luminescent wild-type and *ΔmnmA Mtb* strains were diluted to a multiplicity of infection of 2 bacteria per mouse bone marrow-derived macrophage with or without all-trans retinoic acid (ATRA). Survival was measured by luminescence and normalized to luminescence at time 0. Average values from three independent cultures (n=3) are shown with standard deviations.

Figure 7—source data 1 Mtb growth curve.: https://cdn.elifesciences.org/articles/87146/elife-87146-fig7-data1-v1.xlsx
Download elife-87146-fig7-data1-v1.xlsx

Discussion

Here, we profiled Mtb tRNA modifications with tRNA-seq to provide the first maps of the tRNA modification landscape in this global pathogen. In total, nine modifications, including six modifications without chemical treatment, were identified based on reverse transcription-derived signatures. CMC/alkali treatment and IAA treatment further identified Ψ and m⁷G and sulfur modifications, respectively. Although we did not chemically validate the modifications predicted by tRNA-seq with mass spectrometry, the identification of Mtb homologs of tRNA modifying enzymes strongly bolsters the RT-signature-based predictions. Furthermore, the deletion of truB and mnmA genes in Mtb eliminated the respective modification signatures of pseudouridine and s²U, validating that the enzymes encoded by these genes synthesize these modifications. Finally, the growth defect of the ΔmnmA strain within macrophages but not in a nutrient-rich medium suggests that s²U tRNA modification facilitates Mtb adaptation to the host intracellular environment.

IAA treatment was developed for detecting s⁴U modification in pulse-chase experiments to measure RNA turnover (Herzog et al., 2017). We found that this treatment enhances the RT signatures of additional sulfur modifications as well, including s²C and s²U, indicating that IAA should have general utility in tRNA-seq-based profiling of sulfur modifications, including in studies tracking changes in tRNA sulfuration (Edwards et al., 2022; Laxman et al., 2013).

Multiple modifications do not cause reverse transcriptase errors. Although tRNA-seq provides a simple and rapid method for profiling the landscape of modifications in all tRNA species, this approach does not enable comprehensive identification of modifications. Here, several modifications, including t⁶A, m⁵U, mo⁵U, m⁵C, Cm, Um, and Gm, predicted by the presence of Mtb homolog of tRNA modifying enzymes, such as TsaC2, TsaB, TsaD, TsaE, Trm2, TrcM, TrmR, TrmL, and TrmH, were not detected in tRNA-seq. Additionally, CMC treatment may not determine all Ψ positions, thus, the target of Rv3300 and Rv1540 remain unclear. Since these genes are similar to E. coli rluA, which also targets rRNA, these genes may target rRNAs instead of tRNAs. Mapping and further analysis of these modifications will require different approaches such as tRNA purification and RNA mass spectrometric analysis.

Several strong RT signatures were detected in Mtb tRNAs but not in E. coli (Figure 2A), including G45 in tRNA-Gly2. Most of these signals are strict misincorporations of C or T at A or G position, respectively. Since modifications at a purine position, for example, m¹G and m¹A, generally cause random misincorporation of the other three nucleosides, these Mtb-specific misincorporation signatures are not likely to be derived from modification. Further mass spectrometric analysis of purified tRNAs will be necessary to assess whether these signals are indeed derived from Mtb-specific modifications; such Mtb-specific modifications could be related to Mtb pathogenesis and/or survival in the host.

The role of s²U in Mtb appears to be unusual. s²U is a universally conserved modification observed in all three domains of life at the wobble position in the anticodon. This modification enhances the stacking of s²U with U35 to stabilize the anticodon structure, facilitating codon–anticodon interactions. s²U is also recognized by multiple aminoacyl-tRNA synthetases for efficient amino acylation (Giegé and Eriani, 2023). Although elimination of this sulfur modification causes severe growth phenotypes in most organisms (Kambampati and Lauhon, 2003; Dewez et al., 2008), unexpectedly, the deletion of mnmA in Mtb did not attenuate growth of the pathogen in vitro, suggesting that the Mtb requirement for s²U modification differs from other organisms. The marked growth retardation of ΔmnmA strain within macrophages indicates the specific requirement of this modification within host cells. This modification may be necessary for maintaining general translation efficiency inside host cells and/or facilitate the expression of specific genes that are necessary for survival within macrophages. In fact, lack of mnmA is reported to sensitize E. coli to oxidative stress, raising the possibility that s²U modification promotes Mtb growth under oxidative stress elicited by the host.

Differential codon usage between house-keeping genes and virulence genes could contribute to the differential growth phenotypes observed in vitro and in vivo in the ΔmnmA mutant. Among multiple codons encoding Glu, Gln, and Lys, s²U modification dependent codon usage might be preferentially distributed in genes associated with intracellular growth. For example, Mtb has two tRNA isoacceptors, tRNA-Glu1(s²UUC) and tRNA-Glu2(CUC), to decipher two Glu codons, GAA and GAG. According to the wobble paring rule, GAA is only decoded by tRNA-Glu1(s²UUC), whereas GAG is decoded by both tRNA-Glu1(s²UUC) and tRNA-Glu2(CUC); that is, GAG can be deciphered by an s²U-independent tRNA. Thus, genes required for intracellular growth might be enriched with GAA, an s²U-dependent codon. Similar codon usage differences could be present in Gln and Lys codons deciphered by s²U-containing tRNAs.

In most organisms, s²U is further modified into derivatives containing an additional chemical moiety at position 5 (Björk and Hagervall, 2014; Karlsborn et al., 2014; Asano et al., 2018). However, Mtb does not contain apparent homologs of the tRNA modifying enzymes that introduce the additional modifications to s²U. Thus, Mtb may contain s²U or s²U derivatives synthesized by other types of enzymes. Additional analyses to elucidate the structures of modified s²U in Mtb are warranted.

Another unexpected finding is the presence of a TrmR homolog in Mtb. TrmR in B. subtilis is a methylase that converts a 5-hydroxyuridine (ho⁵U) into a 5-methoxyuridine (mo⁵U) at the tRNA wobble position, suggesting that Mtb possesses mo⁵U. However, Mtb does not encode an apparent homolog for the enzyme that mediates the RNA hydroxylation reaction that yields ho⁵U, the substrate for TrmR. Although we cannot rule out the possibility that Mtb TrmR methylates a substrate other than RNA, it is possible that Mtb utilizes an unknown hydroxylation pathway for synthesizing ho⁵U.

Our findings serve as a valuable starting point for the research community to continue characterizing the physiological roles and mechanisms of Mtb tRNA modifications. Since tRNA-seq offers an efficient and scalable platform for surveying changes in tRNA modifications, this approach will be valuable across growth conditions, and may be extended to growth inside host cells. Finally, further studies elucidating the mechanisms by which tRNA modifications facilitate Mtb growth in host cells should be valuable for designing new therapeutics for tuberculosis.

Materials and methods

Key resources table

Reagent type (species) or resource	Designation	Source or reference	Identifiers	Additional information
Gene (Mycobacterium tuberculosis)	mnmA	NA	Uniprot: Rv3024c; Refseq: NP_217540.1
Gene (Mycobacterium tuberculosis)	truB	NA	Uniprot: Rv2793c; Refseq: NP_217309.1
Strain, strain background (Mycobacterium tuberculosis)	H37Rv	PMID:9634230
Genetic reagent (Mycobacterium tuberculosis)	MtbΔmnmA::zeo	This paper		Mtb H37Rv strain lacking mnmA
Genetic reagent (Mycobacterium tuberculosis)	MtbΔtruB::zeo	This paper		Mtb H37Rv strain lacking truB
Recombinant DNA reagent	pNit-RecET SacBR	NA		For homologous recombination
Chemical compound	Iodoacetamide (IAA)	Sigma
Chemical compound	1-Cyclohexyl-(2-morpholinoethyl) carbodiimide (CMC)	Sigma

Bacterial strains and growth conditions

Request a detailed protocol

Mtb strains were grown from frozen stocks into Middlebrook 7H9 medium supplemented with 0.2% glycerol, 0.05% Tween-80, and ADC (5 g/l bovine serum albumin, 2 g/l dextrose, 3 μg/ml catalase). Cultures were incubated at 37°C. Strains were grown to mid-log phase for all experiments (OD₆₀₀ 0.4–0.6). Growth was measured on a BioTek plate reader for in vitro growth by measuring OD₆₀₀ every 24 hr.

Bacterial strain construction

Supplementary file 4 depicts the strains, plasmids, primers, and recombinant DNA used for this study. Plasmids were built by restriction digest of a parental vector and inserts were prepared by Gibson assembly (Gibson et al., 2009). Plasmids were isolated from E. coli and confirmed via Sanger sequencing carried out by Genewiz, LLC (Massachusetts, USA).

Deletion mutants

Request a detailed protocol

The knockout strain Mtb∆mnmA::zeo (zeocin) was built using double-stranded recombineering in the parental Mtb strain H37Rv. A linear dsDNA fragment was constructed using stitch polymerase chain reaction (PCR) with the primers listed in Supplementary file 4 which consisted of a 500-bp region upstream of mnmA (Rv3024c), 500-bp downstream region, and a lox-zeo-lox fragment. This cassette was transformed into an H37Rv recombineering strain as described (Murphy et al., 2015) and plated on 7H10 + zeocin plates.

Homology search

Request a detailed protocol

Local BLAST was performed to search for Mtb homologs of tRNA modifying enzymes. First, the uniport IDs of tRNA modifying enzymes were obtained from Modomics (Boccaletto et al., 2022), and 12 proteins were manually added to the list (Cho et al., 2023; Kimura et al., 2020; Sakai et al., 2019; Kimura et al., 2022; Takakura et al., 2019; Sakai et al., 2016), including Q47319/TapT, P24188/TrhO, P76403/TrhP, O32034/TrhP1, O32035/TrhP2, P36566/CmoM, and Q87K36/TrcP, O34614/MnmM, Q8N5C7/DTWD1, Q8NBA8/DTWD2, O32036/TrmR, Q9KV41/AcpA. Uniprot ID provides a fasta file of tRNA modifying enzymes from the Uniprot database (UniProt, 2023). A blast database file was generated by ‘makeblastdb’ script using a fasta file of Mtb proteins (H37Rv strain) retrieved from NCBI. Then, the homologs of tRNA modifying enzymes were searched against the Mtb protein database using the fasta file of tRNA modifying enzymes as a query. Output format is defined by the following script: -outfmt ‘7 qacc sacc stitle score qcovs evalue pident’ -evalue 1e-10. The output file was modified by excel (Supplementary file 4).

Phylogenetic analysis of Mtb tRNA modifying enzyme homologs

Request a detailed protocol

Local BLAST was conducted to search for homologs of Mtb tRNA modifying enzymes in 120 manually picked organisms across three domains of life. A custom database was generated by combining the fasta files of organisms’ proteins retrieved from NCBI. Homologs of 18 Mtb tRNA modifying enzymes were searched against the custom protein database using local blast. Log10 E-values were visualized by iTol (Letunic and Bork, 2021) with a phylogenetic tree generated by phyloT (PhyloT, 2015).

tRNA sequencing

Extraction of total RNA

Request a detailed protocol

Strains were grown to mid-log phase with the appropriate antibiotics and inducing agents described above. RNA was collected at the same OD₆₀₀ for each strain (between 0.4 and .6). Cells were left on ice for 20 min, then pelleted by centrifuging at 4000 rpm for 10 min at 4°C. Pellets were resuspended in 0.5–1 ml of TriZol (Life Technologies) and lysed using a BeadBug microtube homogenizer (Millipore Sigma). 200 μl of chloroform was added to each tube, after which samples obtained from Mtb strains were removed from biosafety level 3 precautions. Samples were centrifuged at 15,000 rpm for 15 min at 4°C and the aqueous layer was collected into a fresh tube. To the original tube, 250 μl of sodium acetate buffer (300 mM sodium acetate pH 5.2 and 10 mM ethylenediaminetetraacetic acid (EDTA) pH 8.0) was added, and samples were vortexed at 4°C for 5 min then centrifuged at 15,000 × g for 15 min at 4°C. The aqueous layer was added to the fresh sample-containing tubes. 400 μl chloroform was added, and tubes were briefly vortexed and then centrifuged at 15,000 rpm for 1 min at 4°C. The aqueous phase was collected into a fresh tube and RNA recovered by ethanol precipitation. RNA pellets were resuspended in 10 mM sodium acetate pH 5.2 and stored at −80°C until processed for sequencing. Total RNA samples were alkali treated prior to tRNA extraction to deacylate all tRNAs (1 hr at 37°C in 100 mM Tris–HCl pH 9.0).

Isolation of tRNA fraction

Request a detailed protocol

1–2 μg of total RNA was run on a 10% TBE-UREA gel (Thermo Fisher Scientific) at 250 V for 1 hr. Gels were stained with SYBR Gold (Thermo Fisher Scientific), and tRNA was excised. Excised gels containing tRNA fractions were mashed in RNAse-free tubes, and 300 μl elution buffer (300 mM NaOAc pH 5.5, 1 mM EDTA pH 8.0, 0.10% sodium dodecyl sulfate) was added to each tube. Samples were shaken on a thermoshaker (Eppendorf) for 1–4 hr at 37°C and supernatant was collected using an Ultrafree filter column (Millipore Sigma). tRNA was recovered by isopropanol precipitation.

tRNA dephosphorylation

Request a detailed protocol

tRNA was dephosphorylated using QuickCIP (New England BioLabs) according to manufacturer instructions, and tRNA was collected by phenol–chloroform extraction followed by isopropanol precipitation.

IAA treatment

Request a detailed protocol

IAA treatment was performed as described (Herzog et al., 2017). Briefly, 500 ng of total RNA is combined with 10 mM of IAA, 50 mM NaPO₄ pH 8.0, and 50% dimethyl sulfoxide (DMSO) in a final volume of 50 μl. Reactions were incubated at 50°C for 15 min and quenched with dithiothreitol (DTT).

CMC treatment

Request a detailed protocol

CMC treatment was carried out as described in ref Briefly, 2.5 μg tRNA fraction in 0.5 μl was mixed with 15 μl CMC-BEU buffer with or without CMC (0.34 M or 0 M CMC, 7 M urea, 4 mM EDTA, and 50 mM bicine pH 7.9) and incubated at 37°C for 20 min. Adding 100 μl CMC stop solution (0.3 M NaOAc pH 5.2 and 100 mM EDTA) quenched the reaction. RNA was desalted with PD-10 desalting column (Cytiva) and recovered by ethanol precipitation. RNA was dissolved in 40 μl of 50 mM sodium carbonate buffer (pH 10.4) and incubated at 37°C for 4 hr, followed by ethanol precipitation.

Adapter ligation

Request a detailed protocol

0.5 μL RNase inhibitor was added to 3.5 μl dephosphorylated tRNA (200–250 ng tRNA) and samples were boiled at 80°C for 2 min. Boiled tRNA was mixed with 12 μl PEG buffer mix (10 μl 50% PEG8000, 2 μl 10× buffer B0216S; New England Biolabs). 3 μl of 5′ adenylated linkers (Supplementary file 4) were added (33 pmol/μl) along with 1 μl T4 RNA ligase 2 truncated (New England BioLabs) and incubated at 25°C for 2.5 hr. Samples were recovered by isopropanol precipitation and run on a 10% TBE-Urea PAGE gel for 40 min at 250 V. Ligated products were recovered by gel excision as described above.

Reverse transcription

Request a detailed protocol

Identical quantities of samples with different adapter sequences were pooled for reverse transcription for a total of 200–250 ng tRNA. Reverse transcription was performed by combining 2.1 μl dephosphorylated tRNA with 100 mM Tris–HCl pH 7.5, 0.5 mM EDTA, 1.25 μM RT primer (Supplementary file 4), 450 mM NaCl, 5 mM MgCl₂, 5 mM DTT, 500 nM TGIRT (InGex), and 15% PEG8000 in a final volume of 9 μl. Samples were incubated at 25°C for 30 min, after which 1 μl 10 mM dNTPs (New England BioLabs) were added and reactions incubated at 60°C for 1 hr. 1.15 μl NaOH was added, and samples were boiled for 15 min and run on a 10% TBE Urea PAGE gel at 250 V for 1 hr. Reverse transcription products were excised from gels and cDNA recovered by isopropanol precipitation. Linear single-stranded cDNA was circularized using CircLigase II (Lucigen) in accordance with manufacturer instructions.

PCR of tRNA libraries

Request a detailed protocol

PCR reactions were set up using HF Phusion according to the manufacturer’s instructions using a universal reverse primer (Supplementary file 4) and a different index primer for each pool of samples. PCR reactions were aliquoted into 4 tubes and collected after 6, 8, 10, and 12 cycles. Samples were run on a Native TBE PAGE gel (Thermo Fisher Scientific) at 180 V for 50 min, and amplified products were cut from the same cycle for each sequencing run. Samples were recovered by gel excision and isopropanol precipitation.

Sequencing

Request a detailed protocol

Sequencing was performed on a MiSeq instrument (Illumina) using 150 bp single end reads with a version 3, 150 cycle kit.

Analysis

3′ linker sequences and two nucleotides at the 5′ end were trimmed using cutadapt and fastx-trimmer. Bowtie v1.2.2 was used with default settings to map reads to reference Mtb tRNA sequences retrieved from Mycobrowser (Kapopoulou et al., 2011; Supplementary file 5). Mpileup files were generated using samtools (samtools mpileup -I -A --ff 4 -x -B -q 0 -d 10000000). For analysis of termination frequencies, 5′ end termini of mapped reads were piled up using bedtools genomecov (option, -d -5 -ibam). The number of 5′ termini at each tRNA position was divided by the total number of mapped termini at that position plus all upstream (5′) positions.

Macrophage infection

Request a detailed protocol

Auto-luminescent wild-type and ΔmnmA Mtb strains grown to the same OD₆₀₀ were pelleted by centrifugation and prepared in RPMI media by soft spinning as described (Saito et al., 2017). Briefly, cells were washed, pelleted, resuspended, and centrifuged at 121 × g, with the top half of the centrifuged supernatant used. Suspensions were diluted to a multiplicity of infection of 2 bacteria per mouse bone marrow-derived macrophage by determining the OD₆₀₀. Macrophages were infected for 6 hr, followed by a phosphate-buffered saline wash and addition of RPMI with or without ATRA. ATRA promotes macrophage control of Mtb infection (Babunovic et al., 2022) and was used to assess strain survival in an increasingly restricted macrophage environment. Survival was measured by luminescence in a BioTek plate reader and normalized to luminescence reads at time 0.

Data availability

The sequencing data reported in this paper have been deposited in the NCBI Gene expression omnibus https://www.ncbi.nlm.nih.gov/geo/ (accession code, GSE240200).

The following data sets were generated

1. Tomasi FG
2. Kimura S
3. Rubin EJ
4. Waldor MK
(2023) NCBI Gene Expression Omnibus
ID GSE240200. Profiling tRNA modifications in Mycobacterium tuberculosis.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE240200

References

1. Altschul SF
2. Gish W
3. Miller W
4. Myers EW
5. Lipman DJ
(1990) Basic local alignment search tool
Journal of Molecular Biology 215:403–410.
https://doi.org/10.1016/S0022-2836(05)80360-2
- PubMed
- Google Scholar
1. Anderson J
2. Phan L
3. Cuesta R
4. Carlson BA
5. Pak M
6. Asano K
7. Björk GR
8. Tamame M
9. Hinnebusch AG
(1998) The essential Gcd10p-Gcd14p nuclear complex is required for 1-methyladenosine modification and maturation of initiator methionyl-tRNA
Genes & Development 12:3650–3662.
https://doi.org/10.1101/gad.12.23.3650
- PubMed
- Google Scholar
1. Asano K
2. Suzuki T
3. Saito A
4. Wei F-Y
5. Ikeuchi Y
6. Numata T
7. Tanaka R
8. Yamane Y
9. Yamamoto T
10. Goto T
11. Kishita Y
12. Murayama K
13. Ohtake A
14. Okazaki Y
15. Tomizawa K
16. Sakaguchi Y
17. Suzuki T
(2018) Metabolic and chemical regulation of tRNA modification associated with taurine deficiency and human disease
Nucleic Acids Research 46:1565–1583.
https://doi.org/10.1093/nar/gky068
- PubMed
- Google Scholar
1. Babunovic GH
2. DeJesus MA
3. Bosch B
4. Chase MR
5. Barbier T
6. Dickey AK
7. Bryson BD
8. Rock JM
9. Fortune SM
(2022) CRISPR Interference reveals that All-Trans-Retinoic Acid Promotes Macrophage Control of Mycobacterium tuberculosis by limiting Bacterial access to Cholesterol and Propionyl Coenzyme A
mBio 13:e0368321.
https://doi.org/10.1128/mbio.03683-21
- PubMed
- Google Scholar
1. Björk GR
2. Hagervall TG
(2014) Transfer RNA Modification: presence, synthesis, and function
EcoSal Plus 6:.
https://doi.org/10.1128/ecosalplus.ESP-0007-2013
- PubMed
- Google Scholar
1. Boccaletto P
2. Stefaniak F
3. Ray A
4. Cappannini A
5. Mukherjee S
6. Purta E
7. Kurkowska M
8. Shirvanizadeh N
9. Destefanis E
10. Groza P
11. Avşar G
12. Romitelli A
13. Pir P
14. Dassi E
15. Conticello SG
16. Aguilo F
17. Bujnicki JM
(2022) MODOMICS: a database of RNA modification pathways. 2021 update
Nucleic Acids Research 50:D231–D235.
https://doi.org/10.1093/nar/gkab1083
- Google Scholar
1. Bosch B
2. DeJesus MA
3. Poulton NC
4. Zhang W
5. Engelhart CA
6. Zaveri A
7. Lavalette S
8. Ruecker N
9. Trujillo C
10. Wallach JB
11. Li S
12. Ehrt S
13. Chait BT
14. Schnappinger D
15. Rock JM
(2021) Genome-wide gene expression tuning reveals diverse vulnerabilities of M. tuberculosis
Cell 184:4579–4592.
https://doi.org/10.1016/j.cell.2021.06.033
- PubMed
- Google Scholar
(2014) Pseudouridine profiling reveals regulated mRNA pseudouridylation in yeast and human cells
Nature 515:143–146.
https://doi.org/10.1038/nature13802
- PubMed
- Google Scholar
1. Chionh YH
2. McBee M
3. Babu IR
4. Hia F
5. Lin W
6. Zhao W
7. Cao J
8. Dziergowska A
9. Malkiewicz A
10. Begley TJ
11. Alonso S
12. Dedon PC
(2016) tRNA-mediated codon-biased translation in mycobacterial hypoxic persistence
Nature Communications 7:13302.
https://doi.org/10.1038/ncomms13302
- PubMed
- Google Scholar
1. Cho G
2. Lee J
3. Kim J
(2023) Identification of a novel 5-aminomethyl-2-thiouridine methyltransferase in tRNA modification
Nucleic Acids Research 51:1971–1983.
https://doi.org/10.1093/nar/gkad048
- PubMed
- Google Scholar
1. Dai Q
2. Zhang L-S
3. Sun H-L
4. Pajdzik K
5. Yang L
6. Ye C
7. Ju C-W
8. Liu S
9. Wang Y
10. Zheng Z
11. Zhang L
12. Harada BT
13. Dou X
14. Irkliyenko I
15. Feng X
16. Zhang W
17. Pan T
18. He C
(2023) Quantitative sequencing using BID-seq uncovers abundant pseudouridines in mammalian mRNA at base resolution
Nature Biotechnology 41:344–354.
https://doi.org/10.1038/s41587-022-01505-w
- PubMed
- Google Scholar
1. de Crécy-Lagard V
2. Jaroch M
(2021) Functions of Bacterial tRNA Modifications: from ubiquity to diversity
Trends in Microbiology 29:41–53.
https://doi.org/10.1016/j.tim.2020.06.010
- PubMed
- Google Scholar
1. DeJesus MA
2. Gerrick ER
3. Xu W
4. Park SW
5. Long JE
6. Boutte CC
7. Rubin EJ
8. Schnappinger D
9. Ehrt S
10. Fortune SM
11. Sassetti CM
12. Ioerger TR
(2017) Comprehensive essentiality analysis of the Mycobacterium tuberculosis Genome via saturating Transposon Mutagenesis
mBio 8:e02133-16.
https://doi.org/10.1128/mBio.02133-16
- PubMed
- Google Scholar
1. Delaunay S
2. Pascual G
3. Feng B
4. Klann K
5. Behm M
6. Hotz-Wagenblatt A
7. Richter K
8. Zaoui K
9. Herpel E
10. Münch C
11. Dietmann S
12. Hess J
13. Benitah SA
14. Frye M
(2022) Mitochondrial RNA modifications shape metabolic plasticity in metastasis
Nature 607:593–603.
https://doi.org/10.1038/s41586-022-04898-5
- PubMed
- Google Scholar
1. Dewez M
2. Bauer F
3. Dieu M
4. Raes M
5. Vandenhaute J
6. Hermand D
(2008) The conserved Wobble uridine tRNA thiolase Ctu1-Ctu2 is required to maintain genome integrity
PNAS 105:5459–5464.
https://doi.org/10.1073/pnas.0709404105
- PubMed
- Google Scholar
(2022) Transcriptome-wide mapping reveals a diverse dihydrouridine landscape including mRNA
PLOS Biology 20:e3001622.
https://doi.org/10.1371/journal.pbio.3001622
- PubMed
- Google Scholar
(2022) Sulfur availability impacts accumulation of the 2-Thiouridine tRNA modification in Bacillus subtilis
Journal of Bacteriology 204:e0000922.
https://doi.org/10.1128/jb.00009-22
- PubMed
- Google Scholar
1. El Yacoubi B
2. Lyons B
3. Cruz Y
4. Reddy R
5. Nordin B
6. Agnelli F
7. Williamson JR
8. Schimmel P
9. Swairjo MA
10. de Crécy-Lagard V
(2009) The universal YrdC/Sua5 family is required for the formation of threonylcarbamoyladenosine in tRNA
Nucleic Acids Research 37:2894–2909.
https://doi.org/10.1093/nar/gkp152
- PubMed
- Google Scholar
1. Finet O
2. Yague-Sanz C
3. Krüger LK
4. Tran P
5. Migeot V
6. Louski M
7. Nevers A
8. Rougemaille M
9. Sun J
10. Ernst FGM
11. Wacheul L
12. Wery M
13. Morillon A
14. Dedon P
15. Lafontaine DLJ
16. Hermand D
(2022) Transcription-wide mapping of dihydrouridine reveals that mRNA dihydrouridylation is required for meiotic chromosome segregation
Molecular Cell 82:404–419.
https://doi.org/10.1016/j.molcel.2021.11.003
- PubMed
- Google Scholar
1. Gibson DG
2. Young L
3. Chuang RY
4. Venter JC
5. Hutchison CA
6. Smith HO
(2009) Enzymatic assembly of DNA molecules up to several hundred kilobases
Nature Methods 6:343–345.
https://doi.org/10.1038/nmeth.1318
- PubMed
- Google Scholar
1. Giegé R
2. Eriani G
(2023) The tRNA identity landscape for aminoacylation and beyond
Nucleic Acids Research 51:1528–1570.
https://doi.org/10.1093/nar/gkad007
- PubMed
- Google Scholar
1. Goodarzi H
2. Nguyen HCB
3. Zhang S
4. Dill BD
5. Molina H
6. Tavazoie SF
(2016) Modulated expression of specific tRNAs drives Gene Expression and Cancer Progression
Cell 165:1416–1427.
https://doi.org/10.1016/j.cell.2016.05.046
- PubMed
- Google Scholar
(2010) Insights into the hyperthermostability and unusual region-specificity of archaeal Pyrococcus abyssi tRNA m1A57/58 methyltransferase
Nucleic Acids Research 38:6206–6218.
https://doi.org/10.1093/nar/gkq381
- PubMed
- Google Scholar
1. Herzog VA
2. Reichholf B
3. Neumann T
4. Rescheneder P
5. Bhat P
6. Burkard TR
7. Wlotzka W
8. von Haeseler A
9. Zuber J
10. Ameres SL
(2017) Thiol-linked alkylation of RNA to assess expression dynamics
Nature Methods 14:1198–1204.
https://doi.org/10.1038/nmeth.4435
- PubMed
- Google Scholar
1. Huang HY
2. Hopper AK
(2016) Multiple layers of stress-induced regulation in tRNA Biology
Life 6:16.
https://doi.org/10.3390/life6020016
- PubMed
- Google Scholar
Preprint
1. Jinich A
2. Zaveri A
3. DeJesus MA
4. Flores-Bautista E
5. Smith CM
6. Sassetti CM
7. Rock JM
8. Ehrt S
9. Schnappinger D
10. Ioerger TR
11. Rhee K
(2021) The Mycobacterium Tuberculosis Transposon Sequencing Database (MtbTnDB): A Large-Scale Guide to Genetic Conditional Essentiality
bioRxiv.
https://doi.org/10.1101/2021.03.05.434127
- Google Scholar
1. Jühling F
2. Mörl M
3. Hartmann RK
4. Sprinzl M
5. Stadler PF
6. Pütz J
(2009) tRNAdb 2009: compilation of tRNA sequences and tRNA genes
Nucleic Acids Research 37:D159–D162.
https://doi.org/10.1093/nar/gkn772
- PubMed
- Google Scholar
1. Kambampati R
2. Lauhon CT
(2003) MnmA and IscS are required for in vitro 2-thiouridine biosynthesis in Escherichia coli
Biochemistry 42:1109–1117.
https://doi.org/10.1021/bi026536+
- PubMed
- Google Scholar
(2011) The MycoBrowser portal: a comprehensive and manually annotated resource for mycobacterial genomes
Tuberculosis 91:8–13.
https://doi.org/10.1016/j.tube.2010.09.006
- PubMed
- Google Scholar
1. Karlsborn T
2. Tükenmez H
3. Mahmud AKMF
4. Xu F
5. Xu H
6. Byström AS
(2014) Elongator, a conserved complex required for wobble uridine modifications in eukaryotes
RNA Biology 11:1519–1528.
https://doi.org/10.4161/15476286.2014.992276
- PubMed
- Google Scholar
(2010) Detection of RNA modifications
RNA Biology 7:237–247.
https://doi.org/10.4161/rna.7.2.11468
- PubMed
- Google Scholar
1. Kimura S
2. Waldor MK
(2019) The RNA degradosome promotes tRNA quality control through clearance of hypomodified tRNA
PNAS 116:1394–1403.
https://doi.org/10.1073/pnas.1814130116
- PubMed
- Google Scholar
(2020) Comparative tRNA sequencing and RNA mass spectrometry for surveying tRNA modifications
Nature Chemical Biology 16:964–972.
https://doi.org/10.1038/s41589-020-0558-1
- PubMed
- Google Scholar
(2022) Sequential action of a tRNA base editor in conversion of cytidine to pseudouridine
Nature Communications 13:5994.
https://doi.org/10.1038/s41467-022-33714-x
- PubMed
- Google Scholar
1. Laxman S
2. Sutter BM
3. Wu X
4. Kumar S
5. Guo X
6. Trudgian DC
7. Mirzaei H
8. Tu BP
(2013) Sulfur Amino Acids regulate translational capacity and metabolic Homeostasis through Modulation of tRNA Thiolation
Cell 154:416–429.
https://doi.org/10.1016/j.cell.2013.06.043
- Google Scholar
1. Letunic I
2. Bork P
(2021) Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation
Nucleic Acids Research 49:W293–W296.
https://doi.org/10.1093/nar/gkab301
- PubMed
- Google Scholar
1. Liu Y
2. Zhou J
3. Li X
4. Zhang X
5. Shi J
6. Wang X
7. Li H
8. Miao S
9. Chen H
10. He X
11. Dong L
12. Lee GR
13. Zheng J
14. Liu R-J
15. Su B
16. Ye Y
17. Flavell RA
18. Yi C
19. Wu Y
20. Li H-B
(2022) tRNA-m¹A modification promotes T cell expansion via efficient MYC protein synthesis
Nature Immunology 23:1433–1444.
https://doi.org/10.1038/s41590-022-01301-3
- PubMed
- Google Scholar
1. Marchand V
(2021) Mapping of 7-methylguanosine (m(7)G), 3-methylcytidine (m(3)C), dihydrouridine (D) and 5-hydroxycytidine (ho(5)C) RNA modifications by AlkAniline-Seq
Methods Enzymol 658:25–47.
https://doi.org/10.1016/bs.mie.2021.06.001
- Google Scholar
1. Masuda I
2. Yamaki Y
3. Detroja R
4. Tagore S
5. Moore H
6. Maharjan S
7. Nakano Y
8. Christian T
9. Matsubara R
10. Lowe TM
11. Frenkel-Morgenstern M
12. Hou YM
(2022) tRNA methylation resolves codon usage bias at the limit of cell viability
Cell Reports 41:111539.
https://doi.org/10.1016/j.celrep.2022.111539
- PubMed
- Google Scholar
1. Motorin Y
2. Helm M
(2019) Methods for RNA modification mapping using deep sequencing: established and New Emerging Technologies
Genes 10:35.
https://doi.org/10.3390/genes10010035
- PubMed
- Google Scholar
(2015) Mycobacterial recombineering
Methods in Molecular Biology 1285:177–199.
https://doi.org/10.1007/978-1-4939-2450-9_10
- PubMed
- Google Scholar
1. Nedialkova DD
2. Leidel SA
(2015) Optimization of Codon translation rates via tRNA modifications maintains proteome integrity
Cell 161:1606–1618.
https://doi.org/10.1016/j.cell.2015.05.022
- PubMed
- Google Scholar
(2022) tRNA dysregulation and disease
Nature Reviews. Genetics 23:651–664.
https://doi.org/10.1038/s41576-022-00501-9
- PubMed
- Google Scholar
Software
1. PhyloT
(2015) A tree generator, version V2
PhyloT.

https://phylot.biobyte.de/
1. Saito K
2. Warrier T
3. Somersan-Karakaya S
4. Kaminski L
5. Mi J
6. Jiang X
7. Park S
8. Shigyo K
9. Gold B
10. Roberts J
11. Weber E
12. Jacobs WR
13. Nathan CF
(2017) Rifamycin action on RNA polymerase in antibiotic-tolerant Mycobacterium tuberculosis results in differentially detectable populations
PNAS 114:E4832–E4840.
https://doi.org/10.1073/pnas.1705385114
- PubMed
- Google Scholar
1. Sakai Y
2. Miyauchi K
3. Kimura S
4. Suzuki T
(2016) Biogenesis and growth phase-dependent alteration of 5-methoxycarbonylmethoxyuridine in tRNA anticodons
Nucleic Acids Research 44:509–523.
https://doi.org/10.1093/nar/gkv1470
- PubMed
- Google Scholar
(2019) Dual pathways of tRNA hydroxylation ensure efficient translation by expanding decoding capability
Nature Communications 10:2858.
https://doi.org/10.1038/s41467-019-10750-8
- PubMed
- Google Scholar
1. Schwartz MH
2. Wang H
3. Pan JN
4. Clark WC
5. Cui S
6. Eckwahl MJ
7. Pan DW
8. Parisien M
9. Owens SM
10. Cheng BL
11. Martinez K
12. Xu J
13. Chang EB
14. Pan T
15. Eren AM
(2018) Microbiome characterization by high-throughput transfer RNA sequencing and modification analysis
Nature Communications 9:5353.
https://doi.org/10.1038/s41467-018-07675-z
- PubMed
- Google Scholar
1. Shepherd J
2. Ibba M
(2015) Bacterial transfer RNAs
FEMS Microbiology Reviews 39:280–300.
https://doi.org/10.1093/femsre/fuv004
- PubMed
- Google Scholar
1. Soma A
2. Ikeuchi Y
3. Kanemasa S
4. Kobayashi K
5. Ogasawara N
6. Ote T
7. Kato J
8. Watanabe K
9. Sekine Y
10. Suzuki T
(2003) An RNA-modifying enzyme that governs both the codon and amino acid specificities of isoleucine tRNA
Molecular Cell 12:689–698.
https://doi.org/10.1016/s1097-2765(03)00346-0
- PubMed
- Google Scholar
1. Suzuki T
(2021) The expanding world of tRNA modifications and their disease relevance
Nature Reviews. Molecular Cell Biology 22:375–392.
https://doi.org/10.1038/s41580-021-00342-0
- PubMed
- Google Scholar
(2019) Biogenesis and functions of aminocarboxypropyluridine in tRNA
Nature Communications 10:5542.
https://doi.org/10.1038/s41467-019-13525-3
- PubMed
- Google Scholar
(2018) Cells alter their tRNA abundance to selectively regulate protein synthesis during stress conditions
Science Signaling 11:eaat6409.
https://doi.org/10.1126/scisignal.aat6409
- PubMed
- Google Scholar
1. UniProt C
(2023) UniProt: the universal protein knowledgebase in 2023
Nucleic Acids Research 51:D523–D531.
https://doi.org/10.1093/nar/gkac1052
- Google Scholar
1. van Kessel JC
2. Hatfull GF
(2007) Recombineering in Mycobacterium tuberculosis
Nature Methods 4:147–152.
https://doi.org/10.1038/nmeth996
- PubMed
- Google Scholar
(2004) Mycobacterium tuberculosis Rv2118c codes for a single-component homotetrameric m1A58 tRNA methyltransferase
Nucleic Acids Research 32:1018–1027.
https://doi.org/10.1093/nar/gkh207
- PubMed
- Google Scholar
Software
1. W.H.O
(2021)
Global tuberculosis report

World Health Organisation.
(2002) tadA, an essential tRNA-specific adenosine deaminase from Escherichia coli
The EMBO Journal 21:3841–3851.
https://doi.org/10.1093/emboj/cdf362
- PubMed
- Google Scholar
1. Zhang YJ
2. Reddy MC
3. Ioerger TR
4. Rothchild AC
5. Dartois V
6. Schuster BM
7. Trauner A
8. Wallis D
9. Galaviz S
10. Huttenhower C
11. Sacchettini JC
12. Behar SM
13. Rubin EJ
(2013) Tryptophan biosynthesis protects mycobacteria from CD4 T-cell-mediated killing
Cell 155:1296–1308.
https://doi.org/10.1016/j.cell.2013.10.045
- PubMed
- Google Scholar
1. Zhang W
2. Foo M
3. Eren AM
4. Pan T
(2022) tRNA modification dynamics from individual organisms to metaepitranscriptomics of microbiomes
Molecular Cell 82:891–906.
https://doi.org/10.1016/j.molcel.2021.12.007
- PubMed
- Google Scholar
1. Zheng G
2. Qin Y
3. Clark WC
4. Dai Q
5. Yi C
6. He C
7. Lambowitz AM
8. Pan T
(2015) Efficient and quantitative high-throughput tRNA sequencing
Nature Methods 12:835–837.
https://doi.org/10.1038/nmeth.3478
- PubMed
- Google Scholar

Peer review

Reviewer #1 (Public Review):

Tomasi et al. performed a combination of bioinformatic, next-generation tRNA sequencing experiments to predict the set of tRNA modifications and their corresponding genes in the tRNAs of the pathogenic bacteria Mycobacterium tuberculosis. Long known to be important for translation accuracy and efficiency, tRNA modifications are now emerging as having regulatory roles. However, the basic knowledge of the position and nature of the modifications present in a given organism is very sparse beyond a handful of model organisms. Studies that can generate the tRNA modification maps in different organisms along the tree of life are good starting points for further studies. The focus here on a major human pathogen that is studied by a large community raises the general interest of the study. Finally, deletion of the gene mnmA responsible for the insertion of s2U at position 34 revealed defects in growth in macrophage but in test tubes suggesting regulatory roles that will warrant further studies. The conclusions of the paper are mostly supported by the data but the partial nature of the bioinformatic analysis and absence of Mass-Spectrometry data make it incomplete. The authors do not take advantage of the Mass spec data that is published for Mycobacterium bovis (PMID: 27834374) to discuss what they find.

https://doi.org/10.7554/eLife.87146.3.sa1

Reviewer #2 (Public Review):

In this study, Tomasi et al identify a series of tRNA modifying enzymes from Mtb, show their function in the relevant tRNA modifications and by using at least one deleted strain for MnmA, they show the relevance of tRNA modification in intra-host survival and postulate their potential role in pathogenesis.

Conceptually it is a wonderful study, given that tRNA modifications are so fundamental to all life forms, showing their role in Mtb growth in the host is significant. However, the authors have not thoroughly analyzed the phenotype. The growth defect aspect or impact on pathogenesis needs to be adequately addressed.

- The authors show that ΔmnmA grows equally well in the in vitro cultures as the WT. However, they show attenuated growth in the macrophages. Is it because Glu1_TTC and Gln1-TTG tRNAs are not the preferred tRNAs for incorporation of Glu and Gln, respectively? And for some reason, they get preferred over the alternate tRNAs during infection? What dictates this selectivity?

- As such the growth defect shown in macrophages would be more convincing if the authors also show the phenotype of complementation with WT mnmA.

An important consideration here is the universal nature of these modifications across the life forms. Any strategy to utilize these enzymes as the potential therapeutic candidate would have to factor in this important aspect.

https://doi.org/10.7554/eLife.87146.3.sa2

Reviewer #3 (Public Review):

The work presented in the manuscript tries to identify tRNA modifications present in Mycobacterium tuberculosis (Mtb) using reverse transcription-derived error signatures with tRNA-seq. The study identified enzyme homologs and correlates them with presence of respective tRNA modifications in Mtb. The study used several chemical treatments (IAA and alkali treatment) to further enhance the reverse transcription signals and confirms the presence of modifications in the bases. tRNA modifications by two enzymes TruB and MnmA were established by doing tRNA-seq of respective deletion mutants. Ultimately, authors show that MnmA-dependent tRNA modification is important for intracellular growth of Mtb. Overall, this report identifies multiple tRNA modifications and discuss their implication in Mtb infection.

https://doi.org/10.7554/eLife.87146.3.sa3

Author response

The following is the authors’ response to the original reviews.

We thank the editors and reviewers for their thoughtful consideration of our manuscript. Here, we addressed the reviewers’ points.

Reviewer #1 (Public Review):

Tomasi et al. performed a combination of bioinformatic, next-generation tRNA sequencing experiments to predict the set of tRNA modifications and their corresponding genes in the tRNAs of the pathogenic bacteria Mycobacterium tuberculosis. Long known to be important for translation accuracy and efficiency, tRNA modifications are now emerging as having regulatory roles. However, the basic knowledge of the position and nature of the modifications present in a given organism is very sparse beyond a handful of model organisms. Studies that can generate the tRNA modification maps in different organisms along the tree of life are good starting points for further studies. The focus here on a major human pathogen that is studied by a large community raises the general interest of the study. Finally, deletion of the gene mnmA responsible for the insertion of s2U at position 34 revealed defects in in growth in macrophage but in test tubes suggesting regulatory roles that will warrant further studies. The conclusions of the paper are mostly supported by the data but the partial nature of the bioinformatic analysis and absence of Mass-Spectrometry data make it incomplete. The authors do not take advantage of the Mass spec data that is published for Mycobacterium bovis (PMID: 27834374) to discuss what they find.

1. The authors say they took a list of proteins involved in tRNA modifications from Modomics and added manually a few but we do not know the exact set of proteins that were used to search the M. mycobacterium genome.

Thank you for pointing out this issue. We added the complete list of proteins used for the BLAST query as Supplemental Table 1.

1. The absence of mnmGE genes in TB suggested that the xcm5U derivatives are absent. These are present in M. bovis (PMID: 27834374). Are the MnmEG gene found in M. bovis? If yes, then the authors should perform a phylogenetic distribution analysis in the Mycobacterial clade to see when they disappeared. If they are not present in M. bovis then maybe a non-orthologous set of enzymes do the same reaction and then the authors really do not know what modification is present or not at U34 without LC-MS. The exact same argument can be given for the xmo5U derivatives that are also found in M.bovis but not predicted by the authors in M. tuberculosis.

The reviewer raises a valid point. In M. bovis mnm5U and cmo5U derivatives were observed in LC-MS analysis. However, we did not identify candidate genes known to be involved in the biogenesis of mnm5U and cmo5U in the Mycobacteriaceae, including M. bovis and Mtb, suggesting that if these modifications are indeed present, they are not synthesized through canonical biogenesis pathways in this family. There are several examples where the same modification is generated by distinct modification enzymes (Kimura, 2021). These observations raise the interesting possibility that in the Mycobacteriaceae and most species in actinomycetota (except for Bifidobacterium, Corynebacterium and Rhodococcus species), major wobble modifications are generated by biosynthesis pathways that are distinct from those employed by well-characterized organisms. Future studies will examine this hypothesis.

1. Why is the Psi32 predicted by the authors because of the presence of the Rv3300c/Psu9 gene not detected by CMC-treated tRNA seq while the other Psi residues are? Members of this family can modify both rRNA and tRNA. So the presence of the gene does not guarantee the presence of the modification in tRNAs

Thank you very much for the careful read. We did not include RluA in the list of query proteins because it is not classified as a tRNA modification enzyme in Modomics. Additionally, the CMC-coupled tRNA-seq is imperfect for detection of all pseudouridylated positions. Due to this limitation, we only assigned modifications that are both predicted by the presence of putative biosynthetic enzymes and RT-derived signatures. As the reviewer points out, we cannot rule out that this homolog targets only rRNAs. We clarified this possibility in the revised manuscript by adding the following sentence: “Additionally, CMC treatment may not identify Ψ at all positions, thus, the targets of Rv3300 and Rv1540 remain unclear. Since these genes are similar to E. coli rluA, which also targets rRNA, these genes may target rRNAs instead of tRNAs” (lines 298-300)

In the revised manuscript, RluA was added to the BLAST query for creating Fig. 2. Interestingly, Rv3300c is more similar to Pus9 than RluA, while Rv1540 is the Mtb gene most similar to E. coli RluA suggesting that these two genes encode pseudouridylases that target different species of tRNAs/rRNAs.

1. What are tsaBED not essential but tsaC (called sua5 by the authors) essential?

Thank you for pointing out this interesting observation. We are also curious about differences in the essentiality among t6A biogenesis genes. We speculate that TsaC has critical roles in cell viability other than t6A synthesis. TsaC synthesizes threonylcarbamoyl-AMP as an intermediate for t6A biogenesis. Thus, it is possible that this intermediate has a role in other essential cellular activities besides t6A biogenesis. Further study of these factors in Mtb could reveal interesting crosstalk between modification synthesis and other cellular activities.

Reviewer #2 (Public Review):

In this study, Tomasi et al identify a series of tRNA modifying enzymes from Mtb, show their function in the relevant tRNA modifications and by using at least one deleted strain for MnmA, they show the relevance of tRNA modification in intra-host survival and postulate their potential role in pathogenesis.

Conceptually it is a wonderful study, given that tRNA modifications are so fundamental to all life forms, showing their role in Mtb growth in the host is significant. However, the authors have not thoroughly analyzed the phenotype. The growth defect aspect or impact on pathogenesis needs to be adequately addressed.

The authors show that ΔmnmA grows equally well in the in vitro cultures as the WT. However, they show attenuated growth in the macrophages. Is it because Glu1_TTC and Gln1-TTG tRNAs are not the preferred tRNAs for incorporation of Glu and Gln, respectively? And for some reason, they get preferred over the alternate tRNAs during infection? What dictates this selectivity?

Thank you very much for raising this excellent point. As the reviewer suggests, the attenuation of ΔmnmA Mtb growth inside of macrophages could be caused by disparate codon usage between genes required for in vitro growth and intracellular growth. Among multiple codons encoding Glu, Gln, or Lys, s2U modification-dependent codons might be preferentially distributed in genes associated with intracellular growth. For example, Mtb has two tRNA isoacceptors, Glu1_TTC and Glu2_CTC, to decipher two Glu codons, GAA and GAG. According to the wobble pairing rule, GAA is only decoded by Glu1_TTC, whereas GAG is decoded by both Glu1_TTC and Glu2_CTC; i.e., GAG can be deciphered by an s2U-independent tRNA. Thus, genes required for intracellular growth might be enriched with GAA, an s2U-dependent codon. Similar codon usage differences could be present in Gln and Lys codons deciphered by s2U-containing tRNAs. In the revised manuscript, we included a new paragraph in the discussion explaining the possibility that differences in codon usage could contribute to the intracellular fitness defect of the ΔmnmA Mtb mutant (lines 323-332).

As such the growth defect shown in macrophages would be more convincing if the authors also show the phenotype of complementation with WT mnmA.

The reviewer raises a valid point. We note however, that Rv3023c, a putative transposase, is downstream of MnmA and unlike MnmA, Rv3023c appears to be dispensable for in vivo growth, according to the Tn-seq database (reference 44 and 45). Therefore, it is likely that the intracellular growth defect is caused by loss of mnmA.

An important consideration here is the universal nature of these modifications across the life forms. Any strategy to utilize these enzymes as the potential therapeutic candidate would have to factor in this important aspect.

This is a valid point. Targeting a pathogen-specific system enables avoidance of the adverse side effects caused by many therapeutic reagents. There are a couple of Mtb modification enzymes that are specific to bacteria and critical for Mtb fitness (e.g., TilS). These enzymes represent ideal potential therapeutic targets to impede Mtb intracellular growth.

Reviewer #3 (Public Review):

The work presented in the manuscript tries to identify tRNA modifications present in Mycobacterium tuberculosis (Mtb) using reverse transcription-derived error signatures with tRNA-seq. The study identified enzyme homologs and correlates them with presence of respective tRNA modifications in Mtb. The study used several chemical treatments (IAA and alkali treatment) to further enhance the reverse transcription signals and confirms the presence of modifications in the bases. tRNA modifications by two enzymes TruB and MnmA were established by doing tRNA-seq of respective deletion mutants. Ultimately, authors show that MnmA-dependent tRNA modification is important for intracellular growth of Mtb. Overall, this report identifies multiple tRNA modifications and discuss their implication in Mtb infection.

Important points to be considered:

The presence of tRNA-based modifications is well characterised across life forms including genus Mycobacterium (Mycobacterium tuberculosis: Varshney et al, NAR, 2004; Mycobacterium bovis: Chionh et al, Nat Commun, 2016; Mycobacterium abscessus: Thomas et al, NAR, 2020). These modifications are shown to be essential for pathogenesis of multiple organisms. A comparison of tRNA modification and their respective enzymes with host organism as well as other mycobacterium strains is required. This can be discussed in detail to understand the role of common as well as specific tRNA modifications implicated in pathogenesis.

The reviewer raises a fair point. However, with the exception of Chionh et al., the other studies cited here are not genome-wide characterization of tRNA modification. Re-analysis showed that the distribution of the tRNA modifying enzymes are very similar across mycobacterium strains, e.g., Mycobacterium smegmatis, Mycobacterium tuberculosis, and Mycobacterium abscessus, suggesting that modifications related to pathogenesis in Mtb may have different physiological roles in other Mycobacterium species. We included the distribution of tRNA modification enzymes across multiple mycobacterium species in a revised Fig. 1.

Authors state in line 293 "Several strong signatures were detected in Mtb tRNAs but not in E. coli". Authors can elaborate more on the unique features identified and their relevance in Mtb infection in the discussion or result section.

Thank you for the suggestion. However, the identity of these RT signatures and the relevance of these modifications for Mtb pathogenicity remains speculative at this point.

Deletion of MnmA is shown to be essential for E. coli growth under oxidative stress (Zhao et al, NAR, 2021). In similar lines, MnmA deleted Mtb suffers to grow in macrophage. Is oxidative stress in macrophage responsible for slow Mtb growth?

This is an excellent hypothesis which we have added to the revised manuscript (lines 320-322). “In fact, the absence of mnmA is reported to sensitize E. coli to oxidative stress, raising the possibility that s2U modification promotes Mtb growth under oxidative stress elicited by the host.”

Authors state in line 311-312 "Mtb does not contain apparent homologs of the tRNA modifying enzymes that introduce the additional modifications to s2U". This can be characterised further to rule out the possibility of other enzyme specifically employed by Mtb to introduce additional modification.

The reviewer raises a valid point. As discussed above (Reviewer #1, pt 2), Mtb may employ distinct enzymes to generate certain tRNA modifications. Future mass spec-based analyses of Mtb tRNAs will be carried out to identify the precise chemical structure of the sulfurated uridine, and subsequent studies will attempt to determine the enzymes that account for the biogenesis of these modifications.

Kimura, S. (2021). Distinct evolutionary pathways for the synthesis and function of tRNA modifications. Brief Funct Genomics, 20(2), 125-134. doi:10.1093/bfgp/elaa027

Reviewer #1 (Recommendations For The Authors):

Additional data and Analyses

The Modomics database is far from complete so it would be more rigorous to give the full set of genes that was used to do the searches as supplemental data.

Thank you for the suggestion. We added the list of the query genes as Supplemental Table 1.Minor points to be fixed

1. The authors name the psi32 synthase Rv3300c Pus9 when it is a member of the RluA family. It is not clear why the yeast/eukaryotic name was used.

We included enzymes from diverse species in our query, including eukaryotic genes. Indeed, we found that Rv3300c showed the lowest E-value among our query genes, therefore, we name Rv3300c as Pus9.

1. The sua5 gene name was used it should be tsaC to follow the accepted nomenclature.

We renamed Sua5 to TsaC2.

1. The statement lines 203-296 was totally unclear. I did not understand what the authors were trying to say at all.

This paragraph described how sequence context can result in different reverse transcription-derived signatures from dihydrouridine (D). We added a schematics describing this paragraph as Supplemental Fig. 6.

1. In reference, names with special characters should be fixed such as Börk.

We fixed the names with special characters.

Reviewer #2 (Recommendations For The Authors):

The authors state that at least some of tRNA modifying enzymes, while redundant for growth in vitro, may play a role during growth inside the macrophages, mostly due to the diverse stresses they could encounter.

We added a sentence, “In fact, the absence of mnmA is reported to sensitize E. coli to oxidative stress, suggesting that s2U modification is required for Mtb growth under oxidative stress elicited by the host” in the discussion.

• Ideally, authors could have tested the impact of diverse intracellular stresses that Mtb encounters, like redox stress, nitrosative, pH or nutritional stress, to check whether any of these stresses cause in vitro growth defects in ΔmnmA strain.

Thank you for the suggestion. This point will be addressed in future experiments.

This would be a wonderful way to show that under stress, the essentiality of tRNA modification enzymes changes.

Reviewer #3 (Recommendations For The Authors):

• In general, the clarity of the presentation can be improved.

• Authors state that "MiaA, is non-essential in E. coli, but apparently essential in Mtb". While MiaA is shown to be critical for the fitness and virulence of extraintestinal pathogenic E. coli (Fleming et al, NAR, 2022). This can be clarified.

We rephrase as follows: “Unexpectedly, one modifying enzyme, MiaA, is non-essential in E. coli grown in nutrient-rich medium, but apparently …”

• Line numbers 130-132 is a repetition of line numbers 103-105

We repeated these sentences because the same claim was deduced from different experiments, i.e., BLAST search and tRNAseq.

• Line number 228: The presence of U at position 55 in the tRNAs can be included in the text for a better understanding.

We changed the text as following: “… Termination signatures derived from position 55, which is exclusively uridine in all tRNA species, increased in most tRNA species, suggesting …”

• A detailed pictorial depiction on comparing the modifications and enzymes from E. coli and Mtb can be included for easy understanding.

We created an E. coli tRNA modification map in the same format as Figure 2C and added it to the revised manuscript as a new Supplementary Fig. 1.

https://doi.org/10.7554/eLife.87146.3.sa4

Article and author information

Author details

Francesca G Tomasi

Department of Immunology and Infectious Diseases Harvard T. H. Chan School of Public Health, Boston, United States

Contribution
Conceptualization, Data curation, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

Contributed equally with
Satoshi Kimura

Competing interests
No competing interests declared
Satoshi Kimura
1. Division of Infectious Diseases, Brigham and Women's Hospital, Boston, United States
2. Department of Microbiology, Harvard Medical School, Boston, United States
3. Howard Hughes Medical Institute, Boston, United States
Present address
Department of Microbiology and Immunology, Cornell University, Ithaca, United States

Contribution
Conceptualization, Data curation, Supervision, Investigation, Visualization, Methodology, Writing – review and editing

Contributed equally with
Francesca G Tomasi

For correspondence
sk2956@cornell.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-3555-5877
Eric J Rubin

Department of Immunology and Infectious Diseases Harvard T. H. Chan School of Public Health, Boston, United States

Contribution
Conceptualization, Supervision, Funding acquisition, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-5120-962X
Matthew K Waldor
1. Department of Immunology and Infectious Diseases Harvard T. H. Chan School of Public Health, Boston, United States
2. Division of Infectious Diseases, Brigham and Women's Hospital, Boston, United States
3. Department of Microbiology, Harvard Medical School, Boston, United States
4. Howard Hughes Medical Institute, Boston, United States
Contribution
Conceptualization, Supervision, Funding acquisition, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-1843-7000

Funding

National Institute of Allergy and Infectious Diseases (R01AI-042347)

Matthew K Waldor

National Institute of Allergy and Infectious Diseases (P01AI-095208)

Eric J Rubin

Harvard Medical School (Dean's Innovation Award)

Eric J Rubin

Howard Hughes Medical Institute (MKW)

Matthew K Waldor

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

We appreciate all the members of Waldor lab for fruitful discussion and comments on the manuscript. We also thank Gregory Babunovic for his valuable assistance setting up the macrophage infections used here. This work is supported by NIH/NIAID grants to MKW (R01AI-042347) and EJR (P01AI-095208), a Dean’s Innovation Award from Harvard Medical School to EJR and the Howard Hughes Medical Institute (HHMI) to MKW.

Senior Editor

Bavesh D Kana, University of the Witwatersrand, South Africa

Reviewing Editor

Rajan Sankaranarayanan, Centre for Cellular and Molecular Biology, India

Version history

Preprint posted: February 20, 2023 (view preprint)
Sent for peer review: February 27, 2023
Preprint posted: April 12, 2023 (view preprint)
Preprint posted: August 8, 2023 (view preprint)
Version of Record published: September 27, 2023 (version 1)

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.87146. This DOI represents all versions, and will always resolve to the latest one.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

516

Page views
74

Downloads
0

Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Francesca G Tomasi
Satoshi Kimura
Eric J Rubin
Matthew K Waldor

(2023)

A tRNA modification in Mycobacterium tuberculosis facilitates optimal intracellular growth

eLife 12:RP87146.

https://doi.org/10.7554/eLife.87146.3

Share this article

Cite this article

Phylogenetic distribution of Mtb tRNA modifying enzyme homologs.

Figure 1—source data 1

Heatmap of misincorporation and early termination frequency in sequencing of tRNAs from wild-ype Mtb.

Figure 2—source data 1

Iodoacetamide (IAA) treatment promotes detection of sulfur modifications on tRNAs by enhancing termination signals.

Figure 3—source data 1

Heatmap and plot of early termination frequency from sequencing of tRNAs from wild-type and MtbΔmnmA with and without RNA alkylation.

Figure 4—source data 1

1-Cyclohexyl-(2-morpholinoethyl) carbodiimide (CMC) and alkali treatment facilitate detection of ψ, D, and m7G modifications in E. coli.

Figure 5—source data 1

Heatmap of early termination frequency from sequencing of tRNAs isolated from wild-type (WT) and MtbΔtruB following 1-cyclohexyl-(2-morpholinoethyl) carbodiimide (CMC) treatment.

Figure 6—source data 1

MtbΔmnmA is attenuated in a macrophage infection model.

Figure 7—source data 1

Author details

Francesca G Tomasi

Contribution

Contributed equally with

Competing interests

Satoshi Kimura

Present address

Contribution

Contributed equally with

For correspondence

Competing interests

Eric J Rubin

Contribution

Competing interests

Matthew K Waldor

Contribution

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading

1-Cyclohexyl-(2-morpholinoethyl) carbodiimide (CMC) and alkali treatment facilitate detection of ψ, D, and m⁷G modifications in E. coli.