m6A Epitranscriptome Research Articles

Interpretable deep cross networks unveiled common signatures of dysregulated epitranscriptomes across 12 cancer types

Cancer is a complex and multifaceted group of diseases characterized by uncontrolled cell growth that leads to the formation of malignant tumors. Recent studies suggest that N6-methyladenosine (m6A) RNA methylation plays pivotal roles in cancer pathology by influencing various cellular processes. However, the degree to which these mechanisms are shared across different cancer types remains unclear. In this study, we analyzed an expansive array of 167 m6A epitranscriptome profiles covering 12 distinct cancer types and their originating normal tissues. We trained 12 distinct, cancer type-specific interpretable deep cross network models, which successfully distinguish between specific pairs of normal and cancer m6A contexts using integrated information from both the sequences and curated genomic knowledge. Interestingly, cross-cancer type testing indicated the existence of shared genomic patterns across various cancers at the epitranscriptome level. A pan-cancer model was subsequently developed to identify these shared patterns that could not be observed in a single cancer type. Our analysis uncovered, for the first time, a common epitranscriptome signature shared across multiple cancer types, particularly associated with RNA hybridization process and aberrant splicing. This highlights the importance of a comprehensive understanding of the pan-cancer epitranscriptome and holding potential implications in the development of RNA methylation-based therapeutics for various cancers.

Hypoxic stabilization of RIPOR3 mRNA via METTL3-mediated m6A methylation drives breast cancer progression and metastasis.

Dysregulated N6-methyladenosine (m6A) modification has been associated with breast cancer pathogenesis. Hypoxia which characterizes solid tumors is known to reprogram the m6A epitranscriptome, but the underlying mechanisms of how this process contributes to breast cancer progression remain poorly understood. Through integrative analyses of m6A-RIP sequencing and RNA sequencing databases, we reveal a cluster of mRNAs with upregulated m6A methylation and expression under hypoxia, that are enriched by many oncogenic pathways, including PI3K-Akt signaling. Furthermore, we identify the mRNA, RIPOR3, as a target of METTL3-mediated m6A methylation in response to hypoxia. We find that m6A methylation stabilizes RIPOR3, increasing its protein expression in a METTL3 catalytic activity-dependent manner, and consequently driving breast tumor growth and metastasis. RIPOR3 is found to be overexpressed in breast cancer cell lines and tumor tissues from breast cancer patients, in whom elevated RIPOR3 is associated with a worse prognosis. Mechanistically, we show that RIPOR3 interacts with EGFR and is essential for the PI3K-Akt pathway activation. In conclusion, we identify RIPOR3 as a hypoxia-stabilized oncogenic driver via METTL3-mediated m6A methylation, thus provide a potential therapeutic target for breast cancer.

M6ATM: a deep learning framework for demystifying the m6A epitranscriptome with Nanopore long-read RNA-seq data.

N6-methyladenosine (m6A) is one of the most abundant and well-known modifications in messenger RNAs since its discovery in the 1970s. Recent studies have demonstrated that m6A is involved in various biological processes, such as alternative splicing and RNA degradation, playing an important role in a variety of diseases. To better understand the role of m6A, transcriptome-wide m6A profiling data are indispensable. In recent years, the Oxford Nanopore Technology Direct RNA Sequencing (DRS) platform has shown promise for RNA modification detection based on current disruptions measured in transcripts. However, decoding current intensity data into modification profiles remains a challenging task. Here, we introduce the m6A Transcriptome-wide Mapper (m6ATM), a novel Python-based computational pipeline that applies deep neural networks to predict m6A sites at a single-base resolution using DRS data. The m6ATM model architecture incorporates a WaveNet encoder and a dual-stream multiple-instance learning model to extract features from specific target sites and characterize the m6A epitranscriptome. For validation, m6ATM achieved an accuracy of 80% to 98% across in vitro transcription datasets containing varying m6A modification ratios and outperformed other tools in benchmarking with human cell line data. Moreover, we demonstrated the versatility of m6ATM in providing reliable stoichiometric information and used it to pinpoint PEG10 as a potential m6A target transcript in liver cancer cells. In conclusion, m6ATM is a high-performance m6A detection tool, and our results pave the way for future advancements in epitranscriptomic research.

