In the Murine and Bovine Maternal Mammary Gland Signal Transducer and Activator of Transcription 3 is Activated in Clusters of Epithelial Cells around the Day of Birth.

Signal transducers and activators of transcription (STAT) proteins regulate mammary development. Here we investigate the expression of phosphorylated STAT3 (pSTAT3) in the mouse and cow around the day of birth. We present localised colocation analysis, applicable to other mammary studies requiring identification of spatially congregated events. We demonstrate that pSTAT3-positive events are multifocally clustered in a non-random and statistically significant fashion. Arginase-1 expressing cells, consistent with macrophages, exhibit distinct clustering within the periparturient mammary gland. These findings represent a new facet of mammary STAT3 biology, and point to the presence of mammary sub-microenvironments.

ABCD4 is associated with mammary gland development in mammals.

BACKGROUND: Mammary gland development is a critical process in mammals, crucial for their reproductive success and offspring nourishment. However, the functional roles of key candidate genes associated with teat number, including ABCD4, VRTN, PROX2, and DLST, in this developmental process remain elusive. To address this gap in knowledge, we conducted an in-depth investigation into the dynamic expression patterns, functional implications, and regulatory networks of these candidate genes during mouse mammary gland development. RESULTS: In this study, the spatial and temporal patterns of key genes were characterized in mammary gland development. Using time-series single-cell data, we uncovered differences in the expression of A bcd4, Vrtn, Prox2, and Dlst in cell population of the mammary gland during embryonic and adult stages, while Vrtn was not detected in any cells. We found that only overexpression and knockdown of Abcd4 could inhibit proliferation and promote apoptosis of HC11 mammary epithelial cells, whereas Prox2 and Dlst had no significant effect on these cells. Using RNA-seq and qPCR, further analysis revealed that Abcd4 can induce widespread changes in the expression levels of genes involved in mammary gland development, such as Igfbp3, Ccl5, Tlr2, and Prlr, which were primarily associated with the MAPK, JAK-STAT, and PI3K-AKT pathways by functional enrichment. CONCLUSIONS: These findings revealed ABCD4 as a candidate gene pivotal for regulating mammary gland development and lactation during pregnancy by influencing PRLR expression.

Uncoupling protein 1-driven Cre (Ucp1-Cre) is expressed in the epithelial cells of mammary glands and various non-adipose tissues.

OBJECTIVE: Uncoupling protein 1 (UCP1), a mitochondrial protein responsible for nonshivering thermogenesis in adipose tissue, serves as a distinct marker for thermogenic brown and beige adipocytes. Ucp1-Cre mice are thus widely used to genetically manipulate these thermogenic adipocytes. However, evidence suggests that UCP1 may also be expressed in non-adipocyte cell types. In this study, we investigated the presence of UCP1 expression in different mouse tissues that have not been previously reported. METHODS: We employed Ucp1-Cre mice crossed with Cre-inducible transgenic reporter Nuclear tagging and Translating Ribosome Affinity Purification (NuTRAP) mice to investigate Ucp1-Cre expression in various tissues of adult female mice and developing embryos. Tamoxifen-inducible Ucp1-CreERT2 mice crossed with NuTRAP mice were used to assess active Ucp1 expression in adult mice. Immunostaining, RNA analysis, and single-cell/nucleus RNA-seq (sc/snRNA-seq) data analysis were performed to determine the expression of endogenous UCP1 and Ucp1-Cre-driven reporter expression. We also investigated the impact of UCP1 deficiency on mammary gland development and function using Ucp1-knockout (KO) mice. RESULTS: Ucp1-Cre expression was observed in the mammary glands within the inguinal white adipose tissue of female Ucp1-Cre; NuTRAP mice. Ucp1-Cre was activated during embryonic development in various tissues, including mammary glands, as well as in the brain, kidneys, eyes, and ears, specifically in epithelial cells in these organs. However, Ucp1-CreERT2 showed no or only partial activation in these tissues of adult mice, indicating the potential for low or transient expression of endogenous Ucp1. While sc/snRNA-seq data suggest potential expression of UCP1 in mammary epithelial cells in adult mice and humans, Ucp1-KO female mice displayed normal mammary gland development and function. CONCLUSIONS: Our findings reveal widespread Ucp1-Cre expression in various non-adipose tissue types, starting during early development. These results highlight the importance of exercising caution when interpreting data and devising experiments involving Ucp1-Cre mice.

