MODOMICS: An Operational Guide to the Use of the RNA Modification Pathways Database

Post-transcriptional Modifications of the Large Ribosome Subunit Assembly Intermediates in E. coli Expressing Helicase-Inactive DbpA Variant

RNA post-transcriptional modifications are ubiquitous across all organisms and serve as fundamental regulators of cellular homeostasis, growth, and stress adaptation. Techniques for the simultaneous detection of multiple RNA modifications in a high-throughput, single-nucleotide-resolution manner are largely absent in the field, and developing such techniques is of paramount importance. We used the Escherichia coli ribosome as a model system to develop novel techniques for RNA post-transcriptional modification detection, leveraging its extensive and diverse array of modifications. For modification detection, we performed reverse transcriptase reactions in the presence of Mn2⁺ and quantified the reverse transcriptase deletions and misincorporations at modification positions using Illumina next-generation sequencing. We simultaneously detected the following modifications in ribosomal RNA (rRNA): 1-methylguanosine (m1G), 2-methylguanosine (m2G), 3-methylpseudouridine, N6,N6-dimethyladenosine, and 3-methyluridine, without chemical treatment. Furthermore, subjecting the rRNA samples to 1-cyclohexyl-3-(2-morpholinoethyl) carbodiimide metho-p-toluenesulfonate followed by alkaline conditions allowed us to simultaneously detect pseudouridine, 7-methylguanosine (m7G), 5-hydroxycytidine (OH5C), 2-methyladenosine, and dihydrouridine (D). Finally, subjecting the rRNA samples to KMnO4 followed by alkaline conditions allowed us to simultaneously detect m7G, OH5C, and D. Our results reveal that m1G, m2G, m7G, and D are incorporated prior to the accumulation of the 27S, 35S, and 45S large subunit intermediates in cells expressing the helicase-inactive R331A DbpA construct. These intermediates belong to three distinct stages and pathways of large subunit ribosome assembly. Therefore, our results identify the time points in three pathways at which m1G, m2G, m7G, and D are incorporated into the large ribosome subunit and provide a framework for broader studies on RNA modification dynamics.

The Expanding Universe of Extensions and Tails: Ribosomal Proteins and Histones in RNA and DNA Complex Signaling and Dynamics

This short review bridges two biological fields: ribosomes and nucleosomes-two nucleoprotein assemblies that, along with many viruses, share proteins featuring long filamentous segments at their N- or C-termini. A central hypothesis is that these extensions and tails perform analogous functions in both systems. The evolution of these structures appears closely tied to the emergence of regulatory networks and signaling pathways, facilitating increasingly complex roles for ribosomes and nucleosome alike. This review begins by summarizing the structures and functions of ribosomes and nucleosomes, followed by a detailed comparison highlighting their similarities and differences, particularly in light of recent findings on the roles of ribosomal proteins in signaling and ribosome dynamics. The analysis seeks to uncover whether these systems operate based on shared principles and mechanisms. The nucleosome-ribosome analogy may offer valuable insights into unresolved questions in both fields. For instance, new structural insights from ribosomes might shed light on potential motifs formed by histone tails. From an evolutionary perspective, this study revisits the origins of signaling and regulation in ancient nucleoprotein assemblies, suggesting that tails and extensions may represent remnants of the earliest network systems governing signaling and dynamic control.

PCBP1/2 and TDP43 Function as NAT10 Adaptors to Mediate mRNA ac4C Formation in Mammalian Cells

Massive numbers of modified bases in mRNAs sculpt the epitranscriptome and play vital roles in RNA metabolism. The only known acetylated RNA modification, N-4-acetylcytidine (ac4C), is highly conserved across cell types and among species. Although the GCN5-related acetyltransferase 10 (NAT10) functions as an ac4C writer, the mechanism underlying the acetylation process is largely unknown. In this study, the NAT10/PCBP/TDP43 complex mediated mRNA ac4C formation in mammalian cells is identified. RNA-binding proteins (RBPs) are identified, affiliated with two different families, poly(rC)-binding protein 1/2 (PCBP1/2) and TAR DNA binding protein 43 (TDP43), as NAT10 adaptors for mRNA tethering and substrate selection. Knockdown of the adaptors resulted in decreased mRNA acetylation abundance in HEK293T cells and ablated cytidine-rich ac4C motifs. The adaptors also affect the ac4C sites by recruiting NAT10 to their binding sequences. The presence of the NAT10/PCBP/TDP43 complex in mouse testes highlights its potential physiological functions in vivo. These findings reveal the composition of the mRNA ac4C writer complex in mammalian cells and expand the knowledge of mRNA acetylation and ac4C site preferences.

Domain-knowledge enabled ensemble learning of 5-formylcytosine (f5C) modification sites

5-formylcytidine (f5C) is a unique post-transcriptional RNA modification found in mRNA and tRNA at the wobble site, playing a crucial role in mitochondrial protein synthesis and potentially contributing to the regulation of translation. Recent studies have unveiled that the f5C modifications may drive mitochondrial mRNA translation to power cancer metastasis. Accurate identification of f5C sites is essential for further unraveling their molecular functions and regulatory mechanisms, but there are currently no computational methods available for predicting their locations. In this study, we introduce an innovative ensemble approach, successfully enabling the computational recognition of Saccharomyces cerevisiae f5C. We conducted a comprehensive model selection process that involved multiple basic machine learning and deep learning algorithms such as recurrent neural networks, convolutional neural networks and Transformer-based models. Initially trained only on sequence information, these individual models achieved an AUROC ranging from 0.7104 to 0.7492. Through the integration of 32 novel domain-derived genomic features, the performance of individual models has significantly improved to an AUROC between 0.7309 and 0.8076. To further enhance accuracy and robustness, we then constructed the ensembles of these individual models with different combinations. The best performance attained by our ensemble models reached an AUROC of 0.8391. Shapley additive explanations were conducted to explain the significant contributions of genomic features, providing insights into the putative distribution of f5C across various topological regions and potentially paving the way for revealing their functional relevance within distinct genomic contexts. A freely accessible web server that allows real-time analysis of user-uploaded sites can be accessed at: www.rnamd.org/Resf5C-Pred.