New insights into the dynamics of m6A epitranscriptome: hybrid-seq identifies novel mRNAs of the m6A writers METTL3/14.

Background: N6-methyladenosine (m6A), a prevalent mRNA modification, is dynamically regulated by methyltransferases, including METTL3 and METTL14.Materials & methods: In the current study, we employed a custom hybrid-seq method to identify novel METTL3/14 transcripts, explore their protein-coding capacities and predict the putative role of the METTL isoforms.Results: Demultiplexing of the hybrid-seq barcoded datasets unraveled the expression patterns of the newly identified mRNAs in major malignancies as well as in non-malignant cells, providing a deeper understanding of the methylation pathways. Open reading framequery revealed novel METTL3/14 isoforms, broadening our perspective for the structural diversity within METTL family.Conclusion: Our findings offer significant insights into the intricate transcriptional landscape of METTL3/14, shedding light on the regulatory mechanisms underlying methylation in mRNAs.

The N6-methyladenosine Epitranscriptomic Landscape of Lung Adenocarcinoma.

Comprehensive m6A epitranscriptome profiling of primary tumors remains largely uncharted. Here, we profiled the m6A epitranscriptome of 10 non-neoplastic lung (NL) tissues and 51 lung adenocarcinoma (LUAD) tumors, integrating the corresponding transcriptome, proteome and extensive clinical annotations. We identified distinct clusters and genes that were exclusively linked to disease progression through m6A modifications. In comparison with NL tissues, we identified 430 transcripts to be hypo-methylated and 222 to be hyper-methylated in tumors. Among these genes, EML4 emerged as a novel metastatic driver, displaying significant hyper-methylation in tumors. m6A modification promoted the translation of EML4, leading to its widespread overexpression in primary tumors. Functionally, EML4 modulated cytoskeleton dynamics through interacting with ARPC1A, enhancing lamellipodia formation, cellular motility, local invasion, and metastasis. Clinically, high EML4 protein abundance correlated with features of metastasis. METTL3 small molecule inhibitor markedly diminished both EML4 m6A and protein abundance, and efficiently suppressed lung metastases in vivo.

Antibiotic-Induced Gut Microbiota Dysbiosis Modulates Host Transcriptome and m6A Epitranscriptome via Bile Acid Metabolism.

Gut microbiota can influence host gene expression and physiology through metabolites. Besides, the presence or absence of gut microbiome can reprogram host transcriptome and epitranscriptome as represented by N6-methyladenosine (m6A), the most abundant mammalian mRNA modification. However, which and how gut microbiota-derived metabolites reprogram host transcriptome and m6A epitranscriptome remain poorly understood. Here, investigation is conducted into how gut microbiota-derived metabolites impact host transcriptome and m6A epitranscriptome using multiple mouse models and multi-omics approaches. Various antibiotics-induced dysbiotic mice are established, followed by fecal microbiota transplantation (FMT) into germ-free mice, and the results show that bile acid metabolism is significantly altered along with the abundance change in bile acid-producing microbiota. Unbalanced gut microbiota and bile acids drastically change the host transcriptome and the m6A epitranscriptome in multiple tissues. Mechanistically, the expression of m6A writer proteins is regulated in animals treated with antibiotics and in cultured cells treated with bile acids, indicating a direct link between bile acid metabolism and m6A biology. Collectively, these results demonstrate that antibiotic-induced gut dysbiosis regulates the landscape of host transcriptome and m6A epitranscriptome via bile acid metabolism pathway. This work provides novel insights into the interplay between microbial metabolites and host gene expression.

Small-molecule inhibition of the METTL3/METTL14 complex suppresses neuroblastoma tumor growth and promotes differentiation

The N6-methyladenosine (m6A) RNA modification is an important regulator of gene expression. m6A is deposited by a methyltransferase complex that includes methyltransferase-like 3 (METTL3) and methyltransferase-like 14 (METTL14). High levels of METTL3/METTL14 drive the growth of many types of adult cancer, and METTL3/METTL14 inhibitors are emerging as new anticancer agents. However, little is known about the m6A epitranscriptome or the role of the METTL3/METTL14 complex in neuroblastoma, a common pediatric cancer. Here, we show that METTL3 knockdown or pharmacologic inhibition with the small molecule STM2457 leads to reduced neuroblastoma cell proliferation and increased differentiation. These changes in neuroblastoma phenotype are associated with decreased m6A deposition on transcripts involved in nervous system development and neuronal differentiation, with increased stability of target mRNAs. In preclinical studies, STM2457 treatment suppresses the growth of neuroblastoma tumors invivo. Together, these results support the potential of METTL3/METTL14 complex inhibition as a therapeutic strategy against neuroblastoma.

