MODOMICS: a database of RNA modification pathways. 2021 update

Coordination of RNA modifications in the brain and beyond

Gene expression regulation is a critical process throughout the body, especially in the nervous system. One mechanism by which biological systems regulate gene expression is via enzyme-mediated RNA modifications, also known as epitranscriptomic regulation. RNA modifications, which have been found on nearly all RNA species across all domains of life, are chemically diverse covalent modifications of RNA nucleotides and represent a robust and rapid mechanism for the regulation of gene expression. Although numerous studies have been conducted regarding the impact that single modifications in single RNA molecules have on gene expression, emerging evidence highlights potential crosstalk between and coordination of modifications across RNA species. These potential coordination axes of RNA modifications have emerged as a new direction in the field of epitranscriptomic research. In this review, we will highlight several examples of gene regulation via RNA modification in the nervous system, followed by a summary of the current state of the field of RNA modification coordination axes. In doing so, we aim to inspire the field to gain a deeper understanding of the roles of RNA modifications and coordination of these modifications in the nervous system.

Characterization of epitranscriptome reader proteins experimentally and in silico: Current knowledge and future perspectives beyond the YTH domain

To date, over 150 chemical modifications to the four canonical RNA bases have been discovered, known collectively as the epitranscriptome. Many of these modifications have been implicated in a variety of cellular processes and disease states. Additional work has been done to identify proteins known as "readers" that selectively interact with RNAs that contain specific chemical modifications. Protein interactomes with N6-methyladenosine (m⁶A), N1-methyladenosine (m¹A), N5-methylcytosine (m⁵C), and 8-oxo-7,8-dihydroguanosine (8-oxoG) have been determined, mainly through experimental advances in proteomics techniques. However, relatively few proteins have been confirmed to bind directly to RNA containing these modifications. Furthermore, for many of these protein readers, the exact binding mechanisms as well as the exclusivity for recognition of modified RNA species remain elusive, leading to questions regarding their roles within different cellular processes. In the case of the YT-521B homology (YTH) family of proteins, both experimental and in silico techniques have been leveraged to provide valuable biophysical insights into the mechanisms of m⁶A recognition at atomic resolution. To date, the YTH family is one of the best characterized classes of readers. Here, we review current knowledge about epitranscriptome recognition of the YTH domain proteins from previously published experimental and computational studies. We additionally outline knowledge gaps for proteins beyond the well-studied human YTH domains and the current in silico techniques and resources that can enable investigation of protein interactions with modified RNA outside of the YTH-m⁶A context.

Small RNA modifications: regulatory molecules and potential applications

Small RNAs (also referred to as small noncoding RNAs, sncRNA) are defined as polymeric ribonucleic acid molecules that are less than 200 nucleotides in length and serve a variety of essential functions within cells. Small RNA species include microRNA (miRNA), PIWI-interacting RNA (piRNA), small interfering RNA (siRNA), tRNA-derived small RNA (tsRNA), etc. Current evidence suggest that small RNAs can also have diverse modifications to their nucleotide composition that affect their stability as well as their capacity for nuclear export, and these modifications are relevant to their capacity to drive molecular signaling processes relevant to biogenesis, cell proliferation and differentiation. In this review, we highlight the molecular characteristics and cellular functions of small RNA and their modifications, as well as current techniques for their reliable detection. We also discuss how small RNA modifications may be relevant to the clinical applications for the diagnosis and treatment of human health conditions such as cancer.

mRNA Regulation by RNA Modifications

Chemical modifications on mRNA represent a critical layer of gene expression regulation. Research in this area has continued to accelerate over the last decade, as more modifications are being characterized with increasing depth and breadth. mRNA modifications have been demonstrated to influence nearly every step from the early phases of transcript synthesis in the nucleus through to their decay in the cytoplasm, but in many cases, the molecular mechanisms involved in these processes remain mysterious. Here, we highlight recent work that has elucidated the roles of mRNA modifications throughout the mRNA life cycle, describe gaps in our understanding and remaining open questions, and offer some forward-looking perspective on future directions in the field.

PARP-dependent and NAT10-independent acetylation of N4-cytidine in RNA appears in UV-damaged chromatin

RNA modifications have been known for many years, but their function has not been fully elucidated yet. For instance, the regulatory role of acetylation on N4-cytidine (ac4C) in RNA can be explored not only in terms of RNA stability and mRNA translation but also in DNA repair. Here, we observe a high level of ac4C RNA at DNA lesions in interphase cells and irradiated cells in telophase. Ac4C RNA appears in the damaged genome from 2 to 45 min after microirradiation. However, RNA cytidine acetyltransferase NAT10 did not accumulate to damaged sites, and NAT10 depletion did not affect the pronounced recruitment of ac4C RNA to DNA lesions. This process was not dependent on the G1, S, and G2 cell cycle phases. In addition, we observed that the PARP inhibitor, olaparib, prevents the recruitment of ac4C RNA to damaged chromatin. Our data imply that the acetylation of N4-cytidine, especially in small RNAs, has an important role in mediating DNA damage repair. Ac4C RNA likely causes de-condensation of chromatin in the vicinity of DNA lesions, making it accessible for other DNA repair factors involved in the DNA damage response. Alternatively, RNA modifications, including ac4C, could be direct markers of damaged RNAs.

A tRNA modification in Mycobacterium tuberculosis facilitates optimal intracellular growth

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

Influence of 2'-Modifications (O-Methylation, Fluorination, and Stereochemical Inversion) on the Base Pairing Energies of Protonated Cytidine Nucleoside Analogue Base Pairs: Implications for the Stabilities of i-Motif Structures

Naturally occurring and chemically engineered modifications are among the most powerful strategies explored for fine-tuning the conformational characteristics and intrinsic stability of nucleic acids topologies. Modifications at the 2'-position of the ribose or 2'-deoxyribose moieties differentiate nucleic acid structures and have a significant impact on their electronic properties and base-pairing interactions. 2'-O-Methylation, a common post-transcriptional modification of tRNA, is directly involved in modulating specific anticodon-codon base-pairing interactions. 2'-Fluorinated and arabino nucleosides possess novel and beneficial medicinal properties and find use as therapeutics for treating viral diseases and cancer. However, the potential to deploy 2'-modified cytidine chemistries for tuning i-motif stability is largely unknown. To address this knowledge gap, the effects of 2'-modifications including O-methylation, fluorination, and stereochemical inversion on the base-pairing interactions of protonated cytidine nucleoside analogue base pairs, the core stabilizing interactions of i-motif structures, are examined using complementary threshold collision-induced dissociation techniques and computational methods. The 2'-modified cytidine nucleoside analogues investigated here include 2'-O-methylcytidine, 2'-fluoro-2'-deoxycytidine, arabinofuranosylcytosine, 2'-fluoro-arabinofuranosylcytosine, and 2',2'-difluoro-2'-deoxycytidine. All five 2'-modifications examined here are found to enhance the base-pairing interactions relative to the canonical DNA and RNA cytidine nucleosides with the greatest enhancements arising from 2'-O-methylation and 2',2'-difluorination, suggesting that these modifications should well be tolerated in the narrow grooves of i-motif conformations.

