tRNA stabilization by modified nucleotides

Conservation of Archaeal C/D Box sRNA-Guided RNA Modifications

Post-transcriptional modifications fulfill many important roles during ribosomal RNA maturation in all three domains of life. Ribose 2'-O-methylations constitute the most abundant chemical rRNA modification and are, for example, involved in RNA folding and stabilization. In archaea, these modification sites are determined by variable sets of C/D box sRNAs that guide the activity of the rRNA 2'-O-methyltransferase fibrillarin. Each C/D box sRNA contains two guide sequences that can act in coordination to bridge rRNA sequences. Here, we will review the landscape of archaeal C/D box sRNA genes and their target sites. One focus is placed on the apparent accelerated evolution of guide sequences and the varied pairing of the two individual guides, which results in different rRNA modification patterns and RNA chaperone activities.

The expanding world of tRNA modifications and their disease relevance

Transfer RNA (tRNA) is an adapter molecule that links a specific codon in mRNA with its corresponding amino acid during protein synthesis. tRNAs are enzymatically modified post-transcriptionally. A wide variety of tRNA modifications are found in the tRNA anticodon, which are crucial for precise codon recognition and reading frame maintenance, thereby ensuring accurate and efficient protein synthesis. In addition, tRNA-body regions are also frequently modified and thus stabilized in the cell. Over the past two decades, 16 novel tRNA modifications were discovered in various organisms, and the chemical space of tRNA modification continues to expand. Recent studies have revealed that tRNA modifications can be dynamically altered in response to levels of cellular metabolites and environmental stresses. Importantly, we now understand that deficiencies in tRNA modification can have pathological consequences, which are termed 'RNA modopathies'. Dysregulation of tRNA modification is involved in mitochondrial diseases, neurological disorders and cancer.

Posttranscriptional modifications at the 37th position in the anticodon stem-loop of tRNA: structural insights from MD simulations

Transfer RNA (tRNA) is the most diversely modified RNA. Although the strictly conserved purine position 37 in the anticodon stem-loop undergoes modifications that are phylogenetically distributed, we do not yet fully understand the roles of these modifications. Therefore, molecular dynamics simulations are used to provide molecular-level details for how such modifications impact the structure and function of tRNA. A focus is placed on three hypermodified base families that include the parent i6A, t6A, and yW modifications, as well as derivatives. Our data reveal that the hypermodifications exhibit significant conformational flexibility in tRNA, which can be modulated by additional chemical functionalization. Although the overall structure of the tRNA anticodon stem remains intact regardless of the modification considered, the anticodon loop must rearrange to accommodate the bulky, dynamic hypermodifications, which includes changes in the nucleotide glycosidic and backbone conformations, and enhanced or completely new nucleobase-nucleobase interactions compared to unmodified tRNA or tRNA containing smaller (m1G) modifications at the 37th position. Importantly, the extent of the changes in the anticodon loop is influenced by the addition of small functional groups to parent modifications, implying each substituent can further fine-tune tRNA structure. Although the dominant conformation of the ASL is achieved in different ways for each modification, the molecular features of all modified tRNA drive the ASL domain to adopt the functional open-loop conformation. Importantly, the impact of the hypermodifications is preserved in different sequence contexts. These findings highlight the likely role of regulating mRNA structure and translation.

Human transfer RNA modopathies: diseases caused by aberrations in transfer RNA modifications

tRNA molecules are post-transcriptionally modified by tRNA modification enzymes. Although composed of different chemistries, more than 40 types of human tRNA modifications play pivotal roles in protein synthesis by regulating tRNA structure and stability as well as decoding genetic information on mRNA. Many tRNA modifications are conserved among all three kingdoms of life, and aberrations in various human tRNA modification enzymes cause life-threatening diseases. Here, we describe the class of diseases and disorders caused by aberrations in tRNA modifications as 'tRNA modopathies'. Aberrations in over 50 tRNA modification enzymes are associated with tRNA modopathies, which most frequently manifest as dysfunctions of the brain and/or kidney, mitochondrial diseases, and cancer. However, the molecular mechanisms that link aberrant tRNA modifications to human diseases are largely unknown. In this review, we provide a comprehensive compilation of human tRNA modification functions, tRNA modification enzyme genes, and tRNA modopathies, and we summarize the elucidated pathogenic mechanisms underlying several tRNA modopathies. We will also discuss important questions that need to be addressed in order to understand the molecular pathogenesis of tRNA modopathies.

N6-methyladenosine (m6A) is an endogenous A3 adenosine receptor ligand

About 150 post-transcriptional RNA modifications have been identified in all kingdoms of life. During RNA catabolism, most modified nucleosides are resistant to degradation and are released into the extracellular space. In this study, we explored the physiological role of these extracellular modified nucleosides and found that N⁶-methyladenosine (m⁶A), widely recognized as an epigenetic mark in RNA, acts as a ligand for the human adenosine A3 receptor, for which it has greater affinity than unmodified adenosine. We used structural modeling to define the amino acids required for specific binding of m⁶A to the human A3 receptor. We also demonstrated that m⁶A was dynamically released in response to cytotoxic stimuli and facilitated type I allergy in vivo. Our findings implicate m⁶A as a signaling molecule capable of activating G protein-coupled receptors (GPCRs) and triggering pathophysiological responses, a previously unreported property of RNA modifications.

Cell culture NAIL-MS allows insight into human tRNA and rRNA modification dynamics in vivo

Recently, studies about RNA modification dynamics in human RNAs are among the most controversially discussed. As a main reason, we identified the unavailability of a technique which allows the investigation of the temporal processing of RNA transcripts. Here, we present nucleic acid isotope labeling coupled mass spectrometry (NAIL-MS) for efficient, monoisotopic stable isotope labeling in both RNA and DNA in standard cell culture. We design pulse chase experiments and study the temporal placement of modified nucleosides in tRNAPhe and 18S rRNA. In existing RNAs, we observe a time-dependent constant loss of modified nucleosides which is masked by post-transcriptional methylation mechanisms and thus undetectable without NAIL-MS. During alkylation stress, NAIL-MS reveals an adaptation of tRNA modifications in new transcripts but not existing ones. Overall, we present a fast and reliable stable isotope labeling strategy which allows in-depth study of RNA modification dynamics in human cell culture.

tRNA Modifications as a Readout of S and Fe-S Metabolism

Iron-Sulfur (Fe-S) clusters function as core prosthetic groups known to modulate the activity of metalloenzymes, act as trafficking vehicles for biological iron and sulfur, and participate in several intersecting metabolic pathways. The formation of these clusters is initiated by a class of enzymes called cysteine desulfurases, whose primary function is to shuttle sulfur from the amino acid L-cysteine to a variety of sulfur transfer proteins involved in Fe-S cluster synthesis as well as in the synthesis of other thiocofactors. Thus, sulfur and Fe-S cluster metabolism are connected through shared enzyme intermediates, and defects in their associated pathways cause a myriad of pleiotropic phenotypes, which are difficult to dissect. Post-transcriptionally modified transfer RNA (tRNA) represents a large class of analytes whose synthesis often requires the coordinated participation of sulfur transfer and Fe-S enzymes. Therefore, these molecules can be used as biologically relevant readouts for cellular Fe and S status. Methods employing LC-MS technology provide a valuable experimental tool to determine the relative levels of tRNA modification in biological samples and, consequently, to assess genetic, nutritional, and environmental factors modulating reactions dependent on Fe-S clusters. Herein, we describe a robust method for extracting RNA and analytically evaluating the degree of Fe-S-dependent and -independent tRNA modifications via an LC-MS platform.

