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

Methylated guanosine and uridine modifications in S. cerevisiae mRNAs modulate translation elongation

Chemical modifications to protein encoding messenger RNAs (mRNAs) influence their localization, translation, and stability within cells. Over 15 different types of mRNA modifications have been observed by sequencing and liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) approaches. While LC-MS/MS is arguably the most essential tool available for studying analogous protein post-translational modifications, the high-throughput discovery and quantitative characterization of mRNA modifications by LC-MS/MS has been hampered by the difficulty of obtaining sufficient quantities of pure mRNA and limited sensitivities for modified nucleosides. We have overcome these challenges by improving the mRNA purification and LC-MS/MS pipelines. The methodologies we developed result in no detectable non-coding RNA modifications signals in our purified mRNA samples, quantify 50 ribonucleosides in a single analysis, and provide the lowest limit of detection reported for ribonucleoside modification LC-MS/MS analyses. These advancements enabled the detection and quantification of 13 S. cerevisiae mRNA ribonucleoside modifications and reveal the presence of four new S. cerevisiae mRNA modifications at low to moderate levels (1-methyguanosine, N2-methylguanosine, N2,N2-dimethylguanosine, and 5-methyluridine). We identified four enzymes that incorporate these modifications into S. cerevisiae mRNAs (Trm10, Trm11, Trm1, and Trm2, respectively), though our results suggest that guanosine and uridine nucleobases are also non-enzymatically methylated at low levels. Regardless of whether they are incorporated in a programmed manner or as the result of RNA damage, we reasoned that the ribosome will encounter the modifications that we detect in cells. To evaluate this possibility, we used a reconstituted translation system to investigate the consequences of modifications on translation elongation. Our findings demonstrate that the introduction of 1-methyguanosine, N2-methylguanosine and 5-methyluridine into mRNA codons impedes amino acid addition in a position dependent manner. This work expands the repertoire of nucleoside modifications that the ribosome must decode in S. cerevisiae. Additionally, it highlights the challenge of predicting the effect of discrete modified mRNA sites on translation de novo because individual modifications influence translation differently depending on mRNA sequence context.

WDR4 gene polymorphisms and Wilms tumor susceptibility in Chinese children: A five-center case-control study

Wilms tumor is the most common embryonal renal malignancy in children. WDR4 is an indispensable noncatalytic subunit of the RNA N7-methylguanosine (m7G) methyltransferase complex and plays an essential role in tumorigenesis. However, the relationship between polymorphisms in the WDR4 gene and susceptibility to Wilms tumor remains to be fully investigated. We performed a large case-control study involving 414 patients and 1199 cancer-free controls to investigate whether single nucleotide polymorphisms (SNPs) in the WDR4 gene are associated with Wilms tumor susceptibility. WDR4 gene polymorphisms (rs2156315 C > T, rs2156316 C > G, rs6586250 C > T, rs15736 G > A, and rs2248490 C > G) were genotyped using the TaqMan assay. In addition, unconditioned logistic regression analysis was performed, odds ratios (ORs) and 95% confidence intervals (CIs) were used to assess the association between WDR4 gene SNPs and Wilms tumor susceptibility as well as the strength of the associations. We found that only the rs6586250 C>T polymorphism was significantly associated with an increased risk of Wilms tumor (adjusted OR=2.99, 95% CI = 1.28-6.97, P = 0.011 for the rs6586250 TT genotype; adjusted OR=3.08, 95% CI = 1.33-7.17, P = 0.009 for the rs6586250 CC/CT genotype). Furthermore, the stratification analysis revealed that patients with the rs6586250 TT genotype and carriers with 1-5 risk genotypes exhibited statistically significant associations with increased Wilms tumor risk in specific subgroups. However, the rs2156315 CT/TT genotype was identified as having a protective effect against Wilms tumor in the age >18 months subgroup compared with the rs2156315 CC genotype. In brief, our study demonstrated that the rs6586250 C > T polymorphism of the WDR4 gene was significantly associated with Wilms tumor. This finding may contribute to the understanding of the genetic mechanism of Wilms tumor.

RNA Modification in the Immune System

Characterization of RNA modifications has identified their distribution features and molecular functions. Dynamic changes in RNA modification on various forms of RNA are essential for the development and function of the immune system. In this review, we discuss the value of innovative RNA modification profiling technologies to uncover the function of these diverse, dynamic RNA modifications in various immune cells within healthy and diseased contexts. Further, we explore our current understanding of the mechanisms whereby aberrant RNA modifications modulate the immune milieu of the tumor microenvironment and point out outstanding research questions.

The diverse role of RNA methylation in esophageal cancer

Esophageal cancer is one of the major life-threatening diseases in the world. RNA methylation is the most common post-transcriptional modification and a wide-ranging regulatory system controlling gene expression. Numerous studies have revealed that dysregulation of RNA methylation is critical for cancer development and progression. However, the diverse role of RNA methylation and its regulators in esophageal cancer remains to be elucidated and summarized. In this review, we focus on the regulation of major RNA methylation, including m ⁶A, m ⁵C, and m ⁷G, as well as the expression patterns and clinical implications of its regulators in esophageal cancer. We systematically summarize how these RNA modifications affect the "life cycle" of target RNAs, including mRNA, microRNA, long non-coding RNA, and tRNA. The downstream signaling pathways associated with RNA methylation during the development and treatment of esophageal cancer are also discussed in detail. Further studies on how these modifications function together in the microenvironment of esophageal cancer will draw a clearer picture of the clinical application of novel and specific therapeutic strategies.

Emerging roles of tRNA in cancer

Transfer RNAs (tRNAs) play pivotal roles in the transmission of genetic information, and abnormality of tRNAs directly leads to translation disorders and causes diseases, including cancer. The complex modifications enable tRNA to execute its delicate biological function. Alteration of appropriate modifications may affect the stability of tRNA, impair its ability to carry amino acids, and disrupt the pairing between anticodons and codons. Studies confirmed that dysregulation of tRNA modifications plays an important role in carcinogenesis. Furthermore, when the stability of tRNA is impaired, tRNAs are cleaved into small tRNA fragments (tRFs) by specific RNases. Though tRFs have been found to play vital regulatory roles in tumorigenesis, its formation process is far from clear. Understanding improper tRNA modifications and abnormal formation of tRFs in cancer is conducive to uncovering the role of metabolic process of tRNA under pathological conditions, which may open up new avenues for cancer prevention and treatment.

An Update of Epigenetic Drugs for the Treatment of Cancers and Brain Diseases: A Comprehensive Review