Human SAMD9 is a poxvirus-activatable anticodon nuclease inhibiting codon-specific protein synthesis

As a defense strategy against viruses or competitors, some microbes use anticodon nucleases (ACNases) to deplete essential tRNAs, effectively halting global protein synthesis. However, this mechanism has not been observed in multicellular eukaryotes. Here, we report that human SAMD9 is an ACNase that specifically cleaves phenylalanine tRNA (tRNAPhe), resulting in codon-specific ribosomal pausing and stress signaling. While SAMD9 ACNase activity is normally latent in cells, it can be activated by poxvirus infection or rendered constitutively active by SAMD9 mutations associated with various human disorders, revealing tRNAPhe depletion as an antiviral mechanism and a pathogenic condition in SAMD9 disorders. We identified the N-terminal effector domain of SAMD9 as the ACNase, with substrate specificity primarily determined by a eukaryotic tRNAPhe-specific 2'-O-methylation at the wobble position, making virtually all eukaryotic tRNAPhe susceptible to SAMD9 cleavage. Notably, the structure and substrate specificity of SAMD9 ACNase differ from known microbial ACNases, suggesting convergent evolution of a common immune defense strategy targeting tRNAs.

Epitranscriptional Regulation: From the Perspectives of Cardiovascular Bioengineering

The central dogma of gene expression involves DNA transcription to RNA and RNA translation into protein. As key intermediaries and modifiers, RNAs undergo various forms of modifications such as methylation, pseudouridylation, deamination, and hydroxylation. These modifications, termed epitranscriptional regulations, lead to functional changes in RNAs. Recent studies have demonstrated crucial roles for RNA modifications in gene translation, DNA damage response, and cell fate regulation. Epitranscriptional modifications play an essential role in development, mechanosensing, atherogenesis, and regeneration in the cardiovascular (CV) system, and their elucidation is critically important to understanding the molecular mechanisms underlying CV physiology and pathophysiology. This review aims at providing biomedical engineers with an overview of the epitranscriptome landscape, related key concepts, recent findings in epitranscriptional regulations, and tools for epitranscriptome analysis. The potential applications of this important field in biomedical engineering research are discussed.

Underexplored reciprocity between genome-wide methylation status and long non-coding RNA expression reflected in breast cancer research: potential impacts for the disease management in the framework of 3P medicine

Breast cancer (BC) is the most common female malignancy reaching a pandemic scale worldwide. A comprehensive interplay between genetic alterations and shifted epigenetic regions synergistically leads to disease development and progression into metastatic BC. DNA and histones methylations, as the most studied epigenetic modifications, represent frequent and early events in the process of carcinogenesis. To this end, long non-coding RNAs (lncRNAs) are recognized as potent epigenetic modulators in pathomechanisms of BC by contributing to the regulation of DNA, RNA, and histones' methylation. In turn, the methylation status of DNA, RNA, and histones can affect the level of lncRNAs expression demonstrating the reciprocity of mechanisms involved. Furthermore, lncRNAs might undergo methylation in response to actual medical conditions such as tumor development and treated malignancies. The reciprocity between genome-wide methylation status and long non-coding RNA expression levels in BC remains largely unexplored. Since the bio/medical research in the area is, per evidence, strongly fragmented, the relevance of this reciprocity for BC development and progression has not yet been systematically analyzed. Contextually, the article aims at:consolidating the accumulated knowledge on both-the genome-wide methylation status and corresponding lncRNA expression patterns in BC andhighlighting the potential benefits of this consolidated multi-professional approach for advanced BC management. Based on a big data analysis and machine learning for individualized data interpretation, the proposed approach demonstrates a great potential to promote predictive diagnostics and targeted prevention in the cost-effective primary healthcare (sub-optimal health conditions and protection against the health-to-disease transition) as well as advanced treatment algorithms tailored to the individualized patient profiles in secondary BC care (effective protection against metastatic disease). Clinically relevant examples are provided, including mitochondrial health control and epigenetic regulatory mechanisms involved.

Nanopore sequencing of internal 2'-PO4 modifications installed by RNA repair

Ligation by plant and fungal RNA ligases yields an internal 2'-phosphate group on each RNA ligation product. In budding yeast, this covalent mark occurs at the splice junction of two targets of ligation: intron-containing tRNAs and the messenger RNA HAC1 The repertoire of RNA molecules repaired by RNA ligation has not been explored due to a lack of unbiased approaches for identifying RNA ligation products. Here, we define several unique signals produced by 2'-phosphorylated RNAs during nanopore sequencing. A 2'-phosphate at the splice junction of HAC1 mRNA inhibits 5' → 3' degradation, enabling detection of decay intermediates in yeast RNA repair mutants by nanopore sequencing. During direct RNA sequencing, intact 2'-phosphorylated RNAs on HAC1 and tRNAs produce diagnostic changes in nanopore current properties and base calling features, including stalls produced as the modified RNA translocates through the nanopore motor protein. These approaches enable directed studies to identify novel RNA repair events in other contexts.

Disruption of the mouse liver epitranscriptome by long-term aroclor 1260 exposure

Chronic environmental exposure to polychlorinated biphenyls (PCBs) is associated with non-alcoholic fatty liver disease (NAFLD) and exacerbated by a high fat diet (HFD). Here, chronic (34 wks.) exposure of low fat diet (LFD)-fed male mice to Aroclor 1260 (Ar1260), a non-dioxin-like (NDL) mixture of PCBs, resulted in steatohepatitis and NAFLD. Twelve hepatic RNA modifications were altered with Ar1260 exposure including reduced abundance of 2'-O-methyladenosine (Am) and N(6)-methyladenosine (m6A), in contrast to increased Am in the livers of HFD-fed, Ar1260-exposed mice reported previously. Differences in 13 RNA modifications between LFD- and HFD- fed mice, suggest that diet regulates the liver epitranscriptome. Integrated network analysis of epitranscriptomic modifications identified a NRF2 (Nfe2l2) pathway in the chronic, LFD, Ar1260-exposed livers and an NFATC4 (Nfatc4) pathway for LFD- vs. HFD-fed mice. Changes in protein abundance were validated. The results demonstrate that diet and Ar1260 exposure alter the liver epitranscriptome in pathways associated with NAFLD.

