Single-Cell Chromatin Analysis of Mammary Gland Development Reveals Cell-State Transcriptional Regulators and Lineage Relationships

Comprehensive single-cell aging atlas of healthy mammary tissues reveals shared epigenomic and transcriptomic signatures of aging and cancer

Aging is the greatest risk factor for breast cancer; however, how age-related cellular and molecular events impact cancer initiation is unknown. In this study, we investigated how aging rewires transcriptomic and epigenomic programs of mouse mammary glands at single-cell resolution, yielding a comprehensive resource for aging and cancer biology. Aged epithelial cells exhibit epigenetic and transcriptional changes in metabolic, pro-inflammatory and cancer-associated genes. Aged stromal cells downregulate fibroblast marker genes and upregulate markers of senescence and cancer-associated fibroblasts. Among immune cells, distinct T cell subsets (Gzmk+, memory CD4+, γδ) and M2-like macrophages expand with age. Spatial transcriptomics reveals co-localization of aged immune and epithelial cells in situ. Lastly, we found transcriptional signatures of aging mammary cells in human breast tumors, suggesting possible links between aging and cancer. Together, these data uncover that epithelial, immune and stromal cells shift in proportions and cell identity, potentially impacting cell plasticity, aged microenvironment and neoplasia risk.

Single-Cell Analysis in the Mouse and Human Mammary Gland

The mammary gland is a complex organ, host to a rich array of different cell types. As the only organ to complete its development in adulthood, it delicately balances both cell intrinsic and external signalling from hormones, growth factors and other stimulants. The gland can undergo vast proliferation, restructuring and functional maturation during pregnancy and undo these gross changes to a nearly identical resting state during involution. The adaptive nature of the mammary gland underpins its function but also increases its susceptibility to cancer. While already characterised at a macro scale, understanding the complexities of mammary gland morphogenesis, development and tumorigenesis requires interrogation of cellular and molecular mechanisms. As outlined below, single-cell analysis is a key approach for this, allowing us to unbiasedly explore and characterise the functions and properties of individual cells from the genome to the proteome. Here, we introduce key single-cell analysis methods and give brief introductions to their respective workflows. We then discuss the structure, cell types and development of the mammary gland from birth, puberty and through pregnancy, as well as cancer formation. Additionally, we highlight the benefits and caveats of implementing single-cell methodologies and mouse models for studying critical time points of human development and disease. Finally, we highlight some limitations and future directions of single-cell techniques. This chapter provides a starting point for users hoping to further their understanding of mammary gland development and its link to cancer as explained by single-cell analysis studies.

Collagen signaling and matrix stiffness regulate multipotency in glandular epithelial stem cells in mice

Glandular epithelia, including mammary gland (MG) and prostate, are composed of luminal and basal cells. During embryonic development, glandular epithelia arise from multipotent stem cells (SCs) that are replaced after birth by unipotent basal and unipotent luminal SCs. Different conditions, such as basal cell transplantation, luminal cell ablation, and oncogene expression can reinduce adult basal SC (BaSCs) multipotency in different glandular epithelia. The mechanisms regulating the reactivation of multipotency are incompletely understood. Here, we have found that Collagen I expression is commonly upregulated in BaSCs across the different multipotent conditions. Increasing collagen concentration or stiffness of the extracellular matrix (ECM) promotes BaSC multipotency in MG and prostate organoids. Single cell RNA-seq of MG organoids in stiff conditions have uncovered the importance of β1 integrin/FAK/AP-1 axis in the regulation of BaSC multipotency. Altogether our study uncovers the key role of Collagen signaling and ECM stiffness in the regulation of multipotency in glandular epithelia.

Brca1 haploinsufficiency promotes early tumor onset and epigenetic alterations in a mouse model of hereditary breast cancer

Germline BRCA1 mutation carriers face a high breast cancer risk; however, the underlying mechanisms for this risk are not completely understood. Using a new genetically engineered mouse model of germline Brca1 heterozygosity, we demonstrate that early tumor onset in a Brca1 heterozygous background cannot be fully explained by the conventional 'two-hit' hypothesis, suggesting the existence of inherent tumor-promoting alterations in the Brca1 heterozygous state. Single-cell RNA sequencing and assay for transposase-accessible chromatin with sequencing analyses uncover a unique set of differentially accessible chromatin regions in ostensibly normal Brca1 heterozygous mammary epithelial cells, distinct from wild-type cells and partially mimicking the chromatin and RNA-level changes in tumor cells. Transcription factor analyses identify loss of ELF5 and gain of AP-1 sites in these epigenetically primed regions; in vivo experiments further implicate AP-1 and Wnt10a as strong promoters of Brca1-related breast cancer. These findings reveal a previously unappreciated epigenetic effect of Brca1 haploinsufficiency in accelerating tumorigenesis, advancing our mechanistic understanding and informing potential therapeutic strategies.

The molecular chronology of mammary epithelial cell fate switching

The adult mammary gland is maintained by lineage-restricted progenitor cells through pregnancy, lactation, involution, and menopause. Injury resolution, transplantation-associated mammary gland reconstitution, and tumorigenesis are unique exceptions, wherein mammary basal cells gain the ability to reprogram to a luminal state. Here, we leverage newly developed cell-identity reporter mouse strains, and time-resolved single-cell epigenetic and transcriptomic analyses to decipher the molecular programs underlying basal-to-luminal fate switching in vivo. We demonstrate that basal cells rapidly reprogram toward plastic cycling intermediates that appear to hijack molecular programs we find in bipotent fetal mammary stem cells and puberty-associatiated cap cells. Loss of basal-cell specifiers early in dedifferentiation coincides with activation of Notch and BMP, among others. Pharmacologic blockade of each pathway disrupts basal-to-luminal transdifferentiation. Our studies provide a comprehensive map and resource for understanding the coordinated molecular changes enabling terminally differentiated epithelial cells to transition between cell lineages and highlights the stunning rapidity by which epigenetic reprogramming can occur in response to disruption of tissue structure.

Differential chromatin accessibility and transcriptional dynamics define breast cancer subtypes and their lineages

Breast cancer (BC) is defined by distinct molecular subtypes with different cells of origin. The transcriptional networks that characterize the subtype-specific tumor-normal lineages are not established. In this work, we applied bulk, single-cell and single-nucleus multi-omic techniques as well as spatial transcriptomics and multiplex imaging on 61 samples from 37 patients with BC to show characteristic links in gene expression and chromatin accessibility between BC subtypes and their putative cells of origin. Regulatory network analysis of transcription factors underscored the importance of BHLHE40 in luminal BC and luminal mature cells and KLF5 in basal-like tumors and luminal progenitor cells. Furthermore, we identify key genes defining the basal-like (SOX6 and KCNQ3) and luminal A/B (FAM155A and LRP1B) lineages. Exhausted CTLA4-expressing CD8+ T cells were enriched in basal-like BC, suggesting an altered means of immune dysfunction. These findings demonstrate analysis of paired transcription and chromatin accessibility at the single-cell level is a powerful tool for investigating cancer lineage and highlight transcriptional networks that define basal and luminal BC lineages.

