Cytosolic DNA is characteristic of chromosomally unstable metastatic cancer cells, resulting in constitutive activation of the cGAS–STING innate immune pathway. How tumors co-opt inflammatory signaling while evading immune surveillance remains unknown. Here, we show that the ectonucleotidase ENPP1 promotes metastasis by selectively degrading extracellular cGAMP, an immune-stimulatory metabolite whose breakdown products include the immune suppressor adenosine. ENPP1 loss suppresses metastasis, restores tumor immune infiltration, and potentiates response to immune checkpoint blockade in a manner dependent on tumor cGAS and host STING. Conversely, overexpression of wild-type ENPP1, but not an enzymatically weakened mutant, promotes migration and metastasis, in part through the generation of extracellular adenosine, and renders otherwise sensitive tumors completely resistant to immunotherapy. In human cancers, ENPP1 expression correlates with reduced immune cell infiltration, increased metastasis, and resistance to anti–PD-1/PD-L1 treatment. Thus, cGAMP hydrolysis by ENPP1 enables chromosomally unstable tumors to transmute cGAS activation into an immune-suppressive pathway.

Significance:

Chromosomal instability promotes metastasis by generating chronic tumor inflammation. ENPP1 facilitates metastasis and enables tumor cells to tolerate inflammation by hydrolyzing the immunotransmitter cGAMP, preventing its transfer from cancer cells to immune cells.

This article is highlighted in the In This Issue feature, p. 995

Chromosomal instability (CIN) is a hallmark of human cancer, and it is associated with metastasis, immune evasion, and therapeutic resistance (1–5). In addition to the generation of chromosome copy-number heterogeneity, which serves as a substrate for natural selection, CIN also promotes tumor progression by inducing chronic inflammatory signaling, leading to increased cancer cell migration and invasion (1, 6). Chromosome segregation errors lead to the formation of micronuclei (7, 8). Micronuclear envelopes are highly rupture-prone, often exposing genomic double-stranded DNA (dsDNA) to the cytosol (1, 9–12). Cytosolic dsDNA is sensed by cGAS, which, upon binding to its substrate, catalyzes the formation of the cyclic dinucleotide cGAMP (13). A potent immune-stimulatory molecule, cGAMP promotes inflammatory signaling in a manner dependent on its downstream effector STING (14, 15).

Given the pervasive nature of CIN in human cancer (4), tumor cells must cope with the presence of persistent inflammatory signaling arising from cGAS sensing of cytosolic dsDNA. The activation of cGAS–STING has cell-autonomous and cell-nonautonomous consequences, and therefore cancer cells must mitigate the effects of this inflammatory pathway at multiple levels. One mechanism by which chromosomally unstable cancer cells have evolved to cope with chronic cGAS–STING activation is through silencing of downstream type I IFN signaling while selecting for NFκB-dependent activity to spread to distant organs (1). In line with this, an analysis of STING (Tmem173) expression in The Cancer Genome Atlas (TCGA) database found that tumor STING primarily correlates with NFκB-dependent transcriptional programs, such as the senescence-associated secretory phenotype, rather than IFN-stimulated genes (16). The switch from type I IFN to NFκB-predominant signaling downstream of STING has been proposed to enable cancer cells to simultaneously evade immune surveillance—arising from IFN signaling—while activating noncanonical NFκB-dependent migratory programs, culminating in metastatic progression (1, 6).

In addition to its cell-intrinsic effects, cGAMP is readily exported to the extracellular space where it can promote antitumor immune responses by activating STING in host cells present in the tumor microenvironment (17–19). Unlike cancer cells, host cells respond to STING activation by inducing a robust type I IFN signaling central to a productive cell-mediated immunity. How tumor cells with CIN evolve to eschew the deleterious effects of paracrine cGAMP signaling remains poorly understood. Understanding the adaptive mechanisms employed by cancer cells to evade immune surveillance in response to chronic inflammatory signaling represents an attractive therapeutic opportunity to selectively target tumor cells with CIN, by unmasking them to the immune system, while sparing normal cells devoid of cytosolic dsDNA.

ENPP1 Is Upregulated in Cells with CIN

To investigate the status of cGAS–STING signaling in cancer cells with CIN, we used a human triple-negative breast cancer (TNBC) cell line, MDA-MB-231, that was engineered to exhibit different rates of CIN through overexpression of the kinesin-13 proteins KIF2B or MCAK, or the dominant-negative mutant isoform of MCAK (dnMCAK; refs. 1, 20). We have previously shown that in these otherwise isogenic cell lines, expression of dnMCAK promotes increased chromosome missegregation, leading to the formation of micronuclei, chronic activation of cGAS–STING signaling, and increased metastasis (1). In addition, we used three syngeneic metastasis-competent mouse models of TNBC (4T1 and E0771) and colorectal cancer (CT26). All three models exhibited evidence for CIN, including the presence of chromosome missegregation during anaphase and a preponderance of micronuclei with robust cGAS staining indicative of cytosolic exposure of genomic dsDNA (Supplementary Fig. S1A and S1B). To test if cGAS localization to micronuclei also led to pathway activation, we measured cGAMP levels in total cell lysates of 4T1 cells and upon CRISPR/Cas9-mediated knockout (KO) of Cgas. Loss of cGAS resulted in a significant reduction in the levels of the cyclic dinucleotide, in line with constitutive activation of the pathway in chromosomally unstable cells (Supplementary Fig. S1C and S1D). Furthermore, cGAMP levels were nearly 15-fold higher in conditioned media after 24 hours as compared with cell lysates when both were normalized to cell counts (Supplementary Fig. S1D), suggesting that cGAMP is readily exported from cancer cells, as previously proposed (17–19).

To determine how chromosomally unstable cells adapt to ongoing cGAMP production, we performed pairwise differential expression analysis of otherwise isogenic CINhi (highly metastatic) and CINlo (poorly metastatic) MDA-MB-231 cells. Among the large number of differentially expressed genes, ENPP1 stood out because of its reported role as a negative regulator of cGAMP (21). An ectonucleotidase with a single transmembrane domain, ENPP1 localizes to the plasma membrane with its catalytic site facing the extracellular space where it has been proposed to selectively hydrolyze the extracellular pool of cGAMP (19). Both ENPP1 messenger and protein levels were markedly increased in CINhi cells compared with their CINlo counterparts (Supplementary Fig. S1E and S1F, Log2 fold change = 1.23, FDRq = 8.4 × 10−4). Staining of MDA-MB-231 CINhi cells using an anti-ENPP1 antibody revealed strong membrane localization that was abolished upon shRNA-mediated depletion (Fig. 1A). A similar pattern of cell membrane staining was seen in orthotopically transplanted tumors, where specificity was validated using shRNA-mediated depletion (Fig. 1B; Supplementary Fig. S2A and S2B).

Figure 1.

ENPP1 promotes metastasis of chromosomally unstable tumors. A, Representative immunofluorescence images of control and ENPP1-depleted MDA-MB-231 CINhi cells stained with DAPI (DNA) and anti-ENPP1 antibody; scale bar, 5 μm. B, IHC of an orthotopically transplanted MBA-MB-231 tumor using anti-ENPP1 antibody. C,Enpp1 mRNA expression in various stages of lung adenocarcinoma progression; bars represent mean ± SEM. DTCs, disseminated tumor cells. D, Extracellular-to-intracellular cGAMP ratio in 4T1, CT26, and E0771 cells; bars represent median, n = 10 independent experiments; *, P < 0.05; **, P < 0.01, two-sided Mann–Whitney test. E, Overall survival of animals that were orthotopically transplanted by control and Enpp1-KO 4T1 tumors followed by tumor resection 7 days later (n = 15 animals per condition), and significance tested using log-rank test. F, Left, quantification of surface lung metastases after tail-vein injection of control and Enpp1-KO 4T1 cells; bars represent median, n = 13–15 animals per condition; ****, P < 0.0001, two-sided Mann–Whitney test. Right, representative hematoxylin and eosin–stained lungs from animals injected with control and Enpp1-KO 4T1 cells; scale bar, 3 mm. G, Surface lung metastases after tail-vein injection of eGFP- and eGFP-ENPP1–expressing CT26 cells; bars represent median, n = 15 animals per condition; ****, P < 0.0001, two-sided Mann–Whitney test. FPKM, fragments per kilobase of transcript per million mapped reads.

Figure 1.

ENPP1 promotes metastasis of chromosomally unstable tumors. A, Representative immunofluorescence images of control and ENPP1-depleted MDA-MB-231 CINhi cells stained with DAPI (DNA) and anti-ENPP1 antibody; scale bar, 5 μm. B, IHC of an orthotopically transplanted MBA-MB-231 tumor using anti-ENPP1 antibody. C,Enpp1 mRNA expression in various stages of lung adenocarcinoma progression; bars represent mean ± SEM. DTCs, disseminated tumor cells. D, Extracellular-to-intracellular cGAMP ratio in 4T1, CT26, and E0771 cells; bars represent median, n = 10 independent experiments; *, P < 0.05; **, P < 0.01, two-sided Mann–Whitney test. E, Overall survival of animals that were orthotopically transplanted by control and Enpp1-KO 4T1 tumors followed by tumor resection 7 days later (n = 15 animals per condition), and significance tested using log-rank test. F, Left, quantification of surface lung metastases after tail-vein injection of control and Enpp1-KO 4T1 cells; bars represent median, n = 13–15 animals per condition; ****, P < 0.0001, two-sided Mann–Whitney test. Right, representative hematoxylin and eosin–stained lungs from animals injected with control and Enpp1-KO 4T1 cells; scale bar, 3 mm. G, Surface lung metastases after tail-vein injection of eGFP- and eGFP-ENPP1–expressing CT26 cells; bars represent median, n = 15 animals per condition; ****, P < 0.0001, two-sided Mann–Whitney test. FPKM, fragments per kilobase of transcript per million mapped reads.

Close modal

We next surveyed ENPP1 expression across mouse cancer cell lines and found that 4T1 had the highest mRNA expression levels when compared with CT26 and E0771. Interestingly, E0771.LMB, a more metastatic E0771 derivative (22), had significantly increased levels of Enpp1 mRNA (Supplementary Fig. S2C), suggesting that ENPP1 might be highly expressed in metastatic cancer cells which also frequently exhibit high rates of chromosome missegregation (1). In line with this, Enpp1 mRNA was significantly elevated in 4T1 cells derived from lung metastases compared with the parental cell line (Supplementary Fig. S2D). We next analyzed Enpp1 expression in the various stages of tumorigenesis in a genetically engineered mouse model of lung adenocarcinoma driven by oncogenic KRASG12D and loss of Trp53 (23). In this model, gene expression of barcoded cells was analyzed in the normal lung, benign hyperplasia, primary tumors with various metastatic proclivities, disseminated tumor cells, and overt metastases. Strikingly, mRNA levels of Enpp1 exhibited a stepwise increase during the progression from normal tissue to primary tumors to metastases. Furthermore, primary tumors that seeded metastases had higher Enpp1 expression compared with their nonmetastatic counterparts (Fig. 1C). ENPP1 protein expression mirrored this trend in orthotopically transplanted TNBC tumors, with increased levels observed selectively in tumor cells that have invaded nearby intramammary lymph nodes (Fig. 1B).

