Evaluation of commercial RNA extraction kits for long-read metatranscriptomics in soil.

Metatranscriptomic analysis of the soil microbiome has the potential to reveal molecular mechanisms that drive soil processes regulated by the microbial community. Therefore, RNA samples must be of sufficient yield and quality to robustly quantify differential gene expression. While short-read sequencing technology is often favoured for metatranscriptomics, long-read sequencing has the potential to provide several benefits over short-read technologies. The ability to resolve complete transcripts on a portable sequencing platform for a relatively low capital expenditure makes Oxford Nanopore Technology an attractive prospect for addressing many of the challenges of soil metatranscriptomics. To fully enable long-read metatranscriptomic analysis of the functional molecular pathways expressed in these diverse habitats, RNA purification methods from soil must be optimised for long-read sequencing. Here we compare RNA samples purified using five commercially available extraction kits designed for use with soil. We found that the Qiagen RNeasy PowerSoil Total RNA Kit performed the best across RNA yield, quality and purity and was robust across different soil types. We found that sufficient sequencing depth can be achieved to characterise the active community for total RNA samples using Oxford Nanopore Technology, and discuss its current limitations for differential gene expression analysis in soil studies.

RluA is the major mRNA pseudouridine synthase in Escherichia coli.

Pseudouridine (Ψ) is an ubiquitous RNA modification, present in the tRNAs and rRNAs of species across all domains of life. Conserved pseudouridine synthases modify the mRNAs of diverse eukaryotes, but the modification has yet to be identified in bacterial mRNAs. Here, we report the discovery of pseudouridines in mRNA from E. coli. By testing the mRNA modification capacity of all 11 known pseudouridine synthases, we identify RluA as the predominant mRNA-modifying enzyme. RluA, a known tRNA and 23S rRNA pseudouridine synthase, modifies at least 31 of the 44 high-confidence sites we identified in E. coli mRNAs. Using RNA structure probing data to inform secondary structures, we show that the target sites of RluA occur in a common sequence and structural motif comprised of a ΨURAA sequence located in the loop of a short hairpin. This recognition element is shared with previously identified target sites of RluA in tRNAs and rRNA. Overall, our work identifies pseudouridine in key mRNAs and suggests the capacity of Ψ to regulate the transcripts that contain it.

Mechanistic analysis of Riboswitch Ligand interactions provides insights into pharmacological control over gene expression.

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Exploring the transcription start sites and other genomic features facilitates the accurate identification and annotation of small RNAs across multiple stress conditions in Mycobacterium tuberculosis.

Mycobacterium tuberculosis (MTB) is a pathogen that is known for its ability to persist in harsh environments and cause chronic infections. Understanding the regulatory networks of MTB is crucial for developing effective treatments. Small regulatory RNAs (sRNAs) play important roles in gene expression regulation in all kingdoms of life, and their classification based solely on genomic location can be imprecise due to the computational-based prediction of protein-coding genes in bacteria, which often neglects segments of mRNA such as 5'UTRs, 3'UTRs, and intercistronic regions of operons. To address this issue, our study simultaneously discovered genomic features such as TSSs, UTRs, and operons together with sRNAs in the M. tuberculosis H37Rv strain (ATCC 27294) across multiple stress conditions. Our analysis identified 1,376 sRNA candidates and 8,173 TSSs in MTB, providing valuable insights into its complex regulatory landscape. TSS mapping enabled us to classify these sRNAs into more specific categories, including promoter-associated sRNAs, 5'UTR-derived sRNAs, 3'UTR-derived sRNAs, true intergenic sRNAs, and antisense sRNAs. Three of these sRNA candidates were experimentally validated using 3'-RACE-PCR: predictedRNA_0240, predictedRNA_0325, and predictedRNA_0578. Future characterization and validation are necessary to fully elucidate the functions and roles of these sRNAs in MTB. Our study is the first to simultaneously unravel TSSs and sRNAs in MTB and demonstrate that the identification of other genomic features, such as TSSs, UTRs, and operons, allows for more accurate and specific classification of sRNAs.

