Functions of Bacterial tRNA Modifications: From Ubiquity to Diversity

Analysis of bacterial transcriptome and epitranscriptome using nanopore direct RNA sequencing

Bacterial gene expression is a complex process involving extensive regulatory mechanisms. Along with growing interests in this field, Nanopore Direct RNA Sequencing (DRS) provides a promising platform for rapid and comprehensive characterization of bacterial RNA biology. However, the DRS of bacterial RNA is currently deficient in the yield of mRNA-mapping reads and has yet to be exploited for transcriptome-wide RNA modification mapping. Here, we showed that pre-processing of bacterial total RNA (size selection followed by ribosomal RNA depletion and polyadenylation) guaranteed high throughputs of sequencing data and considerably increased the amount of mRNA reads. This way, complex transcriptome architectures were reconstructed for Escherichia coli and Staphylococcus aureus and extended the boundaries of 225 known E. coli operons and 89 defined S. aureus operons. Utilizing unmodified in vitro-transcribed (IVT) RNA libraries as a negative control, several Nanopore-based computational tools globally detected putative modification sites in the E. coli and S. aureus transcriptomes. Combined with Next-Generation Sequencing-based N6-methyladenosine (m6A) detection methods, 75 high-confidence m6A candidates were identified in the E. coli protein-coding transcripts, while none were detected in S. aureus. Altogether, we demonstrated the potential of Nanopore DRS in systematic and convenient transcriptome and epitranscriptome analysis.

Development, validation and application of an LC-MS/MS method quantifying free forms of the micronutrients queuine and queuosine in human plasma using a surrogate matrix approach

Queuosine (Q) is a hypermodified 7-deaza-guanosine nucleoside exclusively synthesized by bacteria. This micronutrient and its respective nucleobase form queuine (q) are salvaged by humans either from gut microflora or digested food. Depletion of Q-tRNA in human or mouse cells causes protein misfolding that triggers endoplasmic reticular stress and the activation of the unfolded protein responses. In vivo, this reduces the neuronal architecture of the mouse brain affecting learning and memory. Herein, a sensitive method for quantifying free q and Q in human blood was developed, optimised and validated. After evaluating q/Q extraction efficiency in several different solid-phase sorbents, Bond Elut PBA (phenylboronic acid) cartridges were found to have the highest extraction recovery for q (82%) and Q (71%) from pooled human plasma. PBS with 4% BSA was used as surrogate matrix for method development and validation. An LC-MS/MS method was validated across the concentration range of 0.0003-1 µM for both q and Q, showing excellent linearity (r2 = 0.997 (q) and r2 = 0.998 (Q)), limit of quantification (0.0003 µM), accuracy (100.39-125.71%) and precision (CV% < 15.68%). In a sampling of healthy volunteers (n = 44), there was no significant difference in q levels between male (n = 14; mean = 0.0068 µM) and female (n = 30; mean = 0.0080 µM) participants (p = 0.50). Q was not detected in human plasma. This validated method can now be used to further substantiate the role of q/Q in nutrition, physiology and pathology.

Temperature-Dependent tRNA Modifications in Bacillales

Transfer RNA (tRNA) modifications are essential for the temperature adaptation of thermophilic and psychrophilic organisms as they control the rigidity and flexibility of transcripts. To further understand how specific tRNA modifications are adjusted to maintain functionality in response to temperature fluctuations, we investigated whether tRNA modifications represent an adaptation of bacteria to different growth temperatures (minimal, optimal, and maximal), focusing on closely related psychrophilic (P. halocryophilus and E. sibiricum), mesophilic (B. subtilis), and thermophilic (G. stearothermophilus) Bacillales. Utilizing an RNA sequencing approach combined with chemical pre-treatment of tRNA samples, we systematically profiled dihydrouridine (D), 4-thiouridine (s4U), 7-methyl-guanosine (m7G), and pseudouridine (Ψ) modifications at single-nucleotide resolution. Despite their close relationship, each bacterium exhibited a unique tRNA modification profile. Our findings revealed increased tRNA modifications in the thermophilic bacterium at its optimal growth temperature, particularly showing elevated levels of s4U8 and Ψ55 modifications compared to non-thermophilic bacteria, indicating a temperature-dependent regulation that may contribute to thermotolerance. Furthermore, we observed higher levels of D modifications in psychrophilic and mesophilic bacteria, indicating an adaptive strategy for cold environments by enhancing local flexibility in tRNAs. Our method demonstrated high effectiveness in identifying tRNA modifications compared to an established tool, highlighting its potential for precise tRNA profiling studies.

Impact of the chemical modification of tRNAs anticodon loop on the variability and evolution of codon usage in proteobacteria

Despite the highly conserved nature of the genetic code, the frequency of usage of each codon can vary significantly. The evolution of codon usage is shaped by two main evolutionary forces: mutational bias and selection pressures. These pressures can be driven by environmental factors, but also by the need for efficient translation, which depends heavily on the concentration of transfer RNAs (tRNAs) within the cell. The data presented here supports the proposal that tRNA modifications play a key role in shaping the overall preference of codon usage in proteobacteria. Interestingly, some codons, such as CGA and AGG (encoding arginine), exhibit a surprisingly low level of variation in their frequency of usage, even across genomes with differing GC content. These findings suggest that the evolution of GC content in proteobacterial genomes might be primarily driven by changes in the usage of a specific subset of codons, whose usage is itself influenced by tRNA modifications.

