A human breast atlas integrating single-cell proteomics and transcriptomics

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.

Luminal breast epithelial cells from wildtype and BRCA mutation carriers harbor copy number alterations commonly associated with breast cancer

Cancer-associated mutations have been documented in normal tissues, but the prevalence and nature of somatic copy number alterations and their role in tumor initiation and evolution is not well understood. Here, using single cell DNA sequencing, we describe the landscape of CNAs in >42,000 breast epithelial cells from women with normal or high risk of developing breast cancer. Accumulation of individual cells with one or two of a specific subset of CNAs (e.g. 1q gain and 16q, 22q, 7q, and 10q loss) is detectable in almost all breast tissues and, in those from BRCA1 or BRCA2 mutations carriers, occurs prior to loss of heterozygosity (LOH) of the wildtype alleles. These CNAs, which are among the most common associated with ductal carcinoma in situ (DCIS) and malignant breast tumors, are enriched almost exclusively in luminal cells not basal myoepithelial cells. Allele-specific analysis of the enriched CNAs reveals that each allele was independently altered, demonstrating convergent evolution of these CNAs in an individual breast. Tissues from BRCA1 or BRCA2 mutation carriers contain a small percentage of cells with extreme aneuploidy, featuring loss of TP53 , LOH of BRCA1 or BRCA2 , and multiple breast cancer-associated CNAs in addition to one or more of the common CNAs in 1q, 10q or 16q. Notably, cells with intermediate levels of CNAs are not detected, arguing against a stepwise gradual accumulation of CNAs. Overall, our findings demonstrate that chromosomal alterations in normal breast epithelium partially mirror those of established cancer genomes and are chromosome- and cell lineage-specific.

Multimodal joint deconvolution and integrative signature selection in proteomics

Deconvolution is an efficient approach for detecting cell-type-specific (cs) transcriptomic signals without cellular segmentation. However, this type of methods may require a reference profile from the same molecular source and tissue type. Here, we present a method to dissect bulk proteome by leveraging tissue-matched transcriptome and proteome without using a proteomics reference panel. Our method also selects the proteins contributing to the cellular heterogeneity shared between bulk transcriptome and proteome. The deconvoluted result enables downstream analyses such as cs-protein Quantitative Trait Loci (cspQTL) mapping. We benchmarked the performance of this multimodal deconvolution approach through CITE-seq pseudo bulk data, a simulation study, and the bulk multi-omics data from human brain normal tissues and breast cancer tumors, individually, showing robust and accurate cell abundance quantification across different datasets. This algorithm is implemented in a tool MICSQTL that also provides cspQTL and multi-omics integrative visualization, available at https://bioconductor.org/packages/MICSQTL .

A single-cell atlas enables mapping of homeostatic cellular shifts in the adult human breast

Here we use single-cell RNA sequencing to compile a human breast cell atlas assembled from 55 donors that had undergone reduction mammoplasties or risk reduction mastectomies. From more than 800,000 cells we identified 41 cell subclusters across the epithelial, immune and stromal compartments. The contribution of these different clusters varied according to the natural history of the tissue. Age, parity and germline mutations, known to modulate the risk of developing breast cancer, affected the homeostatic cellular state of the breast in different ways. We found that immune cells from BRCA1 or BRCA2 carriers had a distinct gene expression signature indicative of potential immune exhaustion, which was validated by immunohistochemistry. This suggests that immune-escape mechanisms could manifest in non-cancerous tissues very early during tumor initiation. This atlas is a rich resource that can be used to inform novel approaches for early detection and prevention of breast cancer.

Integrated analysis reveals FLI1 regulates the tumor immune microenvironment via its cell-type-specific expression and transcriptional regulation of distinct target genes of immune cells in breast cancer

BACKGROUND

Immunotherapy is a practical therapeutic approach in breast cancer (BRCA), and the role of FLI1 in immune regulation has gradually been unveiled. However, the specific role of FLI1 in BRCA was conflicted; thus, additional convincing evidence is needed.

METHODS

We explored the upstream regulation of FLI1 expression via summary data-based Mendelian randomization (SMR) analysis and ncRNA network construction centering on FLI1 using BRCA genome-wide association study (GWAS) summary data with expression quantitative trait loci (eQTLs) and DNA methylation quantitative trait loci (mQTLs) from the blood and a series of in silico analyses, respectively. We illuminated the downstream function of FLI1 in immune regulation by integrating a series of analyses of single-cell RNA sequence data (scRNA-seq).

RESULTS

We verified a causal pathway from FLI1 methylation to FLI1 gene expression to BRCA onset and demonstrated that FLI1 was downregulated in BRCA. FLI1, a transcription factor, served as myeloid and T cells' communication regulator by targeting immune-related ligands and receptor transcription in BRCA tissues. We constructed a ceRNA network centering on FLI1 that consisted of three LncRNAs (CKMT2-AS1, PSMA3-AS1, and DIO3OS) and a miRNA (hsa-miR-324-5p), and the expression of FLI1 was positively related to a series of immune-related markers, including immune cell infiltration, biomarkers of immune cells, and immune checkpoints.

CONCLUSION

