Emerging Roles of m6A RNA Methylation Regulators in Gynecological Cancer

Synchronous profiling of mRNA N6-methyladenosine modifications and mRNA expression in high-grade serous ovarian cancer: a pilot study

This study aimed to synchronously determine epitranscriptome-wide RNA N6-methyladenosine (m6A) modifications and mRNA expression profile in high grade serous ovarian cancer (HGSOC). The methylated RNA immunoprecipitation sequencing (MeRIP-seq) was used to comprehensively examine the m6A modification profile and the RNA-sequencing (RNA-seq) was performed to analyze the mRNA expression profile in HGSOC and normal fallopian tube (FT) tissues. Go and KEGG analyses were carried out in the enrichment of those differentially methylated and expressed genes. MeRIP-seq data showed 53,794 m6A methylated peaks related to 19,938 genes in the HGSOC group and 51,818 m6A peaks representing 19,681 genes in the FT group. RNA-seq results revealed 2321 upregulated and 2486 downregulated genes in HGSOC. Conjoint analysis of MeRIP-seq and RNA-seq data identified differentially expressed genes in which 659 were hypermethylated (330 up- and 329 down-regulated) and 897 were hypomethylated (475 up- and 422 down-regulated). Functional enrichment analysis indicated that these differentially modulated genes are involved in pathways related to cancer development. Among methylation regulators, the m6A eraser (FTO) expression was significantly lower, but the m6A readers (IGF2BP2 and IGF2BP3) were higher in HGSOC, which was validated by the subsequent real-time PCR assay. Exploration through public databases further corroborated their possible clinical application of certain methylation regulators and differentially expressed genes. For the first time, our study screens the epitranscriptome-wide m6A modification and expression profiles of their modulated genes and signaling pathways in HGSOC. Our findings provide an alternative direction in exploring the molecular mechanisms of ovarian pathogenesis and potential biomarkers in the diagnosis and predicting the prognosis of the disease.

RBM15 Knockdown Impairs the Malignancy of Cervical Cancer by Mediating m6A Modification of Decorin

Cervical cancer (CC) is considered to be the most prevalent female malignancies across the globe and a prime cause of mortality among women. RNA-binding motif protein 15 (RBM15) has been elucidated to participate in tumorigenesis in various cancers by regulating RNA N6-methyladenosine (m6A) methylation. However, its significance and detailed molecular mechanisms remain uncertain in CC. Using CGA database and qRT-PCR, the RBM15 expression was found to be elevated in CC tissues. After performing EdU, wound healing, Transwell migration, and xenograft tumor assays, RBM15 knockdown inhibited the malignant properties of CC cells along with the tumor development of CC cells in vivo. Moreover, qRT-PCR, MeRIP, and western blotting experiments were also confirmed that decorin (DCN) downregulated in CC was a direct substrate of RBM15 m6A methylation, and RBM15 knockdown could enhance DCN expression in CC cells. The anti-tumor effects of RBM15 knockdown could be abolished by DCN silencing. Overall, RBM15 knockdown lowered the tumorigenesis of CC both in vitro and in vivo, and it does so via mediating m6A modification of DCN mRNA in CC cells.

New understandings of the genetic regulatory relationship between non-coding RNAs and m6A modification

One of the most frequent epigenetic modifications of RNA in eukaryotes is N6 methyladenosine (m6A), which is mostly present in messenger RNAs. Through the influence of several RNA processing stages, m6A modification is a crucial approach for controlling gene expression, especially in cancer progression. It is universally acknowledged that numerous non-coding RNAs (ncRNAs), such as microRNAs, circular RNAs, long non-coding RNAs, and piRNAs, are also significantly affected by m6A modification, and the complex genetic regulatory relationship between m6A and ncRNAs plays a pivotal role in the development of cancer. The connection between m6A modifications and ncRNAs offers an opportunity to explore the oncogene potential regulatory mechanisms and suggests that m6A modifications and ncRNAs could be vital biomarkers for multiple cancers. In this review, we discuss the mechanisms of interaction between m6A methylation and ncRNAs in cancer, and we also summarize diagnostic and prognostic biomarkers for clinical cancer detection. Furthermore, our article includes some methodologies for identifying m6A sites when assessing biomarker potential.

