Monoallelic Becn1 loss results in accelerated mammary fat pad filling and excessive ductal side-branching

In agreement with earlier results using a different Becn1^+/− transgenic line,⁴ we also detected increased proliferation in mammary ducts and terminal end buds (TEBs) from 5-wk-old Becn1^+/− mice⁵ (Fig. S1A and S1B); apoptosis rates were similar (Fig. S1C and S1D). Consistent with the partially defective autophagic response observed in Becn1^+/− iMMECs under metabolic stress,⁶ lower levels of MAP1LC3B/LC3B (microtubule-associated protein 1 light chain 3 β) and accumulation of SQSTM1/p62 (sequestosome 1) were observed in MGs from Becn1^+/− mice (Fig. S1E and S1F).

Systematic examination of MGs from cohoused Becn1^+/+ and Becn1^+/− littermates revealed that 6.5-wk-old Becn1^+/− mice exhibited accelerated mammary fat-pad filling and increased ductal side-branching (Fig. 2A and B). The latter phenotype persisted and was further accentuated in MGs from 6-mo-old nulliparous Becn1^+/− mice (Fig. 2C, upper panel); at that age, Becn1^+/− mice also exhibited mammary hyperplasias (Fig. 2C, lower panel). The difference in ductal side-branching in MGs from Becn1^+/+ and Becn1^+/− mice became less pronounced by 12 mo of age (Fig. 2D), possibly indicating that increased ductal side-branching may be a hormonally regulated phenotype.

Open in a separate window

Figure 2.

Monoallelic Becn1 loss results in accelerated mammary fat pad filling and excessive ductal side-branching. (A–D) Analysis of mouse mammary epithelium from Becn1^+/+ (left) and Becn1^+/− (right) mice (n = 3 to 5): (A) MG whole mounts from 6.5-wk-old mice, with circle identifying an unfilled (left) vs. filled (right) fat pad and (B) (top) quantification of number of side-branches per primary branch per field and (bottom) quantification of percentage of mammary fat pad 4 filled with epithelium. Three MG specimen per genotype from cohoused Becn1^+/+ and Becn1^+/− littermates were evaluated. Results are presented as means ± SDs. (C) MG whole mounts (top) and hematoxylin and eosin staining of MG sections (bottom) from 6-mo-old mice, with mammary hyperplasia in separate insert. (D) MG whole mounts from 12-mo-old mice. Scale bar: (A, C,and D) 5 mm for whole mounts; (A, C, and D) 2 mm for whole mount enlargements. Scale bar: (C) 100 μm for H&Es; 200 μm for hyperplasia. *P < 0.05 by a 2-tailed Student t test.

Monoallelic Becn1 loss results in functionally enriched mammary stem and progenitor cell populations

Faster mammary fat pad filling and excessive ductal side-branching result from increased mammary stem and progenitor cell populations/activity, respectively.²² We utilized fluorescence-activated cell sorting (FACS) to isolate PI⁻ LIN⁻ CD24⁺ CD29^hi, i.e. mammary stem cell (MaSC)-containing, and PI⁻ LIN⁻ CD24⁺ CD29^lo, i.e. mammary progenitor cell (MaPC)-containing, populations²² from cohoused Becn1^+/+ and Becn1^+/− littermates. Although the relative proportions of MaSC-, MaPC-, and CD61 (cluster differentiation 61)⁺ luminal progenitor cell-containing populations^22-24 were similar between the 2 genotypes (Fig. S2A and S2B), MGs from Becn1^+/− mice consistently demonstrated higher MaSC and MaPC activities by in vitro and in vivo functional assays (Fig. 3; Fig. S2C and S2D). Specifically, CD24⁺ CD29^hi Becn1^+/− MECs exhibited increased MaSC activity, as measured by primary and limiting dilution mammosphere (Fig. 3A, left and right panel, respectively) and 3D colony-formation (Fig. 3B) assays.^23,25 Furthermore, CD24⁺ CD29^lo Becn1^+/− MECs exhibited increased MaPC activity (Fig. 3C), as measured by number of colonies formed in 2D-matrigel assays.^23,25

Open in a separate window

Figure 3.

Monoallelic Becn1 loss results in functionally enriched mammary stem and progenitor cell populations. (A–C) Freshly isolated Becn1^+/+ (blue) and Becn1^+/− (red) MECs were used, experiments were performed 3 independent times and outcomes are presented as means ± SDs. (A) Primary mammosphere assays performed using Pi⁻ LIN(CD31,CD45,LY76)⁻ CD24⁺ CD29^hi MECs plated at 20,000 cells per well (left) and in limiting dilution of 500 to 25 CD24⁺ CD29^hi MECs plated (right) reveal increased MaSC activity in Becn1^+/− MECs. (B) 3D colony-formation by CD24⁺ CD29^hi Becn1^+/− compared with Becn1^+/+ MECs plated at 20,000 cells per well shows increased MaSC activity. (C) Colony formation by CD24⁺ CD29^lo Becn1^+/− compared with Becn1^+/+ MECs plated at 1,000 cells per well reveals increased colony forming ability in 2D-matrigel conditions. (D) Increased repopulation frequency and mammary fat pad filling are seen following transplantation of 250 CD24⁺ CD29^hi Becn1^+/− MECs isolated from 5-wk-old mice. Becn1^+/+ (left) and Becn1^+/− (right), contralateral transplantations in wild-type recipient mouse. Circles represent fat pads and the black color represents the percentage of the mammary fat pad that is filled in. (E) Representative whole mount images of outgrowths from contralateral CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MEC transplantation are shown and reveal increased ductal side branching in Becn1^+/− samples. *P < 0.05 by a 2-tailed Student t test. Scale bar: (E) 5 mm for whole mounts; 2 mm for whole mount enlargements.

The in vivo stem/progenitor cell activity of CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MECs was evaluated by contralateral transplantation of these cells into cleared 4th mammary fat pads of wild-type mice.²² CD24⁺ CD29^hi Becn1^+/− MECs consistently resulted in higher number of successful (i.e. giving rise to mammary outgrowths) transplantations (Fig. 3D; Fig. S2C and S2D). Six weeks post-transplantation, comparison of similarly filled fat pads revealed that CD24⁺ CD29^hi Becn1^+/− MECs gave rise to outgrowths with more extensive side-branching than CD24⁺ CD29^hi Becn1^+/+ MECs (Fig. 3E). This phenotype was similar to that observed in MGs from Becn1^+/− mice (Fig. 2), suggesting that transplantation of MaSC-containing Becn1^+/− MECs in wild-type mice sufficiently reproduced the side-branching phenotype observed in native MGs from Becn1^+/− mice, where all cell types (MECs and stroma) are missing one Becn1 allele.