Low-methylation-induced or ncRNA-mediated downregulation of FLI1 is associated with poor prognosis, and FLI1 might regulate the tumor immune microenvironment via a cell-type-specific target genes manner in BRCA.

The androgen receptor interacts with GATA3 to transcriptionally regulate a luminal epithelial cell phenotype in breast cancer

BACKGROUND

The androgen receptor (AR) is a tumor suppressor in estrogen receptor (ER) positive breast cancer, a role sustained in some ER negative breast cancers. Key factors dictating AR genomic activity in a breast context are largely unknown. Herein, we employ an unbiased chromatin immunoprecipitation-based proteomic technique to identify endogenous AR interacting co-regulatory proteins in ER positive and negative models of breast cancer to gain new insight into mechanisms of AR signaling in this disease.

RESULTS

The DNA-binding factor GATA3 is identified and validated as a novel AR interacting protein in breast cancer cells irrespective of ER status. AR activation by the natural ligand 5α-dihydrotestosterone (DHT) increases nuclear AR-GATA3 interactions, resulting in AR-dependent enrichment of GATA3 chromatin binding at a sub-set of genomic loci. Silencing GATA3 reduces but does not prevent AR DNA binding and transactivation of genes associated with AR/GATA3 co-occupied loci, indicating a co-regulatory role for GATA3 in AR signaling. DHT-induced AR/GATA3 binding coincides with upregulation of luminal differentiation genes, including EHF and KDM4B, established master regulators of a breast epithelial cell lineage. These findings are validated in a patient-derived xenograft model of breast cancer. Interaction between AR and GATA3 is also associated with AR-mediated growth inhibition in ER positive and ER negative breast cancer.

CONCLUSIONS

AR and GATA3 interact to transcriptionally regulate luminal epithelial cell differentiation in breast cancer regardless of ER status. This interaction facilitates the tumor suppressor function of AR and mechanistically explains why AR expression is associated with less proliferative, more differentiated breast tumors and better overall survival in breast cancer.

Single-Cell Transcription Mapping of Murine and Human Mammary Organoids Responses to Female Hormones

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

Onco-Breastomics: An Eco-Evo-Devo Holistic Approach

Known as a diverse collection of neoplastic diseases, breast cancer (BC) can be hyperbolically characterized as a dynamic pseudo-organ, a living organism able to build a complex, open, hierarchically organized, self-sustainable, and self-renewable tumor system, a population, a species, a local community, a biocenosis, or an evolving dynamical ecosystem (i.e., immune or metabolic ecosystem) that emphasizes both developmental continuity and spatio-temporal change. Moreover, a cancer cell community, also known as an oncobiota, has been described as non-sexually reproducing species, as well as a migratory or invasive species that expresses intelligent behavior, or an endangered or parasite species that fights to survive, to optimize its features inside the host's ecosystem, or that is able to exploit or to disrupt its host circadian cycle for improving the own proliferation and spreading. BC tumorigenesis has also been compared with the early embryo and placenta development that may suggest new strategies for research and therapy. Furthermore, BC has also been characterized as an environmental disease or as an ecological disorder. Many mechanisms of cancer progression have been explained by principles of ecology, developmental biology, and evolutionary paradigms. Many authors have discussed ecological, developmental, and evolutionary strategies for more successful anti-cancer therapies, or for understanding the ecological, developmental, and evolutionary bases of BC exploitable vulnerabilities. Herein, we used the integrated framework of three well known ecological theories: the Bronfenbrenner's theory of human development, the Vannote's River Continuum Concept (RCC), and the Ecological Evolutionary Developmental Biology (Eco-Evo-Devo) theory, to explain and understand several eco-evo-devo-based principles that govern BC progression. Multi-omics fields, taken together as onco-breastomics, offer better opportunities to integrate, analyze, and interpret large amounts of complex heterogeneous data, such as various and big-omics data obtained by multiple investigative modalities, for understanding the eco-evo-devo-based principles that drive BC progression and treatment. These integrative eco-evo-devo theories can help clinicians better diagnose and treat BC, for example, by using non-invasive biomarkers in liquid-biopsies that have emerged from integrated omics-based data that accurately reflect the biomolecular landscape of the primary tumor in order to avoid mutilating preventive surgery, like bilateral mastectomy. From the perspective of preventive, personalized, and participatory medicine, these hypotheses may help patients to think about this disease as a process governed by natural rules, to understand the possible causes of the disease, and to gain control on their own health.

ELF5: A Molecular Clock for Breast Aging and Cancer Susceptibility

Breast cancer is predominantly an age-related disease, with aging serving as the most significant risk factor, compounded by germline mutations in high-risk genes like BRCA1/2. Aging induces architectural changes in breast tissue, particularly affecting luminal epithelial cells by diminishing lineage-specific molecular profiles and adopting myoepithelial-like characteristics. ELF5 is an important transcription factor for both normal breast and breast cancer development. This review focuses on the role of ELF5 in normal breast development, its altered expression throughout aging, and its implications in cancer. It discusses the lineage-specific expression of ELF5, its regulatory mechanisms, and its potential as a biomarker for breast-specific biological age and cancer risk.

SingPro: a knowledge base providing single-cell proteomic data