Northern Blotting Technique for Detection and Expression Analysis of mRNAs and Small RNAs

Northern analysis is a conventional but gold standard method for detection and quantification of gene expression changes. It not only detects the presence of a transcript but also indicates size and relative comparison of transcript abundance on a single membrane. In recent years it has been aptly adapted to validate and study the size and expression of small noncoding RNAs. Here, we describe protocols employed in our laboratory for conventional northern analysis with total RNA/mRNA to study gene expression and validation of small noncoding RNAs using low molecular weight fraction of RNAs. A brief account on the recent advancements for improving the sensitivity and efficiency of northern blot detection is also included in this chapter.

Unnatural bases for recognition of noncoding nucleic acid interfaces

The notion of using synthetic heterocycles instead of the native bases to interface with DNA and RNA has been explored for nearly 60 years. Unnatural bases compatible with the DNA/RNA coding interface have the potential to expand the genetic code and co-opt the machinery of biology to access new macromolecular function; accordingly, this body of research is core to synthetic biology. While much of the literature on artificial bases focuses on code expansion, there is a significant and growing effort on docking synthetic heterocycles to noncoding nucleic acid interfaces; this approach seeks to illuminate major processes of nucleic acids, including regulation of transcription, translation, transport, and transcript lifetimes. These major avenues of research at the coding and noncoding interfaces have in common fundamental principles in molecular recognition. Herein, we provide an overview of foundational literature in biophysics of base recognition and unnatural bases in coding to provide context for the developing area of targeting noncoding nucleic acid interfaces with synthetic bases, with a focus on systems developed through iterative design and biophysical study.

TET-mediated 5-methylcytosine oxidation in tRNA promotes translation

Oxidation of 5-methylcytosine (5mC) in DNA by the ten-eleven translocation (TET) family of enzymes is indispensable for gene regulation in mammals. More recently, evidence has emerged to support a biological function for TET-mediated m5C oxidation in messenger RNA. Here, we describe a previously uncharacterized role of TET-mediated m5C oxidation in transfer RNA (tRNA). We found that the TET-mediated oxidation product 5-hydroxylmethylcytosine (hm5C) is specifically enriched in tRNA inside cells and that the oxidation activity of TET2 on m5C in tRNAs can be readily observed in vitro. We further observed that hm5C levels in tRNA were significantly decreased in Tet2 KO mouse embryonic stem cells (mESCs) in comparison with wild-type mESCs. Reciprocally, induced expression of the catalytic domain of TET2 led to an obvious increase in hm5C and a decrease in m5C in tRNAs relative to uninduced cells. Strikingly, we also show that TET2-mediated m5C oxidation in tRNA promotes translation in vitro. These results suggest TET2 may influence translation through impacting tRNA methylation and reveal an unexpected role for TET enzymes in regulating multiple nodes of the central dogma.

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.

Comparative patterns of modified nucleotides in individual tRNA species from a mesophilic and two thermophilic archaea

To improve and complete our knowledge of archaeal tRNA modification patterns, we have identified and compared the modification pattern (type and location) in tRNAs of three very different archaeal species, Methanococcus maripaludis (a mesophilic methanogen), Pyrococcus furiosus (a hyperthermophile thermococcale), and Sulfolobus acidocaldarius (an acidophilic thermophilic sulfolobale). Most abundant isoacceptor tRNAs (79 in total) for each of the 20 amino acids were isolated by two-dimensional gel electrophoresis followed by in-gel RNase digestions. The resulting oligonucleotide fragments were separated by nanoLC and their nucleotide content analyzed by mass spectrometry (MS/MS). Analysis of total modified nucleosides obtained from complete digestion of bulk tRNAs was also performed. Distinct base- and/or ribose-methylations, cytidine acetylations, and thiolated pyrimidines were identified, some at new positions in tRNAs. Novel, some tentatively identified, modifications were also found. The least diversified modification landscape is observed in the mesophilic Methanococcus maripaludis and the most complex one in Sulfolobus acidocaldarius Notable observations are the frequent occurrence of ac⁴C nucleotides in thermophilic archaeal tRNAs, the presence of m⁷G at positions 1 and 10 in Pyrococcus furiosus tRNAs, and the use of wyosine derivatives at position 37 of tRNAs, especially those decoding U1- and C1-starting codons. These results complete those already obtained by others with sets of archaeal tRNAs from Methanocaldococcus jannaschii and Haloferax volcanii.

TET2 chemically modifies tRNAs and regulates tRNA fragment levels

The ten-eleven translocation 2 (TET2) protein, which oxidizes 5-methylcytosine in DNA, can also bind RNA; however, the targets and function of TET2-RNA interactions in vivo are not fully understood. Using stringent affinity tags introduced at the Tet2 locus, we purified and sequenced TET2-crosslinked RNAs from mouse embryonic stem cells (mESCs) and found a high enrichment for tRNAs. RNA immunoprecipitation with an antibody against 5-hydroxymethylcytosine (hm5C) recovered tRNAs that overlapped with those bound to TET2 in cells. Mass spectrometry (MS) analyses revealed that TET2 is necessary and sufficient for the deposition of the hm5C modification on tRNA. Tet2 knockout in mESCs affected the levels of several small noncoding RNAs originating from TET2-bound tRNAs that were enriched by hm5C immunoprecipitation. Thus, our results suggest a new function of TET2 in promoting the conversion of 5-methylcytosine to hm5C on tRNA and regulating the processing or stability of different classes of tRNA fragments.

Structural and functional insights into Archaeoglobus fulgidus m2G10 tRNA methyltransferase Trm11 and its Trm112 activator

tRNAs play a central role during the translation process and are heavily post-transcriptionally modified to ensure optimal and faithful mRNA decoding. These epitranscriptomics marks are added by largely conserved proteins and defects in the function of some of these enzymes are responsible for neurodevelopmental disorders and cancers. Here, we focus on the Trm11 enzyme, which forms N2-methylguanosine (m2G) at position 10 of several tRNAs in both archaea and eukaryotes. While eukaryotic Trm11 enzyme is only active as a complex with Trm112, an allosteric activator of methyltransferases modifying factors (RNAs and proteins) involved in mRNA translation, former studies have shown that some archaeal Trm11 proteins are active on their own. As these studies were performed on Trm11 enzymes originating from archaeal organisms lacking TRM112 gene, we have characterized Trm11 (AfTrm11) from the Archaeoglobus fulgidus archaeon, which genome encodes for a Trm112 protein (AfTrm112). We show that AfTrm11 interacts directly with AfTrm112 similarly to eukaryotic enzymes and that although AfTrm11 is active as a single protein, its enzymatic activity is strongly enhanced by AfTrm112. We finally describe the first crystal structures of the AfTrm11-Trm112 complex and of Trm11, alone or bound to the methyltransferase inhibitor sinefungin.

Characterization of UVA-Induced Alterations to Transfer RNA Sequences

Ultraviolet radiation (UVR) adversely affects the integrity of DNA, RNA, and their nucleoside modifications. By employing liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based RNA modification mapping approaches, we identified the transfer RNA (tRNA) regions most vulnerable to photooxidation. Photooxidative damage to the anticodon and variable loop regions was consistently observed in both modified and unmodified sequences of tRNA upon UVA (λ 370 nm) exposure. The extent of oxidative damage measured in terms of oxidized guanosine, however, was higher in unmodified RNA compared to its modified version, suggesting an auxiliary role for nucleoside modifications. The type of oxidation product formed in the anticodon stem-loop region varied with the modification type, status, and whether the tRNA was inside or outside the cell during exposure. Oligonucleotide-based characterization of tRNA following UVA exposure also revealed the presence of novel photoproducts and stable intermediates not observed by nucleoside analysis alone. This approach provides sequence-specific information revealing potential hotspots for UVA-induced damage in tRNAs.

