Transfer RNAs (tRNAs) were first discovered in the 1950s (Kresge et al., 2005) and have since been found to play a vital role in the central dogma of molecular biology in all living systems (Crick, 1970). During the 1960s, tRNAs were also reported in viruses of bacteria (bacteriophages or phages) (Weiss et al., 1968). We now know that phage-encoded tRNAs are widespread, especially among virulent phages (Bailly-Bechet et al., 2007). Multiple hypotheses have been proposed for the role of these phage-encoded tRNAs. The most established being codon compensation, where codons rarely used by the host but necessary to the phage are supplemented by the tRNAs encoded by the phage (Bailly-Bechet et al., 2007). Why phages are pushed toward these alternative codons is generally believed to be a side effect of differences in the GC content of phage and host (Bailly-Bechet et al., 2007; Lucks et al., 2008; Limor-Waisberg et al., 2011). A recent study by Yang et al., 2021 may have hinted at an additional factor: phage tRNAs represent a means to counteract the depletion of host tRNAs that occurs as an early response to phage infection. The host employs several mechanisms to deplete its tRNA pool, such as downregulating the expression of its tRNAs, modifying tRNAs to make them unusable for translation, and most notably cleaving the tRNAs using anticodon nucleases (Thompson and Parker, 2009; Yang et al., 2021; Bailly-Bechet et al., 2007; Amitsur et al., 1989; Wolfram-Schauerte et al., 2022). Exactly what activates host tRNA cleavage is often unknown; an exception is anticodon nuclease PrrC, which cleaves tRNA-Lys(ttt) when triggered after sensing phage T4-encoded protein Stp (Kaufmann, 2011). In response, phage T4 encodes a tRNA ligase that repairs the cleaved tRNA-Lys (Kaufmann, 2011). Recently, a phage-encoded tRNA was found to replenish the depleted host tRNA caused by phage defense system Retron Ec78 , thereby preventing the inhibition of phage propagation (Azam et al., 2023). However, it remains unclear how phage tRNAs avoid being degraded by the same mechanism that results in the depletion of host tRNAs during phage infection. Here, we have analyzed phage-encoded tRNAs and hypothesize that the tRNAs encoded by phages are insensitive to tRNA anticodon nuclease activity, preventing depletion of the tRNA pool and translation stalling during phage infection.

To investigate our hypothesis that phage tRNAs are insensitive to tRNA anticodon nucleases, we analyzed the tRNAs encoded by a large and well-characterized dataset of tRNA-rich bacteriophages (33 tRNAs per phage on average) that infect mycobacteria: mycobacteriophage cluster C1 (Russell and Hatfull, 2017; Figure 1A and B). We specifically selected this mycobacteriophage dataset because the bacterial host encodes a range of well-characterized tRNA nucleases (tRNAses), such as VapC, MazF, and RelE (Winther et al., 2016; Chauhan et al., 2022; Cruz et al., 2015; Cintrón et al., 2019; Barth et al., 2021; Pedersen et al., 2003). A subset of these tRNAses target the tRNA anticodon loop and are activated upon various stress responses, including phage infection (Calcuttawala et al., 2022). Upon activation, these anticodon nucleases cleave specific tRNAs in conserved regions within the anticodon loop to inactivate the tRNAs and thereby regulate protein translation of the host (Winther et al., 2016). The cleavage region within the tRNA anticodon loop is sequence-dependent and highly specific to the type of tRNA. Mutations in the recognition and cleavage sequence within the anticodon loop have been found to result in insensitivity to these anticodon nucleases (Winther et al., 2016; Cruz et al., 2015). To check phage tRNAs for mutations that are known to cause insensitivity to anticodon nucleases (Winther et al., 2016; Cruz et al., 2015; Chauhan et al., 2022), we compared them to those of their host. We found that all 10 mycobacteriophage-encoded tRNAs that are targeted by anticodon nucleases contained anticodon loop mutations known to affect cleavage by VapC (Figure 1C, Supplementary file 1a, and Supplementary file 1b). These findings support the idea that phage-encoded tRNAs are insensitive to cleavage by anticodon nucleases (Figure 1C). We propose that these phage tRNAs serve as a mechanism to counteract the depletion of tRNAs by anticodon nucleases during phage infection, thereby allowing the phage to translate its proteins and successfully complete the infection cycle (Figure 1D).

Figure 1 with 1 supplement see all

Download asset Open asset