Deciphering breast cancer dynamics: insights from single-cell and spatial profiling in the multi-omics era

As one of the most common tumors in women, the pathogenesis and tumor heterogeneity of breast cancer have long been the focal point of research, with the emergence of tumor metastasis and drug resistance posing persistent clinical challenges. The emergence of single-cell sequencing (SCS) technology has introduced novel approaches for gaining comprehensive insights into the biological behavior of malignant tumors. SCS is a high-throughput technology that has rapidly developed in the past decade, providing high-throughput molecular insights at the individual cell level. Furthermore, the advent of multitemporal point sampling and spatial omics also greatly enhances our understanding of cellular dynamics at both temporal and spatial levels. The paper provides a comprehensive overview of the historical development of SCS, and highlights the most recent advancements in utilizing SCS and spatial omics for breast cancer research. The findings from these studies will serve as valuable references for future advancements in basic research, clinical diagnosis, and treatment of breast cancer.

FGD5 in basal cells induces CXCL14 secretion that initiates a feedback loop to promote murine mammary epithelial growth and differentiation

The interactions of environmental compartments with epithelial cells are essential for mammary gland development and homeostasis. Currently, the direct crosstalk between the endothelial niche and mammary epithelial cells remains poorly understood. Here, we show that faciogenital dysplasia 5 (FGD5) is enriched in mammary basal cells (BCs) and mediates critical interactions between basal and endothelial cells (ECs) in the mammary gland. Conditional deletion of Fgd5 reduced, whereas conditional knockin of Fgd5 increased, the engraftment and expansion of BCs, regulating ductal morphogenesis in the mammary gland. Mechanistically, murine mammary BC-expressed FGD5 inhibited the transcriptional activity of activating transcription factor 3 (ATF3), leading to subsequent transcriptional activation and secretion of CXCL14. Furthermore, activation of CXCL14/CXCR4/ERK signaling in primary murine mammary stromal ECs enhanced the expression of HIF-1α-regulated hedgehog ligands, which initiated a positive feedback loop to promote the function of BCs. Collectively, these findings identify functionally important interactions between BCs and the endothelial niche that occur through the FGD5/CXCL14/hedgehog axis.

Ehf controls mammary alveolar lineage differentiation and is a putative suppressor of breast tumorigenesis

The transcription factor EHF is highly expressed in the lactating mammary gland, but its role in mammary development and tumorigenesis is not fully understood. Utilizing a mouse model of Ehf deletion, herein, we demonstrate that loss of Ehf impairs mammary lobuloalveolar differentiation at late pregnancy, indicated by significantly reduced levels of milk genes and milk lipids, fewer differentiated alveolar cells, and an accumulation of alveolar progenitor cells. Further, deletion of Ehf increased proliferative capacity and attenuated prolactin-induced alveolar differentiation in mammary organoids. Ehf deletion also increased tumor incidence in the MMTV-PyMT mammary tumor model and increased the proliferative capacity of mammary tumor organoids, while low EHF expression was associated with higher tumor grade and poorer outcome in luminal A and basal human breast cancers. Collectively, these findings establish EHF as a non-redundant regulator of mammary alveolar differentiation and a putative suppressor of mammary tumorigenesis.

Epigenetic activation of SOX11 is associated with recurrence and progression of ductal carcinoma in situ to invasive breast cancer

BACKGROUND

Risk of recurrence and progression of ductal carcinoma in situ (DCIS) to invasive cancer remains uncertain, emphasizing the need for developing predictive biomarkers of aggressive DCIS.

METHODS

Human cell lines and mouse models of disease progression were analyzed for candidate risk predictive biomarkers identified and validated in two independent DCIS cohorts.

RESULTS

RNA profiling of normal mammary and DCIS tissues (n = 48) revealed that elevated SOX11 expression correlates with MKI67, EZH2, and DCIS recurrence score. The 21T human cell line model of DCIS progression to invasive cancer and two mouse models developing mammary intraepithelial neoplasia confirmed the findings. AKT activation correlated with chromatin accessibility and EZH2 enrichment upregulating SOX11 expression. AKT and HER2 inhibitors decreased SOX11 expression along with diminished mammosphere formation. SOX11 was upregulated in HER2+ and basal-like subtypes (P < 0.001). Longitudinal DCIS cohort (n = 194) revealed shorter recurrence-free survival in SOX11+ than SOX11- patients (P = 0.0056 in all DCIS; P < 0.0001 in HER2+ subtype) associated with increased risk of ipsilateral breast event/IBE (HR = 1.9, 95%CI = 1.2-2.9; P = 0.003).

DISCUSSION

Epigenetic activation of SOX11 drives recurrence of DCIS and progression to invasive cancer, suggesting SOX11 as a predictive biomarker of IBE.

Spatially distinct epithelial and mesenchymal cell subsets along progressive lineage restriction in the branching embryonic mammary gland

How cells coordinate morphogenetic cues and fate specification during development remains a fundamental question in organogenesis. The mammary gland arises from multipotent stem cells (MaSCs), which are progressively replaced by unipotent progenitors by birth. However, the lack of specific markers for early fate specification has prevented the delineation of the features and spatial localization of MaSC-derived lineage-committed progenitors. Here, using single-cell RNA sequencing from E13.5 to birth, we produced an atlas of matched mouse mammary epithelium and mesenchyme and reconstructed the differentiation trajectories of MaSCs toward basal and luminal fate. We show that murine MaSCs exhibit lineage commitment just prior to the first sprouting events of mammary branching morphogenesis at E15.5. We identify early molecular markers for committed and multipotent MaSCs and define their spatial distribution within the developing tissue. Furthermore, we show that the mammary embryonic mesenchyme is composed of two spatially restricted cell populations, and that dermal mesenchyme-produced FGF10 is essential for embryonic mammary branching morphogenesis. Altogether, our data elucidate the spatiotemporal signals underlying lineage specification of multipotent MaSCs, and uncover the signals from mesenchymal cells that guide mammary branching morphogenesis.

Scalable and unbiased sequence-informed embedding of single-cell ATAC-seq data with CellSpace

Standard scATAC sequencing (scATAC-seq) analysis pipelines represent cells as sparse numeric vectors relative to an atlas of peaks or genomic tiles and consequently ignore genomic sequence information at accessible loci. Here we present CellSpace, an efficient and scalable sequence-informed embedding algorithm for scATAC-seq that learns a mapping of DNA k-mers and cells to the same space, to address this limitation. We show that CellSpace captures meaningful latent structure in scATAC-seq datasets, including cell subpopulations and developmental hierarchies, and can score transcription factor activities in single cells based on proximity to binding motifs embedded in the same space. Importantly, CellSpace implicitly mitigates batch effects arising from multiple samples, donors or assays, even when individual datasets are processed relative to different peak atlases. Thus, CellSpace provides a powerful tool for integrating and interpreting large-scale scATAC-seq compendia.

Single-cell analyses reveal evolution mimicry during the specification of breast cancer subtype

