MeTDiff: A Novel Differential RNA Methylation Analysis for MeRIP-Seq Data

N6-Methyladenosine (m6A) transcriptome methylation is an exciting new research area that just captures the attention of research community. We present in this paper, MeTDiff, a novel computational tool for predicting differential m6A methylation sites from Methylated RNA immunoprecipitation sequencing (MeRIP-Seq) data. Compared with the existing algorithm exomePeak, the advantages of MeTDiff are that it explicitly models the reads variation in data and also devices a more power likelihood ratio test for differential methylation site prediction. Comprehensive evaluation of MeTDiff's performance using both simulated and real datasets showed that MeTDiff is much more robust and achieved much higher sensitivity and specificity over exomePeak. The R package “MeTDiff” and additional details are available at: https://github.com/compgenomics/MeTDiff.

[1]  D. Bates,et al.  Newton-Raphson and EM Algorithms for Linear Mixed-Effects Models for Repeated-Measures Data , 1988 .

[2]  Yufei Huang,et al.  HEPeak: an HMM-based exome peak-finding package for RNA epigenome sequencing data , 2015, BMC Genomics.

[3]  Chengqi Yi,et al.  N6-Methyladenosine in Nuclear RNA is a Major Substrate of the Obesity-Associated FTO , 2011, Nature chemical biology.

[4]  A. Krainer,et al.  THE SPLICING FACTOR SRSF1 REGULATES APOPTOSIS AND PROLIFERATION TO PROMOTE MAMMARY EPITHELIAL CELL TRANSFORMATION , 2011, Nature Structural &Molecular Biology.

[5]  O. Elemento,et al.  Comprehensive Analysis of mRNA Methylation Reveals Enrichment in 3′ UTRs and near Stop Codons , 2012, Cell.

[6]  Amy J. Hawkins,et al.  SRSF1 Regulates the Alternative Splicing of Caspase 9 Via A Novel Intronic Splicing Enhancer Affecting the Chemotherapeutic Sensitivity of Non–Small Cell Lung Cancer Cells , 2011, Molecular Cancer Research.

[7]  Clifford A. Meyer,et al.  Model-based Analysis of ChIP-Seq (MACS) , 2008, Genome Biology.

[8]  Simon Hess,et al.  The fat mass and obesity associated gene (Fto) regulates activity of the dopaminergic midbrain circuitry , 2013, Nature Neuroscience.

[9]  K. Robertson.  DNA methylation and human disease , 2005, Nature Reviews Genetics.

[10]  Yang Wang,et al.  N6-methyladenosine modification destabilizes developmental regulators in embryonic stem cells , 2014, Nature Cell Biology.

[11]  A. Bird,et al.  DNA methylation landscapes: provocative insights from epigenomics , 2008, Nature Reviews Genetics.

[12]  Minoru Yoshida,et al.  RNA-Methylation-Dependent RNA Processing Controls the Speed of the Circadian Clock , 2013, Cell.

[13]  R J Roberts,et al.  Sequence specificity of the human mRNA N6-adenosine methylase in vitro. , 1990, Nucleic acids research.

[14]  Yong Zhang,et al.  Identifying ChIP-seq enrichment using MACS , 2012, Nature Protocols.

[15]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[16]  Xiaodong Cui,et al.  Exome-based analysis for RNA epigenome sequencing data , 2013, Bioinform..

[17]  M. Kupiec,et al.  Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq , 2012, Nature.

[18]  E. S. Pearson,et al.  On the Problem of the Most Efficient Tests of Statistical Hypotheses , 1933 .

[19]  K. Zhao,et al.  ChIP-Seq: technical considerations for obtaining high-quality data , 2011, Nature Immunology.

[20]  Cole Trapnell,et al.  Computational methods for transcriptome annotation and quantification using RNA-seq , 2011, Nature Methods.

[21]  Schraga Schwartz,et al.  High-Resolution Mapping Reveals a Conserved, Widespread, Dynamic mRNA Methylation Program in Yeast Meiosis , 2013, Cell.

[22]  Yufei Huang,et al.  A protocol for RNA methylation differential analysis with MeRIP-Seq data and exomePeak R/Bioconductor package. , 2014, Methods.

[23]  Helga Thorvaldsdóttir,et al.  Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration , 2012, Briefings Bioinform..

[24]  Timothy L. Bailey,et al.  Gene expression Advance Access publication May 4, 2011 DREME: motif discovery in transcription factor ChIP-seq data , 2011 .

[25]  Gideon Rechavi,et al.  Transcriptome-wide mapping of N6-methyladenosine by m6A-seq based on immunocapturing and massively parallel sequencing , 2013, Nature Protocols.

[26]  Marc D. Perry,et al.  ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia , 2012, Genome research.

[27]  Philip Machanick,et al.  MEME-ChIP: motif analysis of large DNA datasets , 2011, Bioinform..

[28]  Samie R. Jaffrey,et al.  The dynamic epitranscriptome: N6-methyladenosine and gene expression control , 2014, Nature Reviews Molecular Cell Biology.

[29]  P. Park.  ChIP–seq: advantages and challenges of a maturing technology , 2009, Nature Reviews Genetics.

[30]  Cole Trapnell,et al.  TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions , 2013, Genome Biology.

[31]  Miao Yu,et al.  A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation , 2013, Nature chemical biology.