ReMeDy: DNAメチル化調節異常の領域ベース検出のための柔軟な統計的フレームワーク

Region-based epigenome-wide association studies have demonstrated improved statistical power and bio-logical interpretability compared with probe-wise analyses of DNA methylation data. However, most existing region-based methods characterize methylation dysregulation primarily through changes in mean methylation levels associated with a phenotype of interest. Substantial evidence indicates that phenotype-associated methylation alterations may also manifest through changes in methylation variability or through joint shifts in mean and variability. Despite this, no existing statistical framework jointly models mean-variance methylation changes in a region-based manner. We propose ReMeDy, a flexible statistical framework that uses a hierarchical likelihood approach within a generalized linear model setting to identify differentially methylated regions, variably methylated regions, and regions exhibiting joint differential and variable methylation at a genome-wide scale. Unlike existing models, ReMeDy operates directly on biologically defined co-methylated regions, allowing it to naturally capture spatial correlation inherent in DNA methylation array data, while avoiding reliance on heuristic, user-defined tuning parameters such as smoothing spans and kernel bandwidths that can substantially influence results and introduce subjectivity. Through extensive simulation studies and comprehensive benchmarking against popular models, we demonstrate that ReMeDy maintains false discovery and type-I error rates at nominal levels while achieving consistently higher statistical power across a wide range of realistic scenarios. Application to population-level DNA methylation data further shows that ReMeDy identifies biologically meaningful regions and pathways implicated in complex human diseases that are not captured by conventional mean-based analyses alone. ReMeDy is implemented as an open-source R package and is freely available at https://github.com/SChatLab/ReMeDy .