ENPP1 Promotes Cancer Metastasis

To directly test the role of ENPP1 in metastasis, we performed CRISPR/Cas9 KO of Enpp1 in 4T1 cells (Supplementary Fig. S2E). We also overexpressed wild-type (WT) ENPP1 or an enzymatically weakened mutant isoform containing a threonine-to-alanine substitution in the catalytic domain (T238A; ref. 24) in CT26 and E0771 cells which express low baseline levels of this enzyme (Supplementary Fig. S2C). As expected, loss of ENPP1 led to a significant increase in the extracellular-to-intracellular cGAMP ratio (Fig. 1D). Conversely, overexpression of WT ENPP1, but not the enzymatically weakened mutant, led to a reduction in the extracellular-to-intracellular cGAMP ratio in CT26 and E0771 cells (Fig. 1D). Enpp1 KO did not affect cellular proliferation in vitro or primary tumor growth in vivo when 4T1 cells were orthotopically transplanted in the mammary fat pad (Supplementary Fig. S2F and S2G). We then transplanted parental and Enpp1-KO 4T1 cells into BALB/c hosts, either through tail-vein inoculation or orthotopic transplantation followed by primary tumor excision. Loss of ENPP1 led to significantly longer overall survival and a marked reduction in local tumor recurrence and metastasis regardless of whether cells were introduced directly into the tail vein or orthotopically transplanted followed by surgical excision of the primary tumor (Fig. 1E and F; Supplementary Fig. S2H–S2K). Conversely, overexpression of WT ENPP1 led to a significant increase in the number of surface lung metastases upon tail-vein inoculation of CT26 cells (Fig. 1G).

To further examine whether ENPP1 disrupts paracrine tumor-to-host cGAMP transfer during metastatic progression, we overexpressed WT ENPP1 or ENPP1T328A in E0771 and quantified metastatic dissemination using bioluminescence imaging (BLI). Only WT ENPP1—but not ENPP1T328A—led to increased metastatic dissemination (Fig. 2A). Importantly, the role of ENPP1 in metastasis was dependent on host STING, as both control and WT ENPP1–overexpressing cells had similar metastatic proclivity when transplanted into MPYS−/− (Tmem173−/−) hosts (Fig. 2A). Collectively, these results suggest that ENPP1 promotes metastatic progression through extracellular cGAMP hydrolysis, preventing protective STING activation in host cells.

Figure 2.

ENPP1 promotes extracellular adenosine production. A, Left, total BLI of WT or Tmem173−/− animals inoculated with E0771 cells expressing WT or enzymatically weakened ENPP1 (T328A); bars represent median, n = 13–15 mice per group for the WT animals and 11–12 for the Tmem173−/− animals; *, P < 0.05, Welch t test. B, Schematic showing the generation of adenosine from extracellular cGAMP and ATP hydrolysis. C, Normalized adenosine concentration (per 107 cells after 16-hour incubation in serum-free media) in conditioned media of control, Cgas-KO, and Enpp1-KO 4T1 cells; bars represent mean ± SEM, n = 4 independent experiments; *, P < 0.05; **, P < 0.01 two-sided t test. D, Percent wound remaining after 24 hours in control, Cgas-KO, and Enpp1-KO 4T1 cells treated with cGAMP or cGAMP and the adenosine receptor blocker PSB1115. E,Nt5e and Entpd1 mRNA expression in various stages of lung adenocarcinoma progression; bars, mean ± SEM. FPKM, fragments per kilobase of transcript per million mapped reads. F, Surface lung metastases after tail-vein injection of control, Enpp1-KO, Nt5e-KO, and Enpp1/Nt5e double KO 4T1 cells; bars, median, n = 15 animals per condition; ***, P < 0.001, two-sided Mann–Whitney test.

Figure 2.

ENPP1 promotes extracellular adenosine production. A, Left, total BLI of WT or Tmem173−/− animals inoculated with E0771 cells expressing WT or enzymatically weakened ENPP1 (T328A); bars represent median, n = 13–15 mice per group for the WT animals and 11–12 for the Tmem173−/− animals; *, P < 0.05, Welch t test. B, Schematic showing the generation of adenosine from extracellular cGAMP and ATP hydrolysis. C, Normalized adenosine concentration (per 107 cells after 16-hour incubation in serum-free media) in conditioned media of control, Cgas-KO, and Enpp1-KO 4T1 cells; bars represent mean ± SEM, n = 4 independent experiments; *, P < 0.05; **, P < 0.01 two-sided t test. D, Percent wound remaining after 24 hours in control, Cgas-KO, and Enpp1-KO 4T1 cells treated with cGAMP or cGAMP and the adenosine receptor blocker PSB1115. E,Nt5e and Entpd1 mRNA expression in various stages of lung adenocarcinoma progression; bars, mean ± SEM. FPKM, fragments per kilobase of transcript per million mapped reads. F, Surface lung metastases after tail-vein injection of control, Enpp1-KO, Nt5e-KO, and Enpp1/Nt5e double KO 4T1 cells; bars, median, n = 15 animals per condition; ***, P < 0.001, two-sided Mann–Whitney test.

Close modal

Extracellular cGAMP Hydrolysis by ENPP1 Generates Adenosine

We next explored the fate of tumor-derived extracellular cGAMP and asked whether the breakdown products of this metabolite might contribute to the production of extracellular adenosine, an immune-suppressive and tumor-promoting metabolite (25). cGAMP hydrolysis by ENPP1 leads to the formation of AMP and GMP. AMP can be subsequently hydrolyzed into adenosine by NT5E (also known as CD73; Fig. 2B). Measuring adenosine in conditioned media is technically challenging given the presence of enzymes that either degrade this nucleoside [adenosine deaminase (ADA)] or promote its cellular reuptake (Supplementary Fig. S3A). To overcome these challenges, we added serum-free media to 4T1 cells in the presence of erythro-9-(2-hydroxy-3-nonyl)adenine (EHNA), an ADA inhibitor, along with dipyridamole and 6-S-[(4-Nitrophenyl)methyl]-6-thioinosine (NBMPR), which prevent cellular reuptake of adenosine (Supplementary Fig. S3A; ref. 26). Extracellular adenosine levels—as assessed by LC/MS in conditioned media—were reduced by up to 40% upon KO of either Cgas or Enpp1 (Fig. 2C). Using an orthogonal approach, we added exogenous cGAMP to 4T1 cells and used a fluorescence-based method to detect hydrogen peroxide (H2O2) resulting from the oxidation of hypoxanthine, a breakdown product of adenosine (Supplementary Fig. S3A). By comparing fluorescence in the presence and absence of EHNA, we were able to assess relative contribution from adenosine degradation toward H2O2 production and observed a concentration-dependent increase in H2O2 production after the addition of exogenous cGAMP (Supplementary Fig. S3B), suggesting that this cyclic dinucleotide can be readily converted into adenosine in the extracellular environment.

Through its ability to bind extracellular adenosine receptors in both tumor and immune cells, adenosine promotes cancer cell migration and imparts potent immune-suppressive effects (25, 27). Interestingly, KO of either Cgas or Enpp1 in 4T1 cells led to a significant reduction in cellular migration, whereas exogenous addition of cGAMP to the conditioned media rescued migration only in Cgas-KO—but not Enpp1-KO—tumor cells (Fig. 2D). The effect of cGAMP was dependent on activation of the extracellular adenosine receptors and was abolished upon the addition of PSB1115, an inhibitor of the adenosine A2B receptor on cancer cells (Fig. 2D). Conversely, overexpression of WT ENPP1—but not ENPP1T328A—in E0771 or CT26 cells led to increased migration, an effect that was abolished upon treatment of the conditioned media with ADA (Supplementary Fig. S3C and S3D).

In addition to cGAMP hydrolysis by ENPP1, ATP hydrolysis by either ENPP1 or ENTPD1 (also known as CD39) is considered to be a major source of extracellular AMP. Interestingly, in the lung adenocarcinoma tumorigenesis model, expression of mouse Nt5e mirrored that of Enpp1 in that it progressively increased from normal tissues to primary tumors to metastases (Fig. 2E). On the contrary, Entpd1 expression followed the opposite trend with the lowest expression levels observed in metastatic lesions (Fig. 2E). These opposing trends suggest that although ATP hydrolysis might represent a major source of extracellular adenosine in primary tumors, the relative contribution from cGAMP hydrolysis as an adenosine source increases along with metastatic progression. In line with this finding, KO of either Enpp1 or Nt5e in 4T1 cells led to a significant reduction in the number of lung metastases in a manner commensurate with combined loss of both enzymes (Fig. 2F; Supplementary Fig. S3E).

We had recently shown that tumor cell–intrinsic STING activation by intracellular cGAMP can also promote cellular migration and metastasis (1). To test the relative contributions of tumor cell STING activation and extracellular cGAMP hydrolysis by ENPP1, we assessed metastatic potential of control, Enpp1-KO, Tmem173-KO, and Enpp1/Tmem173 double KO 4T1 cells by comparing animal survival after tail-vein inoculation. Loss of either ENPP1 or STING in tumor cells led to reduced metastasis and lifespan extension, and their combined KO led to an additive effect (Supplementary Fig. S3F). Collectively, this suggests that intracellular cGAMP-dependent STING activation and extracellular cGAMP hydrolysis by ENPP1 independently contribute to metastatic progression. Furthermore, these results also indicate that the impact of ENPP1 on metastasis is mediated through activation of host—but not tumor cell—STING.

ENPP1 Promotes Tumor Immune Evasion

We next examined the effect of ENPP1 loss on tumor immune infiltration using shRNA-mediated depletion or CRISPR/Cas9 KO in CINhi MDA-MB-231 orthotopic xenografts and 4T1 metastatic allografts, respectively. Loss of ENPP1 led to increased tumor necrosis and enhanced infiltration of natural killer (NK) cells in MDA-MB-231 tumors (Supplementary Fig. S4A and S4B), in line with previous reports demonstrating a role for cGAMP transfer in activating NK cells (17). In the 4T1 model, metastatic lesions formed from Enpp1-KO 4T1 cells exhibited significant infiltration by CD45+ cells and an approximately 3- to 5-fold enrichment with CD8+ T cells compared with WT counterparts (Fig. 3A and B). Flow cytometry–based immune profiling of dissociated lungs revealed a significant increase in CD45+ cells, CD4+ T cells, and granulocytic CD11b+Ly6G+ cells as compared with controls (Fig. 3C; Supplementary Fig. S4C). Unlike our IHC-based results, we did not observe an absolute enrichment for CD8+ T cells in the injected lungs using flow cytometry; however, there was a significant increase in PD-1+ subpopulations of CD3+CD8+ and CD3+CD4+ cells (Fig. 3C). The overall preponderance of granulocytic cells was notable, given that Enpp1-KO tumors had higher levels of GM-CSF as measured using ELISA-based assays (Fig. 3D). Collectively, these findings suggest that in addition to lymphocytes, granulocytic cells may also play a role in restricting metastatic colonization of Enpp1-KO cells, in line with previous reports showing an antitumor and proinflammatory effect of CD11b+Ly6G+ cells (28–30).

