Invasive lobular carcinoma (ILC) of the breast typically presents with clinical biomarkers consistent with a favorable response to endocrine therapies, and over 90 % of ILC cases express the estrogen receptor (ER). However, a subset of ILC cases may be resistant to endocrine therapies, suggesting that ER biology is unique in ILC. Using ILC cell lines, we previously demonstrated that ER regulates a distinct gene expression program in ILC cells, and we hypothesized that these ER-driven pathways modulate the endocrine response in ILC. One potential novel pathway is via the Wnt ligand WNT4, a critical signaling molecule in mammary gland development regulated by the progesterone receptor.

Scaffold attachment factors SAFB1 and SAFB2 are multifunctional proteins that share >70% sequence similarity. SAFB1-knockout (SAFB1(-/-)) mice display a high degree of lethality, severe growth retardation, and infertility in male mice. To assess the in vivo role of SAFB2, and to identify unique functions of the two paralogs, we generated SAFB2(-/-) mice. In stark contrast to SAFB1(-/-), SAFB2(-/-) offspring were born at expected Mendelian ratios and did not show any obvious defects in growth or fertility. Generation of paralog-specific antibodies allowed extensive expression analysis of SAFB1 and SAFB2 in mouse tissues, showing high expression of both SAFB1 and SAFB2 in the immune system, and in hormonally controlled tissues, with especially high expression of SAFB2 in the male reproductive tract. Further analysis showed a significantly increased testis weight in SAFB2(-/-) mice, which was associated with an increased number of Sertoli cells. Our data suggest that this is at least in part caused by alterations in androgen-receptor function and expression upon deletion of SAFB2. Thus, despite a high degree of sequence similarity, SAFB1(-/-) and SAFB2(-/-) mice do not totally phenocopy each other. SAFB2(-/-) mice are viable, and do not show any major defects, and our data suggest a role for SAFB2 in the differentiation and activity of Sertoli cells that deserves further study.

Estrogen receptor (ER) activity is critical for the development and progression of the majority of breast cancers. It is known that ER is differentially bound to DNA leading to transcriptomic and phenotypic changes in different breast cancer models. We investigated whether single nucleotide variants (SNVs) in ER binding sites (regSNVs) contribute to ER action through changes in the ER cistrome, thereby affecting disease progression. Here we developed a computational pipeline to identify SNVs in ER binding sites using chromatin immunoprecipitation sequencing (ChIP-seq) data from ER+ breast cancer models.

Werner syndrome (WS) is a rare disorder characterized by the premature onset of several pathologies associated with aging. The gene responsible for WS codes for a RecQ-type DNA helicase and is believed to be involved in different aspects of DNA repair, replication, and transcription. We recently identified the Scaffold attachment factor B1 (SAFB1) as a potential interactants in human cells. SAFB1 is a multifunctional protein that binds both nucleic acids and is involved in the attachment of chromatin to the nuclear matrix, transcription, and stress response. Mice lacking SAFB1 exhibit developmental abnormalities in their lungs, high incidence of perinatal lethality, and adults develop different types of tumors. Mouse embryonic fibroblasts from Safb1-null animals are immortalized in culture. In this study, mice with a mutation in the helicase domain of the Wrn gene were crossed to Safb1-null mice. Double homozygous mutant mice exhibited increased apoptosis, a lower cell proliferation rate in their lungs and a higher incidence of perinatal death compared to Safb1-null mice. Few double homozygous mutants survived weaning and died before the age of six months. Finally, mouse embryonic fibroblasts lacking a functional Wrn helicase inhibited the immortalization of Safb1-null cells. These results indicate that an intact Wrn protein is required for immortalization and tumorigenesis in Safb1-null mice.