Multicellular, IVT-derived, unmodified human transcriptome for nanopore-direct RNA analysis

Nanopore direct RNA sequencing (DRS) enables measurements of RNA modifications. Modification-free transcripts are a practical and targeted control for DRS, providing a baseline measurement for canonical nucleotides within a matched and biologically-derived sequence context. However, these controls can be challenging to generate and carry nanopore-specific nuances that can impact analyses. We produced DRS datasets using modification-free transcripts from in vitro transcription of cDNA from six immortalized human cell lines. We characterized variation across cell lines and demonstrated how these may be interpreted. These data will serve as a versatile control and resource to the community for RNA modification analyses of human transcripts.

RNA m6a Methylation Regulator Expression in Castration-Resistant Prostate Cancer Progression and Its Genetic Associations

N6-methyladenosine (m6A) methylation, a prevalent epitranscriptomic modification, plays a crucial role in regulating mRNA expression, stability, and translation in mammals. M6A regulators have gained attention for their potential implications in tumorigenesis and clinical applications, such as cancer diagnosis and therapeutics. The existing literature predominantly addresses m6A regulators in the context of primary prostate cancer (PCa). However, a notable gap in the knowledge emerges regarding the dynamic expression patterns of these regulators as PCa progresses towards the castration-resistant stage (CRPC). Employing sequential window acquisition of all theoretical mass spectra (SWATH-MS) and RNAseq analysis, we comprehensively profiled the expression of 27 m6A regulators in hormone/androgen-dependent and -independent PCa cell lines, revealing distinct clustering between tumor and adjacent normal prostate tissues. High-grade PCa tumors demonstrated the upregulation of METTL3, RBM15B, and HNRNAPA2B1 and the downregulation of ZC3H13, NUDT21, and FTO. Notably, we identified six m6A regulators associated with PCa survival. Additionally, association analysis of the PCa-associated risk loci in the cancer genome atlas program (TCGA) data unveiled genetic variations near the WTAP, HNRNPA2B1, and FTO genes as significant expression quantitative trait loci. In summary, our study unraveled abnormalities in m6A regulator expression in PCa progression, elucidating their association with PCa risk loci. Considering the heterogeneity within the PCa phenotypes and treatment responses, our findings suggest that prognostic stratification based on m6A regulator expression could enhance PCa diagnosis and prognosis.

N-acetyltransferase NAT10 controls cell fates via connecting mRNA cytidine acetylation to chromatin signaling

Cell fate transition involves dynamic changes of gene regulatory network and chromatin landscape, requiring multiple levels of regulation, yet the cross-talk between epitranscriptomic modification and chromatin signaling remains largely unknown. Here, we uncover that suppression of N-acetyltransferase 10 (NAT10), the writer for mRNA N4-acetylcytidine (ac4C) modification, can notably affect human embryonic stem cell (hESC) lineage differentiation and pluripotent reprogramming. With integrative analysis, we identify that NAT10-mediated ac4C modification regulates the protein levels of a subset of its targets, which are strongly enriched for fate-instructive chromatin regulators, and among them, histone chaperone ANP32B is experimentally verified and functionally relevant. Furthermore, NAT10-ac4C-ANP32B axis can modulate the chromatin landscape of their downstream genes (e.g., key regulators of Wnt and TGFβ pathways). Collectively, we show that NAT10 is an essential regulator of cellular plasticity, and its catalyzed mRNA cytidine acetylation represents a critical layer of epitranscriptomic modulation and uncover a previously unrecognized, direct cross-talk between epitranscriptomic modification and chromatin signaling during cell fate transitions.

Quantitative analysis of tRNA abundance and modifications by nanopore RNA sequencing

Transfer RNAs (tRNAs) play a central role in protein translation. Studying them has been difficult in part because a simple method to simultaneously quantify their abundance and chemical modifications is lacking. Here we introduce Nano-tRNAseq, a nanopore-based approach to sequence native tRNA populations that provides quantitative estimates of both tRNA abundances and modification dynamics in a single experiment. We show that default nanopore sequencing settings discard the vast majority of tRNA reads, leading to poor sequencing yields and biased representations of tRNA abundances based on their transcript length. Re-processing of raw nanopore current intensity signals leads to a 12-fold increase in the number of recovered tRNA reads and enables recapitulation of accurate tRNA abundances. We then apply Nano-tRNAseq to Saccharomyces cerevisiae tRNA populations, revealing crosstalks and interdependencies between different tRNA modification types within the same molecule and changes in tRNA populations in response to oxidative stress.

ECT12, an YTH-domain protein, is a potential mRNA m6A reader that affects abiotic stress responses by modulating mRNA stability in Arabidopsis

N6-methyladenosine (m6A), the most abundant modification found in eukaryotic mRNAs, is interpreted by m6A "readers," thus playing a crucial role in regulating RNA metabolism. The YT521-B homology-domain (YTHD) proteins, also known as EVOLUTIONARILY CONSERVED C-TERMINAL REGION (ECT), are recognized as m6A reader proteins in plants and animals. Among the 13 potential YTHD family proteins in Arabidopsis thaliana, the functions of only a few members are known. In this study, we determined the function of ECT12 (YTH11) as a potential m6A reader that plays a crucial role in response to abiotic stresses. The loss-of-function ect12 mutants showed no noticeable developmental defects under normal conditions but displayed hypersensitivity to salt or dehydration stress. The salt- or dehydration-hypersensitive phenotypes were correlated with altered levels of several m6A-modified stress-responsive transcripts. Notably, the increased or decreased transcript levels were associated with each transcript's reduced or enhanced decay, respectively. Electrophoretic mobility shift and RNA-immunoprecipitation assays showed that ECT12 binds to m6A-modified RNAs both in vitro and in planta, suggesting its role as an m6A reader. Collectively, these results indicate that the potential m6A reader ECT12 regulates the stability of m6A-modified RNA transcripts, thereby facilitating the response of Arabidopsis to abiotic stresses.