Epigenetics has long been recognized as a significant field in biology and is defined as the investigation of any alteration in gene expression patterns that is not attributed to changes in the DNA sequences. Epigenetic marks, including histone modifications, non-coding RNAs, and DNA methylation, play crucial roles in gene regulation. Numerous studies in humans have been carried out on single-nucleotide resolution of DNA methylation, the CpG island, new histone modifications, and genome-wide nucleosome positioning. These studies indicate that epigenetic mutations and aberrant placement of these epigenetic marks play a critical role in causing the disease. Consequently, significant development has occurred in biomedical research in identifying epigenetic mechanisms, their interactions, and changes in health and disease conditions. The purpose of this review article is to provide comprehensive information about the different types of diseases caused by alterations in epigenetic factors such as DNA methylation and histone acetylation or methylation. Recent studies reported that epigenetics could influence the evolution of human cancer via aberrant methylation of gene promoter regions, which is associated with reduced gene function. Furthermore, DNA methyltransferases (DNMTs) in the DNA methylation process as well as histone acetyltransferases (HATs)/histone deacetylases (HDACs) and histone methyltransferases (HMTs)/demethylases (HDMs) in histone modifications play important roles both in the catalysis and inhibition of target gene transcription and in many other DNA processes such as repair, replication, and recombination. Dysfunction in these enzymes leads to epigenetic disorders and, as a result, various diseases such as cancers and brain diseases. Consequently, the knowledge of how to modify aberrant DNA methylation as well as aberrant histone acetylation or methylation via inhibitors by using epigenetic drugs can be a suitable therapeutic approach for a number of diseases. Using the synergistic effects of DNA methylation and histone modification inhibitors, it is hoped that many epigenetic defects will be treated in the future. Numerous studies have demonstrated a link between epigenetic marks and their effects on brain and cancer diseases. Designing appropriate drugs could provide novel strategies for the management of these diseases in the near future.

When N7-methyladenosine modification meets cancer: Emerging frontiers and promising therapeutic opportunities

N⁷-methylguanosine (m⁷G) methylation, one of the most common RNA modifications in eukaryotes, has recently gained considerable attention. The biological functions of m⁷G modification in RNAs, including tRNA, rRNA, mRNA, and miRNA, remain largely unknown in human diseases. Owing to rapid advances in high-throughput technologies, increasing evidence suggests that m⁷G modification plays a critical role in cancer initiation and progression. As m⁷G modification and hallmarks of cancer are inextricably linked together, targeting m⁷G regulators may provide new possibilities for future cancer diagnoses and potential intervention targets. This review summarizes various detection methods for m⁷G modification, recent advances in m⁷G modification and tumor biology regarding their interplay and regulatory mechanisms. We conclude with an outlook on the future of diagnosing and treating m⁷G-related diseases.

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

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

Emerging role of mRNA methylation in regulating the hallmarks of cancer

The dynamic chemical modifications of DNA, RNA, and proteins can transform normal cells into malignant ones. While the DNA and protein modifications in cancer have been described extensively in the literature, there are fewer reports about the role of RNA modifications in cancer. There are over 100 forms of RNA modifications and one of these, mRNA methylation, plays a critical role in the malignant properties of the cells. mRNA methylation is a reversible modification responsible for regulating protein expression at the post-transcriptional level. Despite being discovered in the 1970s, a complete understanding of the different proteins involved and the mechanism behind mRNA methylation remains largely unknown. However, these mRNA methylations have been shown to foster cancer hallmarks via specific cellular targets inside the cell. In this review, we provide a brief overview of mRNA methylation and its emerging role in regulating the various hallmarks of cancer.

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

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

Quantitative sequencing using BID-seq uncovers abundant pseudouridines in mammalian mRNA at base resolution

Functional characterization of pseudouridine (Ψ) in mammalian mRNA has been hampered by the lack of a quantitative method that maps Ψ in the whole transcriptome. We report bisulfite-induced deletion sequencing (BID-seq), which uses a bisulfite-mediated reaction to convert pseudouridine stoichiometrically into deletion upon reverse transcription without cytosine deamination. BID-seq enables detection of abundant Ψ sites with stoichiometry information in several human cell lines and 12 different mouse tissues using 10-20 ng input RNA. We uncover consensus sequences for Ψ in mammalian mRNA and assign different 'writer' proteins to individual Ψ deposition. Our results reveal a transcript stabilization role of Ψ sites installed by TRUB1 in human cancer cells. We also detect the presence of Ψ within stop codons of mammalian mRNA and confirm the role of Ψ in promoting stop codon readthrough in vivo. BID-seq will enable future investigations of the roles of Ψ in diverse biological processes.

Expression pattern and prognostic value of key regulators for N7-methylguanosine RNA modification in prostate cancer

Alterations in the regulators of RNA methylation modifications, such as N7-methylguanosine (m7G), have been implicated in a variety of diseases. Therefore, the analysis and identification of disease-related m7G modification regulators will accelerate advances in understanding disease pathogenesis. However, the implications of alterations in the regulators of m7G modifications remain poorly understood in prostate adenocarcinoma. In the present study, we analyze the expression patterns of 29 m7G RNA modification regulators in prostate adenocarcinoma using The Cancer Genome Atlas (TCGA) and perform consistent clustering analysis of differentially expressed genes (DEGs). We find that 18 m7G-related genes are differentially expressed in tumor and normal tissues. In different cluster subgroups, DEGs are mainly enriched in tumorigenesis and tumor development. Furthermore, immune analyses demonstrate that patients in cluster 1 have significantly higher scores for stromal and immune cells, such as B cells, T cells, and macrophages. Then, a TCGA-related risk model is developed and successfully validated using a Gene Expression Omnibus external dataset. Two genes ( EIF4A1 and NCBP2) are determined to be prognostically significant. Most importantly, we construct tissue microarrays from 26 tumor specimens and 20 normal specimens, and further confirm that EIF4A1 and NCBP2 are associated with tumor progression and Gleason score. Therefore, we conclude that the m7G RNA methylation regulators may be involved in the poor prognosis of patients with prostate adenocarcinoma. The results of this study may provide support for exploring the underlying molecular mechanisms of m7G regulators, especially EIF4A1 and NCBP2.

Determining RNA Natural Modifications and Nucleoside Analog-Labeled Sites by a Chemical/Enzyme-Induced Base Mutation Principle

The natural chemical modifications of messenger RNA (mRNA) in living organisms have shown essential roles in both physiology and pathology. The mapping of mRNA modifications is critical for interpreting their biological functions. In another dimension, the synthesized nucleoside analogs can enable chemical labeling of cellular mRNA through a metabolic pathway, which facilitates the study of RNA dynamics in a pulse-chase manner. In this regard, the sequencing tools for mapping both natural modifications and nucleoside tags on mRNA at single base resolution are highly necessary. In this work, we review the progress of chemical sequencing technology for determining both a variety of naturally occurring base modifications mainly on mRNA and a few on transfer RNA and metabolically incorporated artificial base analogs on mRNA, and further discuss the problems and prospects in the field.

m6A-SAC-seq for quantitative whole transcriptome m6A profiling

N6-methyladenosine (m6A) is the most abundant mRNA modification in mammalian cells, regulating many physiological processes. Here we describe a method for base-resolution, quantitative m6A sequencing in the whole transcriptome. The enzyme and small-molecule cofactor used in this protocol are prepared by recombinant protein expression and organic synthesis, respectively. Then the library can be prepared from various types of RNA samples using a ligation-based strategy, with m6A modifications being labeled by the enzyme and cofactor. Detailed instructions on ensuing data analysis are also included in this protocol. The method generates highly reproducible results, uncovering 31,233-129,263 sites using as little as 2 ng of poly A+ RNA. These identified sites correspond well with previous m6A profiling results, covering over 65% of peaks detected by the antibody-based approaches. Compared with other currently available methods, this method can be applied to various types of biological samples, including fresh and frozen tissues as well as formalin-fixed paraffin-embedded samples, providing a quantitative method to uncover new insights into m6A biology. The protocol requires basic expertise in molecular biology, recombinant protein expression and organic synthesis. The whole protocol can be done in 15 days, with the library preparation taking 5 days.