Vital roles of m5C RNA modification in cancer and immune cell biology

RNA modification plays an important role in epigenetics at the posttranscriptional level, and 5-methylcytosine (m5C) has attracted increasing attention in recent years due to the improvement in RNA m5C site detection methods. By influencing transcription, transportation and translation, m5C modification of mRNA, tRNA, rRNA, lncRNA and other RNAs has been proven to affect gene expression and metabolism and is associated with a wide range of diseases, including malignant cancers. RNA m5C modifications also substantially impact the tumor microenvironment (TME) by targeting different groups of immune cells, including B cells, T cells, macrophages, granulocytes, NK cells, dendritic cells and mast cells. Alterations in immune cell expression, infiltration and activation are highly linked to tumor malignancy and patient prognosis. This review provides a novel and holistic examination of m5C-mediated cancer development by examining the exact mechanisms underlying the oncogenicity of m5C RNA modification and summarizing the biological effects of m5C RNA modification on tumor cells as well as immune cells. Understanding methylation-related tumorigenesis can provide useful insights for the diagnosis as well as the treatment of cancer.

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

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

MRT-ModSeq - Rapid detection of RNA modifications with MarathonRT

Chemical modifications are essential regulatory elements that modulate the behavior and function of cellular RNAs. Despite recent advances in sequencing-based RNA modification mapping, methods combining accuracy and speed are still lacking. Here, we introduce MRT- ModSeq for rapid, simultaneous detection of multiple RNA modifications using MarathonRT. MRT-ModSeq employs distinct divalent cofactors to generate 2-D mutational profiles that are highly dependent on nucleotide identity and modification type. As a proof of concept, we use the MRT fingerprints of well-studied rRNAs to implement a general workflow for detecting RNA modifications. MRT-ModSeq rapidly detects positions of diverse modifications across a RNA transcript, enabling assignment of m1acp3Y, m1A, m3U, m7G and 2'-OMe locations through mutation-rate filtering and machine learning. m1A sites in sparsely modified targets, such as MALAT1 and PRUNE1 could also be detected. MRT-ModSeq can be trained on natural and synthetic transcripts to expedite detection of diverse RNA modification subtypes across targets of interest.

Evaluation and development of deep neural networks for RNA 5-Methyluridine classifications using autoBioSeqpy

Post-transcriptionally RNA modifications, also known as the epitranscriptome, play crucial roles in the regulation of gene expression during development. Recently, deep learning (DL) has been employed for RNA modification site prediction and has shown promising results. However, due to the lack of relevant studies, it is unclear which DL architecture is best suited for some pyrimidine modifications, such as 5-methyluridine (m⁵U). To fill this knowledge gap, we first performed a comparative evaluation of various commonly used DL models for epigenetic studies with the help of autoBioSeqpy. We identified optimal architectural variations for m⁵U site classification, optimizing the layer depth and neuron width. Second, we used this knowledge to develop Deepm5U, an improved convolutional-recurrent neural network that accurately predicts m⁵U sites from RNA sequences. We successfully applied Deepm5U to transcriptomewide m⁵U profiling data across different sequencing technologies and cell types. Third, we showed that the techniques for interpreting deep neural networks, including LayerUMAP and DeepSHAP, can provide important insights into the internal operation and behavior of models. Overall, we offered practical guidance for the development, benchmark, and analysis of deep learning models when designing new algorithms for RNA modifications.

Detection, regulation, and functions of RNA N6-methyladenosine modification in plants

N⁶-Methyladenosine (m⁶A) is the most abundant internal chemical modification in eukaryotic mRNA and plays important roles in gene expression regulation, including transcriptional and post-transcriptional regulation. m⁶A is a reversible modification that is installed, removed, and recognized by methyltransferases (writers), demethylases (erasers), and m⁶A-binding proteins (readers), respectively. Recently, the breadth of research on m⁶A in plants has expanded, and the vital roles of m⁶A in plant development, biotic and abiotic stress responses, and crop trait improvement have been investigated. In this review, we discuss recent developments in research on m⁶A and highlight the detection methods, distribution, regulatory proteins, and molecular and biological functions of m⁶A in plants. We also offer some perspectives on future investigations, providing direction for subsequent research on m⁶A in plants.

Identification of Up47 in three thermophilic archaea, one mesophilic archaeon, and one hyperthermophilic bacterium

Analysis of the profile of the tRNA modifications in several Archaea allowed us to observe a novel modified uridine in the V-loop of several tRNAs from two species: Pyrococcus furiosus and Sulfolobus acidocaldarius Recently, Ohira and colleagues characterized 2'-phosphouridine (Up) at position 47 in tRNAs of thermophilic Sulfurisphaera tokodaii, as well as in several other archaea and thermophilic bacteria. From the presence of the gene arkI corresponding to the RNA kinase responsible for Up47 formation, they also concluded that Up47 should be present in tRNAs of other thermophilic Archaea Reanalysis of our earlier data confirms that the unidentified residue in tRNAs of both P. furiosus and S. acidocaldarius is indeed 2'-phosphouridine followed by m5C48. Moreover, we find this modification in several tRNAs of other Archaea and of the hyperthermophilic bacterium Aquifex aeolicus.

An old friend with a new face: tRNA-derived small RNAs with big regulatory potential in cancer biology

Transfer RNAs (tRNAs) are small non-coding RNAs (sncRNAs) essential for protein translation. Emerging evidence suggests that tRNAs can also be processed into smaller fragments, tRNA-derived small RNAs (tsRNAs), a novel class of sncRNAs with powerful applications and high biological relevance to cancer. tsRNAs biogenesis is heterogeneous and involves different ribonucleases, such as Angiogenin and Dicer. For many years, tsRNAs were thought to be just degradation products. However, accumulating evidence shows their roles in gene expression: either directly via destabilising the mRNA or the ribosomal machinery, or indirectly via regulating the expression of ribosomal components. Furthermore, tsRNAs participate in various biological processes linked to cancer, including apoptosis, cell cycle, immune response, and retroviral insertion into the human genome. It is emerging that tsRNAs have significant therapeutic potential. Endogenous tsRNAs can be used as cancer biomarkers, while synthetic tsRNAs and antisense oligonucleotides can be employed to regulate gene expression. In this review, we are recapitulating the regulatory roles of tsRNAs, with a focus on cancer biology.

