Transcriptome-wide Mapping of Internal N7-Methylguanosine Methylome in Mammalian mRNA

RNA-related DNA damage and repair: The role of N7-methylguanosine in the cell nucleus exposed to UV light

Background

Chemical modifications in mRNAs, tRNAs, rRNAs, and non-coding RNAs stabilize these nucleic acids and regulate their function. In addition to regulating the translation of genetic information from mRNA to proteins, it has been revealed that modifications in RNAs regulate repair processes in the genome.

Methods

Using local laser microirradiation, confocal microscopy, dot blots, and mass spectrometry we studied the role of N7-methylguanosine (m7G), which is co-transcriptionally installed in RNA.

Results

Here, we show that after UVC and UVA irradiation, the level of m7G RNA is increased initially in the cytoplasm, and after local laser microirradiation, m7G RNA is highly abundant in UVA-damaged chromatin. This process is poly(ADP-ribose) polymerase (PARP)-dependent, but not accompanied by changes in the level of m7G-writers, including methyltransferases RNMT, METTL1, and WBSCR22. We also observed that METTL1 deficiency does not affect the recruitment of m7G RNA to microirradiated chromatin. Analyzing the levels of mRNA, let-7e, and miR-203a in both the cytoplasm and the cell nucleus, we revealed that UVC irradiation changed the level of mRNA, and significantly increased the pool of both let-7e and miR-203a, which correlated with radiation-induced m7G RNA increase in the cytoplasm.

Conclusions

Irradiation by UV light increases the m7G RNA pool in the cytoplasm and in the microirradiated genome. Thus, epigenetically modified RNAslikely contribute to DNA damage responses or m7G signals the presence of RNA damage.

Prognostic value of RNA methylation-related genes in gastric adenocarcinoma based on bioinformatics

Background

Gastric cancer (GC) is a malignant tumor that originates from the epithelium of the gastric mucosa and has a poor prognosis. Stomach adenocarcinoma (STAD) covers 95% of total gastric cancer. This study aimed to identify the prognostic value of RNA methylation-related genes in gastric cancer.

Methods

In this study, The Cancer Genome Atlas (TCGA)-STAD and GSE84426 cohorts were downloaded from public databases. Patients were classified by consistent cluster analysis based on prognosis-related differentially expressed RNA methylation genes Prognostic genes were obtained by differential expression, univariate Cox and least absolute shrinkage and selection operator (LASSO) analyses. The prognostic model was established and validated in the training set, test set and validation set respectively. Independent prognostic analysis was implemented. Finally, the expression of prognostic genes was affirmed by reverse transcription quantitative PCR (RT-qPCR).

Results

In total, four prognostic genes (ACTA2, SAPCD2, PDK4 and APOD) related to RNA methylation were identified and enrolled into the risk signature. The STAD patients were divided into high- and low-risk groups based on the medium value of the risk score, and patients in the high-risk group had a poor prognosis. In addition, the RNA methylation-relevant risk signature was validated in the test and validation sets, and was authenticated as a reliable independent prognostic predictor. The nomogram was constructed based on the independent predictors to predict the 1/3/5-year survival probability of STAD patients. The gene set enrichment analysis (GSEA) result suggested that the poor prognosis in the high-risk subgroup may be related to immune-related pathways. Finally, the experimental results indicated that the expression trends of RNA methylation-relevant prognostic genes in gastric cancer cells were in agreement with the result of bioinformatics.

Conclusion

Our study established a novel RNA methylation-related risk signature for STAD, which was of considerable significance for improving prognosis of STAD patients and offering theoretical support for clinical therapy.

General Principles and Limitations for Detection of RNA Modifications by Sequencing

Among the many analytical methods applied to RNA modifications, a particularly pronounced surge has occurred in the past decade in the field of modification mapping. The occurrence of modifications such as m6A in mRNA, albeit known since the 1980s, became amenable to transcriptome-wide analyses through the advent of next-generation sequencing techniques in a rather sudden manner. The term "mapping" here refers to detection of RNA modifications in a sequence context, which has a dramatic impact on the interpretation of biological functions. As a consequence, an impressive number of mapping techniques were published, most in the perspective of what now has become known as "epitranscriptomics". While more and more different modifications were reported to occur in mRNA, conflicting reports and controversial results pointed to a number of technical and theoretical problems rooted in analytics, statistics, and reagents. Rather than finding the proverbial needle in a haystack, the tasks were to determine how many needles of what color in what size of a haystack one was looking at.As the authors of this Account, we think it important to outline the limitations of different mapping methods since many life scientists freshly entering the field confuse the accuracy and precision of modification mapping with that of normal sequencing, which already features numerous caveats by itself. Indeed, we propose here to qualify a specific mapping method by the size of the transcriptome that can be meaningfully analyzed with it.We here focus on high throughput sequencing by Illumina technology, referred to as RNA-Seq. We noted with interest the development of methods for modification detection by other high throughput sequencing platforms that act directly on RNA, e.g., PacBio SMRT and nanopore sequencing, but those are not considered here.In contrast to approaches relying on direct RNA sequencing, current Illumina RNA-Seq protocols require prior conversion of RNA into DNA. This conversion relies on reverse transcription (RT) to create cDNA; thereafter, the cDNA undergoes a sequencing-by-synthesis type of analysis. Thus, a particular behavior of RNA modified nucleotides during the RT-step is a prerequisite for their detection (and quantification) by deep sequencing, and RT properties have great influence on the detection efficiency and reliability. Moreover, the RT-step requires annealing of a synthetic primer, a prerequisite with a crucial impact on library preparation. Thus, all RNA-Seq protocols must feature steps for the introduction of primers, primer landing sites, or adapters on both the RNA 3'- and 5'-ends.

BID-seq for transcriptome-wide quantitative sequencing of mRNA pseudouridine at base resolution

