Main

Effective clearance of apoptotic cells is important in maintaining tissue homeostasis. Weaning-induced mammary gland involution is a unique model for studying apoptotic cell clearance, as 80–90% of the milk-producing mammary epithelium undergoes apoptosis to return the gland to a non-secretory state.1 Professional phagocytes, such as macrophages, are recruited into the involuting mammary gland; however, they are thought to have a limited role in the clearance of dying secretory cells, as in mice, peak macrophage infiltration occurs after the majority of apoptotic cell removal.2 Rather, the neighboring mammary epithelial cells themselves appear to be the primary cell type responsible for apoptotic cell clearance during involution.2 Rapid and efficient apoptotic cell clearance is essential, as persistence of apoptotic cells can result in the release of cell fragments into the local environment and subsequent autoimmunity.3 Importantly, impaired apoptotic cell clearance in the postpartum mammary gland results in local inflammation, fibrosis and epithelial cell hyperplasia.4, 5

Although there is increasing evidence that phagocytosis by mammary epithelial cells has a crucial role in maintaining tissue homeostasis in the involuting murine mammary gland, little is known about how mammary epithelial cells become phagocytic during postpartum involution. One of the key changes in the mammary epithelium that may contribute to acquisition of a phagocytic phenotype is reorganization of epithelial cell junctions. During lactation, tight junctions between mammary epithelial cells become highly impermeable, which assures localization of milk within the mammary ducts.6 With weaning, this impermeability is lost,6 consistent with tight junction reorganization. Furthermore, reorganization of adherens junctions is also observed upon the switch from lactation to involution.7 Given that professional phagocytes such as macrophages do not exist in monolayers with cell cell junctions, disruption of epithelial cell junctions at the onset of mammary gland involution may be required for mammary epithelial cells to become phagocytic.

One candidate cytokine for promoting epithelial cell junction reorganization and phagocytosis is transforming growth factor-β (TGF-β). Binding of TGF-β to the TGF-β type II receptor (TβRII) activates canonical signaling through a signaling cascade involving the TGF-β type I receptor, receptor-associated Smads (Smad2/3) and Smad4. TGF-β protein and mRNA levels are significantly increased in the mammary gland on the switch from lactation to involution, with increased expression persisting through at least 9 days post weaning.8, 9 Of the three TGF-β isoforms (TGF-β1, -β2 and -β3), TGF-β3 increases the greatest upon the lactation-to-involution switch.8, 9, 10, 11 Overexpressing TGF-β3 or depleting Smad3 or TβRII in the mammary epithelium reveals a necessary role for TGF-β in promoting apoptosis early during involution.10, 12, 13, 14 However, sustained TGF-β expression throughout the postpartum involution window suggests additional roles for TGF-β that extend beyond apoptosis induction, including influencing extracellular matrix remodeling and immune cell composition.8, 10, 12, 13, 14, 15 TGF-β is known to increase the phagocytic capacity of retinal pigment epithelial cells, fibroblasts and macrophages,16, 17, 18 although a role for TGF-β in mediating apoptotic cell clearance by phagocytic mammary epithelial cells has not been explored. Furthermore, TGF-β is implicated in tight junction disruption in the mammary gland and has known roles in adherens junction disassembly, making it an intriguing target to investigate in the promotion of a phagocytic phenotype in mammary epithelial cells.6, 19

Currently, it is unknown whether the mammary epithelium has a role in apoptotic cell clearance in species other than mice. Therefore, we evaluated rat and human involution mammary tissue for apoptotic cell clearance by the mammary epithelium. Further, as addressing the role of TGF-β in promoting phagocytosis by mammary epithelial cells during gland involution is challenging due to impaired cell death in the absence of TGF-β signaling,12, 13, 14 we developed an in vitro model to investigate the role of TGF-β3 in mammary epithelial cell junction reorganization and phagocytosis.

We demonstrate engulfment of apoptotic cells by mammary epithelial cells during weaning-induced involution in both rats and women, supportive of phagocytosis being a conserved feature of mammary epithelium during postpartum involution. Using our murine mammary epithelial culture model that mimics the high junctional resistance of the lactating gland, we show that TGF-β3 promotes phagocytic capability and identify a potential role for cell–cell junction disruption in epithelial cell phagocytosis. Furthermore, we identify a previously unreported role for the intramembrane protease γ-secretase in the promotion of phagocytosis by TGF-β3. In light of the negative consequences of impaired apoptotic cell clearance during postpartum involution,4, 5 our data provide insight into how mammary epithelial cells may contribute to apoptotic cell clearance during this time.

Results

In mice, mammary epithelial cells have been identified as the primary phagocyte responsible for apoptotic cell clearance in the postpartum involuting mammary gland.2, 4 To determine whether apoptotic cell clearance is a conserved function of mammary epithelium between species, dual fluorescent staining for apoptotic cell markers and the epithelium was performed using postpartum involuting mammary tissue from rats and women. In both species, increased apoptotic cell clearance by the alveolar mammary epithelium was observed during involution (Figures 1a–c), indicating conservation of phagocytosis by mammary epithelial cells in rat and human. Furthermore, during involution, human mammary epithelial cells upregulated the expression of the mannose receptor (Figure 1d), a scavenger receptor involved in engulfment.24 Mannose receptor expression by mammary epithelial cells during involution was transient and was detected in 100% of cases within 1 month of lactation cessation with gradual decline to 0% detection by 12 months postpartum (Figure 1e). This unique pattern of mannose receptor expression and the presence of apoptotic cells within the mammary epithelium are the first demonstration of a ‘clearance’ phenotype in human breast epithelium.

Figure 1
figure 1

