Skip to main content
  • Research article
  • Open access
  • Published:

Parity induces differentiation and reduces Wnt/Notch signaling ratio and proliferation potential of basal stem/progenitor cells isolated from mouse mammary epithelium

Abstract

Introduction

Early pregnancy has a strong protective effect against breast cancer in humans and rodents, but the underlying mechanism is unknown. Because breast cancers are thought to arise from specific cell subpopulations of mammary epithelia, we studied the effect of parity on the transcriptome and the differentiation/proliferation potential of specific luminal and basal mammary cells in mice.

Methods

Mammary epithelial cell subpopulations (luminal Sca1-, luminal Sca1+, basal stem/progenitor, and basal myoepithelial cells) were isolated by flow cytometry from parous and age-matched virgin mice and examined by using a combination of unbiased genomics, bioinformatics, in vitro colony formation, and in vivo limiting dilution transplantation assays. Specific findings were further investigated with immunohistochemistry in entire glands of parous and age-matched virgin mice.

Results

Transcriptome analysis revealed an upregulation of differentiation genes and a marked decrease in the Wnt/Notch signaling ratio in basal stem/progenitor cells of parous mice. Separate bioinformatics analyses showed reduced activity for the canonical Wnt transcription factor LEF1/TCF7 and increased activity for the Wnt repressor TCF3. This finding was specific for basal stem/progenitor cells and was associated with downregulation of potentially carcinogenic pathways and a reduction in the proliferation potential of this cell subpopulation in vitro and in vivo. As a possible mechanism for decreased Wnt signaling in basal stem/progenitor cells, we found a more than threefold reduction in the expression of the secreted Wnt ligand Wnt4 in total mammary cells from parous mice, which corresponded to a similar decrease in the proportion of Wnt4-secreting and estrogen/progesterone receptor-positive cells. Because recombinant Wnt4 rescued the proliferation defect of basal stem/progenitor cells in vitro, reduced Wnt4 secretion appears to be causally related to parity-induced alterations of basal stem/progenitor cell properties in mice.

Conclusions

By revealing that parity induces differentiation and downregulates the Wnt/Notch signaling ratio and the in vitro and in vivo proliferation potential of basal stem/progenitor cells in mice, our study sheds light on the long-term consequences of an early pregnancy. Furthermore, it opens the door to future studies assessing whether inhibitors of the Wnt pathway may be used to mimic the parity-induced protective effect against breast cancer.

Introduction

Pregnancy is the most significant modifiable factor known for breast cancer risk in women. Although an initial increase in risk occurs immediately after parturition in women older than 25 years, the overall lifetime risk of breast cancer decreases after pregnancy [1, 2]. This protective effect is > 50% if a full-term pregnancy has occurred before the age of 20 years [1]. Similarly, pregnancy and pregnancy-mimicking hormones have a strong protective effect against mammary tumors in rodents. This is true both for carcinogen-induced mammary tumors [3] and for genetically engineered mouse models of breast cancer [4].

The cellular and molecular mechanisms underlying the breast cancer-protective effect of early pregnancy remain unclear. Frequently raised hypotheses involve cell nonautonomous mechanisms such as systemic changes in circulating hormones and/or changes in the stromal composition of the mammary gland [5, 6], and cell autonomous processes such as changes in the differentiation state of mammary epithelial cells [7]. Furthermore, numerous parity-induced changes in gene expression have been identified in genome-wide expression profiles of entire lobular breast tissues of women or entire mammary glands of rats and mice [8–10]. However, it is not known to what extent these tissue studies reflect alterations in gene-expression profiles of distinct mammary epithelial cell subpopulations. Hence, given that breast cancers arise from specific subpopulations of mammary epithelial cells [11], investigations of early parity-induced gene-expression changes in distinct mammary epithelial cell subpopulations are warranted.

The mammary epithelium is hierarchically organized into differentiated luminal and basal (myoepithelial) cells, luminal and basal progenitor cells, and mammary stem cells [12, 13]. Whereas the latter were originally thought to lie exclusively in the basal compartment and to be multipotent (able to form both luminal and basal epithelial cells), recent lineage-tracing experiments indicated the existence of unipotent basal and luminal mammary stem cells and identified multipotent mammary stem cells solely in the embryonic and possibly in the pregnant gland [14, 15]. Distinct mammary epithelial cell subpopulations, including luminal progenitor and basal stem/progenitor cells can be isolated with fluorescence-activated cell sorting (FACS) by using specific cell-surface markers from both parous and virgin mice [16–21]. Whereas progenitor cells in general can be characterized in vitro by their colony-forming potential [16, 22, 23], the basal stem/progenitor cell subpopulation has the additional capacity to repopulate deepithelialized mouse mammary fat pads in vivo [16–18]. Although previous studies in total mammary epithelial cells indicated either no change or a decrease in the mammary repopulating capacity after parity [20, 21], the consequences of parity on the transcriptome and functionality of specific mammary epithelial cell subpopulations have not been investigated.

Therefore, we examined in this study whether pregnancy alters the gene-expression profiles ("gene signature") and the differentiation/proliferation potentials of the various mammary epithelial cell subpopulations. The results indicate that early parity decreases Wnt4 expression in luminal epithelial cells, leading to a reduction in the Wnt/Notch signaling ratio specifically in basal stem/progenitor cells. As expected, the decrease in the Wnt/Notch signaling ratio is associated with a concomitant strong prodifferentiation and antiproliferation phenotype in basal stem/progenitor cells. Because a decrease in Wnt signaling is known to have an anticarcinogenic effect [24, 25], the findings support the hypothesis that a reduction in the Wnt/Notch signaling ratio in basal mammary stem/progenitor cells plays a role in the mitigating effect of early pregnancy on breast tumorigenesis.

Methods

Animals and animal experimentation

All experiments were conducted in genetically homogenous FVB/NHanHsd mice purchased from Harlan Laboratories. The mice were bred and maintained in the animal facility of the Friedrich Miescher Institute, according to the Swiss guidelines on animal experimentation. All experiments were performed under permit 2159-2, in accordance with the animal-welfare ordinance and approved by the cantonal veterinary office of Basel Stadt, Switzerland. For the early-pregnancy protocol, mice were time-mated when 42 days old and allowed to lactate for 21 days. The postweaning period until cell harvest was 40 days, unless stated otherwise. To control for the estrus cycle in the transcriptome analyses, at least five mice were grouped for gland harvesting. For immunohistochemical analyses and determination of blood progesterone levels, mice in estrus were used, as assessed by the presence of a vaginal plug after an overnight mating. Age-matched virgin control mice were maintained under the same conditions as parous mice.

Whole mounts

Whole mounts were prepared by fixing the glands on glass slides with methacarnoy solution (60% methanol, 30% chloroform, 10% glacial acetic acid) for 4 hours at room temperature. The mounts were hydrated by sequential incubation in ethanol solutions of decreasing concentration: 100%, (overnight), 70%, 50%, and 30% (15 minutes each), distilled water (2 × 5 minutes), and stained overnight with an aqueous solution of 2% carmine (Sigma, Buchs, Switzerland) and 5% aluminum potassium sulfate (Sigma, Buchs, Switzerland). The mounts were dehydrated in ethanol solutions (70%, 90%, 95%, and 2 × 100%, for 15 minutes each) and cleared with xylene overnight. Images were captured with an Epson Expression 1600 Pro scanner.

Mammary cell preparation

The fourth mammary glands from virgin and parous mice were collected after lymph node removal and pooled. Mammary epithelial organoids were prepared as described [26]. Adipocytes were removed by repeated centrifugations (300 g). Red blood cells (RBCs) were eliminated by incubation with 8.3 g/L ammonium chloride (pH 7.5) for 5 minutes. The number of fibroblasts in gland extracts was reduced by their attachment to polystyrene cell-culture flasks (Corning, Buchs, Switzerland) during a 45-minute incubation step in Dulbecco modified Eagle medium (Invitrogen, Zug, Switzerland) with 10% FCS (Sigma, Buchs, Switzerland) at 37°C/5% CO2 [26]. The epithelial organoids were directly processed to single-cell suspensions by digestion in Hyclone HyQTase (ThermoScientific, Lausanne, Switzerland) with gentle pipetting for 15 minutes at 37°C. The cell suspension was filtered through a 40-μm cell strainer (BD Falcon, Basel, Switzerland), and the final cell suspension pelleted at 650 g for 4 minutes.

Cell labeling and flow cytometry

Cells were labeled as previously described [26] by using the antibodies PE-Cy7-CD45 (1:33), FITC-CD24 (1:40), PE-CD49f (1:40), and APC-Sca1 (1:40). Detailed antibody information is given later. DAPI (0.2%, Invitrogen, Zug, Switzerland) was added 10 minutes before cell sorting (1:250). FACS was carried out on a MoFlo cell sorter (Becton Dickinson, Basel, Switzerland). Cells were gated based on their forward- and side-scatter profiles (FS Area and SS Area). A time-of-flight approach (pulse width) was used to exclude doublets and higher-order cell clumps. Dead cells (DAPI bright) and immune cells (CD45+) were gated out (see Additional file 1). The gate for basal stem/progenitor cells was set at the top 5% of CD49f-expressing cells, as described [19, 21] (see Figure 1). Routine examination of the sorted mammary epithelial cell subpopulations showed a degree of purification higher than 95%.

Figure 1
figure 1

The CD24/Sca1 and CD49fHigh/CD24 flow-cytometry profiles of parous and age-matched virgin mice are similar. (A) Schematic illustration of the cell-isolation strategy and representative flow-cytometry pseudocolor plots of mammary cells from age-matched virgin control mice. After depletion of CD45+ white blood cells, luminal and basal mammary epithelial cells were separated on the basis of CD24 and Sca1 expression. Further separation of basal cells into myoepithelial and basal stem/progenitor cell subpopulations was based on the expression of CD24 and CD49f. The isolated mammary epithelial cell subpopulations included luminal Sca1+ (CD24+HighSca1+) cells, luminal Sca1- (CD24+HighSca1-) cells, basal CD49fHigh (CD24+LowSca1-CD49fHigh) or basal stem/progenitor cells, and basal myoepithelial (CD24+LowSca1-CD49fLow) cells. (B) Outline of the mouse mating, parturition, weaning, and involution protocol. (C) Representative flow-cytometry pseudocolor plots of mammary cells from parous mice. The gates applied were the same as those for age-matched virgin controls. (D) Bar graph showing the distribution of mammary epithelial cell subpopulations comparing cells from parous with age-matched virgin control mice. Data represent the mean ± SEM of seven cell-isolation experiments with a minimum of 10 mice per experiment. The proportion of luminal Sca1+ cells was reduced by approximately 50% in parous mice (P = 0.02 with a two-tailed unpaired Student t test).