Pseudouridine (Ψ) is an abundant RNA modification that is present in and affects the functions of diverse non-coding RNA species, including rRNA, tRNA and small nuclear RNA. Ψ also exists in mammalian mRNA and probably exhibits functional roles; however, functional investigations of mRNA Ψ modifications in mammals have been hampered by the lack of a quantitative method that detects Ψ at base precision. We have recently developed bisulfite-induced deletion sequencing (BID-seq), which provides the community with a quantitative method to map RNA Ψ distribution transcriptome-wide at single-base resolution. Here, we describe an optimized BID-seq protocol for mapping Ψ distribution across cellular mRNAs, which includes fast steps in both library preparation and data analysis. This protocol generates highly reproducible results by inducing high deletion ratios at Ψ modification within diverse sequence contexts, and meanwhile displayed almost zero background deletions at unmodified uridines. When used for transcriptome-wide Ψ profiling in mouse embryonic stem cells, the current protocol uncovered 8,407 Ψ sites from as little as 10 ng of polyA+ RNA input. This optimized BID-seq workflow takes 5 days to complete and includes four main sections: RNA preparation, library construction, next-generation sequencing (NGS) and data analysis. Library construction can be completed by researchers who have basic knowledge and skills in molecular biology and genetics. In addition to the experimental protocol, we provide BID-pipe ( https://github.com/y9c/pseudoU-BIDseq ), a user-friendly data analysis pipeline for Ψ site detection and modification stoichiometry quantification, requiring only basic bioinformatic and computational skills to uncover Ψ signatures from BID-seq data.

M7G-related tumor immunity: novel insights of RNA modification and potential therapeutic targets

RNA modifications play a pivotal role in regulating cellular biology by exerting influence over distribution features and molecular functions at the post-transcriptional level. Among these modifications, N7-methylguanosine (m7G) stands out as one of the most prevalent. Over recent years, significant attention has been directed towards understanding the implications of m7G modification. This modification is present in diverse RNA molecules, including transfer RNAs, messenger RNAs, ribosomal RNAs, and other noncoding RNAs. Its regulation occurs through a series of specific methyltransferases and m7G-binding proteins. Notably, m7G modification has been implicated in various diseases, prominently across multiple cancer types. Earlier studies have elucidated the significance of m7G modification in the context of immune biology regulation within the tumor microenvironment. This comprehensive review culminates in a synthesis of findings related to the modulation of immune cells infiltration, encompassing T cells, B cells, and various innate immune cells, all orchestrated by m7G modification. Furthermore, the interplay between m7G modification and its regulatory proteins can profoundly affect the efficacy of diverse adjuvant therapeutics, thereby potentially serving as a pivotal biomarker and therapeutic target for combinatory interventions in diverse cancer types.

METTL1 mediated tRNA m7G modification promotes leukaemogenesis of AML via tRNA regulated translational control

BACKGROUND

RNA modifications have been proven to play fundamental roles in regulating cellular biology process. Recently, maladjusted N7-methylguanosine (m7G) modification and its modifiers METTL1/WDR4 have been confirmed an oncogene role in multiple cancers. However, the functions and molecular mechanisms of METTL1/WDR4 in acute myeloid leukemia (AML) remain to be determined.

METHODS

METTL1/WDR4 expression levels were quantified using qRT-PCR, western blot analysis on AML clinical samples, and bioinformatics analysis on publicly available AML datasets. CCK-8 assays and cell count assays were performed to determine cell proliferation. Flow cytometry assays were conducted to assess cell cycle and apoptosis rates. Multiple techniques were used for mechanism studies in vitro assays, such as northern blotting, liquid chromatography-coupled mass spectrometry (LC-MS/MS), tRNA stability analysis, transcriptome sequencing, small non-coding RNA sequencing, quantitative proteomics, and protein synthesis measurements.

RESULTS

METTL1/WDR4 are significantly elevated in AML patients and associated with poor prognosis. METTL1 knockdown resulted in reduced cell proliferation and increased apoptosis in AML cells. Mechanically, METTL1 knockdown leads to significant decrease of m7G modification abundance on tRNA, which further destabilizes tRNAs and facilitates the biogenesis of tsRNAs in AML cells. In addition, profiling of nascent proteins revealed that METTL1 knockdown and transfection of total tRNAs that were isolated from METTL1 knockdown AML cells decreased global translation efficiency in AML cells.

CONCLUSIONS

Taken together, our study demonstrates the important role of METTL1/WDR4 in AML leukaemogenesis, which provides a promising target candidate for AML therapy.

Genetic variants of m7G modification genes influence neuroblastoma susceptibility

Objective

Neuroblastoma is a life-threatening pediatric solid tumor whose etiology remains unclear. N7-methylguanosine (m7G) is one of the most important epigenetic modifications of RNA, which plays a crucial role in tumorigenesis. The m7G-mediated genes METTL1 and WDR4 also have been reported to be dysregulated in various cancers. However, the implications of METTL1 and WDR4 in neuroblastoma have not been clarified.

Methods

Given the oncogenic potential of m7G modification, we performed a case-control study to assess the association of METTL1 and WDR4 genes polymorphisms with neuroblastoma risk in a Chinese population consisting of 402 cases and 473 controls. Odds ratios (ORs) and 95 % confidence intervals (CIs) were applied to evaluate the associations between studied polymorphisms and neuroblastoma risk. The adjusted odds ratio (AOR) was adjusted for age and gender.

Results

Overall, four polymorphisms were significantly associated with neuroblastoma risk, including METTL1 rs2291617 (recessive model: adjusted OR = 1.59, 95 % CI = 1.08-2.34, P = 0.019), WDR4 rs2156316 (dominant model: adjusted OR = 0.74, 95 % CI = 0.57-0.97, P = 0.028), WDR4 rs6586250 (dominant model: adjusted OR = 0.59, 95 % CI = 0.42-0.84, P = 0.004) and WDR4 rs15736 (dominant model: adjusted OR = 0.60, 95 % CI = 0.42-0.85, P = 0.004). Stratified analysis showed stronger correlations between significant polymorphisms and neuroblastoma risk among subgroups divided by age, gender, tumor origin, and clinical stage. Furthermore, expression quantitative trait loci (eQTL) analysis revealed that significant polymorphisms were associated with the expression of the adjacent genes.

Conclusions

Our study indicated that four polymorphisms in m7G-mediated genes contribute to neuroblastoma susceptibility in the eastern Chinese population. However, our findings should be verified further by large-scale and well-designed studies.

RMBase v3.0: decode the landscape, mechanisms and functions of RNA modifications

Although over 170 chemical modifications have been identified, their prevalence, mechanism and function remain largely unknown. To enable integrated analysis of diverse RNA modification profiles, we have developed RMBase v3.0 (http://bioinformaticsscience.cn/rmbase/), a comprehensive platform consisting of eight modules. These modules facilitate the exploration of transcriptome-wide landscape, biogenesis, interactome and functions of RNA modifications. By mining thousands of epitranscriptome datasets with novel pipelines, the 'RNA Modifications' module reveals the map of 73 RNA modifications of 62 species. the 'Genes' module allows to retrieve RNA modification profiles and clusters by gene and transcript. The 'Mechanisms' module explores 23 382 enzyme-catalyzed or snoRNA-guided modified sites to elucidate their biogenesis mechanisms. The 'Co-localization' module systematically formulates potential correlations between 14 histone modifications and 6 RNA modifications in various cell-lines. The 'RMP' module investigates the differential expression profiles of 146 RNA-modifying proteins (RMPs) in 18 types of cancers. The 'Interactome' integrates the interactional relationships between 73 RNA modifications with RBP binding events, miRNA targets and SNPs. The 'Motif' illuminates the enriched motifs for 11 types of RNA modifications identified from epitranscriptome datasets. The 'Tools' introduces a novel web-based 'modGeneTool' for annotating modifications. Overall, RMBase v3.0 provides various resources and tools for studying RNA modifications.

m7GHub V2.0: an updated database for decoding the N7-methylguanosine (m7G) epitranscriptome

With recent progress in mapping N7-methylguanosine (m7G) RNA methylation sites, tens of thousands of experimentally validated m7G sites have been discovered in various species, shedding light on the significant role of m7G modification in regulating numerous biological processes including disease pathogenesis. An integrated resource that enables the sharing, annotation and customized analysis of m7G data will greatly facilitate m7G studies under various physiological contexts. We previously developed the m7GHub database to host mRNA m7G sites identified in the human transcriptome. Here, we present m7GHub v.2.0, an updated resource for a comprehensive collection of m7G modifications in various types of RNA across multiple species: an m7GDB database containing 430 898 putative m7G sites identified in 23 species, collected from both widely applied next-generation sequencing (NGS) and the emerging Oxford Nanopore direct RNA sequencing (ONT) techniques; an m7GDiseaseDB hosting 156 206 m7G-associated variants (involving addition or removal of an m7G site), including 3238 disease-relevant m7G-SNPs that may function through epitranscriptome disturbance; and two enhanced analysis modules to perform interactive analyses on the collections of m7G sites (m7GFinder) and functional variants (m7GSNPer). We expect that m7Ghub v.2.0 should serve as a valuable centralized resource for studying m7G modification. It is freely accessible at: www.rnamd.org/m7GHub2.

Base-Resolution Sequencing Methods for Whole-Transcriptome Quantification of mRNA Modifications

RNA molecules are not merely a combination of four bases of A, C, G, and U. Chemical modifications occur in almost all RNA species and play diverse roles in gene expression regulation. The abundant cellular RNAs, such as ribosomal RNA (rRNA) and transfer RNA (tRNA), are known to have the highest density of RNA modifications, which exert critical functions in rRNA and tRNA biogenesis, stability, and subsequent translation. In recent years, modifications on low-abundance RNA species in mammalian cells, such as messenger RNA (mRNA), regulatory noncoding RNA (ncRNA), and chromatin-associated RNA (caRNA), have been shown to contain multiple different chemical modifications with functional significance. As the most abundant mRNA modification in mammals, N6-methyladenosine (m6A) affects nearly every stage of mRNA processing and metabolism, with the antibody-based m6A-MeRIP-seq (methylated RNA immunoprecipitation sequencing) followed by high-throughput sequencing widely employed in mapping m6A distribution transcriptome-wide in diverse biological systems. In addition to m6A, other chemical modifications such as pseudouridine (Ψ), 2'-O-methylation (Nm), 5-methylcytidine (m5C), internal N7-methylguanosine (m7G), N1-methyladenosine (m1A), N4-acetylcytidine (ac4C), etc. also exist in polyA-tailed RNA in mammalian cells, requiring effective mapping approaches for whole-transcriptome profiling of these non-m6A mRNA modifications. Like m6A, the antibody-based enrichment followed by sequencing has been the primary method to study distributions of these modifications. Methods to more quantitatively map these modifications would dramatically improve our understanding of distributions and modification density of these chemical marks on RNA, thereby bettering informing functional implications. In this Account, aimed at both single-base resolution and modification fraction quantification, we summarize our recent advances in developing a series of chemistry- or biochemistry-based methods to quantitatively map RNA modifications, including m6A, Ψ, m5C, m1A, 2'-O-methylation (Nm), and internal m7G, in mammalian mRNA at base resolution. These new methods, including m6A-SAC-seq, eTAM-seq, BID-seq, UBS-seq, DAMM-seq, m1A-quant-seq, Nm-Mut-seq, and m7G-quant-seq, promise to conduct base-resolution mapping of most major mRNA modifications with low RNA input and uncover dynamic changes in modification stoichiometry during biological and physiological processes, facilitating future investigations on these RNA modifications in regulating cellular gene expression and as potential biomarkers for clinical diagnosis and prognosis. These quantitative sequencing methods allow the mapping of most mRNA modifications with limited input sample requirements. The same modifications on diverse RNA species, such as caRNA, ncRNA, nuclear nascent RNA, mitochondrial RNA, cell-free RNA (cfRNA), etc., could be sequenced using the same methods.

Emerging roles of RNA ac4C modification and NAT10 in mammalian development and human diseases

RNA ac4C modification is a novel and rare chemical modification observed in mRNA. Traditional biochemical studies had primarily associated ac4C modification with tRNA and rRNA until in 2018, Arango D et al. first reported the presence of ac4C modification on mRNA and demonstrated its critical role in mRNA stability and translation regulation. Furthermore, they established that the ac4C modification on mRNA is mediated by the classical N-acetyltransferase NAT10. Subsequent studies have underscored the essential implications of NAT10 and mRNA ac4C modification across both physiological and pathological regulatory processes. In this review, we aimed to explore the discovery history of RNA ac4C modification, its detection methods, and its regulatory mechanisms in disease and physiological development. We offer a forward-looking examination and discourse concerning the employment of RNA ac4C modification as a prospective therapeutic strategy across diverse diseases. Furthermore, we comprehensively summarize the functions and mechanisms of NAT10 in gene expression regulation and pathogenesis independent of RNA ac4C modification.

M7G modification of FTH1 and pri-miR-26a regulates ferroptosis and chemotherapy resistance in osteosarcoma

Doxorubicin and platinum are widely used in the frontline treatment of osteosarcoma, but resistance to chemotherapy limits its curative effect. Here, we have identified that METTL1 mediated N7-Methyladenosine (m7G) low expressed in osteosarcoma tissues, plays a critical oncogenic role, and enhances osteosarcoma chemosensitivity in osteosarcoma. Mechanistically, AlkAniline-Seq data revealed that Ferritin heavy chain (FTH1), the main component of ferritin, which is crucial for iron homeostasis and the inhibition of lipid peroxidation, is one of the top 10 genes with the most significant change in m7G methylation sites mediated by METTL1 in human osteosarcoma cells. Interestingly, METTL1 significantly increased the expression of FTH1 at the mRNA level but was remarkably suppressed at the protein level. We then identified primary (pri)-miR-26a and pri-miR-98 in the Top 20 m7G-methylated pri-miRNAs with highly conserved species. Further results confirmed that METTL1 enhances cell ferroptosis by targeting FTH1 and primary (pri)-miR-26a, promoting their maturity by enhancing RNA stability dependent on m7G methylation. The increase of mature miR-26a-5p that resulted from METTL1 overexpression could further target FTH1 mRNA and eliminate FTH1 translation efficiency. Moreover, the reduction of FTH1 translation dramatically increases cell ferroptosis and promotes the sensitivity of osteosarcoma cells to chemotherapy drugs. Collectively, our study demonstrates the METTL1/pri-miR-26a/FTH1 axis signaling in osteosarcoma and highlights the functional importance of METTL1 and m7G methylation in the progression and chemotherapy resistance of osteosarcoma, suggesting that reprogramming RNA m7G methylation as a potential and promising strategy for osteosarcoma treatment.

Global and single-nucleotide resolution detection of 7-methylguanosine in RNA

RNA modifications, including N-7-methylguanosine (m7G), are pivotal in governing RNA stability and gene expression regulation. The accurate detection of internal m7G modifications is of paramount significance, given recent associations between altered m7G deposition and elevated expression of the methyltransferase METTL1 in various human cancers. The development of robust m7G detection techniques has posed a significant challenge in the field of epitranscriptomics. In this study, we introduce two methodologies for the global and accurate identification of m7G modifications in human RNA. We introduce borohydride reduction sequencing (Bo-Seq), which provides base resolution mapping of m7G modifications. Bo-Seq achieves exceptional performance through the optimization of RNA depurination and scission, involving the strategic use of high concentrations of NaBH4, neutral pH and the addition of 7-methylguanosine monophosphate (m7GMP) during the reducing reaction. Notably, compared to NaBH4-based methods, Bo-Seq enhances the m7G detection performance, and simplifies the detection process, eliminating the necessity for intricate chemical steps and reducing the protocol duration. In addition, we present an antibody-based approach, which enables the assessment of m7G relative levels across RNA molecules and biological samples, however it should be used with caution due to limitations associated with variations in antibody quality between batches. In summary, our novel approaches address the pressing need for reliable and accessible methods to detect RNA m7G methylation in human cells. These advancements hold the potential to catalyse future investigations in the critical field of epitranscriptomics, shedding light on the complex regulatory roles of m7G in gene expression and its implications in cancer biology.

Binding asymmetry and conformational studies of the AtGSDA dimer

Guanosine deaminase (GSDA) is an important deaminase that converts guanosine to xanthosine, a key intermediate in nitrogen recycling in plants. We previously solved complex structures of Arabidopsis thaliana GSDA bound by various ligands and examined its catalytic mechanism. Here, we report cocrystal structures of AtGSDA bound by inactive guanosine derivatives, which bind relatively weakly to the enzyme and mostly have poor binding geometries. The two protomers display unequal binding performances, and molecular dynamics simulation identified diverse conformations during the enzyme-ligand interactions. Moreover, intersubunit, tripartite salt bridges show conformational differences between the two protomers, possibly acting as "gating" systems for substrate binding and product release. Our structural and biochemical studies provide a comprehensive understanding of the enzymatic behavior of this intriguing enzyme.

N7-methylguanosine-related miRNAs predict hepatocellular carcinoma prognosis and immune therapy

N7-methylguanosine (m7G) modification has been notably linked with the development of many tumors. However, no investigations have been conducted on whether m7G-related miRNA (m7G-miRNA) is a prognostic index of hepatocellular carcinoma (HCC). Therefore, this investigation aimed to establish a predictive m7G-miRNA signature for efficient HCC prognosis and elucidate the associated immune cell infiltration (ICI) and functions in the tumor microenvironment. RNA sequencing and clinical data on 375 HCC and 50 healthy tissue samples were acquired from The Cancer Genome Atlas database. The m7G-miRNA regulators methyltransferase-like 1 and WD repeat domain 4 were acquired from the TargetScan database. Univariate Cox regression analysis was conducted on the 63 differentially expressed m7G-miRNAs identified. A prognostic signature that consisted of seven miRNAs was identified. According to their risk scores, individuals with HCC were divided into high-risk (HR) and low-risk (LR) cohorts. A Kaplan-Meier test revealed that survival in the HR HCC patients was poorer than in the LR cohort (p < 0.001). The area under the receiver operating characteristic curves of 1-, 3-, and 5-year overall survival were 0.706, 0.695, and 0.715, respectively. A nomogram of sex, risk score, age, and stage indicated the HCC patients' overall survival. Furthermore, it was indicated that the HR and LR patients had different degrees of ICI and immune function. A pathway enrichment analysis revealed the association of several immunity-linked pathways with the risk model. In conclusion, the signature established has great prognostic value and could be used as a new immunotherapy target for individuals with HCC.

RNA modification in cardiovascular disease: implications for therapeutic interventions

Cardiovascular disease (CVD) is the leading cause of death in the world, with a high incidence and a youth-oriented tendency. RNA modification is ubiquitous and indispensable in cell, maintaining cell homeostasis and function by dynamically regulating gene expression. Accumulating evidence has revealed the role of aberrant gene expression in CVD caused by dysregulated RNA modification. In this review, we focus on nine common RNA modifications: N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C), N7-methylguanosine (m7G), N4-acetylcytosine (ac4C), pseudouridine (Ψ), uridylation, adenosine-to-inosine (A-to-I) RNA editing, and modifications of U34 on tRNA wobble. We summarize the key regulators of RNA modification and their effects on gene expression, such as RNA splicing, maturation, transport, stability, and translation. Then, based on the classification of CVD, the mechanisms by which the disease occurs and progresses through RNA modifications are discussed. Potential therapeutic strategies, such as gene therapy, are reviewed based on these mechanisms. Herein, some of the CVD (such as stroke and peripheral vascular disease) are not included due to the limited availability of literature. Finally, the prospective applications and challenges of RNA modification in CVD are discussed for the purpose of facilitating clinical translation. Moreover, we look forward to more studies exploring the mechanisms and roles of RNA modification in CVD in the future, as there are substantial uncultivated areas to be explored.

RNA modification regulator DDC in endometrial cancer affects the tumor microenvironment and patient prognosis

Uterine corpus endometrial carcinoma (UCEC) is infiltrated by immune cells, which are involved in the growth and proliferation of malignant tumors and resistance to immunotherapy. This study suggested that RNA modification regulators played an important role in the development and prognosis of UCEC. Many studies confirmed that RNA modification played an essential role in tumor immune regulation, and abnormal RNA modification contributed to tumorigenesis and cancer progression. Based on the RNA modification regulatory factors, the UCEC samples from TCGA (The Cancer Genome Atlas) were classified into two clusters, namely Cluster A and Cluster B, using unsupervised consensus clustering. We obtained DEG (differentially expressed genes) between the two clusters, and constructed a risk model of RNA modification-related genes using DEGs. Cluster A had lower RNA modification regulatory factors, richer immune cell infiltration, and better prognosis. The differentially expressed genes between the two clusters were obtained, and these genes were used for modeling. This model divided patients with UCEC into two groups. The low-risk group had better immune infiltration, and the ROC (receiver operating characteristic) curve showed that this model had good predictive efficacy. The low-risk group had a better response to immunotherapy by immune checkpoint prediction. We obtained the key gene L-dopa decarboxylase (DDC) through the intersection of LASSO model genes and GEO dataset GSE17025. We evaluated the potential biological functions of DDC. The differences in the expression of DDC were verified by immunohistochemistry. We evaluated the relationship between DDC and immune cell infiltration and verified this difference using immunofluorescence. Cluster A with low expression of RNA modification regulators has better prognosis and richer immune cell infiltration, therefore, we believed that RNA modification regulators in UCEC were closely related to the tumor microenvironment. Also, the risk score could well predict the prognosis of patients and guide immunotherapy, which might benefit patients with UCEC.

Exploring the brain epitranscriptome: perspectives from the NSAS summit

Increasing evidence reinforces the essential function of RNA modifications in development and diseases, especially in the nervous system. RNA modifications impact various processes in the brain, including neurodevelopment, neurogenesis, neuroplasticity, learning and memory, neural regeneration, neurodegeneration, and brain tumorigenesis, leading to the emergence of a new field termed neuroepitranscriptomics. Deficiency in machineries modulating RNA modifications has been implicated in a range of brain disorders from microcephaly, intellectual disability, seizures, and psychiatric disorders to brain cancers such as glioblastoma. The inaugural NSAS Challenge Workshop on Brain Epitranscriptomics hosted in Crans-Montana, Switzerland in 2023 assembled a group of experts from the field, to discuss the current state of the field and provide novel translational perspectives. A summary of the discussions at the workshop is presented here to simulate broader engagement from the general neuroscience field.

The regulation of antiviral innate immunity through non-m6A RNA modifications

The post-transcriptional RNA modifications impact the dynamic regulation of gene expression in diverse biological and physiological processes. Host RNA modifications play an indispensable role in regulating innate immune responses against virus infection in mammals. Meanwhile, the viral RNAs can be deposited with RNA modifications to interfere with the host immune responses. The N6-methyladenosine (m6A) has boosted the recent emergence of RNA epigenetics, due to its high abundance and a transcriptome-wide widespread distribution in mammalian cells, proven to impact antiviral innate immunity. However, the other types of RNA modifications are also involved in regulating antiviral responses, and the functional roles of these non-m6A RNA modifications have not been comprehensively summarized. In this Review, we conclude the regulatory roles of 2'-O-methylation (Nm), 5-methylcytidine (m5C), adenosine-inosine editing (A-to-I editing), pseudouridine (Ψ), N1-methyladenosine (m1A), N7-methylguanosine (m7G), N6,2'-O-dimethyladenosine (m6Am), and N4-acetylcytidine (ac4C) in antiviral innate immunity. We provide a systematic introduction to the biogenesis and functions of these non-m6A RNA modifications in viral RNA, host RNA, and during virus-host interactions, emphasizing the biological functions of RNA modification regulators in antiviral responses. Furthermore, we discussed the recent research progress in the development of antiviral drugs through non-m6A RNA modifications. Collectively, this Review conveys knowledge and inspiration to researchers in multiple disciplines, highlighting the challenges and future directions in RNA epitranscriptome, immunology, and virology.

Epitranscriptomics as a New Layer of Regulation of Gene Expression in Skeletal Muscle: Known Functions and Future Perspectives

Epitranscriptomics refers to post-transcriptional regulation of gene expression via RNA modifications and editing that affect RNA functions. Many kinds of modifications of mRNA have been described, among which are N6-methyladenosine (m6A), N1-methyladenosine (m1A), 7-methylguanosine (m7G), pseudouridine (Ψ), and 5-methylcytidine (m5C). They alter mRNA structure and consequently stability, localization and translation efficiency. Perturbation of the epitranscriptome is associated with human diseases, thus opening the opportunity for potential manipulations as a therapeutic approach. In this review, we aim to provide an overview of the functional roles of epitranscriptomic marks in the skeletal muscle system, in particular in embryonic myogenesis, muscle cell differentiation and muscle homeostasis processes. Further, we explored high-throughput epitranscriptome sequencing data to identify RNA chemical modifications in muscle-specific genes and we discuss the possible functional role and the potential therapeutic applications.

Shining a spotlight on m6A and the vital role of RNA modification in endometrial cancer: a review

RNA modifications are mostly dynamically reversible post-transcriptional modifications, of which m6A is the most prevalent in eukaryotic mRNAs. A growing number of studies indicate that RNA modification can finely tune gene expression and modulate RNA metabolic homeostasis, which in turn affects the self-renewal, proliferation, apoptosis, migration, and invasion of tumor cells. Endometrial carcinoma (EC) is the most common gynecologic tumor in developed countries. Although it can be diagnosed early in the onset and have a preferable prognosis, some cases might develop and become metastatic or recurrent, with a worse prognosis. Fortunately, immunotherapy and targeted therapy are promising methods of treating endometrial cancer patients. Gene modifications may also contribute to these treatments, as is especially the case with recent developments of new targeted therapeutic genes and diagnostic biomarkers for EC, even though current findings on the relationship between RNA modification and EC are still very limited, especially m6A. For example, what is the elaborate mechanism by which RNA modification affects EC progression? Taking m6A modification as an example, what is the conversion mode of methylation and demethylation for RNAs, and how to achieve selective recognition of specific RNA? Understanding how they cope with various stimuli as part of in vivo and in vitro biological development, disease or tumor occurrence and development, and other processes is valuable and RNA modifications provide a distinctive insight into genetic information. The roles of these processes in coping with various stimuli, biological development, disease, or tumor development in vivo and in vitro are self-evident and may become a new direction for cancer in the future. In this review, we summarize the category, characteristics, and therapeutic precis of RNA modification, m6A in particular, with the purpose of seeking the systematic regulation axis related to RNA modification to provide a better solution for the treatment of EC.

Combined signature of N7-methylguanosine regulators with their related genes and the tumor microenvironment: a prognostic and therapeutic biomarker for breast cancer

Background

Identifying predictive markers for breast cancer (BC) prognosis and immunotherapeutic responses remains challenging. Recent findings indicate that N7-methylguanosine (m7G) modification and the tumor microenvironment (TME) are critical for BC tumorigenesis and metastasis, suggesting that integrating m7G modifications and TME cell characteristics could improve the predictive accuracy for prognosis and immunotherapeutic responses.

Methods

We utilized bulk RNA-sequencing data from The Cancer Genome Atlas Breast Cancer Cohort and the GSE42568 and GSE146558 datasets to identify BC-specific m7G-modification regulators and associated genes. We used multiple m7G databases and RNA interference to validate the relationships between BC-specific m7G-modification regulators (METTL1 and WDR4) and related genes. Single-cell RNA-sequencing data from GSE176078 confirmed the association between m7G modifications and TME cells. We constructed an m7G-TME classifier, validated the results using an independent BC cohort (GSE20685; n = 327), investigated the clinical significance of BC-specific m7G-modifying regulators by reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis, and performed tissue-microarray assays on 192 BC samples.

Results

Immunohistochemistry and RT-qPCR results indicated that METTL1 and WDR4 overexpression in BC correlated with poor patient prognosis. Moreover, single-cell analysis revealed relationships between m7G modification and TME cells, indicating their potential as indicators of BC prognosis and treatment responses. The m7G-TME classifier enabled patient subgrouping and revealed significantly better survival and treatment responses in the m7Glow+TMEhigh group. Significant differences in tumor biological functions and immunophenotypes occurred among the different subgroups.

Conclusions

The m7G-TME classifier offers a promising tool for predicting prognosis and immunotherapeutic responses in BC, which could support personalized therapeutic strategies.

Regulation and functions of non-m6A mRNA modifications

Nucleobase modifications are prevalent in eukaryotic mRNA and their discovery has resulted in the emergence of epitranscriptomics as a research field. The most abundant internal (non-cap) mRNA modification is N6-methyladenosine (m6A), the study of which has revolutionized our understanding of post-transcriptional gene regulation. In addition, numerous other mRNA modifications are gaining great attention because of their major roles in RNA metabolism, immunity, development and disease. In this Review, we focus on the regulation and function of non-m6A modifications in eukaryotic mRNA, including pseudouridine (Ψ), N6,2'-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), inosine, 5-methylcytidine (m5C), N4-acetylcytidine (ac4C), 2'-O-methylated nucleotide (Nm) and internal N7-methylguanosine (m7G). We highlight their regulation, distribution, stoichiometry and known roles in mRNA metabolism, such as mRNA stability, translation, splicing and export. We also discuss their biological consequences in physiological and pathological processes. In addition, we cover research techniques to further study the non-m6A mRNA modifications and discuss their potential future applications.

Mettl1-mediated internal m7G methylation of Sptbn2 mRNA elicits neurogenesis and anti-alzheimer's disease

BACKGROUND

N7-methylguanosine (m7G) is one of the most conserved modifications in nucleosides impacting mRNA export, splicing, and translation. However, the precise function and molecular mechanism of internal mRNA m7G methylation in adult hippocampal neurogenesis and neurogenesis-related Alzheimer's disease (AD) remain unknown.

RESULTS

We profiled the dynamic Mettl1/Wdr4 expressions and m7G modification during neuronal differentiation of neural stem cells (NSCs) in vitro and in vivo. Adult hippocampal neurogenesis and its molecular mechanisms were examined by morphology, biochemical methods and biological sequencing. The translation efficiency of mRNA was detected by polysome profiling. The stability of Sptbn2 mRNA was constructed by RNA stability assay. APPswe/PS1ΔE9 (APP/PS1) double transgenic mice were used as model of AD. Morris water maze was used to detect the cognitive function.

METHODS

We found that m7G methyltransferase complex Mettl1/Wdr4 as well as m7G was significantly elevated in neurons. Functionally, silencing Mettl1 in neural stem cells (NSCs) markedly decreased m7G modification, neuronal genesis and proliferation in addition to increasing gliogenesis, while forced expression of Mettl1 facilitated neuronal differentiation and proliferation. Mechanistically, the m7G modification of Sptbn2 mRNA by Mettl1 enhanced its stability and translation, which promoted neurogenesis. Importantly, genetic defciency of Mettl1 reduced hippocampal neurogenesis and spatial memory in the adult mice. Furthermore, Mettl1 overexpression in the hippocampus of APP/PS1 mice rescued neurogenesis and behavioral defects.

CONCLUSION

Our findings unravel the pivotal role of internal mRNA m7G modification in Sptbn2-mediated neurogenesis, and highlight Mettl3 regulation of neurogenesis as a novel therapeutic target in AD treatment.

Fibroblast-specific knockout of METTL1 attenuates myocardial infarction-induced cardiac fibrosis

Cardiac fibrosis, a common pathology in inherited and acquired heart diseases, necessitates the identification of diagnostic and therapeutic targets. Methyltransferase Like 1 (METTL1), an enzyme responsible for RNA modification by methylating guanosine to form m7G, is an emerging area of research in understanding cellular processes and disease pathogenesis. Dysregulation of m7G modification has been implicated in various diseases. However, the role of METTL1 in cardiac fibrosis remains unclear. This study aimed to investigate the role of METTL1 in myocardial infarction-induced heart failure and cardiac fibrosis. Our findings demonstrate that elevated METTL1-mediated RNA m7G methylation is observed in cardiac fibrosis tissues and TGF-β1-induced cardiac fibroblast proliferation and myofibroblast transformation. Furthermore, fibroblast-specific knockout of METTL1 attenuated myocardial infarction-induced heart failure and cardiac fibrosis. Additionally, METTL1 knockout decreased m7G methylated fibrotic genes and impaired their translation efficiency. These results suggest a novel pro-fibrosis role of METTL1-mediated RNA m7G methylation, highlighting its potential as a therapeutic target in cardiac fibrosis.

Biological roles of RNA m7G modification and its implications in cancer

M7G modification, known as one of the common post-transcriptional modifications of RNA, is present in many different types of RNAs. With the accurate identification of m7G modifications within RNAs, their functional roles in the regulation of gene expression and different physiological functions have been revealed. In addition, there is growing evidence that m7G modifications are crucial in the emergence of cancer. Here, we review the most recent findings regarding the detection techniques, distribution, biological functions and Regulators of m7G. We also summarize the connections between m7G modifications and cancer development, drug resistance, and tumor microenvironment as well as we discuss the research's future directions and trends.

Base-resolution quantitative DAMM-seq for mapping RNA methylations in tRNA and mitochondrial polycistronic RNA

The human AlkB family proteins, such as FTO and ALKBH5, are known to mediate RNA m6A demethylation. However, although ALKBH7 localizes in mitochondria and affects metabolism, the detailed biological function and mechanism have remained unknown for years. We developed Demethylation-Assisted Multiple Methylation sequencing (DAMM-seq) to simultaneously detect N1-methyladenosine (m1A), N3-methylcytidine (m3C), N1-methylguanosine (m1G) and N2,N2-dimethylguanosine (m22G) methylations in both steady-state RNA and nascent RNA, and discovered that human ALKBH7 demethylates m22G and m1A within mt-Ile and mt-Leu1 pre-tRNA regions, respectively, in mitochondrial polycistronic RNA. DAMM-seq quantitatively and sensitively monitors the methylation stoichiometry change at pre-tRNA junctions within nascent mt-RNA, revealing the target region where ALKBH7 regulates RNA processing and local structural switch of polycistronic mt-RNAs. A new RNA demethylase in human cells was characterized through the base-resolution quantification of multiple RNA methylations in nascent mt-RNA, resolving the long-standing question about the functional substrate of ALKBH7.

Sequencing methods and functional decoding of mRNA modifications

More than 160 types of post-transcriptional RNA modifications have been reported; there is substantial variation in modification type, abundance, site, and function across species, tissues, and RNA type. The recent development of high-throughput detection technology has enabled identification of diverse dynamic and reversible RNA modifications, including N6,2'-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), 5-methylcytosine (m5C), N6-methyladenosine (m6A), pseudouridine (Ψ), and inosine (I). In this review, we focus on eukaryotic mRNA modifications. We summarize their biogenesis, regulatory mechanisms, and biological functions, as well as high-throughput methods for detection of mRNA modifications. We also discuss challenges that must be addressed in mRNA modification research.

Prognostic and immune implications of a novel 7-methylguanosine-related microRNA signature in breast invasive carcinoma: from exploration to validation

OBJECTIVES

This study aims to develop and validate a prognostic signature based on 7-methylguanosine-related (M7G-related) miRNAs for predicting prognosis and immune implications in breast invasive carcinoma (BRCA).

MATERIALS AND METHODS

M7G-related miRNA data of BRCA were obtained from The Cancer Genome Atlas (TCGA). Least absolute shrinkage and selection operator (LASSO)-penalized, univariate, and multivariate Cox regression analyses were used to construct the prognostic signature. Furthermore, the predictive validity was verified using Kaplan-Meier (KM) survival risk and receiver operating characteristic (ROC) plots. Internal random sampling verification was used to simplify and validate the signature. RT-qPCR was used to quantify the expression level of transcriptional profiles. The independent prognostic role of the risk score was validated using univariate and multivariate regression. Single-sample Gene Set Enrichment Analysis (ssGSEA) was used for functional and immune enrichment analysis.

RESULTS

A total of 18 M7G-related miRNAs were identified to construct the prognostic signature in BRCA. The low-risk group exhibited significantly higher overall survival than the high-risk group in the KM survival plot (P < 0.001). The area under the curve (AUC) for 1-, 3-, and 5-year survivals in the ROC curve were 0.737, 0.724, and 0.702, respectively. The survival significance in the training and testing cohorts was confirmed by random sampling verification. The most prominent miRNAs in the signature were the miR-7, miR-139, miR-10b, and miR-4728. Furthermore, immune scores for B, mast, and Th1 cells varied between risk groups. Our research demonstrated that CD52 was the most positively correlated gene with immune cells and functions in BRCA.

CONCLUSION

Our study presents a comprehensive and systematic analysis of M7G-related miRNAs to construct a prognostic signature in BRCA. The signature demonstrated excellent prognostic validity, with the risk score as an independent prognostic factor. These results provide critical evidence for further investigation of M7G miRNAs and offer new insights for BRCA patients in the context of effective immunotherapy.

RNA epigenetic modifications in ovarian cancer: The changes, chances, and challenges

Ovarian cancer (OC) is the most common female cancer worldwide. Patients with OC have high mortality because of its complex and poorly understood pathogenesis. RNA epigenetic modifications, such as m6 A, m1 A, and m5 C, are closely associated with the occurrence and development of OC. RNA modifications can affect the stability of mRNA transcripts, nuclear export of RNAs, translation efficiency, and decoding accuracy. However, there are few overviews that summarize the link between m6 A RNA modification and OC. Here, we discuss the molecular and cellular functions of different RNA modifications and how their regulation contributes to the pathogenesis of OC. By improving our understanding of the role of RNA modifications in the etiology of OC, we provide new perspectives for their use in OC diagnosis and treatment. This article is categorized under: RNA Processing > RNA Editing and Modification RNA in Disease and Development > RNA in Disease.

Regulation of the epigenome through RNA modifications

Chemical modifications of nucleotides expand the complexity and functional properties of genomes and transcriptomes. A handful of modifications in DNA bases are part of the epigenome, wherein DNA methylation regulates chromatin structure, transcription, and co-transcriptional RNA processing. In contrast, more than 150 chemical modifications of RNA constitute the epitranscriptome. Ribonucleoside modifications comprise a diverse repertoire of chemical groups, including methylation, acetylation, deamination, isomerization, and oxidation. Such RNA modifications regulate all steps of RNA metabolism, including folding, processing, stability, transport, translation, and RNA's intermolecular interactions. Initially thought to influence all aspects of the post-transcriptional regulation of gene expression exclusively, recent findings uncovered a crosstalk between the epitranscriptome and the epigenome. In other words, RNA modifications feedback to the epigenome to transcriptionally regulate gene expression. The epitranscriptome achieves this feat by directly or indirectly affecting chromatin structure and nuclear organization. This review highlights how chemical modifications in chromatin-associated RNAs (caRNAs) and messenger RNAs (mRNAs) encoding factors involved in transcription, chromatin structure, histone modifications, and nuclear organization affect gene expression transcriptionally.

TET2-mediated mRNA demethylation regulates leukemia stem cell homing and self-renewal

TET2 is recurrently mutated in acute myeloid leukemia (AML) and its deficiency promotes leukemogenesis (driven by aggressive oncogenic mutations) and enhances leukemia stem cell (LSC) self-renewal. However, the underlying cellular/molecular mechanisms have yet to be fully understood. Here, we show that Tet2 deficiency significantly facilitates leukemogenesis in various AML models (mediated by aggressive or less aggressive mutations) through promoting homing of LSCs into bone marrow (BM) niche to increase their self-renewal/proliferation. TET2 deficiency in AML blast cells increases expression of Tetraspanin 13 (TSPAN13) and thereby activates the CXCR4/CXCL12 signaling, leading to increased homing/migration of LSCs into BM niche. Mechanistically, TET2 deficiency results in the accumulation of methyl-5-cytosine (m5C) modification in TSPAN13 mRNA; YBX1 specifically recognizes the m5C modification and increases the stability and expression of TSPAN13 transcripts. Collectively, our studies reveal the functional importance of TET2 in leukemogenesis, leukemic blast cell migration/homing, and LSC self-renewal as an mRNA m5C demethylase.

Identification and validation of a m7G-related lncRNA signature for predicting the prognosis and therapy response in hepatocellular carcinoma

BACKGROUND

N7-methylguanosine (m7G) is one of the most common RNA posttranscriptional modifications; however, its potential role in hepatocellular carcinoma (HCC) remains unknown. We developed a prediction signature based on m7G-related long noncoding RNAs (lncRNAs) to predict HCC prognosis and provide a reference for immunotherapy and chemotherapy.

METHODS

RNA-seq data from The Cancer Genome Atlas (TCGA) database and relevant clinical data were used. Univariate and multivariate Cox regression analyses were conducted to identify m7G-related lncRNAs with prognostic value to build a predictive signature. We evaluated the prognostic value and clinical relevance of this signature and explored the correlation between the predictive signature and the chemotherapy treatment response of HCC. Moreover, an in vitro study to validate the function of CASC19 was performed.

RESULTS

Six m7G-related lncRNAs were identified to create a signature. This signature was considered an independent risk factor for the prognosis of patients with HCC. TIDE analyses showed that the high-risk group might be more sensitive to immunotherapy. ssGSEA indicated that the predictive signature was strongly related to the immune activities of HCC. HCC in high-risk patients was more sensitive to the common chemotherapy drugs bleomycin, doxorubicin, gemcitabine, and lenalidomide. In vitro knockdown of CASC19 inhibited the proliferation, migration and invasion of HCC cells.

CONCLUSION

We established a 6 m7G-related lncRNA signature that may assist in predicting the prognosis and response to chemotherapy and immunotherapy of HCC.

Comprehensive analysis of circRNAs for N7-methylguanosine methylation modification in human oral squamous cell carcinoma

N7-methylguanosine (m7G) modification is closely related to the occurrence of tumors. However, the m7G modification of circRNAs in oral squamous cell carcinoma (OSCC) remains to be investigated. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) was used to measure the methylation levels of m7G and identify m7G sites in circRNAs in human OSCC and normal tissues. The host genes of differentially methylated and differentially expressed circRNAs were analyzed by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, and circRNA-miRNA-mRNA networks were predicted using the miRanda and miRDB databases. The analysis identified 2348 m7G peaks in 624 circRNAs in OSCC tissues. In addition, the source of m7G-methylated circRNAs in OSCC was mainly the sense overlap region compared with normal tissues. The most conserved m7G motif in OSCC tissues was CCUGU, whereas the most conserved motif in normal tissues was RCCUG (R = G/A). Importantly, GO enrichment and KEGG pathway analysis showed that the host genes of differentially methylated and differentially expressed circRNAs were involved in many cellular biological functions. Furthermore, the significantly differentially expressed circRNAs were analyzed to predict the circRNA-miRNA-mRNA networks. This study revealed the whole profile of circRNAs of differential m7G methylation in OSCC and suggests that m7G-modified circRNAs may impact the development of OSCC.

