For studying breast epithelia differentiation, DMRs were called from WGBS, requiring at least 5 aligned reads to call differentially-methylated CpG, and at least 3 differentiallymethylated CpGs within a distance of 200 bp of each other44. For studying tissue environment vs. developmental origin, DMRs were called from MeDIP and MRE data using the M&M algorithm55.