M6ACali: machine learning-powered calibration for accurate m6A detection in MeRIP-Seq.

We present m6ACali, a novel machine-learning framework aimed at enhancing the accuracy of N6-methyladenosine (m6A) epitranscriptome profiling by reducing the impact of non-specific antibody enrichment in MeRIP-Seq. The calibration model serves as a genomic feature-based classifier that refines the identification of m6A sites, distinguishing those genuinely present from those that can be detected in in-vitro transcribed (IVT) control experiments. We find that m6ACali effectively identifies non-specific binding peaks reported by exomePeak2 and MACS2 in novel MeRIP-Seq datasets without the need for paired IVT controls. The model interpretation revealed that off-target antibody binding sites commonly occur at short exons and short mRNAs, originating from high read coverage regions that share the motif sequence with true m6A sites. We also reveal that the ML strategy can efficiently adjust differentially methylated peaks and other antibody-dependent, base-resolution m6A detection techniques. As a result, m6ACali offers a promising method for the universal enhancement of m6A profiles generated by MeRIP-Seq experiments, elevating the benchmark for omics-level m6A data integration.

Reading the epitranscriptome of the human malaria parasite

Epigenetic machinery has emerged as a central player in gene regulation and chromatin organization in Plasmodium spp. Epigenetic modifications on histones and their role in antigenic variation in P. falciparum are widely studied. Recent discoveries on nucleic acid methylome are exciting and provide a new dimension to the apicomplexan protozoan parasite's gene regulatory process. Reports have confirmed that N6-methyl adenosine (m6A) methylation plays a crucial role in the translational plasticity of the human malaria parasite during its development in RBC. The YTH domain (YT521-B Homology) protein in P. falciparum binds to m6A epitranscriptome modifications on the mRNA and regulates protein translation. The binding of the PfYTH domain protein to the m6A-modified mRNA is mediated through a binding pocket formed by aromatic amino acids. The P. falciparum genome encodes two members of YTH domain proteins, i.e., YTH1 and YTH2, and both have distinct roles in dictating the epitranscriptome in human malaria parasites. This review highlights recent advancements in the functions and mechanisms of YTH domain protein's role in translational plasticity in the various developmental stages of the parasite.

The RNA m6A writer METTL3 in tumor microenvironment: emerging roles and therapeutic implications.

The tumor microenvironment (TME) is a heterogeneous ecosystem comprising cancer cells, immune cells, stromal cells, and various non-cellular components, all of which play critical roles in controlling tumor progression and response to immunotherapies. Methyltransferase-like 3 (METTL3), the core component of N 6-methyladenosine (m6A) writer, is frequently associated with abnormalities in the m6A epitranscriptome in different cancer types, impacting both cancer cells and the surrounding TME. While the impact of METTL3 on cancer cells has been extensively reviewed, its roles in TME and anti-cancer immunity have not been comprehensively summarized. This review aims to systematically summarize the functions of METTL3 in TME, particularly its effects on tumor-infiltrating immune cells. We also elaborate on the underlying m6A-dependent mechanism. Additionally, we discuss ongoing endeavors towards developing METTL3 inhibitors, as well as the potential of targeting METTL3 to bolster the efficacy of immunotherapy.

RNA N6-Methyladenosine Affects Copper-Induced Oxidative Stress Response in Arabidopsis thaliana.