Figure 3.

ENPP1 reduces tumor immune infiltration. A, Representative IHC of control and ENPP1-KO TNBC lung metastases stained using an anti-CD45 antibody. B, The number of metastasis-infiltrating CD8+ T cells (left) and representative IHC of control Enpp1-KO TNBC lung metastases stained using anti-CD8 antibody (right); bars, median, n = 13–31 metastases; ****, P < 0.0001, two-sided Mann–Whitney test. C, Percentage of CD45+, CD11b+Ly6G+, CD4+, and CD8+ cells out of the total cells as well as the percentage of PD-1+ cells out of the CD3+CD4+ and CD3+CD8+ cells obtained from dissociated lungs after injection with control or Enpp1-KO 4T1 cells, n = 5 animals per group. ***, P < 0.001, two-sided Mann-Whitney test. D, GM-CSF levels measured in orthotopically transplanted control and ENPP1-KO tumors; bars, median, n = 15 tumors per condition; **, P < 0.01, two-sided Mann–Whitney test. E, Percentage of CD8+ T cells, CD4+ T cells (and the PD-1+ and CD44+ fractions of thereof), and NK cells obtained from dissociated subcutaneously transplanted control and ENPP1-expressing CT26 tumors, n = 5 animals per group, bars represent median; *, P < 0.05. Treg, regulatory T cell.

Figure 3.

ENPP1 reduces tumor immune infiltration. A, Representative IHC of control and ENPP1-KO TNBC lung metastases stained using an anti-CD45 antibody. B, The number of metastasis-infiltrating CD8+ T cells (left) and representative IHC of control Enpp1-KO TNBC lung metastases stained using anti-CD8 antibody (right); bars, median, n = 13–31 metastases; ****, P < 0.0001, two-sided Mann–Whitney test. C, Percentage of CD45+, CD11b+Ly6G+, CD4+, and CD8+ cells out of the total cells as well as the percentage of PD-1+ cells out of the CD3+CD4+ and CD3+CD8+ cells obtained from dissociated lungs after injection with control or Enpp1-KO 4T1 cells, n = 5 animals per group. ***, P < 0.001, two-sided Mann-Whitney test. D, GM-CSF levels measured in orthotopically transplanted control and ENPP1-KO tumors; bars, median, n = 15 tumors per condition; **, P < 0.01, two-sided Mann–Whitney test. E, Percentage of CD8+ T cells, CD4+ T cells (and the PD-1+ and CD44+ fractions of thereof), and NK cells obtained from dissociated subcutaneously transplanted control and ENPP1-expressing CT26 tumors, n = 5 animals per group, bars represent median; *, P < 0.05. Treg, regulatory T cell.

Close modal

We next assessed the impact of WT ENPP1 overexpression on subcutaneously transplanted CT26 tumors. Expectedly, exogenous expression of Enpp1 led to reduced CD8+ T cells and NK cells as well as the proportion of PD-1+ CD8+ and CD4+ T cells. In line with these findings, there was a decrease in the proportion of CD44+ T cells, suggesting reduced T-cell activation (Fig. 3E; Supplementary Fig. S5A). The fraction of FOXP3+ T regulatory cells remained constant with a significant reduction in the CD8+:FOXP3+ ratio noted, consistent with an immunosuppressive response (Fig. 3E; Supplementary Fig. S5A).

To determine whether the increased immune infiltration upon ENPP1 loss was dependent on tumor cell–derived cGAMP, we performed population-level depletion of Cgas using CRISPR KO and found a trend toward reduced CD45+ cell and CD8+ T-cell infiltration when cGAS was codepleted in Enpp1-KO 4T1 cells (Supplementary Fig. S5B–S5D). We posit that the lack of complete rescue might be due to the residual fraction of cells with functional cGAS or alternative sources of cGAMP in the tumor microenvironment. Nonetheless, these data suggest that ENPP1 dampens proinflammatory tumor immune infiltration through extracellular cGAMP hydrolysis.

ENPP1 Inhibition Potentiates Response to Immune Checkpoint Blockade Therapy

We then asked whether targeting ENPP1 might represent a selective therapeutic vulnerability to sensitize otherwise resistant chromosomally unstable tumors to immune checkpoint blockade (ICB) therapy. Interestingly, baseline Enpp1 mRNA expression levels in the three mouse cancer cell lines (Supplementary Fig. S2C) mirrored their previously reported sensitivities to ICB therapy, with CT26 and E0771 being considered responsive to ICB treatment in stark contrast to the highly resistant 4T1 model (18, 31). We postulated that Enpp1 KO would render 4T1 tumors responsive to ICB therapy, whereas its overexpression would confer resistance to otherwise sensitive CT26 and E0771 tumors (Fig. 4A; Supplementary Fig. S6A). Luciferase-expressing 4T1 cells were orthotopically transplanted into the mammary fat pad of BALB/c mice, and primary tumor growth was assessed over the span of 25 days (Fig. 4B; Supplementary Fig. S6B and S6C). Animals were treated with combined ICB [anti–PD-1 (aPD-1) and anti-CTLA4 (aCTLA4)] starting at day 6 after tumor cell inoculation for four doses followed by maintenance aCTLA4 treatment every 3 days for four additional doses. Enpp1-KO tumors, derived from two independent KO lines, exhibited reduced tumor growth rates compared with their WT counterparts when both were treated with combined ICB therapy, leading to significantly prolonged survival of the former (Fig. 4B and C; Supplementary Fig. S6C). Importantly, Cgas KO in Enpp1-KO cells diminished the responsiveness of 4T1 tumors, leading to significantly shorter survival (Fig. 4C). Notably, loss of cGAS did not lead to a full rescue of tumor response seen upon Enpp1 KO, suggesting that the hydrolysis of either non–tumor-derived cGAMP or ATP might contribute to the immune evasion phenotype mediated by ENPP1.

Figure 4.

ENPP1 promotes resistance to ICB therapy. A, Schematic diagram of immunotherapy experiments. B, Growth curves of control, Enpp1-KO, Cgas-KO, and Enpp1/Cgas double KO orthotopically transplanted tumors 4T1 upon treatment with combined ICB or corresponding isotype controls; data points, mean ± SEM, n = 15 animals per group; ***, P < 0.001; ****, P < 0.0001, two-sided t test. C, Survival of animals after orthotopic transplantation with control, Enpp1-KO, Cgas-KO, or Enpp1/Cgas double KO 4T1 cells treated with combined ICB or corresponding isotype controls; significance tested using log-rank test; *, P < 0.05; ***, P < 0.001, n = 15 animals per group. D, Survival of BALB/c mice injected with eGFP- or eGFP-ENPP1–expressing CT26 cells, treated with combined ICB or isotype controls, n = 15 animals per group; significance tested using log-rank test; ***, P < 0.001. E, Survival of WT or Tmem173−/− C57BL/6 mice orthotopically transplanted with eGFP- or eGFP-ENPP1–expressing E0771 tumors treated with combined ICB or isotype control antibodies, n = 10 and 4–5 animals per group for the WT and Tmem173−/− C57BL/6 mice, respectively; significance tested using log-rank test, *, P < 0.05; ***, P < 0.001. ns, not significant.

Figure 4.

ENPP1 promotes resistance to ICB therapy. A, Schematic diagram of immunotherapy experiments. B, Growth curves of control, Enpp1-KO, Cgas-KO, and Enpp1/Cgas double KO orthotopically transplanted tumors 4T1 upon treatment with combined ICB or corresponding isotype controls; data points, mean ± SEM, n = 15 animals per group; ***, P < 0.001; ****, P < 0.0001, two-sided t test. C, Survival of animals after orthotopic transplantation with control, Enpp1-KO, Cgas-KO, or Enpp1/Cgas double KO 4T1 cells treated with combined ICB or corresponding isotype controls; significance tested using log-rank test; *, P < 0.05; ***, P < 0.001, n = 15 animals per group. D, Survival of BALB/c mice injected with eGFP- or eGFP-ENPP1–expressing CT26 cells, treated with combined ICB or isotype controls, n = 15 animals per group; significance tested using log-rank test; ***, P < 0.001. E, Survival of WT or Tmem173−/− C57BL/6 mice orthotopically transplanted with eGFP- or eGFP-ENPP1–expressing E0771 tumors treated with combined ICB or isotype control antibodies, n = 10 and 4–5 animals per group for the WT and Tmem173−/− C57BL/6 mice, respectively; significance tested using log-rank test, *, P < 0.05; ***, P < 0.001. ns, not significant.

Close modal

We next asked whether overexpression of ENPP1 would confer ICB therapy resistance in otherwise sensitive CT26 and E0771 tumors (Fig. 4A; Supplementary Fig. S6A). CT26-bearing mice were treated with combined ICB starting at day 6 for a total of five doses. Strikingly, not only did eGFP-ENPP1 expression lead to increased metastasis and reduced survival of isotype control–treated mice, it also rendered this model completely resistant to combined ICB (Fig. 4D). Conversely, eGFP-expressing CT26 tumors were responsive to combined ICB, with 60% of animals surviving for more than 140 days. Similarly, overexpression of eGFP-ENPP1 in orthotopically transplanted E0771 tumors led to their resistance upon three treatments of aPD-1 antibody, wherein 50% of animals bearing eGFP-expressing E0771 tumors underwent a durable complete response compared with 0% of their eGFP-ENPP1–expressing tumor-bearing counterparts (Fig. 4E; Supplementary Fig. S6D). Importantly, the difference in response between eGFP- and eGFP-ENPP1–expressing tumors was abolished when they were transplanted in MPYS−/− (Tmem173−/−) hosts (Fig. 4E; Supplementary Fig. S6D). Collectively, these results suggest that ENPP1 inhibition represents an attractive therapeutic strategy to potentiate the response of chromosomally unstable cancers cells to ICB therapy.

ENPP1 Is Associated with Metastasis in Human Cancer