m5CPred-SVM: a novel method for predicting m5C sites of RNA

BACKGROUND

As one of the most common post-transcriptional modifications (PTCM) in RNA, 5-cytosine-methylation plays important roles in many biological functions such as RNA metabolism and cell fate decision. Through accurate identification of 5-methylcytosine (m5C) sites on RNA, researchers can better understand the exact role of 5-cytosine-methylation in these biological functions. In recent years, computational methods of predicting m5C sites have attracted lots of interests because of its efficiency and low-cost. However, both the accuracy and efficiency of these methods are not satisfactory yet and need further improvement.

RESULTS

In this work, we have developed a new computational method, m5CPred-SVM, to identify m5C sites in three species, H. sapiens, M. musculus and A. thaliana. To build this model, we first collected benchmark datasets following three recently published methods. Then, six types of sequence-based features were generated based on RNA segments and the sequential forward feature selection strategy was used to obtain the optimal feature subset. After that, the performance of models based on different learning algorithms were compared, and the model based on the support vector machine provided the highest prediction accuracy. Finally, our proposed method, m5CPred-SVM was compared with several existing methods, and the result showed that m5CPred-SVM offered substantially higher prediction accuracy than previously published methods. It is expected that our method, m5CPred-SVM, can become a useful tool for accurate identification of m5C sites.

CONCLUSION

In this study, by introducing position-specific propensity related features, we built a new model, m5CPred-SVM, to predict RNA m5C sites of three different species. The result shows that our model outperformed the existing state-of-art models. Our model is available for users through a web server at https://zhulab.ahu.edu.cn/m5CPred-SVM .

Environmental epitranscriptomics

Chemical modifications of RNA molecules have gained increasing attention since evidence emerged for their substantive roles in a range of biological processes, such as the stability and translation of mRNA transcripts. More than 150 modifications have been identified in different organisms to date, collectively known as the 'epitranscriptome', with 6-methyladenosine (m⁶A), 5-methylcytidine (m⁵C), pseudouridine and N1-methyladenosine (m¹A) the most extensively investigated. Although we are just beginning to elucidate the roles of these modifications in cellular functions, there is already evidence for their dysregulation in diseases such as cancer and neurodevelopmental disorders. There is currently more limited knowledge regarding how environmental exposures affect the epitranscriptome and how this may mediate disease risk, but evidence is beginning to emerge. Here, we review the current evidence for the impact of environmental exposures such as benzo[a]pyrene, bisphenol A, pesticides, metals and nanoparticles upon RNA modifications and the expression of their 'writers' (methyl transferases), 'erasers' (demethylases) and 'readers'. We discuss future directions of the field and identify areas of particular promise and consider the technical challenges that are faced.

Functional Characterization of a Putative RNA Demethylase ALKBH6 in Arabidopsis Growth and Abiotic Stress Responses

RNA methylation and demethylation, which is mediated by RNA methyltransferases (referred to as “writers”) and demethylases (referred to as “erasers”), respectively, are emerging as a key regulatory process in plant development and stress responses. Although several studies have shown that AlkB homolog (ALKBH) proteins are potential RNA demethylases, the function of most ALKBHs is yet to be determined. The Arabidopsis thaliana genome contains thirteen genes encoding ALKBH proteins, the functions of which are largely unknown. In this study, we characterized the function of a potential eraser protein, ALKBH6 (At4g20350), during seed germination and seedling growth in Arabidopsis under abiotic stresses. The seeds of T-DNA insertion alkbh6 knockdown mutants germinated faster than the wild-type seeds under cold, salt, or abscisic acid (ABA) treatment conditions but not under dehydration stress conditions. Although no differences in seedling and root growth were observed between the alkbh6 mutant and wild-type under normal conditions, the alkbh6 mutant showed a much lower survival rate than the wild-type under salt, drought, or heat stress. Cotyledon greening of the alkbh6 mutants was much higher than that of the wild-type upon ABA application. Moreover, the transcript levels of ABA signaling-related genes, including ABI3 and ABI4, were down-regulated in the alkbh6 mutant compared to wild-type plants. Importantly, the ALKBH6 protein had an ability to bind to both m⁶A-labeled and m⁵C-labeled RNAs. Collectively, these results indicate that the potential eraser ALKBH6 plays important roles in seed germination, seedling growth, and survival of Arabidopsis under abiotic stresses.

Naturally occurring modified ribonucleosides

The chemical identity of RNA molecules beyond the four standard ribonucleosides has fascinated scientists since pseudouridine was characterized as the "fifth" ribonucleotide in 1951. Since then, the ever-increasing number and complexity of modified ribonucleosides have been found in viruses and throughout all three domains of life. Such modifications can be as simple as methylations, hydroxylations, or thiolations, complex as ring closures, glycosylations, acylations, or aminoacylations, or unusual as the incorporation of selenium. While initially found in transfer and ribosomal RNAs, modifications also exist in messenger RNAs and noncoding RNAs. Modifications have profound cellular outcomes at various levels, such as altering RNA structure or being essential for cell survival or organism viability. The aberrant presence or absence of RNA modifications can lead to human disease, ranging from cancer to various metabolic and developmental illnesses such as Hoyeraal-Hreidarsson syndrome, Bowen-Conradi syndrome, or Williams-Beuren syndrome. In this review article, we summarize the characterization of all 143 currently known modified ribonucleosides by describing their taxonomic distributions, the enzymes that generate the modifications, and any implications in cellular processes, RNA structure, and disease. We also highlight areas of active research, such as specific RNAs that contain a particular type of modification as well as methodologies used to identify novel RNA modifications. This article is categorized under: RNA Processing > RNA Editing and Modification.

Correlation between gene expression levels under drought stress and synonymous codon usage in rice plant by in-silico study

We studied the correlation of synonymous codon usage (SCU) on gene expression levels under drought stress in rice. Sixty genes related to drought stress (with high, intermediate and low expression) were selected from rice meta-analysis data and various codon usage indices such as the effective number of codon usage (ENC), codon adaptation index (CAI) and relative synonymous codon usage (RSCU) were calculated. We found that in genes highly expressing under drought 1) GC content was higher, 2) ENC value was lower, 3) the preferred codons of some amino acids changed and 4) the RSCU ratio of GC-end codons relative to AT-end codons for 18 amino acids increased significantly compared with those in other genes. We introduce ARSCU as the Average ratio of RSCUs of GC-end codons to AT-end codons in each gene that could significantly separate high-expression genes under drought from low-expression genes. ARSCU is calculated using the program ARSCU-Calculator developed by our group to help predicting expression level of rice genes under drought. An index above ARSCU threshold is expected to indicate that the gene under study may belong to the "high expression group under drought". This information may be applied for codon optimization of genes for rice genetic engineering. To validate these findings, we further used 60 other genes (randomly selected subset of 43233 genes studied for their response to drought stress). ARSCU value was able to predict the level of expression at 88.33% of the cases. Using third set of 60 genes selected amongst high expressing genes not related to drought, only 31.65% of the genes showed ARSCU value of higher than the set threshold. This indicates that the phenomenon we described in this report may be unique for drought related genes. To justify the observed correlation between CUB and high expressing genes under drought, possible role of tRNA post transcriptional modification and tRFs was hypothesized as possible underlying biological mechanism.