The modification landscape of Pseudomonas aeruginosa tRNAs

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

tRNA modification profiling reveals epitranscriptome regulatory networks in Pseudomonas aeruginosa

Transfer RNA (tRNA) modifications have emerged as critical posttranscriptional regulators of gene expression affecting diverse biological and disease processes. While there is extensive knowledge about the enzymes installing the dozens of post-transcriptional tRNA modifications - the tRNA epitranscriptome - very little is known about how metabolic, signaling, and other networks integrate to regulate tRNA modification levels. Here we took a comprehensive first step at understanding epitranscriptome regulatory networks by developing a high-throughput tRNA isolation and mass spectrometry-based modification profiling platform and applying it to a Pseudomonas aeruginosa transposon insertion mutant library comprising 5,746 strains. Analysis of >200,000 tRNA modification data points validated the annotations of predicted tRNA modification genes, uncovered novel tRNA-modifying enzymes, and revealed tRNA modification regulatory networks in P. aeruginosa. Platform adaptation for RNA-seq library preparation would complement epitranscriptome studies, while application to human cell and mouse tissue demonstrates its utility for biomarker and drug discovery and development.

RNA modifying enzymes shape tRNA biogenesis and function

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

A tRNA modification pattern that facilitates interpretation of the genetic code

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

Functional redundancy in tRNA dihydrouridylation

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

tRNA modification reprogramming contributes to artemisinin resistance in Plasmodium falciparum

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

Alternate routes to mnm5s2U synthesis in Gram-positive bacteria

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

Prophage terminase with tRNase activity sensitizes Salmonella enterica to oxidative stress

Phage viruses shape the evolution and virulence of their bacterial hosts. The Salmonella enterica genome encodes several stress-inducible prophages. The Gifsy-1 prophage terminase protein, whose canonical function is to process phage DNA for packaging in the virus head, unexpectedly acts as a transfer ribonuclease (tRNase) under oxidative stress, cleaving the anticodon loop of tRNALeu. The ensuing RNA fragmentation compromises bacterial translation, intracellular survival, and recovery from oxidative stress in the vertebrate host. S. enterica adapts to this transfer RNA (tRNA) fragmentation by transcribing the RNA repair Rtc system. The counterintuitive translational arrest provided by tRNA cleavage may subvert prophage mobilization and give the host an opportunity for repair as a way of maintaining bacterial genome integrity and ultimately survival in animals.

Bacterial nucleoid is a riddle wrapped in a mystery inside an enigma

Bacterial chromosome, the nucleoid, is traditionally modeled as a rosette of DNA mega-loops, organized around proteinaceous central scaffold by nucleoid-associated proteins (NAPs), and mixed with the cytoplasm by transcription and translation. Electron microscopy of fixed cells confirms dispersal of the cloud-like nucleoid within the ribosome-filled cytoplasm. Here, I discuss evidence that the nucleoid in live cells forms DNA phase separate from riboprotein phase, the "riboid." I argue that the nucleoid-riboid interphase, where DNA interacts with NAPs, transcribing RNA polymerases, nascent transcripts, and ssRNA chaperones, forms the transcription zone. An active part of phase separation, transcription zone enforces segregation of the centrally positioned information phase (the nucleoid) from the surrounding action phase (the riboid), where translation happens, protein accumulates, and metabolism occurs. I speculate that HU NAP mostly tiles up the nucleoid periphery-facilitating DNA mobility but also supporting transcription in the interphase. Besides extruding plectonemically supercoiled DNA mega-loops, condensins could compact them into solenoids of uniform rings, while HU could support rigidity and rotation of these DNA rings. The two-phase cytoplasm arrangement allows the bacterial cell to organize the central dogma activities, where (from the cell center to its periphery) DNA replicates and segregates, DNA is transcribed, nascent mRNA is handed over to ribosomes, mRNA is translated into proteins, and finally, the used mRNA is recycled into nucleotides at the inner membrane. The resulting information-action conveyor, with one activity naturally leading to the next one, explains the efficiency of prokaryotic cell design-even though its main intracellular transportation mode is free diffusion.

Beyond the Anticodon: tRNA Core Modifications and Their Impact on Structure, Translation and Stress Adaptation

Transfer RNAs (tRNAs) are heavily decorated with post-transcriptional chemical modifications. Approximately 100 different modifications have been identified in tRNAs, and each tRNA typically contains 5-15 modifications that are incorporated at specific sites along the tRNA sequence. These modifications may be classified into two groups according to their position in the three-dimensional tRNA structure, i.e., modifications in the tRNA core and modifications in the anticodon-loop (ACL) region. Since many modified nucleotides in the tRNA core are involved in the formation of tertiary interactions implicated in tRNA folding, these modifications are key to tRNA stability and resistance to RNA decay pathways. In comparison to the extensively studied ACL modifications, tRNA core modifications have generally received less attention, although they have been shown to play important roles beyond tRNA stability. Here, we review and place in perspective selected data on tRNA core modifications. We present their impact on tRNA structure and stability and report how these changes manifest themselves at the functional level in translation, fitness and stress adaptation.

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

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