Association of RNA m7G Modification Gene Polymorphisms with Pediatric Glioma Risk

Glioma stemming from glial cells of the central nervous system (CNS) is one of the leading causes of cancer death in childhood. The genetic predisposition of glioma is not fully understood. METTL1-WDR4 methyltransferase complex is implicated in tumorigenesis by catalyzing N7-methylguanosine (m⁷G) modification of RNA. This study is aimed at determining the association of glioma risk with three polymorphisms (rs2291617, rs10877013, and rs10877012) in METTL1 and five polymorphisms (rs2156315 rs2156316, rs6586250, rs15736, and rs2248490) in WDR4 gene in children of Chinese Han. We enrolled 314 cases and 380 controls from three independent hospitals. Genotypes of these polymorphisms were determined using the TaqMan assay. We found the WDR4 gene rs15736 was significantly associated with reduced glioma risk (GA/AA vs. GG: adjusted odds ratio = 0.63, 95%confidence interval = 0.42 - 0.94, P = 0.023) out of the eight studied polymorphisms. Stratified analyses showed that the association of rs15736 with the risk of glioma remained significant in children aged 60 months or older, girls, the subgroups with astrocytic tumors, or grade I + II glioma. We also found the combined effects of five WDR4 gene polymorphisms on glioma risk. Finally, expression quantitative trait locus (eQTL) analyses elucidated that the rs15736 polymorphism was related to the expression level of WDR4 and neighboring gene cystathionine-beta-synthase (CBS). Our finding provided evidence of a causal association between WDR4 gene polymorphisms and glioma susceptibility in Chinese Han children.

BID-seq: The Quantitative and Base-Resolution Sequencing Method for RNA Pseudouridine

Quantitative and base-resolution sequencing methods are critical to investigations of the biological functions of diverse RNA modifications. These methods may also be employed for clinical studies and clinical applications in the future. In this In Focus article, we introduce and discuss the development of Bisulfite-Induced Deletion sequencing (BID-seq) for quantitatively detecting mRNA pseudouridine (Ψ) modifications at base resolution.

Internal m7G methylation: A novel epitranscriptomic contributor in brain development and diseases

In recent years, N7-methylguanosine (m7G) methylation, originally considered as messenger RNA (mRNA) 5' caps modifications, has been identified at defined internal positions within multiple types of RNAs, including transfer RNAs, ribosomal RNAs, miRNA, and mRNAs. Scientists have put substantial efforts to discover m7G methyltransferases and methylated sites in RNAs to unveil the essential roles of m7G modifications in the regulation of gene expression and determine the association of m7G dysregulation in various diseases, including neurological disorders. Here, we review recent findings regarding the distribution, abundance, biogenesis, modifiers, and functions of m7G modifications. We also provide an up-to-date summary of m7G detection and profile mapping techniques, databases for validated and predicted m7G RNA sites, and web servers for m7G methylation prediction. Furthermore, we discuss the pathological roles of METTL1/WDR-driven m7G methylation in neurological disorders. Last, we outline a roadmap for future directions and trends of m7G modification research, particularly in the central nervous system.

RMDisease V2.0: an updated database of genetic variants that affect RNA modifications with disease and trait implication

Recent advances in epitranscriptomics have unveiled functional associations between RNA modifications (RMs) and multiple human diseases, but distinguishing the functional or disease-related single nucleotide variants (SNVs) from the majority of 'silent' variants remains a major challenge. We previously developed the RMDisease database for unveiling the association between genetic variants and RMs concerning human disease pathogenesis. In this work, we present RMDisease v2.0, an updated database with expanded coverage. Using deep learning models and from 873 819 experimentally validated RM sites, we identified a total of 1 366 252 RM-associated variants that may affect (add or remove an RM site) 16 different types of RNA modifications (m6A, m5C, m1A, m5U, Ψ, m6Am, m7G, A-to-I, ac4C, Am, Cm, Um, Gm, hm5C, D and f5C) in 20 organisms (human, mouse, rat, zebrafish, maize, fruit fly, yeast, fission yeast, Arabidopsis, rice, chicken, goat, sheep, pig, cow, rhesus monkey, tomato, chimpanzee, green monkey and SARS-CoV-2). Among them, 14 749 disease- and 2441 trait-associated genetic variants may function via the perturbation of epitranscriptomic markers. RMDisease v2.0 should serve as a useful resource for studying the genetic drivers of phenotypes that lie within the epitranscriptome layer circuitry, and is freely accessible at: www.rnamd.org/rmdisease2.

Recent Advances of RNA m6A Modifications in Cancer Immunoediting and Immunotherapy

Cancer immunotherapy, which modulates immune responses against tumors using immune-checkpoint inhibitors or adoptive cell transfer, has emerged as a novel and promising therapy for tumors. However, only a minority of patients demonstrate durable responses, while the majority of patients are resistant to immunotherapy. The immune system can paradoxically constrain and promote tumor development and progression. This process is referred to as cancer immunoediting. The mechanisms of resistance to immunotherapy seem to be that cancer cells undergo immunoediting to evade recognition and elimination by the immune system. RNA modifications, specifically N6-methyladenosine (m6A) methylation, have emerged as a key regulator of various post-transcriptional gene regulatory processes, such as RNA export, splicing, stability, and degradation, which play unappreciated roles in various physiological and pathological processes, including immune system development and cancer pathogenesis. Therefore, a deeper understanding of the mechanisms by which RNA modifications impact the cancer immunoediting process can provide insight into the mechanisms of resistance to immunotherapies and the strategies that can be used to overcome such resistance. In this chapter, we briefly introduce the background of cancer immunoediting and immunotherapy. We also review and discuss the roles and mechanisms of RNA m6A modifications in fine-tuning the innate and adaptive immune responses, as well as in regulating tumor escape from immunosurveillance. Finally, we summarize the current strategies targeting m6A regulators for cancer immunotherapy.

N6-methyladenosine functions and its role in skin cancer

N6-methyladenosine (m6A) methylation is the most abundant mammalian mRNA modification. m6A regulates RNA processing, splicing, nucleation, translation and stability by transferring, removing and recognizing m6A methylation sites, which are critical for cancer initiation, progression, metabolism and metastasis. m6A is involved in pathophysiological tumour development by altering m6A modification and expression levels in tumour oncogenes and suppressor genes. Skin cancers are by far the most common malignancies in humans, with well over a million cases diagnosed each year. Skin cancers are grouped into two main categories: melanoma and non-melanoma skin cancers (NMSC), based on cell origin and clinical behaviour. In this review, we summarize m6A methylation functions in different skin cancers, and discuss how m6A methylation is involved in disease development and progression. Moreover, we review potential prognostic biomarkers and molecular targets for early skin cancer diagnosis and treatment.

tRNA Modifications and Modifying Enzymes in Disease, the Potential Therapeutic Targets

tRNA is one of the most conserved and abundant RNA species, which plays a key role during protein translation. tRNA molecules are post-transcriptionally modified by tRNA modifying enzymes. Since high-throughput sequencing technology has developed rapidly, tRNA modification types have been discovered in many research fields. In tRNA, numerous types of tRNA modifications and modifying enzymes have been implicated in biological functions and human diseases. In our review, we talk about the relevant biological functions of tRNA modifications, including tRNA stability, protein translation, cell cycle, oxidative stress, and immunity. We also explore how tRNA modifications contribute to the progression of human diseases. Based on previous studies, we discuss some emerging techniques for assessing tRNA modifications to aid in discovering different types of tRNA modifications.

