We apply for the first time genome-wide expression–methylation quantitative trait loci (emQTL) analysis between DNA methylation and gene expression, and discover that DNA methylation at enhancers and ERα, FOXA1, and GATA3 binding regions is a breast cancer subtype-specific phenotypic feature. This evidence concerns the gene ESR1 and breast cancer.