In vitro colony formation assay and quantification

Freshly sorted cells of each subpopulation (500 cells) were plated onto irradiated 3T3L1 feeder cells (5,000 cells) in 24-well Primaria plates (Becton Dickinson, Basel, Switzerland). The cells were cultured in DMEM/Ham F12 (Invitrogen, Zug, Switzerland) with 10% FCS, 100 IU/ml penicillin, 100 μg/ml streptomycin (Invitrogen, Zug, Switzerland), 5 μg/ml bovine pancreatic insulin (Sigma, cell culture-tested solution, Buchs, Switzerland), 10 ng/ml mouse EGF (BD Biosciences, Basel, Switzerland), and 10 ng/ml cholera toxin (Sigma, Buchs, Switzerland). After 24 hours, the medium was renewed, and 4 days later, the colonies were fixed with acetone/methanol (1:1), washed, and rehydrated with PBS and 0.05% NaN3. In separate rescue experiments with selected mammary epithelial cells from parous mice, the incubation medium was supplemented with recombinant mouse Wnt4 (R&D Systems, Abingdon, UK) at 500 ng/ml. For quantification, colonies were immunofluorescently stained with Krt18 and Krt14 antibodies and with Hoechst 33342. The stained colonies were imaged with a MacroFluo Z16 microscope (Leica, Heerbrugg, Switzerland) at 2× magnification, and a Z1 microscope (Zeiss, Feldbach, Switzerland) at 5× magnification. The feeder cells served as negative controls for Krt18 and Krt14 staining. The number of colonies (that is, clusters of more than three cells [23]) per well was determined manually. Colonies defined as Krt18/Krt14 double-positive were double positive over a minimum of 20% of colony area.

Mammary fat pad transplantation

Freshly sorted basal stem/progenitor cells were resuspended in PBS with 50% Matrigel (BD Biosciences, Basel, Switzerland) and injected (30 μl) in limiting dilution numbers into the cleared fourth mammary glands of 3-week-old syngeneic female mice. Glands from recipient mice were harvested 8 to 9 weeks after transplantation, processed as whole mounts (see earlier), and scanned with an Epson Expression 1600 Pro scanner. The number of outgrowths was counted, with an outgrowth defined as an epithelial structure with ducts starting from a central point and with lobules and/or terminal end buds (see Additional file 2) [27].

Microarray analyses

Microarray analyses were performed on unsorted total mammary cell suspensions (see Figure 2) and on FACS-sorted mammary cell subpopulations (see Figures 3 and 4, as well as Table 1). For unsorted mammary cells, RNA was extracted by using the RNeasy Plus Mini Kit (Qiagen, Hilden, Germany). Genomic DNA was removed by using gDNA Eliminator Mini Spin Columns (Qiagen, Hilden, Germany). The RNA concentration was measured with a Nanodrop 1000 machine, and RNA quality determined with an Agilent 2100 bioanalyzer and RNA Nano Chips. Aliquots of 100 ng of isolated total RNA were amplified by using the Ambion WT Expression kit (Ambion, Zug, Switzerland). For FACS-sorted mammary epithelial cell subpopulations, total RNA was isolated by using the Arcturus PicoPure Frozen RNA Isolation Kit (Life Technologies, Zug, Switzerland). Genomic DNA was removed by using RNase-Free DNase (Qiagen, Hilden, Germany), the RNA concentration determined by using the RiboGreen Assay, and RNA quality assessed by using an Agilent 2100 bioanalyzer and RNA Pico Chips. Total RNA was used as the input for synthesis of amplified cDNA with the NuGen Ovation Pico WTA System (NuGen Inc., Leek, The Netherlands).

Figure 2
figure 2

Parity-induced gene signature in total mammary cell suspensions. (A) Heat-map and cluster analysis of differential gene expression in total isolated mammary cells from age-matched virgin control and parous mice. Gene expression is presented as normalized Z-scores, defined as Z = (x-μ)/sd to allow visualization. A cut-off of P < 0.05 was applied. If multiple probe sets existed for the same gene, the probe set with the largest change in expression was selected. Four independent experiments were performed with 10 mice (five virgins; five parous) per experiment. (B) qPCR validation of pregnancy-induced gene-expression changes in total mammary cell suspensions. Fold changes are shown in relation to expression in cell suspensions from age-matched virgin control mice. Ct values were normalized to the reference gene Hprt. As control for luminal epithelial cell number, qPCR for Krt8 was performed. Data represent the mean ± SEM of four independent experiments with 10 mice (five virgins; five parous) per experiment.

Figure 3
figure 3

Parity leads to differentiation and decreases the Wnt/Notch signaling ratio in basal stem/progenitor cells. (A) Bar graph depicting the number of gene-expression changes in FACS-sorted mammary stromal and epithelial cell subpopulations from parous mice compared with age-matched virgin control mice by using a cut-off of fold change > 1.5 and an adjusted P value < 0.05. By far the most gene-expression changes were observed in basal stem/progenitor cells. Three independent experiments were performed with 10 mice (five virgins; five parous) per experiment. (B) Schematic illustration of prominent gene-expression changes in FACS-sorted basal stem/progenitor cells from parous as compared with age-matched virgin control mice. Fold changes are shown in parentheses with upregulated genes denoted as positive (+), and downregulated genes, as negative (-). Differentiation genes were upregulated (blue), Wnt target genes were downregulated (green), Wnt inhibitor Sfrp1 was upregulated (green) and overall Notch signaling (orange) was increased in basal stem/progenitor cells from parous mice. (C) qPCR validation of the changes in gene expression in basal stem/progenitor cells of parous mice. All classic Wnt target genes were downregulated, including Lgr5, Axin2, and versican (Vcan), whereas the more ubiquitously regulated target Myc was unchanged. In all cases, fold changes are shown relative to cells from age-matched virgin control mice. Ct values were normalized to the reference genes Hprt and Ubc [63]. Data represent the mean ± SEM of three independent experiments with 10 mice (five virgins; five parous) per experiment. (D) Representative images and quantification of immunostaining for the Wnt target gene versican in mammary gland sections from age-matched virgin and parous mice in estrus. Scale bar, 50 μm. Quantitative data represent the mean ± SEM from 60 randomly selected images from three virgin and three parous mice.

Figure 4
figure 4

Effects of parity on Wnt transcription-factor activities and nuclear β-catenin in basal mammary epithelial cells. (A) Wnt transcription factor LEF1/TCF7 motif activity in basal stem/progenitor cells from parous as compared with virgin control mice, as predicted by MARA [32]. The binding motif of the LEF1/TCF7 transcription factor is shown in color. (B) Representative images of immunostaining for β-catenin in mammary gland sections from age-matched virgin and parous mice in estrus. Arrow, basal mammary epithelial cells with nuclear β-catenin. Arrowhead, basal mammary epithelial cells lacking nuclear β-catenin. Scale bar, 10 μm. (C) Bar graph representing the quantification of nuclear β-catenin in mammary gland sections of virgin and parous mice. Data represent the mean ± SD (virgin mice: n = 3; parous mice: n = 3). P = 0.004 with the two-tailed unpaired Student t test. (D) Transcription-factor activities in basal stem/progenitor cells, predicted on the basis of target gene expression by Ingenuity IPA [40]. Activity is reported as Z scores (positive Z score, upregulation; negative Z score, downregulation) by using a cut-off of linear fold change > 1.5 and P value < 0.05. (E) BioFunctions most strongly downregulated (blue) or upregulated (orange) in basal stem/progenitor cells, as calculated by Ingenuity IPA [40]. The color of the fields represents the Z score; the size of the fields represents the number of genes. A stringent cut-off of linear fold change > 2 and P value < 0.01 was used to minimize the number of false positives.

Table 1 Wnt signaling is decreased and Notch signaling is upregulated in basal stem/progenitor cells

Resulting double-stranded cDNA was fragmented and labeled by using the Affymetrix GeneChip WT Terminal Labeling kit (Affymetrix, High Wycombe, UK). Affymetrix Gene Chip Mouse gene 1.0 ST microarrays were hybridized according to the GeneChip Whole Transcript (WT) Sense Target Labeling Assay Manual (Affymetrix, High Wycombe, UK) with a hybridization time of 16 hours. The Affymetrix Fluidics protocol FS450_0007 was used for washing. Scanning was performed with Affymetrix GCC Scan Control Software v. 3.0.0.1214 on a GeneChip Scanner 3000 7 G with autoloader. Arrays were normalized, and probeset-level expression values calculated with R/Bioconductor's (v2.14) "affy" package by using the rma() function. Differential gene expression between experimental and control samples was determined by using linear modeling as implemented in the R/Bioconductor package "limma." For general analysis of gene expression in total mammary cell suspensions and in FACS-sorted mammary cell subpopulations, we used the cut-off linear fold change > 1.5, adjusted P value < 0.05, and average linear expression between conditions greater than 4. To determine the 10 most up- or downregulated genes in FACS-sorted mammary epithelial cell subpopulations, a cut-off of linear fold change > 2.0, adjusted P value < 0.001, and average linear expression between conditions ave > 4 was used. For FACS-sorted mammary epithelial cell subpopulations both resulting lists of differential genes were imported into Ingenuity IPA (Ingenuity, content version 12710793) for pathway analysis. Gene-set enrichments in FACS-sorted mammary epithelial cell subpopulations were determined by using the JAVA application from the Broad Institute v2.0 [28, 29] and gene sets v2.5 and v3.0, as well as custom gene sets (see Additional file 3).

The microarray data from this publication have been submitted to the NCBI Gene Expression Omnibus [30] and are deposited as GSE40875 (mouse mammary cell subtypes), GSE40876 (total mammary epithelial cells in mice), and GSE40877 (both total mammary epithelial cells and mammary cell subtypes in mice).

Quantitative PCR

RNA was isolated as described earlier and converted to cDNA by using the WT-Ovation Exon Module Version 1.0 (NuGen Inc., Leek, The Netherlands). Real-time PCR was performed on the unamplified cDNA corresponding to the specified number of cells or on 25 ng of amplified cDNA (by using the NuGen Ovation Pico WTA System (NuGen Inc., Leek, The Netherlands) for amplification). The TaqMan probe-based system was applied in combination with the TaqMan Universal PCR Master Mix (Applied Biosystems, Zug, Switzerland). The probe IDs are given later. Cycling was performed with 7500 Fast and Step OnePlus Real-Time PCR Systems (both from Applied Biosystems, Zug, Switzerland).

Probe IDs for quantitative PCR