Characterization of m6A modifiers and RNA modifications in uterine fibroids

Uterine leiomyoma or fibroids are the most common prevalent noncancerous tumors of the uterine muscle layer. Common symptoms associated with fibroids include pelvic pain, heavy menstrual bleeding, anemia, and pelvic pressure. These tumors are a leading cause of gynecological care but lack long-term therapy as the origin and development of fibroids are not well understood. Several next-generation sequencing technologies have been performed to identify the underlying genetic and epigenetic basis of fibroids. However, there remains a systemic gap in our understanding of molecular and biological process that define uterine fibroids. Recent epitranscriptomics studies have unraveled RNA modifications that are associated with all forms of RNA and are thought to influence both normal physiological functions and the progression of diseases. We quantified RNA expression profiles by analyzing publicly available RNA-seq data for 15 known epigenetic mediators to identify their expression profile in uterine fibroids compared to myometrium. To validate our findings, we performed RT-qPCR on a separate cohort of uterine fibroids targeting these modifiers confirming our RNA-seq data. We then examined protein profiles of key m6A modifiers in fibroids and their matched myometrium. In concordance with our RNA expression profiles, no significant differences were observed in these proteins in uterine fibroids compared to myometrium. To determine abundance of RNA modifications, mRNA and small RNA from fibroids and matched myometrium were analyzed by UHPLC MS/MS. In addition to the prevalent N6-methyladenosine (m6A), we identified 11 other known modifiers but did not identify any aberrant expression in fibroids. We then mined a previously published dataset and identified differential expression of m6A modifiers that were specific to fibroid genetic sub-type. Our analysis also identified m6A consensus motifs on genes previously identified to be dysregulated in uterine fibroids. Overall, using state-of-the-art mass spectrometry, RNA expression and protein profiles, we characterized and identified differentially expressed m6A modifiers in relation to driver mutations. Despite the use of several different approaches, we identified limited differential expression of RNA modifiers and associated modifications in uterine fibroids. However, considering the highly heterogenous genomic and cellular nature of fibroids, and the possible contribution of single molecule m6A modifications to fibroid pathology, there is a need for greater in-depth characterization of m6A marks and modifiers in a larger and varied patient cohort.

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.

Breast cancer prognosis and immunological characteristics are predicted using the m6A/m5C/m1A/m7G-related long noncoding RNA signature

According to statistics, breast cancer (BC) has replaced lung cancer as the most common cancer in the world. Therefore, specific detection markers and therapeutic targets need to be explored as a way to improve the survival rate of BC patients. We first identified m6A/m5C/m1A/m7G-related long noncoding RNAs (MRlncRNAs) and developed a model of 16 MRlncRNAs. Kaplan-Meier survival analysis was applied to assess the prognostic power of the model, while univariate Cox analysis and multivariate Cox analysis were used to assess the prognostic value of the constructed model. Then, we constructed a nomogram to illustrate whether the predicted results were in good agreement with the actual outcomes. We tried to use the model to distinguish the difference in sensitivity to immunotherapy between the two groups and performed some analyses such as immune infiltration analysis, ssGSEA and IC50 prediction. To explore the novel anti-tumor drug response, we reclassified the patients into two clusters. Next, we assessed their response to clinical treatment by the R package pRRophetic, which is determined by the IC50 of each BC patient. We finally identified 11 MRlncRNAs and based on them, a risk model was constructed. In this model, we found good agreement between calibration plots and prognosis prediction. The AUC of ROC curves was 0.751, 0.734, and 0.769 for 1-year, 2-year, and 3-year overall survival (OS), respectively. The results showed that the IC50 was significantly different between the risk groups, suggesting that the risk groups can be used as a guide for systemic treatment. We regrouped patients into two clusters based on 11 MRlncRNAs expression. Next, we conducted immune scores for 2 clusters, which showed that cluster 1 had higher stromal scores, immune scores and higher estimated (microenvironment) scores, demonstrating that TME of cluster 1 was different from cluster 2. The results of this study support that MRlncRNAs can predict tumor prognosis and help differentiate patients with different sensitivities to immunotherapy as a basis for individualized treatment for BC patients.

The role of RNA methyltransferase METTL3 in gynecologic cancers: Results and mechanisms