Integrated Metabolomic and Transcriptomic Analysis of Modified Nucleosides for Biomarker Discovery in Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) accounts for ~75% of kidney cancers. The biallelic inactivation of the von Hippel-Lindau tumor suppressor gene (VHL) is the truncal driver mutation of most cases of ccRCC. Cancer cells are metabolically reprogrammed and excrete modified nucleosides in larger amounts due to their increased RNA turnover. Modified nucleosides occur in RNAs and cannot be recycled by salvage pathways. Their potential as biomarkers has been demonstrated for breast or pancreatic cancer. To assess their suitability as biomarkers in ccRCC, we used an established murine ccRCC model, harboring Vhl, Trp53 and Rb1 (VPR) knockouts. Cell culture media of this ccRCC model and primary murine proximal tubular epithelial cells (PECs) were investigated by HPLC coupled to triple-quadrupole mass spectrometry using multiple-reaction monitoring. VPR cell lines were significantly distinguishable from PEC cell lines and excreted higher amounts of modified nucleosides such as pseudouridine, 5-methylcytidine or 2'-O-methylcytidine. The method's reliability was confirmed in serum-starved VPR cells. RNA-sequencing revealed the upregulation of specific enzymes responsible for the formation of those modified nucleosides in the ccRCC model. These enzymes included Nsun2, Nsun5, Pus1, Pus7, Naf1 and Fbl. In this study, we identified potential biomarkers for ccRCC for validation in clinical trials.

Systematic comparison of tools used for m6A mapping from nanopore direct RNA sequencing

N6-methyladenosine (m6A) has been increasingly recognized as a new and important regulator of gene expression. To date, transcriptome-wide m6A detection primarily relies on well-established methods using next-generation sequencing (NGS) platform. However, direct RNA sequencing (DRS) using the Oxford Nanopore Technologies (ONT) platform has recently emerged as a promising alternative method to study m6A. While multiple computational tools are being developed to facilitate the direct detection of nucleotide modifications, little is known about the capabilities and limitations of these tools. Here, we systematically compare ten tools used for mapping m6A from ONT DRS data. We find that most tools present a trade-off between precision and recall, and integrating results from multiple tools greatly improve performance. Using a negative control could improve precision by subtracting certain intrinsic bias. We also observed variation in detection capabilities and quantitative information among motifs, and identified sequencing depth and m6A stoichiometry as potential factors affecting performance. Our study provides insight into the computational tools currently used for mapping m6A based on ONT DRS data and highlights the potential for further improving these tools, which may serve as the basis for future research.

Interplay of RNA 2'-O-methylations with viral replication

Viral RNAs (vRNAs) are decorated by post-transcriptional modifications, including methylation of nucleotides. Methylations regulate biological functions linked to the sequence, structure, and protein interactome of RNA. Several RNA viruses were found to harbor 2'-O-methylations, affecting the ribose moiety of RNA. This mark was initially shown to target the first and second nucleotides of the 5'-end cap structure of mRNA. More recently, nucleotides within vRNA were also reported to carry 2'-O-methylations. The consequences of such methylations are still puzzling since they were associated with both proviral and antiviral effects. Here, we focus on the mechanisms governing vRNA 2'-O-methylation and we explore the possible roles of this epitranscriptomic modification for viral replication.

Ribonucleosides from tRNA in hyperglycemic mammalian cells and diabetic murine cardiac models

AIMS

Cardiomyopathy is a diabetic comorbidity with few molecular targets. To address this, we evaluated transfer RNA (tRNA) modifications in the diabetic heart because tRNA modifications have been implicated in diabetic etiologies.

MAIN METHODS

tRNA was isolated from aorta, apex, and atrial tissue of healthy and diabetic murine hearts and related hyperglycemic cell models. tRNA modifications and canonical ribonucleosides were quantified by liquid-chromatography tandem mass spectrometry (LC-MS/MS) using stable isotope dilution. Correlations between ribonucleosides and diabetic comorbidity pathology were assessed using statistical analyses.

KEY FINDINGS

Total tRNA ribonucleoside levels were analyzed from cell types and healthy and diabetic murine heart tissue. Each heart structure had characteristic ribonucleoside profiles and quantities. Several ribonucleosides were observed as significantly different in hyperglycemic cells and diabetic tissues. In hyperglycemic models, ribonucleosides N4-acetylcytidine (ac4C), 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2U), 5-methylcytidine (m5C), and N1-methylguanosine (m1G) were anomalous. Specific tRNA modifications known to be on murine tRNAIni(CAU) were higher in diabetic heart tissue which suggests that tRNA modifications could be regulating translation in diabetes.

SIGNIFICANCE

We identified tRNA ribonucleosides and tRNA species associated with hyperglycemia and diabetic etiology.

RNA levers and switches controlling viral gene expression

RNA viruses are diverse and abundant pathogens that are responsible for numerous human diseases. RNA viruses possess relatively compact genomes and have therefore evolved multiple mechanisms to maximize their coding capacities, often by encoding overlapping reading frames. These reading frames are then decoded by mechanisms such as alternative splicing and ribosomal frameshifting to produce multiple distinct proteins. These solutions are enabled by the ability of the RNA genome to fold into 3D structures that can mimic cellular RNAs, hijack host proteins, and expose or occlude regulatory protein-binding motifs to ultimately control key process in the viral life cycle. We highlight recent findings focusing on less conventional mechanisms of gene expression and new discoveries on the role of RNA structures.

How does precursor RNA structure influence RNA processing and gene expression?

RNA is a fundamental biomolecule that has many purposes within cells. Due to its single-stranded and flexible nature, RNA naturally folds into complex and dynamic structures. Recent technological and computational advances have produced an explosion of RNA structural data. Many RNA structures have regulatory and functional properties. Studying the structure of nascent RNAs is particularly challenging due to their low abundance and long length, but their structures are important because they can influence RNA processing. Precursor RNA processing is a nexus of pathways that determines mature isoform composition and that controls gene expression. In this review, we examine what is known about human nascent RNA structure and the influence of RNA structure on processing of precursor RNAs. These known structures provide examples of how other nascent RNAs may be structured and show how novel RNA structures may influence RNA processing including splicing and polyadenylation. RNA structures can be targeted therapeutically to treat disease.