We next sought to interrogate the role of ENPP1 in human cancers by analyzing Enpp1 mRNA and protein expression in a large number of tumors from various tissues of origin. Enpp1 mRNA was investigated in tumors found in the TCGA, an independent set of primary and metastatic tumors, two separate sarcoma cohorts, and in tumor-derived organoids. ENPP1 protein expression was also analyzed in three independent breast cancer cohorts, including two estrogen receptor–negative (ER) cohorts (n = 223 and 91) and one ER-positive (ER+) cohort (n = 115), as well as in mucosal melanoma primary and metastatic tumors (n = 24).

Enpp1 mRNA expression was highly variable across cancer types found in the TCGA, with the highest expression levels observed in sarcomas and liver, breast, and thyroid cancers (Supplementary Fig. S7A). Elevated Enpp1 mRNA was associated with reduced overall survival in multiple tumor types including breast cancer, irrespective of its hormone receptor status (Supplementary Fig. S7B–S7D). To determine if Enpp1 expression was associated with metastatic progression, we first compared Enpp1 expression levels in a large number of primary and metastatic tumor samples as well as in a collection of tumor-derived organoids. In both cases, Enpp1 mRNA was higher in metastases compared with primary tumors (Fig. 5A; Supplementary Fig. S8A). When metastatic tumors were stratified by tissue site, we found liver and brain metastases to contain the highest expression levels of Enpp1 (Fig. 5A). We next surveyed ENPP1 protein expression in primary and metastatic mucosal melanoma tumors. Unlike cutaneous melanoma, mucosal melanoma is characterized by elevated CIN, reduced tumor mutational burden, and increased resistance to ICB (32, 33). In these tumors, membrane ENPP1 expression was seen in both tumor cells and the stroma, and this pattern was evenly distributed across primary tumor samples. Conversely, metastases displayed significantly increased cancer cell–specific ENPP1 staining (Fig. 5B). Tumor cell–intrinsic ENPP1 protein expression was most remarkable in lymph node metastases where cancer cell clusters displayed strong ENPP1 expression in an otherwise immune cell replete microenvironment (Fig. 5C and D).

Figure 5.

ENPP1 expression is associated with metastasis in human cancer. A, ENPP1 expression across primary and metastatic tumors, stratified by the site of metastasis, n = 180 tumors for primary tumors and 331 tumors for metastases; bars, median; *, P < 0.05; **, P < 0.01; ***, P < 0.001. B, Percentage of patients with mucosal melanoma with tumor-specific or stromal-specific ENPP1 staining patterns in primary as well as metastatic mucosal melanoma human tumor samples. *, P < 0.05, χ2 test. C and D, Representative immunofluorescence images of low (C) and high (D) magnification images of lymph node metastases from mucosal melanoma stained using DAPI (DNA) and anti-ENPP1 antibody showing selective membrane staining of ENPP1 on metastatic cancer cells. Scale bar, 1 mm (C) and 50 μm (D). E, Representative images of human TNBCs stained using anti-ENPP1 antibody; scale bar, 100 μm. F, Distant metastasis-free survival (DMFS) in patients with TNBC stratified based on their Enpp1 and Cgas expression; n = 159, significance tested using log-rank test.

Figure 5.

ENPP1 expression is associated with metastasis in human cancer. A, ENPP1 expression across primary and metastatic tumors, stratified by the site of metastasis, n = 180 tumors for primary tumors and 331 tumors for metastases; bars, median; *, P < 0.05; **, P < 0.01; ***, P < 0.001. B, Percentage of patients with mucosal melanoma with tumor-specific or stromal-specific ENPP1 staining patterns in primary as well as metastatic mucosal melanoma human tumor samples. *, P < 0.05, χ2 test. C and D, Representative immunofluorescence images of low (C) and high (D) magnification images of lymph node metastases from mucosal melanoma stained using DAPI (DNA) and anti-ENPP1 antibody showing selective membrane staining of ENPP1 on metastatic cancer cells. Scale bar, 1 mm (C) and 50 μm (D). E, Representative images of human TNBCs stained using anti-ENPP1 antibody; scale bar, 100 μm. F, Distant metastasis-free survival (DMFS) in patients with TNBC stratified based on their Enpp1 and Cgas expression; n = 159, significance tested using log-rank test.

Close modal

To investigate the impact of ENPP1 protein expression on metastasis, we analyzed a total of 429 primary breast tumors from three independent cohorts for which there were long-term clinical follow-up data available. Similar to our findings in mucosal melanoma, we observed three distinct patterns of ENPP1 protein expression: tumor cell–dominant, stroma-dominant, and negative (Fig. 5E). Overall, 64% of primary TNBCs exhibited moderate or strong ENPP1 staining in either tumor cells or the stroma—a distribution that was consistent across the two ER cohorts (Supplementary Fig. S8B). On the other hand, 90% of ER+ tumors exhibited elevated ENPP1 protein expression. Notably, the tissue distribution and expression patterns varied between the two breast cancer subtypes, with ER tumors displaying both stromal and tumor cell–specific expression compared with their ER+ counterparts, which had a proclivity for tumor cell–specific staining (Supplementary Fig. S8B). Irrespective of the expression patterns, however, moderate-to-strong ENPP1 staining in the tumor was associated with poor prognosis, as evidenced by reduced overall survival, distant metastasis-free survival, and recurrence-free survival (Supplementary Fig. S8C–S8E). We next reasoned that if the association between ENPP1 expression and prognosis was related to its function as a negative regulator of cGAS-STING signaling, then its expression levels should be discriminatory only in tumors with high cGAS expression and activity in micronuclei. Staining using anti-cGAS antibodies revealed predominant staining at micronuclei in human tumors (example shown in Supplementary Fig. S9A and S9B). Indeed, ENPP1 protein expression was associated with reduced distant metastasis-free survival only in tumors with a preponderance of cGAS-positive micronuclei, and it had no significant association with metastasis in tumors with sparse cGAS-positive micronuclei (Fig. 5F). Collectively, these data are in agreement with our in vivo experimental results and further support the role of ENPP1 as an important determinant of cancer progression through its suppression of CIN-induced inflammatory signaling.

ENPP1 Is Associated with Immune Suppression in Human Cancer

We next correlated ENPP1 protein levels with tumor-infiltrating lymphocytes (TIL) and CD8+ T-cell density across breast cancers and found an inverse correlation between ENPP1 IHC expression intensity and lymphocytic infiltration (Fig. 6A and B; Supplementary Fig. S9C and S9D). Similar patterns were seen across the TCGA breast tumor cohort. We segregated 1,079 breast tumors into four subsets based on their relative Cgas and Enpp1 expression levels and used the CIBERSORT method to infer the prevalence of immune cell subsets from tissue expression profiles (34). Expectedly, Enpp1 expression was minimally associated with the immune cell fraction in tumors with low Cgas expression, whereas in those with high Cgas mRNA, it was inversely correlated with the overall leukocyte fraction as well as with the proportion of CD8+ T cells, CD4+ T cells, and proinflammatory macrophages (Fig. 6C). Furthermore, PD-L1 expression was highest in tumors with high Cgas and low Enpp1 expression. Gene set enrichment analysis (GSEA) comparing CgashiEnpp1hi with CgashiEnpp1lo breast tumors revealed upregulation of inflammatory pathways related to allograft rejection, type I IFN, and IFNγ-associated responses in the latter subset of tumors (Supplementary Fig. S9E). These findings suggest that Enpp1-to-Cgas ratio might be more predictive of tumor immune infiltration compared with Enpp1 expression levels alone. We orthogonally validated this assumption in sarcomas and mucosal melanoma tumors. In sarcomas, ENPP1-to-CGAS expression ratio was more strongly associated with the cytotoxic lymphocyte score compared with ENPP1 expression levels alone (Supplementary Fig. S9F). In mucosal melanomas, tumors with numerous cGAS-positive micronuclei and low ENPP1 expression exhibited increased CD8+ T-cell density, whereas those with elevated ENPP1 expression in the setting of widespread cGAS-positive micronuclei exhibited significantly reduced CD8+ T-cell infiltration (Supplementary Fig. S10A and S10B).

Figure 6.

Enpp1 expression is associated with reduced lymphocytic infiltration in human cancer. A, Percentage of TILs in breast tumors stratified based on their ENPP1 expression; bars, mean ± SEM; ***, P < 0.001, two-tailed t test. B, Representative images of human breast cancers stained using anti-ENPP1 or anti-CD8 antibodies. Scale bar, 100 μm. C, Tumor immune infiltration inferred using the CIBERSORT method on breast tumors found in the TCGA; box plots represent median, lower, and upper quartiles; error bars represent 10th and 90th percentiles, n = 1,079 tumors; ****, P < 0.0001, two-sided Mann–Whitney test. D, Percent objective response rate to anti–PD-1/PD-L1 therapy as a function of Enpp1 expression by cancer type for tumor histologies with high levels of Cgas expression. E, Schematic illustrating the consequence of ENPP1 activity (right) or its absence (left) on cancer metastasis and immune evasion.

Figure 6.

Enpp1 expression is associated with reduced lymphocytic infiltration in human cancer. A, Percentage of TILs in breast tumors stratified based on their ENPP1 expression; bars, mean ± SEM; ***, P < 0.001, two-tailed t test. B, Representative images of human breast cancers stained using anti-ENPP1 or anti-CD8 antibodies. Scale bar, 100 μm. C, Tumor immune infiltration inferred using the CIBERSORT method on breast tumors found in the TCGA; box plots represent median, lower, and upper quartiles; error bars represent 10th and 90th percentiles, n = 1,079 tumors; ****, P < 0.0001, two-sided Mann–Whitney test. D, Percent objective response rate to anti–PD-1/PD-L1 therapy as a function of Enpp1 expression by cancer type for tumor histologies with high levels of Cgas expression. E, Schematic illustrating the consequence of ENPP1 activity (right) or its absence (left) on cancer metastasis and immune evasion.

Close modal

In line with its role modulating tumor immune responses, we found that Enpp1 expression within a given cancer type negatively correlates with its overall response rate to aPD-1/anti–PD-L1 (aPD-L1) therapy (35). This inverse association was again restricted to tumor types characterized by elevated overall levels of Cgas expression (Fig. 6D; Supplementary Fig. S10C). We next analyzed the mRNA expression levels of CGAS and Enpp1 in 228 bladder cancers treated with aPD-L1 therapy and a smaller cohort of 52 TNBC tumors treated with aPD-1 (36, 37). Based on our TCGA analysis, these two cancer types exhibit relatively distinct Enpp1 expression levels, representing opposite end of the spectrum. Nonetheless, there was an overall positive correlation between CGAS and Enpp1 expression in bladder tumors where Enpp1 levels were significantly lower in the CGAShi subset of tumors that responded to aPD-L1 therapy. Expectedly, a low Enpp1-to-Cgas expression ratio was significantly associated with tumor response across both the bladder cancer and TNBC cohorts (Supplementary Fig. S10D and S10E).

