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.

A Novel Analytical Strategy To Identify Fusion Transcripts Between Repetitive Elements And Protein Coding-Exons Using Rna-Seq

Long Read, Isoform Aware Sequencing of Mouse Nucleus Accumbens after Chronic Cocaine Treatment.

De Novo Identification of Expressed Cancer Somatic Mutations from Single-Cell RNA Sequencing Data

Chronic Cocaine-Regulated Epigenomic Changes in Mouse Nucleus Accumbens.

Droplet-based Transcriptome Profiling of Individual Synapses.

Fusion Transcripts And Transcribed Retrotransposed Loci Discovered Through Comprehensive Transcriptome Analysis Using Paired-End Ditags (Pets)

1148P Identification and Validation of RET Fusions in Lung Adenocarcinoma Through DNA and RNA Sequencing

Optimized Exon-Exon Junction Library and its Application on Rodents' Brain Transcriptome Analysis

Transposable Element (TE) insertion predictions from RNAseq inputs and TE impact on RNA splicing and gene expression in Drosophila brain transcriptomes

Transposable Element (TE) insertion predictions from RNAseq inputs and TE impact on RNA splicing and gene expression in Drosophila brain transcriptomes.

Decoding RNA Metabolism by RNA-linked CRISPR Screening in Human Cells

Ab Initio Reconstruction of Cell Type–specific Transcriptomes in Mouse Reveals the Conserved Multi-Exonic Structure of Lincrnas

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

Enhancer RNAs predict enhancer–gene regulatory links and are critical for enhancer function in neuronal systems

Identification and analysis of mouse non-coding RNA using transcriptome data

Comprehensive Evaluation of Fusion Transcript Detection Algorithms and a Meta-Caller to Combine Top Performing Methods in Paired-End RNA-seq Data.

A single-nucleus transcriptomic atlas of medium spiny neurons in the rat nucleus accumbens

Transcriptome-Wide Profiling of Nascent RNA in Neurons with Enriched H3K27ac Signal Elevates eRNA Identification Efficiency

Integrated Analysis of a Compendium of RNA-Seq Datasets for Splicing Factors.

Detection of Alu Exonization Events in Human Frontal Cortex From RNA-Seq Data

High-throughput identification of RNA localization elements reveals a regulatory role for A/G rich sequences