Aberrant DNA methylation is an early event in breast carcinogenesis and plays a critical role in regulating gene expression. Here, we perform genome-wide expression-methylation Quantitative Trait Loci (emQTL) analysis through the integration of DNA methylation and gene expression to identify disease-driving pathways under epigenetic control. By grouping the emQTLs using biclustering we identify associations representing important biological processes associated with breast cancer pathogenesis including regulation of proliferation and tumor-infiltrating fibroblasts. We report genome-wide loss of enhancer methylation at binding sites of proliferation-driving transcription factors including CEBP-β, FOSL1, and FOSL2 with concomitant high expression of proliferation-related genes in aggressive breast tumors as we confirm with scRNA-seq. The identified emQTL-CpGs and genes were found connected through chromatin loops, indicating that proliferation in breast tumors is under epigenetic regulation by DNA methylation. Interestingly, the associations between enhancer methylation and proliferation-related gene expression were also observed within known subtypes of breast cancer, suggesting a common role of epigenetic regulation of proliferation. Taken together, we show that proliferation in breast cancer is linked to loss of methylation at specific enhancers and transcription factor binding and gene activation through chromatin looping.

Background Abnormal DNA methylation is observed as an early event in breast carcinogenesis. However, how such alterations arise is still poorly understood. microRNAs (miRNAs) regulate gene expression at the post-transcriptional level and play key roles in various biological processes. Here, we integrate miRNA expression and DNA methylation at CpGs to study how miRNAs may affect the breast cancer methylome and how DNA methylation may regulate miRNA expression. Methods miRNA expression and DNA methylation data from two breast cancer cohorts, Oslo2 (n = 297) and The Cancer Genome Atlas (n = 439), were integrated through a correlation approach that we term miRNA-methylation Quantitative Trait Loci (mimQTL) analysis. Hierarchical clustering was used to identify clusters of miRNAs and CpGs that were further characterized through analysis of mRNA/protein expression, clinicopathological features, in silico deconvolution, chromatin state and accessibility, transcription factor binding, and long-range interaction data. Results Clustering of the significant mimQTLs identified distinct groups of miRNAs and CpGs that reflect important biological processes associated with breast cancer pathogenesis. Notably, two major miRNA clusters were related to immune or fibroblast infiltration, hence identifying miRNAs associated with cells of the tumor microenvironment, while another large cluster was related to estrogen receptor (ER) signaling. Studying the chromatin landscape surrounding CpGs associated with the estrogen signaling cluster, we found that miRNAs from this cluster are likely to be regulated through DNA methylation of enhancers bound by FOXA1, GATA2, and ER-alpha. Further, at the hub of the estrogen cluster, we identified hsa-miR-29c-5p as negatively correlated with the mRNA and protein expression of DNA methyltransferase DNMT3A, a key enzyme regulating DNA methylation. We found deregulation of hsa-miR-29c-5p already present in pre-invasive breast lesions and postulate that hsa-miR-29c-5p may trigger early event abnormal DNA methylation in ER-positive breast cancer. Conclusions We describe how miRNA expression and DNA methylation interact and associate with distinct breast cancer phenotypes.

Background:Abnormal DNA methylation is observed as an early event in breast carcinogenesis. However, how such alterations arise is still poorly understood. microRNAs (miRNAs) regulate gene expression at the post-transcriptional level and have been shown to play key roles in various biological processes. Here, we integrate miRNA expression and DNA methylation at CpGs to study how miRNAs may affect the breast cancer methylome and how DNA methylation may regulate miRNA expression. ResultsmiRNA expression and DNA methylation data from two breast cancer cohorts were subjected to genome-wide correlation analysis. Clustering of the miRNA expression-DNA methylation association pairs significant in both cohorts identified distinct clusters of miRNAs and CpGs. These clusters recapitulated important biological processes associated with breast cancer pathogenesis. Notably, two major clusters were related to immune or fibroblast infiltration, hence identifying miRNAs associated with cells of the tumor microenvironment, while another large cluster was related to estrogen receptor (ER) signaling. Studying the chromatin landscape surrounding the CpGs associated with the estrogensignaling cluster, we found that miRNAs from this cluster are likely to be regulated through DNA methylation of enhancers bound by FOXA1, GATA2 and ER-alpha. Further, at the hub of the estrogen-cluster, we identified hsa-miR-29c-5p as negatively correlated with the mRNA and protein expression of the DNA methyltransferase DNMT3A, a key enzyme regulating DNA methylation. We found deregulation of hsa-miR-29c-5p already in pre-invasive breast lesions and postulate that hsa-miR-29c-5p may trigger early event abnormal DNA methylation in ER positive breast cancer. ConclusionsWe describe how miRNA expression and DNA methylation interact and associate with distinct breast cancer phenotypes.

Breast cancer is a highly heterogeneous disease driven by multiple factors including genetic and epigenetic alterations. DNA methylation patterns have been shown to be altered on a genome-wide scale and previous studies have highlighted the critical role of aberrant DNA methylation on gene expression and breast cancer pathogenesis. Here, we perform genome-wide expression-methylation Quantitative Trait Loci (emQTL), a method for integration of CpG methylation and gene expression to identify disease-driving genes under epigenetic control. By grouping these emQTLs by biclustering we identify associations representing important biological processes associated with breast cancer pathogenesis such as proliferation and tumor infiltrating fibroblasts. We report hypomethylation at enhancers carrying transcription factor binding sites of key proliferation-driving transcription factors such as CEBP-β, FOSL1, and FOSL2, with concomitant high expression of cell cycle- and proliferation-related genes in aggressive breast tumors. The identified CpGs and genes were found to be connected through chromatin loops, together indicating that proliferation in aggressive breast tumors is under epigenetic regulation by DNA methylation. Interestingly, there was a significant correlation between proliferation-related DNA methylation and gene expression also within subtypes of breast cancer, thereby showing that varying proliferation may be explained by epigenetic profiles across breast cancer subtypes. Indeed, the identified proliferation gene signature was prognostic both in the Luminal A and Luminal B subtypes. Taken together, we show that proliferation in breast cancer is linked to hypomethylation at specific enhancers and transcription factor binding mediated through chromatin loops.