OxyBS: estimation of 5-methylcytosine and 5-hydroxymethylcytosine from tandem-treated oxidative bisulfite and bisulfite DNA

The use of sodium bisulfite (BS) treatment followed by hybridization to an Illumina Infinium BeadChip (HumanMethylation450 and MethylationEPIC) is a common method for interrogating 5-methylcytosine (5mC) at single nucleotide resolution. However, standard treatment of DNA with BS does not allow disambiguation of 5mC from an additional cytosine modification, 5-hydroxymethylcytosine (5hmC). Recently, it has been demonstrated that paired BS and oxidative bisulfite (oxBS) treatment on the same sample followed by hybridization to an Infinium microarray permits the differentiation of 5hmC from 5mC. Nevertheless, estimation of 5hmC and 5mC from tandem-treated arrays has been shown to produce irregular estimates of cytosine modifications. We present a novel method using maximum likelihood estimation to accurately estimate the parameters of unmethylated cytosine (5C), 5mC and 5hmC from Infinium microarray data given the signal intensities from the oxBS and BS replicates. OxyBS is an R package available on CRAN. Andres.Houseman@oregonstate.edu Supplementary data are available at Bioinformatics online.