m7G-quant-seq: Quantitative Detection of RNA Internal N7-Methylguanosine

Methods for the precise detection and quantification of RNA modifications are critical to uncover functional roles of diverse RNA modifications. The internal m⁷G modification in mammalian cytoplasmic tRNAs is known to affect tRNA function and impact embryonic stem cell self-renewal, tumorigenesis, cancer progression, and other cellular processes. Here, we introduce m⁷G-quant-seq, a quantitative method that accurately detects internal m⁷G sites in human cytoplasmic tRNAs at single-base resolution. The efficient chemical reduction and mild depurination can almost completely convert internal m⁷G sites into RNA abasic sites (AP sites). We demonstrate that RNA abasic sites induce a mixed variation pattern during reverse transcription, including G → A or C or T mutations as well as deletions. We calculated the total variation ratio to quantify the m⁷G modification fraction at each methylated site. The calibration curves of all relevant motif contexts allow us to more quantitatively determine the m⁷G methylation level. We detected internal m⁷G sites in 22 human cytoplasmic tRNAs from HeLa and HEK293T cells and successfully estimated the corresponding m⁷G methylation stoichiometry. m⁷G-quant-seq could be applied to monitor the tRNA m⁷G methylation level change in diverse biological processes.

The role of post-transcriptional modifications during development

While the existence of post-transcriptional modifications of RNA nucleotides has been known for decades, in most RNA species the exact positions of these modifications and their physiological function have been elusive until recently. Technological advances, such as high-throughput next-generation sequencing (NGS) methods and nanopore-based mapping technologies, have made it possible to map the position of these modifications with single nucleotide accuracy, and genetic screens have uncovered the “writer”, “reader” and “eraser” proteins that help to install, interpret and remove such modifications, respectively. These discoveries led to intensive research programmes with the aim of uncovering the roles of these modifications during diverse biological processes. In this review, we assess novel discoveries related to the role of post-transcriptional modifications during animal development, highlighting how these discoveries can affect multiple aspects of development from fertilization to differentiation in many species.

Analysis of m7G methylation modification patterns and pulmonary vascular immune microenvironment in pulmonary arterial hypertension

Background

M⁷G methylation modification plays an important role in cardiovascular disease development. Dysregulation of the immune microenvironment is closely related to the pathogenesis of PAH. However, it is unclear whether m⁷G methylation is involved in the progress of PAH by affecting the immune microenvironment.

Methods

The gene expression profile of PAH was obtained from the GEO database, and the m⁷G regulatory factors were analyzed for differences. Machine learning algorithms were used to screen characteristic genes, including the least absolute shrinkage and selection operator, random forest, and support vector machine recursive feature elimination analysis. Constructed a nomogram model, and receiver operating characteristic was used to evaluate the diagnosis of disease characteristic genes value. Next, we used an unsupervised clustering method to perform consistent clustering analysis on m⁷G differential genes. Used the ssGSEA algorithm to estimate the relationship between the m⁷G regulator in PAH and immune cell infiltration and analyze the correlation with disease-characteristic genes. Finally, the listed drugs were evaluated through the screened signature genes.

Results

We identified 15 kinds of m⁷G differential genes. CYFIP1, EIF4E, and IFIT5 were identified as signature genes by the machine learning algorithm. Meanwhile, two m⁷G molecular subtypes were identified by consensus clustering (cluster A/B). In addition, immune cell infiltration analysis showed that activated CD4 T cells, regulatory T cells, and type 2 T helper cells were upregulated in m⁷G cluster B, CD56 dim natural killer cells, MDSC, and monocyte were upregulated in the m⁷G cluster A. It might be helpful to select Calpain inhibitor I and Everolimus for the treatment of PAH.

Conclusion

Our study identified CYFIP1, EIF4E, and IFIT5 as novel diagnostic biomarkers in PAH. Furthermore, their association with immune cell infiltration may facilitate the development of immune therapy in PAH.

Constructing and validating of m7G-related genes prognostic signature for hepatocellular carcinoma and immune infiltration: potential biomarkers for predicting the overall survival

Background

To investigate the prognostic significance of N7-methylguanosine (m7G) regulators and immune infiltration in liver hepatocellular carcinoma (LIHC).

Methods

The research measured predictive m7G genes in LIHC samples from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) datasets. Data on the stemness index based on mRNA expression (mRNAsi), gene mutations, and corresponding clinical characteristics were obtained from TCGA and ICGC. Lasso regression was used to construct the prediction model to assess the m7G prognostic signals in LIHC. Based on these genes, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to identify key biological functions and pathways. The correlation between m7G RNA methylation regulators and the prognosis and immune infiltration of LIHC was evaluated.

Results

There were 21 m7G-related differentially expressed genes (DEGs) in LIHC and healthy tissues, and LIHC patients could be divided into two categories by consensus clustering of these DEGs. A five-gene predictive approach was employed using least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Patients in the low-risk group showed a significantly higher survival rate compared with those in the high-risk group (P=0.001). Validations using the ICGC database. Also, univariate and multivariate Cox regression analyses suggested that the risk score produced by the predictive model is an independent predictor for LIHC [hazard ratio (HR): 1.848, 95% confidence interval (CI): 1.286-2.656; HR: 2.597, 95% CI: 1.358-4.965]. The ROC curves of the ICGC cohort revealed that the five-gene prediction model performed well [area under the curve (AUC) =0.642 at 1 year, AUC =0.686 at 2 years, and AUC =0.667 at 3 years]. Immuno-oncology scoring revealed that in the high-risk group, among 16 immune cells, the expressions of neutrophils and natural killer (NK) cells were low and that of regulatory T-cells (Tregs) was high.

Conclusions

LIHC occurrence and progression are linked to m7G-related genes. Corresponding prognostic models help forecast the prognosis of LIHC patients. m7G-related genes and associated immune cell infiltration in the TME may serve as potential therapeutic targets in LIHC, which requires further trials. In addition, the m7G-related gene signature offers a viable alternative to predict LIHC, and these m7G-related genes show a prospective research area for LIHC targeted treatment in the future.

METTL1 drives tumor progression of bladder cancer via degrading ATF3 mRNA in an m7G-modified miR-760-dependent manner

7-methylguanosine (m⁷G) modification is recently found to conservatively exist in RNA internal position besides mRNA caps and mediates the various RNA metabolisms. As the core confirmed transmethylase of m⁷G modification, METTL1 has been reported in certain human cancers. However, the role of internal m⁷G at miRNAs and its core writer METTL1 in bladder cancer (BCa) remains to be elucidated. Here, we demonstrated that METTL1 was indispensable for BCa proliferation and metastasis in vitro and in vivo. By combining miRNA sequencing, m⁷G methylated RNA immunoprecipitation (MeRIP) and RIP, we identified METTL1 promoted the processing of miR-760 in an m⁷G-dependent manner. Transcription sequencing suggested that METTL1 indirectly degrades tumor suppressor ATF3 mRNA mediated by miR-760. Together, we concluded a regulatory axis composed of METTL1/m⁷G/miR-760/ATF3 in regulating BCa progression and provided potential therapeutic targets for BCa.

Epitranscriptome: Review of Top 25 Most-Studied RNA Modifications

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