Stacking-ac4C: an ensemble model using mixed features for identifying n4-acetylcytidine in mRNA

N4-acetylcytidine (ac4C) is a modification of cytidine at the nitrogen-4 position, playing a significant role in the translation process of mRNA. However, the precise mechanism and details of how ac4C modifies translated mRNA remain unclear. Since identifying ac4C sites using conventional experimental methods is both labor-intensive and time-consuming, there is an urgent need for a method that can promptly recognize ac4C sites. In this paper, we propose a comprehensive ensemble learning model, the Stacking-based heterogeneous integrated ac4C model, engineered explicitly to identify ac4C sites. This innovative model integrates three distinct feature extraction methodologies: Kmer, electron-ion interaction pseudo-potential values (PseEIIP), and pseudo-K-tuple nucleotide composition (PseKNC). The model also incorporates the robust Cluster Centroids algorithm to enhance its performance in dealing with imbalanced data and alleviate underfitting issues. Our independent testing experiments indicate that our proposed model improves the Mcc by 15.61% and the ROC by 5.97% compared to existing models. To test our model's adaptability, we also utilized a balanced dataset assembled by the authors of iRNA-ac4C. Our model showed an increase in Sn of 4.1%, an increase in Acc of nearly 1%, and ROC improvement of 0.35% on this balanced dataset. The code for our model is freely accessible at https://github.com/louliliang/ST-ac4C.git, allowing users to quickly build their model without dealing with complicated mathematical equations.

RNA Post-transcriptional Modifications of an Early-Stage Large-Subunit Ribosomal Intermediate

Protein production by ribosomes is fundamental to life, and proper assembly of the ribosome is required for protein production. The RNA, which is post-transcriptionally modified, provides the platform for ribosome assembly. Thus, a complete understanding of ribosome assembly requires the determination of the RNA post-transcriptional modifications in all of the ribosome assembly intermediates and on each pathway. There are 26 RNA post-transcriptional modifications in 23S RNA of the mature Escherichia coli (E. coli) large ribosomal subunit. The levels of these modifications have been investigated extensively only for a small number of large subunit intermediates and under a limited number of cellular and environmental conditions. In this study, we determined the level of incorporations of 2-methyl adenosine, 3-methyl pseudouridine, 5-hydroxycytosine, and seven pseudouridines in an early-stage E. coli large-subunit assembly intermediate with a sedimentation coefficient of 27S. The 27S intermediate is one of three large subunit intermediates accumulated in E. coli cells lacking the DEAD-box RNA helicase DbpA and expressing the helicase inactive R331A DbpA construct. The majority of the investigated modifications are incorporated into the 27S large subunit intermediate to similar levels to those in the mature 50S large subunit, indicating that these early modifications or the enzymes that incorporate them play important roles in the initial events of large subunit ribosome assembly.

Establishment of a prognostic model toward lung squamous cell carcinoma based on m7G-related genes in the cancer genome atlas

Lung squamous cell carcinoma (LUSC) is a non-small cell lung cancer with a poor prognosis owing to late diagnosis. New molecular markers are urgently needed to improve the diagnosis and prognosis of LUSC. 7-Methylguanosine (m7G) modifications, a tRNA modification, are common in eubacteria, eukaryotes, and a few archaea. These modifications promote the turnover and stability of some mRNAs to prevent mRNA decay, improve translation efficiency, and reduce ribosomal pausing but are associated with poor survival in human cancer cells. However, expression of m7G-related genes in LUSC and their association with prognosis remain unclear. In the present study, we identified nine differentially expressed genes related to prognosis by comparing the expression profiles of tumor tissues (502 LUSC reports) with normal tissues (49 adjacent nontumor lung tissue reports). The genes included six upregulated genes (KLK7, LCE3E, AREG, KLK6, ZBED2, and MAPK4) and three downregulated genes (ADH1C, NTS, and ERLIN2). Based on these nine genes, patients with LUSC were classified into low- and high-risk groups to analyze the trends in prognosis. We found that the nine m7G-related genes play important roles in immune regulation, hormone regulation, and drug sensitivity through pathways including antigen processing and presentation, adherent plaques, extracellular matrix receptor interactions, drug metabolism of cytochrome P-450, and metabolism of cytochrome P-450 to xenobiotics; the functions of these genes are likely accomplished in part by m6A modifications. The effect of m7G-related genes on the diagnosis and prognosis of LUSC was further indicated by population analysis.NEW & NOTEWORTHY Based on the differential expression of 7-methylguanosine (m7G) modification-associated genes between normal and lung squamous cell carcinoma (LUSC) tissues, and considering the performance of our m7G-related gene risk profiles as independent risk factors in predicting overall survival, we conclude that m7G modification is closely linked to the development of LUSC. In addition, this study offers a new genetic marker for predicting the prognosis of patients with LUSC and presents a crucial theoretical foundation for future investigations on the relationship between m7G modification-related genes, immunity, and drug sensitivity in LUSC.

Transcriptome-wide m6A methylation in natural yellow leaf of Catalpa fargesii