Is bovine somatotropin an alternative strategy to overcome the detrimental effects of high-gain diets on prepubertal Holstein × Gyr heifers?

Feeding high-gain diets and an inadequate energy and protein ratio during pre-puberty may lead to impaired growth and mammary gland development of heifers. Thus, frequent application of bovine somatotropin (bST) may prevent future losses in productivity, improve mammary development and animal performance. We aimed to evaluate the effects of bST on digestibility, performance, blood metabolites, mammary gland development, and carcass composition of high-performance prepubertal Holstein × Gyr heifers. Thirty-four Holstein × Gyr heifers with an average initial body weight of 218 ± 49 kg and 14 ± 4 months of age were submitted to an 84-day trial evaluating the effects of no bST or bST injections. Treatments were randomly assigned to each animal within one of the tree blocks. The bST did not influence digestibility or performance parameters. Regarding blood results, IGF1 concentration presented an interaction between treatment and day, where bST heifers had the highest IGF1 concentration. Heifers receiving bST also showed increased ribeye area; however, only an experimental day effect for backfat thickness was observed, with greater accumulation of carcass fat on day 84. Heifers receiving bST had lower pixels/mm² on parenchyma, characteristic of greater parenchymal tissue. Moreover, heifers on bST treatment also had reduced pixels/mm2, characteristic of reduced fat pad tissue. Lastly, bST injections did not influence liver and muscle gene expression, nor most genes evaluated in mammary gland tissue, except for IGFBP3 expression, which was greater for bST heifers. In summary, we confirm the efficacy of bST injections to overcome the detrimental effects of high-gain diets on mammary gland growth and to improve lean carcass gain of prepubertal Holstein × Gyr heifers.

Stabilization of Epithelial ß-Catenin Compromises Mammary Cell Fate Acquisition and Branching Morphogenesis.

The Wnt/ß-catenin pathway plays a critical role in cell fate specification, morphogenesis, and stem cell activation across diverse tissues, including the skin. In mammals, the embryonic surface epithelium gives rise to the epidermis as well as the associated appendages including hair follicles and mammary glands, both of which depend on epithelial Wnt/ß-catenin activity for initiation of their development. Later on, Wnts are thought to enhance mammary gland growth and branching, whereas in hair follicles, they are essential for hair shaft formation. In this study, we report a strong downregulation of epithelial Wnt/ß-catenin activity as the mammary bud progresses to branching. We show that forced activation of epithelial ß-catenin severely compromises embryonic mammary gland branching. However, the phenotype of conditional Lef1-deficient embryos implies that a low level of Wnt/ß-catenin activity is necessary for mammary cell survival. Transcriptomic profiling suggests that sustained high ß-catenin activity leads to maintenance of mammary bud gene signature at the expense of outgrowth/branching gene signature. In addition, it leads to upregulation of epidermal differentiation genes. Strikingly, we find a partial switch to hair follicle fate early on upon stabilization of ß-catenin, suggesting that the level of epithelial Wnt/ß-catenin signaling activity may contribute to the choice between skin appendage identities.

miR-30a-3p Regulates Autophagy in the Involution of Mice Mammary Glands.