RNA modifications in cancer

Currently, more than 170 modifications have been identified on RNA. Among these RNA modifications, various methylations account for two-thirds of total cases and exist on almost all RNAs. Roles of RNA modifications in cancer are garnering increasing interest. The research on m6A RNA methylation in cancer is in full swing at present. However, there are still many other popular RNA modifications involved in the regulation of gene expression post-transcriptionally besides m6A RNA methylation. In this review, we focus on several important RNA modifications including m1A, m5C, m7G, 2'-O-Me, Ψ and A-to-I editing in cancer, which will provide a new perspective on tumourigenesis by peeking into the complex regulatory network of epigenetic RNA modifications, transcript processing, and protein translation.

Advances in mRNA therapeutics for cancer immunotherapy: From modification to delivery

RNA vaccines have demonstrated their ability to solve the issues posed by the COVID-19 pandemic. This success has led to the renaissance of research into mRNA and their nanoformulations as potential therapeutic modalities for various diseases. The potential of mRNA as a template for synthesizing proteins and protein fragments for cancer immunotherapy is now being explored. Despite the promise, the use of mRNA in cancer immunotherapy is limited by challenges, such as low stability against extracellular RNases, poor delivery efficiency to the target organs and cells, short circulatory half-life, variable expression levels and duration. This review highlights recent advances in chemical modification and advanced delivery systems that are helping to address these challenges and unlock the biological and pharmacological potential of mRNA therapeutics in cancer immunotherapy. The review concludes by discussing future perspectives for mRNA-based cancer immunotherapy, which holds great promise as a next-generation therapeutic modality.