These were the probes: Areg (Mm00437583_m1), Axin2 (Mm00443610_m1), B2M (Mm00437762_m1), Cdk8 (Mm01223097_m1), Cel (Mm00486975_m1), Csn1s1 (Mm00514430_m1), Csn2 (Mm04207885_m1), Dsc2 (Mm00516355_m1), Dusp1 (Mm00457274_g1), Esr1 (Mm00433149_m1), Fst (Mm00514982_m1), Gata3 (Mm00484683_m1), Hprt (Mm00446968_m1), Id3 (Mm00492575_m1), Igfbp3 (Mm_01187817_m1), Jag2 (Mm01325629_m1), Krt8 (Mm00835759_m1), Krt14 (Mm00516876_m1), Krt15 (Mm00492972_m1), Krt19 (Mm00492980_m1), Lgr5 (Mm00438890_m1), Ltf (Mm00434787_m1), Maml2 (Mm00620617_m1), Myc (Mm00487804_m1), Pgr (Mm00435628_m1), Sfrp1 (Mm00489161_m1), Ubc (Mm01201237_m1), Vcan (Mm01283063_m1), Wnt4 (Mm01194003_m1).

Immunofluorescent staining

For single-cell staining, freshly sorted cells were allowed to air-dry on poly-L-lysine-coated slides and stored at -20°C. The dried cells were blocked with PBS, 2.5% goat serum, and 0.05% NaN3 for 60 minutes under UV light. The UV treatment was used to attenuate residual fluorescence from bound FACS antibodies. Primary antibody staining was performed overnight at 4°C by using Krt18 antibody (1:1,000) and Krt14 antibody (1:500) as luminal and basal cell markers, respectively. Secondary antibody staining was carried out for 60 minutes at room temperature by using anti-guinea pig Ig-Alexa488 (1:500) and anti-rabbit Ig-Alexa546 (1:1,000). Hoechst 33342 (0.5 μg/ml; Invitrogen, Zug, Switzerland) staining was performed for 10 minutes at room temperature. The stained cells were mounted with ProLong Gold antifade reagent (Invitrogen, Zug, Switzerland) and imaged with a Z1 microscope (Zeiss, Feldbach, Switzerland) at 63× magnification. No primary antibody was added as a negative control.

Immunohistochemistry

The fourth mammary glands of parous or age-matched virgin control mice were collected 40 days after weaning from mice in estrus. The glands were fixed in 4% PFA and embedded in paraffin. For immunostaining with ERα and PR, the sections were dewaxed and subjected to antigen retrieval by boiling in 10 mM citrate buffer for 10 minutes. Subsequently, the sections were cooled to room temperature, quenched for 10 minutes with PBS and 3% H2O2, and blocked for 30 to 60 minutes with PBS and 2.5% NGS. ERα and PR primary antibody staining was performed overnight at 4°C at a 1:1,000 (ERα) and a 1:200 (PR) dilution. Secondary antibody staining was carried out for 30 to 60 minutes at room temperature with biotinylated anti-rabbit IgG. Immunohistochemistry for versican, β-catenin, and p21 was performed on the Ventana DiscoveryXT instrument (Roche Diagnostics, Rotkreuz, Switzerland) by using the Research IHC Dap Map XT procedure. In brief, dewaxing was performed in the machine, and slides were pretreated with mildCC1 (versican and p21) or standardCC1 (β-catenin) (Roche Diagnostics, Rotkreuz, Switzerland). Primary antibodies were incubated at 37°C for 1 hour at the following dilutions: versican, 1:50; p21, 1:50; and β-catenin, 1:500. As secondary antibody, Immpress, an anti-rabbit HRP conjugated polymer, was used. All slides were counterstained with Hematoxylin II (Roche Diagnostics, Rotkreuz, Switzerland) and Bluing Reagent (Roche Diagnostics, Rotkreuz, Switzerland). Images were captured at 20-fold (ERα or PR) or 40-fold (versican, β-catenin, p21) magnification by using an Eclipse E600 microscope (Nikon, Egg, Switzerland). For quantification of ERα and PR positivity, at least 2,000 epithelial cells per mouse were counted. For quantification of nuclear β-catenin and p21, at least 500 epithelial cells per mouse were counted. The quantification of versican was performed with the MATLAB software by using color segmentation based on Mahalanobis distance to determine the pixels with a particular RGB-color distribution [31].

Determination of blood progesterone concentrations

Blood was collected from the right atrium of mice in estrus in EDTA-covered tubes (Sarstedt, Numbrecht-Rommelsdorf, Germany), and plasma was extracted by centrifugation at 1,500 g for 15 minutes at 4°C. Progesterone concentrations were assessed with ELISA, as specified by the manufacturer's guidelines (DRG, catalog no. EIA-1561, Marburg, Germany).

Antibodies

For flow cytometry, the following anti-mouse antibodies were used: PE-Cy7-CD45 (clone 30-F11), FITC-CD24 (clone M1/69), PE-CD49f (clone GoH3) (all from BD Pharmingen, Basel, Switzerland), and APC-Sca1 (clone E13-161.7, from Biolegend, San Diego, USA).

For immunofluorescence staining, the primary anti-mouse antibodies used were guinea pig keratin 18 (Krt18; Fitzgerald; catalog no. GP11, North Acton, USA) and rabbit keratin 14 (Krt14; ThermoScientific, catalog no. RB-9020, Lausanne, Switzerland). The secondary antibodies used were anti-guinea pig Ig-Alexa488 and anti-rabbit Ig-Alexa546 (Molecular Probe, Invitrogen, Zug, Switzerland).

For immunohistochemistry staining, the primary antibodies included estrogen receptor alpha (ERα; Santa Cruz Biotechnology; catalog no. SC-542, Dallas, USA), progesterone receptor (PR; Clone SP2, ThermoScientific; catalog no. RM-9102, Lausanne, Switzerland), versican (Millipore/Chemicon; catalog no. ab1033, Billerica, USA), β-catenin (Cell Signaling; catalog no. 9587, Danvers, USA), and p21 (Abcam; catalog no. ab2961, Cambridge, UK). The secondary antibodies used were biotinylated anti-rabbit IgG (Vector Labs; catalog no. BA-1000, Petersborough, UK) and anti-rabbit HRP conjugated polymer Immpress (Vector Labs; catalog no. MP-7401, Petersborough, UK).

Statistics

The limited dilution transplantation data were analyzed statistically as published previously [18]. The two-tailed unpaired Student t test was used to determine statistical significance of comparisons.

Motif Activity Response Analysis (MARA)

The MARA model [32] combines knowledge of gene-expression levels (measured by microarray) with transcription-factor binding sites to answer the question of which transcription factors are driving expression changes in mammary stem/progenitor cells in parous as compared with age-matched virgin control mice. Specifically, log-expression levels of all genes present on the microarray were modeled as linear combinations of transcription factor activities. The coefficients of these combinations were determined by the number of transcription-factor binding sites in the proximal promoter regions. For each transcription-factor binding motif m and each sample (microarray) s, we estimated the activity A ms with the corresponding error. Furthermore, we quantified the significance of activity change of each binding motif in parous as compared with virgin control mice.

Results

Early pregnancy decreases luminal Sca1+ cells, but does not change the proportions of the other mammary epithelial cell subpopulations

To investigate the influence of early parity on the proportions of mammary epithelial cell subtypes, we first established FACS profiles of epithelial cell subpopulations in virgin FVB control mice [21, 26]. Luminal CD24+High Sca1+ cells, luminal CD24+High Sca1- cells, basal CD24+Low Sca1- CD49fHigh cells, and basal CD24+Low Sca1- CD49fLow myoepithelial cells were isolated from these mice (Figure 1A). Use of the established cell markers keratin 18 (Krt18) and keratin 14 (Krt14) confirmed the luminal and basal origin of the isolated cell subpopulations, and qPCR for CD49f and Sca1 affirmed the purity of the isolated cell subpopulations (see Additional file 4). Because basal CD49fHigh cells are considered to be enriched for basal mammary stem/progenitor cells in virgin mice [19, 21], the term "basal CD49fHigh cells" is used synonymously with "basal stem/progenitor cells" throughout this article. To extend these analyses to parous mice, we confirmed that involution was complete 28 and 40 days after weaning (see Additional file 5). To allow a margin of safety, 40 days after weaning was used for all subsequent cell-isolation experiments in a standardized parturition protocol with mating at 42 days (Figure 1B). The FACS profiles of epithelial cell subpopulations from parous mice and age-matched virgin control mice were similar (Figure 1C, D; Additional file 1) with the exception of luminal Sca1+ cells, which decreased by about 50% in parous mice (P = 0.02). Of note, luminal Sca1+ cells have been shown to be enriched for hormone receptor-positive cells [17]. This was verified by qPCR (see Additional file 6). These data demonstrate that the adopted experimental procedure permits the isolation of all epithelial cell subpopulations from parous mice, including basal stem/progenitor cells, at levels adequate for transcriptomic and functional analyses.

Parity upregulates differentiation genes in all cell subpopulations and decreases the Wnt/Notch signaling ratio in the basal stem/progenitor cell subpopulation

Next we investigated the effects of early pregnancy on the gene-expression profiles of the isolated mammary epithelial cell subpopulations. To control for the effect of the cell-isolation procedure on gene expression, we performed first a transcriptome and cluster analysis in non-FACS-sorted total mammary cell suspensions from age-matched virgin and parous mice (Figure 2A). The analysis showed that pregnancy induces an upregulation of many immunoglobulin and differentiation genes (for example, whey acidic protein Wap, and carboxyl ester lipase Cel) and a downregulation of growth factors (for example, amphiregulin Areg, betacellulin Btc, tumor-associated calcium signal transducer 2 Tacstd2), and extracellular matrix (ECM) elements (for example, laminin, gamma 2 Lamc2, desmocollin 2 Dsc2). These data are consistent with the published pregnancy-induced gene signature determined in snap-frozen rodent mammary glands [9, 10]. Thus, the gene signature was not lost during the isolation procedure, confirming the validity of our experimental system. Furthermore and most interestingly, a novel 3.4-fold downregulation of the Wnt signaling protein Wnt4 was observed in cells from parous mice (Figure 2). These microarray data were validated by qPCR for four genes, including Wnt4. Importantly, the expression of the luminal marker Krt8 was not altered on parity, demonstrating that the observed changes in gene expression were independent of unspecific alterations in total luminal cell numbers (Figure 2B).

Because the number of isolated basal stem/progenitor cells for gene-profiling analysis was limited compared with the other cell subpopulations (that is, myoepithelial, luminal Sca1-, and luminal Sca1+ cells), it was important to evaluate the influence of different cell numbers on transcriptome analysis. As shown in Additional Material, we found that it is valid to use cell numbers in the range of 2,000 to 50,000 cells for comparison of transcriptomes from different mammary cell subpopulations (Additional file 6). Furthermore, qPCR analysis of nonamplified and amplified cDNA for known basal (Krt14) and luminal (Krt8; Krt19) cell marker genes confirmed that the amplification process of the microarray analysis was unbiased (Additional file 6).