In addition, we observe that mycobacteriophages rarely encode tRNAs that are cleaved within the anticodon itself (Supplementary file 1b), suggesting that these anticodons are avoided by the phage and that no tRNAs evolved to be insensitive to cleavage. Specifically, this avoidance is seen for the serine-coding tRNAs that are cleaved at the GA site within the anticodon (tRNA-Ser(gga), tRNA-Ser(tga), tRNA-Ser(cga), and tRNA-Ser(aga)) (Winther et al., 2016; Supplementary file 1b). To compensate for this, the phage encodes a serine tRNA (tRNA-Ser(gct)) that is not targeted by nucleases. We observed the same avoidance for known targets of RelE (Pedersen et al., 2003), including stop codon (cta), tRNA-Leu(tta), and tRNA-Gln(cga) (Supplementary file 1b). Interestingly, when comparing the codon frequency of phage genes, we observed no differences in the codon frequency between codons for which the phage encodes a tRNA and those for which it does not (Welch’s two-sample t-test, t = 1.0471, df = 41.591, p-value=0.3011) (Supplementary file 1b). Moreover, we did not observe a difference in the codon frequency between phage and host genes (paired t-test, p-value=0.999) (Supplementary file 1b). We also found that in only 2 out of 23 instances the preferred codon (i.e. the most frequently encoded codon per amino acid) of the phage did not match that of the host (tRNA-Val(cac) and tRNA-Ala(cgc)) (Supplementary file 1b). Together, these observations suggest that the phage-encoded tRNAs were likely not selected for codon compensation. Overall, our findings support the hypothesis that phage tRNAs in mycobacteria evolved to resist anticodon nucleases in order to overcome host tRNA-depletion strategies that limit phage propagation. To investigate whether this hypothesis could extend more generally outside mycobacteria, we examined other species with anticodon nucleases, including Shigella flexneri (VapC-LT2) (Winther and Gerdes, 2011), Escherichia coli (VapC, PrrC, Colicin D, and Colicin E5) (Winther and Gerdes, 2011; Kaufmann, 2011; Amitsur et al., 1987; Ogawa et al., 2006), and Salmonella enterica (VapC/MvpT) (Winther and Gerdes, 2011). It is important to point out that the following analyses are based solely on the known cleavage site, and that the recognition sequence of these anticodon nucleases is unknown, thus limiting the ability to identify possible insensitivity-causing mutations in the recognition site of the anticodon loop. Despite this limitation, we observed mutations in the anticodon loop near or at the cleavage site of anticodon loop-targeted phage-encoded tRNAs in nearly all instances (9 out of 11) (Figure 1—figure supplement 1 and Supplementary file 1c). When the anticodon itself was the target of the anticodon nuclease, as is the case for VapC-LT2 in S. flexneri, we found that the phage generally avoided encoding the target tRNA (tRNA-fMet), except for two instances. One of these is the phage-encoded tRNA-Lys(ttt) in coliphages. In this instance, the phage encodes a tRNA-Lys ligase (rnl1 and rnl2) that repairs tRNA-Lys(ttt) after it has been cleaved by PrrC (Kaufmann, 2011). The other instance involves phage-encoded tRNAs targeted by E. coli Colicin E5. The cleavage activity of Colicin E5 depends on modifying the wobble position (Ogawa et al., 2006), which might be absent in these phage tRNAs possibly caused by mutations in the anticodon loop that we observed for three out of the four targeted tRNAs (Figure 1—figure supplement 1c). Alternatively, these mutations might affect the recognition by Colicin E5. In summary, our findings in mycobacteria are consistent with almost all (9 out of 11) currently known instances of targeted phage tRNAs in enterobacteria (Supplementary file 1c). The discrepancies in these cases can be accounted for by the presence of tRNA ligases, reliance on tRNA modifications of the cleavage site, or mutations within the uncharacterized recognition site. Thus, we speculate that our hypothesis may be extended beyond mycobacteria and enterobacteria, given the ubiquity of virus-encoded tRNAs and host tRNAses (Ogawa et al., 2006; Cavard and Lazdunski, 1979; Jones et al., 2017).