M7G-Related lncRNAs predict prognosis and regulate the immune microenvironment in lung squamous cell carcinoma

Background

N7-Methylguanosine (m7G) and long non-coding RNAs (lncRNAs) have been widely studied in cancer and have been found to be useful for assessing tumor progression. However, the role of m7G-related lncRNAs in lung squamous cell carcinoma (LUSC) remains unclear. Thus, it is crucial to identify m7G-associated lncRNAs with definitive prognostic value. This study aimed to investigate the prognostic value, correlation with tumor mutation burden, and impact on the tumor immune microenvironment of m7G-related lncRNAs in LUSC. 

Methods

LUSC transcriptome data and clinical data were downloaded from The Cancer Genome Atlas, and an m7G-related lncRNA-mRNA co-expression network was constructed using Pearson’s correlation analysis. Cox regression analyses were used to determine a risk model for m7G-associated lncRNAs with prognostic value. The risk signature was verified using the Kaplan–Meier method, receiver operating characteristic curve analysis, and principal component analysis. A nomogram based on risk scores and clinical characteristics was then developed. Gene set enrichment analysis was used for functional annotation to analyze the risk signature. The association among the risk signature, tumor mutational burden, and tumor-infiltrating immune cells was then analyzed. RT-qPCR was used to investigate the expression of 6 m7G-related lncRNAs in LUSC cells. The cytological function of SRP14-AS1 was verified by wound-healing assay and transwell assay.

Results

A total of 293 m7G-related lncRNAs were identifed, 27 candidate m7G-related lncRNAs were signifcantly associated with overall survival (OS). Six of these lncRNAs (CYP4F26P, LINC02178, MIR22HG, SRP14-AS1, TMEM99, PTCSC2) were selected for establishment of the risk model. The OS of patients in the low-risk group was higher than that of patients in the high-risk group (p < 0.001). Multivariate cox regression analysis indicated that the model could be an independent prognostic factor for LUSC (HR = 1.859; 95% CI 1.452–2.380, p < 0.001). The ROC curve analysis revealed that the AUCs for OS in the 3-, and 5-year were 0.682, 0.657, respectively. GSEA analysis revealed that the risk model was closely related to immune-related pathways. Compared with normal lung epithelial cells, four m7G-related lncRNAs were higher expressed in cancer cells and two were lower expressed, among which knockdown of SRP14-AS1 promoted the proliferation and migration of LUSC cells.

Conclusion

A risk model based on six m7G-related lncRNAs with prognostic value may be a promising prognostic tool in LUSC and guide individualized patient treatment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-10232-z.

BRPCA: Bounded Robust Principal Component Analysis to Incorporate Similarity Network for N7-Methylguanosine(m7G) Site-Disease Association Prediction

Recent studies have revealed that N7-methylguanosine(m⁷G) plays a pivotal role in various biological processes and disease pathogenesis. To date, transcriptome-wide m⁷G modification sites have been identified by high-throughput sequencing approaches, and some related information has been recorded in a few biological databases. However, the mechanism of site action in disease remains uncharted. Wet experiments can help identify true m⁷G sites with high confidence, but it is time-consuming to find the true ones in such a large number of sites, which will also cost too much. Thus, computational methods are emergently needed to predict the associations between m⁷G sites and various diseases, thus help to uncover potential active sites for specific diseases. In this article, we proposed a bounded robust principal component analysis (BRPCA) method to predict unknown m⁷G-disease association based on similarity information. Importantly, BRPCA tolerates the noise and redundancy existing in association and similarity information. Moreover, a suitable bounded constraint is incorporated into BRPCA to ensure that the predicted association scores locate in a meaningful interval. The extensive experiments demonstrate the superiority and robustness of the BRPCA.

RNA binding proteins in MLL-rearranged leukemia

RNA binding proteins (RBPs) have recently emerged as important post-transcriptional gene expression regulators in both normal development and disease. RBPs influence the fate of mRNAs through multiple mechanisms of action such as RNA modifications, alternative splicing, and miR-mediated regulation. This complex and, often, combinatorial regulation by RBPs critically impacts the expression of oncogenic transcripts and, thus, the activation of pathways that drive oncogenesis. Here, we focus on the major features of RBPs, their mechanisms of action, and discuss the current progress in investigating the function of important RBPs in MLL-rearranged leukemia.

The m7G Modification Level and Immune Infiltration Characteristics in Patients with COVID-19

Purpose

The 7-methylguanosine (m7G)-related genes were used to identify the clinical severity and prognosis of patients with coronavirus disease 2019 (COVID-19) and to identify possible therapeutic targets.

Patients and Methods

The GSE157103 dataset provides the transcriptional spectrum and clinical information required to analyze the expression of m7G-related genes and the disease subtypes. R language was applied for immune infiltration analysis, functional enrichment analysis, and nomogram model construction.

Results

Most m7G-related genes were up-regulated in COVID-19 and were closely related to immune cell infiltration. Disease subtypes were grouped using a clustering algorithm. It was found that the m7G-cluster B was associated with higher immune infiltration, lower mechanical ventilation, lower intensive care unit (ICU) status, higher ventilator-free days, and lower m7G scores. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that differentially expressed genes (DEGs) between m7G-cluster A and B were enriched in viral infection and immune-related aspects, including COVID-19 infection; Th17, Th1, and Th2 cell differentiation, and human T-cell leukemia virus 1 infection. Finally, through machine learning, six disease characteristic genes, NUDT4B, IFIT5, LARP1, EIF4E, LSM1, and NUDT4, were screened and used to develop a nomogram model to estimate disease risk.

Conclusion

The expression of most m7G genes was higher in COVID-19 patients compared with that in non-COVID-19 patients. The m7G-cluster B showed higher immune infiltration and milder symptoms. The predictive nomogram based on the six m7G genes can be used to accurately assess risk.

A novel m7G methylation–related signature associated with chromosome homeostasis in patients with lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a malignant tumor of the respiratory system with poor prognosis. Recent studies have revealed that N7-methylguanosine (m7G) methylation is a widespread modification occurring in RNA. But the expression of m7G methylation–related genes in LUAD and their correlations with prognosis are still unclear. In this study, we found 12 m7G methylation–related regulators with differential expression between LUAD and normal lung tissues. According to differentially expressed genes (DEGs), all LUAD cases were separated into two subtypes. The prognostic value of each m7G methylation–related gene for survival was evaluated to construct a multigene signature using The Cancer Genome Atlas (TCGA) cohort. Finally, an m7G methylation–related prognostic signature based on three genes was built to classify LUAD patients into two risk groups. Patients in the high-risk group showed significantly reduced overall survival (OS) when compared with patients in the low-risk group (p &lt; 0.05). The receiver operating characteristic (ROC) curve analysis confirmed the predictive capacity of the signature. The Gene Ontology (GO) functional annotation analysis disclosed that chromosome homeostasis plays an important role in this process. The gene set enrichment analysis (ssGSEA) implied that the immune status was decreased in the high-risk group. To sum up, m7G methylation–related genes play a vital role in tumor immunity and the related signature is a reliable predictor for LUAD prognosis.

Comprehensive pan-cancer analysis of N7-methylguanosine regulators: Expression features and potential implications in prognosis and immunotherapy

