Search

Background: Increased concentrations of choline-containing compounds are frequently observed in breast carcinomas, and may serve as biomarkers for both diagnostic and treatment monitoring purposes. However, underlying mechanisms for the abnormal choline metabolism are poorly understood.Methods: The concentrations of choline-derived metabolites were determined in xenografted primary human breast carcinomas, representing basal-like and luminal-like subtypes. Quantification of metabolites in fresh frozen tissue was performed using high-resolution magic angle spinning magnetic resonance spectroscopy (HR MAS MRS). The expression of genes involved in phosphatidylcholine (PtdCho) metabolism was retrieved from whole genome expression microarray analyses. The metabolite profiles from xenografts were compared with profiles from human breast cancer, sampled from patients with estrogen/progesterone receptor positive (ER+/PgR+) or triple negative (ER-/PgR-/HER2-) breast cancer.Results: In basal-like xenografts, glycerophosphocholine (GPC) concentrations were higher than phosphocholine (PCho) concentrations, whereas this pattern was reversed in luminal-like xenografts. These differences may be explained by lower choline kinase (CHKA, CHKB) expression as well as higher PtdCho degradation mediated by higher expression of phospholipase A2 group 4A (PLA2G4A) and phospholipase B1 (PLB1) in the basal-like model. The glycine concentration was higher in the basal-like model. Although glycine could be derived from energy metabolism pathways, the gene expression data suggested a metabolic shift from PtdCho synthesis to glycine formation in basal-like xenografts. In agreement with results from the xenograft models, tissue samples from triple negative breast carcinomas had higher GPC/PCho ratio than samples from ER+/PgR+ carcinomas, suggesting that the choline metabolism in the experimental models is representative for luminal-like and basal-like human breast cancer.Conclusions: The differences in choline metabolite concentrations corresponded well with differences in gene expression, demonstrating distinct metabolic profiles in the xenograft models representing basal-like and luminal-like breast cancer. The same characteristics of choline metabolite profiles were also observed in patient material from ER+/PgR+ and triple-negative breast cancer, suggesting that the xenografts are relevant model systems for studies of choline metabolism in luminal-like and basal-like breast cancer. [ABSTRACT FROM AUTHOR]