Phagocytosis is a conserved function of mammary epithelium in rat and human. (a) Apoptotic cells within mammary epithelial cells in rat (top) and human (bottom) tissue as demonstrated by dual staining for the epithelial cell marker E-cad (green) and the apoptotic cell markers TUNEL or cleaved caspase-3 (CC3) (red); arrows, apoptotic cells within epithelial cells; scale bar=10 μm. (b) Quantification of TUNEL-positive nuclei within live mammary epithelial cells per 100 × field in tissue isolated from lactating (Lac) or involution (Inv) day 4 rats; ***P<0.001, two-tailed, unpaired t-test. (c) Quantification of TUNEL-positive cells within live mammary epithelial cells in lactational and involution lobules in breast tissue from women; **P<0.01, two-tailed, unpaired t-test (N=5 cases). (d) Mannose receptor (MR) expression (brown) by mammary epithelial cells in human breast tissue from involuting (Inv), but not nulliparous (Nullip) women; scale bar=50 μm. (e) Percentage of human cases with mannose receptor-positive mammary epithelial cells; Nullip, nulliparous; Preg, pregnant; MECs, mammary epithelial cells; N, number of cases/reproductive category

How mammary epithelial cells transition from a predominantly milk-secreting cell during lactation to function as a phagocyte during postpartum involution is unexplored. Based on reports of TGF-β enhancing phagocytic ability in both professional and non-professional phagocytes,16, 17, 18 we hypothesize that TGF-β promotes epithelial cell phagocytosis during involution. Supportive of this hypothesis, we observed increased TGF-β protein in postpartum involuting rat mammary tissue (Figure 2a), in agreement with previously reported murine data.8, 9, 10, 11 TGF-β was between 150 and 250 kDa, consistent with TGF-β being in the latent complex. However, during involution, strong nuclear phospho-Smad2 (p-Smad2) was observed, demonstrating increased active TGF-β during involution. Further, nuclear p-Smad2 was found within the mammary epithelium, but not within the stroma (Figure 2b), identifying mammary epithelial cells as highly TGF-β responsive.

Figure 2
figure 2

Increased TGF-β signaling in the rat mammary gland during postpartum involution. (a) Immunoblot analysis (immunoblot, left panel; quantification, right panel) for TGF-β using pooled rat mammary tissue lysate (N=6 rats/stage) demonstrates increased TGF-β between 150 and 250 kDa, consistent with latent TGF-β, during postpartum involution; Lac, lactation; Inv2–10, involution days 2–10. (b) Nuclear p-Smad2 (IHC staining, left panel; quantification, right panel), as a measure of TGF-β signaling, increases in mammary epithelial cells during postpartum involution in the rat; dashed line=representative alveoli, alveolar lumens pseudocolored white to remove nonspecific staining of milk; scale bar=20 μm, **P<0.01 versus lactation, Dunnett’s multiple comparison’s test (~14 500 total nuclei quantified)

To test whether TGF-β promoted phagocytosis by mammary epithelial cells, an in vitro mammary epithelium model was developed. Non-transformed murine mammary epithelial EpH4 cells were cultured at high density on transwell filters to form a polarized monolayer with functional tight and adherens junctions, demonstrated by zona occludens-1 (ZO-1) and E-cadherin (E-cad) localization at cell–cell contacts, respectively (Figure 3a). In this model, the monolayer acquires high transepithelial electrical resistance (TEER), which is a functional readout of junction impermeability25 (Figure 3b). This mammary epithelium model mimics aspects of lactational epithelium in which tight junctions are well organized and highly impermeable.26 Following tight monolayer formation, EpH4 cells were treated for 24 h with TGF-β3, which was chosen over other TGF-β isoforms based on its increased abundance during early involution compared with TGF-β1 (Figure 3c and previous studies8, 9, 10, 11) and our observation that TGF-β2 associates with lactation, but not involution (Supplementary Figure S2). After TGF-β3 treatment, EpH4 cell monolayers were incubated with an enriched population of apoptotic cells, consisting of ~64% apoptotic, 18% live and 18% necrotic cells (Supplementary Figure S3A), and positive for the apoptosis marker cleaved caspase-3 (Supplementary Figure S3B). Following TGF-β3 treatment, an ~4-fold increase in apoptotic EpH4 cells was observed intracellularly within the live mammary epithelial cell monolayer (Figures 3d and e). Labeled cells located within the live epithelial cell monolayer had fragmented nuclei or lacked detectable DNA. We did not find evidence for epithelial cell engulfment of live cells, as phagocytized cells with intact nuclei were not observed. These data are consistent with increased engulfment of apoptotic EpH4 cells by the epithelial cell monolayer following TGF-β3 treatment. TGF-β3 treatment also increased engulfment of carboxylated latex beads, a surrogate for apoptotic cells,27, 28 which accumulate within the cytoplasm of the epithelial cells due to lack of degradation27, 28 (Figures 3f and g). Further, increased phagocytic ability following TGF-β3 treatment was dose responsive (Figure 3h). Importantly, as TGF-β can induce apoptosis in non-transformed mammary epithelial cells,10, 29 we confirmed that induction of apoptosis was not observed in EpH4 cell monolayers following 1 ng/ml TGF-β3 treatment (Supplementary Figure S4). For subsequent phagocytosis assays, we used the 1 ng/ml TGF-β3 dose.

Figure 3
figure 3

