Search

BackgroundEarly changes in breast intratumor heterogeneity during neoadjuvant chemotherapy may reflect the tumor's ability to adapt and evade treatment. We investigated the combination of precision medicine predictors of genomic and MRI data towards improved prediction of recurrence free survival (RFS).MethodsA total of 100 women from the ACRIN 6657/I-SPY 1 trial were retrospectively analyzed. We estimated MammaPrint, PAM50 ROR-S, and p53 mutation scores from publicly available gene expression data and generated four, voxel-wise 3-D radiomic kinetic maps from DCE-MR images at both pre- and early-treatment time points. Within the primary lesion from each kinetic map, features of change in radiomic heterogeneity were summarized into 6 principal components.ResultsWe identify two imaging phenotypes of change in intratumor heterogeneity (p

I-SPY2 is an adaptively randomized phase 2 clinical trial evaluating novel agents in combination with standard-of-care paclitaxel followed by doxorubicin and cyclophosphamide in the neoadjuvant treatment of breast cancer. Ganitumab is a monoclonal antibody designed to bind and inhibit function of the type I insulin-like growth factor receptor (IGF-1R). Ganitumab was tested in combination with metformin and paclitaxel (PGM) followed by AC compared to standard-of-care alone. While pathologic complete response (pCR) rates were numerically higher in the PGM treatment arm for hormone receptor-negative, HER2-negative breast cancer (32% versus 21%), this small increase did not meet I-SPY's prespecified threshold for graduation. PGM was associated with increased hyperglycemia and elevated hemoglobin A1c (HbA1c), despite the use of metformin in combination with ganitumab. We evaluated several putative predictive biomarkers of ganitumab response (e.g., IGF-1 ligand score, IGF-1R signature, IGFBP5 expression, baseline HbA1c). None were specific predictors of response to PGM, although several signatures were associated with pCR in both arms. Any further development of anti-IGF-1R therapy will require better control of anti-IGF-1R drug-induced hyperglycemia and the development of more predictive biomarkers.

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 PER3VNRT 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.

In recent years, neoadjuvant treatment trials have shown that breast cancer subtypes identified on the basis of genomic and/or molecular signatures exhibit different response rates and recurrence outcomes, with the implication that subtype-specific treatment approaches are needed. Estrogen receptor-positive (ER+) breast cancers present a unique set of challenges for determining optimal neoadjuvant treatment approaches. There is increased recognition that not all ER+ breast cancers benefit from chemotherapy, and that there may be a subset of ER+ breast cancers that can be treated effectively using endocrine therapies alone. With this uncertainty, there is a need to improve the assessment and to optimize the treatment of ER+ breast cancers. While pathology-based markers offer a snapshot of tumor response to neoadjuvant therapy, non-invasive imaging of the ER disease in response to treatment would provide broader insights into tumor heterogeneity, ER biology, and the timing of surrogate endpoint measurements. In this review, we provide an overview of the current landscape of breast imaging in neoadjuvant studies and highlight the technological advances in each imaging modality. We then further examine some potential imaging markers for neoadjuvant treatment response in ER+ breast cancers.

PurposeThe phosphatidylinositol 3-kinase/Akt/mammalian target of rapamycin is a key pathway of survival and therapeutic resistance in breast cancer. We evaluated the pan-Akt inhibitor MK-2206 in combination with standard therapy in patients with high-risk early-stage breast cancer.Patients and methodsI-SPY 2 is a multicenter, phase II, open-label, adaptively randomized neoadjuvant platform trial that screens experimental therapies and efficiently identifies potential predictive biomarker signatures. Patients are categorized by human epidermal growth factor receptor 2 (HER2), hormone receptor (HR), and MammaPrint statuses in a 2 × 2 × 2 layout. Patients within each of these 8 biomarker subtypes are adaptively randomly assigned to one of several experimental therapies, including MK-2206, or control. Therapies are evaluated for 10 biomarker signatures, each of which is a combination of these subtypes. The primary end point is pathologic complete response (pCR). A therapy graduates with one or more of these signatures if and when it has an 85% Bayesian predictive probability of success in a hypothetical phase III trial, adjusting for biomarker covariates. Patients in the current report received standard taxane- and anthracycline-based neoadjuvant therapy without (control) or with oral MK-2206 135 mg/week.ResultsMK-2206 graduated with 94 patients and 57 concurrently randomly assigned controls in 3 graduation signatures: HR-negative/HER2-positive, HR-negative, and HER2-positive. Respective Bayesian mean covariate-adjusted pCR rates and percentage probability that MK-2206 is superior to control were 0.48:0.29 (97%), 0.62:0.36 (99%), and 0.46:0.26 (94%). In exploratory analyses, MK-2206 evinced a numerical improvement in event-free survival in its graduating signatures. The most significant grade 3-4 toxicity was rash (14% maculopapular, 8.6% acneiform).ConclusionThe Akt inhibitor MK-2206 combined with standard neoadjuvant therapy resulted in higher estimated pCR rates in HR-negative and HER2-positive breast cancer. Although MK-2206 is not being further developed at this time, this class of agents remains of clinical interest.

