for Illumina HumanMethylation450 (HM450) BeadChip analysis. The total amount of DNA from each breast milk sample was bisulfite
treated with the Zymo EZ DNA methylation kit (Zymo Research, Irvine, CA) and 1µl aliquots were used for MethyLight quality
control (QC) analyses to determine the completeness of conversion and the amount of converted DNA available for the HM450 assay
[26 (link)]. The remaining bisulfite converted DNA samples were further processed using the
Illumina FFPE Restoration Solution (Illumina, San Diego, CA) as specified by the manufacturer. The Restoration Solution repairs
degraded DNAs for use in genome-scale genotyping and DNA methylation assay platforms. The entire restored sample was then used as
a substrate for the Illumina HM450 BeadArrays, as recommended by the manufacturer and described previously [27 (link)]. BeadArrays were scanned using Illumina iScan readers and the raw signal intensities were extracted
from the *.IDAT files and normalized using the R package sesame [28 (link), 29 (link)], a recently developed R package that masks problematic probes (i.e. probes for which DNA
methylation is invalid because they overlap SNPs or repeats).