Extracurricular Functions of tRNA Modifications in Microorganisms

Transfer RNAs (tRNAs) are essential adaptors that mediate translation of the genetic code. These molecules undergo a variety of post-transcriptional modifications, which expand their chemical reactivity while influencing their structure, stability, and functionality. Chemical modifications to tRNA ensure translational competency and promote cellular viability. Hence, the placement and prevalence of tRNA modifications affects the efficiency of aminoacyl tRNA synthetase (aaRS) reactions, interactions with the ribosome, and transient pairing with messenger RNA (mRNA). The synthesis and abundance of tRNA modifications respond directly and indirectly to a range of environmental and nutritional factors involved in the maintenance of metabolic homeostasis. The dynamic landscape of the tRNA epitranscriptome suggests a role for tRNA modifications as markers of cellular status and regulators of translational capacity. This review discusses the non-canonical roles that tRNA modifications play in central metabolic processes and how their levels are modulated in response to a range of cellular demands.

Transcriptome-Wide 5-Methylcytosine Functional Profiling of Long Non-Coding RNA in Hepatocellular Carcinoma

Background

Growing evidence indicates that methylation status is associated with the pathogenesis of numerous types of cancers. Among these, hepatocellular carcinoma (HCC) is a deadly disease threatening global human health. Although 5-methylcytosine (m5C) has been identified as an important regulatory modification, its distribution in solid tumors, including HCC, remains unclear. The present study aimed to explore the distribution of m5C in HCC.

Materials and Methods

Six pairs of human HCC tissues and adjacent non-tumor tissues were collected to analyze the transcriptome-wide m5C methylation of long non-coding RNA (lncRNA). RNA MeRIP-seq was performed to identify m5C peaks on lncRNA and differences in m5C distribution between HCC and adjacent tissues. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathway analyses were explored to predict the possible roles of m5C.

Results

Using m5C peak sequencing, we observed that a sequence motif was necessary for m5C methylation in HCC lncRNA. Unsupervised hierarchical cluster analysis confirmed that lncRNA m5C methylation occurred more frequently in HCC than adjacent non-tumor tissues. RNA sequencing data demonstrated that more genes were up-regulated by methylation in HCC, while methylation down-regulated more genes in adjacent non-tumor tissues. GO and KEGG pathway analyses revealed that genes having a significant correlation with m5C sites in lncRNA were involved in HCC signaling pathways.

Conclusion

Our results revealed the substantially different amounts and distributions of m5C in HCC compared to adjacent non-tumor tissue. We further predicted the cellular functions in HCC that m5C may participate in to provide evidence implicating m5C lncRNA epigenetic regulation in the tumorigenesis and progression in HCC.

Production and purification of endogenously modified tRNA-derived small RNAs

During particular stress conditions, transfer RNAs (tRNAs) become substrates of stress-induced endonucleases, resulting in the production of distinct tRNA-derived small RNAs (tsRNAs). These small RNAs have been implicated in a wide range of biological processes, but how isoacceptor and even isodecoder-specific tsRNAs act at the molecular level is still poorly understood. Importantly, stress-induced tRNA cleavage affects only a few tRNAs of a given isoacceptor or isodecoder, raising the question as to how such limited molecule numbers could exert measurable biological impact. While the molecular function of individual tsRNAs is likely mediated through association with other molecules, addressing the interactome of specific tsRNAs has only been attempted by using synthetic RNA sequences. Since tRNAs carry post-transcriptional modifications, tsRNAs are likely modified but the extent of their modifications remains largely unknown. Here, we developed a biochemical framework for the production and purification of specific tsRNAs using human cells. Preparative scale purification of tsRNAs from biological sources should facilitate experimentally addressing as to how exactly these small RNAs mediate the multitude of reported molecular functions.

Unraveling the RNA modification code with mass spectrometry

The discovery and analysis of modifications on proteins and nucleic acids has provided functional information that has rapidly accelerated the field of epigenetics. While protein post-translational modifications (PTMs), especially on histones, have been highlighted as critical components of epigenetics, the post-transcriptional modification of RNA has been a subject of more recently emergent interest. Multiple RNA modifications have been known to be present in tRNA and rRNA since the 1960s, but the exploration of mRNA, small RNA, and inducible tRNA modifications remains nascent. Sequencing-based methods have been essential to the field by creating the first epitranscriptome maps of m6A, m5C, hm5C, pseudouridine, and inosine; however, these methods possess significant limitations. Here, we discuss the past, present, and future of the application of mass spectrometry (MS) to the study of RNA modifications.

The Arabidopsis TRM61/TRM6 complex is a bona fide tRNA N1-methyladenosine methyltransferase

tRNA molecules, which contain the most abundant post-transcriptional modifications, are crucial for proper gene expression and protein biosynthesis. Methylation at N1 of adenosine 58 (A58) is critical for maintaining the stability of initiator methionyl-tRNA (tRNAiMet) in bacterial, archaeal, and eukaryotic tRNAs. However, although research has been conducted in yeast and mammals, it remains unclear how A58 in plant tRNAs is modified and involved in development. In this study, we identify the nucleus-localized complex AtTRM61/AtTRM6 in Arabidopsis as tRNA m1A58 methyltransferase. Deficiency or a lack of either AtTRM61 or AtTRM6 leads to embryo arrest and seed abortion. The tRNA m1A level decreases in conditionally complemented Attrm61/LEC1pro::AtTRM61 plants and this is accompanied by reduced levels of tRNAiMet, indicating the importance of the tRNA m1A modification for tRNAiMet stability. Taken together, our results demonstrate that tRNA m1A58 modification is necessary for tRNAiMet stability and is required for embryo development in Arabidopsis.

Structure of tRNA methyltransferase complex of Trm7 and Trm734 reveals a novel binding interface for tRNA recognition

The complex between Trm7 and Trm734 (Trm7–Trm734) from Saccharomyces cerevisiae catalyzes 2′-O-methylation at position 34 in tRNA. We report biochemical and structural studies of the Trm7–Trm734 complex. Purified recombinant Trm7–Trm734 preferentially methylates tRNA^Phe transcript variants possessing two of three factors (Cm32, m¹G37 and pyrimidine34). Therefore, tRNA^Phe, tRNA^Trp and tRNA^Leu are specifically methylated by Trm7–Trm734. We have solved the crystal structures of the apo and S-adenosyl-L-methionine bound forms of Trm7–Trm734. Small angle X-ray scattering reveals that Trm7–Trm734 exists as a hetero-dimer in solution. Trm7 possesses a Rossmann-fold catalytic domain, while Trm734 consists of three WD40 β-propeller domains (termed BPA, BPB and BPC). BPA and BPC form a unique V-shaped cleft, which docks to Trm7. The C-terminal region of Trm7 is required for binding to Trm734. The D-arm of substrate tRNA is required for methylation by Trm7–Trm734. If the D-arm in tRNA^Phe is docked onto the positively charged area of BPB in Trm734, the anticodon-loop is located near the catalytic pocket of Trm7. This model suggests that Trm734 is required for correct positioning of tRNA for methylation. Additionally, a point-mutation in Trm7, which is observed in FTSJ1 (human Trm7 ortholog) of nosyndromic X-linked intellectual disability patients, decreases the methylation activity.