Single-cell proteomics (SCP) has emerged as a powerful tool for detecting cellular heterogeneity, offering unprecedented insights into biological mechanisms that are masked in bulk cell populations. With the rapid advancements in AI-based time trajectory analysis and cell subpopulation identification, there exists a pressing need for a database that not only provides SCP raw data but also explicitly describes experimental details and protein expression profiles. However, no such database has been available yet. In this study, a database, entitled 'SingPro', specializing in single-cell proteomics was thus developed. It was unique in (a) systematically providing the SCP raw data for both mass spectrometry-based and flow cytometry-based studies and (b) explicitly describing experimental detail for SCP study and expression profile of any studied protein. Anticipating a robust interest from the research community, this database is poised to become an invaluable repository for OMICs-based biomedical studies. Access to SingPro is unrestricted and does not mandate a login at: http://idrblab.org/singpro/.

Signaling Pathway Alterations Driven by BRCA1 and BRCA2 Germline Mutations are Sufficient to Initiate Breast Tumorigenesis by the PIK3CAH1047R Oncogene

Single-cell transcriptomics studies have begun to identify breast epithelial cell and stromal cell specific transcriptome differences between BRCA1/2 mutation carriers and non-carriers. We generated a single-cell transcriptome atlas of breast tissues from BRCA1, BRCA2 mutation carriers and compared this single-cell atlas of mutation carriers with our previously described single-cell breast atlas of healthy non-carriers. We observed that BRCA1 but not BRCA2 mutations altered the ratio between basal (basal-myoepithelial), luminal progenitor (luminal adaptive secretory precursor, LASP), and mature luminal (luminal hormone sensing) cells in breast tissues. A unique subcluster of cells within LASP cells is underrepresented in case of BRCA1 and BRCA2 mutation carriers compared with non-carriers. Both BRCA1 and BRCA2 mutations specifically altered transcriptomes in epithelial cells which are an integral part of NFκB, LARP1, and MYC signaling. Signaling pathway alterations in epithelial cells unique to BRCA1 mutations included STAT3, BRD4, SMARCA4, HIF2A/EPAS1, and Inhibin A signaling. BRCA2 mutations were associated with upregulation of IL6, PDK1, FOXO3, and TNFSF11 signaling. These signaling pathway alterations are sufficient to alter sensitivity of BRCA1/BRCA2-mutant breast epithelial cells to transformation as epithelial cells from BRCA1 mutation carriers overexpressing hTERT + PIK3CAH1047R generated adenocarcinomas, whereas similarly modified mutant BRCA2 cells generated basal carcinomas in NSG mice. Thus, our studies provide a high-resolution transcriptome atlas of breast epithelial cells of BRCA1 and BRCA2 mutation carriers and reveal their susceptibility to PIK3CA mutation-driven transformation.

SIGNIFICANCE

This study provides a single-cell atlas of breast tissues of BRCA1/2 mutation carriers and demonstrates that aberrant signaling due to BRCA1/2 mutations is sufficient to initiate breast cancer by mutant PIK3CA.

Mechanosensitive hormone signaling promotes mammary progenitor expansion and breast cancer risk

Tissue stem-progenitor cell frequency has been implicated in tumor risk and progression, but tissue-specific factors linking these associations remain ill-defined. We observed that stiff breast tissue from women with high mammographic density, who exhibit increased lifetime risk for breast cancer, associates with abundant stem-progenitor epithelial cells. Using genetically engineered mouse models of elevated integrin mechanosignaling and collagen density, syngeneic manipulations, and spheroid models, we determined that a stiff matrix and high mechanosignaling increase mammary epithelial stem-progenitor cell frequency and enhance tumor initiation in vivo. Augmented tissue mechanics expand stemness by potentiating extracellular signal-related kinase (ERK) activity to foster progesterone receptor-dependent RANK signaling. Consistently, we detected elevated phosphorylated ERK and progesterone receptors and increased levels of RANK signaling in stiff breast tissue from women with high mammographic density. The findings link fibrosis and mechanosignaling to stem-progenitor cell frequency and breast cancer risk and causally implicate epidermal growth factor receptor-ERK-dependent hormone signaling in this phenotype.

Metabolic switch regulates lineage plasticity and induces synthetic lethality in triple-negative breast cancer

Metabolic reprogramming is key for cancer development, yet the mechanism that sustains triple-negative breast cancer (TNBC) cell growth despite deficient pyruvate kinase M2 (PKM2) and tumor glycolysis remains to be determined. Here, we find that deficiency in tumor glycolysis activates a metabolic switch from glycolysis to fatty acid β-oxidation (FAO) to fuel TNBC growth. We show that, in TNBC cells, PKM2 directly interacts with histone methyltransferase EZH2 to coordinately mediate epigenetic silencing of a carnitine transporter, SLC16A9. Inhibition of PKM2 leads to impaired EZH2 recruitment to SLC16A9, and in turn de-represses SLC16A9 expression to increase intracellular carnitine influx, programming TNBC cells to an FAO-dependent and luminal-like cell state. Together, these findings reveal a new metabolic switch that drives TNBC from a metabolically heterogeneous-lineage plastic cell state to an FAO-dependent-lineage committed cell state, where dual targeting of EZH2 and FAO induces potent synthetic lethality in TNBC.