N6-methyladenosine (m6A) methylation is the most prevalent mRNA modification in eukaryotes, and it is defined as the methylation of nitrogen atoms on the six adenine (A) bases of RNA in the presence of methyltransferases. Methyltransferase-like 3 (Mettl3), one of the components of m6A methyltransferase, plays a decisive catalytic role in m6A methylation. Recent studies have confirmed that m6A is associated with a wide spectrum of biological processes and it significantly affects disease progression and prognosis of patients with gynecologic tumors, in which the role of Mettl3 cannot be ignored. Mettl3 is involved in numerous pathophysiological functions, such as embryonic development, fat accumulation, and tumor progression. Moreover, Mettl3 may serve as a potential target for treating gynecologic malignancies, thus, it may benefit the patients and prolong survival. However, there is a need to further study the role and mechanism of Mettl3 in gynecologic malignancies. This paper reviews the recent progression on Mettl3 in gynecologic malignancies, hoping to provide a reference for further research.

Comprehensive analysis of prognosis and immune function of CD70-CD27 signaling axis in pan-cancer

The immune checkpoint molecule CD70 and its receptor CD27 constitute the signal transduction axis, which is abnormally expressed in many solid tumors and is crucial for T cell co-stimulation and immune escape. Tumor cells regulate CD27 expression in the tumor microenvironment by expressing CD70, which promotes immune escape. Although current research evidence suggests a link between CD70 and tumors, no pan-cancer analysis is available. Using the Cancer Genome Atlas, Gene Expression Omnibus datasets, and online databases, we first explored the potential carcinogenic role of the CD70-CD27 signaling axis in human malignancies. Furthermore, qRT-PCR, Western blot, immunohistochemistry, and a T cell-mediated tumor cell killing assay were used to assess the biological function of the CD70-CD27 signaling axis. CD70 expression is upregulated in most cancers and has an obvious correlation with the prognosis of tumor patients. The expression of CD70 and CD27 is associated with the level of regulatory T cell (Treg) infiltration. In addition, T cell receptor signaling pathways, PI3K-AKT, NF-κB, and TNF signaling pathways are also involved in CD70-mediated immune escape. CD70 mainly regulates tumor immune escape by regulating T cell-mediated tumor killing, with Tregs possibly being its primary T cell subset. Our first pan-cancer study provides a relatively comprehensive understanding of the carcinogenic role of the CD70-CD27 signaling axis in different tumors.

Role of N6-Methyladenosine Methylation Regulators in the Drug Therapy of Digestive System Tumours

Digestive system tumours, including stomach, colon, esophagus, liver and pancreatic tumours, are serious diseases affecting human health. Although surgical treatment and postoperative chemoradiotherapy effectively improve patient survival, current diagnostic and therapeutic strategies for digestive system tumours lack sensitivity and specificity. Moreover, the tumour's tolerance to drug therapy is enhanced owing to tumour cell heterogeneity. Thus, primary or acquired treatment resistance is currently the main hindrance to chemotherapy efficiency. N6-methyladenosine (m6A) has various biological functions in RNA modification. m6A modification, a key regulator of transcription expression, regulates RNA metabolism and biological processes through the interaction of m6A methyltransferase ("writers") and demethylase ("erasers") with the binding protein decoding m6A methylation ("readers"). Additionally, m6A modification regulates the occurrence and development of tumours and is a potential driving factor of tumour drug resistance. This review systematically summarises the regulatory mechanisms of m6A modification in the drug therapy of digestive system malignancies. Furthermore, it clarifies the related mechanisms and therapeutic prospects of m6A modification in the resistence of digestive system malignancies to drug therapy.

The Role of m6A RNA Methylation in Cancer: Implication for Nature Products Anti-Cancer Research

