Progression of epithelial cancers predominantly proceeds by collective invasion of cell groups with coordinated cell–cell junctions and multicellular cytoskeletal activity. Collectively invading breast cancer cells express the gap junction protein connexin-43 (Cx43), yet whether Cx43 regulates collective invasion remains unclear. We here show that Cx43 mediates gap-junctional coupling between collectively invading breast cancer cells and, via hemichannels, adenosine nucleotide/nucleoside release into the extracellular space. Using molecular interference and rescue strategies, we identify that Cx43 hemichannel function, but not intercellular communication, induces leader cell activity and collective migration through the engagement of the adenosine receptor 1 (ADORA1) and AKT signaling. Accordingly, pharmacological inhibition of ADORA1 or AKT signaling caused leader cell collapse and halted collective invasion. ADORA1 inhibition further reduced local invasion of orthotopic mammary tumors in vivo, and joint up-regulation of Cx43 and ADORA1 in breast cancer patients correlated with decreased relapse-free survival. This identifies autocrine purinergic signaling, through Cx43 hemichannels, as a critical pathway in leader cell function and collective invasion.

Collective invasion of multicellular groups with intact cell–cell junctions is an important mechanism to local invasion and metastasis in epithelial cancers (Gaggioli et al., 2007; Friedl and Gilmour, 2009; Cheung and Ewald, 2016). Collective invasion is initiated and maintained by leader cells that extend actin-rich protrusions, engage with the ECM to exert traction and ECM remodeling, and maintain cell–cell coupling with follower cells (Cheung et al., 2013; Westcott et al., 2015). Mechanical connection and cell-to-cell signaling between moving cells are provided by adherens junctions, through E-, N-, and other cadherins (Theveneau and Mayor, 2012). Cadherins coordinate the actomyosin cytoskeleton through catenins and vinculin (Hartsock and Nelson, 2008; Han and de Rooij, 2016) and secure supracellular force balance, front–rear and apicobasal polarity, and multicellular branching (Tambe et al., 2011).

In nonneoplastic collective processes during morphogenesis and regeneration, moving cell sheets are further connected by gap junctions (Ashton et al., 1999; Huang et al., 1998; Kotini et al., 2018; Marins et al., 2009). Gap junctions consist of connexins (Cxs) oligomerized into hemichannels that engage across cell membranes (Goodenough and Paul, 2009). The resulting transmembrane connections mediate gap-junctional transfer of ions and small molecules (<1 kD), including Ca2+, phosphoinosites, and nucleotides (Goodenough and Paul, 2009). Cxs mediate multicellular contractility of cardiomyocytes and the myoepithelial layer of mammary ducts (Kumai et al., 2000; Mroue et al., 2015) and, by unknown mechanisms, collective movement of neural and endothelial cells during morphogenesis (Ashton et al., 1999; Huang et al., 1998; Marins et al., 2009). The mechanisms by which Cxs support coordinated cytoskeletal contractility and multicellular dynamics vary. Direct cell-to-cell signaling occurs through gap-junctional intercellular transfer of second messengers, including IP3 and cAMP, energy equivalents (ATP and glucose), and Ca2+ wave propagation (Boitano et al., 1992; Goldberg et al., 2004; Howe, 2004). Cxs further contribute to gene expression of cadherins (Kotini et al., 2018) and/or the release of chemotactic factors to induce cell polarity and migration (Barletta et al., 2012; Haynes et al., 2006; Kaczmarek et al., 2005). In morphogenesis, Cxs are indispensable for coordinated tissue growth (Sinyuk et al., 2018). By regulating cell differentiation, Cxs further counteract neoplastic transformation (Bazzoun et al., 2019; Fostok et al., 2019; Saunders et al., 2001; Zhang et al., 2003), and forced expression of Cxs in transformed cells reduces tumor cell growth and inhibits invasion by reverting the epithelial-to-mesenchymal transition (Kazan et al., 2019; McLachlan et al., 2006). On the other hand, Cxs remain expressed in several solid tumors, and their expression increases in metastases (Bos et al., 2009; Elzarrad et al., 2008; Kanczuga-Koda et al., 2006; Stoletov et al., 2013). Interference with Cx43 expression or channel function inhibits cancer cell migration in vitro (Ogawa et al., 2012) and reduces metastatic seeding by reducing binding of circulating tumor cells to vascular endothelial cells (el-Sabban and Pauli, 1991; Elzarrad et al., 2008; Stoletov et al., 2013). Because Cxs exert a range of functions, and these may vary between tumor types and experimental conditions, an integrating concept on how Cxs either suppress or enhance neoplastic progression is lacking (Aasen et al., 2019; Naus and Laird, 2010).

Whereas a role of Cxs in collective processes in morphogenesis has been established, the contribution of gap-junctional communication in collective tumor cell invasion remains unclear (Friedl and Gilmour, 2009). We here revisited Cx43 expression and function in breast cancer models of collective invasion. Using pharmacological inhibition, molecular interference, and rescue strategies, we show that Cx43 mediates gap-junctional intercellular communication (GJIC) between collectively invading cells and, via hemichannels, the release of purine derivatives into the extracellular space. Cx43 hemichannel-mediated release of adenosine (ADO) nucleotides initiate an autocrine loop inducing leader cell function and collective invasion, but not single-cell motility, through ADO receptor 1 (ADORA1) and AKT/protein kinase B signaling.

Breast cancer cells express Cx43 during collective invasion

To address whether Cxs are expressed in invasive ductal breast carcinoma infiltrating the fibrous tumor stroma and adipose tissue, Cx43 expression was detected in patient samples by multispectral imaging (Fig. S1, a and b; and Mascaux et al., 2019). Nonneoplastic ducts consisted of luminal epithelial cells with negative to low Cx43 expression and Cx43-positive myoepithelial cells (Fig. 1, a and b; and Fig. S1 c). In tumors, however, multicellular epithelial strands or roundish clusters with cytokeratin expression located in the fibrous or adipose tissue showed increased Cx43 expression in 9 out of 12 samples, and expression was increased compared with normal luminal epithelium (Fig. 1 a, arrowheads; Fig. 1 c; and Fig. S1, c and d). The fraction of Cx43bright cells, as quantified by single-cell in situ cytometry in Cx43-positive samples (Fig. 1 c, red label), varied from 5% to 75% (Fig. S1 d). This indicates that invasive ductal carcinomas contain subregions of high Cx43 expression in the invasion zones. Using meta-analysis of the Györffy cohort of breast cancer patients (Györffy et al., 2010), high Cx43 gene expression was associated with decreased relapse-free survival (RFS) in patients with estrogen receptor (ER)–negative subtypes, including the basal subtype (Fig. 1 d), and Her2-positive tumors (Fig. S1 e). No inverse association was detected when patients were classified based on tumor grade or on the ER-positive subtype (Fig. S1 e).

To address whether Cx43 supports invasion, the ER-negative metastatic breast cancer cells 4T1, which express E-cadherin, and mouse mammary tumor (MMT) cells, which lack E-cadherin but expresses N-cadherin (Ilina et al., 2018), were embedded as multicellular spheroids in 3D collagen lattices and monitored for invasion over 24–48 h (Videos 1 and 2). 4T1 and MMT cells invaded as collective strands and expressed high levels of Cx43 and Cx45 mRNA and low levels of other Cxs after isolation from collagen matrix culture (Fig. S1 f). During collective invasion, Cx43 was localized along cell–cell junctions (Fig. 1, e and f, white arrowheads) and in nonjunctional domains, along cell protrusions (Fig. 1, e and f, red arrowheads; and Fig. 1 g). Junctional and nonjunctional Cx43 localization was further observed during invasion into the mammary fat pad in vivo (Fig. 1 h). In contrast to Cx43, Cx45 was absent along cell–cell junctions but localized in the cytosol (Fig. S1 g). Thus, Cx43 is expressed at the surface of collectively invading breast cancer cells in vitro and in vivo.

Cx43 mediates GJIC during collective invasion

Collective invasion in 3D spheroid culture was initiated by single or few leader cells, which extended actin-rich protrusions to the front and maintained cadherin-based junctions with following cells (Fig. S2, a and b). To examine whether gap junctions are active during collective invasion, intercellular dye transfer was measured in 3D spheroid invasion cultures using gap FRAP (Video 3 and Fig. 2 a). Within minutes after photobleaching, leader and follower cells recovered the calcein signal by 20–30% (Fig. 2, b and c), with unbleached neighboring cells acting as donors (Video 3; and Fig. S2, c and d). As internal control, detached single cells lacked dye transfer after photobleaching (Fig. 2, b and c; Fig. S2, c and d, purple cell outline; Video 3). Fluorescence recovery was reduced by 70–85% after inhibition of Cx channel function by carbenoxolone (CBX; Fig. 2, b–d; and Video 3). Fluorescence recovery was further impaired when Cx43 expression was stably down-regulated (by >80%) using shRNA (Fig. 2, e and f; and Fig. S2, e–h).

To test whether all cells maintain GJIC or only specific cell subsets, dye transfer between freshly interacting breast cancer cells was measured (parachute assay; Abbaci et al., 2008; Fig. 2 g). Within a 12-h observation period, the majority of control 4T1 cells (90–95%) received dye from the calcein-loaded monolayer (Fig. 2 h). The extent of dye transfer was strongly reduced by 67–84% by CBX treatment or Cx43 down-regulation (Figs. 2 i and S2 i). These data indicate that collectively invading breast cancer cells interact by GJIC through Cx43.

Cx43 hemichannels release purine nucleotides

Besides gap junctions, Cxs can form transmembrane hemichannels (Ye et al., 2003), which release small molecules from the cytosol into the extracellular space, including glutamate, prostaglandins, and the purine nucleotides ATP and ADP (Baroja-Mazo et al., 2013; Eltzschig et al., 2006; Retamal et al., 2007). Extracellular ATP and ADP and their degradation product, ADO, are important energy equivalents. In addition, extracellular nucleotides/nucleosides are effective extracellular signaling molecules, which promote the directional migration in endothelial cells (Kaczmarek et al., 2005), microglia (Haynes et al., 2006), and neutrophils (Barletta et al., 2012), as well as in individually moving cancer cells (transwell assay; Zhou et al., 2015). In addition to its localization at cell–cell contacts between invading cells, Cx43-positive foci were present in protrusions of leader cells (Fig. 1, e, f, and h, red arrowheads; and Fig. 1 g). To address whether breast cancer cells release purine nucleotides during invasion, supernatants from 3D collagen cultures were analyzed for ADO, AMP, ADP, and ATP by HPLC. Whereas cell-free media contained no purine metabolites, 3D spheroid invasion cultures after 24 h contained high levels of extracellular ATP (0.1–0.25 µg/ml) and moderate levels of ADO, AMP, and ADP (10–50 ng/ml; Fig. 2 j).

