Abstract:Despite the prevalent studies of DNA/Chromatin related epigenetics, such as, histone modifications and DNA methylation, RNA epigenetics did not receive deserved attention due to the lack of high throughput approach for profiling epitranscriptome. Recently, a new affinity-based sequencing approach MeRIPseq was developed and applied to survey the global mRNA N6-methyladenosine (m6A) in mammalian cells. As a marriage of ChIPseq and RNAseq, MeRIPseq has the potential to study, for the first time, the transcriptome-wide distribution of different types of post-transcriptional RNA modifications. Yet, this technology introduced new computational challenges that have not been adequately addressed. We have previously developed a MATLAB-based package `exomePeak' for detection of RNA methylation sites from MeRIPseq data. Here, we extend the features of exomePeak by including a novel computational framework that enables differential analysis to unveil the dynamics in RNA epigenetic regulations. The novel differential analysis monitors the percentage of modified RNA molecules among the total transcribed RNAs, which directly reflects the impact of RNA epigenetic regulations. In contrast, current available software packages developed for sequencing-based differential analysis such as DESeq or edgeR monitors the changes in the absolute amount of molecules, and, if applied to MeRIPseq data, might be dominated by transcriptional gene differential expression. The algorithm is implemented as an R-package `exomePeak' and freely available. It takes directly the aligned BAM files as input, statistically supports biological replicates, corrects PCR artifacts, and outputs exome-based results in BED format, which is compatible with all major genome browsers for convenient visualization and manipulation. Examples are also provided to depict how exomePeak R-package is integrated with exiting tools for MeRIPseq based peak calling and differential analysis. Particularly, the rationales behind each processing step as well as the specific method used, the best practice, and possible alternative strategies are briefly discussed. The algorithm was applied to the human HepG2 cell MeRIPseq data sets and detects more than 16000 RNA m6A sites, many of which are differentially methylated under ultraviolet radiation. The challenges and potentials of MeRIPseq in epitranscriptome studies are discussed in the end.

Differential analysis of RNA methylome with improved spatial resolution

Spatially Enhanced Differential RNA Methylation Analysis from Affinity-Based Sequencing Data with Hidden Markov Model.

Differential RNA Methylation Analysis for MeRIP-seq Data under General Experimental Design

Evaluation of epitranscriptome-wide N6-methyladenosine differential analysis methods

Unveiling the Dynamics in RNA Epigenetic Regulations.

A Protocol for RNA Methylation Differential Analysis with MeRIP-Seq Data and Exomepeak R/Bioconductor Package.

DRME: Count-based Differential RNA Methylation Analysis at Small Sample Size Scenario

Detecting m6A methylation regions from Methylated RNA Immunoprecipitation Sequencing

Mapping Messenger RNA Methylations at Single Base Resolution

Modeling Of Replicates Variances For Detecting Rna Methylation Site In Merip-Seq Data

HMM-Fisher: identifying differential methylation using a hidden Markov model and Fisher’s exact test

QNB: differential RNA methylation analysis for count-based small-sample sequencing data with a quad-negative binomial model

Recent Progress And Challenges In High Throughput Rna Methylation Sequencing Data Analysis

Identification of Differentially Methylated Single-Nucleotide M6a Sites by Incorporating Site-Specific Antibody Specificity

Identification of differentially methylated single-nucleotide m A sites by incorporating site-specific antibody specificity

MeRIP-PF: an Easy-to-use Pipeline for High-Resolution Peak-Finding in MeRIP-Seq Data.

Detection of differentially methylated regions from bisulfite-seq data by hidden Markov models incorporating genome-wide methylation level distributions

A Comprehensive Evaluation of Alignment Software for Reduced Representation Bisulfite Sequencing Data

Detecting differentially methylated mRNA from MeRIP-Seq with likelihood ratio test.

Qdmr: A Quantitative Method For Identification Of Differentially Methylated Regions By Entropy

A comprehensive evaluation of computational tools to identify differential methylation regions using RRBS data