N6-methyladenosine (m6A) RNA methylation is identified as the most common, abundant and reversible RNA epigenetic modification in messenger RNA (mRNA) and non-coding RNA, especially within eukaryotic messenger RNAs (mRNAs), which post-transcriptionally directs many important processes of RNA. It has also been demonstrated that m6A modification plays a pivotal role in the occurrence and development of tumors by regulating RNA splicing, localization, translation, stabilization and decay. Growing number of studies have indicated that natural products have outstanding anti-cancer effects of their unique advantages of high efficiency and minimal side effects. However, at present, there are very few research articles to study and explore the relationship between natural products and m6A RNA modification in tumorigenesis. m6A is dynamically deposited, removed, and recognized by m6A methyltransferases (METTL3/14, METTL16, WTAP, RBM15/15B, VIRMA, CBLL1, and ZC3H13, called as “writers”), demethylases (FTO and ALKBH5, called as “erasers”), and m6A-specific binding proteins (YTHDF1/2/3, YTHDC1/2, IGH2BP1/2/3, hnRNPs, eIF3, and FMR1, called as “readers”), respectively. In this review, we summarize the biological function of m6A modification, the role of m6A and the related signaling pathway in cancer, such as AKT, NF-kB, MAPK, ERK, Wnt/β-catenin, STAT, p53, Notch signaling pathway, and so on. Furthermore, we reviewed the current research on nature products in anti-tumor, and further to get a better understanding of the anti-tumor mechanism, thus provide an implication for nature products with anti-cancer research by regulating m6A modification in the future.

The Methylation Game: Epigenetic and Epitranscriptomic Dynamics of 5-Methylcytosine

DNA and RNA methylation dynamics have been linked to a variety of cellular processes such as development, differentiation, and the maintenance of genome integrity. The correct deposition and removal of methylated cytosine and its oxidized analogues is pivotal for cellular homeostasis, rapid responses to exogenous stimuli, and regulated gene expression. Uncoordinated expression of DNA/RNA methyltransferases and demethylase enzymes has been linked to genome instability and consequently to cancer progression. Furthermore, accumulating evidence indicates that post-transcriptional DNA/RNA modifications are important features in DNA/RNA function, regulating the timely recruitment of modification-specific reader proteins. Understanding the biological processes that lead to tumorigenesis or somatic reprogramming has attracted a lot of attention from the scientific community. This work has revealed extensive crosstalk between epigenetic and epitranscriptomic pathways, adding a new layer of complexity to our understanding of cellular programming and responses to environmental cues. One of the key modifications, m5C, has been identified as a contributor to regulation of the DNA damage response (DDR). However, the various mechanisms of dynamic m5C deposition and removal, and the role m5C plays within the cell, remains to be fully understood.

An m6A/m5C/m1A/m7G-Related Long Non-coding RNA Signature to Predict Prognosis and Immune Features of Glioma

Background: Gliomas are the most common and fatal malignant type of tumor of the central nervous system. RNA post-transcriptional modifications, as a frontier and hotspot in the field of epigenetics, have attracted increased attention in recent years. Among such modifications, methylation is most abundant, and encompasses N6-methyladenosine (m6A), 5-methylcytosine (m5C), N1 methyladenosine (m1A), and 7-methylguanosine (m7G) methylation.

Methods: RNA-sequencing data from healthy tissue and low-grade glioma samples were downloaded from of The Cancer Genome Atlas database along with clinical information and mutation data from glioblastoma tumor samples. Forty-nine m6A/m5C/m1A/m7G-related genes were identified and an m6A/m5C/m1A/m7G-lncRNA signature of co-expressed long non-coding RNAs selected. Least absolute shrinkage and selection operator Cox regression analysis was used to identify 12 m6A/m5C/m1A/m7G-related lncRNAs associated with the prognostic characteristics of glioma and their correlation with immune function and drug sensitivity analyzed. Furthermore, the Chinese Glioma Genome Atlas dataset was used for model validation.

Results: A total of 12 m6A/m5C/m1A/m7G-related genes (AL080276.2, AC092111.1, SOX21-AS1, DNAJC9-AS1, AC025171.1, AL356019.2, AC017104.1, AC099850.3, UNC5B-AS1, AC006064.2, AC010319.4, and AC016822.1) were used to construct a survival and prognosis model, which had good independent prediction ability for patients with glioma. Patients were divided into low and high m6A/m5C/m1A/m7G-LS groups, the latter of which had poor prognosis. In addition, the m6A/m5C/m1A/m7G-LS enabled improved interpretation of the results of enrichment analysis, as well as informing immunotherapy response and drug sensitivity of patients with glioma in different subgroups.

Conclusion: In this study we constructed an m6A/m5C/m1A/m7G-LS and established a nomogram model, which can accurately predict the prognosis of patients with glioma and provides direction toward promising immunotherapy strategies for the future.