The modification landscape of P. aeruginosa tRNAs

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

Anticodon stem-loop tRNA modifications influence codon decoding and frame maintenance during translation

RNAs are central to protein synthesis, with ribosomal RNA, transfer RNAs and messenger RNAs comprising the core components of the translation machinery. In addition to the four canonical bases (uracil, cytosine, adenine, and guanine) these RNAs contain an array of enzymatically incorporated chemical modifications. Transfer RNAs (tRNAs) are responsible for ferrying amino acids to the ribosome, and are among the most abundant and highly modified RNAs in the cell across all domains of life. On average, tRNA molecules contain 13 post-transcriptionally modified nucleosides that stabilize their structure and enhance function. There is an extensive chemical diversity of tRNA modifications, with over 90 distinct varieties of modifications reported within tRNA sequences. Some modifications are crucial for tRNAs to adopt their L-shaped tertiary structure, while others promote tRNA interactions with components of the protein synthesis machinery. In particular, modifications in the anticodon stem-loop (ASL), located near the site of tRNA:mRNA interaction, can play key roles in ensuring protein homeostasis and accurate translation. There is an abundance of evidence indicating the importance of ASL modifications for cellular health, and in vitro biochemical and biophysical studies suggest that individual ASL modifications can differentially influence discrete steps in the translation pathway. This review examines the molecular level consequences of tRNA ASL modifications in mRNA codon recognition and reading frame maintenance to ensure the rapid and accurate translation of proteins.

The enigmatic epitranscriptome of bacteriophages: putative RNA modifications in viral infections

RNA modifications play essential roles in modulating RNA function, stability, and fate across all kingdoms of life. The entirety of the RNA modifications within a cell is defined as the epitranscriptome. While eukaryotic RNA modifications are intensively studied, understanding bacterial RNA modifications remains limited, and knowledge about bacteriophage RNA modifications is almost nonexistent. In this review, we shed light on known mechanisms of bacterial RNA modifications and propose how this knowledge might be extended to bacteriophages. We build hypotheses on enzymes potentially responsible for regulating the epitranscriptome of bacteriophages and their host. This review highlights the exciting prospects of uncovering the unexplored field of bacteriophage epitranscriptomics and its potential role to shape bacteriophage-host interactions.

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

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

Quantitative analysis of tRNA abundance and modifications by nanopore RNA sequencing

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

Precision run-on sequencing (PRO-seq) for microbiome transcriptomics

Bacteria respond to environmental stimuli through precise regulation of transcription initiation and elongation. Bulk RNA sequencing primarily characterizes mature transcripts, so to identify actively transcribed loci we need to capture RNA polymerase (RNAP) complexed with nascent RNA. However, such capture methods have only previously been applied to culturable, genetically tractable organisms such as E. coli and B. subtilis. Here we apply precision run-on sequencing (PRO-seq) to profile nascent transcription in cultured E. coli and diverse uncultured bacteria. We demonstrate that PRO-seq can characterize the transcription of small, structured, or post-transcriptionally modified RNAs, which are often absent from bulk RNA-seq libraries. Applying PRO-seq to the human microbiome highlights taxon-specific RNAP pause motifs and pause-site distributions across non-coding RNA loci that reflect structure-coincident pausing. We also uncover concurrent transcription and cleavage of CRISPR guide RNAs and transfer RNAs. We demonstrate the utility of PRO-seq for exploring transcriptional dynamics in diverse microbial communities.

The specificity landscape of bacterial ribonuclease P

Developing quantitative models of substrate specificity for RNA processing enzymes is a key step toward understanding their biology and guiding applications in biotechnology and biomedicine. Optimally, models to predict relative rate constants for alternative substrates should integrate an understanding of structures of the enzyme bound to "fast" and "slow" substrates, large datasets of rate constants for alternative substrates, and transcriptomic data identifying in vivo processing sites. Such data are either available or emerging for bacterial ribonucleoprotein RNase P a widespread and essential tRNA 5' processing endonuclease, thus making it a valuable model system for investigating principles of biological specificity. Indeed, the well-established structure and kinetics of bacterial RNase P enabled the development of high throughput measurements of rate constants for tRNA variants and provided the necessary framework for quantitative specificity modeling. Several studies document the importance of conformational changes in the precursor tRNA substrate as well as the RNA and protein subunits of bacterial RNase P during binding, although the functional roles and dynamics are still being resolved. Recently, results from cryo-EM studies of E. coli RNase P with alternative precursor tRNAs are revealing prospective mechanistic relationships between conformational changes and substrate specificity. Yet, extensive uncharted territory remains, including leveraging these advances for drug discovery, achieving a complete accounting of RNase P substrates, and understanding how the cellular context contributes to RNA processing specificity in vivo.

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

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