Treatment with the anti-angiogenic drug bevacizumab in addition to chemotherapy has shown efficacy for breast cancer in some clinical trials, but better biomarkers are needed to optimally select patients for treatment. Here, we present an omics approach where DNA methylation profiles are integrated with gene expression and results from proteomic data in breast cancer patients to predict response to therapy and pinpoint response-related epigenetic events. Fresh-frozen tumor biopsies taken before, during, and after treatment from human epidermal growth factor receptor 2 negative non-metastatic patients receiving neoadjuvant chemotherapy with or without bevacizumab were subjected to molecular profiling. Here, we report that DNA methylation at enhancer CpGs related to cell cycle regulation can predict response to chemotherapy and bevacizumab for the estrogen receptor positive subset of patients (AUC = 0.874), and we validated this observation in an independent patient cohort with a similar treatment regimen (AUC = 0.762). Combining the DNA methylation scores with the scores from a previously published protein signature resulted in a slight increase in the prediction performance (AUC = 0.784). We also show that tumors receiving the combination treatment underwent more extensive epigenetic alterations. Finally, we performed an integrative expression-methylation quantitative trait loci analysis on alterations in DNA methylation and gene expression levels, showing that the epigenetic alterations that occur during treatment are different between responders and non-responders and that these differences may be explained by the proliferation-epithelial-to-mesenchymal transition axis through the activity of grainyhead like transcription factor 2. Using tumor purity computed from copy number data, we developed a method for estimating cancer cell-specific methylation to confirm that the association to response reflects DNA methylation in cancer cells. Taken together, these results support the potential for clinical benefit of the addition of bevacizumab to chemotherapy when administered to the correct patients., (© 2024 The Authors. Molecular Oncology published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Disease recurrence and drug resistance are major challenges in the clinical management of patients with colorectal cancer liver metastases (CLM), and because tumors are generally microsatellite stable (MSS), responses to immune therapies are poor. The mesenchymal phenotype is overrepresented in treatment-resistant cancers and is associated with an immunosuppressed microenvironment. The aim of this work was to molecularly identify and characterize a mesenchymal subgroup of MSS CLM to identify novel therapeutic approaches. We here generated a mesenchymal gene expression signature by analysis of resection specimens from 38 patients with CLM using ranked expression level of the epithelial-to-mesenchymal transition-related transcription factor PRRX1 . Downstream pathway analysis based on the resulting gene signature was performed and independent, publicly available datasets were used to validate the findings. A subgroup comprising 16% of the analyzed CLM samples were classified as mesenchymal, or belonging to the PRRX1 high group. Analysis of the PRRX1 signature genes revealed a distinct immunosuppressive phenotype with high expression of immune checkpoints HAVCR2/TIM-3 and VISTA, in addition to the M2 macrophage marker CD163. The findings were convincingly validated in datasets from three external CLM cohorts. Upregulation of immune checkpoints HAVCR2/TIM-3 and VISTA in the PRRX1 high subgroup is a novel finding, and suggests immune evasion beyond the PD-1/PD-L1 axis, which may contribute to poor response to PD-1/PD-L1-directed immune therapy in MSS colorectal cancer. Importantly, these checkpoints represent potential novel opportunities for immune-based therapy approaches in a subset of MSS CLM., Significance: CLM is an important cause of colorectal cancer mortality where the majority of patients have yet to benefit from immunotherapies. In this study of gene expression profiling analyses, we uncovered novel immune checkpoint targets in a subgroup of patients with MSS CLMs harboring a mesenchymal phenotype., Competing Interests: V. Nygaard reports grants from Research Council of Norway, South-Eastern Norway Health Authority, and Sister Institution Network Fund during the conduct of the study. Å.A. Fretland reports personal fees from Bayer and Olympus outside the submitted work. M.H. Haugen reports grants from South-Eastern Healt Authorities Norway and MD Anderson Sister Institution Network Fund (SINF) during the conduct of the study. No disclosures were reported by the other authors., (© 2023 The Authors; Published by the American Association for Cancer Research.)

Most cancer alterations occur in the noncoding portion of the human genome, where regulatory regions control gene expression. The discovery of noncoding mutations altering the cells' regulatory programs has been limited to few examples with high recurrence or high functional impact. Here, we show that transcription factor binding sites (TFBSs) have similar mutation loads to those in protein-coding exons. By combining cancer somatic mutations in TFBSs and expression data for protein-coding and miRNA genes, we evaluate the combined effects of transcriptional and post-transcriptional alterations on the regulatory programs in cancers. The analysis of seven TCGA cohorts culminates with the identification of protein-coding and miRNA genes linked to mutations at TFBSs that are associated with a cascading trans-effect deregulation on the cells' regulatory programs. Our analyses of cis-regulatory mutations associated with miRNAs recurrently predict 12 mature miRNAs (derived from 7 precursors) associated with the deregulation of their target gene networks. The predictions are enriched for cancer-associated protein-coding and miRNA genes and highlight cis-regulatory mutations associated with the dysregulation of key pathways associated with carcinogenesis. By combining transcriptional and post-transcriptional regulation of gene expression, our method predicts cis-regulatory mutations related to the dysregulation of key gene regulatory networks in cancer patients., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Background: Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear., Methods: Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes., Results: Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions., Conclusion: This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction., (© 2021. The Author(s).)

Basal-like breast cancer (BLBC) is the most aggressive subtype of breast tumors with poor prognosis and limited molecular-targeted therapy options. We show that BLBC cells have a high Cys demand and reprogrammed Cys metabolism. Patient-derived BLBC tumors from four different cohorts exhibited elevated expression of the transsulfuration enzyme cystathione β-synthetase (CBS). CBS silencing (shCBS) made BLBC cells less invasive, proliferate slower, more vulnerable to oxidative stress and cystine (CySSCy) deprivation, prone to ferroptosis, and less responsive to HIF1-α activation under hypoxia. shCBS xenograft tumors grew slower than controls and exhibited impaired angiogenesis and larger necrotic areas. Sulfur metabolite profiling suggested that realigned sulfide/persulfide-inducing functions of CBS are important in BLBC tumor progression. Supporting this, the exclusion of serine, a substrate of CBS for producing Cys but not for producing sulfide/persulfide, did not exacerbate CySSCy deprivation-induced ferroptosis in shCBS BLBC cells. Impaired Tyr phosphorylation was detected in shCBS cells and xenografts, likely due to persulfidation-inhibited phosphatase functions. Overexpression of cystathione γ-lyase (CSE), which can also contribute to cellular sulfide/persulfide production, compensated for the loss of CBS activities, and treatment of shCBS xenografts with a CSE inhibitor further blocked tumor growth. Glutathione and protein-Cys levels were not diminished in shCBS cells or xenografts, but levels of Cys persulfidation and the persulfide-catabolizing enzyme ETHE1 were suppressed. Finally, expression of enzymes of the oxidizing Cys catabolism pathway was diminished, but expression of the persulfide-producing CARS2 was elevated in human BLBC tumors. Hence, the persulfide-producing pathways are major targetable determinants of BLBC pathology that could be therapeutically exploited., Competing Interests: The authors declare no competing interest.

Polymorphisms in the PER3 gene have been associated with several human disease phenotypes, including sleep disorders and cancer. In particular, the long allele of a variable number of tandem repeat (VNTR) polymorphism has been previously linked to an increased risk of breast cancer. Here we carried out a combined germline and somatic genetic analysis of the role of the PER3 VNRT polymorphism in breast cancer. The combined data from 8284 individuals showed a non-significant trend towards increased breast cancer risk in the 5-repeat allele homozygous carriers (OR = 1.17, 95% CI: 0.97-1.42). We observed allelic imbalance at the PER3 locus in matched blood and tumor DNA samples, showing a significant retention of the long variant (risk) allele in tumor samples, and a preferential loss of the short repetition allele (p = 0.0005). Gene co-expression analysis in healthy and tumoral breast tissue samples uncovered significant associations between PER3 expression levels with those from genes which belong to several cancer-associated pathways. Finally, relapse-free survival (RFS) analysis showed that low expression levels of PER3 were linked to a significant lower RSF in luminal A (p = 3 × 10 -12 ) but not in the rest of breast cancer subtypes., (© 2021. The Author(s).)

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.

Despite significant advancements in breast cancer (BC) research, clinicians lack robust serological protein markers for accurate diagnostics and tumor stratification. Tumor interstitial fluid (TIF) accumulates aberrantly externalized proteins within the local tumor space, which can potentially gain access to the circulatory system. As such, TIF may represent a valuable starting point for identifying relevant tumor-specific serological biomarkers. The aim of the study was to perform comprehensive proteomic profiling of TIF to identify proteins associated with BC tumor status and subtype. A liquid chromatography tandem mass spectrometry (LC-MS/MS) analysis of 35 TIFs of three main subtypes: luminal (19), Her2 (4), and triple-negative (TNBC) (12) resulted in the identification of > 8800 proteins. Unsupervised hierarchical clustering segregated the TIF proteome into two major clusters, luminal and TNBC/Her2 subgroups. High-grade tumors enriched with tumor infiltrating lymphocytes (TILs) were also stratified from low-grade tumors. A consensus analysis approach, including differential abundance analysis, selection operator regression, and random forest returned a minimal set of 24 proteins associated with BC subtypes, receptor status, and TIL scoring. Among them, a panel of 10 proteins, AGR3, BCAM, CELSR1, MIEN1, NAT1, PIP4K2B, SEC23B, THTPA, TMEM51, and ULBP2, was found to stratify the tumor subtype-specific TIFs. In particular, upregulation of BCAM and CELSR1 differentiates luminal subtypes, while upregulation of MIEN1 differentiates Her2 subtypes. Immunohistochemistry analysis showed a direct correlation between protein abundance in TIFs and intratumor expression levels for all 10 proteins. Sensitivity and specificity were estimated for this protein panel by using an independent, comprehensive breast tumor proteome dataset. The results of this analysis strongly support our data, with eight of the proteins potentially representing biomarkers for stratification of BC subtypes. Five of the most representative proteomics databases currently available were also used to estimate the potential for these selected proteins to serve as putative serological markers., (© 2020 The Authors. Molecular Oncology published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Purpose: Antiangiogenic therapy using bevacizumab has proven effective for a number of cancers; however, in breast cancer (BC), there is an unmet need to identify patients who benefit from such treatment., Patients and Methods: In the NeoAva phase II clinical trial, patients (N = 132) with large (≥ 25 mm) human epidermal growth factor receptor 2 (HER2)-negative primary tumors were randomly assigned 1:1 to treatment with neoadjuvant chemotherapy (CTx) alone or in combination with bevacizumab (Bev plus CTx). The ratio of the tumor size after relative to before treatment was calculated into a continuous response scale. Tumor biopsies taken prior to neoadjuvant treatment were analyzed by reverse-phase protein arrays (RPPA) for expression levels of 210 BC-relevant (phospho-) proteins. Lasso regression was used to derive a predictor of tumor shrinkage from the expression of selected proteins prior to treatment., Results: We identified a nine-protein signature score named vascular endothelial growth factor inhibition response predictor (ViRP) for use in the Bev plus CTx treatment arm able to predict with accuracy pathologic complete response (pCR) (area under the curve [AUC] = 0.85; 95% CI, 0.74 to 0.97) and low residual cancer burden (RCB 0/I) (AUC = 0.80; 95% CI, 0.68 to 0.93). The ViRP score was significantly lower in patients with pCR ( P < .001) and in patients with low RCB ( P < .001). The ViRP score was internally validated on mRNA data and the resultant surrogate mRNA ViRP score significantly separated the pCR patients ( P = .016). Similarly, the mRNA ViRP score was validated ( P < .001) in an independent phase II clinical trial (PROMIX)., Conclusion: Our ViRP score, integrating the expression of nine proteins and validated on mRNA data both internally and in an independent clinical trial, may be used to increase the likelihood of benefit from treatment with bevacizumab combined with chemotherapy in patients with HER2-negative BC., Competing Interests: The following represents disclosure information provided by authors of this manuscript. All relationships are considered compensated unless otherwise noted. Relationships are self-held unless noted. I = Immediate Family Member, Inst = My Institution. Relationships may not relate to the subject matter of this manuscript. For more information about ASCO’s conflict of interest policy, please refer to www.asco.org/rwc or ascopubs.org/cci/author-center. Open Payments is a public database containing information reported by companies about payments made to US-licensed physicians (Open Payments). Mads H. HaugenPatents, Royalties, Other Intellectual Property: Patent application 62/969770 to USPTOOle Christian LingjærdeConsulting or Advisory Role: NovartisThomas HatschekConsulting or Advisory Role: Roche, Pfizer, Pierre Fabre Research Funding: Roche, Pfizer Travel, Accommodations, Expenses: RocheAnne-Lise Børresen-DaleEmployment: Arctic Pharma AS, PubGene Stock and Other Ownership Interests: Arctic Pharma ASGordon B. MillsStock and Other Ownership Interests: Catena, SignalChem, Tarveda Therapeutics, ImmunoMET Honoraria: Nuevolution: AstraZeneca, Tarveda Therapeutics, Tesaro, Symphogen, PDX Pharmacy, ImmunoMET, Lilly Consulting or Advisory Role: AstraZeneca, SignalChem, Tarveda Therapeutics, Symphogen, Takeda/Millennium, PDX Pharmacy, ImmunoMET, Lilly, Turbine, ION Pharma, Zentalis Research Funding: Adelson Medical Research Foundation, AstraZeneca, NanoString Technologies, Breast Cancer Research Foundation, Karus Therapeutics, Pfizer, Prospect Creek Foundation, Tarveda Therapeutics, Ions Pharmaceuticals, ImmunoMET Patents, Royalties, Other Intellectual Property: HRD assay to Myriad Genetics, DSP technology patent with Nanostring Travel, Accommodations, Expenses: AstraZeneca, Pfizer, Symphogen, Chrysalis Biomedical Advisors, ImmunoMET, Michigan Primary Care ConsortiumGunhild M. MælandsmoPatents, Royalties, Other Intellectual Property: Patent application submitted for a nine-protein/gene panel predicting response to anti VEGF therapies in combination with chemotherapyOlav EngebraatenPatents, Royalties, Other Intellectual Property: Patent application pending for a biomarker for antibody drug conjugates, Patent application submitted for a nine-protein/gene panel predicting response to anti VEGF therapies in combination with chemotherapy No other potential conflicts of interest were reported., (© 2021 by American Society of Clinical Oncology.)

Antiangiogenic drugs are potentially a useful supplement to neoadjuvant chemotherapy for a subgroup of patients with human epidermal growth factor receptor 2 (HER2) negative breast cancer, but reliable biomarkers for improved response are lacking. Here, we report on a randomized phase II clinical trial to study the added effect of bevacizumab in neoadjuvant chemotherapy with FEC100 (5-fluorouracil, epirubicin and cyclophosphamide) and taxanes (n = 132 patients). Gene expression from the tumors was obtained before neoadjuvant treatment, and treatment response was evaluated by residual cancer burden (RCB) at time of surgery. Bevacizumab increased the proportion of complete responders (RCB class 0) from 5% to 20% among patients with estrogen receptor (ER) positive tumors (P = .02). Treatment with bevacizumab was associated with improved 8-year disease-free survival (P = .03) among the good responders (RCB class 0 or I). Patients treated with paclitaxel (n = 45) responded better than those treated with docetaxel (n = 21; P = .03). Improved treatment response was associated with higher proliferation rate and an immune phenotype characterized by high presence of classically activated M1 macrophages, activated NK cells and memory activated CD4 T cells. Treatment with bevacizumab increased the number of adverse events, including hemorrhage, hypertension, infection and febrile neutropenia, but despite this, the ECOG status was not affected., (© 2020 The Authors. International Journal of Cancer published by John Wiley & Sons Ltd on behalf of UICC.)

Background: In precision cancer medicine, the challenge is to prioritize DNA driver events, account for resistance markers, and procure sufficient information for treatment that maintains patient safety. The MetAction project, exploring how tumor molecular vulnerabilities predict therapy response, first established the required workflow for DNA sequencing and data interpretation (2014-2015). Here, we employed it to identify molecularly matched therapy and recorded outcome in end-stage cancer (2016-2019). Material and methods: Metastatic tissue from 26 patients (16 colorectal cancer cases) was sequenced by the Oncomine assay. The study tumor boards interpreted called variants with respect to sensitivity or resistance to matched therapy and recommended single-agent or combination treatment if considered tolerable. The primary endpoint was the rate of progression-free survival 1.3-fold longer than for the most recent systemic therapy. The objective response rate and overall survival were secondary endpoints. Results: Both common and rare actionable alterations were identified. Thirteen patients were found eligible for therapy following review of tumor sensitivity and resistance variants and patient tolerability. The interventions were inhibitors of ALK/ROS1-, BRAF-, EGFR-, FGFR-, mTOR-, PARP-, or PD-1-mediated signaling for 2-3 cases each. Among 10 patients who received treatment until radiologic evaluation, 6 (46% of the eligible cases) met the primary endpoint. Four colorectal cancer patients (15% of the total study cohort) had objective response. The only serious adverse event was a transient colitis, which appeared in 1 of the 2 patients given PD-1 inhibitor with complete response. Apart from those two, overall survival was similar for patients who did and did not receive study treatment. Conclusions: The systematic MetAction approach may point forward to a refined framework for how to interpret the complexity of sensitivity versus resistance and patient safety that resides in tumor sequence data, for the possibly improved outcome of precision cancer medicine in future studies. ClinicalTrials.gov, identifier: NCT02142036.

Somatic copy number alterations are a frequent sign of genome instability in cancer. A precise characterization of the genome architecture would reveal underlying instability mechanisms and provide an instrument for outcome prediction and treatment guidance. Here we show that the local spatial behavior of copy number profiles conveys important information about this architecture. Six filters were defined to characterize regional traits in copy number profiles, and the resulting Copy Aberration Regional Mapping Analysis (CARMA) algorithm was applied to tumors in four breast cancer cohorts (n = 2919). The derived motifs represent a layer of information that complements established molecular classifications of breast cancer. A score reflecting presence or absence of motifs provided a highly significant independent prognostic predictor. Results were consistent between cohorts. The nonsite-specific occurrence of the detected patterns suggests that CARMA captures underlying replication and repair defects and could have a future potential in treatment stratification.

Identifying the underlying genetic drivers of the heritability of breast cancer prognosis remains elusive. We adapt a network-based approach to handle underpowered complex datasets to provide new insights into the potential function of germline variants in breast cancer prognosis. This network-based analysis studies ~7.3 million variants in 84,457 breast cancer patients in relation to breast cancer survival and confirms the results on 12,381 independent patients. Aggregating the prognostic effects of genetic variants across multiple genes, we identify four gene modules associated with survival in estrogen receptor (ER)-negative and one in ER-positive disease. The modules show biological enrichment for cancer-related processes such as G-alpha signaling, circadian clock, angiogenesis, and Rho-GTPases in apoptosis.

Gene expression profiling of tumours is an important source of information for cancer patient stratification. Detecting subtle alterations of gene expression remains a challenge, however. Here, we propose a novel tool for high-sensitivity detection of differential pathway activity in tumours. For a pathway defined by a collection of genes, the samples are projected onto a low-dimensional manifold in the subspace spanned by those genes. For each sample, a score is next found by calculating the distance between each projected sample and the projection of a subgroup of reference samples. Depending on the aim of the analysis and the available data, the reference samples may represent e.g. normal tissue or tumour samples with a particular genotype or phenotype. The proposed tool, PathTracer, is demonstrated on gene expression data from 1952 invasive breast cancer samples, 10 DCIS, 9 benign samples and 144 tumour adjacent normal breast tissue samples. PathTracer scores are shown to predict survival, clinical subtypes, cellular proliferation and genomic instability. Furthermore, predictions are shown to outperform those obtained with other comparable methods.

Background: Anti-tumorigenic vs pro-tumorigenic roles of estrogen receptor-beta (ESR2) in breast cancer remain unsettled. We investigated the potential of TP53 status to be a determinant of the bi-faceted role of ESR2 and associated therapeutic implications for triple negative breast cancer (TNBC)., Methods: ESR2-TP53 interaction was analyzed with multiple assays including the in situ proximity ligation assay. Transcriptional effects on TP53-target genes and cell proliferation in response to knocking down or overexpressing ESR2 were determined. Patient survival according to ESR2 expression levels and TP53 mutation status was analyzed in the basal-like TNBC subgroup in the Molecular Taxonomy of Breast Cancer International Consortium (n = 308) and Roswell Park Comprehensive Cancer Center (n = 46) patient cohorts by univariate Cox regression and log-rank test. All statistical tests are two-sided., Results: ESR2 interaction with wild-type and mutant TP53 caused pro-proliferative and anti-proliferative effects, respectively. Depleting ESR2 in cells expressing wild-type TP53 resulted in increased expression of TP53-target genes CDKN1A (control group mean [SD] = 1 [0.13] vs ESR2 depletion group mean [SD] = 2.08 [0.24], P = .003) and BBC3 (control group mean [SD] = 1 [0.06] vs ESR2 depleted group mean [SD] = 1.92 [0.25], P = .003); however, expression of CDKN1A (control group mean [SD] = 1 [0.21] vs ESR2 depleted group mean [SD] = 0.56 [0.12], P = .02) and BBC3 (control group mean [SD] = 1 [0.03] vs ESR2 depleted group mean [SD] = 0.55 [0.09], P = .008) was decreased in cells expressing mutant TP53. Overexpressing ESR2 had opposite effects. Tamoxifen increased ESR2-mutant TP53 interaction, leading to reactivation of TP73 and apoptosis. High levels of ESR2 expression in mutant TP53-expressing basal-like tumors is associated with better prognosis (Molecular Taxonomy of Breast Cancer International Consortium cohort: log-rank P = .001; hazard ratio = 0.26, 95% confidence interval = 0.08 to 0.84, univariate Cox P = .02)., Conclusions: TP53 status is a determinant of the functional duality of ESR2. Our study suggests that ESR2-mutant TP53 combination prognosticates survival in TNBC revealing a novel strategy to stratify TNBC for therapeutic intervention potentially by repurposing tamoxifen., (© The Author(s) 2019. Published by Oxford University Press.)

One of the hallmarks of cancer is sustained angiogenesis. Favorable results have been reported in some breast cancer (BC) patients receiving antiangiogenic therapy with bevacizumab (Bev) in combination with chemotherapy, and further knowledge on how Bev can be optimally combined with conventional treatment to increase efficacy is strongly needed. In this randomized, neoadjuvant phase II clinical trial, 132 patients with HER2-negative, nonmetastatic BC were treated with Bev in combination with sequential chemotherapy. Biopsies were sampled before treatment, after 12 weeks with anthracycline and after taxane therapy at week 25. MicroRNA (miRNA) expression profiling was performed on biopsies from each time point. Altogether, 241 biopsies were analyzed with the aim of identifying miRNA-based biomarkers of response to therapy. Results from the miRNA analyses were reported for the ER-positive cohort, which were previously demonstrated to benefit from antiangiogenic therapy in this study. For both treatment arms of this cohort, significantly different expression was observed for 217 miRNAs between objective responding and nonresponding patients before treatment initiation. These miRNAs have been linked to regulation of epithelial-mesenchymal transition, metastasis, and tumor growth, among other processes. Bev in combination with chemotherapy resulted in similar miRNA changes to chemotherapy alone. However, the deregulation of miRNA expression occurred earlier in the Bev arm. In both arms, tumor suppressor miRNAs were found upregulated after treatment, while oncogenic miRNAs were downregulated in the Bev arm. Patients responding to Bev showed a strong correlation between deregulated miRNAs and decreased proliferation score during the course of treatment, with downregulation of miR-4465 as the strongest indicator of reduced proliferation. Integrative analyses at miRNA-, gene-, and protein expression further indicated a longitudinal decrease in proliferation. Altogether, the results indicate that proliferation might represent a predictive factor for increased Bev sensitivity, which may aid in the identification of patients who could potentially benefit from Bev., (© 2019 The Authors. Published by FEBS Press and John Wiley & Sons Ltd.)

The usefulness of mechanistic models to disentangle complex multiscale cancer processes, such as treatment response, has been widely acknowledged. However, a major barrier for multiscale models to predict treatment outcomes in individual patients lies in their initialization and parametrization, which needs to reflect individual cancer characteristics accurately. In this study, we use multitype measurements acquired routinely on a single breast tumor, including histopathology, MRI, and molecular profiling, to personalize parts of a complex multiscale model of breast cancer treated with chemotherapeutic and antiangiogenic agents. The model accounts for drug pharmacokinetics and pharmacodynamics. We developed an open-source computer program that simulates cross-sections of tumors under 12-week therapy regimens and used it to individually reproduce and elucidate treatment outcomes of 4 patients. Two of the tumors did not respond to therapy, and model simulations were used to suggest alternative regimens with improved outcomes dependent on the tumor's individual characteristics. It was determined that more frequent and lower doses of chemotherapy reduce tumor burden in a low proliferative tumor while lower doses of antiangiogenic agents improve drug penetration in a poorly perfused tumor. Furthermore, using this model, we were able to correctly predict the outcome in another patient after 12 weeks of treatment. In summary, our model bridges multitype clinical data to shed light on individual treatment outcomes. SIGNIFICANCE: Mathematical modeling is used to validate possible mechanisms of tumor growth, resistance, and treatment outcome., (©2019 American Association for Cancer Research.)

Most patients whose large bowel cancer has spread to other organs do not respond to immune therapy. We detected a rare gene mutation, termed 9p24.1 copy-number gain (CNG), in an otherwise incurable colorectal cancer that provoked an immune therapy response. We identified this gene mutation by gene-panel sequencing of DNA from a liver metastasis biopsy from a patient who had disease refractory to standard therapies. Following immune checkpoint blockade (ICB) with pembrolizumab (anti-PD-1), the patient experienced conversion of the tumor phenotype from one with epithelial features to that of an inflamed microenvironment, detected by high-resolution RNA sequencing. Circulating tumor DNA disappeared over the first weeks of therapy. As assessed by standard radiographic measurement, the patient had a partial response that was durable. This patient's response may support the use of histology-agnostic ICB in solid tumors that carry the rare 9p24.1 CNG., (©2019 American Association for Cancer Research.)

Global loss of DNA methylation and CpG island (CGI) hypermethylation are key epigenomic aberrations in cancer. Global loss manifests itself in partially methylated domains (PMDs) which extend up to megabases. However, the distribution of PMDs within and between tumor types, and their effects on key functional genomic elements including CGIs are poorly defined. We comprehensively show that loss of methylation in PMDs occurs in a large fraction of the genome and represents the prime source of DNA methylation variation. PMDs are hypervariable in methylation level, size and distribution, and display elevated mutation rates. They impose intermediate DNA methylation levels incognizant of functional genomic elements including CGIs, underpinning a CGI methylator phenotype (CIMP). Repression effects on tumor suppressor genes are negligible as they are generally excluded from PMDs. The genomic distribution of PMDs reports tissue-of-origin and may represent tissue-specific silent regions which tolerate instability at the epigenetic, transcriptomic and genetic level.

In the preceding decades, molecular characterization has revolutionized breast cancer (BC) research and therapeutic approaches. Presented herein, an unbiased analysis of breast tumor proteomes, inclusive of 9995 proteins quantified across all tumors, for the first time recapitulates BC subtypes. Additionally, poor-prognosis basal-like and luminal B tumors are further subdivided by immune component infiltration, suggesting the current classification is incomplete. Proteome-based networks distinguish functional protein modules for breast tumor groups, with co-expression of EGFR and MET marking ductal carcinoma in situ regions of normal-like tumors and lending to a more accurate classification of this poorly defined subtype. Genes included within prognostic mRNA panels have significantly higher than average mRNA-protein correlations, and gene copy number alterations are dampened at the protein-level; underscoring the value of proteome quantification for prognostication and phenotypic classification. Furthermore, protein products mapping to non-coding genomic regions are identified; highlighting a potential new class of tumor-specific immunotherapeutic targets.

Background: We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry., Methods: Meta-analyses included summary estimates based on Cox models of twelve datasets using ~10.4 million variants for 96,661 women with breast cancer and 7697 events (breast cancer-specific deaths). Oestrogen receptor (ER)-specific analyses were based on 64,171 ER-positive (4116) and 16,172 ER-negative (2125) patients. We evaluated the probability of a signal to be a true positive using the Bayesian false discovery probability (BFDP)., Results: We did not find any variant associated with breast cancer-specific mortality at P < 5 × 10 -8 . For ER-positive disease, the most significantly associated variant was chr7:rs4717568 (BFDP = 7%, P = 1.28 × 10 -7 , hazard ratio [HR] = 0.88, 95% confidence interval [CI] = 0.84-0.92); the closest gene is AUTS2. For ER-negative disease, the most significant variant was chr7:rs67918676 (BFDP = 11%, P = 1.38 × 10 -7 , HR = 1.27, 95% CI = 1.16-1.39); located within a long intergenic non-coding RNA gene (AC004009.3), close to the HOXA gene cluster., Conclusions: We uncovered germline variants on chromosome 7 at BFDP < 15% close to genes for which there is biological evidence related to breast cancer outcome. However, the paucity of variants associated with mortality at genome-wide significance underpins the challenge in providing genetic-based individualised prognostic information for breast cancer patients.

Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative ( R < 0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2 , CREBBP , and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer., (© 2019 Smid et al.; Published by Cold Spring Harbor Laboratory Press.)

In the Methods section of this Article, 'greater than' should have been 'less than' in the sentence 'Putative regions of clustered rearrangements were identified as having an average inter-rearrangement distance that was at least 10 times greater than the whole-genome average for the individual sample. '. The Article has not been corrected.

The original version of this Article omitted a declaration from the competing interests statement, which should have included the following: 'K.P.W. is President of Tempus Lab, Inc., Chicago, IL, USA'. This has now been corrected in both the PDF and HTML versions of the Article.

Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: HER2 positive Breast Cancers (BC) have aggressive behavior and poor prognosis. Previously, we have identified miR-342-5p as an upstream regulator of HER2 signaling, as well as inhibitor of HER2 positive BC cell line growth., Objective: Here, we aimed to further investigate the molecular mechanisms behind miR-342-5pinduced HER2 pathway deregulation., Method: Two HER2 amplified breast cancer cell lines were transiently transfected with miR-342-5p mimic or negative control, and gene expression was analyzed by Agilent microarrays. Three clinical datasets with BC patients were used to identify correlations between candidate genes and miR-342- 5p, and associations with survival., Results: Pathway analyses of all deregulated genes revealed a significant suppression of the HER2 downstream pathways ERK/MAPK and SAPK/JNK, whereas the miR-342-5p predicted target genes were enriched for pathways associated with cell motility.Biological functions linked to mitochondrial stability were ranked among the top toxicological functions in both gene lists. Among the most deregulated genes, Cytochrome B5 Reductase 3 (CYB5R3) and Rap Guanine Nucleotide Exchange Factor 6 (RAPGEF6) significantly anticorrelated and correlated, respectively, with miR-342-5p in all three clinical BC datasets. Low CYB5R3 levels and high RAPGEF6 levels were significantly associated with survival, although this was not directly associated with HER2 expression., Conclusion: Our data suggest that miR-342-5p overexpression in HER2 positive BC cell lines elicits broad effects on HER2 downstream signaling, cell motility and mitochondrial stability. Together these effects may render cells less proliferative and more sensitive to cellular stress., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

The original version of this Article contained an error in the author affiliations. The affiliation of Kevin P. White with Tempus Labs, Inc., Chicago, IL, USA was inadvertently omitted.This has now been corrected in both the PDF and HTML versions of the Article.

Background: Chemotherapeutic agents such as anthracyclines and taxanes are commonly used in the neoadjuvant setting. Bevacizumab is an antibody which binds to vascular endothelial growth factor A (VEGFA) and inhibits its receptor interaction, thus obstructing the formation of new blood vessels., Methods: A phase II randomized clinical trial of 123 patients with Her2-negative breast cancer was conducted, with patients treated with neoadjuvant chemotherapy (fluorouracil (5FU)/epirubicin/cyclophosphamide (FEC) and taxane), with or without bevacizumab. Serial biopsies were obtained at time of diagnosis, after 12 weeks of treatment with FEC ± bevacizumab, and after 25 weeks of treatment with taxane ± bevacizumab. A time course study was designed to investigate the genomic landscape at the three time points when tumor DNA alterations, tumor percentage, genomic instability, and tumor clonality were assessed. Substantial differences were observed with some tumors changing mainly between diagnosis and at 12 weeks, others between 12 and 25 weeks, and still others changing in both time periods., Results: In both treatment arms, good responders (GR) and non-responders (NR) displayed significant difference in genomic instability index (GII) at time of diagnosis. In the combination arm, copy number alterations at 25 loci at the time of diagnosis were significantly different between the GR and NR. An inverse aberration pattern was also observed between the two extreme response groups at 6p22-p12 for patients in the combination arm. Signs of subclonal reduction were observed, with some aberrations disappearing and others being retained during treatment. Increase in subclonal amplification was observed at 6p21.1, a locus which contains the VEGFA gene for the protein which are targeted by the study drug bevacizumab. Of the 13 pre-treatment samples that had a gain at VEGFA, 12 were responders. Significant decrease of frequency of subclones carrying gains at 17q21.32-q22 was observed at 12 weeks, with the peak occurring at TMEM100, an ALK1 receptor signaling-dependent gene essential for vasculogenesis. This implies that cells bearing amplifications of VEGFA and TMEM100 are particularly sensitive to this treatment regime., Conclusions: Taken together, these results suggest that heterogeneity and subclonal architecture influence the response to targeted treatment in combination with chemotherapy, with possible implications for clinical decision-making and monitoring of treatment efficacy., Trial Registration: NCT00773695 . Registered 15 October 2008.

Cancers elicit an immune response by modifying the microenvironment. The immune system plays a pivotal role in cancer recognition and eradication. While the potential clinical value of infiltrating lymphocytes at the tumor site has been assessed in breast cancer, circulating cytokines - the molecules coordinating and fine-tuning immune response - are still poorly characterized. Using two breast cancer cohorts (MicMa, n = 131, DCTB, n = 28) and the multiplex Luminex platform, we measured the levels of 27 cytokines in the serum of breast cancer patients prior to treatment. We investigated the cytokine levels in relation to clinicopathological characteristics and in perspective of the tumor infiltrating immune cells predicted from the bulk mRNA expression data. Unsupervised clustering analysis of the serum cytokine levels in the MicMa cohort identified a cluster of pro-inflammatory, pro-angiogenic, and Th2-related cytokines which was associated with poor prognosis. Notably high levels of platelet derived growth factor BB (PDGF) reflected a more aggressive tumor phenotype and larger tumor size. A significant positive correlation between serum levels of interferon gamma-induced protein 10 (IP10) and its mRNA expression at the tumor site suggested that tumor-IP10-production may outflow to the bloodstream. High IP10 serum levels were associated with a worse prognosis. Finally, we found serum levels of both PDGF and IP10 associated with enrichment scores of specific tumor infiltrating immune cells. Our study suggests that monitoring cytokine circulating levels in breast cancer could be used to characterize breast cancers and the immune composition of their microenvironment through readily available biological material.

Targeted therapy for patients with HER2-positive (HER2+) breast cancer has improved overall survival, but many patients still suffer relapse and death from the disease. Intratumor heterogeneity of both estrogen receptor (ER) and HER2 expression has been proposed to play a key role in treatment failure, but little work has been done to comprehensively study this heterogeneity at the single-cell level. In this study, we explored the clinical impact of intratumor heterogeneity of ER protein expression, HER2 protein expression, and HER2 gene copy number alterations. Using combined immunofluorescence and in situ hybridization on tissue sections followed by a validated computational approach, we analyzed more than 13 000 single tumor cells across 37 HER2+ breast tumors. The samples were taken both before and after neoadjuvant chemotherapy plus HER2-targeted treatment, enabling us to study tumor evolution as well. We found that intratumor heterogeneity for HER2 copy number varied substantially between patient samples. Highly heterogeneous tumors were associated with significantly shorter disease-free survival and fewer long-term survivors. Patients for which HER2 characteristics did not change during treatment had a significantly worse outcome. This work shows the impact of intratumor heterogeneity in molecular diagnostics for treatment selection in HER2+ breast cancer patients and the power of computational scoring methods to evaluate in situ molecular markers in tissue biopsies., (© 2018 The Authors. Published by FEBS Press and John Wiley & Sons Ltd.)

The tumor microenvironment (TME) may influence both cancer progression and therapeutic response. In breast cancer, particularly in the aggressive triple-negative/basal-like subgroup, patient outcome is strongly associated with the tumor's inflammatory profile. Tumor-associated macrophages (TAMs) are among the most abundant immune cells in the TME, shown to be linked to poor prognosis and therapeutic resistance. In this study, we investigated the effect of the metastasis- and inflammation-associated microenvironmental factor S100A4 on breast cancer cells (BCCs) of different subtypes and explored their further interactions with myeloid cells. We demonstrated that extracellular S100A4 activates BCCs, particularly the basal-like subtype, to elevate secretion of pro-inflammatory cytokines. The secreted factors promoted conversion of monocytes to TAM-like cells that exhibited protumorigenic activities: stimulated epithelial-mesenchymal transition, proliferation, chemoresistance, and motility in cancer cells. In conclusion, we have shown that extracellular S100A4 instigates a tumor-supportive microenvironment, involving a network of cytokines and TAM-like cells, which was particularly characteristic for basal-like BCCs and potentiated their aggressive properties. The S100A4-BCC-TAM interaction cascade could be an important contributor to the aggressive behavior of this subtype and should be further explored for therapeutic targeting., (© 2018 The Authors. Published by FEBS Press and John Wiley & Sons Ltd.)

A high concentration of circulating vascular endothelial growth factor (VEGF) in cancer patients is associated with an aggressive tumor phenotype. Here, serum levels of 27 cytokines and blood cell counts were assessed in breast cancer patients receiving neoadjuvant chemotherapy with or without bevacizumab (Bev) in a randomized cohort of 132 patients with non-metastatic HER2-negative tumors. Cytokine levels were determined prior to treatment and at various time-points. The cytotoxic chemotherapy regimen of fluorouracil, epirubicin, and cyclophosphamide (FEC) had a profound impact on both circulating white blood cells and circulating cytokine levels. At the end of FEC treatment, the global decrease in cytokine levels correlated with the drop in white blood cell counts and was significantly greater in the patients of the Bev arm for cytokines, such as VEGF-A, IL-12, IP-10 and IL-10. Among patients who received Bev, those with pathological complete response (pCR) exhibited significantly lower levels of VEGF-A, IFN-γ, TNF-α and IL-4 than patients without pCR. This effect was not observed in the chemotherapy-only arm. Certain circulating cytokine profiles were found to correlate with different immune cell types at the tumor site. For the Bev arm patients, the serum cytokine levels correlated with higher levels of cytotoxic T cells at the end of the therapy regimen, which was indicative of treatment response. The higher response rate for Bev-treated patients and stronger correlations between serum cytokine levels and infiltrating CD8T cells merits further investigation.

Particular N-glycan structures are known to be associated with breast malignancies by coordinating various regulatory events within the tumor and corresponding microenvironment, thus implying that N-glycan patterns may be used for cancer stratification and as predictive or prognostic biomarkers. However, the association between N-glycans secreted by breast tumor and corresponding clinical relevance remain to be elucidated. We profiled N-glycans by HILIC UPLC across a discovery dataset composed of tumor interstitial fluids (TIF, n = 85), paired normal interstitial fluids (NIF, n = 54) and serum samples (n = 28) followed by independent evaluation, with the ultimate goal of identifying tumor-related N-glycan patterns in blood of patients with breast cancer. The segregation of N-linked oligosaccharides revealed 33 compositions, which exhibited differential abundances between TIF and NIF. TIFs were depleted of bisecting N-glycans, which are known to play essential roles in tumor suppression. An increased level of simple high mannose N-glycans in TIF strongly correlated with the presence of tumor infiltrating lymphocytes within tumor. At the same time, a low level of highly complex N-glycans in TIF inversely correlated with the presence of infiltrating lymphocytes within tumor. Survival analysis showed that patients exhibiting increased TIF abundance of GP24 had better outcomes, whereas low levels of GP10, GP23, GP38, and coreF were associated with poor prognosis. Levels of GP1, GP8, GP9, GP14, GP23, GP28, GP37, GP38, and coreF were significantly correlated between TIF and paired serum samples. Cross-validation analysis using an independent serum dataset supported the observed correlation between TIF and serum, for five of nine N-glycan groups: GP8, GP9, GP14, GP23, and coreF. Collectively, our results imply that profiling of N-glycans from proximal breast tumor fluids is a promising strategy for determining tumor-derived glyco-signature(s) in the blood. N-glycans structures validated in our study may serve as novel biomarkers to improve the diagnostic and prognostic stratification of patients with breast cancer., (© 2018 The Authors. Published by FEBS Press and John Wiley & Sons Ltd.)

After publication of the original article [1] the authors found that the article contained an incorrect version of Fig. 4. This does not affect the results and conclusions of the article.

Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10 -8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.