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.

MeT-DB: a Database of Transcriptome Methylation in Mammalian Cells.

MeT-DB V2.0: Elucidating Context-Specific Functions of N6-methyl-adenosine Methyltranscriptome.

Mapping Messenger RNA Methylations at Single Base Resolution

MethSMRT: an integrative database for DNA N6-methyladenine and N4-methylcytosine generated by single-molecular real-time sequencing

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

CVm6A: A Visualization and Exploration Database for m⁶As in Cell Lines.

m6A-Atlas v2.0: updated resources for unraveling the N 6-methyladenosine (m6A) epitranscriptome among multiple species

CVm6A: A Visualization and Exploration Database for m<sup>6</sup>As in Cell Lines

MethyRNA: a Web Server for Identification of N6-methyladenosine Sites.

MeDReaders: a Database for Transcription Factors That Bind to Methylated DNA.

MethBank: a database integrating next-generation sequencing single-base-resolution DNA methylation programming data.

m<SUP>6</SUP>A-Atlas: a comprehensive knowledgebase for unraveling the <i>N</i><SUP>6</SUP>-methyladenosine (m<SUP>6</SUP>A) epitranscriptome

m6Am-seq reveals the dynamic m6Am methylation in the human transcriptome

MethBank 3.0: a database of DNA methylomes across a variety of species

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

ImmuMethy, a database of DNA methylation plasticity at a single cytosine resolution in human blood and immune cells.

Atlas of quantitative single-base-resolution N6-methyl-adenine methylomes

scMethBank: a database for single-cell whole genome DNA methylation maps

CanMethdb: a database for genome-wide DNA methylation annotation in cancers

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

Landscape and Regulation of m6A and m6Am Methylome across Human and Mouse Tissues