M7G methylated core genes (METTL1 and WDR4) and associated RNA risk signatures are associated with prognosis and immune escape in HCC

N7 methylguanosine (m7G) has a crucial role the development of hepatocellular carcinoma (HCC). This study aimed to investigate the impact of the m7G methylation core genes (METTL1 and WDR4) and associated RNA risk signatures on HCC. we found m7G methylation core genes (METTL1 and WDR4) were upregulated in four HCC cell lines, and downregulation of METTL1 and WDR4 attenuated HCC cell proliferation, migration, and invasion. Moreover, METTL1 and WDR4 are upregulated in HCC tissues, and that there is a significant positive correlation between them. METTL1 and WDR4 were identified as independent prognostic markers for HCC by employing overall survival (OS), disease-specific survival (DSS), Progression Free Interval survival (PFI), and univariate/multivariate Cox analyses. We identified 1479 coding RNAs (mRNAs) and 232 long non-coding RNAs (lncRNAs) associated with METTL1 / WDR4 by using weighted coexpression network analysis (WGCNA) and co-clustering analysis. The least absolute shrinkage and selection operator (lasso) were used to constructing mRNA and lncRNA risk signatures associated with the METTL1 / WDR4. These risk were independent poor prognostic factors in HCC. Furthermore, we found that METTL1 / WDR4 expression and mRNA / lncRNA risk scores were closely associated with TP53 mutations. Clinicopathological features correlation results showed that METTL1 / WDR4 expression and mRNA / lncRNA risk score were associated with the stage and invasion depth (T) of HCC. To predict the overall survival of HCC individuals, we constructed a nomogram with METTL1/WDR4 expression, mRNA/lncRNA risk score, and clinicopathological features. In addition, we combined single-cell sequencing datasets and immune escape-related checkpoints to construct an immune escape-related protein-protein interaction(PPI) network. In conclusion, M7G methylated core genes (METTL1 and WDR4) and associated RNA risk signatures are associated with prognosis and immune escape in HCC.

