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.

DRUM: Inference of Disease-Associated M6a RNA Methylation Sites from a Multi-Layer Heterogeneous Network

Predict disease-related RNA methylation sites from the methylation-expression association by using hypergeometric test

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

Large-scale Functional Annotation of Individual RNA Methylation Sites by Mining Complex Biological Networks

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

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

Predicting Disease-Associated N7-Methylguanosine (m7g) Sites Via Random Walk on Heterogeneous Network.

Prediction of RNA Methylation Status from Gene Expression Data Using Classification and Regression Methods

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

RMDGCN: Prediction of RNA Methylation and Disease Associations Based on Graph Convolutional Network with Attention Mechanism.

EDLm6APred: ensemble deep learning approach for mRNA m6A site prediction

Defining disease-related modules based on weighted miRNA synergistic network

M6acomet: Large-Scale Functional Prediction of Individual M6a RNA Methylation Sites from an RNA Co-Methylation Network.

Inferring Latent MicroRNA-Disease Associations on a Gene-Mediated Tripartite Heterogeneous Multiplexing Network.

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

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

Detecting m6A methylation regions from Methylated RNA Immunoprecipitation Sequencing

A Framework for Integrating Multiple Biological Networks to Predict Microrna-Disease Associations

Significance of RNA N6-methyladenosine regulators in the diagnosis and subtype classification of coronary heart disease using the Gene Expression Omnibus database

ConsRM: collection and large-scale prediction of the evolutionarily conserved RNA methylation sites, with implications for the functional epitranscriptome

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