The importance of pseudouridylation: human disorders related to the fifth nucleoside

Pseudouridylation is one of the most abundant RNA modifications in eukaryotes, making pseudouridine known as the "fifth nucleoside." This highly conserved alteration affects all non-coding and coding RNA types. Its role and importance have been increasingly widely researched, especially considering that its absence or damage leads to serious hereditary diseases. Here, we summarize the human genetic disorders described to date that are related to the participants of the pseudouridylation process.

Rewired m6A epitranscriptomic networks link mutant p53 to neoplastic transformation

N6-methyladenosine (m6A), one of the most prevalent mRNA modifications in eukaryotes, plays a critical role in modulating both biological and pathological processes. However, it is unknown whether mutant p53 neomorphic oncogenic functions exploit dysregulation of m6A epitranscriptomic networks. Here, we investigate Li-Fraumeni syndrome (LFS)-associated neoplastic transformation driven by mutant p53 in iPSC-derived astrocytes, the cell-of-origin of gliomas. We find that mutant p53 but not wild-type (WT) p53 physically interacts with SVIL to recruit the H3K4me3 methyltransferase MLL1 to activate the expression of m6A reader YTHDF2, culminating in an oncogenic phenotype. Aberrant YTHDF2 upregulation markedly hampers expression of multiple m6A-marked tumor-suppressing transcripts, including CDKN2B and SPOCK2, and induces oncogenic reprogramming. Mutant p53 neoplastic behaviors are significantly impaired by genetic depletion of YTHDF2 or by pharmacological inhibition using MLL1 complex inhibitors. Our study reveals how mutant p53 hijacks epigenetic and epitranscriptomic machinery to initiate gliomagenesis and suggests potential treatment strategies for LFS gliomas.

Potential Roles of m6A and FTO in Synaptic Connectivity and Major Depressive Disorder

RNA modifications known as epitranscriptomics have emerged as a novel layer of transcriptomic regulation. Like the well-studied epigenetic modifications characterized in DNA and on histone-tails, they have been shown to regulate activity-dependent gene expression and play a vital role in shaping synaptic connections in response to external stimuli. Among the hundreds of known RNA modifications, N6-methyladenosine (m6A) is the most abundant mRNA modification in eukaryotes. Through recognition of its binding proteins, m6A can regulate various aspects of mRNA metabolism and is essential for maintaining higher brain functions. Indeed, m6A is highly enriched in synapses and is involved in neuronal plasticity, learning and memory, and adult neurogenesis. m6A can also respond to environmental stimuli, suggesting an important role in linking molecular and behavioral stress. This review summarizes key findings from fields related to major depressive disorder (MDD) including stress and learning and memory, which suggest that activity-dependent m6A changes may, directly and indirectly, contribute to synaptic connectivity changes underlying MDD. Furthermore, we will highlight the roles of m6A and FTO, a m6A eraser, in the context of depressive-like behaviors. Although we have only begun to explore m6A in the context of MDD and psychiatry, elucidating a link between m6A and MDD presents a novel molecular mechanism underlying MDD pathogenesis.

Chemical Space Virtual Screening against Hard-to-Drug RNA Methyltransferases DNMT2 and NSUN6

Targeting RNA methyltransferases with small molecules as inhibitors or tool compounds is an emerging field of interest in epitranscriptomics and medicinal chemistry. For two challenging RNA methyltransferases that introduce the 5-methylcytosine (m⁵C) modification in different tRNAs, namely DNMT2 and NSUN6, an ultra-large commercially available chemical space was virtually screened by physicochemical property filtering, molecular docking, and clustering to identify new ligands for those enzymes. Novel chemotypes binding to DNMT2 and NSUN6 with affinities down to K_D,app = 37 µM and K_D,app = 12 µM, respectively, were identified using a microscale thermophoresis (MST) binding assay. These compounds represent the first molecules with a distinct structure from the cofactor SAM and have the potential to be developed into activity-based probes for these enzymes. Additionally, the challenges and strategies of chemical space docking screens with special emphasis on library focusing and diversification are discussed.

Functional analysis of tRNA modification enzymes using mutational profiling

Transfer RNA (tRNA) delivers amino acids to the ribosome and functions as an essential adapter molecule for decoding codons on the messenger RNA (mRNA) during protein synthesis. Before attaining their proper activity, tRNAs undergo multiple post-transcriptional modifications with highly diversified roles such as stabilization of the tRNA structure, recognition of aminoacyl tRNA synthetases, precise codon-anticodon recognition, support of viral replication and onset of immune responses. The synthesis of the majority of modified nucleosides is catalyzed by a site-specific tRNA modification enzyme. This chapter provides a detailed protocol for using mutational profiling to analyze the enzymatic function of a tRNA methyltransferase in a high-throughput manner. In a previous study, we took tRNA m1A22 methyltransferase TrmK from Geobacillus stearothermophilus as a model tRNA methyltransferase and applied this protocol to gain mechanistic insights into how TrmK recognizes the substrate tRNAs. In theory, this protocol can be used unaltered for studying enzymes that catalyze modifications at the Watson-Crick face such as 1-methyladenosine (m1A), 3-methylcytosine (m3C), 3-methyluridine (m3U), 1-methylguanosine (m1G), and N2,N2-dimethylguanosine (m22G).

Epitranscriptomics in the development, functions, and disorders of cancer stem cells

Biomolecular modifications play an important role in the development of life, and previous studies have investigated the role of DNA and proteins. In the last decade, with the development of sequencing technology, the veil of epitranscriptomics has been gradually lifted. Transcriptomics focuses on RNA modifications that affect gene expression at the transcriptional level. With further research, scientists have found that changes in RNA modification proteins are closely linked to cancer tumorigenesis, progression, metastasis, and drug resistance. Cancer stem cells (CSCs) are considered powerful drivers of tumorigenesis and key factors for therapeutic resistance. In this article, we focus on describing RNA modifications associated with CSCs and summarize the associated research progress. The aim of this review is to identify new directions for cancer diagnosis and targeted therapy.

Construction of a Single-Molecule Biosensor for Antibody-Free Detection of Locus-Specific N6-Methyladenosine in Cancer Cells and Tissues