Identification of aberrant luminal progenitors and mTORC1 as a potential breast cancer prevention target in BRCA2 mutation carriers

Inheritance of a BRCA2 pathogenic variant conveys a substantial life-time risk of breast cancer. Identification of the cell(s)-of-origin of BRCA2-mutant breast cancer and targetable perturbations that contribute to transformation remains an unmet need for these individuals who frequently undergo prophylactic mastectomy. Using preneoplastic specimens from age-matched, premenopausal females, here we show broad dysregulation across the luminal compartment in BRCA2mut/+ tissue, including expansion of aberrant ERBB3lo luminal progenitor and mature cells, and the presence of atypical oestrogen receptor (ER)-positive lesions. Transcriptional profiling and functional assays revealed perturbed proteostasis and translation in ERBB3lo progenitors in BRCA2mut/+ breast tissue, independent of ageing. Similar molecular perturbations marked tumours bearing BRCA2-truncating mutations. ERBB3lo progenitors could generate both ER+ and ER- cells, potentially serving as cells-of-origin for ER-positive or triple-negative cancers. Short-term treatment with an mTORC1 inhibitor substantially curtailed tumorigenesis in a preclinical model of BRCA2-deficient breast cancer, thus uncovering a potential prevention strategy for BRCA2 mutation carriers.

Redox and metabolic reprogramming in breast cancer and cancer-associated adipose tissue

Redox and metabolic processes are tightly coupled in both physiological and pathological conditions. In cancer, their integration occurs at multiple levels and is characterized by synchronized reprogramming both in the tumor tissue and its specific but heterogeneous microenvironment. In breast cancer, the principal microenvironment is the cancer-associated adipose tissue (CAAT). Understanding how the redox-metabolic reprogramming becomes coordinated in human breast cancer is imperative both for cancer prevention and for the establishment of new therapeutic approaches. This review aims to provide an overview of the current knowledge of the redox profiles and regulation of intermediary metabolism in breast cancer while considering the tumor and CAAT of breast cancer as a unique Warburg's pseudo-organ. As cancer is now recognized as a systemic metabolic disease, we have paid particular attention to the cell-specific redox-metabolic reprogramming and the roles of estrogen receptors and circadian rhythms, as well as their crosstalk in the development, growth, progression, and prognosis of breast cancer.

Systematic benchmarking of imaging spatial transcriptomics platforms in FFPE tissues

Emerging imaging spatial transcriptomics (iST) platforms and coupled analytical methods can recover cell-to-cell interactions, groups of spatially covarying genes, and gene signatures associated with pathological features, and are thus particularly well-suited for applications in formalin fixed paraffin embedded (FFPE) tissues. Here, we benchmarked the performance of three commercial iST platforms on serial sections from tissue microarrays (TMAs) containing 23 tumor and normal tissue types for both relative technical and biological performance. On matched genes, we found that 10x Xenium shows higher transcript counts per gene without sacrificing specificity, but that all three platforms concord to orthogonal RNA-seq datasets and can perform spatially resolved cell typing, albeit with different false discovery rates, cell segmentation error frequencies, and with varying degrees of sub-clustering for downstream biological analyses. Taken together, our analyses provide a comprehensive benchmark to guide the choice of iST method as researchers design studies with precious samples in this rapidly evolving field.

Cellular Plasticity in Mammary Gland Development and Breast Cancer

Cellular plasticity is a phenomenon where cells adopt different identities during development and tissue homeostasis as a response to physiological and pathological conditions. This review provides a general introduction to processes by which cells change their identity as well as the current definition of cellular plasticity in the field of mammary gland biology. Following a synopsis of the evolving model of the hierarchical development of mammary epithelial cell lineages, we discuss changes in cell identity during normal mammary gland development with particular emphasis on the effect of the gestation cycle on the emergence of new cellular states. Next, we summarize known mechanisms that promote the plasticity of epithelial lineages in the normal mammary gland and highlight the importance of the microenvironment and extracellular matrix. A discourse of cellular reprogramming during the early stages of mammary tumorigenesis that follows focuses on the origin of basal-like breast cancers from luminal progenitors and oncogenic signaling networks that orchestrate diverse developmental trajectories of transforming epithelial cells. In addition to the epithelial-to-mesenchymal transition, we highlight events of cellular reprogramming during breast cancer progression in the context of intrinsic molecular subtype switching and the genesis of the claudin-low breast cancer subtype, which represents the far end of the spectrum of epithelial cell plasticity. In the final section, we will discuss recent advances in the design of genetically engineered models to gain insight into the dynamic processes that promote cellular plasticity during mammary gland development and tumorigenesis in vivo.

A Cellular atlas of the human fallopian tube reveals the metamorphosis of secretory epithelial cells during the menstrual cycle and menopause