Structure of a tRNA-specific deaminase with compromised deamination activity

Nucleotide 34 in tRNA is extensively modified to ensure translational fidelity and efficacy in cells. The deamination of adenosine at this site catalyzed by the enzyme TadA gives rise to inosine (I), which serves as a typical example of the wobble hypothesis due to its diverse basepairing capability. However, recent studies have shown that tRNAArgACG in Mycoplasma capricolum contains unmodified adenosine, in order to decode the CGG codon. The structural basis behind the poorly performing enzyme M. capricolum TadA (McTadA) is largely unclear. Here we present the structures of the WT and a mutant form of McTadA determined at high resolutions. Through structural comparison between McTadA and other active TadA enzymes as well as modeling efforts, we found that McTadA presents multiple structural conflicts with RNA substrates and thus offered support to previous studies from a structural perspective. These clashes would potentially lead to reduced substrate binding affinity of McTadA, consistent with our in vitro deamination activity and binding assays. To rescue the deamination activity of McTadA, we carried out two rounds of protein engineering through structure-guided design. The unsuccessful attempts of the activity restoration could be attributed to the altered dimer interface and stereo hindrance from the non-catalytic subunit of McTadA, which could be the inevitable outcome of the natural evolution. Our study provides structural insight into an alternative decoding and evolutionary strategy by a compromised TadA enzyme at a molecular level.

The large repertoire of 2'-O-methylation guided by C/D snoRNAs on Trypanosoma brucei rRNA

The parasite Trypanosoma brucei cycles between insect and mammalian hosts, and is the causative agent of sleeping sickness. Here, we performed genome-wide mapping of 2'-O-methylations (Nms) on trypanosome rRNA using three high-throughput sequencing methods; RibOxi-seq, RiboMeth-seq and 2'-OMe-seq. This is the first study using three genome-wide mapping approaches on rRNA from the same species showing the discrepancy among the methods. RibOxi-seq detects all the sites, but RiboMeth-seq is the only method to evaluate the level of a single Nm site. The sequencing revealed at least ninety-nine Nms guided by eighty-five snoRNAs among these thirty-eight Nms are trypanosome specific sites. We present the sequence and target of the C/D snoRNAs guiding on rRNA. This is the highest number of Nms detected to date on rRNA of a single cell parasite. Based on RiboMeth-seq, several Nm sites were found to be differentially regulated at the two stages of the parasite life cycle, the insect procyclic form (PCF) versus the bloodstream form (BSF) in the mammalian host.

Reductive Evolution and Diversification of C5-Uracil Methylation in the Nucleic Acids of Mollicutes

The C5-methylation of uracil to form 5-methyluracil (m5U) is a ubiquitous base modification of nucleic acids. Four enzyme families have converged to catalyze this methylation using different chemical solutions. Here, we investigate the evolution of 5-methyluracil synthase families in Mollicutes, a class of bacteria that has undergone extensive genome erosion. Many mollicutes have lost some of the m5U methyltransferases present in their common ancestor. Cases of duplication and subsequent shift of function are also described. For example, most members of the Spiroplasma subgroup use the ancestral tetrahydrofolate-dependent TrmFO enzyme to catalyze the formation of m5U54 in tRNA, while a TrmFO paralog (termed RlmFO) is responsible for m5U1939 formation in 23S rRNA. RlmFO has replaced the S-adenosyl-L-methionine (SAM)-enzyme RlmD that adds the same modification in the ancestor and which is still present in mollicutes from the Hominis subgroup. Another paralog of this family, the TrmFO-like protein, has a yet unidentified function that differs from the TrmFO and RlmFO homologs. Despite having evolved towards minimal genomes, the mollicutes possess a repertoire of m5U-modifying enzymes that is highly dynamic and has undergone horizontal transfer.

Identification of the 3-amino-3-carboxypropyl (acp) transferase enzyme responsible for acp3U formation at position 47 in Escherichia coli tRNAs

tRNAs from all domains of life contain modified nucleotides. However, even for the experimentally most thoroughly characterized model organism Escherichia coli not all tRNA modification enzymes are known. In particular, no enzyme has been found yet for introducing the acp³U modification at position 47 in the variable loop of eight E. coli tRNAs. Here we identify the so far functionally uncharacterized YfiP protein as the SAM-dependent 3-amino-3-carboxypropyl transferase catalyzing this modification and thereby extend the list of known tRNA modification enzymes in E. coli. Similar to the Tsr3 enzymes that introduce acp modifications at U or m¹Ψ nucleotides in rRNAs this protein contains a DTW domain suggesting that acp transfer reactions to RNA nucleotides are a general function of DTW domain containing proteins. The introduction of the acp³U-47 modification in E. coli tRNAs is promoted by the presence of the m⁷G-46 modification as well as by growth in rich medium. However, a deletion of the enzymes responsible for the modifications at position 46 and 47 in the variable loop of E. coli tRNAs did not lead to a clearly discernible phenotype suggesting that these two modifications play only a minor role in ensuring the proper function of tRNAs in E. coli.

Interchangeable parts: The evolutionarily dynamic tRNA population in plant mitochondria

Transfer RNAs (tRNAs) remain one of the very few classes of genes still encoded in the mitochondrial genome. These key components of the protein translation system must interact with a large enzymatic network of nuclear-encoded gene products to maintain mitochondrial function. Plants have an evolutionarily dynamic mitochondrial tRNA population, including ongoing tRNA gene loss and replacement by both horizontal gene transfer from diverse sources and import of nuclear-expressed tRNAs from the cytosol. Thus, plant mitochondria represent an excellent model for understanding how anciently divergent genes can act as "interchangeable parts" during the evolution of complex molecular systems. In particular, understanding the integration of the mitochondrial translation system with elements of the corresponding machinery used in cytosolic protein synthesis is a key area for eukaryotic cellular evolution. Here, we review the increasingly detailed phylogenetic data about the evolutionary history of mitochondrial tRNA gene loss, transfer, and functional replacement that has created extreme variation in mitochondrial tRNA populations across plant species. We describe emerging tRNA-seq methods with promise for refining our understanding of the expression and subcellular localization of tRNAs. Finally, we summarize current evidence and identify open questions related to coevolutionary changes in nuclear-encoded enzymes that have accompanied turnover in mitochondrial tRNA populations.

RF-PseU: A Random Forest Predictor for RNA Pseudouridine Sites

One of the ubiquitous chemical modifications in RNA, pseudouridine modification is crucial for various cellular biological and physiological processes. To gain more insight into the functional mechanisms involved, it is of fundamental importance to precisely identify pseudouridine sites in RNA. Several useful machine learning approaches have become available recently, with the increasing progress of next-generation sequencing technology; however, existing methods cannot predict sites with high accuracy. Thus, a more accurate predictor is required. In this study, a random forest-based predictor named RF-PseU is proposed for prediction of pseudouridylation sites. To optimize feature representation and obtain a better model, the light gradient boosting machine algorithm and incremental feature selection strategy were used to select the optimum feature space vector for training the random forest model RF-PseU. Compared with previous state-of-the-art predictors, the results on the same benchmark data sets of three species demonstrate that RF-PseU performs better overall. The integrated average leave-one-out cross-validation and independent testing accuracy scores were 71.4% and 74.7%, respectively, representing increments of 3.63% and 4.77% versus the best existing predictor. Moreover, the final RF-PseU model for prediction was built on leave-one-out cross-validation and provides a reliable and robust tool for identifying pseudouridine sites. A web server with a user-friendly interface is accessible at http://148.70.81.170:10228/rfpseu.