In the subsequent transcriptome analysis of FACS-sorted mammary epithelial cell subpopulations from age-matched virgin control and parous mice, all mammary epithelial cell subpopulations, except immune cell-depleted stromal cells, showed parity-induced changes in gene expression with by far the most prominent effects being scored for basal stem/progenitor cells (Figure 3A). Furthermore, although differentiation genes were upregulated in all epithelial cell subpopulations, the strongest prodifferentiation effects were seen in basal stem/progenitor cells from parous mice. For example, casein alpha 1 (Csn1s1), casein beta (Csn2), and lactotransferrin (Ltf) were upregulated 41-, 19-, and 14-fold, respectively (Figure 3B). Apart from these differentiation genes, the Notch co-activator Maml2 was found among the 10 most-upregulated genes (by fold change and/or by P value) in basal stem/progenitor cells from parous mice. In contrast, the 10 most-downregulated genes included the Wnt target and co-receptor Lgr5, the Wnt target and epithelial stem cell marker keratin 15 (Krt15) [33], and the Wnt targets versican (Vcan) and Igfbp3 [34] (Figure 3B). Extension of the analysis to all data for signaling-pathway genes revealed further Wnt target genes, which were downregulated (for example, CD44 and follistatin (Fst) [35, 36]), and the Wnt inhibitor Sfrp1, which was upregulated in basal stem/progenitor cells from parous mice. Moreover, in the same epithelial cell subpopulation of parous mice, the Notch ligand Jag2 and the Notch target genes Gata3, Id3, and Dusp1 were upregulated, whereas the Notch inhibitor Itch and the Maml2 RBP-J complex inhibitors Cdk8 and Pak1 [37] were downregulated (Figure 3B). These gene-expression profiling data were validated with qPCR for several selected genes (Figure 3C). Furthermore, one of the classic Wnt target genes (versican) was examined on the protein level and found to be strongly downregulated in the basal compartment of mammary glands from parous mice (Figure 3D). Thus, the data demonstrate strong upregulation of differentiation genes, a downregulation of Wnt target genes, and an increase in Notch signaling in the basal stem/progenitor-cell subpopulation of parous mice.

Further verification of these conclusions came from Motif Activity Response Analysis (MARA) of transcription factor activities (see Methods) [38]. The canonical Wnt transcription factor LEF1/TCF7 was shown to have significantly decreased activity in basal stem/progenitor cells from parous mice (Figure 4A), which is consistent with the observed decrease in canonical Wnt signaling in the previous gene-expression profile analysis.

Final confirmation of parity-induced downregulation of Wnt signaling in mammary glands was provided by immunohistochemical staining of mammary gland sections from age-matched virgin control and parous mice for β-catenin. As illustrated in Figure 4B, nuclear β-catenin was observed in basal but not in luminal mammary epithelial cells of virgin mice. In parous mice, the proportion of basal mammary epithelial cells positive for nuclear β-catenin was significantly decreased (Figure 4C). Because nuclear β-catenin is necessary for Wnt target gene expression, this finding represents an additional verification of parity-induced downregulation of Wnt signaling (Figure 3).

The parity-induced decrease in Wnt/Notch signaling ratio is specific for the basal stem/progenitor cell subpopulation

To assess whether the observed decrease in canonical Wnt and increase in Notch signaling were specific for basal stem/progenitor cells, we next investigated the enrichment of Wnt/Notch signaling genes over all genes altered in the various FACS-sorted mammary epithelial cell subpopulations from parous as compared with virgin control mice. Such analysis showed that canonical Wnt signaling was significantly downregulated in basal stem/progenitor but not in luminal Sca1- or luminal Sca1+ cells from parous mice (Table 1; see Additional file 7). In myoepithelial cells, a trend toward a decrease in Wnt signaling was not significant when applying a very stringent cut-off (see Table 1 legend) and was probably due to contamination of the myoepithelial cell subpopulation with basal stem and/or progenitor cells. Similarly, Notch signaling was found to be significantly upregulated in basal stem/progenitor cells but in no other mammary epithelial cell subpopulation (Table 1). Whereas downregulation of Wnt and upregulation of Notch signaling were specific for basal stem/progenitor cells, similar enrichment analyses for genes involved in other signaling pathways revealed that the previously reported upregulation of the p53-p21 pathway [39] occurred in all epithelial cell subpopulations tested (Additional file 7).

Further bioinformatics analysis of the data with different software [40] provided a second line of verification of the specific decrease in Wnt signaling in basal stem/progenitor cells after parity. Performing transcription-factor activity analyses based on target gene expression, we found TCF3, an inhibitor of the canonical Wnt signaling pathway [41], to be the transcription factor with the highest z-score (Z = 3.521) and the protooncogene MYC to have a very low z-score (Z = -4.108) in basal stem/progenitor cells but not in other mammary epithelial cell subpopulations of parous mice (Figure 4D, z-score defined as (x-μ)/sd and used as a measure for transcription factor activity based on the expression levels of target genes). These findings confirm the results obtained in the MARA analysis. Furthermore, inhibition of MYC leads to upregulation of Sfrp1, which in turn inhibits canonical Wnt signaling [42]. Upregulation of Sfrp1 was observed in the microarray analysis and validated with qPCR (Figure 3B, C), thus directly reflecting the expected effects of the bioinformatic predictions.

Downregulation of canonical Wnt signaling and MYC activity would be expected to decrease the propensity for cancer, and indeed, in an analysis of biofunctions, a marked and consistent decrease in cancer-associated functions, was observed for basal stem/progenitor cells (Figure 4E), but no other mammary epithelial cell subpopulation. This potential anticancer phenotype of basal stem/progenitor cells was underscored by gene-enrichment analyses on all available pathway gene sets (see Additional file 8), which showed a strong downregulation of proliferation- and tumorigenesis-associated gene sets.

Parity decreases the in vitro clonogenic potential to the greatest extent in the basal stem/progenitor cell subpopulation

Because decreased Wnt signaling and increased Notch signaling have been shown to decrease in vitro and in vivo proliferation of basal stem/progenitor cells [43, 44], we next assessed the in vitro colony-formation capacities of mammary epithelial cell subpopulations from parous and age-matched virgin mice (Figure 5A). In virgin control mice, luminal Sca1- cells had the highest colony-formation capacity, with an average of 107 colonies per well (Figure 5B). This strong clonogenic potential suggests a pronounced progenitor identity of luminal Sca1- cells and is consistent with previous observations in younger virgin mice [16, 17]. A high colony-formation capacity was also observed for the basal stem/progenitor cells of virgin mice (63 colonies per well) (Figure 5B), which is consistent with the notion that CD49fHigh cells contain a high proportion of basal progenitor cells as well as putative mammary stem cells [45]. With the exception of luminal Sca1- cells, the colony-formation capacities of all epithelial cell subpopulations were lower in parous mice than in age-matched virgins (Figure 5B). Thereby, by far the most pronounced difference was observed in basal stem/progenitor cells (Figure 5B). A substantial decrease in the colony-formation capacity was also seen for the myoepithelial cell subpopulation, which also contains basal progenitor cells. Of note, basal stem/progenitor cells from parous mice did not die but remained as quiescent single cells or divided only once during 5 days of culture (Figure 5C).

Figure 5
figure 5

Parity reduces the progenitor potential of mammary epithelial cell subpopulations. (A) Representative images of individual wells with colonies formed by the specified cell subpopulations from age-matched virgin and parous mice. Scale bar, 4 mm. (B) Bar graph comparing the colony-forming capacities of myoepithelial cells, basal stem/progenitor cells, luminal Sca1-, and luminal Sca1+ cells of age-matched virgin and parous mice. Data are from three independent experiments and represent the mean ± SEM of colonies per well; 18 wells were assessed per cell type. *P < 0.015, NS: not significant (P = 0.08), by using two-tailed unpaired Student t test. (C) Representative images for the immunophenotyping of 5-day-old colonies grown from myoepithelial cells, basal stem/progenitor cells, luminal Sca1-, and luminal Sca1+ cells. The colonies were stained for luminal Krt18 (green) and basal Krt14 (red) expression. Hoechst 33342 (blue) was used to distinguish nuclei and to label feeder cells (negative control). Scale bar, 200 μm.

The reduced progenitor potential of basal mammary epithelial cell subtypes and especially of basal stem/progenitor cells from parous mice was further confirmed by phenotypic analysis of colonies by colony size (cell number) and immunophenotyping of luminal and basal markers Krt18 and Krt14, respectively. The largest difference in colony size was observed for the basal stem/progenitor cell subpopulation, where parity induced a decrease in colonies of ≥ 20 cells from 66% to 9% (Table 2A; Figure 5C). Considerably smaller parity-induced reductions in colony size were observed for the other epithelial cell subpopulations, although a substantial decrease was also seen for basal myoepithelial cells (Table 2A). With regard to Krt18/Krt14 double positivity, parity induced a reduction in Krt18/Krt14 double-positive colonies derived from basal stem/progenitor cells from 26% to 0 (Table 2A). In contrast, although Krt18/Krt14 double-positivity decreased by 27%, 57% of colonies derived from luminal Sca1- cells of parous mice maintained double positivity (Table 2A). No effects of parity on Krt18/Krt14 double positivity were observed for colonies derived from luminal Sca1+ or basal myoepithelial cells (Table 2A), whereby myoepithelial cells did not give rise to double-positive colonies even when originating from virgin control mice. Thus, consistent with the basal mammary stem/progenitor cell subpopulation-specific reduction in the Wnt/Notch signaling ratio, parity decreased colony-formation capacity, in vitro proliferation potential (colony number and size), and Krt18/Krt14 double positivity most prominently in the basal stem/progenitor cell subpopulation.

Table 2 Functional characterization of isolated epithelial cell subpopulations

Parity decreases the in vivo reconstitution efficiency of the basal stem/progenitor cell subpopulation

To test their proliferation potential in vivo, we transplanted basal stem/progenitor cells into deepithelialized mammary glands ("cleared fat pads") [12]. It already was demonstrated that basal CD49fHigh cells have the highest mammary gland reconstitution ability of all mammary cell subpopulations in virgin mice [16, 19]. Transplantation of FACS-isolated basal stem/progenitor cells into cleared fat pads demonstrated a significant decrease in the number of large outgrowths (≥ 10% and ≥ 25% of fat pad filled), indicating a decrease in the in vivo proliferation potential. Interestingly, when assessing also for rudimentary outgrowths (≥ 3% of fat pad filled), no significant difference in the number of outgrowths was observed between parous and virgin donors (Table 2B). Apart from the change in size, no qualitative differences were apparent between outgrowths from virgin and parous donors. In both cases, ductal as well as lobular structures were formed. Hence, early parity led to a reduction in in vivo mammary repopulating efficiency of FACS-isolated basal stem/progenitor cells, whereas their ability to differentiate into different mammary epithelial structures was maintained.