Dual chemical labeling enables nucleotide-resolution mapping of DNA abasic sites and common alkylation damage in human mitochondrial DNA

Mitochondrial DNA (mtDNA) modifications play an emerging role in innate immunity and inflammatory diseases. Nonetheless, relatively little is known regarding the locations of mtDNA modifications. Such information is critically important for deciphering their roles in mtDNA instability, mtDNA-mediated immune and inflammatory responses, and mitochondrial disorders. The affinity probe-based enrichment of lesion-containing DNA represents a key strategy for sequencing DNA modifications. Existing methods are limited in the enrichment specificity of abasic (AP) sites, a prevalent DNA modification and repair intermediate. Herein, we devise a novel approach, termed dual chemical labeling-assisted sequencing (DCL-seq), for mapping AP sites. DCL-seq features two designer compounds for enriching and mapping AP sites specifically at single-nucleotide resolution. For proof of principle, we mapped AP sites in mtDNA from HeLa cells under different biological conditions. The resulting AP site maps coincide with mtDNA regions with low TFAM (mitochondrial transcription factor A) coverage and with potential G-quadruplex-forming sequences. In addition, we demonstrated the broader applicability of the method in sequencing other DNA modifications in mtDNA, such as N7-methyl-2'-deoxyguanosine and N3-methyl-2'-deoxyadenosine, when coupled with a lesion-specific repair enzyme. Together, DCL-seq holds the promise to sequence multiple DNA modifications in various biological samples.