N⁶-Methyladenosine (m⁶A) has emerged as a key post-transcriptional regulator in mRNA metabolism, and its dysregulation is associated with multiple human diseases. Herein, we construct a single-molecule fluorescent biosensor for antibody-free detection of locus-specific m⁶A in cancer cells and tissues. A 5'-biotinylated capture probe and a 3'-hydroxylated assistant probe are designed for the recognition of specific m⁶A-mRNA. The m⁶A-sensitive endoribonuclease MazF can identify and cleave the unmethylated mRNA, and the retained intact m⁶A-mRNA can hybridize with assistant probes and capture probes to achieve sandwich hybrids. The sandwich hybrids are immobilized on magnetic beads (MBs) to initiate the terminal deoxynucleotidyl transferase (TdT)-assisted polymerization, facilitating the continuous incorporation of Cy5-dATP to form long Cy5-polyA tails for the production of an on-bead amplified fluorescence signal. After magnetic separation and exonuclease digestion, numerous Cy5 fluorophores are released and subsequently measured by single-molecule detection. Especially, this biosensor is implemented simply and isothermally without the involvement of either radiolabeling or m⁶A-specific antibody. Moreover, this biosensor shows ultrahigh sensitivity with a detection limit of 2.24 × 10^-17 M, and it can discriminate a 0.01% m⁶A level from a large pool of coexisting counterparts. Furthermore, this biosensor can be used for monitoring cellular m⁶A-mRNA expression and differentiating the m⁶A level in the breast cancer patient tissues from that in the healthy person tissues, providing a new avenue for clinical diagnosis and epitranscriptomic research.

Two-step conversion of uridine and cytidine to variously C5-C functionalized analogs

C5-substituted pyrimidine nucleosides are an important class of molecules that have practical use as biological probes and pharmaceuticals. Herein we report an operationally simple protocol for C5-functionalization of uridine and cytidine via transformation of underexploited 5-trifluoromethyluridine or 5-trifluoromethylcytidine, respectively. The unique reactivity of the CF₃ group in the aromatic ring allowed the direct incorporation of several distinct C5-C "carbon substituents": carboxyl, nitrile, ester, amide, and amidine.

M6A-BERT-Stacking: A Tissue-Specific Predictor for Identifying RNA N6-Methyladenosine Sites Based on BERT and Stacking Strategy

As the most abundant RNA methylation modification, N6-methyladenosine (m6A) could regulate asymmetric and symmetric division of hematopoietic stem cells and play an important role in various diseases. Therefore, the precise identification of m6A sites around the genomes of different species is a critical step to further revealing their biological functions and influence on these diseases. However, the traditional wet-lab experimental methods for identifying m6A sites are often laborious and expensive. In this study, we proposed an ensemble deep learning model called m6A-BERT-Stacking, a powerful predictor for the detection of m6A sites in various tissues of three species. First, we utilized two encoding methods, i.e., di ribonucleotide index of RNA (DiNUCindex_RNA) and k-mer word segmentation, to extract RNA sequence features. Second, two encoding matrices together with the original sequences were respectively input into three different deep learning models in parallel to train three sub-models, namely residual networks with convolutional block attention module (Resnet-CBAM), bidirectional long short-term memory with attention (BiLSTM-Attention), and pre-trained bidirectional encoder representations from transformers model for DNA-language (DNABERT). Finally, the outputs of all sub-models were ensembled based on the stacking strategy to obtain the final prediction of m6A sites through the fully connected layer. The experimental results demonstrated that m6A-BERT-Stacking outperformed most of the existing methods based on the same independent datasets.

N6-methyladenosine reader YTHDF family in biological processes: Structures, roles, and mechanisms

As the most abundant and conserved internal modification in eukaryote RNAs, N6-methyladenosine (m6A) is involved in a wide range of physiological and pathological processes. The YT521-B homology (YTH) domain-containing family proteins (YTHDFs), including YTHDF1, YTHDF2, and YTHDF3, are a class of cytoplasmic m6A-binding proteins defined by the vertebrate YTH domain, and exert extensive functions in regulating RNA destiny. Distinct expression patterns of the YTHDF family in specific cell types or developmental stages result in prominent differences in multiple biological processes, such as embryonic development, stem cell fate, fat metabolism, neuromodulation, cardiovascular effect, infection, immunity, and tumorigenesis. The YTHDF family mediates tumor proliferation, metastasis, metabolism, drug resistance, and immunity, and possesses the potential of predictive and therapeutic biomarkers. Here, we mainly summary the structures, roles, and mechanisms of the YTHDF family in physiological and pathological processes, especially in multiple cancers, as well as their current limitations and future considerations. This will provide novel angles for deciphering m6A regulation in a biological system.

A Split CRISPR/Cas13b System for Conditional RNA Regulation and Editing

The CRISPR/Cas13b system has been demonstrated as a robust tool for versatile RNA studies and relevant applications. New strategies enabling precise control of Cas13b/dCas13b activities and minimal interference with native RNA activities will further facilitate the understanding and regulation of RNA functions. Here, we engineered a split Cas13b system that can be conditionally activated and deactivated under the induction of abscisic acid (ABA), which achieved the downregulation of endogenous RNAs in dosage- and time-dependent manners. Furthermore, an ABA inducible split dCas13b system was generated to achieve temporally controlled deposition of m6A at specific sites on cellular RNAs through conditional assembly and disassembly of split dCas13b fusion proteins. We also showed that the activities of split Cas13b/dCas13b systems can be modulated by light via using a photoactivatable ABA derivative. Overall, these split Cas13b/dCas13b platforms expand the existing repertoire of the CRISPR and RNA regulation toolkit to achieve targeted manipulation of RNAs in native cellular environments with minimal functional disruption to these endogenous RNAs.

In-source fragmentation of nucleosides in electrospray ionization towards more sensitive and accurate nucleoside analysis

Nucleosides have been found to suffer in-source fragmentation (ISF) in electrospray ionization mass spectrometry, which leads to reduced sensitivity and ambiguous identification. In this work, a combination of theoretical calculations and nuclear magnetic resonance analysis revealed the key role of protonation at N3 near the glycosidic bond during ISF. Therefore, an ultrasensitive liquid chromatography-tandem mass spectrometry system for 5-formylcytosine detection was developed with 300 fold signal enhancement. Also, we established a MS1-only platform for nucleoside profiling and successfully identified sixteen nucleosides in the total RNA of MCF-7 cells. Taking ISF into account, we can realize analysis with higher sensitivity and less ambiguity, not only for nucleosides, but for other molecules with similar protonation and fragmentation behaviors.

Recent advances in the plant epitranscriptome

