Abstract
The adult human breast is comprised of an intricate network of epithelial ducts and lobules that are embedded in connective and adipose tissue1,2,3. Although most previous studies have focused on the breast epithelial system4,5,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.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
$199.00 per year
only $3.90 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
The HBCA website can be viewed at http://www.breastatlas.org. The data are available at the Gene Expression Omnibus (GSE195665). The data are also available from CZI in the CELLxGENE database (https://cellxgene.cziscience.com/collections/4195ab4c-20bd-4cd3-8b3d-65601277e731).
Code availability
The scripts associated with the analysis are available at GitHub (https://github.com/navinlabcode/HumanBreastCellAtlas).
References
Hassiotou, F. & Geddes, D. Anatomy of the human mammary gland: current status of knowledge. Clin. Anat. 26, 29–48 (2013).
Russo, J., Rivera, R. & Russo, I. H. Influence of age and parity on the development of the human breast. Breast Cancer Res. Treat. 23, 211–218 (1992).
Gusterson, B. A. & Stein, T. Human breast development. Semin. Cell Dev. Biol. 23, 567–573 (2012).
Nguyen, Q. H. et al. Profiling human breast epithelial cells using single cell RNA sequencing identifies cell diversity. Nat. Commun. 9, 2028 (2018).
Bach, K. et al. Differentiation dynamics of mammary epithelial cells revealed by single-cell RNA sequencing. Nat. Commun. 8, 2128 (2017).
Bhat-Nakshatri, P. et al. A single-cell atlas of the healthy breast tissues reveals clinically relevant clusters of breast epithelial cells. Cell Rep. Med. 2, 100219 (2021).
Russo, J. & Russo, I. H. Toward a physiological approach to breast cancer prevention. Cancer Epidemiol. Biomarkers Prev. 3, 353–364 (1994).
Crawford, Y. G. et al. Histologically normal human mammary epithelia with silenced p16(INK4a) overexpress COX-2, promoting a premalignant program. Cancer Cell 5, 263–273 (2004).
Wu, S. Z. et al. A single-cell and spatially resolved atlas of human breast cancers. Nat. Genet. 53, 1334–1347 (2021).
Shackleton, M. et al. Generation of a functional mammary gland from a single stem cell. Nature 439, 84–88 (2006).
Visvader, J. E. & Stingl, J. Mammary stem cells and the differentiation hierarchy: current status and perspectives. Genes Dev. 28, 1143–1158 (2014).
Ewald, A. J., Brenot, A., Duong, M., Chan, B. S. & Werb, Z. Collective epithelial migration and cell rearrangements drive mammary branching morphogenesis. Dev. Cell 14, 570–581 (2008).
Zwick, R. K. et al. Adipocyte hypertrophy and lipid dynamics underlie mammary gland remodeling after lactation. Nat. Commun. 9, 3592 (2018).
Pal, B. et al. A single-cell RNA expression atlas of normal, preneoplastic and tumorigenic states in the human breast. EMBO J. 40, e107333 (2021).
Gray, G. K. et al. A human breast atlas integrating single-cell proteomics and transcriptomics. Dev. Cell https://doi.org/10.1016/j.devcel.2022.05.003 (2022).
Twigger, A. J. et al. Transcriptional changes in the mammary gland during lactation revealed by single cell sequencing of cells from human milk. Nat. Commun. 13, 562 (2022).
Azizi, E. et al. Single-cell map of diverse immune phenotypes in the breast tumor microenvironment. Cell 174, 1293–1308 (2018).
Pullan, S. et al. Requirement of basement membrane for the suppression of programmed cell death in mammary epithelium. J. Cell Sci. 109, 631–642 (1996).
Polyak, K. & Kalluri, R. The role of the microenvironment in mammary gland development and cancer. Cold Spring Harb. Perspect. Biol. 2, a003244 (2010).
Hanasoge Somasundara, A. V. et al. Parity-induced changes to mammary epithelial cells control NKT cell expansion and mammary oncogenesis. Cell Rep. 37, 110099 (2021).
Rozenblatt-Rosen, O., Stubbington, M. J. T., Regev, A. & Teichmann, S. A. The Human Cell Atlas: from vision to reality. Nature 550, 451–453 (2017).
Hovey, R. C. & Aimo, L. Diverse and active roles for adipocytes during mammary gland growth and function. J. Mammary Gland Biol. Neoplasia 15, 279–290 (2010).
Emont, M. P. et al. A single-cell atlas of human and mouse white adipose tissue. Nature 603, 926–933 (2022).
Aibar, S. et al. SCENIC: single-cell regulatory network inference and clustering. Nat. Methods 14, 1083–1086 (2017).
Pellacani, D., Tan, S., Lefort, S. & Eaves, C. J. Transcriptional regulation of normal human mammary cell heterogeneity and its perturbation in breast cancer. EMBO J. 38, e100330 (2019).
Stahl, P. L. et al. Visualization and analysis of gene expression in tissue sections by spatial transcriptomics. Science 353, 78–82 (2016).
Goltsev, Y. et al. Deep profiling of mouse splenic architecture with CODEX multiplexed imaging. Cell 174, 968–981 (2018).
Menendez, J. A. & Lupu, R. Fatty acid synthase regulates estrogen receptor-α signaling in breast cancer cells. Oncogenesis 6, e299 (2017).
Riese, D. J. 2nd & Cullum, R. L. Epiregulin: roles in normal physiology and cancer. Semin. Cell Dev. Biol. 28, 49–56 (2014).
Kannan, N. et al. The luminal progenitor compartment of the normal human mammary gland constitutes a unique site of telomere dysfunction. Stem Cell Rep. 1, 28–37 (2013).
Zheng, L. et al. Pan-cancer single-cell landscape of tumor-infiltrating T cells. Science 374, abe6474 (2021).
Cheng, S. et al. A pan-cancer single-cell transcriptional atlas of tumor infiltrating myeloid cells. Cell 184, 792–809 (2021).
Hu, Q. et al. Atlas of breast cancer infiltrated B-lymphocytes revealed by paired single-cell RNA-sequencing and antigen receptor profiling. Nat. Commun. 12, 2186 (2021).
Wherry, E. J. & Kurachi, M. Molecular and cellular insights into T cell exhaustion. Nat. Rev. Immunol. 15, 486–499 (2015).
Milner, J. J. et al. Runx3 programs CD8+ T cell residency in non-lymphoid tissues and tumours. Nature 552, 253–257 (2017).
Colleluori, G., Perugini, J., Barbatelli, G. & Cinti, S. Mammary gland adipocytes in lactation cycle, obesity and breast cancer. Rev. Endocr. Metab. Disord. 22, 241–255 (2021).
Kalucka, J. et al. Single-cell transcriptome atlas of murine endothelial cells. Cell 180, 764–779 (2020).
Schupp, J. C. et al. Integrated single-cell atlas of endothelial cells of the human lung. Circulation 144, 286–302 (2021).
Takeda, A. et al. Single-cell survey of human lymphatics unveils marked endothelial cell heterogeneity and mechanisms of homing for neutrophils. Immunity 51, 561–572 e565 (2019).
Sweeney, M. D., Ayyadurai, S. & Zlokovic, B. V. Pericytes of the neurovascular unit: key functions and signaling pathways. Nat. Neurosci. 19, 771–783 (2016).
Armulik, A., Genove, G. & Betsholtz, C. Pericytes: developmental, physiological, and pathological perspectives, problems, and promises. Dev. Cell 21, 193–215 (2011).
Brozovich, F. V. et al. Mechanisms of vascular smooth muscle contraction and the basis for pharmacologic treatment of smooth muscle disorders. Pharmacol. Rev. 68, 476–532 (2016).
Visvader, J. E. Keeping abreast of the mammary epithelial hierarchy and breast tumorigenesis. Genes Dev. 23, 2563–2577 (2009).
Schmid, P. et al. Pembrolizumab for early triple-negative breast cancer. N. Engl. J. Med. 382, 810–821 (2020).
Choi, H. Y. et al. Preoperative axillary lymph node evaluation in breast cancer: current issues and literature review. Ultrasound Q. 33, 6–14 (2017).
Kothari, C., Diorio, C. & Durocher, F. The importance of breast adipose tissue in breast cancer. Int. J. Mol. Sci. https://doi.org/10.3390/ijms21165760 (2020).
Garbe, J. C. et al. Accumulation of multipotent progenitors with a basal differentiation bias during aging of human mammary epithelia. Cancer Res. 72, 3687–3701 (2012).
Pelissier Vatter, F. A. et al. High-dimensional phenotyping identifies age-emergent cells in human mammary epithelia. Cell Rep. 23, 1205–1219 (2018).
Fraser, I. S., Critchley, H. O., Munro, M. G. & Broder, M. Can we achieve international agreement on terminologies and definitions used to describe abnormalities of menstrual bleeding? Hum. Reprod. 22, 635–643 (2007).
Murrow, L. M. et al. Mapping hormone-regulated cell-cell interaction networks in the human breast at single-cell resolution. Cell Syst. 13, 644–664 (2022).
Navin, N. et al. Tumour evolution inferred by single-cell sequencing. Nature 472, 90–94 (2011).
Baslan, T. et al. Genome-wide copy number analysis of single cells. Nat. Protoc. 7, 1024–1041 (2012).
Lander, E. S. et al. Initial sequencing and analysis of the human genome. Nature 409, 860–921 (2001).
Schneider, V. A. et al. Evaluation of GRCh38 and de novo haploid genome assemblies demonstrates the enduring quality of the reference assembly. Genome Res. 27, 849–864 (2017).
Stuart, T. et al. Comprehensive integration of single-cell data. Cell 177, 1888–1902 (2019).
Aran, D. et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat. Immunol. 20, 163–172 (2019).
Tirosh, I. et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science 352, 189–196 (2016).
Sergushichev, A. A. An algorithm for fast preranked gene set enrichment analysis using cumulative statistic calculation. Preprint at bioRxiv https://doi.org/10.1101/060012 (2016).
Yu, G., Wang, L. G., Han, Y. & He, Q. Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–287 (2012).
Bankhead, P. et al. QuPath: open source software for digital pathology image analysis. Sci Rep. 7, 16878 (2017).
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012).
A, L. & M, W. Classification and regression by randomForest. R News 3, 18–22 (2002).
Wei, R. et al. Spatial charting of single-cell transcriptomes in tissues. Nat. Biotechnol. https://doi.org/10.1038/s41587-022-01233-1 (2022).
Stringer, C., Wang, T., Michaelos, M. & Pachitariu, M. Cellpose: a generalist algorithm for cellular segmentation. Nat. Methods 18, 100–106 (2021).
Satija, R., Farrell, J. A., Gennert, D., Schier, A. F. & Regev, A. Spatial reconstruction of single-cell gene expression data. Nat. Biotechnol. 33, 495–502 (2015).
Efremova, M., Vento-Tormo, M., Teichmann, S. A. & Vento-Tormo, R. CellPhoneDB: inferring cell–cell communication from combined expression of multi-subunit ligand–receptor complexes. Nat. Protoc. 15, 1484–1506 (2020).
Acknowledgements
This work was supported by major funding from the Chan-Zuckerberg Initiative (CZI) SEED Network Grant (CZF2019-002432); and grants to N. Navin from the NIH National Cancer Institute (RO1CA240526, RO1CA236864, 1R01CA234496, F30CA243419), the CPRIT Single Cell Genomics Center (RP180684), the American Cancer Society (132551-RSG-18-194-01-DDC). N. Navin is an AAAS Fellow, AAAS Wachtel Scholar, Damon-Runyon Rachleff Innovator, Andrew Sabin Fellow and Jack & Beverly Randall Innovator. T.K. is funded by the NCI T32 Translational Genomics Fellowship and Rosalie B. Hite fellowship. This study was supported by the MD Anderson Sequencing Core Facility Grant (CA016672). M.P. is supported by a fellowship from the CIRM Training Grant (EDUC4-12822). The work was also supported by a Damon-Runyon Quantitative Biology Postdoctoral Fellow to R.W. We thank B. Marshall, N. Tavares and J. Cool for their guidance and support; J. Waters, S. Stingley and L. Vann from for their support on this project; J. Wiley, A. Wood, A. Alexander and A. Contreras for clinical support; A. Longworth, S. Mallya and M. Curran for their advice; Y. Lin and R. Ye at MD Anderson for help with experiments; J. Zamanian and J. Yu-Sheng Chien for assistance with data depositing; and all of the women who participated in the HBCA and donated their breast tissue to this project. We thank Enable Medicine for their help with CODEX data generation. This publication is part of the HCA (www.humancellatlas.org).
Author information
Authors and Affiliations
Contributions
scRNA-seq experiments were performed by K.N., Q.H.N., S.B., T.K., E.S., J.L., M.P., T.T. and S.M. Spatial genomics experiments were performed by S.B., E.S., B.N., N.K., O.B., B.B.C. and N. Nikulina. Single-cell data analysis was performed by T.K., R.W., S. He., K.B., M.P., Y.G., M.H., A.K.C., B.N., N.K., K.C. and T.L. Spatial data analysis was performed by R.W., S. He and A. Thennavan. Tissue samples and clinical coordination was performed by O.B., B.B.C., H.C., A.K.C., M.T., B.M., H.J., J.M., R.E., D.F.T., C.N., E.L., R.P., S.W., S.M., A. Thompson, B.L. and S. Hanson. Tissue pathological analysis was performed by H.C., C.N. and A. Thennavan. Project management and manuscript writing was performed by B.L., D.A.L., N. Navin and K.K. N. Navin and K.K. are the coordinators for the Breast Atlas Bionetwork that is part of the HCA Project.
Corresponding authors
Ethics declarations
Competing interests
N. Navin previously served on the scientific advisory board for Resolve Biosciences (2020–2022) but did not receive any compensation. O.B., B.B.C. and N. Nikulina are employees at Akoya Biosciences. B.N. and N.K. are employees of Resolve Biosciences. The other authors declare no competing interests.
Peer review
Peer review information
Nature thanks Itai Yanai and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Frequency of Major Breast Cell Types Across Women and Sample Types.
a, Experimental workflow for breast tissue processing for scRNA-seq, showing different conditions used for digestion times and trypsin treatments. b, Pie chart showing ethnic backgrounds of women who provided tissue samples for the breast atlas. c, Cell type frequencies for different tissue sources (reduction mammoplasties - RM, prophylactic mastectomies - PM and contralateral mastectomies - CM), cells vs nuclei single cell RNA-seq protocols, and experimental dissociation protocol. d, Major cell type frequencies of matched left and right breasts from 22 women (left) and averages across all left and all right breast tissues (right). e, Stacked barplot showing the variation of cell type frequencies across the 126 women in scRNA-seq data. The top annotation bar shows the experimental workflow used (short/medium/long). f, Top regulons identified with SCENIC for each cell type cluster from the snRNA-seq data. g, Top regulons identified with SCENIC for each cell type cluster from the scRNA-seq data. h, Multi-dimensional scaling and Procrustes analysis to determine the concordance of left and right breast cell type frequencies and Pearson correlations for 22 women with matched breast tissue samples. P-value was calculated based on a two-sided test.
Extended Data Fig. 2 Ligand-receptor Interaction Analysis Between Cell Types.
Ligand-receptor interaction plots predicted from scRNA-seq data using CellPhoneDB between the major breast cell types. a, Interaction plot between the epithelial (Basal, LumHR and LumSec) cell types. b, Interaction plot between the epithelial and immune (B-cells, T-cells, Myeloid cells) cell types. c, Interaction plot between the epithelial and stromal (Fibroblasts, Perivascular and Vascular endothelial cells) cell types. d, Interaction plot within the stromal (Fibroblasts, Myeloid, Lymphatic, Vascular and Perivascular cells) cell types. e, Interaction plot between adipocytes and stromal cell types.
Extended Data Fig. 3 Spatial Transcriptomic Analysis of Breast Cell Types.
a, Integrated UMAP and unbiased clustering of ST data from 10 breast samples, showing 9 ST clusters. b, Histopathological images, and spatial distribution of ST clusters in the ST data from the breast tissues. c, Concordance of ST clusters and the scRNA-seq clusters of the major cell types using Fisher’s exact test. d, Pearson correlation analysis of marker gene expression levels between the ST clusters and the scRNA-seq data for different cell types. All p-values were calculated based on two-sided tests.
Extended Data Fig. 4 Spatial Analysis of Breast Cell Types with CODEX and smFISH.
a, Cell segmentation results of smFISH (Resolve) data across 12 tissue samples profiled from 5 different women. Cells were annotated based on combinations of markers for each cell type as described in Supplementary Table 6. b, Densities of cell types across three topographic areas using 12 tissues profiled by smFISH (Resolve). c, Heatmap of the top 5 targeted maker genes for each cell type in the smFISH (Resolve) data from 12 combined tissue samples. d, Cell segmentation results of CODEX data from 8 different women. Cells were annotated based on combinations or single protein markers to identify different cell types. e, Densities of cell types across three topographic areas from 8 different women by CODEX. f, Heatmap showing protein levels for markers that were used to identify different cell types in the CODEX data. (D: ducts, L: lobules and C: connective regions).
Extended Data Fig. 5 Analysis of Single Cell and Spatial Epithelial Data.
a, UMAPs of snRNA-seq data showing the expression of hormone receptor genes. b, Epithelial cell state frequencies across the 126 women in scRNA-seq data, where the top annotation bar represents the dissociation protocol. c, UMAP feature plots showing the expression of previously reported stem cell marker genes in the scRNA-seq epithelial dataset. d, Ligand-receptor interactions within the epithelial cell states predicted with CellPhoneDB. e, Cell cycle scoring of S-phase for different epithelial cell states detected in the scRNA-seq data. f, Cell cycle scoring for S-phase in the epithelial cell type clusters detected in the snRNA-seq data. g, smFISH (Resolve) data showing the expression of the MKI67 proliferation marker in the epithelial cells of the ducts and lobules from 4 different breast tissues. h, UMAP of different LumSec cell states and ELF5, LTF signature scores, respectively. i, Histopathological image of adjacent H&E section showing the anatomic annotation of ducts and lobules (left) and smFISH MERFISH (right panel) from P101 showing the spatial distribution of different LumSec cell states across different regions. j, Stacked barplot showing the distribution of different LumSec cell proportions in ducts and lobules across 3 MERFISH samples. k, Histopathological image (left panel) and smFISH MERFISH (right panel) from P101 showing the spatial distribution of the LumHR-SCGB population in a specific region of epithelium.
Extended Data Fig. 6 Spatial analysis of epithelial cells in ductal and lobular structures.
a, Spatial transcriptomic analysis showing clusters labelled as duct or lobule/TDLU from 3 breast tissues (P10, P35 and P47). b, smFISH (Resolve) data (P46-S1 and P46-S4) showing a subset of Keratin markers (left) and hormone receptor genes (right) and their localization to different breast tissue regions annotated as either duct or lobule/TDLU. c, CODEX data from P131 showing KRT5 in ducts and KRT19 in lobules/TDLU regions, with enlarged panels of the right. d, CODEX analysis from P130 of ductal and lobular/TDLU regions, showing differences for KRT14 levels in ducts and lobules. e, CODEX data from P131 showing protein levels of KRT8 and progesterone receptor (PR) in epithelial cells in the ducts and lobular/TDLU regions.
Extended Data Fig. 7 Immune cell subtypes in the breast and their variation in women.
a, H&E staining of plasma B-cells, T-cells, mast cell and macrophages (arrows) in human breast tissues. b, Stacked barplot showing the cell type frequencies of T, B and myeloid cells across 126 women in scRNA-seq data. Top annotation bar represents different tissue dissociation protocols that were utilized. c-e, Stacked barplots showing the cell state frequencies of T, B and myeloid cells across 126 women in scRNA-seq data respectively. f, Dot plot showing expression of checkpoint/exhaustion markers in NK and T cell states from the scRNA-seq data of 126 women. g, Ligand-receptor interaction analysis predicted with CellPhoneDB between the fibroblasts cell states and macrophage cell states.
Extended Data Fig. 8 Spatial analysis of immune cells in human breast tissues.
a, CODEX data from patient P130 and P131 showing localization of different immune cells with epithelial marker KRT19 and vascular marker CD31. Yellow and white arrows indicate CD4 Tregs and DCs, respectively. b, Frequency of T cells with the RUNX3 tissue residency marker in CODEX data. c, CODEX data (P130) showing immune cells in ductal, lobular and connective regions. d, Stacked barplots of CODEX data showing the density of immune cell types in each spatial region in 8 women. e, smFISH (Resolve) data (P46-S1) showing RNA localization of T, B and myeloid cells. f, Segmented smFISH (Resolve) data (P46-S1) showing cell localization of T, B and myeloid cells. g, smFISH (Resolve) data (P46-S1 and P47-S1) showing immune cell localization of B, T and myeloid cells across ducts, lobules and connective regions. h, Stacked barplots of smFISH (Resolve) data showing the density and proportion of immune cell types in different spatial regions. i, Adjacent histopathological tissue section (left) and segmented smFISH MERFISH data (right) from patient P91 showing the spatial distribution of m1, m2 macrophages and cDC2 populations in different regions of human breast tissue. j, Stacked barplot showing the density of m1, m2 macrophages and cDC2 populations in different regions across 3 smFISH MERFISH samples. k, Adjacent histopathological tissue section (left) and segmented smFISH MERFISH data (right) from patient P96 showing the spatial distribution of different B-cell states in different regions of human breast tissue. l, Stacked barplot showing the density of different B-cell states in different regions across three smFISH MERFISH samples.
Extended Data Fig. 9 Fibroblast cell states in the human breast.
a, Stacked barplot showing the fibroblasts cell state frequencies across 126 women in scRNA-seq data with top annotation bar representing the tissue dissociation protocol. b, Gene ontology enrichment analysis showing top enriched biological process gene sets associated with each cell state (Pos: positive; Neg: negative; Reg: regulation; RSTK: receptor protein serine/threonine kinase; TGF: transforming growth factor; IGF: insulin-like growth factor). c, CODEX data from P132 showing fibroblasts marked by VIM in the connective tissue (I) and interlobular (II) regions. d, smFISH (Resolve) data showing fibroblast markers in areas of connective tissue regions (I) and epithelial regions (II) from two women (P47-S1 and P69-S3). e, smFISH (Resolve) data (P35-S1) indicating spatial proximity regions with epithelial-proximal (Epi-prox), epithelial-middle (Epi-mid) and epithelial-distant (Epi-Dist) regions for 4 marker genes. f, Percentages of 4 markers that are proximal, middle or distant to the epithelial cells, quantified from the smFISH (Resolve) data. g, RNAscope in situ hybridization of breast tissues using an MMP3 probe in combination with anti-Vimentin and anti-PanCK immunofluorescent staining, with enlarged panel (right). h, Ligand-receptor interactions between fibroblasts, adipocytes and myeloid cell states predicted using CellPhoneDB.
Extended Data Fig. 10 Endothelial cell diversity in the Human Breast.
a, Stacked barplot showing the endothelial cell state frequencies across the 126 women in scRNA-seq data, with top annotation bar showing the tissue dissociation protocol. b, Dot plot of gene ontology enrichment results for 4 lymphatic cell states. c, Heatmap showing top gene expression for vascular and lymphatic endothelial clusters detected in the ST data. d, smFISH (Resolve) data showing veins (ACKR1) and capillaries (RBP7), as well as a canonical vascular marker (VWF) in two different HBCA samples (P46-S3 and P69-S3). e, Adjacent H&E tissue section with pathological annotations (left panel) and segmented smFISH MERFISH data (right panel) from P101 showing the spatial distribution of vascular endothelial cell states in different regions of human breast tissue. f, Stacked barplot showing the density of vascular endothelial states in different regions across 3 smFISH MERFISH samples.
Extended Data Fig. 11 Perivascular cells in Human Breast Tissues.
a, Stacked barplot showing the perivascular cell state frequencies across the 126 women in scRNA-seq data with top annotation bars indicating the tissue dissociation protocol. b, UMAPs of pericytes and vascular smooth muscle cells (VSMCs) and feature plots of the VSMCs marker genes (SYNM and ACTG2). c-e, smFISH (Resolve) data showing expression of pericyte marker RGS5, together with vascular marker VWF and fibroblast marker COL1A1 in lobular and ductal regions from 2 different breast tissue samples (P47-S1 and P46-S3). f, CODEX results from P131 showing vascular cells (anti-CD31) and pericytes (anti-LIF) in a TDLU region. g, smFISH MERFISH from P96 showing the spatial distribution of vascular endothelial cell states (left panel) and perivascular cell states (right panel) in different regions of human breast tissue. h, smFISH (MERFISH) data showing arteries (SOX17) and VSMCs (ATCG2 and SYNM) in breast tissue.
Extended Data Fig. 12 Metadata correlations with breast cell types and states.
a, Boxplots showing the major cell type frequencies across ethnicity status in the n = 69 women using Wilcoxon rank sum test (top). Significant associations of cell states with ethnicity status using Fisher’s exact test (bottom). b, Boxplots showing the major cell type frequencies across pre- and post-menopause status in the n = 71 women using Wilcoxon rank sum test (top). Significant associations of cell states with menopause status using Fisher’s exact test (bottom). c, Boxplots showing the major cell type frequencies across different age groups using Wilcoxon rank sum test, young (<50 years) and old (>50 years) for n = 76 women (top). Significant associations of cell states with age groups using Fisher’s exact test (bottom). d, Boxplots showing the major cell type frequencies across different breast density (high, low) groups in the n = 16 women using Wilcoxon rank sum test (top). Significant associations of cell states with breast density using Fisher’s exact test (bottom). e, Boxplots showing the major cell type frequencies across different BMI status in 73 women using Wilcoxon rank sum test, overweight (BMI >= 25 and < 30) and obese (BMI >= 30). f, Boxplots showing the major cell type frequencies across different parity status (nulliparous, parous) status in the n=64 women using Wilcoxon rank sum test. All p-values were calculated based on two-sided tests. Boxplots show the median with interquartile ranges (25–75%), while whiskers extend to 1.5× the interquartile range from the box.
Extended Data Fig. 13 Summary of the Major Cell Types and States in Breast Tissues.
This illustration summarizes all of the breast cell types and cell states that were identified in the HBCA study. a, Summary of cell lineages from cell types to cell states. b, Mapping of cell types and cell states to the four major spatial regions (Adipose, Connective, Ductal, Lobular) that were supported by the spatial technologies. Not all cell states were assigned to specific spatial regions, in cases where the data did not support their assignment. Individual figures were created with BioRender.com.
Supplementary information
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Kumar, T., Nee, K., Wei, R. et al. A spatially resolved single-cell genomic atlas of the adult human breast. Nature 620, 181–191 (2023). https://doi.org/10.1038/s41586-023-06252-9
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-023-06252-9
This article is cited by
-
A single-cell atlas enables mapping of homeostatic cellular shifts in the adult human breast
Nature Genetics (2024)
-
White adipocyte dysfunction and obesity-associated pathologies in humans
Nature Reviews Molecular Cell Biology (2024)
-
A cell transcriptomic profile provides insights into adipocytes of porcine mammary gland across development
Journal of Animal Science and Biotechnology (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
