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.

Identification of additional regulatory RNPs that impact rRNA and U6 snRNA methylation

Rna N-6-Methyladenosine Methylation In Post-Transcriptional Gene Expression Regulation

Nucleotide Modifications of Eukaryotic Rrnas: the World of Small Nucleolar Rna Guides Revisited

Beyond rRNA and snRNA: tRNA as a 2'-O-methylation target for nucleolar and Cajal body box C/D RNPs

Intron-Encoded, Antisense Small Nucleolar Rnas: the Characterization of Nine Novel Species Points to Their Direct Role As Guides for the 2'-O-ribose Methylation of Rrnas

RNA helicase-mediated regulation of snoRNP dynamics on pre-ribosomes and rRNA 2′-O-methylation

Comprehensive map of ribosomal 2′-O-methylation and C/D box snoRNAs in Drosophila melanogaster

Tuning the ribosome: The influence of rRNA modification on eukaryotic ribosome biogenesis and function

Conserved Composition of Mammalian Box H/ACA and Box C/D Small Nucleolar Ribonucleoprotein Particles and Their Interaction with the Common Factor Nopp140.

Regulation of translation by site-specific ribosomal RNA methylation

Alteration of 28S rRNA 2′-O-methylation by etoposide correlates with decreased SMN phosphorylation and reduced Drosha levels

A snoRNA-tRNA modification network governs codon-biased cellular states

2'-O-methylation in RNA: progress, challenges, and future directions

Turning Uridines around: Role of rRNA Pseudouridylation in Ribosome Biogenesis and Ribosomal Function

A genome-wide analysis of C/D and H/ACA-like small nucleolar RNAs in Trypanosoma brucei reveals a trypanosome-specific pattern of rRNA modification.

RNA modification in Cajal bodies

Ribosomal RNA Pseudouridylation: Will Newly Available Methods Finally Define the Contribution of This Modification to Human Ribosome Plasticity?

The Arabidopsis 2′-O-Ribose-Methylation and Pseudouridylation Landscape of rRNA in Comparison to Human and Yeast

Activity‐Based Profiling of RNA Modifying Enzymes

Modification of messenger RNA by 2′-O-methylation regulates gene expression in vivo

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