1Department of Pharmacology,
2Department of Genetics and Genome Sciences and Department of General Medical Sciences-Oncology,
© The Author(s) 2019. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
The existence of mammary epithelial stem cell (MaSC) populations capable of mediating mammary gland development and homeostasis has been established for over a decade. A combination of lineage tracing and mammary gland transplantation studies has affirmed that MaSCs and their downstream progenitors are organized in a hierarchal manner; however, these techniques have failed to illuminate the complete spectrum of epithelial intermediate populations or their spatial and temporal relationships. The advent of single cell sequencing technology has allowed for characterization of highly heterogeneous tissues at high resolution. In the last two years, the remarkable advances in single cell RNA sequencing technologies have been leveraged to address the heterogeneity of the mammary epithelium. These studies have afforded fresh insights into the transcriptional differentiation hierarchy and its chronology. Importantly, these data have led to a major conceptual shift in which the rigid boundaries separating stem, progenitor, and differentiated epithelial populations have been deconstructed, resulting in a new more fluid and flexible model of epithelial differentiation. The emerging view of the mammary epithelial hierarchy has important implications for mammary development, carcinogenesis, and metastasis, providing novel insights into the underlying cellular states that may promote malignant phenotypes.
Mammary gland, mammary epithelial cells, breast cancer, stem cell
Many of the processes that occur during mammary gland morphogenesis are appropriated during the formation and malignant progression of breast cancer. These include the regeneration and hormonally-driven proliferative expansion of the mammary epithelium, remodeling of the extracellular matrix and invasion into surrounding stroma, aversion of apoptosis, and recruitment of vasculature[1,2]. Thus, understanding the endocrine, cellular, and molecular mechanisms that control mammary gland development should provide critical insights into the processes that also underlie breast cancer growth and metastases and reveal novel approaches for preventing and treating aggressive disease. This review focuses on recent advances in our understanding of the development of the mammary epithelial cell hierarchy, revealing the abundant progression of cell states, many of which could undergo malignant transformation and give rise to the heterogeneous disease known as “breast cancer”.
The mammary gland is a uniquely dynamic organ that undergoes extensive remodeling throughout life in response to cytokine signaling and hormonal cues[3,4]. The primitive mammary epithelium is first detected around embryonic day 11.5 (E11.5) in the mouse. Between embryonic day 16 and embryonic day 18 (E18), this rudimentary structure invades the mammary fat pad, forming a primitive ductal structure. During puberty, the mammary epithelium expands to form a branched, ductal network that invades the stroma until it reaches the edge of the mammary fat pad. This ductal structure is maintained into adulthood where it remains highly flexible in response to oscillating levels of ovarian hormones during the estrous cycle, undergoing extensive proliferation and apoptosis. In pregnancy, a unique hormonal milieu triggers the formation of secretory alveolar structures to form at the terminal end of each duct. Alveolar units function to produce milk proteins during lactation. Once lactation has ceased, the gland undergoes a process of involution in which massive architectural remodeling occurs, resulting in an epithelial structure that resembles the pre-pregnant state. Although morphologically similar, genetically engineered mouse models permitted tracing of “terminally” differentiated lactating cells and revealed that these cells contribute to the remodeled post-involuted gland and are molecularly distinct from their virgin counterparts[6,7]. The events of expansion, differentiation, and repopulation of the mammary tree reoccur with each subsequent pregnancy indicating extensive plasticity of the cells that comprise the gland. Closer examination reveals that the mammary gland consists of two main epithelial lineages, the inner luminal lineage and the outer myoepithelial lineage. The luminal lineage can be further subdivided into mature ductal and milk-producing alveolar cells. Myoepithelial cells are commonly referred to as “basal” since they are located on the basal side of the luminal cells. These cells are highly contractile, and function to facilitate the movement of milk though the ducts.
The ability of the gland to undergo successive rounds of reconstruction in response to multiple pregnancies implicates the existence of a mammary stem cell (MaSC) population capable of generating all mature cell types within the gland. Studies utilizing transplantation and lineage tracing techniques further established that mammary epithelial differentiation is organized in a hierarchical, branched pathway wherein multipotent MaSCs give rise to a series of downstream intermediates[8-10]. With each differentiation step, cells become increasingly lineage-restricted moving from multi- to bi- to uni-potent progenitor states until terminally differentiating into mature epithelial populations.
Discoveries over the last 50 years have been integral to our understanding of the mammary epithelial cell hierarchy; however, the advent of single-cell (sc) sequencing technologies coupled with microfluidic cell sorting has recently uncovered intricacies in lineage relationships and heterogeneity of the mammary epithelium. Recent evidence has demonstrated that each cell is unique in its transcriptomic signature[11-14]. This revelation has redefined the differentiation cascade as a continuous trajectory of heterogeneous populations passing through gradual transcriptional states. These studies suggest that historical, rigid, models of the epithelial hierarchy should be replaced by a model that incorporates the diverse transcriptional profiles and paths represented throughout development, underscoring the immense diversity and flexibility that likely underlies mammary gland development and gives rise to breast cancer heterogeneity. Here we discuss the recent studies utilizing scRNA-sequencing of mammary epithelium and how these new findings have shifted our current perspectives of mammary gland development and breast cancer.
In 1959, DeOme and colleagues developed the mammary transplantation assay. Using this technique, transplantation of murine ductal fragments[16,17] or dissociated suspensions of normal mammary epithelium into the cleared fat pad of a recipient mouse resulted in the regeneration of a complete ductal structure. Smith and Medina then reported that morphologically distinct cells with regenerating capacity could be found throughout the entire mammary tree and persisted throughout pregnancy, lactation, and involution. Both the frequency and regenerative ability of these cells was not affected by age or reproductive history[19,20]. Later studies demonstrated that mammary gland reconstitution was achievable upon transplantation of a single cell. This was first reported by Kordon and Smith, using serial transplantation and limiting dilution assays of mammary epithelium labeled with unique viral insertions generated by mouse mammary tumor virus (MMTV). Similarly, by labeling cells with a LacZ transgene, Shackleton et al. confirmed this finding and demonstrated that stem cell populations were enriched in the basal epithelial fraction (defined by Lin-CD29hiCD24+ cell surface markers). Additionally, cells from reconstituted glands harbored the same repopulating activity of the original and were capable of fully differentiating into milk-producing alveolar cells during pregnancy. These studies implicated the existence of a long lived multipotent mammary stem cell population (MaSC), capable of giving rise to all downstream lineages within the gland.
Following the disclosure of MaSCs, the search for molecular markers specific to the stem cell population became a major focus within the field. Transplantation of luminal vs. basal epithelial populations revealed that cells exhibiting stem-like properties were located within the basal fraction and could be enriched by FACS sorting using CD49fhiCD29hiCD24+/modEpCAM+Sca1low expression as markers[22-25]. Additionally, expression of s-SHIP, LGR5[27-29], Lrp5/6, and Axin2 enriched for populations with repopulating activity upon transplantation; however, stem cell activity was not exclusive to these subsets and the degree to which these markers may enrich for overlapping stem cell populations was not clear. Difficulty surrounding the identification of MaSC markers was due, in part, to the rarity of MaSCs in the adult gland. The frequency of repopulating cells was found to be higher during embryogenesis (~7% of total cells by E18 compared to ~2%-5% frequency in the adult gland)[22,23,32,33]; however, enriched populations of fetal (f)MaSCs displayed altered gene expression profiles from their adult counterparts further complicating the isolation and characterization of a multipotent population.
Transplantation studies demonstrated the existence of a stem cell entity that was capable of multi-lineage differentiation within the adult gland. Yet, various studies employing lineage tracing analyses led to confusion surrounding these findings. While some groups provided evidence to support multipotent cells, others reported that only unipotent cells remain after birth. Early lineage tracing experiments by Van Keymeulen et al., used the promoters of the keratin 14, 5, and 8 (K14, K5, K8) genes, to drive tracer expression and track both luminal (K8) and myoepithelial (K14, K5) lineages. This analysis demonstrated that prior to birth, K14+ cells were multipotent; however, after birth, these cells became unipotent, only giving rise to myoepithelial cells. Cells labeled by the K8 promoter revealed a second unipotent population responsible for maintaining the luminal lineage after birth. These unipotent populations supported epithelial expansion during puberty and pregnancy underscoring their long lived capacity within the gland. Additionally, lineage tracing at the single cell level revealed that unipotent populations were dispersed throughout the ductal tree, estimating that within each major duct, ~20 luminal and ~15 basal unipotent cells contribute to ductal expansion during puberty. Furthermore, by employing clonal analyses at saturation using DOX-inducible systems to trace luminal (K8rtTA-CreTetO/Rosa-TdTomato) and myoepithelial (K14rtTA-CreTetO/Rosa-TdTomato) lineages, Wuidart et al. also reported unipotent progenitors as solely responsible for maintenance of adult lineages. Evidence of luminal unipotent progenitors that are restricted to the alveolar fate has also been reported. In this case, the “parity-induced mammary epithelial cell” (PI-MEC) population consists of alveolar-like cells that originate during a female’s first pregnancy. Unlike the majority of alveolar cells, PI-MECs do not undergo apoptosis during involution. Rather, these cells remain as a residual population within the parous gland and serve as alveolar precursors during subsequent pregnancies[6,37].
Contrary to these results, some genetic fate-mapping studies have supported the existence of multi- and bi-potent mammary epithelial cell populations within the adult. For example, tracing of the Wnt target protein C receptor (Procr) illustrated the presence of a multipotent cell within the basal layer of the adult gland that was capable of giving rise to all mature lineages. Furthermore, studies have also reported bipotent progenitors in the adult gland. By combining a multi-color Cre reporter with three dimensional imaging techniques, Rios and colleagues traced luminal (ELF5+) and basal (Keratin 5+) cells during puberty and adulthood. This analysis resulted in the tracing of discreet clonal patches of luminal or basal cells as well as patches containing both lineages, indicating a common cellular origin. Similarly, an inducible Cre reporter mouse line, driven by the Axin2 promoter, labeled both bipotent and unipotent cells within the adult gland. Further complicating these studies, Axin2+ cells can undergo cell fate switching. This phenomenon has also been reported in mammary epithelial cells that are positive for Lgr5. Interestingly, Axin2+ unipotent cells are basally-restricted during adulthood but become bipotent during pregnancy, suggesting that cell fate switching may occur in response to hormonal cues[10,38]. Pregnancy-induced bipotent mammary epithelial cells were also recently reported in a study by Song et al.. Lineage tracing of K8+ luminal cells revealed bipotent luminal cells that gave rise to basal progeny during pregnancy or upon stimulation with estrogen or progesterone. Additionally, transplantation of luminal-derived basal cells into cleared fat pads of recipient mice, demonstrated their capacity to generate new mammary gland structures with normal morphology.
The conflicting results generated from transplantation studies with those from lineage tracing has been a source of controversy in the field. The primary question still remains: does a multipotent stem cell exist after birth as indicated by transplantation assays; or are progenitor populations solely responsible for postnatal development? Inherent flaws within transplantation and lineage tracing methods likely contribute to the conflicting evidence within the field.
Transplantation assays have been highly scrutinized since this method includes the dissociation of mammary epithelial cells from their normal tissue architecture followed by their reintroduction into a new mammary tissue environment that has undergone invasive surgical resection. The drastic micro-environmental changes that epithelial cells must endure for this procedure has raised concerns that this assay may impart multipotent potential onto cells that would otherwise be restricted to specific differentiation outcomes. For example, Wnt/-catenin-responsive cells only give rise to myoepithelial cells during puberty and pregnancy, however; upon transplantation these cells can regenerate both luminal and myoepithelial lineages. Additionally, in lineage tracing experiments, K14+ cells are restricted to the myoepithelial lineage after birth; however, transplantation of K14-labeled cells from 4-week old mice results in the formation of new ductal structures containing both myoepithelial and luminal cells. The ability of transplantation to alter stem cell fate has also been demonstrated in other systems, including hair follicle development and hematopoiesis. It is plausible that differences in the microenvironment such as stromal, hormonal, and inflammatory cues, may impart phenotypic changes in mammary epithelial cell populations, activating primitive precursor pathways in their lineage. This may be especially relevant in the case of mammary transplants as the surgical resection needed to clear the innate epithelia likely elicits a “wound healing” response in the local environment of the transplant. Thus, transplantation assays may not accurately represent stem cell biology as it would intrinsically occur in vivo.
Although lineage tracing allows for the visualization of cellular differentiation within the context of the normal tissue environment, it is not without its own set of limitations. The ability to detect rare heterogeneous populations of cells using genetic lineage tracing is highly dependent on recombination efficiency. High variability in the extent of recombination not only exists between mouse models (i.e., K5-Cre versus K14-Cre) but also between individual mice from the same strain. In addition, labeling efficiency varies significantly throughout the ductal tree of a single mammary gland. Several explanations for this variability have been postulated[8,10,42]. Inadequate strength and specificity of the promoter driving recombination could result in inconsistent labeling. This may be particularly relevant to quiescent stem cells that are transcriptionally uncommitted to a certain lineage. Additionally, inducible creERT2 models are reliant on accurate dosing of Tamoxifen, a selective estrogen receptor (ER) modulator, to induce sufficient recombination. Since mammary gland development is highly dependent on ER signaling, tamoxifen administration, even at low does, can alter normal mammary gland biology and lead to false readouts.
Taken together, the results from both transplantation and lineage tracing studies have led to a highly compartmentalized view of the mammary epithelial hierarchy in which fetal MaSCs give rise to a common bipotent progenitor which differentiates into lineage-restricted and unipotent luminal and myoepithelial progenitors. Unipotent progenitors then give rise to mature myoepithelial and luminal lineages (for detailed reviews see Visvader and Stingl, 2014, and Yang et al., 2017) [Figure 1]. Although transplantation and lineage tracing methods have been fundamental to our understanding of mammary epithelial biology, it has become increasingly evident that these methods are insufficient to fully characterize rare or highly heterogeneous cell populations. Gene expression signatures on bulk isolates of the major epithelial compartments including luminal progenitors, mature luminal, and MaSC-enriched basal epithelial cells, have been generated. However, the utility of these signatures to identify critical molecular regulators and markers of epithelial subsets has been limited due to the heterogeneity of each epithelial population in which only a small fraction of cells expresses the same set of transcripts (as in the case of stem populations within the basal fraction)[22,32,46]. As a single cell can give rise to an entire mammary gland, it is imperative that the molecular states and instructive signals occurring at the single cell level be examined to fully grasp how individual cells traverse the differentiation cascade as well as discern how specific cell states can re-emerge in breast cancer and contribute to metastasis.
Figure 1. Historical model of the mammary epithelial cell differentiation hierarchy. Multipotent mammary stem cells (MaSCs) have the ability to self-renew and give rise to all downstream lineages. Bipotent stem cells give rise to committed unipotent progenitor populations that maintain basal/myoepithelial and luminal lineages. These unipotent cells are responsible for maintaining mature basal/myoepithelial and luminal lineages in the postnatal gland. Two luminal progenitor populations have been reported, one giving rise to alveolar cells during pregnancy and the other responsible for maintaining the mature ductal cells
Single-cell sequencing has emerged as a powerful tool capable of uncovering the transcriptional heterogeneity of diverse cell types within a tissue. Although transcriptome profiling of single cells was reported over 20 years ago[48,49], only recently has it become a widely utilized technique. Technological advances in areas such as microfluidics and in situ barcoding have allowed for high throughput, cost effective sequencing of hundreds to thousands of cells within a given sample. Furthermore, the ability to analyze large numbers of cells at high resolution has mandated the development of innovative approaches for downstream data analysis. Software tools have recently become available to infer evolutionary history and lineage relationships within heterogeneous cell populations.
Recently, several independent laboratories have utilized scRNA-seq to identify epithelial subpopulations and their differentiation trajectory within the mammary gland[11-14]. Together these studies have examined the complete mammary developmental cycle encompassing embryonic through adult stages and pregnancy through post-weaning. The transcriptional landscape generated from these studies has provided critical insight into epithelial differentiation, transforming the way in which we view the mammary epithelial hierarchy. The new model emerging from these studies emphasizes the fluid and gradual progression of epithelial differentiation and reveals a greater degree of cellular heterogeneity than had previously been acknowledged.
Prior to birth, sc-RNA sequencing has shown that the epithelium consists of a continuum of primitive cellular states in which unique subpopulations could not be detected. This supports previous findings generated from gene expression analysis of bulk fMaSCs, demonstrating that fMaSCs are comprised of a single CD24+CD49f+ population. During embryogenesis, the majority of mammary epithelial cells exhibit dual expression of both luminal and myoepithelial signature genes[12,32,52]. This hybrid gene expression signature is associated with multipotent potential and a loss of dual lineage expression as cells become restricted to unipotent potential.
At birth, many mammary epithelial cells maintain a hybrid gene expression profile. This dual lineage signature may include rare bipotent progenitor populations previously observed by lineage tracing studies[28,38]. However, the majority of this population likely represents cells spanning various transitional stages of lineage restriction. Cells become increasingly more lineage restricted throughout puberty; however, the precise timing of the major bifurcation event between luminal and basal cells remains unclear. In work by Pal et al., scRNA-sequencing of mammary epithelium from 2 (pre-puberty), 5 (puberty), and 10 (adult) week-old mice demonstrated that before puberty, the mammary epithelial cell populations are largely homogeneous, expressing a basal (myoepithelial)-like gene signature. Upon puberty, a fraction of cells upregulate expression of luminal genes such as Epcam, and Keratins 8, 18, and 19, suggesting that bifurcation initiates at this time. However, Giraddi et al., specifically assessed changes in the early postnatal epithelium by examining glands from mice at day zero (PN0) and day 4 (PN4) after birth. In contrast to the report by Pal and colleagues, they found that two distinct cellular populations, one containing myoepithelial cells and the other containing luminal precursors, begin to emerge as early as PN4. While these studies confirm that the mammary epithelium remains in a primitive state until postnatal stages during which lineage restriction commences, the specific timing of those commitments remains unclear.
Sc-RNA sequencing of adult stages clearly demonstrate that mature epithelial cells are transcriptionally segregated into basal (myoepithelial) and luminal linages[11,12,53]. Within these overarching phenotypes, up to 15 transcriptionally distinct clusters (11 luminal and 4 basal) have been reported across nulliparous and reproductive stages. However, the majority of these subpopulations cannot be characterized by expression of a single gene, indicating that cells remain highly heterogeneous and many of these clusters may represent transient states of differentiation. The luminal lineage within the nulliparous gland undergoes multiple differentiation events. Luminal progenitors give rise to a series of intermediate states that bifurcate to generate a hormone responsive lineage (ER+) and a secretory alveolar lineage (ER-)[11,12,14]. Upon pregnancy, ER- cells undergo massive expansion to produce the alveolar cells needed for lactation. The expanded population expresses genes associated with proliferation, fatty-acid transport, and lipid biosynthetic and homeostatic processes. Recent lineage tracing studies have proposed that ER+ and ER- luminal lineages arise from independent progenitor populations[54,55]. In work by Wang et al., PROM1+ progenitors were found to give rise to ER+ luminal cells whereas SOX9+ progenitors give rise to ER- cells. Thus, as heterogeneous luminal populations undergo gradual differentiation, the acquisition of SOX9 or PROM1 expression may be a major determinant of further luminal cell fate specification.
The sc-RNA-sequencing studies reported to date have not detected a single transcriptionally distinct multipotent stem cell population within the adult murine gland[11,12,53]. However, dormant populations that possess stem potential under specific environmental cues may exist. In a study by Sun et al. mammary epithelial stem cell markers including Bcl11b, Tspan8, and Procr were detected in the adult gland; however, expression on these genes was mutually exclusive, indicating each labels an independent population. Further analysis of a novel population marked by Cdh5+, revealed increased stem activity within this group of cells as demonstrated by both mammosphere and in vivo transplantation assays.
In human breast, bipotent progenitor populations within the adult gland have been identified in cells isolated from reduction mammoplasties. Using Monocle software analysis of sc-RNA sequencing data collected from seven individuals, Nguyen et al., mapped both luminal and myoepithelial lineages to a common cellular origin. These cells exhibited a basal gene expression signature as well as expression of ZEB1. This contrasts with the inability to identify subpopulations of bipotent progenitors using sc-RNA sequencing of mouse mammary glands. Discrepancies between mouse and human studies may be due to differences in the methodologies employed to generate scRNA-seq profiles [Table 1]. These include epithelial cell isolation methods, read depth, and data analysis and interpretation. For example, in work by Giraddi et al. (2018) and Bach et al. (2017), sc-RNA sequencing clustering results were based on epithelial cells that were positive for EPCAM expression. Thus, there remains the possibility that an EPCAM-stem cell population was overlooked in the murine studies. Additionally, cells double positive for the major luminal and basal cell markers, cytokeratin 8 and 14 respectively, are detected within the basal epithelial layer of human adult breast tissue. This is in contrast with mouse basal epithelium where double positive cells are not detected after birth. This difference in dual cytokeratin expression coupled with the interpretation that cells expressing both luminal and basal markers represent bipotent progenitors may contribute to the differences between mouse and human scRNA-seq findings. Furthermore, evidence from transplantation assays suggests that stem populations are widely distributed within the gland[16,17,19]; however, it is unknown if specific subpopulations are spatially restricted. In this regard, it has been reported that a group of highly quiescent MaSCs, marked by Tspan8 expression, are located within proximal regions of the adult murine gland, closest to the nipple. Thus, discrepancies between data sets may also arise depending on area of the mammary gland from which the sequenced cells were collected.
Comparison of scRNA-seq methodologies used to examine the mammary epithelial hierarchy
|Author||Species and developmental time points examined||Cell isolation method||scRNA-seq Platform||Number of cells sequenced||Sequencing system||Read coverage (x) or Number of reads||Genes detected per cell|
|Bach et al. 2017||Mouse:
Day 14.5 Gestation
Day 6 Lactation
11 days post-involution
|10 × Genomics Chromium||25,010||Illumina HiSeq 2,500||~ 162 million||2,118|
|Pal et al. 2017||Mouse:
Pre-puberty (2 weeks)
Puberty (5 weeks)
Adult (10 weeks)
|10 × Genomics Chromium||3,308||NextSeq 500||Not reported||> 1,500|
|Fluidigm C1||460||Illumina HiSeq 2000||2 million||8000|
|Giraddi et al. 2018||Mouse:
Embryonic day 16
Embryonic day 18
Post-natal day 4
|10 × Genomics Chromium||6,424||Illumina HiSeq 2,500||5000-50,000||500-2000|
|Fluidigm C1||262||Illumina HiSeq 2,500||> 1.5 million||4000-9000|
|Sun et al. 2018||Mouse:
Day 12 Gestation
|Fluidigm C1||318||Illumina HiSeq 2,500||> 1 million||12,688|
|Nguyen et al. 2018||Human:
Mammoplasties from age matched post-pubertal and pre-pubertal females
|CD49fhiEpCAM+ (Basal)||Fluidigm C1||868||Illumina HiSeq 2,500||~ 1.6 million||~ 4,500|
|CD49f+EpCAMhi (luminal)||10 × Genomics Chromium||24,646||Illumina HiSeq 4000||~ 60,000|| > 500
Luminal progenitors that are first identified as mixed lineage cells gradually gain luminal orientation during postnatal development. By puberty, a distinct luminal progenitor population can be detected. However, scRNA-sequencing has further uncovered novel subcategories within the greater luminal progenitor population. For example, prior to puberty, CD55 is a marker of rare basal cells that give rise to a subset of luminal progenitors during adulthood. CD55+ cells within the luminal fraction exhibit increased colony forming ability compared to CD55- cells. In addition, populations positive for both CD55 and the progenitor marker, CD14 possess an even greater colony-forming ability. Luminal progenitors acquire additional complexity during pregnancy. Single-cell analysis of luminal progenitors in the nulliparous adult compared to post-involution glands revealed permanent transcriptomic alterations within a subset of post-involution cells. Following involution, luminal progenitors retain expression of genes associated with lactation and milk production. It is likely that this population overlaps with the previously described PI-MEC population.
The ability of mammary epithelial cells to express both luminal and myoepithelial signature genes at various stages of maturation and reproduction suggests that a fraction of cells may maintain lineage plasticity throughout life[12,14,52,53]. Studies investigating chromatin accessibility within the mammary epithelium at fetal and adult stages have demonstrated that bi-lineage gene expression is achievable due to epigenetic regulation of genes associated with multiple lineages[12,58]. In a recent study by Dravis et al., gene expression profiles and open chromatin regions were compared between fetal (f)MECs, luminal progenitors, and mature lineages using RNA-seq and ATAC-seq, respectively. fMaSCs expressed both canonical luminal and basal genes and gene loci were associated with open chromatin, reflecting their immature, multi-lineage state (consistent with Giraddi et al.). While, the chromatin architecture of adult myoepithelial cells resembled that of fMaSCs in which luminal and myoepithelial gene loci were associated with open chromatin; only genes associated with the myoepithelial lineage were expressed. These results indicate that although myoepithelial cells express lineage-restricted genes, they retain the flexibility to adopt a luminal identity similar to multipotent fMaSCs. This may explain the multipotent capacity of basal cells observed upon transplantation[22,23]. Several epigenetic factors have been shown to play a role in defining mammary epithelial cell fate. These include the histone methylation reader PYGO2, the polycomb-repressive complex-1 member, BMI1[60,61], and the polycomb-repressive complex 3 member, EZH2. BMI1 and EZH2 also play key roles in metastasis of breast cancers, likely due to their activation of more primitive cell fates[63-66]. Yet how these players collaborate to orchestrate the epigenetic continuum that drives fate determination and heterogeneity is not fully understood.
Differential chromatin accessibility of both myoepithelial and luminal genes in basal cells may explain the ability of transcriptional drivers to induce cell fate switching in committed cells. In a recent example, overexpression of Notch1 in Smooth muscle actin (SMA)+ or K5+ myoepithelial cells, was sufficient to commit cells to a luminal fate. Complementary to this, overexpression of p63 in luminal cells was sufficient to reprogram them into a myoepithelial state. Although the mechanism of p63-mediated cell fate switching is unknown, chromatin remodeling is likely required as mature luminal cells exhibit repressed chromatin at basal gene loci. SOX10 has also emerged as a major transcriptional regulator of epithelial cell fate. Within fMaSCs, SOX10 motifs are enriched at accessible chromatin regions flanking highly expressed genes. Furthermore, tumors expressing high SOX10 levels exhibit neural crest-like features and high SOX10 expression is correlated with the aggressive, basal-like breast cancer subtype. As embryonic neural crest cells are multipotent and highly mobile, SOX10 reprograming may lead to primitive cellular states that could contribute to aggressive tumor phenotypes.
Cumulative evidence from the recent papers utilizing sc-RNA-sequencing has revealed that the historical models by which the mammary epithelial hierarchy has been traditionally represented do not accurately portray the complexity of the system[11-14,53]. Previous models in which cellular states are depicted as discrete populations that differentiate along restricted paths is an oversimplification [Figure 1]. Instead, recent results indicate that heterogeneous cell populations gradually advance, and likely also retreat, through a differentiation trajectory[11,12,53]. This is supported by the work from several groups, reporting a vast array of epithelial transcriptional profiles throughout the stages of mammary gland development[11-14,53]. Similar findings have redefined the hierarchical visualization of hematopoiesis. Using the new visual representation of hematopoietic differentiation suggested by Laurenti and Gottgens as an example, we propose a similar paradigm to represent the mammary epithelial hierarchy [Figure 2]. As depicted in this model, groups of heterogeneous epithelial cells traverse through the differentiation landscape, passing through cellular states that have been previously defined including fMaSCs, luminal and myoepithelial progenitors, and mature luminal and myoepithelial lineages. The newly reported heterogeneity of these populations implies that each cell may take a slightly different transcriptional path from the next as it undergoes differentiation. Additionally, the array of cellular states detected within the gland at any one time suggests that variation in temporal dynamics of differentiation may exist between individual cells. Moreover, this model implies significant forward and reverse plasticity of cell states that could be impacted by the microenvironment. This could partially explain historical difficulty in attempting to isolate and characterize specific subpopulations; however, further studies are needed to test these predictions.
Figure 2. A progression model of mammary epithelial cell fate specification. Mammary epithelial differentiation occurs as cells pass through a continuous cascade of transcriptional states, with the majority of cells becoming increasingly more lineage-restricted throughout life and in response to pregnancy. Progressive lineage restriction is represented by gradual changes of the colors of cells (red, orange, yellow, green and blue). As cells differentiate, they pass through previously defined phenotypic compartments [fetal mammary stem cells (fMaSCs), Luminal and myoepithelial (myo-epi) progenitors, mature luminal (ML) ductal and alveolar cells and mature myoepithelial cells (rectangles)]. Grayscale circles (ranging from white to black) on the horizontal lines represent the various epigenetic and transcriptomic states found within each phenotypic compartment with the number and different shading patterns representing the degree of heterogeneity within each major phenotypic group. Arching and angled lines connecting those circles reflect the various transcriptional and epigenetic trajectories that cells pass through before eventually terminating into mature, differentiated lineages (myoepithelial, luminal, and alveolar). Cells at various points along the differentiation cascade may serve as tumor initiators, leading to inter- and intra-tumoral heterogeneity. Figure adapted from Giraddi et al.
Breast cancer is an amalgam of diseases that exhibit both inter- and intra-tumoral heterogeneity. Gene expression profiling has led to the classification of five overarching subtypes, including luminal A, luminal B, HER2+, basal-like, and claudin-low[71-74]. However, within each subtype, tumors can exhibit further variability in gene expression, molecular function, and drug susceptibility conveying distinct patient outcomes. The heterogenic nature of breast cancer is thought to arise from the combination of cellular origin, genetic and epigenetic changes, and environmental context[1,10,70,75].
Molecular profiling studies have revealed that the gene expression patterns of cancer subtypes align with those of normal mammary epithelial cell lineages[76,77]. This suggests that tumor subtypes may originate from distinct mammary epithelial subpopulations. It is widely speculated that epithelial stem/progenitor populations may serve as tumor initiating cells since their longevity and ability to self-renew affords the accumulation of genetic mutations. This has been supported by gene expression profiling studies. For example, the gene expression profile of claudin-low tumors most closely resembles that of MaSC/basal epithelial cells, whereas basal-like tumors more closely resemble luminal progenitors[45,72,76,78]. Further evidence supporting luminal progenitors as the cell of origin for the basal-like subtype has been shown in tumors harboring BRCA-1 mutations[76,77,79]. More recently, sc-RNA-sequencing of human breast epithelial cells demonstrated that cells expressing luminal progenitor genes correlated most highly with basal-like breast cancers from the METABRIC dataset. Despite these findings, not all cancer subtypes have been correlated with a cell of origin and the inability to isolate pure populations of mammary stem/progenitor cells has impeded functional testing. Moreover, it is also feasible that cancers acquire the transcriptional profiles of specific cell states due to the mutations that drive tumor development or progression and that the cells of origin may not be consistent with final tumor deposition. In this case, identifying the contextual cues that promote normal epithelial hierarchical progression would likely reveal pathways that are also activated or suppressed in tumors that promote their phenotypes.
The new view of the mammary epithelial hierarchy suggests that the continued search for a tumor cell of origin based on molecular profiles that have historically been used to characterize the major epithelial populations (mature luminal, myoepithelial, luminal progenitor ext.) is unlikely to reveal discreet cells from which different breast cancer subtypes emanate. The complexity of the normal epithelium suggests a similar heterogeneity is present in the primary and secondary tumors that arise from it. Indeed, a study comparing sc-RNA-seq profiles generated from triple-negative patient tumors revealed that a subset of cells within each tumor expressed profiles resembling a spectrum of epithelial differentiation statuses. It now can be presumed that cells residing anywhere along the continuous differentiation trajectory may acquire tumor-initiating mutations [Figure 2] and that the location of the cell within the differentiation hierarchy will determine its transcriptional status and influence whether acquired mutations will be oncogenic. This is most consistent with the inability to identify a single lineage that contributes to HER2/Neu-induced tumors because amplification or overexpression of the ERBB2 gene may occur in an array of cell states rather than a single cell of origin for all tumors[81,82]. Similar findings have been reported in leukemia, in which progenitor cells rely on the activation of the β-catenin pathway for oncogenic transformation. Although it is noted that some populations may be more susceptible to mutagenic events than others[84,85], the possibility that phenotypically similar tumors arise from cells residing in the exact same cellular state is slim. This notion is supported by the transcriptional heterogeneity observed across breast cancer patient tumors within the same subtype. Several genetic mutations have been associated with the stratification of breast cancer subtypes. These include mutations in GATA3, RUNX1, or NCOR1 which are associated with the Luminal A subtype; whereas mutations in the tumor suppressor retinoblastoma, RB1, or CDH1 have been associated with luminal B tumors[10,87,88]. Although these mutations are associated with the acquisition of a specific breast cancer subtype, it is likely the combination of the preexisting epigenetic and transcriptional status of the cell along with the transforming mutations that explain tumor heterogeneity and give rise to the variable therapeutic responses observed between patients.
Accumulating evidence suggests that a rare population of stem-like progenitor cells are responsible for tumor initiation, progression, metastasis, and therapeutic resistance. These cells are referred to as cancer stem cells (CSCs) or “tumor initiating cells” and share many characteristics with normal mammary epithelial stem cells including the ability to self-renew and differentiate. Additionally, CSCs share gene expression profiles that closely resemble those of primitive mammary epithelial cells, including the expression of stem and EMT associated genes. This has led to many to posit that CSCs arise from mammary stem/progenitor cells, however definitive evidence for the origin of CSCs in breast cancer is lacking. The extensive heterogeneity within the mammary epithelium that has been revealed by scRNA-seq makes it reasonable to presume that many cellular states have the potential to give rise to CSCs. Furthermore, the unique transcriptional and epigenetic context of an epithelial cell prior to transformation likely determines the aggressive nature of the arising tumor, including its ability to metastasize. Indeed, many studies have identified transcriptional regulators of normal mammary development that also encourage aggressive behaviors during cancer. This is the case for transcriptional drivers of stem cell and EMT programs including SOX9, SLUG, and SOX10. In a study by Guo et al., expression of both SOX9 and SLUG was found to be necessary for maintenance of the MaSC phenotype. Moreover, coexpression of SOX9 and SLUG in breast cancer cell lines was sufficient to induce stem cell-like properties, including activation of an EMT program and metastasis-seeding abilities. Similar results were observed for SOX10. Binding of SOX10 within regions of open chromatin that flank genes associated with EMT or with embryonic neural crest cells (NCC) was found to promote migratory and invasive behavior within mammary tumor cells. Thus, it is possible that epithelial cells expressing high levels of transcriptional drivers of EMT and stem-like states at the time of transformation may be predisposed to malignant phenotypes. Likewise, the chromatin accessibility of genes associated with processes such as EMT and NCC may ‘prime’ cells for aggressive behaviors associated with CSC phenotypes.
The formation of CSCs and their ability to seed new tumors at distant sites is also highly dependent on the surrounding microenvironment[93-95]. Notably, contextual stimuli such as exposure to TGFβ, can initiate stem cell and EMT programs in normal mammary epithelial cells that then promote their transition to a CSC state. The resulting CSCs remain highly susceptible to surrounding cues. This is especially relevant in the context of metastasis, in which Paget’s “seed and soil” hypothesis proposes that a tumor-permissive environment is required for disseminating tumor cells to engraft and metastasize into new tissue sites[93,96]. As breast cancer arises from the normal epithelium, the contextual cues regulating normal epithelial behavior are likely mirrored in breast cancer cells. Therefore, understanding the reciprocal interactions between epithelial cells and their niche should also inform the mechanisms underlying tumor initiation and metastasis.
In this regard, the recent studies using sc-RNA-sequencing outside the context of disruptive transplantation assays, have revealed that epithelial populations may be more highly influenced by the surrounding environmental cues than previously appreciated. As discussed earlier, previous studies using transplantation assays demonstrated the multipotent potential of rare mammary epithelial cells within the adult gland[22,23]. However, lineage tracing analyses were unable to replicate these results and recent sc-RNA-sequencing has further supported the lack of multipotent cells during adulthood. Thus, it has become increasingly apparent that introduction of isolated cells into wounded stromal environment may elicit cellular phenotypes that are not supported under normal physiological conditions. A similar phenomenon is likely to occur with breast cancer initiation and progression. While acquisition of an “oncogenic” mutation may occur in varied cells along the differentiation hierarchy, its impact on tumorigenesis and progression may rely upon the stromal microenvironment, particularly an environment that is reflective of wound repair such as that which occurs with high breast density, obesity, and in post-lactational remodeling. The same is also likely to be true with metastasizing cancer cells. Instructional interactions between the epithelium and surrounding stroma are well established. This was most exquisitely demonstrated in extensive studies by Smith and colleagues revealing that the local microenvironment of the mammary gland is sufficient to reprogram progenitor cells from completely different tissues (including testicular and embryonic stem cells) into mammary progenitors that were capable of forming normal mammary gland structures[100-104]. Furthermore, alterations in the microenvironment may influence cellular transformation. In work by Barcellos-Hoff and Ravani, irradiating the cleared fat pads of recipient mice, lead to alterations in growth factor activity and extracellular matrix composition. Following irradiation, implantation of non-tumorigenic mammary epithelial cells resulted in tumor formation more readily in irradiated mice. Conversely, epithelial tumor cells injected into non-irradiated mammary fat pads reverted to form normal mammary gland ducts rather than tumors.
The effects of contextual cues on the epigenetic and transcriptional status of major transcriptional regulators has also been recently highlighted. Although SOX10 is highly expressed in cells isolated from MMTV-PyMT tumors, this expression was lost when cells were grown in 2D culture conditions. The loss of SOX10 expression was associated with a switch from open to closed chromatin at the SOX10 locus, emphasizing the role of environmental cues on epigenetic status. This is not a new concept, as changes in the epigenetic landscape are known to occur in response to hormonal exposure[62,107].
The scRNA-sequencing studies reviewed here have provided a detailed transcriptional map of mammary epithelial differentiation. However, these data indicate that previously underappreciated complexity undoubtedly underlies the regulation of the heterogeneous transcriptional networks observed. The unique transcriptional status within each epithelial cell is likely correlated with a unique epigenetic profile and it is anticipated that these profiles have a considerable level of plasticity that is dictated by the microenvironment and, in the case of tumors, mutational status. Studies examining the chromatin accessibility within bulk isolations of epithelial populations (MaSCs, basal, luminal and luminal progenitors) have demonstrated that each population is associated with distinctive active and repressed regions, many of which are associated with functional consequences and specific tumor subtypes[12,58,107]. The emerging ability to simultaneously assess both the transcriptional and epigenetic status at a single cell level should allow for clearer understanding of the lineage relationships in the normal gland and in tumors as well as reveal whether tumor cells can readily transition through states as has been observed in the normal breast. It is important to note, however, that while single cell sequencing technologies can provide high resolution maps of cellular states, they cannot make definitive predictive connections to the resulting cellular fates. The use of newly developed technologies, such as “CellTagging” that allow for early transcriptional statuses to be linked to cellular fates will allow for a clearer understanding of epithelial lineage relationships within the mammary gland. Furthermore, studies layering single cell technologies that assess transcriptional, epigenetic, and proteomic data within the normal gland, primary tumor, and associated metastatic lesions should provide important information regarding the extent of plasticity of primary and metastatic tumor cells that can be therapeutically leveraged to promote differentiation of breast cancers into less malignant states.
Conceptualized and drafted the manuscript: Anstine LJ
Revised the manuscript for important intellectual content: Keri RAvailability of data and materials
Not applicable.Financial support and sponsorship
None.Conflicts of interest
Both authors declared that there are no conflicts of interest.Ethical approval and consent to participate
Not applicable.Consent for publication
© The Author(s) 2019.
1. Wiseman BS, Werb Z. Stromal effects on mammary gland development and breast cancer. Science 2002;296:1046-9.DOIPubMedPMC
2. Hu GF, Li LZ, Xu W. Extracellular matrix in mammary gland development and breast cancer progression. Frontiers in Laboratory Medicine 2017;1:36-9.DOI
3. Watson CJ, Khaled WT. Mammary development in the embryo and adult: a journey of morphogenesis and commitment. Development 2008;135:995-1003.DOIPubMed
4. Gjorevski N, Nelson CM. Integrated morphodynamic signalling of the mammary gland. Nat Rev Mol Cell Biol 2011;12:581-93.DOIPubMed
5. Richert MM, Schwertfeger KL, Ryder JW, Anderson SM. An atlas of mouse mammary gland development. J Mammary Gland Biol Neoplasia 2000;5:227-41.DOIPubMed
6. Wagner KU, Boulanger CA, Henry MD, Sgagias M, Hennighausen L, et al. An adjunct mammary epithelial cell population in parous females: its role in functional adaptation and tissue renewal. Development 2002;129:1377-86.PubMed
7. Smith GH. Experimental mammary epithelial morphogenesis in an in vivo model: evidence for distinct cellular progenitors of the ductal and lobular phenotype. Breast Cancer Res Treat 1996;39:21-31.DOIPubMed
8. Fu N, Lindeman GJ, Visvader JE. The mammary stem cell hierarchy. Curr Top Dev Biol 2014;107:133-60.DOIPubMed
9. Rodilla V, Fre S. Cellular plasticity of mammary epithelial cells underlies heterogeneity of breast cancer. Biomedicines 2018;6:E103.DOIPubMedPMC
10. Visvader JE, Stingl J. Mammary stem cells and the differentiation hierarchy: current status and perspectives. Genes Dev 2014;28:1143-58.DOIPubMedPMC
11. Bach K, Pensa S, Grzelak M, Hadfield J, Adams DJ, et al. Differentiation dynamics of mammary epithelial cells revealed by single-cell RNA sequencing. Nat Commun 2017;8:2128.DOIPubMedPMC
12. Giraddi RR, Chung CY, Heinz RE, Balcioglu O, Novotny M, et al. Single-cell transcriptomes distinguish stem cell state changes and lineage specification programs in early mammary gland development. Cell Rep 2018;24:1653-1666.e7.DOIPubMedPMC
13. Sun H, Miao Z, Zhang X, Chan UI, Su SM, et al. Single-cell RNA-Seq reveals cell heterogeneity and hierarchy within mouse mammary epithelia. J Biol Chem 2018;293:8315-29.DOIPubMedPMC
14. Nguyen QH, Pervolarakis N, Blake K, Ma D, Davis RT, et al. Profiling human breast epithelial cells using single cell RNA sequencing identifies cell diversity. Nat Commun 2018;9:2028.DOIPubMedPMC
15. Deome KB, Faulkin LJ Jr, 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-20.PubMed
16. Daniel CW, De Ome KB, Young JT, Blair PB, Faulkin LJ Jr. The in vivo life span of normal and preneoplastic mouse mammary glands: a serial transplantation study. Proc Natl Acad Sci U S A 1968;61:53-60.DOIPubMedPMC
17. Hoshino K. Morphogenesis and growth potentiality of mammary glands in mice. I. Transplantability and growth potentiality of mammary tissue of virgin mice. J Natl Cancer Inst 1962;29:835-51.DOIPubMed
18. Smith GH, Gallahan D, Zwiebel JA, Freeman SM, Bassin RH, et al. Long-term in vivo expression of genes introduced by retrovirus-mediated transfer into mammary epithelial cells. J Virol 1991;65:6365-70.PubMedPMC
19. Smith GH, Medina D. A morphologically distinct candidate for an epithelial stem cell in mouse mammary gland. J Cell Sci 1988;90:173-83.PubMed
20. Young LJ, Medina D, DeOme KB, Daniel CW. The influence of host and tissue age on life span and growth rate of serially transplanted mouse mammary gland. Exp Gerontol 1971;6:49-56.DOIPubMed
21. Kordon EC, Smith GH. An entire functional mammary gland may comprise the progeny from a single cell. Development 1998;125:1921-30.PubMed
22. Shackleton M, Vaillant F, Simpson KJ, Stingl J, Smyth GK, et al. Generation of a functional mammary gland from a single stem cell. Nature 2006;439:84-8.DOIPubMed
23. Stingl J, Eirew P, Ricketson I, Shackleton M, Vaillant F, et al. Purification and unique properties of mammary epithelial stem cells. Nature 2006;439:993-7.DOIPubMed
24. Shehata M, Teschendorff A, Sharp G, Novcic N, Russell IA, et al. Phenotypic and functional characterisation of the luminal cell hierarchy of the mammary gland. Breast Cancer Res 2012;14:R134.DOIPubMedPMC
25. 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.DOIPubMedPMC
26. Bai L, Rohrschneider LR. s-SHIP promoter expression marks activated stem cells in developing mouse mammary tissue. Genes Dev 2010;24:1882-92.DOIPubMedPMC
27. Plaks V, Brenot A, Lawson DA, Linnemann JR, Van Kappel EC, et al. Lgr5-expressing cells are sufficient and necessary for postnatal mammary gland organogenesis. Cell Rep 2013;3:70-8.DOIPubMedPMC
28. Rios AC, Fu NY, Lindeman GJ, Visvader JE. In situ identification of bipotent stem cells in the mammary gland. Nature 2014;506:322-7.DOIPubMed
29. de Visser KE, Ciampricotti M, Michalak EM, Tan DW, Speksnijder EN, et al. Developmental stage-specific contribution of LGR5(+) cells to basal and luminal epithelial lineages in the postnatal mammary gland. J Pathol 2012;228:300-9.DOI
30. Badders NM, Goel S, Clark RJ, Klos KS, Kim S, et al. The Wnt receptor, Lrp5, is expressed by mouse mammary stem cells and is required to maintain the basal lineage. PLoS One 2009;4:e6594.DOIPubMedPMC
31. 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-77.DOIPubMedPMC
32. Spike BT, Engle DD, Lin JC, Cheung SK, La J, et al. A mammary stem cell population identified and characterized in late embryogenesis reveals similarities to human breast cancer. Cell Stem Cell 2012;10:183-97.DOIPubMedPMC
33. Wang D, Cai C, Dong X, Yu QC, Zhang XO, et al. Identification of multipotent mammary stem cells by protein C receptor expression. Nature 2015;517:81-4.DOIPubMed
34. Van Keymeulen A, Rocha AS, Ousset M, Beck B, et al. Distinct stem cells contribute to mammary gland development and maintenance. Nature 2011;479:189-93.DOIPubMed
35. Davis FM, Lloyd-Lewis B, Harris OB, Kozar S, Winton DJ, et al. Single-cell lineage tracing in the mammary gland reveals stochastic clonal dispersion of stem/progenitor cell progeny. Nat Commun 2016;7:13053.DOIPubMedPMC
36. Wuidart A, Ousset M, Rulands S, Simons BD, Van Keymeulen A, et al. Quantitative lineage tracing strategies to resolve multipotency in tissue-specific stem cells. Genes Dev 2016;30:1261-77.DOIPubMedPMC
37. Boulanger CA, Wagner KU, Smith GH. Parity-induced mouse mammary epithelial cells are pluripotent, self-renewing and sensitive to TGF-beta1 expression. Oncogene 2005;24:552-60.DOIPubMed
38. 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.DOIPubMed
39. Song W, Wang R, Jiang W, Yin Q, Peng G, et al. Hormones induce the formation of luminal-derived basal cells in the mammary gland. Cell Res 2019; doi: 10.1038/s41422-018-0137-0.DOIPubMed
40. Morris RJ, Liu Y, Marles L, Yang Z, Trempus C, et al. Capturing and profiling adult hair follicle stem cells. Nat Biotechnol 2004;22:411-7.DOIPubMed
41. Rodriguez-Fraticelli AE, Wolock SL, Weinreb CS, Panero R, Patel SH, et al. Clonal analysis of lineage fate in native haematopoiesis. Nature 2018;553:212-6.DOIPubMedPMC
42. Rios AC, Fu NY, Cursons J, Lindeman GJ, Visvader JE. The complexities and caveats of lineage tracing in the mammary gland. Breast Cancer Res 2016;18:116.DOIPubMedPMC
43. Shehata M, van Amerongen R, Zeeman AL, Giraddi RR, Stingl J. The influence of tamoxifen on normal mouse mammary gland homeostasis. Breast Cancer Res 2014;16:411.DOIPubMedPMC
44. Yang X, Wang H, Jiao B. Mammary gland stem cells and their application in breast cancer. Oncotarget 2017;8:10675-91.DOIPubMedPMC
45. Lim E, Wu D, Pal B, Bouras T, Asselin-Labat ML, et al. Transcriptome analyses of mouse and human mammary cell subpopulations reveal multiple conserved genes and pathways. Breast Cancer Res 2010;12:R21.DOIPubMedPMC
46. Trejo CL, Luna G, Dravis C, Spike BT, Wahl GM. Lgr5 is a marker for fetal mammary stem cells, but is not essential for stem cell activity or tumorigenesis. NPJ Breast Cancer 2017;3:16.DOIPubMedPMC
47. Shi XS, Chakraborty P, Chaudhuri A. Unmasking tumor heterogeneity and clonal evolution by single-cell analysis. J Cancer Metastasis Treat 2018;4:47.DOI
48. Eberwine J, Yeh H, Miyashiro K, Cao Y, Nair S, et al. Analysis of gene expression in single live neurons. Proc Natl Acad Sci U S A 1992;89:3010-4.DOIPubMedPMC
49. Cano-Gauci DF, Lualdi JC, Ouellette AJ, Brady G, Iscove NN, et al. In vitro cDNA amplification from individual intestinal crypts: a novel approach to the study of differential gene expression along the crypt-villus axis. Exp Cell Res 1993;208:344-9.DOIPubMed
50. Svensson V, Vento-Tormo R, Teichmann SA. Exponential scaling of single-cell RNA-seq in the past decade. Nat Protoc 2018;13:599-604.DOIPubMed
51. Kuipers J, Jahn K, Beerenwinkel N. Advances in understanding tumour evolution through single-cell sequencing. Biochim Biophys Acta Rev Cancer 2017;1867:127-38.DOIPubMedPMC
52. Wuidart A, Sifrim A, Fioramonti M, Matsumura S, Brisebarre A, et al. Early lineage segregation of multipotent embryonic mammary gland progenitors. Nat Cell Biol 2018;20:666-76.DOIPubMedPMC
53. Pal B, Chen Y, Vaillant F, Jamieson P, Gordon L, et al. Construction of developmental lineage relationships in the mouse mammary gland by single-cell RNA profiling. Nat Commun 2017;8:1627.DOIPubMedPMC
54. Wang C, Christin JR, Oktay MH, Guo W. Lineage-biased stem cells maintain estrogen-receptor-positive and -negative mouse mammary luminal lineages. Cell Rep 2017;18:2825-35.DOIPubMedPMC
55. Van Keymeulen A, Fioramonti M, Centonze A, Bouvencourt G, Achouri Y, et al. Lineage-restricted mammary stem cells sustain the development, homeostasis, and regeneration of the estrogen receptor positive lineage. Cell Rep 2017;20:1525-32.DOIPubMedPMC
56. Cai S, Kalisky T, Sahoo D, Dalerba P, Feng W, et al. A quiescent Bcl11b high stem cell population is required for maintenance of the mammary gland. Cell Stem Cell 2017;20:247-60.e5.DOIPubMedPMC
57. Fu NY, Rios AC, Pal B, Law CW, Jamieson P, et al. Identification of quiescent and spatially restricted mammary stem cells that are hormone responsive. Nat Cell Biol 2017;19:164-76.DOIPubMed
58. Dravis C, Chung CY, Lytle NK, Herrera-Valdez J, Luna G, et al. Epigenetic and transcriptomic profiling of mammary gland development and tumor models disclose regulators of cell state plasticity. Cancer Cell 2018;34:466-82.e6.DOIPubMedPMC
59. Gu B, Sun P, Yuan Y, Moraes RC, Li A, et al. Pygo2 expands mammary progenitor cells by facilitating histone H3 K4 methylation. J Cell Biol 2009;185:811-26.DOIPubMedPMC
60. Liu S, Dontu G, Mantle ID, Patel S, Ahn NS, et al. Hedgehog signaling and Bmi-1 regulate self-renewal of normal and malignant human mammary stem cells. Cancer Res 2006;66:6063-71.DOIPubMedPMC
61. Pietersen AM, Evers B, Prasad AA, Tanger E, Cornelissen-Steijger P, et al. Bmi1 regulates stem cells and proliferation and differentiation of committed cells in mammary epithelium. Curr Biol 2008;18:1094-9.DOIPubMed
62. Pal B, Bouras T, Shi W, Vaillant F, Sheridan JM, et al. Global changes in the mammary epigenome are induced by hormonal cues and coordinated by Ezh2. Cell Rep 2013;3:411-26.DOIPubMed
63. Hoenerhoff MJ, Chu I, Barkan D, Liu ZY, Datta S, et al. BMI1 cooperates with H-RAS to induce an aggressive breast cancer phenotype with brain metastases. Oncogene 2009;28:3022-32.DOIPubMedPMC
64. Paranjape AN, Balaji SA, Mandal T, Krushik EV, Nagaraj P, et al. Bmi1 regulates self-renewal and epithelial to mesenchymal transition in breast cancer cells through Nanog. BMC Cancer 2014;14:785.DOIPubMedPMC
65. Moore HM, Gonzalez ME, Toy KA, Cimino-Mathews A, Argani P, et al. EZH2 inhibition decreases p38 signaling and suppresses breast cancer motility and metastasis. Breast Cancer Res Treat 2013;138:741-52.DOIPubMedPMC
66. Crea F, Fornaro L, Bocci G, Sun L, Farrar WL, et al. EZH2 inhibition: targeting the crossroad of tumor invasion and angiogenesis. Cancer Metastasis Rev 2012;31:753-61.DOIPubMed
67. Lilja AM, Rodilla V, Huyghe M, Hannezo E, Landragin C, et al. Clonal analysis of Notch1-expressing cells reveals the existence of unipotent stem cells that retain long-term plasticity in the embryonic mammary gland. Nat Cell Biol 2018;20:677-87.DOIPubMed
68. Mayor R, Theveneau E. The neural crest. Development 2013;140:2247-51.DOIPubMed
69. Laurenti E, Gottgens B. From haematopoietic stem cells to complex differentiation landscapes. Nature 2018;553:418-26.DOIPubMedPMC
70. Visvader JE. Cells of origin in cancer. Nature 2011;469:314-22.DOIPubMed
71. Prat A, Perou CM. Deconstructing the molecular portraits of breast cancer. Mol Oncol 2011;5:5-23.DOIPubMedPMC
72. Prat A, Parker JS, Karginova O, Fan C, Livasy C, et al. Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res 2010;12:R68.DOIPubMedPMC
73. Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A 2001;98:10869-74.DOIPubMedPMC
74. Neve RM, Chin K, Fridlyand J, Yeh J, Baehner FL, et al. A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell 2006;10:515-27.DOIPubMedPMC
75. Easwaran H, Tsai HC, Baylin SB. Cancer epigenetics: tumor heterogeneity, plasticity of stem-like states, and drug resistance. Mol Cell 2014;54:716-27.DOIPubMedPMC
76. Lim E, Vaillant F, Wu D, Forrest NC, Pal B, et al. Aberrant luminal progenitors as the candidate target population for basal tumor development in BRCA1 mutation carriers. Nat Med 2009;15:907-13.DOIPubMed
77. Molyneux G, Geyer FC, Magnay FA, McCarthy A, Kendrick H, et al. BRCA1 basal-like breast cancers originate from luminal epithelial progenitors and not from basal stem cells. Cell Stem Cell 2010;7:403-17.DOIPubMed
78. Herschkowitz JI, Simin K, Weigman VJ, Mikaelian I, Usary J, et al. Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome Biol 2007;8:R76.DOIPubMedPMC
79. Proia TA, Keller PJ, Gupta PB, Klebba I, Jones AD, et al. Genetic predisposition directs breast cancer phenotype by dictating progenitor cell fate. Cell Stem Cell 2011;8:149-63.DOIPubMedPMC
80. Karaayvaz M, Cristea S, Gillespie SM, Patel AP, Mylvaganam R, et al. Unravelling subclonal heterogeneity and aggressive disease states in TNBC through single-cell RNA-seq. Nat Commun 2018;9:3588.DOIPubMedPMC
81. Ferrari A, Sertier AS, Vincent-Salomon A, Pivot X, Pauporte I, et al. A phenotypic and mechanistic perspective on heterogeneity of HER2-positive breast cancers. Mol Cell Oncol 2016;3:e1232186.DOIPubMedPMC
82. Vaillant F, Asselin-Labat ML, Shackleton M, Forrest NC, Lindeman GJ, et al. The mammary progenitor marker CD61/beta3 integrin identifies cancer stem cells in mouse models of mammary tumorigenesis. Cancer Res 2008;68:7711-7.DOIPubMed
83. Wang Y, Krivtsov AV, Sinha AU, North TE, Goessling W, et al. The Wnt/beta-catenin pathway is required for the development of leukemia stem cells in AML. Science 2010;327:1650-3.DOIPubMedPMC
84. Thompson EG, Fares H, Dixon K. BRCA1 requirement for the fidelity of plasmid DNA double-strand break repair in cultured breast epithelial cells. Environ Mol Mutagen 2012;53:32-43.DOIPubMed
85. Magee JA, Piskounova E, Morrison SJ. Cancer stem cells: impact, heterogeneity, and uncertainty. Cancer Cell 2012;21:283-96.DOIPubMedPMC
86. Chung W, Eum HH, Lee HO, Lee KM, Lee HB, et al. Single-cell RNA-seq enables comprehensive tumour and immune cell profiling in primary breast cancer. Nat Commun 2017;8:15081.DOIPubMedPMC
87. Banerji S, Cibulskis K, Rangel-Escareno C, Brown KK, Carter SL, et al. Sequence analysis of mutations and translocations across breast cancer subtypes. Nature 2012;486:405-9.DOIPubMedPMC
88. Ellis MJ, Perou CM. The genomic landscape of breast cancer as a therapeutic roadmap. Cancer Discov 2013;3:27-34.DOIPubMedPMC
89. Shipitsin M, Polyak K. The cancer stem cell hypothesis: in search of definitions, markers, and relevance. Lab Invest 2008;88:459-63.DOIPubMedPMC
90. Mani SA, Guo W, Liao MJ, Eaton EN, Ayyanan A, et al. The epithelial-mesenchymal transition generates cells with properties of stem cells. Cell 2008;133:704-15.DOIPubMedPMC
91. Sin WC, Lim CL. Breast cancer stem cells-from origins to targeted therapy. Stem Cell Investig 2017;4:96.DOIPubMedPMC
92. Guo W, Keckesova Z, Donaher JL, Shibue T, Tischler V, et al. Slug and Sox9 cooperatively determine the mammary stem cell state. Cell 2012;148:1015-28.DOIPubMedPMC
93. Psaila B, Lyden D. The metastatic niche: adapting the foreign soil. Nat Rev Cancer 2009;9:285-93.DOIPubMedPMC
94. Chaffer CL, Marjanovic ND, Lee T, Bell G, Kleer CG, et al. Poised chromatin at the ZEB1 promoter enables breast cancer cell plasticity and enhances tumorigenicity. Cell 2013;154:61-74.DOIPubMedPMC
95. Lau EY, Ho NP, Lee TK. Cancer stem cells and their microenvironment: biology and therapeutic implications. Stem Cells Int 2017;2017:3714190.DOIPubMedPMC
96. Paget S. The distribution of secondary growths in cancer of the breast. 1889. Cancer Metastasis Rev 1989;8:98-101.DOIPubMed
97. Huo CW, Hill P, Chew G, Neeson PJ, Halse H, et al. High mammographic density in women is associated with protumor inflammation. Breast Cancer Res 2018;20:92.DOIPubMedPMC
98. Mentoor I, Engelbrecht AM, van Jaarsveld PJ, Nell T. Chemoresistance: intricate interplay between breast tumor cells and adipocytes in the tumor microenvironment. Front Endocrinol (Lausanne) 2018;9:758.DOIPubMedPMC
99. Martinson HA, Jindal S, Durand-Rougely C, Borges VF, Schedin P. Wound healing-like immune program facilitates postpartum mammary gland involution and tumor progression. Int J Cancer 2015;136:1803-13.DOIPubMedPMC
100. Bruno RD, Smith GH. Reprogramming non-mammary and cancer cells in the developing mouse mammary gland. Semin Cell Dev Biol 2012;23:591-8.DOIPubMedPMC
101. Bruno RD, Boulanger CA, Rosenfield SM, Anderson LH, Lydon JP, et al. Paracrine-rescued lobulogenesis in chimeric outgrowths comprising progesterone-receptor-null mammary epithelium and redirected wild-type testicular cells. J Cell Sci 2014;127:27-32.DOIPubMedPMC
102. Boulanger CA, Bruno RD, Rosu-Myles M, Smith GH. The mouse mammary microenvironment redirects mesoderm-derived bone marrow cells to a mammary epithelial progenitor cell fate. Stem Cells Dev 2012;21:948-54.DOIPubMedPMC
103. Bruno RD, Fleming JM, George AL, Boulanger CA, Schedin P, et al. Mammary extracellular matrix directs differentiation of testicular and embryonic stem cells to form functional mammary glands in vivo. Sci Rep 2017;7:40196.DOIPubMedPMC
104. Boulanger CA, Bruno RD, Mack DL, Gonzales M, Castro NP, et al. Embryonic stem cells are redirected to non-tumorigenic epithelial cell fate by interaction with the mammary microenvironment. PLoS One 2013;8:e62019.DOIPubMedPMC
105. Barcellos-Hoff MH, Ravani SA. Irradiated mammary gland stroma promotes the expression of tumorigenic potential by unirradiated epithelial cells. Cancer Res 2000;60:1254-60.PubMed
106. Maffini MV, Calabro JM, Soto AM, Sonnenschein C. Stromal regulation of neoplastic development: age-dependent normalization of neoplastic mammary cells by mammary stroma. Am J Pathol 2005;167:1405-10.DOIPubMedPMC
107. Casey AE, Sinha A, Singhania R, Livingstone J, Waterhouse P, et al. Mammary molecular portraits reveal lineage-specific features and progenitor cell vulnerabilities. J Cell Biol 2018;217:2951-74.DOIPubMedPMC
108. Biddy BA, Kong WJ, Kamimoto K, Guo C, Waye SE, et al. Single-cell mapping of lineage and identity in direct reprogramming. Nature 2018;564:219-24.DOIPubMed
Ana Armada et al., Cancer Drug Resistance, 2019
Aiko Sueta et al., Cancer Drug Resistance, 2019
Silvia Dei et al., Cancer Drug Resistance, 2019
Pooja Murthy et al., Cancer Drug Resistance, 2019
Jose E. Belizario et al., Cancer Drug Resistance, 2019
Miki Mori et al., Cancer Drug Resistance, 2019
Philippa D. Darbre et al., Journal of Cancer Metastasis and Treatment, 2019
Dillon D. McCourt et al., Journal of Unexplored Medical Data, 2019
Saertje Verkoeijen et al., Journal of Cancer Metastasis and Treatment, 2019
Logeswari Ponnusamy et al., Cancer Drug Resistance, 2019