Epigenomics of mammary gland development
Breast Cancer Research volume 20, Article number: 100 (2018)
Differentiation of stem cells into highly specialised cells requires gene expression changes brought about by remodelling of the chromatin architecture. During this lineage-commitment process, the majority of DNA needs to be packaged into inactive heterochromatin, allowing only a subset of regulatory elements to remain open and functionally required genes to be expressed. Epigenetic mechanisms such as DNA methylation, post-translational modifications to histone tails, and nucleosome positioning all potentially contribute to the changes in higher order chromatin structure during differentiation. The mammary gland is a particularly useful model to study these complex epigenetic processes since the majority of its development is postnatal, the gland is easily accessible, and development occurs in a highly reproducible manner. Inappropriate epigenetic remodelling can also drive tumourigenesis; thus, insights into epigenetic remodelling during mammary gland development advance our understanding of breast cancer aetiology. We review the current literature surrounding DNA methylation and histone modifications in the developing mammary gland and its implications for breast cancer.
Lineage commitment in the mammary gland
The mammary gland is a dynamic tissue with rapid changes in tissue architecture occurring throughout the lifetime of the mammal in response to hormonal cues (reviewed in [1, 2]). The gland is comprised of an epithelial ductal tree embedded within a stromal fat pad comprised of a variety of cell types including adipocytes, fibroblasts, immune cells, lymphatic cells, and vascular cells that interact with each other to maintain a functional organ . At birth, the gland contains a rudimentary ductal structure. The presence of oestrogen at puberty causes the ducts to undergo branching morphogenesis, generating a ductal tree that invades the stromal fat pad. Ductal elongation is driven by proliferation of cap cells located at the tips of the terminal end buds (TEBs) [1, 2]. During pregnancy and lactation, progesterone and prolactin cause extensive secondary and tertiary side branching and the formation of alveolar units that produce and secrete milk. Weaning of the offspring initiates the process of involution, which essentially remodels the mammary gland back to the virgin state [1, 2].
The mammary ductal epithelium is comprised of two main cell lineages: the inner luminal population containing ductal and alveolar cells, and the outer basal population containing myoepithelial cells (Fig. 1). The basal population is enriched for cells capable of self-renewal and multi-lineage differentiation upon serial transplantation into cleared mammary fat pads. These cells, known as mammary stem cells (MaSCs) [3, 4], lack unique cell surface markers, complicating their purification from the bulk basal cell population (referred to as ‘basal’). The luminal compartment contains proliferative luminal progenitor and mature luminal cells. During pregnancy and lactation, luminal progenitors differentiate into alveolar cells via alveolar progenitors. Transplantation experiments support a model whereby a bipotent MaSC at the apex of a differentiation hierarchy gives rise to both myoepithelial and luminal lineages . The existence of bipotent MaSCs under physiological conditions is debated, with various lineage-tracing experiments yielding irreconcilable results. Some groups have found that adult basal cells give rise to both mature luminal and basal cells [5, 6]. Other groups have found that bipotent MaSCs only exist during embryonic development and that basal and luminal lineages of the adult gland are maintained by distinct pools of unipotent stem cells [7,8,9,10] (Fig. 1). The different conclusions may be due to reliance on different genetic reporters that mark distinct cell populations with discrete differentiation potentials.
It is increasingly accepted that cellular differentiation is not unidirectional and that ‘terminally differentiated’ cells may exhibit plasticity under certain conditions of stress, injury, or experimental stimuli . Indeed, differentiated myoepithelial and luminal cells have been shown to adopt stem-like properties when cultured ex vivo . The molecular mechanisms underlying this cellular plasticity are largely unknown.
Epigenetics and development
Extensive changes in gene expression are required for a stem cell to undergo lineage commitment and functional differentiation during development. Changes in gene expression are associated with heritable epigenetic modifications to DNA and chromatin without changes to the DNA sequence. There are several layers of epigenetic regulation involved in the moderation of gene expression, including DNA methylation, post-translational modification to histone tails, chromatin remodelling, and higher order chromosome organisation (Fig. 2). DNA methylation is central to transcriptional repression, genomic imprinting, X-chromosome inactivation, and suppression of repetitive genomic elements [13, 14]. The basic chromatin subunit is the nucleosome, which is made up of ~ 147 nucleotides of DNA wrapped around a core histone octamer made up of two of each of the histone proteins H2A, H2B, H3, and H4  (Fig. 2). Protruding N-terminal histone tails are subject to covalent post-translational modifications including acetylation, phosphorylation, ubiquitination, and methylation. The resulting ‘histone code’ ultimately influences gene transcription through multiple mechanisms (reviewed in ). The organisation and positioning of nucleosomes also determines which regions of the genome are active or inactive. Chromatin remodellers use ATP hydrolysis to move, destabilise, eject, or restructure nucleosomes to change chromatin accessibility and gene transcription .
Epigenetic regulation of differentiation in embryonic stem (ES) cells and the haematopoietic system is well characterised [18, 19]; however, far less is known about the epigenetic mechanisms underlying self-renewal and differentiation of tissue-specific epithelial stem and progenitor cells. Lineage-specific epigenetic programmes can become deregulated and result in oncogenesis . Understanding these processes under physiological circumstances provides insights into how these complex programmes are altered during carcinogenesis.
The majority of CpG dinucleotides in the mammalian genome are methylated with the exception of CG-dense regions located around transcriptional start sites, known as CpG islands. These occur in approximately 70% of annotated gene promoters . During development, DNA methyltransferase enzymes (DNMT1, DNMT3A, and DNMT3B) deposit and maintain methylation in a small subset of CpG island promoters causing transcriptional repression . Conversely, methylation in gene bodies correlates with transcriptional activation  (Fig. 2). Cellular differentiation is accompanied by progressive gain of CpG island methylation resulting in silencing of developmental and non-related lineage-specific genes . Demethylation of tissue-specific genes also defines cellular identity and is mediated by ten-eleven translocation (Tet) methylcytosine dioxygenases (Tet1, Tet2, Tet3) [23, 24].
DNA methylation patterns define mammary epithelial cell identity
Mammary epithelial identity is, in part, shaped by the DNA methylation landscape, which varies between cell types [25, 26]. For example, methylation profiles from mammary epithelial cells are more similar to skin cells than blood cells, in keeping with the ectodermal origin of mammary and skin epithelial cells .
There have been numerous efforts to analyse DNA methylation patterns of stem/progenitor cells and differentiated luminal and myoepithelial cells in mice and humans [26,27,28,29,30]. Stem/progenitor cells isolated from human mammary glands are hypomethylated compared with differentiated luminal and myoepithelial cell types. This suggests that DNA methylation increases as the cells undergo lineage restriction, complementing work in other stem cell populations . Genes that are hypomethylated and highly expressed in stem/progenitor cells include transcription factors with known roles in stem cell maintenance (e.g. HOXA1 and TCF7L1) . Conversely, genes with promoter methylation (silenced) in stem/progenitor cells and gene body methylation (activated) in mature luminal cells include luminal-driving transcription factors (e.g. GATA3). This implies that DNA methylation is important in regulating the expression of lineage-specific transcription factors . Comprehensive DNA methylation profiling of mouse mammary epithelial subpopulations has revealed similar patterns; when compared to basal cells, luminal cells have hypermethylation and silencing of stem and basal cell-related genes (e.g. Angptl2 and Krt5) and hypomethylation and activation of epithelial differentiation genes (e.g. Elf5, Cldn4, and Krt8) [25, 26, 30]. Thus, DNA methylation is important for controlling expression of transcription factors that define basal cells and luminal populations.
Murine studies have been particularly useful in determining the effects of ageing and pregnancy on DNA methylation in mammary epithelial cells [25, 26]. Pregnancy induces changes to the DNA methylome of both basal and luminal populations that persist throughout life [25, 26]. Genes involved in lactation and involution become hypomethylated in response to the pregnancy and are primed to respond robustly to subsequent pregnancies . An epigenetic memory of pregnancy may explain the protective effect of early pregnancy on breast cancer risk in humans later in life .
DNMT1 maintains mammary stem/progenitor cells
DNMT1 has high affinity for hemimethylated DNA and is responsible for restoring the original methylation pattern present before DNA replication . Its activity is required for the maintenance of adult stem cells including MaSCs  (Table 1). Expression of DNMT1 is similar in basal and luminal cells and increases during pregnancy . Mammary-specific deletion of DNMT1 in mice severely effects TEB development and ductal elongation, results in fewer proliferative Ki67+ mammary epithelial cells, and mammary cells have reduced mammosphere-forming capacity when cultured ex vivo . Correspondingly, inhibition of DNMT activity by administering 5-azacitidine (5-AzaC) to mice decreases Cyclin D1 expression and an overall reduction in mammary cell numbers . Taken together, DNMT1-mediated DNA methylation is essential for maintenance of stem/progenitor cells in the mammary gland.
Histone modifications occur on lysine and arginine residues and regulate DNA accessibility, or act as protein docking sites for the initiation of downstream biological processes, including chromatin compaction, transcriptional regulation, and DNA repair . Histone acetylation is associated with transcriptional activation , while histone methylation can either activate (trimethylation of lysine 4 on histone 3; designated as H3K4Me3) or repress (H3K27Me3 or H3K9Me3) transcription [35, 36]. Bivalent promoters, containing both active H3K4Me3 and inactive H3K27Me3 marks, are often found at the promoters of key developmental genes and signal a repressed yet poised state, which allows for rapid activation or silencing of genes during differentiation . These were originally discovered in ES cells and have since been identified in adult stem cells. Histones flanking active enhancer regions are marked by H3K4Me1 and H3K27ac modifications . Histone modifications are most effectively mapped using chromatin immunoprecipitation followed by DNA sequencing (ChIP-Seq).
Rewriting the histone code through mammary gland development
Mammary subpopulations and regulatory elements governing transcription have been defined through the integration of ChIP-Seq data with transcriptional signatures [28, 39,40,41,42]. The diverse chromatin states that distinguish between mammary epithelial subpopulations are as extensive as those that distinguish epithelial subpopulations from developmentally unrelated stromal cells . The largest variations in chromatin state occurred in enhancer regions, although there were also significant variations in promoter regions . This shows that the chromatin state of cell-specific regulatory elements is a key determinant of cell type, even within the same epithelial lineage.
Histone methylation changes correlate with gene expression changes during lineage restriction. In the mouse mammary gland, luminal progenitor-defining genes  have higher H3K4Me3 coverage and lower H3K27Me3 coverage in luminal progenitor cells compared with basal cells . Genes within the mature luminal signature  show the same pattern when comparing mature luminal cells with luminal progenitors . Histone modifications are also altered upon pregnancy. Luminal cells differentiate into milk-secreting alveolar cells in preparation for lactation, which results in decreased repressive marks and increased active marks in key luminal differentiation and milk-production genes. Similar results were found in sorted human mammary epithelial populations [28, 40]. Interestingly, the genes repressed in luminal and basal subsets are often present in large regions enriched for H3K27Me3 marks (K27 blocs), which may allow for cell-type specific co-ordinated gene silencing . These studies demonstrate that distinct histone methylation profiles influence gene expression changes that direct basal to luminal progenitor differentiation and the maturation of luminal progenitors to luminal and alveolar cells (Fig. 1). Histone modifications also regulate gene expression within the heterogeneous basal population. A small population of quiescent MaSCs, marked by co-expression of LGR5 and TSPAN8, have recently been purified from the mouse mammary epithelium . These cells share similarities with quiescent stem cells in other tissues and have distinct H3K4Me3 and H3K27Me3 landscapes compared with the other basal cells .
Analysis of bivalent promoters in the mammary epithelial subpopulations has led to some interesting discoveries. Luminal progenitor cells have intermediate promoter features between basal and mature luminal cells. For example, promoters of genes involved in basal functions (such as extracellular matrix organisation) are marked by H3K4Me3 in basal cells, are bivalent in luminal progenitors, and are marked with H3K27Me3 in mature luminal cells . This corresponds with a decrease in gene expression . This is consistent with the model that luminal progenitors are derived from basal cells. However, it is not possible to discern whether luminal commitment of basal cells happens in the adult gland or if it occurs during embryonic development and is maintained in the adult gland. While each subset contains bivalent domains, the highest number occurs in the terminally differentiated mature luminal subset [39, 40]. This challenges the dogma that bivalent domains operate predominantly in stem cells to restrict lineage-specific gene expression . Instead, bivalent domains may be a more general phenomenon in cells of all differentiation states. Key developmental transcription factors within basal and luminal cells contain bivalent promoters; this may keep these genes poised, enabling rapid response to environmental stimuli [28, 40]. For example, in differentiated myoepithelial cells, the epithelial-mesenchymal transition (EMT) transcription factor ZEB1 is held in a bivalent state. Stimulation with transforming growth factor (TGF)β results in the removal of the repressive H3K27Me3 mark while maintaining the active H3K4Me3 mark. This leads to rapid transcription of ZEB1 and de-differentiation of the myoepithelial cells into stem-like cells in culture .
The following section will summarise the current literature surrounding the effector proteins mediating the epigenetic changes to the histone code.
Polycomb group proteins maintain stem/progenitor cells
Polycomb group (PcG) proteins are epigenetic repressors that participate in the establishment and maintenance of cell identity. PcG proteins bind and repress genes that drive differentiation in embryonic and somatic stem cells. Differentiation is accompanied by loss of PcG binding and increased activation of PcG target genes . In mammals there are two PcG chromatin-modifying complexes, Polycomb repressive complexes 1 and 2 (PRC1 and PRC2). These complexes work in a co-ordinated fashion to mediate repression . The PRC2 complex is comprised of SUZ12, EED, RBBP4 or RBBP7, and EZH2 or EZH1 (Table 1). Initial recruitment of PRC2 to the chromatin depends on the DNA methylation state, pre-existing histone modifications, and recruitment by sequence-specific transcription factors . EZH2 catalyses trimethylation of H3K27, which leads PRC1 recruitment. PRC1 is comprised of one member from each of the following paralog groups: PCGF, CBX, PHC, SCML, and RING  (Table 1). PRC1 ubiquitinates H2AK119, which represses transcription by condensing the chromatin and pausing RNA polymerase II .
PcG proteins are important for mammary gland development. As discussed above, H3K27 methylation is associated with gene expression changes that accompany mammary lineage commitment. EZH2, the catalytic subunit of PRC2, has been implicated in co-ordinating histone methylation changes during differentiation. Changes in EZH2 activity are regulated by progesterone during pregnancy and are mirrored by changes in global H3K27Me3 levels, coupling hormonal cues to changes in the epigenetic landscape . Over-expression of EZH2 leads to multi-layered ducts and luminal cell hyperplasia, suggesting that EZH2 drives luminal expansion . Conversely, EZH2 knock-out mice have delayed ductal elongation, and cells derived from these mice have a lower re-populating capacity in vivo and lower clonogenic activity in vitro. Doxycycline-inducible knock-out of EZH2 depletes the luminal progenitor pool, strengthening the role for EZH2 in maintaining luminal progenitors . EZH2 supresses cell cycle inhibitors (e.g. Ink4a and Cdkn1a) and genes involved in epidermal differentiation, suggesting that EZH2 plays a critical role in progenitor cell proliferation and preventing activation of extraneous differentiation programmes . Paradoxically, loss of EZH2 in the mammary epithelium did not alter H3K27Me3 ChIP-Seq profiles. However, because this was conducted on unsorted mammary tissue it is possible that stromal cells with intact EZH2 expression masked the true epigenetic modifications in the mammary epithelial cells. Another possibility is that EZH1, a weaker histone methyltransferase, compensates for EZH2 . EZH2-deficient mice produce milk yet cannot support their pups. This is likely due to a combination of impaired progenitor proliferation leading to reduced alveolar unit density and impaired alveolar differentiation [39, 50].
BMI1, a member of the PCGF paralog group of the PRC1 complex, has also been implicated in mammary gland development. Like EZH2 knock-out mice, BMI1 knock-out mice have severely stunted ductal trees during pubertal development. Mammary epithelial cells from these mice have a 14-fold reduction in re-populating capacity upon serial transplantation, demonstrating the role of BMI1 in MaSC self-renewal . BMI1 is also important for human MaSC self-renewal, where over-expression and knock-down of BMI1 increases and decreases mammosphere-forming capacity, respectively . BMI1 also restricts alveolar differentiation; over-expression of BMI1 blocks alveolar differentiation and loss of BMI1 causes premature alveolar development . The opposing roles of PcG proteins EZH2 and BMI1 in promoting and restricting alveolar differentiation, respectively, is an unexpected result that warrants further investigation.
Lysine demethylases JARID1B, JMJ2B, and KDM6A drive luminal fate commitment
While histone methylation is well studied, enzymes that remove these marks, histone demethylases (known as KDMs), have only been identified more recently. Before their discovery, it was thought that histone methylation was an irreversible modification. It is now evident that histone demethylases play pivotal roles in modifying histones to determine whether a cell maintains multipotency or differentiates. As well as modifying the histone tails directly, certain family members can recruit PcG proteins to further modify chromatin . There are six families of histone demethylase proteins, KDM1–6, each with multiple members that have distinct substrate specificity . Members of the KDM2–6 families contain a Jumonji (or JmjC) domain, which uses a demethylation mechanism distinct from KDM1. So far, three JmjC domain-containing proteins have been identified as important regulators of mammary gland development: JARID1B (KDM5B, PLU-1), JMJ2B (KDM4B), and KDM6A (UTX) (Table 1).
JARID1B removes tri- and di-methylation of H3K4 and is thought to repress transcription [55, 56]. Complete and functional JARID1B knock-out mice have defects in pubertal mammary gland development, including a reduced number of TEBs, less side branching, and impaired ductal elongation [57, 58]. JARID1B mRNA is expressed in both murine basal and luminal lineages (Holliday et al., unpublished data) and seems to be important for rapid epithelial cell proliferation during puberty but not in alveolar development during pregnancy since JARID1B knock-out mice are able to produce milk . The mammary developmental defects in these mice are partially due to a deficiency in systemic oestrogen levels; however, JARID1B also has a mammary cell-intrinsic function . Gene expression analysis on primary mammary epithelial cells and cell lines with perturbed JARID1B expression revealed that JARID1B promotes luminal lineage-specific gene expression and represses basal-specific genes [56, 58, 59]. Key luminal lineage commitment genes (Elf5, Esr1, Pgr, Prlr, and Stat5a) are downregulated in the JARID1B-deficient mammary epithelial cells and breast cancer cell lines [58, 59]. In some but not all cases, increased gene expression was a consequence of JARID1B binding directly to chromatin, overlapping with active H3K4Me2/3 modifications. This is counter-intuitive given its histone demethylase activity on H3K4. It has been suggested that the PHD domain of JARID1B mediates recognition and binding to H3K4Me3 marks, leading to fine-tuning of H3K4 methylation [59, 60]. Zou et al.  observed co-binding of JARID1B and the luminal transcription factor GATA3 at the Foxa1 and Stat5a promoters and loss of GATA3 binding in JARID1B knock-out mammary epithelial cells, arguing that JARID1B and GATA3 can act co-operatively to mediate transcription. JARID1B has also been shown to interact with oestrogen receptor (ER)α in the COS-7 cell line . However, this has not been validated in mammary epithelial cells.
Trimethylation of H3K9 and H3K27 is associated with an inactive chromatin state. JMJD3B and KDM6A demethylate H3K9 and H3K27, respectively, and are therefore transcriptional activators . Knock-out of either of these proteins results in defects in pubertal mammary gland development [61, 62]. JMJD2B interacts with ERα and oestrogen stimulation causes JMJD2B and ERα to localise to chromatin and demethylate H3K9 at ERα target genes . KDM6A knock-out luminal mammary epithelial cells have a gene expression signature more similar to wild-type basal cells than wild-type luminal cells . Like JMJD2B, KDM6A may also be a co-factor for luminal transcription factors since ChIP-Seq analysis on whole mammary glands reveals that KDM6A could bind to promoters and enhancers of ERα, progesterone receptor (PR), and ELF5 target genes . Paradoxically, H3K27Me3 marks were unchanged upon KDM6A knock-out, suggesting that KDM6A has histone demethylase-independent functions, or that KDM6B can compensate for KDM6A loss .
Pygo2 maintains mammary stem cells
The chromatin binding protein Pygo2 is part of the Pygopus family of proteins, which contain a highly conserved PHD domain. Pygo2 is an essential component of the Wnt/β-catenin signalling pathway , fundamental for maintaining stem cell self-renewal in many tissues including the mammary gland . Unlike the PcG and the KDM proteins, Pygo2 does not directly modify chromatin; instead it recognises and binds active H3K4Me3 modifications via its PHD domain. Basal-specific deletion of Pygo2 results in a two-fold reduction of the basal population, decreased the re-populating capacity, and a basal cell gene expression profile which more closely resembles a luminal signature than a MaSC/basal signature [64, 65]. Components of the Notch signalling pathway, a key driver of luminal cell differentiation , are upregulated in Pygo2-deficient basal cells suggesting that Pygo2 normally acts to repress Notch signalling. Pygo2 is required for recruiting β-catenin to the Notch3 locus and maintaining the Notch3 gene in a bivalent state, such that loss of Pygo2 permits Notch-mediated luminal differentiation . Taken together, these studies highlight the role of Pygo2 as a Wnt/β-catenin co-factor that maintains the basal fate by suppression of Notch signalling.
Relevance to breast cancer
There is increasing evidence that different breast cancer subtypes arise from distinct developmental stages along the differentiation hierarchy and retain characteristics of their cell of origin . The epigenetic processes that determine cell fate in normal cells are often hijacked by cancer cells [20, 68]. A comprehensive review of epigenetic perturbation in breast cancer is beyond the scope of this review, although several key examples of developmental epigenetic mechanisms gone awry in breast cancer are discussed below.
The genomes of cancer cells are globally hypomethylated compared with normal cells, resulting in genomic instability . Cancer cells also harbour selective hypermethylation of promoter CpG islands in tumour suppressor genes [69, 70]. Consistent with the role of DNMT1 in maintaining MaSCs, it is also required for cancer stem cell (CSC) maintenance in the MMTV-Neu-Tg mouse mammary tumour model of HER2+ breast cancer [33, 71]. DNMT1 is highly expressed in breast CSCs where it methylates and silences several tumour suppressor genes including Isl1, whose gene product inhibits ERα-mediated transcription . Treatment of mice with the clinical DNMT inhibitor 5-AzaC reduces the CSC pool and significantly improves survival, especially when combined with a histone deacetylase (HDAC) inhibitor .
In addition to DNA methylation, cancer cells also display perturbation of histone modifiers . For example, in non-malignant cells, the core PRC2 component, EZH2, is required for the proliferation of progenitor cells [39, 51]. EZH2 is over-expressed in high-grade basal-like breast cancer where it likely plays a similar role . Similarly, the PRC1 protein BMI1, whose role is to sustain self-renewal of normal mammary stem cells [52, 53], also has proto-oncogenic functions in breast cancer. BMI1 is over-expressed in aggressive basal-like breast cancers where it promotes EMT and self-renewal and is thought to confer drug resistance [73, 74]. JARID1B normally promotes luminal differentiation  and is frequently amplified and over-expressed in luminal breast cancer . In this context, JARID1B is a highly active luminal lineage-specific proto-oncogene that correlates with poor patient outcome . These examples highlight the need to understand the context-dependent regulation of these epigenetic factors in normal development.
Unlike genetic aberrations, epigenetic modifications are reversible, making epigenetic therapy an attractive avenue for treatment. Bromodomain and extra-terminal (BET) domain inhibitors are currently under clinical trial for the treatment of various cancers . These drugs target proteins with a bromodomain, which recognises acetylated histone marks. Bromodomain-containing protein 4 (BRD4) co-operates with FOXO1 to promote the expression of basal genes in MCF10A normal mammary epithelial cells . BET domain inhibitors are generally more efficacious in triple-negative breast cancer (TNBC) cell lines compared with luminal and HER2+ breast cancer cell lines . Interestingly, these drugs seem to push the TNBC cells to a more differentiated luminal state .
Breast tumours are not only made up of epithelial cells, but also contain a heterogeneous microenvironment encompassing multiple different cell types including cancer-associated fibroblasts, dendritic cells, macrophages, and lymphocytes . Epigenetic drugs do not only target the epithelial cells intrinsically but can also target the cross-talk between the epithelial cells and stromal cells. It has recently been shown that epigenetic therapies can sensitise cancer cells to the host immune system and boost the effects of immunotherapies such as check-point inhibitors (reviewed in ). This is particularly relevant to breast cancer, where treatment with check-point inhibitors has been shown to have limited efficacy compared with other cancer types .
Studying the context-dependent role of epigenetic modifiers in development has shaped our understanding of their role in different breast cancer subtypes. This research has uncovered exciting potential novel therapeutic targets and this list is continuing to grow.
In response to microenvironmental cues, multiple layers of epigenetic regulation work in concert to maintain the undifferentiated MaSC state, or to direct the differentiation into specialised myoepithelial, luminal, and alveolar cells (Fig. 1). This review has focused on DNA methylation and the histone code but did not include nucleosome positioning and the three-dimensional organisation of chromatin within the nucleus. Systematic studies of genome-wide chromatin remodelling through mammary gland development and cell fate decisions are lacking and present an exciting area for further investigation.
Perturbation of epigenetic mechanisms can lead to the onset of different subtypes of cancer in a highly lineage-specific manner. A rigorous understanding of the epigenetic processes governing normal mammary development is central to our understanding of breast cancer aetiology and also for employing epigenetic therapies, which are becoming more commonly used in cancer treatment.
There is an increasing appreciation that differentiation occurs as a continuous spectrum, rather than proceeding through stable populations of cells with discrete identities . This presents a challenge when studying bulk cell populations, as do the majority of studies described in this review. Technical advancement in this space has made it possible to perform single-cell ChIP-Seq, RNA-Seq, and the newly developed ATAC-Seq to decipher regions of open and closed chromatin . Future studies will employ these cutting-edge technologies to generate chromatin maps and gene expression profiles to better understand the epigenomic and transcriptomic events that accompany lineage commitment at cellular resolution.
Assay for transposase-accessible chromatin
Bromodomain and extra-terminal
Cancer stem cell
Mammary stem cell
Polycomb repressive complex
Terminal end bud
Triple-negative breast cancer
Visvader JE, Stingl J. Mammary stem cells and the differentiation hierarchy: current status and perspectives. Genes Dev. 2014;28:1143–58.
Macias H, Hinck L. Mammary gland development. Wiley Interdiscip Rev Dev Biol. 2012;1:533–57.
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–7.
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–8.
Wang D, Cai C, Dong X, Yu QC, Zhang XO, Yang L, Zeng YA. Identification of multipotent mammary stem cells by protein C receptor expression. Nature. 2015;517:81–4.
Rios AC, Fu NY, Lindeman GJ, Visvader JE. In situ identification of bipotent stem cells in the mammary gland. Nature. 2014;506:322–7.
Van Keymeulen A, Rocha AS, Ousset M, Beck B, Bouvencourt G, Rock J, Sharma N, Dekoninck S, Blanpain C. Distinct stem cells contribute to mammary gland development and maintenance. Nature. 2011;479:189–93.
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.
Prater MD, Petit V, Alasdair Russell I, Giraddi RR, Shehata M, Menon S, Schulte R, Kalajzic I, Rath N, Olson MF, et al. Mammary stem cells have myoepithelial cell properties. Nat Cell Biol. 2014;16:942–50. 941-947
Wuidart A, Ousset M, Rulands S, Simons BD, Van Keymeulen A, Blanpain C. Quantitative lineage tracing strategies to resolve multipotency in tissue-specific stem cells. Genes Dev. 2016;30:1261–77.
Galliot B, Ghila L. Cell plasticity in homeostasis and regeneration. Mol Reprod Dev. 2010;77:837–55.
Chaffer CL, Marjanovic ND, Lee T, Bell G, Kleer CG, Reinhardt F, D'Alessio AC, Young RA, Weinberg RA. Poised chromatin at the ZEB1 promoter enables breast cancer cell plasticity and enhances tumorigenicity. Cell. 2013;154:61–74.
Holliday R, Pugh JE. DNA modification mechanisms and gene activity during development. Science. 1975;187:226–32.
Smith ZD, Meissner A. DNA methylation: roles in mammalian development. Nat Rev Genet. 2013;14:204–20.
Luger K, Mader AW, Richmond RK, Sargent DF, Richmond TJ. Crystal structure of the nucleosome core particle at 2.8 A resolution. Nature. 1997;389:251–60.
Bannister AJ, Kouzarides T. Regulation of chromatin by histone modifications. Cell Res. 2011;21:381–95.
Skulte KA, Phan L, Clark SJ, Taberlay PC. Chromatin remodeler mutations in human cancers: epigenetic implications. Epigenomics. 2014;6:397–414.
Liang G, Zhang Y. Embryonic stem cell and induced pluripotent stem cell: an epigenetic perspective. Cell Res. 2013;23:49–69.
Sharma S, Gurudutta G. Epigenetic regulation of hematopoietic stem cells. Int J Stem Cells. 2016;9:36–43.
Sharma S, Kelly TK, Jones PA. Epigenetics in cancer. Carcinogenesis. 2010;31:27–36.
Bestor TH. The DNA methyltransferases of mammals. Hum Mol Genet. 2000;9:2395–402.
Laurent L, Wong E, Li G, Huynh T, Tsirigos A, Ong CT, Low HM, Kin Sung KW, Rigoutsos I, Loring J, et al. Dynamic changes in the human methylome during differentiation. Genome Res. 2010;20:320–31.
Huang K, Fan G. DNA methylation in cell differentiation and reprogramming: an emerging systematic view. Regen Med. 2010;5:531–44.
Bogdanovic O. Tet proteins: master regulators of vertebrate body plan formation? Epigenomics. 2017;9:93–6.
Huh SJ, Clement K, Jee D, Merlini A, Choudhury S, Maruyama R, Yoo R, Chytil A, Boyle P, Ran FA, et al. Age- and pregnancy-associated DNA methylation changes in mammary epithelial cells. Stem Cell Rep. 2015;4:297–311.
Dos Santos CO, Dolzhenko E, Hodges E, Smith AD, Hannon GJ. An epigenetic memory of pregnancy in the mouse mammary gland. Cell Rep. 2015;11:1102–9.
Bloushtain-Qimron N, Yao J, Snyder EL, Shipitsin M, Campbell LL, Mani SA, Hu M, Chen H, Ustyansky V, Antosiewicz JE, et al. Cell type-specific DNA methylation patterns in the human breast. Proc Natl Acad Sci U S A. 2008;105:14076–81.
Maruyama R, Choudhury S, Kowalczyk A, Bessarabova M, Beresford-Smith B, Conway T, Kaspi A, Wu Z, Nikolskaya T, Merino VF, et al. Epigenetic regulation of cell type-specific expression patterns in the human mammary epithelium. PLoS Genet. 2011;7:e1001369.
Gascard P, Bilenky M, Sigaroudinia M, Zhao J, Li L, Carles A, Delaney A, Tam A, Kamoh B, Cho S, et al. Epigenetic and transcriptional determinants of the human breast. Nat Commun. 2015;6:6351.
Lee HJ, Hinshelwood RA, Bouras T, Gallego-Ortega D, Valdes-Mora F, Blazek K, Visvader JE, Clark SJ, Ormandy CJ. Lineage specific methylation of the Elf5 promoter in mammary epithelial cells. Stem Cells. 2011;29:1611–9.
Berdasco M, Esteller M. DNA methylation in stem cell renewal and multipotency. Stem Cell Res Ther. 2011;2:42.
Britt K, Ashworth A, Smalley M. Pregnancy and the risk of breast cancer. Endocr Relat Cancer. 2007;14:907–33.
Pathania R, Ramachandran S, Elangovan S, Padia R, Yang P, Cinghu S, Veeranan-Karmegam R, Arjunan P, Gnana-Prakasam JP, Sadanand F, et al. DNMT1 is essential for mammary and cancer stem cell maintenance and tumorigenesis. Nat Commun. 2015;6:6910.
Hebbes TR, Thorne AW, Crane-Robinson C. A direct link between core histone acetylation and transcriptionally active chromatin. EMBO J. 1988;7:1395–402.
Liang G, Lin JC, Wei V, Yoo C, Cheng JC, Nguyen CT, Weisenberger DJ, Egger G, Takai D, Gonzales FA, et al. Distinct localization of histone H3 acetylation and H3-K4 methylation to the transcription start sites in the human genome. Proc Natl Acad Sci U S A. 2004;101:7357–62.
Black JC, Van Rechem C, Whetstine JR. Histone lysine methylation dynamics: establishment, regulation, and biological impact. Mol Cell. 2012;48:491–507.
Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, Fry B, Meissner A, Wernig M, Plath K, et al. A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006;125:315–26.
Calo E, Wysocka J. Modification of enhancer chromatin: what, how, and why? Mol Cell. 2013;49:825–37.
Pal B, Bouras T, Shi W, Vaillant F, Sheridan JM, Fu N, Breslin K, Jiang K, Ritchie ME, Young M, et al. Global changes in the mammary epigenome are induced by hormonal cues and coordinated by Ezh2. Cell Rep. 2013;3:411–26.
Pellacani D, Bilenky M, Kannan N, Heravi-Moussavi A, Knapp DJ, Gakkhar S, Moksa M, Carles A, Moore R, Mungall AJ, et al. Analysis of normal human mammary epigenomes reveals cell-specific active enhancer states and associated transcription factor networks. Cell Rep. 2016;17:2060–74.
Locke WJ, Zotenko E, Stirzaker C, Robinson MD, Hinshelwood RA, Stone A, Reddel RR, Huschtscha LI, Clark SJ. Coordinated epigenetic remodelling of transcriptional networks occurs during early breast carcinogenesis. Clin Epigenetics. 2015;7:52.
Lim E, Wu D, Pal B, Bouras T, Asselin-Labat ML, Vaillant F, Yagita H, Lindeman GJ, Smyth GK, Visvader JE. Transcriptome analyses of mouse and human mammary cell subpopulations reveal multiple conserved genes and pathways. Breast Cancer Res. 2010;12:R21.
Lim E, Vaillant F, Wu D, Forrest NC, Pal B, Hart AH, Asselin-Labat ML, Gyorki DE, Ward T, Partanen A, et al. Aberrant luminal progenitors as the candidate target population for basal tumor development in BRCA1 mutation carriers. Nat Med. 2009;15:907–13.
Fu NY, Rios AC, Pal B, Law CW, Jamieson P, Liu R, Vaillant F, Jackling F, Liu KH, Smyth GK, et al. Identification of quiescent and spatially restricted mammary stem cells that are hormone responsive. Nat Cell Biol. 2017;19:164–76.
Sauvageau M, Sauvageau G. Polycomb group proteins: multi-faceted regulators of somatic stem cells and cancer. Cell Stem Cell. 2010;7:299–313.
van Kruijsbergen I, Hontelez S, Veenstra GJ. Recruiting polycomb to chromatin. Int J Biochem Cell Biol. 2015;67:177–87.
Schwartz YB, Pirrotta V. A new world of Polycombs: unexpected partnerships and emerging functions. Nat Rev Genet. 2013;14:853–64.
Zhou W, Zhu P, Wang J, Pascual G, Ohgi KA, Lozach J, Glass CK, Rosenfeld MG. Histone H2A monoubiquitination represses transcription by inhibiting RNA polymerase II transcriptional elongation. Mol Cell. 2008;29:69–80.
Li X, Gonzalez ME, Toy K, Filzen T, Merajver SD, Kleer CG. Targeted overexpression of EZH2 in the mammary gland disrupts ductal morphogenesis and causes epithelial hyperplasia. Am J Pathol. 2009;175:1246–54.
Michalak EM, Nacerddine K, Pietersen A, Beuger V, Pawlitzky I, Cornelissen-Steijger P, Wientjens E, Tanger E, Seibler J, van Lohuizen M, et al. Polycomb group gene Ezh2 regulates mammary gland morphogenesis and maintains the luminal progenitor pool. Stem Cells. 2013;31:1910–20.
Yoo KH, Oh S, Kang K, Hensel T, Robinson GW, Hennighausen L. Loss of EZH2 results in precocious mammary gland development and activation of STAT5-dependent genes. Nucleic Acids Res. 2015;43:8774–89.
Pietersen AM, Evers B, Prasad AA, Tanger E, Cornelissen-Steijger P, Jonkers J, van Lohuizen M. Bmi1 regulates stem cells and proliferation and differentiation of committed cells in mammary epithelium. Curr Biol. 2008;18:1094–9.
Liu S, Dontu G, Mantle ID, Patel S, Ahn NS, Jackson KW, Suri P, Wicha MS. Hedgehog signaling and Bmi-1 regulate self-renewal of normal and malignant human mammary stem cells. Cancer Res. 2006;66:6063–71.
Nottke A, Colaiacovo MP, Shi Y. Developmental roles of the histone lysine demethylases. Development. 2009;136:879–89.
Yamane K, Tateishi K, Klose RJ, Fang J, Fabrizio LA, Erdjument-Bromage H, Taylor-Papadimitriou J, Tempst P, Zhang Y. PLU-1 is an H3K4 demethylase involved in transcriptional repression and breast cancer cell proliferation. Mol Cell. 2007;25:801–12.
Scibetta AG, Santangelo S, Coleman J, Hall D, Chaplin T, Copier J, Catchpole S, Burchell J, Taylor-Papadimitriou J. Functional analysis of the transcription repressor PLU-1/JARID1B. Mol Cell Biol. 2007;27:7220–35.
Catchpole S, Spencer-Dene B, Hall D, Santangelo S, Rosewell I, Guenatri M, Beatson R, Scibetta AG, Burchell JM, Taylor-Papadimitriou J. PLU-1/JARID1B/KDM5B is required for embryonic survival and contributes to cell proliferation in the mammary gland and in ER+ breast cancer cells. Int J Oncol. 2011;38:1267–77.
Zou MR, Cao J, Liu Z, Huh SJ, Polyak K, Yan Q. Histone demethylase jumonji AT-rich interactive domain 1B (JARID1B) controls mammary gland development by regulating key developmental and lineage specification genes. J Biol Chem. 2014;289:17620–33.
Yamamoto S, Wu Z, Russnes HG, Takagi S, Peluffo G, Vaske C, Zhao X, Moen Vollan HK, Maruyama R, Ekram MB, et al. JARID1B is a luminal lineage-driving oncogene in breast cancer. Cancer Cell. 2014;25:762–77.
Klein BJ, Piao L, Xi Y, Rincon-Arano H, Rothbart SB, Peng D, Wen H, Larson C, Zhang X, Zheng X, et al. The histone-H3K4-specific demethylase KDM5B binds to its substrate and product through distinct PHD fingers. Cell Rep. 2014;6:325–35.
Kawazu M, Saso K, Tong KI, McQuire T, Goto K, Son DO, Wakeham A, Miyagishi M, Mak TW, Okada H. Histone demethylase JMJD2B functions as a co-factor of estrogen receptor in breast cancer proliferation and mammary gland development. PLoS One. 2011;6:e17830.
Yoo KH, Oh S, Kang K, Wang C, Robinson GW, Ge K, Hennighausen L. Histone demethylase KDM6A controls the mammary luminal lineage through enzyme-independent mechanisms. Mol Cell Biol. 2016;36:2108–20.
Belenkaya TY, Han C, Standley HJ, Lin X, Houston DW, Heasman J, Lin X. Pygopus encodes a nuclear protein essential for wingless/Wnt signaling. Development. 2002;129:4089–101.
Gu B, Sun P, Yuan Y, Moraes RC, Li A, Teng A, Agrawal A, Rheaume C, Bilanchone V, Veltmaat JM, et al. Pygo2 expands mammary progenitor cells by facilitating histone H3 K4 methylation. J Cell Biol. 2009;185:811–26.
Gu B, Watanabe K, Sun P, Fallahi M, Dai X. Chromatin effector Pygo2 mediates Wnt-notch crosstalk to suppress luminal/alveolar potential of mammary stem and basal cells. Cell Stem Cell. 2013;13:48–61.
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–41.
Visvader JE. Cells of origin in cancer. Nature. 2011;469:314–22.
Locke WJ, Clark SJ. Epigenome remodelling in breast cancer: insights from an early in vitro model of carcinogenesis. Breast Cancer Res. 2012;14:215.
Robertson KD. DNA methylation and human disease. Nat Rev Genet. 2005;6:597–610.
Jones PA, Issa JP, Baylin S. Targeting the cancer epigenome for therapy. Nat Rev Genet. 2016;17:630–41.
Pathania R, Ramachandran S, Mariappan G, Thakur P, Shi H, Choi JH, Manicassamy S, Kolhe R, Prasad PD, Sharma S, et al. Combined inhibition of DNMT and HDAC blocks the tumorigenicity of cancer stem-like cells and attenuates mammary tumor growth. Cancer Res. 2016;76:3224–35.
Kleer CG, Cao Q, Varambally S, Shen R, Ota I, Tomlins SA, Ghosh D, Sewalt RG, Otte AP, Hayes DF, et al. EZH2 is a marker of aggressive breast cancer and promotes neoplastic transformation of breast epithelial cells. Proc Natl Acad Sci U S A. 2003;100:11606–11.
Guo BH, Feng Y, Zhang R, Xu LH, Li MZ, Kung HF, Song LB, Zeng MS. Bmi-1 promotes invasion and metastasis, and its elevated expression is correlated with an advanced stage of breast cancer. Mol Cancer. 2011;10:10.
Paranjape AN, Balaji SA, Mandal T, Krushik EV, Nagaraj P, Mukherjee G, Rangarajan A. Bmi1 regulates self-renewal and epithelial to mesenchymal transition in breast cancer cells through Nanog. BMC Cancer. 2014;14:785.
Noguchi-Yachide T. BET bromodomain as a target of epigenetic therapy. Chem Pharm Bull. 2016;64:540–7.
Nagarajan S, Bedi U, Budida A, Hamdan FH, Mishra VK, Najafova Z, Xie W, Alawi M, Indenbirken D, Knapp S, et al. BRD4 promotes p63 and GRHL3 expression downstream of FOXO in mammary epithelial cells. Nucleic Acids Res. 2016;45:3130–45
Shu S, Lin CY, He HH, Witwicki RM, Tabassum DP, Roberts JM, Janiszewska M, Huh SJ, Liang Y, Ryan J, et al. Response and resistance to BET bromodomain inhibitors in triple-negative breast cancer. Nature. 2016;529:413–7.
Soysal SD, Tzankov A, Muenst SE. Role of the tumor microenvironment in breast cancer. Pathobiology. 2015;82:142–52.
Dunn J, Rao S. Epigenetics and immunotherapy: the current state of play. Mol Immunol. 2017;87:227–39.
Emens LA. Breast cancer immunotherapy: facts and hopes. Clin Cancer Res. 2018;24:511–20.
Bach K, Pensa S, Grzelak M, Hadfield J, Adams DJ, Marioni JC, Khaled WT. Differentiation dynamics of mammary epithelial cells revealed by single-cell RNA sequencing. Nat Commun. 2017;8:2128.
Clark SJ, Lee HJ, Smallwood SA, Kelsey G, Reik W. Single-cell epigenomics: powerful new methods for understanding gene regulation and cell identity. Genome Biol. 2016;17:72.
The author would like to thank Dr. Kate Patterson for illustrating the template for Fig. 2.
This work was supported by funding from the NHMRC in the form of Australian Postgraduate Awards (to HH and LAB) and a Project Grant (to AS). The McMurtrie Family and the estate of the late R.T. Hall also provided support for this work.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Holliday, H., Baker, L.A., Junankar, S.R. et al. Epigenomics of mammary gland development. Breast Cancer Res 20, 100 (2018). https://doi.org/10.1186/s13058-018-1031-x
- Mammary development
- Breast cancer
- Stem cell