Our work reveals an adaptive mechanism by which chromosomally unstable tumors co-opt cancer cell–intrinsic cGAS-STING signaling without eliciting antitumor immune surveillance (Fig. 6E). By virtue of their constant exposure to cytosolic dsDNA in micronuclei, cancer cells with CIN must address the consequences of cGAMP leakage into the extracellular space and its potential uptake by cells in the tumor microenvironment. By acquiring the ability to degrade cGAMP selectively in the extracellular environment, tumor cells can maintain relatively high levels of this metabolite in the intracellular compartment where it promotes metastatic progression (1), while minimizing antitumor paracrine STING activation in neighboring immune cells (Fig. 6E).

Previous work has linked ENPP1 to the ability of tumor cells to disseminate especially in the context of bone metastasis (38), yet the precise mechanisms underlying this relationship had remained poorly understood. One possible mechanism by which ENPP1 would facilitate tumor spread to the bone is through its contribution to pyrophosphate metabolism, promoting bone remodeling (39). Our data, however, indicate that the role of ENPP1 in tumor progression extends beyond osseous metastases, owing to its ability to hydrolyze cGAMP, and therefore suppresses the host's ability to control metastatic progression through activation of protective STING signaling in the tumor microenvironment.

Extracellular cGAMP hydrolysis by ENPP1 generates AMP, a substrate for adenosine production, thereby transforming an immune-stimulatory pathway into an immune-suppressive mechanism that promotes tumor progression (Fig. 6E). Our findings suggest that cGAMP represents a significant source of extracellular adenosine. Furthermore, the stepwise increase in ENPP1 levels—and concomitant decrease of CD39—during the evolution from primary tumors to metastasis suggests dynamic changes in the extracellular sources of adenosine, with ATP representing a significant source in primary tumors and the fractional contribution of cGAMP as an adenosine source increasing during tumor progression. Targeting extracellular adenosine production and signaling is currently being investigated at the preclinical and clinical stages (25). ENPP1 inhibition would achieve the dual purpose of reducing the extracellular levels of an immune suppressor while simultaneously increasing extracellular levels of the immunostimulatory metabolite cGAMP. These findings highlight an important STING-independent function for tumor cGAS and suggest that, in the presence of ENPP1, high tumor cGAS activity might in fact be paradoxically immune-suppressive, enabling tolerance for CIN and pervasive cytosolic dsDNA in advanced cancers.

Through extensive assessment of ENPP1 mRNA and protein expression levels across human tumors, our work positions ENPP1 into the broader clinical context and makes the case for the development of ENPP1 inhibitors for the treatment of advanced and chromosomally unstable cancers (19, 40, 41). Interestingly, cancer types with elevated ENPP1 expression are generally thought to be less responsive to ICB therapy, raising the possibility that extracellular purine metabolism might represent an important innate immune checkpoint that must be overcome for the full activation of the adaptive immune response against cancer. Indeed, our work suggests that ENPP1 inhibition is a viable mechanism to sensitize otherwise resistant tumors to ICB therapy. Interestingly, the widespread stromal staining patterns of ENPP1 in human cancers—reminiscent of fibroblast expression—suggest that this mechanism of immune evasion might arise not only from tumor cells but also from cells in the tumor microenvironment. Given its low expression levels in normal tissues, it will be important to dissect tumor-derived factors that promote induction of ENPP1 in the stroma. Nonetheless, our data suggest that in metastatic cancers, ENPP1 staining is biased toward a cancer cell–intrinsic pattern, raising the possibility that tumor cells that acquire the ability to transmute cGAMP-mediated immune activation into immune suppression have a selective advantage to spread to distant organs.

Therapies that activate STING (also known as STING agonists) have been the focus of intense investigation given their ability to elicit antitumor immunity through type I IFN signaling (42). Inhibition of ENPP1 is distinct from direct pharmacologic activation of STING in a number of important ways. First, ENPP1 tilts the relative balance of STING activation away from cancer cells, where it promotes metastatic progression (1), and toward host cells where it potentiates antitumor immunity. STING agonists indiscriminately activate STING in both cancer cells and the host, promoting dichotomous outcomes. Second, inhibition of cGAMP hydrolysis by ENPP1 would primarily affect cGAMP concentrations at the microscopic scales relevant to paracrine tumor cell–host cell interactions. This is particularly relevant given the short half-lives of extracellular cellular cGAMP and adenosine (21). Furthermore, this critical spatial consideration is likely to minimize any potential side effects that might be observed during the systemic administration of STING agonists, thus offering a larger therapeutic window. Third, ENPP1 is selectively upregulated in metastatic and chromosomally unstable tumor cells, and a systemic ENPP1 inhibitor would interfere with the ability of disseminated tumor cells to evade immune surveillance arising from CIN, bypassing the need for technically challenging intratumoral administration that is typical of STING agonists. In summary, our work highlights the therapeutic utility of selectively targeting cancer cell dependencies on CIN and the mechanism by which they have evolved to tolerate it.

Cell Culture

4T1 (ATCC; catalog no. CRL-2539), CT26 (ATCC; catalog no. CRL-2638), and B16F10 (ATCC; catalog no. CRL-6475) cell lines were purchased from the ATCC, and E0771 was a gift from Alexander Rudensky. Cells were cultured in DMEM (B16F10 and E0771) or RPMI (4T1 and CT26) supplemented with 10% FBS and 2 mM l-glutamine in the presence of penicillin (50 U/mL) and streptomycin (50 μg/mL). All cells were found to be negative for Mycoplasma upon repeated testing every 2 months using the MycoAlert Mycoplasma Detection Kit (Lonza; catalog no. LT07-318). Cells were used within three to five passages. Details of cell line generation using CRISPR/Cas9 KO and shRNA knockdowns are included in the Supplementary Methods section, and gRNA and shRNA sequences are listed in Supplementary Table S1.

Immunofluorescence and Immunoblotting

Detailed protocols for immunoblotting and immunofluorescence are described in the Supplementary Methods, and antibodies used in these protocols are listed in Supplementary Tables S2 and S3, respectively.

cGAMP Quantification

For intracellular and extracellular cGAMP quantification in cancer cell lines, cancer cells were seeded in 15 cm culture dishes. When culture plates were 80% to 90% confluent, media were changed to serum-free phenol red–free RPMI (Corning). Sixteen hours following media exchange, the conditioned media were removed and centrifuged at ≥ 600 × g at 4°C for 15 minutes. Supernatant was assayed directly. All the steps were performed on ice. Cells were washed with PBS twice and then trypsinized for 5 minutes at 37°C, and cell counts were measured. Cells were then centrifuged at ≥ 600 × g at 4°C for 15 minutes. Whole cell lysates were generated by lysing the cell pellet in LP2 lysis buffer [Tris HCl pH 7.7, 20 mM, NaCl 100 mM, NaF 10 mM, beta-glycerophosphate 20 mmol/L, MgCl2 5 mM, Triton X-100 0.1% (v/v), Glycerol 5% (v/v)]. The homogenate was then subjected to centrifugation at 10,000 × g for 15 minutes. cGAMP ELISA was performed according to the manufacturer's protocol using DetectX Direct 2′,3′-Cyclic GAMP Enzyme Immunoassay Kit (Arbo Assay).

Hematoxylin and Eosin (H&E) Staining and Immune Phenotyping of Lung Metastases

All antibodies used in IHC are listed in Supplementary Table S4. Lungs were excised from euthanized mice and submerged in 4% paraformaldehyde (PFA) overnight at 4°C and then were transferred to 70% ethanol. Tissue embedding, slide sectioning, and H&E staining were performed by the Molecular Cytology Core Facility at Memorial Sloan Kettering Cancer Center (MSKCC). IHC for CD8 and CD45 staining was performed using anti-CD8 (Cell Signaling Technology no. 98941) and anti-CD45 (Biosciences 550539) by the Laboratory of Comparative Pathology at MSKCC. For immune profiling using flow cytometry, animals were sacrificed 18 days after tail-vein injection with control and Enpp1-KO 4T1 cells. Lungs were perfused through the right ventricle with 10 to 15 mL of PBS. The lungs were removed, and the large airways, thymus, and lymph nodes were dissected from the peripheral lung tissue. The peripheral lung tissue was minced and transferred into 50 mL falcon tubes and processed in digestion buffer by mouse tumor dissociation kit (Miltenyi), according to the manufacturer's instructions. Homogenized lungs were passed through 40-μm nylon mesh to obtain a single-cell suspension. The remaining red blood cells were lysed using BD Pharm Lyse (BD Biosciences). Cells were stained with viability dye LIVE/DEAD Fixable Blue Dead Cell Stain Kit (Invitrogen), followed by incubation with FcBlock (Invitrogen), and stained with a mixture of fluorochrome-conjugated antibodies (see Supplementary Table S5 for a list of antibodies, clones, fluorochromes, and manufacturers). Data were acquired on a BD LSR II flow cytometer using BD FACS Diva software (BD Biosciences); compensation and data analysis were performed using FCS express 7 software. Unstained biological controls and single-color controls were used. Cell populations were identified using sequential gating strategy (Supplementary Fig. S4C).

Adenosine Measurements

4T1 cells were seeded in 10 cm culture dishes in quadruplicates. When culture plates reached 80% to 90% confluence, 7 mL serum-free phenol red–free RPMI (Corning) with and without inhibitors (EHNA 100 μM, NBMPR 100 μM, and dipyridamole 40 μM) was added to plates. Conditioned media were collected after 16-hour incubation. Conditioned media were centrifuged at 10,000 × g for 10 minutes at 4°C. Cells were harvested and cell counts were recorded for back calculations. Direct quantification of adenosine in flash-frozen conditioned media was performed by Charles River Laboratories Inc. Adenosine concentrations were determined by high-performance liquid chromatography with MS-MS detection in multiple reaction monitoring (MRM) mode. In brief, 4 μL of internal standard solution containing 10 nM Adenosine-13C5 was added to 10 μL of undiluted experimental sample. Ten microliters was injected into an Infinity 1290 LC system (Agilent) by an automated sample injector (SIL-20AD, Shimadzu). Analytes were separated by LC using a linear gradient of mobile phase B at a flow rate of 0.200 mL/min on a reversed-phase Atlantis T3 C18 column (2.1 × 150 mm, 3.0 μm particle size; Waters) held at a temperature of 40°C. Mobile phase A consisted of 5 mM ammonium formate in ultrapure water. Mobile phase B was methanol. Acquisitions were achieved in the positive ionization mode using a QTrap 5500 (Applied Biosystems) equipped with a Turbo Ion Spray interface. The ion spray voltage was set at 5.0 kV, and the probe temperature was 500°C. The collision gas (nitrogen) pressure was kept at the medium setting level. The following MRM transitions were used for quantification: m/z 268.2/136.1 for adenosine. Data were calibrated and quantified using the Analyst data system (Applied Biosystems, version 1.5.2). Indirect adenosine measurements in conditioned media after cGAMP addition were performed using the adenosine assay Kit (Cell Biolabs) according to a modified manufacturer's protocol; for each sample, we measured fluorescence intensity at 600 nm with and without the ADA inhibitor EHNA (Supplementary Fig. S3A and S3B).