2'-O-ribose methylation of transfer RNA promotes recovery from oxidative stress in Saccharomyces cerevisiae

Chemical modifications that regulate protein expression at the translational level are emerging as vital components of the cellular stress response. Transfer RNAs (tRNAs) are significant targets for methyl-based modifications, which are catalyzed by tRNA methyltransferases (Trms). Here, Saccharomyces cerevisiae served as a model eukaryote system to investigate the role of 2'-O-ribose tRNA methylation in the cell's response to oxidative stress. Using 2'-O-ribose deletion mutants for trms 3, 7, 13, and 44, in acute and chronic exposure settings, we demonstrate a broad cell sensitivity to oxidative stress-inducing toxicants (i.e., hydrogen peroxide, rotenone, and acetic acid). A global analysis of hydrogen peroxide-induced tRNA modifications shows a complex profile of decreased, or undetectable, 2'-O-ribose modification events in 2’-O-ribose trm mutant strains, providing a critical link between this type of modification event and Trm status post-exposure. Based on the pronounced oxidative stress sensitivity observed for trm7 mutants, we used a bioinformatic tool to identify transcripts as candidates for regulation by Trm7-catalyzed modifications (i.e., enriched in UUC codons decoded by tRNA^PheGmAA). This screen identified transcripts linked to diverse biological processes that promote cellular recovery after oxidative stress exposure, including DNA repair, chromatin remodeling, and nutrient acquisition (i.e., CRT10, HIR3, HXT2, and GNP1); moreover, these mutants were also oxidative stress-sensitive. Together, these results solidify a role for TRM3, 7, 13, and 44, in the cellular response to oxidative stress, and implicate 2'-O-ribose tRNA modification as an epitranscriptomic strategy for oxidative stress recovery.

Improved RNA modification mapping of cellular non-coding RNAs using C- and U-specific RNases

Locating ribonucleoside modifications within an RNA sequence requires digestion of the RNA into oligoribonucleotides of amenable size for subsequent analysis by LC-MS (liquid chromatography-mass spectrometry). This approach, widely referred to as RNA modification mapping, is facilitated through ribonucleases (RNases) such as T1 (guanosine-specific), U2 (purine-selective) and A (pyrimidine-specific) among others. Sequence coverage by these enzymes depends on positioning of the recognized nucleobase (such as guanine or purine or pyrimidine) in the sequence and its ribonucleotide composition. Using E. coli transfer RNA (tRNA) and ribosomal RNA (rRNA) as model samples, we demonstrate the ability of complementary nucleobase-specific ribonucleases cusativin (C-specific) and MC1 (U-specific) to generate digestion products that facilitate confident mapping of modifications in regions such as G-rich and pyrimidine-rich segments of RNA, and to distinguish C to U sequence differences. These enzymes also increase the number of oligonucleotide digestion products that are unique to a specific RNA sequence. Further, with these additional RNases, multiple modifications can be localized with high confidence in a single set of experiments with minimal dependence on the individual tRNA abundance in a mixture. The sequence overlaps observed with these complementary digestion products and that of RNase T1 improved sequence coverage to 75% or above. A similar level of sequence coverage was also observed for the 2904 nt long 23S rRNA indicating their utility has no dependence on RNA size. Wide-scale adoption of these additional modification mapping tools could help expedite the characterization of modified RNA sequences to understand their structural and functional role in various living systems.

Genome-wide analysis revealed novel molecular features and evolution of Anti-codons in cyanobacterial tRNAs

Transfer RNAs (tRNAs) play important roles to decode the genetic information contained in mRNA in the process of translation. The tRNA molecules possess conserved nucleotides at specific position to regulate the unique function. However, several nucleotides at different position of the tRNA undergo modification to maintain proper stability and function. The major modifications include the presence of pseudouridine (Ψ) residue instead of uridine and the presence of m⁵-methylation sites. We found that, Ψ13 is conserved in D-stem, whereas Ψ38 & Ψ39 were conserved in the anti-codon loop (AL) and anti-codon arm (ACA), respectively. Furthermore, Ψ55 found to be conserved in the Ψ loop. Although, fourteen possible methylation sites can be found in the tRNA, cyanobacterial tRNAs were found to possess conserved G9, m³C₃₂, C₃₆, A₃₇, m⁵C₃₈ and U₅₄ methylation sites. The presence of multiple conserved methylation sites might be responsible for providing necessary stability to the tRNA. The evolutionary study revealed, tRNA^Met and tRNA^Ile were evolved earlier than other tRNA isotypes and their evolution is date back to at least 4000 million years ago. The presence of novel pseudouridination and m⁵-methylation sites in the cyanobacterial tRNAs are of particular interest for basic biology. Further experimental study can delineate their functional significance in protein translation.

TrmB, a tRNA m7G46 methyltransferase, plays a role in hydrogen peroxide resistance and positively modulates the translation of katA and katB mRNAs in Pseudomonas aeruginosa

Cellular response to oxidative stress is a crucial mechanism that promotes the survival of Pseudomonas aeruginosa during infection. However, the translational regulation of oxidative stress response remains largely unknown. Here, we reveal a tRNA modification-mediated translational response to H₂O₂ in P. aeruginosa. We demonstrated that the P. aeruginosa trmB gene encodes a tRNA guanine (46)-N7-methyltransferase that catalyzes the formation of m⁷G46 in the tRNA variable loop. Twenty-three tRNA substrates of TrmB with a guanosine residue at position 46 were identified, including 11 novel tRNA substrates. We showed that loss of trmB had a strong negative effect on the translation of Phe- and Asp-enriched mRNAs. The trmB-mediated m⁷G modification modulated the expression of the catalase genes katA and katB, which are enriched with Phe/Asp codons at the translational level. In response to H₂O₂ exposure, the level of m⁷G modification increased, consistent with the increased translation efficiency of Phe- and Asp-enriched mRNAs. Inactivation of trmB led to decreased KatA and KatB protein abundance and decreased catalase activity, resulting in H₂O₂-sensitive phenotype. Taken together, our observations reveal a novel role of m⁷G46 tRNA modification in oxidative stress response through translational regulation of Phe- and Asp-enriched genes, such as katA and katB.

Biogenesis and functions of aminocarboxypropyluridine in tRNA

Transfer (t)RNAs contain a wide variety of post-transcriptional modifications, which play critical roles in tRNA stability and functions. 3-(3-amino-3-carboxypropyl)uridine (acp3U) is a highly conserved modification found in variable- and D-loops of tRNAs. Biogenesis and functions of acp3U have not been extensively investigated. Using a reverse-genetic approach supported by comparative genomics, we find here that the Escherichia coli yfiP gene, which we rename tapT (tRNA aminocarboxypropyltransferase), is responsible for acp3U formation in tRNA. Recombinant TapT synthesizes acp3U at position 47 of tRNAs in the presence of S-adenosylmethionine. Biochemical experiments reveal that acp3U47 confers thermal stability on tRNA. Curiously, the ΔtapT strain exhibits genome instability under continuous heat stress. We also find that the human homologs of tapT, DTWD1 and DTWD2, are responsible for acp3U formation at positions 20 and 20a of tRNAs, respectively. Double knockout cells of DTWD1 and DTWD2 exhibit growth retardation, indicating that acp3U is physiologically important in mammals.

Interplay between substrate recognition, 5' end tRNA processing and methylation activity of human mitochondrial RNase P