Elevated levels of estrogen are a risk factor for breast cancer. In addition to inducing DNA damage, estrogens can enhance cell proliferation as well as modulate fatty acid metabolism that collectively contributes to mammary tumorigenesis. Sulforaphane (SFN) is an isothiocyanate derived from broccoli that is currently under evaluation in multiple clinical trials for prevention of several diseases, including cancer. Previous studies showed that SFN suppressed DNA damage and lipogenesis pathways. Therefore, we hypothesized that administering SFN to animals that are co-exposed to 17β-estradiol (E2) would prevent mammary tumor formation. In our study, 4-6 week old female August Copenhagen Irish rats were implanted with slow-release E2 pellets (3 mg x 3 times) and gavaged 3x/week with either vehicle or 100 μmol/kg SFN for 56 weeks. SFN-treated rats were protected significantly against mammary tumor formation compared to vehicle controls. Mammary glands of SFN-treated rats showed decreased DNA damage while serum free fatty acids and triglyceride species were 1.5 to 2-fold lower in SFN-treated rats. Further characterization also showed that SFN diminished expression of enzymes involved in mammary gland lipogenesis. This study indicated that SFN protects against breast cancer development through multiple potential mechanisms in a clinically relevant hormonal carcinogenesis model.

Mutations in ESR1 that confer constitutive estrogen receptor alpha (ER) activity in the absence of ligand are acquired by ≥40% of metastatic breast cancers (MBC) resistant to adjuvant aromatase inhibitor (AI) therapy. To identify targetable vulnerabilities in MBC, we examined steroid hormone receptors and tumor-infiltrating immune cells in metastatic lesions with or without ER mutations. ER and progesterone receptor (PR) were significantly lower in metastases with wild-type (WT) ER compared with those with mutant ER, suggesting that metastases that evade AI therapy by mechanism(s) other than acquiring ER mutations lose dependency on ER and PR. Metastases with mutant ER had significantly higher T regulatory and Th cells, total macrophages, and programmed death ligand-1 (PD-L1)-positive immune-suppressive macrophages than those with WT ER. Breast cancer cells with CRISPR-Cas9-edited ER (D538G, Y537S, or WT) and patient-derived xenografts harboring mutant or WT ER revealed genes and proteins elevated in mutant ER cells, including androgen receptor (AR), chitinase-3-like protein 1 (CHI3L1), and IFN-stimulated genes (ISG). Targeting these proteins blunted the selective advantage of ER-mutant tumor cells to survive estrogen deprivation, anchorage independence, and invasion. Thus, patients with mutant ER MBC might respond to standard-of-care fulvestrant or other selective ER degraders when combined with AR or CHI3L1 inhibition, perhaps with the addition of immunotherapy. SIGNIFICANCE: Targetable alterations in MBC, including AR, CHI3L1, and ISG, arise following estrogen-deprivation, and ER-mutant metastases may respond to immunotherapies due to elevated PD-L1+ macrophages.See related article by Arnesen et al., p. 539.

Metastatic breast cancer is a major cause of cancer-related deaths in woman. Brain metastasis is a common and devastating site of relapse for several breast cancer molecular subtypes, including oestrogen receptor-positive disease, with life expectancy of less than a year. While efforts have been devoted to developing therapeutics for extra-cranial metastasis, drug penetration of blood-brain barrier (BBB) remains a major clinical challenge. Defining molecular alterations in breast cancer brain metastasis enables the identification of novel actionable targets.

Constitutively active estrogen receptor α (ER/ESR1) mutations have been identified in approximately one-third of ER+ metastatic breast cancers. Although these mutations are known as mediators of endocrine resistance, their potential role in promoting metastatic disease has not yet been mechanistically addressed. In this study, we show the presence of ESR1 mutations exclusively in distant but not local recurrences in five independent breast cancer cohorts. In concordance with transcriptomic profiling of ESR1-mutant tumors, genome-edited ESR1 Y537S and D538G-mutant cell models exhibited a reprogrammed cell adhesive gene network via alterations in desmosome/gap junction genes and the TIMP3/MMP axis, which functionally conferred enhanced cell-cell contacts while decreasing cell-extracellular matrix adhesion. In vivo studies showed ESR1-mutant cells were associated with larger multicellular circulating tumor cell (CTC) clusters with increased compactness compared with ESR1 wild-type CTCs. These preclinical findings translated to clinical observations, where CTC clusters were enriched in patients with ESR1-mutated metastatic breast cancer. Conversely, context-dependent migratory phenotypes revealed cotargeting of Wnt and ER as a vulnerability in a D538G cell model. Mechanistically, mutant ESR1 exhibited noncanonical regulation of several metastatic pathways, including secondary transcriptional regulation and de novo FOXA1-driven chromatin remodeling. Collectively, these data provide evidence for ESR1 mutation-modulated metastasis and suggest future therapeutic strategies for targeting ESR1-mutant breast cancer.