Animal Metastasis Studies

Animal experiments were performed in accordance with protocols approved by the MSKCC Institutional Animal Care and Use Committee (IACUC). For survival experiments in 4T1 experiments, power analysis indicated that 15 mice per group would be sufficient to detect a difference at relative HRs of <0.25 or >4.0 with 80% power and 95% confidence, given a median survival of 58 days in the control group and a total follow-up period of 180 days also accounting for accidental animal death during procedures. There was no need to randomize animals. Investigators were not blinded to group allocation. For tail-vein injections, 2.5 × 104 4T1 or 5 × 104 CT26 cells were injected into the tail vein of 6- to 7-week-old BALB/c mice. Metastasis was primarily assessed through overall survival. Overall survival endpoint was met when the mice died or met the criteria for euthanasia under the IACUC protocol. Surface lung metastases were assessed at endpoint by direct visual examination after euthanasia at which point lungs were perfused and fixed in 4% PFA (4T1 experiments) or stained using india ink (CT26 experiments). Furthermore, lung metastasis after injection of 4T1 cells was qualitatively assessed using routine H&E staining as shown in Fig. 5E. Metastatic dissemination in Supplementary Fig. S2J was determined using BLI. Mice were injected with d-luciferin (150 mg/kg) and subjected to BLI using tan IVIS Spectrum Xenogen instrument (Caliper Life Sciences) to image locoregional recurrence as well as distant metastases. BLI images were analyzed using Living Image Software v.2.50. For orthotopic tumor implantation, 1.25 × 105 4T1 cells in 50 μL PBS were mixed 1:1 with Matrigel (BD Biosciences) and injected into the fourth mammary fat pad. Only one tumor was implanted per animal. Primary tumors were surgically excised on day 7 after implantation, and metastatic dissemination was assessed by monitoring overall survival or on day 30 through quantification of surface lung metastases upon euthanasia. In the E0771 metastasis model, 2.5 × 105 tdTomato-Luciferase–expressing E0771 cells were injected into the tail vein of 7- to 12-week old C56BL/6 or MPYS−/− (Tmem173−/−, The Jackson Laboratory stock number 025805) mice. Metastatic dissemination was accessed by BLI.

RNA-Sequencing Analysis of TCGA Tumors

RNA-sequencing (RNA-seq) data for human tumor samples from TCGA patients were obtained from https://gdc.cancer.gov/about-data/publications/pancanatlas (43). The data are upper-quartile–normalized RNA-seq by expectation minimization for batch-corrected mRNA gene expression and are from 33 different cancer types. Overall leukocyte fractions and CIBERSORT immune fractions for the TCGA Breast Cancer (BRCA) patients were obtained from https://gdc.cancer.gov/node/998 (44). The absolute abundance of the CIBERSORT immune cell types was obtained by multiplying the leukocyte fraction by the CIBERSORT immune fractions. The expression values for ENNP1 and CGAS from the TCGA RNA-seq data were utilized to categorize tumors into the four groups: Enpp1loCgaslo, Enpp1hiCgaslo, Enpp1loCgashi, and Enpp1hiCgashi. The median expression value per cancer type was used to categorize tumors into Enpp1lo and Enpp1hi groups. Tumors with expression values less than or equal to the median for a given cancer type were considered Enpp1lo, whereas tumors with expression values above the median were considered Enpp1hi. The bottom tertile expression value per cancer type was used to categorize tumors into Cgaslo and Cgashi groups. Tumors with expression values less than or equal to the bottom tertile (<33%) of CGAS expression in a given cancer type were categorized as CGASlo, whereas tumors with expression values greater than the bottom tertile (>33%) were categorized as Cgashi. The Wilcoxon rank-sum test was used to compare the relative abundance of CIBERSORT immune cell types between different Cgas/Enpp1 expression subgroups. For pathway enrichment analysis, the DESeq2 R package (45) was used to identify differentially expressed genes between the Enpp1loCgashi and Enpp1hiCgashi groups within the TCGA BRCA cohort. The GSEA method (46) was used to perform a pathway enrichment analysis between the Enpp1loCgashi and Enpp1hiCgashi groups. A preranked gene list from DESeq2 was created and sorted by the following: sign of the log fold change × -log(adjusted P value). The sorted preranked list was run in GSEA with the Hallmark gene set database that was downloaded from the Molecular Signatures Database (46). Survival analysis across TCGA tumor types was performed using KMPlot (http://www.kmplot.com) using autoselection for best cutoff between the 25th and 75th percentiles.

Animal Immunotherapy Experiments

To assess the role of ENPP1 in the primary tumor growth upon the ICB, we adopted the 4T1 orthotopic mammary fat pad implantation model. First, 4T1 (4T1-Luc) cells and 4T1-Luc Enpp1-KO cells were generated by stably integrating the Lentivirus pLVX vector expressing the tdTomato-Luciferase fusion gene in the 4T1 and 4T1 Enpp1-KO cells, respectively. Fifteen 7-week-old mice were used for each of the arms, including four combinations of two cell lines (4T1-Luc and 4T1-Luc ENPP1 KO) and two conditions (ICB and the isotype control treatment). 4T1-Luc cells or 4T1-Luc Enpp1-KO cells (1.25 × 105) in PBS:Matrigel (1:1) mix were injected into the mammary fat pad of BALB/c mice. Two hundred microgram rat anti-mouse PD-1 IgG2a antibody (aPD-1) and 100 μg mouse anti-mouse CTLA4 IgG2b antibody (aCTLA4) or their corresponding isotype control antibodies were delivered intraperitoneally in 100 mL of PBS to mice every 3 days starting at day 6 after implantation. After four doses of combined ICB, maintenance aCTLA4 treatment and the corresponding isotype control were given every 3 days. The length (L) and width (W) of the tumor were measured using calipers. The tumor size was calculated according to the following formula: L × W2/2. For experiment in Fig. 4CE, endpoint was determined when primary tumor reached the size of 2,000 mm3. For the CT26 model, 5 × 104 eGFP- or eGFP-ENPP1–expressing CT26 cells were delivered intravenously to 7-week-old BALB/c mice. Treatment with aPD-1/aCTLA4 antibodies and their corresponding isotype control antibodies was initiated intraperitoneally starting on day 6 and given every 3 days for five total doses. Animals were monitored for overall survival. For the E0771 model, 5 × 105 eGFP- or eGFP-ENPP1–expressing E0771-Luc cells in PBS:Matrigel (1:1) mix were injected into the mammary fat pad of C57BL/6 WT mice or MPYS−/− (Tmem173−/−, The Jackson Laboratory stock number 025805) at the age of 7 weeks. Treatment with 200 μg of aPD-1 or its corresponding isotype control antibody was given on days 6, 10, and 13.

Data Availability

Tumor DNA and RNA-seq data used in this article are publicly available and cited as appropriate in the text and Methods section. No new code was used in this article.