BackgroundBreast cancer patients with estrogen receptor (ER)-positive disease have a continuous long-term risk for fatal breast cancer, but the biological factors influencing this risk are unknown. We aimed to determine whether high intratumor heterogeneity of ER predicts an increased long-term risk (25 years) of fatal breast cancer.MethodsThe STO-3 trial enrolled 1780 postmenopausal lymph node-negative breast cancer patients randomly assigned to receive adjuvant tamoxifen vs not. The fraction of cancer cells for each ER intensity level was scored by breast cancer pathologists, and intratumor heterogeneity of ER was calculated using Rao's quadratic entropy and categorized into high and low heterogeneity using a predefined cutoff at the second tertile (67%). Long-term breast cancer-specific survival analyses by intra-tumor heterogeneity of ER were performed using Kaplan-Meier and multivariable Cox proportional hazard modeling adjusting for patient and tumor characteristics.ResultsA statistically significant difference in long-term survival by high vs low intratumor heterogeneity of ER was seen for all ER-positive patients (P < .001) and for patients with luminal A subtype tumors (P = .01). In multivariable analyses, patients with high intratumor heterogeneity of ER had a twofold increased long-term risk as compared with patients with low intratumor heterogeneity (ER-positive: hazard ratio [HR] = 1.98, 95% confidence interval [CI] = 1.31 to 3.00; luminal A subtype tumors: HR = 2.43, 95% CI = 1.18 to 4.99).ConclusionsPatients with high intratumor heterogeneity of ER had an increased long-term risk of fatal breast cancer. Interestingly, a similar long-term risk increase was seen in patients with luminal A subtype tumors. Our findings suggest that intratumor heterogeneity of ER is an independent long-term prognosticator with potential to change clinical management, especially for patients with luminal A tumors.

This corrects the article DOI: 10.1038/ncomms5999.

To address cancer as a multifaceted adaptive system, the increasing momentum for cross-disciplinary connectivity between cancer biologists, physical scientists, mathematicians, chemists, biomedical engineers, computer scientists, clinicians, and advocates is fueling the emergence of new scientific frontiers, principles, and opportunities within physical sciences and oncology. In parallel to highlighting the advances, challenges, and acceptance of advocates as credible contributors, we offer recommendations for addressing real world hurdles in advancing equitable partnerships among advocacy stakeholders.

This corrects the article DOI: 10.1038/ncomms5999.

The many improvements in breast cancer therapy in recent years have so lowered rates of recurrence that it is now difficult or impossible to conduct adequately powered adjuvant clinical trials. Given the many new drugs and potential synergistic combinations, the neoadjuvant approach has been used to test benefit of drug combinations in clinical trials of primary breast cancer. A recent FDA-led meta-analysis showed that pathologic complete response (pCR) predicts disease-free survival (DFS) within patients who have specific breast cancer subtypes. This meta-analysis motivated the FDA's draft guidance for using pCR as a surrogate endpoint in accelerated drug approval. Using pCR as a registration endpoint was challenged at ASCO 2014 Annual Meeting with the presentation of ALTTO, an adjuvant trial in HER2-positive breast cancer that showed a nonsignificant reduction in DFS hazard rate for adding lapatinib, a HER-family tyrosine kinase inhibitor, to trastuzumab and chemotherapy. This conclusion seemed to be inconsistent with the results of NeoALTTO, a neoadjuvant trial that found a statistical improvement in pCR rate for the identical lapatinib-containing regimen. We address differences in the two trials that may account for discordant conclusions. However, we use the FDA meta-analysis to show that there is no discordance at all between the observed pCR difference in NeoALTTO and the observed HR in ALTTO. This underscores the importance of appropriately modeling the two endpoints when designing clinical trials. The I-SPY 2/3 neoadjuvant trials exemplify this approach.