Background: The stem or progenitor antecedents confer developmental plasticity and unique cell identities to cancer cells via genetic and epigenetic programs. A comprehensive characterization and mapping of the cell-of-origin of breast cancer using novel technologies to unveil novel subtype-specific therapeutic targets is still absent. Methods: We integrated 195,144 high-quality cells from normal breast tissues and 406,501 high-quality cells from primary breast cancer samples to create a large-scale single-cell atlas of human normal and cancerous breasts. Potential heterogeneous origin of malignant cells was explored by contrasting cancer cells against reference normal epithelial cells. Multi-omics analyses and both in vitro and in vivo experiments were performed to screen and validate potential subtype-specific treatment targets. Novel biomarkers of identified immune and stromal cell subpopulations were validated by immunohistochemistry in our cohort. Results: Tumor stratification based on cancer cell-of-origin patterns correlated with clinical outcomes, genomic aberrations and diverse microenvironment constitutions. We found that the luminal progenitor (LP) subtype was robustly associated with poor prognosis, genomic instability and dysfunctional immune microenvironment. However, the LP subtype patients were sensitive to neoadjuvant chemotherapy (NAC), PARP inhibitors (PARPi) and immunotherapy. The LP subtype-specific target PLK1 was investigated by both in vitro and in vivo experiments. Besides, large-scale single-cell profiling of breast cancer inspired us to identify a range of clinically relevant immune and stromal cell subpopulations, including subsets of innate lymphoid cells (ILCs), macrophages and endothelial cells. Conclusion: The present single-cell study revealed the cellular repertoire and cell-of-origin patterns of breast cancer. Combining single-cell and bulk transcriptome data, we elucidated the evolution mimicry from normal to malignant subtypes and expounded the LP subtype with vital clinical implications. Novel immune and stromal cell subpopulations of breast cancer identified in our study could be potential therapeutic targets. Taken together, Our findings lay the foundation for the precise prognostic and therapeutic stratification of breast cancer.

Research progress of SWI/SNF complex in breast cancer

In the past decade, numerous epigenetic mechanisms have been discovered to be associated with cancer. The mammalian SWI/SNF complex is an ATP-dependent chromatin remodeling complex whose mutations are associated with various malignancies including breast cancer. As the SWI/SNF complex has become one of the most commonly mutated complexes in cancer, targeting epigenetic mutations acquired during breast cancer progress is a potential means of improving clinical efficacy in treatment strategies. This article reviews the composition of the SWI/SNF complex, its main roles and research progress in breast cancer, and links these findings to the latest discoveries in cancer epigenomics to discuss the potential mechanisms and therapeutic potential of SWI/SNF in breast cancer.

Cancer as a Disease of Development Gone Awry

In the 160 years since Rudolf Virchow first postulated that neoplasia arises by the same law that regulates embryonic development, scientists have come to recognize the striking overlap between the molecular and cellular programs used by cancers and embryos. Advances in cancer biology and molecular techniques have further highlighted the similarities between carcinogenesis and embryogenesis, where cellular growth, differentiation, motility, and intercellular cross talk are mediated by common drivers and regulatory networks. This review highlights the many connections linking cancer biology and developmental biology to provide a deeper understanding of how a tissue's developmental history may both enable and constrain cancer cell evolution.

Targeting ABCA12-controlled ceramide homeostasis inhibits breast cancer stem cell function and chemoresistance

Cancer stem cells (CSCs) drive tumor growth, metastasis, and chemoresistance. While emerging evidence suggests that CSCs have a unique dependency on lipid metabolism, the functions and regulation of distinct lipid species in CSCs remain poorly understood. Here, we developed a stem cell factor SOX9-based reporter for isolating CSCs in primary tumors and metastases of spontaneous mammary tumor models. Transcriptomic analyses uncover that SOX9high CSCs up-regulate the ABCA12 lipid transporter. ABCA12 down-regulation impairs cancer stemness and chemoresistance. Lipidomic analyses reveal that ABCA12 maintains cancer stemness and chemoresistance by reducing intracellular ceramide abundance, identifying a CSC-associated function of ABCA subfamily transporter. Ceramide suppresses cancer stemness by inhibiting the YAP-SOX9 signaling pathway in CSCs. Increasing ceramide levels in tumors enhances their sensitivity to chemotherapy and prevents the enrichment of SOX9high CSCs. In addition, SOX9high and ABCA12high cancer cells contribute to chemoresistance in human patient-derived xenografts. These findings identify a CSC-suppressing lipid metabolism pathway that can be exploited to inhibit CSCs and overcome chemoresistance.

Three-dimensional genome architecture coordinates key regulators of lineage specification in mammary epithelial cells

Although lineage-specific genes have been identified in the mammary gland, little is known about the contribution of the 3D genome organization to gene regulation in the epithelium. Here, we describe the chromatin landscape of the three major epithelial subsets through integration of long- and short-range chromatin interactions, accessibility, histone modifications, and gene expression. While basal genes display exquisite lineage specificity via distal enhancers, luminal-specific genes show widespread promoter priming in basal cells. Cell specificity in luminal progenitors is largely mediated through extensive chromatin interactions with super-enhancers in gene-body regions in addition to interactions with polycomb silencer elements. Moreover, lineage-specific transcription factors appear to be controlled through cell-specific chromatin interactivity. Finally, chromatin accessibility rather than interactivity emerged as a defining feature of the activation of quiescent basal stem cells. This work provides a comprehensive resource for understanding the role of higher-order chromatin interactions in cell-fate specification and differentiation in the adult mouse mammary gland.

Differential chromatin accessibility and transcriptional dynamics define breast cancer subtypes and their lineages

Breast cancer is a heterogeneous disease, and treatment is guided by biomarker profiles representing distinct molecular subtypes. Breast cancer arises from the breast ductal epithelium, and experimental data suggests breast cancer subtypes have different cells of origin within that lineage. The precise cells of origin for each subtype and the transcriptional networks that characterize these tumor-normal lineages are not established. In this work, we applied bulk, single-cell (sc), and single-nucleus (sn) multi-omic techniques as well as spatial transcriptomics and multiplex imaging on 61 samples from 37 breast cancer patients to show characteristic links in gene expression and chromatin accessibility between breast cancer subtypes and their putative cells of origin. We applied the PAM50 subtyping algorithm in tandem with bulk RNA-seq and snRNA-seq to reliably subtype even low-purity tumor samples and confirm promoter accessibility using snATAC. Trajectory analysis of chromatin accessibility and differentially accessible motifs clearly connected progenitor populations with breast cancer subtypes supporting the cell of origin for basal-like and luminal A and B tumors. Regulatory network analysis of transcription factors underscored the importance of BHLHE40 in luminal breast cancer and luminal mature cells, and KLF5 in basal-like tumors and luminal progenitor cells. Furthermore, we identify key genes defining the basal-like ( PRKCA , SOX6 , RGS6 , KCNQ3 ) and luminal A/B ( FAM155A , LRP1B ) lineages, with expression in both precursor and cancer cells and further upregulation in tumors. Exhausted CTLA4-expressing CD8+ T cells were enriched in basal-like breast cancer, suggesting altered means of immune dysfunction among breast cancer subtypes. We used spatial transcriptomics and multiplex imaging to provide spatial detail for key markers of benign and malignant cell types and immune cell colocation. These findings demonstrate analysis of paired transcription and chromatin accessibility at the single cell level is a powerful tool for investigating breast cancer lineage development and highlight transcriptional networks that define basal and luminal breast cancer lineages.