N6-methyladenosine (m6A) is the most abundant internal modification in eukaryotic messenger RNA, and involved in various biological processes in plants. However, the distribution features and functions of mRNA m6A methylation have been poorly explored in woody perennial plants. In this study, a new natural variety with yellow-green leaves, named Maiyuanjinqiu, was screened from the seedlings of Catalpa fargesii. Based on the preliminary experiment, the m6A methylation levels in the leaves of Maiyuanjinqiu were significantly higher than those in C. fargesii. Furthermore, a parallel analysis of m6A-seq and RNA-seq was carried out in different leaf color sectors. The result showed that m6A modification were mostly identified around the 3'-untranslated regions (3'-UTR), which was slightly negatively correlated with the mRNA abundance. KEGG and GO analyses showed that m6A methylation genes were associated with photosynthesis, pigments biosynthesis and metabolism, oxidation-reduction and response to stress, etc. The overall increase of m6A methylation levels in yellow-green leaves might be associated with the decreased the expression of RNA demethylase gene CfALKBH5. The silencing of CfALKBH5 caused a chlorotic phenotype and increased m6A methylation level, which further confirmed our hypothesis. Our results suggested that mRNA m6A methylation could be considered as a vital epigenomic mark and contribute to the natural variations in plants.

FLK is an mRNA m6A reader that regulates floral transition by modulating the stability and splicing of FLC in Arabidopsis

N⁶-methyladenosine (m⁶A), which is added, removed, and interpreted by m⁶A writers, erasers, and readers, respectively, is the most abundant modification in eukaryotic mRNAs. The m⁶A marks play a pivotal role in the regulation of floral transition in plants. FLOWERING LOCUS K (FLK), an RNA-binding protein harboring K-homology (KH) motifs, is known to regulate floral transition by repressing the levels of a key floral repressor FLOWERING LOCUS C (FLC) in Arabidopsis. However, the molecular mechanism underlying FLK-mediated FLC regulation remains unclear. In this study, we identified FLK as a novel mRNA m⁶A reader protein that directly binds the m⁶A site in the 3'-untranslated region of FLC transcripts to repressing FLC levels by reducing its stability and splicing. Importantly, FLK binding of FLC transcripts was abolished in vir-1, an m⁶A writer mutant, and the late-flowering phenotype of the flk mutant could not be rescued by genetic complementation using the mutant FLK^m gene, in which the m⁶A reader encoding function was eliminated, indicating that FLK binds and regulates FLC expression in an m⁶A-dependent manner. Collectively, our study has addressed a long-standing question of how FLK regulates FLC transcript levels and established a molecular link between the FLK-mediated recognition of m⁶A modifications on FLC transcripts and floral transition in Arabidopsis.

Eliminating METTL1-mediated accumulation of PMN-MDSCs prevents hepatocellular carcinoma recurrence after radiofrequency ablation

BACKGROUND AND AIMS

Radiofrequency ablation (RFA) is an important curative therapy in hepatocellular carcinoma (HCC), but recurrence rate remains as high as all the other HCC therapeutic modalities. Methyltransferase 1 (METTL1), an enzyme for m 7 G tRNA modification, was reported to promote HCC development. Here, we assessed the role of METTL1 in shaping the immunosuppressive tumor microenvironment after insufficient RFA (iRFA).

APPROACH AND RESULTS

By immunohistochemistry and multiplex immunofluorescence (mIF) staining, we showed that METTL1 expression was enhanced in post-RFA recurrent HCC, accompanied by increased CD11b + CD15 + polymorphonuclear-myeloid-derived suppressor cells (PMN-MDSCs) and decreased CD8 + T cells. Mechanistically, heat-mediated METTL1 upregulation enhanced TGF-β2 translation to form the immunosuppressive environment by induction of myeloid-derived suppressor cell. Liver-specific overexpression or knockdown of Mettl1 significantly affected the accumulation of PMN-MDSCs and subsequently affected CD8 + T cell infiltration. Complete RFA successfully eliminated the tumor, whereas iRFA-treated mice exhibited enhanced tumor growth and metastasis with increased PMN-MDSC accumulation and decreased CD8 + T cells compared to sham surgery. Interrupting METTL1-TGF-β2-PMN-MDSC axis by anti-Ly6G antibody, or knockdown of hepatoma-intrinsic Mettl1 or Tgfb2 , or TGF-β signaling blockade significantly mitigated tumor progression induced by iRFA and restored CD8 + T cell population.

CONCLUSIONS

Our study sheds light on the pivotal role of METTL1 in modulating an immunosuppressive microenvironment and demonstrated that interrupting METTL1-TGF-β2-PMN-MDSC axis could be a therapeutic strategy to restore antitumor immunity and prevent HCC recurrence after RFA treatment, meriting further clinical studies.

Roles of oncogenes in esophageal squamous cell carcinoma and their therapeutic potentials

Esophageal squamous cell carcinoma (ESCC) is the most common type of esophageal cancer (EC) in Asia. It is a malignant digestive tract tumor with abundant gene mutations. Due to the lack of specific diagnostic markers and early cancer screening markers, most patients are diagnosed at an advanced stage. Genetic and epigenetic changes are closely related to the occurrence and development of ESCC. Here, We review the activation of proto-oncogenes into oncogenes through gene mutation and gene amplification in ESCC from a genetic and epigenetic genome perspective, We also discuss the specific regulatory mechanisms through which these oncogenes mainly affect the biological function and occurrence and development of ESCC through specific regulatory mechanisms. In addition, we summarize the clinical application value of these oncogenes is summarized, and it provides a feasible direction for clinical use as potential therapeutic and diagnostic markers.

m6A readers, writers, erasers, and the m6A epitranscriptome in breast cancer