Chemical modifications of RNAs, known as the epitranscriptome, are emerging as widespread regulatory mechanisms underlying gene regulation. The field of epitranscriptomics advances recently due to improved transcriptome-wide sequencing strategies for mapping RNA modifications and intensive characterization of writers, erasers, and readers that deposit, remove, and recognize RNA modifications, respectively. Herein, we review recent advances in characterizing plant epitranscriptome and its regulatory mechanisms in post-transcriptional gene regulation and diverse physiological processes, with main emphasis on N⁶-methyladenosine (m⁶A) and 5-methylcytosine (m⁵C). We also discuss the potential and challenges for utilization of epitranscriptome editing in crop improvement.

A New Bacterial Adenosine-Derived Nucleoside as an Example of RNA Modification Damage

The fields of RNA modification and RNA damage both exhibit a plethora of non-canonical nucleoside structures. While RNA modifications have evolved to improve RNA function, the term RNA damage implies detrimental effects. Based on stable isotope labelling and mass spectrometry, we report the identification and characterisation of 2-methylthio-1,N6-ethenoadenosine (ms² ϵA), which is related to 1,N6-ethenoadenine, a lesion resulting from exposure of nucleic acids to alkylating chemicals in vivo. In contrast, a sophisticated isoprene labelling scheme revealed that ms² ϵA biogenesis involves cleavage of a prenyl moiety in the known transfer RNA (tRNA) modification 2-methylthio-N6-isopentenyladenosine (ms² i⁶ A). The relative abundance of ms² ϵA in tRNAs from translating ribosomes suggests reduced function in comparison to its parent RNA modification, establishing the nature of the new structure in a newly perceived overlap of the two previously separate fields, namely an RNA modification damage.

N6-methyladenosine of Spi2a attenuates inflammation and sepsis-associated myocardial dysfunction in mice

Bacteria-triggered sepsis is characterized by systemic, uncontrolled inflammation in affected individuals. Controlling the excessive production of pro-inflammatory cytokines and subsequent organ dysfunction in sepsis remains challenging. Here, we demonstrate that Spi2a upregulation in lipopolysaccharide (LPS)-stimulated bone marrow-derived macrophages reduces the production of pro-inflammatory cytokines and myocardial impairment. In addition, exposure to LPS upregulates the lysine acetyltransferase, KAT2B, to promote METTL14 protein stability through acetylation at K398, leading to the increased m⁶A methylation of Spi2a in macrophages. m⁶A-methylated Spi2a directly binds to IKKβ to impair IKK complex formation and inactivate the NF-κB pathway. The loss of m⁶A methylation in macrophages aggravates cytokine production and myocardial damage in mice under septic conditions, whereas forced expression of Spi2a reverses this phenotype. In septic patients, the mRNA expression levels of the human orthologue SERPINA3 negatively correlates with those of the cytokines, TNF, IL-6, IL-1β and IFNγ. Altogether, these findings suggest that m⁶A methylation of Spi2a negatively regulates macrophage activation in the context of sepsis.

Epitranscriptomics in metabolic disease

While epigenetic modifications of DNA and histones play main roles in gene transcription regulation, recently discovered post-transcriptional RNA modifications, known as epitranscriptomic modifications, have been found to have a profound impact on gene expression by regulating RNA stability, localization and decoding efficiency. Importantly, genetic variations or environmental perturbations of epitranscriptome modifiers (that is, writers, erasers and readers) are associated with obesity and metabolic diseases, such as type 2 diabetes. The epitranscriptome is closely coupled to epigenetic signalling, adding complexity to our understanding of gene expression in both health and disease. Moreover, the epitranscriptome in the parental generation can affect organismal phenotypes in the next generation. In this Review, we discuss the relationship between epitranscriptomic modifications and metabolic diseases, their relationship with the epigenome and possible therapeutic strategies.

The Role of m6A Modifications in B-Cell Development and B-Cell-Related Diseases

B cells are a class of professional antigen-presenting cells that produce antibodies to mediate humoral immune response and participate in immune regulation. m⁶A modification is the most common RNA modification in mRNA; it involves almost all aspects of RNA metabolism and can affect RNA splicing, translation, stability, etc. This review focuses on the B-cell maturation process as well as the role of three m⁶A modification-related regulators-writer, eraser, and reader-in B-cell development and B-cell-related diseases. The identification of genes and modifiers that contribute to immune deficiency may shed light on regulatory requirements for normal B-cell development and the underlying mechanism of some common diseases.

METTL3 regulates breast cancer-associated alternative splicing switches

Alternative splicing (AS) enables differential inclusion of exons from a given transcript, thereby contributing to the transcriptome and proteome diversity. Aberrant AS patterns play major roles in the development of different pathologies, including breast cancer. N⁶-methyladenosine (m⁶A), the most abundant internal modification of eukaryotic mRNA, influences tumor progression and metastasis of breast cancer, and it has been recently linked to AS regulation. Here, we identify a specific AS signature associated with breast tumorigenesis in vitro. We characterize for the first time the role of METTL3 in modulating breast cancer-associated AS programs, expanding the role of the m⁶A-methyltransferase in tumorigenesis. Specifically, we find that both m⁶A deposition in splice site boundaries and in splicing and transcription factor transcripts, such as MYC, direct AS switches of specific breast cancer-associated transcripts. Finally, we show that five of the AS events validated in vitro are associated with a poor overall survival rate for patients with breast cancer, suggesting the use of these AS events as a novel potential prognostic biomarker.

Identification of a novel 5-aminomethyl-2-thiouridine methyltransferase in tRNA modification

The uridine at the 34th position of tRNA, which is able to base pair with the 3'-end codon on mRNA, is usually modified to influence many aspects of decoding properties during translation. Derivatives of 5-methyluridine (xm5U), which include methylaminomethyl (mnm-) or carboxymethylaminomethyl (cmnm-) groups at C5 of uracil base, are widely conserved at the 34th position of many prokaryotic tRNAs. In Gram-negative bacteria such as Escherichia coli, a bifunctional MnmC is involved in the last two reactions of the biosynthesis of mnm5(s2)U, in which the enzyme first converts cmnm5(s2)U to 5-aminomethyl-(2-thio)uridine (nm5(s2)U) and subsequently installs the methyl group to complete the formation of mnm5(s2)U. Although mnm5s2U has been identified in tRNAs of Gram-positive bacteria and plants as well, their genomes do not contain an mnmC ortholog and the gene(s) responsible for this modification is unknown. We discovered that MnmM, previously known as YtqB, is the methyltransferase that converts nm5s2U to mnm5s2U in Bacillus subtilis through comparative genomics, gene complementation experiments, and in vitro assays. Furthermore, we determined X-ray crystal structures of MnmM complexed with anticodon stem loop of tRNAGln. The structures provide the molecular basis underlying the importance of U33-nm5s2U34-U35 as the key determinant for the specificity of MnmM.