During the type-setting of the final version of the article [1] some of the additional files were swapped. The correct files are republished in this Erratum.

We recently identified a novel susceptibility variant, rs865686, for estrogen-receptor positive breast cancer at 9q31.2. Here, we report a fine-mapping analysis of the 9q31.2 susceptibility locus using 43 160 cases and 42 600 controls of European ancestry ascertained from 52 studies and a further 5795 cases and 6624 controls of Asian ancestry from nine studies. Single nucleotide polymorphism (SNP) rs676256 was most strongly associated with risk in Europeans (odds ratios [OR] = 0.90 [0.88-0.92]; P-value = 1.58 × 10(-25)). This SNP is one of a cluster of highly correlated variants, including rs865686, that spans ∼14.5 kb. We identified two additional independent association signals demarcated by SNPs rs10816625 (OR = 1.12 [1.08-1.17]; P-value = 7.89 × 10(-09)) and rs13294895 (OR = 1.09 [1.06-1.12]; P-value = 2.97 × 10(-11)). SNP rs10816625, but not rs13294895, was also associated with risk of breast cancer in Asian individuals (OR = 1.12 [1.06-1.18]; P-value = 2.77 × 10(-05)). Functional genomic annotation using data derived from breast cancer cell-line models indicates that these SNPs localise to putative enhancer elements that bind known drivers of hormone-dependent breast cancer, including ER-α, FOXA1 and GATA-3. In vitro analyses indicate that rs10816625 and rs13294895 have allele-specific effects on enhancer activity and suggest chromatin interactions with the KLF4 gene locus. These results demonstrate the power of dense genotyping in large studies to identify independent susceptibility variants. Analysis of associations using subjects with different ancestry, combined with bioinformatic and genomic characterisation, can provide strong evidence for the likely causative alleles and their functional basis.

Previous studies have suggested that polymorphisms in CASP8 on chromosome 2 are associated with breast cancer risk. To clarify the role of CASP8 in breast cancer susceptibility, we carried out dense genotyping of this region in the Breast Cancer Association Consortium (BCAC). Single-nucleotide polymorphisms (SNPs) spanning a 1 Mb region around CASP8 were genotyped in 46 450 breast cancer cases and 42 600 controls of European origin from 41 studies participating in the BCAC as part of a custom genotyping array experiment (iCOGS). Missing genotypes and SNPs were imputed and, after quality exclusions, 501 typed and 1232 imputed SNPs were included in logistic regression models adjusting for study and ancestry principal components. The SNPs retained in the final model were investigated further in data from nine genome-wide association studies (GWAS) comprising in total 10 052 case and 12 575 control subjects. The most significant association signal observed in European subjects was for the imputed intronic SNP rs1830298 in ALS2CR12 (telomeric to CASP8), with per allele odds ratio and 95% confidence interval [OR (95% confidence interval, CI)] for the minor allele of 1.05 (1.03-1.07), P = 1 × 10(-5). Three additional independent signals from intronic SNPs were identified, in CASP8 (rs36043647), ALS2CR11 (rs59278883) and CFLAR (rs7558475). The association with rs1830298 was replicated in the imputed results from the combined GWAS (P = 3 × 10(-6)), yielding a combined OR (95% CI) of 1.06 (1.04-1.08), P = 1 × 10(-9). Analyses of gene expression associations in peripheral blood and normal breast tissue indicate that CASP8 might be the target gene, suggesting a mechanism involving apoptosis.