Inhibition of Cx channel-function by CBX and Cx43 knockdown both reduced the extracellular levels of ADO (by 26–42%), AMP (by 40–41%), ADP (by 60–64%), and ATP (by 28–52%; Figs. 2 k and S2 j). This indicates that breast cancer cells release purine nucleosides and nucleotides through Cx43 hemichannels in 3D invasion culture.

Cx43 is required for leader cell function

To determine whether Cx43 function, including GJIC and hemichannel activity, is required for collective cell invasion, 3D spheroid cultures were treated with the pharmacological channel blockers CBX or 18 α-glycyrrhetinic acid (18aGA) at doses that did not compromise growth (Fig. S3, a and b). Cultures were treated either before (prevention) or after collective invasion has started (intervention). Both compounds dose-dependently reduced strand initiation and impaired the emergence of leader cells compared with the inactive homologue glycyrrhizic acid ammonium salt (GLZ) or vehicle control (Fig. 3 a, arrowheads, Fig. 3, b and c; Fig. S3, c and d, arrowheads; and Fig. 3 e). Consequently, the invasion speed and cumulative length of collective invasion strands after 24 h were reduced by up to 10-fold (Fig. 3 b; and Fig. S3, f–h).

Likewise, stable down-regulation of Cx43 by shRNA compromised both the number and length of collectively invading strands (Fig. 3, d–g). To rule out adaptation or off-target effects of stable Cx43 down-regulation, transient Cx43 down-regulation using individual RNAi probes was used, leading to reduced invasive strand initiation and elongation (Fig. 3, h and i). Neither pharmacological nor RNAi with GJIC altered cell–cell cohesion or shape. Invasion strands, emerging despite interference, retained intact morphology and cortical F-actin organization (Fig. 3 c, arrowheads; and Fig. 3 d).

The decreased initiation and prolongation of invasion strands after interference with Cx43 suggested that the activity of leader cells was compromised. Leader cells are critical in initiating and maintaining collective invasion (Cheung et al., 2013; Zhang et al., 2019). To test whether leader cell function depends on Cx43, fluorescent mosaic spheroids of cells expressing nontargeting or Cx43 shRNA at equal ratio were allowed to invade. Cells with down-regulated Cx43 were 50% less likely to acquire leader cell position (Fig. 3, j and k). CBX treatment reduced the speed of leader cells as assessed by cell tracking (Fig. S3 f and Video 4). Besides decreasing leader cell initiation, intermittent inhibition of GJIC by CBX in spheroids with established collective invasion affected leader cell function. CBX caused near-instantaneously collapsing protrusive tips, followed by speed reduction, and, ultimately, strand retraction (Fig. 3 l and Video 5). Leader cell functions were restored after CBX washout (Video 5), indicating that Cx channels were required for the initiation and maintenance of collective invasion.

We next aimed to understand which of the two channel functions of Cx43, GJIC or hemichannel activity, is required for collective invasion. Therefore, we bypassed hemichannel function by exogenous addition of ADO in Cx-inhibited cultures and tested whether invasion was restored. Exogenous ADO was added to 4T1 and MMT spheroid cultures in the presence of CBX or after Cx43 down-regulation. Increasing concentrations of ADO partially (35–50%) restored the frequency of leader cells compared with untreated cells (Fig. 4, a–c; and Fig. S4 a). ADO substitution further restored defective leader cell function caused by down-regulation of Cx43 in mosaic spheroids (Fig. 4 d, green cells; and Fig. 4 e). These data suggest that nucleotide/nucleoside release via hemichannels is essential, while GJIC is dispensable for collective invasion. The restored leader cell function by exogenous ADO thus indicates an autocrine Cx43-hemichannel-dependent purinergic signaling loop to secure collective migration.

Autocrine purinergic receptor signaling maintains collective cell invasion

The signaling function of extracellular ADO depends on its binding to P1 purinergic receptors (Barletta et al., 2012). P1 receptors belong to the G-protein–coupled family of receptors, with the members ADORA1, ADORA2a, ADORA2b, and ADORA3 (Antonioli et al., 2013). 4T1 and MMT cells isolated from 3D invasion culture expressed ADORA1 and ADORA2b mRNA but no other ADO or ADP receptor (Fig. S4 b). We thus tested whether ADORA1 or ADORA2b mediate collective invasion using pharmacological interference. PSB36 and 8-cyclopentyl-1,3-dipropylxanthine (DPCPX), both selective ADORA1 antagonists with established activity in vitro and in vivo (Abo-Salem et al., 2004; Töpfer et al., 2008), caused dose-dependent inhibition of collective invasion in 4T1 and MMT cells (Fig. 5, a–c; Fig. S4 c; and Fig. S4 d). PSB36 did not affect the mitotic frequency of the cells (Fig. S4 e). Unlike ADORA1, interference with ADORA2b using the antagonist PSB1115 showed no effect on the emergence of leader cells of 4T1 and MMT cells (Fig. S4 f).

PSB36 abolished the effects of exogenous ADO on rescuing collective invasion (Fig. 5, d and e) and prevented leader cells to develop pointed anterior protrusions (Video 6). As a consequence of leader cell collapse, invasion strands retracted, and this was reversible after washout of the inhibitor (Video 6). Thus, leader cell function in collective invasion critically depends on ADO signaling and ADORA1.

ADORA1 localization during collective invasion in vitro was predominantly cytosolic, and its expression levels were not affected by CBX treatment (Fig. S4, g and h). The localization of ADORA1 is in line with previously reported cytosolic distribution of ADORA1 and other chemokine receptors (G-protein–coupled receptors) such as in nerve terminals and gastric cancer tissues (Iwasa et al., 2009; Rebola et al., 2003).

To address whether Cx43/ADORA1 regulates collective movement of nontumorigenic cells, we embedded spheroids of the normal mammary cell line NMuMG (Nemir et al., 2000) in 3D collagen. Normal Mus muculus mammary gland (NMuMG) cells are deficient in both junctional expression of Cxs (including Cx43) and GJIC (Sheng et al., 2007). Consistent with defective GJIC, Cx43 was enriched in the cytosols of NMuMG leader cells guiding collective invasion strands (Fig. S5 a, arrowheads). Treatment with the channel blocker CBX or ADORA1 inhibitor PSB36 reduced invasion of NMuMG (Fig. S5, b–e), suggesting that Cx43 hemichannels and extracellular ADO signaling rather than GJIC support collective movement. Thus, the autocrine nucleotide loop is not restricted to cancer cells but may represent a general collective steering mechanism.

We next tested whether single-cell invasion is also regulated by Cx hemichannels and ADORA1. 4T1 cells were individualized and embedded in 3D collagen and their movement followed by cell tracking over a time course of 40 h. Although individual cells maintained the release of nucleotides and nucleosides into the culture supernatant (Fig. S5 f), only 20–30% of individualized 4T1 cells developed elongated morphology, and single cells only migrated on the spot with an average speed of 0.009 µm/min, which is at least six times below the average speed of collectively invading 4T1 cells (0.06 µm/min; Fig. S5 g). Treatment with the pan-Cx inhibitor CBX or the ADORA1 inhibitor PSB36 did not affect the number of protrusive cells (Fig. S5, h and i). Thus, the hemichannel–ADORA loop does not support cell protrusion of individual breast cancer cells.

AKT phosphorylation in Cx/ADORA-mediated invasion

In 2D sheet migration of nontransformed epithelial cells, leader cell function depends on phosphoinositide 3-kinase and AKT signaling (Yamaguchi et al., 2015). Since AKT activation is a major event downstream of ADO signaling via ADORA1 (Park et al., 2010; Umapathy et al., 2013), we investigated the involvement of AKT in Cx/ADORA1-mediated invasion. In 3D invasion culture, the active form of AKT (AKT-pSer473) was abundant near the protrusion tips of leaders, but not in follower cells (Fig. 6 a; and Fig. S5, j and k). The focal AKT-Ser473 intensity in leader cells was reduced after addition of CBX (Fig. 6, a–d). Likewise, PSB36 reduced AKT-pSer473 levels (Fig. S5 l), whereas AKT-pSer473 levels increased within 1 h of ADO addition (Fig. S5 m), in line with previous reports showing time-dependent regulation of AKT phosphorylation in response to purinergic signaling (Jacques-Silva et al., 2004). To determine whether AKT activation is required for invasion, we treated 4T1 and MMT 3D invasion cultures with the AKT inhibitor IV. The emergence of leader cells and strand formation in both 4T1 and MMT spheroids were reduced by >84%, compared with the DMSO controls (Fig. 6, e and f). These data suggest that the hemichannel/ADO loop regulates collective invasion by activating AKT in leader cells.

ADORA1 in breast cancer invasion and progression in vivo

To address whether ADORA1 signaling supports collective breast cancer cell invasion in vivo, multicellular tumoroids were implanted into the mammary fat pad and the effect of systemic administration of PSB36 on progressing invasion was monitored by intravital multiphoton microscopy through a body window for 4 d (Ilina et al., 2018). Untreated tumors invaded into the mammary fat pad via collective strands with elongated leading tip cells (Fig. 7 a, arrowheads), as described previously (Ilina et al., 2018). PSB36 inhibited the extent of invasion with fewer strands with tips and reduced strand length (Fig. 7 a, arrowheads; and Fig. 7 b). PSB36 did not affect tumor growth, as assessed by the number of nuclei per spheroid (Fig. 7 c), and this resulted in more compact tumors with increased packing density of cells. These data establish ADORA1 in promoting collective breast cancer cell invasion in vivo.

In tumor lesions in vivo, extracellular ADO levels are elevated to micromolar ranges (Blay et al., 1997; Kumar, 2013), and ADORA1 mRNA and protein expression were reported to be increased in whole-tissue lysates of primary tumors compared with normal tissue (Khoo et al., 1996; Mirza et al., 2005). To verify that ADORA1 was also present in the invasion zone of ductal breast cancer, clinical samples were stained for ADORA1, and the collective invasion patterns located in the fibrous and adipose tissues were analyzed by multispectral microscopy (Fig. 7 d). Multicellular strands and nests containing cytokeratin-positive epithelial cancer cells expressed ADORA1 (in six out of eight samples) with increased levels and heterogeneity compared with the luminal epithelial cells of normal ducts (Fig. 7, d and e; and Fig. S6 a). In ADORA1-positive samples, 14–68% of the cancer cells within the invasion patterns were ADORA1 bright (Fig. S6 b) with a predominant cytosolic distribution, consistent with the ADORA1 distribution detected in 3D invasion culture in vitro.