Single-cell and spatial analyses reveal a tradeoff between murine mammary proliferation and lineage programs associated with endocrine cues

Although distinct epithelial cell types have been distinguished in glandular tissues such as the mammary gland, the extent of heterogeneity within each cell type and the degree of endocrine control of this diversity across development are incompletely understood. By combining mass cytometry and cyclic immunofluorescence, we define a rich array of murine mammary epithelial cell subtypes associated with puberty, the estrous cycle, and sex. These subtypes are differentially proliferative and spatially segregate distinctly in adult versus pubescent glands. Further, we identify systematic suppression of lineage programs at the protein and RNA levels as a common feature of mammary epithelial expansion during puberty, the estrous cycle, and gestation and uncover a pervasive enrichment of ribosomal protein genes in luminal cells elicited specifically during progesterone-dominant expansionary periods. Collectively, these data expand our knowledge of murine mammary epithelial heterogeneity and connect endocrine-driven epithelial expansion with lineage suppression.

Comprehensive single cell aging atlas of mammary tissues reveals shared epigenomic and transcriptomic signatures of aging and cancer

Aging is the greatest risk factor for breast cancer; however, how age-related cellular and molecular events impact cancer initiation is unknown. We investigate how aging rewires transcriptomic and epigenomic programs of mouse mammary glands at single cell resolution, yielding a comprehensive resource for aging and cancer biology. Aged epithelial cells exhibit epigenetic and transcriptional changes in metabolic, pro-inflammatory, or cancer-associated genes. Aged stromal cells downregulate fibroblast marker genes and upregulate markers of senescence and cancer-associated fibroblasts. Among immune cells, distinct T cell subsets (Gzmk+, memory CD4+, γδ) and M2-like macrophages expand with age. Spatial transcriptomics reveal co-localization of aged immune and epithelial cells in situ. Lastly, transcriptional signatures of aging mammary cells are found in human breast tumors, suggesting mechanistic links between aging and cancer. Together, these data uncover that epithelial, immune, and stromal cells shift in proportions and cell identity, potentially impacting cell plasticity, aged microenvironment, and neoplasia risk.

Implications of tumour heterogeneity on cancer evolution and therapy resistance: lessons from breast cancer

Tumour heterogeneity is pervasive amongst many cancers and leads to disease progression, and therapy resistance. In this review, using breast cancer as an exemplar, we focus on the recent advances in understanding the interplay between tumour cells and their microenvironment using single cell sequencing and digital spatial profiling technologies. Further, we discuss the utility of lineage tracing methodologies in pre-clinical models of breast cancer, and how these are being used to unravel new therapeutic vulnerabilities and reveal biomarkers of breast cancer progression. © 2023 The Pathological Society of Great Britain and Ireland.

Imagine beyond: recent breakthroughs and next challenges in mammary gland biology and breast cancer research

On 8 December 2022 the organizing committee of the European Network for Breast Development and Cancer labs (ENBDC) held its fifth annual Think Tank meeting in Amsterdam, the Netherlands. Here, we embraced the opportunity to look back to identify the most prominent breakthroughs of the past ten years and to reflect on the main challenges that lie ahead for our field in the years to come. The outcomes of these discussions are presented in this position paper, in the hope that it will serve as a summary of the current state of affairs in mammary gland biology and breast cancer research for early career researchers and other newcomers in the field, and as inspiration for scientists and clinicians to move the field forward.

Single-cell chromatin accessibility of developing murine pancreas identifies cell state-specific gene regulatory programs

Numerous studies have characterized the existence of cell subtypes, along with their corresponding transcriptional profiles, within the developing mouse pancreas. The upstream mechanisms that initiate and maintain gene expression programs across cell states, however, remain largely unknown. Here, we generate single-nucleus ATAC-Sequencing data of developing murine pancreas and perform an integrated, multi-omic analysis of both chromatin accessibility and RNA expression to describe the chromatin landscape of the developing pancreas at both E14.5 and E17.5 at single-cell resolution. We identify candidate transcription factors regulating cell fate and construct gene regulatory networks of active transcription factor binding to regulatory regions of downstream target genes. This work serves as a valuable resource for the field of pancreatic biology in general and contributes to our understanding of lineage plasticity among endocrine cell types. In addition, these data identify which epigenetic states should be represented in the differentiation of stem cells to the pancreatic beta cell fate to best recapitulate in vitro the gene regulatory networks that are critical for progression along the beta cell lineage in vivo.

Chromatin Organization and Transcriptional Programming of Breast Cancer Cell Identity

The advent of sequencing technologies for assessing chromosome conformations has provided a wealth of information on the organization of the 3-dimensional genome and its role in cancer progression. It is now known that changes in chromatin folding and accessibility can promote aberrant activation or repression of transcriptional programs that can drive tumorigenesis and progression in diverse cancers. This includes breast cancer, which comprises several distinct subtypes defined by their unique transcriptomes that dictate treatment response and patient outcomes. Of these, basal-like breast cancer is an aggressive subtype controlled by a pluripotency-enforcing transcriptome. Meanwhile, the more differentiated luminal subtype of breast cancer is driven by an estrogen receptor-dominated transcriptome that underlies its responsiveness to antihormone therapies and conveys improved patient outcomes. Despite the clear differences in molecular signatures, the genesis of each subtype from normal mammary epithelial cells remains unclear. Recent technical advances have revealed key distinctions in chromatin folding and organization between subtypes that could underlie their transcriptomic and, hence, phenotypic differences. These studies also suggest that proteins controlling particular chromatin states may be useful targets for treating aggressive disease. In this review, we explore the current state of understanding of chromatin architecture in breast cancer subtypes and its potential role in defining their phenotypic characteristics.

Single-cell chromatin landscapes of mouse skin development

The coat of mammals is produced by hair follicles, and hair follicle is an important and complex accessory organ of skin. As a complex physiological regulation process, hair follicle morphogenesis is regulated by a series of signal pathway factors, involves the interaction of multiple cell types and begins in the early embryonic stage. However, its transcriptional regulatory mechanism is unclear. We have therefore utilized single-cell ATAC sequencing to obtain the chromatin accessibility landscapes of 6,928, 6,961 and 7,374 high-quality cells from the dorsal skins of E13.5, E16.5 and P0 mice (Mus musculus), respectively. Based on marker gene activity clustering, we defined 6, 8 and 5 distinct cell types in E13.5, E16.5 and P0 stages, respectively. Furtherly, we integrated the fibroblasts and keratinocytes clusters, performed further analysis and re-clustered. The single cell map of the chromatin open area was drawn from each cell type and the mechanism of cell transcription regulation was explored. Collectively, our data provide a reference for deeply exploring the epigenetic regulation mechanism of mouse hair follicles development.