The fallopian tube, connecting the uterus with the ovary, is a dynamic organ that undergoes cyclical changes and is the site of several diseases, including serous cancer. Here, we use single-cell technologies to construct a comprehensive cell map of healthy pre-menopausal fallopian tubes, capturing the impact of the menstrual cycle and menopause on different fallopian tube cells at the molecular level. The comparative analysis between pre- and post-menopausal fallopian tubes reveals substantial shifts in cellular abundance and gene expression patterns, highlighting the physiological changes associated with menopause. Further investigations into menstrual cycle phases illuminate distinct molecular states in secretory epithelial cells caused by hormonal fluctuations. The markers we identified characterizing secretory epithelial cells provide a valuable tool for classifying ovarian cancer subtypes.

Graphical summary

Graphical summary of results. During the proliferative phase (estrogen high ) of the menstrual cycle, SE2 cells (OVGP1 + ) dominate the fallopian tube (FT) epithelium, while SE1 cells (OVGP1 - ) dominate the epithelium during the secretory phase. Though estrogen levels decrease during menopause, SE post-cells (OVGP1 + , CXCL2 + ) make up most of the FT epithelium.

FXYD3 functionally demarcates an ancestral breast cancer stem cell subpopulation with features of drug-tolerant persisters

The heterogeneity of cancer stem cells (CSCs) within tumors presents a challenge in therapeutic targeting. To decipher the cellular plasticity that fuels phenotypic heterogeneity, we undertook single-cell transcriptomics analysis in triple-negative breast cancer (TNBC) to identify subpopulations in CSCs. We found a subpopulation of CSCs with ancestral features that is marked by FXYD domain-containing ion transport regulator 3 (FXYD3), a component of the Na+/K+ pump. Accordingly, FXYD3+ CSCs evolve and proliferate, while displaying traits of alveolar progenitors that are normally induced during pregnancy. Clinically, FXYD3+ CSCs were persistent during neoadjuvant chemotherapy, hence linking them to drug-tolerant persisters (DTPs) and identifying them as crucial therapeutic targets. Importantly, FXYD3+ CSCs were sensitive to senolytic Na+/K+ pump inhibitors, such as cardiac glycosides. Together, our data indicate that FXYD3+ CSCs with ancestral features are drivers of plasticity and chemoresistance in TNBC. Targeting the Na+/K+ pump could be an effective strategy to eliminate CSCs with ancestral and DTP features that could improve TNBC prognosis.

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.

Advancements in Human Breast Organoid Culture: Modeling Complex Tissue Structures and Developmental Insights

Organoids have been widely used for studying tissue growth and modeling diseases, but achieving physiologically relevant architecture, size, and function has remained a challenge. Here, we develop a next-generation organotypic culture method that enables the formation of a highly patterned, complex, branched tissue that is spatially organized to accurately recapitulate the morphology, scale, cellular, transcriptional, and tissue-level heterogeneity of human breast tissue. Hormone responsiveness of organoids is also a feature allowing for examination of androgen therapy or post-menopausal changes to breast tissue development and regeneration. Live imaging allows for studying stem cell dynamics during organoid formation and is adaptable to a high throughput setting. Real-time imaging of organoid formation reveals activation of latent epithelial organogenesis programs and inductive cellular dynamics that drive formation of a miniature breast tissue along with its mesenchyme akin to tissue stroma. By advancing human breast organoid technology, this model can elucidate cell- and tissue-level consequences to hormonal changes and therapy. In addition, this method can lead to new insights into the cellular, molecular, and tissue-level processes involved in organogenesis and regeneration, as well as disease.

Single-cell transcription mapping of murine and human mammary organoids responses to female hormones

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

Stromal heterogeneity may explain increased incidence of metaplastic breast cancer in women of African descent

The biologic basis of genetic ancestry-dependent variability in disease incidence and outcome is just beginning to be explored. We recently reported enrichment of a population of ZEB1-expressing cells located adjacent to ductal epithelial cells in normal breasts of women of African ancestry compared to those of European ancestry. In this study, we demonstrate that these cells have properties of fibroadipogenic/mesenchymal stromal cells that express PROCR and PDGFRα and transdifferentiate into adipogenic and osteogenic lineages. PROCR + /ZEB1 + /PDGFRα+ (PZP) cells are enriched in normal breast tissues of women of African compared to European ancestry. PZP: epithelial cell communication results in luminal epithelial cells acquiring basal cell characteristics and IL-6-dependent increase in STAT3 phosphorylation. Furthermore, level of phospho-STAT3 is higher in normal and cancerous breast tissues of women of African ancestry. PZP cells transformed with HRasG12V ± SV40-T/t antigens generate metaplastic carcinoma suggesting that these cells are one of the cells-of-origin of metaplastic breast cancers.

A transcriptional response to replication stress selectively expands a subset of Brca2-mutant mammary epithelial cells

Germline BRCA2 mutation carriers frequently develop luminal-like breast cancers, but it remains unclear how BRCA2 mutations affect mammary epithelial subpopulations. Here, we report that monoallelic Brca2mut/WT mammary organoids subjected to replication stress activate a transcriptional response that selectively expands Brca2mut/WT luminal cells lacking hormone receptor expression (HR-). While CyTOF analyses reveal comparable epithelial compositions among wildtype and Brca2mut/WT mammary glands, Brca2mut/WT HR- luminal cells exhibit greater organoid formation and preferentially survive and expand under replication stress. ScRNA-seq analysis corroborates the expansion of HR- luminal cells which express elevated transcript levels of Tetraspanin-8 (Tspan8) and Thrsp, plus pathways implicated in replication stress survival including Type I interferon responses. Notably, CRISPR/Cas9-mediated deletion of Tspan8 or Thrsp prevents Brca2mut/WT HR- luminal cell expansion. Our findings indicate that Brca2mut/WT cells activate a transcriptional response after replication stress that preferentially favours outgrowth of HR- luminal cells through the expression of interferon-responsive and mammary alveolar genes.