QKI shuttles internal m7G-modified transcripts into stress granules and modulates mRNA metabolism

N7-methylguanosine (m7G) modification, routinely occurring at mRNA 5' cap or within tRNAs/rRNAs, also exists internally in messenger RNAs (mRNAs). Although m7G-cap is essential for pre-mRNA processing and protein synthesis, the exact role of mRNA internal m7G modification remains elusive. Here, we report that mRNA internal m7G is selectively recognized by Quaking proteins (QKIs). By transcriptome-wide profiling/mapping of internal m7G methylome and QKI-binding sites, we identified more than 1,000 high-confidence m7G-modified and QKI-bound mRNA targets with a conserved "GANGAN (N = A/C/U/G)" motif. Strikingly, QKI7 interacts (via C terminus) with the stress granule (SG) core protein G3BP1 and shuttles internal m7G-modified transcripts into SGs to regulate mRNA stability and translation under stress conditions. Specifically, QKI7 attenuates the translation efficiency of essential genes in Hippo signaling pathways to sensitize cancer cells to chemotherapy. Collectively, we characterized QKIs as mRNA internal m7G-binding proteins that modulate target mRNA metabolism and cellular drug resistance.

Quantitative base-resolution sequencing technology for mapping pseudouridines in mammalian mRNA

Posttranscriptional RNA modifications occur in almost all types of RNA in all life forms. As an abundant RNA modification in mammals, pseudouridine (Ψ) regulates diverse biological functions of different RNA species such as ribosomal RNA (rRNA), transfer RNA (tRNA), small nuclear RNA (snRNA), etc. However, the functional investigation of mRNA pseudouridine (Ψ) has been hampered by the lack of a quantitative method that can efficiently map Ψ transcriptome-wide. We developed bisulfite-induced deletion sequencing (BID-seq), with the optimized bisulfite-based chemical reaction to convert pseudouridine selectively and completely into Ψ-BS adduct without cytosine deamination. The Ψ-BS adduct can be further read out as deletion signatures during reverse transcription. The deletion ratios induced by Ψ sites were used for estimating the modification stoichiometry at each modified site. BID-seq starts with 10-20 ng polyA+ RNA and detects thousands of mRNA Ψ sites with stoichiometry information in cell lines and tissues. We uncovered consensus motifs for Ψ in mammalian mRNA and assigned specific 'writer' proteins to individual Ψ deposition. BID-seq also confirmed the presence of Ψ within stop codons of mammalian mRNA. BID-seq set the stage for future investigations of Ψ functions in diverse biological processes.