Loss of Epigenetic Regulation Disrupts Lineage Integrity, Induces Aberrant Alveogenesis, and Promotes Breast Cancer

Systematically investigating the scores of genes mutated in cancer and discerning disease drivers from inconsequential bystanders is a prerequisite for precision medicine but remains challenging. Here, we developed a somatic CRISPR/Cas9 mutagenesis screen to study 215 recurrent "long-tail" breast cancer genes, which revealed epigenetic regulation as a major tumor-suppressive mechanism. We report that components of the BAP1 and COMPASS-like complexes, including KMT2C/D, KDM6A, BAP1, and ASXL1/2 ("EpiDrivers"), cooperate with PIK3CAH1047R to transform mouse and human breast epithelial cells. Mechanistically, we find that activation of PIK3CAH1047R and concomitant EpiDriver loss triggered an alveolar-like lineage conversion of basal mammary epithelial cells and accelerated formation of luminal-like tumors, suggesting a basal origin for luminal tumors. EpiDriver mutations are found in ∼39% of human breast cancers, and ∼50% of ductal carcinoma in situ express casein, suggesting that lineage infidelity and alveogenic mimicry may significantly contribute to early steps of breast cancer etiology.

SIGNIFICANCE

Infrequently mutated genes comprise most of the mutational burden in breast tumors but are poorly understood. In vivo CRISPR screening identified functional tumor suppressors that converged on epigenetic regulation. Loss of epigenetic regulators accelerated tumorigenesis and revealed lineage infidelity and aberrant expression of alveogenesis genes as potential early events in tumorigenesis. This article is highlighted in the In This Issue feature, p. 2711.

Breast cancer plasticity is restricted by a LATS1-NCOR1 repressive axis

Breast cancer, the most frequent cancer in women, is generally classified into several distinct histological and molecular subtypes. However, single-cell technologies have revealed remarkable cellular and functional heterogeneity across subtypes and even within individual breast tumors. Much of this heterogeneity is attributable to dynamic alterations in the epigenetic landscape of the cancer cells, which promote phenotypic plasticity. Such plasticity, including transition from luminal to basal-like cell identity, can promote disease aggressiveness. We now report that the tumor suppressor LATS1, whose expression is often downregulated in human breast cancer, helps maintain luminal breast cancer cell identity by reducing the chromatin accessibility of genes that are characteristic of a "basal-like" state, preventing their spurious activation. This is achieved via interaction of LATS1 with the NCOR1 nuclear corepressor and recruitment of HDAC1, driving histone H3K27 deacetylation near NCOR1-repressed "basal-like" genes. Consequently, decreased expression of LATS1 elevates the expression of such genes and facilitates slippage towards a more basal-like phenotypic identity. We propose that by enforcing rigorous silencing of repressed genes, the LATS1-NCOR1 axis maintains luminal cell identity and restricts breast cancer progression.

Spatially resolved gene regulatory and disease-related vulnerability map of the adult Macaque cortex

Single cell approaches have increased our knowledge about the cell type composition of the non-human primate (NHP), but a detailed characterization of area-specific regulatory features remains outstanding. We generated single-cell transcriptomic and chromatin accessibility (single-cell ATAC) data of 358,237 cells from prefrontal cortex (PFC), primary motor cortex (M1) and primary visual cortex (V1) of adult female cynomolgus monkey brain, and integrated this dataset with Stereo-seq (spatial enhanced resolution omics-sequencing) of the corresponding cortical areas to assign topographic information to molecular states. We identified area-specific chromatin accessible sites and their targeted genes, including the cell type-specific transcriptional regulatory network associated with excitatory neurons heterogeneity. We reveal calcium ion transport and axon guidance genes related to specialized functions of PFC and M1, identified the similarities and differences between adult macaque and human oligodendrocyte trajectories, and mapped the genetic variants and gene perturbations of human diseases to NHP cortical cells. This resource establishes a transcriptomic and chromatin accessibility combinatory regulatory landscape at a single-cell and spatially resolved resolution in NHP cortex.

Correlation of liver and kidney indicators with foetal vital organ function

This study aimed to investigate the correlation between indicators of liver and kidney function and foetal vital organ function. One hundred and eighty-five pregnant women who underwent cordocentesis and whose foetuses were diagnosed with abnormal foetal organ function were enrolled. The indicators of liver and kidney function were compared between foetuses with abnormal vital organ function and healthy foetuses. There was a significant difference between foetuses with and those without normal cardiovascular systems in terms of total protein, albumin, total bile acid, and creatinine levels (P < .05). A significant difference in aspartate aminotransferase (AST) and lactate dehydrogenase (LDH) levels was observed in foetuses with and those without normal foetal urinary systems (P < .05). A difference between foetuses with normal and those without normal musculoskeletal systems was noted when comparing LDH levels. Further, there was a significant difference in gestational age and AST, alanine aminotransferase, albumin, total bilirubin, alkaline phosphatase, LDH, adenosine dehydrogenase, fibronectin, and creatinine levels between foetuses with normal versus abnormal blood systems (P < .05). Thus, hepatic and renal function indicators may be associated with abnormal foetal vital organ function.Impact statementWhat is already known on this subject? Foetal cardiac function is currently evaluated using colour Doppler ultrasound and magnetic resonance imaging in clinical practice, but there are few predictive indicators of the function of other vital organs. It is difficult to determine whether children have abnormalities in the urinary system, digestive system, nervous system, or other vital organs.What do the results of this study add? In this study, it was found that total protein, albumin, total bile acid, creatinine, aspartate aminotransferase, lactate dehydrogenase, fibronectin, alanine aminotransferase, total bilirubin, alkaline phosphatase, adenosine dehydrogenase, and other liver and kidney function indicators may be associated with foetal vital organ dysfunction. However, the forecast range of specific indicators must be further improved upon.What are the implications of these findings for clinical practice and/or further research? This study provides an additional reference for predicting foetal cardiac function.

Newly identified breast luminal progenitor and gestational stem cell populations likely give rise to HER2-overexpressing and basal-like breast cancers

Breast Cancer (BrC) is a common malignancy with genetically diverse subtypes. There is evidence that specific BrC subtypes originate from particular normal mammary cell populations. However, the cell populations that give rise to most BrC subtypes are unidentified. Several human breast scRNAseq datasets are available. In this research, we utilized a robust human scRNAseq dataset to identify population-specific marker genes and then identified the expression of these marker genes in specific BrC subtypes. In humans, several BrC subtypes, HER2-enriched, basal-like, and triple-negative (TN), are more common in women who have had children. This observation suggests that cell populations that originate during pregnancy give rise to these BrCs. The current human datasets have few normal parous samples, so we supplemented this research with mouse datasets, which contain mammary cells from various developmental stages. This research identified two novel normal breast cell populations that may be the origin of the basal-like and HER2-overexpressing subtypes, respectively. A stem cell-like population, SC, that expresses gestation-specific genes has similar gene expression patterns to basal-like BrCs. A novel luminal progenitor cell population and HER2-overexpressing BrCs are marked by S100A7, S100A8, and S100A9 expression. We bolstered our findings by examining SC gene expression in TN BrC scRNAseq datasets and S100A7-A9 gene expression in BrC cell lines. We discovered that several potential cancer stem cell populations highly express most of the SC genes in TN BrCs and confirmed S100A8 and A9 overexpression in a HER2-overexpressing BrC cell line. In summary, normal SC and the novel luminal progenitor cell population likely give rise to basal-like and HER2-overexpressing BrCs, respectively. Characterizing these normal cell populations may facilitate a better understanding of specific BrCs subtypes.

