{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,11,8]],"date-time":"2025-11-08T17:58:25Z","timestamp":1762624705630,"version":"3.41.2"},"reference-count":39,"publisher":"Oxford University Press (OUP)","issue":"1","license":[{"start":{"date-parts":[[2021,9,29]],"date-time":"2021-09-29T00:00:00Z","timestamp":1632873600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"funder":[{"DOI":"10.13039\/100000002","name":"National Institutes of Health","doi-asserted-by":"publisher","award":["R35 GM124998"],"award-info":[{"award-number":["R35 GM124998"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2022,1,17]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:sec>\n                  <jats:title>Motivation<\/jats:title>\n                  <jats:p>m6A methylation is a highly prevalent post-transcriptional modification in eukaryotes. MeRIP-seq or m6A-seq, which comprises immunoprecipitation of methylation fragments , is the most common method for measuring methylation signals. Existing computational tools for analyzing MeRIP-seq data sets and identifying differentially methylated genes\/regions are not most optimal. They either ignore the sparsity or dependence structure of the methylation signals within a gene\/region. Modeling the methylation signals using univariate distributions could also lead to high type I error rates and low sensitivity. In this paper, we propose using mean vector testing (MVT) procedures for testing differential methylation of RNA at the gene level. MVTs use a distribution-free test statistic with proven ability to control type I error even for extremely small sample sizes. We performed a comprehensive simulation study comparing the MVTs to existing MeRIP-seq data analysis tools. Comparative analysis of existing MeRIP-seq data sets is presented to illustrate the advantage of using MVTs.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Results<\/jats:title>\n                  <jats:p>Mean vector testing procedures are observed to control type I error rate and achieve high power for detecting differential RNA methylation using m6A-seq data. Results from two data sets indicate that the genes detected identified as having different m6A methylation patterns have high functional relevance to the study conditions.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Availability<\/jats:title>\n                  <jats:p>The dimer software package for differential RNA methylation analysis is freely available at https:\/\/github.com\/ouyang-lab\/DIMER.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Supplementary information<\/jats:title>\n                  <jats:p>Supplementary data are available at Briefings in Bioinformatics online.<\/jats:p>\n               <\/jats:sec>","DOI":"10.1093\/bib\/bbab309","type":"journal-article","created":{"date-parts":[[2021,9,21]],"date-time":"2021-09-21T11:23:21Z","timestamp":1632223401000},"source":"Crossref","is-referenced-by-count":5,"title":["Differential RNA methylation using multivariate statistical methods"],"prefix":"10.1093","volume":"23","author":[{"given":"Deepak Nag","family":"Ayyala","sequence":"first","affiliation":[{"name":"Division of Biostatistics and Data Science, Department of Population Health Sciences, Medical College of Georgia, Augusta University, Augusta, GA 30912, USA"}]},{"given":"Jianan","family":"Lin","sequence":"additional","affiliation":[{"name":"Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, MA 01003, USA"}]},{"given":"Zhengqing","family":"Ouyang","sequence":"additional","affiliation":[{"name":"The Jackson Laboratory for Genomic Medicine, Farmington CT, 06032, USA"}]}],"member":"286","published-online":{"date-parts":[[2021,9,29]]},"reference":[{"issue":"2","key":"2022011921010885300_ref1","first-page":"162","article-title":"Two generalizations of the binomial distribution","volume":"27","author":"Altham","year":"1978","journal-title":"J R Stat Soc Ser C Appl Stat"},{"issue":"10","key":"2022011921010885300_ref2","doi-asserted-by":"crossref","first-page":"R106","DOI":"10.1186\/gb-2010-11-10-r106","article-title":"Differential expression analysis for sequence count data","volume":"11","author":"Anders","year":"2010","journal-title":"Genome Biol"},{"issue":"6","key":"2022011921010885300_ref3","doi-asserted-by":"crossref","first-page":"926","DOI":"10.1093\/bib\/bbv089","article-title":"Statistical methods for detecting differentially methylated regions based on MethylCap-seq data","volume":"17","author":"Ayyala","year":"2016","journal-title":"Brief Bioinform"},{"key":"2022011921010885300_ref4","doi-asserted-by":"crossref","first-page":"195","DOI":"10.1093\/nar\/gkq1028","article-title":"The RNA modification database, RNAMDB: 2011 update","volume":"39","author":"Cantara","year":"2011","journal-title":"Nucleic Acids Res"},{"issue":"2","key":"2022011921010885300_ref5","doi-asserted-by":"crossref","first-page":"808","DOI":"10.1214\/09-AOS716","article-title":"A two-sample test for high-dimensional data with applications to gene-set testing","volume":"38","author":"Chen","year":"2010","journal-title":"Ann Stat"},{"key":"2022011921010885300_ref6","doi-asserted-by":"crossref","first-page":"26","DOI":"10.1016\/j.ab.2015.08.021","article-title":"iRNA-Methyl: identifying N6-methyladenosine sites using pseudo nucleotide composition","volume":"490","author":"Chen","year":"2015","journal-title":"Anal Biochem"},{"issue":"13859","key":"2022011921010885300_ref7","article-title":"Identification and analysis of the N6-methyladenosine in the Saccharomyces cerevisiae transcriptome","volume":"5","author":"Chen","year":"2015","journal-title":"Sci Rep"},{"issue":"33","key":"2022011921010885300_ref8","doi-asserted-by":"crossref","first-page":"20117","DOI":"10.1016\/S0021-9258(17)30477-5","article-title":"Sequence specificity of mRNA N6-Adenosine methyltransferase","volume":"265","author":"Csepany","year":"1990","journal-title":"J Biol Chem"},{"issue":"11","key":"2022011921010885300_ref9","doi-asserted-by":"crossref","first-page":"2622","DOI":"10.1016\/j.celrep.2017.02.059","article-title":"m6A RNA methylation regulates the self-renewal and tumorigenesis of glioblastoma stem cells","volume":"18","author":"Cui","year":"2017","journal-title":"Cell Rep"},{"key":"2022011921010885300_ref10","doi-asserted-by":"crossref","first-page":"S2","DOI":"10.1186\/1471-2164-16-S4-S2","article-title":"HEPeak: an HMM-based exome peak-finding package for RNA epigenome sequencing data","volume":"16","author":"Cui","year":"2015","journal-title":"BMC Genomics"},{"issue":"c","key":"2022011921010885300_ref11","first-page":"1","article-title":"MeTDiff: a novel differential RNA methylation analysis for MeRIP-Seq data","volume":"5963","author":"Cui","year":"2015","journal-title":"IEEE\/ACM Trans Comput Biol Bioinform"},{"issue":"24","key":"2022011921010885300_ref12","doi-asserted-by":"crossref","first-page":"3423","DOI":"10.1093\/bioinformatics\/btr539","article-title":"Pybedtools: a flexible Python library for manipulating genomic datasets and annotations","volume":"27","author":"Dale","year":"2011","journal-title":"Bioinformatics"},{"issue":"10","key":"2022011921010885300_ref13","doi-asserted-by":"crossref","first-page":"3971","DOI":"10.1073\/pnas.71.10.3971","article-title":"Identification of methylated nucleosides in messenger RNA from Novikoff hepatoma cells","volume":"71","author":"Desrosiers","year":"1974","journal-title":"Proc Natl Acad Sci U S A"},{"issue":"7397","key":"2022011921010885300_ref14","doi-asserted-by":"crossref","first-page":"201","DOI":"10.1038\/nature11112","article-title":"Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq","volume":"485","author":"Dominissini","year":"2012","journal-title":"Nature"},{"issue":"5","key":"2022011921010885300_ref15","doi-asserted-by":"crossref","first-page":"293","DOI":"10.1038\/nrg3724","article-title":"Gene expression regulation mediated through reversible m6A RNA methylation","volume":"15","author":"Fu","year":"2014","journal-title":"Nat Rev Genet"},{"issue":"6225","key":"2022011921010885300_ref16","doi-asserted-by":"crossref","first-page":"1002","DOI":"10.1126\/science.1261417","article-title":"m6A mRNA methylation facilitates resolution of na\u00efve pluripotency toward differentiation","volume":"347","author":"Geula","year":"2015","journal-title":"Science"},{"issue":"19","key":"2022011921010885300_ref17","doi-asserted-by":"crossref","first-page":"5735","DOI":"10.1093\/nar\/18.19.5735","article-title":"Sequence specificity of the human mRNA N6-adenosine methylase in vitro","volume":"18","author":"Harper","year":"1990","journal-title":"Nucleic Acids Res"},{"issue":"4","key":"2022011921010885300_ref18","doi-asserted-by":"crossref","first-page":"740","DOI":"10.1016\/j.cell.2013.10.028","article-title":"m6A mRNA methylation: a new circadian pacesetter","volume":"155","author":"Hastings","year":"2013","journal-title":"Cell"},{"issue":"1","key":"2022011921010885300_ref19","doi-asserted-by":"crossref","first-page":"8","DOI":"10.1016\/0003-9861(68)90553-5","article-title":"Methylated bases of ribosomal ribonucleic acid from HeLa cells","volume":"126","author":"Iwanami","year":"1968","journal-title":"Arch Biochem Biophys"},{"issue":"4","key":"2022011921010885300_ref20","doi-asserted-by":"crossref","first-page":"523","DOI":"10.1093\/bioinformatics\/btt703","article-title":"Causal analysis approaches in ingenuity pathway analysis","volume":"30","author":"Kr\u00e4mer","year":"2014","journal-title":"Bioinformatics"},{"issue":"5","key":"2022011921010885300_ref21","doi-asserted-by":"crossref","first-page":"666","DOI":"10.1016\/j.chom.2016.10.002","article-title":"Dynamics of human and viral RNA methylation during Zika virus infection","volume":"20","author":"Lichinchi","year":"2016","journal-title":"Cell Host Microbe"},{"issue":"8","key":"2022011921010885300_ref22","doi-asserted-by":"crossref","first-page":"767","DOI":"10.1038\/nmeth.3453","article-title":"Single-nucleotide-resolution mapping of m6A and m6Am throughout the transcriptome","volume":"12","author":"Bastian Linder","year":"2015","journal-title":"Nat Methods"},{"key":"2022011921010885300_ref23","doi-asserted-by":"crossref","DOI":"10.1038\/nchembio.1432","article-title":"A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation","volume":"10","author":"Liu","year":"2014","journal-title":"Nat Chem Biol"},{"key":"2022011921010885300_ref24","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1016\/j.ab.2016.01.014","article-title":"DRME: count-based differential RNA methylation analysis at small sample size scenario","volume":"499","author":"Liu","year":"2016","journal-title":"Anal Biochem"},{"issue":"1","key":"2022011921010885300_ref25","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1016\/0047-259X(77)90027-6","article-title":"Characterization of the multivariate poisson distribution","volume":"7","author":"Lukacs","year":"1977","journal-title":"J Multivar Anal"},{"issue":"11","key":"2022011921010885300_ref26","doi-asserted-by":"crossref","first-page":"e95","DOI":"10.1093\/nar\/gku351","article-title":"De novo detection of differentially bound regions for ChIP-seq data using peaks and windows: controlling error rates correctly","volume":"42","author":"Lun","year":"2014","journal-title":"Nucleic Acids Res"},{"issue":"12","key":"2022011921010885300_ref27","doi-asserted-by":"crossref","first-page":"1565","DOI":"10.1093\/bioinformatics\/btt171","article-title":"Exome-based analysis for RNA epigenome sequencing data","volume":"29","author":"Meng","year":"2013","journal-title":"Bioinformatics"},{"issue":"7","key":"2022011921010885300_ref28","doi-asserted-by":"crossref","first-page":"1635","DOI":"10.1016\/j.cell.2012.05.003","article-title":"Comprehensive analysis of mRNA methylation reveals enrichment in 3\u2019 UTRs and near stop codons","volume":"149","author":"Meyer","year":"2012","journal-title":"Cell"},{"issue":"10","key":"2022011921010885300_ref29","doi-asserted-by":"crossref","first-page":"951","DOI":"10.1056\/NEJMoa1600651","article-title":"Zika virus associated with microcephaly","volume":"374","author":"Mlakar","year":"2016","journal-title":"N Engl J Med"},{"issue":"5","key":"2022011921010885300_ref30","doi-asserted-by":"crossref","first-page":"929","DOI":"10.1016\/j.jspi.2012.11.001","article-title":"A test for the mean vector in large dimension and small samples","volume":"143","author":"Park","year":"2013","journal-title":"J Stat Plan Inf"},{"key":"2022011921010885300_ref31","doi-asserted-by":"crossref","first-page":"177","DOI":"10.1038\/cr.2014.3","article-title":"Mammalian WTAP is a regulatory subunit of the RNA N6-methyladenosine methyltransferase","volume":"24","author":"Ping","year":"2014","journal-title":"Cell Res"},{"issue":"6","key":"2022011921010885300_ref32","doi-asserted-by":"crossref","first-page":"841","DOI":"10.1093\/bioinformatics\/btq033","article-title":"BEDTools: a flexible suite of utilities for comparing genomic features","volume":"26","author":"Quinlan","year":"2010","journal-title":"Bioinformatics"},{"issue":"2","key":"2022011921010885300_ref33","doi-asserted-by":"crossref","first-page":"264","DOI":"10.1016\/0005-2787(69)90078-1","article-title":"Isolation and characterization of N6-methyladenosine from Escherichia coli valine transfer RNA","volume":"190","author":"Saneyoshi","year":"1969","journal-title":"Biochim Biophys Acta Nucleic Acids Protein Synth"},{"key":"2022011921010885300_ref34","doi-asserted-by":"crossref","first-page":"349","DOI":"10.1016\/j.jmva.2012.08.014","article-title":"A two sample test in high dimensional data","volume":"114","author":"Srivastava","year":"2013","journal-title":"J Multivar Anal"},{"issue":"2","key":"2022011921010885300_ref35","doi-asserted-by":"crossref","first-page":"256","DOI":"10.1016\/0047-259X(90)90084-U","article-title":"Some representations of the multivariate Bernoulli and binomial distributions","volume":"32","author":"Teugels","year":"1990","journal-title":"J Multivar Anal"},{"key":"2022011921010885300_ref36","doi-asserted-by":"crossref","first-page":"1502","DOI":"10.1038\/npp.2017.31","article-title":"The role of the RNA demethylase FTO (fat mass and obesity-associated) and mRNA methylation in hippocampal memory formation","volume":"42","author":"Walters","year":"2017","journal-title":"Neuropsychopharmacology"},{"issue":"9","key":"2022011921010885300_ref37","doi-asserted-by":"crossref","first-page":"R137","DOI":"10.1186\/gb-2008-9-9-r137","article-title":"Model-based analysis of ChIP-Seq (MACS)","volume":"9","author":"Zhang","year":"2008","journal-title":"Genome Biol"},{"issue":"5","key":"2022011921010885300_ref38","doi-asserted-by":"crossref","first-page":"1278","DOI":"10.1105\/tpc.108.058883","article-title":"MTA is an Arabidopsis messenger RNA adenosine methylase and interacts with a homolog of a sex-specific splicing factor","volume":"20","author":"Zhong","year":"2008","journal-title":"Plant Cell Online"},{"issue":"10","key":"2022011921010885300_ref39","doi-asserted-by":"crossref","first-page":"e91","DOI":"10.1093\/nar\/gkw104","article-title":"SRAMP: prediction of mammalian N6-methyladenosine (m6A) sites based on sequence-derived features","volume":"44","author":"Zhou","year":"2016","journal-title":"Nucleic Acids Res"}],"container-title":["Briefings in Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bib\/article-pdf\/23\/1\/bbab309\/42230040\/bbab309.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bib\/article-pdf\/23\/1\/bbab309\/42230040\/bbab309.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2022,1,19]],"date-time":"2022-01-19T21:01:42Z","timestamp":1642626102000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bib\/article\/doi\/10.1093\/bib\/bbab309\/6377721"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2021,9,29]]},"references-count":39,"journal-issue":{"issue":"1","published-print":{"date-parts":[[2022,1,17]]}},"URL":"https:\/\/doi.org\/10.1093\/bib\/bbab309","relation":{},"ISSN":["1467-5463","1477-4054"],"issn-type":[{"type":"print","value":"1467-5463"},{"type":"electronic","value":"1477-4054"}],"subject":[],"published-other":{"date-parts":[[2022,1]]},"published":{"date-parts":[[2021,9,29]]},"article-number":"bbab309"}}