Since transplantations of CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MECs were performed in wild-type recipient mice, we examined the autophagy status of resulting outgrowths and found that similar to MGs from Becn1^+/− mice (Fig. S1E and S1F), CD24⁺ CD29^hi Becn1^+/− MEC-generated outgrowths exhibited LC3B expression decrease and SQSTM1 accumulation (Fig. S2E and S2F). Furthermore, similar to MGs from Becn1^+/− mice (Fig. S1A and S1B), proliferation rates in TEBs and mammary ducts were significantly higher in outgrowths derived from CD24⁺ CD29^hi Becn1^+/− MEC transplantation (Fig. 4A and B). In contrast to similar apoptosis rates observed in native glands (Fig. S1C and S1D), in CD24⁺ CD29^hi Becn1^+/− MEC-generated outgrowths, apoptosis was suppressed (Fig. 4C and D), possibly due to enrichment in more apoptosis-resistant MaSCs and MaPCs.

Open in a separate window

Figure 4.

Outgrowths from CD24⁺ CD29^hi Becn1^+/− MEC transplantation display increased proliferation and TNFSF-NFKB pathway activation. (A–D, and F) Examination of outgrowths from contralateral transplantations of CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MECs, and (E) native MGs from 6.5-wk-old cohoused Becn1^+/+ and Becn1^+/− littermates. (A) Representative images show more MKI67⁺ cells in ducts (top) and TEBs (bottom) of outgrowths from CD24⁺ CD29^hi Becn1^+/− MEC transplantation; (B) quantification of MKI67⁺ cells in mammary ducts and TEBs is shown as means ± SDs (n = 3 mammary gland specimens per genotype). (C) Representative images show fewer cleaved CASP3-staining cells in ducts (top) and TEBs (bottom) of outgrowths from CD24⁺ CD29^hi Becn1^+/− MEC transplantation; (D) quantification of cleaved CASP3⁺ cells in mammary ducts and TEBs is shown as means ± SDs (n = 3 mammary gland specimens per genotype). Increased TNFSF11 (red) expression is seen in (E) native MGs from 6.5-wk-old Becn1^+/− mice and (F) outgrowths from contralateral CD24⁺ CD29^hi Becn1^+/− MEC transplantation (top panel). Increased TNFRSF11A (middle panel) and RELA (bottom panel) levels are also seen in outgrowths from CD24⁺ CD29^hi Becn1^+/− MEC transplantation. Scale bar: (A, C, E, and F, middle and bottom panels) 30 μm ; (F, top panel) 40 μm.

MEC subpopulations exhibit different levels of autophagic activity

Given that monoallelic Becn1 loss resulted in deregulation of the mammary hierarchy, thus implicating BECN1, and possibly autophagy in general, in mammary differentiation, we isolated MEC populations from wild-type mice, as described above, and examined basal autophagy in these populations by electron microscopy (EM)¹⁷ (Fig. S3). The CD24⁺ CD29^hi population exhibited paucity in mitochondria, autophagosomes, and other cytoplasmic organelles (Fig. S3, left panel), whereas the CD24⁺ CD29^lo population showed high autophagic activity, as illustrated by the presence of numerous autophagosomes (Fig. S3, middle panel). Interestingly, the autophagosomal cargo was different between CD24⁺ CD29^hi and CD24⁺ CD29^lo MECs; in the former (MaSC-containing) population, autophagosomes appeared largely carrying mitochondria, while in the latter (progenitor cell-containing) population, diverse cargo was observed in autophagosomes. Lastly, CD24⁻ MECs contained fewer autophagosomes than the CD24⁺ CD29^lo population, but more than the CD24⁺ CD29^hi compartment, and more mitochondria than both other MEC populations examined (Fig. S3, right panel). High autophagy induction in CD24⁺ CD29^lo MECs, though not a proof, suggests that autophagy may play a significant role in mammary differentiation, which is also accompanied by enrichment in mitochondria.

Krt14-driven biallelic Atg7 deletion delays mammary fat pad filling, without altering ductal side-branching

In an effort to investigate whether the increased ductal side-branching and enriched MaSC and progenitor cell activities observed in MGs from Becn1^+/− mice were secondary to defective autophagy associated with monoallelic Becn1 loss or to a nonautophagy-related BECN1 function, we generated mice with conditional biallelic deletion of the essential autophagy regulator Atg7²⁶ utilizing Krt14 promoter-driven expression of CRE recombinase.²⁷ In agreement with prior studies showing that Atg7 is not happloinsufficient for autophagy deficiency or other phenotypes,²⁶ we found that MGs from Atg7^+/+ and Atg7^F/+;Krt14-Cre mice displayed similar morphological characteristics (data not shown).

In MGs from atg7^F/F;Krt14-Cre mice, Atg7 deletion and defective autophagy status were confirmed (Fig. S4A and S4B). In contrast to our observations in developing MGs from Becn1^+/− mice (Fig. 2A and B), whole MG mounts revealed no difference in ductal side-branching between Atg7^+/+ and atg7^F/F;Krt14-Cre mice and, more impressively, showed delayed—rather than accelerated—fat pad filling in 6.5-wk-old atg7^F/F;Krt14-Cre mice (Fig. S4C and S4D). While proliferation rates were similar, increased apoptosis was detected in MGs from atg7^F/F;Krt14-Cre compared with Atg7^+/+ mice (Fig. S4E and S4F). Furthermore, KRT6, KRT8, and KRT14 levels were comparable in MGs from atg7^F/F;Krt14-Cre and Atg7^+/+ mice (Fig. S4G and S4H). Although functional in vitro and in vivo studies for MaSC and progenitor cell activities in Atg7^+/+ and Atg7-null CD24⁺ CD29^hi and CD24⁺ CD29^lo MECs were not performed, decreased mammary fat pad filling in pubertal MGs from atg7^F/F;Krt14-Cre mice suggests, at least partial, impairment of MaSC or progenitor cell function in association with MEC-specific biallelic Atg7 deletion, and points to a nonautophagy-related role of BECN1 in mammary differentiation. Alternatively, and together with the evidence for high autophagic activity during mammary epithelial differentiation (Fig. S3) and the mammary phenotypes observed in MGs from Becn1^+/− mice (Figs. 1, ​,2,2, and ​and3),3), it raises the possibility of a dose-dependent requirement for functional autophagy in MaSC and progenitor cell maintenance. If this is the case, partial autophagy deficiency may affect mammary cell hierarchy differently than a more severe autophagy defect. Further investigation of these 2 not necessarily mutually exclusive, mechanistic possibilities will be needed.