A human breast atlas integrating single-cell proteomics and transcriptomics

The breast is a dynamic organ whose response to physiological and pathophysiological conditions alters its disease susceptibility, yet the specific effects of these clinical variables on cell state remain poorly annotated. We present a unified, high-resolution breast atlas by integrating single-cell RNA-seq, mass cytometry, and cyclic immunofluorescence, encompassing a myriad of states. We define cell subtypes within the alveolar, hormone-sensing, and basal epithelial lineages, delineating associations of several subtypes with cancer risk factors, including age, parity, and BRCA2 germline mutation. Of particular interest is a subset of alveolar cells termed basal-luminal (BL) cells, which exhibit poor transcriptional lineage fidelity, accumulate with age, and carry a gene signature associated with basal-like breast cancer. We further utilize a medium-depletion approach to identify molecular factors regulating cell-subtype proportion in organoids. Together, these data are a rich resource to elucidate diverse mammary cell states.

Alveolar cells in the mammary gland: lineage commitment and cell death

The mammary gland provides a spectacular example of physiological cell death whereby the cells that produce milk during lactation are removed swiftly, efficiently, and without inducing inflammation upon the cessation of lactation. The milk-producing cells arise primarily during pregnancy and comprise the alveolar lineage that is specified by signalling pathways and factors that are activated in response to pregnancy hormones. There are at least two alveolar sub-lineages, one of which is marked by the presence of binucleate cells that are especially susceptible to programmed cell death during involution. This process of post-lactational regression, or involution, is carefully orchestrated and occurs in two phases, the first results in a rapid switch in cell fate with the secretory epithelial cells becoming phagocytes whereupon they destroy dead and dying cells from milk. This reversible phase is followed by the second phase that is marked by an influx of immune cells and a remodelling of the gland to replace the alveolar cells with re-differentiated adipocytes, resulting in a return to the pre-pregnant state in preparation for any subsequent pregnancies. The mouse mammary gland provides an excellent experimental tool with which to investigate lineage commitment and the mechanisms of programmed cell death that occur in a normal physiological process. Importantly, involution has highlighted a role for lysoptosis, a mechanism of cell death that is mediated by lysosomal cathepsins and their endogenous inhibitors, serpins. In this review, I discuss alveolar lineage commitment during pregnancy and the programmed cell death pathways that destroy these cells during involution.

The Mammary Gland: Basic Structure and Molecular Signaling during Development

The mammary gland is a compound, branched tubuloalveolar structure and a major characteristic of mammals. The mammary gland has evolved from epidermal apocrine glands, the skin glands as an accessory reproductive organ to support postnatal survival of offspring by producing milk as a source of nutrition. The mammary gland development begins during embryogenesis as a rudimentary structure that grows into an elementary branched ductal tree and is embedded in one end of a larger mammary fat pad at birth. At the onset of ovarian function at puberty, the rudimentary ductal system undergoes dramatic morphogenetic change with ductal elongation and branching. During pregnancy, the alveolar differentiation and tertiary branching are completed, and during lactation, the mature milk-producing glands eventually develop. The early stages of mammary development are hormonal independent, whereas during puberty and pregnancy, mammary gland development is hormonal dependent. We highlight the current understanding of molecular regulators involved during different stages of mammary gland development.

Gpr125 is a unifying hallmark of multiple mammary progenitors coupled to tumor latency

Gpr125 is an orphan G-protein coupled receptor, with homology to cell adhesion and axonal guidance factors, that is implicated in planar polarity and control of cell movements. By lineage tracing we demonstrate that Gpr125 is a highly specific marker of bipotent mammary stem cells in the embryo and of multiple long-lived unipotent basal mammary progenitors in perinatal and postnatal glands. Nipple-proximal Gpr125+ cells express a transcriptomic profile indicative of chemo-repulsion and cell movement, whereas Gpr125+ cells concentrated at invasive ductal tips display a hybrid epithelial-mesenchymal phenotype and are equipped to bind chemokine and growth factors and secrete a promigratory matrix. Gpr125 progenitors acquire bipotency in the context of transplantation and cancer and are greatly expanded and massed at the pushing margins of short latency MMTV-Wnt1 tumors. High Gpr125 expression identifies patients with particularly poor outcome within the basal breast cancer subtype highlighting its potential utility as a factor to stratify risk.

Gpr125 has emerged as a specific marker of mammary stem cells and basal progenitors. Here they show that Gpr125 cells congregate at ductal tips during morphogenesis and amass at tumor margins, and that high Gpr125 predicts early tumor onset and poor outcome in basal breast cancer.

Mechanism of activation and the rewired network: New drug design concepts

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

Single-Cell Analysis Unveils the Role of the Tumor Immune Microenvironment and Notch Signaling in Dormant Minimal Residual Disease

Tumor dormancy is a stage in which residual cancer cells remain inactive, but regrowth of dormant cancer cells contributes to recurrence. The complex ecosystem in cancer that promotes cell survival and the factors that eventually overcome growth constraints and result in proliferation remain to be fully elucidated. Doing so may provide new insights and help identify novel strategies to prolong cancer dormancy and prevent disease recurrence. To dissect the molecular pathways and the microenvironments involved in regulation of dormancy, we utilized a novel immunocompetent transgenic model to study minimal residual disease and relapse. This model revealed a significant reorganization of cancer cell structures, stroma, and immune cells, with cancer cells showing dormant cell signatures. Single-cell RNA sequencing uncovered remodeling of myeloid and lymphoid compartments. In addition, the Jagged-1/Notch signaling pathway was shown to regulate many aspects of tumorigenesis, including stem cell development, epithelial-to-mesenchymal transition, and immune cell homeostasis during minimal residual disease. Treatment with an anti-Jagged-1 antibody inhibited the Jagged-1/Notch signaling pathway in tumor cells and the microenvironment, delaying tumor recurrence. These findings uncover a cascade of regulatory changes in the microenvironment during dormancy and identify a therapeutic strategy to undercut these changes.

SIGNIFICANCE

Single-cell RNA-sequencing analysis reveals dormancy-associated changes in immune and stromal cells and demonstrates a rationale to pursue Jagged-1/Notch pathway inhibition as a viable therapeutic strategy to reduce disease recurrence.

Coordinate control of basal epithelial cell fate and stem cell maintenance by core EMT transcription factor Zeb1

