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.

BDBB: A Novel Beta-distribution-based Biclustering Algorithm for Revealing Local Co-methylation Patterns in Epi-transcriptome Profiling Data

BBM: A novel beta-binomial-distribution-based biclustering algorithm for mining m 6 A co-methylation patterns

Genome-scale sequence data processing and epigenetic analysis of DNA methylation]

A novel family of beta mixture models for the differential analysis of DNA methylation data: An application to prostate cancer

REW-ISA V2: A Biclustering Method Fusing Homologous Information for Analyzing and Mining Epi-Transcriptome Data

A family of mixture models for beta valued DNA methylation data

A Hybrid Method Based on Beta-Binomial Bayesian Hierarchical Model and Combination of Ranking Method to Detect Differential Methylation Regions in Bisulfite Sequencing Data

A Novel Computational Method For Detecting Dna Methylation Sites With Dna Sequence Information And Physicochemical Properties

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

A Classification Approach for DNA Methylation Profiling with Bisulfite Next-Generation Sequencing Data

REW-ISA: unveiling local functional blocks in epi-transcriptome profiling data via an RNA expression-weighted iterative signature algorithm

Identification of Methylation States of DNA Regions for Illumina Methylation BeadChip

Gbdmr: identifying differentially methylated CpG regions in the human genome via generalized beta regressions

Higher order methylation features for clustering and prediction in epigenomic studies

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

A Bayesian framework for genome-wide inference of DNA methylation levels

HMM-DM: identifying differentially methylated regions using a hidden Markov model

Functional Grouping of Yeast Genes Via Biclustering Microarray Data

High resolution profiling of human exon methylation by liquid hybridization capture-based bisulfite sequencing

A New Strategy of Cooperativity of Biclustering and Hierarchical Clustering: a Case of Analyzing Yeast Genomic Microarray Datasets.

A Framework for Analyzing DNA Methylation Data from Illumina Infinium HumanMethylation450 BeadChip