GWAS have identified a breast cancer susceptibility locus on 2q35. Here we report the fine mapping of this locus using data from 101,943 subjects from 50 case-control studies. We genotype 276 SNPs using the 'iCOGS' genotyping array and impute genotypes for a further 1,284 using 1000 Genomes Project data. All but two, strongly correlated SNPs (rs4442975 G/T and rs6721996 G/A) are excluded as candidate causal variants at odds against >100:1. The best functional candidate, rs4442975, is associated with oestrogen receptor positive (ER+) disease with an odds ratio (OR) in Europeans of 0.85 (95% confidence interval=0.84-0.87; P=1.7 × 10(-43)) per t-allele. This SNP flanks a transcriptional enhancer that physically interacts with the promoter of IGFBP5 (encoding insulin-like growth factor-binding protein 5) and displays allele-specific gene expression, FOXA1 binding and chromatin looping. Evidence suggests that the g-allele confers increased breast cancer susceptibility through relative downregulation of IGFBP5, a gene with known roles in breast cell biology.

IntroductionWe have previously shown that a tag single nucleotide polymorphism (rs10235235), which maps to the CYP3A locus (7q22.1), was associated with a reduction in premenopausal urinary estrone glucuronide levels and a modest reduction in risk of breast cancer in women age ≤50 years.MethodsWe further investigated the association of rs10235235 with breast cancer risk in a large case control study of 47,346 cases and 47,570 controls from 52 studies participating in the Breast Cancer Association Consortium. Genotyping of rs10235235 was conducted using a custom Illumina Infinium array. Stratified analyses were conducted to determine whether this association was modified by age at diagnosis, ethnicity, age at menarche or tumor characteristics.ResultsWe confirmed the association of rs10235235 with breast cancer risk for women of European ancestry but found no evidence that this association differed with age at diagnosis. Heterozygote and homozygote odds ratios (ORs) were OR = 0.98 (95% CI 0.94, 1.01; P = 0.2) and OR = 0.80 (95% CI 0.69, 0.93; P = 0.004), respectively (P(trend) = 0.02). There was no evidence of effect modification by tumor characteristics. rs10235235 was, however, associated with age at menarche in controls (P(trend) = 0.005) but not cases (P(trend) = 0.97). Consequently the association between rs10235235 and breast cancer risk differed according to age at menarche (P(het) = 0.02); the rare allele of rs10235235 was associated with a reduction in breast cancer risk for women who had their menarche age ≥15 years (OR(het) = 0.84, 95% CI 0.75, 0.94; OR(hom) = 0.81, 95% CI 0.51, 1.30; P(trend) = 0.002) but not for those who had their menarche age ≤11 years (OR(het) = 1.06, 95% CI 0.95, 1.19, OR(hom) = 1.07, 95% CI 0.67, 1.72; P(trend) = 0.29).ConclusionsTo our knowledge rs10235235 is the first single nucleotide polymorphism to be associated with both breast cancer risk and age at menarche consistent with the well-documented association between later age at menarche and a reduction in breast cancer risk. These associations are likely mediated via an effect on circulating hormone levels.

Abstract Background First-generation molecular profiles for human breast cancers have enabled the identification of features that can predict therapeutic response; however, little is known about how the various data types can best be combined to yield optimal predictors. Collections of breast cancer cell lines mirror many aspects of breast cancer molecular pathobiology, and measurements of their omic and biological therapeutic responses are well-suited for development of strategies to identify the most predictive molecular feature sets. Results We used least squares-support vector machines and random forest algorithms to identify molecular features associated with responses of a collection of 70 breast cancer cell lines to 90 experimental or approved therapeutic agents. The datasets analyzed included measurements of copy number aberrations, mutations, gene and isoform expression, promoter methylation and protein expression. Transcriptional subtype contributed strongly to response predictors for 25% of compounds, and adding other molecular data types improved prediction for 65%. No single molecular dataset consistently out-performed the others, suggesting that therapeutic response is mediated at multiple levels in the genome. Response predictors were developed and applied to TCGA data, and were found to be present in subsets of those patient samples. Conclusions These results suggest that matching patients to treatments based on transcriptional subtype will improve response rates, and inclusion of additional features from other profiling data types may provide additional benefit. Further, we suggest a systems biology strategy for guiding clinical trials so that patient cohorts most likely to respond to new therapies may be more efficiently identified.

