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.

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

Single-base mapping of m<sup>6</sup>A by an antibody-independent method

Single-base mapping of m6A by an antibody-independent method

Differential analysis of RNA methylome with improved spatial resolution

Evaluation of epitranscriptome-wide N6-methyladenosine differential analysis methods

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

Detecting m6A methylation regions from Methylated RNA Immunoprecipitation Sequencing

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

TargetM6A: Identifying N6-Methyladenosine Sites from RNA Sequences Via Position-Specific Nucleotide Propensities and a Support Vector Machine

m6ACali: machine learning-powered calibration for accurate m6A detection in MeRIP-Seq

Endogenous labeling empowers accurate detection of m A from single long reads of direct RNA sequencing

Mass-Spectrometry-Based Assay at Single-Base Resolution for Simultaneously Detecting M6a and M6am in RNA

Detection Methods of Epitranscriptomic Mark N6-methyladenosine.

Comprehensive analysis across mammalian tissues and cells decipher the underlying mechanism of m6A specificity

Precise Antibody-Independent M6a Identification Via 4Sedttp-Involved and FTO-Assisted Strategy at Single-Nucleotide Resolution.

Refined RIP-seq protocol for epitranscriptome analysis with low input materials

m6AConquer: a Data Resource for Unified Quantification and Integration of m6A Detection Techniques

Computational Identification of N6-methyladenosine Sites in Multiple Tissues of Mammals.

Unveiling the Dynamics in RNA Epigenetic Regulations.

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