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.

MOTIVATION: Estimating a phenotype distribution conditional on a set of discrete-valued predictors is a commonly encountered task. For example, interest may be in how the density of a quantitative trait varies with single nucleotide polymorphisms and patient characteristics. The subset of important predictors is not usually known in advance. This becomes more challenging with a high-dimensional predictor set when there is the possibility of interaction. RESULTS: We demonstrate a novel non-parametric Bayes method based on a tensor factorization of predictor-dependent weights for Gaussian kernels. The method uses multistage predictor selection for dimension reduction, providing succinct models for the phenotype distribution. The resulting conditional density morphs flexibly with the selected predictors. In a simulation study and an application to molecular epidemiology data, we demonstrate advantages over commonly used methods.

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/

SUMMARY: Proteomics technology has shown promise in identifying biomarkers for disease, toxicant exposure and stress. We show by example that the genetic algorithm/k-nearest neighbors method, developed for mining high-dimensional microarray gene expression data, is also capable of mining surface enhanced laser desorption/ionization-time-of-flight proteomics data. AVAILABILITY: The source code of the program and documentation on how to use it are freely available to non-commercial users at http://dir.niehs.nih.gov/dirbb/lifiles/softlic.htm