Importance

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

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

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

Unrealized targets in the discovery of antibiotics for Gram-negative bacterial infections

Advances in areas that include genomics, systems biology, protein structure determination and artificial intelligence provide new opportunities for target-based antibacterial drug discovery. The selection of a 'good' new target for direct-acting antibacterial compounds is the first decision, for which multiple criteria must be explored, integrated and re-evaluated as drug discovery programmes progress. Criteria include essentiality of the target for bacterial survival, its conservation across different strains of the same species, bacterial species and growth conditions (which determines the spectrum of activity of a potential antibiotic) and the level of homology with human genes (which influences the potential for selective inhibition). Additionally, a bacterial target should have the potential to bind to drug-like molecules, and its subcellular location will govern the need for inhibitors to penetrate one or two bacterial membranes, which is a key challenge in targeting Gram-negative bacteria. The risk of the emergence of target-based drug resistance for drugs with single targets also requires consideration. This Review describes promising but as-yet-unrealized targets for antibacterial drugs against Gram-negative bacteria and examples of cognate inhibitors, and highlights lessons learned from past drug discovery programmes.

Integrative Approach to Probe Alternative Redox Mechanisms in RNA Modifications

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

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

IMPORTANCE

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

Different modification pathways for m1A58 incorporation in yeast elongator and initiator tRNAs

As essential components of the protein synthesis machinery, tRNAs undergo a tightly controlled biogenesis process, which include the incorporation of numerous posttranscriptional modifications. Defects in these tRNA maturation steps may lead to the degradation of hypomodified tRNAs by the rapid tRNA decay (RTD) and nuclear surveillance pathways. We previously identified m1A58 as a late modification introduced after modifications Ψ55 and T54 in yeast elongator tRNAPhe. However, previous reports suggested that m1A58 is introduced early during the tRNA modification process, in particular on primary transcripts of initiator tRNAiMet, which prevents its degradation by RNA decay pathways. Here, aiming to reconcile this apparent inconsistency on the temporality of m1A58 incorporation, we examined its introduction into yeast elongator and initiator tRNAs. We used specifically modified tRNAs to report on the molecular aspects controlling the Ψ55 → T54 → m1A58 modification circuit in elongator tRNAs. We also show that m1A58 is efficiently introduced on unmodified tRNAiMet, and does not depend on prior modifications. Finally, we show that m1A58 has major effects on the structural properties of initiator tRNAiMet, so that the tRNA elbow structure is only properly assembled when this modification is present. This observation provides a structural explanation for the degradation of hypomodified tRNAiMet lacking m1A58 by the nuclear surveillance and RTD pathways.

A tRNA modification in Mycobacterium tuberculosis facilitates optimal intracellular growth

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

RNA folding studies inside peptide-rich droplets reveal roles of modified nucleosides at the origin of life

Compartmentalization of RNA in biopolymer-rich membraneless organelles is now understood to be pervasive and critical for the function of extant biology and has been proposed as a prebiotically plausible way to accumulate RNA. However, compartment-RNA interactions that drive encapsulation have the potential to influence RNA structure and function in compartment- and RNA sequence-dependent ways. Here, we detail next-generation sequencing (NGS) experiments performed in membraneless compartments called complex coacervates to characterize the fold of many different transfer RNAs (tRNAs) simultaneously under the potentially denaturing conditions of these compartments. Notably, we find that natural modifications favor the native fold of tRNAs in these compartments. This suggests that covalent RNA modifications could have played a critical role in metabolic processes at the origin of life.

Discovering RNA modification enzymes using a comparative genomics approach

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

A tRNA modification in Mycobacterium tuberculosis facilitates optimal intracellular growth

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

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

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

A-to-I RNA Editing in Klebsiella pneumoniae Regulates Quorum Sensing and Affects Cell Growth and Virulence

Millions of adenosine (A) to inosine (I) RNA editing events are reported and well-studied in eukaryotes; however, many features and functions remain unclear in prokaryotes. By combining PacBio Sequel, Illumina whole-genome sequencing, and RNA Sequencing data of two Klebsiella pneumoniae strains with different virulence, a total of 13 RNA editing events are identified. The RNA editing event of badR is focused, which shows a significant difference in editing levels in the two K. pneumoniae strains and is predicted to be a transcription factor. A hard-coded Cys is mutated on DNA to simulate the effect of complete editing of badR. Transcriptome analysis identifies the cellular quorum sensing (QS) pathway as the most dramatic change, demonstrating the dynamic regulation of RNA editing on badR related to coordinated collective behavior. Indeed, a significant difference in autoinducer 2 activity and cell growth is detected when the cells reach the stationary phase. Additionally, the mutant strain shows significantly lower virulence than the WT strain in the Galleria mellonella infection model. Furthermore, RNA editing regulation of badR is highly conserved across K. pneumoniae strains. Overall, this work provides new insights into posttranscriptional regulation in bacteria.