Recently, post-transcriptional regulation of mRNA mediated by N6-methyladenosine (m6A) has been found to have profound effects on transcriptome regulation during plant responses to various abiotic stresses. However, whether this RNA modification can affect an oxidative stress response in plants has not been studied. To assess the role of m6A modifications during copper-induced oxidative stress responses, m6A-IP-seq was performed in Arabidopsis seedlings exposed to high levels of copper sulfate. This analysis revealed large-scale shifts in this modification on the transcripts most relevant for oxidative stress. This altered epitranscriptomic mark is known to influence transcript abundance and translation; therefore we scrutinized these possibilities. We found an increased abundance of copper-enriched m6A-containing transcripts. Similarly, we also found increased ribosome occupancy of copper-enriched m6A-containing transcripts, specifically those encoding proteins involved with stress responses relevant to oxidative stressors. Furthermore, the significance of the m6A epitranscriptome on plant oxidative stress tolerance was uncovered by assessing germination and seedling development of the mta (N6-methyladenosine RNA methyltransferase A mutant complemented with ABI3:MTA) mutant exposed to high copper treatment. These analyses suggested hypersensitivity of the mta mutant compared to the wild-type plants in response to copper-induced oxidative stress. Overall, our findings suggest an important role for m6A in the oxidative stress response of Arabidopsis.

M6A epitranscriptome analysis reveals differentially methylated transcripts that drive early chemoresistance in bladder cancer.

N 6-Methyladenosine (m6A) RNA modifications dynamically regulate messenger RNA processing, differentiationand cell fate. Given these functions, we hypothesized that m6A modifications play a role in the transition to chemoresistance. To test this, we took an agnostic discovery approach anchored directly to chemoresistance rather than to any particular m6A effector protein. Specifically, we used methyl-RNA immunoprecipitation followed by sequencing (MeRIP-seq) in parallel with RNA sequencing to identify gene transcripts that were both differentially methylated and differentially expressed between cisplatin-sensitive and cisplatin-resistant bladder cancer (BC) cells. We filtered and prioritized these genes using clinical and functional database tools, andthen validated several of the top candidates via targeted quantitative polymerase chain reaction (qPCR) and MeRIP-PCR. In cisplatin-resistant cells, SLC7A11transcripts had decreased methylation associated with decreased m6A reader YTHDF3 binding, prolonged RNA stability, and increased RNA and protein levels, leading to reduced ferroptosis and increased survival. Consistent with this, cisplatin-sensitive BC cell lines and patient-derived organoids exposed to cisplatin for as little as 48 h exhibited similar mechanisms of SLC7A11 upregulation and chemoresistance, trends that were also reflected in public cancer survival databases. Collectively, these findings highlight epitranscriptomic plasticity as a mechanism of rapid chemoresistance and a potential therapeutic target.

A deregulated m6A writer complex axis driven by BRD4 confers an epitranscriptomic vulnerability in combined DNA repair–targeted therapy

Aberrant transcripts expression of the m6A methyltransferase complex (MTC) is widely found across human cancers, suggesting a dysregulated signaling cascade which integrates m6A epitranscriptome to drive tumorigenesis. However, the responsible transcriptional machinery directing the expression of distinct MTC subunits remains unclear. Here, we identified an unappreciated interplay between the histone acetyl-lysine reader BRD4 and the m6A writer complex across human cancers. BRD4 directly stimulates transcripts expression of seven MTC subunits, allowing the maintenance of the nuclear writer complex integrity. Upon BET inhibition, this BRD4-MTC signaling cascade accounts for global m6A reduction and the subsequent dynamic alteration of BRD4-dependent transcriptome, resulting in impaired DNA damage response that involves activation of homologous recombination (HR) repair and repression of apoptosis. We further demonstrated that the combined synergy upon BET/PARP inhibition largely relies on disrupted m6A modification of HR and apoptotic genes, counteracting PARP inhibitor (PARPi) resistance in patient-derived xenograft models. Our study revealed a widespread active cross-talk between BRD4-dependent epigenetic and MTC-mediated epitranscriptomic networks, which provides a unique therapeutic vulnerability that can be leveraged in combined DNA repair-targeted therapy.

Exon junction complex modulates the formation of the m6A epitranscriptome.

N6-methyladenosine (m6A) is one of the most abundant chemical modifications in mRNA and plays essential roles in diverse physiological and pathological processes. m6A is highly enriched near stop codons and in long internal exons of mRNA, but the mechanism leading to this specific distribution has been unclear. Recently, three papers have solved this major problem by revealing that exon junction complexes (EJCs) act as m6A suppressors and shape the formation of the m6A epitranscriptome. Here, we briefly introduce the m6A pathway, elaborate the roles of EJC on the formation of m6A modification based on these results, and describe the effect of exon-intron structure on mRNA stability via m6A, which will help us better understand the latest progress in the m6A RNA modification field.