Although immunotherapy has made great strides in cancer therapy, its effectiveness varies widely among individual patients as well as tumor types, and there is an urgent need to develop biomarkers for effectively assessing immunotherapy response. In recent years, RNA methylation regulators have demonstrated to be novel potential biomarkers for prognosis as well as immunotherapy of cancers, such as N6-methyladenine (m6A) and 5-methylcytosine (m5C). N7-methylguanosine (m7G) is a prevalent RNA modification in eukaryotes, but the relationship between m7G regulators and prognosis as well as tumor immune microenvironment is still unclear. In this study, a pan-cancer analysis of 26 m7G regulators across 17 cancer types was conducted based on the bioinformatics approach. On the one hand, a comprehensive analysis of expression features, genetic variations and epigenetic regulation of m7G regulators was carried out, and we found that the expression tendency of m7G regulators were different among tumors and their aberrant expression in cancers could be affected by single nucleotide variation (SNV), copy number variation (CNV), DNA methylation and microRNA (miRNA) separately or simultaneously. On the other hand, the m7Gscore was modeled based on single sample gene set enrichment analysis (ssGSEA) for evaluating the relationships between m7G regulators and cancer clinical features, hallmark pathways, tumor immune microenvironment, immunotherapy response as well as pharmacotherapy sensitivity, and we illustrated that the m7Gscore exhibited tight correlations with prognosis, several immune features, immunotherapy response and drug sensitivity in most cancers. In conclusion, our pan-cancer analysis revealed that m7G regulators may exert critical roles in the tumor progression and immune microenvironment, and have the potential as biomarkers for predicting prognosis, immunotherapy response as well as candidate drug compounds for cancer patients.

Detection technologies for RNA modifications

To date, more than 170 chemical modifications have been characterized in RNA, providing a new layer of gene expression regulation termed the 'epitranscriptome'. RNA modification detection methods and tools advance the functional studies of the epitranscriptome. According to the detection throughput and principles, existing RNA modification detection technologies can be categorized into four classes, including quantification methods, locus-specific detection methods, next-generation sequencing-based detection technologies and nanopore direct RNA sequencing-based technologies. In this review, we summarize the current knowledge about these RNA modification detection technologies and discuss the challenges for the existing detection tools, providing information for a comprehensive understanding of the epitranscriptome.

Role of N7-methylguanosine (m7G) in cancer

N7-methylguanosine (m7G) is a significant RNA modification occurring during epigenetic regulation. An increasing number of investigations have proved that the effect of m7G in suppressing cancer deserves more attention. Herein, we explore potential valuable targets based on present studies of m7G-related cancer to improve disease therapy and prognoses.

Comprehensive analysis of m7G modification patterns based on potential m7G regulators and tumor microenvironment infiltration characterization in lung adenocarcinoma

Background: The non-negligible role of epigenetic modifications in cancer development and tumor microenvironment (TME) has been demonstrated in recent studies. Nonetheless, the potential regulatory role of N7-methylguanosine (m⁷G) modification in shaping and impacting the TME remains unclear.

Methods: A comprehensive analysis was performed to explore the m⁷G modification patterns based on 24 potential m⁷G regulators in 817 lung adenocarcinoma (LUAD) patients, and the TME landscape in distinct m⁷G modification patterns were evaluated. The m⁷G score was established based on principal component analysis (PCA) to quantify m⁷G modification patterns and evaluate the TME cell infiltrating characteristics of individual tumors. Further, correlation analyses of m7Gscore with response to chemotherapy and immunotherapy were performed.

Results: We identified three distinct m⁷G modification patterns with the biological pathway enrichment and TME cell infiltrating characteristics corresponded to immune-desert, immune-inflamed and immune-excluded phenotype, respectively. We further demonstrated the m⁷Gscore could predict the TME infiltrating characteristics, tumor mutation burden (TMB), response to immunotherapy and chemotherapy, as well as prognosis of individual tumors. High m⁷Gscore was associated with increased component of immune cell infiltration, low TMB and survival advantage, while low m⁷Gscore was linked to decreased immune cell infiltration and increased TMB. Additionally, patients with lower m⁷Gscore demonstrated significant therapeutic advantages.

Conclusion: This study demonstrated the regulatory mechanisms of m⁷G modification on TME formation and regulation of lung adenocarcinoma. Identification of individual tumor m⁷G modification patterns will contribute to the understanding of TME characterization and guiding more effective immunotherapy strategies.

Integrative pan-cancer analysis and clinical characterization of the N7-methylguanosine (m7G) RNA modification regulators in human cancers

Background: RNA modification is one of the epigenetic mechanisms that regulates post-transcriptional gene expression, and abnormal RNA modifications have been reported to play important roles in tumorigenesis. N7-methylguanosine (m7G) is an essential modification at the 5′ cap of human mRNA. However, a systematic and pan-cancer analysis of the clinical relevance of m7G related regulatory genes is still lacking.

Methods: We used univariate Cox model and Kaplan-Meier analysis to generate the forest plot of OS, PFI, DSS and identified the correlation between the altered expression of m7G regulators and patient survival in 33 cancer types from the TCGA and GTEx databases. Then, the “estimate” R-package, ssGSEA and CIBERSORT were used to depict the pan-cancer immune landscape. Through Spearman’s correlation test, we analyzed the correlation between m7G regulators and the tumor microenvironment (TME), immune subtype, and drug sensitivity of the tumors, which was further validated in NSCLC. We also assessed the changes in the expression of m7G related regulatory genes in NSCLC with regards to the genetic and transcriptional aspects and evaluated the correlation of METTL1 and WDR4 expression with TMB, MSI and immunotherapy in pan-cancer.

Results: High expression of most of the m7G regulators was significantly associated with worse prognosis. Correlation analyses revealed that the expression of majority of the m7G regulators was correlated with tumor immune infiltration and tumor stem cell scores. Drug sensitivity analysis showed that the expression of CYFP1,2 was closely related to drug sensitivity for various anticancer agents (p &lt; 0.001). Analysis of the pan-cancer immune subtype revealed significant differences in the expression of m7G regulators between different immune subtypes (p &lt; 0.001). Additionally, the types and proportions of mutations in METTL1 and WDR4 and their relevance to immunotherapy were further described.

Conclusion: Our study is the first to evaluate the correlation between the altered expression of m7G regulators and patient survival, the degree of immune infiltration, TME and drug sensitivity in pan-cancer datasets.

RNA modifications: importance in immune cell biology and related diseases

RNA modifications have become hot topics recently. By influencing RNA processes, including generation, transportation, function, and metabolization, they act as critical regulators of cell biology. The immune cell abnormality in human diseases is also a research focus and progressing rapidly these years. Studies have demonstrated that RNA modifications participate in the multiple biological processes of immune cells, including development, differentiation, activation, migration, and polarization, thereby modulating the immune responses and are involved in some immune related diseases. In this review, we present existing knowledge of the biological functions and underlying mechanisms of RNA modifications, including N⁶-methyladenosine (m⁶A), 5-methylcytosine (m⁵C), N¹-methyladenosine (m¹A), N⁷-methylguanosine (m⁷G), N⁴-acetylcytosine (ac⁴C), pseudouridine (Ψ), uridylation, and adenosine-to-inosine (A-to-I) RNA editing, and summarize their critical roles in immune cell biology. Via regulating the biological processes of immune cells, RNA modifications can participate in the pathogenesis of immune related diseases, such as cancers, infection, inflammatory and autoimmune diseases. We further highlight the challenges and future directions based on the existing knowledge. All in all, this review will provide helpful knowledge as well as novel ideas for the researchers in this area.