Epitranscriptomic modification of RNA regulates human development, health, and disease. The true diversity of the transcriptome in breast cancer including chemical modification of transcribed RNA (epitranscriptomics) is not well understood due to limitations of technology and bioinformatic analysis. N-6-methyladenosine (m6A) is the most abundant epitranscriptomic modification of mRNA and regulates splicing, stability, translation, and intracellular localization of transcripts depending on m6A association with reader RNA-binding proteins. m6A methylation is catalyzed by the METTL3 complex and removed by specific m6A demethylase ALKBH5, with the role of FTO as an 'eraser' uncertain. In this review, we provide an overview of epitranscriptomics related to mRNA and focus on m6A in mRNA and its detection. We summarize current knowledge on altered levels of writers, readers, and erasers of m6A and their roles in breast cancer and their association with prognosis. We summarize studies identifying m6A peaks and sites in genes in breast cancer cells.

Nucleic acid adductomics - The next generation of adductomics towards assessing environmental health risks

This Discussion article aims to explore the potential for a new generation of assay to emerge from cellular and urinary DNA adductomics which brings together DNA-RNA- and, to some extent, protein adductomics, to better understand the role of the exposome in environmental health. Components of the exposome have been linked to an increased risk of various, major diseases, and to identify the precise nature, and size, of risk, in this complex mixture of exposures, powerful tools are needed. Modification of nucleic acids (NA) is a key consequence of environmental exposures, and a goal of cellular DNA adductomics is to evaluate the totality of DNA modifications in the genome, on the basis that this will be most informative. Consequently, an approach which encompasses modifications of all nucleic acids (NA) would be potentially yet more informative. This article focuses on NA adductomics, which brings together the assessment of both DNA and RNA modifications, including modified (2'-deoxy)ribonucleosides (2'-dN/rN), modified nucleobases (nB), plus: DNA-DNA, RNA-RNA, DNA-RNA, DNA-protein, and RNA-protein crosslinks (DDCL, RRCL, DRCL, DPCL, and RPCL, respectively). We discuss the need for NA adductomics, plus the pros and cons of cellular vs. urinary NA adductomics, and present some evidence for the feasibility of this approach. We propose that NA adductomics provides a more comprehensive approach to the study of nucleic acid modifications, which will facilitate a range of advances, including the identification of novel, unexpected modifications e.g., RNA-RNA, and DNA-RNA crosslinks; key modifications associated with mutagenesis; agent-specific mechanisms; and adductome signatures of key environmental agents, leading to the dissection of the exposome, and its role in human health/disease, across the life course.

RNA Modifications in Cancer Metabolism and Tumor Microenvironment

RNA modifications have recently been recognized as essential posttranscriptional regulators of gene expression in eukaryotes. Investigations over the past decade have revealed that RNA chemical modifications have profound effects on tumor initiation, progression, refractory, and recurrence. Tumor cells are notorious for their robust plasticity in response to the stressful microenvironment and undergo metabolic adaptations to sustain rapid cell proliferation, which is termed as metabolic reprogramming. Meanwhile, cancer-associated metabolic reprogramming leads to substantial alterations of intracellular and extracellular metabolites, which further reshapes the tumor microenvironment (TME). Moreover, cancer cells compete with tumor-infiltrating immune cells for the limited nutrients to maintain their proliferation and function in the TME. In this chapter, we review recent interesting findings on the engagement of epitranscriptomic pathways, especially the ones associated with N6-methyladenosine (m6A), in the regulation of cancer metabolism and the surrounding microenvironment. We also discuss the promising therapeutic approaches targeting RNA modifications for anti-tumor therapy.

Epigenetic and epitranscriptomic mechanisms of chromium carcinogenesis

Hexavalent chromium [Cr(VI)], a Group I carcinogen classified by the International Agency for Research on Cancer (IARC), represents one of the most common occupational and environmental pollutants. The findings from human epidemiological and laboratory animal studies show that long-term exposure to Cr(VI) causes lung cancer and other cancer. Although Cr(VI) is a well-recognized carcinogen, the mechanism of Cr(VI) carcinogenesis has not been well understood. Due to the fact that Cr(VI) undergoes a series of metabolic reductions once entering cells to generate reactive Cr metabolites and reactive oxygen species (ROS) causing genotoxicity, Cr(VI) is generally considered as a genotoxic carcinogen. However, more and more studies have demonstrated that acute or chronic Cr(VI) exposure also causes epigenetic dysregulations including changing DNA methylation, histone posttranslational modifications and regulatory non-coding RNA (microRNA and long non-coding RNA) expressions. Moreover, emerging evidence shows that Cr(VI) exposure is also capable of altering cellular epitranscriptome. Given the increasingly recognized importance of epigenetic and epitranscriptomic dysregulations in cancer initiation and progression, it is believed that Cr(VI) exposure-caused epigenetic and epitranscriptomic changes could play important roles in Cr(VI) carcinogenesis. The goal of this chapter is to review the epigenetic and epitranscriptomic effects of Cr(VI) exposure and discuss their roles in Cr(VI) carcinogenesis. Better understanding the mechanism of Cr(VI) carcinogenesis may identify new molecular targets for more efficient prevention and treatment of cancer resulting from Cr(VI) exposure.

ADAR2 Protein Is Associated with Overall Survival in GBM Patients and Its Decrease Triggers the Anchorage-Independent Cell Growth Signature

Background: Epitranscriptomic mechanisms, such as A-to-I RNA editing mediated by ADAR deaminases, contribute to cancer heterogeneity and patients’ stratification. ADAR enzymes can change the sequence, structure, and expression of several RNAs, affecting cancer cell behavior. In glioblastoma, an overall decrease in ADAR2 RNA level/activity has been reported. However, no data on ADAR2 protein levels in GBM patient tissues are available; and most data are based on ADARs overexpression experiments. Methods: We performed IHC analysis on GBM tissues and correlated ADAR2 levels and patients’ overall survival. We silenced ADAR2 in GBM cells, studied cell behavior, and performed a gene expression/editing analysis. Results: GBM tissues do not all show a low/no ADAR2 level, as expected by previous studies. Although, different amounts of ADAR2 protein were observed in different patients, with a low level correlating with a poor patient outcome. Indeed, reducing the endogenous ADAR2 protein in GBM cells promotes cell proliferation and migration and changes the cell’s program to an anchorage-independent growth mode. In addition, deep-seq data and bioinformatics analysis indicated multiple RNAs are differently expressed/edited upon siADAR2. Conclusion: ADAR2 protein is an important deaminase in GBM and its amount correlates with patient prognosis.

