Resistance to endocrine therapy remains a major clinical problem in breast cancer. Genetic studies highlight the potential role of estrogen receptor-α (ESR1) mutations, which show increased prevalence in the metastatic, endocrine-resistant setting. No naturally occurring ESR1 mutations have been reported in in vitro models of BC either before or after the acquisition of endocrine resistance making functional consequences difficult to study. We report the first discovery of naturally occurring ESR1 Y537C and ESR1 Y537S mutations in MCF7 and SUM44 ESR1-positive cell lines after acquisition of resistance to long-term-estrogen-deprivation (LTED) and subsequent resistance to fulvestrant (ICIR). Mutations were enriched with time, impacted on ESR1 binding to the genome and altered the ESR1 interactome. The results highlight the importance and functional consequence of these mutations and provide an important resource for studying endocrine resistance.

BackgroundTherapies targeting estrogenic stimulation in estrogen receptor-positive (ER+) breast cancer (BC) reduce mortality, but resistance remains a major clinical problem. Molecular studies have shown few high-frequency mutations to be associated with endocrine resistance. In contrast, expression profiling of primary ER+ BC samples has identified several promising signatures/networks for targeting.MethodsTo identify common adaptive mechanisms associated with resistance to aromatase inhibitors (AIs), we assessed changes in global gene expression during adaptation to long-term estrogen deprivation (LTED) in a panel of ER+ BC cell lines cultured in 2D on plastic (MCF7, T47D, HCC1428, SUM44 and ZR75.1) or in 3D on collagen (MCF7) to model the stromal compartment. Furthermore, dimethyl labelling followed by LC-MS/MS was used to assess global changes in protein abundance. The role of target genes/proteins on proliferation, ER-mediated transcription and recruitment of ER to target gene promoters was analysed.ResultsThe cholesterol biosynthesis pathway was the common upregulated pathway in the ER+ LTED but not the ER– LTED cell lines, suggesting a potential mechanism dependent on continued ER expression. Targeting the individual genes of the cholesterol biosynthesis pathway with siRNAs caused a 30–50 % drop in proliferation. Further analysis showed increased expression of 25-hydroxycholesterol (HC) in the MCF7 LTED cells. Exogenous 25-HC or 27-HC increased ER-mediated transcription and expression of the endogenous estrogen-regulated gene TFF1 in ER+ LTED cells but not in the ER– LTED cells. Additionally, recruitment of the ER and CREB-binding protein (CBP) to the TFF1 and GREB1 promoters was increased upon treatment with 25-HC and 27-HC. In-silico analysis of two independent studies of primary ER+ BC patients treated with neoadjuvant AIs showed that increased expression of MSMO1, EBP, LBR and SQLE enzymes, required for cholesterol synthesis and increased in our in-vitro models, was significantly associated with poor response to endocrine therapy.ConclusionTaken together, these data provide support for the role of cholesterol biosynthesis enzymes and the cholesterol metabolites, 25-HC and 27-HC, in a novel mechanism of resistance to endocrine therapy in ER+ BC that has potential as a therapeutic target.Electronic supplementary materialThe online version of this article (doi:10.1186/s13058-016-0713-5) contains supplementary material, which is available to authorized users.

Genome-wide association studies (GWAS) have identified approximately 100 breast cancer risk loci. Translating these findings into a greater understanding of the mechanisms that influence disease risk requires identification of the genes or non-coding RNAs that mediate these associations. Here, we use Capture Hi-C (CHi-C) to annotate 63 loci; we identify 110 putative target genes at 33 loci. To assess the support for these target genes in other data sources we test for associations between levels of expression and SNP genotype (eQTLs), disease-specific survival (DSS), and compare them with somatically mutated cancer genes. 22 putative target genes are eQTLs, 32 are associated with DSS and 14 are somatically mutated in breast, or other, cancers. Identifying the target genes at GWAS risk loci will lead to a greater understanding of the mechanisms that influence breast cancer risk and prognosis.