MODOMICS is a database of RNA modifications that provides comprehensive information concerning the chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. In the current database version, we included the following new features and data: extended mass spectrometry and liquid chromatography data for modified nucleosides; links between human tRNA sequences and MINTbase - a framework for the interactive exploration of mitochondrial and nuclear tRNA fragments; new, machine-friendly system of unified abbreviations for modified nucleoside names; sets of modified tRNA sequences for two bacterial species, updated collection of mammalian tRNA modifications, 19 newly identified modified ribonucleosides and 66 functionally characterized proteins involved in RNA modification. Data from MODOMICS have been linked to the RNAcentral database of RNA sequences. MODOMICS is available at http://modomics.genesilico.pl.

The analytical identification of positional isomers (e.g., 3-, N⁴-, 5-methylcytidine) within the > 160 different post-transcriptional modifications found in RNA can be challenging. Conventional liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) approaches rely on chromatographic separation for accurate identification because the collision-induced dissociation (CID) mass spectra of these isomers nearly exclusively yield identical nucleobase ions (BH₂⁺) from the same molecular ion (MH⁺). Here, we have explored higher-energy collisional dissociation (HCD) as an alternative fragmentation technique to generate more informative product ions that can be used to differentiate positional isomers. LC-MS/MS of modified nucleosides characterized using HCD led to the creation of structure- and HCD energy-specific fragmentation patterns that generated unique fingerprints, which can be used to identify individual positional isomers even when they cannot be separated chromatographically. While particularly useful for identifying positional isomers, the fingerprinting capabilities enabled by HCD also offer the potential to generate HPLC-independent spectral libraries for the rapid analysis of modified ribonucleosides.

Extensive knowledge of both the nature and position of tRNA modifications in all cellular tRNAs has been limited to two bacteria, Escherichia coli and Mycoplasma capricolum. Bacillus subtilis sp subtilis strain 168 is the model Gram-positive bacteria and the list of the genes involved in tRNA modifications in this organism is far from complete. Mass spectrometry analysis of bulk tRNA extracted from B. subtilis, combined with next generation sequencing technologies and comparative genomic analyses, led to the identification of 41 tRNA modification genes with associated confidence scores. Many differences were found in this model Gram-positive bacteria when compared to E. coli. In general, B. subtilis tRNAs are less modified than those in E. coli, even if some modifications, such as m¹A22 or ms²t⁶A, are only found in the model Gram-positive bacteria. Many examples of non-orthologous displacements and of variations in the most complex pathways are described. Paralog issues make uncertain direct annotation transfer from E. coli to B. subtilis based on homology only without further experimental validation. This difficulty was shown with the identification of the B. subtilis enzyme that introduces ψ at positions 31/32 of the tRNAs. This work presents the most up to date list of tRNA modification genes in B. subtilis, identifies the gaps in knowledge, and lays the foundation for further work to decipher the physiological role of tRNA modifications in this important model organism and other bacteria.

The reduction of epoxyqueuosine (oQ) is the last step in the synthesis of the tRNA modification queuosine (Q). While the epoxyqueuosine reductase (EC 1.17.99.6) enzymatic activity was first described 30 years ago, the encoding gene queG was only identified in Escherichia coli in 2011. Interestingly, queG is absent from a large number of sequenced genomes that harbor Q synthesis or salvage genes, suggesting the existence of an alternative epoxyqueuosine reductase in these organisms. By analyzing phylogenetic distributions, physical gene clustering, and fusions, members of the Domain of Unknown Function 208 (DUF208) family were predicted to encode for an alternative epoxyqueuosine reductase. This prediction was validated with genetic methods. The Q modification is present in Lactobacillus salivarius, an organism missing queG but harboring the duf208 gene. Acinetobacter baylyi ADP1 is one of the few organisms that harbor both QueG and DUF208, and deletion of both corresponding genes was required to observe the absence of Q and the accumulation of oQ in tRNA. Finally, the conversion oQ to Q was restored in an E. coli queG mutant by complementation with plasmids harboring duf208 genes from different bacteria. Members of the DUF208 family are not homologous to QueG enzymes, and thus, duf208 is a non-orthologous replacement of queG. We propose to name DUF208 encoding genes as queH. While QueH contains conserved cysteines that could be involved in the coordination of a Fe/S center in a similar fashion to what has been identified in QueG, no cobalamin was identified associated with recombinant QueH protein.