Dihydrouridine in the Transcriptome: New Life for This Ancient RNA Chemical Modification

Until recently, post-transcriptional modifications of RNA were largely restricted to noncoding RNA species. However, this belief seems to have quickly dissipated with the growing number of new modifications found in mRNA that were originally thought to be primarily tRNA-specific, such as dihydrouridine. Recently, transcriptomic profiling, metabolic labeling, and proteomics have identified unexpected dihydrouridylation of mRNAs, greatly expanding the catalog of novel mRNA modifications. These data also implicated dihydrouridylation in meiotic chromosome segregation, protein translation rates, and cell proliferation. Dihydrouridylation of tRNAs and mRNAs are introduced by flavin-dependent dihydrouridine synthases. In this review, we will briefly outline the current knowledge on the distribution of dihydrouridines in the transcriptome, their chemical labeling, and highlight structural and mechanistic aspects regarding the dihydrouridine synthases enzyme family. A special emphasis on important research directions to be addressed will also be discussed. This new entry of dihydrouridine into mRNA modifications has definitely added a new layer of information that controls protein synthesis.

NAT10-mediated N4-acetylcytidine modification is required for meiosis entry and progression in male germ cells

Post-transcriptional RNA modifications critically regulate various biological processes. N⁴-acetylcytidine (ac⁴C) is an epi-transcriptome, which is highly conserved in all species. However, the in vivo physiological functions and regulatory mechanisms of ac⁴C remain poorly understood, particularly in mammals. In this study, we demonstrate that the only known ac⁴C writer, N-acetyltransferase 10 (NAT10), plays an essential role in male reproduction. We identified the occurrence of ac⁴C in the mRNAs of mouse tissues and showed that ac⁴C undergoes dynamic changes during spermatogenesis. Germ cell-specific ablation of Nat10 severely inhibits meiotic entry and leads to defects in homologous chromosome synapsis, meiotic recombination and repair of DNA double-strand breaks during meiosis. Transcriptomic profiling revealed dysregulation of functional genes in meiotic prophase I after Nat10 deletion. These findings highlight the crucial physiological functions of ac⁴C modifications in male spermatogenesis and expand our understanding of its role in the regulation of specific physiological processes in vivo.

Identification and mapping of post-transcriptional modifications on the HIV-1 antisense transcript Ast in human cells

The human immunodeficiency virus type 1 (HIV-1) encodes multiple RNA molecules. Transcripts that originate from the proviral 5' long terminal repeat (LTR) function as messenger RNAs for the expression of 16 different mature viral proteins. In addition, HIV-1 expresses an antisense transcript (Ast) from the 3'LTR, which has both protein-coding and noncoding properties. While the mechanisms that regulate the coding and noncoding activities of Ast remain unknown, post-transcriptional modifications are known to influence RNA stability, interaction with protein partners, and translation capacity. Here, we report the nucleoside modification profile of Ast obtained through liquid chromatography coupled with mass spectrometry (LC-MS) analysis. The epitranscriptome includes a limited set of modified nucleosides but predominantly ribose methylations. A number of these modifications were mapped to specific positions of the sequence through RNA modification mapping procedures. The presence of modifications on Ast is consistent with the RNA-modifying enzymes interacting with Ast The identification and mapping of Ast post-transcriptional modifications is expected to elucidate the mechanisms through which this versatile molecule can carry out diverse activities in different cell compartments. Manipulation of post-transcriptional modifications on the Ast RNA may have therapeutic implications.

Chronic Hexavalent Chromium Exposure Upregulates the RNA Methyltransferase METTL3 Expression to Promote Cell Transformation, Cancer Stem Cell-Like Property, and Tumorigenesis

Hexavalent chromium [Cr(VI)] is a common environmental carcinogen causing lung cancer in humans. This study investigates the mechanism of Cr(VI) carcinogenesis focusing on the role of the epitranscriptomic dysregulation. The epitranscriptomic effect of Cr(VI) was determined in Cr(VI)-transformed human bronchial epithelial cells, chromate-exposed mouse and human lungs. The epitranscriptomic effect and its role in Cr(VI)-induced cell transformation, cancer stem cell (CSC)-like property, and tumorigenesis were determined by microarray analysis, soft agar colony formation, suspension spheroid formation, and mouse xenograft tumorigenesis assays. It was found that chronic Cr(VI) exposure causes epitranscriptomic dysregulations as evidenced by the increased levels of total RNA N6-methyladenosine (m6A) modification and the RNA m6A methyltransferase like-3 (METTL3) in Cr(VI)-transformed cells and chromate exposure-caused mouse and human lung tumors. Knockdown of METTL3 expression in Cr(VI)-transformed cells significantly reduces their m6A levels and transformed phenotypes and tumorigenicity in mice. Moreover, knockdown of METTL3 expression in parental nontransformed cells significantly reduces the capability of chronic Cr(VI) exposure to induce cell transformation and CSC-like property. Together, this study reveals that chronic Cr(VI) exposure is capable of altering cellular epitranscriptome by increasing the m6A RNA modification via upregulating the RNA methyltransferase METTL3 expression, which plays an important role in Cr(VI)-induced cell transformation, CSC-like property, and tumorigenesis.

m6A Regulator-Mediated Methylation Modification Patterns and Tumor Microenvironment Cell-Infiltration Characterization in Head and Neck Cancer