N6-methyladenosine is required for efficient RNA synthesis of Ebola virus and other haemorrhagic fever viruses

ABSTRACT N6-methyladenosine (m6A) is one of the most abundant modifications of cellular RNA, where it serves various functions. m6A methylation of many viral RNA species has also been described; however, little is known about the m6A epitranscriptome of haemorrhagic fever-causing viruses like Ebola virus (EBOV). Here, we analysed the importance of the methyltransferase METTL3 for the life cycle of this virus. We found that METTL3 interacts with the EBOV nucleoprotein and the transcriptional activator VP30 to support viral RNA synthesis, and that METTL3 is recruited into EBOV inclusions bodies, where viral RNA synthesis occurs. Analysis of the m6A methylation pattern of EBOV mRNAs showed that they are methylated by METTL3. Further studies revealed that METTL3 interaction with the viral nucleoprotein, as well as its importance for RNA synthesis and protein expression, is also observed for other haemorrhagic fever viruses such as Junín virus (JUNV) and Crimean-Congo haemorrhagic fever virus (CCHFV). The negative effects on viral RNA synthesis due to loss of m6A methylation are independent of innate immune sensing, as METTL3 knockout did not affect type I interferon induction in response to viral RNA synthesis or infection. Our results suggest a novel function for m6A that is conserved among diverse haemorrhagic fever-causing viruses (i.e. EBOV, JUNV and CCHFV), making METTL3 a promising target for broadly-acting antivirals.

Roles of N6-methyladenosine epitranscriptome in non-alcoholic fatty liver disease and hepatocellular carcinoma.

Non-alcoholic fatty liver disease (NAFLD) is a typical chronic liver disease connected to a high risk of developing hepatocellular carcinoma (HCC). The development of NAFLD and HCC has been associated with changes in epigenetics, such as histone modifications and micro RNA (miRNA)-mediated processes. Recently, in the realm of epitranscriptomics, RNA alterations have become important regulators. N6-methyladenosine (m6A) is the most common and crucial alteration for controlling mRNA stability, splicing, and translation. It is particularly important for controlling liver disease progression and hepatic function. This review aims to conclude recent research on the functions of m6A epitranscriptome in the molecular mechanisms behind NAFLD and HCC development, with special attention to the effects of m6A alteration on how HCC develops and its possible roles in the progression of NAFLD to HCC. Additionally, the review discusses the possible effects of m6A alteration on the treatment and diagnostic of NAFLD and HCC. It is crucial to remember that m6A modification is a reversible action controlled via the coordinated functions of the proteins that write and delete, enabling quick adaptability to environmental changes. The review also discusses m6A-binding proteins' function in mRNA alternative splicing, translation, and degradation and their ability to modulate mRNA stability and processing. Understanding RNA modification regulation and its part in the emergence of HCC and NAFLD may provide new avenues for diagnosing and treating these diseases.

GPX8 deficiency–induced oxidative stress reprogrammed m6A epitranscriptome of oral cancer cells

ABSTRACT Glutathione peroxidase 8 (GPX8) is a key regulator of redox homoeostasis. Whether its antioxidant activity participates in the regulation of m6A modification is a crucial issue, which has important application value in cancer treatment. In this study, MeRIP-seq was used to explore the characteristics of transcriptome-wide m6A modification in GPX8-deficient oral cancer cells. Oxidative stress caused by the lack of GPX8 resulted in 1,279 hyper- and 2,287 hypo-methylated m6A peaks and 2,036 differentially expressed genes in GPX8-KO cells. Twenty-eight differentially expressed genes were related to the cell response to oxidative stress, and half of them changed their m6A modification. In GPX8-KO cells, m6A regulators IGF2BP2 and IGF2BP3 were upregulated, while FTO, RBM15, VIRMA, ZC3H13, and YTHDC2 were downregulated. After H2O2 treatment, the expression changes of RBM15, IGF2BP2, and IGF2BP3 were further enhanced. These data indicated that GPX8-mediated redox homoeostasis regulated m6A modification, thereby affecting the expression and function of downstream genes. This study highlights the possible significance of GPX8 and the corresponding m6A regulatory or regulated genes as novel targets for antioxidant intervention in cancer therapy.

M6A epitranscriptomic regulation of tissue homeostasis during primate aging.