Novel roles of METTL1/WDR4 in tumor via m7G methylation

As one of the prevalent posttranscriptional modifications of RNA, N⁷-methylguanosine (m⁷G) plays essential roles in RNA processing, metabolism, and function, mainly regulated by the methyltransferase-like 1 (METTL1) and WD repeat domain 4 (WDR4) complex. Emerging evidence suggests that the METTL1/WDR4 complex promoted or inhibited the processes of many tumors, including head and neck, lung, liver, colon, bladder cancer, and teratoma, dependent on close m⁷G methylation modification of tRNA or microRNA (miRNA). Therefore, METTL1 and m⁷G modification can be used as biomarkers or potential intervention targets, providing new possibilities for early diagnosis and treatment of tumors. This review will mainly focus on the mechanisms of METTL1/WDR4 via m⁷G in tumorigenesis and the corresponding detection methods.

METTL1 promotes neuroblastoma development through m7G tRNA modification and selective oncogenic gene translation

Background

Neuroblastoma (NBL) is the most common extra-cranial solid tumour in childhood, with prognosis ranging from spontaneous remission to high risk for rapid and fatal progression. Despite existing therapy approaches, the 5-year event-free survival (EFS) for patients with advanced NBL remains below 30%, emphasizing urgent necessary for novel therapeutic strategies. Studies have shown that epigenetic disorders play an essential role in the pathogenesis of NBL. However, the function and mechanism of N7-methylguanosine (m⁷G) methyltransferase in NBL remains unknown.

Methods

The expression levels of m⁷G tRNA methyltransferase Methyltransferase-like 1 (METTL1) were analyzed by querying the Gene Expression Omnibus (GEO) database and further confirmed by immunohistochemistry (IHC) assay. Kaplan-Meier, univariate and multivariate cox hazard analysis were performed to reveal the prognostic role of METTL1. Cell function assays were performed to evaluate how METTL1 works in proliferation, apoptosis and migration in cell lines and xenograft mouse models. The role of METTL1 on mRNA translation activity of NBL cells was measured using puromycin intake assay and polysome profiling assay. The m⁷G modified tRNAs were identified by tRNA reduction and cleavage sequencing (TRAC-seq). Ribosome nascent-chain complex-bound mRNA sequencing (RNC-seq) was utilized to identify the variation of gene translation efficiency (TE). Analyzed the codon frequency decoded by m⁷G tRNA to clarify the translation regulation and mechanism of m⁷G modification in NBL.

Results

This study found that METTL1 were significantly up-regulated in advanced NBL, which acted as an independent risk factor and predicted poor prognosis. Further in NBL cell lines and BALB/c-nu female mice, we found METTL1 played a crucial role in promoting NBL progression. Furthermore, m⁷G profiling and translation analysis revealed downregulation of METTL1 would inhibit puromycin intake efficiency of NBL cells, indicating that METTL1 did count crucially in regulation of NBL cell translation. With all tRNAs with m⁷G modification identified in NBL cells, knockdown of METTL1 would significantly reduce the levels of both m⁷G modification and m⁷G tRNAs expressions. Result of RNC-seq shew there were 339 overlapped genes with impaired translation in NBL cells upon METTL1 knockdown. Further analysis revealed these genes contained higher frequency of codons decoded by m⁷G-modified tRNAs and were enriched in oncogenic pathways.

Conclusion

This study revealed the critical role and mechanism of METTL1-mediated tRNA m⁷G modification in regulating NBL progression, providing new insights for developing therapeutic approaches for NBL patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40364-022-00414-z.

The role of RNA modification in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a highly mortal type of primary liver cancer. Abnormal epigenetic modifications are present in HCC, and RNA modification is dynamic and reversible and is a key post-transcriptional regulator. With the in-depth study of post-transcriptional modifications, RNA modifications are aberrantly expressed in human cancers. Moreover, the regulators of RNA modifications can be used as potential targets for cancer therapy. In RNA modifications, N6-methyladenosine (m6A), N7-methylguanosine (m7G), and 5-methylcytosine (m5C) and their regulators have important regulatory roles in HCC progression and represent potential novel biomarkers for the confirmation of diagnosis and treatment of HCC. This review focuses on RNA modifications in HCC and the roles and mechanisms of m6A, m7G, m5C, N1-methyladenosine (m1A), N3-methylcytosine (m3C), and pseudouridine (ψ) on its development and maintenance. The potential therapeutic strategies of RNA modifications are elaborated for HCC.

Identification of N7-methylguanosine-related IncRNA signature as a potential predictive biomarker for colon adenocarcinoma

N⁷-Methylguanosine (m7G) is an RNA modification serving as a key part of colon cancer development. Thus, a comprehensive analysis was executed to explore prognostic roles and associations with the immune status of the m7G-related lncRNA (m7G-RNAs) in colon adenocarcinoma (COAD). Identification of m7G-RNAs was achieved via Pearson’s correlation analysis of lncRNAs in the TCGA-COAD dataset and m7G regulators. A prognostic signature was developed via LASSO analyses. ESTIMATE, CIBERSORT, and ssGSEA algorithms were utilized to assess immune infiltration between different risk groups. Survival analysis suggested the high-risk group possesses poor outcomes compared with the low-risk group. According to the ROC curves, the m7G-RNAs signature exhibited a reliable capability of prediction (AUCs at 1, 3, and 5 years were 0.770, 0.766, and 0.849, respectively). Multivariate hazard analysis proved that the signature was an independent predictive indicator for OS. Moreover, the risk score was related to infiltration levels of naïve B cells, CD4⁺ memory T cells, and resting NK cells. The result revealed the prognostic value of m7G modification in COAD and provided a novel perspective on personalized immunotherapy strategies.

Identification of N7-methylguanosine related signature for prognosis and immunotherapy efficacy prediction in lung adenocarcinoma

Background

Lung adenocarcinoma (LUAD) is one of the most frequent causes of tumor-related mortality worldwide. Recently, the role of N7-methylguanosine (m7G) in tumors has begun to receive attention, but no investigation on the impact of m7G on LUAD. This study aims to elucidate the significance of m7G on the prognosis and immunotherapy in LUAD.

Methods

Consensus clustering was employed to determine the molecular subtype according to m7G-related regulators extracted from The Cancer Genome Atlas (TCGA) database. Survival, clinicopathological features and tumor mutational burden (TMB) analysis were applied to research molecular characteristics of each subtype. Subsequently, “limma” package was used to screen differentially expressed genes (DEGs) between subtypes. In the TCGA train cohort (n = 245), a prognostic signature was established by univariate Cox regression, lasso regression and multivariate Cox regression analysis according to DEGs and survival analysis was employed to assess the prognosis. Then the prognostic value of the signature was verified by TCGA test cohort (n = 245), TCGA entire cohort (n = 490) and GSE31210 cohort (n = 226). Moreover, the association among immune infiltration, clinical features and the signature was investigated. The immune checkpoints, TMB and tumor immune dysfunction and exclusion (TIDE) were applied to predict the immunotherapy response.

Results