Novel insight into RNA modifications in tumor immunity: Promising targets to prevent tumor immune escape

An immunosuppressive state is a typical feature of the tumor microenvironment. Despite the dramatic success of immune checkpoint inhibitor (ICI) therapy in preventing tumor cell escape from immune surveillance, primary and acquired resistance have limited its clinical use. Notably, recent clinical trials have shown that epigenetic drugs can significantly improve the outcome of ICI therapy in various cancers, indicating the importance of epigenetic modifications in immune regulation of tumors. Recently, RNA modifications (N⁶-methyladenosine [m⁶A], N¹-methyladenosine [m¹A], 5-methylcytosine [m⁵C], etc.), novel hotspot areas of epigenetic research, have been shown to play crucial roles in protumor and antitumor immunity. In this review, we provide a comprehensive understanding of how m⁶A, m¹A, and m⁵C function in tumor immunity by directly regulating different immune cells as well as indirectly regulating tumor cells through different mechanisms, including modulating the expression of immune checkpoints, inducing metabolic reprogramming, and affecting the secretion of immune-related factors. Finally, we discuss the current status of strategies targeting RNA modifications to prevent tumor immune escape, highlighting their potential.

TROP2 translation mediated by dual m6A/m7G RNA modifications promotes bladder cancer development

RNA modifications, including adenine methylation (m6A) of mRNA and guanine methylation (m7G) of tRNA, are crucial for the biological function of RNA. However, the mechanism underlying the translation of specific genes synergistically mediated by dual m6A/m7G RNA modifications in bladder cancer (BCa) remains unclear. We demonstrated that m6A methyltransferase METTL3-mediated programmable m6A modification of oncogene trophoblast cell surface protein 2 (TROP2) mRNA promoted its translation during malignant transformation of bladder epithelial cells. m7G methyltransferase METTL1 enhanced TROP2 translation by mediating m7G modification of certain tRNAs. TROP2 protein inhibition decreased the proliferation and invasion of BCa cells in vitro and in vivo. Moreover, synergistical knockout of METTL3/METTL1 inhibited BCa cell proliferation, migration, and invasion; however, TROP2 overexpression partially abrogated its effect. Furthermore, TROP2 expression was significantly positively correlated with the expression levels of METTL3 and METTL1 in BCa patients. Overall, our results revealed that METTL3/METTL1-mediated dual m6A/m7G RNA modifications enhanced TROP2 translation and promoted BCa development, indicating a novel RNA epigenetic mechanism in BCa.