A 3'UTR-derived small RNA represses pneumolysin synthesis and facilitates pneumococcal brain invasion.

Pneumolysin (Ply) of Streptococcus pneumoniae (pneumococcus) at relatively high and low levels facilitates pneumococcal invasion into the lung and brain, respectively; however, the regulatory mechanisms of Ply expression are poorly understood. Here, we find that a small RNA plyT, processed from the 3'UTR of the ply operon, is expressed higher in anaerobically- than in statically-cultured pneumococcus D39. Using bioinformatic, biochemical and genetic approaches, we reveal that PlyT inhibits Ply synthesis and hemolytic activities by pairing with an RBS-embedded intergenic region of the ply operon. The RNA-binding protein SPD_1558 facilitates the pairing. Importantly, PlyT inhibition of Ply synthesis is stronger in anaerobic culture and leads to lower Ply abundance. Deletion of plyT decreases the number of pneumococci in the infected mouse brain and reduces the virulence, demonstrating that PlyT-regulated lower Ply in oxygen-void microenvironments, such as the blood, is important for pneumococcus to cross the blood-brain barrier and invade the brain. PlyT-mediated repression of Ply synthesis at anoxic niches is also verified in pneumococcal serotype 4 and 14 strains; moreover, the ply operon with a 3'UTR-embedded plyT, and the pairing sequences of IGR and plyT are highly conserved among pneumococcal strains, implying PlyT-regulated Ply synthesis might be widely employed by pneumococcus.

Antiterminator LoaP loads onto RNA to chase a runaway RNA polymerase.

In this issue of Structure, Elghondakly et al.1 present the crystal structure of Thermoanaerobacter pseudethanolicus antiterminator LoaP, a member of a ubiquitous family of NusG transcription factors, bound to its target, a dfn RNA hairpin. LoaP uses RNA as a recognition determinant, which is unique among NusG paralogs and makes unusual contacts in the major groove of the RNA.

Rli51 Attenuates Transcription of the Listeria Pathogenicity Island 1 Gene mpl and Functions as a Trans-Acting sRNA in Intracellular Bacteria.

Listeria pathogenicity island 1 (LIPI-1) is a genetic region containing a cluster of genes essential for virulence of the bacterial pathogen Listeria monocytogenes. Main virulence factors in LIPI-1 include long 5' untranslated regions (5'UTRs), among which is Rli51, a small RNA (sRNA) in the 5'UTR of the Zn-metalloprotease-coding mpl. So far, Rli51 function and molecular mechanisms have remained obscure. Here, we show that Rli51 exhibits a dual mechanism of regulation, functioning as a cis- and as a trans-acting sRNA. Under nutrient-rich conditions, rli51-mpl transcription is prematurely terminated, releasing a short 121-nucleotide-long sRNA. Rli51 is predicted to function as a transcription attenuator that can fold into either a terminator or a thermodynamically more stable antiterminator. We show that the sRNA Rli21/RliI binds to a single-stranded RNA loop in Rli51, which is essential to mediate premature transcription termination, suggesting that sRNA binding could stabilize the terminator fold. During intracellular infection, rli51 transcription is increased, which generates a higher abundance of the short Rli51 sRNA and allows for transcriptional read-through into mpl. Comparative intracellular bacterial transcriptomics in rli51-null mutants and the wild-type reference strain EGD-e suggests that Rli51 upregulates iron-scavenging proteins and downregulates virulence factors from LIPI-1. MS2 affinity purification confirmed that Rli51 binds transcripts of the heme-binding protein Lmo2186 and Lmo0937 in vivo. These results prove that Rli51 functions as a trans-acting sRNA in intracellular bacteria. Our research shows a growth condition-dependent mechanism of regulation for Rli51, preventing unintended mpl transcription in extracellular bacteria and regulating genes important for virulence in intracellular bacteria.

Hfq-binding small RNA PqsS regulates Pseudomonas aeruginosa pqs quorum sensing system and virulence.