Decrease in the proportion of Wnt4-secreting cells after early parity can explain decreased Wnt signaling and reduced proliferation capacity in basal mammary stem/progenitor cells

Having observed most of the expected functional consequences of a decrease in the Wnt/Notch signaling ratio in basal stem/progenitor cells from parous mice, we finally examined the possible cause of parity-induced reduction in canonical Wnt signaling and proliferation capacity in basal stem/progenitor cells. Because parity induced a greater than threefold decrease in Wnt ligand Wnt4 gene expression (Figure 2A, B), and Wnt4 is known to be secreted in response to progesterone by hormone-sensing luminal cells [46], thus inducing canonical Wnt signaling in mammary stem/progenitor cells [47], a parity-induced decrease in estrogen/progesterone-sensitive luminal cells could explain the overall decrease in Wnt signaling in mammary stem cells. This hypothesis is supported by the reduction in the proportion of luminal Sca1+ cells isolated from parous mice (Figure 1D) and by the demonstration that luminal Sca1+ cells are hormone receptor positive (Additional file 6) [17]. Furthermore, immunohistochemical analysis of mammary gland sections for estrogen receptor alpha (ERα) and its target progesterone receptor (PR) showed a twofold decrease in ERα- and a threefold decrease in PR-positive cells in parous compared with age-matched virgin control mice (Figure 6A, B). These results were additionally verified by qPCR in total mammary cell suspensions (Figure 6C). Notably, expression of the luminal marker Krt8 was similar in cell suspensions from parous and age-matched virgin control mice, supporting the conclusion of a specific decrease in hormone receptor-positive cells rather than a general cell loss after pregnancy. Furthermore, parity-induced reduction in progesterone-stimulated Wnt4 expression was independent of blood progesterone, because average blood progesterone concentrations were similar in parous mice and age-matched virgin control mice in estrus (see Additional file 9). Finally, supplementation of the culture medium with recombinant Wnt4 stimulated in vitro proliferation capacity of basal myoepithelial and basal stem/progenitor cells from parous mice by +138% ± 22% and +140% ± 17%, respectively (Figure 7A). In contrast, no significant effects of recombinant Wnt4 on the colony-formation capacity of luminal Sca1- (+3.7% ± 3.0%) and luminal Sca1+ cells from parous mice (-5.1% ± 2.9%) were observed. These data strongly suggest a causal relation between reduced number of luminal progesterone receptor-positive/Wnt4-secreting cells and decreased Wnt/Notch signaling and proliferation potential of basal stem/progenitor cells after early pregnancy (Figure 7B).

Figure 6
figure 6

Early pregnancy decreases the proportion of progesterone receptor (PR) and estrogen receptor α (ERα)-positive cells. (A) Representative images of immunostaining for PR and ERα in mammary gland sections from age-matched virgin and parous mice in estrus. Scale bar, 50 μm. (B) Bar graph comparing the relative frequency of estrogen- and progesterone receptor-positive cells between mammary glands of virgin and parous mice. Data represent the mean ± SD (virgin mice, n = 6; parous mice, n = 5). For PR, P = 3.70E-07; for ERα, P = 0.003, by using two-tailed unpaired Student t test. (C) qPCR for progesterone receptor (Pgr), estrogen receptor alpha (Esr1), and the luminal marker keratin 8 (Krt8) genes in total mammary cell suspensions. Data are expressed as the mean ± SEM from four groups of a minimum of five parous and five age-matched virgin control mice.

Figure 7
figure 7

Effect of Wnt4 on the proliferation capacity of basal stem/progenitor cells after early pregnancy in mice. (A) Recombinant Wnt4 rescues the parity-induced in vitro proliferation defect in basal mammary epithelial cells. Selected mammary epithelial cells from parous mice were cultured in the absence or presence of recombinant Wnt4. Three independent experiments were performed. Data represent the mean ± SEM of colonies per well, with six to nine wells assessed per cell type. *P ≤ 0.02 (two-tailed unpaired Student t test). (B) Mechanistic model illustrating the parity-induced decrease in hormone-sensing and Wnt4-secreting luminal cells on the Wnt/Notch signaling pathways and the proliferation/differentiation potential in basal stem/progenitor cells.

In conclusion, all collected data indicate that early parity induces a decrease in luminal hormone-sensing ERα- and PR-positive cells, which leads to decreased Wnt4 expression levels and to reduced Wnt signaling in basal stem/progenitor cells. Consistent with the decrease in the Wnt/Notch signaling ratio, differentiation is promoted and proliferation inhibited in basal stem/progenitor cells of parous mice.

Discussion

This study demonstrates that a history of early pregnancy changes the gene-expression profiles and functional properties of mammary epithelial cell subpopulations in a cell subtype-specific fashion. Most important, the following parity-induced alterations were observed in mice: (a) an induction of differentiation and downregulation of the Wnt/Notch signaling ratio in basal stem/progenitor cells; (b) a decrease in the in vitro and in vivo proliferation potential of isolated basal stem/progenitor cells; (c) a selective downregulation of potentially tumorigenic biofunctions in the basal stem/progenitor cell subpopulation; (d) a reduction in estrogen- and progesterone-responsive and Wnt4-secreting luminal cells; and (e) a rescue of the proliferation defect in basal stem/progenitor cells in vitro by recombinant Wnt4. The finding of a decreased Wnt/Notch signaling ratio provides direct experimental evidence for the hypothesis that early pregnancy changes the "genomic signature" of mammary stem/progenitor cells [7], causing their differentiation and reducing their proliferation potential. Furthermore, the data indicate a novel causal relation between parity-induced reduction in hormone-sensing and Wnt4-secreting luminal cells and altered biofunctions in basal stem/progenitor cells.

The basal CD49fHigh cells, as isolated in this study, are a subfraction of basal epithelial cells [19, 21]. Basal CD49fHigh cells have been demonstrated previously to be enriched in mammary repopulating units (MRUs) (known as mammary stem cells (MaSCs)) [16, 19] and to correspond to the Lin- CD24+ CD29High epithelial cell subpopulation isolated by an alternative method [18, 19, 45]. However, the isolated basal CD49fHigh epithelial cell subpopulation represents a heterogeneous cell fraction containing, in addition to MaSCs, basal progenitor cells and possibly mature myoepithelial cells [45]. Progenitor cells can be characterized in vitro by their colony-formation capacity [22, 23], whereas MaSCs have traditionally been defined by their in vivo regenerative capacity [16–18]. The observed effects of parity on the in vitro and in vivo proliferation capacity of the CD49fHigh cell subpopulation (Table 2) suggests that both basal progenitor cells and basal MaSCs are the target of an early pregnancy within the mammary epithelium. Although the dramatic decrease in the in vitro proliferation capacity of the CD49fHigh cell subpopulation (Figure 5) indicates a predominant effect of an early pregnancy on basal progenitor cells, the additional reduction in large in vivo outgrowths (≥ 10% of fat pad filled) and the prevalence of rudimentary outgrowths (≥ 3% of fat pad filled) after parity (Table 2) suggest that isolated basal MaSCs are also affected by an early pregnancy. Because basal mammary stem and progenitor cells are closely related and likely to be interdependent in their proliferation potentials, our data do not permit a definite discrimination between basal stem and progenitor cells as primary targets of pregnancy. Therefore, we adhered to the combined term basal stem/progenitor cells throughout this study.

We found the p53-p21 pathway to be enriched to a similar extent in all mammary epithelial cell subpopulations (Additional file 7) [39], and hence, although parity-caused induction of the p53-p21 pathway may explain the relatively modest decrease in in vitro colony-formation potential of the luminal Sca1+ cell subpopulation (Figure 5), it may contribute but cannot account for the almost complete proliferation block in basal stem/progenitor cells. The most prominent parity-induced alterations in gene expression in basal stem/progenitor cells were downregulation of the Wnt-signaling pathway, upregulation of the Notch-signaling pathway, and upregulation of differentiation genes. Decreased Wnt signaling in basal stem/progenitor cells from parous mice was verified on the protein level by measuring versican and nuclear β-catenin expression (Figures 3D and 4B, C). Wnt signaling has been shown to promote long-term expansion of cultured Lin- CD24+ CD29High cells and to provide a competitive advantage in mammary gland reconstitution assays [43]. The latter is especially true for the expression of the classic Wnt target Lgr5 [48], which was found in our study to be downregulated in the basal stem/progenitor cell subpopulation after parity. Notably, Wnt signaling inhibition was demonstrated to have an antiproliferative effect in CD29High cells [43]. Furthermore, Notch signaling was observed to reduce in vitro and in vivo proliferation of CD29High cells while promoting their differentiation [44, 49]. Because the differentiation processes of stem and progenitor cells in many organs and in several model systems are dependent on the Wnt/Notch signaling ratio [50], an overall reduction in Wnt/Notch signaling ratio would be expected to have a dramatic antiproliferation and prodifferentiation effect in mammary basal stem/progenitor cells. This is exactly what we observed in the CD49fHigh cell subpopulation of parous mice. Thereby, the overall conclusion of our study is strengthened by the fact that all assays used (that is, transcriptome analysis, bioinformatics transcription factor, and gene-enrichment analyses, in vitro colony-forming assay, in vivo transplantation assay, and immunohistochemistry) pointed into the same direction. Thus, analysis of specific mammary epithelial cell subpopulations allowed the discovery of a decrease in the Wnt/Notch signaling ratio, which so far is the only plausible explanation for the observed differentiation burst and dramatic proliferation block experienced by basal stem/progenitor cells of parous mice.

As possible explanations for the parity-induced decrease in Wnt signaling, we found a marked increase in the activity of the Wnt repressor TCF3 in basal stem/progenitor cells and a more than threefold reduction in expression of the secreted Wnt ligand Wnt4 in total mammary cells from parous mice. The latter corresponded to a similar decrease in Wnt4-secreting [46] and estrogen and progesterone receptor-positive luminal cells. Notably, a similar decrease in progesterone receptor α-positive cells after parity has also been observed in human breast epithelium [51]. Hence, mechanistically, early parity decreases the hormone responsiveness of the mammary gland in mice by decreasing the number of estrogen/progesterone receptor-positive luminal cells. This reduces the paracrine signaling cascade mediated by Wnt4, inducing TCF3-dependent repression [52] and/or primary downregulation [47] of canonical Wnt signaling and secondary (reactive) upregulation of Notch signaling in basal stem/progenitor cells.