Maintenance of undifferentiated, long-lived, and often quiescent stem cells in the basal compartment is important for homeostasis and regeneration of multiple epithelial tissues, but the molecular mechanisms that coordinately control basal cell fate and stem cell quiescence are elusive. Here, we report an epithelium-intrinsic requirement for Zeb1, a core transcriptional inducer of epithelial-to-mesenchymal transition, for mammary epithelial ductal side branching and for basal cell regenerative capacity. Our findings uncover an evolutionarily conserved role of Zeb1 in promoting basal cell fate over luminal differentiation. We show that Zeb1 loss results in increased basal cell proliferation at the expense of quiescence and self-renewal. Moreover, Zeb1 cooperates with YAP to activate Axin2 expression, and inhibition of Wnt signaling partially restores stem cell function to Zeb1-deficient basal cells. Thus, Zeb1 is a transcriptional regulator that maintains both basal cell fate and stem cell quiescence, and it functions in part through suppressing Wnt signaling.

C/EBPB-dependent adaptation to palmitic acid promotes tumor formation in hormone receptor negative breast cancer

Epidemiological studies have established a positive association between obesity and the incidence of postmenopausal breast cancer. Moreover, it is known that obesity promotes stem cell-like properties of breast cancer cells. However, the cancer cell-autonomous mechanisms underlying this correlation are not well defined. Here we demonstrate that obesity-associated tumor formation is driven by cellular adaptation rather than expansion of pre-existing clones within the cancer cell population. While there is no correlation with specific mutations, cellular adaptation to obesity is governed by palmitic acid (PA) and leads to enhanced tumor formation capacity of breast cancer cells. This process is governed epigenetically through increased chromatin occupancy of the transcription factor CCAAT/enhancer-binding protein beta (C/EBPB). Obesity-induced epigenetic activation of C/EBPB regulates cancer stem-like properties by modulating the expression of key downstream regulators including CLDN1 and LCN2. Collectively, our findings demonstrate that obesity drives cellular adaptation to PA drives tumor initiation in the obese setting through activation of a C/EBPB dependent transcriptional network.

Obesity is linked to cancer risk in post-menopausal breast cancer. At the molecular level this is governed by cellular adaption to palmitic acid through epigenetic activation of a C/EBPB-dependent transcriptional network that drives tumor formation.

Single-Cell Transcriptomic and Epigenetic Analyses of Mouse Mammary Development Starting with the Embryo

Cancers are caricatures of normal development. Yet, for most organs we are only beginning to learn about the molecular events underlying the embryonic antecedents of organogenesis and when differentiation into the cell types found in the adult actually begins. Here, we will focus on the powerful single-cell RNA sequencing and Assay for Transposase Accessible DNA by DNA sequencing (ATAC-seq) that we and others have been using to decipher the key regulators and signal transduction pathways involved in normal mammary development. We will first describe the techniques we use to identify, dissect, and isolate embryonic mammary rudiments and their constituent cells. We then describe the methods we have employed to perform single-cell RNA-seq and single-nucleus ATAC-seq using the small number of cells obtainable from mouse embryos. Finally, we will discuss the bioinformatic techniques we have used to interpret the vast amount of data obtained with these methods.

Parity-induced changes to mammary epithelial cells control NKT cell expansion and mammary oncogenesis

Pregnancy reprograms mammary epithelial cells (MECs) to control their responses to pregnancy hormone re-exposure and carcinoma progression. However, the influence of pregnancy on the mammary microenvironment is less clear. Here, we used single-cell RNA sequencing to profile the composition of epithelial and non-epithelial cells in mammary tissue from nulliparous and parous female mice. Our analysis indicates an expansion of γδ natural killer T-like immune cells (NKTs) following pregnancy and upregulation of immune signaling molecules in post-pregnancy MECs. We show that expansion of NKTs following pregnancy is due to elevated expression of the antigen-presenting molecule CD1d on MECs. Loss of CD1d expression on post-pregnancy MECs, or overall lack of activated NKTs, results in mammary oncogenesis. Collectively, our findings illustrate how pregnancy-induced changes modulate the communication between MECs and the immune microenvironment and establish a causal link between pregnancy, the immune microenvironment, and mammary oncogenesis.

Single-cell evaluation reveals shifts in the tumor-immune niches that shape and maintain aggressive lesions in the breast

There is an unmet clinical need for stratification of breast lesions as indolent or aggressive to tailor treatment. Here, single-cell transcriptomics and multiparametric imaging applied to a mouse model of breast cancer reveals that the aggressive tumor niche is characterized by an expanded basal-like population, specialization of tumor subpopulations, and mixed-lineage tumor cells potentially serving as a transition state between luminal and basal phenotypes. Despite vast tumor cell-intrinsic differences, aggressive and indolent tumor cells are functionally indistinguishable once isolated from their local niche, suggesting a role for non-tumor collaborators in determining aggressiveness. Aggressive lesions harbor fewer total but more suppressed-like T cells, and elevated tumor-promoting neutrophils and IL-17 signaling, disruption of which increase tumor latency and reduce the number of aggressive lesions. Our study provides insight into tumor-immune features distinguishing indolent from aggressive lesions, identifies heterogeneous populations comprising these lesions, and supports a role for IL-17 signaling in aggressive progression.

Progesterone receptor antagonists reverse stem cell expansion and the paracrine effectors of progesterone action in the mouse mammary gland

Background

The ovarian hormones estrogen and progesterone (EP) are implicated in breast cancer causation. A specific consequence of progesterone exposure is the expansion of the mammary stem cell (MSC) and luminal progenitor (LP) compartments. We hypothesized that this effect, and its molecular facilitators, could be abrogated by progesterone receptor (PR) antagonists administered in a mouse model.

Methods

Ovariectomized FVB mice were randomized to 14 days of treatment: sham, EP, EP + telapristone (EP + TPA), EP + mifepristone (EP + MFP). Mice were then sacrificed, mammary glands harvested, and mammary epithelial cell lineages separated by flow cytometry using cell surface markers. RNA from each lineage was sequenced and differential gene expression was analyzed using DESeq. Quantitative PCR was performed to confirm the candidate genes discovered in RNA seq. ANOVA with Tukey post hoc analysis was performed to compare relative expression. Alternative splicing events were examined using the rMATs multivariate analysis tool.

Results

Significant increases in the MSC and luminal mature (LM) cell fractions were observed following EP treatment compared to control (p < 0.01 and p < 0.05, respectively), whereas the LP fraction was significantly reduced (p < 0.05). These hormone-induced effects were reversed upon exposure to TPA and MFP (p < 0.01 for both). Gene Ontology analysis of RNA-sequencing data showed EP-induced enrichment of several pathways, with the largest effect on Wnt signaling in MSC, significantly repressed by PR inhibitors. In LP cells, significant induction of Wnt4 and Rankl, and Wnt pathway intermediates Lrp2 and Axin2 (confirmed by qRTPCR) were reversed by TPA and MFP (p < 0.0001). Downstream signaling intermediates of these pathways (Lrp5, Mmp7) showed similar effects. Expression of markers of epithelial-mesenchymal transition (Cdh1, Cdh3) and the induction of EMT regulators (Zeb1, Zeb2, Gli3, Snai1, and Ptch2) were significantly responsive to progesterone. EP treatment was associated with large-scale alternative splicing events, with an enrichment of motifs associated with Srsf, Esrp, and Rbfox families. Exon skipping was observed in Cdh1, Enah, and Brd4.