Invasive lobular carcinoma (ILC) is a histological subtype of breast cancer with distinct molecular and clinical features from the more common subtype invasive ductal carcinoma (IDC). ILC cells exhibit anchorage-independent growth in ultra-low attachment (ULA) suspension cultures, which is largely attributed to the loss of E-cadherin. In addition to anoikis resistance, herein we show that human ILC cell lines exhibit enhanced cell proliferation in ULA cultures as compared to IDC cells. Proteomic comparison of ILC and IDC cell lines identified induction of PI3K/Akt and p90-RSK pathways specifically in ULA culture in ILC cells. Further transcriptional profiling uncovered unique upregulation of the inhibitors of differentiation family transcription factors ID1 and ID3 in ILC ULA culture, the knockdown of which diminished the anchorage-independent growth of ILC cell lines through cell cycle arrest. We find that ID1 and ID3 expression is higher in human ILC tumors as compared to IDC, correlated with worse prognosis uniquely in patients with ILC and associated with upregulation of angiogenesis and matrisome-related genes. Altogether, our comprehensive study of anchorage independence in human ILC cell lines provides mechanistic insights and clinical implications for metastatic dissemination of ILC and implicates ID1 and ID3 as novel drivers and therapeutic targets for lobular breast cancer.

Breast cancer is categorized by the molecular and histologic presentation of the tumor, with the major histologic subtypes being No Special Type (NST) and Invasive Lobular Carcinoma (ILC). ILC are characterized by growth in a single file discohesive manner with stromal infiltration attributed to their hallmark pathognomonic loss of E-cadherin ( CDH1 ). Few ILC cell line models are available to researchers. Here we report the successful establishment and characterization of a novel ILC cell line, WCRC-25, from a metastatic pleural effusion from a postmenopausal Caucasian woman with metastatic ILC. WCRC-25 is an ER-negative luminal epithelial ILC cell line with both luminal and Her2-like features. It exhibits anchorage independent growth and haptotactic migration towards Collagen I. Sequencing revealed a CDH1 Q706* truncating mutation, together with mutations in FOXA1, CTCF, BRCA2 and TP53 , which were also seen in a series of metastatic lesions from the patient. Copy number analyses revealed amplification and deletion of genes frequently altered in ILC while optical genome mapping revealed novel structural rearrangements. RNA-seq analysis comparing the primary tumor, metastases and the cell line revealed signatures for cell cycle progression and receptor tyrosine kinase signaling. To assess targetability, we treated WCRC-25 with AZD5363 and Alpelisib confirming WCRC-25 as susceptible to PI3K/AKT signaling inhibition as predicted by our RNA sequencing analysis. In conclusion, we report WCRC-25 as a novel ILC cell line with promise as a valuable research tool to advance our understanding of ILC and its therapeutic vulnerabilities.

Invasive lobular carcinoma (ILC) is the second most common histological subtype of breast cancer following invasive ductal carcinoma (IDC). ILC differs from IDC in a number of histological and clinical features, such as single strand growth, difficulty in detection, and frequent late recurrences. To understand the molecular pathways involved in the clinical characteristics of ILC, we compared the gene expression profiles of luminal A ILC and luminal A IDC using data from TCGA and utilized samples from METABRIC as a validation data set. Top pathways that were significantly enriched in ILC were related to immune response. ILC exhibited a higher activity of almost all types of immune cells based on cell type-specific signatures compared to IDC. Conversely, pathways that were less enriched in ILC were related to protein translation and metabolism, which we functionally validated in cell lines. The higher immune activity uncovered in our study highlights the currently unexplored potential of a response to immunotherapy in a subset of patients with ILC. Furthermore, the lower rates of protein translation and metabolism - known features of tumor dormancy - may play a role in the late recurrences of ILC and lower detection rate in mammography and PET scanning.

