A human breast atlas integrating single-cell proteomics and transcriptomics

Knock-out of CD73 delays the onset of HR-negative breast cancer by reprogramming lipid metabolism and is associated with increased tumor mutational burden

OBJECTIVE

CD73 (ecto-5'-nucleotidase, NT5E), a cell-surface enzyme converting 5'-AMP to adenosine, is crucial for cancer progression. However, its role in the tumorigenesis process remains mostly obscure. We aimed to demonstrate CD73's role in breast cancer (BC) tumorigenesis through metabolic rewiring of fatty acid metabolism, a process recently indicated to be regulated by BC major prognostic markers, hormone receptors (HR) for estrogen (ER), and progesterone (PR).

METHODS

A murine model of chemically induced mammary gland tumorigenesis was applied to analyze CD73 knock-out (KO)-induced changes at the transcriptome (RNA-seq), proteome (IHC, WB), and lipidome (GC-EI-MS) levels. CD73 KO-induced changes were correlated with scRNA-seq and bulk RNA-seq data for human breast tissues and BCs from public collections and confirmed at the proteome level with IHC or WB analysis of BC tissue microarrays and cell lines.

RESULTS

CD73 KO delayed the onset of HR/PR-negative mammary tumors in a murine model. This delay correlated with increased expression of genes related to biosynthesis and β-oxidation of fatty acids (FAs) in the CD73 KO group at the initiation stage. STRING analysis based on RNA-seq data indicated an interplay between CD73 KO, up-regulated expression of PR-coding gene, and DEGs involved in FA metabolism, with PPARγ, a main regulator of FA synthesis, as a main connective node. In epithelial cells of mammary glands, PPARγ expression correlated with CD73 at the RNA level. With cancer progression, CD73 KO increased the levels of PUFAn3/6 (polyunsaturated omega 3/6 FAs), known ligands of PPARγ and target for lipid peroxidation, which may lead to oxidative DNA damage. It correlated with the downregulation of genes involved in cellular stress response (Mlh1, Gsta3), PR-or CD73-dependent changes in the intracellular ROS levels and expression or activation of proteins involved in DNA repair or oxidative stress response in mammary tumor or human BC cell lines, increased tumor mutational burden (TMB) and genomic instability markers in CD73 low HR-negative human BCs, and the prolonged onset of tumors in the CD73 KO HR/PR-negative group.

CONCLUSIONS

CD73 has a significant role in tumorigenesis driving the reprogramming of lipid metabolism through the regulatory loop with PR and PPARγ in epithelial cells of mammary glands. Low CD73 expression/CD73 KO might enhance mutational burden by disrupting this regulatory loop, delaying the onset of HR-negative tumors. Our results support combining therapy targeting the CD73-adenosine axis and tumor lipidome against HR-negative tumors, especially at their earliest developmental stage.

Volumetric analysis of the terminal ductal lobular unit architecture and cell phenotypes in the human breast

The major lactiferous ducts of the human breast branch out and end at terminal ductal lobular units (TDLUs). Despite their functional and clinical importance, the three-dimensional (3D) architecture of TDLUs has remained undetermined. Our quantitative and volumetric imaging of healthy human breast tissue demonstrates that highly branched TDLUs, which exhibit increased proliferation, are uncommon in the resting tissue regardless of donor age, parity, or hormonal contraception. Overall, TDLUs have a consistent shape and branch parameters, and they contain a main subtree that dominates in bifurcation events and exhibits a more duct-like keratin expression pattern. Simulation of TDLU branching morphogenesis in three dimensions suggests that evolutionarily conserved mechanisms regulate mammary gland branching in humans and mice despite their anatomical differences. In all, our data provide structural insight into 3D anatomy and branching of the human breast and exemplify the power of volumetric imaging in gaining a deeper understanding of breast biology.

Loss of STING impairs lactogenic differentiation

Heightened energetic and nutrient demand during lactogenic differentiation of the mammary gland elicits upregulation of various stress responses to support cellular homeostasis. Here, we identify the stimulator of interferon genes (STING) as an immune supporter of the functional development of mouse mammary epithelial cells (MECs). An in vitro model of MEC differentiation revealed that STING is activated in a cGAS-independent manner to produce both type I interferons and proinflammatory cytokines in response to the accumulation of mitochondrial reactive oxygen species. Induction of STING activity was found to be dependent on the breast tumor suppressor gene single-minded 2 (SIM2). Using mouse models of lactation, we discovered that loss of STING activity results in early involution of #3 mammary glands, severely impairing lactational performance. Our data suggest that STING is required for successful functional differentiation of the mammary gland and bestows a differential lactogenic phenotype between #3 mammary glands and the traditionally explored inguinal 4|9 pair. These findings affirm unique development of mammary gland pairs that is essential to consider in future investigations into normal development and breast cancer initiation.

Tumor editing suppresses innate and adaptive antitumor immunity and is reversed by inhibiting DNA methylation

Cancer cells edit gene expression to evade immunosurveillance. However, genome-wide studies of gene editing during early tumorigenesis are lacking. Here we used single-cell RNA sequencing in a breast cancer genetically engineered mouse model (GEMM) to identify edited genes without bias. Late tumors repressed antitumor immunity genes, reducing infiltrating immune cells and tumor-immune cell communications. Innate immune genes, especially interferon-stimulated genes, dominated the list of downregulated tumor genes, while genes that regulate cell-intrinsic malignancy were mostly unedited. Naive and activated CD8+ T cells in early tumors were replaced with exhausted or precursor-exhausted cells in late tumors. Repression of immune genes was reversed by inhibiting DNA methylation using low-dose decitabine, which suppressed tumor growth and restored immune control, increasing the number, functionality and memory of tumor-infiltrating lymphocytes and reducing the number of myeloid suppressor cells. Decitabine induced important interferon, pyroptosis and necroptosis genes, inflammatory cell death and immune control in GEMM and implanted breast and melanoma tumors.

Midkine as a driver of age-related changes and increase in mammary tumorigenesis

Aging is a pivotal risk factor for cancer, yet the underlying mechanisms remain poorly defined. Here, we explore age-related changes in the rat mammary gland by single-cell multiomics. Our findings include increased epithelial proliferation, loss of luminal identity, and decreased naive B and T cells with age. We discover a luminal progenitor population unique to old rats with profiles reflecting precancerous changes and identify midkine (Mdk) as a gene upregulated with age and a regulator of age-related luminal progenitors. Midkine treatment of young rats mimics age-related changes via activating PI3K-AKT-SREBF1 pathway and promotes nitroso-N-methylurea-induced mammary tumorigenesis. Midkine levels increase with age in human blood and mammary epithelium, and higher MDK in normal breast tissue is associated with higher breast cancer risk in younger women. Our findings reveal a link between aging and susceptibility to tumor initiation and identify midkine as a mediator of age-dependent increase in breast tumorigenesis.

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.

Mcam stabilizes a luminal progenitor-like breast cancer cell state via Ck2 control and Src/Akt/Stat3 attenuation

Cell state control is crucial for normal tissue development and cancer cell mimicry of stem/progenitor states, contributing to tumor heterogeneity, therapy resistance, and progression. Here, we demonstrate that the cell surface glycoprotein Mcam maintains the tumorigenic luminal progenitor (LP)-like epithelial cell state, leading to Basal-like mammary cancers. In the Py230 mouse mammary carcinoma model, Mcam knockdown (KD) destabilized the LP state by deregulating the Ck2/Stat3 axis, causing a switch to alveolar and basal states, loss of an estrogen-sensing subpopulation, and resistance to tamoxifen-an effect reversed by Ck2 and Stat3 inhibitors. In vivo, Mcam KD blocked generation of Basal-like tumors and Sox10+Krt14+ cells. In human tumors, MCAM loss was largely exclusive of the Basal-like subtype, linked instead to proliferative Luminal subtypes, including often endocrine-resistant Luminal B cancers. This study has implications for developing therapies targeting MCAM, CK2, and STAT3 and their likely effective contexts.

A triple hormone receptor ER, AR, and VDR signature is a robust prognosis predictor in breast cancer

BACKGROUND

Despite evidence indicating the dominance of cell-of-origin signatures in molecular tumor patterns, translating these genome-wide patterns into actionable insights has been challenging. This study introduces breast cancer cell-of-origin signatures that offer significant prognostic value across all breast cancer subtypes and various clinical cohorts, compared to previously developed genomic signatures.

METHODS

We previously reported that triple hormone receptor (THR) co-expression patterns of androgen (AR), estrogen (ER), and vitamin D (VDR) receptors are maintained at the protein level in human breast cancers. Here, we developed corresponding mRNA signatures (THR-50 and THR-70) based on these patterns to categorize breast tumors by their THR expression levels. The THR mRNA signatures were evaluated across 56 breast cancer datasets (5040 patients) using Kaplan-Meier survival analysis, Cox proportional hazard regression, and unsupervised clustering.

RESULTS

The THR signatures effectively predict both overall and progression-free survival across all evaluated datasets, independent of subtype, grade, or treatment status, suggesting improvement over existing prognostic signatures. Furthermore, they delineate three distinct ER-positive breast cancer subtypes with significant survival in differences-expanding on the conventional two subtypes. Additionally, coupling THR-70 with an immune signature identifies a predominantly ER-negative breast cancer subgroup with a highly favorable prognosis, comparable to ER-positive cases, as well as an ER-negative subgroup with notably poor outcome, characterized by a 15-fold shorter survival.

CONCLUSIONS

The THR cell-of-origin signature introduces a novel dimension to breast cancer biology, potentially serving as a robust foundation for integrating additional prognostic biomarkers. These signatures offer utility as a prognostic index for stratifying existing breast cancer subtypes and for de novo classification of breast cancer cases. Moreover, THR signatures may also hold promise in predicting hormone treatment responses targeting AR and/or VDR.

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.

Single-Cell RNA-Sequencing: Opening New Horizons for Breast Cancer Research

Breast cancer is the most prevalent malignant tumor among women with high heterogeneity. Traditional techniques frequently struggle to comprehensively capture the intricacy and variety of cellular states and interactions within breast cancer. As global precision medicine rapidly advances, single-cell RNA sequencing (scRNA-seq) has become a highly effective technique, revolutionizing breast cancer research by offering unprecedented insights into the cellular heterogeneity and complexity of breast cancer. This cutting-edge technology facilitates the analysis of gene expression profiles at the single-cell level, uncovering diverse cell types and states within the tumor microenvironment. By dissecting the cellular composition and transcriptional signatures of breast cancer cells, scRNA-seq provides new perspectives for understanding the mechanisms behind tumor therapy, drug resistance and metastasis in breast cancer. In this review, we summarized the working principle and workflow of scRNA-seq and emphasized the major applications and discoveries of scRNA-seq in breast cancer research, highlighting its impact on our comprehension of breast cancer biology and its potential for guiding personalized treatment strategies.

Single-nucleus chromatin accessibility and transcriptomic map of breast tissues of women of diverse genetic ancestry

Single-nucleus analysis allows robust cell-type classification and helps to establish relationships between chromatin accessibility and cell-type-specific gene expression. Here, using samples from 92 women of several genetic ancestries, we developed a comprehensive chromatin accessibility and gene expression atlas of the breast tissue. Integrated analysis revealed ten distinct cell types, including three major epithelial subtypes (luminal hormone sensing, luminal adaptive secretory precursor (LASP) and basal-myoepithelial), two endothelial and adipocyte subtypes, fibroblasts, T cells, and macrophages. In addition to the known cell identity genes FOXA1 (luminal hormone sensing), EHF and ELF5 (LASP), TP63 and KRT14 (basal-myoepithelial), epithelial subtypes displayed several uncharacterized markers and inferred gene regulatory networks. By integrating breast epithelial cell gene expression signatures with spatial transcriptomics, we identified gene expression and signaling differences between lobular and ductal epithelial cells and age-associated changes in signaling networks. LASP cells and fibroblasts showed genetic ancestry-dependent variability. An estrogen receptor-positive subpopulation of LASP cells with alveolar progenitor cell state was enriched in women of Indigenous American ancestry. Fibroblasts from breast tissues of women of African and European ancestry clustered differently, with accompanying gene expression differences. Collectively, these data provide a vital resource for further exploring genetic ancestry-dependent variability in healthy breast biology.

Computational Analyses Reveal Deregulated Clock Genes Associated with Breast Cancer Development in Night Shift Workers

Breast cancer (BC) is the leading cause of cancer death among women worldwide. Women employed in shift jobs face heightened BC risk due to prolonged exposure to night shift work (NSW), classified as potentially carcinogenic by the International Agency for Research on Cancer (IARC). This risk is linked to disruptions in circadian rhythms governed by clock genes at the cellular level. However, the molecular mechanisms are unclear. This study aimed to assess clock genes as potential BC biomarkers among women exposed to long-term NSW. Clock gene expression was analysed in paired BC and normal breast tissues within Nurses' Health Studies I and II GEO datasets. Validation was performed on additional gene expression datasets from healthy night shift workers and women with varying BC susceptibility, as well as single-cell sequencing datasets. Post-transcriptional regulators of clock genes were identified through miRNA analyses. Significant alterations in clock gene expression in BC compared to normal tissues were found. BHLHE40, CIART, CLOCK, PDPK1, and TIMELESS were over-expressed, while HLF, NFIL3, NPAS3, PER1, PER3, SIM1, and TEF were under-expressed. The downregulation of PER1 and TEF and upregulation of CLOCK correlated with increased BC risk in healthy women. Also, twenty-six miRNAs, including miR-10a, miR-21, miR-107, and miR-34, were identified as potential post-transcriptional regulators influenced by NSW. In conclusion, a panel of clock genes and circadian miRNAs are suggested as BC susceptibility biomarkers among night shift workers, supporting implications for risk stratification and early detection strategies.

Dissection of triple-negative breast cancer microenvironment and identification of potential therapeutic drugs using single-cell RNA sequencing analysis

Breast cancer remains a leading cause of mortality in women worldwide. Triple-negative breast cancer (TNBC) is a particularly aggressive subtype characterized by rapid progression, poor prognosis, and lack of clear therapeutic targets. In the clinic, delineation of tumor heterogeneity and development of effective drugs continue to pose considerable challenges. Within the scope of our study, high heterogeneity inherent to breast cancer was uncovered based on the landscape constructed from both tumor and healthy breast tissue samples. Notably, TNBC exhibited significant specificity regarding cell proliferation, differentiation, and disease progression. Significant associations between tumor grade, prognosis, and TNBC oncogenes were established via pseudotime trajectory analysis. Consequently, we further performed comprehensive characterization of the TNBC microenvironment. A crucial epithelial subcluster, E8, was identified as highly malignant and strongly associated with tumor cell proliferation in TNBC. Additionally, epithelial-mesenchymal transition (EMT)-associated fibroblast and M2 macrophage subclusters exerted an influence on E8 through cellular interactions, contributing to tumor growth. Characteristic genes in these three cluster cells could therefore serve as potential therapeutic targets for TNBC. The collective findings provided valuable insights that assisted in the screening of a series of therapeutic drugs, such as pelitinib. We further confirmed the anti-cancer effect of pelitinib in an orthotopic 4T1 tumor-bearing mouse model. Overall, our study sheds light on the unique characteristics of TNBC at single-cell resolution and the crucial cell types associated with tumor cell proliferation that may serve as potent tools in the development of effective anti-cancer drugs.

Spatial molecular profiling of mixed invasive ductal and lobular breast cancers reveals heterogeneity in intrinsic molecular subtypes, oncogenic signatures, and mutations

Mixed invasive ductal and lobular carcinoma (MDLC) is a rare histologic subtype of breast cancer displaying both E-cadherin positive ductal and E-cadherin negative lobular morphologies within the same tumor, posing challenges with regard to anticipated clinical management. It remains unclear whether these distinct morphologies also have distinct biology and risk of recurrence. Our spatially resolved transcriptomic, genomic, and single-cell profiling revealed clinically significant differences between ductal and lobular tumor regions including distinct intrinsic subtype heterogeneity - e.g., MDLC with triple-negative breast cancer (TNBC) or basal ductal and estrogen receptor positive (ER+) luminal lobular regions, distinct enrichment of cell cycle arrest/senescence and oncogenic (ER and MYC) signatures, genetic and epigenetic CDH1 inactivation in lobular but not ductal regions, and single-cell ductal and lobular subpopulations with unique oncogenic signatures further highlighting intraregional heterogeneity. Altogether, we demonstrated that the intratumoral morphological/histological heterogeneity within MDLC is underpinned by intrinsic subtype and oncogenic heterogeneity which may result in prognostic uncertainty and therapeutic dilemma.

Organoid Cultures for the Study of Mammary Biology and Breast Cancer: The Promise and Challenges

During the last decade, biomedical research has experienced a resurgence in the use of three-dimensional culture models for studies of normal and cancer biology. This resurgence has been driven by the development of models in which primary cells are grown in tissue-mimicking media and extracellular matrices to create organoid or organotypic cultures that more faithfully replicate the complex architecture and physiology of normal tissues and tumors. In addition, patient-derived tumor organoids preserve the three-dimensional organization and characteristics of the patient tumors ex vivo, becoming excellent preclinical models to supplement studies of tumor xenografts transplanted into immunocompromised mice. In this perspective, we provide an overview of how organoids are being used to investigate normal mammary biology and as preclinical models of breast cancer and discuss improvements that would enhance their utility and relevance to the field.

Breast organoid suspension cultures maintain long-term estrogen receptor expression and responsiveness

Organoid cultures offer a powerful technology to investigate many different aspects of development, physiology, and pathology of diverse tissues. Unlike standard tissue culture of primary breast epithelial cells, breast organoids preserve the epithelial lineages and architecture of the normal tissue. However, existing organoid culture methods are tedious, difficult to scale, and do not robustly retain estrogen receptor (ER) expression and responsiveness in long-term culture. Here, we describe a modified culture method to generate and maintain organoids as suspension cultures in reconstituted basement membrane (™Matrigel). This method improves organoid growth and uniformity compared to the conventional Matrigel dome embedding method, while maintaining the fidelity of the three major epithelial lineages. Using this adopted method, we are able to culture and passage purified hormone sensing (HS) cells that retain ER responsiveness upon estrogen stimulation in long-term culture. This culture system presents a valuable platform to study the events involved in initiation and evolution of ER-positive breast cancer.

Spatial molecular profiling of mixed invasive ductal-lobular breast cancers reveals heterogeneity in intrinsic molecular subtypes, oncogenic signatures, and mutations

Mixed invasive ductal and lobular carcinoma (MDLC) is a rare histologic subtype of breast cancer displaying both E-cadherin positive ductal and E-cadherin negative lobular morphologies within the same tumor, posing challenges with regard to anticipated clinical management. It remains unclear whether these distinct morphologies also have distinct biology and risk of recurrence. Our spatially-resolved transcriptomic, genomic, and single-cell profiling revealed clinically significant differences between ductal and lobular tumor regions including distinct intrinsic subtype heterogeneity (e.g., MDLC with TNBC/basal ductal and ER+/luminal lobular regions), distinct enrichment of senescence/dormancy and oncogenic (ER and MYC) signatures, genetic and epigenetic CDH1 inactivation in lobular, but not ductal regions, and single-cell ductal and lobular sub-populations with unique oncogenic signatures further highlighting intra-regional heterogeneity. Altogether, we demonstrated that the intra-tumoral morphological/histological heterogeneity within MDLC is underpinned by intrinsic subtype and oncogenic heterogeneity which may result in prognostic uncertainty and therapeutic dilemma.

Significance

MDLC displays both ductal and lobular tumor regions. Our multi-omic profiling approach revealed that these morphologically distinct tumor regions harbor distinct intrinsic subtypes and oncogenic features that may cause prognostic uncertainty and therapeutic dilemma. Thus histopathological/molecular profiling of individual tumor regions may guide clinical decision making and benefit patients with MDLC, particularly in the advanced setting where there is increased reliance on next generation sequencing.

Investigating phenotypic plasticity due to toxicants with exposure disparities in primary human breast cells in vitro

Introduction

Breast cancer is the second most diagnosed cancer, as well as the primary cause of cancer death in women worldwide. Of the different breast cancer subtypes, triple-negative breast cancer (TNBC) is particularly aggressive and is associated with poor prognosis. Black women are two to three times more likely to be diagnosed with TNBCs than white women. Recent experimental evidence suggests that basal-like TNBCs may derive from luminal cells which acquire basal characteristics through phenotypic plasticity, a newly recognized hallmark of cancer. Whether chemical exposures can promote phenotypic plasticity in breast cells is poorly understood.

Methods

To investigate further, we developed a high-content immunocytochemistry assay using normal human breast cells to test whether chemical exposures can impact luminal/basal plasticity by unbiased quantification of keratin 14 (KRT14), a basal-myoepithelial marker; keratin 8 (KRT8), a luminal-epithelial marker; and Hoechst 33342, a DNA marker. Six cell lines established from healthy tissue from donors to the Susan G. Komen Normal Tissue Bank were exposed for 48 hours to three different concentrations (0.1μM, 1μM, and 10μM) of eight ubiquitous chemicals (arsenic, BPA, BPS, cadmium, copper, DDE, lead, and PFNA), with documented exposure disparities in US Black women, in triplicate. Automated fluorescence image quantification was performed using Cell Profiler software, and a random-forest classifier was trained to classify individual cells as KRT8 positive, KRT14 positive, or hybrid (both KRT8 and KRT14 positive) using Cell Profiler Analyst.

Results and discussion

Results demonstrated significant concentration-dependent increases in hybrid populations in response to BPA, BPS, DDE, and PFNA. The increase in hybrid populations expressing both KRT14 and KRT8 is indicative of a phenotypically plastic progenitor-like population in line with known theories of carcinogenesis. Furthermore, BPA, BPS, DDE, and copper produced significant increases in cell proliferation, which could be indicative of a more malignant phenotype. These results further elucidate the relationship between chemical exposure and breast phenotypic plasticity and highlight potential environmental factors that may impact TNBC risk.

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.

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.