Mapping and Quantification of tRNA 2'-O-Methylation by RiboMethSeq

A snoRNA-tRNA modification network governs codon-biased cellular states

The dynamic balance between tRNA supply and codon usage demand is a fundamental principle in the cellular translation economy. However, the regulation and functional consequences of this balance remain unclear. Here, we use PARIS2 interactome capture, structure modeling, conservation analysis, RNA-protein interaction analysis, and modification mapping to reveal the targets of hundreds of snoRNAs, many of which were previously considered orphans. We identify a snoRNA-tRNA interaction network that is required for global tRNA modifications, including 2'-O-methylation and others. Loss of Fibrillarin, the snoRNA-guided 2'-O-methyltransferase, induces global upregulation of tRNA fragments, a large group of regulatory RNAs. In particular, the snoRNAs D97/D133 guide the 2'-O-methylation of multiple tRNAs, especially for the amino acid methionine (Met), a protein-intrinsic antioxidant. Loss of D97/D133 snoRNAs in human HEK293 cells reduced target tRNA levels and induced codon adaptation of the transcriptome and translatome. Both single and double knockouts of D97 and D133 in HEK293 cells suppress Met-enriched proliferation-related gene expression programs, including, translation, splicing, and mitochondrial energy metabolism, and promote Met-depleted programs related to development, differentiation, and morphogenesis. In a mouse embryonic stem cell model of development, knockdown and knockout of D97/D133 promote differentiation to mesoderm and endoderm fates, such as cardiomyocytes, without compromising pluripotency, consistent with the enhanced development-related gene expression programs in human cells. This work solves a decades-old mystery about orphan snoRNAs and reveals a function of snoRNAs in controlling the codon-biased dichotomous cellular states of proliferation and development.

Quantification of substoichiometric modification reveals global tsRNA hypomodification, preferences for angiogenin-mediated tRNA cleavage, and idiosyncratic epitranscriptomes of human neuronal cell-lines

Modification of tRNA is an integral part of the epitranscriptome with a particularly pronounced potential to generate diversity in RNA expression. Eukaryotic tRNA contains modifications in up to 20% of their nucleotides, but not all sites are always fully modified. Combinations and permutations of partially modified sites in tRNAs can generate a plethora of tRNA isoforms, termed modivariants. Here, we investigate the stoichiometry of incompletely modified sites in tRNAs from human cell lines for their information content. Using a panel of RNA modification mapping methods, we assess the stoichiometry of sites that contain the modifications 5-methylcytidine (m5C), 2'-O-ribose methylation (Nm), 3-methylcytidine (m3C), 7-methylguanosine (m7G), and Dihydrouridine (D). We discovered that up to 75% of sites can be incompletely modified and that the differential modification status of a cellular tRNA population holds information that allows to discriminate e.g. different cell lines. As a further aspect, we investigated potential causal connectivity between tRNA modification and its processing into tRNA fragments (tiRNAs and tRFs). Upon exposure of cultured living cells to cell-penetrating angiogenin, the modification patterns of the corresponding RNA populations was changed. Importantly, we also found that tsRNAs were significantly less modified than their parent tRNAs at numerous sites, suggesting that tsRNAs might derive chiefly from hypomodified tRNAs.

Epitranscriptome: Review of Top 25 Most-Studied RNA Modifications

The alphabet of building blocks for RNA molecules is much larger than the standard four nucleotides. The diversity is achieved by the post-transcriptional biochemical modification of these nucleotides into distinct chemical entities that are structurally and functionally different from their unmodified counterparts. Some of these modifications are constituent and critical for RNA functions, while others serve as dynamic markings to regulate the fate of specific RNA molecules. Together, these modifications form the epitranscriptome, an essential layer of cellular biochemistry. As of the time of writing this review, more than 300 distinct RNA modifications from all three life domains have been identified. However, only a few of the most well-established modifications are included in most reviews on this topic. To provide a complete overview of the current state of research on the epitranscriptome, we analyzed the extent of the available information for all known RNA modifications. We selected 25 modifications to describe in detail. Summarizing our findings, we describe the current status of research on most RNA modifications and identify further developments in this field.

Computational methods for RNA modification detection from nanopore direct RNA sequencing data

The covalent modification of RNA molecules is a pervasive feature of all classes of RNAs and has fundamental roles in the regulation of several cellular processes. Mapping the location of RNA modifications transcriptome-wide is key to unveiling their role and dynamic behaviour, but technical limitations have often hampered these efforts. Nanopore direct RNA sequencing is a third-generation sequencing technology that allows the sequencing of native RNA molecules, thus providing a direct way to detect modifications at single-molecule resolution. Despite recent advances, the analysis of nanopore sequencing data for RNA modification detection is still a complex task that presents many challenges. Many works have addressed this task using different approaches, resulting in a large number of tools with different features and performances. Here we review the diverse approaches proposed so far and outline the principles underlying currently available algorithms.