Monoallelic Becn1 loss results in TNFRSF11A-NFKB pathway activation in mammary tissues

TNFSF11/TNF11/RANKL (tumor necrosis factor [ligand] superfamily, member 11) has been implicated in progesterone-induced MaSC expansion²⁸ and mice overexpressing TNFSF11 or its receptor, TNFRSF11A/RANK (tumor necrosis factor receptor superfamily, member 11a, NFKB activator), under the mouse mammary tumor virus (MMTV) promoter,^29,30 exhibit mammary phenotypes similar to Becn1^+/− mice, including excessive ductal side-branching (Fig. 2), increased proliferation (Fig. S1A and S1B), and mammary hyperplasias (Fig. 2C), but no mammary tumors in nulliparous mice.^4,15 We hypothesized that ectopic TNFRSF11A-TNFSF11 signaling may be mediating the increased mammary stem and progenitor cell activities in MGs from Becn1^+/− mice. TNFSF11 was indeed detected in MGs from pubertal Becn1^+/− mice (Fig. 4E) and robustly in outgrowths from CD24⁺ CD29^hi Becn1^+/− MEC transplantations (Fig. 4F). Consistent with aberrant TNFRSF11A pathway activation, outgrowths from CD24⁺ CD29^hi Becn1^+/− MEC transplantation also displayed increased TNFRSF11A expression (Fig. 4F) and likely activation of the downstream effector NFKB axis,³¹ as indicated by increased expression, though not apparent nuclear localization, of RELA (Fig. 4F). In agreement with previous studies,^28-30,32-34 this data suggests that aberrant TNFRSF11A pathway activation may play a role in the MaSC and progenitor cell expansion observed in association with monoallelic Becn1 loss. However, further studies will need to be conducted to determine whether TNFRSF11A signaling is responsible for the mammary phenotypes observed in Becn1^+/− mice.

Locally increased progesterone signaling in MGs from Becn1^+/− mice

Given that progesterone is required for mammary ductal side-branching³⁵ and TNFSF11 is a mediator of both progesterone signaling^28,32 and progestin-induced mammary tumorigenesis,³⁴ we examined MGs from cohoused Becn1^+/+ and Becn1^+/− littermate mice for evidence of progesterone signaling. As shown (Fig. 5A and B), monoallelic Becn1 loss resulted in higher proportion of PGR (progesterone receptor)-positive MECs. Transplantation of CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MECs into contralateral fat pads of a wild-type recipient resulted in mammary outgrowths of both Becn1 genotypes in the same mouse; similar to native MGs, and despite uniform systemic hormone levels, the proportion of nuclear PGR-positive cells was again higher in outgrowths generated by Becn1^+/− MECs (Fig. 5C and D), indicating that monoallelic Becn1 loss results in locally-regulated increased progesterone signaling.

Open in a separate window

Figure 5.

Augmented progesterone signaling, cytokines, and oxidative stress detected upon monoallelic Becn1 loss. (A and C) Nuclear PGR (red) and KRT6 (green) cell staining in native MGs from cohoused Becn1^+/+ and Becn1^+/− MEC and outgrowths from contralateral CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MEC transplantation. (A) Representative images of PGR⁺ and KRT6⁺ cells and (B) quantification of nuclear PGR⁺ cells (left) and KRT6⁺ cells (right) in MGs from 5-wk old mice reveals increased numbers of PGR⁺ and KRT6⁺ cells in Becn1^+/− mice. (C) Representative images of nuclear PGR⁺ cells (top panel) and KRT6⁺ cells (bottom panel) in outgrowths from contralateral CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MEC transplantation with (D) quantification of PGR⁺ cells (top panel) and KRT6⁺ cells (bottom panel). Results are presented as means ± SDs (n = 3 mammary gland specimens per genotype). (E) Quantification of ductal side-branching in MGs from hormonally stimulated following ovariectomization, where Sham = ovariectomized, Vehicle = ovariectomized + 100 μL sesame oil daily, and E + P = ovariectomized + 5 μg 17B-estradiol + 5 mg progesterone in 100 μL of sesame oil daily. (F–H) Examination of outgrowths from contralateral transplantations of CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MECs. (F) Increased IL1B (top panel) and TNF (bottom panel) seen in outgrowths from Becn1^+/− MEC transplantation. (G) Similar F4/80 staining detected in outgrowths from Becn1^+/+ and Becn1^+/− MEC transplantation. (H) Increased 8-O-dG is seen in outgrowths from Becn1^+/− MEC transplantation. *P < 0.05, **P < 0.01, and ***P < 0.001 by a 2-tailed Student t test. Scale bar: (A, C, F, G, and H) 30 μm; (E) 50 μm.

To investigate the role of hormonal signaling in the mammary hyper side-branching phenotype observed in association with monoallelic Becn1 loss, Becn1^+/+ and Becn1^+/− mice were ovariectomized at 3 to 4 wk of age. In the absence of endogenous hormone production, MGs from Becn1^+/+ and Becn1^+/− mice displayed similar ductal side-branching (Fig. 5E), whereas hormonal stimulation of ovariectomized mice with exogenous 17β-estradiol and progesterone³² rescued the hyper ductal side-branching phenotype in Becn1^+/− vs. Becn1^+/+ mice (Fig. 5E). This data implies that hormonal signaling is required for the increased ductal side-branching seen in MGs from Becn1^+/− mice, and supports our observations that this phenotype is associated with progesterone axis activation (Fig. 2C and D) and is likely hormonally regulated.

In progesterone-driven TNFRSF11A signaling, exogenous progesterone activates PGR-positive MECs to secrete TNFSF11, which in turn induces proliferation of neighboring TNFRSF11A-expressing mammary stem and progenitor cells and subsequent generation of more PGR-positive MECs, thus creating a positive feedback loop.^28,32 In macrophages, TNFSF11 production and activation of TNFRSF11A-NFKB signaling is driven by proinflammatory cytokines.³⁶ In our studies, where contralaterally transplanted Becn1^+/+ and Becn1^+/− MECs are exposed to the same systemic hormones, outgrowths generated by CD24⁺ CD29^hi Becn1^+/− MECs displayed increased IL1B (interleukin 1, β) and TNF/TNFα (tumor necrosis factor) levels (Fig. 5F), which were not accompanied by enhanced macrophage infiltration (Fig. 5G). Together these studies indicate that IL1B and TNF were likely secreted by and acting autonomously on Becn1^+/− MECs.