Background: Recently, RNA modifications have emerged as essential epigenetic regulators of gene expression. However, the mechanism of how RNA N⁶-methyladenosine (m6A) modification interacts with tumor microenvironment (TME) infiltration remains obscure.

Methods: A total of 876 head and neck cancer samples considering 21 m6A regulators were included and analyzed to determine the m6A modification patterns. These modification patterns were then correlated with TME immune cell-infiltrating characteristics. A scoring system, the m6Ascore, was constructed using principal component analysis algorithms to quantify m6A modification of tumors.

Results: Three m6A modification patterns were identified, with TME infiltrating characteristics highly consistent with tumors with three distinct immune phenotypes, including immune-inflamed, immune-exclude, and immune-desert phenotypes. It was demonstrated that the identification of the m6A modification patterns via m6Ascore could predict tumor progression, subtypes, TME stromal activity, variation of relevant genes, and patient prognosis. Low m6Ascore, identified to be an inflamed phenotype, is found to be associated with low stroma activity and tumor mutation burden, high survival probability, increased tumor neoantigen burden, and enhanced response to anti-PD-1/L1 immunotherapy. The therapeutic advantages and clinical benefits of patients with low m6Ascore were further verified in two immunotherapy cohorts.

Conclusion: This study identified the significant role that the m6A modification played in the formation of TME characteristics. A more comprehensive understanding of the m6A modification patterns and their correlation with TME infiltration will contribute to the discovery of immunotherapy strategies with better efficacy.

m6A: An Emerging Role in Programmed Cell Death

Programmed cell death is an active extinction process, including autophagy, ferroptosis, pyroptosis, apoptosis, and necroptosis. m6A is a reversible RNA modification which undergoes methylation under the action of methylases (writers), and is demethylated under the action of demethylases (erasers). The RNA base site at which m6A is modified is recognized by specialized enzymes (readers) which regulate downstream RNA translation, decay, and stability. m6A affects many aspects of mRNA metabolism, and also plays an important role in promoting the maturation of miRNA, the translation and degradation of circRNA, and the stability of lncRNA. The regulatory factors including writers, erasers and readers promote or inhibit programmed cell death via up-regulating or down-regulating downstream targets in a m6A-dependent manner to participate in the process of disease. In this review, we summarize the functions of m6A with particular reference to its role in programmed cell death.

N6-Methyladenosine-Related Gene Signature Associated With Monocyte Infiltration Is Clinically Significant in Gestational Diabetes Mellitus

OBJECTIVE

The objective of this study was to reveal the potential crosstalk between immune infiltration and N⁶- methyladenosine (m⁶A) modification in the placentas of patients with gestational diabetes mellitus (GDM), and to construct a model for the diagnosis of GDM.

METHODS

We analyzed imbalanced immune infiltration and differentially expressed m⁶A-related genes (DMRGs) in the placentas of patients with GDM, based on the GSE70493 dataset. An immune-related DMRG signature, with significant classifying power and diagnostic value, was identified using a least absolute shrinkage and selection operator (LASSO) regression. Based on the selected DMRGs, we developed and validated a nomogram model using GSE70493 and GSE92772 as the training and validation sets, respectively.

RESULTS

Infiltration of monocytes was higher in GDM placentas than in control samples, while the infiltration of macrophages (M1 and M2) in GDM placentas was lower than in controls. A total of 14 DMRGs were strongly associated with monocyte infiltration, seven of which were significant in distinguishing patients with GDM from normal controls. These genes were CD81, CFH, FABP5, GBP1, GNG11, IL1RL1, and SLAMF6. The calibration curve, decision curve, clinical impact curve, and receiver operating characteristic curve showed that the nomogram recognized GDM with high accuracy in both the training and validation sets.

CONCLUSIONS

Our results provide clues that crosstalk between m⁶A modification and immune infiltration may have implications in terms of novel biomarkers and therapeutic targets for GDM.

S-adenosylmethionine-dependent methyltransferase inhibitor DZNep blocks transcription and translation of SARS-CoV-2 genome with a low tendency to select for drug-resistant viral variants

We report the in vitro antiviral activity of DZNep (3-Deazaneplanocin A; an inhibitor of S-adenosylmethionine-dependent methyltransferase) against SARS-CoV-2, besides demonstrating its protective efficacy against lethal infection of infectious bronchitis virus (IBV, a member of the Coronaviridae family). DZNep treatment resulted in reduced synthesis of SARS-CoV-2 RNA and proteins without affecting other steps of viral life cycle. We demonstrated that deposition of N6-methyl adenosine (m6A) in SARS-CoV-2 RNA in the infected cells recruits heterogeneous nuclear ribonucleoprotein A1 (hnRNPA1), an RNA binding protein which serves as a m6A reader. DZNep inhibited the recruitment of hnRNPA1 at m6A-modified SARS-CoV-2 RNA which eventually suppressed the synthesis of the viral genome. In addition, m6A-marked RNA and hnRNPA1 interaction was also shown to regulate early translation to replication switch of SARS-CoV-2 genome. Furthermore, abrogation of methylation by DZNep also resulted in defective synthesis of the 5’ cap of viral RNA, thereby resulting in its failure to interact with eIF4E (a cap-binding protein), eventually leading to a decreased synthesis of viral proteins. Most importantly, DZNep-resistant mutants could not be observed upon long-term sequential passage of SARS-CoV-2 in cell culture. In summary, we report the novel role of methylation in the life cycle of SARS-CoV-2 and propose that targeting the methylome using DZNep could be of significant therapeutic value against SARS-CoV-2 infection.

Beyond sequencing: machine learning algorithms extract biology hidden in Nanopore signal data