We finally tested whether inhibition of ADORA1 also reduces spontaneous metastasis to the lungs. 4T1 cells were implanted into the mammary fat pad, and the mice were treated daily with DMSO or PSB36. PSB36 therapy had no effect on the number or size of metastases to the lungs 30–36 d after tumor implantation (Fig. S6, c–e). Thus, although PSB36 inhibited local invasion in vivo, it had no effect on spontaneous metastasis in the 4T1 model. ADORA1 expression was inversely correlated with RFS in ER-negative patients with high Cx43 expression, but not ADORA1 alone (Fig. S6, f and g). This indicates that Cx43/ADORA1 expression correlates with poor outcome in ER-negative cancer, including the basal subtype.

We here identify an autocrine mechanism of leader cell activation and maintenance through Cx43 hemichannels, which release nucleotides that induce AKT activation in leader cells and collective invasion of breast cancer cells via a signaling cascade involving ADORA1. ADORA1 drives local invasion and dissemination from implanted tumors in mice, and ER-negative breast cancer patients with high Cx43/ADORA1 expression suffer from worsened prognosis. Autocrine nucleotide signaling thus provides a cell-autonomous mechanism sustaining collective invasion in breast cancer and may deliver a rational basis for targeted intervention using hemichannel and ADORA inhibitors.

The data indicate a previously unappreciated role for nucleotide release via Cx43 hemichannels and autocrine ADO signaling in inducing leader cells that guide collective invasion of breast cancer cells. The concept of hemichannel-mediated nucleotide release in our study is supported by (1) focal Cx43 localization at cytoplasmic extensions of leader cells and lateral sites of follower cells, (2) dependence of extracellular nucleotide levels in 3D tumor spheroid culture on Cx43 channel function, (3) the decreased probability of cells with low Cx43 expression to reach leader cell position, (4) the rapid loss of leader cell activity after Cx channel inhibition or antagonization of ADORA1, and (5) reversion of invasion inhibition by CBX or Cx43 down-regulation by exogenous ADO addition. These data strongly suggest Cx43 hemichannel-mediated nucleotide release as major pathway to collective invasion.

The release of nucleotides by Cx hemichannels reveals a tumor-cell–autonomous, autocrine mechanism of induction and continuation of collective invasion. Hemichannel/purinergic signaling loops are known to regulate polarized protrusions of microglial cells during response to brain injury (Davalos et al., 2005) and single-cell migration of neutrophils during chemotaxis (Wang et al., 2017). In contrast to microglia and leukocytes, the autocrine nucleotide loop was not sufficient to induce polarization and migration of individualized breast cancer cells. Several parameters may account for this differential hemichannel function in epithelial contexts. Collective, but not individual-cell, polarity is reinforced by cell–cell junctions, including cadherins and Cxs, and receptor–ligand intercellular interactions, such as notch/delta and ephrin/Eph signaling (Friedl and Mayor, 2017). Intercellular junctions and interactions mediate mechanical and biochemical transductions between cells that enhance promigratory signals such as ECM contact and paracrine factors. As example, cadherin-dependent mechanocoupling promotes front–rear asymmetry by directing protrusion orientation and ECM binding toward the cell-free edge (Khalil and de Rooij, 2019; Mayor and Carmona-Fontaine, 2010; Plutoni et al., 2016). Thus, to induce leader cell activity and collective invasion, purinergic signaling may cooperate with the cell–cell junction machinery.

Our data further indicate that Cx43 hemichannel activity is essential but GJIC is dispensable for leader cell function and collective invasion. Whether, beyond the short-term effects detected here, GJIC contributes to mid- and long-term sustainability of collective invasion remains to be addressed. By allowing cytoplasmic continuity between invading cells, Cx43 intercellular- and hemichannels may act in series and contribute to ADO delivery and availability to the invasion front for release into the extracellular space.

In the extracellular space, ATP and ADP are degraded into ADO by ectonucleotidases, including CD39 and CD73 (de Andrade Mello et al., 2017), and ADO levels in the high-micromolar range are present in the tumor microenvironment (Blay et al., 1997); therefore ADO released in 3D invasion cultures and in vivo acts as both, metabolite and signaling effector of cancer progression (Antonioli et al., 2013). The ADO levels released by 4T1 and MMT cells during in vitro culture were one to five orders lower than ADO levels detected in the extracellular space of the tumor microenvironment (10–100 µM; Vaupel and Mayer, 2016; Vaupel and Multhoff, 2016). It is difficult to compare ADO levels present in the tumor tissue to the levels detected in spheroid invasion culture. Extracellular ADO in intact tissues is contributed by both tumor and stromal cells (de Lourdes Mora-García et al., 2016) over long periods of time. The resulting steady-state level of ADO reached in vivo may thus be higher compared with levels reached in invasion cultures over 24 h, with cancer cells as the sole source of ADO. Leader cell induction by auto- or paracrine nucleotides depends on downstream signaling function through ADORAs to induce polarized cell–substrate interaction, and this is shared with nucleotide activity inducing single-cell migration in other contexts. AKT signaling, which is induced by G-protein–coupled receptors (ADORA1, ADORA2a, ADORA2b, and ADORA3), mediates polarization of moving cells downstream of ADORA signaling (Chen et al., 2006; Gao et al., 2001; Othman et al., 2003; Stagg et al., 2010; Umapathy et al., 2013; Wen et al., 2011). Active AKT was increased in leader, but not follower, cells and required both Cx channel and ADORA1 function. In single-cell migration, individualized keratinocyte migration along electrical gradient depends on the Cx/G-protein–coupled receptor axis (Riding and Pullar, 2015). In brain-damage responses, astrocyte processes protrude in a Cx-ADORA–dependent manner to interact and close the microdefect (Davalos et al., 2005; Riding and Pullar, 2015). Although both ADORA1 and ADORA2b are expressed, only ADORA1 mediated ADO-dependent invasion of 4T1 and MMT cells. Possibly, the reported 70-fold higher affinity of ADO to ADORA1 (Fredholm et al., 2001) accounts for its dominance in leader cell activation. However, in difference to ADORA-stimulated single-cell migration, autocrine stimulation of tip cells favors their engagement toward ECM while remaining coupled to follower cells, thus mediating the maintenance and elongation of collective invasion strands.

While ADORA1 function was clearly involved in local tissue invasion, and collective invasion in breast carcinoma samples is positively associated with clinical progression to metastasis (Khalil et al., 2017), spontaneous metastasis to the lungs was not affected by pharmacological ADORA1 antagonization. The 4T1 in vivo model resulted in 100% penetrance of spontaneous micro- and macrometastases, and the used dosing of PSB36 inhibited invasion in mammary fat pads, suggesting that the drug is efficient in vivo. The lack of inhibition of metastases may indicate that leader cell function maintained by ADORA1 is of lower relevance in secondary organ colonization and metastatic outgrowth. A similar disconnect between local invasion and metastatic ability in distant organs was recently shown after deletion of E-cadherin, which caused increased local invasion in breast cancer models but strongly compromised metastatic ability due to a survival defect (Padmanaban et al., 2019). Interfering with ADORA1 pathways may thus be justified to limit locally invasive disease, which cannot be effectively treated by surgery (Bakst et al., 2019). Beyond enhancing invasion, hemichannels and ADORAs also contribute to tumor growth and immunomodulation (Schalper et al., 2014) and orchestrate an integrated program to enhance local cancer progression. Purinergic signaling thus establishes a proinvasive microenvironment through multipronged mechanisms, rendering pharmacological interference of nucleotide signaling through the Cx/ADORA axis in patient subsets as an attractive route to reprogram the tumor stroma and neoplastic invasion.

Antibodies and reagents

The following antibodies were used: rabbit anti-human Cx43 (3512; CST), chicken anti-human vimentin (ab24525; Abcam), mouse anti-mouse β-catenin (clone 14/β-Catenin; BD Biosciences), rabbit anti-mouse AKT (9272; CST), mouse anti-human AKT (2920; CST); rabbit anti-human phospho-AKT Ser473 (4060; CST); mouse anti-human Pan Keratin (4545; CST); rabbit anti-rat ADORA1 (ab82477; Abcam); rabbit anti-human ADORA1 (ab124780; Abcam); and secondary polyclonal Alexa Fluor 488/647–conjugated goat anti-mouse, anti-rabbit, or anti-chicken antibody. For F-actin visualization, Alexa Fluor 568–conjugated phalloidin (Invitrogen) was used. ADO (A4036; Sigma-Aldrich) was dissolved in water at a concentration of 1 mg/ml and stored in aliquots at −20°C. The following inhibitors were used: gap junction channel inhibitors CBX disodium salt, 18aGA, and its inactive homologue, GLZ (Sigma-Aldrich; Kenny et al., 2002; Marins et al., 2009); the ADORA1 selective and competitive antagonists 1-butyl-3-(3-hydroxypropyl)-8-(3-noradamantyl) xanthine (PSB36) and 8-cyclopentyl-1,3-dipropylxanthine (Sigma-Aldrich); and the ADORA2b selective antagonist 4-(2,3,6,7-tetrahydro-2, 6-dioxo-1-propyl-1H-purin-8-yl)-benzenesulfonic acid (PSB1115; Tocris; Lin et al., 2010). CBX was dissolved in water (stock concentration: 163 mM); all other inhibitors and antagonists were dissolved in DMSO (stock concentration 100–200 mM). All drugs were aliquoted (5–20 µl) and stored at −20°C.

Cell lines and culture

The following metastatic breast cancer cell lines were used: 4T1 (CRL-2539 ATCC); MMT wild-type and MMT dual-color cells (H2b-GFP and cytosolic RFP; Tsuji et al., 2006; a generous gift from R. Hoffman, AntiCancer Inc., San Diego, CA); and NMuMG cells (a generous gift from M. Vooijs, Maastricht University Medical Center, Netherlands; Owens et al., 1974). The identity of 4T1, MMT and NMuMG cells was verified by short tandem repeat DNA profiling (IDEXX BioResearch). No mammalian interspecies contamination was detected. Cells were routinely tested for mycoplasma contamination (MycoAlert; Lonza). All cell lines were negative for mycoplasma contamination. Cells were maintained (37°C, 5% CO2, humidified atmosphere) in RPMI media (Invitrogen) supplemented with 10% fetal bovine serum (Sigma-Aldrich), penicillin (100 U/ml) and streptomycin (100 µg/ml; both PAA), and l-glutamine (2 mM; Invitrogen). For generation of 4T1 and MMT cells with stable Cx43 RNA down-regulation, mission lentiviral vector targeting mouse GJA1 (TRCN0000068473, sequence: 5′-CCG​GCC​CAC​CTT​TGT​GTC​TTC​CAT​ACT​CGA​GTA​TGG​AAG​ACA​CAA​AGG​TGG​GTT​TTT​G-3′) and nontargeting control sequence (5′-CCG​GCA​ACA​AGA​TGA​AGA​GCA​CCA​ACT​CGA​GTT​GGT​GCT​CTT​CAT​CTT​GTT​GTT​TTT-3′) was used (Sigma-Aldrich). For tumor implantation into the mammary fat pad in vivo, a stable 4T1/H2B-mCherry/LifeAct-GFP cells were used, as described previously (Ilina et al., 2018). Fluorescence-expressing cells did not deviate from the wild-type counterparts in parameters such as growth and invasion (data not shown).