Expanding the Epitranscriptomic RNA Sequencing and Modification Mapping Mass Spectrometry Toolbox with Field Asymmetric Waveform Ion Mobility and Electrochemical Elution Liquid Chromatography

Post-transcriptional modifications of RNA strongly influence the RNA structure and function. Recent advances in RNA sequencing and mass spectrometry (MS) methods have identified over 140 of these modifications on a wide variety of RNA species. Most next-generation sequencing approaches can only map one RNA modification at a time, and while MS can assign multiple modifications simultaneously in an unbiased manner, MS cannot accurately catalog and assign RNA modifications in complex biological samples due to limitations in the fragment length and coverage depth. Thus, a facile method to identify novel RNA modifications while simultaneously locating them in the context of their RNA sequences is still lacking. We combined two orthogonal modes of RNA ion separation before MS identification: high-field asymmetric ion mobility separation (FAIMS) and electrochemically modulated liquid chromatography (EMLC). FAIMS RNA MS increases both coverage and throughput, while EMLC LC-MS orthogonally separates RNA molecules of different lengths and charges. The combination of the two methods offers a broadly applicable platform to improve the length and depth of MS-based RNA sequencing while providing contextual access to the analysis of RNA modifications.

Ultradeep characterisation of translational sequence determinants refutes rare-codon hypothesis and unveils quadruplet base pairing of initiator tRNA and transcript

Translation is a key determinant of gene expression and an important biotechnological engineering target. In bacteria, 5'-untranslated region (5'-UTR) and coding sequence (CDS) are well-known mRNA parts controlling translation and thus cellular protein levels. However, the complex interaction of 5'-UTR and CDS has so far only been studied for few sequences leading to non-generalisable and partly contradictory conclusions. Herein, we systematically assess the dynamic translation from over 1.2 million 5'-UTR-CDS pairs in Escherichia coli to investigate their collective effect using a new method for ultradeep sequence-function mapping. This allows us to disentangle and precisely quantify effects of various sequence determinants of translation. We find that 5'-UTR and CDS individually account for 53% and 20% of variance in translation, respectively, and show conclusively that, contrary to a common hypothesis, tRNA abundance does not explain expression changes between CDSs with different synonymous codons. Moreover, the obtained large-scale data provide clear experimental evidence for a base-pairing interaction between initiator tRNA and mRNA beyond the anticodon-codon interaction, an effect that is often masked for individual sequences and therefore inaccessible to low-throughput approaches. Our study highlights the indispensability of ultradeep sequence-function mapping to accurately determine the contribution of parts and phenomena involved in gene regulation.

Adaptation to Overflow Metabolism by Mutations That Impair tRNA Modification in Experimentally Evolved Bacteria

When microbes grow in foreign nutritional environments, selection may enrich mutations in unexpected pathways connecting growth and homeostasis. An evolution experiment designed to identify beneficial mutations in Burkholderia cenocepacia captured six independent nonsynonymous substitutions in the essential gene tilS, which modifies tRNA^Ile2 by adding a lysine to the anticodon for faithful AUA recognition. Further, five additional mutants acquired mutations in tRNA^Ile2, which strongly suggests that disrupting the TilS-tRNA^Ile2 interaction was subject to strong positive selection. Mutated TilS incurred greatly reduced enzymatic function but retained capacity for tRNA^Ile2 binding. However, both mutant sets outcompeted the wild type by decreasing the lag phase duration by ~3.5 h. We hypothesized that lysine demand could underlie fitness in the experimental conditions. As predicted, supplemental lysine complemented the ancestral fitness deficit, but so did the additions of several other amino acids. Mutant fitness advantages were also specific to rapid growth on galactose using oxidative overflow metabolism that generates redox imbalance, not resources favoring more balanced metabolism. Remarkably, 13 tilS mutations also evolved in the long-term evolution experiment with Escherichia coli, including four fixed mutations. These results suggest that TilS or unknown binding partners contribute to improved growth under conditions of rapid sugar oxidation at the predicted expense of translational accuracy. IMPORTANCE There is growing evidence that the fundamental components of protein translation can play multiple roles in maintaining cellular homeostasis. Enzymes that interact with transfer RNAs not only ensure faithful decoding of the genetic code but also help signal the metabolic state by reacting to imbalances in essential building blocks like free amino acids and cofactors. Here, we present evidence of a secondary function for the essential enzyme TilS, whose only prior known function is to modify tRNA^Ile(CAU) to ensure accurate translation. Multiple nonsynonymous substitutions in tilS, as well as its cognate tRNA, were selected in evolution experiments favoring rapid, redox-imbalanced growth. These mutations alone decreased lag phase and created a competitive advantage, but at the expense of most primary enzyme function. These results imply that TilS interacts with other factors related to the timing of exponential growth and that tRNA-modifying enzymes may serve multiple roles in monitoring metabolic health.

tRNAs Are Stable After All: Pitfalls in Quantification of tRNA from Starved Escherichia coli Cultures Exposed by Validation of RNA Purification Methods

