- Research article
- Open Access
Gene expression in murine mammary epithelial stem cell-like cells shows similarities to human breast cancer gene expression
Breast Cancer Research volume 11, Article number: R26 (2009)
Mammary stem cells are bipotential and suggested to be the origin of breast cancer development, but are elusive and vaguely characterized. Breast tumors can be divided into subgroups, each one requiring specific treatment. To determine a possible association between mammary stem cells and breast cancer, a detailed characterization of the transcriptome in mammary stem cells is essential.
We have used a murine mammary epithelial stem-like cell line (HC11) and made a thorough investigation of global gene-expression changes during stepwise differentiation using dual-color comparative microarray technique. Subsequently, we have performed a cross-species comparison to reveal conserved gene expression between stem cells and subtype-specific and prognosis gene signatures, and correlated gene expression to in vivo mammary gland development.
Our analysis of mammary stem-like and stepwise cell differentiation, and an in-depth description of our findings in a breast cancer perspective provide a unique map of the transcriptomic changes and a number of novel mammary stem cell markers. We correlate the alterations to in vivo mammary gland differentiation, and describe novel changes in nuclear receptor gene expression. Interestingly, our comparisons show that specific subtypes of breast cancers with poor prognosis and metastasizing capabilities show resemblance to stem-like gene expression.
The transcriptional characterization of these mammary stem-like cells and their differentiation-induced gene expression patterns is here made widely accessible and provides a basis for research on mammary stem-like cells. Our comparisons suggest that some tumors are more stem-like than others, with a corresponding worse prognosis. This information would, if established, be important for treatment decisions. We also suggest several marker candidates valuable to investigate further.
The mammary gland exhibits unique developmental features during puberty, pregnancy and lactation. For each round of pregnancy, the mammary gland undergoes sequential cycles of proliferation, differentiation and apoptosis, and alveolar ducts form and grow, differentiate to produce milk and, after lactation, cease, revert and regress to the pre-pregnancy state. These and several other observations point towards the presence of stem cells as the basis for the capacity for alveolar renewal in each pregnancy (reviewed in ). These stem cells would be the origin of the epithelial compartment, where committed precursor cells become restricted to either a myoepithelial or luminal (ductal or alveolar) fate .
The cellular identity and markers of a rapidly cycling population of normal adult mammary stem cells have been suggested . Some current models explain cancer development by the stem cell and clonal evolution hypothesis . A bipotential mammary epithelial stem cell and/or a luminal progenitor cell would, via acquired genetic alterations, epigenetic changes and paracrine signals from surrounding cells, abandon its controlled stem cell self-renewal and develop into a breast tumor cell with uncontrolled growth [3, 4]. Depending on the cell of origin (stem cell or progenitor cell) and acquired changes, different forms of breast cancer can develop. Tumors are further hypothesized to be sustained by the subpopulation of cancer stem cells , which should be specifically targeted for long-term successful therapeutic intervention.
The only factor known to consistently decrease lifetime breast cancer risk is early childbirth and breastfeeding . This is speculated to be due to decreased circulating hormones, increased differentiation of the mammary cells and/or a decrease of mammary stem cells in the gland . Exposure to estrogen is a risk factor for the development of breast cancer. The effects of estrogens are mediated via two estrogen receptor (ER) isoforms: ERα and ERβ. Activation of ERα drives the proliferation in most of ERα-positive breast tumors, whereas ERβ can behave as an antagonist to ERα at the transcriptional level  and may be of protective value in breast cancer. Epidemiological studies also suggest that the breast is at particular risk to acquire deleterious genetic changes before or during puberty, which is thought to be a period of stem cell expansion. At this time, exposure to phytoestrogens, ligands of ERβ, has also been shown to be protective .
The overwhelming majority of breast cancers have luminal epithelial characteristics. Five molecular subgroups of breast cancer have been defined that can be distinguished according to their gene expression profiles: basal like, luminal A, luminal C, Her2+ and normal breast like . Two nuclear receptors, ERα and its downstream target the progesterone receptor, are two of the most important prognostic markers of breast cancer, whose expression is an indication that anti-estrogenic therapy can be successful. Most luminal A and luminal C tumors express these markers, but luminal A tumors have a considerably better prognosis than luminal C tumors. Most basal-like, normal breast-like and Her2+ subtypes do not express ERα; of these subtypes, Her2+ and basal like have a worse prognosis . There is a need for improved prognostic markers and tailor-made treatment for each of these groups.
The normal mammary epithelial cell line HC11 originates from mid-pregnancy BALB/c mice , and resembles mammary stem cells or progenitor cells remarkably. HC11 cells exhibit both self-renewal and ability for pluripotency: these cells can be cultured for an unlimited number of passages in a proliferating stem cell-like phase, they can be made to functionally differentiate and express milk proteins in vitro  and they can reconstitute the ductal epithelium of a cleared mammary fat pad with myoepithelial, alveolar and ductal luminal cells in vivo . This cell line can be cultured without requirements for complex exogenously added, without extracellular matrix or without co-cultivation with other cell types for the prolactin-dependent in vitro induction of differentiation .
Few studies have been performed to investigate the underlying mechanisms when proliferating mammary epithelial cells differentiate using such pure cell systems of normal stem-like cells. Desrivières and colleagues have performed initial studies at the proteome level [13, 14] in the HC11 cells, and different mammary epithelial subpopulations (basal/myoepithelial, luminal ERα-positive and luminal ERα-negative) were recently characterized at the transcriptome level ; however, a broad description of mammary stem cell or stem cell-like transcriptomes and their change during differentiation is lacking. In vitro models have the benefit of representing homogeneous cell populations, giving enough material for thorough, repeated and high-quality analysis, allowing manipulations to better understand development, characteristics and responses and minimizing animal use.
In the present study, the undifferentiated stage of the HC11 cells is referred to as stem like, since the cells have the capability of both self-renewal and pluripotency, but it is not entirely clear whether the HC11 cells are complete stem cells or are in the process of becoming progenitor cells. In order to investigate the transcriptional events behind this stem cell-like state, to provide a map of mammary epithelial cell differentiation and to determine a possible connection between mammary stem cells and breast cancer, we have performed a complete analysis of all 36,000 known genes and transcripts using dual-color microarray of the stem cell-like stage and different differentiation stages of these cells. We have compared this stem cell-like transcriptome with in vivo mammary gland development and with gene expression signatures of different types of human breast cancers and their metastases, and have analyzed the biological consequences in the context of stem cell and breast cancer interaction in detail. We report potential stem-like markers, and show their specificity at the protein level.
Materials and methods
Cell culture and induction of differentiation
HC11 mammary epithelial cells were routinely grown in complete medium (RPMI 1640, 10% FBS, L-glutamine, 5 μg/ml insulin, 10 ng/ml epidermal growth factor, and 50 μg/ml gentamicin; all from Sigma, St Louis, MO, USA), and proliferating cells were obtained under these growth conditions. When cells reached confluence, the complete medium was changed to medium without epidermal growth factor (RPMI 1640, 2% FBS, 5 μg/ml insulin, and 50 μg/ml gentamicin), and pre-differentiated (competent) cells were obtained after 48 hours of growth in this medium. To induce differentiation of competent cells, cultures were treated for 72 hours with medium without epidermal growth factor containing 100 nM dexamethasone and 1 μg/ml ovine prolactin.
Animals and mammary gland tissue
Female Balb/c mice were fed ad libitum and were kept under a 12-hour light, 12-hour dark cycle. Mammary glands from 2-month old virgin, 10-day pregnant and 6-day lactation mice were excised and frozen in liquid nitrogen for RNA extraction. All animal experimentation was conducted in accordance with accepted standards of humane animal care as outlined by the Stockholm South Ethical Committee of the Swedish National Board for Laboratory Animals.
RNA was extracted with TRIzol (Invitrogen, Carlsbad, CA, USA) according to standard protocol, and purified using Qiagen RNeasy spin columns (Qiagen, Hilden, Germany) with on-column DNase I digestion.
Microarray experiment and analysis
Raw data and detailed protocols for the microarray analysis are available from the ArrayExpress data repository [E-MEXP-969]. In short, Operon's long-oligonucleotide spotted arrays covering the whole known genome represented by 36,000 genes and gene transcripts were used (printed by the microarray facility at KTH, Royal Institute of Technology, Stockholm, Sweden). Triplicate cultures of stem-like, pre-differentiated and differentiated stages of HC11 were analyzed in 10 hybridizations using 20 μg total RNA per sample (20 samples); each microarray comparison was replicated and dye-swapped with different sets of cultures to compensate for unequal dye incorporation and variance between cultures and differentiations. Hybridizations were performed in a loop design, where stem-like stage was hybridized to pre-differentiated stage, pre-differentiated stage to fully differentiated stage and, in addition, stem-like stage to fully differentiated stage. Scanning was carried out at 10 μm resolution using the G2565BA DNA microarray scanner (Agilent Technologies, Santa Clara, CA, USA), for which the photomultiplier tube was set to 100.
The obtained TIFF images were analyzed using the GenePix Pro 6 software (Molecular Devices Corp., Union City, CA, USA), in which TIFF files generated for Cy3 channels and Cy5 channels were superimposed upon each other. Spot identification, manual examination of the surface of the array and flagging of spots/regions with poor quality were performed in GenePix. Result files (gpr files) produced by GenePix were imported into the R environment for statistical computing and programming, where data processing and identification of differentially expressed genes were carried out using the Bioconductor package bundle, Limma, aroma package and the kth-package. Unreliable spots with abnormal physical properties were removed using several filters. After filtering, the slides were normalized with print-tip loess (local regression) normalization. To identify differentially expressed genes, a parametric empirical Bayes approach implemented in Limma was used . This test statistic will assign a score (B-score) to each gene. The B-score was used to rank the genes so that the gene with the highest score has the highest probability of being differentially expressed. When differences were being investigated, two criteria had to be fulfilled for a gene to be regarded as differentially expressed: the gene had to have a B-score of more than 0 and an |M-value| of more than 0.5 (an M-value is the second logarithm of the fold change).
Genes close to cut-off (B > 0.0, |2logFC| > 0.4 and P < 0.005) could be confirmed by real-time PCR in three independent differentiation cultures. Confirmations of 22 genes were performed.
Analysis of over-represented themes and classification into gene ontology functional groups was performed using Expression Analysis Systematic Explorer . The complete mouse transcriptome was used to calculate expected frequencies of over-represented themes. Gene ontology groups were considered over-represented when the calculated Expression Analysis Systematic Explorer score (modified Fischer Exact probability test) was ≤ 0.2.
cDNA was synthesized with Superscript III (Invitrogen) and random hexamers using 1 μg DNase-I-treated and purified RNA. Real-time PCR was performed with SYBR-Green Mastermix (Applied Biosystems) according to the manufacturer's protocol, using 10 ng cDNA per 10 μl reaction in an ABI PRISM 7500 (Applied Biosystems, Foster City, CA, USA) under the following conditions: 95°C for 10 minutes, followed by 40 cycles at 95°C for 15 seconds and 60°C for 40 seconds. All runs were performed in triplicate from triplicate cultures, and specific amplification was checked with melting curve analysis. Primer sequences of intron-spanning fragments will be provided on request. The ΔCt formula was used to determine fold-change differences, and 18S was used as the reference gene.
HC11 cells were grown in 10 cm plates and were differentiated as described above. Cells were washed with PBS, collected in a 1.5 ml tube, and pelleted by centrifugation at 4°C for 2 minutes. The cell pellets underwent one cycle of freeze–thaw and were resuspended in lysis buffer (1% NP-40, 50 mM Tris–HCl, pH 7.5, 140 mM NaCl, 2 mM ethylenediamine tetraacetic acid, 1 mM phenylmethanesulphonylfluoride, 1 mM Na3VO4 and protease inhibitor cocktail; Roche, Mannheim, Germany). The lysate was kept on ice for 20 minutes and was centrifuged at 14,000 rpm for 10 minutes at 4°C. The protein concentration was determined with Bradford reagent (BioRad, Hercules, CA, USA). Whole-cell extracts (40 μg protein) were resolved on SDS-PAGE and transferred onto a Polyvinylidene Fluoride membrane.
The membranes were blocked with 5% (w/v) milk protein dissolved in PBS and were incubated overnight at 4°C with the following primary antibodies: rabbit anti-ADAMTS1 (GTX11557; Genetex, Irvine, CA, USA), chicken anti-DUSP6 (GTX19115; Genetex), rabbit anti-Birc5 (also known as survivin, #2808; Cell Signaling, Houston, TX, USA), rabbit anti-Melk (#2274; Cell Signaling), and mouse anti-COUP TF2 (ab41859; AbCam, Cambridge, UK) and rabbit anti-tubulin (#2125; Cell Signaling) was used as loading control. The antibodies were considered specific based on previous data and detection of few bands in the blots, and in all cases the molecular weight of the corresponding bands was calculated using Quantity One software (BioRad). The secondary antibodies were coupled to horseradish peroxidase (Sigma). The luminescent signal was detected with the enhanced chemiluminescence kit (Amersham, Buckinghamshire, UK).
Immunocytochemistry and immunohistochemistry
HC11 cells were cultured on eight-well glass chamber slides (BD Biosciences, Meylan Cedex, France). Cells were fixed in 10% buffered formalin and permeabilized in 0.5% Triton–PBS for 30 minutes. Unspecific binding was blocked by incubation in blocking solution for 1 hour (10% FBS in 0.1% Tween–PBS). Incubation with the primary antibodies proceeded overnight at room temperature. The primary antibodies used were the same as for western blots and also anti-mouse CD44 (#553991; BD Pharmigen, San Diego, CA, USA), mouse anti-keratin 5/8 (MA1-35858; Affinity Bioreagents, Rockford, IL, USA) and mouse anti-Rab4 (ab13252; AbCam). After three washes with PBS, the corresponding secondary antibodies – anti-mouse-Cy3 (Sigma), anti-rabbit or anti-chicken Alexa 488 (Molecular Probes, Invitrogen, Eugene, OR, USA) and anti-mouse Alexa 568 (Molecular Probes) – were added (1/500 in blocking solution) for 1 hour. Cells were washed four times with PBS and nuclei were stained with ToPro-3 (1/1,000; Molecular Probes). To assess whether the staining was due to binding of the primary antibody, a group of samples was stained in parallel but only with the secondary antibody. Images were captured using a laser-scanning confocal microscope (Zeiss 510, Stockholm, Sweden). Images were acquired under the same settings; when edited with Adobe Photoshop 6.0, the same adjustments were applied to all images.
Results and discussion
Mammary stem-like cell transcriptome
To explore the largely unknown transcriptome of mammary stem cells we compared the murine mammary stem-like cells with their pre-differentiated and fully differentiated stage. In this way, all changes of genes and important networks for commitment of cell fate, implementation of differentiation programs, specification and survival were revealed. The most highly changed genes, at stem cell-like exit, are presented in Tables 1 and 2, gene enrichment analysis of over-represented and under-represented processes is listed in Table 3, and a schematic overview of changed processes is depicted in Figure 1. All regulated genes are available in Additional file 1, and real-time confirmations of 22 genes are available as Additional files 2 to 4. Confirmations at the protein level are depicted in Figures 2 to 5.
Our analysis confirmed previous findings, including: an increase of Stat1 and Stat5a during differentiation, which also mimics the in vivo control of proliferation, differentiation and survival during mammary gland differentiation ; an increase of Ctgf expressed, also found in the in vivo mouse mammary gland during pregnancy and lactation ; an increase of Lamp1, a differentiation marker for HC11 ; and an increase of Igfbp5, whose expression is stimulated during cellular differentiation by lactogenic hormones . Further, a decrease of Igfbp2 , Hnrpd , Cyclin D1 and Myc  during the differentiation of these cells was confirmed. In addition to the above verifications, which in themselves demonstrate a high reproducibility of the differentiation process and of the analysis, a total of 2,251 genes was observed to be changed in the first step of commitment to differentiation, and a further 1,010 alterations during the subsequent stage of functional differentiation. The majority of these alterations are novel, and the implications of many of them are at this time unknown.
When cells left the stem-like stage and entered the pre-differentiation stage, gene enrichment analysis of processes – as defined by the Gene Ontology Consortium – showed that differentiation (particularly epithelial differentiation), skeletal development, cell adhesion, regulation of apoptosis and several types of metabolisms (coenzyme, lipid, carbohydrate, oxygen and sulfur metabolism) increased (Figure 1, box A). At the same time, there was a robust decrease in expression of mitotic cell cycle-associated genes, metabolism of DNA, RNA, nucleotide and proteins, respectively, and angiogenesis (Figure 1, box B). Most genes that were altered in the first transition remained at this level in the second transition, such as upregulated genes within vesicle-mediated transport, cell adhesion and epithelial cell differentiation (Figure 1, box G). Expression of other genes continued to increase: for example, skeletal development and positive regulation of transcription (Figure 1, box E). There was also a significant increase of protein modifications and additional cell adhesion gene expression specific for the second transition (Figure 1, box F). The strong slow-down of proliferation in the first differentiation transition partially continued in the second differentiation step (Figure 1, box J – also signified by a strong decrease in proliferation marker Mki67), whereas most of the decreased cell cycle genes did not change further in the second transition (Figure 1, box H). Organelle organization, biogenesis, organogenesis and regulation of proteolysis and peptidolysis started to decrease in the second transition (Figure 1, box I). A temporary increase of 58 genes, pronounced within programmed cell death genes (Figure 1, box C), and a temporary decrease of protein metabolism genes (Figure 1, box D) were specific changes that only occurred in the pre-differentiated stage.
Among the most strongly altered genes (Tables 1 and 2), transcript levels of amino acid phosphorylation proteins (including Melk) and dephosphorylation proteins (including Dusp6) were highly elevated at the stem-like cell stage and decreased considerably during differentiation. Melk and Dusp6 expression was further corroborated at the protein level (Figure 2a, b), where the decrease upon differentiation was evident. A large group of different histone transcripts were also found to be strongly downregulated during differentiation – for example, Hist1h1a confirmed by real-time PCR (Additional file 2), Hist1h4h, Hist1h3c and H2afz (Table 1). When differentiation starts, gene expression of proteinase inhibitor activity (Expi, confirmed by real-time PCR; Additional file 2), proteolysis proteins (such as cathepsin D and cathepsin A) and transcription-related proteins (Ehf, confirmed by real-time PCR (Additional file 2), Plagl1, Stat1, Stat3, Stat5a and Stat6) increased. From these data we can conclude that the HC11 stem-like cells have a high activity of cell cycle, protein phosphorylation and angiogenic activities coupled with low adhesion, apoptosis, transcriptional activity and differentiation.
Mammary stem cell characteristics in the HC11 stem-like cells
The exact features of mammary stem cells are not fully known, and there are disagreeing reports of mammary stem cell characteristics. Also the exact stem cell characteristics of the cell line used here are unspecified, but exhibit several trademarks of stem cells. The HC11 cells are immortal and have the capability of self-renewal and pluripotency, and maintain properties that allow them to differentiate in vitro in response to lactogenic hormones. Their repopulation potential in vivo may not be perfectly reproducible, however, and the expression of stem cell markers is not well established. Further, these cells are proliferating, which our current dogma suggests stem cells of the mammary gland are not. These cells may be in the process of becoming progenitor cells, or they may possibly constitute a mixture of stem and progenitor cells. Regardless, it is of great interest to define the potential stem-like gene expression pattern and/or markers for these cells.
To evaluate our material for stem cell characteristics, we compared the gene expression of the stem-like stage with findings reported or suggested by others. A bipotential human stem cell is hypothesized to be Cd44+ , an adhesion molecule with roles in signaling, migration and homing. We found that Cd44 was highly expressed at the stem-like cell stage and decreased extensively during differentiation at the mRNA level. Cd44 at the protein level also decreased with differentiation (Figure 3a), although occasional cells remained positive also in the differentiated stage. Furthermore, a receptor with similar properties to Cd44, Hmmr, was strongly overexpressed in the stem-like stage. A strong correlation with other suggested mammary stem cell marker genes was observed; that is, Brca1 , Krt6, Krt5  and Melk, a suggested stem cell gene in hematopoietic , neural  and possibly epidermal stem cells – the latter two genes were observed and corroborated also at protein level (Figure 3b and Figure 2a, respectively).
A novel mechanism for the control of stem cell proliferation in embryonic and neural stem cells involving histone H2afx was recently described . We found this gene to be overexpressed in our stem-like cells. This histone gene has also been shown to exhibit copy number changes in sporadic breast cancer . Vim (vimentin), a suggested stem cell marker in mesenchymal cells , also showed a decrease during differentiation (corroborated at protein level; Additional file 6). Other possible stem cell markers, however, showed an opposite expression in these cells: Lrp5 (a reported cell surface marker for somatic mammary stem cells), Musashi homolog 2 (Msi2, a neuronal stem cell marker) and both Kit and Kitl (stem cell markers in hematopoietic stem cells) were all upregulated at the transcript level during differentiation of HC11 mammary stem cell-like cells.
Murine mammary stem cells have further been selected using protein markers CD45/Ptprc-, Ter119/Ly76-, CD31/Pecam1-, Sca-1/Ly6alow, CD24/Cd24amed and CD49f/Itga6high by Stingl and coworkers . In our HC11 material, however, these markers were not significantly changed at the transcript level. One reason for these apparent differences may be that the protein levels at the cell surface do not always follow the mRNA levels, so even if the transcript for a specific gene does not change, other mechanisms can affect both protein localization and stability. Another reason could be related to where in the process between pure stem cells towards progenitor cells our stem-like cells and/or the literature reported cells are residing. The transcription factor Etv4 (Pea3), suggested to function in multipotential mammary progenitors to regulate their lineage-specific differentiation potential by Kurpios and colleagues , was at the highest expression in the stem-like stage.
HC11 cells have the capacity to differentiate in vivo into both myoepithelial and luminal (ductal and alveolar) epithelial cells. Both markers of myoepithelial lineage (Mme/Cd10) and luminal epithelial lineage (Krt18) increased strongly during differentiation. Further, we compared the expression patterns of the HC11 cell differentiation stages with the three different mammary epithelial cell subpopulations – basal/myoepithelial, luminal ERα-positive and luminal ERα-negative – from virgin mouse mammary gland, characterized at the transcriptome level in a study by Kendrick and colleagues . Genes specific for each of the three subgroups were increasing in expression during differentiation of the HC11 cells, further proving that all lineages are represented in the differentiated stage. We found that whereas nearly all luminal specific genes that were changed during the differentiation of HC11 cells were upregulated (79 out of 82 for luminal ER-positive and 79 out of 84 for luminal ER-negative), about one-half of the genes specific for basal/myoepithelial lineage (99 genes) were upregulated during differentiation, whereas the remaining 91 genes decreased. The reason for this could be that, because the basal cell layer also contains the mammary epithelial stem cell compartment , and genes in the stem cells should be downregulated during differentiation, genes in differentiated cells in the basal/myoepithelial lineage should be correspondingly upregulated during differentiation. The luminal cells should be mostly represented by differentiated cells, and genes in these cells would be expected to be upregulated during HC11 differentiation.
During midgestational mammary development in vivo, several signal transducers and activators of transcription are known to be increased by prolactin; in the HC11 cells, Stat1 (Table 2), Stat3, Stat5a and Stat6 as well as the prolactin-induced protein increased, indicating that the cells retain this mammary epithelial cell-differentiation attribute. Furthermore, the milk protein β-casein was expressed in differentiated cells (Additional file 7).
Activity of the Wnt, Hedgehog and Notch pathways are other hallmarks of stem cell characteristics. In the HC11 cells, drastic changes of expression of Wnt members during differentiation was not observed; however, several known positively regulated Wnt target genes  were strongly downregulated during differentiation (Ccnd1, Birc5, Cd44 and Myc, all corroborated at mRNA and/or protein level (Figures 3 and 4, Additional file 2), and Cyr61, Fosl1, Cd87/Plaur, Met, Fst, Emp1, Abcb1b, Ptgs2, Abcb1b, Runx2, Gja1, IL6, Mycbp), and targets known to be downregulated by Wnt were correspondingly upregulated (Cdh1, confirmed at mRNA level (Additional file 2), Sox9 and Postn), in agreement with a higher Wnt activity at the stem-like cell stage. This normal downregulation of the Wnt pathway during differentiation correlates to a previous report by Shackleton and colleagues showing that normal differentiation could be inhibited by the over activation of the Wnt pathway . Several members of the Hedgehog pathway (including Snail and Prkca) and the Notch pathway (including Jag1, Jag2, Hr, Lfng, Hes1) were also changed during differentiation.
A large study of heterogeneous collections of gene expression data generated from 83 mouse stem cell-related samples defined four super-families of stem cell markers associated with differentiation: serine proteinase inhibitors (serpins), cytochrome P450 family, Rab family GTPases, and nuclear receptors . In HC11 mammary cells, differentially expressed genes signified all of these four groups. Two serpins (Serpine1 and Serpine2) were highly reduced as the stem-like cells underwent differentiation. Serpine1 is also involved in regulation of angiogenesis in breast cancer, and Serpine2 in cell differentiation. In a previous study of differentiating hematopoietic stem cells we also identified two serpins as strongly reduced (Serpin a3g and a3n) , of which expression of Serpina3g has been shown to prevent stem cells from differentiation . We further found an increase of Group 2 genes (two cytochrome P450 members, Cyp2f2 and Cyp4x1, in the HC11 cells) and several Rab family GTPases in Group 3 (increase of Rab4a (confirmed at protein level in Figure 5), Rab1, Rab3a, Rab5b, Rab15, Rab18, Rab25, and related gene family members Rhoj, Rhoq, Rhou; and decrease of two Rab family GTPases Rab12 and Rab32) as the mammary stem-like cells differentiated.
Group 4 (nuclear receptors) was also differentially expressed in HC11 mammary stem-like cells. Ten nuclear receptors changed their expression when the mammary stem-like cells differentiated. COUP-TFII, COUP-TFI, FXRβ, NGFIB, NURR1 and ERβ decreased (Additional file 5). COUP-TFII and COUP-TFI influence proliferation of breast cancer cells [36, 37] and are implicated in metastasis . In the fly the common ancestral gene (svp/NR2F3) regulates stem cell identity of neuroblasts . Figure 2c clearly shows the downregulation of COUP-TFII also at the protein level when the cells start differentiating. Further, RORα, VDR, EAR2 and ERα increased (Additional file 4). RORα is frequently inactivated in breast cancers, and VDR is indicated to be protective against breast cancer . Both COUP-TFII and RORα are among 426 selected markers of stem cells described . Further, we noted that both ERα and ERβ were among the nuclear receptors that changed during differentiation. Both are involved, in opposing manners, in mammary development and breast cancer .
Taken together, a number of stem cell-related changes as well as changes indicative of a mixture of both myoepithelial and luminal cell fates are in line with the stem cell characteristics of these cells.
Mammary stem-like cells show resemblance with breast cancer signatures
To investigate whether there are similarities between human breast cancer gene expression and the stem cell-like expression of the murine HC11 cells, supporting the hypothesis that breast cancer primarily arises from mammary stem/progenitor cells , we compared our material of mammary stem-like gene expression with gene-profiling signatures of breast tumors. These profiles have a prognostic value equal to or better than clinicopathologic variables , and a 70-gene signature is able to distinguish sporadic breast cancer tumors with poor prognosis from those with favorable prognosis . Comparison showed that the poor prognosis signature overlapped with the expression profile of the mammary stem cell-like genes. In contrast, none of the genes indicating favorable prognosis were differentially expressed in the stem-like cells (Table 4).
Among signatures distinguishing subclasses of breast carcinomas , our stem-like signature showed an irrelevant overlap with subclasses nominated as normal breast like, as luminal epithelial containing ER, and as Her2+. For subtype luminal C and basal-like tumors, however, there was a considerable overlap of genes (75% and 60%, respectively; see Table 5). Patients with these two subclasses of breast tumors also show among the lowest survival . The basal-like phenotype has been suggested to resemble normal mammary stem cells [1, 44], which is here demonstrated at the gene expression level. Both the HC11 cells and the poor-prognosis and basal-like tumors are characterized by high proliferation; nonetheless, many of the shared gene signatures are not directly linked to proliferation but to adhesion (for example, Tnc, Ly6e, Cdh3), protein phosphorylation (Melk), transcription (Id1), development (Ext1) and signal transduction (Ect2, Gpr126) (Tables 4 to 6).
Our results are in line with the recent report of poorly differentiated aggressive human tumors showing an embryonic stem cell-like gene expression signature . Here a core set of nine embryonic transcription regulators was found to be overexpressed at the mRNA level in many poorly differentiated tumors, and in the HC11 cells we observed four of these transcription factors to be overexpressed at the stem-like cell stage (Mybl2, Hmga1, Hmgb3, Tead4). In conclusion, the comparisons presented here show that breast cancer subtypes defined by, for example, Sorlie and colleagues  can be further subdivided according to stem cell-like resemblance, and comparison with both the 70-gene signature, and Sorlie's classification reveals that stem cell-like expression infers worse prognosis. We speculate that subtypes with a higher degree of stem cell-like gene expression may have a higher fraction of cancer stem cells, yielding a more aggressive cancer. Specific markers to determine whether a tumor is stem cell-like, of which we here suggest several candidates, could be important for diagnosis and treatment decisions.
Further, in the transcriptome analysis we observed a relation between mammary stem-like cell differentiation and regulation of skeletal development genes (including osteoblastic stem cell markers Spark and Spp1 ). In addition, both Il6, which functions as a differentiation regulator of preosteoblast cells , and the corresponding downstream osteoblast-specific differentiation marker Runx2 decreased when stem-like cells entered differentiation, and Ocil, a negative regulator of osteoclast differentiation, showed a robust increase at both differentiation stages. This finding may indicate why breast tumors have a preference for skeletal metastases, and these genes may have a potential as metastasis markers. Indeed, other genes implicated in bone metastasis of breast cancer cells (Ctfg, Fst and Dusp1  and Adamts1 [47, 48]) were altered during the differentiation of mammary stem-like cells. Expression of Adamts1 at the protein level is also shown in Figure 2d as elevated at the stem-like cell stage. In addition, lung metastasis gene expression also has an apparent parallel to stem cell-like gene expression (Table 6); lung metastasis signature genes ) change considerably during differentiation, most of them being elevated at the stem cell-like stage – for example, the cytokine angiopoietin-like 4, shown to prime breast cancer cells for lung metastasis . The correlation of stem cell-like gene expression to metastasis signatures may in part explain the above correlation to poor prognosis.
Our approach using cross-species comparisons of murine mammary stem cell-like expression and human tumor gene expression to unlock evolutionarily conserved breast cancer–stem cell networks has provided highly concordant observations. Furthermore, this approach recently gained support, as cross-species comparisons were shown to be a powerful means of identifying essential connections .
Mammary stem-like cell differentiation compared with in vivomammary gland
We were interested to see whether the changes observed during stem-like cell differentiation showed any resemblance to the in vivo mammary gland differentiation, keeping in mind that in the mammary gland the stem cells only constitute a minority of all cells, and that their gene expression is likely to be masked by changes in other cells as well as by changes in the relative proportions of different cell types. The cellular three-dimensional structure, interaction with stroma and in vitro versus in vivo signaling, also makes the two systems very different. We compared mRNA levels of HC11 cells with 2-month-old virgin, pregnant and lactating mammary glands, where the proliferating stage could be compared with the actively proliferating pregnant mammary gland. Genes that were overexpressed at the stem cell-like stage (Birc5, Areg, Ereg, Cyclin D, Lif and Hist1h1a) all had increased expression in the pregnant glands but decreased their expression in the lactating glands (data available in Additional file 2). For the genes whose expression was low at the stem cell-like stage but was upregulated as the cells differentiated, several were also expressed at a low level in the virgin gland and were upregulated in pregnant and/or lactating gland (Expi, Ecad, Perp, Ehf, Nfat), whereas two genes (Msi2h and Mmp15) showed an opposite regulation and decreased during in vivo gland differentiation (Additional file 2). Although this comparison is relatively simplistic, our data indicate that stem cell-like proliferation is highest in the pregnant mammary gland whereas genes robustly expressed at terminal differentiation of HC11 cells are also highly expressed in the lactating mammary gland.
Among the nuclear receptors, we found that COUP-TFII and NURR1 decreased during differentiation of the HC11 cells and in the transition from pregnant to lactating mammary gland. ERβ and VDR, both suggested to be protective in breast cancer [8, 40, 51], were further induced in differentiated (lactating) mammary gland compared with virgin gland, whereas ERα – which is often increased in breast cancers – was reduced. Comparison of a published study of in vivo mouse mammary glands, investigating 10-week nulliparous mammary glands and 18-day pregnant glands , with our data on HC11 cellular differentiation shows that as many as 279 genes changed in HC11 stem-like cells were also changed in the pregnant mammary gland, further showing that there is numerous correlations between the in vitro and in vivo systems in terms of differentiation.
The aim of the present study was to characterize the differentiation process of the stem cell-like HC11 cell line and to define the transcriptome of proliferating undifferentiated mammary epithelial stem cell-like cells, in relation to their differentiated counterparts. This provides a basis for research on mammary stem cells; both known and novel stem cell gene expression characteristics were found. Characterized mammary stem cell markers are highly needed, and several potentially suitable targets detected in our study will be valuable to investigate further.
We explored whether there is a link between mammary stem-like cell gene expression and that of breast cancer. The present study was performed using an in vitro system of mammary stem-like cell differentiation. In vitro systems have drawbacks, and rarely fully resemble the in vivo situation, but can nonetheless yield significant information. We found an interesting correlation between the pattern of stem-like cell expression and that of human breast cancer with poor prognosis, metastasis and tumor subtypes. This may indicate that some breast tumors have a high ratio of cancer stem cells, and may require specific and aggressive treatment. An amount of novel gene expression data is presented, with implications for stem cell and mammary gland development biology and breast cancer research. A scheme is provided where key differences between differentiation steps can be dissected. We conclude that HC11 cells are relevant for studying mammary stem cells, their differentiation and their relationship to breast cancer.
fetal bovine serum
polymerase chain reaction.
Vaillant F, Asselin-Labat ML, Shackleton M, Lindeman GJ, Visvader JE: The emerging picture of the mouse mammary stem cell. Stem Cell Rev. 2007, 3: 114-123. 10.1007/s12015-007-0018-2.
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.
Polyak K: Breast cancer: origins and evolution. J Clin Invest. 2007, 117: 3155-3163. 10.1172/JCI33295.
Cobaleda C, Cruz J, Gonzalez-Sarmiento R, Sanchez-Garcia I, Perez-Losada J: The emerging picture of human breast cancer as a stem cell-based disease. Stem Cell Rev. 2008, 4: 67-79. 10.1007/s12015-008-9012-6.
Visvader JE, Lindeman GJ: Cancer stem cells in solid tumours: accumulating evidence and unresolved questions. Nat Rev Cancer. 2008, 8: 755-768. 10.1038/nrc2499.
Kelsey JL, Gammon MD, John EM: Reproductive factors and breast cancer. Epidemiol Rev. 1993, 15: 36-47.
Britt K, Ashworth A, Smalley M: Pregnancy and the risk of breast cancer. Endocr Relat Cancer. 2007, 14: 907-933. 10.1677/ERC-07-0137.
Williams C, Edvardsson K, Lewandowski SA, Strom A, Gustafsson J-A: A genome-wide study of the repressive effects of estrogen receptor beta on estrogen receptor alpha signaling in breast cancer cells. Oncogene. 2008, 27: 1019-1032. 10.1038/sj.onc.1210712.
Warri A, Saarinen NM, Makela S, Hilakivi-Clarke L: The role of early life genistein exposures in modifying breast cancer risk. Br J Cancer. 2008, 98: 1485-1493. 10.1038/sj.bjc.6604321.
Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, Rijn van de M, Jeffrey SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Lonning PE, Borresen-Dale A-L: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci. 2001, 98: 10869-10874. 10.1073/pnas.191367098.
Ball RK, Friis RR, Schoenenberger CA, Doppler W, Groner B: Prolactin regulation of beta-casein gene expression and of a cytosolic 120-kd protein in a cloned mouse mammary epithelial cell line. EMBO J. 1988, 7: 2089-2095.
Humphreys R, Rosen J: Stably transfected HC11 cells provide an in vitro and in vivo model system for studying Wnt gene function. Cell Growth Differ. 1997, 8: 839-849.
Desrivières S, Prinz T, Castro-Palomino Laria N, Meyer M, Boehm G, Bauer U, Schafer J, Neumann T, Shemanko C, Groner B: Comparative proteomic analysis of proliferating and functionally differentiated mammary epithelial cells. Mol Cell Proteomics. 2003, 2: 1039-1054. 10.1074/mcp.M300032-MCP200.
Desrivières S, Kuhn K, Müller J, Gläser M, Castro-Palomino LN, Korder J, Sonnentag M, Neumann T, Schwarz J, Schäfer J, Hamon C, Groner B, Prinz T: Comparison of the nuclear proteomes of mammary epithelial cells at different stages of functional differentiation. Proteomics. 2007, 7: 2019-2037. 10.1002/pmic.200600994.
Kendrick H, Regan J, Magnay F-A, Grigoriadis A, Mitsopoulos C, Zvelebil M, Smalley M: Transcriptome analysis of mammary epithelial subpopulations identifies novel determinants of lineage commitment and cell fate. BMC Genomics. 2008, 9: 591-10.1186/1471-2164-9-591.
Smyth G: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: article 3-
Hosack DA, Dennis JG, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4: R70-10.1186/gb-2003-4-10-r70.
Petersen H, Haldosén L-A: EGF modulates expression of STAT5 in mammary epithelial cells. Exp Cell Res. 1998, 243: 347-358. 10.1006/excr.1998.4160.
Weihan W, Bethanie M, Traci G, Cynthia CJ, Nicholas K, Cutler ML: Glucocorticoid induced expression of connective tissue growth factor contributes to lactogenic differentiation of mouse mammary epithelial cells. J Cell Physiol. 2008, 214: 38-46. 10.1002/jcp.21159.
Cella N, Cornejo-Uribe RR, Montes GS, Hynes NE, Chammas R: The lysosomal-associated membrane protein LAMP-1 is a novel differentiation marker for HC11 mouse mammary epithelial cells. Differentiation. 1996, 61: 113-120. 10.1046/j.1432-0436.1996.6120113.x.
Phillips K, Park M, Quarrie L, Boutinaud M, Lochrie J, Flint D, Allan G, Beattie J: Hormonal control of IGF-binding protein (IGFBP)-5 and IGFBP-2 secretion during differentiation of the HC11 mouse mammary epithelial cell line. J Mol Endocrinol. 2003, 31: 197-208. 10.1677/jme.0.0310197.
Nagaoka K, Tanaka T, Imakawa K, Sakai S: Involvement of RNA binding proteins AUF1 in mammary gland differentiation. Exp Cell Res. 2007, 313: 2937-2945. 10.1016/j.yexcr.2007.04.017.
Grolli S, Accornero P, Ramoni R, Donofrio G, Whitelaw CBA: Expression of c-myc is down-regulated as mouse mammary epithelial cells become confluent. Biochem Biophys Res Commun. 1997, 239: 566-569. 10.1006/bbrc.1997.7498.
Liu S, Ginestier C, Charafe-Jauffret E, Foco H, Kleer CG, Merajver SD, Dontu G, Wicha MS: BRCA1 regulates human mammary stem/progenitor cell fate. Proc Natl Acad Sci. 2008, 105: 1680-1685. 10.1073/pnas.0711613105.
Boecker W, Moll R, Dervan P, Buerger H, Poremba C, Diallo RI, Herbst H, Schmidt A, Lerch MM, Buchwalow IB: Usual ductal hyperplasia of the breast is a committed stem (progenitor) cell lesion. J Pathol. 2002, 198: 458-467. 10.1002/path.1241.
Saito R, Tabata Y, Muto A, Arai K-i, Watanabe S: Melk-like kinase plays a role in hematopoiesis in the zebra fish. Mol Cell Biol. 2005, 25: 6682-6693. 10.1128/MCB.25.15.6682-6693.2005.
Nakano I, Paucar AA, Bajpai R, Dougherty JD, Zewail A, Kelly TK, Kim KJ, Ou J, Groszer M, Imura T, Freije WA, Nelson SF, Sofroniew MV, Wu H, Liu X, Terskikh AV, Geschwind DH, Kornblum HI: Maternal embryonic leucine zipper kinase (MELK) regulates multipotent neural progenitor proliferation. J Cell Biol. 2005, 170: 413-427. 10.1083/jcb.200412115.
Andäng M, Hjerling-Leffler J, Moliner A, Lundgren TK, Castelo-Branco G, Nanou E, Pozas E, Bryja V, Halliez S, Nishimaru H, Wilbertz J, Arenas E, Koltzenburg M, Charnay P, El Manira A, Ibañez CF, Ernfors P: Histone H2AX-dependent GABA(A) receptor regulation of stem cell proliferation. Nature. 2008, 451: 460-464. 10.1038/nature06488.
Srivastava N, Gochhait S, Gupta P, Bamezai RN: Copy number alterations of the H2AFX gene in sporadic breast cancer patients. Cancer Genet Cytogenet. 2008, 180: 121-128. 10.1016/j.cancergencyto.2007.09.024.
Krzyzanowski P, Andrade-Navarro M: Identification of novel stem cell markers using gap analysis of gene expression data. Genome Biol. 2007, 8: R193-10.1186/gb-2007-8-9-r193.
Kurpios NA, MacNeil L, Shepherd TG, Gludish DW, Giacomelli AO, Hassell JA: The Pea3 Ets transcription factor regulates differentiation of multipotent progenitor cells during mammary gland development. Dev Biol. 2009, 325: 106-121. 10.1016/j.ydbio.2008.09.033.
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.
The Wnt Homepage. [http://www.stanford.edu/~rnusse/wntwindow.html]
Richter K, Wirta V, Dahl L, Bruce S, Lundeberg J, Carlsson L, Williams C: Global gene expression analyses of hematopoietic stem cell-like cell lines with inducible Lhx2 expression. BMC Genomics. 2006, 7: 75-10.1186/1471-2164-7-75.
Hampson IN, Hampson L, Pinkoski M, Cross M, Heyworth CM, Bleackley RC, Atkinson E, Dexter TM: Identification of a serpin specifically expressed in multipotent and bipotent hematopoietic progenitor cells and in activated T cells. Blood. 1997, 89: 108-118.
Nakshatri H, Mendonca MS, Bhat-Nakshatri P, Patel NM, Goulet RJ, Cornetta K: The orphan receptor COUP-TFII regulates G2/M progression of breast cancer cells by modulating the expression/activity of p21WAF1/CIP1, cyclin D1, and cdk2. Biochem Biophys Res Commun. 2000, 270: 1144-1153. 10.1006/bbrc.2000.2562.
Le Dily F, Métivier R, Guéguen MM, Le Péron C, Flouriot G, Tas P, Pakdel F: COUP-TFI modulates estrogen signaling and influences proliferation, survival and migration of breast cancer cells. Breast Cancer Res Treat. 2008, 110: 69-83. 10.1007/s10549-007-9693-6.
Minn AJ, Gupta GP, Siegel PM, Bos PD, Shu W, Giri DD, Viale A, Olshen AB, Gerald WL, Massague J: Genes that mediate breast cancer metastasis to lung. Nature. 2005, 436: 518-524. 10.1038/nature03799.
Kanai MI, Okabe M, Hiromi Y: Seven-up controls switching of transcription factors that specify temporal identities of Drosophila neuroblasts. Dev Cell. 2005, 8: 203-213. 10.1016/j.devcel.2004.12.014.
Cui Y, Rohan TE: Vitamin D, calcium, and breast cancer risk: a review. Cancer Epidemiol Biomarkers Prev. 2006, 15: 1427-1437. 10.1158/1055-9965.EPI-06-0075.
Wicha MS, Liu S, Dontu G: Cancer stem cells: an old idea – a paradigm shift. Cancer Res. 2006, 66: 1883-1890. 10.1158/0008-5472.CAN-05-3153.
Gong Y, Symmans WF, Pusztai L: Gene-expression microarrays provide new prognostic and predictive tests for breast cancer. Pharmacogenomics. 2007, 8: 1359-1368. 10.2217/146224126.96.36.1999.
van't Veer LJ, Dai H, Vijver van de MJ, He YD, Hart AAM, Mao M, Peterse HL, Kooy van der K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415: 530-536. 10.1038/415530a.
Yehiely F, Moyano JV, Evans JR, Nielsen TO, Cryns VL: Deconstructing the molecular portrait of basal-like breast cancer. Trends Mol Med. 2006, 12: 537-544. 10.1016/j.molmed.2006.09.004.
Ben-Porath I, Thomson MW, Carey VJ, Ge R, Bell GW, Regev A, Weinberg RA: An embryonic stem cell-like gene expression signature in poorly differentiated aggressive human tumors. Nat Genet. 2008, 40: 499-507. 10.1038/ng.127.
Li Y, Bäckesjö C-M, Haldosén L-A, Lindgren U: IL-6 receptor expression and IL-6 effects change during osteoblast differentiation. Cytokine. 2008, 43: 165-173. 10.1016/j.cyto.2008.05.007.
Kang Y, Siegel PM, Shu W, Drobnjak M, Kakonen SM, Cordón-Cardo C, Guise TA, Massagué J: A multigenic program mediating breast cancer metastasis to bone. Cancer Cell. 2003, 3: 537-549. 10.1016/S1535-6108(03)00132-6.
Liu YJ, Xu Y, Yu Q: Full-length ADAMTS-1 and the ADAMTS-1 fragments display pro- and antimetastatic activity, respectively. Oncogene. 2006, 25: 2452-2467. 10.1038/sj.onc.1209287.
Padua D, Zhang XH, Wang Q, Nadal C, Gerald WL, Gomis RR, Massagué J: TGFβ primes breast tumors for lung metastasis seeding through angiopoietin-like 4. Cell. 2008, 133: 66-77. 10.1016/j.cell.2008.01.046.
Bennett C, Green J: Unlocking the power of cross-species genomic analyses: identification of evolutionarily conserved breast cancer networks and validation of preclinical models. Breast Cancer Res. 2008, 10: 213-10.1186/bcr2125.
Palmieri C, Cheng G, Saji S, Zelada-Hedman M, Warri A, Weihua Z, Van Noorden S, Wahlstrom T, Coombes R, Warner M, Gustafsson J: Estrogen receptor beta in breast cancer. Endocr Relat Cancer. 2002, 9: 1-13. 10.1677/erc.0.0090001.
Wang M, Master S, Chodosh L: Computational expression deconvolution in a complex mammalian organ. BMC Bioinformatics. 2006, 7: 328-10.1186/1471-2105-7-328.
The authors thank MSc Maria Berling for valuable assistance. The study was supported by grants from the Swedish Cancer Society, the Lars Hierta Memorial Foundation, the Magnus Bergvall Foundation, and the David and Britt Haydens Foundation.
The authors declare that they have no competing interests.
CW conceived of the study and participated in the design and coordination, carried out the transcriptome studies, real-time PCR analysis, bioinformatic analysis and comparisons, and drafted the manuscript. LH participated in the design, carried out the cell cultures, mammary gland tissue dissection and immunostainings, and helped to draft the manuscript. KE performed laboratory work and helped to draft the manuscript. L-AH participated in the design of the study and helped to draft the manuscript. J-ÅG helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional data file 4: Adobe file containing a figure that shows the real-time PCR confirmations of differentially expressed genes: expression of nuclear receptors and/or related genes increasing during HC11 mammary stem-like differentiation and correlating changes in vivo mammary glands. (PDF 17 KB)
Additional data file 5: Adobe file containing a figure that shows the real-time PCR confirmations of differentially expressed genes: expression of nuclear receptors and/or related genes decreasing during HC11 mammary stem-like differentiation and correlating changes in vivo mammary glands. (PDF 18 KB)
About this article
Cite this article
Williams, C., Helguero, L., Edvardsson, K. et al. Gene expression in murine mammary epithelial stem cell-like cells shows similarities to human breast cancer gene expression. Breast Cancer Res 11, R26 (2009). https://doi.org/10.1186/bcr2256
- Breast Cancer
- Mammary Gland
- Stem Cell Marker
- Sporadic Breast Cancer
- Mammary Stem Cell