Cancer spheroid and single-cell culture

Spheroids from 4T1, MMT, and NMuMG cells used in short-term in vitro migration assays were generated using the hanging drop method (Venhuizen et al., 2019). In brief, cells from subconfluent culture were detached with EDTA (1 mM) and trypsin (0.075%; Invitrogen), washed in PBS, suspended in medium/methylcellulose (2.4%; Sigma), and maintained as hanging droplets (25 µl), each containing 1,000 cells. Hanging drops were kept in 37°C, 5% CO2 humidified atmosphere for 24 h (Del Duca et al., 2004). Spheroids were harvested and placed into DMEM media without FCS and without penicillin/streptomycin in a six-well plate.

After cell aggregation, spheroids were washed in PBS and incorporated into 3D type I collagen lattices consisting of nonpepsinized rat-tail collagen (BD Biosciences) at a final concentration of 4 or 6 mg/ml for 4T1 or MMT cells (Ilina et al., 2018). In brief, spheroids were washed in complete medium, mixed with collagen solution, and prevented from sinking during collagen polymerization (37°C, 5% CO2) by rotation of the sample. Only spheroids invading into the 3D phase of the collagen matrix were analyzed, whereas spheroids moving as 2D sheets underneath or on top of the culture were excluded.

For single-cell invasion, cells from subconfluent culture were detached with EDTA (2 mM), washed and embedded in 3D nonpepsinized rat-tail collagen (4 mg/ml; BD Biosciences) at a density of 20,000 cells/100 μl collagen.

Transient down-regulation of Cx43 by RNAi

4T1 and MMT cells were seeded in a 12-well plate at a density of 60 × 103/well, in the presence of siRNA duplex (Dharmacon, Thermo Fisher Scientific). ON-Target plus nontargeting siRNA (D-001810-10) and two Cx43-specific targeting sequences (J-051694-05 and J-051694-06) were used at a final concentration of 50 nM using Dharmafect transfection reagent 4 (Thermo Fisher Scientific). After 12 h, RNAi containing supernatant was removed and replaced by penicillin/streptomycin free media. After a recovery period of at least 6 h, cells were used for spheroid generation. For assessment of knockdown efficiency, whole-cell lysates were analyzed after 24, 48, and 72 h after removal of siRNA by Western blot.

Resazurin cytotoxicity assay

Mitochondrial activity was assessed by resazurin labeling (Czekanska, 2011). 4T1 and MMT spheroids embedded in 3D collagen were incubated with resazurin-containing media (100 µg/ml) for 4h (at 37°C, 5% CO2, humidified atmosphere). Gels were gently washed with PBS, and cell-bound fluorescence was measured (excitation/emission 560/590 nm) using an ELISA reader.

Real-time quantitative PCR