Conclusions

PR inhibition reverses known tumorigenic pathways in the mammary gland and suppresses a previously unknown effect of progesterone on RNA splicing events. In total, our results strengthen the case for reconsideration of PR inhibitors for breast cancer prevention.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01455-2.

Breast cancer heterogeneity through the lens of single-cell analysis and spatial pathologies

Breast cancer ecosystems are composed of complex cell types, including tumor, stromal and immune cells, each of which can assume diverse phenotypes. Both the heterogeneous composition and spatially distinct tumor microenvironment impact breast cancer progression, treatment response and therapeutic resistance. Thus, a deeper understanding of breast cancer heterogeneity may help facilitate the development of novel therapies and improve outcomes for patients. The advent of paradigm shifting single-cell analysis and spatial pathologies allows for a comprehensive analysis of the tumor ecosystem as well as the interactions between its components at unprecedented resolution. In this review, we discuss the insights gained through single-cell analysis and spatial pathologies on breast cancer heterogeneity.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01445-4.

Mammary cell gene expression atlas links epithelial cell remodeling events to breast carcinogenesis

The female mammary epithelium undergoes reorganization during development, pregnancy, and menopause, linking higher risk with breast cancer development. To characterize these periods of complex remodeling, here we report integrated 50 K mouse and 24 K human mammary epithelial cell atlases obtained by single-cell RNA sequencing, which covers most lifetime stages. Our results indicate a putative trajectory that originates from embryonic mammary stem cells which differentiates into three epithelial lineages (basal, luminal hormone-sensing, and luminal alveolar), presumably arising from unipotent progenitors in postnatal glands. The lineage-specific genes infer cells of origin of breast cancer using The Cancer Genome Atlas data and single-cell RNA sequencing of human breast cancer, as well as the association of gland reorganization to different breast cancer subtypes. This comprehensive mammary cell gene expression atlas ( https://mouse-mammary-epithelium-integrated.cells.ucsc.edu ) presents insights into the impact of the internal and external stimuli on the mammary epithelium at an advanced resolution.

The molecular basis of mammary gland development and epithelial differentiation

Our understanding of the molecular events underpinning the development of mammalian organ systems has been increasing rapidly in recent years. With the advent of new and improved next-generation sequencing methods, we are now able to dig deeper than ever before into the genomic and epigenomic events that play critical roles in determining the fates of stem and progenitor cells during the development of an embryo into an adult. In this review, we detail and discuss the genes and pathways that are involved in mammary gland development, from embryogenesis, through maturation into an adult gland, to the role of pregnancy signals in directing the terminal maturation of the mammary gland into a milk producing organ that can nurture the offspring. We also provide an overview of the latest research in the single-cell genomics of mammary gland development, which may help us to understand the lineage commitment of mammary stem cells (MaSCs) into luminal or basal epithelial cells that constitute the mammary gland. Finally, we summarize the use of 3D organoid cultures as a model system to study the molecular events during mammary gland development. Our increased investigation of the molecular requirements for normal mammary gland development will advance the discovery of targets to predict breast cancer risk and the development of new breast cancer therapies.

Mammary gland development from a single cell 'omics view

Understanding the complexity and heterogeneity of mammary cell subpopulations is vital to delineate the mechanisms behind breast cancer development, progression and prevention. Increasingly sophisticated tools for investigating these cell subtypes has led to the development of a greater understanding of these cell subtypes, complex interplay of certain subtypes and their developmental potential. Of note, increasing accessibility and affordability of single cell technologies has led to a plethora of studies being published containing data from mammary cell subtypes and their differentiation potential in both mice and human data sets. Here, we review the different types of single cell technologies and how they have been used to improve our understanding of mammary gland development.

RNA sequencing of long-term label-retaining colon cancer stem cells identifies novel regulators of quiescence

Recent data suggest that therapy-resistant quiescent cancer stem cells (qCSCs) are the source of relapse in colon cancer. Here, using colon cancer patient-derived organoids and xenografts, we identify rare long-term label-retaining qCSCs that can re-enter the cell cycle to generate new tumors. RNA sequencing analyses demonstrated that these cells display the molecular hallmarks of quiescent tissue stem cells, including expression of p53 signaling genes, and are enriched for transcripts common to damage-induced quiescent revival stem cells of the regenerating intestine. In addition, we identify negative regulators of cell cycle, downstream of p53, that we show are indicators of poor prognosis and may be targeted for qCSC abolition in both p53 wild-type and mutant tumors. These data support the temporal inhibition of downstream targets of p53 signaling, in combination with standard-of-care treatments, for the elimination of qCSCs and prevention of relapse in colon cancer.

•

Colon tumors contain therapy-resistant quiescent cancer stem cells (qCSCs)

•

qCSC gene expression mirrors that of quiescent stem cells of the regenerating gut

•

qCSCs are enriched for p53 signaling genes

•

qCSC elimination may be achieved by inhibiting downstream targets of p53 signaling

Sox11 regulates mammary tumour-initiating and metastatic capacity in Brca1-deficient mouse mammary tumour cells

Little is known about the role of Sox11 in the regulation of mammary progenitor cells. Sox11 is expressed by mammary bud epithelial cells during embryonic mammary gland development and is not detected in mammary epithelial cells after birth. As Sox11 is an oncofetal gene, we investigated the effects of reducing Sox11 levels in embryonic mammary progenitor cells and found that Sox11 regulates proliferative state, stem cell activity and lineage marker expression. We also investigated the effect of reducing Sox11 levels in two transplantable Brca1-deficient oestrogen receptor-negative mouse mammary tumour cell lines, to assess whether Sox11 regulates similar functions in tumour progenitor cells. When Sox11 levels were reduced in one Brca1-deficient mammary tumour cell line that expressed both epithelial and mesenchymal markers, similar effects on proliferation, stem cell activity and expression of lineage markers to those seen in the embryonic mammary progenitor cells were observed. Orthotopic grafting of mammary tumour cells with reduced Sox11 levels led to alterations in tumour-initiating capacity, latency, expression of lineage markers and metastatic burden. Our results support a model in which tumours expressing higher levels of Sox11 have more stem and tumour-initiating cells, and are less proliferative, whereas tumours expressing lower levels of Sox11 become more proliferative and capable of morphogenetic/metastatic growth, similar to what occurs during embryonic mammary developmental progression.

Summary:Brca1^−/− mammary tumours expressing Sox11 at high levels have more stem- and tumour-initiating cells, and are less proliferative, whereas tumours expressing Sox11 at lower levels become more proliferative and are capable of morphogenetic/metastatic growth.