The mammary gland undergoes intensive remodeling during the lactation cycle, and the involution process of mammary gland contains extensive epithelial cells involved in the process of autophagy. Our studies of mice mammary glands suggest that miR-30a-3p expression was low during involution compared with its high expression in the mammary glands of lactating mice. Then, we revealed that miR-30a-3p negatively regulated autophagy by autophagy related 12 (Atg12) in mouse mammary gland epithelial cells (MMECs). Restoring ATG12, knocking down autophagy related 5 (Atg5), starvation, and Rapamycin were used to further confirm this conclusion. Overexpression of miR-30a-3p inhibited autophagy and altered mammary structure in the involution of the mammary glands of mice, which was indicative of alteration in mammary remodeling. Taken together, these results elucidated the molecular mechanisms of miR-30a-3p as a key induction mediator of autophagy by targeting Atg12 within the transition period between lactation and involution in mammary glands.

Distinct Requirements for Adaptor Proteins NCK1 and NCK2 in Mammary Gland Development.

The adaptor proteins NCK1 and NCK2 are well-established signalling nodes that regulate diverse biological processes including cell proliferation and actin dynamics in many tissue types. Here we have investigated the distribution and function of Nck1 and Nck2 in the developing mouse mammary gland. Using publicly available single-cell RNA sequencing data, we uncovered distinct expression profiles between the two paralogs. Nck1 showed widespread expression in luminal, basal, stromal and endothelial cells, while Nck2 was restricted to luminal and basal cells, with prominent enrichment in hormone-sensing luminal subtypes. Next, using mice with global knockout of Nck1 or Nck2, we assessed mammary gland development during and after puberty (5, 8 and 12 weeks of age). Mice lacking Nck1 or Nck2 displayed significant defects in ductal outgrowth and branching at 5 weeks compared to controls, and the defects persisted in Nck2 knockout mice at 8 weeks before normalizing at 12 weeks. These defects were accompanied by an increase in epithelial cell proliferation at 5 weeks and a decrease at 8 weeks in both Nck1 and Nck2 knockout mice. We also profiled expression of several key genes associated with mammary gland development at these timepoints and detected temporal changes in transcript levels of hormone receptors as well as effectors of cell proliferation and migration in Nck1 and Nck2 knockout mice, in line with the distinct phenotypes observed at 5 and 8 weeks. Together these studies reveal a requirement for NCK proteins in mammary gland morphogenesis, and suggest that deregulation of Nck expression could drive breast cancer progression and metastasis.

ERBB Receptors and Their Ligands in the Developing Mammary Glands of Different Species: Fifteen Characters in Search of an Author.

The ERBB tyrosine kinase receptors and their ligands belong to a complex family that has diverse biological effects and expression profiles in the developing mammary glands, where its members play an essential role in translating hormone signals into local effects. While our understanding of these processes stems mostly from mouse models, there is the potential for differences in how this family functions in the mammary glands of other species, particularly in light of their unique histomorphological features. Herein we review the postnatal distribution and function of ERBB receptors and their ligands in the mammary glands of rodents and humans, as well as for livestock and companion animals. Our analysis highlights the diverse biology for this family and its members across species, the regulation of their expression, and how their roles and functions might be modulated by varying stromal composition and hormone interactions. Given that ERBB receptors and their ligands have the potential to influence processes ranging from normal mammary development to diseased states such as cancer and/or mastitis, both in human and veterinary medicine, a more complete understanding of their biological functions should help to direct future research and the identification of new therapeutic targets.

Ultrasound characterization of mammary gland development in heifer calves fed at two different levels until weaning.