M. Salto-Tellez reports other support from Roche and Philips, and personal fees from Incyte, MSD/Targos, and Sonrai Analytics outside the submitted work. M.G. Hanna reports other support from Paige outside the submitted work. K. Litchfield reports personal fees from Roche Tissue Diagnostics, Monopteros, and CIC, and grants from CRUK TDL/Ono/LifeArc alliance outside the submitted work; in addition, K. Litchfield has a patent on indel burden and checkpoint inhibitor response pending to CRUK. D. Biswas reports personal fees from NanoString and grants from National Institute for Health Research Biomedical Research Centre, Innovation (i2i) Crick translation scheme supported by the Medical Research Council, Jean Shanks Foundation, and University College London MBPhD programme outside the submitted work; in addition, D. Biswas has a patent for method of predicting survival rates for cancer patients issued. T. Merghoub is a cofounder and holds an equity in IMVAQ Therapeutics; is a consultant of Immunos Therapeutics and Pfizer; has research support from Bristol-Myers Squibb, Surface Oncology, Kyn Therapeutics, Infinity Pharmaceuticals, Inc., Peregrine Pharmaceuticals, Inc., Adaptive Biotechnologies, Leap Therapeutics, Inc., and Aprea; and has patents on applications related to work on oncolytic viral therapy, alpha virus-based vaccine, neoantigen modeling, CD40, GITR, OX40, PD-1, and CTLA4. J.S. Reis-Filho reports personal fees from Paige, Repare Therapeutics, Goldman Sachs, Grupo Oncoclinicas, Roche Tissue Diagnostics, Novartis, Genentech, Roche, and inviCRO outside the submitted work. N. Riaz reports grants from Pfizer and BMS, and grants and personal fees from REPARE Therapeutics outside the submitted work. S.-S.M. Su reports personal fees from Volastra Therapeutics during the conduct of the study and personal fees from Volastra Therapeutics outside the submitted work; in addition, S.-S.M. Su has a patent pending application. N. Vasan reports personal fees from Volastra Therapeutics, Novartis, and Petra Pharmaceuticals, and other from Heligenics, Inc. outside the submitted work. S.N. Powell reports personal fees from AstraZeneca, Varian, and Philips outside the submitted work. J.D. Wolchok reports grants and personal fees from Bristol Myers Squibb during the conduct of the study; and personal fees and other support from Tizona Pharmaceuticals, Imvaq Therapeutics, Beigene, Apricity Therapeutics, Arsenal IO, Georgiamune, Serametrix, and Adaptive Biotech; personal fees from Linneaus, Amgen, Ascentage, Astellas, AstraZeneca, Bayer, Boehringer Ingelheim, Celgene, Chugai, Eli Lilly, Elucida, F Star, Kyowa Hakko Kirin, Merck, Neon Therapeutics, Polynoma, Psioxus, Recepta, Takara Bio, Trieza, Truvax, Sellas, Surface Oncology, Syndax, Syntalogic, and Werewolf Therapeutics; and grants from Sephora outside the submitted work; in addition, J.D. Wolchok has a patent for Alphavirus replicon particles expressing issued, a patent for Newcastle Disease viruses for Cancer Therapy issued, a patent for anti–PD-1 antibody issued and licensed to Agenus, a patent for anti-CTLA4 antibodies issued and licensed to Agenus, a patent for anti-GITR antibodies and methods of use thereof issued and licensed to Agenus/Incyte, a patent for anti-CD40 agonist mAb fused to Monophosphoryl Lipid A (MPL) for cancer therapy pending, a patent for CAR+ T cells targeting differentiation antigens as means to treat cancer pending, a patent for Identifying and Treating Subjects at Risk for Checkpoint Blockade Therapy Associated Colitis pending, a patent for immunosuppressive follicular helper-like T cells modulated by immune checkpoint blockade pending, a patent for phosphatidylserine targeting agents and uses thereof for adoptive T-cell therapies pending, a patent for Xenogeneic DNA Vaccines licensed and with royalties paid from Merial, and a patent for myeloid-derived suppressor cell assay licensed and with royalties paid from Serametrix. O. Elemento reports grants and nonfinancial support from Volastra Therapeutics during the conduct of the study; and other support from One Three Biotech outside the submitted work; O. Elemento is on the scientific advisory board of Freenome, OWKIN, and Champions Oncology. C. Swanton reports grants and personal fees from AstraZeneca, Pfizer, Ono Pharmaceutical, Roche Ventana, and Bristol Myers Squibb; grants from Archer Dx Inc.(collaboration in minimal residual disease sequencing technologies) and Boehringer Ingelheim; personal fees from Novartis, MSD, Illumina, Celgene, GSK, Sarah Canon Research Institute, Medixci, GRAIL, Genentech, Bicycle Therapeutics, and Epic Biosciences; personal fees and other support from Achilles Therapeutics and EPIC Biosciences; other support from ApoGen Biotech; stock option in GRAIL; and being Chief Investigator for the MeRmaiD1 clinical trial during the conduct of the study. C. Swanton also reports grants, personal fees, and other from AstraZeneca; grants from Archer Dx Inc. and Boehringer Ingelheim; grants and personal fees from Pfizer, Ono Pharmaceutical, Roche Ventana, and Bristol Myers Squibb; personal fees from Novartis, MSD, Illumina, Celgene, GSK, Sarah Canon Research Institute, Medixci, Genentech, and Bicycle Therapeutics; personal fees and other support from GRAIL; stock options in Epic Biosciences and Apogen Biotech; and stock options and being cofounder of Achilles Therapeutics outside the submitted work; in addition, C. Swanton has a patent for Immune checkpoint intervention in cancer (PCT/EP2016/071471) issued, a patent for method for treating cancer based on identification of clonal neo-antigens (PCT/EP2016/059401) issued, a patent for methods for lung cancer detection (PCT/US2017/028013) issued, a patent for method of detecting tumor recurrence (PCT/GB2017/053289) issued, a patent for method of treating cancer by targeting insertion/deletion mutations (PCT/GB2018/051893) issued, a patent for method of identifying insertion/deletion mutation targets (PCT/GB2018/051892) issued, a patent for method for determining whether an HLA allele is lost in a tumor (PCT/GB2018/052004) issued, a patent for method for identifying responders to cancer treatment (PCT/GB2018/051912) issued, and a patent for method of predicting survival rates for cancer patients (PCT/GB2020/050221) issued. C. Swanton is an editorial board member for Cell and PLOS Medicine, associate editor and editorial board member for Annals of Oncology, scientific editor for Cancer Discovery, advisory board member for Nature Reviews Clinical Oncology and Cancer Cell, and Royal Society Napier Research Professor (RP150154). This work was supported by the Francis Crick Institute that receives its core funding from Cancer Research UK (FC001169), the UK Medical Research Council (FC001169), and the Wellcome Trust (FC001169). C. Swanton is funded by Cancer Research UK (TRACERx, PEACE, and CRUK Cancer Immunotherapy Catalyst Network), Cancer Research UK Lung Cancer Centre of Excellence, the Rosetrees Trust, Butterfield and Stoneygate Trusts, NovoNordisk Foundation (ID16584), Royal Society Professorship Enhancement Award (RP/EA/180007), the NIHR BRC at University College London Hospitals, the CRUK-UCL Centre, Experimental Cancer Medicine Centre, and the Breast Cancer Research Foundation. This research is supported by a Stand Up To Cancer–LUNGevity–American Lung Association Lung Cancer Interception Dream Team Translational Research Grant (grant number: SU2C-AACR-DT23-17). Stand Up To Cancer (SU2C) is a program of the Entertainment Industry Foundation. Research grants are administered by the American Association for Cancer Research, the Scientific Partner of SU2C. C. Swanton receives funding from the European Research Council (ERC) under the European Union's Seventh Framework Programme (FP7/2007-2013) Consolidator Grant (FP7-THESEUS-617844), European Commission ITN (FP7-PloidyNet 607722), an ERC Advanced Grant (PROTEUS) from the European Research Council under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 835297), and Chromavision from the European Union's Horizon 2020 research and innovation programme (grant agreement No. 665233). A.N. Shoushtari reports grants and personal fees from Bristol-Myers Squibb and Immunocore; grants from Xcovery, Novartis, and Pfizer; and personal fees from Castle Biosciences outside the submitted work. E.E. Parkes reports personal fees from Boehringer Ingelheim outside the submitted work. B. Izar reports personal fees from Merck and Volastra Therapeutics outside the submitted work. S.F. Bakhoum reports personal fees from Volastra Therapeutics Inc. during the conduct of the study; and personal fees from Sanofi and other support from Cancer Research UK and Prostate Cancer Foundation outside the submitted work; in addition, S.F. Bakhoum has a patent for targeting chromosomal instability and downstream cytosolic DNA signaling for cancer treatment licensed to Volastra Therapeutics Inc. and a patent for methods and strategies to target cytosolic dsDNA signaling in chromosomally unstable cancers pending. No disclosures were reported by the other authors.

J. Li: Conceptualization, data curation, formal analysis, investigation. M.A. Duran: Data curation, formal analysis, investigation. N. Dhanota: Data curation, formal analysis, investigation. W.K. Chatila: Data curation, formal analysis, investigation. S.E. Bettigole: Conceptualization, investigation, methodology. J. Kwon: Data curation, formal analysis, investigation. R.K. Sriram: Investigation. M.P. Humphries: Data curation. M. Salto-Tellez: Data curation. J.A. James: Data curation. M.G. Hanna: Formal analysis. J.C. Melms: Investigation. S. Vallabhaneni: Investigation. K. Litchfield: Data curation, formal analysis. I. Usaite: Formal analysis. D. Biswas: Formal analysis. R. Bareja: Formal analysis. H.W. Li: Investigation. M.L. Martin: Data curation. P. Dorsaint: Formal analysis. J.-A. Cavallo: Investigation. P. Li: Methodology. C. Pauli: Formal analysis. L. Gottesdiener: Formal analysis, investigation. B.J. DiPardo: Formal analysis. T.J. Hollmann: Data curation, investigation. T. Merghoub: Methodology. H.Y. Wen: Data curation, methodology. J.S. Reis-Filho: Methodology. N. Riaz: Formal analysis. S.-S.M. Su: Conceptualization. A. Kalbasi: Formal analysis. N. Vasan: Formal analysis. S.N. Powell: Conceptualization, funding acquisition. J.D. Wolchok: Conceptualization, methodology. O. Elemento: Data curation. C. Swanton: Conceptualization, formal analysis. A.N. Shoushtari: Data curation, formal analysis. E.E. Parkes: Conceptualization, supervision, investigation. B. Izar: Conceptualization, supervision, investigation. S.F. Bakhoum: Conceptualization, resources, supervision, funding acquisition, investigation, writing–original draft.

The authors would like to thank Stephen McQuaid and Christine Greene (Northern Ireland Biobank) for technical support, and members of the Bakhoum Laboratory (MSKCC), John Maciejowski (MSKCC), Ashley Laughney (Weill Cornell Medicine), Mathieu Bakhoum (University of California, San Diego), and Christopher Garris (Rockefeller University) for constructive feedback. They would also like to thank the MSKCC genomics, molecule cytology, and RNAi cores. S.F. Bakhoum is supported by the Office of the Director, National Institutes of Health, under Award Number DP5OD026395 High-Risk High-Reward Program, the NCI Breast Cancer SPORE (P50CA247749), the Burroughs Wellcome Fund Career Award for Medical Scientists, the Parker Institute for Immunotherapy at MSKCC, the Josie Robertson Foundation, and the MSKCC core grant P30-CA008748. J.S. Reis-Filho is funded in part by the Breast Cancer Research Foundation. B. Izar is supported by NIH grants K08CA222663 and U54CA225088, the Burroughs Wellcome Fund Career Award for Medical Scientists, the Louis V. Gerstner, Jr. Scholars Program, and Herbert Irving Comprehensive Cancer Center core grant P30CA013696. E.E. Parkes is supported by the Oxford Institute for Radiation Oncology, the Prostate Cancer Foundation, the American Society of Clinical Oncology, and the Academy of Medical Sciences. The Northern Ireland Biobank has received funds from HSC Research Development Division of the Public Health Agency in Northern Ireland and the Friends of the Cancer Centre. J.D. Wolchok is supported by Swim Across America, Ludwig Cancer Research. S.F. Bakhoum, A.N. Shoushtari, and J.D. Wolchok are supported by the Nonna's Garden Foundation.