tRNAs and ribosomal RNAs are often considered stable RNAs. In contrast to this view, we recently proposed that tRNAs are degraded during amino acid starvation and drug-induced transcription inhibition. However, reevaluation of our experimental approach revealed that common RNA extraction methods suffer from alarming extraction and size biases that can lead to gross underestimation of RNA levels in starved Escherichia coli populations. Quantification of tRNAs suffers additional biases due to differing fractions of tRNAs with base modifications in growing versus starved bacteria. Applying an improved methodology, we measured tRNA levels after starvation for amino acids, glucose, phosphate, or ammonium and transcription inhibition by rifampicin. We report that tRNA levels remain largely unaffected in all tested conditions, including several days of starvation. This confirms that tRNAs are remarkably stable RNAs and serves as a cautionary tale about quantification of RNA from cells cultured outside the steady-state growth regime. rRNA, conversely, is extensively degraded during starvation. Thus, E. coli downregulates the translation machinery in response to starvation by reducing the ribosome pool through rRNA degradation, while a high concentration of tRNAs available to supply amino acids to the remaining ribosomes is maintained. IMPORTANCE We show that E. coli tRNAs are remarkably stable during several days of nutrient starvation, although rRNA is degraded extensively under these conditions. The levels of these two major RNA classes are considered to be strongly coregulated at the level of transcription. We demonstrate that E. coli can control the ratio of tRNAs per ribosome under starvation by means of differential degradation rates. The question of tRNA stability in stressed E. coli cells has become subject to debate. Our in-depth analysis of RNA quantification methods reveals hidden technical pitfalls at every step of the analysis, from RNA extraction to target detection and normalization. Most importantly, starved E. coli populations were more resilient to RNA extraction than unstarved populations. The current results underscore that the seemingly trivial task of quantifying an abundant RNA species is not straightforward for cells cultured outside the exponential growth regime.

RNA folding studies inside peptide-rich droplets reveal roles of modified nucleosides at the origin of life

Compartmentalization of RNA in biopolymer-rich membraneless organelles is now understood to be pervasive and critical for the function of extant biology and has been proposed as a prebiotically-plausible way to accumulate RNA. However, compartment-RNA interactions that drive encapsulation have the potential to influence RNA structure and function in compartment- and RNA sequence-dependent ways. Herein, we detail Next-Generation Sequencing (NGS) experiments performed for the first time in membraneless compartments called complex coacervates to characterize the fold of many different transfer RNAs (tRNAs) simultaneously under the potentially denaturing conditions of these compartments. Strikingly, we find that natural modifications favor the native fold of tRNAs in these compartments. This suggests that covalent RNA modifications could have played a critical role in metabolic processes at the origin of life.

One Sentence Summary

We demonstrate that RNA folds into native secondary and tertiary structures in protocell models and that this is favored by covalent modifications, which is critical for the origins of life.

Bacteriophage Infection of the Marine Bacterium Shewanella glacialimarina Induces Dynamic Changes in tRNA Modifications

Viruses are obligate intracellular parasites that, throughout evolution, have adapted numerous strategies to control the translation machinery, including the modulation of post-transcriptional modifications (PTMs) on transfer RNA (tRNA). PTMs are critical translation regulators used to further host immune responses as well as the expression of viral proteins. Yet, we lack critical insight into the temporal dynamics of infection-induced changes to the tRNA modification landscape (i.e., 'modificome'). In this study, we provide the first comprehensive quantitative characterization of the tRNA modificome in the marine bacterium Shewanella glacialimarina during Shewanella phage 1/4 infection. Specifically, we show that PTMs can be grouped into distinct categories based on modification level changes at various infection stages. Furthermore, we observe a preference for the UAC codon in viral transcripts expressed at the late stage of infection, which coincides with an increase in queuosine modification. Queuosine appears exclusively on tRNAs with GUN anticodons, suggesting a correlation between phage codon usage and PTM modification. Importantly, this work provides the basis for further studies into RNA-based regulatory mechanisms employed by bacteriophages to control the prokaryotic translation machinery.

Pseudomonas aeruginosa GidA modulates the expression of catalases at the posttranscriptional level and plays a role in virulence

Pseudomonas aeruginosa gidA, which encodes a putative tRNA-modifying enzyme, is associated with a variety of virulence phenotypes. Here, we demonstrated that P. aeruginosa gidA is responsible for the modifications of uridine in tRNAs in vivo. Loss of gidA was found to have no impact on the mRNA levels of katA and katB, but it decreased KatA and KatB protein levels, resulting in decreased total catalase activity and a hydrogen peroxide-sensitive phenotype. Furthermore, gidA was found to affect flagella-mediated motility and biofilm formation; and it was required for the full virulence of P. aeruginosa in both Caenorhabditis elegans and macrophage models. Together, these observations reveal the posttranscriptional impact of gidA on the oxidative stress response, highlight the complexity of catalase gene expression regulation, and further support the involvement of gidA in the virulence of P. aeruginosa.

The Influence of the Nucleotide Composition of Genes and Gene Regulatory Elements on the Efficiency of Protein Expression in Escherichia coli