Immunotherapy strategies have been emerging as powerful weapons against cancer. Early clinical trials reveal that overall response to immunotherapy is low in breast cancer patients, suggesting that effective strategies to overcome resistance to immunotherapy are urgently needed. In this study, we investigated whether epigenetic reprograming by modulating histone methylation could enhance effector T lymphocyte trafficking and improve therapeutic efficacy of immune checkpoint blockade in breast cancer with focus on triple-negative breast cancer (TNBC) subtype. In silico analysis of The Cancer Genome Atlas (TCGA) data shows that expression of histone lysine-specific demethylase 1 (LSD1) is inversely associated with the levels of cytotoxic T cell-attracting chemokines (C-C motif chemokine ligand 5 (CCL5), C-X-C motif chemokine ligand 9 and 10 (CXCL9, CXCL10)) and programmed death-ligand 1 (PD-L1) in clinical TNBC specimens. Tiling chromatin immunoprecipitation study showed that re-expression of chemokines by LSD1 inhibition is associated with increased H3K4me2 levels at proximal promoter regions. Rescue experiments using concurrent treatment with small interfering RNA or inhibitor of chemokine receptors blocked LSD1 inhibitor-enhanced CD8+ T cell migration, indicating a critical role of key T cell chemokines in LSD1-mediated CD8+ lymphocyte trafficking to the tumor microenvironment. In mice bearing TNBC xenograft tumors, anti-PD-1 antibody alone failed to elicit obvious therapeutic effect. However, combining LSD1 inhibitors with PD-1 antibody significantly suppressed tumor growth and pulmonary metastasis, which was associated with reduced Ki-67 level and augmented CD8+ T cell infiltration in xenograft tumors. Overall, these results suggest that LSD1 inhibition may be an effective adjuvant treatment with immunotherapy as a novel management strategy for poorly immunogenic breast tumors.

MicroRNAs (miRNAs) are post-transcriptional regulators that bind to their target mRNAs through base complementarity. Predicting miRNA targets is a challenging task and various studies showed that existing algorithms suffer from high number of false predictions and low to moderate overlap in their predictions. Until recently, very few algorithms considered the dynamic nature of the interactions, including the effect of less specific interactions, the miRNA expression level, and the effect of combinatorial miRNA binding. Addressing these issues can result in a more accurate miRNA:mRNA modeling with many applications, including efficient miRNA-related SNP evaluation. We present a novel thermodynamic model based on the Fermi-Dirac equation that incorporates miRNA expression in the prediction of target occupancy and we show that it improves the performance of two popular single miRNA target finders. Modeling combinatorial miRNA targeting is a natural extension of this model. Two other algorithms show improved prediction efficiency when combinatorial binding models were considered. ComiR (Combinatorial miRNA targeting), a novel algorithm we developed, incorporates the improved predictions of the four target finders into a single probabilistic score using ensemble learning. Combining target scores of multiple miRNAs using ComiR improves predictions over the naïve method for target combination. ComiR scoring scheme can be used for identification of SNPs affecting miRNA binding. As proof of principle, ComiR identified rs17737058 as disruptive to the miR-488-5p:NCOA1 interaction, which we confirmed in vitro. We also found rs17737058 to be significantly associated with decreased bone mineral density (BMD) in two independent cohorts indicating that the miR-488-5p/NCOA1 regulatory axis is likely critical in maintaining BMD in women. With increasing availability of comprehensive high-throughput datasets from patients ComiR is expected to become an essential tool for miRNA-related studies.

Chromatin immunoprecipitation coupled with massively parallel sequencing (ChIP-seq) is increasingly being applied to study genome-wide binding sites of transcription factors. There is an increasing interest in understanding the mechanism of action of co-regulator proteins, which do not bind DNA directly, but exert their effects by binding to transcription factors such as the estrogen receptor (ER). However, due to the nature of detecting indirect protein-DNA interaction, ChIP-seq signals from co-regulators can be relatively weak and thus biologically meaningful interactions remain difficult to identify.