Advances in targeting RNA modifications for anticancer therapy

Numerous strategies are employed by cancer cells to control gene expression and facilitate tumorigenesis. In the study of epitranscriptomics, a diverse set of modifications to RNA represent a new player of gene regulation in disease and in development. N6-methyladenosine (m6A) is the most common modification on mammalian messenger RNA and tends to be aberrantly placed in cancer. Recognized by a series of reader proteins that dictate the fate of the RNA, m6A-modified RNA could promote tumorigenesis by driving protumor gene expression signatures and altering the immunologic response to tumors. Preclinical evidence suggests m6A writer, reader, and eraser proteins are attractive therapeutic targets. First-in-human studies are currently testing small molecule inhibition against the methyltransferase-like 3 (METTL3)/methyltransferase-like 14 (METTL14) methyltransferase complex. Additional modifications to RNA are adopted by cancers to drive tumor development and are under investigation.

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.

TMEM11 regulates cardiomyocyte proliferation and cardiac repair via METTL1-mediated m7G methylation of ATF5 mRNA

The mitochondrial transmembrane (TMEM) protein family has several essential physiological functions. However, its roles in cardiomyocyte proliferation and cardiac regeneration remain unclear. Here, we detected that TMEM11 inhibits cardiomyocyte proliferation and cardiac regeneration in vitro. TMEM11 deletion enhanced cardiomyocyte proliferation and restored heart function after myocardial injury. In contrast, TMEM11-overexpression inhibited neonatal cardiomyocyte proliferation and regeneration in mouse hearts. TMEM11 directly interacted with METTL1 and enhanced m⁷G methylation of Atf5 mRNA, thereby increasing ATF5 expression. A TMEM11-dependent increase in ATF5 promoted the transcription of Inca1, an inhibitor of cyclin-dependent kinase interacting with cyclin A1, which suppressed cardiomyocyte proliferation. Hence, our findings revealed that TMEM11-mediated m⁷G methylation is involved in the regulation of cardiomyocyte proliferation, and targeting the TMEM11-METTL1-ATF5-INCA1 axis may serve as a novel therapeutic strategy for promoting cardiac repair and regeneration.

Fast and Efficient Design of Deep Neural Networks for Predicting N7-Methylguanosine Sites Using autoBioSeqpy

N⁷-Methylguanosine (m⁷G) is a crucial post-transcriptional RNA modification that plays a pivotal role in regulating gene expression. Accurately identifying m⁷G sites is a fundamental step in understanding the biological functions and regulatory mechanisms associated with this modification. While whole-genome sequencing is the gold standard for RNA modification site detection, it is a time-consuming, expensive, and intricate process. Recently, computational approaches, especially deep learning (DL) techniques, have gained popularity in achieving this objective. Convolutional neural networks and recurrent neural networks are examples of DL algorithms that have emerged as versatile tools for modeling biological sequence data. However, developing an efficient network architecture with superior performance remains a challenging task, requiring significant expertise, time, and effort. To address this, we previously introduced a tool called autoBioSeqpy, which streamlines the design and implementation of DL networks for biological sequence classification. In this study, we utilized autoBioSeqpy to develop, train, evaluate, and fine-tune sequence-level DL models for predicting m⁷G sites. We provided detailed descriptions of these models, along with a step-by-step guide on their execution. The same methodology can be applied to other systems dealing with similar biological questions. The benchmark data and code utilized in this study can be accessed for free at http://github.com/jingry/autoBioSeeqpy/tree/2.0/examples/m7G.

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.

Identification of m7G-Related miRNA Signatures Associated with Prognosis, Oxidative Stress, and Immune Landscape in Lung Adenocarcinoma

The role of N7-methylguanosine(m7G)-related miRNAs in lung adenocarcinoma (LUAD) remains unclear. We used LUAD data from The Cancer Genome Atlas (TCGA) to establish a risk model based on the m7G-related miRNAs, and divided patients into high-risk or low-risk subgroups. A nomogram for predicting overall survival (OS) was then constructed based on the independent risk factors. In addition, we performed a functional enrichment analysis and defined the oxidative stress-related genes, immune landscape as well as a drug response profile in the high-risk and low-risk subgroups. This study incorporated 28 m7G-related miRNAs into the risk model. The data showed a significant difference in the OS between the high-risk and low-risk subgroups. The receiver operating characteristic curve (ROC) predicted that the area under the curve (AUC) of one-year, three-year and five-year OS was 0.781, 0.804 and 0.853, respectively. The C-index of the prognostic nomogram for predicting OS was 0.739. We then analyzed the oxidative stress-related genes and immune landscape in the high-risk and low-risk subgroups. The data demonstrated significant differences in the expression of albumin (ALB), estimated score, immune score, stromal score, immune cell infiltration and functions between the high-risk and low-risk subgroups. In addition, the drug response analysis showed that low-risk subgroups may be more sensitive to tyrosine kinase inhibitor (TKI) and histone deacetylase (HDAC) inhibitors. We successfully developed a novel risk model based on m7G-related miRNAs in this study. The model can predict clinical prognosis and guide therapeutic regimens in patients with LUAD. Our data also provided new insights into the molecular mechanisms of m7G in LUAD.

Epigenetic regulation in the tumor microenvironment: molecular mechanisms and therapeutic targets

Over decades, researchers have focused on the epigenetic control of DNA-templated processes. Histone modification, DNA methylation, chromatin remodeling, RNA modification, and noncoding RNAs modulate many biological processes that are crucial to the development of cancers. Dysregulation of the epigenome drives aberrant transcriptional programs. A growing body of evidence suggests that the mechanisms of epigenetic modification are dysregulated in human cancers and might be excellent targets for tumor treatment. Epigenetics has also been shown to influence tumor immunogenicity and immune cells involved in antitumor responses. Thus, the development and application of epigenetic therapy and cancer immunotherapy and their combinations may have important implications for cancer treatment. Here, we present an up-to-date and thorough description of how epigenetic modifications in tumor cells influence immune cell responses in the tumor microenvironment (TME) and how epigenetics influence immune cells internally to modify the TME. Additionally, we highlight the therapeutic potential of targeting epigenetic regulators for cancer immunotherapy. Harnessing the complex interplay between epigenetics and cancer immunology to develop therapeutics that combine thereof is challenging but could yield significant benefits. The purpose of this review is to assist researchers in understanding how epigenetics impact immune responses in the TME, so that better cancer immunotherapies can be developed.