Recombinant proteins expressed in Escherichia coli are widely used in biochemical research and industrial processes. At the same time, achieving higher protein expression levels and correct protein folding still remains the key problem, since optimization of nutrient media, growth conditions, and methods for induction of protein synthesis do not always lead to the desired result. Often, low protein expression is determined by the sequences of the expressed genes and their regulatory regions. The genetic code is degenerated; 18 out of 20 amino acids are encoded by more than one codon. Choosing between synonymous codons in the coding sequence can significantly affect the level of protein expression and protein folding due to the influence of the gene nucleotide composition on the probability of formation of secondary mRNA structures that affect the ribosome binding at the translation initiation phase, as well as the ribosome movement along the mRNA during elongation, which, in turn, influences the mRNA degradation and the folding of the nascent protein. The nucleotide composition of the mRNA untranslated regions, in particular the promoter and Shine-Dalgarno sequences, also affects the efficiency of mRNA transcription, translation, and degradation. In this review, we describe the genetic principles that determine the efficiency of protein production in Escherichia coli.

Metal retention and replacement in QueD2 protect queuosine-tRNA biosynthesis in metal-starved Acinetobacter baumannii

In response to bacterial infection, the vertebrate host employs the metal-sequestering protein calprotectin (CP) to withhold essential transition metals, notably Zn(II), to inhibit bacterial growth. Previous studies of the impact of CP-imposed transition-metal starvation in A. baumannii identified two enzymes in the de novo biosynthesis pathway of queuosine-transfer ribonucleic acid (Q-tRNA) that become cellularly abundant, one of which is QueD2, a 6-carboxy-5,6,7,8-tetrahydropterin (6-CPH4) synthase that catalyzes the initial, committed step of the pathway. Here, we show that CP strongly disrupts Q incorporation into tRNA. As such, we compare the AbQueD2 "low-zinc" paralog with a housekeeping, obligatory Zn(II)-dependent enzyme QueD. The crystallographic structure of Zn(II)-bound AbQueD2 reveals a distinct catalytic site coordination sphere and assembly state relative to QueD and possesses a dynamic loop, immediately adjacent to the catalytic site that coordinates a second Zn(II) in the structure. One of these loop-coordinating residues is an invariant Cys18, that protects QueD2 from dissociation of the catalytic Zn(II) while maintaining flux through the Q-tRNA biosynthesis pathway in cells. We propose a "metal retention" model where Cys18 introduces coordinative plasticity into the catalytic site which slows metal release, while also enhancing the metal promiscuity such that Fe(II) becomes an active cofactor. These studies reveal a complex, multipronged evolutionary adaptation to cellular Zn(II) limitation in a key Zn(II) metalloenzyme in an important human pathogen.

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

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

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

Impact of queuosine modification of endogenous E. coli tRNAs on sense codon reassignment

The extent to which alteration of endogenous tRNA modifications may be exploited to improve genetic code expansion efforts has not been broadly investigated. Modifications of tRNAs are strongly conserved evolutionarily, but the vast majority of E. coli tRNA modifications are not essential. We identified queuosine (Q), a non-essential, hypermodified guanosine nucleoside found in position 34 of the anticodons of four E. coli tRNAs as a modification that could potentially be utilized to improve sense codon reassignment. One suggested purpose of queuosine modification is to reduce the preference of tRNAs with guanosine (G) at position 34 of the anticodon for decoding cytosine (C) ending codons over uridine (U) ending codons. We hypothesized that introduced orthogonal translation machinery with adenine (A) at position 34 would reassign U-ending codons more effectively in queuosine-deficient E. coli. We evaluated the ability of introduced orthogonal tRNAs with AUN anticodons to reassign three of the four U-ending codons normally decoded by Q34 endogenous tRNAs: histidine CAU, asparagine AAU, and aspartic acid GAU in the presence and absence of queuosine modification. We found that sense codon reassignment efficiencies in queuosine-deficient strains are slightly improved at Asn AAU, equivalent at His CAU, and less efficient at Asp GAU codons. Utilization of orthogonal pair-directed sense codon reassignment to evaluate competition events that do not occur in the standard genetic code suggests that tRNAs with inosine (I, 6-deaminated A) at position 34 compete much more favorably against G34 tRNAs than Q34 tRNAs. Continued evaluation of sense codon reassignment following targeted alterations to endogenous tRNA modifications has the potential to shed new light on the web of interactions that combine to preserve the fidelity of the genetic code as well as identify opportunities for exploitation in systems with expanded genetic codes.

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

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

THUMPD1 bi-allelic variants cause loss of tRNA acetylation and a syndromic neurodevelopmental disorder