Autophagy-deficient tissues commonly harbor a proinflammatory microenvironment,³⁷ mainly attributed to increased reactive oxygen species (ROS) due to impaired removal of damaged mitochondria.³⁷ Consistent with their functional autophagy defect (Fig. S2E and S2F), CD24⁺ CD29^hi Becn1^+/− MEC-generated outgrowths exhibited increased oxidative stress, as evaluated by 8-oxo-7,8-dihydroguanine (8-O-dG) staining (Fig. 5H).

Becn1^+/− mice develop spontaneous mammary tumors following parity

In earlier studies which reported mammary hyperplasias, but no mammary tumors, in Becn1^+/− mice independent of age and parity,^4,15 the mouse genetic background was C57BL/6, which is known to suppress mammary tumorigenesis.³⁸ To ensure that mammary tumor studies in Becn1^+/− mice were not confounded by the C57BL/6 status, we changed the genetic background of Becn1^+/− mice⁵ to FVB/N, which is more amenable to mammary tumorigenesis studies.³⁸

Nulliparous and multiparous (post 2 pregnancy/lactation/involution cycles) Becn1^+/+ and Becn1^+/− female mouse cohorts were generated, cohoused and monitored for mammary tumor formation. Similar to C57BL/6 mice,^4,15 both nulliparous and multiparous Becn1^+/− FVB/N mice developed lymphomas, lung and liver tumors, and mammary hyperplasias with age. Furthermore, mammary ducts from older (24-mo-old) nulliparous Becn1^+/− FVB/N mice showed increased cellular density and irregular epithelial cell layering (Fig. 6A). Similar, but more pronounced, changes of increased epithelial density and disorganized ductal architecture were observed in MGs from multiparous Becn1^+/− multiparous mice at a younger age (14 to 17 mo) (Fig. 6B). In contrast to C57BL/6 mice and nulliparous FVB/N mice, 4 out of 21 mice (19%) of multiparous Becn1^+/− female mice developed mammary tumors between the ages of 14 and 24 mo (Fig. 6C). Tumor #1 displayed squamous and focal glandular differentiation, with keratinizing differentiated squamous and nonkeratinized poorly differentiated regions. Tumor #2 displayed well-differentiated keratinizing squamous differentiation with foci of basaloid squamous and glandular components; both tumors #1 and 2 showed pervasive KRT6- and KRT14-positivity, but no PGR staining. Tumor #3 had characteristics of glandular differentiation with small foci of well-differentiated keratinizing squamous component, and was positive for KRT6 and KRT14, but not PGR. Tumors #1–3 were characterized as invasive adenosquamous carcinomas, a type of metaplastic carcinoma seen in the mammary gland. Tumor #4 displayed a marked cribriform pattern and luminal eosinophilic cylinders resembling the invasive mammary tumor known as an adenoid cystic carcinoma, and was KRT6-negative, but PGR-positive. Despite diversity in histology likely due to tumor heterogeneity,³⁹ most tumors were squamous, KRT6-positive, and basal-like (KRT14-positive) and all showed negative correlation between KRT6 and PGR expression.

Open in a separate window

Figure 6.