The term breast cancer covers many different conditions, whose clinical course ranges from indolent to aggressive. However, current practice in breast cancer prevention and care, and in breast cancer epidemiology, does not take into account the heterogeneity of the disease. A comprehensive understanding of the etiology and progression of different breast cancer subtypes would enable a more patient-centered approach to breast health care: assessing an individual's risk of getting specific subtypes of the disease, providing risk-based screening and prevention recommendations, and, for those diagnosed with the disease, tailored treatment options based on risk and timing of progression and mortality. The Athena Breast Health Network is an initiative of the five University of California medical and cancer centers to prototype this approach and to enable the development of a rapid learning system-connecting risk and outcome information from a heterogeneous patient population in real time and using new knowledge from research to continuously improve the quality of care. The Network is based on integrating clinical and research processes to create a comprehensive approach to accelerating patient-centered breast health care. Since its inception in 2009, the Network has developed a multi-site, transdisciplinary collaboration that enables the learning system. The five-campus collaboration has implemented a shared informatics platform, standardized electronic patient intake questionnaires, and common biospecimen protocols, as well as new clinical programs and multi-center research projects. The Athena Breast Health Network can serve as a model of a rapid learning system that integrates epidemiologic, behavioral, and clinical research with clinical care improvements.

PurposeWe tested the hypotheses that CHEK2*1100delC heterozygosity is associated with increased risk of early death, breast cancer-specific death, and risk of a second breast cancer in women with a first breast cancer.Patients and methodsFrom 22 studies participating in the Breast Cancer Association Consortium, 25,571 white women with invasive breast cancer were genotyped for CHEK2*1100delC and observed for up to 20 years (median, 6.6 years). We examined risk of early death and breast cancer-specific death by estrogen receptor status and risk of a second breast cancer after a first breast cancer in prospective studies.ResultsCHEK2*1100delC heterozygosity was found in 459 patients (1.8%). In women with estrogen receptor-positive breast cancer, multifactorially adjusted hazard ratios for heterozygotes versus noncarriers were 1.43 (95% CI, 1.12 to 1.82; log-rank P = .004) for early death and 1.63 (95% CI, 1.24 to 2.15; log-rank P < .001) for breast cancer-specific death. In all women, hazard ratio for a second breast cancer was 2.77 (95% CI, 2.00 to 3.83; log-rank P < .001) increasing to 3.52 (95% CI, 2.35 to 5.27; log-rank P < .001) in women with estrogen receptor-positive first breast cancer only.ConclusionAmong women with estrogen receptor-positive breast cancer, CHEK2*1100delC heterozygosity was associated with a 1.4-fold risk of early death, a 1.6-fold risk of breast cancer-specific death, and a 3.5-fold risk of a second breast cancer. This is one of the few examples of a genetic factor that influences long-term prognosis being documented in an extensive series of women with breast cancer.

Invasive lobular carcinoma (ILC) has been reported to be less responsive to neoadjuvant chemotherapy (NAC) than invasive ductal carcinoma (IDC). We sought to determine whether ILC histology indeed predicts poor response to NAC by analyzing tumor characteristics such as protein expression, gene expression, and imaging features, and by comparing NAC response rates to those seen in IDC after adjustment for these factors. We combined datasets from two large prospective NAC trials, including in total 676 patients, of which 75 were of lobular histology. Eligible patients had tumors ≥3 cm in diameter or pathologic documentation of positive nodes, and underwent serial biopsies, expression microarray analysis, and MRI imaging. We compared pathologic complete response (pCR) rates and breast conservation surgery (BCS) rates between ILC and IDC, adjusted for clinicopathologic factors. On univariate analysis, ILCs were significantly less likely to have a pCR after NAC than IDCs (11 vs. 25 %, p = 0.01). However, the known differences in tumor characteristics between the two histologic types, including hormone receptor (HR) status, HER2 status, histological grade, and p53 expression, accounted for this difference with the lowest pCR rates among HR+/HER2- tumors in both ILC and IDC (7 and 5 %, respectively). ILC which were HR- and/or HER2+ had a pCR rate of 25 %. Expression subtyping, particularly the NKI 70-gene signature, was correlated with pCR, although the small numbers of ILC in each group precluded significant associations. BCS rate did not differ between IDC and ILC after adjusting for molecular characteristics. We conclude that ILC represents a heterogeneous group of tumors which are less responsive to NAC than IDC. However, this difference is explained by differences in molecular characteristics, particularly HR and HER2, and independent of lobular histology.