Collagen gels containing spheroids were dissolved with Trizol (250 µl/100 µl gel; Invitrogen) and homogenized before chloroform (Merck) phase separation (Farhat, 2012). After centrifugation (15 min, 12 × 103g, 4°C), the upper phase was mixed with 70% vol/vol ethanol and transferred to an RNeasy MinElute spin column for RNA extraction according to the manufacturer’s instructions (Qiagen’s RNeasy kit). After RNA quality control, reverse transcription (SuperScript II Reverse transcription; Thermo Fisher Scientific) was performed. Primers for Cxs and purinergic receptors (Table S1) were designed using HomoloGene (http://www.ncbi.nlm.nih.gov) and validated in silico using Oligoanalyzer software and in vitro using positive controls (RNA isolated from mouse tissues including brain and heart). Quantitative PCR was performed using CFX96 Real-Time PCR Detection System with C1000 Thermocycler (Bio-Rad) and analyzed using the Bio-Rad CFX Manager software (version 2.0). RNA levels for the genes of interest were normalized to the pooled levels of three housekeeping genes (β-actin, GAPDH, and ΥWhaz).

Western blot analysis

Cells from 2D culture were lysed using SDS sample buffer (62.5 mM Tris-HCl, 2% SDS, 10% glycerol, 50 mM DTT, and bromophenol blue). Collagen gels containing spheroids were digested by highly purified type VII collagenase (C0773; Sigma-Aldrich), followed by whole-cell lyses using SDS buffer. Western blotting was performed by SDS-PAGE. Briefly, whole-cell lysates were loaded and separated on a 10% acrylamide gel with Tris-glycine running buffer followed by blotting to polyvinylidene difluoride membrane by wet transfer. Polyvinylidene difluoride membranes were blotted with monoclonal mouse anti-β-tubulin (E7; Department of Cell Biology, Radboud University Medical Center, Nijmegen, Netherlands) as loading control, anti-Cx43 (dilution 1:1,000), anti-AKT (1:1,000), and anti-phospho-AKT Ser473 (1:500) followed by fluorescence detection (Odyssey; LI-COR Biosciences) and densitometric analysis.

HPLC analysis for purine content in media

Supernatants from 3D single-cell or spheroid cultures (24 h) were used for purine analysis. Media contained 80°C heated FCS, which minimizes the degradation of nucleotides and nucleosides (Gendaszewska-Darmach et al., 2003). ADO, AMP, ADP, and ATP levels were determined by HPLC using a modified method from Bhatt et al. (2012). In brief, culture supernatant (1 ml) was mixed with chloroacetaldehyde (250 µl; 6x diluted in 1 M acetate buffer, pH 4.5; Sigma-Aldrich), followed by derivatization (60 min, 70°C, 500 rpm) using Thermomixer comfort (Eppendorf AG) and centrifugation (3 min, RT, 12,100 × g). For HPLC, supernatant (400 µl) was transferred to a HPLC vial and injected. Purines were separated by HPLC (1200 Series; Agilent Technologies) using a Polariz C18-A column (150 × 4.6 mm) with a gradient elution using eluent A (0.1 M K2HPO4, 10 mM tetrabutylammonium hydrogen sulfate, pH 6.5, and 2% MeOH) and eluent B (H2O/ACN/THF 50:49:1). Retention times were 6.5 (ADO), 7.7 (AMP), 11.8 (ADP), 15.6 (ATP), and 15.0 min (cAMP). Quantification was based on peak areas of the samples, and reference standards were measured with fluorescence detection (excitation, 280 nm; emission, 420 nm).

Bright-field microscopy

Single or collective invasion was detected by digital time-lapse bright-field microscopy, as described previously (Wolf et al., 2013). In brief, single cells or 3D spheroid cultures were maintained at 37°C and images were taken with 5-min intervals for up to 72 h. For pharmacological intervention, spheroids were allowed to migrate for 24 h, followed by addition of inhibitor or vehicle (solvent) for 24 h and washout using 1 ml of growth media and follow-up for additional 24 h. Single cells and leader cells were tracked using the manual tracking plugin from ImageJ (1.40v; National Institutes of Health). The number, length, and elongation speed of invading single cells or strands were quantified using ImageJ. For endpoint experiments, bright-field micrographs of spheroid cultures in 3D collagen after 15, 24, and 48 h of incubation were recorded using a CCD camera (Sentech), 10x/0.20 NA air objective (Leica) and frame grabber system (Vistek). The number of pointed and protrusive tips (leader cells) per spheroid was counted manually; length of invading strands were quantified using the image analysis software ImageJ.

Parachute assay and FACS analysis

Receiver 4T1 cells (150 × 103 cells/well) were incubated for 24 h (37°C, 5% CO2 humidified atmosphere) in 48-well plates to reach 95–100% confluency. Donor 4T1 cells were loaded with calcein-AM (1 µM; Life Technologies) for 30 min, washed, and incubated for a further 45 min (Abbaci et al., 2008). Calcein-labeled cells were detached with EDTA (1 mM) and trypsin (0.075%), centrifuged (251 g, 5 min), resuspended in media, and 75 × 103 cells/well were added over the receiver cell culture at a ratio of 1:2 (donor/receiver) and co-cultured (3, 6, or 12 h). For analysis of dye transfer, cells were detached with EDTA (1 mM) and trypsin (0.075%), assayed by flow cytometry (FACS Caliber), and analyzed using FCS Express software.

Immunofluorescence staining and confocal microscopy

3D spheroid samples were fixed with 4% paraformaldehyde at room temperature (RT) for 15 min. For detection of intracellular antigens, gels were washed with PBS, permeabilized using 5–10% normal goat serum, 0.3% Triton X-100 or 0.2% Saponin in PBS (1 h, 20°C), incubated with antibody in PBS containing 0.1% BSA, 0.3% Triton X-100, or 0.2% Saponin. Antibody dilutions were as follows: Cx43 (1:50), AKT (1:50), S473 pAKT (1:50), β-Catenin (1:200), ZO-1 (1:100), and vimentin (1:200). Primary antibody was incubated at 4°C with shaking overnight, followed by at least five washing steps of 15 min each with PBS. Samples were then incubated with secondary antibody (1:400) in PBS buffer for intracellular staining, together with DAPI (5 µg/ml) and phalloidin (1:200; 4 h, RT), washed at least five times (15 min each) and imaged by confocal microscopy (Olympus FV1000; long working distance objectives 20×/NA 0.50, 40×/NA 0.80, or 60×/NA 1.35). Subcellular localization and quantification of Cx43, β-catenin, AKT, and pAKT at cell–cell contacts and cytoplasmic extensions were determined by manually outlining F-actin–positive regions followed by background correction using ImageJ.

For time-lapse confocal microscopy, 3D spheroid cultures consisting of MMT dual-color cells or 4T1/H2B-mCherry within 3D collagen were imaged (37°C, 5% CO2, humidified atmosphere). Time-lapse confocal microscopy (20×/NA 0.5 air objective, LSM 510 system; Carl Zeiss) was performed with time intervals of 10–15 min for 24 h with usually four to seven z-scans comprising all of the strands per spheroid (∼100 µm in depth). The confocal slices were reconstructed as maximum intensity projection and image analysis, including leader cell tracking, was performed using ImageJ. The minimum and maximum pixel intensity values of different samples were equally adjusted per experiment using the color balance tool of ImageJ.

Gap FRAP

Cancer cells grown as monolayers on glass-bottom dishes or as multicellular spheroids in 3D collagen were labeled with calcein-AM (2–3 µM, 30 min) using a modified method described previously (Kuzma-Kuzniarska et al., 2014). Regions of interest were defined for single cells or small regions and bleaching was performed using 100% Argon laser output 1 mW with 200 iterations using a large pinhole (∼8 airy units) and 20-s interval at 37°C, 5% CO2 and humidified atmosphere. For 12 min of follow-up, calcein transfer between cells seeded on 2D coverglass or within collagen gels was monitored using confocal microscopy (LSM 510; 20×/NA 0.5 air; Carl Zeiss; 37°C, 5% CO2). Inhibition of GJIC was achieved using CBX (100 µM) at least 45 min before gap FRAP imaging. Time-resolved bleaching and postbleaching recovery was quantified using ImageJ as the mean fluorescence intensity. For normalizing fluorescence after photobleaching, cells in nonbleached control regions were coregistered and the signal used for background correction. Detached single cells were used as an internal negative control for GJIC-dependent fluorescence recovery.

Breast cancer tissue section immunofluorescence

For multispectral imaging of Cx43 and ADORA1, tissue sections with positive invasion margin were selected from 19 breast cancer patients (Table 1). The use of coded tumor tissue was approved by the institutional review board, according to national law (Nagelkerke et al., 2011). Formalin-fixed paraffin-embedded breast tissue sections (5 µm thickness) were deparaffinized, followed by antigen retrieval using Tris-EDTA buffer (15 min, 95–100°C) followed by a washing step with PBS. Sections were blocked with 5% normal goat serum in 0.2% Tween for 1 h at 20°C. Primary antibody was incubated at 4°C with shaking overnight using the following dilutions: Cx43 (1:100), vimentin (1:200), pan-cytokeratin (1:200), and ADORA1 (1:200). Samples were then washed four times (PBS, 5 min each), incubated with secondary Alexa Fluor 488/546/647–conjugated goat anti-mouse, anti-rabbit, and anti-chicken antibodies, together with DAPI (1 μg/ml), for 1 h at 20°C and again washed four times (PBS, 5 min each). Dried sections were embedded in Mowiol (Sigma-Aldrich) and scanned using the automated Vectra Intelligent Slide Analysis System (version 2.0.8; PerkinElmer).

Multispectral imaging and quantitative digital analysis

Tissue slides were imaged using Vectra Intelligent Slide Analysis System (Version 2.0.8; PerkinElmer), as described (Mascaux et al., 2019). This imaging technique combines imaging with spectroscopy by collecting the spectral range from 400 to 720 nm in 10-nm steps in an automatic manner. To define the spectrum of each fluorophore as well as of tissue autofluorescence, stained single-fluorescence and unstained tissues were recorded with the Nuance Multispectral Imaging System (Version 3.0.2; PerkinElmer). A spectral library was built from the extracted spectra to enable the quantitative separation and spectral unmixing, thus removing cross-talk between fluorophores and interfering tissue autofluorescence. A selection of 10 representative original multispectral images was loaded into an advanced user-trainable morphological image analysis software (InForm Version 2.1; PerkinElmer) which utilizes machine learning for user-trained assignment of morphological and spectral intensity patterns to discriminate tumor cells, stromal regions, and normal ducts. The analysis software utilized the multimarker, multispectral data to extract and calculate single-cell parameters, including nuclear, cytoplasmic, and membrane intensity values used to classify each cell (e.g., as luminal, cancerous, or stroma cell; Fig. S1 b). After training, the algorithm was applied as batch analysis of multiple original multispectral images of different samples. Cx43, pan-cytokeratin, and vimentin expression levels were quantified in normal bilayered ducts, tumor cells located within the fibrous and adipose tissue, and their surrounding stroma. The expression levels of Cx43 were measured in each cell and the positivity threshold was determined by receiver operating characteristic (ROC) analysis whereby Cx43 and ADORA1 values of all the luminal epithelium cells available in different patients were used as negative control. With 99% specificity and 50% sensitivity for Cx43 and 93% specificity and 31% selectivity for ADORA1, the minimum threshold positive staining was calculated to be 0.875 and 2.500, respectively (Figs. 1 c and 7 e, dashed line).

In silico analysis of Cx43 and ADORA1 expression in breast cancers

We used the publicly available database (http://kmplot.com), which contains expression data on several genes, including Cx43 and ADORA1, in a cohort consisting of 3,951 breast cancer patients. The database classifies patients by trichotomizing the expression levels of the gene of interest, into high (T3) or low expression (T1), to extract correlations between gene expression level and patient outcome.

Mice

BALB/c female mice (6–8 wk old) were purchased from Charles River Laboratories. All animal procedures were approved by the Ethical Committee on Animal Experiments at Radboud University, Nijmegen (protocol number 140201) and performed according to guidelines of the Animal Welfare Committee of the Royal Netherlands Academy of Arts and Sciences, the Netherlands.

Implantation of tumor cell spheroids in the mammary fat pad

All surgical procedures were performed under isoflurane inhalation anesthesia (1–2% isoflurane/O2 mixture) and performed under conditions of microsurgery, as described (Ilina et al., 2018). In brief, 4T1 dual-color spheroids (500 cells/spheroid) were implanted adjacent to the third mammary fat pad into a central area of the mammary imaging window by nontraumatic microimplantation. After implantation, a sterile coverglass was placed into the window frame, and the integrity of both implanted tumoroid and stroma was verified by fluorescence microscopy.

Intravital microscopy and image analysis

Tumor growth and cell invasion were monitored for up to 4 d by intravital epifluorescence and multiphoton microscopy, as described (Ilina et al., 2018). In brief, mice were anesthetized with isoflurane, and the imaging window was stably fixed on a temperature-controlled stage (37°C) using a custom holder. Blood vessels were visualized by i.v. injection of 70-kD dextran labeled with Alexa Fluor 750. Invasion was detected using multiphoton microscopy (TrimScope II, LaVision BioTec; 20x objective, NA 0.95). Images were analyzed to quantify both, the number of nuclei and the surface/invasion area of implanted tumors at different time points using ImageJ.

Spontaneous metastasis analysis

The development of spontaneous lung metastasis was obtained from orthotopically implanted mammary tumors (Ilina et al., 2018). In brief, 4T1 cells (105 in 50 µl) were injected into the fourth mammary fat pad of BALB/c 8-wk-old female mice (Charles River Laboratories). DMSO (0.2%) and PSB36 (30 mg/kg) dissolved in PBS (Bilkei-Gorzo et al., 2008) were injected intraperitoneally once per day. Tumor size was recorded weekly by a caliper, and mice were euthanized 30–36 d after cancer cell injection or at the humane endpoint (tumor volume of 2 cm3). Lungs were harvested, fixed in buffered formalin, embedded in paraffin, and sectioned (5 µm thickness). Metastatic foci were detected by immunohistochemistry of cytokeratin-8 and quantified by operator-based blinded scoring (Ilina et al., 2018). In brief, the number of neoplastic events was counted from sequential sections and represented as the number of metastatic foci per slice. Metastasis was classified as micrometastasis (3–15 cytokeratin-8–positive cells) and macrometastasis (>15 tumor cells in nodular topology).

Statistical analysis

Graph Prism (version 5.0) was used for statistical analysis. The two-tailed unpaired Mann–Whitney (two groups) or Kruskal–Wallis test (more than two groups) with Dunns’ correction were used for all invasion analysis in vitro (number of leader cells/spheroid, cumulative invasion distance, and average speed). As an exception, data in Fig. S3 (e and h) were analyzed using ANOVA with Bonferroni correction after confirming normal distribution of the data (Kolmogorov–Smirnov normality test). The two-tailed unpaired Mann–Whitney test was further used for comparing invasion and tumor growth between DMSO and PSB36 treated mice in vivo. Data from Resazurin assay and HPLC analysis were tested for normal distribution using Kolmogorov–Smirnov test followed by ANOVA analysis with Bonferroni correction. Kaplan–Meier survival curves were compared using log-rank test.

Online supplemental material

Fig. S1 shows the expression of Cx43 in collective invasion patterns within breast lesions and Cx expression profile of 4T1 and MMT cells in 3D collagen (related to Fig. 1). Fig. S2 shows that collectively invading cancer cells maintain Cx43-mediated GJIC and extracellular nucleotide release (related to Fig. 2). Fig. S3 shows the dependence of leader cell function on Cx43 expression and channel function (related to Fig. 3). Fig. S4 shows purinergic receptor mRNA expression and the effect of ADORA inhibition on invasion and proliferation (related to Figs. 4 and 5). Fig. S5 shows the requirement of the Cx hemichannel-ADORA1 loop in other collective, but not single-cell, invasion models and downstream AKT activation (related to Fig. 6). Fig. S6 shows ADORA1 in normal human mammary epithelium, its role in preclinical spontaneous metastasis, and its prognostic relevance in clinical samples (related to Fig. 7). Video 1 shows 4T1 protrusive tips during initiation of collective invasion in 3D collagen I by bright-field time-lapse microscopy. Video 2 shows the initiation and elongation of collective invasion of MMT cells (H2b-GFP and cytosolic dsRed) in 3D collagen I by time-lapse confocal microscopy. Video 3 shows calcein intercellular dye transfer between collectively invading strands in the presence and absence of CBX. Video 4 shows MMT leader cells driving invasion (H2b-GFP) in 3D collagen I in the presence and absence of CBX. Video 5 and Video 6 show changes in the dynamics of established strands after intervention with CBX and PSB36 treatment, respectively. Table S1 shows primer sequences.

We gratefully acknowledge Huib Croes, Gert-Jan Bakker, and Dagmar Verweij for technical support and Michael Dustin and Pavlo G. Gritsenko for helpful discussions.

This work was supported by a German Excellence Initiative grant to the Graduate School of Life Sciences, University of Würzburg (to A.A. Khalil), the Netherlands Science Organization (NWO-VICI 918.11.626), the European Research Council (617430-DEEPINSIGHT), and the Cancer Genomics Center, Netherlands (to P. Friedl).

The authors declare no competing financial interests.

Author contributions: A.A. Khalil and P. Friedl designed the experiments. A.A. Khalil, O. Ilina, A. Vasaturo, J.-H. Venhuizen, M. Vullings, V. Venhuizen, and A. Bilos performed the experiments. A.A. Khalil, O. Ilina, A. Vasaturo, J.-H. Venhuizen, M. Vullings, V. Venhuizen, A. Bilos, and P. Friedl analyzed the data. C.G. Figdor provided the platform for multispectral imaging. P.N. Span provided and interpreted clinical samples. A.A. Khalil and P. Friedl wrote the manuscript. All authors read the manuscript.

Aasen
,
T.
,
E.
Leithe
,
S.V.
Graham
,
P.
Kameritsch
,
M.D.
Mayán
,
M.
Mesnil
,
K.
Pogoda
, and
A.
Tabernero
.
2019
.
Connexins in cancer: bridging the gap to the clinic
.
Oncogene
.
38
:
4429
4451
.
Abbaci
,
M.
,
M.
Barberi-Heyob
,
W.
Blondel
,
F.
Guillemin
, and
J.
Didelon
.
2008
.
Advantages and limitations of commonly used methods to assay the molecular permeability of gap junctional intercellular communication
.
Biotechniques
.
45
:
33
52–62
.
Abo-Salem
,
O.M.
,
A.M.
Hayallah
,
A.
Bilkei-Gorzo
,
B.
Filipek
,
A.
Zimmer
, and
C.E.
Müller
.
2004
.
Antinociceptive effects of novel A2B adenosine receptor antagonists
.
J. Pharmacol. Exp. Ther
.
308
:
358
366
.
Antonioli
,
L.
,
C.
Blandizzi
,
P.
Pacher
, and
G.
Haskó
.
2013
.
Immunity, inflammation and cancer: a leading role for adenosine
.
Nat. Rev. Cancer
.
13
:
842
857
.
Ashton
,
A.W.
,
R.
Yokota
,
G.
John
,
S.
Zhao
,
S.O.
Suadicani
,
D.C.
Spray
, and
J.A.
Ware
.
1999
.
Inhibition of endothelial cell migration, intercellular communication, and vascular tube formation by thromboxane A(2)
.
J. Biol. Chem
.
274
:
35562
35570
.
Bakst
,
R.L.
,
C.M.
Glastonbury
,
U.
Parvathaneni
,
N.
Katabi
,
K.S.
Hu
, and
S.S.
Yom
.
2019
.
Perineural Invasion and Perineural Tumor Spread in Head and Neck Cancer
.
Int. J. Radiat. Oncol. Biol. Phys
.
103
:
1109
1124
.
Barletta
,
K.E.
,
K.
Ley
, and
B.
Mehrad
.
2012
.
Regulation of neutrophil function by adenosine
.
Arterioscler. Thromb. Vasc. Biol
.
32
:
856
864
.
Baroja-Mazo
,
A.
,
M.
Barberà-Cremades
, and
P.
Pelegrín
.
2013
.
The participation of plasma membrane hemichannels to purinergic signaling
.
Biochim. Biophys. Acta
.
1828
:
79
93
.
Bazzoun
,
D.
,
H.A.
Adissu
,
L.
Wang
,
A.
Urazaev
,
I.
Tenvooren
,
S.F.
Fostok
,
S.
Chittiboyina
,
J.
Sturgis
,
K.
Hodges
,
G.
Chandramouly
, et al
.
2019
.
Connexin 43 maintains tissue polarity and regulates mitotic spindle orientation in the breast epithelium
.
J. Cell Sci
.
132
. jcs223313.
Bhatt
,
D.P.
,
X.
Chen
,
J.D.
Geiger
, and
T.A.
Rosenberger
.
2012
.
A sensitive HPLC-based method to quantify adenine nucleotides in primary astrocyte cell cultures
.
J. Chromatogr. B Analyt. Technol. Biomed. Life Sci
.
889-890
:
110
115
.
Bilkei-Gorzo
,
A.
,
O.M.
Abo-Salem
,
A.M.
Hayallah
,
K.
Michel
,
C.E.
Müller
, and
A.
Zimmer
.
2008
.
Adenosine receptor subtype-selective antagonists in inflammation and hyperalgesia
.
Naunyn Schmiedebergs Arch. Pharmacol
.
377
:
65
76
.
Blay
,
J.
,
T.D.
White
, and
D.W.
Hoskin
.
1997
.
The extracellular fluid of solid carcinomas contains immunosuppressive concentrations of adenosine
.
Cancer Res
.
57
:
2602
2605
.
Boitano
,
S.
,
E.R.
Dirksen
, and
M.J.
Sanderson
.
1992
.
Intercellular propagation of calcium waves mediated by inositol trisphosphate
.
Science
.
258
:
292
295
.
Bos
,
P.D.
,
X.H.
Zhang
,
C.
Nadal
,
W.
Shu
,
R.R.
Gomis
,
D.X.
Nguyen
,
A.J.
Minn
,
M.J.
van de Vijver
,
W.L.
Gerald
,
J.A.
Foekens
, et al
.
2009
.
Genes that mediate breast cancer metastasis to the brain
.
Nature
.
459
:
1005
1009
.
Chen
,
Y.
,
R.
Corriden
,
Y.
Inoue
,
L.
Yip
,
N.
Hashiguchi
,
A.
Zinkernagel
,
V.
Nizet
,
P.A.
Insel
, and
W.G.
Junger
.
2006
.
ATP release guides neutrophil chemotaxis via P2Y2 and A3 receptors
.
Science
.
314
:
1792
1795
.
Cheung
,
K.J.
, and
A.J.
Ewald
.
2016
.
A collective route to metastasis: Seeding by tumor cell clusters
.
Science
.
352
:
167
169
.
Cheung
,
K.J.
,
E.
Gabrielson
,
Z.
Werb
, and
A.J.
Ewald
.
2013
.
Collective invasion in breast cancer requires a conserved basal epithelial program
.
Cell
.
155
:
1639
1651
.
Czekanska
,
E.M.
.
2011
.
Assessment of cell proliferation with resazurin-based fluorescent dye
.
Methods Mol. Biol
.
740
:
27
32
.
Davalos
,
D.
,
J.
Grutzendler
,
G.
Yang
,
J.V.
Kim
,
Y.
Zuo
,
S.
Jung
,
D.R.
Littman
,
M.L.
Dustin
, and
W.B.
Gan
.
2005
.
ATP mediates rapid microglial response to local brain injury in vivo
.
Nat. Neurosci
.
8
:
752
758
.
de Andrade Mello
,
P.
,
R.
Coutinho-Silva
, and
L.E.B.
Savio
.
2017
.
Multifaceted Effects of Extracellular Adenosine Triphosphate and Adenosine in the Tumor-Host Interaction and Therapeutic Perspectives
.
Front. Immunol
.
8
:
1526
.
de Lourdes Mora-García
,
M.
,
R.
García-Rocha
,
O.
Morales-Ramírez
,
J.J.
Montesinos
,
B.
Weiss-Steider
,
J.
Hernández-Montes
,
L.R.
Ávila-Ibarra
,
C.A.
Don-López
,
M.A.
Velasco-Velázquez
,
V.
Gutiérrez-Serrano
, et al
.
2016
.
Mesenchymal stromal cells derived from cervical cancer produce high amounts of adenosine to suppress cytotoxic T lymphocyte functions
.
J. Transl. Med
.
14
:
302
.
Del Duca
,
D.
,
T.
Werbowetski
, and
R.F.
Del Maestro
.
2004
.
Spheroid preparation from hanging drops: characterization of a model of brain tumor invasion
.
J. Neurooncol
.
67
:
295
303
.
el-Sabban
,
M.E.
, and
B.U.
Pauli
.
1991
.
Cytoplasmic dye transfer between metastatic tumor cells and vascular endothelium
.
J. Cell Biol
.
115
:
1375
1382
.
Eltzschig
,
H.K.
,
T.
Eckle
,
A.
Mager
,
N.
Küper
,
C.
Karcher
,
T.
Weissmüller
,
K.
Boengler
,
R.
Schulz
,
S.C.
Robson
, and
S.P.
Colgan
.
2006
.
ATP release from activated neutrophils occurs via connexin 43 and modulates adenosine-dependent endothelial cell function
.
Circ. Res
.
99
:
1100
1108
.
Elzarrad
,
M.K.
,
A.
Haroon
,
K.
Willecke
,
R.
Dobrowolski
,
M.N.
Gillespie
, and
A.B.
Al-Mehdi
.
2008
.
Connexin-43 upregulation in micrometastases and tumor vasculature and its role in tumor cell attachment to pulmonary endothelium
.
BMC Med
.
6
:
20
.
Farhat
,
Y.
.
2012
.
RNA Extraction from Collagen Gels Using Qiagen’s RNeasy Kit
.
The Protocol Place
.
Fostok
,
S.
,
M.
El-Sibai
,
D.
Bazzoun
,
S.
Lelièvre
, and
R.
Talhouk
.
2019
.
Connexin 43 Loss Triggers Cell Cycle Entry and Invasion in Non-Neoplastic Breast Epithelium: A Role for Noncanonical Wnt Signaling
.
Cancers (Basel)
.
11
:
339
.
Fredholm
,
B.B.
,
E.
Irenius
,
B.
Kull
, and
G.
Schulte
.
2001
.
Comparison of the potency of adenosine as an agonist at human adenosine receptors expressed in Chinese hamster ovary cells
.
Biochem. Pharmacol
.
61
:
443
448
.
Friedl
,
P.
, and
D.
Gilmour
.
2009
.
Collective cell migration in morphogenesis, regeneration and cancer
.
Nat. Rev. Mol. Cell Biol
.
10
:
445
457
.
Friedl
,
P.
, and
R.
Mayor
.
2017
.
Tuning Collective Cell Migration by Cell-Cell Junction Regulation
.
Cold Spring Harb. Perspect. Biol
.
9
. a029199.
Gaggioli
,
C.
,
S.
Hooper
,
C.
Hidalgo-Carcedo
,
R.
Grosse
,
J.F.
Marshall
,
K.
Harrington
, and
E.
Sahai
.
2007
.
Fibroblast-led collective invasion of carcinoma cells with differing roles for RhoGTPases in leading and following cells
.
Nat. Cell Biol
.
9
:
1392
1400
.
Gao
,
Z.
,
B.S.
Li
,
Y.J.
Day
, and
J.
Linden
.
2001
.
A3 adenosine receptor activation triggers phosphorylation of protein kinase B and protects rat basophilic leukemia 2H3 mast cells from apoptosis
.
Mol. Pharmacol
.
59
:
76
82
.
Gendaszewska-Darmach
,
E.
,
M.
Maszewska
,
M.
Zakłos
, and
M.
Koziołkiewicz
.
2003
.
Degradation of extracellular nucleotides and their analogs in HeLa and HUVEC cell cultures
.
Acta Biochim. Pol
.
50
:
973
984
.
Goldberg
,
G.S.
,
V.
Valiunas
, and
P.R.
Brink
.
2004
.
Selective permeability of gap junction channels
.
Biochim. Biophys. Acta
.
1662
:
96
101
.
Goodenough
,
D.A.
, and
D.L.
Paul
.
2009
.
Gap junctions
.
Cold Spring Harb. Perspect. Biol
.
1
. a002576.
Györffy
,
B.
,
A.
Lanczky
,
A.C.
Eklund
,
C.
Denkert
,
J.
Budczies
,
Q.
Li
, and
Z.
Szallasi
.
2010
.
An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients
.
Breast Cancer Res. Treat
.
123
:
725
731
.
Han
,
Mitchell K.L.
, and
Johan
de Rooij
.
2016
.
Converging and Unique Mechanisms of Mechanotransduction at Adhesion Sites
.
Trends in Cell Biology
.
Hartsock
,
A.
, and
W.J.
Nelson
.
2008
.
Adherens and tight junctions: structure, function and connections to the actin cytoskeleton
.
Biochim. Biophys. Acta
.
1778
:
660
669
.
Haynes
,
S.E.
,
G.
Hollopeter
,
G.
Yang
,
D.
Kurpius
,
M.E.
Dailey
,
W.B.
Gan
, and
D.
Julius
.
2006
.
The P2Y12 receptor regulates microglial activation by extracellular nucleotides
.
Nat. Neurosci
.
9
:
1512
1519
.
Howe
,
A.K.
.
2004
.
Regulation of actin-based cell migration by cAMP/PKA
.
Biochim. Biophys. Acta
.
1692
:
159
174
.
Huang
,
G.Y.
,
E.S.
Cooper
,
K.
Waldo
,
M.L.
Kirby
,
N.B.
Gilula
, and
C.W.
Lo
.
1998
.
Gap junction-mediated cell-cell communication modulates mouse neural crest migration
.
J. Cell Biol
.
143
:
1725
1734
.
Ilina
,
O.
,
L.
Campanello
,
P.G.
Gritsenko
,
M.
Vullings
,
C.
Wang
,
P.
Bult
,
W.
Losert
, and
P.
Friedl
.
2018
.
Intravital microscopy of collective invasion plasticity in breast cancer
.
Dis. Model. Mech
.
11
. dmm034330.
Iwasa
,
S.
,
T.
Yanagawa
,
J.
Fan
, and
R.
Katoh
.
2009
.
Expression of CXCR4 and its ligand SDF-1 in intestinal-type gastric cancer is associated with lymph node and liver metastasis
.
Anticancer Res
.
29
:
4751
4758
.
Jacques-Silva
,
M.C.
,
R.
Rodnight
,
G.
Lenz
,
Z.
Liao
,
Q.
Kong
,
M.
Tran
,
Y.
Kang
,
F.A.
Gonzalez
,
G.A.
Weisman
, and
J.T.
Neary
.
2004
.
P2X7 receptors stimulate AKT phosphorylation in astrocytes
.
Br. J. Pharmacol
.
141
:
1106
1117
.
Kaczmarek
,
E.
,
L.
Erb
,
K.
Koziak
,
R.
Jarzyna
,
M.R.
Wink
,
O.
Guckelberger
,
J.K.
Blusztajn
,
V.
Trinkaus-Randall
,
G.A.
Weisman
, and
S.C.
Robson
.
2005
.
Modulation of endothelial cell migration by extracellular nucleotides: involvement of focal adhesion kinase and phosphatidylinositol 3-kinase-mediated pathways
.
Thromb. Haemost
.
93
:
735
742
.
Kanczuga-Koda
,
L.
,
S.
Sulkowski
,
A.
Lenczewski
,
M.
Koda
,
A.
Wincewicz
,
M.
Baltaziak
, and
M.
Sulkowska
.
2006
.
Increased expression of connexins 26 and 43 in lymph node metastases of breast cancer
.
J. Clin. Pathol
.
59
:
429
433
.
Kazan
,
J.M.
,
J.
El-Saghir
,
J.
Saliba
,
A.
Shaito
,
N.
Jalaleddine
,
L.
El-Hajjar
,
S.
Al-Ghadban
,
L.
Yehia
,
K.
Zibara
, and
M.
El-Sabban
.
2019
.
Cx43 Expression Correlates with Breast Cancer Metastasis in MDA-MB-231 Cells In Vitro, In a Mouse Xenograft Model and in Human Breast Cancer Tissues
.
Cancers (Basel)
.
11
:
460
.
Kenny
,
L.C.
,
P.N.
Baker
,
D.A.
Kendall
,
M.D.
Randall
, and
W.R.
Dunn
.
2002
.
The role of gap junctions in mediating endothelium-dependent responses to bradykinin in myometrial small arteries isolated from pregnant women
.
Br. J. Pharmacol
.
136
:
1085
1088
.
Khalil
,
A.A.
, and
J.
de Rooij
.
2019
.
Cadherin mechanotransduction in leader-follower cell specification during collective migration
.
Exp. Cell Res
.
376
:
86
91
.
Khalil
,
A.A.
,
O.
Ilina
,
P.G.
Gritsenko
,
P.
Bult
,
P.N.
Span
, and
P.
Friedl
.
2017
.
Collective invasion in ductal and lobular breast cancer associates with distant metastasis
.
Clin. Exp. Metastasis
.
34
:
421
429
.
Khoo
,
H.E.
,
C.L.
Ho
,
V.J.
Chhatwal
,
S.T.
Chan
,
S.S.
Ngoi
, and
S.M.
Moochhala
.
1996
.
Differential expression of adenosine A1 receptors in colorectal cancer and related mucosa
.
Cancer Lett
.
106
:
17
21
.
Kotini
,
M.
,
E.H.
Barriga
,
J.
Leslie
,
M.
Gentzel
,
V.
Rauschenberger
,
A.
Schambony
, and
R.
Mayor
.
2018
.
Gap junction protein Connexin-43 is a direct transcriptional regulator of N-cadherin in vivo
.
Nat. Commun
.
9
:
3846
.
Kumai
,
M.
,
K.
Nishii
,
K.
Nakamura
,
N.
Takeda
,
M.
Suzuki
, and
Y.
Shibata
.
2000
.
Loss of connexin45 causes a cushion defect in early cardiogenesis
.
Development
.
127
:
3501
3512
.
Kumar
,
V.
.
2013
.
Adenosine as an endogenous immunoregulator in cancer pathogenesis: where to go?
Purinergic Signal
.
9
:
145
165
.
Kuzma-Kuzniarska
,
M.
,
C.
Yapp
,
T.W.
Pearson-Jones
,
A.K.
Jones
, and
P.A.
Hulley
.
2014
.
Functional assessment of gap junctions in monolayer and three-dimensional cultures of human tendon cells using fluorescence recovery after photobleaching
.
J. Biomed. Opt
.
19
:
15001
.
Lin
,
Z.
,
P.
Yin
,
S.
Reierstad
,
M.
O’Halloran
,
V.J.
Coon
,
E.K.
Pearson
,
G.M.
Mutlu
, and
S.E.
Bulun
.
2010
.
Adenosine A1 receptor, a target and regulator of estrogen receptoralpha action, mediates the proliferative effects of estradiol in breast cancer
.
Oncogene
.
29
:
1114
1122
.
Marins
,
M.
,
A.L.
Xavier
,
N.B.
Viana
,
F.S.
Fortes
,
M.M.
Fróes
, and
J.R.
Menezes
.
2009
.
Gap junctions are involved in cell migration in the early postnatal subventricular zone
.
Dev. Neurobiol
.
69
:
715
730
.
Mascaux
,
C.
,
M.
Angelova
,
A.
Vasaturo
,
J.
Beane
,
K.
Hijazi
,
G.
Anthoine
,
B.
Buttard
,
F.
Rothe
,
K.
Willard-Gallo
,
A.
Haller
, et al
.
2019
.
Immune evasion before tumour invasion in early lung squamous carcinogenesis
.
Nature
.
571
:
570
575
.
Mayor
,
R.
, and
C.
Carmona-Fontaine
.
2010
.
Keeping in touch with contact inhibition of locomotion
.
Trends Cell Biol
.
20
:
319
328
.
McLachlan
,
E.
,
Q.
Shao
,
H.L.
Wang
,
S.
Langlois
, and
D.W.
Laird
.
2006
.
Connexins act as tumor suppressors in three-dimensional mammary cell organoids by regulating differentiation and angiogenesis
.
Cancer Res
.
66
:
9886
9894
.
Mirza
,
A.
,
A.
Basso
,
S.
Black
,
M.
Malkowski
,
L.
Kwee
,
J.A.
Pachter
,
J.E.
Lachowicz
,
Y.
Wang
, and
S.
Liu
.
2005
.
RNA interference targeting of A1 receptor-overexpressing breast carcinoma cells leads to diminished rates of cell proliferation and induction of apoptosis
.
Cancer Biol. Ther
.
4
:
1355
1360
.
Mroue
,
R.
,
J.
Inman
,
J.
Mott
,
I.
Budunova
, and
M.J.
Bissell
.
2015
.
Asymmetric expression of connexins between luminal epithelial- and myoepithelial- cells is essential for contractile function of the mammary gland
.
Dev. Biol
.
399
:
15
26
.
Nagelkerke
,
A.
,
S.J.
van Kuijk
,
F.C.
Sweep
,
I.D.
Nagtegaal
,
N.
Hoogerbrugge
,
J.W.
Martens
,
M.A.
Timmermans
,
H.W.
van Laarhoven
,
J.
Bussink
, and
P.N.
Span
.
2011
.
Constitutive expression of γ-H2AX has prognostic relevance in triple negative breast cancer
.
Radiother. Oncol
.
101
:
39
45
.
Naus
,
C.C.
, and
D.W.
Laird
.
2010
.
Implications and challenges of connexin connections to cancer
.
Nat. Rev. Cancer
.
10
:
435
441
.
Nemir
,
M.
,
D.
Bhattacharyya
,
X.
Li
,
K.
Singh
,
A.B.
Mukherjee
, and
B.B.
Mukherjee
.
2000
.
Targeted inhibition of osteopontin expression in the mammary gland causes abnormal morphogenesis and lactation deficiency
.
J. Biol. Chem
.
275
:
969
976
.
Ogawa
,
K.
,
P.
Pitchakarn
,
S.
Suzuki
,
T.
Chewonarin
,
M.
Tang
,
S.
Takahashi
,
A.
Naiki-Ito
,
S.
Sato
,
S.
Takahashi
,
M.
Asamoto
, et al
.
2012
.
Silencing of connexin 43 suppresses invasion, migration and lung metastasis of rat hepatocellular carcinoma cells
.
Cancer Sci
.
103
:
860
867
.
Othman
,
T.
,
H.
Yan
, and
S.A.
Rivkees
.
2003
.
Oligodendrocytes express functional A1 adenosine receptors that stimulate cellular migration
.
Glia
.
44
:
166
172
.
Owens
,
Smith
,
Hackett
, et al
.
1974
.
Epithelial cell cultures from normal glandular tissue of mice
.
J Natl Cancer Inst
.
Padmanaban
,
V.
,
I.
Krol
,
Y.
Suhail
,
B.M.
Szczerba
,
N.
Aceto
,
J.S.
Bader
, and
A.J.
Ewald
.
2019
.
E-cadherin is required for metastasis in multiple models of breast cancer
.
Nature
.
573
:
439
444
.
Park
,
S.W.
,
S.W.
Chen
,
M.
Kim
,
K.M.
Brown
,
V.D.
D’Agati
, and
H.T.
Lee
.
2010
.
Protection against acute kidney injury via A(1) adenosine receptor-mediated Akt activation reduces liver injury after liver ischemia and reperfusion in mice
.
J. Pharmacol. Exp. Ther
.
333
:
736
747
.
Plutoni
,
C.
,
E.
Bazellieres
,
M.
Le Borgne-Rochet
,
F.
Comunale
,
A.
Brugues
,
M.
Séveno
,
D.
Planchon
,
S.
Thuault
,
N.
Morin
,
S.
Bodin
, et al
.
2016
.
P-cadherin promotes collective cell migration via a Cdc42-mediated increase in mechanical forces
.
J. Cell Biol
.
212
:
199
217
.
Rebola
,
N.
,
P.C.
Pinheiro
,
C.R.
Oliveira
,
J.O.
Malva
, and
R.A.
Cunha
.
2003
.
Subcellular localization of adenosine A(1) receptors in nerve terminals and synapses of the rat hippocampus
.
Brain Res
.
987
:
49
58
.
Retamal
,
M.A.
,
N.
Froger
,
N.
Palacios-Prado
,
P.
Ezan
,
P.J.
Sáez
,
J.C.
Sáez
, and
C.
Giaume
.
2007
.
Cx43 hemichannels and gap junction channels in astrocytes are regulated oppositely by proinflammatory cytokines released from activated microglia
.
J. Neurosci
.
27
:
13781
13792
.
Riding
,
A.
, and
C.E.
Pullar
.
2015
.
ATP release and P2Y Receptor Signalling are Essential for Keratinocyte Galvanotaxis
.
J. Cell. Physiol
.
231
:
181
191
.
Saunders
,
M.M.
,
M.J.
Seraj
,
Z.
Li
,
Z.
Zhou
,
C.R.
Winter
,
D.R.
Welch
, and
H.J.
Donahue
.
2001
.
Breast cancer metastatic potential correlates with a breakdown in homospecific and heterospecific gap junctional intercellular communication
.
Cancer Res
.
61
:
1765
1767
.
Schalper
,
K.A.
,
D.
Carvajal-Hausdorf
, and
M.P.
Oyarzo
.
2014
.
Possible role of hemichannels in cancer
.
Front. Physiol
.
5
:
237
.
Sheng
,
W.
,
H.
Dong
,
D.Y.
Lee
,
W.Y.
Lu
, and
B.B.
Yang
.
2007
.
Versican modulates gap junction intercellular communication
.
J. Cell. Physiol
.
211
:
213
219
.
Sinyuk
,
M.
,
E.E.
Mulkearns-Hubert
,
O.
Reizes
, and
J.
Lathia
.
2018
.
Cancer Connectors: Connexins, Gap Junctions, and Communication
.
Front. Oncol
.
8
:
646
.
Stagg
,
J.
,
U.
Divisekera
,
N.
McLaughlin
,
J.
Sharkey
,
S.
Pommey
,
D.
Denoyer
,
K.M.
Dwyer
, and
M.J.
Smyth
.
2010
.
Anti-CD73 antibody therapy inhibits breast tumor growth and metastasis
.
Proc. Natl. Acad. Sci. USA
.
107
:
1547
1552
.
Stoletov
,
K.
,
J.
Strnadel
,
E.
Zardouzian
,
M.
Momiyama
,
F.D.
Park
,
J.A.
Kelber
,
D.P.
Pizzo
,
R.
Hoffman
,
S.R.
VandenBerg
, and
R.L.
Klemke
.
2013
.
Role of connexins in metastatic breast cancer and melanoma brain colonization
.
J. Cell Sci
.
126
:
904
913
.
Tambe
,
D.T.
,
C.C.
Hardin
,
T.E.
Angelini
,
K.
Rajendran
,
C.Y.
Park
,
X.
Serra-Picamal
,
E.H.
Zhou
,
M.H.
Zaman
,
J.P.
Butler
,
D.A.
Weitz
, et al
.
2011
.
Collective cell guidance by cooperative intercellular forces
.
Nat. Mater
.
10
:
469
475
.
Theveneau
,
E.
, and
R.
Mayor
.
2012
.
Cadherins in collective cell migration of mesenchymal cells
.
Curr. Opin. Cell Biol
.
24
:
677
684
.
Töpfer
,
M.
,
C.E.
Burbiel
,
C.E.
Müller
,
J.
Knittel
, and
E.J.
Verspohl
.
2008
.
Modulation of insulin release by adenosine A1 receptor agonists and antagonists in INS-1 cells: the possible contribution of 86Rb+ efflux and 45Ca2+ uptake
.
Cell Biochem. Funct
.
26
:
833
843
.
Tsuji
,
K.
,
K.
Yamauchi
,
M.
Yang
,
P.
Jiang
,
M.
Bouvet
,
H.
Endo
,
Y.
Kanai
,
K.
Yamashita
,
A.R.
Moossa
, and
R.M.
Hoffman
.
2006
.
Dual-color imaging of nuclear-cytoplasmic dynamics, viability, and proliferation of cancer cells in the portal vein area
.
Cancer Res
.
66
:
303
306
.
Umapathy
,
S.N.
,
E.
Kaczmarek
,
N.
Fatteh
,
N.
Burns
,
R.
Lucas
,
K.R.
Stenmark
,
A.D.
Verin
, and
E.V.
Gerasimovskaya
.
2013
.
Adenosine A1 receptors promote vasa vasorum endothelial cell barrier integrity via Gi and Akt-dependent actin cytoskeleton remodeling
.
PLoS One
.
8
. e59733.
Vaupel
,
P.
, and
A.
Mayer
.
2016
.
Hypoxia-Driven Adenosine Accumulation: A Crucial Microenvironmental Factor Promoting Tumor Progression
.
Adv. Exp. Med. Biol
.
876
:
177
183
.
Vaupel
,
P.
, and
G.
Multhoff
.
2016
.
Adenosine can thwart antitumor immune responses elicited by radiotherapy: Therapeutic strategies alleviating protumor ADO activities
.
Strahlenther. Onkol
.
192
:
279
287
.
Venhuizen
,
J.H.
,
S.
Sommer
,
P.N.
Span
,
P.
Friedl
, and
M.M.
Zegers
.
2019
.
Differential expression of p120-catenin 1 and 3 isoforms in epithelial tissues
.
Sci. Rep
.
9
:
90
.
Wang
,
X.
,
W.
Qin
,
X.
Xu
,
Y.
Xiong
,
Y.
Zhang
,
H.
Zhang
, and
B.
Sun
.
2017
.
Endotoxin-induced autocrine ATP signaling inhibits neutrophil chemotaxis through enhancing myosin light chain phosphorylation
.
Proc. Natl. Acad. Sci. USA
.
114
:
4483
4488
.
Wen
,
J.
,
A.
Grenz
,
Y.
Zhang
,
Y.
Dai
,
R.E.
Kellems
,
M.R.
Blackburn
,
H.K.
Eltzschig
, and
Y.
Xia
.
2011
.
A2B adenosine receptor contributes to penile erection via PI3K/AKT signaling cascade-mediated eNOS activation
.
FASEB J
.
25
:
2823
2830
.
Westcott
,
J.M.
,
A.M.
Prechtl
,
E.A.
Maine
,
T.T.
Dang
,
M.A.
Esparza
,
H.
Sun
,
Y.
Zhou
,
Y.
Xie
, and
G.W.
Pearson
.
2015
.
An epigenetically distinct breast cancer cell subpopulation promotes collective invasion
.
J. Clin. Invest
.
125
:
1927
1943
.
Wolf
,
K.
,
M.
Te Lindert
,
M.
Krause
,
S.
Alexander
,
J.
Te Riet
,
A.L.
Willis
,
R.M.
Hoffman
,
C.G.
Figdor
,
S.J.
Weiss
, and
P.
Friedl
.
2013
.
Physical limits of cell migration: control by ECM space and nuclear deformation and tuning by proteolysis and traction force
.
J. Cell Biol
.
201
:
1069
1084
.
Yamaguchi
,
N.
,
T.
Mizutani
,
K.
Kawabata
, and
H.
Haga
.
2015
.
Leader cells regulate collective cell migration via Rac activation in the downstream signaling of integrin β1 and PI3K
.
Sci. Rep
.
5
:
7656
.
Ye
,
Z.C.
,
M.S.
Wyeth
,
S.
Baltan-Tekkok
, and
B.R.
Ransom
.
2003
.
Functional hemichannels in astrocytes: a novel mechanism of glutamate release
.
J. Neurosci
.
23
:
3588
3596
.
Zhang
,
Y.W.
,
K.
Nakayama
,
K.
Nakayama
, and
I.
Morita
.
2003
.
A novel route for connexin 43 to inhibit cell proliferation: negative regulation of S-phase kinase-associated protein (Skp 2)
.
Cancer Res
.
63
:
1623
1630
.
Zhang
,
J.
,
K.F.
Goliwas
,
W.
Wang
,
P.V.
Taufalele
,
F.
Bordeleau
, and
C.A.
Reinhart-King
.
2019
.
Energetic regulation of coordinated leader-follower dynamics during collective invasion of breast cancer cells
.
Proc. Natl. Acad. Sci. USA
.
116
:
7867
7872
.
Zhou
,
J.Z.
,
M.A.
Riquelme
,
X.
Gao
,
L.G.
Ellies
,
L.Z.
Sun
, and
J.X.
Jiang
.
2015
.
Differential impact of adenosine nucleotides released by osteocytes on breast cancer growth and bone metastasis
.
Oncogene
.
34
:
1831
1842
.

Author notes

A.A. Khalil’s present address is Department of Molecular Cancer Research, Center for Molecular Medicine, University Medical Center Utrecht, Utrecht, Netherlands.

O. Ilina’s present address is Department of Tumor Immunology, Radboud University Medical Center, Nijmegen, Netherlands.

A. Vasaturo’s present address is Ultivue, Segrate, Milan, Italy.

This article is distributed under the terms of an Attribution–Noncommercial–Share Alike–No Mirror Sites license for the first six months after the publication date (see http://www.rupress.org/terms/). After six months it is available under a Creative Commons License (Attribution–Noncommercial–Share Alike 4.0 International license, as described at https://creativecommons.org/licenses/by-nc-sa/4.0/).

Supplementary data