Spontaneously arising mammary tumors are detected in multiparous Becn1^+/− FVB/N mice. (A and B) Hematoxylin and eosin staining (H&E) of MG sections from Becn1^+/+ (left) and Becn1^+/− (right) FVB/N mice: (A) Representative images from MGs of nulliparous 21- to 22-mo-old mice reveals increased epithelial cell density and inflammatory cell infiltration (brown, identified with arrow) in Becn1^+/− mice. (B) Representative H&E images of MGs from parous 14- to 17-mo-old mice reveals increased epithelial cell density and disorganized cell layering in MGs from Becn1^+/− mice. (C) Representative images from 4 spontaneously arising mammary tumors in parous Becn1^+/− FVB/N mice (Tumors #1 to 4). Left panels are representative H&Es; middle panels, KRT6 (green) and nuclear PGR (red) staining; and right panels, KRT8 (green) and KRT14 (red) staining. (_*) and (_*) denote keratin pearls and basement membrane/eosinophilic hyaline materials, respectively. Tumor type denotations are identified by: (i) for glandular components, (ii) for poorly to moderately differentiated squamous cell carcinoma, (iii) for well-differentiated squamous cell carcinoma, and (iv) for basaloid squamous cell carcinoma. Scale bar: (A and B) 30 μm and KRT8 and KRT14 (C); (C) 200 μm for H&Es; (C) 50 μm for KRT6 and PGR.

Monoallelic Becn1 loss accelerates WNT1-driven mammary tumorigenesis

Given our findings that monoallelic Becn1 loss promotes KRT6-positive immature mammary cell accumulation (Fig. 1), we investigated whether Becn1 heterozygosity has an impact on mammary tumorigenesis in the context of oncogenic WNT1 activation, where KRT6-positive bipotent MaPCs are considered tumor-initiating cells.^18,20

As shown (Fig. 7A; Fig. S5A), Becn1^+/−;MMTV-Wnt1 mice displayed significantly shorter mammary tumor-free survival (4 vs. 7.2 mo; P 0.004) and overall survival (4.7 vs. 7.5 mo; P < 0.001) compared with Becn1^+/+;MMTV-Wnt1 mice. Consistent with defective autophagy status, mammary tumors from Becn1^+/−;MMTV-Wnt1 mice showed reduced LC3B levels and SQSTM1 aggregates (Fig. S5B). Tumors arising in Becn1^+/+;MMTV-Wnt1 mice displayed known characteristics of MMTV-Wnt1 mammary tumors,⁴⁰ as they were KRT6-positive adenosquamous carcinomas with high KRT8 and low KRT14 levels, and variable ESR (estrogen receptor) and PGR expression (Fig. 7B). Mammary tumors from Becn1^+/−;MMTV-Wnt1 mice were also mostly adenosquamous carcinomas, but exhibited higher KRT6 levels, were more basal-like, as KRT14 expression was very prominent, and did not express ER or PGR (Fig. 7B). Furthermore, while proliferation was similar, apoptosis was suppressed in tumors from Becn1^+/−;MMTV-Wnt1 mice (Fig. S5C and S5D).

Open in a separate window

Figure 7.

Monoallelic Becn1 loss accelerates WNT1-driven tumorigenesis and gives rise to mammary tumors with TNFRSF11A-NFKB1 pathway activation and basal-like characteristics. Basal-like characteristics and TNFRSF11A-NFKB1 signaling activation are detected in spontaneous mammary tumors that arise faster in Becn1^+/−;MMTV-Wnt1 (right) compared with Becn1^+/+;MMTV-Wnt1 mice (left). (A) Kaplan-Meier curves for mammary tumor-free survival in Becn1^+/+;MMTV-Wnt1 (n = 25) and Becn1^+/−;MMTV-Wnt1 (n = 20) mice demonstrates decreased mammary tumor-free survival in Becn1^+/−;MMTV-Wnt1 mice. (B–E) Examination of mammary tumors and premalignant MGs from Becn1^+/+;MMTV-Wnt1 (left) and Becn1^+/−;MMTV-Wnt1 (right) mice. (B) Representative images of H&E staining (top panel), KRT6 (green) and nuclear PGR (red) cell staining (second panel), KRT8 (green) and KRT14 (red) cell staining (third panel), and nuclear ESR1 staining (bottom panel) on tumor sections. (C) Heat map representation of microarray analysis confirms basal-like tumor characteristics along with higher CTNNB1 expression in tumors from Becn1^+/−;MMTV-Wnt1 mice. Statistical significance of upregulation of a given signature on a group of samples was determined using a Fischer exact test for enrichment of samples with significant (P < 0.05) signature expression in the group relative to samples outside the group. (D) Representative images of TNFSF11 (top panel) and RELA (bottom panel) expression in mammary tumors, and (E) TNFSF11 (top panel), TNFRSF11A (center panel), and RELA (bottom panel) expression in premalignant MGs from 6- to 10-wk-old cohoused Becn1^+/+;MMTV-Wnt1 and Becn1^+/−;MMTV-Wnt1 littermates. P < 0.01 determined by a Mantel-Cox test for Kaplan-Meier curves. Scale bar: (B) 200 μm for H&Es; (B, D, and E) 30 μm for all other panels.

Gene expression profiling also showed that mammary tumors from Becn1^+/−;MMTV-Wnt1 mice showed higher expression of KRT6 and other basal keratins, including KRT5, KRT14, and KRT17; they also displayed CTNNB1/β-catenin (catenin [cadherin-associated protein], β 1, 88 kDa) upregulation (Fig. 7C), indicating that monoallelic Becn1 loss rendered WNT1-induced mammary tumors basal-like and further enhanced CTNNB1 signaling already present due to WNT1 activation. An earlier publication reported that TP53 expression is regulated by BECN1 and, as such, may contribute to BECN1-associated tumor suppression.⁴¹ In our studies, mammary tumors from Becn1^+/+;MMTV-Wnt1 and Becn1^+/−;MMTV-Wnt1 mice did not exhibit differences in expression of TP53 transcriptional target genes (data not shown), thus suggesting that in the context of WNT1-driven mammary tumorigenesis, differential TP53 activation is not a likely contributor to mammary tumor acceleration in association with monoallelic Becn1 loss.

Consistent with our finding that monoallelic Becn1 loss results in increased mammary stem and progenitor cell activity in association with TNFRSF11A-NFKB pathway activation (Figs. 1, ​,3,3, and ​4),4), tumors from Becn1^+/−; MMTV-Wnt1 mice exhibited increased TNFSF11 expression and nuclear RELA staining (Fig. 7D). Premalignant MGs from Becn1^+/−;MMTV-Wnt1 mice displayed even higher TNFSF11, TNFRSF11A, and nuclear RELA levels (Fig. 7E), indicating that early activation of TNFRSF11A-NFKB signaling may contribute to promotion of WNT1-driven mammary tumorigenesis by monoallelic Becn1 deletion. However, additional studies must be conducted to further investigate whether aberrant TNFRSF11A pathway activation is responsible for the accelerated mammary tumor development in Becn1^+/−;MMTV-Wnt1 mice.

Monoallelic Becn1 loss and TNFRSF11A-NFKB signaling activation in mammary tissues

Our studies demonstrate that MGs from pubertal and adult Becn1^+/− mice exhibit excessive ductal side-branching in association with abnormal TNFSF11 expression and TNFRSF11A-NFKB axis activation. TNFSF11 is a member of the tumor necrosis factor (TNF) ligand superfamily, primarily known for its role in bone remodeling and metastasis³³ and also implicated in MaPC expansion and progestin-induced tumorigenesis.^34,46 High TNFRSF11A and TNFSF11 levels were also reported in breast cancers lacking ESR1 and PGR expression and in high-grade and highly proliferative breast tumors;⁴⁷ consistent with these results, we also found evidence of TNFRSF11A pathway activation in aggressive basal-like tumors exhibiting low BECN1 expression and an activated WNT pathway gene signature. Furthermore, TNFRSF11A overexpression in breast cancer cell lines has been reported to disrupt mammary cell fate, induce epithelial-to-mesenchymal transition (EMT) and promote mammary tumorigenesis and metastasis.⁴⁷ In our studies, tumors originating in both Becn1^+/+;MMTV-Wnt1 and Becn1^+/−;MMTV-Wnt1 mice retained epithelial keratin expression, and neither expressed mesenchymal markers, such as VIM/vimentin (data not shown) nor gave rise to visible metastases at the time of necropsy. However, since metastatic lesions are not usually noted in MMTV-Wnt1 mice sacrificed when their primary mammary tumors reach experimental endpoints,⁴⁰ specifically designed metastasis studies, involving mammary tumor resection survival surgery and/or intracardial injection of tumor-derived mammary tumor cell lines, will be necessary to determine the role of BECN1 in metastasis and EMT in the context of WNT pathway activation.

Our findings, together with the striking similarities in mammary development and tumor phenotypes between Becn1^+/− and MMTV-Tnfrsf11a (or MMTV-Tnfsf11) mice, suggest that the TNFRSF11A pathway may mediate the impact of monoallelic Becn1 loss on mammary cell hierarchy, morphogenesis and tumorigenesis. However, whether ectopic TNFRSF11A pathway activation is responsible for, or contributes to, MaSC and progenitor expansion and increased tumor susceptibility following parity or WNT pathway activation upon monoallelic Becn1 loss requires further mechanistic investigation.