An overview of the analysed data supporting the findings of this study are available within the paper and in the Supplementary Data. All genomic sequences of the C1 mycobacteriophages were obtained from the publicly available actinobacteriophage database (PhagesDB; link: https://phagesdb.org/subclusters/C1/). Mycobacterium smegmatis MC2-155 (CP000480.1) and Mycobacterium tuberculosis H37Rv (NC_000962.3) were used as the reference for obtaining the host tRNA sequences.

1. 1. Amitsur M
   2. Levitz R
   3. Kaufmann G

   (1987) Bacteriophage T4 anticodon nuclease, polynucleotide kinase and RNA ligase reprocess the host lysine tRNA

   The EMBO Journal 6:2499–2503.

   https://doi.org/10.1002/j.1460-2075.1987.tb02532.x

   • PubMed
   • Google Scholar
2. 1. Amitsur M
   2. Morad I
   3. Kaufmann G

   (1989) In vitro reconstitution of Anticodon Nuclease from components encoded by Phage T4 and Escherichia coli Ctr5X

   The EMBO Journal 8:2411–2415.

   https://doi.org/10.1002/j.1460-2075.1989.tb08371.x

   • PubMed
   • Google Scholar
3. Preprint

   1. Azam AH
   2. Chihara K
   3. Kondo K
   4. Nakamura T
   5. Ojima S
   6. Tamura A
   7. Yamashita W
   8. Cui L
   9. Takahashi Y
   10. Watashi K
   11. Kiga K

   (2023) Viruses Encode TRNA and Anti-Retron to Evade Bacterial Immunity

   bioRxiv.

   https://doi.org/10.1101/2023.03.15.532788

   • Google Scholar
4. 1. Bailly-Bechet M
   2. Vergassola M
   3. Rocha E

   (2007) Causes for the intriguing presence of tRNAs in Phages

   Genome Research 17:1486–1495.

   https://doi.org/10.1101/gr.6649807

   • PubMed
   • Google Scholar
5. 1. Barth VC
   2. Chauhan U
   3. Zeng J
   4. Su X
   5. Zheng H
   6. Husson RN
   7. Woychik NA

   (2021) Mycobacterium tuberculosis Vapc4 toxin engages small Orfs to initiate an integrated oxidative and copper stress response

   PNAS 118:e2022136118.

   https://doi.org/10.1073/pnas.2022136118

   • PubMed
   • Google Scholar
6. 1. Calcuttawala F
   2. Shaw R
   3. Sarbajna A
   4. Dutta M
   5. Sinha S
   6. K. Das Gupta S
   7. Minnick MF

   (2022) Apoptosis like symptoms associated with abortive infection of Mycobacterium Smegmatis by Mycobacteriophage D29

   PLOS ONE 17:e0259480.

   https://doi.org/10.1371/journal.pone.0259480

   • Google Scholar
7. 1. Cavard D
   2. Lazdunski C

   (1979) Interaction of Colicin E4 with specific receptor sites mediates its cleavage into two fragments inactive towards whole cells

   European Journal of Biochemistry 96:525–533.

   https://doi.org/10.1111/j.1432-1033.1979.tb13066.x

   • PubMed
   • Google Scholar
8. Book

   1. Chan PP
   2. Lowe TM

   (2019) tRNAscan-SE: searching for tRNA genes in Genomic sequences

   In: Chan PP, editors. Gene Prediction. Humana. pp. 1–14.

   https://doi.org/10.1007/978-1-4939-9173-0

   • Google Scholar
9. 1. Chauhan U
   2. Barth VC
   3. Woychik NA

   (2022) tRNAfMet Inactivating Mycobacterium tuberculosis Vapbc toxin-Antitoxin systems as therapeutic targets

   Antimicrobial Agents and Chemotherapy 66:e0189621.

   https://doi.org/10.1128/aac.01896-21

   • PubMed
   • Google Scholar
10. 1. Cintrón M
    2. Zeng JM
    3. Barth VC
    4. Cruz JW
    5. Husson RN
    6. Woychik NA

    (2019) Accurate target identification for Mycobacterium tuberculosis Endoribonuclease toxins requires expression in their native host

    Scientific Reports 9:5949.

    https://doi.org/10.1038/s41598-019-41548-9

    • PubMed
    • Google Scholar
11. 1. Crick F

    (1970) Molecular biology in the year 2000

    Nature 228:613–615.

    https://doi.org/10.1038/228613a0

    • PubMed
    • Google Scholar
12. 1. Cruz JW
    2. Sharp JD
    3. Hoffer ED
    4. Maehigashi T
    5. Vvedenskaya IO
    6. Konkimalla A
    7. Husson RN
    8. Nickels BE
    9. Dunham CM
    10. Woychik NA

    (2015) Growth-regulating Mycobacterium tuberculosis Vapc-Mt4 toxin is an Isoacceptor-specific tRNase

    Nature Communications 6:7480.

    https://doi.org/10.1038/ncomms8480

    • PubMed
    • Google Scholar
13. 1. Davidov E
    2. Kaufmann G

    (2008) Rloc: a wobble nucleotide-Excising and zinc-responsive bacterial tRNase

    Molecular Microbiology 69:1560–1574.

    https://doi.org/10.1111/j.1365-2958.2008.06387.x

    • PubMed
    • Google Scholar
14. 1. Jones AM
    2. Garza-Sánchez F
    3. So J
    4. Hayes CS
    5. Low DA

    (2017) Activation of contact-dependent Antibacterial tRNase toxins by translation elongation factors

    PNAS 114:E1951–E1957.

    https://doi.org/10.1073/pnas.1619273114

    • Google Scholar
15. Book

    1. Kaufmann G

    (2011) Rloc: A translation-disabling tRNase implicated in Phage exclusion during recovery from DNA damage

    In: Kaufmann G, editors. DNA Repair. Wiley. pp. 1–21.

    https://doi.org/10.5772/1751

    • Google Scholar
16. 1. Kresge N
    2. Simoni RD
    3. Hill RL

    (2005) The discovery of tRNA by Paul C

    Journal of Biological Chemistry 280:e37–e39.

    https://doi.org/10.1016/S0021-9258(20)79029-0

    • Google Scholar
17. 1. Laslett D
    2. Canback B

    (2004) ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences

    Nucleic Acids Research 32:11–16.

    https://doi.org/10.1093/nar/gkh152

    • PubMed
    • Google Scholar
18. 1. Limor-Waisberg K
    2. Carmi A
    3. Scherz A
    4. Pilpel Y
    5. Furman I

    (2011) Specialization versus adaptation: two strategies employed by Cyanophages to enhance their translation efficiencies

    Nucleic Acids Research 39:6016–6028.

    https://doi.org/10.1093/nar/gkr169

    • PubMed
    • Google Scholar
19. 1. Lin BY
    2. Chan PP
    3. Lowe TM

    (2019) tRNAviz: explore and visualize tRNA sequence features

    Nucleic Acids Research 47:W542–W547.

    https://doi.org/10.1093/nar/gkz438

    • PubMed
    • Google Scholar
20. 1. Lopes APY
    2. Lopes LM
    3. Fraga TR
    4. Chura-Chambi RM
    5. Sanson AL
    6. Cheng E
    7. Nakajima E
    8. Morganti L
    9. Martins EAL
    10. Hoskisson P

    (2014) Vapc from the Leptospiral Vapbc toxin-Antitoxin Module displays Ribonuclease activity on the initiator tRNA

    PLOS ONE 9:e101678.

    https://doi.org/10.1371/journal.pone.0101678

    • PubMed
    • Google Scholar
21. 1. Lucks JB
    2. Nelson DR
    3. Kudla GR
    4. Plotkin JB

    (2008) Genome landscapes and Bacteriophage Codon usage

    PLOS Computational Biology 4:e1000001.

    https://doi.org/10.1371/journal.pcbi.1000001

    • PubMed
    • Google Scholar
22. 1. Miyamoto T
    2. Ota Y
    3. Yokota A
    4. Suyama T
    5. Tsuneda S
    6. Noda N

    (2017) Characterization of a Deinococcus radiodurans Mazf: A UACA-specific RNA endoribonuclease

    MicrobiologyOpen 6:e00501.

    https://doi.org/10.1002/mbo3.501

    • PubMed
    • Google Scholar
23. 1. Ogawa T
    2. Inoue S
    3. Yajima S
    4. Hidaka M
    5. Masaki H

    (2006) Sequence-specific recognition of Colicin E5, a tRNA-targeting Ribonuclease

    Nucleic Acids Research 34:6065–6073.

    https://doi.org/10.1093/nar/gkl629

    • PubMed
    • Google Scholar
24. 1. Pedersen K
    2. Zavialov AV
    3. Pavlov MY
    4. Elf J
    5. Gerdes K
    6. Ehrenberg M

    (2003) The bacterial toxin Rele displays Codon-specific cleavage of mRNAs in the Ribosomal A site

    Cell 112:131–140.

    https://doi.org/10.1016/s0092-8674(02)01248-5

    • PubMed
    • Google Scholar
25. 1. Rice P
    2. Longden I
    3. Bleasby A

    (2000) EMBOSS: the European molecular biology open software suite

    Trends in Genetics 16:276–277.

    https://doi.org/10.1016/s0168-9525(00)02024-2

    • PubMed
    • Google Scholar
26. 1. Russell DA
    2. Hatfull GF

    (2017) Phagesdb: the Actinobacteriophage database

    Bioinformatics 33:784–786.

    https://doi.org/10.1093/bioinformatics/btw711

    • PubMed
    • Google Scholar
27. 1. Thompson DM
    2. Parker R

    (2009) Stressing out over tRNA cleavage

    Cell 138:215–219.

    https://doi.org/10.1016/j.cell.2009.07.001

    • PubMed
    • Google Scholar
28. Preprint

    1. Tynecki P
    2. Guziński A
    3. Kazimierczak J
    4. Jadczuk M
    5. Dastych J
    6. Onisko A

    (2020) PhageAI - Bacteriophage Life Cycle Recognition with Machine Learning and Natural Language Processing

    bioRxiv.

    https://doi.org/10.1101/2020.07.11.198606

    • Google Scholar
29. 1. Weiss SB
    2. Hsu WT
    3. Foft JW
    4. Scherberg NH

    (1968) Transfer RNA coded by the T4 Bacteriophage genome

    PNAS 61:114–121.

    https://doi.org/10.1073/pnas.61.1.114

    • PubMed
    • Google Scholar
30. 1. Winther KS
    2. Gerdes K

    (2011) Enteric virulence associated protein Vapc inhibits translation by cleavage of initiator tRNA

    PNAS 108:7403–7407.

    https://doi.org/10.1073/pnas.1019587108

    • PubMed
    • Google Scholar
31. 1. Winther K
    2. Tree JJ
    3. Tollervey D
    4. Gerdes K

    (2016) Vapcs of Mycobacterium tuberculosis Cleave Rnas essential for translation

    Nucleic Acids Research 44:9860–9871.

    https://doi.org/10.1093/nar/gkw781

    • PubMed
    • Google Scholar
32. 1. Wolfram-Schauerte M
    2. Pozhydaieva N
    3. Viering M
    4. Glatter T
    5. Höfer K

    (2022) Integrated Omics reveal time-resolved insights into T4 Phage infection of E. coli on Proteome and Transcriptome levels

    Viruses 14:2502.

    https://doi.org/10.3390/v14112502

    • PubMed
    • Google Scholar
33. 1. Yang JY
    2. Fang W
    3. Miranda-Sanchez F
    4. Brown JM
    5. Kauffman KM
    6. Acevero CM
    7. Bartel DP
    8. Polz MF
    9. Kelly L

    (2021) Degradation of host Translational machinery drives tRNA acquisition in viruses

    Cell Systems 12:771–779.

    https://doi.org/10.1016/j.cels.2021.05.019

    • PubMed
    • Google Scholar

Author details

1. Daan F van den Berg

   1. Department of Bionanoscience, Delft University of Technology, Delft, Netherlands
   2. Kavli Institute of Nanoscience, Delft, Netherlands

   Contribution

   Conceptualization, Data curation, Software, Formal analysis, Investigation, Methodology, Writing - original draft

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-2217-4074

2. Baltus A van der Steen

   1. Department of Bionanoscience, Delft University of Technology, Delft, Netherlands
   2. Kavli Institute of Nanoscience, Delft, Netherlands

   Contribution

   Investigation

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-8193-4814

3. Ana Rita Costa

   1. Department of Bionanoscience, Delft University of Technology, Delft, Netherlands
   2. Kavli Institute of Nanoscience, Delft, Netherlands
   3. Fagenbank, Delft, Netherlands

   Contribution

   Supervision, Investigation

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-6749-6408

4. Stan JJ Brouns

   1. Department of Bionanoscience, Delft University of Technology, Delft, Netherlands
   2. Kavli Institute of Nanoscience, Delft, Netherlands
   3. Fagenbank, Delft, Netherlands

   Contribution

   Supervision, Funding acquisition, Investigation, Project administration, Writing - review and editing

   For correspondence

   stanbrouns@gmail.com

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-9573-1724

• 3,093

  views

• 336

  downloads

• 17

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.