Silencing mediator of retinoic acid and thyroid hormone receptor (SMRT), also known as nuclear corepressor 2 (NCOR2) is a transcriptional corepressor for multiple members of the nuclear receptor superfamily of transcription factors, including estrogen receptor-α (ERα). In the classical model of corepressor action, SMRT binds to antiestrogen-bound ERα at target promoters and represses ERα transcriptional activity and gene expression. Herein SMRT mRNA and protein expression was examined in a panel of 30 breast cancer cell lines. Expression of both parameters was found to vary considerably amongst lines and the correlation between protein and mRNA expression was very poor (R (2) = 0.0775). Therefore, SMRT protein levels were examined by immunohistochemical staining of a tissue microarray of 866 patients with stage I-II breast cancer. Nuclear and cytoplasmic SMRT were scored separately according to the Allred score. The majority of tumors (67 %) were negative for cytoplasmic SMRT, which when detected was found at very low levels. In contrast, nuclear SMRT was broadly detected. There was no significant difference in time to recurrence (TTR) according to SMRT expression levels in the ERα-positive tamoxifen-treated patients (P = 0.297) but the difference was significant in the untreated patients (P = 0.01). In multivariate analysis, ERα-positive tamoxifen-untreated patients with high nuclear SMRT expression (SMRT 5-8, i.e., 2nd to 4th quartile) had a shorter TTR (HR = 1.94, 95 % CI, 1.24-3.04; P = 0.004) while there was no association with SMRT expression for ERα-positive tamoxifen-treated patients. There was no association between SMRT expression and overall survival for patients, regardless of whether they received tamoxifen. Thus while SMRT protein expression was not predictive of outcome after antiestrogen therapy, it may have value in predicting tumor recurrence in patients not receiving adjuvant tamoxifen therapy.

Invasive lobular breast carcinoma (ILC) accounts for 10% to 15% of breast cancers diagnosed annually. Evidence suggests that some aspects of endocrine treatment response might differ between invasive ductal carcinoma (IDC) and ILC, and that patients with ILC have worse long-term survival. We analyzed The Cancer Genome Atlas dataset and observed lower levels of ESR1 mRNA (P = 0.002) and ERα protein (P = 0.038) in ER+ ILC (n = 137) compared to IDC (n = 554), and further confirmed the mRNA difference in a local UPMC cohort (ILC, n = 143; IDC, n = 877; P < 0.005). In both datasets, the correlation between ESR1 mRNA and ERα protein was weaker in ILC, suggesting differential post-transcriptional regulation of ERα. In vitro, 17β-estradiol (E2) decreased the rate of degradation and increased the half-life of ERα in ILC cell lines, whereas the opposite was observed in IDC cell lines. Further, E2 failed to induce robust ubiquitination of ERα in ILC cells. To determine the potential clinical relevance of these findings, we evaluated the effect of 2 selective estrogen receptor downregulators (SERDs), ICI 182,780 and AZD9496, on ERα turnover and cell growth. While ICI 182,780 and AZD9496 showed similar effects in IDC cells, in ILC cell lines, AZD9496 was not as effective as ICI 182,780 in decreasing ERα stability and E2-induced proliferation. Furthermore, AZD9496 exhibited partial agonist activity in growth assays in ILC cell lines. Our study provides evidence for a distinct ERα regulation by SERDs in ILC cell lines, and therefore it is important to include ILC models into preclinical and clinical testing of novel SERDs.

The molecular events and transcriptional plasticity driving brain metastasis in clinically relevant breast tumor subtypes has not been determined. Here we comprehensively dissect genomic, transcriptomic and clinical data in patient-matched longitudinal tumor samples, and unravel distinct transcriptional programs enriched in brain metastasis. We report on subtype specific hub genes and functional processes, central to disease-affected networks in brain metastasis. Importantly, in luminal brain metastases we identify homologous recombination deficiency operative in transcriptomic and genomic data with recurrent breast mutational signatures A, F and K, associated with mismatch repair defects, TP53 mutations and homologous recombination deficiency (HRD) respectively. Utilizing PARP inhibition in patient-derived brain metastatic tumor explants we functionally validate HRD as a key vulnerability. Here, we demonstrate a functionally relevant HRD evident at genomic and transcriptomic levels pointing to genomic instability in breast cancer brain metastasis which is of potential translational significance.