Which cell type produces TNFSF11 in MGs from Becn1^+/− mice is another topic of interest. TNFSF11 production by MECs is primarily considered the result of progesterone stimulation,^28,32 whereas in macrophages, adipocytes, and breast cancer cells, proinflammatory cytokines can also upregulate TNFSF11.^36,48,49 Our transplantation studies argue in support of TNFSF11 production by MECs in response to proinflammatory cytokine release by MECs themselves: (1) mammary outgrowths were exposed to same systemic progesterone levels, as CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MECs were transplanted in contralateral fat pads of recipient mice. Increased progesterone signaling in association with monoallelic Becn1 loss must therefore result from a locally-induced expansion of PGR-expressing cells; (2) 2 proinflammatory cytokines known to induce TNFSF11 production, namely IL1B and TNF, were detected in CD24⁺ CD29^hi Becn1^+/− MEC-generated outgrowths; (3) CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MEC transplantations were performed in wild-type recipient mice, thus all stromal cells were wild type. The fact that, similar to MGs from Becn1^+/− mice, outgrowths from CD24⁺ CD29^hi Becn1^+/− MEC transplantation also exhibited increased ductal side-branching and TNFSF11 production indicate that MECs are sufficient to induce these mammary phenotypes. Furthermore, macrophage infiltration was similar in CD24⁺ CD29^hi Becn1^+/+ and Becn1^+/− MEC-generated outgrowths; (4) increased oxidative stress, a inducer of proinflammatory cytokine production,³⁷ was documented in CD24⁺ CD29^hi Becn1^+/− MEC-generated outgrowths; (5) autophagy-deficient status, commonly associated with inflammation and cytokine release,^14,37 was maintained in CD24⁺ CD29^hi Becn1^+/− MEC-generated outgrowths. This study provides evidence for an essential role of BECN1 in restricting mammary stem and progenitor cell expansion and limiting ectopic TNFRSF11A pathway activation, 2 pathways known to be deregulated in particular contexts of mammary tumorigenesis.^18,30,34

Shedding light in the complex role of Becn1 in mammary tumorigenesis

BECN1, and autophagy in general, have complex roles in tumorigenesis: allelic BECN1 deletions have been reported in human cancers,² and Becn1^+/−4,5 and other Atg-deficient^50,51 mice are tumor-prone, in support of a tumor suppressor role for autophagy in cancer. However, autophagy also preserves tumor cell survival under stress and likely contributes to tumor progression¹¹ and treatment resistance.⁵² While apoptosis-deficient Becn1^+/− iMMECs, compared with their wild-type counterparts, result in earlier mammary tumor development in nude mice,⁶ monoallelic Becn1 deletion does not alter ERBB2- or PyMT-driven,⁸ while it delays Palb2 loss-associated,¹³ mammary tumorigenesis.

Our studies define specific contexts, either following parity or oncogenic WNT1 activation, in which monoallelic Becn1 deletion promotes mammary tumorigenesis. Spontaneous tumors arising in Becn1^+/− mice after parity display mostly squamous cell characteristics, which is similar to the squamous differentiation seen by Wei and colleagues in lung cancers with low BECN1 expression.⁹ These 2 studies indicate that low BECN1 expression may drive particular lung and breast cancer subsets toward squamous, rather than adenocarcinoma, histologies, thus also raising the possibility of concurrent EGFR activation.

In the case of WNT pathway activation, we propose that WNT1-driven mammary stem and progenitor cell expansion cooperates with Becn1 heterozygosity-induced deregulation of the mammary hierarchy to create more transformation-susceptible immature mammary cells (Fig. 8E). Since monoallelic Becn1 loss results in hormone receptor-negative, basal-like mammary tumors, it is possible that corresponding tumor-initiating cells (TICs) reside within the expanded mammary stem and progenitor cell populations,²³ a hypothesis that is worthy of further investigation. If this is the case, other oncogenes deregulating the mammary hierarchy and resulting in expansion of progenitor cell populations may also cooperate with monoallelic Becn1 loss to promote mammary tumorigenesis.

Mammosphere and colony-forming assays

Primary mammosphere assays were performed as previously described.^23,25 MaSC viability assays under ultra low attachment were performed as previously described, and mammospheres were counted. 20,000 CD24⁺ CD29^hi MECs per well were plated on ultra-low attachment plates (Costar, 3471) for 10 to 15 d in Dulbecco's modified Eagle's medium (DMEM)/F12, EGF (epidermal growth factor; Sigma, E4127), FGF2 [fibroblast growth factor 2 (basic); Invitrogen, PHG0026], B27 (Gibco, 17504-044), and heparin (Sigma, H6279). Results are presented as means ± standard deviation (SD) from 3 independent experiments. Primary mammospheres (2 or more cells) were counted. For limiting dilution mammosphere assays, 25 to 500 CD24⁺ CD29^hi MECs per well were plated in duplicate.

The 3D colony-forming assays were performed as previously described.^23,25 20,000 CD24⁺ CD29^hi MECs per well were plated in reduced growth factor Matrigel (BD, 354230) and supplemented with mammary epithelial growth medium (MEGM) (Lonza, CC3051A), 5% FBS, EGF, and FGF2. Colonies were counted after 9 to 14 d. Results are presented as means ± SDs from 3 independent experiments.

The 2D colony-forming assays were performed as previously described.^23,25 1,000 CD24⁺ CD29^lo cells per well were plated on reduced growth factor Matrigel, supplemented with MEGM, 10% FBS, EGF, FGF2. Colonies were counted after 7 to 9 d. Results are presented as means ± SDs from 3 independent experiments.

Transplantation assays