Ultrasound technologies allow for a non-invasive assessment of mammary gland (MG) development, the differentiation between the tissue types of the MG, and the evaluation of changes in its composition. This study aimed to work out a detailed description of the different stages of MG development that are visually discernible by ultrasonography for providing a template to classify the different structures. With this basis, the qualitative categorization of the developmental stage, as well as further quantitative assessments via pixel densities in the structures of interest, should be facilitated. Ultrasonic images were acquired from all four quarters of 37 German Holstein heifer calves fed either at a high feeding level of milk replacer (MR; 14% solids) at 10 L/day (1.4 kg MR/day; HI, n = 18) or at a restrictive low level of 5.7 L/day (0.8 kg MR/day; RES, n = 19) until linear weaning from week 13 to 14 of life. Ultrasound MG scans were performed first in week 3 of life, fortnightly from week 8-16, and in week 20 of life, in standing position, of each quarter, using a B-mode ultrasound device equipped with a linear probe (18 MHz). The developmental stages of the mammary gland parenchyma (PAR), visible in ultrasound images, obtained over 20 weeks of life, were categorized, described, and drawn by hand. On this basis, a template for classifying the visible categories of mammary PAR development and its surrounding tissue (SURR), and for measuring their pixel brightness was created thus providing an ultrasonographic atlas of the developing bovine MG, describing 11 categories. The ultrasound images were further classified by PAR structure, and pixel brightness was measured in PAR and SURR by using ImageJ Fiji. The difference in pixel brightness between PAR and SURR, the delta (Δ) pixel value was calculated. With increasing age, higher atlas categories of PAR developmental stages were shown. Pixel values, i.e. the brightness of PAR, its SURR, and Δ pixel value changed with age but were neither affected by the feeding group nor by a group × time interaction. With progressing PAR development, its pixel brightness increased from week 10 to 20 of life, i.e., the PAR became more hyperechoic since it spread and grew into its SURR. The atlas can serve as a template for the categorization and qualitative assessment of MG structures and for the quantitative assessment of PAR's development by measuring pixel brightness. With our study, we could show the structural development in PAR as well as in SURR in MG simultaneously in early life and confirm the spreading of PAR into its SURR by ultrasound scanning.

Activation function 1 of progesterone receptor is required for mammary development and regulation of RANKL during pregnancy.

Progesterone receptor (PGR) is a member of the nuclear receptor superfamily of transcription factors. It is critical for mammary stem cells expansion, mammary ductal branching and alveologenesis. The transcriptional activity of PGR is mainly mediated by activation functions AF1 and AF2. Although the discovery of AF1 and AF2 propelled the understanding of the mechanism of gene regulation by nuclear receptors, their physiological roles are still poorly understood. This is largely due to the lack of suitable genetic models. The present study reports gain or loss of AF1 function mutant mouse models in the study of mammary development. The gain of function mutant AF1_QQQ exhibits hyperactivity while the loss of function mutant AF1_FFF shows hypoactivity on mammary development. However, the involvement of AF1 is context dependent. Whereas the AF1_FFF mutation causes significant impairment in mammary development during pregnancy or in response to estrogen and progesterone, it has no effect on mammary development in nulliparous mice. Furthermore, Rankl, but not Wnt4 and Areg is a major target gene of AF1. In conclusion, PGR AF1 is a pivotal ligand-dependent activation domain critical for mammary development during pregnancy and it exerts gene specific effect on PGR regulated genes.

Multidimensional Fluorescence Imaging of Embryonic and Postnatal Mammary Gland Development.

Multidimensional fluorescence imaging represents a powerful approach for studying the dynamic cellular processes underpinning the development, function, and maintenance of the mammary gland. Here, we describe key multidimensional imaging strategies that enable visualization of mammary branching morphogenesis and epithelial cell fate dynamics during postnatal and embryonic mammary gland development. These include 4-dimensional intravital microscopy and ex vivo imaging of embryonic mammary cultures, in addition to methods that facilitate 3-dimensional imaging of the ductal epithelium at single-cell resolution within its native stroma. Collectively, these approaches provide a window into mammary developmental dynamics, and the perturbations underlying tissue dysfunction and disease.

LGL1 binds to Integrin ß1 and inhibits downstream signaling to promote epithelial branching in the mammary gland.