The occurrence order and cross-talk of different tRNA modifications

Chemical modifications expand the composition of RNA molecules from four standard nucleosides to over 160 modified nucleosides, which greatly increase the complexity and utility of RNAs. Transfer RNAs (tRNAs) are the most heavily modified cellular RNA molecules and contain the largest variety of modifications. Modification of tRNAs is pivotal for protein synthesis and also precisely regulates the noncanonical functions of tRNAs. Defects in tRNA modifications lead to numerous human diseases. Up to now, more than 100 types of modifications have been found in tRNAs. Intriguingly, some modifications occur widely on all tRNAs, while others only occur on a subgroup of tRNAs or even only a specific tRNA. The modification frequency of each tRNA is approximately 7% to 25%, with 5-20 modification sites present on each tRNA. The occurrence and modulation of tRNA modifications are specifically noticeable as plenty of interplays among different sites and modifications have been discovered. In particular, tRNA modifications are responsive to environmental changes, indicating their dynamic and highly organized nature. In this review, we summarized the known occurrence order, cross-talk, and cooperativity of tRNA modifications.

Mapping of 7-methylguanosine (m7G), 3-methylcytidine (m3C), dihydrouridine (D) and 5-hydroxycytidine (ho5C) RNA modifications by AlkAniline-Seq

Precise and reliable mapping of modified nucleotides in RNA is a challenging task in epitranscriptomics analysis. Only deep sequencing-based methods are able to provide both, a single-nucleotide resolution and sufficient selectivity and sensitivity. A number of protocols employing specific chemical reagents to distinguish modified RNA nucleotides from canonical parental residues have already proven their performance. We developed a deep-sequencing analytical pipeline for simultaneous detection of several modified nucleotides of different nature (methylation, hydroxylation, reduction) in RNA. The AlkAniline-Seq protocol uses intrinsic fragility of the N-glycosidic bond present in certain modified residues (7-methylguanosine (m⁷G), 3-methylcytidine (m³C), dihydrouridine (D) and 5-hydroxycytidine (ho⁵C)) to induce cleavage under heat combined with alkaline conditions. The resulting RNA abasic site is decomposed by aniline-driven β-elimination and creates a 5'-phosphate (5'-P) at the adjacent N+1 residue. This 5'-P is the crucial entry point for a highly selective ligation of sequencing adapters during the subsequent Illumina library preparation protocol. AlkAniline-Seq protocol has a very low background, and is both highly sensitive and specific. Applications of AlkAniline-Seq include mapping of m⁷G, m³C, D, and ho⁵C in variety of cellular RNAs, including in particular rRNAs and tRNAs.

Current Practice in Bicistronic IRES Reporter Use: A Systematic Review

Bicistronic reporter assays have been instrumental for transgene expression, understanding of internal ribosomal entry site (IRES) translation, and identification of novel cap-independent translational elements (CITE). We observed a large methodological variability in the use of bicistronic reporter assays and data presentation or normalization procedures. Therefore, we systematically searched the literature for bicistronic IRES reporter studies and analyzed methodological details, data visualization, and normalization procedures. Two hundred fifty-seven publications were identified using our search strategy (published 1994-2020). Experimental studies on eukaryotic adherent cell systems and the cell-free translation assay were included for further analysis. We evaluated the following methodological details for 176 full text articles: the bicistronic reporter design, the cell line or type, transfection methods, and time point of analyses post-transfection. For the cell-free translation assay, we focused on methods of in vitro transcription, type of translation lysate, and incubation times and assay temperature. Data can be presented in multiple ways: raw data from individual cistrons, a ratio of the two, or fold changes thereof. In addition, many different control experiments have been suggested when studying IRES-mediated translation. In addition, many different normalization and control experiments have been suggested when studying IRES-mediated translation. Therefore, we also categorized and summarized their use. Our unbiased analyses provide a representative overview of bicistronic IRES reporter use. We identified parameters that were reported inconsistently or incompletely, which could hamper data reproduction and interpretation. On the basis of our analyses, we encourage adhering to a number of practices that should improve transparency of bicistronic reporter data presentation and improve methodological descriptions to facilitate data replication.

Non-Redundant tRNA Reference Sequences for Deep Sequencing Analysis of tRNA Abundance and Epitranscriptomic RNA Modifications

Analysis of RNA by deep-sequencing approaches has found widespread application in modern biology. In addition to measurements of RNA abundance under various physiological conditions, such techniques are now widely used for mapping and quantification of RNA modifications. Transfer RNA (tRNA) molecules are among the frequent targets of such investigation, since they contain multiple modified residues. However, the major challenge in tRNA examination is related to a large number of duplicated and point-mutated genes encoding those RNA molecules. Moreover, the existence of multiple isoacceptors/isodecoders complicates both the analysis and read mapping. Existing databases for tRNA sequencing provide near exhaustive listings of tRNA genes, but the use of such highly redundant reference sequences in RNA-seq analyses leads to a large number of ambiguously mapped sequencing reads. Here we describe a relatively simple computational strategy for semi-automatic collapsing of highly redundant tRNA datasets into a non-redundant collection of reference tRNA sequences. The relevance of the approach was validated by analysis of experimentally obtained tRNA-sequencing datasets for different prokaryotic and eukaryotic model organisms. The data demonstrate that non-redundant tRNA reference sequences allow improving unambiguous mapping of deep sequencing data.