1.
Bakhoum
SF
,
Ngo
B
,
Laughney
AM
,
Cavallo
JA
,
Murphy
CJ
,
Ly
P
, et al
Chromosomal instability drives metastasis through a cytosolic DNA response
.
Nature
2018
;
553
:
467
72
.
2.
Burrell
RA
,
McGranahan
N
,
Bartek
J
,
Swanton
C
. 
The causes and consequences of genetic heterogeneity in cancer evolution
.
Nature
2013
;
501
:
338
45
.
3.
Rosenthal
R
,
Cadieux
EL
,
Salgado
R
,
Bakir
MA
,
Moore
DA
,
Hiley
CT
, et al
Neoantigen-directed immune escape in lung cancer evolution
.
Nature
2019
;
567
:
479
85
.
4.
Watkins
TBK
,
Lim
EL
,
Petkovic
M
,
Elizalde
S
,
Birkbak
NJ
,
Wilson
GA
, et al
Pervasive chromosomal instability and karyotype order in tumour evolution
.
Nature
2020
;
587
:
126
32
.
5.
McGranahan
N
,
Rosenthal
R
,
Hiley
CT
,
Rowan
AJ
,
Watkins
TBK
,
Wilson
GA
, et al
Allele-specific HLA loss and immune escape in lung cancer evolution
.
Cell
2017
;
171
:
1259
71
.
6.
Bakhoum
SF
,
Cantley
LC
. 
The multifaceted role of chromosomal instability in cancer and its microenvironment
.
Cell
2018
;
174
:
1347
60
.
7.
Crasta
K
,
Ganem
NJ
,
Dagher
R
,
Lantermann
AB
,
Ivanova
EV
,
Pan
Y
, et al
DNA breaks and chromosome pulverization from errors in mitosis
.
Nature
2012
;
482
:
53
8
.
8.
Kato
H
,
Sandberg
AA
. 
Chromosome pulverization in human cells with micronuclei
.
J Natl Cancer Inst
1968
;
40
:
165
79
.
9.
Hatch
EM
,
Fischer
AH
,
Deerinck
TJ
,
Hetzer
MW
. 
Catastrophic nuclear envelope collapse in cancer cell micronuclei
.
Cell
2013
;
154
:
47
60
.
10.
Yang
H
,
Wang
H
,
Ren
J
,
Chen
Q
,
Chen
ZJ
. 
cGAS is essential for cellular senescence
.
Proc Natl Acad Sci U S A
2017
;
114
:
E4612
E20
.
11.
Mackenzie
KJ
,
Carroll
P
,
Martin
CA
,
Murina
O
,
Fluteau
A
,
Simpson
DJ
, et al
cGAS surveillance of micronuclei links genome instability to innate immunity
.
Nature
2017
;
548
:
461
5
.
12.
Harding
SM
,
Benci
JL
,
Irianto
J
,
Discher
DE
,
Minn
AJ
,
Greenberg
RA
. 
Mitotic progression following DNA damage enables pattern recognition within micronuclei
.
Nature
2017
;
548
:
466
70
.
13.
Sun
L
,
Wu
J
,
Du
F
,
Chen
X
,
Chen
ZJ
. 
Cyclic GMP-AMP synthase is a cytosolic DNA sensor that activates the type I interferon pathway
.
Science
2013
;
339
:
786
91
.
14.
Ishikawa
H
,
Barber
GN
. 
STING is an endoplasmic reticulum adaptor that facilitates innate immune signalling
.
Nature
2008
;
455
:
674
8
.
15.
Abe
T
,
Barber
GN
. 
Cytosolic-DNA-mediated, STING-dependent proinflammatory gene induction necessitates canonical NF-kappaB activation through TBK1
.
J Virol
2014
;
88
:
5328
41
.
16.
Dou
Z
,
Ghosh
K
,
Vizioli
MG
,
Zhu
J
,
Sen
P
,
Wangensteen
KJ
, et al
Cytoplasmic chromatin triggers inflammation in senescence and cancer
.
Nature
2017
;
550
:
402
6
.
17.
Marcus
A
,
Mao
AJ
,
Lensink-Vasan
M
,
Wang
L
,
Vance
RE
,
Raulet
DH
. 
Tumor-derived cGAMP triggers a STING-mediated interferon response in non-tumor cells to activate the NK cell response
.
Immunity
2018
;
49
:
754
63
.
18.
Schadt
L
,
Sparano
C
,
Schweiger
NA
,
Silina
K
,
Cecconi
V
,
Lucchiari
G
, et al
Cancer-cell-intrinsic cGAS expression mediates tumor immunogenicity
.
Cell Rep
2019
;
29
:
1236
48
.
19.
Carozza
JA
,
Böhnert
V
,
Nguyen
KC
,
Skariah
G
,
Shaw
KE
,
Brown
JA
, et al
Extracellular cGAMP is a cancer-cell-produced immunotransmitter involved in radiation-induced anticancer immunity
.
Nat Cancer
2020
;
1
:
184
96
.
20.
Bakhoum
SF
,
Thompson
SL
,
Manning
AL
,
Compton
DA
. 
Genome stability is ensured by temporal control of kinetochore-microtubule dynamics
.
Nat Cell Biol
2009
;
11
:
27
35
.
21.
Li
L
,
Yin
Q
,
Kuss
P
,
Maliga
Z
,
Millan
JL
,
Wu
H
, et al
Hydrolysis of 2′3′-cGAMP by ENPP1 and design of nonhydrolyzable analogs
.
Nat Chem Biol
2014
;
10
:
1043
8
.
22.
Johnstone
CN
,
Smith
YE
,
Cao
Y
,
Burrows
AD
,
Cross
RS
,
Ling
X
, et al
Functional and molecular characterisation of EO771.LMB tumours, a new C57BL/6-mouse-derived model of spontaneously metastatic mammary cancer
.
Dis Model Mech
2015
;
8
:
237
51
.
23.
Chuang
CH
,
Greenside
PG
,
Rogers
ZN
,
Brady
JJ
,
Yang
D
,
Ma
RK
, et al
Molecular definition of a metastatic lung cancer state reveals a targetable CD109-Janus kinase-Stat axis
.
Nat Med
2017
;
23
:
291
300
.
24.
Kato
K
,
Nishimasu
H
,
Oikawa
D
,
Hirano
S
,
Hirano
H
,
Kasuya
G
, et al
Structural insights into cGAMP degradation by Ecto-nucleotide pyrophosphatase phosphodiesterase 1
.
Nat Commun
2018
;
9
:
4424
.
25.
Vijayan
D
,
Young
A
,
Teng
MWL
,
Smyth
MJ
. 
Targeting immunosuppressive adenosine in cancer
.
Nat Rev Cancer
2017
;
17
:
709
24
.
26.
Lofgren
L
,
Pehrsson
S
,
Hagglund
G
,
Tjellstrom
H
,
Nylander
S
. 
Accurate measurement of endogenous adenosine in human blood
.
PLoS One
2018
;
13
:
e0205707
.
27.
Stagg
J
,
Divisekera
U
,
McLaughlin
N
,
Sharkey
J
,
Pommey
S
,
Denoyer
D
, et al
Anti-CD73 antibody therapy inhibits breast tumor growth and metastasis
.
Proc Natl Acad Sci U S A
2010
;
107
:
1547
52
.
28.
Fridlender
ZG
,
Sun
J
,
Kim
S
,
Kapoor
V
,
Cheng
G
,
Ling
L
, et al
Polarization of tumor-associated neutrophil phenotype by TGF-beta: “N1” versus “N2” TAN
.
Cancer Cell
2009
;
16
:
183
94
.
29.
Fischer
MA
,
Davies
ML
,
Reider
IE
,
Heipertz
EL
,
Epler
MR
,
Sei
JJ
, et al
CD11b(+), Ly6G(+) cells produce type I interferon and exhibit tissue protective properties following peripheral virus infection
.
PLoS Pathog
2011
;
7
:
e1002374
.
30.
Liu
Y
,
O'Leary
CE
,
Wang
LS
,
Bhatti
TR
,
Dai
N
,
Kapoor
V
, et al
CD11b+Ly6G+ cells inhibit tumor growth by suppressing IL-17 production at early stages of tumorigenesis
.
Oncoimmunology
2016
;
5
:
e1061175
.
31.
Kim
K
,
Skora
AD
,
Li
Z
,
Liu
Q
,
Tam
AJ
,
Blosser
RL
, et al
Eradication of metastatic mouse cancers resistant to immune checkpoint blockade by suppression of myeloid-derived cells
.
Proc Natl Acad Sci U S A
2014
;
111
:
11774
9
.
32.
Dabas
N
,
Byrnes
DM
,
Rosa
AM
,
Eller
MS
,
Grichnik
JM
. 
Diagnostic role of chromosomal instability in melanoma
.
J Skin Cancer
2012
;
2012
:
914267
.
33.
D'Angelo
SP
,
Larkin
J
,
Sosman
JA
,
Lebbe
C
,
Brady
B
,
Neyns
B
, et al
Efficacy and safety of nivolumab alone or in combination with ipilimumab in patients with mucosal melanoma: a pooled analysis
.
J Clin Oncol
2017
;
35
:
226
35
.
34.
Newman
AM
,
Liu
CL
,
Green
MR
,
Gentles
AJ
,
Feng
W
,
Xu
Y
, et al
Robust enumeration of cell subsets from tissue expression profiles
.
Nat Methods
2015
;
12
:
453
7
.
35.
Yarchoan
M
,
Hopkins
A
,
Jaffee
EM
. 
Tumor mutational burden and response rate to PD-1 inhibition
.
N Engl J Med
2017
;
377
:
2500
1
.
36.
Mariathasan
S
,
Turley
SJ
,
Nickles
D
,
Castiglioni
A
,
Yuen
K
,
Wang
Y
, et al
TGFbeta attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells
.
Nature
2018
;
554
:
544
8
.
37.
Voorwerk
L
,
Slagter
M
,
Horlings
HM
,
Sikorska
K
,
van de Vijver
KK
,
de Maaker
M
, et al
Immune induction strategies in metastatic triple-negative breast cancer to enhance the sensitivity to PD-1 blockade: the TONIC trial
.
Nat Med
2019
;
25
:
920
8
.
38.
Lau
WM
,
Doucet
M
,
Stadel
R
,
Huang
D
,
Weber
KL
,
Kominsky
SL
. 
Enpp1: a potential facilitator of breast cancer bone metastasis
.
PLoS One
2013
;
8
:
e66752
.
39.
Mackenzie
NC
,
Zhu
D
,
Milne
EM
,
van ‘t Hof
R
,
Martin
A
,
Darryl Quarles
L
, et al
Altered bone development and an increase in FGF-23 expression in Enpp1(-/-) mice
.
PLoS One
2012
;
7
:
e32177
.
40.
Carozza
JA
,
Brown
JA
,
Bohnert
V
,
Fernandez
D
,
AlSaif
Y
,
Mardjuki
RE
, et al
Structure-aided development of small-molecule inhibitors of ENPP1, the extracellular phosphodiesterase of the immunotransmitter cGAMP
.
Cell Chem Biol
2020
;
27
:
1347
58
.
41.
Cogan
D
,
Bakhoum
SF
. 
Re-awakening innate immune signaling in cancer: the development of highly potent ENPP1 inhibitors
.
Cell Chem Biol
2020
;
27
:
1327
8
.
42.
Corrales
L
,
Glickman
LH
,
McWhirter
SM
,
Kanne
DB
,
Sivick
KE
,
Katibah
GE
, et al
Direct activation of STING in the tumor microenvironment leads to potent and systemic tumor regression and immunity
.
Cell Rep
2015
;
11
:
1018
30
.
43.
Hoadley
KA
,
Yau
C
,
Hinoue
T
,
Wolf
DM
,
Lazar
AJ
,
Drill
E
, et al
Cell-of-origin patterns dominate the molecular classification of 10,000 tumors from 33 types of cancer
.
Cell
2018
;
173
:
291
304
.
44.
Thorsson
V
,
Gibbs
DL
,
Brown
SD
,
Wolf
D
,
Bortone
DS
,
Ou Yang
TH
, et al
The immune landscape of cancer
.
Immunity
2018
;
48
:
812
30
.
45.
Love
MI
,
Huber
W
,
Anders
S
. 
Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2
.
Genome Biol
2014
;
15
:
550
.
46.
Subramanian
A
,
Tamayo
P
,
Mootha
VK
,
Mukherjee
S
,
Ebert
BL
,
Gillette
MA
, et al
Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles
.
Proc Natl Acad Sci U S A
2005
;
102
:
15545
50
.