Pseudomonas aeruginosa is a widespread nosocomial pathogen with a significant to cause both severe planktonic acute and biofilm-related chronic infections. Small RNAs (sRNAs) are noncoding regulatory molecules that are stabilized by the RNA chaperone Hfq to trigger various virulence-related signaling pathways. Here, we identified an Hfq-binding sRNA in P. aeruginosa PAO1, PqsS, which promotes bacterial pathogenicity and pseudomonas quinolone signal quorum sensing (pqs QS) system. Specifically, PqsS enhanced acute bacterial infections by inducing host cell death and promoting rhamnolipid-regulated swarming motility. Meanwhile, PqsS reduced chronic infection traits including biofilm formation and antibiotic resistance. Moreover, PqsS repressed pqsL transcript, increasing PQS levels for pqs QS. A PQS-rich environment promoted PqsS expression, thus forming a positive feedback loop. Furthermore, we demonstrated that the PqsS interacts and destabilizes the pqsL mRNA by recruiting RNase E to drive degradation. These findings provide insights for future research on P. aeruginosa pathogenesis and targeted treatment.

Genome-wide sRNA and mRNA transcriptomic profiling insights into carbapenem-resistant Acinetobacter baumannii.

Introduction: Acinetobacter baumannii (AB) is rising as a human pathogen of critical priority worldwide as it is the leading cause of opportunistic infections in healthcare settings and carbapenem-resistant AB is listed as a "super bacterium" or "priority pathogen for drug resistance" by the World Health Organization. Methods: Clinical isolates of A. baumannii were collected and tested for antimicrobial susceptibility. Among them, carbapenem-resistant and carbapenem-sensitive A. baumannii were subjected to prokaryotic transcriptome sequencing. The change of sRNA and mRNA expression was analyzed by bioinformatics and validated by quantitative reverse transcription-PCR. Results: A total of 687 clinical isolates were collected, of which 336 strains of A. baumannii were resistant to carbapenem. Five hundred and six differentially expressed genes and nineteen differentially expressed sRNA candidates were discovered through transcriptomic profile analysis between carbapenem-resistant isolates and carbapenem-sensitive isolates. Possible binding sites were predicted through software for sRNA21 and adeK, sRNA27 and pgaC, sRNA29 and adeB, sRNA36 and katG, indicating a possible targeting relationship. A negative correlation was shown between sRNA21 and adeK (r = -0.581, P = 0.007), sRNA27 and pgaC (r = -0.612, P = 0.004), sRNA29 and adeB (r = -0.516, P = 0.020). Discussion: This study preliminarily screened differentially expressed mRNA and sRNA in carbapenem-resistant A. baumannii, and explored possible targeting relationships, which will help further reveal the resistance mechanism and provide a theoretical basis for the development of drugs targeting sRNA for the prevention and treatment of carbapenem-resistant A. baumannii infection.

Physiological significance of the two isoforms of initiator tRNAs in Escherichia coli.

Escherichia coli possesses four initiator tRNA (i-tRNA) genes, three of which are present together as metZWV and the fourth one as metY. In E. coli B, all four genes (metZWV and metY) encode i-tRNAfMet1, in which the G at position 46 is modified to m7G46 by TrmB (m7G methyltransferase). However, in E. coli K, because of a single-nucleotide polymorphism, metY encodes a variant, i-tRNAfMet2, having an A in place of m7G46. We generated E. coli strains to explore the importance of this polymorphism in i-tRNAs. The strains were sustained either on metYA46 (metY of E. coli K origin encoding i-tRNAfMet2) or its derivative metYG46 (encoding i-tRNAfMet1) in single (chromosomal) or plasmid-borne copies. We show that the strains sustained on i-tRNAfMet1 have a growth fitness advantage over those sustained on i-tRNAfMet2. The growth fitness advantages are more pronounced for the strains sustained on i-tRNAfMet1 in nutrient-rich media than in nutrient-poor media. The growth fitness of the strains correlates well with the relative stabilities of the i-tRNAs in vivo. Furthermore, the atomistic molecular dynamics simulations support the higher stability of i-tRNAfMet1 than that of i-tRNAfMet2. The stability of i-tRNAfMet1 remains unaffected upon the deletion of TrmB. These studies highlight how metYG46 and metYA46 alleles might influence the growth fitness of E. coli under certain nutrient-limiting conditions. IMPORTANCE: Escherichia coli harbors four initiator tRNA (i-tRNA) genes: three of these at metZWV and the fourth one at metY loci. In E. coli B, all four genes encode i-tRNAfMet1. In E. coli K, because of a single-nucleotide polymorphism, metY encodes a variant, i-tRNAfMet2, having an A in place of G at position 46 of i-tRNA sequence in metY. We show that G46 confers stability to i-tRNAfMet1. The strains sustained on i-tRNAfMet1 have a growth fitness advantage over those sustained on i-tRNAfMet2. Strains harboring metYG46 (B mimic) or metYA46 (K mimic) show that while in the nutrient-rich media, the K mimic is outcompeted rapidly; in the nutrient-poor medium, the K mimic is outcompeted less rapidly.