Tumor-associated macrophages (TAMs) exert profound influence over breast cancer progression, promoting immunosuppression, angiogenesis, and metastasis. Neuropilin-2 (NRP2), consisting of the NRP2a and NRP2b isoforms, is a co-receptor for heparin-binding growth factors including VEGF-C and Class 3 Semaphorins. Selective upregulation in response to environmental stimuli and independent signaling pathways endow the NRP2 isoforms with unique functionality, with NRP2b promoting increased Akt signaling via receptor tyrosine kinases including VEGFRs, MET, and PDGFR. Although NRP2 has been shown to regulate macrophage/TAM biology, the role of the individual NRP2a/NRP2b isoforms in TAMs has yet to be evaluated. Using transcriptional profiling and spectral flow cytometry, we show that NRP2 isoform expression was significantly higher in TAMs from murine mammary tumors. NRP2a/NRP2b levels in human breast cancer metastasis were dependent upon the anatomic location of the tumor and significantly correlated with TAM infiltration in both primary and metastatic breast cancers. We define distinct phenotypes of NRP2 isoform-expressing TAMs in mouse models of breast cancer and within malignant pleural effusions from breast cancer patients which were exclusive of neuropilin-1 expression. Genetic depletion of either NRP2 isoform in macrophages resulted in a dramatic reduction of LPS-induced IL-10 production, defects in phagosomal processing of apoptotic breast cancer cells, and increase in cancer cell migration following co-culture. By contrast, depletion of NRP2b, but not NRP2a, inhibited production of IL-6. These results suggest that NRP2 isoforms regulate both shared and unique functionality in macrophages and are associated with distinct TAM subsets in breast cancer.

Model organisms are instrumental substitutes for human studies to expedite basic, translational, and clinical research. Despite their indispensable role in mechanistic investigation and drug development, molecular congruence of animal models to humans has long been questioned and debated. Little effort has been made for an objective quantification and mechanistic exploration of a model organism's resemblance to humans in terms of molecular response under disease or drug treatment. We hereby propose a framework, namely Congruence Analysis for Model Organisms (CAMO), for transcriptomic response analysis by developing threshold-free differential expression analysis, quantitative concordance/discordance scores incorporating data variabilities, pathway-centric downstream investigation, knowledge retrieval by text mining, and topological gene module detection for hypothesis generation. Instead of a genome-wide vague and dichotomous answer of "poorly" or "greatly" mimicking humans, CAMO assists researchers to numerically quantify congruence, to dissect true cross-species differences from unwanted biological or cohort variabilities, and to visually identify molecular mechanisms and pathway subnetworks that are best or least mimicked by model organisms, which altogether provides foundations for hypothesis generation and subsequent translational decisions.

As one of the most successful cancer therapeutic targets, estrogen receptor-α (ER/ESR1) has been extensively studied over the past few decades. Sequencing technological advances have enabled genome-wide analysis of ER action. However, comparison of individual studies is limited by different experimental designs, and few meta-analyses are available. Here, we established the EstroGene database through unified processing of data from 246 experiments including 136 transcriptomic, cistromic, and epigenetic datasets focusing on estradiol (E2)-triggered ER activation across 19 breast cancer cell lines. A user-friendly browser (https://estrogene.org/) was generated for multiomic data visualization involving gene inquiry under user-defined experimental conditions and statistical thresholds. Notably, annotation of metadata associated with public datasets revealed a considerable lack of experimental details. Comparison of independent RNA-seq or ER ChIP-seq data with the same design showed large variability and only strong effects could be consistently detected. Temporal estrogen response metasignatures were defined, and the association of E2 response rate with temporal transcriptional factors, chromatin accessibility, and heterogeneity of ER expression was evaluated. Unexpectedly, harmonizing 146 E2-induced transcriptomic datasets uncovered a subset of genes harboring bidirectional E2 regulation, which was linked to unique transcriptional factors and highly associated with immune surveillance in the clinical setting. Furthermore, the context dependent E2 response programs were characterized in MCF7 and T47D cell lines, the two most frequently used models in the EstroGene database. Collectively, the EstroGene database provides an informative and practical resource to the cancer research community to uniformly evaluate key reproducible features of ER regulomes and unravels modes of ER signaling.