Human mitochondrial ribonuclease P (mtRNase P) is an essential three-protein complex that catalyzes the 5' end maturation of mitochondrial precursor tRNAs (pre-tRNAs). Mitochondrial RNase P Protein 3 (MRPP3), a protein-only RNase P (PRORP), is the nuclease component of the mtRNase P complex and requires a two-protein S-adenosyl-methionine (SAM)-dependent methyltransferase MRPP1/2 subcomplex to function. Dysfunction of mtRNase P is linked to several human mitochondrial diseases, such as mitochondrial myopathies. Despite its central role in mitochondrial RNA processing, little is known about how the protein subunits of mtRNase P function synergistically. Here, we use purified mtRNase P to demonstrate that mtRNase P recognizes, cleaves, and methylates some, but not all, mitochondrial pre-tRNAs in vitro. Additionally, mtRNase P does not process all mitochondrial pre-tRNAs uniformly, suggesting the possibility that some pre-tRNAs require additional factors to be cleaved in vivo. Consistent with this, we found that addition of the TRMT10C (MRPP1) cofactor SAM enhances the ability of mtRNase P to bind and cleave some mitochondrial pre-tRNAs. Furthermore, the presence of MRPP3 can enhance the methylation activity of MRPP1/2. Taken together, our data demonstrate that the subunits of mtRNase P work together to efficiently recognize, process, and methylate human mitochondrial pre-tRNAs.

Serratia liquefaciens FG3 isolated from a metallophyte plant sheds light on the evolution and mechanisms of adaptive traits in extreme environments

Serratia liquefaciens strain FG3 (SlFG3), isolated from the flower of Stachytarpheta glabra in the Brazilian ferruginous fields, has distinctive genomic, adaptive, and biotechnological potential. Herein, using a combination of genomics and molecular approaches, we unlocked the evolution of the adaptive traits acquired by S1FG3, which exhibits the second largest chromosome containing the largest conjugative plasmids described for Serratia. Comparative analysis revealed the presence of 18 genomic islands and 311 unique protein families involved in distinct adaptive features. S1FG3 has a diversified repertoire of genes associated with Nonribosomal peptides (NRPs/PKS), a complete and functional cluster related to cellulose synthesis, and an extensive and functional repertoire of oxidative metabolism genes. In addition, S1FG3 possesses a complete pathway related to protocatecuate and chloroaromatic degradation, and a complete repertoire of genes related to DNA repair and protection that includes mechanisms related to UV light tolerance, redox process resistance, and a laterally acquired capacity to protect DNA using phosphorothioation. These findings summarize that SlFG3 is well-adapted to different biotic and abiotic stress situations imposed by extreme conditions associated with ferruginous fields, unlocking the impact of the lateral gene transfer to adjust the genome for extreme environments, and providing insight into the evolution of prokaryotes.

Reengineering a tRNA Methyltransferase To Covalently Capture New RNA Substrates

Covalent RNA modifications can alter the function and metabolism of RNA transcripts. Altering the RNA substrate specificities of the enzymes that install these modifications can provide powerful tools to study and manipulate RNA. To develop new tools and probe the plasticity of the substrate specificity of one of these enzymes, we examined the engineerability of the uridine-54 tRNA methyltransferase, TrmA. Starting from a mutant that remains covalently bound to its substrate RNA (E358Q, TrmA*), we were able to use both rational design and a high-throughput sequencing assay to examine the RNA substrates of TrmA*. Although rational engineering substantially changed TrmA* specificity, the rationally designed substrate mutants we developed still retained activity with the wild-type protein. Using high-throughput substrate screening of additional TrmA* mutants, we identified a triple mutant of the substrate RNA (C56A;A58G;C60U) that is efficiently trapped by a TrmA* double mutant (E49R;R51E) but not by the wild-type TrmA*. This work establishes a foundation for using protein engineering to reconfigure substrate specificities of RNA-modifying enzymes and covalently trap RNAs with engineered proteins.

RNAm5CPred: Prediction of RNA 5-Methylcytosine Sites Based on Three Different Kinds of Nucleotide Composition