Two novel molecular subtypes (C1 and C2) of LUAD were identified. Compared to C2 subtype, C1 subtype had poorer prognosis and higher TMB. Subsequently, the signature (called the “m7G score”) was constructed according to four key genes (E2F7, FAM83A, PITX3, and HOXA13). The distribution of m7G score were significantly different between two molecular subtypes. The patients with lower m7G score had better prognosis in TCGA train cohort and three verification cohort. The m7G score was intensively related to immune infiltration. Compared with the lower score, the higher m7G score was related to remarkable upregulation of the PD-1 and PD-L1, the higher TMB and the lower TIDE score.

Conclusion

This study established a m7G-related signature for predicting prognosis and immunotherapy in LUAD, which may contribute to the development of new therapeutic strategies for LUAD.

METTL1 gene polymorphisms synergistically confer hepatoblastoma susceptibility

INTRODUCTION

Hepatoblastoma is a rare but devastating pediatric liver malignancy. Overexpressed methyltransferase-like 1 (METTL1) is a methyltransferase that catalyzes essential N7-methylguanosine (m7G) modification of eukaryotic mRNA. Accumulating evidence has revealed the oncogenic potential of METTL1. However, whether METTL1 gene polymorphisms confer susceptibility to hepatoblastoma has not been reported. This study aimed to identify causal relationships between genetic variants of this gene and susceptibility to hepatoblastoma.

MATERIALS AND METHODS

Using the TaqMan assay, we genotyped three METTL1 polymorphisms (rs2291617 G > T, rs10877013 T > C, rs10877012 T > G) in germline DNA samples from 1759 Chinese children of Han ethnicity (313 cases vs. 1446 controls).

RESULTS

None of these polymorphisms were associated with hepatoblastoma risk. However, combination analysis showed that children with 1 to 3 risk genotypes were associated with increased hepatoblastoma risk (adjusted odds ratio = 1.47, 95% confidence interval  1.07-2.02; P = 0.018). Stratified analyses revealed significant effects of combined polymorphisms mainly among young children (< 17 months of age), boys, and those with advanced hepatoblastoma.

CONCLUSION

We identified some potential functional METTL1 gene polymorphisms that work together to increase the risk of hepatoblastoma among Chinese Han children; single polymorphism showed only weak effects. These METTL1 polymorphisms may be promising biomarkers for screening high-risk individuals for hepatoblastoma. These findings are inspiring and deserve to be validated among individuals of different ethnicities.

A novel prognostic 7-methylguanosine signature reflects immune microenvironment and alternative splicing in glioma based on multi-omics analysis

Glioma is the most common type of central nervous system tumor with increasing incidence. 7-methylguanosine (m7G) is one of the diverse RNA modifications that is known to regulate RNA metabolism and its dysregulation was associated with various cancers. However, the expression pattern of m7G regulators and their roles in regulating tumor immune microenvironments (TIMEs) as well as alternative splicing events (ASEs) in glioma has not been reported. In this study, we showed that m7G regulators displayed a close correlation with each other and most of them were differentially expressed between normal and glioma tissues. Two m7G signatures were then constructed to predict the overall survival of both GBM and LGG patients with moderate predictive performance. The risk score calculated from the regression coefficient and expression level of signature genes was proved to be an independent prognostic factor for patients with LGG, thus, a nomogram was established on the risk score and other independent clinical parameters to predict the survival probability of LGG patients. We also investigated the correlation of m7G signatures with TIMEs in terms of immune scores, expression levels of HLA and immune checkpoint genes, immune cell composition, and immune-related functions. While exploring the correlation between signature genes and the ASEs in glioma, we found that EIF4E1B was a key regulator and might play dual roles depending on glioma grade. By incorporating spatial transcriptomic data, we found a cluster of cells featured by high expression of PTN exhibited the highest m7G score and may communicate with adjacent cancer cells via SPP1 and PTN signaling pathways. In conclusion, our work brought novel insights into the roles of m7G modification in TIMEs and ASEs in glioma, suggesting that evaluation of m7G in glioma could predict prognosis. Moreover, our data suggested that blocking SPP1 and PTN pathways might be a strategy for combating glioma.

RNA Modifications in Gastrointestinal Cancer: Current Status and Future Perspectives

Gastrointestinal (GI) cancer, referring to cancers of the digestive system such as colorectal cancer (CRC), gastric cancer (GC), and liver cancer, is a major cause of cancer-related deaths in the world. A series of genetic, epigenetic, and epitranscriptomic changes occur during the development of GI cancer. The identification of these molecular events provides potential diagnostic, prognostic, and therapeutic targets for cancer patients. RNA modification is required in the posttranscriptional regulation of RNA metabolism, including splicing, intracellular transport, degradation, and translation. RNA modifications such as N6-methyladenosine (m6A) and N1-methyladenosine (m1A) are dynamically regulated by three different types of regulators named methyltransferases (writers), RNA binding proteins (readers), and demethylases (erasers). Recent studies have pointed out that abnormal RNA modification contributes to GI tumorigenesis and progression. In this review, we summarize the latest findings on the functional significance of RNA modification in GI cancer and discuss the therapeutic potential of epitranscriptomic inhibitors for cancer treatment.

Prognostic value of 12 m7G methylation-related miRNA markers and their correlation with immune infiltration in breast cancer

RNA guanine-7 methyltransferase (RNMT), in complex with FAM103A1, plays an important role in tumorigenesis and development. The aim of this study was to establish a prognostic model of RNMT and FAM103A1-based upstream microRNAs and explore its correlation with immune cell infiltration in breast cancer (BC) while investigating its potential prognostic value and verify the model by quantitative real-time polymerase chain reaction (qRT-PCR). The miRNA expression data upstream of the m7G methyltransferase complex RNMT/FAM103A1 in BC was obtained from The Cancer Genome Atlas and TargetScan databases. We performed univariate Cox regression, LASSO regression, Kaplan-Meier survival, and principal component analyses, along with risk prognostic modelling. Based on multivariate Cox regression analysis, a total of 12 m7G methyltransferase-related miRNAs were found. The model showed good accuracy for predicting the 1-, 3-,5-, and 10-year survival rates, and the areas under the curve were almost &gt;0.7. To characterize the risk-level model constructed from 12 miRNAs, 12 differentially expressed mRNAs related to prognosis and immune infiltration were obtained. The prognosis of BC patients is well predicted by the risk model we constructed. This model is also closely related to immune infiltration, and new immunotherapy targets can be explored from this field.

The Epitranscriptome in miRNAs: Crosstalk, Detection, and Function in Cancer

The epitranscriptome encompasses all post-transcriptional modifications that occur on RNAs. These modifications can alter the function and regulation of their RNA targets, which, if dysregulated, result in various diseases and cancers. As with other RNAs, miRNAs are highly modified by epitranscriptomic modifications such as m6A methylation, 2'-O-methylation, m5C methylation, m7G methylation, polyuridine, and A-to-I editing. miRNAs are a class of small non-coding RNAs that regulates gene expression at the post-transcriptional level. miRNAs have gathered high clinical interest due to their role in disease, development, and cancer progression. Epitranscriptomic modifications alter the targeting, regulation, and biogenesis of miRNAs, increasing the complexity of miRNA regulation. In addition, emerging studies have revealed crosstalk between these modifications. In this review, we will summarize the epitranscriptomic modifications-focusing on those relevant to miRNAs-examine the recent crosstalk between these modifications, and give a perspective on how this crosstalk expands the complexity of miRNA biology.