Branching morphogenesis is a fundamental process by which organs in invertebrates and vertebrates form branches to expand their surface areas. The current dogma holds that directional cell migration determines where a new branch forms and thus patterns branching. Here, we asked whether mouse Lgl1, a homolog of the Drosophila tumor suppressor Lgl, regulates epithelial polarity in the mammary gland. Surprisingly, mammary glands lacking Lgl1 have normal epithelial polarity, but they form fewer branches. Moreover, we find that Lgl1 null epithelium is unable to directionally migrate, suggesting that migration is not essential for mammary epithelial branching as expected. We show that LGL1 binds to Integrin ß1 and inhibits its downstream signaling, and Integrin ß1 overexpression blocks epithelial migration, thus recapitulating the Lgl1 null phenotype. Altogether, we demonstrate that Lgl1 modulation of Integrin ß1 signaling is essential for directional migration and that epithelial branching in invertebrates and the mammary gland is fundamentally distinct.

Evaluation of early post-natal pig mammary gland development and human breast cancer gene expression.

Breast cancer is the second leading cause of death in women after lung cancer, and only 5% of patients with metastatic breast cancer survive beyond ten years of diagnosis. Considering the heterogeneous subclasses of breast cancer, current cancer models have shortfalls due to copy number variants, and genetic differences of humans and immunocompromised animal models. Preclinical studies indicate stem cell activity in early post-natal mammary development may be reactivated in the human adult as a trigger to initiate cell proliferation leading to breast cancer. The goal of the work reported herein was to compare genetic expression of early development, post-natal pig mammary glands to the literature reported genes implicated in different subclasses of human breast cancer. Differentially expressed genes associated with breast cancer and present in early developing pig samples include NUCB2, ANGPTL4 and ACE. Histological staining confirmed E-cadherin, Vimentin, N-cadherin, and Claudin-1, which are all implicated in malignant cancer. Due to the homology of gene expression patterns in the developing pig mammary gland and reported genes in human breast cancer profiles, this research is worthy of further study to address a potential model using mammary development cues to unravel breast cancer biology.

Gestational and lactational xenoestrogen exposure disrupts morphology and inflammatory aspects in mammary gland of gerbil mothers during involution.

In the mammary gland (MG), the developmental window for gestational/lactational differentiation and growth is highly vulnerable to hormonal disruption. Here we describe that the MG involution process in female gerbil mothers is delayed by bisphenol A (BPA) exposure during gestation and lactation. The process is directly influenced by changes in expression of extracellular matrix proteases MMP-2, MMP-9, and FAP, and the incidence of collagen and elastin is reduced after 7 and 14 days of weaning. A pro-inflammatory environment in the late involution process was confirmed by higher expression of TNF-α, COX-2 and phospho-STAT3 n the MG stroma, allied to increases in the incidence of macrophages and mast cells. These aspects impacted the proliferative pattern of epithelial cells, which decreased on the 14th post-weaning day. These data confirm that the milk production window of susceptibility is vulnerable to the impact of BPA, which promotes a suggestive pro-tumoral microenvironment during mammary involution.

Unique Transcriptomic Changes Underlie Hormonal Interactions During Mammary Histomorphogenesis in Female Pigs.

Successful lactation and the risk for developing breast cancer depend on growth and differentiation of the mammary gland (MG) epithelium that is regulated by ovarian steroids (17ß-estradiol [E] and progesterone [P]) and pituitary-derived prolactin (PRL). Given that the MG of pigs share histomorphogenic features present in the normal human breast, we sought to define the transcriptional responses within the MG of pigs following exposure to all combinations of these hormones. Hormone-ablated female pigs were administered combinations of E, medroxyprogesterone 17-acetate (source of P), and either haloperidol (to induce PRL) or 2-bromo-α-ergocryptine. We subsequently monitored phenotypic changes in the MG including mitosis, receptors for E and P (ESR1 and PGR), level of phosphorylated STAT5 (pSTAT5), and the frequency of terminal ductal lobular unit (TDLU) subtypes; these changes were then associated with all transcriptomic changes. Estrogen altered the expression of approximately 20% of all genes that were mostly associated with mitosis, whereas PRL stimulated elements of fatty acid metabolism and an inflammatory response. Several outcomes, including increased pSTAT5, highlighted the ability of E to enhance PRL action. Regression of transcriptomic changes against several MG phenotypes revealed 1669 genes correlated with proliferation, among which 29 were E inducible. Additional gene expression signatures were associated with TDLU formation and the frequency of ESR1 or PGR. These data provide a link between the hormone-regulated genome and phenome of the MG in a species having a complex histoarchitecture like that in the human breast, and highlight an underexplored synergy between the actions of E and PRL during MG development.