Nanopore sequencing provides signal data corresponding to the nucleotide motifs sequenced. Through machine learning-based methods, these signals are translated into long-read sequences that overcome the read size limit of short-read sequencing. However, analyzing the raw nanopore signal data provides many more opportunities beyond just sequencing genomes and transcriptomes: algorithms that use machine learning approaches to extract biological information from these signals allow the detection of DNA and RNA modifications, the estimation of poly(A) tail length, and the prediction of RNA secondary structures. In this review, we discuss how developments in machine learning methodologies contributed to more accurate basecalling and lower error rates, and how these methods enable new biological discoveries. We argue that direct nanopore sequencing of DNA and RNA provides a new dimensionality for genomics experiments and highlight challenges and future directions for computational approaches to extract the additional information provided by nanopore signal data.

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

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

Human Mitochondrial Ribosomal RNA Modification-Based Classification Contributes to Discriminate the Prognosis and Immunotherapy Response of Glioma Patients

Background

Epigenetic regulations of the tumor microenvironment (TME) and immunotherapy have been investigated in recent years. Nevertheless, the potential value of mitochondrial ribosomal RNA (mt-rRNA) modification in regulation of the TME and immunotherapy remains unknown.

Methods

We comprehensively investigated the mt-rRNA-modification patterns in glioma patients based on nine regulators of mt-rRNA. Subsequently, these modification patterns were correlated systematically with immunologic characteristics and immunotherapy. An “mt-rRNA predictor” was constructed and validated in multiple publicly available cohorts to provide guidance for prognosis prediction and immunotherapy of glioma patients.

Results

Two distinct patterns of mt-rRNA modification were determined based on the evidence that nine regulators of mt-rRNA correlated significantly with most clinicopathologic characteristics, immunomodulators, TME, immune-checkpoint blockers (ICBs), and prognosis. Patients with mt-rRNA subtype II presented significantly poorer overall survival/progression-free survival (OS/PFS), but higher tumor mutational burden (TMB), more somatic mutations, and copy number variation (CNV). These two mt-rRNA subtypes had distinct TME patterns and responses to ICB therapy. An mt-rRNA predictor was constructed and validated in four glioma cohorts. The subtype with high mt-rRNA score, characterized by increased TMB, infiltration of immune cells, and activation of immunity, suggested an immune-activated phenotype, and was also linked to greater sensitivity to immunotherapy using anti-programmed cell death protein 1 (PD-1) but resistance to temozolomide.

Conclusions

Regulators of mt-rRNA modification have indispensable roles in the complexity and diversity of the TME and prognosis. This novel classification based on patterns of mt-rRNA modification could provide an effective prognostic predictor and guide more appropriate immunotherapy/chemotherapy strategies for glioma patients.

NAT10-Mediated N4-Acetylcytidine of RNA Contributes to Post-transcriptional Regulation of Mouse Oocyte Maturation in vitro

N4-acetylcytidine (ac4C), a newly identified epigenetic modification within mRNA, has been characterized as a crucial regulator of mRNA stability and translation efficiency. However, the role of ac4C during oocyte maturation, the process mainly controlled via post-transcriptional mechanisms, has not been explored. N-acetyltransferase 10 (NAT10) is the only known enzyme responsible for ac4C production in mammals and ac4C-binding proteins have not been reported yet. In this study, we have documented decreasing trends of both ac4C and NAT10 expression from immature to mature mouse oocytes. With NAT10 knockdown mediated by small interfering RNA (siRNA) in germinal vesicle (GV)-stage oocytes, ac4C modification was reduced and meiotic maturation in vitro was significantly retarded. Specifically, the rate of first polar body extrusion was significantly decreased with NAT10 knockdown (34.6%) compared to control oocytes without transfection (74.6%) and oocytes transfected with negative control siRNA (72.6%) (p < 0.001), while rates of germinal vesicle breakdown (GVBD) were not significantly different (p = 0.6531). RNA immunoprecipitation and high-throughput sequencing using HEK293T cells revealed that the modulated genes were enriched in biological processes associated with nucleosome assembly, chromatin silencing, chromatin modification and cytoskeletal anchoring. In addition, we identified TBL3 as a potential ac4C-binding protein by a bioinformatics algorithm and RNA pulldown with HEK293T cells, which may mediate downstream cellular activities. Taken together, our results suggest that NAT10-mediated ac4C modification is an important regulatory factor during oocyte maturation in vitro and TBL3 is a potential ac4C-binding protein.

RNA methyltransferase METTL16: Targets and function

The N6-methyladenosine (m6A) RNA methyltransferase METTL16 is an emerging player in the RNA modification landscape of the human cell. Originally thought to be a ribosomal RNA methyltransferase, it has now been shown to bind and methylate the MAT2A messenger RNA (mRNA) and U6 small nuclear RNA (snRNA). It has also been shown to bind the MALAT1 long noncoding RNA and several other RNAs. METTL16's methyltransferase domain contains the Rossmann-like fold of class I methyltransferases and uses S-adenosylmethionine (SAM) as the methyl donor. It has an RNA methylation consensus sequence of UACAGARAA (modified A underlined), and structural requirements for its known RNA interactors. In addition to the methyltransferase domain, METTL16 protein has two other RNA binding domains, one of which resides in a vertebrate conserved region, and a putative nuclear localization signal. The role of METTL16 in the cell is still being explored, however evidence suggests it is essential for most cells. This is currently hypothesized to be due to its role in regulating the splicing of MAT2A mRNA in response to cellular SAM levels. However, one of the more pressing questions remaining is what role METTL16's methylation of U6 snRNA plays in splicing and potentially cellular survival. METTL16 also has several other putative coding and noncoding RNA interactors but the definitive methylation status of those RNAs and the role METTL16 plays in their life cycle is yet to be determined. Overall, METTL16 is an intriguing RNA binding protein and methyltransferase whose important functions in the cell are just beginning to be understood. This article is categorized under: RNA Processing > RNA Editing and Modification RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.