Covalent tRNA modifications play multi-faceted roles in tRNA stability, folding, and recognition, as well as the rate and fidelity of translation, and other cellular processes such as growth, development, and stress responses. Mutations in genes that are known to regulate tRNA modifications lead to a wide array of phenotypes and diseases including numerous cognitive and neurodevelopmental disorders, highlighting the critical role of tRNA modification in human disease. One such gene, THUMPD1, is involved in regulating tRNA N4-acetylcytidine modification (ac4C), and recently was proposed as a candidate gene for autosomal-recessive intellectual disability. Here, we present 13 individuals from 8 families who harbor rare loss-of-function variants in THUMPD1. Common phenotypic findings included global developmental delay, speech delay, moderate to severe intellectual deficiency, behavioral abnormalities such as angry outbursts, facial dysmorphism, and ophthalmological abnormalities. We demonstrate that the bi-allelic variants identified cause loss of function of THUMPD1 and that this defect results in a loss of ac4C modification in small RNAs, and of individually purified tRNA-Ser-CGA. We further corroborate this effect by showing a loss of tRNA acetylation in two CRISPR-Cas9-generated THUMPD1 KO cell lines. In addition, we also show the resultant amino acid substitution that occurs in a missense THUMPD1 allele identified in an individual with compound heterozygous variants results in a marked decrease in THUMPD1 stability and RNA-binding capacity. Taken together, these results suggest that the lack of tRNA acetylation due to THUMPD1 loss of function results in a syndromic form of intellectual disability associated with developmental delay, behavioral abnormalities, hearing loss, and facial dysmorphism.

Eukaryotic tRNA sequences present conserved and amino acid-specific structural signatures

Metazoan organisms have many tRNA genes responsible for decoding amino acids. The set of all tRNA genes can be grouped in sets of common amino acids and isoacceptor tRNAs that are aminoacylated by corresponding aminoacyl-tRNA synthetases. Analysis of tRNA alignments shows that, despite the high number of tRNA genes, specific tRNA sequence motifs are highly conserved across multicellular eukaryotes. The conservation often extends throughout the isoacceptors and isodecoders with, in some cases, two sets of conserved isodecoders. This study is focused on non-Watson–Crick base pairs in the helical stems, especially GoU pairs. Each of the four helical stems may contain one or more conserved GoU pairs. Some are amino acid specific and could represent identity elements for the cognate aminoacyl tRNA synthetases. Other GoU pairs are found in more than a single amino acid and could be critical for native folding of the tRNAs. Interestingly, some GoU pairs are anticodon-specific, and others are found in phylogenetically-specific clades. Although the distribution of conservation likely reflects a balance between accommodating isotype-specific functions as well as those shared by all tRNAs essential for ribosomal translation, such conservations may indicate the existence of specialized tRNAs for specific translation targets, cellular conditions, or alternative functions.

Methyltransferase METTL8 is required for 3-methylcytosine modification in human mitochondrial tRNAs

A subset of eukaryotic tRNAs is methylated in the anticodon loop, forming 3-methylcytosine (m³C) modifications. In mammals, the number of tRNAs containing m³C modifications has been expanded to include mitochondrial (mt) tRNA-Ser-UGA and mt-tRNA-Thr-UGU. However, whereas the enzymes catalyzing m³C formation in nuclear-encoded tRNAs have been identified, the proteins responsible for m³C modification in mt-tRNAs are unknown. Here, we show that m³C formation in human mt-tRNAs is dependent upon the Methyltransferase-Like 8 (METTL8) enzyme. We find that METTL8 is a mitochondria-associated protein that interacts with mitochondrial seryl-tRNA synthetase, as well as with mt-tRNAs containing m³C. We demonstrate that human cells deficient in METTL8 exhibit loss of m³C modification in mt-tRNAs, but not nuclear-encoded tRNAs. Consistent with the mitochondrial import of METTL8, the formation of m³C in METTL8-deficient cells could be rescued by re-expression of wildtype METTL8, but not by a METTL8 variant lacking the N-terminal mitochondrial localization signal. Notably, we found METTL8-deficiency in human cells causes alterations in the native migration pattern of mt-tRNA-Ser-UGA, suggesting a role for m³C in tRNA folding. Altogether, these findings demonstrate that METTL8 is required for m³C formation in mitochondrial tRNAs and uncover a potential function for m³C modification in mitochondrial tRNA structure.

Different modes of spacer acquisition by the Staphylococcus epidermidis type III-A CRISPR-Cas system

CRISPR-Cas systems provide prokaryotic organisms with an adaptive defense mechanism that acquires immunological memories of infections. This is accomplished by integration of short fragments from the genome of invaders such as phages and plasmids, called ‘spacers’, into the CRISPR locus of the host. Depending on their genetic composition, CRISPR-Cas systems can be classified into six types, I-VI, however spacer acquisition has been extensively studied only in type I and II systems. Here, we used an inducible spacer acquisition assay to study this process in the type III-A CRISPR-Cas system of Staphylococcus epidermidis, in the absence of phage selection. Similarly to type I and II spacer acquisition, this type III system uses Cas1 and Cas2 to preferentially integrate spacers from the chromosomal terminus and free dsDNA ends produced after DNA breaks, in a manner that is enhanced by the AddAB DNA repair complex. Surprisingly, a different mode of spacer acquisition from rRNA and tRNA loci, which spans only the transcribed sequences of these genes and is not enhanced by AddAB, was also detected. Therefore, our findings reveal both common mechanistic principles that may be conserved in all CRISPR-Cas systems, as well as unique and intriguing features of type III spacer acquisition.