Transcriptome analyses of mouse and human mammary cell subpopulations reveal multiple conserved genes and pathways
- Elgene Lim†1, 2,
- Di Wu†1, 2,
- Bhupinder Pal†1,
- Toula Bouras1,
- Marie-Liesse Asselin-Labat1,
- François Vaillant1,
- Hideo Yagita3,
- Geoffrey J Lindeman1, 4, 5,
- Gordon K Smyth1, 2 and
- Jane E Visvader1, 2Email author
© Lim et al.; licensee BioMed Central Ltd. 2010
Received: 2 January 2010
Accepted: 26 March 2010
Published: 26 March 2010
Molecular characterization of the normal epithelial cell types that reside in the mammary gland is an important step toward understanding pathways that regulate self-renewal, lineage commitment, and differentiation along the hierarchy. Here we determined the gene expression signatures of four distinct subpopulations isolated from the mouse mammary gland. The epithelial cell signatures were used to interrogate mouse models of mammary tumorigenesis and to compare with their normal human counterpart subsets to identify conserved genes and networks.
RNA was prepared from freshly sorted mouse mammary cell subpopulations (mammary stem cell (MaSC)-enriched, committed luminal progenitor, mature luminal and stromal cell) and used for gene expression profiling analysis on the Illumina platform. Gene signatures were derived and compared with those previously reported for the analogous normal human mammary cell subpopulations. The mouse and human epithelial subset signatures were then subjected to Ingenuity Pathway Analysis (IPA) to identify conserved pathways.
The four mouse mammary cell subpopulations exhibited distinct gene signatures. Comparison of these signatures with the molecular profiles of different mouse models of mammary tumorigenesis revealed that tumors arising in MMTV-Wnt-1 and p53-/- mice were enriched for MaSC-subset genes, whereas the gene profiles of MMTV-Neu and MMTV-PyMT tumors were most concordant with the luminal progenitor cell signature. Comparison of the mouse mammary epithelial cell signatures with their human counterparts revealed substantial conservation of genes, whereas IPA highlighted a number of conserved pathways in the three epithelial subsets.
The conservation of genes and pathways across species further validates the use of the mouse as a model to study mammary gland development and highlights pathways that are likely to govern cell-fate decisions and differentiation. It is noteworthy that many of the conserved genes in the MaSC population have been considered as epithelial-mesenchymal transition (EMT) signature genes. Therefore, the expression of these genes in tumor cells may reflect basal epithelial cell characteristics and not necessarily cells that have undergone an EMT. Comparative analyses of normal mouse epithelial subsets with murine tumor models have implicated distinct cell types in contributing to tumorigenesis in the different models.
The mammary gland comprises a ductal epithelial network embedded in a stromal matrix. The ducts are composed of an inner layer of luminal cells and an outer layer of myoepithelial cells. Pregnancy is accompanied by the expansion and differentiation of alveolar luminal cells, resulting in secretory cells that produce and secrete milk. Although the function of the mammary gland is preserved across species, marked anatomic differences exist between human and mouse mammary tissue. The human mammary gland is characterised by a branching network of ducts that terminate in clusters of small ductules that constitute the terminal ductal lobular units (TDLUs). In contrast, the mouse mammary epithelial tree does not contain TDLUs, although small alveolar buds are formed during each estrous cycle. Moreover, the human breast parenchyma is significantly more fibrous than the mouse stroma, which contains predominantly adipocytes. Despite these architectural differences, accumulating evidence suggests that remarkable parallels are found between the hierarchy of epithelial cells that exist in the mammary glands of humans and mice .
Distinct epithelial subtypes have been prospectively isolated from both mouse [2–5] and human mammary glands [6–10]. Functionally analogous subpopulations have been identified: the MaSC-enriched/bipotent progenitor, committed luminal progenitor and mature luminal cell subsets. In the mouse, MaSCs are found within the basal CD49fhiCD29hiCD24+Sca1- subset (referred to as MaSC-enriched), whereas committed luminal progenitor cells exhibit a CD29loCD24+CD61+ (or Sca-1-CD24+) phenotype, and mature luminal cells display a CD29loCD24+CD61-phenotype [2, 3]. In human mammary tissue, the CD49fhiEpCAM-/lo subpopulation has been demonstrated to be enriched for MaSCs, based on in vivo transplantation either into the mouse mammary fat pad  or under the renal capsule . Luminal progenitor and differentiated cells prospectively isolated from human breast tissue are characterized by CD49fhiEpCAM+ and CD49f-EpCAM+ phenotypes, respectively.
There are similarities as well as species-specific differences in the expression of cell-surface markers on the epithelial subsets. Both the mouse and human MaSC-enriched populations express high levels of CD49f. However, CD24 is a marker of epithelial cells in the mouse mammary gland, but not in human breast tissue, where it specifically marks luminal cells [3–5, 7, 11]. Significantly, both the human and mammary MaSC-enriched populations lack expression of the steroid hormone receptors ERα and PR [7, 12]. Moreover, these MaSCs do not express detectable levels of ERBB2/HER2, reminiscent of the triple-negative receptor phenotype that characterizes many basal cancers .
Understanding the relation between normal epithelial cell types and the different molecular subtypes of breast cancer is fundamental to gaining insight into cell types predisposed to carcinogenesis. At least six distinct subtypes of breast tumors have been defined on the basis of gene expression profiling. These include the luminal A and B, basal-like, claudin-low, HER2/ERBB2-overexpressing, and normal breast-like subtypes . We recently used the emerging human mammary hierarchy as a framework for understanding aberrant cell subsets that may arise during breast oncogenesis . The claudin-low subtype was found to be most closely associated with the gene signature of the MaSC-enriched population, whereas the molecular profiles of the basal-like subtype of breast cancer showed remarkable concordance with the luminal progenitor gene signature. Not surprisingly, the expression profiles of the luminal A and B subtypes were closest to that of mature luminal epithelial cells. Interestingly, the molecular portrait of premalignant tissue from BRCA1 mutation carriers, who usually develop basal-like breast cancers, showed striking similarity to the luminal progenitor signature .
In the context of the mouse mammary gland, transcriptome analyses of epithelial cells have highlighted the differences between basal and luminal cells and revealed a number of potential regulators [5, 15]. Here we performed genome-wide transcriptome analyses of three different mouse epithelial subpopulations and established pathways that are conserved in functionally equivalent subsets in humans by using specific gene signatures. We further used these signatures to interrogate mouse models of mammary tumors, providing insight into cell types that contribute to breast oncogenesis.
Materials and methods
Mice and mammary cell preparations
A single cell suspension of mammary cells was prepared from freshly harvested mammary glands and sorted by flow cytometry, as previously described . Mice were on a pure FVB/N background. All experiments were approved by the WEHI Animal Ethics Committee, and the care of animals was in accordance with institutional guidelines. Experiments using human tissue obtained from the Royal Melbourne Hospital Tissue Bank were approved by the Human Research Ethics Committees of The Walter and Eliza Hall Institute of Medical Research and Melbourne Health.
Antibodies, staining and cell sorting
Unless otherwise specified, antibodies for flow cytometry were obtained from BD Pharmingen. Antibodies against mouse antigens were PE-conjugated antibody to CD24, FITC-conjugated antibody to CD29 (clone HMbeta1-1 from H. Yagita) , biotin-conjugated antibodies to CD31, CD45, and TER119, and APC-conjugated antibody to CD61 (Caltag). Antibodies used for human antigens have previously been described . The Alexa Fluor 647 anti-human CD24 antibody (Biolegend) was used for analysis of human breast epithelial subsets. Antibody staining and cell sorting was as previously described . Data were analyzed by using WEASEL software .
RNA preparation and quantitative RT-PCR analysis
Total RNA was isolated from primary mammary cell subpopulations with the RNeasy Micro kit (Qiagen). Reverse transcription by using oligo(dT) primer and Moloney murine leukemia virus reverse transcriptase (Invitrogen) was according to the manufacturer's protocol. Quantitative RT-PCR was carried out by using a Rotorgene RG-6000 (Corbett Research) and SensiMix (dT) DNA kit (Quantace) under the following conditions: 10 min at 95°C followed by 35 cycles consisting of 15 seconds at 95°C, 20 seconds at 62°C, and 20 seconds at 72°C. Gene expression was determined with the Rotor-Gene software (version 1.7). The primer sequences used are listed in Supplementary Methods in Additional file 1.
Total RNA was purified from sorted cell populations by using the RNeasy Micro kit (Qiagen). RNA quality was assessed with the Agilent Bioanalyzer 2100 (Agilent Technologies) by using the Agilent RNA 6000 Nanokit (Agilent Technologies) according to the manufacturer's protocol. Up to 500 ng of RNA was labeled with the standard Total Prep RNA amplification kit (Ambion), and complementary RNA (1.5 μg) was hybridized to Illumina MouseWG-6 v2.0 BeadChips. After washing, the chips were coupled with Cy3 and scanned by using an Illumina BeadArray Reader. Unnormalized summary probe profiles, with associated probe annotation, were output from BeadStudio.
Microarray data analysis: normal cell subpopulations
Raw intensities were normalized by using the neqc function, which performs normexp background correction and quantile normalization by using control probes . Probes were filtered if not detected in any sample (detection p value, 0.01). The mouse data are deposited as GEO series GSE19446, and the human, as GSE16997.
Microarray data analysis: mouse model tumors
Expression profiles of mouse tumors were downloaded from GEO series GSE3165 . Fifty-six Agilent arrays (Agilent-011978 Mouse Microarray G4121A) profiling mouse tumor models of interest were included in the analysis. The samples and arrays are described by Herschkowitz et al. . Data analysis used the raw Agilent Feature Extraction data files and probe annotation from GEO. Control probes were filtered, and then expression values were normexp background corrected with offset 16 , and then log-ratios were global loess normalized . Two MMTV-Wnt-1 samples and one MMTV-Neu sample were removed as outliers on the basis of unsupervised clustering.
Subpopulation expression signatures
Pairwise comparisons were made between the three epithelial cell populations by using empiric Bayes-moderated t statistics  and array quality weights . Allowance was made for possible correlations between RNA samples drawn from the same pool of mice . The false discovery rate (FDR) was controlled by using the Benjamini and Hochberg algorithm. Probes with FDR < 0.05 and fold-change > 1.5 were judged to be differentially expressed. For each subpopulation (MaSC-enriched, luminal progenitor, and mature luminal), signature probes were defined as those that were significantly differentially expressed in the same direction versus both of the other two cell subpopulations. For stromal cells, the signature probes were defined relative to the three epithelial cell populations.
For each target sample (mouse tumor or normal mouse mammary cell subpopulation), a set of signature scores was computed to measure the transcriptional activity of each mouse cell subpopulation in that sample, by using a method previously described . The signature score is essentially the average log-expression of the signature genes in the target sample, weighted by the direction and magnitude of change of those genes in the mouse subpopulation used to define the signature. Higher scores indicate that the transcriptional signature of the mouse cell subpopulation is found in the target sample.
Conserved signature genes
A larger set of mouse signature genes were defined by using the "nestedF" multiple-testing option of limma with FDR < 0.1. The 1.5 fold-change threshold was maintained. Mouse and human probes were matched by gene symbol by using the Jackson Laboratory orthology report of 13 November 2009 . If multiple probes mapped to the same symbol, the probe with the highest average log-expression was used. Human signature genes were defined as for mouse, with the multiple testing step repeated only for those human genes with orthologues among the mouse signature genes.
Ingenuity pathway analysis
Ingenuity Pathway Analysis (IPA)  was conducted on conserved signature genes. For the MaSC-enriched subpopulation, only the top 300 signature genes were used, to make the numbers comparable for the different subpopulations. The signature sets were overlaid with canonic pathways. Canonic pathways were selected based on known biologic significance of the most highly overlapping pathways, and were displayed by using "subcellular layout". Direct associations between signature genes were drawn by using the "connect" option. The luminal progenitor signature set was too small to generate connections, so direct associations were drawn from the KIT and CYP24A1 genes by using the "grow" option. Genes without connections to other signature genes were removed from the final figures.
Derivation of distinct gene signatures for mouse mammary cell subpopulations
Number of signature probes and genes for each mouse cell subpopulation (stringent criteria)
Mammary epithelial cell subsets
Mouse genes with human orthologues
Mouse gene signatures correlate with specific mouse models of breast cancer
Comparison of the expression profiles of human and mouse subpopulations
To relate the two species more closely, we examined the transcriptional activity scores of the mouse subtype-specific signature genes in each of the human RNA samples. This demonstrated a conserved expression profile for each epithelial cell subtype across the two species. For each subpopulation, the mouse transcriptional signature was consistently highest in the corresponding human subtype for every biologic replicate (Figure 3c). The two luminal subtypes showed intermediate cross-over transcriptional activity with each other, whereas the transcriptional activity of the MaSC-enriched subset was relatively specific (Figure 3c).
Conservation of signature genes between mouse and human mammary subpopulations
Number of mouse signature probes and genes conserved in human mammary epithelial cell subpopulations
Mammary epithelial cell subsets
Mouse genes with human orthologues
Mouse genes conserved in humans
The conserved upregulated genes in the MaSC-enriched population spanned diverse gene ontology groups, including transcription factors (for example, Irx4, Mef2c, Slug, Egr2, Twist2, Tbx2, Id4, p63, and Sox11), cytokeratins (Krt5, 14, 16), and plasma membrane proteins (for example, Lgr6 and the receptors for Oxytocin, Oncostatin M, and Lif). The Notch ligand Jag2 was highly expressed in this subpopulation, and its product may directly signal through Notch receptors expressed on adjacent luminal cells . The Wnt/β-catenin pathway is anticipated to be active in self-renewing MaSCs, compatible with the observation that Fzd8 and Tcf4 are components of the conserved upregulated gene signature in the MaSC-enriched population. The Wnt-pathway inhibitors Wif1 and Dkk3, however, were also found to be abundantly expressed. These antagonists may be expressed and secreted by mature myoepithelial cells present within this population to attenuate Wnt signal transduction in the basally located MaSCs.
For the luminal progenitor signature, Kit (receptor tyrosine kinase), Elf5 (Ets transcription factor), Cyp24A1 (vitamin D metabolizing enzyme), Lbp, and Cxcr4 were highly expressed in both species. Aldh1a3 was also upregulated in luminal progenitor cells versus other cell types, although another isoform Aldh5a1 was identified in the luminal-restricted population isolated by Raouf et al. . Interestingly, virtually all the ALDH activity in human breast epithelium resides within the luminal progenitor population (unpublished data) rather than the more primitive mammary stem cell subset . In mature luminal cells, highly expressed genes included the transcription factors Foxa1, Myb, estrogen receptor (ER), progesterone receptor (PgR), and Tbx3, as well as the prolactin receptor (Prlr) and Rank ligand (Tnfsf11).
Conservation of canonic pathways between mouse and human subpopulations
In this study, we describe a comparative transcriptome analysis of functionally analogous human and mouse mammary cell populations using an Illumina platform. Four prospectively isolated populations were evaluated, corresponding to those enriched for basal/mammary stem cells, committed luminal progenitor, mature luminal epithelial, and stromal cells. Distinct gene signatures were apparent for the mouse subpopulations, reminiscent of that found for human mammary cell subsets . Comparison of the mammary epithelial signatures across human and mouse, combined with Ingenuity pathway analysis, revealed a number of conserved genes and pathways that are likely to regulate key processes during mammary ontogeny.
The MaSC-enriched subset exhibited the largest number of genes conserved across species. This subset comprises stem cells, likely basal progenitor cells, as well as mature myoepithelial cells. These cells share many common cell-surface markers that have impeded efforts to fractionate this population. Multiple transcriptional regulators (Irx4, Mef2C, Slug, Egr2, Twist2, Tbx2) were found to be highly expressed in this basal subset. Interestingly, the leucine-rich repeat-containing G protein-coupled receptor Lgr6 , which belongs to the same subgroup as Lgr5, a stem cell marker of small intestine, colon, and hair follicles , was identified as a component of the MaSC-gene signature. A prominent integrin network also emerged; these proteins play an important role in mediating interactions between basal cells (including MaSCs) and the underlying extracellular matrix. Of relevance, several genes attributed to cells that have undergone an EMT , such as slug, vimentin, and absence of E-cadherin expression, also characterize basal cells in the mammary gland. Therefore, the expression of these genes in breast tumor cells may indicate the acquisition of basal cell characteristics rather than an EMT. The recently described link between Wnt signaling and the EMT  may also reflect an active Wnt pathway in MaSCs or other cells in this basal population.
Kit, Cyp24A1, and Elf5 appear to be defining markers of committed luminal progenitor cells in both species. Interestingly, the tyrosine kinase KIT was reported to be overexpressed in basal breast cancers  and BRCA1-associated basal cancers , suggesting that it may serve as a useful prognostic marker or therapeutic target. Elf5 has been demonstrated to be important for driving alveolar cell differentiation during pregnancy  but may play an earlier role in regulating luminal cell-fate decisions. Interestingly, triple-negative breast cancer patients have been shown to have lower serum vitamin D levels, and Cyp24A1 is known to catabolize both 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D. It is therefore tempting to speculate that higher levels of CYP24A1 might be linked to increased breast cancer risk . Other interesting candidates include CXCR4, a receptor implicated in mediating metastasis of breast cancer cells through its ligand SDF-1 , and CD14 and lipopolysaccharide-binding protein (LBP), which are implicated in Toll-like receptor signaling and LPS-mediated inhibition. In the mature luminal population, active pathways identified by IPA include the Wnt ligands (Wnt4, 5A, 7B), which may act on MaSCs to enhance their self-renewal or proliferation. Expression of the transcriptional regulator Lmo4 was downregulated in the mature luminal subset, consistent with findings that this oncoprotein is important for promoting mammary epithelial cell proliferation and inhibiting differentiation [39, 40]. Conversely, the expression of other transcriptional regulators (ERα, Myb, PR, and Cited1) was significantly upregulated in mature luminal cells.
A high degree of concordance was found between the expression profiles of the basal and mature luminal cell subsets in the mouse mammary gland described here and those previously reported [5, 15], although the luminal progenitor expression profiles proved to be more divergent. For the basal/MaSC-enriched population, conserved pathways such as the Ephrin, Wnt, and extracellular matrix networks were also identified as nodes in interaction mapping of the basal subset by Kendrick et al. . Similarly, the gene expression profiles of the mature luminal subset (reported here) shared substantial overlap with that of the ER+ population described by Kendrick et al. , with the ER/glucocorticoid receptor signalling network emerging as one of the predominant nodes. The expression profile of the luminal progenitor subpopulation, however, exhibited substantial differences from that of the ER-  and Ma-CFC subsets , indicating that they may represent distinct or heterogeneous cell populations. Nevertheless, the Kit and TLR signaling pathways identified here using Ingenuity Pathway analysis were also revealed as distinct modules in the ER- network by ROCK analysis . The gene profiles determined for the same three epithelial subsets isolated from human mammary tissue by Raouf et al.  show similarities but also differences that likely reflect short-term culture of their cells before gene expression studies .
Interrogation of breast cancer subtypes with the gene signatures of normal human epithelial cell subsets has revealed striking relations. Intriguingly, the luminal progenitor gene signature shared marked similarity with the basal subtype of breast cancer and preneoplastic breast tissue from BRCA1 mutation carriers . Moreover, aberrant luminal progenitor cells were detected in BRCA1 mutation carriers, suggesting that they serve as a target population for further oncogenic events . To extend these studies and identify candidate cell types that might contribute to oncogenesis in mouse models, we explored the link between the mouse mammary epithelial hierarchy and models of mammary tumorigenesis. The MaSC-enriched transcriptional signature was highest in MMTV-Wnt-1 and p53-/- tumors, indicating that cells within these tumors exhibit similarities with MaSCs or basal progenitors. Although cancer stem cell populations have been identified in these tumors [41–43], one cannot conclude that these bear resemblance to MaSCs based on expression profiling studies. Rather, the molecular profiles may indicate a cell type that has been expanded during tumor progression. It is notable that preneoplastic tissue from MMTV-Wnt-1 transgenic mice has been shown to harbor an expanded mammary stem cell pool as well as aberrant bipotent progenitor cells, indicating that more than one cell of origin may exist in the Wnt-1 model .
The luminal progenitor signature was highest in MMTV-Neu and MMTV-PyMT tumors. Compatible with this observation for the MMTV-Neu model, FACS analyses of these tumors has indicated a homogeneous population of cells expressing high levels of the luminal progenitor marker CD61+ . Thus, luminal progenitor cells may have undergone expansion in these tumors. The MMTV-Neu strain, however, does not accurately recapitulate HER2-overexpressing cancers arising in women, because MMTV-Neu tumors do not show significant gene overlap with the HER2-positive subtype but are more similar to human "luminal" tumors . Interestingly, the mature luminal signature was highest in MMTV-Myc tumors. The small progenitor subset in the 'mature' population [2, 28], rather than the differentiated luminal cells, is likely to contribute to tumorigenesis in this model.
The high degree of conservation between analogous epithelial subtypes across species supports the use of mice as a model system to study normal mammary gland development and oncogenesis. The conserved pathways pinpoint those that are likely to be involved in cell-fate decisions and lineage differentiation in the basal or luminal epithelial cell lineages. In the context of breast cancer, genes within the conserved signatures, such as those that characterize the more purified luminal progenitor subset (for example, KIT, CYP24A1, ELF5), have the potential to provide novel prognostic markers or therapeutic targets in breast cancer.
epithelial mesenchymal transition
Ingenuity Pathway Analysis
mammary stem cell
mouse mammary tumor virus
standard error of the mean
terminal ductal lobular unit.
We are grateful to K. Stoev and M. Everest for excellent animal husbandry and microarray expression profiling, respectively. Microarray experiments were carried out in the Australian Genome Research Facility, Melbourne. This work was supported by the Victorian Cancer Agency through the Victorian Breast Cancer Research Consortium, National Health and Medical Research Council (NHMRC, Australia) and the Australian Cancer Research Foundation. EL was supported by the NHMRC and National Breast Cancer Foundation; MA, by the Australian Research Council; TB, by the National Breast Cancer Foundation; and GJL, GKS, and JEV by the NHMRC.
- Visvader JE: Keeping abreast of the mammary epithelial hierarchy and breast tumorigenesis. Genes Dev. 2009, 23: 2563-2577. 10.1101/gad.1849509.PubMedPubMed CentralView ArticleGoogle Scholar
- Asselin-Labat ML, Sutherland KD, Barker H, Thomas R, Shackleton M, Forrest NC, Hartley L, Robb L, Grosveld FG, Wees van der J, Lindeman GJ, Visvader JE: Gata-3 is an essential regulator of mammary-gland morphogenesis and luminal-cell differentiation. Nat Cell Biol. 2007, 9: 201-209. 10.1038/ncb1530.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- Eirew P, Stingl J, Raouf A, Turashvili G, Aparicio S, Emerman JT, Eaves CJ: A method for quantifying normal human mammary epithelial stem cells with in vivo regenerative ability. Nat Med. 2008, 14: 1384-1389. 10.1038/nm.1791.PubMedView ArticleGoogle Scholar
- Lim E, Vaillant F, Wu D, Forrest NC, Pal B, Hart AH, Asselin-Labat ML, Gyorki DE, Ward T, Partanen A, Feleppa F, Huschtscha LI, Thorne HJ, Fox SB, Yan M, French JD, Brown MA, Smyth GK, Visvader JE, Lindeman GJ: Aberrant luminal progenitors as the candidate target population for basal tumor development in BRCA1 mutation carriers. Nat Med. 2009, 15: 907-913. 10.1038/nm.2000.PubMedView ArticleGoogle Scholar
- Shipitsin M, Campbell LL, Argani P, Weremowicz S, Bloushtain-Qimron N, Yao J, Nikolskaya T, Serebryiskaya T, Beroukhim R, Hu M, Halushka MK, Sukumar S, Parker LM, Anderson KS, Harris LN, Garber JE, Richardson AL, Schnitt SJ, Nikolsky Y, Gelman RS, Polyak K: Molecular definition of breast tumor heterogeneity. Cancer Cell. 2007, 11: 259-273. 10.1016/j.ccr.2007.01.013.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Villadsen R, Fridriksdottir AJ, Ronnov-Jessen L, Gudjonsson T, Rank F, LaBarge MA, Bissell MJ, Petersen OW: Evidence for a stem cell hierarchy in the adult human breast. J Cell Biol. 2007, 177: 87-101. 10.1083/jcb.200611114.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Asselin-Labat ML, Shackleton M, Stingl J, Vaillant F, Forrest NC, Eaves CJ, Visvader JE, Lindeman GJ: Steroid hormone receptor status of mouse mammary stem cells. J Natl Cancer Inst. 2006, 98: 1011-1014. 10.1093/jnci/djj267.PubMedView ArticleGoogle Scholar
- Carey LA, Dees EC, Sawyer L, Gatti L, Moore DT, Collichio F, Ollila DW, Sartor CI, Graham ML, Perou CM: The triple negative paradox: primary tumor chemosensitivity of breast cancer subtypes. Clin Cancer Res. 2007, 13: 2329-2334. 10.1158/1078-0432.CCR-06-1109.PubMedView ArticleGoogle Scholar
- Herschkowitz JI, Simin K, Weigman VJ, Mikaelian I, Usary J, Hu Z, Rasmussen KE, Jones LP, Assefnia S, Chandrasekharan S, Backlund MG, Yin Y, Khramtsov AI, Bastein R, Quackenbush J, Glazer RI, Brown PH, Green JE, Kopelovich L, Furth PA, Palazzo JP, Olopade OI, Bernard PS, Churchill GA, Van Dyke T, Perou CM: Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome Biol. 2007, 8: R76-10.1186/gb-2007-8-5-r76.PubMedPubMed CentralView ArticleGoogle Scholar
- Kendrick H, Regan JL, Magnay FA, Grigoriadis A, Mitsopoulos C, Zvelebil M, Smalley MJ: 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Noto K, Kato K, Okumura K, Yagita H: Identification and functional characterization of mouse CD29 with a mAb. Int Immunol. 1995, 7: 835-842. 10.1093/intimm/7.5.835.PubMedView ArticleGoogle Scholar
- WEASEL for Flow Cytometry Data Analysis. [http://www.wehi.edu.au/faculty/advanced_research_technologies/flow_cytometry/weasel_for_flow_cytometry_data_analysis/]
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.PubMedPubMed CentralView ArticleGoogle Scholar
- Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article 3Google Scholar
- Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, Smyth GK: A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007, 23: 2700-2707. 10.1093/bioinformatics/btm412.PubMedView ArticleGoogle Scholar
- Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo]
- Smyth GK, Speed T: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.PubMedView ArticleGoogle Scholar
- Ritchie ME, Diyagama D, Neilson J, van Laar R, Dobrovic A, Holloway A, Smyth GK: Empirical array quality weights in the analysis of microarray data. BMC Bioinformatics. 2006, 7: 261-10.1186/1471-2105-7-261.PubMedPubMed CentralView ArticleGoogle Scholar
- Smyth GK, Michaud J, Scott HS: Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics. 2005, 21: 2067-2075. 10.1093/bioinformatics/bti270.PubMedView ArticleGoogle Scholar
- Mouse Genomics Informatics - Mammalian Orthology. [http://www.informatics.jax.org/orthology.shtml]
- Ingenuity Systems. [http://www.ingenuity.com]
- Cardiff RD, Anver MR, Gusterson BA, Hennighausen L, Jensen RA, Merino MJ, Rehm S, Russo J, Tavassoli FA, Wakefield LM, Ward JM, Green JE: The mammary pathology of genetically engineered mice: the consensus report and recommendations from the Annapolis meeting. Oncogene. 2000, 19: 968-988. 10.1038/sj.onc.1203277.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Ginestier C, Hur MH, Charafe-Jauffret E, Monville F, Dutcher J, Brown M, Jacquemier J, Viens P, Kleer CG, Liu S, Schott A, Hayes D, Birnbaum D, Wicha MS, Dontu G: ALDH1 is a marker of normal and malignant human mammary stem cells and a predictor of poor clinical outcome. Cell Stem Cell. 2007, 1: 555-567. 10.1016/j.stem.2007.08.014.PubMedPubMed CentralView ArticleGoogle Scholar
- Charafe-Jauffret E, Ginestier C, Iovino F, Wicinski J, Cervera N, Finetti P, Hur MH, Diebel ME, Monville F, Dutcher J, Brown M, Viens P, Xerri L, Bertucci F, Stassi G, Dontu G, Birnbaum D, Wicha MS: Breast cancer cell lines contain functional cancer stem cells with metastatic capacity and a distinct molecular signature. Cancer Res. 2009, 69: 1302-1313. 10.1158/0008-5472.CAN-08-2741.PubMedPubMed CentralView ArticleGoogle Scholar
- Hsu SY, Kudo M, Chen T, Nakabayashi K, Bhalla A, Spek van der PJ, van Duin M, Hsueh AJ: The three subfamilies of leucine-rich repeat-containing G protein-coupled receptors (LGR): identification of LGR6 and LGR7 and the signaling mechanism for LGR7. Mol Endocrinol. 2000, 14: 1257-1271. 10.1210/me.14.8.1257.PubMedView ArticleGoogle Scholar
- Barker N, van Es JH, Kuipers J, Kujala P, Born van den M, Cozijnsen M, Haegebarth A, Korving J, Begthel H, Peters PJ, Clevers H: Identification of stem cells in small intestine and colon by marker gene Lgr5. Nature. 2007, 449: 1003-1007. 10.1038/nature06196.PubMedView ArticleGoogle Scholar
- Thiery JP, Acloque H, Huang RY, Nieto MA: Epithelial-mesenchymal transitions in development and disease. Cell. 2009, 139: 871-890. 10.1016/j.cell.2009.11.007.PubMedView ArticleGoogle Scholar
- DiMeo TA, Anderson K, Phadke P, Fan C, Perou CM, Naber S, Kuperwasser C: A novel lung metastasis signature links Wnt signaling with cancer cell self-renewal and epithelial-mesenchymal transition in basal-like breast cancer. Cancer Res. 2009, 69: 5364-5373. 10.1158/0008-5472.CAN-08-4135.PubMedPubMed CentralView ArticleGoogle Scholar
- Schneider BP, Winer EP, Foulkes WD, Garber J, Perou CM, Richardson A, Sledge GW, Carey LA: Triple-negative breast cancer: risk factors to potential targets. Clin Cancer Res. 2008, 14: 8010-8018. 10.1158/1078-0432.CCR-08-1208.PubMedView ArticleGoogle Scholar
- Oakes SR, Naylor MJ, Asselin-Labat ML, Blazek KD, Gardiner-Garden M, Hilton HN, Kazlauskas M, Pritchard MA, Chodosh LA, Pfeffer PL, Lindeman GJ, Visvader JE, Ormandy CJ: The Ets transcription factor Elf5 specifies mammary alveolar cell fate. Genes Dev. 2008, 22: 581-586. 10.1101/gad.1614608.PubMedPubMed CentralView ArticleGoogle Scholar
- Holick MF: Vitamin D deficiency. N Engl J Med. 2007, 357: 266-281. 10.1056/NEJMra070553.PubMedView ArticleGoogle Scholar
- Orimo A, Gupta PB, Sgroi DC, Arenzana-Seisdedos F, Delaunay T, Naeem R, Carey VJ, Richardson AL, Weinberg RA: Stromal fibroblasts present in invasive human breast carcinomas promote tumor growth and angiogenesis through elevated SDF-1/CXCL12 secretion. Cell. 2005, 121: 335-348. 10.1016/j.cell.2005.02.034.PubMedView ArticleGoogle Scholar
- Sum EY, Segara D, Duscio B, Bath ML, Field AS, Sutherland RL, Lindeman GJ, Visvader JE: Overexpression of LMO4 induces mammary hyperplasia, promotes cell invasion, and is a predictor of poor outcome in breast cancer. Proc Natl Acad Sci USA. 2005, 102: 7659-7664. 10.1073/pnas.0502990102.PubMedPubMed CentralView ArticleGoogle Scholar
- Visvader JE, Venter D, Hahm K, Santamaria M, Sum EY, O'Reilly L, White D, Williams R, Armes J, Lindeman GJ: The LIM domain gene LMO4 inhibits differentiation of mammary epithelial cells in vitro and is overexpressed in breast cancer. Proc Natl Acad Sci USA. 2001, 98: 14452-14457. 10.1073/pnas.251547698.PubMedPubMed CentralView ArticleGoogle Scholar
- Cho RW, Wang X, Diehn M, Shedden K, Chen GY, Sherlock G, Gurney A, Lewicki J, Clarke MF: Isolation and molecular characterization of cancer stem cells in MMTV-Wnt-1 murine breast tumors. Stem Cells. 2008, 26: 364-371. 10.1634/stemcells.2007-0440.PubMedView ArticleGoogle Scholar
- Vaillant F, Asselin-Labat ML, Shackleton M, Forrest NC, Lindeman GJ, Visvader JE: The mammary progenitor marker CD61/beta3 integrin identifies cancer stem cells in mouse models of mammary tumorigenesis. Cancer Res. 2008, 68: 7711-7717. 10.1158/0008-5472.CAN-08-1949.PubMedView ArticleGoogle Scholar
- Zhang M, Behbod F, Atkinson RL, Landis MD, Kittrell F, Edwards D, Medina D, Tsimelzon A, Hilsenbeck S, Green JE, Michalowska AM, Rosen JM: Identification of tumor-initiating cells in a p53-null mouse model of breast cancer. Cancer Res. 2008, 68: 4674-4682. 10.1158/0008-5472.CAN-07-6353.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.