TGF-β3 treatment leads to increased phagocytosis and cell–cell junction reorganization in mammary epithelial cell monolayers. (a) Localization of E-cad (green) and ZO-1 (red) to cell–cell junctions in EpH4 cell monolayers before TGF-β3 treatment; scale bar=10 μm. (b) Average TEER in EpH4 cell monolayers at 5, 6 and 7 days post plating before treatment with TGF-β3 (N=10 experiments). (c) Analysis of publically available data sets demonstrates increased Tgbf3 gene expression before Tgfb1 in mouse mammary tissue. (d) TGF-β3 (1 ng/ml) treatment induces phagocytosis of apoptotic cells by EpH4 cells; **P<0.01, two-tailed t-test. (e) Engulfment of apoptotic cells labeled with CellTracker red (red) by EpH4 cells (E-cad, green) following 1 ng/ml TGF-β3 treatment; dashed line, indicates level of z-stack; scale bar=10 μm. (f) TGF-β3 (1 ng/ml) induces phagocytosis of carboxylated latex beads by EpH4 cells (N=13 experiments); *P<0.05, one-sample, one-tailed t-test compared with vehicle set to 1. (g) Engulfment of carboxylated latex beads (green spheres) by EpH4 cells following 1 ng/ml TGF-β3 treatment; scale bar=50 μm. (h) Enhanced phagocytosis of carboxylated latex beads with increasing doses of TGF-β3; *P<0.05, **P<0.01, ***P<0.001, one-tailed t-test compared with vehicle (N=3 wells per condition). (i) EpH4 cell monolayers fail to increase resistance following TGF-β3 treatment; *P<0.005, two-tailed t-test (N=7 experiments). (j) TGF-β3 induces relocalization of ZO-1 (red) away from cell–cell junctions; scale bar=10 μm. (k) Increased junction disruption following TGF-β3 treatment (reported as the difference in TEER between vehicle- and TGF-β3-treated cells (ΔTEER) associates with increased phagocytosis. (l) Increased E-cad cytoplasmic localization (green; scale bar=10 μm) and cleavage (m) following TGF-β3 treatment

Following TGF-β3 treatment, EpH4 cell monolayers did not increase TEER (Figure 3i) and had focal disruption of the tight-junction protein ZO-1 at the cell–cell interface (Figure 3j). Loss of tight junction integrity correlated with increased phagocytosis (Figure 3k), suggesting a role for junction disruption in increasing the ability to phagocytose. To further investigate the role of cell junction disruption in epithelial cell phagocytosis, we evaluated adherens junctions. Following TGF-β3 treatment, EpH4 monolayers displayed decreased E-cad junctional staining and increased cytoplasmic localization, as well as increased E-cad cleavage (Figures 3l and m). In addition, gene expression analysis revealed downregulation of multiple claudins, classical cadherns and cell adhesion molecules with TGF-β3 treatment (Table 1), demonstrating global changes in cell junctions with acquisition of the phagocytic phenotype. In transformed epithelium, TGF-β is known to induce cell junction disruption through activation of pathways involved in epithelial-to-mesenchymal transition (EMT);30 therefore, expression of mesenchymal markers and transcription factors involved in EMT were evaluated. Notably, Snai1, Snai2, Twist1, Twist2, Zeb1, Vim and Fn1, genes highly associated with EMT, were not upregulated in EpH4 cell monolayers following TGF-β3 treatment and, in fact, Zeb1 was significantly downregulated (Supplemental Table S3). These data indicate that the junctional changes in our phagocytic model are probably independent of TGF-β-induced EMT.

Table 1 Genes associated with cell junctions that are downregulated following TGF-β3 treatment

To determine whether E-cad reorganization occurs in vivo at the time of the phagocytic switch, as occurs in our in vitro model, E-cad localization was assessed in rat and human mammary tissue. In the nulliparous rat, mammary epithelial cells exhibited weak E-cad cytoplasmic and apical staining, and moderate to strong junctional staining (Figure 4a). During pregnancy, staining was primarily junctional, but at a lower intensity. Localization of E-cad to the lateral cell junctions was greatest during lactation (Figure 4a), consistent with the known junctional impermeability during this time.26 Following lactation, during early involution, E-cad relocalized to the apical cell surface (Figure 4a, Inv2-4, arrows) and in some cells was found to be strongly cytoplasmic (Figure 4a, Inv2-6, asterisk). E-cad junctional disruption was transient, as lateral cell junction localization was restored to nulliparous levels in the fully regressed gland (Figure 4a, Reg). In human breast tissue, E-cad was similarly localized primarily to the lateral cell junctions during lactation (Figure 4b, Lac panel), with apical (Figure 4b, Inv panel, arrow) and cytoplasmic localization during involution (Figure 4b, Inv panel, asterisk). Cumulatively, these data indicate that relocalization of E-cad away from lateral mammary epithelial cell junctions during involution is conserved in rats and humans. Further evidence of E-cad-based junction remodeling during involution was indicated by partial E-cad proteolysis specific to the involution time points (Figure 4c). E-cad cleavage during involution is reported to result in the release of the β-catenin-binding domain.7 To determine whether E-cad cleavage associated with relocalization of β-catenin away from the lateral epithelial cell junctions, β-catenin localization was evaluated in rat mammary tissue. Nuclear accumulation of β-catenin within the mammary epithelial cells was detected as early as involution day 4 and peaked at involution day 6 (Figure 4d). These data are consistent with E-cad cleavage resulting in the release of β-catenin from the epithelial cell membranes during involution.

Figure 4
figure 4

Mammary epithelial cell adherens junctions reorganize during postpartum involution. (a) Representative IHC images of rat mammary tissue E-cad show increased apical (arrow) and cytoplasmic (asterisk) localization during postpartum involution; Vir, virgin; Preg, pregnant; Lac, lactation; Inv2–10, involution days 2–10; Reg, regressed; scale bar=10 μm. (b) Increased apical (arrow) and cytoplasmic (asterisk) E-cad observed in human involution breast tissue; scale bar=20 μm. (c) Increased fragmentation of E-cad in rat mammary tissue during postpartum involution (lysates pooled from six rats/stage); FL, full length; CTF1–3, C-terminal fragment 1–3. (d, left panel) IHC images of β-catenin (β-cat) show increased cytoplasmic (asterisk) and nuclear (arrow head) localization during postpartum involution in the rat mammary gland; scale bar=20 μm. (d, right panel) Quantification of β-catenin IHC images; MEC, mammary epithelial cells, #P<0.05, ##P<0.0001; one-tailed t-test compared with lactation (~31 000 nuclei quantified)

The release of the β-catenin-binding domain from E-cad has been shown to occur through cleavage of E-cad by the intramembrane protease γ-secretase;7 therefore, we investigated a role for γ-secretase in the rat involution model as well as in our in vitro phagocytic mammary epithelial cell model. In the rat mammary gland, γ-secretase activity increased from lactation to involution days 4 and 6 (Figure 5a). The γ-secretase enzyme complex consists of four known subunits—nicastrin, anterior pharynx-defective 1, presenilin (PS) enhancer 2 and PS. PS is the catalytic component of the complex, whereas the other three proteins serve as obligate co-factors.31 Analysis of mouse microarray data sets32, 33 demonstrated increased PS-1 gene (Psen1) expression in the mammary gland during involution (Figure 5b), consistent with increased γ-secretase activity during this time. In vitro, TGF-β3 modestly, but consistently, increased Psen1 mRNA expression by EpH4 cell monolayers (Figure 5c), as well as PS-1 dimers and the C-terminal fragment (Figure 5d), which represent stable pools and active PS-1 protein, respectively.34, 35

Figure 5
figure 5

Evidence for increased γ-secretase activity during involution in vivo and following TGF-β3 treatment in vitro. (a) Increased γ-secretase activity in rat mammary tissue lysate during postpartum involution. Data are reported as fluorescent units normalized to GAPDH, as described in Methods; *P<0.05, ****P<0.001 Dunnett’s multiple comparison test. (b) Increased Psen1 gene expression in the mouse mammary gland during early stages of postpartum involution, using two publically available data sets. Data set 1 error bars=S.D., N=3, P<0.05 versus lactation, Dunnett’s multiple comparison’s test. (c) TGF-β3 induces a modest increase in Psen1 expression in EpH4 cell monolayers; *P<0.05, one-sample, two-tailed t-test compared with 1 (N=3 experiments). (d) Increased PS-1 dimers and C-terminal fragment (CTF) in EpH4 cell monolayers following TGF-β3 treatment (Immunoblot, left panel; quantification, right panel), *P<0.05, two-tailed t-test (N=3 experiments)

We next sought to determine whether γ-secretase activity and β-catenin signaling were required for induction of phagocytosis in the mammary epithelial cells following TGF-β3 treatment. This question was addressed using pharmacologic inhibitors. Inhibition of γ-secretase activity with DAPT36 or L-685,45837 decreased phagocytosis following TGF-β3 treatment (Figures 6a and b), indicating that γ-secretase activity is required for increased phagocytic ability in our model. As γ-secretase cleavage of E-cad can lead to the release of β-catenin from the cell membrane, we next sought to determine whether increased phagocytosis following TGF-β3 treatment was β-catenin-dependent. Blocking β-catenin with the inhibitor iCRT338 reduced phagocytosis downstream of TGF-β3 by ~60% (Figure 6c), indicating a role for β-catenin in the increased ability of EpH4 cells to phagocytose following TGF-β3 treatment. To address whether β-catenin signaling is sufficient to increase phagocytic ability, EpH4 cell monolayers were treated with LiCl at doses that induced β-catenin target genes Axin2 and Tcf7, but not Ccnd1 (Cyclin D1) (Figure 6d), which we speculate is due to EpH4 cells in a mature monolayer being non-proliferative. Induction of β-catenin signaling with LiCl treatment did not enhance phagocytic ability in this model (Figure 6e). These data suggest that although β-catenin signaling is needed for increased phagocytic ability, it is not sufficient. Based on these data, we propose a model in which increased phagocytic ability in the mammary epithelial cell monolayers following TGF-B3 treatment requires both junctional reorganization following γ-secretase cleavage of E-cad as well as β-catenin signaling (Figure 6f).

Figure 6
figure 6

γ-Secretase and β-catenin mediate TGF-β3-induced phagocytosis in EpH4 cells. Inhibition of γ-secretase with DAPT (a) or L-685,458 (b) reduces TGF-β3-induced engulfment of carboxylated beads; two-tailed t-test, *P<0.05 versus TGF-β3+DMSO, n.s.=non-significant. (c) Inhibition of β-catenin with the inhibitor iCRT3 reduces TGF-β3-induced engulfment of carboxylated beads, *P<0.05, one-sample, two-tailed t-test. (N=3 experiments). LiCl treatment induces expression of β-catenin target genes (d), but is insufficient to induce phagocytosis (e). N=3 experiments. (f) Proposed model of phagocytosis induction in mammary epithelial cells. In the absence of TGF-β signaling (left panel), cell–cell tight and adherens junction are intact (left panel). In the presence of TGF-β signaling, there is decreased integrity of cell–cell junctions and γ-secretase cleavage of E-cad, both of which are necessary for induction of phagocytosis. Following E-cad cleavage, β-catenin is freed from the cell membrane, allowing it to move to the nucleus where it activates transcription of genes involved in the acquisition of a phagocytic phenotype and may contribute to downregulation of cell junction-associated genes

Discussion

Deficiencies in apoptotic cell clearance during postpartum mammary gland involution have substantial negative consequences on tissue architecture and function.4, 5, 39 Upon the cessation of lactation, mammary epithelial cells lose their secretory phenotype, and—somewhat surprisingly—acquire phagocytic properties to function as the primary cell responsible for the clearance of dying secretory epithelial cells during postpartum involution.2 Currently, little is known about how mammary epithelial cells obtain phagocytic capacity to clear the gland of the immense number of apoptotic cells generated during involution.2, 4 Here we describe a novel in vitro model for studying the induction of a phagocytic phenotype in mammary epithelial cells that permits the maturation of a polarized epithelium with impermeable tight junctions similar to those found in the mammary gland during lactation.6, 26, 40 TGF-β3 treatment in this model increased phagocytic ability and recapitulates multiple aspects of the lactation to involution switch including TGF-β signaling in the mammary epithelium, loss of tight and adherens junction integrity in an epithelial cell monolayer with previously established cell junctions, increased γ-secretase activity and increased phagocytosis by the epithelium. Our data unveil a previously undescribed role for TGF-β3 in the enhancement of phagocytosis by mammary epithelial cells and implicate a mechanism dependent on γ-secretase and β-catenin. Furthermore, using clinical samples, we demonstrate that human breast epithelial cells function in apoptotic cell clearance during involution. Our in vivo finding that human breast epithelial cells express the mannose receptor specifically during the involution window provides further evidence that the mammary epithelium has a unique, phagocytic phenotype during this time. Mannose receptor is a scavenger receptor involved in the engulfment of glycoproteins, as well as opsinization-independent recognition and engulfment of bacterial, viral and fungal pathogens.24 Mannose receptor expression by the mammary epithelial cells raises the question of whether their phagocytic capacity extends beyond that of apoptotic cell engulfment to include clearance of milk constituents and host defense.24

Based on our in vivo and in vitro data, we predict that upon the switch from lactation to involution, TGF-β3 induces γ-secretase activity, leading to E-cad cleavage, release of β-catenin from the cell membrane, translocation of β-catenin into the nucleus and increased transcription of genes involved in phagocytosis (Figure 6f). In support of this model and consistent with our results, TGF-β has been shown by others to increase the expression levels and activity of PS-1,41 and that PS-1 can cleave E-cadheren and promote β-catenin signaling.42 Our demonstration that two different γ-secretase inhibitors reduced TGF-β3-induced phagocytosis suggests that loss of membranous β-catenin at the lactation to involution switch is dependent on γ-secretase enzyme activity. Our studies using the β-catenin activator LiCl suggest that induction of β-catenin signaling alone may not be sufficient to induce phagocytosis in mammary epithelial cells. However, a limitation to this study is that LiCl functions through inhibition of glycogen synthase kinase-3β, which has a diverse array of proposed substrates in addition to β-catenin,43 and the effect of these other LiCl targets on inhibiting phagocytosis cannot be discounted. Furthermore, in addition to regulation of the β-catenin pathway, E-cad also influences the phosphatidylinositol-4,5-bisphosphate 3-kinase/Akt, Rho GTPase and nuclear factor κB pathways, all of which have been shown to have a role in phagocytosis in other models.44, 45, 46, 47 Genetic ablation and constitutive activation approaches would directly test the roles of E-cad cleavage and β-catenin signaling in the induction of the phagocytic phenotype. However, transient transfections of mature, electrically resistant EpH4 cell monolayers were unsuccessful and permanent genetic manipulation of these adhesion proteins can prevent the formation of mammary epithelial monolayers with functional junctions.49, 50 For example, E-cad cleavage by γ-secretase requires binding of PS-1 to residues 604–615 of the mature processed E-cad.37, 48 Thus, although mutants lacking this region of E-cad cannot be processed by γ-secretase, loss of these residues also prevents proper formation of adherens junctions, precluding the use of such a mutant in the model described here.37, 49 Likewise, expression of the cytoplasmic domain of E-cad functions in a dominant-negative manner to disrupt the normal structure of both adherens and tight junctions, indicating that this approach is also not able to be used to investigate the potential role of the cytoplasmic E-cad fragment generated by γ-secretase cleavage in phagocytosis.50, 51

The work described here raises intriguing questions regarding the role of phagocytic mammary epithelial cells in vivo. The correlation between decreased junctional integrity and increased phagocytic capacity indicates that the loss of cell junctions may be an important step in mammary epithelial cells becoming phagocytic upon the switch from lactation to involution; however, cell junctions are tighter during lactation than any other stage, begging the question of whether or not mammary epithelial cells are normally phagocytic, except for during lactation. Given that the mammary epithelium undergoes small cycles of proliferation and apoptosis in response to ovarian hormone cycling,52 it is possible that mammary epithelial cells also function in apoptotic cell clearance to maintain tissue homeostasis in the non-pregnant, premenopausal state. However, mannose receptor expression by mammary epithelial cells in human breast tissue exclusively during early involution suggests that even if mammary epithelial cells have a basal phagocytic capacity, they obtain a unique phagocytic phenotype during postpartum involution. Murine models using up-the-teat injections of carboxylated latex beads, apoptotic cells and other targets of engulfment will provide further insight into the extent of the phagocytic capacity of the mammary epithelium across reproductive stages.

Further understanding of how mammary epithelial cells become phagocytic is of interest as deficiencies in apoptotic cell clearance during postpartum involution lead to increased macrophage infiltration, collagen deposition, vascularization and epithelial proliferation.4, 5, 39 These stromal changes and hyperproliferative epithelium are similar to the changes observed with ductal carincoma in situ of the breast, leading to the speculation that defective apoptotic cell clearance may provide a catalyst for breast cancer development.53, 54, 55, 56, 57 In addition, the immunosuppressive immune milieu associated with clearance of apoptotic cells by the mammary epithelium has also been proposed to contribute to breast cancer progression.58, 59 Women have a transient increase in breast cancer risk during the postpartum period and breast cancers diagnosed postpartum have poor prognosis.60, 61 In preclinical postpartum breast cancer models, postpartum tumors are associated with increased immune-suppressive cytokines, macrophages and regulatory T cells.62, 63 Apoptotic cell clearance by the macrophages can promote production of the immune-suppressive Th2 cytokines IL-4, IL-10 and TGF-β1,63, 64, 65, 66 supporting the hypothesis that apoptotic cell clearance during postpartum involution promotes an immune-suppressive microenvironment that allows tumor cells to escape immune surveillance and metastasize.58, 59 Consistent with this hypothesis, Standford et al.63 demonstrated decreased immune-suppressive macrophages and a dramatic reduction in metastasis following inhibition of the apoptotic cell receptor MerTK in postpartum mouse mammary tumors. However, the effects of MerTK inhibition on the primary tumor cells were not addressed in this model and cannot be discounted, as MerTK upregulation is observed in metastatic breast cancers67 and MerTK can mediate cancer cell invasion in vitro.68, 69 The potential roles for apoptotic cell clearance in breast cancer promotion and the increased risk for poor-prognosis breast cancers during the postpartum period highlight the need for in vitro models, such as the one described here, which can be used to further understand mechanisms involved in apoptotic cell clearance by the mammary epithelium, subsequent consequences on the local microenvironment and breast cancer development.

Materials and Methods

Human tissue acquisition

Under protocols approved by the University of Colorado Institutional Review Board, 53 breast tissue specimens were obtained from premenopausal women aged 21–45 years, who underwent clinically indicated biopsy (n=2) or surgery (n=51). Histologically normal tissue was identified by pathologic review and analyzed an average of >8 mm away from the cancer. Based on patient histories, cases were stratified into the following reproductive categories: nulliparous (n=7), pregnant (n=3), ≤1 month post-lactation (n=6), >1 month and ≤6 months postpartum (n=6), >6 months and ≤12 months postpartum (n=4), >12 months and ≤18 months postpartum (n=2), >18 months and ≤24 months postpartum (n=7), >24 months and ≤3 years postpartum (n=4), >3 years and ≤6 years postpartum (n=8), and >6 years postpartum (n=6).

Animal husbandry, reproductive staging and tissue preparation

Animal procedures were approved by the University of Colorado Anschutz Medical Campus Institutional Care and Use Committee. Sprague–Dawley rats (Harlan Laboratories, Indianapolis, IN, USA) 70 days of age were bred, normalized to six pups per dam at 2 days post-parturition and mammary tissue collected at the following reproductive stages as described:20 age-matched virgin; pregnant (days 18–20); lactation (days 10–11); involution day 2, 4, 6, 8 and 10 (2–10 days post-wean); and fully regressed (4 weeks post-wean).

Immunohistochemical and TUNEL staining, image acquisition and quantification of breast and mammary tissues

For human and rodent tissue, 4-μm-thick, formalin-fixed paraffin-embedded sections were deparaffinized and rehydrated. For immunohistochemical (IHC) staining, except p-Smad2, sections were subjected to antigen retrieval (125 °C under pressure for 5 min) and sequentially incubated in 3% H2O2 in methanol (room temperature (RT), 5 min), Protein Block (RT, 10 min; Dako, Carpinteria, CA, USA), primary antibody (RT, 1 h), secondary antibody (RT, 30 min), 3,3’-diaminobenzidine (RT, 10 min; Dako) and hematoxylin (RT, 6 min; Dako). See Supplementary Table S1 for specific antigen retrieval and antibody information. p-Smad2 staining was performed as described.21 For dual E-cad and TUNEL (terminal deoxynucleotidyl transferase dUTP nick end labeling) staining of rat tissue, E-cad IHC was performed followed by TUNEL staining using the TACS 2-TdT-Fluor In Situ Apoptosis Detection Kit (Trevigen, Gaithersburg, MD, USA) with the Mn2+ cation at a final concentration of 1 × as per the manufacturer’s instructions. Fluorescent images were acquired using an Olympus IX81 inverted microscope and Slidebook v.4.067 software (Intelligent Imaging Innovations, Denver, CO, USA). For TUNEL staining of human tissue, the TACS 2-TdT-DAB In Situ Apoptotic Detection Kit (Trevigen) was used with the Co2+ cation at a final concentration of 1 × as per the manufacturer’s instructions. For rat tissue, quantification of apoptotic cell clearance by mammary epithelial cells was performed by counting the number of TUNEL+ nuclei within the mammary epithelium in ten 100 × fields per rat (N=7 rats per stage). For human tissue, the number of TUNEL+ nuclei within the mammary epithelium was counted in 7–10 lactation lobules and 7–10 involution lobules within the same case (N=5 cases) in breast tissue of women who were within 6 weeks post-lactation. For each case, the data are reported relative to lactation. Images of mannose receptor-stained tissue were achieved by extracting regions from whole-slide images that were acquired using a ScanScope T3 scanner (Aperio Technologies, Vista, CA, USA) at 0.46 μm per pixel. Owing to heterogeneous mannose receptor expression in the alveolar epithelium across a single case, quantification of mannose receptor expression by mammary epithelial cells was performed by determining whether or not each case had any mannose receptor-positive mammary alveoli. In cases without mannose receptor staining in the mammary epithelial cells, evidence that the staining worked was provided by identifying positive stromal cells. For E-cad, β-catenin and p-Smad2, light microscopy images were acquired using a Zeiss Axioskop (Zeiss, Thornwood, NY, USA) and SPOT RT3 camera (SPOT Imaging Solutions, Sterling Heights, MI, USA). Quantification of nuclear β-catenin was performed by counting the total number of epithelial nuclei and determining the percentage that were β-catenin positive. Mammary epithelium from the lymph node region of the fourth right mammary gland was analyzed from four rats per stage, with five representative fields analyzed per rat.

In vitro phagocytosis assay

Non-transformed murine mammary EpH4 cells were cultured in DMEM with high glucose (Sigma-Aldrich, St. Louis, MO, USA), without L-glutamine or sodium pyruvate, with additives at the following final concentrations: 5% FBS, 2 mM L-glutamine, 10 mM HEPES and 100 units/ml each of penicillin and streptomycin (all additives from Hyclone, Logan, UT, USA). In order to promote junctional integrity before use in phagocytosis assays, cells were maintained as confluent monolayers split 1:3 every 5–7 days with daily media changes. For phagocytosis assays, 70 000 EpH4 cells were plated in complete medium onto 0.4 μm filters (BD Biosciences, San Jose, CA, USA) in a 24-well format with the media changed daily. Cells were cultured for ~7 days, to allow for the formation of a monolayer with intact cell junctions. To monitor integrity of cell junctions, the TEER was measured daily, before changing the media. Once the TEER reached 300 Ω*cm2, complete medium was removed, cells were washed twice with phosphate buffered saline (PBS) and then treated with 1, 3 or 5 ng/ml recombinant human (rh) TGF-β3, 20 or 40 mM LiCl, or vehicle in serum-free medium for 24 h. For inhibition studies, cells were treated with 500 nM L-685,458, 20 μM DAPT or 25 μM iCRT3 simultaneously with 1 ng/ml TGF-β3. TGF-β3 and DAPT were purchased from R&D Systems (Minneapolis, MN, USA), and LiCl, L-685,458 and iCRT3 were purchased from Sigma-Aldrich. After 24 h, medium containing the treatment was removed and cells were washed twice with PBS. EpH4 monolayers were then incubated with the phagocytic target (either 10 μm Fluoresbrite BB carboxylated latex microspheres/beads (Polysciences, Inc., Warrington, PA, USA) or apoptotic cells) for 4–6 h, followed by washing three times with ice-cold PBS to remove bound, but unengulfed, beads or cells. EpH4 cell monolayers contain similar cell numbers across treatment conditions (data not shown); therefore, phagocytic index was quantified by multiplying the number of phagocytes by the average number of targets engulfed per phagocyte, with four fields analyzed per filter and three filters per condition per experiment. Results are reported as fold change compared with vehicle treatment due to variation in baseline phagocytosis levels across experiments.

Generation of apoptotic cells

EpH4 cells at 95% confluence were labeled with 1 μM CellTracker Red CMPTX (Molecular Probes, Life Technologies, Grand Island, NY, USA) at 37 °C for 30 min. Following labeling, apoptosis was induced in EpH4 cells by subculturing 1:2 into apoptosis induction media (DMEM high-glucose base media with no L-glutamine or sodium pyruvate, plus 1% 1 M HEPES and 1% pen/strep) at 37 °C for 16 h. After 16 h, the percent live, apoptotic and necrotic cells in the preparation was quantified by morphological criteria and trypan blue exclusion. Live cells were rounded with no evidence of apoptotic morphology and trypan blue negative, apoptotic cells were those with blebbing or apoptotic nuclei and trypan blue negative, and necrotic cells were those that were trypan blue positive.

Immunocytochemistry of mammary epithelial cells in culture

EpH4 cell monolayers on transwell filters were fixed in 4% paraformaldehyde for 15 min and permeabilized using 0.5% Triton X-100 (Sigma-Aldrich). Fixed cells on filters were then incubated sequentially at RT with 1.1 M glycine for 5 min to quench autofluorescence, 20 mg/ml bovine serum albumin for 1 h to block nonspecific protein interactions and with primary and secondary antibodies (Supplemental Table S2) for 60 and 30 min, respectively. Filters were then removed from transwells with a scalpel and mounted onto slides using Prolong Gold Antifade Reagent with DAPI (Life Technologies).

γ-Secretase activity assay

For measurement of γ-secretase activity in mammary gland tissue lysates, rat mammary tissue with the lymph node removed was pulverized in liquid nitrogen and 50 mg per rat was homogenized in 400 μl of assay buffer (50 mM Tris-HCl pH 6.8, 2 mM EDTA acid, 0.25% CHAPSO). Protein concentrations were determined using Bio-Rad Protein Assay Dye Reagent (Hercules, CA, USA). Ten micrograms of total protein and 8 μM fluorogenic γ-secretase substrate (Merck Millipore, Billerica, MA, USA) were incubated in assay buffer in a final reaction volume of 150 μl for 5 h at 37 °C. After incubation, reactions were centrifuged at 16 100 × g for 15 min at 4 °C. Supernatant was transferred to a black-walled 96-well plate (Thermo Scientific, Waltham, MA, USA) and fluorescence was measured at excitation wavelength of 360/40 nm and emission wavelength of 460/40 nm using a BioTek Synergy 2 (Winooski, VT, USA). To account for differences in cellular content in the mammary tissue from lactation to involution day 6, fluorescence was normalized to glyceraldehyde 3-phosphate dehydrogenase (GAPDH) level in homogenized mammary tissue samples as determined by immunoblotting (Supplementary Figure S1). For immunoblotting, 15 μg of mammary tissue lysate homogenized in γ-secretase activity assay buffer was separated by SDS-polyacrylamide gel electrophoresis (SDS-PAGE) using a 10% acrylamide gel. Immunoblotting using rabbit-anti-GAPDH (Sigma-Aldrich; 1 h, RT) and densitometry were performed as described below.

Immunoblotting and quantification

Pooled rat mammary tissue lysates were prepared as described (n=6 animals/reproductive group).22 EpH4 cell lysates were prepared by collecting cells directly into 100 μl radioimmunoprecipitation assay buffer (10 μM Tris pH 7.4, 150 μM NaCl, 0.1% SDS, 1% sodium deoxycholate and 1% Triton X-100) containing protease and phosphatase inhibitors (1 × Protease Inhibitor Cocktail (Sigma-Aldrich), 200 μg/ml phenylmethylsulfonyl fluoride, 1 mM Na3VO4 and 50 mM NaF). For TGF-β, samples were run under non-reducing conditions by preparing mammary tissue lysate samples in loading buffer containing 10% SDS, 50% glycerol, 0.4 M Tris pH 6.8 and 0.1% w/v bromophenol blue without reducing agents. Thirty microliters of total protein were separated by SDS-PAGE using a 10% acrylamide gel. Membranes were blocked in 5% non-fat dry milk in Tris-buffered saline (10 mM Tris-HCl pH 8.0, 150 mM NaCl) containing 0.1% Tween 20 for 1 h at RT. Immunoblotting was performed using mouse anti-bovine TGF-β1, 2 and 3 (MAB1835, R&D Systems; overnight, 4 °C), mouse anti-human E-cad (610181, BD Biosciences; overnight, 4 °C), rabbit anti-human PS-1 (ab65293, Abcam, Cambridge, MA, USA) or rabbit anti-mouse GAPDH (G9545, Sigma-Aldrich; 1 h, RT), followed by anti-mouse or anti-rabbit horseradish peroxidase-conjugated secondary antibody (Santa Cruz Biotechnologies, Santa Cruz, CA, USA, and Bio-Rad Laboratories, respectively; 30 min, RT) with detection using Enhanced Chemiluminescence Western Blotting Substrate (Thermo Scientific). Densitometry was performed using ImageJ Software (http://imagej.nih.gov/ij/).

RNA isolation, cDNA synthesis and quantitative RT-PCR

Two hundred microliters of TRIzol Reagent (Life Technologies) were added to EpH4 cell monolayers per 24-well transwell filter immediately following removal of medium and samples were collected using a cell scraper modified to fit within the transwell filter. RNA was isolated as per the manufacturer’s protocol and, to facilitate RNA precipitation, 10 μg UltraPure glycogen (Life Technologies) was added per sample during the RNA precipitation. Resultant RNA was resuspended in 20 μl nuclease-free water and RNA concentration determined using a NanoDrop 2000 UV-Vis Spectrophotometer (Thermo Scientific). Reverse-transcriptase reactions were performed using the iScript cDNA Synthesis Kit (Bio-Rad Laboratories) with 1 μg total RNA. For quantitative RT-PCR, 1 μl of reverse transcriptase reactions was combined with SYBR green (Bio-Rad Laboratories) and forward and reverse primers at a concentration of 1 μM each, and brought up to a final reaction volume of 20 μl with nuclease-free water. The following primers were used, from Integrated DNA Technologies (Coralville, Iowa, USA): Psen1, forward: 5′-ACC TGC ACC TTT GTC CTA CTT CCA-3′, reverse: 5′-TGT CAA GTC TCT GCC TGT CAT GCT-3′; Actb, forward: 5′-GCA ACG AGC GGT TCC G-3′, reverse: 5′-CCC AAG AAG GAA GGC TGG A-3′; Ccnd1, forward: 5′-TAC CAT GGA GGG TGG GTT GGA AAT-3′, reverse: 5′-TGC TGC AAA TGG AAC TGC TTC TGG-3′; axin2, forward: 5′-ACG CAC TGA CCG ACG ATT C-3′, reverse: 5′-CCA TGC GGT AAG GAG GGA C-3′; and Tcf7, forward: 5′-AAG GTC ATT GCT GAG TGC ACA C-3′, reverse: 5′-TGC ATG CCA CCT GCG AC-3′. Reactions were run using the following program: 3 min, 95 °C; 10 s, 95 °C; 1 min, 60 °C (42 cycles); 1 min, 95 °C; 1 min, 55 °C. Primer specificity was confirmed via melt-curve analysis. Data were collected using the MyiQ Single Color Real-Time PCR Detection System with iQ5 software (Bio-Rad Laboratories) and the Psen1 gene expression was normalized to actin using the ΔΔCt method.

Microarray gene expression

Affymetrix Mouse Gene 1.1 ST Array Strips (Affymetrix, Santa Clara, CA, USA) were used to measure gene expression using the GeneAtlas system (Affymetrix). RNA was isolated from EpH4 cells using TRIzol Reagent (Life Technologies), as described above, and 250 ng per sample was prepared for microarray analysis and hybridized using the Ambion WT Expression Kit (Life Technologies) and GeneChip WT Terminal Labeling Kit (Affymetrix), as per the manufacturer’s protocol. Raw gene expression data was extracted and normalized by Robust Multiarray Average algorithm23 using Affymetrix Power Tools program. Microarray data has been deposited at the NCBI Gene Expression Omnibus under accession number GES62126.

Statistical analysis

Statistical analyses were performed using GraphPad Prism v.6 software (GraphPad Software, Inc., La Jolla, CA, USA), as indicated.