Neoadjuvant chemotherapy for breast cancer allows individual tumor response to be assessed depending on molecular subtype, and to judge the impact of response to therapy on recurrence-free survival (RFS). The multicenter I-SPY 1 TRIAL evaluated patients with ≥ 3 cm tumors by using early imaging and molecular signatures, with outcomes of pathologic complete response (pCR) and RFS. The current analysis was performed using data from patients who had molecular profiles and did not receive trastuzumab. The various molecular classifiers tested were highly correlated. Categorization of breast cancer by molecular signatures enhanced the ability of pCR to predict improvement in RFS compared to the population as a whole. In multivariate analysis, the molecular signatures that added to the ability of HR and HER2 receptors, clinical stage, and pCR in predicting RFS included 70-gene signature, wound healing signature, p53 mutation signature, and PAM50 risk of recurrence. The low risk signatures were associated with significantly better prognosis, and also identified additional patients with a good prognosis within the no pCR group, primarily in the hormone receptor positive, HER-2 negative subgroup. The I-SPY 1 population is enriched for tumors with a poor prognosis but is still heterogeneous in terms of rates of pCR and RFS. The ability of pCR to predict RFS is better by subset than it is for the whole group. Molecular markers improve prediction of RFS by identifying additional patients with excellent prognosis within the no pCR group.

Abstract Introduction The detection of circulating tumour cells (CTCs) in the peripheral blood and disseminated tumour cells (DTCs) in the bone marrow are promising prognostic tools for risk stratification in early breast cancer. There is, however, a need for further validation of these techniques in larger patient cohorts with adequate follow-up periods. Methods We assayed CTCs and DTCs at primary surgery in 733 stage I or II breast cancer patients with a median follow-up time of 7.6 years. CTCs were detected in samples of peripheral blood mononuclear cells previously stored in liquid-nitrogen using a previously-developed multi-marker quantitative PCR (QPCR)-based assay. DTCs were detected in bone marrow samples by immunocytochemical analysis using anti-cytokeratin antibodies. Results CTCs were detected in 7.9% of patients, while DTCs were found in 11.7%. Both CTC and DTC positivity predicted poor metastasis-free survival (MFS) and breast cancer-specific survival (BCSS); MFS hazard ratio (HR) = 2.4 (P < 0.001)/1.9 (P = 0.006), and BCSS HR = 2.5 (P < 0.001)/2.3 (P = 0.01), for CTC/DTC status, respectively). Multivariate analyses demonstrated that CTC status was an independent prognostic variable for both MFS and BCSS. CTC status also identified a subset of patients with significantly poorer outcome among low-risk node negative patients that did not receive adjuvant systemic therapy (MFS HR 2.3 (P = 0.039), BCSS HR 2.9 (P = 0.017)). Using both tests provided increased prognostic information and indicated different relevance within biologically dissimilar breast cancer subtypes. Conclusions These results support the use of CTC analysis in early breast cancer to generate clinically useful prognostic information.

Ectonucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1) expression correlates with poor prognosis in many cancers, and we previously discovered that ENPP1 is the dominant hydrolase of extracellular cGAMP: a cancer-cell-produced immunotransmitter that activates the anticancer stimulator of interferon genes (STING) pathway. However, ENPP1 has other catalytic activities and the molecular and cellular mechanisms contributing to its tumorigenic effects remain unclear. Here, using single-cell RNA-seq, we show that ENPP1 in both cancer and normal tissues drives primary breast tumor growth and metastasis by dampening extracellular 2'3'-cyclic-GMP-AMP (cGAMP)-STING-mediated antitumoral immunity. ENPP1 loss-of-function in both cancer cells and normal tissues slowed primary tumor growth and abolished metastasis. Selectively abolishing the cGAMP hydrolysis activity of ENPP1 phenocopied ENPP1 knockout in a STING-dependent manner, demonstrating that restoration of paracrine cGAMP-STING signaling is the dominant anti-cancer mechanism of ENPP1 inhibition. Finally, ENPP1 expression in breast tumors deterministically predicated whether patients would remain free of distant metastasis after pembrolizumab (anti-PD-1) treatment followed by surgery. Altogether, ENPP1 blockade represents a strategy to exploit cancer-produced extracellular cGAMP for controlled local activation of STING and is therefore a promising therapeutic approach against breast cancer. [ABSTRACT FROM AUTHOR]