Signal recognition particle RNA is critical for genetic competence and virulence of Streptococcus pneumoniae.

Streptococcus pneumoniae (pneumococcus) causes a wide range of important human infectious diseases, including pneumonia, pneumonia-derived sepsis, otitis media, and meningitis. Pneumococcus produces numerous secreted proteins that are critical for normal physiology and pathogenesis. The membrane targeting and translocation of these secreted proteins are partly mediated by the signal recognition particle (SRP) complex, which consists of 4.5S small cytoplasmic RNA (ScRNA), and the Ffh, and FtsY proteins. Here, we report that pneumococcal ∆scRNA, ∆ffh, and ∆ftsY mutants were significantly impaired in competence induction, competence pili production, exogenous DNA uptake, and genetic transformation. Also, the ∆scRNA mutant was significantly attenuated in the mouse models of bacteremia and pneumonia. Interestingly, unlike the ∆scRNA, both ∆ffh and ∆ftsY mutants had growth defects on Todd-Hewitt Agar, which were alleviated by the provision of free amino acids or serum. Differences in nutritional requirements between ∆ffh and ∆ftsY vs ∆scRNA suggest that Ffh and FtsY may be partially functional in the absence of ScRNA. Finally, the insertase YidC2, which could functionally rescue some SRP mutations in other streptococcal species, was not essential for pneumococcal genetic transformation. Collectively, these results indicate that ScRNA is crucial for the successful development of genetic competence and virulence in pneumococcus. IMPORTANCE: Streptococcus pneumoniae (pneumococcus) causes multiple important infectious diseases in humans. The signal recognition particle (SRP) complex, which comprised 4.5S small cytoplasmic RNA (ScRNA), and the Ffh and FtsY proteins, mediates membrane targeting and translocation of secreted proteins in all organisms. However, the role of SRP and ScRNA has not been characterized during the induction of the competence system for genetic transformation and virulence in pneumococcus. By using a combination of genetic, biochemical, proteomic, and imaging approaches, we demonstrated that the SRP complex plays a significant role in membrane targeting of competence system-regulated effectors important for genetic transformation, virulence during bacteremia and pneumonia infections, and nutritional acquisition.

Rho-dependent transcriptional switches regulate the bacterial response to cold shock.

Binding of the bacterial Rho helicase to nascent transcripts triggers Rho-dependent transcription termination (RDTT) in response to cellular signals that modulate mRNA structure and accessibility of Rho utilization (Rut) sites. Despite the impact of temperature on RNA structure, RDTT was never linked to the bacterial response to temperature shifts. We show that Rho is a central player in the cold-shock response (CSR), challenging the current view that CSR is primarily a posttranscriptional program. We identify Rut sites in 5'-untranslated regions of key CSR genes/operons (cspA, cspB, cspG, and nsrR-rnr-yjfHI) that trigger premature RDTT at 37°C but not at 15°C. High concentrations of RNA chaperone CspA or nucleotide changes in the cspA mRNA leader reduce RDTT efficiency, revealing how RNA restructuring directs Rho to activate CSR genes during the cold shock and to silence them during cold acclimation. These findings establish a paradigm for how RNA thermosensors can modulate gene expression.