MECs were isolated as previously described.²² FACS was used to collect CD24⁺ CD29^hi MECs, i.e., the MaSC-containing population, and transplantations were performed as previously described into cleared fat pads of 3- to 4-wk-old recipient wild-type C57BL/6 female mice,²² in accordance with IACUC-approved protocols. Primary MEC donors were littermates. CD24⁺ CD29^hi MECs were injected in a 200 μL volume (100 μL reduced growth factor Matrigel + 100 μL 1× phosphate-buffered saline (PBS; Gibco, 14190) at cell numbers specified. Mammary outgrowths were allowed to grow for 6 wk following transplantation and were then dissected, mounted on microscope slides, and processed for histological analysis. Three independent transplantation experiments were performed.

Mammary gland whole mounting

Native mammary glands from pubertal mice and mammary outgrowths generated by CD24⁺ CD29^hi MEC transplantation were mounted onto microscope slides and fixed in 10% formalin overnight. Mammary glands from aged mice were mounted on microscope slides and fixed in Carnoy fixative for de-fatting. Mounts were then subjected to graded ethanol to water washes, and stained with carmine solution (Sigma, C1022 prepared as directed, with Aluminum Potassium Sulfate, Sigma, A7167) overnight with shaking. Stained slides were then subjected to graded water to ethanol washes, destained and stored in Histo-Clear II (National Diagnostics, HS-202). Imaging was done on Olympus SZ61 dissecting microscope (Rutgers Cancer Institute of New Jersey) using the ProRes Capture Pro2.5 software.

Stable cell lines and culture conditions

Cells were maintained in F12 (Gibco, 11765), 10% FBS, 1% penicillin/streptomycin, 5 ng/mL EGF, 1 μg/mL hydrocortisone (Sigma, H0888), 5 μg/mL insulin (Sigma, I9278). Cells were plated on glass coverslips, and allowed to attach overnight. For immunofluorescence (IF), cells were fixed with 1:1 methanol:acetone at −20 °C for 10 min and rinsed with 1× PBS. Slides were incubated with 50 μL KRT6 primary antibody 1:200 in 3% bovine serum albumin (BSA; Benchmark, 100-106), 0.01% Tween 20 (BIO-RAD, 170-6531), 1× PBS at 37 °C for 90 min with shaking. After washing with 0.01% Tween 20, 1× PBS, slides were incubated with 50 μL FITC-conjugated anti-rabbit antibody (Jackson ImmunoResearch, 711-095-152) 1:100 in 3% BSA, 0.01% Tween 20, 1× PBS at 37 °C for 30 min with shaking. Cells were stained with DAPI (Invitrogen, P36931) and mounted with ProLong Gold Antifade Mounting Media (Life Technologies, P363930). Slides were imaged using the Olympus 1×51 scope and Metamorph Image program.

Primary MEC isolation and FACS for MEC subpopulations

Primary MECs were isolated as previously described.^22,23,25 MECs were stained as previously described^22,23 and sorted using FACSVantage SE with DiVa (BD) at Princeton University. LIN⁻ refers to cluster differentiation 31 (CD31)⁻, cluster differentiation 45 (CD45)⁻, and erythroid-specific LY76 (lymphocyte antigen 76)⁻ negative cells, i.e. MECs. CD24⁺ CD29^hi MECs refer to a MaSC-containing population, CD24⁺ CD29^lo to a mammary progenitor cell-containing, and CD24⁻ to differentiated MECs, as previously described.^22,23 FACS profiles were analyzed using FlowJo software (TreeStar, Inc.). Five independent experiments were performed.

Primary and stable cell lines

Immortalized mouse mammary epithelial cells were generated as previously described.⁶ MaSC and MaPC colony- and mammosphere-forming assays were performed as previously described.^23,25 MaSC transplantations were performed as previously described.²²

Histology of tissues and tumors

Mammary tissues were processed and imaged as previously described.¹⁵ Following Carmine staining, mammary glands and outgrowths from CD24⁺ CD29^hi MEC transplantation were transferred to 70% ethanol overnight, and processed as above. Hematoxylin and eosin (H&E) staining and imaging were performed by the Histopathological and Imaging Core Facility at Rutgers Cancer Institute of New Jersey.

Antibodies

For flow cytommetry, CD31 (BD558737), CD45 (B553078), LY76/TER119 (BD553672), and CD24-PE (BD553262) antibodies were from BD Biosciences, CD61-APC (MCD6105) from Invitrogen, and CD29-FITC (MCA2298FB) from Serotech were used. For immunohistology, KRT6 (PRB-169P) from Covance, KRT17 (ab53707), TNF/TNFα (ab8871), F4/80 (ab6640), and 8-O-dG (ab48508) from Abcam, KRT14 (gift from Dr Segre), PGR (gift from Dr Clarke), CASP3 (9661) from Cell Signaling Technology, MKI67 (NCL-L-Ki67-MM1) from Leica, ESR1/ERα (sc-542) and RELA/p65 (sc-8008) from Santa Cruz Biotechnology, TNFRSF11A/TNR11 (AF692) and TNFSF11/TNF11 (AF462) from R&D Systems, SQSTM1 (BML-PW9860) from Enzo, LC3 (NB100-2331) and IL1B/IL-1β (NBP1-19775) from Novus Biologics, and KRT8 (TROMA-1) and BrdU (G3G4) from the Developmental Studies Hybridoma Bank were used.

Immunostaining

For immunohistochemistry (IHC), slides were deparrafinized with xylenes, rehydrated with ethanol to water washes, and boiled at 97 °C for 16 min in 1× Antigen Unmasking Solution (Vector, H-3300). For PGR staining, retrieval was done by autoclaving at 121 °C, 16 to 17psi, for 30 min in Vector 1× unmasking solution. For ESR1 staining, retrieval was done by microwaving slides in 1× Vector Unmasking Solution for 20 min. Slides were washed in water, followed by endogenous peroxidase blocking with 3% H₂O₂. Washes in 0.05% Tween 20 in 1× PBS or 0.05% Triton X-100 (Sigma, X100) (for PGR) were performed before blocking for 1 h at room temperature in 5% goat serum albumin for rabbit antibody, 5% BSA for goat antibody, or MOM Blocking Reagent (Vector, MKB-2213) for mouse or rat antibodies. For ESR1 staining, retrieval was done by microwaving slides in 1× Vector Unmasking Solution for 20 min. For rabbit and goat antibodies, primary antibodies were added in blocking reagent overnight at 4 °C. Biotin-conjugated anti-rabbit 1:1000 (Novus, NB 730-B) in 5% goat serum albumin (Jackson ImmunoResearch, 005-000-001)or anti-goat 1:500 (Jackson ImmunoResearch, 711-065-152) in 5% BSA were added to slides for 1 h at room temperature. For rat and mouse antibodies, primary antibodies were added in MOM diluent (prepared as per MKB-2202) for 30 min at room temperature as directed. For rat antibody, biotinylated anti-rat (Dako, E0468) 1:500 was added for 1 h at room temperature. For mouse antibody, biotinylated anti-mouse antibody (MKB-2225) 1:250 in MOM diluent was added at room temperature for 10 min (15 min for PGR). For all antibodies, sections were washed with 0.05% Tween 20 in 1× PBS, then incubated in Vector Elite ABC reagent (PK-7100) for 30 min at room temperature. For ESR1, an additional Streptavidin-HRP (Vector, SA-5704) incubation step was performed for 30 min as room temperature. All sections were incubated in DAB + chromogen substrate (Dako, K3468) for 10 min (15 min for PGR and ESR1), and counterstained with Harris Modified Hematoxylin 1:5 (Fischer Scientific, SH26-500D). Slides were mounted with Permount medium (Fisher Scientific, SP15-100). Imaging was done with an Olympus 1×51 microscope, Jenoptik Laser Optik camera system (Rutgers Cancer Institute of New Jersey), and ProgRes C10^Plus software.

For immunofluorescence (IF), slides were deparrafinized with xylene, rehydrated with ethanol to 1× PBS washes, and boiled at 97 °C for 16 min in 1 M Urea (Bio-Rad, 161-0730) or 1× Vector Unmasking Solution (TNFSF11 only). For PGR staining, retrieval was done by autoclaving slides at 121 °C, 16 to 17 psi, for 30 min in 1× Vector Unmasking Solution. Slides were washed 4 times in water and once in 1× PBS, followed by blocking for 1 h at room temperature in 5% BSA (0.5% Tween-20, 1× PBS) or MOM Blocking Reagent (1× PBS for MKI67 and 0.05% Triton X-100, 1× PBS for PGR). Nonmouse primary antibodies were added in 5% BSA at room temperature for 2 h, except for the TNFSF11 antibody, which was added in 1% BSA, 0.5% Tween 20 overnight at 4 °C. PGR antibody was used in MOM diluent in 0.05% Triton X-100, 1× PBS for 2 h at room temperature. Anti-rat, anti-chicken, and anti-rabbit fluorescent antibodies 1:100 (Jackson ImmunoResearch) were used in 5% BSA for 1 h at room temperature. For TNFSF11, biotinylated anti-goat 1:500 (Jackson ImmunoResearch, 711-065-152) was used for 1 h at room temperature, followed by streptavidin-Texas Red 1:100 (Invitrogen, SA1017) for 30 min. For PGR, biotinylated anti-mouse 1:250 in MOM diluent in 0.05% Triton X-100, 1× PBS was used at room temperature for 15 min, followed by streptavidin-Texas Red 1:100. For KRT6, CD29 and CD24 costaining, slides were first stained with KRT6 and CD29 antibodies, anti-rabbit Alexa-647 1:100 (Invitrogen, A20991) and anti-rat TRITC, then with the CD24-PE. After all antibody incubations, sections were washed with 1× PBS 3 times and mounted with ProLong Gold antifade with DAPI. Confocal imaging was done using the LSM510 Meta Confocal Laser Scanning Microscope at the Keck Collaborative Neuroscience Center at Rutgers University.

Data quantification and statistical analysis

For all in vitro studies statistics, 2-tailed Student t tests were performed using http://graphpad.com/quickcalcs/ttest1/. Immunostaining quantification was performed on at least 3 individual mammary gland specimens, including all TEBs and ducts per slide. For mammary side-branching quantification, results are presented as means with standard deviation (± SD) from 3 mammary gland specimens per genotype from C57BL/6 littermates. Same results were obtained in FVB/N background. For statistics on Kaplan Meier curves, Graph Prism software and Mantel-Cox analysis were utilized.

Mouse mammary tumor gene expression analysis

A minimum of 1 million cells were harvested and lysed directly from culture plates using 2 mL Trizol (Invitrogen, 15596-026). Cell lysates were incubated on ice for 5 min, then vortexed for 1 min. Tumors were dissected from mice and snap frozen in liquid nitrogen. Cell lysates and tumor specimens were given to the Functional Genomics Core Facility at Rutgers CINJ for RNA extraction and gene expression profiling on Affymetrix Mouse Genome 430 2.0 microarray. Using GeneSpring GX 11.5 (Agilent Technologies, Inc., Santa Clara, CA, USA), raw gene expression signals were summarized with MAS5. The data was further normalized using baseline transformation by the median of all samples and log scaling. Further, the normalized expression signals were averaged between biological replicates where applicable. Gene expression data was first filtered by flags (labeled as “present,” “marginal,” or “absent”). Probes were only kept if at least one sample was labeled as “present.” Gene lists were generated by pair-wise comparison. Differentially expressed genes were designated by looking for a significant fold change (> 2.0) between the normalized expression values of 2 experiments. Genes were considered significantly stagnant if the fold change was low (< 1.25).

The generated gene list data were analyzed through the use of Ingenuity Pathways Analysis (Ingenuity® Systems, www.ingenuity.com). The Functional Analysis identified the biological functions and/or diseases that were most significant to the data set. Molecules from the data set that met the fold change cutoff of 2-fold and were associated with biological functions and/or diseases in Ingenuity's Knowledge Base were considered for the analysis. The right-tailed Fisher exact test was used to calculate a P value determining the probability that each biological function and/or disease assigned to that data set is due to chance alone.

Human breast cancer cohorts

The Sabatier⁴³ cohort included 266 breast tumor samples analyzed with Affymetrix U133 Plus 2.0 Arrays [PMID: 20490655]. The data was downloaded from Gene Expression Omnibus, series GSE21653. The Hatzis⁴² cohort included 508 ERBB2-negative breast tumor samples analyzed with Affymetrix U133A Arrays [PMID: 21558518]. The data was downloaded from Gene Expression Omnibus, series GSE25066. Raw CEL files were processed using the justRMA function in R Bioconductor. The samples were classified into BECN1 high or BECN1 low depending on whether BECN1 mRNA expression (probe 208945_s_at) was above or below the average across samples. WNT pathway gene signature was based on KEGG pathway: hsa04310 for WNT signaling pathway-Homo sapiens (http://www.genome.jp/dbget-bin/www_bget? hsa04310). An activated WNT pathway gene signature was attributed to tumors with significant (P < 0.05) upregulation of WNT pathway gene signature expression relative to the rest of tumors within each breast cancer cohort; a Fischer exact test was used for this calculation.

We thank Z Yue and S Jin (Mount Sinai School of Medicine and Rutgers University), Dr Komatsu (Tokyo Metropolitan Institute of Medical Science), Dr Mizushima (University of Tokyo), and Y Li (Baylor College of Medicine) for providing Becn1^+/−, Atg7^F/F, GFP-LC3B, and MMTV-Wnt1 mice, respectively. We thank Dr Segre (National Institute of Health) for the KRT14 antibody and Dr Clarke (Westmead Millennium Institute) for the PGR antibody. We also thank N Goldsmith-Kane for assistance with confocal microscopy, the Histopathological, Imaging, and Functional Genomics Core Facilities at Rutgers Cancer Institute of New Jersey; also, the Flow Cytometry Facility at Princeton University. We are indebted to Ed and Sue Goldstein for their unconditional support and encouragement. This work is dedicated to the memory of Stacy Goldstein.