MOTIVATION: Gene Set Enrichment Analysis (GSEA) is widely used to interpret DNA methylation data by associating differentially methylated sites with biological pathways. However, existing GSEA methods struggle with several challenges in methylation data, including probe dependency, probe number bias, and the complexity of gene-probe mapping. These limitations can lead to biased enrichment results, reduced statistical power, and computational inefficiencies. RESULTS: We introduce gsGene and gsPG, two novel GSEA methods specifically designed for DNA methylation data. gsGene aggregates association signals at gene level while correcting for probe dependency and probe number bias, enabling more biologically meaningful enrichment analysis. gsPG takes a different approach by conducting gene set enrichment using summary statistics for independent probe groups based on gene annotation, mitigating biases from multi-mapping probes. Both methods improve computational efficiency, enhance statistical power, and effectively control type I error rates. Comprehensive evaluations in two large datasets demonstrate superior performance compared to existing methods. Furthermore, we propose a novel beta distribution fitting strategy to improve enrichment P-value estimation, providing a computationally efficient alternative to traditional permutation-based gene set methods. AVAILABILITY AND IMPLEMENTATION: These methods are implemented in the R package dmGsea, which is freely available on GitHub and Bioconductor (DOI: 10.18129/B9.bioc.dmGsea). The package supports Illumina 450K, EPIC, and mouse methylation arrays and can be extended to other omics data with user-provided probe-to-gene mapping annotations.

MOTIVATION: DNA methylation-based predictors of various biological metrics have been widely published and are becoming valuable tools in epidemiologic studies of epigenetics and personalized medicine. However, generating these predictors from original source software and web servers is complex and time consuming. Furthermore, different predictors were often derived based on data from different types of arrays, where array differences and batch effects can make predictors difficult to compare across studies. RESULTS: We integrate these published methods into a single R function to produce 158 previously published predictors for chronological age, biological age, exposures, lifestyle traits and serum protein levels using both classical and principal component-based methods. To mitigate batch and array differences, we also provide a modified RCP method (ref-RCP) that normalize input DNA methylation data to reference data prior to estimation. Evaluations in real datasets show that this approach improves estimate precision and comparability across studies. AVAILABILITY AND IMPLEMENTATION: The function was included in software package ENmix, and is freely available from Bioconductor website (https://www.bioconductor.org/packages/release/bioc/html/ENmix.html).

SUMMARY: ipDMR is an R software tool for identification of differentially methylated regions (DMRs) using auto-correlated P-values for individual CpGs from epigenome-wide association analysis using array or bisulfite sequencing data. It summarizes P-values for adjacent CpGs, identifies association peaks and then extends peaks to find boundaries of DMRs. ipDMR uses BED format files as input and is easy to use. Simulations guided by real data found that ipDMR outperformed current available methods and provided slightly higher true positive rates and much lower false discovery rates. AVAILABILITY AND IMPLEMENTATION: ipDMR is available at https://bioconductor.org/packages/release/bioc/html/ENmix.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MOTIVATION: The Illumina HumanMethylation450 BeadChip has been extensively utilized in epigenome-wide association studies. This array and its successor, the MethylationEPIC array, use two types of probes-Infinium I (type I) and Infinium II (type II)-in order to increase genome coverage but differences in probe chemistries result in different type I and II distributions of methylation values. Ignoring the difference in distributions between the two probe types may bias downstream analysis. RESULTS: Here, we developed a novel method, called Regression on Correlated Probes (RCP), which uses the existing correlation between pairs of nearby type I and II probes to adjust the beta values of all type II probes. We evaluate the effect of this adjustment on reducing probe design type bias, reducing technical variation in duplicate samples, improving accuracy of measurements against known standards, and retention of biological signal. We find that RCP is statistically significantly better than unadjusted data or adjustment with alternative methods including SWAN and BMIQ. AVAILABILITY: We incorporated the method into the R package ENmix, which is freely available from the Bioconductor website (https://www.bioconductor.org/packages/release/bioc/html/ENmix.html). CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MOTIVATION: 5-Methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) are important epigenetic regulators of gene expression. 5mC and 5hmC levels can be computationally inferred at single base resolution using sequencing or array data from paired DNA samples that have undergone bisulfite and oxidative bisulfite conversion. Current estimation methods have been shown to produce irregular estimates of 5hmC level or are extremely computation intensive. RESULTS: We developed an efficient method oxBS-MLE based on binomial modeling of paired bisulfite and oxidative bisulfite data from sequencing or array analysis. Evaluation in several datasets showed that it outperformed alternative methods in estimate accuracy and computation speed. AVAILABILITY AND IMPLEMENTATION: oxBS-MLE is implemented in Bioconductor package ENmix. CONTACT: [email protected] information: Supplementary data are available at Bioinformatics online.

UNLABELLED: Genetic association studies increasingly rely on the use of linkage disequilibrium (LD) tag SNPs to reduce genotyping costs. We developed a software package TAGster to select, evaluate and visualize LD tag SNPs both for single and multiple populations. We implement several strategies to improve the efficiency of current LD tag SNP selection algorithms: (1) we modify the tag SNP selection procedure of Carlson et al. to improve selection efficiency and further generalize it to multiple populations. (2) We propose a redundant SNP elimination step to speed up the exhaustive tag SNP search algorithm proposed by Qin et al. (3) We present an additional multiple population tag SNP selection algorithm based on the framework of Howie et al., but using our modified exhaustive search procedure. We evaluate these methods using resequenced candidate gene data from the Environmental Genome Project and show improvements in both computational and tagging efficiency. AVAILABILITY: The software Package TAGster is freely available at http://www.niehs.nih.gov/research/resources/software/tagster/