In Utero Exposure to trans-10, cis-12 Conjugated Linoleic Acid Modifies Postnatal Development of the Mammary Gland and its Hormone Responsiveness.

We previously showed that dietary trans-10, cis-12 conjugated linoleic acid (10,12 CLA) stimulates estrogen-independent mammary growth in young ovariectomized mice. Here we investigated the effects of in utero or postnatal exposure to cis-9, trans-11 (9,11 CLA) and 10,12 CLA on postnatal development of the mammary gland and its responsiveness to ovarian steroids. In the first experiment we fed dams different CLA prior to and during gestation, then cross fostered female pups onto control fed dams prior to assessing the histomorphology of their mammary glands. Pregnant dams in the second experiment were similarly exposed to CLA, after which their female pups were ovariectomized then treated with 17ß-estradiol (E), progesterone (P) or E + P for 5 days. In a third experiment, mature female mice were fed different CLA for 28 days prior to ovariectomy, then treated with E, P or E + P. Our data indicate that 10,12 CLA modifies the responsiveness of the mammary glands to E or E + P when exposure occurs either in utero, or postnatally. These findings underline the sensitivity of the mammary glands to dietary fatty acids and reinforce the potential for maternal nutrition to impact postnatal development of the mammary glands and their risk for developing cancer.

A matrisome RNA signature from early-pregnancy mouse mammary fibroblasts predicts distant metastasis-free breast cancer survival in humans.

BACKGROUND: During pregnancy, the mouse mammary ductal epithelium branches and grows into the surrounding stroma, requiring extensive extracellular matrix (ECM) and tissue remodelling. It therefore shows parallels to cancer invasion. We hypothesised that similar molecular mechanisms may be utilised in both processes, and that assessment of the stromal changes during pregnancy-associated branching may depict the stromal involvement during human breast cancer progression. METHODS: Immunohistochemistry (IHC) was employed to assess the alterations within the mouse mammary gland extracellular matrix during early pregnancy when lateral branching of the primary ductal epithelium is initiated. Primary mouse mammary fibroblasts from three-day pregnant and age-matched non-pregnant control mice, respectively, were 3D co-cultured with mammary epithelial cells to assess differences in their abilities to induce branching morphogenesis in vitro. Transcriptome analysis was performed to identify the underlying molecular changes. A signature of the human orthologues of the differentially expressed matrisome RNAs was analysed by Kaplan-Meier and multi-variate analysis in two large breast cancer RNA datasets (Gene expression-based Outcome for Breast cancer Online (GOBO) und Kaplan-Meier Plotter), respectively, to test for similarities in expression between early-pregnancy mouse mammary gland development and breast cancer progression. RESULTS: The ECM surrounding the primary ductal network showed significant differences in collagen and basement membrane protein distribution early during pregnancy. Pregnancy-associated fibroblasts (PAFs) significantly enhanced branching initiation compared to age-matched control fibroblast. A combined signature of 64 differentially expressed RNAs, encoding matrisome proteins, was a strong prognostic indicator of distant metastasis-free survival (DMFS) independent of other clinical parameters. The prognostic power could be significantly strengthened by using only a subset of 18 RNAs (LogRank P ≤ 1.00e-13; Hazard ratio (HR) = 2.42 (1.8-3.26); p = 5.61e-09). The prognostic power was confirmed in a second breast cancer dataset, as well as in datasets from ovarian and lung cancer patients. CONCLUSIONS: Our results describe for the first time the early stromal changes that accompany pregnancy-associated branching morphogenesis in mice, specify the early pregnancy-associated molecular alterations in mouse mammary fibroblasts, and identify a matrisome signature as a strong prognostic indicator of human breast cancer progression, with particular strength in oestrogen receptor (ER)-negative breast cancers.

