Abstract
Estrogen Receptor alpha (ERα) is the main driver and prime drug target in luminal breast. ERα chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ERα chromatin action, along with its biological implications. Here, we use a large set of ERα ChIP-seq data from 70 ERα+ breast cancers to explore inter-patient heterogeneity in ERα DNA binding, to reveal a striking inter-tumor heterogeneity of ERα action. Interestingly, commonly-shared ERα sites showed the highest estrogen-driven enhancer activity and were most-engaged in long-range chromatin interactions. In addition, the most-commonly shared ERα-occupied enhancers were enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ERα and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we could confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ERα-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ERα landscape, with the most-common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.