The EstroGene Database Reveals Diverse Temporal, Context-Dependent, and Bidirectional Estrogen Receptor Regulomes in Breast Cancer

As one of the most successful cancer therapeutic targets, estrogen receptor-α (ER/ESR1) has been extensively studied over the past few decades. Sequencing technological advances have enabled genome-wide analysis of ER action. However, comparison of individual studies is limited by different experimental designs, and few meta-analyses are available. Here, we established the EstroGene database through unified processing of data from 246 experiments including 136 transcriptomic, cistromic, and epigenetic datasets focusing on estradiol (E2)-triggered ER activation across 19 breast cancer cell lines. A user-friendly browser (https://estrogene.org/) was generated for multiomic data visualization involving gene inquiry under user-defined experimental conditions and statistical thresholds. Notably, annotation of metadata associated with public datasets revealed a considerable lack of experimental details. Comparison of independent RNA-seq or ER ChIP-seq data with the same design showed large variability and only strong effects could be consistently detected. Temporal estrogen response metasignatures were defined, and the association of E2 response rate with temporal transcriptional factors, chromatin accessibility, and heterogeneity of ER expression was evaluated. Unexpectedly, harmonizing 146 E2-induced transcriptomic datasets uncovered a subset of genes harboring bidirectional E2 regulation, which was linked to unique transcriptional factors and highly associated with immune surveillance in the clinical setting. Furthermore, the context dependent E2 response programs were characterized in MCF7 and T47D cell lines, the two most frequently used models in the EstroGene database. Collectively, the EstroGene database provides an informative and practical resource to the cancer research community to uniformly evaluate key reproducible features of ER regulomes and unravels modes of ER signaling.

SIGNIFICANCE

A resource database integrating 246 publicly available ER profiling datasets facilitates meta-analyses and identifies estrogen response temporal signatures, a bidirectional program, and model-specific biases.

A spatially resolved single-cell genomic atlas of the adult human breast

The adult human breast is comprised of an intricate network of epithelial ducts and lobules that are embedded in connective and adipose tissue1-3. Although most previous studies have focused on the breast epithelial system4-6, many of the non-epithelial cell types remain understudied. Here we constructed the comprehensive Human Breast Cell Atlas (HBCA) at single-cell and spatial resolution. Our single-cell transcriptomics study profiled 714,331 cells from 126 women, and 117,346 nuclei from 20 women, identifying 12 major cell types and 58 biological cell states. These data reveal abundant perivascular, endothelial and immune cell populations, and highly diverse luminal epithelial cell states. Spatial mapping using four different technologies revealed an unexpectedly rich ecosystem of tissue-resident immune cells, as well as distinct molecular differences between ductal and lobular regions. Collectively, these data provide a reference of the adult normal breast tissue for studying mammary biology and diseases such as breast cancer.

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.

Label-free single cell proteomics utilizing ultrafast LC and MS instrumentation: A valuable complementary technique to multiplexing

The ability to map a proteomic fingerprint to transcriptomic data would master the understanding of how gene expression translates into actual phenotype. In contrast to nucleic acid sequencing, in vitro protein amplification is impossible and no single cell proteomic workflow has been established as gold standard yet. Advances in microfluidic sample preparation, multi-dimensional sample separation, sophisticated data acquisition strategies, and intelligent data analysis algorithms have resulted in major improvements to successfully analyze such tiny sample amounts with steadily boosted performance. However, among the broad variation of published approaches, it is commonly accepted that highest possible sensitivity, robustness, and throughput are still the most urgent needs for the field. While many labs have focused on multiplexing to achieve these goals, label-free SCP is a highly promising strategy as well whenever high dynamic range and unbiased accurate quantification are needed. We here focus on recent advances in label-free single-cell mass spectrometry workflows and try to guide our readers to choose the best method or combinations of methods for their specific applications. We further highlight which techniques are most propitious in the future and which applications but also limitations we foresee for the field.

Blood Supplementation Enhances Bartonella henselae Growth and Molecular Detection of Bacterial DNA in Liquid Culture

Bacteria of the genus Bartonella, a member of the Alphaproteobacteria, are fastidious, Gram-negative, aerobic bacilli that comprise numerous species, subspecies, and genotypes. Bartonella henselae, with a worldwide distribution, infects cats, dogs, horses, humans, and other mammals. Diagnostically, direct detection of Bartonella henselae in patient blood specimens by culture or molecular methods is required to confirm infection with this bacterium. Enrichment blood culture combined with quantitative PCR (qPCR) or ddPCR enhances the sensitivity of direct detection. The addition of sheep blood to liquid culture media increased the Bartonella henselae DNA concentration compared to controls, additionally improving PCR direct detection sensitivity. IMPORTANCE This study aims to improve diagnostic detection of Bartonella henselae. Patient samples are combined with enriched bacterial cultures aimed at growing Bartonella henselae for the best possible chance at detection. However, current Bartonella growth methods could be improved. The DNA extraction method used by most laboratories should also be optimized. Sheep blood was added to increase the growth of Bartonella henselae and multiple DNA extraction methods were to be compared to each other.

Spatio-temporal transcriptome dynamics coordinate rapid transition of core crop functions in 'lactating' pigeon

Pigeons (Columba livia) are among a select few avian species that have developed a specialized reproductive mode wherein the parents produce a 'milk' in their crop to feed newborn squabs. Nonetheless, the transcriptomic dynamics and role in the rapid transition of core crop functions during 'lactation' remain largely unexplored. Here, we generated a de novo pigeon genome assembly to construct a high resolution spatio-temporal transcriptomic landscape of the crop epithelium across the entire breeding stage. This multi-omics analysis identified a set of 'lactation'-related genes involved in lipid and protein metabolism, which contribute to the rapid functional transitions in the crop. Analysis of in situ high-throughput chromatin conformation capture (Hi-C) sequencing revealed extensive reorganization of promoter-enhancer interactions linked to the dynamic expression of these 'lactation'-related genes between stages. Moreover, their expression is spatially localized in specific epithelial layers, and can be correlated with phenotypic changes in the crop. These results illustrate the preferential de novo synthesis of 'milk' lipids and proteins in the crop, and provides candidate enhancer loci for further investigation of the regulatory elements controlling pigeon 'lactation'.

Changes in the mammary gland during aging and its links with breast diseases

The functional capacity of organisms declines in the process of aging. In the case of breast tissue, abnormal mammary gland development can lead to dysfunction in milk secretion, a primary function, as well as the onset of various diseases, such as breast cancer. In the process of aging, the terminal duct lobular units (TDLUs) within the breast undergo gradual degeneration, while the proportion of adipose tissue in the breast continues to increase and hormonal levels in the breast change accordingly. Here, we review changes in morphology, internal structure, and cellular composition that occur in the mammary gland during aging. We also explore the emerging mechanisms of breast aging and the relationship between changes during aging and breast-related diseases, as well as potential interventions for delaying mammary gland aging and preventing breast disease.

A spatially resolved single cell genomic atlas of the adult human breast

The adult human breast comprises an intricate network of epithelial ducts and lobules that are embedded in connective and adipose tissue. While previous studies have mainly focused on the breast epithelial system, many of the non-epithelial cell types remain understudied. Here, we constructed a comprehensive Human Breast Cell Atlas (HBCA) at single-cell and spatial resolution. Our single-cell transcriptomics data profiled 535,941 cells from 62 women, and 120,024 nuclei from 20 women, identifying 11 major cell types and 53 cell states. These data revealed abundant pericyte, endothelial and immune cell populations, and highly diverse luminal epithelial cell states. Our spatial mapping using three technologies revealed an unexpectedly rich ecosystem of tissue-resident immune cells in the ducts and lobules, as well as distinct molecular differences between ductal and lobular regions. Collectively, these data provide an unprecedented reference of adult normal breast tissue for studying mammary biology and disease states such as breast cancer.

Preneoplastic stromal cells promote BRCA1-mediated breast tumorigenesis

Women with germline BRCA1 mutations (BRCA1+/mut) have increased risk for hereditary breast cancer. Cancer initiation in BRCA1+/mut is associated with premalignant changes in breast epithelium; however, the role of the epithelium-associated stromal niche during BRCA1-driven tumor initiation remains unclear. Here we show that the premalignant stromal niche promotes epithelial proliferation and mutant BRCA1-driven tumorigenesis in trans. Using single-cell RNA sequencing analysis of human preneoplastic BRCA1+/mut and noncarrier breast tissues, we show distinct changes in epithelial homeostasis including increased proliferation and expansion of basal-luminal intermediate progenitor cells. Additionally, BRCA1+/mut stromal cells show increased expression of pro-proliferative paracrine signals. In particular, we identify pre-cancer-associated fibroblasts (pre-CAFs) that produce protumorigenic factors including matrix metalloproteinase 3 (MMP3), which promotes BRCA1-driven tumorigenesis in vivo. Together, our findings demonstrate that precancerous stroma in BRCA1+/mut may elevate breast cancer risk through the promotion of epithelial proliferation and an accumulation of luminal progenitor cells with altered differentiation.

Implantation of engineered adipocytes that outcompete tumors for resources suppresses cancer progression

Tumors acquire an increased ability to obtain and metabolize nutrients. Here, we engineered and implanted adipocytes to outcompete tumors for nutrients and show that they can substantially reduce cancer progression. Growing cells or xenografts from several cancers (breast, colon, pancreas, prostate) alongside engineered human adipocytes or adipose organoids significantly suppresses cancer progression and reduces hypoxia and angiogenesis. Transplanting modulated adipocyte organoids in pancreatic or breast cancer mouse models nearby or distal from the tumor significantly suppresses its growth. To further showcase therapeutic potential, we demonstrate that co-culturing tumor organoids derived from human breast cancers with engineered patient-derived adipocytes significantly reduces cancer growth. Combined, our results introduce a novel cancer therapeutic approach, termed adipose modulation transplantation (AMT), that can be utilized for a broad range of cancers.

Identification and characterization of a proliferative cell population in estrogen receptor-positive metastatic breast cancer through spatial and single-cell transcriptomics

Background

Intratumor heterogeneity is a hallmark of most solid tumors, including breast cancers. We applied spatial transcriptomics and single-cell RNA-sequencing technologies to profile spatially resolved cell populations within estrogen receptor-positive (ER ⁺ ) metastatic breast cancers and elucidate their importance in estrogen-dependent tumor growth.

Methods

Spatial transcriptomics and single-cell RNA-sequencing were performed on two patient-derived xenografts (PDXs) of "ER-high" metastatic breast cancers with opposite estrogen-mediated growth responses: estrogen-suppressed GS3 (80-100% ER) and estrogen-stimulated SC31 (30-75% ER) models. The analyses included samples treated with and without 17β-estradiol. The findings were validated via scRNA-seq analyses on "ER-low" estrogen-accelerating PDX, GS1 (5% ER). The results from our spatial and single-cell analyses were further supported by the analysis of a publicly available single cell dataset and a protein-based dual immunohistochemical (IHC) evaluation using three important clinical markers [i.e., ER, progesterone receptor (PR), and Ki67]. The translational implication of these results was assessed by clinical outcome analyses on public breast cancer cohorts.

Results

Our novel space-gene-function study revealed a "proliferative" cell population in addition to three major spatially distinct compartments within ER ⁺ metastatic breast cancers. These compartments showed functional diversity (i.e., estrogen-responsive, proliferative, hypoxia-induced, and inflammation-related). The "proliferative ( MKI67 ⁺ )" population, not "estrogen-responsive" compartment, was crucial for estrogen-dependent tumor growth, leading to the acquisition of luminal B features. The cells with induction of typical estrogen-responsive genes such as PGR were not directly linked to estrogen-dependent proliferation. Additionally, the dual IHC analyses demonstrated the distinct contribution of the Ki67 ⁺ proliferative cells toward estrogen-mediated growth and their response to palbociclib, a CDK4/6 inhibitor. The gene signatures developed from the proliferative, hypoxia-induced, and inflammation-related compartments were significantly correlated with worse clinical outcomes, while patients with the high estrogen-responsive scores showed better prognosis, confirming that the estrogen-responsive compartment would not be directly associated with estrogen-dependent tumor progression.

Conclusions

For the first time, our study elucidated a "proliferative" cell population distinctly distributed in ER ⁺ metastatic breast cancers. They contribute differently toward progression of these cancers, and the gene signature in the "proliferative" compartment is an important determinant of luminal cancer subtypes.

Cellular Plasticity and Heterotypic Interactions during Breast Morphogenesis and Cancer Initiation

Simple Summary

This review aims to discuss the structure, function and dynamics of the breast gland and how changes to the function of the breast’s cells can lead to different types of cancer.

Abstract

The human breast gland is a unique organ as most of its development occurs postnatally between menarche and menopause, a period ranging from 30 to 40 years. During this period, the monthly menstruation cycle drives the mammary gland through phases of cell proliferation, differentiation, and apoptosis, facilitated via a closely choreographed interaction between the epithelial cells and the surrounding stroma preparing the gland for pregnancy. If pregnancy occurs, maximal differentiation is reached to prepare for lactation. After lactation, the mammary gland involutes to a pre-pregnant state. These cycles of proliferation, differentiation, and involution necessitate the presence of epithelial stem cells that give rise to progenitor cells which differentiate further into the luminal and myoepithelial lineages that constitute the epithelial compartment and are responsible for the branching structure of the gland. Maintaining homeostasis and the stem cell niche depends strongly on signaling between the stem and progenitor cells and the surrounding stroma. Breast cancer is a slowly progressing disease whose initiation can take decades to progress into an invasive form. Accumulating evidence indicates that stem cells and/or progenitor cells at different stages, rather than terminally differentiated cells are the main cells of origin for most breast cancer subgroups. Stem cells and cancer cells share several similarities such as increased survival and cellular plasticity which is reflected in their ability to switch fate by receiving intrinsic and extrinsic signals. In this review, we discuss the concept of cellular plasticity in normal breast morphogenesis and cancer, and how the stromal environment plays a vital role in cancer initiation and progression.

Cellular and molecular mechanisms of breast cancer susceptibility

There is a plethora of recognized risk factors for breast cancer (BC) with poorly understood or speculative biological mechanisms. The lack of prevention options highlights the importance of understanding the mechanistic basis of cancer susceptibility and finding new targets for breast cancer prevention. Until now, we have understood risk and cancer susceptibility primarily through the application of epidemiology and assessing outcomes in large human cohorts. Relative risks are assigned to various human behaviors and conditions, but in general the associations are weak and there is little understanding of mechanism. Aging is by far the greatest risk factor for BC, and there are specific forms of inherited genetic risk that are well-understood to cause BC. We propose that bringing focus to the biology underlying these forms of risk will illuminate biological mechanisms of BC susceptibility.