Developing ovine mammary terminal duct lobular units have a dynamic mucosal and stromal immune microenvironment.

The human breast and ovine mammary gland undergo striking levels of postnatal development, leading to formation of terminal duct lobular units (TDLUs). Here we interrogate aspects of sheep TDLU growth as a model of breast development and to increase understanding of ovine mammogenesis. The distributions of epithelial nuclear Ki67 positivity differ significantly between younger and older lambs. Ki67 expression is polarised to the leading edge of the developing TDLUs. Intraepithelial ductal macrophages exhibit periodicity and considerably increased density in lambs approaching puberty. Stromal macrophages are more abundant centrally than peripherally. Intraepithelial T lymphocytes are more numerous in older lambs. Stromal hotspots of Ki67 expression colocalize with immune cell aggregates that exhibit distinct organisation consistent with tertiary lymphoid structures. The lamb mammary gland thus exhibits a dynamic mucosal and stromal immune microenvironment and constitutes a valuable model system that provides new insights into postnatal breast development.

Microenvironmental control of cell fate decisions in mammary gland development and cancer.

Cell fate decisions are critical for adequate tissue development, maintenance and regeneration. In the mammary gland, epithelial cell fates are tightly controlled by the microenvironment. Here, we review how cell fate decisions are regulated by components of the microenvironment during mammary gland development and how pathological changes in the microenvironment can alter cell fates, leading to malignancy. Specifically, we describe the current understanding of how mammary cell fate is controlled and directed by three elements: the extracellular matrix, the immune microenvironment, and hormones-and how these elements can converge to create microenvironments that promote a fourth element: DNA damage.

Single cell transcriptome atlas of mouse mammary epithelial cells across development.

BACKGROUND: Heterogeneity within the mouse mammary epithelium and potential lineage relationships have been recently explored by single-cell RNA profiling. To further understand how cellular diversity changes during mammary ontogeny, we profiled single cells from nine different developmental stages spanning late embryogenesis, early postnatal, prepuberty, adult, mid-pregnancy, late-pregnancy, and post-involution, as well as the transcriptomes of micro-dissected terminal end buds (TEBs) and subtending ducts during puberty. METHODS: The single cell transcriptomes of 132,599 mammary epithelial cells from 9 different developmental stages were determined on the 10x Genomics Chromium platform, and integrative analyses were performed to compare specific time points. RESULTS: The mammary rudiment at E18.5 closely aligned with the basal lineage, while prepubertal epithelial cells exhibited lineage segregation but to a less differentiated state than their adult counterparts. Comparison of micro-dissected TEBs versus ducts showed that luminal cells within TEBs harbored intermediate expression profiles. Ductal basal cells exhibited increased chromatin accessibility of luminal genes compared to their TEB counterparts suggesting that lineage-specific chromatin is established within the subtending ducts during puberty. An integrative analysis of five stages spanning the pregnancy cycle revealed distinct stage-specific profiles and the presence of cycling basal, mixed-lineage, and 'late' alveolar intermediates in pregnancy. Moreover, a number of intermediates were uncovered along the basal-luminal progenitor cell axis, suggesting a continuum of alveolar-restricted progenitor states. CONCLUSIONS: This extended single cell transcriptome atlas of mouse mammary epithelial cells provides the most complete coverage for mammary epithelial cells during morphogenesis to date. Together with chromatin accessibility analysis of TEB structures, it represents a valuable framework for understanding developmental decisions within the mouse mammary gland.