- Research article
- Open Access
Parity induces differentiation and reduces Wnt/Notch signaling ratio and proliferation potential of basal stem/progenitor cells isolated from mouse mammary epithelium
Breast Cancer Researchvolume 15, Article number: R36 (2013)
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.
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.
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.
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.
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 . 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  and for genetically engineered mouse models of breast cancer .
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 . 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 , 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.
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 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 . 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 . 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  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%.
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 ) 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) .
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).
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. 18.104.22.1684 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  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).
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).
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.
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 .
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).
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).
The limited dilution transplantation data were analyzed statistically as published previously . The two-tailed unpaired Student t test was used to determine statistical significance of comparisons.
Motif Activity Response Analysis (MARA)
The MARA model  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.
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 . 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) , and the Wnt targets versican (Vcan) and Igfbp3  (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  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) . 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  occurred in all epithelial cell subpopulations tested (Additional file 7).
Further bioinformatics analysis of the data with different software  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 , 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 . 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 . 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).
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.
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") . 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 , thus inducing canonical Wnt signaling in mammary stem/progenitor cells , 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) . 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).
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.
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 , 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 . 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) , 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 . The latter is especially true for the expression of the classic Wnt target Lgr5 , 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 . 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 , 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  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 . 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  and/or primary downregulation  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 . 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.  found no effects of late pregnancy (9 weeks) on mammary-repopulating units (MRUs) in total mammary epithelial cells, Siwko et al.  (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 .
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.  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 . 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 .
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  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.
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 , 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.
false discovery rate
family-wise error rate
gene-set enrichment analysis
motif activity response analysis
mammary stem cell
mammary repopulating unit
normalized enrichment score
red blood cell
white blood cell.
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.
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.
Medina D: Mammary developmental fate and breast cancer risk. Endocr Relat Cancer. 2005, 12: 483-495. 10.1677/erc.1.00804.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Many AM, Brown AM: Mammary stem cells and cancer: roles of Wnt signaling in plain view. Breast Cancer Res. 2010, 12: 313-
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.
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.
GSEA homepage. [http://www.broadinstitute.org/gsea]
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.
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.
Gonzalez RC, Woods RE, Eddins SL: Digital Image Processing using MATLAB. 2009, Knoxville, TN: Gatesmark Publishing
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.
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.
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.
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.
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.
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.
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.
Ingenuity Systems. [http://www.ingenuity.com]
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.
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.
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.
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.
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.
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.
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
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.
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.
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.
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.
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.
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.
Visvader JE: Keeping abreast of the mammary epithelial hierarchy and breast tumorigenesis. Genes Dev. 2009, 23: 2563-2577. 10.1101/gad.1849509.
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.
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.
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.
Medina D, Kittrell FS: p53 function is required for hormone-mediated protection of mouse mammary tumorigenesis. Cancer Res. 2003, 63: 6140-6143.
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.
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.
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.
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
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.
Wnt homepage. [http://www.wnt.stanford.edu]
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.
The authors declare that they have no competing interests.
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
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)
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)
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)
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. . 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)
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 . 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)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.