Characterizing RNA modifications in the central nervous system and single cells by RNA sequencing and liquid chromatography-tandem mass spectrometry techniques

Post-transcriptional modifications to RNA constitute a newly appreciated layer of translation regulation in the central nervous system (CNS). The identity, stoichiometric quantity, and sequence position of these unusual epitranscriptomic marks are central to their function, making analytical methods that are capable of accurate and reproducible measurements paramount to the characterization of the neuro-epitranscriptome. RNA sequencing-based methods and liquid chromatography-tandem mass spectrometry (LC-MS/MS) techniques have been leveraged to provide an early glimpse of the landscape of RNA modifications in bulk CNS tissues. However, recent advances in sample preparation, separations, and detection methods have revealed that individual cells display remarkable heterogeneity in their RNA modification profiles, raising questions about the prevalence and function of cell-specific distributions of post-transcriptionally modified nucleosides in the brain. In this Trends article, we present an overview of RNA sequencing and LC-MS/MS methodologies for the analysis of RNA modifications in the CNS with special emphasis on recent advancements in techniques that facilitate single-cell and subcellular detection.

Identification of RNA helicases with unwinding activity on angiogenin-processed tRNAs

Stress-induced tRNA fragmentation upon environmental insult is a conserved cellular process catalysed by endonucleolytic activities targeting mature tRNAs. The resulting tRNA-derived small RNAs (tsRNAs) have been implicated in various biological processes that impact cell-to-cell signalling, cell survival as well as gene expression regulation during embryonic development. However, how endonuclease-targeted tRNAs give rise to individual and potentially biologically active tsRNAs remains poorly understood. Here, we report on the in vivo identification of proteins associated with stress-induced tsRNAs-containing protein complexes, which, together with a 'tracer tRNA' assay, were used to uncover enzymatic activities that can bind and process specific endonuclease-targeted tRNAs in vitro. Among those, we identified conserved ATP-dependent RNA helicases which can robustly separate tRNAs with endonuclease-mediated 'nicks' in their anticodon loops. These findings shed light on the existence of cellular pathways dedicated to producing individual tsRNAs after stress-induced tRNA hydrolysis, which adds to our understanding as to how tRNA fragmentation and the resulting tsRNAs might exert physiological impact.

Tet-dependent 5-hydroxymethyl-Cytosine modification of mRNA regulates the axon guidance genes robo2 and slit in Drosophila

Modifications of mRNA, especially methylation of adenosine, have recently drawn much attention. The much rarer modification, 5-hydroxymethylation of cytosine (5hmC), is not well understood and is the subject of this study. Vertebrate Tet proteins are 5-methylcytosine (5mC) hydroxylases enzymes catalyzing the transition of 5mC to 5hmC in DNA and have recently been shown to have the same function in messenger RNAs in both vertebrates and in Drosophila. The Tet gene is essential in Drosophila because Tet knock-out animals do not reach adulthood. We describe the identification of Tet-target genes in the embryo and larval brain by determining Tet DNA-binding sites throughout the genome and by mapping the Tet-dependent 5hmrC modifications transcriptome-wide. 5hmrC-modified sites can be found along the entire transcript and are preferentially located at the promoter where they overlap with histone H3K4me3 peaks. The identified mRNAs are frequently involved in neuron and axon development and Tet knock-out led to a reduction of 5hmrC marks on specific mRNAs. Among the Tet-target genes were the robo2 receptor and its slit ligand that function in axon guidance in Drosophila and in vertebrates. Tet knock-out embryos show overlapping phenotypes with robo2 and are sensitized to reduced levels of slit. Both Robo2 and Slit protein levels were markedly reduced in Tet KO larval brains. Our results establish a role for Tet-dependent 5hmrC in facilitating the translation of modified mRNAs, primarily in developing nerve cells.

RNA Modification Detection Using Nanopore Direct RNA Sequencing and nanoDoc2

RNA modifications regulate multiple aspects of cellular function including RNA splicing, translation, export, decay, stability, and phase separation. One of the comprehensive ways to detect such modifications is by the recent advancement of direct RNA sequencing from Oxford Nanopore Technologies (ONT). However, this method obtains a large amount of data with high complexity in the form of raw current signal that poses a new informatics challenge to accurately detect those modifications. Here, we provide nanoDoc2, a software to detect multiple types of RNA modification from nanopore direct RNA sequencing data. The nanoDoc2 includes a novel signal segmentation algorithm based on the trace value-a base probability feature that is added by the Guppy basecalling program from ONT during processing of the raw signal. The core of nanoDoc2 includes a machine learning algorithm in which a 6-mer segmented raw current signal is analyzed by deep one-class classification using a WaveNet-based neural network. As an output, an RNA modification is detected by a statistical score in each candidate position. Herein, we describe the detailed instructions on how to use nanoDoc2 for signal segmentation, train/test the neural network, and finally predict RNA modifications present in nanopore direct RNA sequencing data.

Purine catabolism by enterobacteria

Purines are abundant among organic nitrogen sources and have high nitrogen content. Accordingly, microorganisms have evolved different pathways to catabolize purines and their metabolic products such as allantoin. Enterobacteria from the genera Escherichia, Klebsiella and Salmonella have three such pathways. First, the HPX pathway, found in the genus Klebsiella and very close relatives, catabolizes purines during aerobic growth, extracting all four nitrogen atoms in the process. This pathway includes several known or predicted enzymes not previously observed in other purine catabolic pathways. Second, the ALL pathway, found in strains from all three species, catabolizes allantoin during anaerobic growth in a branched pathway that also includes glyoxylate assimilation. This allantoin fermentation pathway originally was characterized in a gram-positive bacterium, and therefore is widespread. Third, the XDH pathway, found in strains from Escherichia and Klebsiella spp., at present is ill-defined but likely includes enzymes to catabolize purines during anaerobic growth. Critically, this pathway may include an enzyme system for anaerobic urate catabolism, a phenomenon not previously described. Documenting such a pathway would overturn the long-held assumption that urate catabolism requires oxygen. Overall, this broad capability for purine catabolism during either aerobic or anaerobic growth suggests that purines and their metabolites contribute to enterobacterial fitness in a variety of environments.