Translational T-box riboswitches bind tRNA by modulating conformational flexibility.

T-box riboswitches are noncoding RNA elements involved in genetic regulation of most Gram-positive bacteria. They regulate amino acid metabolism by assessing the aminoacylation status of tRNA, subsequently affecting the transcription or translation of downstream amino acid metabolism-related genes. Here we present single-molecule FRET studies of the Mycobacterium tuberculosis IleS T-box riboswitch, a paradigmatic translational T-box. Results support a two-step binding model, where the tRNA anticodon is recognized first, followed by interactions with the NCCA sequence. Furthermore, after anticodon recognition, tRNA can transiently dock into the discriminator domain even in the absence of the tRNA NCCA-discriminator interactions. Establishment of the NCCA-discriminator interactions significantly stabilizes the fully bound state. Collectively, the data suggest high conformational flexibility in translational T-box riboswitches; and supports a conformational selection model for NCCA recognition. These findings provide a kinetic framework to understand how specific RNA elements underpin the binding affinity and specificity required for gene regulation.

Critical factors for precise and efficient RNA cleavage by RNase Y in Staphylococcus aureus.

Cellular processes require precise and specific gene regulation, in which continuous mRNA degradation is a major element. The mRNA degradation mechanisms should be able to degrade a wide range of different RNA substrates with high efficiency, but should at the same time be limited, to avoid killing the cell by elimination of all cellular RNA. RNase Y is a major endoribonuclease found in most Firmicutes, including Bacillus subtilis and Staphylococcus aureus. However, the molecular interactions that direct RNase Y to cleave the correct RNA molecules at the correct position remain unknown. In this work we have identified transcripts that are homologs in S. aureus and B. subtilis, and are RNase Y targets in both bacteria. Two such transcript pairs were used as models to show a functional overlap between the S. aureus and the B. subtilis RNase Y, which highlighted the importance of the nucleotide sequence of the RNA molecule itself in the RNase Y targeting process. Cleavage efficiency is driven by the primary nucleotide sequence immediately downstream of the cleavage site and base-pairing in a secondary structure a few nucleotides downstream. Cleavage positioning is roughly localised by the downstream secondary structure and fine-tuned by the nucleotide immediately upstream of the cleavage. The identified elements were sufficient for RNase Y-dependent cleavage, since the sequence elements from one of the model transcripts were able to convert an exogenous non-target transcript into a target for RNase Y.

Synonymous codon substitutions modulate transcription and translation of a divergent upstream gene by modulating antisense RNA production.

Synonymous codons were originally viewed as interchangeable, with no phenotypic consequences. However, substantial evidence has now demonstrated that synonymous substitutions can perturb a variety of gene expression and protein homeostasis mechanisms, including translational efficiency, translational fidelity, and cotranslational folding of the encoded protein. To date, most studies of synonymous codon-derived perturbations have focused on effects within a single gene. Here, we show that synonymous codon substitutions made far within the coding sequence of Escherichia coli plasmid-encoded chloramphenicol acetyltransferase (cat) can significantly increase expression of the divergent upstream tetracycline resistance gene, tetR. In four out of nine synonymously recoded cat sequences tested, expression of the upstream tetR gene was significantly elevated due to transcription of a long antisense RNA (asRNA) originating from a transcription start site within cat. Surprisingly, transcription of this asRNA readily bypassed the native tet transcriptional repression mechanism. Even more surprisingly, accumulation of the TetR protein correlated with the level of asRNA, rather than total tetR RNA. These effects of synonymous codon substitutions on transcription and translation of a neighboring gene suggest that synonymous codon usage in bacteria may be under selection to both preserve the amino acid sequence of the encoded gene and avoid DNA sequence elements that can significantly perturb expression of neighboring genes. Avoiding such sequences may be especially important in plasmids and prokaryotic genomes, where genes and regulatory elements are often densely packed. Similar considerations may apply to the design of genetic circuits for synthetic biology applications.