5-methylcytosine (m5C) is one of the most common and abundant post-transcriptional modifications (PTCMs) in RNA. Recent studies showed that m5C plays important roles in many biological functions such as RNA metabolism and cell fate decision. Because most experimental methods that determine m5C sites across the transcriptome are time-consuming and expensive, it is urgent to develop accurate computational methods to identify m5C sites effectively. A benchmark dataset is important for developing and evaluating computational methods. In this work, we constructed four different datasets according to the data redundancy and imbalance. Based on these datasets, we generated three different kinds of features, i.e., KNFs (K-nucleotide frequencies), KSNPFs (K-spaced nucleotide pair frequencies), and pseDNC (pseudo-dinucleotide composition), and then used a support vector machine (SVM) to build our models. Based on the imbalanced and nonredundant dataset, Met935, we extensively studied the three kinds of features and determined an optimal combination of the features. Based on the feature combination, we built models on the three different datasets and compared them with state-of-the-art models. According to the predictive results of the stringent jackknife test, the models based on the three features, 4NF, 1SNPF, and pseDNC, are superior or comparable to other methods. To determine the best model between the models based on the imbalanced dataset Met935 and the balanced dataset Met240, we further evaluated the two models on an independent test set Test1157. Our results demonstrate that the model based on the balanced dataset Met240 achieved the highest recall (68.79%) and the highest Matthews correlation coefficient (MCC) (0.154). In addition, the model is also superior to other state-of-the-art methods according to the integrated parameter MCC on the independent test set. Thus, we selected the model based on Met240 as our final model, which was named RNAm5CPred. In addition, a web server for RNAm5CPred (http://zhulab.ahu.edu.cn/RNAm5CPred/) has been provided to facilitate experimental research.

Structural Elucidation of Bisulfite Adducts to Pseudouridine That Result in Deletion Signatures during Reverse Transcription of RNA

The recent report of RBS-Seq to map simultaneously the epitranscriptomic modifications N¹-methyladenosine, 5-methylcytosine, and pseudouridine (Ψ) via bisulfite treatment of RNA provides a key advance to locate these important modifications. The locations of Ψ were found by a deletion signature generated during cDNA synthesis after bisulfite treatment for which the chemical details of the reaction are poorly understood. In the present work, the bisulfite reaction with Ψ was explored to identify six isomers of bisulfite adducted to Ψ. We found four of these adducts involved the heterocyclic ring, similar to the reaction with other pyrimidines. The remaining two adducts were bonded to the 1' carbon, which resulted in opening of the ribose ring. The utilization of complementary 1D- and 2D-NMR, Raman, and electronic circular dichroism spectroscopies led to the assignment of the two ribose adducts being the constitutional isomers of an S- and an O-adduct of bisulfite to the ribose, and these are the final products after heating. A mechanistic proposal is provided to rationalize chemically the formation and stereochemistries of all six isomeric bisulfite adducts to Ψ; conversion of intermediate adducts to the two final products is proposed to involve E2, S_N2', and [2,3]-sigmatropic shift reactions. Lastly, a synthetic RNA template with Ψ at a known location was treated with bisulfite, leading to a deletion signature after reverse transcription, supporting the RBS-Seq report. This classical bisulfite reaction used for epigenomic and epitranscriptomic sequencing diverges from the C nucleoside Ψ to form stable bisulfite end products that yield signatures for next-generation sequencing.

Genome-wide Translation Profiling by Ribosome-Bound tRNA Capture

In the ribosome complex, tRNA is a critical element of mRNA translation. A rich repertoire of cell regulation is hypothesized to occur during the recruitment of specific tRNAs in polypeptide formation. However, this basic question in nascent chain biology remains unaddressed due to the lack of technologies to report the complete tRNA complement inside ribosomes during active translation. Here, we characterize a technique for profiling ribosome-embedded tRNA and their modifications. With this method, we generated a comprehensive survey of the quantity and quality of intra-ribosomal tRNAs. In cells under environmental stress, we show that methionine tRNA inside ribosomes is a robust biomarker for the impairment of translation initiation or elongation steps. Concurrent tRNA/mRNA ribosome profiling revealed a stress-dependent incorporation of damaged and uncharged tRNAs into ribosomes causing translation arrest. Thus, tRNA ribosome profiling can provide insights on translation control mechanisms in diverse biological contexts.

Archaeal NSUN6 catalyzes m5C72 modification on a wide-range of specific tRNAs

Human NOL1/NOP2/Sun RNA methyltransferase family member 6 (hNSun6) generates 5-methylcytosine (m5C) at C72 of four specific tRNAs, and its homologs are present only in higher eukaryotes and hyperthermophilic archaea. Archaeal NSun6 homologs possess conserved catalytic residues, but have distinct differences in their RNA recognition motifs from eukaryotic NSun6s. Until now, the biochemical properties and functions of archaeal NSun6 homologs were unknown. In archaeon Pyrococcus horikoshii OT3, the gene encoding the NSun6 homolog is PH1991. We demonstrated that the PH1991 protein could catalyze m5C72 formation on some specific PhtRNAs in vitro and was thus named as PhNSun6. Remarkably, PhNSun6 has a much wider range of tRNA substrates than hNSun6, which was attributed to its tRNA substrate specificity. The mechanism was further elucidated using biochemical and crystallographic experiments. Structurally, the binding pocket for nucleotide 73 in PhNSun6 is specific to accommodate U73 or G73-containing PhtRNAs. Furthermore, PhNSun6 lacks the eukaryotic NSun6-specific Lys-rich loop, resulting in the non-recognition of D-stem region by PhNSun6. Functionally, the m5C72 modification could slightly promote the thermal stability of PhtRNAs, but did not affect the amino acid accepting activity of PhtRNAs.

Exposure to environmental toxicants reduces global N6-methyladenosine RNA methylation and alters expression of RNA methylation modulator genes

The epitranscriptome comprises more than 100 forms of RNA modifications. Of these, N6-methyladenosine (m⁶A) is the most abundantform of RNA methylation, with roles in modulating mRNA transcript processing and regulation. The aims of the study weretoexamine changes inm⁶A RNA methylation in A549 lung epithelial cells in response to environmental toxicants, anddifferential gene expression of m⁶A modulator genes ('readers', 'writers' and 'erasers') in human subjects exposed toparticulate matter (PM) and in lung cancer tissueusing publicly-available microarray datasets. Global m⁶A methylation levelsweremeasured in total RNA after exposuretotwo carcinogens (PM and sodium arsenite) for 24- and 48-h, and totwo endocrine disruptors (bisphenol A and vinclozolin)for 24-h.Global m⁶A methylation level significantly decreased with exposure to >62 μg/mlPM, >1 μM sodium arsenite, >1  μM bisphenol A (BPA), and0.1  μM vinclozolin. In an analysis of a published dataset derived from a population study, we observed that m⁶A writers (METTL3 and WTAP), erasers (FTO and ALKBH5) and readers (HNRPC) showed significantly higher expression among participants in the high-PM_2.5exposure group compared to those in the low-exposure control group (all p < 0.05). Further, the m⁶A writer METTL3shows reduced expression in lung tumors in comparison to normal lung epithelia (p < 0.0001). Our findings reveal that m⁶A RNA methylation can be modified by exposure to environmental toxicants, and exposure to particulate matter is associated with differential expression level of m⁶A RNA methylation modification machinery.

To be or not to be modified: Miscellaneous aspects influencing nucleotide modifications in tRNAs

Transfer RNAs (tRNAs) are essential components of the cellular protein synthesis machineries, but are also implicated in many roles outside translation. To become functional, tRNAs, initially transcribed as longer precursor tRNAs, undergo a tightly controlled biogenesis process comprising the maturation of their extremities, removal of intronic sequences if present, addition of the 3'-CCA amino-acid accepting sequence, and aminoacylation. In addition, the most impressive feature of tRNA biogenesis consists in the incorporation of a large number of posttranscriptional chemical modifications along its sequence. The chemical nature of these modifications is highly diverse, with more than hundred different modifications identified in tRNAs to date. All functions of tRNAs in cells are controlled and modulated by modifications, making the understanding of the mechanisms that determine and influence nucleotide modifications in tRNAs an essential point in tRNA biology. This review describes the different aspects that determine whether a certain position in a tRNA molecule is modified or not. We describe how sequence and structural determinants, as well as the presence of prior modifications control modification processes. We also describe how environmental factors and cellular stresses influence the level and/or the nature of certain modifications introduced in tRNAs, and report situations where these dynamic modulations of tRNA modification levels are regulated by active demodification processes. © 2019 IUBMB Life, 71(8):1126-1140, 2019.

Oligonucleotide-mediated tRNA sequestration enables one-pot sense codon reassignment in vitro

Sense codon reassignment to unnatural amino acids (uAAs) represents a powerful approach for introducing novel properties into polypeptides. The main obstacle to this approach is competition between the native isoacceptor tRNA(s) and orthogonal tRNA(s) for the reassigned codon. While several chromatographic and enzymatic procedures for selective deactivation of tRNA isoacceptors in cell-free translation systems exist, they are complex and not scalable. We designed a set of tRNA antisense oligonucleotides composed of either deoxy-, ribo- or 2′-O-methyl ribonucleotides and tested their ability to efficiently complex tRNAs of choice. Methylated oligonucleotides targeting sequence between the anticodon and variable loop of tRNA^SerGCU displayed subnanomolar binding affinity with slow dissociation kinetics. Such oligonucleotides efficiently and selectively sequestered native tRNA^SerGCU directly in translation-competent Escherichia coli S30 lysate, thereby, abrogating its translational activity and liberating the AGU/AGC codons. Expression of eGFP protein from the template harboring a single reassignable AGU codon in tRNA^SerGCU-depleted E. coli lysate allowed its homogeneous modification with n-propargyl-l-lysine or p-azido-l-phenylalanine. The strategy developed here is generic, as demonstrated by sequestration of tRNA^ArgCCU isoacceptor in E. coli translation system. Furthermore, this method is likely to be species-independent and was successfully applied to the eukaryotic Leishmania tarentolae in vitro translation system. This approach represents a new direction in genetic code reassignment with numerous practical applications.

Absolute Quantification of Noncoding RNA by Microscale Thermophoresis

Accurate quantification of the copy numbers of noncoding RNA has recently emerged as an urgent problem, with impact on fields such as RNA modification research, tissue differentiation, and others. Herein, we present a hybridization-based approach that uses microscale thermophoresis (MST) as a very fast and highly precise readout to quantify, for example, single tRNA species with a turnaround time of about one hour. We developed MST to quantify the effect of tRNA toxins and of heat stress and RNA modification on single tRNA species. A comparative analysis also revealed significant differences to RNA-Seq-based quantification approaches, strongly suggesting a bias due to tRNA modifications in the latter. Further applications include the quantification of rRNA as well as of polyA levels in cellular RNA.