As a final consequence, proliferation is repressed and basal stem/progenitor cells differentiate. This mechanistic model is supported by the ability of recombinant Wnt4 to rescue the proliferation defect of basal stem/progenitor cells from parous mice in vitro (Figure 7).

Our findings in specific mammary epithelial cell subpopulations are in part consistent with and in part contradictory to studies in entire breasts/mammary glands, total mammary cells, or total mammary epithelial cells. With respect to the transcriptome analysis, our studies in total mammary cells agree with previous reports in the entire mammary glands [9, 10]. However, in intact mammary glands or total mammary cells, the additional presence of stromal and dominant epithelial cell subtypes might mask the detection of key signaling-pathway changes. Indeed, our study demonstrates that isolation of specific mammary epithelial cell subpopulations is a prerequisite for the detection of a decrease in Wnt/Notch signaling ratio in basal stem/progenitor cells. A similar masking effect by stromal and dominant epithelial cell subtypes (for example, strong clonogenic luminal Sca1- cells) might also explain why a previous study with a similar early-pregnancy protocol did not observe a parity-induced reduction in in vitro proliferation of total mixed mammary cells [20]. Controversial results have also been reported with respect to the effect of parity on the in vivo mammary-repopulating capacity. Hence, although Britt et al. [21] found no effects of late pregnancy (9 weeks) on mammary-repopulating units (MRUs) in total mammary epithelial cells, Siwko et al. [20] (early-pregnancy protocol) observed a parity-induced reduction in the mammary-repopulating capacity of total mixed mammary cells. With a similar early-pregnancy protocol and the same cut-off for mammary gland outgrowth (≥ 10% of fat pad filled), our findings in isolated mammary basal stem/progenitor cells appear consistent with the observations of Siwko et al. However, given that our studies were performed with isolated mammary basal stem/progenitor cell subpopulations, our findings are not directly comparable with and neither confirm nor contradict previous studies using total mammary (epithelial) cells [20, 21]. Furthermore, the fact that the number of smaller outgrowths (≥ 3% of fat pad filled) was unchanged after parity suggests that MRUs survive after pregnancy despite their reduced reconstitution efficiency. This conclusion is in line with the recent demonstration that Wnt-responsive mammary epithelial stem cells persist after parity in Axin2 reporter mice [15].

It is intriguing to speculate that marked growth inhibitory effects and the downregulation of canonical Wnt signaling in basal stem/progenitor cells account, at least in part, for the cancer-protective effect of early pregnancy. Increases in canonical Wnt signaling have been linked repeatedly to oncogenesis [53, 54]. Moreover, downregulation of the Wnt inhibitory protein Sfrp1 and overexpression of the Wnt target versican have been associated with carcinogenesis [55, 56]. In the transcriptome analysis reported here, Sfrp1 was upregulated and versican downregulated, thus supporting a parity-induced anticarcinogenic effect. Moreover, Li et al. [57] showed that transgenes encoding components of the Wnt signaling pathway preferentially induce mammary cancers from progenitor cells. Hence, the contribution of decreased canonical Wnt signaling in basal stem/progenitor cells to the cancer-protective effect of early pregnancy may be in conjunction with other tumor-suppressing mechanisms, such as parity-induced induction of p53 [58, 59]. Thereby, the increase in the TCF3 repressor activity in basal stem/progenitor cells is expected to elevate the threshold further for the activation of tumorigenic Wnt signaling [52]. In addition to decreased Wnt signaling in basal stem/progenitor cells, decreased ERα- and PR-positive cells could also be a mechanism for the breast cancer-protective effect of an early pregnancy. This is especially relevant, given the specific protective effect of pregnancy against ER/PR-positive tumors [60].

Furthermore, mammary epithelial cell differentiation per se has been suggested to exert a breast cancer-protective effect. This has been challenged, however, by the observation that differentiation-causing agents such as placental lactogen and perphenazine failed to protect against carcinogenesis in rodents [3, 61]. Also, the hypothesis of a potential breast cancer-protective effect of mammary stromal cells [62] is not supported by our study, because stromal cells exhibited by far the fewest parity-induced gene-expression changes. However, our stromal cell subpopulation was not homogenous and devoid of immune cells, which may have masked some parity-induced alterations. In any case, the parity-induced downregulation of the Wnt/Notch signaling ratio in basal stem/progenitor cells represents a possible important mechanism for the breast cancer-protective effect of early pregnancy.

Conclusions

This study identified downregulation of the Wnt/Notch signaling ratio in basal stem/progenitor cells as the dominant early parity-induced alteration of gene expression in mice. This change in gene expression is specific for basal mammary stem/progenitor cells, is associated with proliferation defects in vitro and in vivo, and is probably caused by an early parity-induced decrease in hormone-sensitive and Wnt4-secreting luminal cells. Importantly, because a similar reduction in progesterone receptor α-positive luminal cells has been reported in women [51], parity-induced alterations in Wnt/Notch signaling pathways may also occur in human basal stem/progenitor cells. Testing whether Wnt inhibitors mimic early parity-induced breast cancer protection warrants further investigation.

Abbreviations

ER:

estrogen receptor

FDR:

false discovery rate

FWER:

family-wise error rate

GSEA:

gene-set enrichment analysis

MARA:

motif activity response analysis

MaSC:

mammary stem cell

MRU:

mammary repopulating unit

NES:

normalized enrichment score

PR:

progesterone receptor

RBC:

red blood cell

WBC:

white blood cell.