Effective RNA isolation method for gram-positive and acid-fast bacteria: metamorphosed from conventional RNA isolation techniques.

The RNA-based study provides an excellent indication of an organism's gene expression profile. Obtaining high-yield and high-purity RNA from Gram-positive and acid-fast bacteria is difficult without high-end kits and facilities. We optimised effective and simple protocol for RNA isolation that is a combination of enzymatic, physical and chemical treatment to disrupt cells. We successfully isolated high quality intact total RNA with yields ranging from 23.13 ± 0.40 to 61.51 ± 0.27 µg and the 260/280 purity ratio of 1.95 ± 0.01 to 2.05 ± 0.01 from Staphylococcus aureus, Staphylococcus epidermidis, Enterococcus faecalis, and Mycobacterium smegmatis. These results represents a significantly enhanced yield and purity compared to other combination of techniques which we performed. Compared to previous studies the yield obtained by this method is high for the studied organisms. Furthermore the yielded RNA was successfully used for downstream applications such as quantitative real time PCR. The described method can be easily optimised and used for various bacteria.

Unraveling RNA contribution to the molecular origins of bacterial surface-enhanced Raman spectroscopy (SERS) signals.

Surface-enhanced Raman spectroscopy (SERS) is widely utilized in bacterial analyses, with the dominant SERS peaks attributed to purine metabolites released during sample preparation. Although adenosine triphosphate (ATP) and nucleic acids are potential molecular origins of these metabolites, research on their exact contributions remains limited. This study explored purine metabolite release from E. coli and RNA integrity following various sample preparation methods. Standard water washing generated dominant SERS signals within 10 s, a duration shorter than the anticipated RNA half-lives under starvation. Evaluating RNA integrity indicated that the most abundant ribosomal RNA species remained intact for hours post-washing, whereas messenger RNA and transfer RNA species degraded gradually. This suggests that bacterial SERS signatures observed after the typical washing step could originate from only a small fraction of endogenous purine-containing molecules. In contrast, acid depurination led to degradation of most RNA species, releasing about 40 times more purine derivatives than water washing. Mild heating also instigated the RNA degradation and released more purine derivatives than water washing. Notably, differences were also evident in the dominant SERS signals following these treatments. This work provides insights into SERS-based studies of purine metabolites released by bacteria and future development of methodologies.

What goes up must come down: off switches for regulatory RNAs.

Small RNAs base pair with and regulate mRNA translation and stability. For both bacterial small regulatory RNAs and eukaryotic microRNAs, association with partner proteins is critical for the stability and function of the regulatory RNAs. We review the mechanisms for degradation of these RNAs: displacement of the regulatory RNA from its protein partner (in bacteria) or destruction of the protein and its associated microRNAs (in eukaryotes). These mechanisms can allow specific destruction of a regulatory RNA via pairing with a decay trigger RNA or function as global off switches by disrupting the stability or function of the protein partner.

Induction of bacterial expression at the mRNA level by light.

Vital organismal processes, including development, differentiation and adaptation, involve altered gene expression. Although expression is frequently controlled at the transcriptional stage, various regulation mechanisms operate at downstream levels. Here, we leverage the photoreceptor NmPAL to optogenetically induce RNA refolding and the translation of bacterial mRNAs. Blue-light-triggered NmPAL binding disrupts a cis-repressed mRNA state, thereby relieves obstruction of translation initiation, and upregulates gene expression. Iterative probing and optimization of the circuit, dubbed riboptoregulator, enhanced induction to 30-fold. Given action at the mRNA level, the riboptoregulator can differentially regulate individual structural genes within polycistronic operons. Moreover, it is orthogonal to and can be wed with other gene-regulatory circuits for nuanced and more stringent gene-expression control. We thus advance the pAurora2 circuit that combines transcriptional and translational mechanisms to optogenetically increase bacterial gene expression by >1000-fold. The riboptoregulator strategy stands to upgrade numerous regulatory circuits and widely applies to expression control in microbial biotechnology, synthetic biology and materials science.

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.