How N6-methyladenosine (m6A), the most abundant mRNA modification, contributes to primate tissue homeostasis and physiological aging remains elusive. Here, we characterize the m6A epitranscriptome across the liver, heart and skeletal muscle in young and old nonhuman primates. Our data reveal a positive correlation between m6A modifications and gene expression homeostasis across tissues as well as tissue-type-specific aging-associated m6A dynamics. Among these tissues, skeletal muscle is the most susceptible to m6A loss in aging and shows a reduction in the m6A methyltransferase METTL3. We further show that METTL3 deficiency in human pluripotent stem cell-derived myotubes leads to senescence and apoptosis, and identify NPNT as a key element downstream of METTL3 involved in myotube homeostasis, whose expression and m6A levels are both decreased in senescent myotubes. Our study provides a resource for elucidating m6A-mediated mechanisms of tissue aging and reveals a METTL3-m6A-NPNT axis counteracting aging-associated skeletal muscle degeneration.

DNER and GNL2 are differentially m6A methylated in periodontitis in comparison with periodontal health revealed by m6A microarray of human gingival tissue and transcriptomic analysis

This study aims to investigate the differences in the epigenomic patterns of N6-methyladenosine (m6A) methylation in gingival tissues between patients with periodontitis (PD) and healthy controls, identifying potential biomarkers. As a multifactorial disease, PD involves multiple genetic and environmental effects. The m6A modification is the most prevalent internal mRNA modification and linked to various inflammatory diseases. However, the m6A modification pattern and m6A-related signatures in PD remain unclear. An m6A microarray of human gingival tissues was conducted in eight subjects: four diagnosed with PD and four healthy controls. Microarray analysis was performed to identify the differentially m6A methylated mRNAs (DMGs) and the differentially expressed mRNAs (DEGs). The differentially methylated and expressed mRNAs (DMEGs) were subjected to functional enrichment analysis by Metascape. The weighted gene co-expression network analysis (WGCNA) algorithm, the least absolute shrinkage and selection operator (LASSO) regression, and univariate logistic regression were performed to identify potential biomarkers. The cell type localization of the target genes was determined using single-cell RNA-seq (scRNA-seq) analysis. The m6A methylation level and gene expression of hub genes were subsequently verified by m6A methylated RNA immunoprecipitation (MeRIP) and quantitative real-time PCR (qRT-PCR). In total, 458 DMGs, 750 DEGs, and 279 DMEGs were identified based on our microarray. Pathway analyses conducted for the DMEGs revealed that biological functions were mainly involved in the regulation of stem cell differentiation, ossification, circadian rhythm, and insulin secretion pathways. Besides, the genes involved in crucial biological processes were mainly expressed in fibroblast and epithelial cells. Furthermore, the m6A methylation and expression levels of two hub biomarkers (DNER and GNL2) were validated. The current study exhibited a distinct m6A epitranscriptome, identified and verified two PD-related biomarkers (DNER and GNL2), which may provide novel insights into revealing the new molecular mechanisms and latent targets of PD.

Conserved reduction of m6A RNA modifications during aging and neurodegeneration is linked to changes in synaptic transcripts

N6-methyladenosine (m6A) regulates mRNA metabolism. While it has been implicated in the development of the mammalian brain and in cognition, the role of m6A in synaptic plasticity, especially during cognitive decline, is not fully understood. In this study, we employed methylated RNA immunoprecipitation sequencing to obtain the m6A epitranscriptome of the hippocampal subregions CA1, CA3, and the dentate gyrus and the anterior cingulate cortex (ACC) in young and aged mice. We observed a decrease in m6A levels in aged animals. Comparative analysis of cingulate cortex (CC) brain tissue from cognitively intact human subjects and Alzheimer's disease (AD) patients showed decreased m6A RNA methylation in AD patients. m6A changes common to brains of aged mice and AD patients were found in transcripts linked to synaptic function including calcium/calmodulin-dependent protein kinase 2 (CAMKII) and AMPA-selective glutamate receptor 1 (Glua1). We used proximity ligation assays to show that reduced m6A levels result in decreased synaptic protein synthesis as exemplified by CAMKII and GLUA1. Moreover, reduced m6A levels impaired synaptic function. Our results suggest that m6A RNA methylation controls synaptic protein synthesis and may play a role in cognitive decline associated with aging and AD.