References

  1. MacMahon B, Cole P, Lin TM, Lowe CR, Mirra AP, Ravnihar B, Salber EJ, Valaoras VG, Yuasa S: Age at first birth and breast cancer risk. Bull World Health Org. 1970, 43: 209-221.

    CAS  PubMed  PubMed Central  Google Scholar 

  2. Albrektsen G, Heuch I, Hansen S, Kvale G: Breast cancer risk by age at birth, time since birth and time intervals between births: exploring interaction effects. Br J Cancer. 2005, 92: 167-175. 10.1038/sj.bjc.6602302.

    Article  CAS  PubMed  Google Scholar 

  3. Medina D: Mammary developmental fate and breast cancer risk. Endocr Relat Cancer. 2005, 12: 483-495. 10.1677/erc.1.00804.

    Article  CAS  PubMed  Google Scholar 

  4. Rajkumar L, Kittrell FS, Guzman RC, Brown PH, Nandi S, Medina D: Hormone-induced protection of mammary tumorigenesis in genetically engineered mouse models. Breast Cancer Res. 2007, 9: R12-10.1186/bcr1645.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Thordarson G, Jin E, Guzman RC, Swanson SM, Nandi S, Talamantes F: Refractoriness to mammary tumorigenesis in parous rats: is it caused by persistent changes in the hormonal environment or permanent biochemical alterations in the mammary epithelia?. Carcinogenesis. 1995, 16: 2847-2853. 10.1093/carcin/16.11.2847.

    Article  CAS  PubMed  Google Scholar 

  6. Schedin P, Mitrenga T, McDaniel S, Kaeck M: Mammary ECM composition and function are altered by reproductive state. Mol Carcinog. 2004, 41: 207-220. 10.1002/mc.20058.

    Article  CAS  PubMed  Google Scholar 

  7. Russo J, Moral R, Balogh GA, Mailo D, Russo IH: The protective role of pregnancy in breast cancer. Breast Cancer Res. 2005, 7: 131-142. 10.1186/bcr1029.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Russo J, Balogh GA, Russo IH: Full-term pregnancy induces a specific genomic signature in the human breast. Cancer Epidemiol Biomarkers Prev. 2008, 17: 51-66. 10.1158/1055-9965.EPI-07-0678.

    Article  CAS  PubMed  Google Scholar 

  9. D'Cruz CM, Moody SE, Master SR, Hartman JL, Keiper EA, Imielinski MB, Cox JD, Wang JY, Ha SI, Keister BA, Chodosh LA: Persistent parity-induced changes in growth factors, TGF-beta3, and differentiation in the rodent mammary gland. Mol Endocrinol. 2002, 16: 2034-2051. 10.1210/me.2002-0073.

    Article  PubMed  Google Scholar 

  10. Blakely CM, Stoddard AJ, Belka GK, Dugan KD, Notarfrancesco KL, Moody SE, D'Cruz CM, Chodosh LA: Hormone-induced protection against mammary tumorigenesis is conserved in multiple rat strains and identifies a core gene expression signature induced by pregnancy. Cancer Res. 2006, 66: 6421-6431. 10.1158/0008-5472.CAN-05-4235.

    Article  CAS  PubMed  Google Scholar 

  11. Lindeman GJ, Visvader JE: Insights into the cell of origin in breast cancer and breast cancer stem cells. Asia Pac J Clin Oncol. 2010, 6: 89-97. 10.1111/j.1743-7563.2010.01279.x.

    Article  PubMed  Google Scholar 

  12. Deome KB, Faulkin LJ, Bern HA, Blair PB: Development of mammary tumors from hyperplastic alveolar nodules transplanted into gland-free mammary fat pads of female C3H mice. Cancer Res. 1959, 19: 515-520.

    CAS  PubMed  Google Scholar 

  13. Bruno RD, Smith GH: Functional characterization of stem cell activity in the mouse mammary gland. Stem Cell Rev. 2011, 7: 238-247. 10.1007/s12015-010-9191-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Van Keymeulen A, Rocha AS, Ousset M, Beck B, Bouvencourt G, Rock J, Sharma N, Dekoninck S, Blanpain C: Distinct stem cells contribute to mammary gland development and maintenance. Nature. 2011, 479: 189-193. 10.1038/nature10573.

    Article  CAS  PubMed  Google Scholar 

  15. van Amerongen R, Bowman AN, Nusse R: Developmental stage and time dictate the fate of Wnt/beta-catenin-responsive stem cells in the mammary gland. Cell Stem Cell. 2012, 11: 387-400. 10.1016/j.stem.2012.05.023.

    Article  CAS  PubMed  Google Scholar 

  16. Stingl J, Eirew P, Ricketson I, Shackleton M, Vaillant F, Choi D, Li HI, Eaves CJ: Purification and unique properties of mammary epithelial stem cells. Nature. 2006, 439: 993-997.

    CAS  PubMed  Google Scholar 

  17. Sleeman KE, Kendrick H, Robertson D, Isacke CM, Ashworth A, Smalley MJ: Dissociation of estrogen receptor expression and in vivo stem cell activity in the mammary gland. J Cell Biol. 2007, 176: 19-26. 10.1083/jcb.200604065.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Shackleton M, Vaillant F, Simpson KJ, Stingl J, Smyth GK, Asselin-Labat ML, Wu L, Lindeman GJ, Visvader JE: Generation of a functional mammary gland from a single stem cell. Nature. 2006, 439: 84-88. 10.1038/nature04372.

    Article  CAS  PubMed  Google Scholar 

  19. Smalley MJ, Kendrick H, Sheridan JM, Regan JL, Prater MD, Lindeman GJ, Watson CJ, Visvader JE, Stingl J: Isolation of mouse mammary epithelial subpopulations: a comparison of leading methods. J Mammary Gland Biol Neoplasia. 2012, 17: 91-97. 10.1007/s10911-012-9257-1.

    Article  PubMed  Google Scholar 

  20. Siwko SK, Dong J, Lewis MT, Liu H, Hilsenbeck SG, Li Y: Evidence that an early pregnancy causes a persistent decrease in the number of functional mammary epithelial stem cells: implications for pregnancy-induced protection against breast cancer. Stem Cells. 2008, 26: 3205-3209. 10.1634/stemcells.2008-0103.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Britt KL, Kendrick H, Regan JL, Molyneux G, Magnay FA, Ashworth A, Smalley MJ: Pregnancy in the mature adult mouse does not alter the proportion of mammary epithelial stem/progenitor cells. Breast Cancer Res. 2009, 11: R20-10.1186/bcr2245.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Smalley MJ, Titley J, O'Hare MJ: Clonal characterization of mouse mammary luminal epithelial and myoepithelial cells separated by fluorescence-activated cell sorting. Vitro Cell Dev Biol Anim. 1998, 34: 711-721. 10.1007/s11626-998-0067-0.

    Article  CAS  Google Scholar 

  23. Stingl J, Eaves CJ, Zandieh I, Emerman JT: Characterization of bipotent mammary epithelial progenitor cells in normal adult human breast tissue. Breast Cancer Res Treat. 2001, 67: 93-109. 10.1023/A:1010615124301.

    Article  CAS  PubMed  Google Scholar 

  24. Roarty K, Rosen JM: Wnt and mammary stem cells: hormones cannot fly wingless. Curr Opin Pharmacol. 2010, 10: 643-649. 10.1016/j.coph.2010.07.004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Many AM, Brown AM: Mammary stem cells and cancer: roles of Wnt signaling in plain view. Breast Cancer Res. 2010, 12: 313-

    PubMed  PubMed Central  Google Scholar 

  26. Sleeman KE, Kendrick H, Ashworth A, Isacke CM, Smalley MJ: CD24 staining of mouse mammary gland cells defines luminal epithelial, myoepithelial/basal and non-epithelial cells. Breast Cancer Res. 2006, 8: R7-10.1186/bcr1371.

    Article  PubMed  Google Scholar 

  27. Smalley MJ: Isolation, culture and analysis of mouse mammary epithelial cells. Methods Mol Biol. 2010, 633: 139-170. 10.1007/978-1-59745-019-5_11.

    Article  CAS  PubMed  Google Scholar 

  28. GSEA homepage. [http://www.broadinstitute.org/gsea]

  29. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Edgar R, Domrachev M, Lash AE: Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Gonzalez RC, Woods RE, Eddins SL: Digital Image Processing using MATLAB. 2009, Knoxville, TN: Gatesmark Publishing

    Google Scholar 

  32. ISMARA. [http://www.mara.unibas.ch]

  33. Palmer HG, Anjos-Afonso F, Carmeliet G, Takeda H, Watt FM: The vitamin D receptor is a Wnt effector that controls hair follicle differentiation and specifies tumor type in adult epidermis. PLoS One. 2008, 3: e1483-10.1371/journal.pone.0001483.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Oikonomopoulos A, Sereti KI, Conyers F, Bauer M, Liao A, Guan J, Crapps D, Han JK, Dong H, Bayomy AF, Fine GC, Westerman K, Biechele TL, Moon RT, Force T, Liao R: Wnt signaling exerts an antiproliferative effect on adult cardiac progenitor cells through IGFBP3. Circ Res. 2011, 109: 1363-1374. 10.1161/CIRCRESAHA.111.250282.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Yao HH, Matzuk MM, Jorgez CJ, Menke DB, Page DC, Swain A, Capel B: Follistatin operates downstream of Wnt4 in mammalian ovary organogenesis. Dev Dyn. 2004, 230: 210-215. 10.1002/dvdy.20042.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Willert J, Epping M, Pollack JR, Brown PO, Nusse R: A transcriptional response to Wnt protein in human embryonic carcinoma cells. BMC Dev Biol. 2002, 2: 8-10.1186/1471-213X-2-8.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Borggrefe T, Oswald F: The Notch signaling pathway: transcriptional regulation at Notch target genes. Cell Mol Life Sci. 2009, 66: 1631-1646. 10.1007/s00018-009-8668-7.

    Article  CAS  PubMed  Google Scholar 

  38. FANTOM Consortium, Suzuki H, Forrest AR, van Nimwegen E, Daub CO, Balwierz PJ, Irvine KM, Lassmann T, Ravasi T, Hasegawa Y, de Hoon MJ, Katayama S, Schroder K, Carninci P, Tomaru Y, Kanamori-Katayama M, Kubosaki A, Akalin A, Ando Y, Arner E, Asada M, Asahara H, Bailey T, Bajic VB, Bauer D, Beckhouse AG, Bertin N, Björkegren J, Brombacher F, Bulger E, et al: The transcriptional network that controls growth arrest and differentiation in a human myeloid leukemia cell line. Nat Genet. 2009, 41: 553-562. 10.1038/ng.375.

    Article  Google Scholar 

  39. Sivaraman L, Conneely OM, Medina D, O'Malley BW: p53 is a potential mediator of pregnancy and hormone-induced resistance to mammary carcinogenesis. Proc Natl Acad Sci USA. 2001, 98: 12379-12384. 10.1073/pnas.221459098.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Ingenuity Systems. [http://www.ingenuity.com]

  41. Yi F, Pereira L, Hoffman JA, Shy BR, Yuen CM, Liu DR, Merrill BJ: Opposing effects of Tcf3 and Tcf1 control Wnt stimulation of embryonic stem cell self-renewal. Nat Cell Biol. 2011, 13: 762-770. 10.1038/ncb2283.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Cowling VH, Cole MD: Turning the tables: Myc activates Wnt in breast cancer. Cell Cycle. 2007, 6: 2625-2627. 10.4161/cc.6.21.4880.

    Article  CAS  PubMed  Google Scholar 

  43. Zeng YA, Nusse R: Wnt proteins are self-renewal factors for mammary stem cells and promote their long-term expansion in culture. Cell Stem Cell. 2010, 6: 568-577. 10.1016/j.stem.2010.03.020.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Bouras T, Pal B, Vaillant F, Harburg G, Asselin-Labat ML, Oakes SR, Lindeman GJ, Visvader JE: Notch signaling regulates mammary stem cell function and luminal cell-fate commitment. Cell Stem Cell. 2008, 3: 429-441. 10.1016/j.stem.2008.08.001.

    Article  CAS  PubMed  Google Scholar 

  45. Visvader JE, Smith GH: Murine mammary epithelial stem cells: discovery, function, and current status. Cold Spring Harb Perspect Biol. 2011, 3: a004879-10.1101/cshperspect.a004879.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Brisken C, Heineman A, Chavarria T, Elenbaas B, Tan J, Dey SK, McMahon JA, McMahon AP, Weinberg RA: Essential function of Wnt-4 in mammary gland development downstream of progesterone signaling. Genes Dev. 2000, 14: 650-654.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Rajaram RD, Brisken C: Role of canonical Wnt signaling in mammary gland development. PhD Thesis. ISREC, School of Life Sciences Ecole Polytechnique Federal (EPFL). 2009

    Google Scholar 

  48. Plaks V, Brenot A, Lawson DA, Linnemann JR, Van Kappel EC, Wong KC, de Sauvage F, Klein OD, Werb Z: Lgr5-expressing cells are sufficient and necessary for postnatal mammary gland organogenesis. Cell Rep. 2013, 3: 70-78. 10.1016/j.celrep.2012.12.017.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Raouf A, Zhao Y, To K, Stingl J, Delaney A, Barbara M, Iscove N, Jones S, McKinney S, Emerman J, Aparicio S, Marra M, Eaves C: Transcriptome analysis of the normal human mammary cell commitment and differentiation process. Cell Stem Cell. 2008, 3: 109-118. 10.1016/j.stem.2008.05.018.

    Article  CAS  PubMed  Google Scholar 

  50. Munoz-Descalzo S, de Navascues J, Arias AM: Wnt-Notch signalling: an integrated mechanism regulating transitions between cell states. Bioessays. 2012, 34: 110-118. 10.1002/bies.201100102.

    Article  CAS  PubMed  Google Scholar 

  51. Taylor D, Pearce CL, Hovanessian-Larsen L, Downey S, Spicer DV, Bartow S, Pike MC, Wu AH, Hawes D: Progesterone and estrogen receptors in pregnant and premenopausal non-pregnant normal human breast. Breast Cancer Res Treat. 2009, 118: 161-168. 10.1007/s10549-009-0322-4.

    Article  CAS  PubMed  Google Scholar 

  52. Wu CI, Hoffman JA, Shy BR, Ford EM, Fuchs E, Nguyen H, Merrill BJ: Function of Wnt/beta-catenin in counteracting Tcf3 repression through the Tcf3-beta-catenin interaction. Development. 2012, 139: 2118-2129. 10.1242/dev.076067.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Visvader JE, Lindeman GJ: Cancer stem cells: current status and evolving complexities. Cell Stem Cell. 2012, 10: 717-728. 10.1016/j.stem.2012.05.007.

    Article  CAS  PubMed  Google Scholar 

  54. Visvader JE: Keeping abreast of the mammary epithelial hierarchy and breast tumorigenesis. Genes Dev. 2009, 23: 2563-2577. 10.1101/gad.1849509.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Gauger KJ, Hugh JM, Troester MA, Schneider SS: Down-regulation of sfrp1 in a mammary epithelial cell line promotes the development of a cd44high/cd24low population which is invasive and resistant to anoikis. Cancer Cell Int. 2009, 9: 11-10.1186/1475-2867-9-11.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Du WW, Yang BB, Shatseva TA, Yang BL, Deng Z, Shan SW, Lee DY, Seth A, Yee AJ: Versican G3 promotes mouse mammary tumor cell growth, migration, and metastasis by influencing EGF receptor signaling. PLoS One. 2010, 5: e13828-10.1371/journal.pone.0013828.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Li Y, Welm B, Podsypanina K, Huang S, Chamorro M, Zhang X, Rowlands T, Egeblad M, Cowin P, Werb Z, Tan LK, Rosen JM, Varmus HE: Evidence that transgenes encoding components of the Wnt signaling pathway preferentially induce mammary cancers from progenitor cells. Proc Natl Acad Sci USA. 2003, 100: 15853-15858. 10.1073/pnas.2136825100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Medina D, Kittrell FS: p53 function is required for hormone-mediated protection of mouse mammary tumorigenesis. Cancer Res. 2003, 63: 6140-6143.

    CAS  PubMed  Google Scholar 

  59. Dunphy KA, Blackburn AC, Yan H, O'Connell LR, Jerry DJ: Estrogen and progesterone induce persistent increases in p53-dependent apoptosis and suppress mammary tumors in BALB/c-Trp53+/- mice. Breast Cancer Res. 2008, 10: R43-10.1186/bcr2094.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Colditz GA, Rosner BA, Chen WY, Holmes MD, Hankinson SE: Risk factors for breast cancer according to estrogen and progesterone receptor status. J Natl Cancer Inst. 2004, 96: 218-228. 10.1093/jnci/djh025.

    Article  CAS  PubMed  Google Scholar 

  61. Guzman RC, Yang J, Rajkumar L, Thordarson G, Chen X, Nandi S: Hormonal prevention of breast cancer: mimicking the protective effect of pregnancy. Proc Natl Acad Sci USA. 1999, 96: 2520-2525. 10.1073/pnas.96.5.2520.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Smith GH, Salomon DS, Vonderhar BK: Re: "Evidence that an early pregnancy causes a persistent decrease in the number of functional mammary epithelial stem cells-implications for pregnancy-induced protection against breast cancer" by Siwko et al . Stem Cells. 2009, 27: 1223-10.1002/stem.19. author reply 1224-1225

    Article  PubMed  Google Scholar 

  63. Han LQ, Yang GY, Zhu HS, Wang YY, Wang LF, Guo YJ, Lu WF, Li HJ, Wang YL: Selection and use of reference genes in mouse mammary glands. Genet Mol Res. 2010, 9: 449-456. 10.4238/vol9-1gmr724.

    Article  CAS  PubMed  Google Scholar 

  64. Wnt homepage. [http://www.wnt.stanford.edu]

Download references

Acknowledgements

We thank S. Gasser, N. Hynes, M. Smalley, and P. Meier-Abt for helpful comments on the manuscript, C. Rochlitz, D. Schuebeler and D. Medina for helpful discussions, M. Stadler for bioinformatics support, and U. Mueller, S. Bichet, and R. Thierry for technical assistance, as well as further members of the Bentires-Alj lab for advice and discussions. Research in the lab of MB-A is supported by the Novartis Research Foundation, the European Research Council (ERC starting grant 243211-PTPsBDC), the Swiss Cancer League, and the Krebsliga Beider Basel. FM-A is supported by an SNF fellowship.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mohamed Bentires-Alj.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

FM-A conceived of the hypothesis, prepared the mice, designed and performed the experiments, and wrote the manuscript. EM introduced FM-A to the techniques and assisted with some of the experiments. TR oversaw the microarray studies and assisted in the analysis of the microarray data. HB provided technical assistance with the immunohistochemistry. SD provided technical and intellectual support. DM provided technical assistance with the transplantation experiments, and IK performed fat-pad clearings and cell injections. PB and EvN performed the MARA analysis and provided intellectual support. MB-A conceived of the hypothesis, directed the project, and revised the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

13058_2012_3224_MOESM1_ESM.PDF

Additional file 1: Flow-cytometric separation of mammary epithelial cells from virgin control and parous mice. Representative flow-cytometry pseudocolor plots depicting the first steps in the gating strategy used to eliminate doublets, cell clumps, dead cells (DAPI bright) and white blood cells (CD45+) during the procedure for isolating mammary epithelial cell subpopulations from virgin control (A) and parous (B) mice. Subsequent isolation steps are shown in Figure 1. (C) Bar graph showing the proportion of mammary epithelial cells relative to the total white blood cell depleted (CD45-) mammary cells. The apparent modest decrease in total epithelial cells from parous mice was not significant (P = 0.05 by using two-tailed unpaired Student t test). (PDF 589 KB)

13058_2012_3224_MOESM2_ESM.PDF

Additional file 2: Basal stem/progenitor cells from parous mice show reduced in vivoproliferation potential. Examples of an outgrowth (left) and of no outgrowth (right) from basal stem/progenitor cells of age-matched virgin control mice and of parous mice, respectively. (PDF 2 MB)

13058_2012_3224_MOESM3_ESM.PDF

Additional file 3: Custom gene set of Wnt target genes. The gene set is composed of the Wnt targets that have been reported to be upregulated on canonical Wnt signaling in mammalian systems [64]. (PDF 9 KB)

13058_2012_3224_MOESM4_ESM.PDF

Additional file 4: Verification of luminal/basal origin and purity of isolated mammary epithelial cell subpopulations. (A) Immunofluorescent staining of isolated mammary epithelial cells with the luminal marker keratin 18 (Krt18) and the basal marker keratin 14 (Krt14). Basal myoepithelial cells were negative for Krt18 and positive for Krt14 in > 95% of total cells. Conversely, luminal Sca1- and luminal Sca1+ cells were positive for Krt18 and negative for Krt14 in > 95% of total cells. These data confirm the basal and luminal origin of the isolated cell subpopulations. Basal stem/progenitor cells were positive for Krt14 and Krt18 in > 95% and about 20% of total cells, respectively. Data are representative of three independent experiments. Scale bar, 50 μm. (B/C) qPCR for CD49f and Sca1 in FACS-sorted mammary epithelial cell subpopulations. Fold changes are shown relative to myoepithelial cells. Data are expressed as the mean ± SEM of three independent experiments. (PDF 285 KB)

13058_2012_3224_MOESM5_ESM.PDF

Additional file 5: Control for complete involution. Representative images of whole mounts of mammary glands from virgins and parous mice 28 and 40 days after weaning. Mammary glands were completely involuted at 28 days and certainly at 40 days after weaning. (PDF 3 MB)

13058_2012_3224_MOESM6_ESM.PDF

Additional file 6: Influence of cell number on transcriptome analysis and validation of the amplification method. (A) Pairwise correlation plot of transcriptome data derived from 2,000 and 50,000 myoepithelial and luminal Sca1- cells, and from 2,000 basal CD49fHigh stem/progenitor cells isolated from 11-week-old virgin mice (n = 6). Individual arrays were pairwise correlated by using the unfiltered data as input. Pearson correlation coefficients were calculated and mapped onto a gray scale from black (low values) to white (high values). Higher cell numbers resulted in higher reproducibility, as assessed by high Pearson correlation coefficients. However, although lower cell numbers resulted in lower Pearson correlation coefficients, all arrays of one cell subpopulation were clearly discernible from other subpopulations, irrespective of the cell number used. Thus, in the range of 2,000 to 50,000 cells, cell-subpopulation identity was more determining for cluster analyses than cell number. (B) qPCR on amplified cDNA of mammary epithelial cell subpopulations. Data were normalized to the reference gene B2M and are shown relative to 50,000 cells of myoepithelial cells. The basal marker keratin 14 (Krt14) was expressed by myoepithelial cells and basal CD49fHigh stem/progenitor cells, but not by luminal Sca1- and luminal Sca1+ cells. Conversely, the luminal markers keratin 8 (Krt8) and keratin 19 (Krt19) were expressed by luminal Sca1- and luminal Sca1+ cells, but not by myoepithelial and not by basal CD49fHigh stem/progenitor cells. As expected, the estrogen receptor alpha (Esr1) was expressed by luminal Sca1+ cells only. Data represent the means of duplicates. (C) qPCR on unamplified cDNA. Data were processed and analyzed as in (B). Changes in expression levels of the luminal and basal markers Krt19 and Krt14, respectively, were similar to those of amplified cDNA, indicating that the amplification process was unbiased. (PDF 247 KB)

13058_2012_3224_MOESM7_ESM.PDF

Additional file 7: The decrease in Wnt signaling is specific for basal stem/progenitor cells, whereas the p53-p21 pathway is upregulated to the same degree in all mammary epithelial cell subpopulations from parous mice. (A) Enrichment plots of Wnt target gene-set enrichment analysis [28, 29] for isolated mammary epithelial cell subpopulations. The enrichment score is plotted against the ranked gene list, calculated by subtracting the gene expression levels of cells from age-matched virgins and parous mice. The gene set contained all canonical Wnt target genes reported in mammalian systems (Additional file 3). Thus, positive enrichment scores indicate an upregulation, and negative enrichment scores, a downregulation of Wnt signaling in cells of parous mice. Statistical analysis indicated specific downregulation of Wnt signaling in basal stem/progenitor cells from parous mice (Table 1). The apparent downregulation of Wnt signaling in myoepithelial cells, which are contaminated with basal stem/progenitor cells, was not significant. (B) Bar plot of gene-set enrichment analysis (GSEA)-calculated and normalized enrichment scores for the p53-p21 pathway previously identified by Sivaraman et al. [39]. The p53-p21 gene set was rendered by GSEA as significantly upregulated in all isolated mammary epithelial cell subpopulations from parous mice when testing for all signaling-pathway gene sets contained in v2.5 and v3.0. The nominal P value, the false discovery rate (FDR), and the family-wise error rate (FWER) for this pathway were < 0.01 for all mammary epithelial cell subpopulations; 1,000 permutations were performed with the permutation type "gene set." (C) Representative images of immunostaining for p21 and bar graph comparing the relative frequency of p21-positive epithelial cells in mammary gland sections from age-matched virgin and parous mice in estrus. Data represent the mean ± SD (virgin mice: n = 3; parous mice: n = 3). P = 0.0007, by using two-tailed unpaired Student t test. Scale bar, 25 μm. (PDF 4 MB)

13058_2012_3224_MOESM8_ESM.PDF

Additional file 8: Twenty most significantly downregulated pathways in basal stem/progenitor cells after parity. The list was calculated by using v3.0 of GSEA [28]. The 1,000 permutations were performed by using the permutation type "gene set." In all other cases, the default settings were used. (PDF 18 KB)

13058_2012_3224_MOESM9_ESM.PDF

Additional file 9: Blood progesterone concentrations in parous and age-matched virgin control mice in estrus do not change significantly. Plasma progesterone levels were measured with ELISA. (PDF 82 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Meier-Abt, F., Milani, E., Roloff, T. et al. Parity induces differentiation and reduces Wnt/Notch signaling ratio and proliferation potential of basal stem/progenitor cells isolated from mouse mammary epithelium. Breast Cancer Res 15, R36 (2013). https://doi.org/10.1186/bcr3419

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/bcr3419

Keywords