Bacterial tRNA 2'-O-methylation is dynamically regulated under stress conditions and modulates innate immune response

RNA modifications are a well-recognized way of gene expression regulation at the post-transcriptional level. Despite the importance of this level of regulation, current knowledge on modulation of tRNA modification status in response to stress conditions is far from being complete. While it is widely accepted that tRNA modifications are rather dynamic, such variations are mostly assessed in terms of total tRNA, with only a few instances where changes could be traced to single isoacceptor species. Using Escherichia coli as a model system, we explored stress-induced modulation of 2'-O-methylations in tRNAs by RiboMethSeq. This analysis and orthogonal analytical measurements by LC-MS show substantial, but not uniform, increase of the Gm18 level in selected tRNAs under mild bacteriostatic antibiotic stress, while other Nm modifications remain relatively constant. The absence of Gm18 modification in tRNAs leads to moderate alterations in E. coli mRNA transcriptome, but does not affect polysomal association of mRNAs. Interestingly, the subset of motility/chemiotaxis genes is significantly overexpressed in ΔTrmH mutant, this corroborates with increased swarming motility of the mutant strain. The stress-induced increase of tRNA Gm18 level, in turn, reduced immunostimulation properties of bacterial tRNAs, which is concordant with the previous observation that Gm18 is a suppressor of Toll-like receptor 7 (TLR7)-mediated interferon release. This documents an effect of stress induced modulation of tRNA modification that acts outside protein translation.

From canonical to modified nucleotides: balancing translation and metabolism

Every type of nucleic acid in cells may undergo some kind of post-replicative or post-transcriptional chemical modification. Recent evidence has highlighted their importance in biology and their chemical complexity. In the following pages, we will describe new discoveries of modifications, with a focus on tRNA and mRNA. We will highlight current challenges and advances in modification detection and we will discuss how changes in nucleotide post-transcriptional modifications may affect cell homeostasis leading to malfunction. Although, RNA modifications prevail in all forms of life, the present review will focus on eukaryotic systems, where the great degree of intracellular compartmentalization provides barriers and filters for the level at which a given RNA is modified and will of course affect its fate and function. Additionally, although we will mention rRNA modification and modifications of the mRNA 5'-CAP structure, this will only be discussed in passing, as many substantive reviews have been written on these subjects. Here we will not spend much time describing all the possible modifications that have been observed; truly a daunting task. For reference, Bujnicki and coworkers have created MODOMICS, a useful repository for all types of modifications and their associated enzymes. Instead we will discuss a few examples, which illustrate our arguments on the connection of modifications, metabolism and ultimately translation. The fact remains, a full understanding of the long reach of nucleic acid modifications in cells requires both a global and targeted study of unprecedented scale, which at the moment may well be limited only by technology.

2'-O-methylation within prokaryotic and eukaryotic tRNA inhibits innate immune activation by endosomal Toll-like receptors but does not affect recognition of whole organisms

Bacterial RNA has emerged as an important activator of innate immune responses by stimulating Toll-like receptors TLR7 and TLR8 in humans. Guanosine 2'-O-methylation at position 18 (Gm18) in bacterial tRNA was shown to antagonize tRNA-induced TLR7/8 activation, suggesting a potential role of Gm18 as an immune escape mechanism. This modification also occurs in eukaryotic tRNA, yet a physiological immune function remained to be tested. We therefore set out to investigate the immune modulatory role of Gm18 in both prokaryotic and eukaryotic microorganisms, Escherichia coli and Saccharomyces cerevisiae, and in human cells. Using RiboMethSeq analysis we show that mutation of trmH in E. coli, trm3 in S. cereviase, and CRISPR/Cas9-induced knockout of TARBP1 in H. sapiens results in loss of Gm18 within tRNA. Lack of Gm18 across the kingdoms resulted in increased immunostimulation of peripheral blood mononuclear cells when activated by tRNA preparations. In E. coli, lack of 2'-O-methyltransferase trmH also enhanced immune stimulatory properties by whole cellular RNA. In contrast, lack of Gm18 in yeasts and human cells did not affect immunostimulation by whole RNA preparations. When using live E. coli bacteria, lack of trmH did not affect overall immune stimulation although we detected a defined TLR8/RNA-dependent gene expression signature upon E. coli infection. Together, these results demonstrate that Gm18 is a global immune inhibitory tRNA modification across the kingdoms and contributes to tRNA recognition by innate immune cells, but as an individual modification has insufficient potency to modulate recognition of the investigated microorganisms.