Abstract:Biochemical modifications to mRNA, especially N6-methyladenosine (m(6)A) and 5-methylcytosine (m(5)C), have been recently shown to be associated with crucial biological functions. Despite the intriguing advancements, little is known so far about the dynamic landscape of RNA methylome across different cell types and how the epitranscriptome is regulated at the system level by enzymes, i.e., RNA methyltransferases and demethylases. To investigate this issue, a meta-analysis of m(6)A MeRIP-Seq datasets collected from 10 different experimental conditions (cell type/tissue or treatment) is performed, and the combinatorial epitranscriptome, which consists of 42758 m(6)A sites, is extracted and divided into 3 clusters, in which the methylation sites are likely to be hyper- or hypo-methylated simultaneously (or co-methylated), indicating the sharing of a common methylation regulator. Four different clustering approaches are used, including K-means, hierarchical clustering (HC), Bayesian factor regression model (BFRM) and nonnegative matrix factorization (NMF) to unveil the co-methylation patterns. To validate whether the patterns are corresponding to enzymatic regulators, i.e., RNA methyltransferases or demethylases, the target sites of a known m(6)A regulator, fat mass and obesity-associated protein (FTO), are identified from an independent mouse MeRIP-Seq dataset and lifted to human. Our study shows that 3 out of the 4 clustering approaches used can successfully identify a group of methylation sites overlapping with FTO target sites at a significance level of 0.05 (after multiple hypothesis adjustment), among which, the result of NMF is the most significant (p-value 2.81 x 10(-06)). We defined a new approach evaluating the consistency between two clustering results which shows that clustering results of different methods are highly correlated strongly indicating the existence of co-methylation patterns. Consistent with recent studies, a number of cancer and neuronal disease-related bimolecular functions are enriched in the identified clusters, which are biological functions that can be regulated at the epitranscriptional level, indicating the pharmaceutical prospect of RNA N6-methyladenosine-related studies. This result successfully reveals the linkage between the global RNA co-methylation patterns embedded in the epitranscriptomic data under multiple experimental conditions and the latent enzymatic regulators, suggesting a promising direction towards a more comprehensive understanding of the epitranscriptome.

FunDMDeep-m6A: Identification and Prioritization of Functional Differential M6a Methylation Genes

M6a-Driver: Identifying Context-Specific Mrna M6a Methylation-Driven Gene Interaction Networks.

Global Analysis of N6-methyladenosine Functions and Its Disease Association Using Deep Learning and Network-Based Methods

DLm6Am: A Deep-Learning-Based Tool for Identifying N6,2'-O-Dimethyladenosine Sites in RNA Sequences

EMDL_m6Am: identifying N6,2'-O-dimethyladenosine sites based on stacking ensemble deep learning

Funm6AViewer: a Web Server and R Package for Functional Analysis of Context-Specific M6a RNA Methylation

m6Aminer: Predicting the m6Am Sites on mRNA by Fusing Multiple Sequence-Derived Features into a CatBoost-Based Classifier

EDLm6APred: ensemble deep learning approach for mRNA m6A site prediction

Identification of Functionally Methylated Regions Based on Discriminant Analysis Through Integrating Methylation and Gene Expression Data

BLAM6A-Merge: Leveraging Attention Mechanisms and Feature Fusion Strategies to Improve the Identification of RNA N6-methyladenosine Sites

Decomposition of Rna Methylome Reveals Co-Methylation Patterns Induced by Latent Enzymatic Regulators of the Epitranscriptome

Maexpress-Reader: Prediction of M6a Regulated Expression Genes by Integrating M6a Sites and Reader Binding Information in Specific- Context

Deepm6A-MT: A deep learning-based method for identifying RNA N6-methyladenosine sites in multiple tissues

Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm

M6acancer-Net: Identification of M6a-Mediated Cancer Driver Genes from Gene-Site Heterogeneous Network

TargetM6A: Identifying N6-Methyladenosine Sites from RNA Sequences Via Position-Specific Nucleotide Propensities and a Support Vector Machine

Identification of Differentially Methylated Sites with Weak Methylation Effects.

m <SUP> 6 </SUP> Aexpress-BHM: predicting m<SUP>6</SUP>A regulation of gene expression in multiple-groups context by a Bayesian hierarchical mixture model

<i>m<SUP>6</SUP>A-express</i>: uncovering complex and condition-specific m<SUP>6</SUP>A regulation of gene expression

TS-m6A-DL: Tissue-specific Identification of N6-methyladenosine Sites Using a Universal Deep Learning Model.

Landscape and Regulation of M 6A and M 6am Methylome Across Human and Mouse Tissues